Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

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

Thermodynamics of MgO Atomic Layer Deposition Surface Reactions

The realities of atomic layer deposition (ALD) surface reactions often deviate from the simple ligand exchange frequently used to illustrate the technique. A detailed understanding of these reactions is necessary to develop greater surface synthetic control including chemical selectivity for patterning or to target defects. Here, the thermodynamics of surface reactions relevant to MgO ALD were investigated using pyroelectric calorimetry to measure the time-resolved heat generation. These reactions show exothermic heat generation of 0.12 mJ/cm 2 and 0.15 mJ/cm 2 for alternating Mg(CpEt) 2 and H 2 O reactions, respectively. The total reaction heat closely matches the standard reaction enthalpy for bulk MgO, supplemented with first-principles molecular calculations. First-principles models further reveal that while simple ligand exchange is favorable during surface reactions, the required increase in Mg-coordination number from two in the precursor to six in bulk MgO requires additional coverage-dependent surface reactions, which depend on the availability and stability of proximal surface hydroxyls.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Discrimination of Hexane Isomers by Temperature Swing Adsorption in a Rigid Aluminum Metal–Organic Framework

The efficient separation of alkane isomers with similar physicochemical properties remains a persistent challenge for the petrochemical industry. Adsorptive separation using metal− organic frameworks (MOFs) offers an energy-efficient alternative to conventional distillation. Herein, we report temperature swing discrimination of hexane isomers with different degrees of branching using MIL-120, a rigid aluminum pyromellitate-based MOF. MIL-120 features uniform one-dimensional channels with an aperture of ∼5.5 Å. At 30 °C, it selectively adsorbs linear and monobranched hexanes while excluding the dibranched isomer. Upon heating to 120 °C, both mono- and dibranched isomers are completely excluded, whereas linear hexane remains strongly adsorbed. Breakthrough experiments validate the temperature swing separation performance. Adsorption heat analysis combined with ab initio calculations provides a quantitative measure of distinct differences in adsorption enthalpies, binding energies, and diffusion barriers responsible for the observed separation efficiency, highlighting the potential of this MOF for efficient separation of alkane isomers via temperature swing adsorption.

Adsorption↗
DOE OSTI2022

Hybrid Nanocrystals of Small Molecules and Chemically Disordered Polymers

Organic crystals formed by small molecules can be highly functional but are often brittle or insoluble structures with limited possibilities for use or processing from a liquid phase. A possible solution is the nanoscale integration of polymers into organic crystals without sacrificing long-range order and therefore function. This enables the organic crystals to benefit from the advantageous mechanical and chemical properties of the polymeric component. We report here on a strategy in which small molecules co-crystallize with side chains of chemically disordered polymers to create hybrid nanostructures containing a highly ordered lattice. Synchrotron X-ray scattering, absorption spectroscopy, and coarse-grained molecular dynamics simulations reveal that the polymer backbones form an “exo-crystalline” layer of disordered chains that wrap around the nanostructures becoming a handle for interesting properties. The morphology of this “hybrid bonding polymer” nanostructure is dictated by the competition between the polymers’ entropy and the enthalpy of the lattice allowing for control over the aspect ratio of the nanocrystal by changing the degree of polymer integration. Further, we observed that nanostructures with an exo-crystalline layer of polymer exhibit enhanced fracture strength, self-healing capacity, and dispersion in water, which benefits their use as light-harvesting assemblies in photocatalysis. Guided by computation, future work could also explore these hybrid nanostructures as building blocks for functional materials.

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

Enhancing Room-Temperature Hydrogen Storage Properties in Ti–V–Cr–Nb–Mo/Fe Multielement Alloys via Substitutional Tuning

Building upon previously reported high-entropy alloys (HEA) in the Ti−V−Cr−Nb−Mo/Fe compositional space, this study presents a novel bcc alloy containing 6 at. % of the earthabundant element Fe, exhibiting enhanced hydrogen storage properties at room temperature (25 °C). The optimized bcc Ti 23 V 30 Nb 10 Cr 31 Fe 6 alloy rapidly absorbs hydrogen at room temperature with a maximum gravimetric capacity of 3.4 wt % and an enthalpy of hydride formation of −34 kJ/mol H 2 . This value is among the lowest reported to date for bcc HEAs. The reversible capacity of this alloy remains stable at 2.0 wt % upon 20 absorption/desorption cycles at room temperature, which is among the highest reported values for this material class. Collectively, these advances highlight the potential of Fecontaining multielement alloys as efficient and economical candidates for next-generation hydrogen storage systems.

