Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HEXAGONAL CONFIGURATION”

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 OSTI2020

< c+a > dislocations in {101 - 2} twins in Mg: A kinematic and energetic requirement

Consideration of published TEM and molecular dynamics studies of interactions between matrix dislocations and {101¯2} twin boundaries (TB) in hexagonal close packed metals inspires a hypothesis regarding the formation of dislocations with < c+a > Burgers vectors inside the twins. It was previously reported that when basal < a > dislocations from the matrix are swept by an advancing TB, they are converted into partial dislocations bounding basal I 1 stacking faults (SFs) within the twin. In the present work, we show that this dislocation configuration near TB renders it energetically favorable for a pair of partial dislocations to constrict into a dislocation within the twin. New in situ TEM observations confirm this hypothesis by revealing long < c+a > dislocation dipoles terminated by the TB at one end and by I1 SFs within the twin at the other end. Finally, the resulting near-screw < c+a > dislocation dipoles are observed to be glissile, enabling self-annihilation in some circumstances. Thus, rather than requiring a high energy dislocation pile-up to press two < a > dislocations through a TB into a single < c+a > dislocation, it is demonstrated that < c+a > dislocations naturally form, within {101¯2} twins in Mg, as TB migrates into a matrix, which inevitably contains dislocations with < a > Burgers vectors. Furthermore, these observations help to explain how the twins themselves may subsequently deform and why they may be particularly vulnerable to plastic, shear localization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2022

Magnetoelastic properties of multiferroic hexagonal ErMnO 3

Here, the strength and dynamics of magnetoelastic coupling through the paramagnetic (PM) – antiferromagnetic (AFM) – ferrimagnetic (FIM) transitions in multiferroic hexagonal ErMnO 3 have been investigated by Resonant Ultrasound Spectroscopy. Elastic stiffening by up to 2% below the PM – AFM transition at 80 K arises from biquadratic coupling between strain and the magnetic order parameter with relaxation times longer than ~10 -6 s for the response of spins to changes in strain. In contrast with YMnO 3 , the PM – AFM transition in ErMnO 3 is accompanied by a peak in acoustic loss immediately below the Néel point which is interpreted in terms of strain relaxation accompanying ordering of spins of Er 3+ at 4b sites. Changes in the magnetic ordering scheme at the AFM – FIM transition near 3 K are accompanied by elastic softening of ~0.03 %. During poling of the low temperature ferrimagnetic structure round magnetic hysteresis loops, small changes in elastic stiffness which arise due to the contribution of piezomagnetic and/or piezoelectric moduli are detected. Contributions of piezoelectric moduli to acoustic resonance frequencies also permit changes in the configuration of ferroelectric domains to be detected in response both to cycling through this transition and to application of a magnetic field. A peak in acoustic loss in the vicinity of 250 K is attributed to strain-mediated pinning/freezing of some aspect of the domain microstructure with an activation energy of ~0.25–0.3 eV. A return to the original elastic properties on heating to temperatures above ~250 K is interpreted in terms of backswitching of domains to the configuration they had at the start. These observations confirm the existence of subtle variations in magnetoelastic coupling behaviour relating to both the magnetic order parameters and magnetic domain structures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
DOE OSTI2024

Light-wave-controlled Haldane model in monolayer hexagonal boron nitride

In recent years, the stacking and twisting of atom-thin structures with matching crystal symmetry has provided a unique way to create new superlattice structures in which new properties emerge. In parallel, control over the temporal characteristics of strong light fields has allowed researchers to manipulate coherent electron transport in such atom-thin structures on sublaser-cycle timescales. Here we demonstrate a tailored light-wave-driven analogue to twisted layer stacking. Tailoring the spatial symmetry of the light waveform to that of the lattice of a hexagonal boron nitride monolayer and then twisting this waveform result in optical control of time-reversal symmetry breaking and the realization of the topological Haldane model in a laser-dressed two-dimensional insulating crystal. Further, the parameters of the effective Haldane-type Hamiltonian can be controlled by rotating the light waveform, thus enabling ultrafast switching between band structure configurations and allowing unprecedented control over the magnitude, location and curvature of the bandgap. This results in an asymmetric population between complementary quantum valleys that leads to a measurable valley Hall current, which can be detected by optical harmonic polarimetry. The universality and robustness of our scheme paves the way to valley-selective bandgap engineering on the fly and unlocks the possibility of creating few-femtosecond switches with quantum degrees of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Low-symmetry vacancy-related spin qubit in hexagonal boron nitride

