Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice Field Theory”

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

Unitary coupled-channel three-body amplitude with pions and kaons

Three-body dynamics above threshold is required for the reliable extraction of many amplitudes and resonances from experiment and lattice QCD. The S-matrix principle of unitarity can be used to construct dynamical coupled-channel approaches in which three particles scatter off each other, rearranging two-body subsystems by particle exchange. This paper reports the development of a three-body coupled-channel, amplitude including pions and kaons. The unequal-mass amplitude contains two-body S- and P-wave subsystems (“isobars”) of all isospins, 𝐼 = 0,1/2,1,3/2,2 , and it also allows for transitions within a given isobar. The 𝑓 0 ⁡(500)⁢(𝜎),𝑓 0 ⁡(980),𝜌⁡(700),𝐾$^{*}_{0}$⁡(700)⁢(𝜅), and 𝐾*⁡(892) resonances are included, apart from repulsive isobars. Different methods to evaluate the amplitude for physical momenta are discussed. Production amplitudes for 𝑎 1 quantum numbers are shown as a proof of principle for the numerical implementation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dual Mechanism for Transient Capacitance Anomaly in Improper Ferroelectrics

Negative capacitance (NC) effects in ferroelectrics can potentially break fundamental limits of power dissipation known as “Boltzmann tyranny.” However, the origin of transient NC of ferroelectrics, which is attributed to two different mechanisms involving free-energy landscape and nucleation, is under intense debate. Here, we report the coexistence of transient NC and an S-shaped anomaly during the switching of ferroelectric hexagonal ferrites capacitor in an RC circuit. The early-stage NC arises from the nucleation process, while the late-stage S-shaped anomaly corresponds to a nascent NC associated with the free-energy landscape. The entire waveform can be reproduced using a hybrid model that simultaneously incorporates these two mechanisms. These results highlight the multivariable free-energy landscape of hexagonal ferrites that enables an abrupt change of the internal field and demonstrate that the two mechanisms are not mutually exclusive, resolving the long-standing debate. In conclusion, the behavior of the S-shaped anomaly also provides a pathway to extract parameters of free-energy landscape and switching dynamics.

Capacitance↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Mechanistic origins of excitonic properties in 2D perovskites: Implications for exciton engineering

As the field of 2D halide perovskites (HPs) matures, state-of-the-art techniques to measure important properties, such as the band gap (E g ) and exciton binding energy (E b ), continue to produce inconsistent values. Here, we tackle this long-standing problem by obtaining direct measurements of E g and E b for 31 unique HP structures. The E b values are lower than in previous literature reports and lower than expected from standard theory that assumes excitons are screened by optical-frequency dielectric constants. These low E b values are shown to be a consequence of unique screening effects, such as superlattice screening and phonon screening. Here, we find a strikingly strong correlation between E b and E g and provide design principles to a priori tune E g and E b to their optimal values. As such, this work offers a blueprint for E g -E b engineering of low-dimensional semiconductors as an even more useful replacement for simply band-gap engineering.

14 SOLAR ENERGY↗
DOE OSTI2021

Topological Mott transition in a two band model of spinless fermions with on-site Coulomb repulsion

Highlights: • The effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling lead to the 2D topological semimetal. • The ground state phase diagram is calculated. • The topological semimetal phase is characterized by the zero energy Majorana states. In the framework of mean field approach, we study topological Mott transition in a two band model of spinless fermions on a square lattice at half filling. We consider the combined effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling. The ground state phase diagram is calculated as a function of the strength of the spin–orbit Rashba coupling and the Coulomb repulsion. The spin–orbit Rashba coupling leads to a distinct phase of matter, the topological semimetal. We study a new type of phase transition between the non-topological insulator and topological semimetal states. Topological phase state is characterized by the zero energy Majorana states, which there are in defined region of the wave vectors and are localized at the boundaries of the sample. The region of existence of the zero energy Majorana states tends to zero at the point of the Mott phase transition. The zero energy Majorana states are dispersionless (they can be considered as flat bands), the Chern number and Hall conductance are equal to zero (note in two dimensional model).

36 MATERIALS SCIENCE↗
DOE OSTI2020

The electroweak phase transition: a collider target

