Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SPACE ORIENTATION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Diagenesis is key to unlocking outcrop fracture data suitable for quantitative extrapolation to geothermal targets

Exceptionally large, well-exposed sandstone outcrops in New York provide insights into folds, deformation bands, and fractures that could influence permeability, heat exchange, and stimulation outcomes of geothermal reservoir targets. Cambrian Potsdam Sandstone with <5% porosity contains decimeter-scale open, angular-limbed monoclines <0.5 km apart with associated low-porosity mm-wide cataclastic deformation bands. Crossing and abutting relationships among sub-vertical opening-mode fractures show four chronological Sets A–D, striking NNW, NE, NW, and ENE, respectively. Fracture lengths and heights range from millimeters to tens of meters. Sets A and C macro-fractures, and possibly B and D, contain quartz deposits. All sets have abundant associated quartz cemented microfractures that also record set orientations and crosscutting relations. Quartz cement deposits—evidence of diagenesis—are the key to identifying attributes of outcrop fractures suitable for extrapolation to geothermal targets in sandstones because they show which fractures formed in the subsurface. Set A fluid inclusion homogenization temperatures (120°C–129°C) are compatible with fracture at >3 km depth. Fractures are stiff and those ≥0.05 mm (Set C) and ≥0.1 mm (Set A) are open and potentially conducive to flow. Sets A and D are abundant in outcrops with close fracture spacing—0.18 m and 0.68 m, respectively—and define a rectangular connectivity network dominated by crossing and abutting X and Y nodes. Set A aperture distributions follow a power law with slope –0.8 up to 0.15 mm; other sets have lognormal distributions. Set A and D microfractures are weakly clustered, while macro-fractures commonly have 1D anticlustered (regular or periodic) arrangements at shorter length scales (<0.2 m). Sub-horizontal fractures are barren and may have formed near the surface. Fracture heights, lengths, and spatial arrangements show good trace connectivity but low open connectivity. For geothermal applications, outcrop results predict low initial well-test permeabilities owing to quartz disconnecting open fractures, but stimulation of closely spaced microfractures and partly open macro-fractures could yield high surface area for heat exchange. Quantitative extrapolation of key fracture attributes like abundance, orientation, spatial arrangement, length, and open fracture connectivity is possible from outcrops to fractured reservoirs if differing thermal histories and diagenesis are accounted for.

02 PETROLEUM↗
DOE OSTI2025

Inter-cofactor protein remodeling rewires short-circuited transmembrane electron transfer

Intraprotein electron transfer (ET) requires explicit local control of the environment of cofactors to influence their intermolecular distances, relative orientations, and redox properties. Efficient, longer-range ET often utilizes molecular orbitals of aromatic residues present in the intervening space. Here, revitalization of a vestigial ET pathway in the bacterial photosynthetic reaction center is achieved by scanning with tryptophans to uncover markedly improved routes of electron conduction in a key stabilizing step spanning 15 Å between tetrapyrrole and quinone cofactors. This ET event is maximally enhanced by pairing one or more tryptophans with a threonine to influence quinone binding and/or redox potential. Synergistic effects of these substitutions increase the yield of that ET step to ~95%. Joining these substitutions with mutant residues that improve initial ET steps dramatically enhances transmembrane charge separation via this redesigned version of a pathway that is quantitatively inactive in the native protein-cofactor complex.

Biophysical chemistry↗
DOE OSTI2026

One Earth Energy Static and Dynamic Reservoir Modeling

This report presents the static and dynamic reservoir modeling conducted for the CarbonSAFE Phase III Illinois Storage Corridor project to assess the feasibility of commercial-scale CO 2 storage in the Mt. Simon Sandstone at the One Earth Energy (OEE) site in McLean County, Illinois. Three-dimensional geocellular models of the Mt. Simon storage complex were developed in Petrel ® by integrating petrophysical log data, core analyses, and seismic surveys from the OEE #1 stratigraphic test well and two nearby wells, with multiple model versions created as new data became available. Dynamic reservoir simulations, performed using Landmark's Nexus software, progressed through three phases (preliminary, sensitivity, and UIC Class VI permit studies) evaluating injection scenarios across varying rates, well orientations, permeability models, and multi-well configurations. Results demonstrate that commercial-scale storage is feasible: three injection wells spaced approximately one mile apart can store a total of 90 million tonnes of CO 2 over 20 years, producing a combined plume with an equivalent radius of 3.2 miles and a maximum pressure-front-defined Area of Review of 178 mi 2 at the end of injection that diminishes to 34 mi 2 after 50 years of post-injection monitoring. Sensitivity analyses indicate that a 20% change in porosity or permeability yields approximately a 7% change in AoR radius, and that perforating the high-permeability arkosic zone minimizes the pressure front compared to injection in the upper Mt. Simon Sandstone.