08 HYDROGEN↗
DOE OSTI2022

Process Model-Based Validation of the Intensification of Biomass Fast Pyrolysis in a Fluidized Bed via Autothermal Operation

A model for the fluidized-bed pyrolysis of biomass is extended to enable the simulation of intensified autothermal operation. In this system, partial oxidation of char and pyrolysate species replaces an external heat source to supply the enthalpy of pyrolysis, greatly increasing the throughput of a given fast pyrolysis reactor. Oxidation reactions are compiled from CRECK and other literature sources and extended to cover the species found in pyrolysis products, including a derivation of the catalytic effect of ash on char combustion from experimental studies at Iowa State University (ISU). Results indicate a roughly 3-fold increase in biomass throughput for a given reactor volume at the pilot scale with minimal loss of valuable products, in agreement with published data from the 3 × 10 -2 m 3 reactor operated by ISU. Additionally, a proposal for a 250 ton/day commercial-scale biomass pyrolysis reactor is analyzed, showing that the same size reactor operated by heating the fluidizing gas externally would have only one-tenth the capacity; if operated by heating the sand externally, the sand would need to circulate at an impractical 0.7 bed volume per minute rate to maintain a 250 ton/day biomass throughput.

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

Unexpected “Spontaneous” Evolution of Catalytic, MOF-Supported Single Cu(II) Cations to Catalytic, MOF-Supported Cu(0) Nanoparticles

A desirable feature of metal–organic frameworks (MOFs) is their well-defined structural periodicity and the presence of well-defined catalyst grafting sites (e.g., reactive –OH and –OH 2 groups) that can support single-site heterogeneous catalysts. However, one should not overlook the potential role of residual organic moieties, specifically formate ions that can occupy the catalyst anchoring sites during MOF synthesis. Here we show how these residual formate species in a Zr-based MOF, NU-1000, critically alter the structure, redox capability, and catalytic activity of postsynthetically incorporated Cu(II) ions. Single-crystal X-ray diffraction measurements established that there are two structurally distinct types of Cu(II) ions in NU-1000: one type with residual formate and one without. In NU-1000 with formate, Cu(II) solely binds to the node via the formate-unoccupied, bridging μ 3 –OH, whereas in the formate-free case, it displaces protons from two node hydroxo ligands and resides close to the terminal –OH 2 . Under an inert atmosphere, node-bound formate facilitates the unanticipated reduction of isolated Cu(II) to nanoparticulate Cu(0)—a behavior which is essentially absent in the formate-free analogue because no other sacrificial reductant is present. When the two MOFs were tested as benzyl alcohol oxidation catalysts, we observed that residual formate boosts the catalytic turnover frequency. Density functional calculations showed that node-bound formate acts as a sacrificial two-electron donor and assists in reducing Cu(II) to Cu(0) by a nonradical pathway. The negative Gibbs free energy of reaction (ΔG) and enthalpy of reaction (ΔH) indicate that the reduction is thermodynamically favorable. Here, the work presented here highlights how the often-neglected residual formate prevalent in nearly all zirconium-based MOFs can significantly modulate the properties of supported catalysts.

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

Ga + -Chabazite Zeolite: A Highly Selective Catalyst for Nonoxidative Propane Dehydrogenation

