Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reflection coefficient”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Coarse Woody Debris Decomposition Assessment Tool: Model validation and application

Coarse woody debris (CWD) is a significant component of the forest biomass pool; hence a model is warranted to predict CWD decomposition and its role in forest carbon (C) and nutrient cycling under varying management and climatic conditions. A process-based model, CWDDAT (Coarse Woody Debris Decomposition Assessment Tool) was calibrated and validated using data from the FACE (Free Air Carbon Dioxide Enrichment) Wood Decomposition Experiment utilizing pine ( Pinus taeda ), aspen ( Populous tremuloides ) and birch ( Betula papyrifera ) on nine Experimental Forests (EF) covering a range of climate, hydrology, and soil conditions across the continental USA. The model predictions were evaluated against measured FACE log mass loss over 6 years. Four widely applied metrics of model performance demonstrated that the CWDDAT model can accurately predict CWD decomposition. The R 2 (squared Pearson’s correlation coefficient) between the simulation and measurement was 0.80 for the model calibration and 0.82 for the model validation ( P <0.01). The predicted mean mass loss from all logs was 5.4% lower than the measured mass loss and 1.4% lower than the calculated loss. The model was also used to assess the decomposition of mixed pine-hardwood CWD produced by Hurricane Hugo in 1989 on the Santee Experimental Forest in South Carolina, USA. The simulation reflected rapid CWD decomposition of the forest in this subtropical setting. The predicted dissolved organic carbon (DOC) derived from the CWD decomposition and incorporated into the mineral soil averaged 1.01 g C m -2 y -1 over the 30 years. The main agents for CWD mass loss were fungi (72.0%) and termites (24.5%), the remainder was attributed to a mix of other wood decomposers. These findings demonstrate the applicability of CWDDAT for large-scale assessments of CWD dynamics, and fine-scale considerations regarding the fate of CWD carbon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Advancing Mass Timber Buildings: Novel Methods Improve Thermal Assessment and Material Use

For nearly a century, thermal demand calculations for buildings have relied on simplified models developed to match the technical constraints of their era. The first standards, introduced in Germany and Austria in 1929, established climate zones and material conductivity coefficients that, with only incremental updates, still underpin many current assessments. Yet, methods such as Hot box testing, originally designed for lightweight insulation, continue to be applied for mass timber buildings, overlooking thermodynamic characteristics confer real-world advantages. Recent research at Oak Ridge National Laboratory incorporates updated methodologies, aligned with ASHRAE Standard 55 (ASHRAE, 2023) accounting for factors such as thermal inertia, inner surface temperatures, emissivity, solar gains, and dynamic outdoor conditions. These factors better reflect observed heating and cooling loads and highlight opportunities for efficient use of materials in mass timber construction. This work provides a framework for designing comfortable, resilient, and resource-efficient buildings while aligning with performance expectations in energy codes.