Determining the thermal history of electroweak symmetry breaking (EWSB) is an important challenge for particle physics and cosmology. Lattice simulations indicate that EWSB in the Standard Model (SM) occurs through a crossover transition, while the presence of new physics beyond the SM could alter this thermal history. The occurrence of a first order EWSB transition would be particularly interesting, providing the needed pre-conditions for generation of the cosmic matter-antimatter asymmetry and sources for potentially observable gravitational radiation. I provide simple, generic arguments that if such an alternate thermal history exists, the new particles involved cannot be too heavy with respect to the SM electroweak temperature, nor can they interact too feebly with the SM Higgs boson. These arguments do not rely on the decoupling limit. I derive corresponding quantitative expectations for masses and interaction strengths which imply that their effects could in principle be observed (or ruled out) by prospective next generation high energy colliders. The simple, generic arguments provide a quantitative, parametric understanding of results obtained in a wide range of explicit model studies; relate them explicitly to the electroweak temperature; and delineate broad contours of collider phenomenology pertaining to a non-standard history of EWSB.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
DOE OSTI2025

Theory of terahertz pulse transmission through ferroelectric nanomembranes

An analytical model is developed to predict the temporal evolution of the lattice polarization in ferroelectric nanomembranes upon the excitation by a terahertz (THz) electromagnetic pulse of an arbitrary waveform and the concurrent transmission of the THz pulse in both linear and nonlinear regimes. It involves the use of the perturbation method to solve the equation of motion for the lattice polarization in both unclamped and strained ferroelectric nanomembranes within the framework of Landau-Ginzburg-Devonshire theory. The model is applicable to perovskite oxides such as BaTiO 3 and SrTiO 3 , wurtzite Al 1−𝑥 ⁢Sc 𝑥 ⁢N, and trigonal LiNbO 3 . Our analytical model provides a theoretical basis for determining the thermodynamic and kinetic parameters of ferroelectric materials through a THz transmission experiment. The calculation results also suggest an approach to reversing the chirality of a circularly polarized THz pulse by harnessing the resonant polarization-photon coupling in ferroelectrics. This capability of chirality reversal, along with the high tunability from a strain applied along any arbitrarily oriented in-plane axis, provides new opportunities for THz wave modulation without relying on complex metasurface designs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Parton distributions and lattice-QCD calculations: Toward 3D structure

The strong force which binds hadrons is described by the theory of quantum chromodynamics (QCD). Determining the character and manifestations of QCD is one of the most important and challenging outstanding issues necessary for a comprehensive understanding of the structure of hadrons. Within the context of the QCD parton picture, the parton distribution functions (PDFs) have been remarkably successful in describing a wide variety of processes. However, these PDFs have generally been confined to the description of collinear partons within the hadron. New experiments and facilities provide the opportunity to additionally explore the transverse structure of hadrons which is described by generalized parton distributions (GPDs) and transverse-momentum-dependent parton distribution functions (TMD PDFs). In our previous report Lin et al. (2018), we compared and contrasted the two main approaches used to determine the collinear PDFs: the first based on perturbative QCD factorization theorems, and the second based on lattice-QCD calculations. In the present report, we provide an update of recent progress on the collinear PDFs, and also expand the scope to encompass the generalized PDFs (GPDs and TMD PDFs). We review the current state of the various calculations, and consider what new data might be available in the near future. We also examine how a shared effort can foster dialog between the PDF and lattice-QCD communities, and yield improvements for these generalized PDFs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Snowmass 2021 Whitepaper: Proton Structure at the Precision Frontier

An overwhelming number of theoretical predictions for hadron colliders require parton distribution functions (PDFs), which are an important ingredient of theory infrastructure for the next generation of high-energy experiments. This whitepaper summarizes the status and future prospects for determination of high-precision PDFs applicable in a wide range of energies and experiments, in particular in precision tests of the Standard Model and in new physics searches at the high-luminosity Large Hadron Collider and Electron–Ion Collider. We discuss the envisioned advancements in experimental measurements, QCD theory, global analysis methodology, and computing that are necessary to bring unpolarized PDFs in the nucleon to the N2LO and N3LO accuracy in the QCD coupling strength. Special attention is given to the new tasks that emerge in the era of the precision PDF analysis, such as those focusing on the robust control of systematic factors both in experimental measurements and theoretical computations. Various synergies between experimental and theoretical studies of the hadron structure are explored, including opportunities for studying PDFs for nuclear and meson targets, PDFs with electroweak contributions or dependence on the transverse momentum, for incisive comparisons between phenomenological models for the PDFs and computations on discrete lattice, and for cross-fertilization with machine learning/AI approaches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Electronic structure and two-band superconductivity in unconventional high- T c cuprates Ba 2 CuO 3+δ