09 BIOMASS FUELS↗
DOE OSTI2024

Large magnetoresistance and first-order phase transition in antiferromagnetic single-crystalline EuAg 4 Sb 2

Here, we present the results of a thorough investigation of the physical properties of EuAg 4 Sb 2 single crystals using magnetization, heat capacity, and electrical resistivity measurements. High-quality single crystals, which crystallize in a trigonal structure with space group 𝑅⁢$\bar{3}$𝑚, were grown using a conventional flux method. Temperature-dependent magnetization measurements along different crystallographic orientations confirm two antiferromagnetic phase transitions around 𝑇 𝑁⁢1 = 10.5K and 𝑇 𝑁⁢2 = 7.5K. Isothermal magnetization data exhibit several metamagnetic transitions below these transition temperatures. Antiferromagnetic phase transitions in EuAg 4 Sb 2 are further confirmed by two sharp peaks in the temperature-dependent heat capacity data at 𝑇 𝑁⁢1 and 𝑇 𝑁⁢2 , which shift to lower temperature in the presence of an external magnetic field. Our systematic heat capacity measurements utilizing a long-pulse and single-slope analysis technique allow us to detect a first-order phase transition in EuAg 4 Sb 2 at 7.5 K. The temperature-dependent electrical resistivity data also manifest two features associated with magnetic order. The magnetoresistance exhibits a broad hump due to a field-induced metamagnetic transition. Remarkably, the magnetoresistance keeps increasing without showing any tendency to saturate as the applied magnetic field increases, and it reaches ∼20 000% at 1.6 K and 60 T. At high magnetic fields, several magnetic quantum oscillations are observed, indicating a complex Fermi surface. A large negative magnetoresistance of about −55% is also observed near 𝑇 𝑁⁢1 . Moreover, the 𝐻−𝑇 phase diagram constructed using magnetization, heat capacity, and magnetotransport data indicates complex magnetic behavior in EuAg 4 Sb 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

XY-like incommensurate magnetic order in Ce 2 ⁢SnS 5

We report the synthesis of single crystals of Ce 2 ⁢SnS 5 through a two-stage chemical vapor transport method. The Ce 2 ⁢SnS 5 system is a member of the orthorhombic Pbam (No. 55) space group and realizes a distorted trigonal tricapped prism (TTP) crystal field around each cerium site. We characterized the sample through orientation-dependent magnetization and heat capacity measurements to probe the magnetic anisotropy in the system characteristic of XY-like anisotropic Heisenberg model behavior. Ce 2 ⁢SnS 5 furthermore enters a zero-field ordered phase under 𝑇 𝑁 =2.4 K; powder neutron diffraction measurements reveal incommensurate magnetic order near 𝑇 𝑁 . Furthermore, the system then locks into a commensurate, two-𝑞 magnetic structure below approximately 1.2 K. This commensurate structure belongs to the Shubnikov group 𝑃⁢𝑏′⁢𝑎′⁢𝑚′ ⁡(MSG 55.359) and realizes the propagation vectors $\overrightarrow{𝑞}$ = (1/3, 0, 0) and $\overrightarrow{𝑞}$ = (0, 0, 0).

Antiferromagnetism↗
DOE OSTI2025

Control over Banded Morphologies and Circular Dichroism in Chiral Halide Perovskites

Chiral halide perovskites (c-HPs) merge the chirality of organic cations with the semiconducting properties of metal halide frameworks, creating a family of chiral semiconductors with tunable chiroptoelectronic behavior. Here, we describe the impact of periodic banded morphologies of textured c-HP ( R/S -NEA) 2 PbBr 4 films (NEA = 1- (1-naphthyl)ethylammonium) on their chiroptical behavior. Due to the interplay between the crystalline and glassy phases, the c-HP film growth is driven by rhythmic precipitation, producing a distinctive controllable radial banded pattern with the ( R/S -NEA) 2 PbBr 4 inorganic planes oriented parallel to the substrate. The banded morphology can be controlled, as evidenced by the growth temperature dictating the ridge-to-ridge spacing as well as the density of banded regions. The resulting circular dichroism (CD) spectral shape, intensity, and polarity vary in a seemingly random manner across processing conditions. However, these spectral features can be explained by considering key features of the banded morphology, such as refraction of the incident light due to surface morphology, birefringence, and stacked, rotated crystallites. These effects cannot be canceled by averaging front and back CD spectra of c-HP films, and our model incorporating these effects reproduces all observed CD spectra remarkably well. The control over the c-HP morphology and prediction capabilities of our CD modeling leads to further understanding of this class of semiconductors and the possibility of exploiting structural features for light polarization control akin to enhanced metamaterials.

