Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximation”

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

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A new “gold standard”: Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the “gold standard” CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.

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

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

Chemistry↗
DOE OSTI2024

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi–Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Chemistry↗
DOE OSTI2024

Experimental validation of a collision-radiation dataset for molecular hydrogen in plasmas

Quantitative spectroscopy of molecular hydrogen has generated substantial demand, leading to the accumulation of diverse elementary process data encompassing radiative transitions, electron-impact transitions, predissociations, and quenching. However, their rates currently available are still sparse, and there are inconsistencies among those proposed by different authors. In this study, we demonstrate an experimental validation of such a molecular dataset by composing a collisional-radiative model (CRM) for molecular hydrogen and comparing experimentally obtained vibronic populations across multiple levels. From the population kinetics of molecular hydrogen, the importance of each elementary process in various parameter space is studied. In low-density plasmas (electron density ne≲1017 m−3) the excitation rates from the ground states and radiative decay rates, both of which have been reported previously, determine the excited state population. The inconsistency in the excitation rates affects the population distribution the most significantly in this parameter space. However, in higher density plasmas (ne≳1018 m−3), the excitation rates from excited states become important, which have never been reported in the literature, and may need to be approximated in some way. In order to validate these molecular datasets and approximated rates, we carried out experimental observations for two different hydrogen plasmas; a low-density radio frequency heated plasma (ne≈1016 m−3) and the Large Helical Device (LHD) divertor plasma (ne≳1018 m−3). The visible emission lines from EF1Σg+, HH¯1Σg+, D1Πu±, GK1Σg+, I1Πg±, J1Δg±, h3Σg+, e3Σu+, d3Πu±,g3Σg+, i3Πg±, and j3Δg± states were observed simultaneously and their population distributions were obtained from their intensities. We compared the observed population distributions with the CRM prediction, in particular the CRM with the rates compiled by Janev et al., Miles et al., and those calculated with the molecular convergent close-coupling (MCCC) method. The MCCC prediction gives the best agreement with the experiment, particularly for the emission from the low-density plasma. However, the population distribution in the LHD divertor shows a worse agreement with the CRM than those from low-density plasma, indicating the necessity of the precise excitation rates from excited states. We also found that the rates for the electron attachment is inconsistent with experimental results. This requires further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exciting DeePMD: Learning excited-state energies, forces, and non-adiabatic couplings

We extend the DeePMD neural network architecture to predict electronic structure properties necessary to perform non-adiabatic dynamics simulations. While learning the excited state energies and forces follows a straightforward extension of the DeePMD approach for ground-state energies and forces, how to learn the map between the non-adiabatic coupling vectors (NACV) and the local chemical environment descriptors of DeePMD is less trivial. Most implementations of machine-learning-based non-adiabatic dynamics inherently approximate the NACVs, with an underlying assumption that the energy-difference-scaled NACVs are conservative fields. We overcome this approximation, implementing the method recently introduced by Richardson [J. Chem. Phys. 158, 011102 (2023)], which learns the symmetric dyad of the energy-difference-scaled NACV. Furthermore, the efficiency and accuracy of our neural network architecture are demonstrated through the example of the methaniminium cation CH 2 NH 2 + .

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

Rovibrationally resolved Rayleigh and Raman scattering cross sections for molecular hydrogen

Accurate Rayleigh and Raman scattering cross sections, tensor components, depolarization ratios, and reversal coefficients for all rovibrational transitions within the X1Σg+ ground electronic state of H2 have been calculated. Raman spectra have been generated using these data. A method for calculating Raman scattering cross sections is formulated that is valid below the ionization threshold and in the region containing resonances, which explicitly accounts for all bound and dissociative vibrational levels of the bound intermediate electronic states and approximately accounts for the ionization continuum. A representative set of cross sections is presented for incident photon energies below 15 eV and compared with existing results in the literature where possible. Convergence of our results with an increasing number of bound intermediate electronic states is demonstrated. The accuracy of the Placzek–Teller approximation is discussed. The effect of accounting for the intermediate ionization continuum is investigated. Local thermal equilibrium cross sections are calculated for Rayleigh and Raman scattering. This work represents the most accurate and complete treatment of Raman scattering for molecular hydrogen to date. A total of 9582 Rayleigh and Raman scattering cross sections have been generated and are openly available on Zenodo under an open-source Creative Commons Attribution license at https://zenodo.org/doi/10.5281/zenodo.13441471.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Linear gyrokinetic simulations of toroidal Alfvén eigenmodes in the Mega-Amp Spherical Tokamak

