Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TWO-DIMENSIONAL SYSTEMS”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analytical Model for Atomic Relaxation in Twisted Moiré Materials

By virtue of being atomically thin, the electronic properties of heterostructures built from two-dimensional materials are strongly influenced by atomic relaxation. The atomic layers behave as flexible membranes rather than rigid crystals. Here we develop an analytical theory of lattice relaxation in twisted moiré materials. We obtain analytical results for the lattice displacements and corresponding pseudo gauge fields, as a function of twist angle. We benchmark our results for twisted bilayer graphene and twisted WSe 2 bilayers using large-scale molecular dynamics simulations. Our single-parameter theory is valid in graphene bilayers for twist angles 𝜃 ≳ 0.7°, and in twisted WSe 2 for 𝜃 ≳ 1.6°. Furthermore, we also investigate how relaxation alters the electronic structure in twisted bilayer graphene, providing a simple extension to the continuum model to account for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

In Situ Imaging Reveals Efficient Charge Separation in Monolayer MoS 2 –WS 2 Type-II Heterojunctions

Covalently bonded in-plane two-dimensional (2D) transition metal dichalcogenide (TMD) heterojunctions with atomically sharp interfaces hold great promise for photocatalytic applications in solar energy conversion and environmental remediation; however, their spatially resolved charge distribution and transport, particularly under operando conditions, remain poorly understood. Here, we employ photoscanning electrochemical microscopy (photo-SECM) to directly visualize photoinduced charge separation in monolayer MoS 2 –WS 2 in-plane heterojunctions. Spatial separation of photogenerated carriers is observed, with electrons accumulating in MoS 2 and holes in WS 2 , leading to strongly asymmetric interfacial kinetics: Fc + reduction proceeds rapidly on MoS 2 (0.6 cm s –1 ), whereas Fc oxidation on WS 2 is significantly slower (0.008 cm s –1 ). High-resolution surface photovoltage microscopy (SPVM) enables a quantitative comparison of charge-separation capacity across architectures. The in-plane MoS 2 –WS 2 heterojunction shows the largest photovoltage contrast (−35 mV in MoS 2 , 20 mV in WS 2 ), exceeding the vertical heterojunction (−18 mV in MoS 2 , 11 mV in WS 2 ) and the individual monolayers (−12 mV for MoS 2 , – 1 mV for WS 2 ), establishing the following trend: in-plane > vertical > monolayers. Ultraviolet photoelectron spectroscopy (UPS) indicates that this directional charge separation is driven by intrinsic type-II band alignment, while photoluminescence (PL) imaging shows that the interface acts as a recombination center that limits efficient carrier extraction. These results provide direct experimental evidence of type-II-driven charge separation in in-plane heterojunctions and offer critical insights for interface design in high-efficiency photocatalytic and optoelectronic systems.

electrical properties↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optical control of integer and fractional Chern insulators

Optical control of topology, particularly in the presence of electron correlations, is an interesting topic with broad scientific and technological impact. Twisted MoTe 2 bilayer (tMoTe 2 ) is a zero-field fractional Chern insulator (FCI), exhibiting the fractionally quantized anomalous Hall effect. As the chirality of the edge states and sign of the Chern number are determined by the underlying ferromagnetic polarization, manipulation of ferromagnetism would realize control of the Chern insulator (CI)/FCI states. Here, in this work, we demonstrate control of ferromagnetic polarization, and thus the CI and FCI states, by circularly polarized optical pumping in tMoTe 2 . At low excitation power, we achieve on-demand preparation of ferromagnetic polarization by optical training, that is, electrically tuning the system from non-ferromagnetic to desirable ferromagnetic states under helicity-selective optical pumping. With increased excitation power, we further realize direct optical switching of ferromagnetic polarization at a temperature far below the Curie temperature. Both optical training and direct switching are most effective near CI and FCI states, which we attribute to a gap-enhanced valley polarization of optically pumped holes. The magnetization can be dynamically switched by modulating the helicity of optical excitation. Spatially resolved measurements further demonstrate optical writing of ferromagnetic, and thus CI (or FCI) domains. Our work realizes precise optical control of a topological quantum many-body system with potential applications in topological spintronics, quantum memories and creation of exotic edge states by programmable patterning of integer and fractionally quantized anomalous Hall domains.