We report that the recently discovered cuprate superconductor Ba 2 CuO 3+δ exhibits a high T c ≃73 K at δ ≃ 0.2. The polycrystal grown under high pressure has a structure similar to La 2 CuO 4 but with dramatically different lattice parameters due to the CuO 6 octahedron compression. The crystal field in the compressed Ba 2 CuO 4 leads to an inverted Cu 3d e g complex with the d x 2 -y 2 orbital sitting below the d 3z 2 - r 2 and an electronic structure highly unusual compared to the conventional cuprates. We construct a two-orbital Hubbard model for the Cu d 9 state at hole doping x = 2δ and study the orbital-dependent strong correlation and superconductivity. For the undoped case at x = 0 , we found that strong correlation drives an orbital-polarized Mott-insulating state with the spin-1/2 moment of the localized d 3z 2 - r 2 orbital. In contrast to the single-band cuprates where superconductivity is suppressed in the overdoped regime, hole doping the two-orbital Mott insulator leads to orbital-dependent correlations and the robust spin and orbital exchange interactions produce a high-T c antiphase d-wave superconductor even in the heavily doped regime at x = 0.4 . We conjecture that Ba 2 CuO 3+δ realizes mixtures of such heavily hole-doped superconducting Ba 2 CuO 4 and disordered Ba 2 CuO 3 chains in a single-layer or predominately separated bilayer structure. Our findings suggest that unconventional cuprates with liberated orbitals as doped two-band Mott insulators can be a direction for realizing high-T c superconductivity with enhanced transition temperature T c .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Determining the gluonic gravitational form factors of the proton

The proton is one of the main building blocks of all visible matter in the Universe. Among its intrinsic properties are its electric charge, mass and spin. These properties emerge from the complex dynamics of its fundamental constituents—quarks and gluons—described by the theory of quantum chromodynamics. The electric charge and spin of protons, which are shared among the quarks, have been investigated previously using electron scattering. An example is the highly precise measurement of the electric charge radius of the proton. By contrast, little is known about the inner mass density of the proton, which is dominated by the energy carried by gluons. Gluons are hard to access using electron scattering because they do not carry an electromagnetic charge. Here we investigated the gravitational density of gluons using a small colour dipole, through the threshold photoproduction of the J/ψ particle. We determined the gluonic gravitational form factors of the proton from our measurement. Here we used a variety of models and determined, in all cases, a mass radius that is notably smaller than the electric charge radius. In some, but not all cases, depending on the model, the determined radius agrees well with first-principle predictions from lattice quantum chromodynamics. This work paves the way for a deeper understanding of the salient role of gluons in providing gravitational mass to visible matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Separation of infrared and bulk in thermal QCD

A new thermal regime of QCD, featuring decoupled scale-invariant infrared glue, has been proposed to exist both in pure-glue (N f =0) and “real-world” (N f =2+1 at physical quark masses) QCD. In this IR phase, elementary degrees of freedom flood the infrared, forming a distinct component independent from the bulk. This behavior necessitates non-analyticities in the theory. In pure-glue QCD, such non-analyticities have been shown to arise via Anderson-like mobility edges in Dirac spectra (λ IR = 0, ±λ A ≠ 0), as manifested in the dimension function d IR (λ). Here, we present the first evidence, based on lattice QCD calculation at a=0.105 fm, that this mechanism is also at work in real-world QCD, thus supporting the existence of the proposed IR regime in nature. An important aspect of our results is that, while at T = 234 MeV we find a dimensional jump between zero modes and lowest near-zero modes very close to unity (d IR = 3 to d IR ≃ 2), similar to the IR phase of pure-glue QCD, at T = 187 MeV we observe a continuous λ-dependence. This suggests that thermal states just above the chiral crossover are non-analytically (in T) connected to thermal state at T = 234 MeV, supporting the key original proposition that the transition into the IR regime occurs at a temperature strictly above the chiral crossover.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

