Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Charging Behavior”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2021

Revealing the competition between charge density wave and superconductivity in CsV 3 Sb 5 through uniaxial strain

In this paper we report the effect of uniaxial strain ε applied along the crystalline a axis on the newly discovered kagome superconductor CsV 3 Sb 5 . At ambient conditions, CsV 3 Sb 5 shows a charge-density wave (CDW) transition at T CDW = 94.5K and superconducts below T c = 3.34K. In our paper, when the uniaxial strain ε is varied from –0.90% to 0.90%,T c monotonically increases by ~33% from 3.0 to 4.0 K, giving rise to the empirical relation T c (ε) = 3.4 + 0.56ε + 0.12ε 2 . On the other hand, for ε changing from –0.76% to 1.26%,TCDW decreases monotonically by ~10% from 97.5 to 87.5 K with T CDW (ε) = 94.5 – 4.72ε – 0.60ε 2 . The opposite response of T c and T CDW to the uniaxial strain suggests strong competition between these two orders. Comparison with hydrostatic pressure measurements indicate that it is the change in the c axis that is responsible for these behaviors of the CDW and superconducting transitions, and that the explicit breaking of the sixfold rotational symmetry by strain has a negligible effect. Combined with our first-principles calculations and phenomenological analysis, we conclude that the enhancement in Tc with decreasing c is caused primarily by the suppression of T CDW , rather than strain-induced modifications in the bare superconducting parameters. We propose that the sensitivity of T CDW with respect to the changes in the c axis arises from the impact of the latter on the trilinear coupling between the $M$$^{+}_{1}$ and the $L$$^{–}_{2}$ phonon modes associated with the CDW. Altogether, our paper reveals that the c-axis lattice parameter, which can be controlled by both pressure and uniaxial strain, is a powerful tuning knob for the phase diagram of CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2026

Pressure‐Induced Structural and Magnetic Evolution in Layered Antiferromagnet YbMn 2 Sb 2

Electronic states under pressure exhibit unconventional spin and charge dynamics that provide a powerful route to uncover exotic phases in quantum materials. Here, we present the structural, magnetic, and electronic evolution of YbMn 2 Sb 2 under pressure. Single-crystal X-ray diffraction reveals a pressure-induced structural transition from the space group trigonal P$\overline{3}$m1 to the monoclinic P2 1 /m phase near 3.5 GPa, which remains stable up to 10 GPa. Magnetization measurements display an anomalously weak net magnetic moment and the absence of Curie–Weiss behavior up to 400 K, suggesting the formation of short-range Mn moment pairs that cancel macroscopically and subsequently evolve into long-range order upon cooling. Temperature-dependent resistivity shows semiconducting behavior with a transition at ∼119 K at ambient pressure, while pressure induces a dramatic suppression of resistance and the emergence of metallic-like temperature dependence, stabilized beyond 5 GPa. This pressure-driven semiconductor-metal transition is consistent with our density functional theory calculations, confirming the closing of the band gap under compression. Neutron diffraction under pressure identifies an incommensurate magnetic structure with antiparallel correlations between paired spins. Together, these results demonstrate how pressure-driven structural tuning and competing exchange interactions stabilize unconventional magnetic states in this low-dimensional magnetic semiconductor.

antiferromagnetic semiconductor↗
DOE OSTI2020

Hole antidoping of oxides

In standard doping, adding charge carrier to a compound results in a shift of the Fermi level towards the conduction band for electron doping and towards the valence band for hole doping. In this paper, we discuss the curious case of antidoping, where the direction of band movements in response to doping is reversed. Specifically, $\textit{p}$-type antidoping moves the previously occupied bands to the principal conduction band resulting in an increase of band gap energy and reduction of electronic conductivity. We find that this is a generic behavior for a class of materials: early transition and rare-earth metal (e.g., Ti, Ce) oxides where the sum of composition-weighed formal oxidation states is positive; such compounds tend to form the well-known electron-trapped intermediate bands localized on the reduced cation orbitals. What is less known is that doping by a hole annihilates a single trapped electron on a cation. The latter thus becomes electronically inequivalent with respect to the normal cation in the undoped lattice, thus representing a symmetry-breaking effect. We give specific theoretical predictions for target compounds where hole antidoping might be observed experimentally: Magnéli-like phases (i.e., CeO 2–x and TiO 2–x ) and ternary compounds (i.e., Ba 2 Ti 6 O 13 and Ba 4 Ti 12 O 27 ), and note that this unique behavior opens the possibility of unconventional control of materials conductivity by doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mechanisms and site requirements for NO and NH 3 oxidation on Cu/SSZ-13