ferromagnetism↗
DOE OSTI2026

Electronic origin of reorganization energy in interfacial electron transfer

Electron transfer (ET) reactions underpin energy conversion and chemical transformations in both biological and a biological systems. The efficiency of any ET process relies on achieving a desired ET rate within an optimal driving force range. Marcus theory provides a microscopic framework for understanding the activation free energy—and therefore the rate—of ET in terms of a key parameter: the reorganization energy. For electrified solid–liquid interfaces, it has long been conventionally understood that only factors in the electrolyte phase are responsible for determining the reorganization energy and that the electronic density of states (DOS) of the electrode only serves to dictate the number of thermally accessible channels for ET. Here we show instead that the electrode DOS plays a central role in governing the reorganization energy, far outweighing its conventionally assumed role. Using atomically layered heterostructures, we tune the DOS of graphene and measure outer-sphere ET kinetics. We find the ensuing variation in ET rate arises from strong modulation in a reorganization energy associated with image potential localization in the electrode. Here we redefine the traditional paradigm of heterogeneous ET kinetics, revealing a deeper role of the electrode electronic structure in interfacial reactivity.

Electrochemistry↗
DOE OSTI2026

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2026

Ion Transport in Charged Membranes: Linking Electric-Field-Driven Mechanisms to Pore Size via Perturbation Analysis

Ion-exchange membranes are a critical component in electrochemical systems. Nevertheless, the understanding and modeling of ion transport within these porous structures have been limited by particular complexity reductions, either ignoring the dimensionality of their porous network architecture or imposing geometric assumptions (i.e., overlapping double layers). Before addressing this morphology-transport gap, a framework that relates the driving forces of transport to the geometry of a single pore is required. In this work, our modeling domain consists of a two-dimensional single pore with charged walls, connecting two identical electrolyte reservoirs. Using the Poisson-Nernst-Planck equations and regular perturbation theory, we decouple the electric fields and analyze the driving forces of ion transport, specifically electromigration and induced electroosmosis within the pore. These processes are described as analytical functions of the interaction aspect ratio, ?, defined as the ratio of the pore radius to the Debye length. Using this parameter, our study (i) describes the interplay between electromigrative and electroosmotic mechanisms that set ionic conductivity, (ii) identifies a dimensionless group of intrinsic electrolyte properties that indicates the predominant driving force, and (iii) provides a qualitative, confinement-dependent perspective on selectivity in ion-conducting membranes.

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

Pseudocapacitive Titanium Oxynitride Nanowires for Ultrahigh Capacitance Supercapacitors

High-quality, multifunctional two-dimensional (2D) titanium oxynitide (TiNO) thin films and one-dimensional (1D) TiNO nanowires have been synthesized using a pulsed laser deposition, a simple, fast, and congruent evaporation method. First-principles calculations as a function of surface orientation and termination indicate that surface oxidation of TiNO nanowires can stabilize the (110) orientation observed experimentally. The specific capacitance value for the TiNO nanowire samples (2725 mF/cm 2 ) has been found to be nearly six times more than that of the TiNO thin film samples (400 mF/cm 2 ), which is attributed to the high packing density of TiNO nanowires over a given area. The nanowire samples have also been found to exhibit a significantly higher energy density (1.35 μWh/cm 2 ) than the TiNO thin-film samples (0.33 μWh/cm 2 ). Thus, the TiNO material system in thin-film and nanowire forms has been demonstrated to be a promising candidate for use as an electrode material in supercapacitors and other charge-storage applications.