New state of matter between the hadronic phase and the quark-gluon plasma?

Lattice-quantum chromodynamics (QCD) simulations and theoretical arguments hint at the existence of an intermediate phase of strongly interacting matter between a confined hadron gas and a deconfined quark-gluon plasma (QGP). We qualitatively and semiquantitatively explore and differentiate the phase structures in the temperature window from the QCD pseudocritical temperature 𝑇 c ≃ 160 MeV to the pure gluonic deconfinement temperature 𝑇 d ≃ 285 MeV. We propose a three-regime picture using a hadron resonance gas description augmented with the exponential spectrum of strings, corresponding to highly excited mesons and glueballs, based on the analysis of a large number 𝑁 c of colors. We estimate the entropy density from our model to confirm that the lattice-QCD data are bracketed with three regimes, i.e., a hadron gas, a QGP, and a new phase for 𝑇 c ≲ 𝑇 ≲ 𝑇 d . In this new phase, which we name a spaghetti of quarks with glueballs (SQGBs), thermal degrees of freedom of quarks are liberated, yet gluons remain confined in glueballs. Since the Hagedorn temperature 𝑇 H ∼ 285 MeV is universal in the meson and the glueball sectors, in the infinite-𝑁 c limit, the phase diagram in the plane of the baryon chemical potential and the temperature is reduced to one with the confined and deconfined phases and quarkyonic matter at high density. At large but finite 𝑁 c , an SQGB window may open between these phases. We point out that the SQGB has interesting similarities with quarkyonic matter and that this matter in the large-𝑁 c limit is confined as measured by the interaction between heavy quarks, but behaves in other respects like a quasifree gas of quarks. As a result of the extrapolation to 𝑁 c = 3, we present a revised phase diagram with the SQGB phase bounded by thermal crossovers. Finally, we give a quantitative analysis of chiral-symmetry restoration in the SQGB phase.

Color confinement↗
DOE OSTI2020

Probing the Local Site Disorder and Distortion in Pyrochlore High-Entropy Oxides