Two series of Cu/SSZ-13 catalysts were synthesized via aqueous solution and solid-state ion exchange using SSZ-13 supports of varying Si/Al ratios. The isolated and multinuclear Cu content of these catalysts were determined by H 2 temperature programmed reduction (H 2 -TPR). Multinuclear Cu in these catalysts, including in situ Cu- dimers formed from ZCu II OH coupling and permanent CuO clusters, are active species for dry NO oxidation. NH 3 oxidation on these catalysts follows an internal SCR (i-SCR) mechanism, i.e., a portion of NH 3 is first oxidized to NO, then NO is selectively reduced by the remaining NH 3 to N 2 . NH 3 oxidation displays distinct kinetic behavior below ~300 °C and above ~400 °C. At low temperature the results indicate that NH 3 -solvated mobile Cu-ions are the active centers. CuO clusters, when present, also contribute to the low temperature activity by catalyzing NH 3 oxidation to NO. At high temperature, in situ Cu-dimers and CuO clusters catalyze NH 3 oxidation to NO, and isolated Cu-ions catalyze SCR to realize the cascade turnovers. For both NO and NH 3 oxidation, Cu-dimers balanced by framework charges of close proximity appear to be more active than Cu-dimers balanced by distant framework charges. However, the former Cu-dimers are less stable than the latter and tend to split into monomers in the presence of vicinal Brønsted acid sites. Via density functional theory (DFT) calculations, the i-SCR mechanism for low temperature NH 3 oxidation, i.e., the energetic favorability for the involvement of the NO intermediate, is justified. Furthermore, the DFT results also agree with experimental data that the formation of Cu-dimers from ZCu II OH dimerization is essential for NH 3 oxidation at high temperature.

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

New precious metal containing normal spinels: LiRhRu 1-x Ir x O 4 , LiFeIr 1-x Ru x O 4 , and LiCoIr 1-x Ru x O 4

Large spin-orbit-coupled cations in geometrically frustrated crystal structures have the most suitable setting for exploring novel exotic states of matter. Spinel oxides (AM 2 O 4 ) are well-known examples of geometrically frustrated systems. In this study, we report for the first time the synthesis of compositions LiRhRu 1-x Ir x O 4 (x = 0–0.5), LiFeIr 1-x Ru x O 4 (x = 0–0.5), and LiCoIr 1-x Ru x O 4 (x = 0–0.3) containing precious metal cations on edge-sharing octahedral M-sites, and systematically investigate their magnetic and electrical properties. 57 Fe Mössbauer spectroscopy revealed that iron is trivalent in all LiFeIr 1-x Ru x O 4 solid solutions. Magnetic measurements indicate deviations from theoretical spin-only magnetic moment values, indicating the influence of spin-orbit coupling owing to the presence of 4d and 5d block elements. The LiFeIr 1-x Ru x O 4 series shows spin-glass-like freezing behavior with T g ≈ 20 K, and a small frustration index (f ≈ 1-2), indicating that the frustration originates from site disorder. LiRhRu 1-x Ir x O 4 and LiCoIr 1-x Ru x O 4 exhibit strongly geometrically frustrated magnetism. Electrical resistivity measurements as a function of temperature indicate that all phases are semiconducting. Seebeck coefficient measurements show that LiRhRu 1-x Ir x O 4 and LiFeIr 1-x Ru x O 4 are p-type semiconductors with holes as the major charge carriers. A sign reversal of the Seebeck coefficient indicates both holes and electrons as carriers for LiCoIr 1-x RuxO 4 (x = 0–0.2), but only holes as major carriers for x = 0.3. Here, the Seebeck coefficient and power factor increase drastically in the LiRhRu 1-x Ir x O 4 solid solution with Ir substitution, reaching a maximum of ≈ +125 μV/K and ≈2.3×10 -6 W/mK 2 at ∼650 K for x = 0.5.