capacitors↗
DOE OSTI2026

Imaging a terahertz superfluid plasmon in a two-dimensional superconductor

The superconducting gap defines the fundamental energy scale for the emergence of dissipationless transport and collective phenomena in a superconductor. In layered high-temperature cuprate superconductors, in which the Cooper pairs are confined to weakly coupled two-dimensional (2D) copper–oxygen (CuO 2 ) planes, terahertz (THz) spectroscopy at subgap millielectronvolt (meV) energies has provided crucial insights into the collective superfluid response perpendicular to the superconducting layers. However, within the CuO 2 planes, the collective superfluid response manifests as plasmonic charge oscillations at energies far exceeding the superconducting gap, obscured by strong dissipation. Here, in this study, we present spectroscopic evidence of a below-gap, 2D superfluid plasmon in few-layer Bi 2 Sr 2 CaCu 2 O 8+x and spatially resolve its deeply subdiffractive THz electrodynamics. By placing the superconductor in the near field of a spintronic THz emitter, we reveal this distinct resonance—absent in bulk samples and observed only in the superconducting phase—and determine its plasmonic nature by mapping the geometric anisotropy and dispersion. Crucially, these measurements offer a direct view of the momentum-dependent and frequency-dependent superconducting transition in two dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2026

MATEY: multiscale adaptive transformer models for spatiotemporal physical systems

Accurate representation of the multiscale features in spatiotemporal physical systems using vision transformer architectures requires extremely long, computationally prohibitive token sequences. To address this issue, we propose two novel adaptive tokenization schemes that dynamically adjust patch sizes based on local features: one ensures convergent behavior to uniform patch refinement, while the other offers better computational efficiency. Moreover, we present a set of spatiotemporal attention schemes, where the temporal or axial spatial dimensions are decoupled, to evaluate their baseline computational and data efficiencies and to determine whether adaptive tokenization can improve this performance. We assess the performance of the proposed multiscale adaptive model, MATEY, in a sequence of experiments. Compared to a full spatiotemporal attention scheme or a scheme that decouples only the temporal dimension, we find that fully decoupled axial attention is less efficient and expressive, requiring more training time and model parameters to achieve the same accuracy. The experiments on the adaptive tokenization schemes show that, compared to a uniformly refined model, the proposed schemes achieve comparable or improved accuracy at a much lower cost in the tested two-dimensional settings. While the asymptotic analysis suggests the potential for favorable scaling, empirical validation at substantially longer sequence lengths remains to be performed in future work. Finally, we demonstrate in two fine-tuning tasks featuring different physics that models pretrained on PDEBench data outperform the ones trained from scratch, especially in the low data regime with frozen attention.

adaptive tokenization↗
DOE OSTI2026

Josephson dynamics in two-dimensional ring-shaped condensates

We investigate Josephson transport in a fully closed, two-dimensional superfluid circuit formed by a ring-shaped 87 Rb Bose-Einstein condensate that contains two optical barriers acting as movable weak links. Translating these barriers at controlled speeds imposes a steady bias current, enabling direct mapping of the current-chemical-potential (𝐼−Δ⁢𝜇) characteristics. For narrow junctions (𝑤 ≈ 1µ⁢m), the circuit exhibits a pronounced dc branch that terminates at a critical current 𝐼 𝑐 = 9⁢(1) × 10 3 s −1 ; above this threshold, the system switches to an ac, resistive regime. Classical-field simulations that include the moving barriers quantitatively reproduce both the nonlinear 𝐼−Δ⁢𝜇 curve and the measured 𝐼 𝑐 , validating the underlying microscopic picture. Analysis of the ensuing phase dynamics shows that dissipation is mediated by the nucleation and traversal of vortex-antivortex pairs through the junctions, while the bulk condensate remains globally phase locked—direct evidence of the ring's topological constraint enforcing quantized circulation. These results establish a cold-atom analog of a superconducting quantum interference device in which Josephson dynamics can be resolved at the single-vortex level, providing a versatile platform for atomtronic circuit elements, nonreciprocal Josephson devices, and on-chip Sagnac interferometers for multiaxis rotation sensing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Diagonal isometric form for tensor network states in two dimensions