Linear gyrokinetic (GK) simulations using the Gyrokinetic Toroidal Code (GTC) have been performed to investigate Toroidicity-driven Alfvén Eigenmodes (TAEs) driven by the neutral beam injection (NBI) induced fast ions in the Mega-Amp Spherical Tokamak (MAST) to identify the non-perturbative and kinetic effects of thermal plasma. A specific TAE in MAST discharge 26887, with an on-axis NBI power of approximately 1.5 MW and plasma current around 800 kA, exhibited frequency chirping, and the tangential soft x-ray camera array resolved the radial mode structure peaked near |q|=1.5. Various excitation methods were used in the GTC linear simulations, illustrating this code's capability to realistically represent the mechanisms and behaviors of fast ion-driven TAEs in spherical tokamaks. The radial structures from these GK simulations closely match measurements and calculations performed using the NOVA ideal MHD code, though with the frequencies approximately 10 kHz lower, likely due to various kinetic and non-perturbative effects. The simulations measured the damping rates due to continuum damping, radiative damping, and ion Landau damping, revealing that ion Landau damping has the most significant contribution to the total damping rate of the TAE. A comparison of growth rates of TAEs excited by fast ion Maxwellian and slowing-down distributions shows that the TAEs excited by a fast ion anisotropic pitch distribution (as part of the slowing-down distributions) are more unstable compared to those excited by a Maxwellian distribution with an equivalent fast ion beta. This shows that the use of fast ion anisotropy alters the number of fast ions to be in shear Alfvén resonance, and hence, it can greatly affect the stability of TAEs. These tests can be performed with the GTC but impossible with ideal MHD simulations, highlighting the necessity of kinetic simulations such as the GTC for a precise prediction of the TAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Excess Density as a Descriptor for Electrolyte Solvent Design

Electrolytes mediate interactions between the cathode and anode and determine performance characteristics of batteries. Mixtures of multiple solvents are often used in electrolytes to achieve desired properties, such as viscosity, dielectric constant, boiling point, and melting point. Conventionally, multi-component electrolyte properties are approximated with linear mixing, but in practice, significant deviations are observed. Excess quantities can provide insights into the molecular behavior of the mixture and could form the basis for designing high-performance electrolytes. Here we investigate the excess density of commonly used Li-ion battery solvents such as cyclic carbonates, linear carbonates, ethers, and nitriles with molecular dynamics simulations. We additionally investigate electrolytes consisting of these solvents and a salt. The results smoothly vary with mole percent and are fit to permutation-invariant Redlich-Kister polynomials. Mixtures of similar solvents, such as cyclic-cyclic carbonate mixtures, tend to have excess properties that are lower in magnitude compared to mixtures of dissimilar substances, such as carbonate-nitrile mixtures. We perform experimental testing using our automated test stand, Clio, to provide validation to the observed simulation trends. We quantify the structure similarity using SOAP fingerprints to create a descriptor for excess density, enabling the design of electrolyte properties. To a first approximation, this will allow us to estimate the deviation of a mixture from ideal behavior based solely upon the structural dissimilarity of the components.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

An efficient surrogate model of secondary electron formation and evolution