Pickett, Robert [International Mass Timber Allianc↗
DOE OSTI2020

Microbial potentiometric sensor technology for real-time detecting and monitoring of toxic metals in aquatic matrices

Considering that toxic metals can affect metabolic processes in microorganisms adversely, it can be hypothesized that these metals in water matrices would induce a decrease in metabolic activity of the biofilm microorganisms populating the surface of a sensing electrode, which could be registered as a change in the open-circuit potential (OCP) generated by the biofilm microorganisms. The goal of this study was to test this hypothesis and demonstrate the underlying principle that microbial potentiometric sensor (MPS) technology could be used for long-term and real-time monitoring and detection of rapid changes in metal concentrations in realistic aquatic environments. To address the goal, four objective were addressed: (1) a batch reactor with three graphite-based MPS electrodes was fabricated; (2) a set of single-ion solutions and one multiple ion solution were prepared reflecting realistic concentrations of metals found in electroplating wastewaters; (3) the responses of the MPS to the simultaneous presence of multiple toxic metal ions in a single solution were measured; and (4) the changes of the MPS signals to the presence of individual metal ion solutions were examined. While the hypothesis was validated, the study also revealed that the MPS was sufficiently sensitive to not only detect, but also quantify, toxic metal ion concentrations in aqueous solutions. The coefficients of determination, which were R 2 > 0.995, and responsiveness of <1 μmol/L for some toxic metal cations, strongly support the performance of MPS technology in the echelons of expensive analytical tools capable detecting and measuring trace elements.The magnitude of the MPS response was toxic metal specific. When the molar concertation normalizes the inhibition portion of the signal area, the assessed sensitivity order was: Se>Cd>Pb>Ag>Ni> Zn. The study provides valuable information for enforcement agents, environmental professionals, and wastewater treatment operators, so toxic metal pollution and its detrimental impacts can be prevented and mitigated.

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

Elastocaloric signatures of symmetric and antisymmetric strain-tuning of quadrupolar and magnetic phases in DyB 2 C 2

The adiabatic elastocaloric effect measures the temperature change of a given system with strain and provides a thermodynamic probe of the entropic landscape in the temperature-strain space. Here, we demonstrate that the DC bias strain-dependence of AC elastocaloric effect allows decomposition of the latter into symmetric (rotation-symmetry-preserving) and antisymmetric (rotation-symmetry-breaking) strain channels, using a tetragonal f -electron intermetallic DyB 2 C 2 —whose antiferroquadrupolar order breaks local fourfold rotational symmetries while globally remaining tetragonal—as a showcase example. We capture the strain evolution of its quadrupolar and magnetic phase transitions using both singularities in the elastocaloric coefficient and its jumps at the transitions, and the latter we show follows a modified Ehrenfest relation. We find that antisymmetric strain couples to the underlying order parameter in a biquadratic (linear-quadratic) manner in the antiferroquadrupolar (canted antiferromagnetic) phase, which are attributed to a preserved (broken) global tetragonal symmetry, respectively. The broken tetragonal symmetry in the magnetic phase is further evidenced by elastocaloric strain-hysteresis and optical birefringence. Additionally, within the staggered quadrupolar order, the observed elastocaloric response reflects a quadratic increase of entropy with antisymmetric strain, analogous to the role magnetic field plays for Ising antiferromagnetic orders by promoting pseudospin flips. Our results demonstrate AC elastocaloric effect as a compact and incisive thermodynamic probe into the coupling between electronic degrees of freedom and strain in free energy, which holds the potential for investigating and understanding the symmetry of a wide variety of ordered phases in broader classes of quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Atomic view of photosynthetic metabolite permeability pathways and confinement in synthetic carboxysome shells

Carboxysomes are protein microcompartments found in cyanobacteria, whose shell encapsulates rubisco at the heart of carbon fixation in the Calvin cycle. Carboxysomes are thought to locally concentrate CO 2 in the shell interior to improve rubisco efficiency through selective metabolite permeability, creating a concentrated catalytic center. However, permeability coefficients have not previously been determined for these gases, or for Calvin-cycle intermediates such as bicarbonate (HCO$^–_3$), 3-phosphoglycerate, or ribulose-1,5-bisphosphate. Starting from a high-resolution cryogenic electron microscopy structure of a synthetic β-carboxysome shell, we perform unbiased all-atom molecular dynamics to track metabolite permeability across the shell. The synthetic carboxysome shell structure, lacking the bacterial microcompartment trimer proteins and encapsulation peptides, is found to have similar permeability coefficients for multiple metabolites, and is not selectively permeable to HCO$^–_3$ relative to CO 2 . To resolve how these comparable permeabilities can be reconciled with the clear role of the carboxysome in the CO 2 -concentrating mechanism in cyanobacteria, complementary atomic-resolution Brownian Dynamics simulations estimate the mean first passage time for CO 2 assimilation in a crowded model carboxysome. Despite a relatively high CO 2 permeability of approximately 10 -2 cm/s across the carboxysome shell, the shell proteins reflect enough CO 2 back toward rubisco that 2,650 CO 2 molecules can be fixed by rubisco for every 1 CO 2 molecule that escapes under typical conditions. The permeabilities determined from all-atom molecular simulation are key inputs into flux modeling, and the insight gained into carbon fixation can facilitate the engineering of carboxysomes and other bacterial microcompartments for multiple applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Shockwave Arrival Times from Operation Redwing and Operation Upshot-Knothole

Back in the 1940s, 50s, and 60s, 210 nuclear tests were conducted in the atmosphere. As part of the testing program, EG&G performed shockwave arrival-time measurements for different shots in different operations using pressure transducers at different distances from ground zero. They then developed correlations unique to each set of data. They did not attempt to correlate all data onto one curve. As a fireball moves, material is picked up and pulled in, which is called entrainment. This extra material increases the mass of the fireball, and as the energy is spread out over a larger and larger mass, the specific internal energy of the fireball gases drops, which causes the temperature to drop as well. When the energy in the solid angle of a spherical shockwave hits the surface, it is reflected. Shortly thereafter, it catches up with the outward-moving shockwave and re-combines with it, causing it to expand at a slightly faster rate. As the shockwave grows, it approaches a hemispherical shape. The purpose of this project is to demonstrate that the shockwave arrival-time data can, in fact, be correlated onto a single function when the shockwave time-dependent solution is modified to include a geometric factor, θ, and is corrected for an entrainment coefficient.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

CFD modeling and simulation for corrosive wear of refractory in molten slag

The development of high-wear resistant refractories having minimal production costs is facilitated by characterizing the wear mechanisms associated with their corrosive wear. Static cup testing is a commonly used method for comparing the corrosion resistance performance of two or more refractory materials. Although the static cup test conditions are not as severe as dynamic tests, this study shows that the thermal gradient present within the system during heating and cooling stages serves to generate movement of the slag leading to mechanical wear. The thermal gradient within the refractory, and between the slag and the refractory, occurs during the ramp stage of the test and lasts until the soaking stage is reached bringing the system to a thermal equilibrium. Using computational fluid dynamics (CFD) capabilities embedded within ANSYS software, this study modelled and quantified the convection currents within the slag and associated shear stresses generated on the refractory walls due to the thermal gradient. A traditional ladle furnace was employed as a case study to verify the results of the studied CFD model. The corrosion rate of the refractory lining was found to depend on the mass transfer coefficient of the refractory dissolution into the slag, and a velocity term which governs the extent of corrosion at any given location. This velocity term is a function of slag viscosity, as well as the concentration gradient and/or temperature gradient at the triple points. In this study, wall shear stress was used as a reliable proxy for identifying high-velocity regions prone to excessive corrosive wear. Elevated wall shear stress near the slag/air and slag/molten steel interfaces align with observed corrosion grooves, which reflects the intensified corrosive wear at these locations.

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗
DOE OSTI2021

High-specific-power flexible transition metal dichalcogenide solar cells

Abstract Semiconducting transition metal dichalcogenides (TMDs) are promising for flexible high-specific-power photovoltaics due to their ultrahigh optical absorption coefficients, desirable band gaps and self-passivated surfaces. However, challenges such as Fermi-level pinning at the metal contact–TMD interface and the inapplicability of traditional doping schemes have prevented most TMD solar cells from exceeding 2% power conversion efficiency (PCE). In addition, fabrication on flexible substrates tends to contaminate or damage TMD interfaces, further reducing performance. Here, we address these fundamental issues by employing: (1) transparent graphene contacts to mitigate Fermi-level pinning, (2) MoO x capping for doping, passivation and anti-reflection, and (3) a clean, non-damaging direct transfer method to realize devices on lightweight flexible polyimide substrates. These lead to record PCE of 5.1% and record specific power of 4.4 W g −1 for flexible TMD (WSe 2 ) solar cells, the latter on par with prevailing thin-film solar technologies cadmium telluride, copper indium gallium selenide, amorphous silicon and III-Vs. We further project that TMD solar cells could achieve specific power up to 46 W g −1 , creating unprecedented opportunities in a broad range of industries from aerospace to wearable and implantable electronics.

42 ENGINEERING↗
DOE OSTI2022

Characteristics and mathematical modeling of apple slice drying in an electrohydrodynamic system with a needle‐plate electrode

Abstract The drying of apple slices in an electrohydrodynamic (EHD) system with a needle‐plate electrode was characterized and mathematically modeled at voltages of 30, 40, and 50 kV and with distances between needles of 5, 6, and 7 cm. The results showed that the drying rate, effective water diffusion coefficient, and rehydration ratio of the apple slices, and specific energy consumption of drying, all increased with increasing voltage. Initially, there was a period when the drying rate was approximately constant in the low‐voltage EHD, but not in the high voltage conditions. At the same voltage, the drying rate was highest when the distance between needles was 6 cm. The voltage and distance between needles had a significant ( p < .05) interactive effect on the effective water diffusion coefficient and specific energy consumption, but not on the rehydration ratio. General mathematical models for drying can be used to represent the drying behavior of apple slices. The Midilli and Kucuk model and the Logarithmic model had the best fits, with correlation coefficients of 0.998803 and 0.998267, respectively. Predictive functions for the moisture ratio of apple slices based on voltage and distance between needles were established and verified. The results indicated that the Logarithmic model could more accurately predict the moisture ratio of apple slices dried using the EHD method. Practical Applications As a new drying technology, electrohydrodynamic (EHD) drying has unequaled advantages. According to the results, using the appropriate voltage and distance between needles improves not only drying efficiency, but also the rehydration ratio of apple slices, which reflects the quality of dried products. The Logarithmic model can accurately predict the EHD drying behavior of apple slices, which is of great significance in optimizing actual production.

Li, Bowen↗
DOE OSTI2023

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗
DOE OSTI2025

Molecular diversity of dissolved organic matter reflects macroecological patterns in river networks

Deciphering dissolved organic matter (DOM) molecular complexity is crucial for understanding ecosystem function. Using the continental-scale Worldwide Hydrobiogeochemistry Observation Network for Dynamic Rivers Systems (WHONDRS) Fourier-transform ion cyclotron resonance mass spectrometry (FTICR-MS) dataset, we reveal fundamental scaling patterns of DOM chemodiversity with watershed characteristics. Analysis of 54 river sites shows local and regional watershed features significantly influence DOM chemodiversity (2500–8718 unique formulae), exhibiting consistent scaling patterns across compound classes and a novel latitudinal gradient (decreasing diversity with increasing latitude). Scaling relationships for DOM composition vary by compound class. Crucially, the scaling parameters (B, baseline chemodiversity; Z, sensitivity) are linearly interrelated. This B–Z relationship is most robust for potentially bio-labile carbohydrates (coefficient of determination R 2 ≈ 0.85), diminishing for recalcitrant, plant-derived molecules (such as lignin), and indicates (potential) biolability-dependent coupling between baseline diversity and environmental responsiveness. These quantitative scaling relationships, with scaling exponents ranging from − 2.1 to 2.2 across compound classes, enable prediction of DOM composition across watersheds, offering a framework to understand ecosystem responses to environmental change. This research bridges biogeochemistry and ecology, providing tools to anticipate molecular transformations across scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Identification and quantification of distinct active sites in Hf-Beta zeolites for transfer hydrogenation catalysis

Despite the significant progress made in characterizing different framework heteroatom sites that exist in Lewis acidic zeolites with probe molecule adsorption and spectroscopy, methods to reliably quantify site counts remain indefinite and have been primarily limited to Sn and Ti Lewis acid sites. Here, methods to quantify framework Lewis acidic Hf 4+ sites in zeolite Beta (Hf-Beta) with two Lewis base titrants (pyridine, deuterated acetonitrile) were developed using infrared (IR) spectroscopy. Lewis acid site counts for Hf-Beta zeolites were validated by measuring integrated molar extinction coefficients (IMECs; ε, cm μmol –1 ) on Sn-Beta zeolites using identical Lewis base titrants to benchmark site counts with established literature procedures to quantify Lewis acid sites in Beta zeolites from IR spectra. Highlighting the importance of benchmarking active site counts against well-established experimental protocols, IMECs of CD 3 CN bound to open (ε(Sn; 2316 cm –1 ): 1.80 ± 0.25) and closed (ε(Sn; 2308 cm –1 ): 3.76 ± 0.33) Sn sites were ~1.8x larger than those previously reported while total Lewis acid site counts agreed with those measured by pyridine (ε(Sn; 1451 cm –1 ): 1.58 ± 0.16) on six different Sn-Beta zeolites. IMECs measured for IR peaks reflecting pyridine bound to Lewis acidic Hf sites (ε(Hf; 1448 cm –1 ): 1.54 ± 0.21) and CD 3 CN bound to open (ε(Hf; 2313 cm –1 ): 2.40 ± 0.22) and closed (ε(Hf; 2307 cm –1 ): 3.55 ± 0.41) Hf sites, gave similar counts for the total number of Lewis acidic sites across six Hf-Beta zeolites (Si/Hf = 100–413). Consistent with previous reports with Sn-Beta catalysts where open Sn sites are responsible for catalytic turnover, apparent first and zero-order MPVO rate constants (0.01–1 M cyclohexanone in 2-butanol; per total Hf, 373 K) correlated with the total number of open Hf sites, per total Hf, but not with the total number of closed Hf sites or total Lewis acid site counts. Measured initial MPVO rates (0.1 M cyclohexanone in 2-butanol, per open Hf, 373 K) were ~25x higher on hydrophobic Hf-Beta-F than on hydrophilic Hf-Beta-OH zeolites. Overall, the apparent first-order MPVO rate constants (2-butanol solvent, per open Hf, 373 K) were ~6x higher on Hf-Beta-F than on Hf-Beta-OH zeolites. The characterization methods reported here enable normalization of MPVO turnover rates on Sn- and Hf-Beta zeolites by their number of open sites. Finally, this enables performing quantitative rate comparisons across Lewis acid zeolites of varying active site identity, solvation, and pore topology used in liquid-phase catalysis.

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

Free fermion entanglement with a semitransparent interface: the effect of graybody factors on entanglement islands

We study the entanglement entropy of free fermions in 2d in the presence of a partially transmitting interface that splits Minkowski space into two half-spaces. We focus on the case of a single interval that straddles the defect, and compute its entanglement entropy in three limits: Perturbing away from the fully transmitting and fully reflecting cases, and perturbing in the amount of asymmetry of the interval about the defect. Using these results within the setup of the Poincaré patch of AdS _2 2 statically coupled to a zero temperature flat space bath, we calculate the effect of a partially transmitting AdS _2 2 boundary on the location of the entanglement island region. The partially transmitting boundary is a toy model for black hole graybody factors. Our results indicate that the entanglement island region behaves in a monotonic fashion as a function of the transmission/reflection coefficient at the interface.

Kruthoff, Jorrit↗
DOE OSTI2020

Equation of State of CO 2 Shock Compressed to 1 TPa

Equation-of-state (pressure, density, temperature, internal energy) and reflectivity measurements on shock-compressed CO 2 at and above the insulating-to-conducting transition reveal new insight into the chemistry of simple molecular systems in the warm-dense-matter regime. CO 2 samples were precompressed in diamond-anvil cells to tune the initial densities from 1.35 g/cm 3 (liquid) to 1.74 g/cm 3 (solid) at room temperature and were then shock compressed up to 1 TPa and 93 000 K. Lastly, variation in initial density was leveraged to infer thermodynamic derivatives including specific heat and Gruneisen coefficient, exposing a complex bonded and moderately ionized state at the most extreme conditions studied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A numerical study of the spectral radiative properties of packed bed with mixed bauxite and silica spheres

Bauxite and silica particles have gained increasing attention for applications in the field of concentrated solar power. In this work, a Monte Carlo ray-tracing simulation is performed to predict the radiative properties (absorptance, reflectance, and transmittance) of packed beds with mixed bauxite and silica spherical particles at wavelengths of 0.5, 2.6, 9.0, and 9.35 µm. These wavelengths are representative for the visible, near-infrared, and the mid-infrared regions that are important for solar and thermal radiation, respectively. A repeating unit-column approach is used to mathematically represent the particle bed. The effects of particle mixing ratio, volume fraction, and wavelength on the predicted radiative properties are examined. The obtained radiative properties are inputted to an inverse method to retrieve the effective absorption and scattering coefficients as well as the scattering albedo, which may be later used in a continuous-scale radiative heat transfer analysis. Furthermore, the independent scattering model is used to obtain the absorption coefficient and scattering albedo based on the absorption and scattering cross sections predicted by a Monte Carlo algorithm for a single particle. It is shown that the independent scattering model underpredicts the scattering coefficient for opaque particles but overpredicts the scattering coefficient for semitransparent particles for sufficiently high particle volume fractions. Here, the radiative properties calculated from the independent scattering model are compared to the full Monte Carlo simulation of the particle bed to examine the influence of particle mixing on dependent scattering.

42 ENGINEERING↗
DOE OSTI2020

Goldsim Modeling of Vadose Zone Transport for E-Area Naval Reactor Component Disposal Areas: Model Description and Benchmarking

This report documents the development of a GoldSim® model of flow and radionuclide transport to the water table through the Naval Reactor Components Disposal Area (NRCDA) waste disposal sites and underlying vadose zones. The model is designed to be used for Monte Carlo uncertainty analysis in support of the E-Area Performance Assessment (PA). This report describes the model and shows results obtained from benchmarking the model to best-estimate deterministic results obtained using a PORFLOW model of NRCDA vadose zone transport. The PORFLOW model is three-dimensional while the GoldSim model is a simplified one-dimensional treatment. Nevertheless, the GoldSim model was able to accurately reproduce PORFLOW results with some adjustment to the nominal dispersion coefficient and vadose zone flow area used as “tuning” parameters. An example of the results obtained comparing GoldSim and PORFLOW calculation of releases of I-129, Tc-99, C-14 and Ni-59 from waste disposal containers at the 643-26E site is shown in Figure 1 below. For all of the test cases evaluated, GoldSim predicted peak concentrations within 6% of the PORFLOW values and peak times agreed within 8% with the majority of the results in better agreement. The close agreement between the two models provides confidence that GoldSim will give results accurately reflecting the behavior of releases from the NRCDA under off-normal operating conditions for sensitivity and uncertainty analysis

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗
DOE OSTI2026

Unexpected symmetry breaking in ferroelectric wurtzite thin films on silicon observed by optical second harmonic generation

Optical second harmonic generation (SHG) exists as a popular tool for probing materials with broken inversion symmetry in the physical and biological sciences. SHG polarimetry can reveal material anisotropy with high sensitivity, including point group and phase transitions. Here, we probe ferroelectric wurtzite films with a nominal 6 mm symmetry under a normal reflection geometry using SHG microscopy and discover an unexpected symmetry breaking. Symmetry considerations would normally forbid the detection of the SHG signal when light propagates along the polar 6-fold rotation axis. Yet a uniquely anisotropic SHG response is observed in this geometry in Al 1-x B x N, Al 1-x Sc x N, Zn 1-x Mg x O, and AlN/Al 1-x Sc x N heterostructures grown on silicon that can be modeled by an average monoclinic symmetry of point group m. A significant enhancement of the SHG signal corresponding to up to a 5.3× increase in the SHG intensity (hence ∼2.3 × in effective SHG tensor coefficient) is observed at antiparallel polar domain walls, suggesting local cooperative alignment of symmetry-breaking structural distortions. Namely, it is found that the monoclinic mirror plane is oriented predominantly perpendicular to the walls. Such increases in domain wall SHG thus reveal that subtle symmetry breaking can be a pathway to large property enhancements.

36 MATERIALS SCIENCE↗