Isometric tensor network states (isoTNS) generalize the isometric form of the one-dimensional matrix product states to tensor networks in two and higher dimensions. Here, we introduce an alternative isometric form for isoTNS by incorporating auxiliary tensors to represent the orthogonality hypersurface. We implement the time evolving block decimation algorithm on this new isometric form and benchmark the method by computing ground states and the real time evolution of the transverse field Ising model in two dimensions on large square lattices of up to 1250 sites. Our results demonstrate that isoTNS can efficiently capture the entanglement structure of two-dimensional area law states. The short-time dynamics is also accurately reproduced even at the critical point. Our isoTNS formulation further allows for a natural extension to different lattice geometries, such as the honeycomb or kagome lattice.

2-dimensional systems↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Dynamics of ballistic photocurrents driven by Coulomb scattering in a two-dimensional material

First-principles real-time time-dependent density-functional theory (rt-TDDFT) calculations reveal the existence of ballistic photocurrents generated by Coulomb scattering, which has not previously been considered as a mechanism for the bulk photovoltaic effect. With monolayer GeS as an example, it is predicted that ballistic currents can be comparable to shift currents under experimentally accessible conditions.

2-dimensional systems↗
DOE OSTI2026

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2026

Molecular Inventory and Comparative Organic Profiling of A Homogenized Aggregate Sample From Asteroid (101955) Bennu

Evaluating the distribution and diversity of solvent-soluble organic molecules in extraterrestrial materials provide information about the physicochemical environments where prebiotic chemistry occurred and is essential for assessing the inventory of prebiotic compounds available to the early Earth. In this work, we conducted a comprehensive molecular characterization of solvent-soluble organic compounds in a homogenized aggregate sample (ID: OREX-800107-128) of unsorted regolith from asteroid (101955) Bennu returned by NASA’s OSIRIS-REx mission. Using two-dimensional gas chromatography coupled with high-resolution time-of-flight mass spectrometry (GC×GC-HRMS), we identified 87 organic compounds across several molecular classes, including alkanes, polycyclic aromatic hydrocarbons, and nitrogen-, oxygen-, and sulfur-bearing species. Comparison with the organic inventory of CI (Orgueil), C2-ungrouped (Tagish Lake and Tarda), and CM2 (Murchison) carbonaceous chondrites revealed that Bennu sample OREX-800107-128 is most similar to Tagish Lake and Tarda yet compositionally distinct from all of these meteorites. This Bennu sample exhibits a relatively high abundance of low-molecular-weight aromatics (e.g., toluene), pyridine derivatives, and sulfur species such as dimethyl sulfite, consistent with episodic low-temperature aqueous alteration. The presence of these organic compounds suggests that this material experienced chemically diverse processes in the parent body conducive to complex organic synthesis, including aldol condensation reactions responsible for N- and O-bearing compounds. We also compared the molecular distributions of this Bennu sample and previously analyzed samples from (162173) Ryugu to assess potential common origins between these two carbonaceous asteroids. Bennu shows more chemically diverse and nitrogen-rich organic materials than Ryugu. While both asteroids share broadly similar primitive compositions, the observed chemical differences indicate distinct evolutionary histories within a likely common formation environment though whether Bennu and Ryugu derive from a single compositionally heterogeneous parent body or from multiple closely related progenitors remains unresolved. Further detailed comparative analyses of material from Bennu, Ryugu, and specific carbonaceous meteorites (e.g., Tagish Lake and Tarda) will provide new insights into how the N-rich molecular profile observed in Bennu samples was formed and evolved in the outer Solar System.

sample return↗