Ga-chabazite zeolites (Ga-CHA) have been found to efficiently catalyze propane dehydrogenation with high propylene selectivity (96%). In situ FTIR spectroscopy and pulse titrations are employed to determine that upon reduction, surface Ga 2 O 3 is reduced and diffuses into the zeolite pores, displacing the Brønsted acid sites (BAS) and forming extra-framework Ga + sites. This isolated Ga + site reacts reversibly with H 2 to form GaHx (2034 cm -1 ) with an enthalpy of formation of ~ -51.2 kJ·mol -1 , a result supported by Density functional theory (DFT) calculations. The initial C 3 H 6 dehydrogenation rates decrease rapidly (40%) during the first 100 min and then decline slowly afterward, while the C 3 H 6 selectivity is stable at ~ 96%. The reduction in the reaction rate is correlated with the formation of polycyclic aromatics inside the zeolite (using UV-vis spectroscopy) indicating that the accumulation of polycyclic aromatics is the main cause of the deactivation. The carbon species formed can be easily oxidized at 600 °C with complete recovery of the PDH catalytic properties. The correlations between GaH x vs. Ga/Al ratio, and PDH rates vs. Ga/Al ratio show that extra-framework Ga + is the active center catalyzing propane dehydrogenation. The higher reaction rate on Ga + than In + in CHA zeolites, by a factor of 43, is the result of differences in the stabilization of the transition state due to the higher stability of Ga 3+ vs. In 3+ . The uniformity of the Ga + sites in this material makes it an excellent model for the molecular understanding of metal cation exchanged hydrocarbon interactions in zeolites.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Influence of 1-Butene Adsorption on the Dimerization Activity of Single Metal Cations on UiO-66 Nodes

Grafting metal cations to missing linker defect sites in zirconium-based metal–organic frameworks, such as UiO-66, produces a uniquely well-defined and homotopic catalytically active site. We present here the synthesis and characterization of a group of UiO-66-supported metal catalysts, M-UiO-66 (M = Ni, Co, Cu, and Cr), for the catalytic dimerization of alkenes. The hydrogen–deuterium exchange via deuterium oxide adsorption followed by infrared spectroscopy showed that the last molecular water ligand desorbs from the sites after evacuation at 300 °C leading to M(OH)-UiO-66 structures. Adsorption of 1-butene is studied using calorimetry and density functional theory techniques to characterize the interactions of the alkene with metal cation sites that are found active for alkene oligomerization. For the most active Ni-UiO-66, the removal of molecular water from the active site significantly increases the 1-butene adsorption enthalpy and almost doubles the catalytic activity for 1-butene dimerization in comparison to the presence of water ligands. Other M-UiO-66 (M = Co, Cu, and Cr) exhibit 1–3 orders of magnitude lower catalytic activities compared to Ni-UiO-66. The catalytic activities correlate linearly with the Gibbs free energy of 1-butene adsorption. Furthermore, density functional theory calculations probing the Cossee–Arlman mechanism for all metals support the differences in activity, providing a molecular level understanding of the metal site as the active center for 1-butene dimerization.

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

Toward Hydrogen Isotope Separations through Strong Hydrogen Adsorption at Open Copper(I) Sites in an Ultramicroporous Metal-Organic Framework