This work extends the adjoint-deep learning framework for runaway electron (RE) evolution, developed by McDevitt et al. [Phys. Plasmas 32, 042503 (2025)], to account for large-angle collisions. By incorporating large-angle collisions, the framework allows the avalanche of REs to be captured, an essential component of RE dynamics. This extension is accomplished by using a Rosenbluth–Putvinski approximation to estimate the distribution of secondary electrons generated by large-angle collisions. By evolving both the primary and multiple generations of secondary electrons, the present formulation can capture both the detailed temporal evolution of a RE population beginning from an arbitrary initial momentum space distribution, along with providing approximations to the saturated growth and decay rates of the RE population. Predictions of the adjoint-deep learning framework are verified against a traditional RE solver, with good agreement present across a broad range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Grad-Shafranov model for compact quasisymmetric stellarators

A Grad-Shafranov equation (GSE) valid for compact quasisymmetric stellarators is derived by an asymptotic expansion around a vacuum field carried to first order. We obtain an equation for the existence of flux surfaces leading up to the GSE. The flux surface label must simultaneously satisfy the existence equation and the GSE, which generally leads to an overdetermined problem. We show how the overdetermined problem can be resolved within our model for a class of hybrid devices similar to that studied by Henneberg and Plunk [Phys. Rev. Res. 6, L022052 (2024)]. We are also able to solve the existence equation for flux surfaces analytically in the most general case by introducing a special coordinate system. This will enable us to carry out an optimization seeking to minimize the error in our GSE while obeying the flux surface existence equation, which will allow us to find solutions outside the class of hybrid devices. This will allow for a coarse-grained approximate search in the space of quasisymmetric equilibria that should be faster than a conventional stellarator optimization. Nevertheless, it would still be necessary to fine-tune the approximate solutions using conventional tools to obtain a more precise optimized equilibrium.

Nikulsin, N. (ORCID:0000000318611777)↗
DOE OSTI2025

The Langdon effect in laser plasmas: Absorption and conduction

A plasma heated by inverse bremsstrahlung absorption of laser light develops a non-Maxwellian electron distribution function, called the Langdon effect [A. B. Langdon, Phys. Rev. Lett. 44, 575 (1980)]. These non-Maxwellian distributions are sufficiently long-lived to impact the absorption processes itself as well as the transport of heat by electrons. The theory of the Langdon effect in a homogeneous plasma is reviewed to clarify some aspects of Langdon's derivation as well as to confirm that the widely used super-Gaussian approximation works fairly well to describe the shape of the distribution function and reduction of the absorption rate. The Langdon effect on thermal conduction in an inhomogeneous plasma is developed by considering perturbations in a homogeneous absorbing plasma, which develops a heat flux due to both temperature and density gradients. A practical theory of the heat flux is developed by fitting the results of Vlasov–Fokker–Planck simulations, which avoids several approximations that compromised the usefulness of past theoretical predictions, most critically, the effect of electron–electron collisions on the fluxes. The present fits parameterize the coefficients of the temperature gradient (thermal conductivity) and the density gradient for a plasma of any ionization state and for any laser intensity where the theory of the Langdon effect remains locally valid. It is expected that this generalized theory of heat flow in an absorbing plasma will improve the predictive capability of radiation-hydrodynamics simulations of laser-produced plasmas, especially those formed in inertial confinement fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

High-frequency electrical behavior in V 3 O 5 thin films

Vanadium oxides are known for their metal–insulator transition (MIT), with V 3 O 5 being notable for its transition temperature exceeding room temperature. At about 430 K, this material shows a change in crystal symmetry accompanied with one order of magnitude increase in its electrical conductivity and alterations in its optical properties. Although the property changes during the MIT in V 3 O 5 are less pronounced than those observed in VO 2 , its transition temperature is 90 K higher, making it appealing for applications requiring elevated temperatures. In this article, the high-frequency characteristics were determined in a V 3 O 5 two-terminal device in the range from 5 to 35 GHz. The S-parameters showed that the return loss at room temperature was close to −1.5 dB, and the isolation between ports was approximately −50 dB. At temperatures above the metal–insulator transition, the isolation decreased to around −40 dB at 35 GHz. For S11 and S22, similar behavior was observed at room temperature, with a notable change in the S-parameter phase of the device. This behavior suggests that V 3 O 5 may function well as a capacitor because the considerable change in phase could control the flow of electrical signals in devices. This property also may be used for matching purposes, especially considering its response to temperature changes. Additionally, conductivity calculation from S-parameters shows a decrease of approximately two orders of magnitude at 500 K and one order of magnitude at 300 K compared to DC values. These findings highlight V 3 O 5 potential for integration into radio frequency devices that demand consistent performance in high-temperature environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluating the gradients of localized diabatic state energies and couplings at minimum cost