Abstract Point defect qubits in semiconductors have demonstrated their outstanding capabilities for high spatial resolution sensing generating broad multidisciplinary interest. Hexagonal boron nitride (hBN) hosting point defect qubits have recently opened up new horizons for quantum sensing by implementing sensing foils. The sensitivity of point defect sensors in hBN is currently limited by the linewidth of the magnetic resonance signal, which is broadened due to strong hyperfine couplings. Here, we report on a vacancy-related spin qubit with an inherently low symmetry configuration, the VB2 center, giving rise to a reduced magnetic resonance linewidth at zero magnetic fields. The VB2 center is also equipped with a classical memory that can be utilized for storing population information. Using scanning transmission electron microscopy imaging, we confirm the existence of the VB2 configuration in free-standing monolayer hBN.

Chemistry↗
DOE OSTI2025

Ultrahigh-pressure crystallographic passage towards metallic hydrogen

The structural evolution of molecular hydrogen H 2 under multi-megabar compression and its relation to atomic metallic hydrogen is a key unsolved problem in condensed-matter physics. Although dozens of crystal structures have been proposed by theory, only one, the simple hexagonal-close-packed (hcp) structure of only spherical disordered H 2 , has been previously confirmed in experiments. Through advancing nano-focused synchrotron X-ray probes, here we report the observation of the transition from hcp H 2 to a post-hcp structure with a six-fold larger supercell at pressures above 212 GPa, indicating the change of spherical H 2 to various ordered configurations. Theoretical calculations based on our XRD results found a time-averaged structure model in the space group $P\bar{6}2c$ with alternating layers of spherically disordered H 2 and new graphene-like layers consisting of H 2 trimers (H 6 ) formed by the association of three H 2 molecules. Here, this supercell has not been reported by any previous theoretical study for the post-hcp phase, but is close to a number of theoretical models with mixed-layer structures. The evidence of a structural transition beyond hcp establishes the trend of H 2 molecular association towards polymerization at extreme pressures, giving clues about the nature of the molecular-to-atomic transition of metallic hydrogen. Considering the spectroscopic behaviours that show strong vibrational and bending peaks of H 2 up to 400 GPa, it would be prudent to speculate the continuation of hydrogen molecular polymerization up to its metallization.