14 SOLAR ENERGY↗
DOE OSTI2025

Wire Plane Transparency: Commissioning the Time Projection Chambers for the Short-Baseline Near Detector

The Short Baseline Near Detector (SBND) serves as the near detector for Fermilab's Short Baseline Neutrino (SBN) programme. It is a 112-ton Liquid Argon Time Projection Chamber (LArTPC) designed to study neutrino-argon interactions and search for new physics phenomena such as sterile neutrinos. Situated just 110 m from the Booster Neutrino Beam (BNB), SBND just began its first physics run and is expected to record over two million neutrino interactions annually. Commissioning the time projection chambers of SBND is crucial for optimising detector performance. Each of the two anodes is composed of three charge-sensing wire planes which detect drifting ionisation electrons from charged particles traversing the detector. The wire planes, arranged with different 2D orientations, facilitate the 3D reconstruction of neutrino interactions. Electron transport induces signals on the first two planes, with the signal being collected by the rear plane. The first two induction planes require transparency to ensure all drift electrons reach the collection plane. This transparency depends on the bias voltages applied to the wires and the spacing between them. Fine-tuning these bias voltages ensures the optimisation of wire plane transparency. In this poster, I will demonstrate how transparency can be assessed by analysing the waveforms from cosmic-ray muon tracks traversing the detector at specific orientations.

McCusker, Bethany [Unlisted, US]↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Geothermal-integrated thermally anisotropic building envelope for energy and peak-demand reduction

Buildings consume large amounts of energy for heating and cooling, while peak electricity demand places significant stress on the power grid. This paper presents a reduced-order co-simulation framework and load-oriented supervisory control strategy for a geothermal-integrated thermally anisotropic building envelope with a ground loop (TABE+GL). In TABE+GL, a hydronic loop embedded in the building envelope is directly coupled with a geothermal ground loop, allowing for bidirectional heat exchange between the envelope, the ground, and the indoor environment. A hybrid co-simulation framework was established by coupling a reduced-order resistor–capacitor (RC) thermal network model with EnergyPlus augmented with GHEDesigner modules. The RC model generated feasible heat flux options under three operating modes, and EnergyPlus predicted sensible loads, energy use, and pump energy demand. At each simulation step, a supervisory control algorithm selected the optimal loop configuration and duty factor that maximizes useful TABE geothermal utilization without exceeding the predicted sensible load, thereby avoiding overheating or cooling. Case studies were conducted for Los Angeles, California, Charleston, South Carolina, and Denver, Colorado. Results showed that the proposed framework reduced HVAC electricity consumption by 43%–67%, natural gas use for space heating by 11%–38%, and peak electricity demand by 43%–88%. These results highlight the potential of combining reduced-order envelope modeling, direct geothermal coupling, and load-oriented supervisory control to improve whole building energy performance and reduce peak demand across diverse weather conditions.

Howard, Daniel [Southern Adventist University]↗
DOE OSTI2025

A Novel Approach for Computing Rigid Body Motion Using Linear Accelerations

Here, a novel approach is presented for computing general rigid body motion based on a few known linear accelerations. This method utilizes linear acceleration data obtained from three distinct points on the body, all within a body-fixed reference frame. The only requirement is that the three chosen points must not be collinear. A system of differential-algebraic equations is derived, combining principles of rigid body kinematics with theory of the rotation group SO(3). These equations provide a framework for numerically computing various motion parameters, including angular velocity, angular acceleration, body orientation, velocity field, acceleration field, and displacement field. By numerically solving this system of equations, we can fully characterize rigid body motion in three-dimensional space. A numerical example is provided to demonstrate the practical implementation and efficacy of the proposed technique, illustrating its potential for accurate motion computation in various applications.

42 ENGINEERING↗
DOE OSTI2025

Hot Hydrogen Exposure of U x Zr 1-x C y Nuclear Fuel: The Effect of Additional Carbon