Metal-organic frameworks with coordinatively unsaturated metal sites (open metal sites) capable of engaging in orbital interactions with pi-acidic gases are of interest for enabling ambient-temperature gas separations, such as hydrogen isotope separations. In view of the weakly pi-acidic nature of H2, we sought to strengthen pi-backbonding-mediated H2 adsorption through pore confinement effects. Toward that end, we synthesized and characterized the ultramicroporous metal-organic framework CuxZn5-xCl4-yHz(bbta)3 (CuIZn-MFU-4; H2bbta = 1H,5H-benzo(1,2-d:4,5-d')bistriazole), featuring pi-basic trigonal pyramidal CuI sites that reside within 7 A of one another at their closest. Gas adsorption measurements reveal an H2 adsorption enthalpy of -38 kJ/mol, exceeding that of the larger-pore analog (CuIZn-MFU-4l; -33 kJ/mol) and representing the strongest H2 adsorption yet achieved in a metal-organic framework. The stronger H2 adsorption in CuIZn-MFU-4 is attributed to a combination of pore confinement effects and the increased ..sigma..-accepting nature of the CuI sites caused by a more electron-withdrawing bbta2- linker, as supported by structural, spectroscopic, and computational evidence. With the strongest H2 adsorption, equilibrium isotope effects in CuIZn-MFU-4 lead to a D2/H2 selectivity (as estimated by ideal adsorbed solution theory) of 1.35 even at 298 K, approaching the values reported below 200 K for conventional porous materials.

08 HYDROGEN↗
DOE OSTI2020

Inorganic Halide Double Perovskites with Optoelectronic Properties Modulated by Sublattice Mixing

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs 2 BB'Cl 6 ) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (t) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs 2 BB'Cl 6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy <0.05 eV/atom). Of these 261 likely synthesizable compounds, 47 contain no toxic elements and have direct or nearly direct (within 100 meV) band gaps between 1 and 3 eV, as computed with hybrid density functional theory (HSE06). Within this set, we identify the triple-alkali perovskites Cs 2 [Alk] + [TM] 3+ Cl 6 , where Alk is a group 1 alkali cation and TM is a transition-metal cation, as a class of Cs 2 BB'Cl 6 double perovskites with remarkable optical properties, including large and tunable exciton binding energies as computed by the GW-Bethe-Salpeter equation (GW-BSE) method. We attribute the unusual electronic structure of these compounds to the mixing of the Alk-Cl and TM-Cl sublattices, leading to materials with small band gaps, large exciton binding energies, and absorption spectra that are strongly influenced by the identity of the transition metal. The role of the double-perovskite structure in enabling these unique properties is probed through an analysis of the electronic structures and chemical bonding of these compounds in comparison with other transition-metal and alkali transition-metal halides.

14 SOLAR ENERGY↗
DOE OSTI2020

Correlating Mechanical Sensitivity with Spin Transition in the Explosive Spin Crossover Complex [Fe(Htrz) 3 ] n [ClO 4 ] 2 n

Spin crossover complexes are known to undergo bond length, volume, and enthalpy changes during spin transition. In an explosive spin crossover complex, these changes could affect the mechanical and initiation sensitivity of the explosive and lead to the development of a new class of sensitivity switchable materials. To explore this relationship, the well-known spin crossover compound [Fe(Htrz) 3 ] n [ClO 4 ] 2n (1) was re-evaluated for its explosive properties, and its mechanical impact sensitivity was correlated to spin transition. A variable temperature impact test was developed and used to evaluate the impact sensitivity of 1 in the low spin (LS, S = 0), thermally accessed high spin (HS, S = 2), and mixed LS and HS states. For comparison, the structurally similar Ni compound, [Ni(Htrz) 3 ] n [ClO 4 ] 2n (2), which does not undergo a spin transition at accessible temperatures, was synthesized and characterized, and its explosive properties and variable temperature impact sensitivity measured. These results reveal a correlation between impact sensitivity and spin transition, where 1 exhibits lower impact sensitivity in the LS state and increases in sensitivity upon transition to the HS state. Here, density functional theory was used to predict structural changes that occur upon spin transition that correlate to the change in sensitivity. This demonstrates, for the first time, an explosive spin crossover compound (ExSCO) that exhibits switchable impact sensitivity with a fully reversible internal switching mechanism.

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

Isopotential Electron Titration of Ammonia Charge Transfer on Metal Catalysts

Electron transfer between adsorbates and surfaces determines the binding strength and reactivity of chemical moieties on materials designed for separations and catalysis. To quantify electron exchange, the extent of charge transfer resulting from ammonia adsorption on a Ru surface was measured by isopotential electron titration (IET) on a Ru catalytic condenser, where isopotential conditions were maintained between Ru and silicon separated by an insulating HfO 2 layer during gas phase ammonia adsorption. Charge transfer upon ammonia adsorption on a Ru catalytic condenser increased from 40 to 1200 nC/cm 2 at 75 and 225 °C, respectively. Charge transfer measurements provided a direct estimate of ammonia adsorption thermodynamics on Ru without knowing surface coverages a priori, revealing an adsorption enthalpy of −53 ± 10 kJ/mol and entropy of −61 ± 26 J/mol·K. Combining experimentally-measured charge transfer with kinetic Monte Carlo simulations informed adsorbate surface coverages determined that 0.058 electrons were transferred to the Ru surface for each molecular ammonia adsorption event (δ NH 3 = 0.058 ± 0.005 e – /NH 3 *), consistent with calculated Bader charges. The ability to measure the extent of charge transfer for adsorbed species provides a fundamental descriptor to understand existing and new chemically functional surfaces, providing a foundational method for the emerging field of thermochemical surface coulometry.

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

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗
DOE OSTI2022

Reconciling and Improving Formulations for Thermodynamics and Conservation Principles in Earth System Models (ESMs)

Abstract This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Learning grain boundary segregation energy spectra in polycrystals

The segregation of solute atoms at grain boundaries (GBs) can profoundly impact the structural properties of metallic alloys, and induce effects that range from strengthening to embrittlement. And, though known to be anisotropic, there is a limited understanding of the variation of solute segregation tendencies across the full, multidimensional GB space, which is critically important in polycrystals where much of that space is represented. Here we develop a machine learning framework that can accurately predict the segregation tendency—quantified by the segregation enthalpy spectrum—of solute atoms at GB sites in polycrystals, based solely on the undecorated (pre-segregation) local atomic environment of such sites. We proceed to use the learning framework to scan across the alloy space, and build an extensive database of segregation energy spectra for more than 250 metal-based binary alloys. The resulting machine learning models and segregation database are key to unlocking the full potential of GB segregation as an alloy design tool, and enable the design of microstructures that maximize the useful impacts of segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Confinement effects and acid strength in zeolites

Chemical reactivity and sorption in zeolites are coupled to confinement and—to a lesser extent—to the acid strength of Brønsted acid sites (BAS). In presence of water the zeolite Brønsted acid sites eventually convert into hydronium ions. The gradual transition from zeolite Brønsted acid sites to hydronium ions in zeolites of varying pore size is examined by ab initio molecular dynamics combined with enhanced sampling based on Well-Tempered Metadynamics and a recently developed set of collective variables. While at low water content (1–2 water/BAS) the acidic protons prefer to be shared between zeolites and water, higher water contents (n > 2) invariably lead to solvation of the protons within a localized water cluster adjacent to the BAS. At low water loadings the standard free energy of the formed complexes is dominated by enthalpy and is associated with the acid strength of the BAS and the space around the site. Conversely, the entropy increases linearly with the concentration of waters in the pores, favors proton solvation and is independent of the pore size/shape.

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

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease

Emerging SARS-CoV-2 variants continue to threaten the effectiveness of COVID-19 vaccines, and small-molecule antivirals can provide an important therapeutic treatment option. The viral main protease (Mpro) is critical for virus replication and thus is considered an attractive drug target. We performed the design and characterization of three covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2 created by splicing components of hepatitis C protease inhibitors boceprevir and narlaprevir, and known SARS-CoV-1 protease inhibitors. A joint X-ray/neutron structure of the M pro /BBH-1 complex demonstrates that a Cys145 thiolate reaction with the inhibitor’s keto-warhead creates a negatively charged oxyanion. Protonation states of the ionizable residues in the M pro active site adapt to the inhibitor, which appears to be an intrinsic property of M pro . Structural comparisons of the hybrid inhibitors with PF-07321332 reveal unconventional F···O interactions of PF-07321332 with Mpro which may explain its more favorable enthalpy of binding. BBH-1, BBH-2 and NBH-2 exhibit comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Releasing chemical energy in spatiallyprogrammed ferroelectrics

Chemical energy ferroelectrics are generally solid macromolecules showing spontaneous polarization and chemical bonding energy. These materials still suffer drawbacks, including the limited control of energy release rate, and thermal decomposition energy well below total chemical energy. To overcome these drawbacks, we report the integrated molecular ferroelectric and energetic material from machine learning-directed additive manufacturing coupled with the ice-templating assembly. The resultant aligned porous architecture shows a low density of 0.35 g cm -3 , polarization-controlled energy release, and an anisotropic thermal conductivity ratio of 15. Thermal analysis suggests that the chlorine radicals react with macromolecules enabling a large exothermic enthalpy of reaction (6180 kJ kg -1 ). In addition, the estimated detonation velocity of molecular ferroelectrics can be tuned from 6.69 ± 0.21 to 7.79 ± 0.25 km s -1 by switching the polarization state. These results provide a pathway toward spatially programmed energetic ferroelectrics for controlled energy release rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