Charge carriers↗
DOE OSTI2025

Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide (SiC) [Zhao, J. Nature 2024, 625 (7993), 60−65, 10.1038/s41586-023-06811-0] provided an important step toward integration of the graphene-based system into active components in postsilicon micro- and nanoelectronics. However, the exact atomic-scale structure and complex bonding configurations of the first epitaxial graphene carbon layer (C buffer ) remain an open problem. Our recent report [Kolmer, M. Communications Physics 2024, 7 (1), 16, 10.1038/s42005-023-01515-3] has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the C buffer –SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Bias voltage and epitaxial graphene thickness-dependent characterization of the collective C buffer –SiC interface showed that “Si” vacancy sites beneath C buffer are stable under STM electric fields. Moreover, the vacancies introduce localized electronic states below the Fermi level, thereby enhancing the charge-transfer phenomenon across the interface.

Thupakula, Umamahesh [Ames Laboratory (AMES), Ames↗
DOE OSTI2023

Ultrafast terahertz emission from emerging symmetry-broken materials

Abstract Nonlinear optical spectroscopies are powerful tools for investigating both static material properties and light-induced dynamics. Terahertz (THz) emission spectroscopy has emerged in the past several decades as a versatile method for directly tracking the ultrafast evolution of physical properties, quasiparticle distributions, and order parameters within bulk materials and nanoscale interfaces. Ultrafast optically-induced THz radiation is often analyzed mechanistically in terms of relative contributions from nonlinear polarization, magnetization, and various transient free charge currents. While this offers material-specific insights, more fundamental symmetry considerations enable the generalization of measured nonlinear tensors to much broader classes of systems. We thus frame the present discussion in terms of underlying broken symmetries, which enable THz emission by defining a system directionality in space and/or time, as well as more detailed point group symmetries that determine the nonlinear response tensors. Within this framework, we survey a selection of recent studies that utilize THz emission spectroscopy to uncover basic properties and complex behaviors of emerging materials, including strongly correlated, magnetic, multiferroic, and topological systems. We then turn to low-dimensional systems to explore the role of designer nanoscale structuring and corresponding symmetries that enable or enhance THz emission. This serves as a promising route for probing nanoscale physics and ultrafast light-matter interactions, as well as facilitating advances in integrated THz systems. Furthermore, the interplay between intrinsic and extrinsic material symmetries, in addition to hybrid structuring, may stimulate the discovery of exotic properties and phenomena beyond existing material paradigms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Out-of-equilibrium charge redistribution in a copper-oxide based superconductor by time-resolved X-ray photoelectron spectroscopy

Charge-transfer excitations are of paramount importance for understanding the electronic structure of copper-oxide based high-temperature superconductors. In this study, we investigate the response of a Bi 2 Sr 2 CaCu 2 O 8+δ crystal to the charge redistribution induced by an infrared ultrashort pulse. Element-selective time-resolved core-level photoelectron spectroscopy with a high energy resolution allows disentangling the dynamics of oxygen ions with different coordination and bonds thanks to their different chemical shifts. Our experiment shows that the O 1s component arising from the Cu–O planes is significantly perturbed by the infrared light pulse. Conversely, the apical oxygen, also coordinated with Sr ions in the Sr-O planes, remains unaffected. This result highlights the peculiar behavior of the electronic structure of the Cu–O planes. It also unlocks the way to study the out-of-equilibrium electronic structure of copper-oxide-based high-temperature superconductors by identifying the O 1s core-level emission originating from the oxygen ions in the Cu–O planes. This ability could be critical to gain information about the strongly-correlated electron ultrafast dynamical mechanisms in the Cu–O plane in the normal and superconducting phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Active learning of polarizable nanoparticle phase diagrams for the guided design of triggerable self-assembling superlattices

Polarizable nanoparticles are of interest in materials science because of their rich and complex phase behavior that can be used to engineer nanostructured materials with long-range crystalline order. To understand and rationally navigate the design space of polarizable nanoparticles for self-assembling highly ordered superlattices, we developed a coarse-grained computational model to describe the nanoparticle-nanoparticle interactions in implicit solvent and employ the computationally efficient image method to model many-body polarization interactions. We conducted high-throughput virtual screening over a five-dimensional particle design space spanned by temperature, particle size, particle charge, particle dielectric, and solvent dielectric using enhanced sampling molecular dynamics calculations within an active learning framework to efficiently map out the regions of thermodynamic stability of the self-assembled aggregates. We validate our predictions in comparisons against small angle x-ray scattering measurements of gold nanoparticles surface functionalized with metal chalcogenide ligands. Lastly, we use our validated phase maps to computationally design switchable nanostructured materials capable of triggered assembly and disassembly as a function of temperature and solvent dielectric with potential applications as sensors, smart windows, optoelectronic devices, and in medical diagnostics.

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

Emergent tetragonality in a fundamentally orthorhombic material

Symmetry plays a key role in determining the physical properties of materials. By Neumann’s principle, the properties of a material remain invariant under the symmetry operations of the space group to which the material belongs. Continuous phase transitions are associated with a spontaneous reduction in symmetry. Less common are examples where proximity to a continuous phase transition leads to an increase in symmetry. We find signatures of an emergent tetragonal symmetry close to a charge density wave (CDW) bicritical point in a fundamentally orthorhombic material, ErTe 3 , for which the two distinct CDW phase transitions are tuned via anisotropic strain. We first establish that tension along the a axis favors an abrupt rotation of the CDW wave vector from the c to a axis and infer the presence of a bicritical point where the two continuous phase transitions meet. We then observe a divergence of the nematic elastoresistivity approaching this putative bicritical point, indicating an emergent tetragonality in the critical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reversal of spin-polarization near the Fermi level of the Rashba semiconductor BiTeCl

Abstract Spin–orbit coupling forms the physical basis for quantum materials with non-trivial topology and potential spintronics applications. The Rashba interaction is a textbook model of spin–orbit interactions, with charge carriers undergoing linear, isotropic spin-splitting in momentum space. Recently, non-centrosymmetric semiconductors in the family BiTe X ( X = Cl, Br, I) have been identified as exemplary Rashba materials due to the strong splitting of their bulk bands, yet a detailed investigation of their spin textures, and their relationships to local crystal symmetry, is currently lacking. We perform high-efficiency spin-resolved photoemission spectroscopy to directly image the spin texture of surface states of BiTeCl, and we find dramatic deviations from idealized behavior, including a reversal of the spin-polarization near the Fermi level. We show that this behavior can be described by higher-order contributions to the canonical Rashba model with the surface states localized to individual trilayers of the crystal. Due to the prominence of these effects near the Fermi level, they should have a strong impact on the spin-dependent transport of carriers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

How well do one-electron self-interaction-correction methods perform for systems with fractional electrons?

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z σ ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of $H$$^{+}_{2}$ and $H$$^{+}_{2}$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li + to F – at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.

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

Extracting the speed of sound in quark–gluon plasma with ultrarelativistic lead–lead collisions at the LHC

Abstract Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark–gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we determine the speed of sound in an extended volume of quark–gluon plasma using lead–lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb −1 . The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 ± 0.002 ( stat ) ± 0.016 ( syst ) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 ± 8 ( syst ) MeV . The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Thermal evolution of ferroelectric behavior in epitaxial Hf 0.5 Zr 0.5 O 2

In this paper, we report a cryogenic-temperature study on the evolution of the ferroelectric properties of epitaxial Hf 0.5 Zr 0.5 O 2 thin films on silicon. Wake-up, endurance, and fatigue of these films are found to be intricately correlated, strongly hysteretic, and dependent on available thermal energy. Field-dependent measurements reveal a decrease in polarization with temperature, which has been determined not to be an intrinsic change of the material property, rather a demonstration of the increase in the coercive bias of the material. Our findings suggest that a deficiency in thermal energy suppresses the mobility of defects presumed to be oxygen vacancies during wake-up and trapped injected charge during fatigue, which is responsible for polarization evolution during cycling. This permits accelerated wake-up and fatigue effects at high temperatures where thermal energy is abundant but delays these effects at cryogenic temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Piezoelectric response and failure behavior of cement paste under external loading

We study the elastic behavior and rupture mechanism of Calcium Silicate Hydrates (C-S-H) in response to tensile loading. For the three tensile axes, initial defects are spotted in the interlayer regions. However, the crack propagation is rather different depending on the axis of tensile. The cup-and-cone failure surface shows the ductile feature of C-S-H in the x and y directions; while the devastating rupture mode confirms the brittle feature of C-S-H in the z direction. Piezoelectricity was investigated using the charge fluctuation method. The polarization variation in response to small strain was also investigated and piezoelectric coefficients were derived. Furthermore, dependence of elastic and piezoelectric properties over external pressure was analyzed. Results reveal an improvement of the mechanical properties and the ultimate strength especially in the z direction; whereas the non-piezoelectricity remains unchanged despite the rearrangement of atoms induced by the external applied pressure.

36 MATERIALS SCIENCE↗