Nuclear thermal propulsion (NTP) using hydrogen as propellant in a solid-state fission reactor to reach temperatures of up to 3000K can achieve significantly higher specific impulse than chemical propulsion. A promising fuel that is again considered for NTP is U x Zr 1-x C y . To reduce the overall mass loss as well as the uranium fuel loss, a range of experiments using additional carbon in hot hydrogen (2600 K, 6 SLPM flow rate) were performed for up to 5-6 h. Both a CH 4 addition (0.1 and 0.2 vol%) in the hydrogen stream and the use of sacrificial graphite pieces upstream of the samples were tested, with sample compositions from 5 at% UC to 30 at% UC. The results showed a strong reduction of mass loss, up to a factor of 3, as well as a reduction of surface uranium losses in the presence of additional carbon in the atmosphere. Even with additional carbon in the atmosphere, material with 30 at% UC was not stable, but material with 20 at% UC was stable for up to 6h. The sample surfaces showed texturing after processing, possibly by grain growth on the surface. This led to a separation into grains with high uranium content close to the (100)-orientation, and with low uranium content in grains close to the (111)-orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗
DOE OSTI2025

Effects of Strain and Strain Rate on Dynamic Grain Growth and Subgrain Evolution During Plastic Deformation of an Interstitial-Free Steel at 850 ° C

Here, the effects of strain and strain rate on dynamic grain growth (DGG) and subgrain evolution are reported for an interstitial-free steel deformed at 850 ° C. Microstructures produced during tension tests at true-strain rates of 10 -4 and to 10 -3 s -1 true strains ranging from 0.02 to 0.2 were preserved following deformation. These were characterized using electron backscatter diffraction (EBSD), including the application of spherical harmonic transform indexing to produce high-angular-resolution EBSD (HR-EBSD) data. HR-EBSD data resolved the small misorientation angles of subgrain boundaries while imaging much larger data fields than possible with previously available techniques. The resulting data confirmed that steady-state flow stress is inversely proportional to the average subgrain size and that subgrain boundary misorientation angle increases with strain. The following new observations are reported. The rate of DGG increased with respect to time but decreased with respect to strain as strain rate increased. This behavior is rationalized through a simple model using separate rate parameters for the effects of time and strain. Subgrain size was not constant during steady-state deformation, but decreased slowly with increasing strain. Subgrain size distributions and subgrain boundary misorientation angle distributions were measured, and both remained approximately log-normal during steady-state deformation. Subgrain evolution demonstrated no dependence on parent grain size, crystallographic orientation, or Taylor factor. These new data suggest that steady-state flow stress is more likely controlled by the dislocation density internal to subgrains than by the spacing between subgrain boundaries.

dynamic grain growth↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
DOE OSTI2024

Demonstration and Automation of Reflected Target Optical Measurement for Heliostats

Accurate optical surfaces are a primary driver of concentrated solar power plant performance. Errors in pointing and tracking mirrors, the canting of individual mirror facets, and the surface slope of the mirror itself can be caused by errors during assembly, transportation, wind loading, gravity, and many other sources. The tools that exist to measure these error sources today largely rely on fringe deflectometry (SOFAST, QDec, others), or photogrammetry with targets attached to the mirror surface. Since 2022, NREL has been developing a measurement method called the Reflected Target Non-intrusive Assessment (ReTNA) system. This system differs from most established methods in that we perform deflectometry with a pattern of coded targets, identified in space with photogrammetry. Reflected target systems have several advantages over traditional fringe deflectometry systems. Firstly, they can be operated in bright or ambient lighting, a challenge for fringe systems that use a projector and screen. Reflected target systems also can use a much lighter and less expensive target than projector-based systems. Lastly, 2D slope measurement can be solved from a single image, which leads to several advantages for accommodating faster measurements and smaller sized targets. These advantages make ReTNA particularly well-suited for applications where there are space or lighting constraints, like performing heliostat quality assurance on an assembly line. It's also useful when a lightweight, flexible system is needed, like for heliostat developers to quickly measure a new heliostat design at different orientations, to observe gravitational effects on the mirror surface shape. In the last year, significant improvements were made to this tool to make it more useful for these applications. These improvements were focused around validation of the ReTNA measurement system, and automation of the setup and measurement process. First, we present an improved ReTNA layout, for use on the heliostat assembly line. Next, we detail the various changes to the ReTNA software and computer vision methods to automate data collection in this new setup, and lessons learned from this process. The goal with this new setup is to perform a full heliostat surface characterization without removing the mirror from the assembly line. Lastly, we share results from several ReTNA validation studies undertaken over the last year. These include repeated ReTNA measurement on demonstration mirror facets, comparisons with other optical measurement tools, and some studies aimed at quantifying the uncertainty of ReTNA measurement under various constraints (mirror-target spacing, camera resolution, etc.). These results are compared with 2024 HelioCon performance targets, and our planned next steps for the ReTNA measurement system are presented.