We calculate the gradients of Boys diabatic state energies and couplings when the electronic vector space is generated by configuration interaction singles. Our approach follows the Lagrangian approach of Paz and Glover (rather than direct differentiation of the adiabatic-to-diabatic approaches that have been published previously). The result is that we achieve a dramatic increase in savings over previous approaches, and the present approach should be immediately useful to scientists focused on electronic relaxation, especially chemists studying electron transfer who wish to go beyond the Condon approximation. Here, a future extension to time-dependent density functional theory in the Tamm–Dancoff approximation is clear.

Chemical physics↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2025

Effect of likelihood misspecification in Gaussian process-driven autonomous experimentation

In recent years, several groups have designed Autonomous Experiment (AE) models with the aim of using them as an alternative method for neutron scattering scanning. In an AE, Gaussian processes (GPs) are most frequently used due to their interpretability, their non-parametric nature, their universal approximation, and their closed-form predictive distribution. GPs have two key components, namely, the model for the likelihood of a neutron count knowing the underlying dynamic structure factor and the acquisition function. In this paper, we investigate the impact, on the quality of an AE, of the likelihood and acquisition function choices, in energy scans and (Q, ω) ones, with respect to the signal-over-noise ratio. While we hypothesized that the quality of GP predictions would decrease when the normal to Poisson likelihood approximation breaks down at low count rates, we found that the use of the correct Poisson likelihood does not improve the quality of the data collected, as well as yields very poor results in (Q, ω) scans at low count rates. In fact, the best results are obtained with a combination of normal likelihood, including the observation noise, and the change in variance acquisition function. In addition, we find that the performance, or quality of the predictive distribution, is a misleading measure of efficiency, that is, of the quality of the data collected.

Perryman, David Elliott [Inst. Laue-Langevin (ILL)↗
DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

Octave-spanning terahertz quarter-wave plates based on over-coupled Fabry–Pérot resonances in reflective metal–dielectric–metal metasurfaces

Compact devices for broadband polarization control in the terahertz (THz) regime remain challenging to realize due to the intrinsic phase dispersion associated with birefringent materials and resonant structures. Here, we demonstrate achromatic THz quarter-wave plates with near-unity efficiency and operating under a wide range of incidence angles, based on over-coupled metal–dielectric–metal reflective metasurfaces. The devices operate as single-port anisotropic Fabry–Pérot cavities, in which the phase dispersion of the over-coupled resonances is engineered to maintain an approximately constant relative phase delay between orthogonal field components over a broad frequency range. Four complementary metasurface designs operating at an incidence angle of 45° collectively cover the 0.25–3 THz frequency range accessible to a typical THz time-domain spectroscopy system. Each device exhibits an approximately octave-wide operational bandwidth, with axial ratios below 3 dB and polarization-conversion efficiencies exceeding 80% across most of the operating band. Systematic structural optimization suppresses coupling to diffraction and guided-wave modes, thereby further extending the usable bandwidth while preserving the required phase relationship between orthogonal polarizations. The demonstrated metasurfaces are compatible with wafer-scale fabrication processes, and experimental results show excellent agreement with simulations. These findings establish over-coupled reflective metasurfaces as a robust, scalable, and versatile platform for broadband THz polarization control.

36 MATERIALS SCIENCE↗