High-entropy oxides (HEOs) have attracted great interest in diverse fields because of their inherent opportunities to tailor and combine materials functionalities. The control of local order/disorder in the class is by extension a grand challenge toward realizing their vast potential. In this paper, we report the first examples of pyrochlore HEOs with five M-site cations, for Nd 2 M 2 O 7 , in which the local structure has been investigated by neutron diffraction and pair distribution function (PDF) analysis. The average structure of the pyrochlores is found to be orthorhombic Imma, in agreement with radius-ratio rules governing the structural archetype. The computed PDFs from density functional theory relaxed special quasirandom structure models are compared with real space PDFs in this work to evaluate M-site order/disorder. Reverse Monte Carlo combined with ab initio molecular dynamics and Metropolis Monte Carlo simulations demonstrates that Nd 2 (Ta 0.2 Sc 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7 is synthesized with its M-site local to nanoscale order highly randomized/disordered, while Nd 2 (Ti 0.2 Nb 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7+x exhibits a strong distortion of the TiO 6 octahedron and small degree of Ti chemical short-range order (SRO) on the subnanometer scale. Calculations suggest that this may be intrinsic, energetically favored SRO rather than due to sample processing. These results offer an important demonstration that the engineered variation of participating ions in HEOs, even among those with very similar radii, provides richly diverse opportunities to control local order/disorder motifs—and therefore materials properties for future designs. This work also hints at the exquisite level of detail that may be needed in computational and experimental data analysis to guide structure–property tuning in the emerging HEO materials class.

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

Classification of electronic nematicity in three-dimensional crystals and quasicrystals

Electronic nematic order has been reported in a rich landscape of materials, encompassing not only a range of intertwined correlated and topological phenomena but also different underlying lattice symmetries. Motivated by these findings, we investigate the behavior of electronic nematicity as the spherical symmetry of three-dimensional (3D) space is systematically lowered by the lattice environment. Here, we consider all 32 crystallographic point groups as well as four major classes of quasicrystalline point groups, given the recent observations of electronic phases of interest in quasicrystalline materials and artificial twisted quasicrystals. Valuable insights are gained by establishing a mapping between the five-component charge-quadrupolar nematic order parameter of the electronic fluid and the 3D tensorial order parameter of nematic liquid crystals. We find that a uniaxial nematic state is only generically realized in polyhedral point groups (icosahedral and cubic), with the nematic director pointing along different sets of rotational symmetry axes. Interestingly, icosahedral point groups are the only ones in which the five nematic order parameter components transform as the same irreducible representation, making them the closest analog of 3D isotropic nematics. In axial point groups, one of the nematic components is always condensed, whereas the other four components decompose into an in-plane and an out-of-plane nematic doublet, resulting in biaxial nematic ground states. Because these two nematic doublets behave as Z q clock order parameters, this allows us to identify the types of crystals and quasicrystals that can host interesting electronic nematic phenomena enabled by the critical properties of the q ≥ 4 clock model, such as emergent continuous nematic fluctuations in 3D, critical phases with quasi-long-range nematic order in two dimensions (2D), and Ashkin-Teller nematicity in 2D.

3-dimensional systems↗
DOE OSTI2025

Bottomonium properties in the quark-gluon plasma from a lattice-QCD informed 𝑇-matrix approach

Recent computations of bottomonium correlation functions with extended sources in lattice-discretized quantum chromodynamics (lQCD) provide new insights into heavy-quark dynamics at distance scales which are of the order of the inverse temperature. We analyze these results employing the thermodynamic 𝑇-matrix framework, in a continued effort to interpret lQCD data for quarkonium correlation functions in a nonperturbative and self-consistently solved quantum-many-body approach to a strongly coupled quark-gluon plasma (QGP). Its key inputs are the in-medium driving kernel (potential) of the scattering equation and an interference function which implements 3-body effects in the quarkonium coupling to the thermal medium. A simultaneous description of lQCD results for the bottomonium correlators with extended operators and the previously analyzed Wilson line correlators only requires minor refinements of the potential but calls for stronger interference effects at larger separation of the bottom quark and antiquark. We then analyze the poles of the self-consistent 𝑇 matrices on the real axis to assess the survival of the various bound states. Here, we estimate the pertinent temperatures where the poles disappear for the various bottomonium states and discuss the relation to the corresponding peaks in the bottomonium spectral functions. We also recalculate the spatial diffusion coefficient of the QGP and find it to be similar to that in our previous study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Chemistry and Interfacial Structure Promoting Quasi-van der Waals Epitaxial Growth of WS 2 Nanosheets on Sapphire for Prospective Application in Field-Effect Transistors

How do chemical and structural modifications to the supporting crystal surface affect the subsequent van der Waals (vdW) or quasi(Q)-vdW epitaxial growth of 2D nanocrystals? Developing an atomic-scale picture of such an interfacial system is crucial for understanding its impact on the physical and chemical properties of the supported 2D materials. The elucidation of the interfacial structure and chemistry needed to promote the Q-vdW epitaxial growth of 2D tungsten disulfide (WS 2 ) nanocrystals contributes to the growth mechanism understanding, thus pushing forward the integration of such atomically thin semiconductors toward real field-effect transistor applications. In addition to an atomic-force microscopy top view, we showcase a combination of X-ray techniques for a top-to-bottom investigation of the complexities of the buried interface structures. Furthermore, this approach uses X-ray photoelectron spectroscopy, X-ray standing wave excited X-ray fluorescence, and crystal truncation rod scattering to produce a highly resolved chemical-state-specific 3D atomic map for the extended interface structure of WS 2 /α-Al 2 O 3 (001). Employing these detailed analysis methods, along with density functional theory to further refine the picoscale structure, we demonstrate how two different types of interface engineering during the pregrowth stage lead to significant differences in the chemical and structural modifications to the terminal surface of c-face sapphire, which in turn leads to substantial differences in the submonolayer growth of supported WS 2 2D nanocrystals in terms of lateral domain sizes, epitaxial registry, vdW gaps, and stability.

36 MATERIALS SCIENCE↗