Ji, Cheng [Center for High Pressure Science and Te↗
DOE OSTI2021

Antiferromagnetic ordering and possible lattice response to dynamic uranium valence in U 3 O 8

Determining the correct electronic structure of U 3 O 8 remains a formidable experimental and theoretical challenge. In the low-temperature phase, two crystallographic U sites are separated into a distinct 2U(V)+1U(VI) oxidation configuration. At low temperatures, the U(V) sites form a distorted honeycomb lattice, but the U(VI) sit on a triangular sublattice, suggesting potential for magnetic frustration effects. The spin configuration of the unpaired f electrons on the U(V) sites is likely antiferromagnetic (AFM) from susceptibility measurements, but this has not been confirmed. Here in this paper, we present a neutron scattering investigation of the structure and dynamics of U 3 O 8 from 1.7 to 600 K. We confirm static AFM ordering onset at between 22 and 25 K, which is present down to at least 1.7 K with AFM peaks corresponding to [0.5 1 1] and [0.5 2 2] in the orthorhombic phase. These measurements rule out static AFM order along the a axis of the Amm2 phase, a configuration previously suggested by theory. Above 100 K a quasielastic scattering channel opens that we speculate arises from a lattice relaxation response to thermally activated electron hopping. This term does not conform to a magnetic form factor, so it is not related to spin relaxations. If correct, this mechanism stabilizes a continuous valence transition from 2U(V)+1U(VI) in the low-temperature (T<600 K) orthorhombic phase to the hexagonal phase that contains only one degenerate U site, wherein the U valence can be dynamically stabilized between U(V)↔U(VI) by phonon-assisted electron hopping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

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

A novel stress-induced martensitic transformation in a single-phase refractory high-entropy alloy

High-entropy alloys (HEAs) provide a new perspective to design metastable alloys with the stress-induced martensitic transformation (SIMT) for overcoming the strength-ductility trade-off. Here, we report a novel SIMT, orthorhombic to hexagonal close-packed martensite, in a single orthorhombic refractory HEA (Ti 16 Zr 35 Hf 35 Ta 14 RHEA), showing a good yield strength-ductility matching. The analysis of the elastic distortion energy (ΔE els ) of Ti 16 Zr 35 Hf 35 Ta 14 and several other RHEAs reveals that severe lattice distortion is a key factor which causes this SIMT. Combined the “d-electron alloy design” approach with the ΔE els , the phase configuration and SIMT path in RHEAs can be well predicted. Our work brings new insights between the lattice distortion and SIMT of RHEAs, benefiting the metastable alloy development.

36 MATERIALS SCIENCE↗
DOE OSTI2024

IRIS Reimagined: Advancements in Intelligent Runtime System for Task-Based Programming

Task-based programming models are gaining traction in scientific computing. IRIS is a portable runtime system that exploits multiple heterogeneous programming systems and can discover available resources and manage multiple diverse programming systems (e.g., CUDA, Hexagon, HIP, Level Zero, OpenCL, and OpenMP) simultaneously. It accounts for the constraints of task dependencies and provides customizable scheduling policies to map those tasks to heterogeneous devices. In this paper, we present new capabilities added to IRIS to improve its portability for heterogeneous programming, build-friendliness, and performance efficiency. The new additions include vendor-specific kernel support, a runtime system with a foreign function interface to eliminate writing wrapper or boilerplate code for heterogeneous kernels, an easy-to-use and configurable CMake-based build environment, automatic and efficient data transfers and orchestration, and the Hunter and DAGGER toolchains to evaluate IRIS’s task scheduling algorithms.

Miniskar, Narasinga Rao↗
DOE OSTI2023

Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau

With the highest albedo of the land surface, snow plays a vital role in Earth's surface energy budget and water cycle. Snow albedo is primarily controlled by snow grain properties (e.g., size and shape) and light-absorbing particles (LAPs) such as black carbon (BC) and dust. The mixing state of LAPs in snow also has impacts on LAP-induced snow albedo reduction and surface radiative forcing (RF). However, most land surface models assume that snow grain shape is spherical and LAPs are externally mixed with the snow grains. This study improves the snow radiative transfer model in the Energy Exascale Earth System Model version 2.0 (E3SM v2.0) Land Model (ELM v2.0) by considering non-spherical snow grain shapes (i.e., spheroid, hexagonal plate, and Koch snowflake) and internal mixing of dust–snow, and it systematically evaluates the impacts on the surface energy budget and water cycle over the Tibetan Plateau (TP). A series of ELM simulations with different treatments of snow grain shape, mixing state of BC–snow and dust–snow, and sub-grid topographic effects (TOP) on solar radiation are performed. Compared with two remote sensing snow products derived from the Moderate Resolution Imaging Spectroradiometer, the control ELM simulation (ELM_Control) with the default configurations of spherical snow grain shape, internal mixing of BC–snow, external mixing of dust–snow, and without TOP as well as the ELM simulation with new model features (ELM_New) can both capture the overall snow distribution reasonably. Additionally, ELM_New overall shows smaller biases in snow cover fraction than ELM_Control in spring when snowmelt is important for water management. The estimated LAP-induced RF in ELM_New ranges from 0 to 19.3 W m -2 with the area-weighted average value of 1.5 W m -2 that is comparable to the reported values in existing studies. The Koch snowflake shape, among other non-spherical shapes, shows the largest difference from the spherical shape in spring when snow processes related to the surface energy budget and water cycle have high importance. The impacts of the mixing state of LAP in snow are smaller than the shape effects and depend on snow grain shape. Compared to external mixing, internal mixing of LAP–snow can lead to larger snow albedo reduction and snowmelt, which further affect the surface energy budget and water cycle. The individual contributions of non-spherical snow shape, mixing state of LAP–snow, and local topography impacts on the snow and surface fluxes have different signs and magnitudes, and their combined effects may be negative or positive due to complex and nonlinear interactions among the factors. Overall, the changes in net solar radiation in spring due to individual and combined effects range from -28.6 to 16.9 W m -2 and -29.7 to 12.2 W m -2 , respectively. This study advances understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offers guidance for improving snow simulations and RF estimates in Earth system models under climate change.

58 GEOSCIENCES↗
DOE OSTI2022

DFT Investigation of Ammonia Formation via a Langmuir–Hinshelwood Mechanism on Mo-Terminated δ-MoN(0001)

In this work, we employed density functional theory to elucidate the energetics associated with elementary steps along a Langmuir-Hinshelwood mechanism for the Haber-Bosch synthesis of ammonia from N 2 and H 2 on a hexagonal, Mo-terminated molybdenum nitride surface. Using nudged elastic band calculations, we determined the energy barriers involved in the reaction processes. An active site consisting of four nearest-neighbor Mo atoms, previously identified as an active site on similar surfaces, was chosen to investigate the reaction processes. Using this approach, we calculate a barrier of ~0.5 eV for the dissociation of N 2 . The superior activity of the dissociation of the strong N 2 bonds is rationalized based on the unique geometric and electronic configurations present at these active sites. Despite the favorable energetics for nitrogen dissociation, the energy cost for hydrogenation of NH x (0 ≤ x ≤ 2) species is shown to be energetically limiting for the formation of ammonia through the Langmuir-Hinshelwood mechanism at these sites, with elementary step activation barriers calculated to be as large as ~2 eV. A comparison to Haber-Bosch results derived from a similar γ-Mo 2 N model system suggests the relative independence of surface chemistry and bulk stoichiometry for rhombic Mo 4 active sites present on molybdenum nitrides.

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

Magnetism and interlayer bonding in pores of Bernal-stacked hexagonal boron nitride

When single-layer h-BN is subjected to a high-energy electron beam, triangular pores with nitrogen edges are formed. Because of the broken sp 2 bonds, these pores are known to possess magnetic states. Here, we report on the magnetism and electronic structure of triangular pores as a function of their size. Moreover, in the Bernal-stacked h-BN (AB-h-BN), multilayer pores with parallel edges can be created, which is not possible in the commonly fabricated multilayer AA'-h-BN. Given that these pores can be manufactured in a well-controlled fashion using an electron beam, it is important to understand the interactions of pores in neighboring layers. We find that in certain configurations, the edges of the neighboring pores remain open and retain their magnetism, and in others, they form interlayer bonds. We present a comprehensive report on these configurations for small nanopores. We find that at low temperatures, these pores have near degenerate magnetic configurations, and may be utilized in magnetoresistance and spintronics applications. In the process of forming larger multilayer nanopores, interlayer bonds can form, reducing the magnetization. Yet, unbonded parallel multilayer edges remain available at all sizes. Understanding these pores is also helpful in a multitude of applications such as DNA sequencing and quantum emission.

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

Magnetic multilayer edges in Bernal-stacked hexagonal boron nitride

Single-layer h-BN is known to have edges with unique magnetism; however, in the commonly fabricated multilayer AA'-h-BN, edge relaxations occur that create interlayer bonds and eliminate the unpaired electrons at the edge. Recently, a robust method of growing the unconventional Bernal-stacked h-BN (AB-h-BN) has been reported. Here, we use theoretical approaches to investigate the nitrogen-terminated zigzag edges in AB-h-BN that can be formed in a controlled fashion using a high-energy electron beam. We find that these "open"edges remain intact in bilayer and multilayer AB-h-BN, enabling researchers potentially to investigate these edge states experimentally. We also investigate the thermodynamics of the spin configurations at the edge by constructing a lattice model that is based on parameters extracted from a set of first-principles calculations. We find that the edge spins in neighboring layers interact very weakly, resulting in a sequence of independent spin chains in multilayer samples. By solving this model using Monte Carlo simulations, we can determine nm-scale correlation lengths at liquid-N 2 temperatures and lower. At low temperatures, these edges may be utilized in magnetoresistance and spintronics applications.

1-dimensional spin chains↗
DOE OSTI2023

CFD Benchmark of Pressure Drop in a 61-Pin Wire-Wrapped Assembly with Blocked Channels Using NekRS

Thermal hydraulic behaviors of wire-wrapped rod bundles continue to be a subject of intense research. One of the leading next-generation designs, the sodium fast reactor, relies on a hexagonal assembly of wire-wrapped fuel pins. The issue of potential channel blockage has been raised as a safety concern due to the tightly packed arrangement of the fuel pins. This has led to several recent experimental and computational studies working to quantify the potential impact on the fluid flow and heat transfer behaviors of such blockages. The objective of the present study is to benchmark the high-fidelity NekRS CFD solver in predicting pressure drop for large blockages against available experimental data. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at various low to moderate Reynolds numbers. The NekRS solver has been shown to yield exponentially decreasing spatial discretization errors with increasing polynomial order. All simulated results agreed well with measured data, which indicates that the overall methodology is adequate and consistent. The results of this benchmark study demonstrate the accuracy of NekRS for sodium fast reactor hydrodynamic simulations, increasing the confidence in its use for design, licensing, and analysis activities.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Ensuring reliable connectivity to cellular-connected UAVs with up-tilted antennas and interference coordination

To integrate unmanned aerial vehicles (UAVs) in future large-scale deployments, a new wireless communication paradigm, namely, the cellular-connected UAV has recently attracted interest. However, the line-of-sight dominant air-to-ground channels along with the antenna pattern of the cellular ground base stations (GBSs) introduce critical interference issues in cellular-connected UAV communications. In particular, the complex antenna pattern and the ground reflection (GR) from the down-tilted antennas create both coverage holes and patchy coverage for the UAVs in the sky, which leads to unreliable connectivity from the underlying cellular network. To overcome these challenges, in this paper, we propose a new cellular architecture that employs an extra set of co-channel antennas oriented towards the sky to support UAVs on top of the existing down-tilted antennas for ground user equipment (GUE). To model the GR stemming from the down-tilted antennas, we propose a path-loss model, which takes both antenna radiation pattern and configuration into account. Next, we formulate an optimization problem to maximize the minimum signal-to-interference ratio (SIR) of the UAVs by tuning the up-tilt (UT) angles of the up-tilted antennas. Since this is an NP-hard problem, we propose a genetic algorithm (GA) based heuristic method to optimize the UT angles of these antennas. After obtaining the optimal UT angles, we integrate the 3GPP Release-10 specified enhanced inter-cell interference coordination (eICIC) to reduce the interference stemming from the down-tilted antennas. Our simulation results based on the hexagonal cell layout show that the proposed interference mitigation method can ensure higher minimum SIRs for the UAVs over baseline methods while creating minimal impact on the SIR of GUEs.

3GPP↗
DOE OSTI2021

Anisotropic moiré optical transitions in twisted monolayer/bilayer phosphorene heterostructures

Moiré superlattices of van der Waals heterostructures provide a powerful way to engineer electronic structures of two-dimensional materials. Many novel quantum phenomena have emerged in graphene and transition metal dichalcogenide moiré systems. Twisted phosphorene offers another attractive system to explore moiré physics because phosphorene features an anisotropic rectangular lattice, different from isotropic hexagonal lattices previously reported. Here we report emerging anisotropic moiré optical transitions in twisted monolayer/bilayer phosphorenes. The optical resonances in phosphorene moiré superlattice depend sensitively on twist angle and are completely different from those in the constitute monolayer and bilayer phosphorene even for a twist angle as large as 19°. Our calculations reveal that the Γ-point direct bandgap and the rectangular lattice of phosphorene give rise to the remarkably strong moiré physics in large-twist-angle phosphorene heterostructures. This work highlights fresh opportunities to explore moiré physics in phosphorene and other van der Waals heterostructures with different lattice configurations.

36 MATERIALS SCIENCE↗