CSP↗
DOE OSTI2025

Morphological Heterogeneity Impact of Film Solid-State Cathode on Utilization and Fracture Dynamics

Structural heterogeneity in solid-state batteries can impact material utilization and fracture mechanisms. Crystallographically oriented lithium cobalt oxide film cathodes serve as a model electrode system for exploring how void distribution contributes to stress relief and build up during cycling. Real- and reciprocal-space operando and ex situ synchrotron based experiments are utilized to understand structural changes across multiple length scales contribute to stress generation and fracture. Nanotomography uncovers a depth-dependent porosity variation in the pristine electrode and highlights preferential fracture in regions of lower porosity during delithiation. Energy-dispersive X-ray diffraction and 3D X-ray absorption near-edge spectroscopy (XANES) reveal the underutilization of cathode material in these regions. 3D XANES also confirms preferential delithiation near the sub-grain boundaries. Chemo-mechanical modeling coupled with site-specific mechanical characterization demonstrate how stress accumulation in dense regions of the electrode leads to fracture and underutilization of active material. In conclusion, our findings reveal the importance of materials design to alleviate stress in small-volume changing cathodes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Realization of Non‐Equilibrium Wurtzite Structure in Heterovalent Ternary MgSiN 2 Film Grown by Reactive Sputtering

The piezoelectric and ferroelectric applications of heterovalent ternary materials are not well explored. Epitaxial MgSiN 2 films are grown at 600 °C on (111)Pt//(001)Al 2 O 3 substrates by the reactive sputtering method using metallic Mg and Si under the N 2 atmosphere. Detailed X-ray diffraction measurements and transmission electron microscopy observations revealed that the epitaxially grown films on the substrates have a hexagonal wurtzite structure with c-axis out-of-plane orientation. The random occupation of this structure by Mg and Si differs from that of the previously reported structure in which these two cations periodically occupy the cationic sites. However, the lattice spacings closely approximate those that are previously reported, irrespective of the ordering, and they are almost comparable with those of (Al 0.8 Sc 0.2 )N. The wide bandgap of >5.0 eV in deposited MgSiN 2 is compatible with that of AlN and suggests durability against the application of strong external electric fields, possibly to induce polarization switching. In addition, MgSiN 2 is shown to have piezoelectric properties with an effective d 33 value of 2.3 pm V −1 for the first time. This work demonstrates the compositional expansion of hexagonal wurtzite to heterovalent ternary nitrides for novel piezoelectric materials, whose ferroelectricity is expected.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Emergent Above-Gap Photoluminescence in Molecularly Engineered Hybrid Bilayer Crystals

Bilayer crystals, built by stacking two-dimensional (2D) covalent monolayers, give rise to coupled excitonic states whose properties are constrained by fixed lattice symmetry and orientation. Replacing one covalent monolayer with a 2D molecular crystal––held together by non-covalent forces––overcomes this limitation, as molecular functional groups afford tunable in-plane lattice geometry, intermolecular spacing, and interlayer coupling, providing a powerful knob for exciton engineering. Here, in this study, we report four-atom-thick hybrid bilayer crystals (HBCs) synthesized by directly growing single-crystalline PDI molecular crystal atop WS 2 monolayers, which exhibit a robust photoluminescence (PL) peak 120 meV above the WS 2 optical band gap alongside a belowgap emission. Both peaks display strong polarization anisotropy—nearing unity for the above-gap emission—and maintain a perfectly linear power-law dependence up to an excitation density of ~10 7 mW/cm 2 , indicative of coexisting localized and delocalized excitonic states. Substituting PDI with a PTCDA monolayer on WS 2 fully quenches PL, demonstrating molecular control over excitonic emission. Lattice scale ab initio GW and GW-BSE calculations reveal a significantly hybridized bilayer band structure in PDI/WS 2 that supports interlayer excitonic species both above and below the WS 2 gap with strong polarization anisotropy, in excellent agreement with experiment. Our work introduces a molecule-based bilayer platform for the bottom-up design and control of excitonic phenomena in atomically thin optoelectronic and quantum materials.

2D materials↗