Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Pareto solutions in multicriteria optimization under uncertainty

We present and analyze several definitions of Pareto optimality for multicriteria optimization or decision problems with uncertainty primarily in their objective function values. In comparison to related notions of Pareto robustness, we first provide a full characterization of an alternative efficient set hierarchy that is based on six different ordering relations both with respect to the multiple objectives and a possibly finite, countably infinite or uncountable number of scenarios. We then establish several scalarization results for the generation of the corresponding efficient points using generalized weighted-sum and epsilon-constraint techniques. Finally, we leverage these scalarization results to also derive more general conditions for the existence of efficient points in each of the corresponding optimality classes, under suitable assumptions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Private, public, and bottled drinking water: Shared contaminant-mixture exposures and effects challenge

Background: Humans are primary drivers of environmental–contaminant exposures worldwide, including in drinking-water (DW). In the United States, point-of-use DW (POU–DW) is supplied via private tapwater (TW), public-supply TW, and bottled water (BW). Differences in management, monitoring, and messaging and lack of directly–intercomparable exposure data influence the actual and perceived quality and safety of different DW supplies and directly impact consumer decision–making. Objectives: The purpose of this paper is to provide a meta-analysis (quantitative synthesis) of POU–DW contaminant–mixture exposures and corresponding potential human–health effects of private-TW, public-TW, and BW by aggregating exposure results and harmonizing apical–health–benchmark–weighted and bioactivity–weighted effects predictions across previous studies by this research group. Discussion: Simultaneous exposures to multiple inorganic and organic contaminants of known or suspected human-health concern are common across all three DW supplies, with substantial variability observed in each and no systematic difference in predicted cumulative risk between supplies. Differences in contaminant or contaminant–class exposures, with important implications for DW–quality improvements, were observed and attributed to corresponding differences in regulation and compliance monitoring. Conclusion: The results indicate that human-health risks from contaminant exposures are common to and comparable in all three DW–supplies, including BW. Importantly, this study’s target analytical coverage, which exceeds that currently feasible for water purveyors or homeowners, nevertheless is a substantial underestimation of the breadth of contaminant mixtures in the environment and potentially present in DW. Thus, the results emphasize the need for improved understanding of the adverse human-health implications of long-term exposures to low–level inorganic–/organic–contaminant mixtures across all three distribution pipelines and do not support commercial messaging of BW as a systematically safer alternative to public-TW. Regardless of the supply, increased public engagement in source-water protection and drinking–water treatment is necessary to reduce risks associated with long-term DW–contaminant exposures, especially in vulnerable populations, and to reduce environmental waste and plastics contamination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗
DOE OSTI2020

Superior performance biodiesel from biomass-derived fusel alcohols and low grade oils: Fatty acid fusel esters (FAFE)

We demonstrated production of a superior performance biodiesel referred to here as fatty acid fusel alcohol esters (FAFE) – by reacting fusel alcohols (isobutanol, 3-methyl-1-butanol, and (S)-(-)-2-methyl-1-butanol) with oil (glyceryl trioleate) using lipase from Aspergillus oryzae. Reaction conditions corresponding to a molar ratio of 5:1 (fusel alcohols to oil), enzyme loading of 2% w/w, reaction temperature of 35 °C, shaking speed of 250 rpm, and reaction time of 24 h achieved >97% conversion to FAFE. Further, FAFE obtained from reacting a fusel alcohol mixture with corn oil were evaluated for use as a fuel for diesel engines. FAFE mixtures showed superior combustion and cold-flow properties, with the derived cetane numbers up to 4.8 points higher, cloud points up to -6 °C lower, and the heat of combustion up to 2.1% higher than the corresponding FAME samples, depending on the fusel mixture used. This represents a significant improvement for all three metrics, which are typically anti-correlated. Finally, FAFE provides a new opportunity for expanded usage of biodiesel by addressing feedstock limitations, fuel performance, and low temperature tolerance.

09 BIOMASS FUELS↗
DOE OSTI2020

A deep insight on the coal ash-to-slag transformation behavior during the entrained flow gasification process

Recent research provided deep insight on the coal ash-to-slag transformation characterization during the entrained flow gasification process, with experimentation on a 40 kg/hr (dry basis) coal-fed opposed multi-burner (OMB) entrained flow gasifier and simulation via FactSage™ software. A general mechanism is presented to relate the gasifier design temperature, ash fluid temperature, and operating temperature with the degree of the slag polymerization. The change of the high temperature zone, the corresponding particle residence time in the high temperature zone, and syngas composition have obvious effects on the slag mineral transformation behavior. Mineral types formed on the wall of the gasifier chamber were mainly anorthite (CaAl 2 Si 2 O 8 ), aluminum oxide (Al 2 O 3 ), and calcium sulfide (CaS). These minerals transformed to anorthite and diopside (CaMgSi 2 O 6 ) at the slag hole zone, while the minerals at the lock hopper were anorthite, orthoclase (KAlSi 3 O 8 ), quartz (SiO 2 ), gypsum (CaSO 4 ), calcite (CaCO 3 ), and halite (NaCl). FactSage™ predicted minerals as anorthite, diopside, orthoclase, and albite (NaAlSi 3 O 8 ), etc., where the slag temperature was below the ash fluid temperature and when the ratios of CO/CO 2 and (CO + H 2 )/CO 2 were lower than 1.0 and 2.0, respectively. By simulation, residual carbon was found to be the dominant factor over syngas composition to cause mineral transformation, and this was verified experimentally. The Ca-based crystals, typically anorthite, was shifted to diopside, near the slag hole zone, and a linear relationship was found between the content ratios of diopside/(anorthite + diopside), CaO/SiO 2 , and (CaO + MgO)/SiO 2 . A dimensionless number, θ, was defined to characterize the changing chemical composition and the degree of slag polymerization, with temperature deviation from the design condition. Three zones of θ were identified and related to the deviation between the actual gasification condition from the design condition. Finally, a low slag polymerization degree corresponded with a higher temperature deviation between the actual condition and design condition, and this proved that increased residual carbon content and changing iron valence state increased the mineral types when the slag temperature was below the ash fluid temperature.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Branches, quivers, and ideals for knot complements

We generalize the F K invariant, i.e. Z^ for the complement of a knot K in the 3-sphere, the knots-quivers correspondence, and A-polynomials of knots, and find several interconnections between them. We associate an F K invariant to any branch of the A-polynomial of K and we work out explicit expressions for several simple knots. We show that these F K invariants can be written in the form of a quiver generating series, in analogy with the knots-quivers correspondence. We discuss various methods to obtain such quiver representations, among others using R-matrices. We generalize the quantum a-deformed A-polynomial to an ideal that contains the recursion relation in the group rank, i.e. in the parameter a, and describe its classical limit in terms of the Coulomb branch of a 3d-5d theory. We also provide t-deformed versions. Furthermore, we study how the quiver formulation for closed 3-manifolds obtained by surgery leads to the superpotential of 3d N = 2 theory T[M 3 ] and to the data of the associated modular tensor category MTC[M 3 ].

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A note on the instability and pattern formation of shrinkage cracks in viscoplastic soils

In this work we present a theoretical study on the conditions for the onset of cracks, as well as the corresponding pattern formation, in saturated viscoplastic soils under isotropic loading (extension). The type of stress applied is left unspecified, to cover a variety of loadings including shrinkage due to desiccation, isotropic thermal expansion, mechanical loading and so forth. By treating the saturated soil as rigid viscoplastic, we obtain a 2D extension of the Cnoidal Waves equations (Veveakis and Regenauer-Lieb, 2015). By numerically solving the corresponding boundary value problem, we retrieve conditions for the onset of cracking instability in 2D loading, and identify the characteristic spacing between cracks to be a length scale combining all the hydro-mechanical parameters of the problem. Finally, we show that in a rectangular slab of clay under isotropic extension, patterns of triangular, rectangular and hexagonal cracks can tessellate the domain, with the hexagonal pattern being the energetically favored, as it minimizes the free energy of the system.

58 GEOSCIENCES↗
DOE OSTI2023

Measured CO 2 sorption isotherms with 25 Bakken Petroleum System rock samples from the Lower and Upper Shales, Middle Bakken, and Three Forks formations

Isotherms were measured with 25 Bakken Petroleum System (BPS) rock samples using a magnetic suspension balance at reservoir conditions of 110 °C and pressures up to 345 bar. Samples from four wells represented the Middle Bakken (MB) and Three Forks (TF) production zones, and the Upper Bakken Shale (UBS) and Lower Bakken Shale (LBS) source shales. UBS and LBS CO 2 sorption capacities at 345 bar averaged 12–14 g/kg rock compared to 1.3 to 2.5 g/kg rock for MB and TF. Increasing pressure from 150 to 345 bar resulted in only small increases in CO 2 sorption for all lithofacies. For the UBS and LBS source shales, the higher average total organic carbon (TOC) values (14.7 and 11.4 wt.%) and clay values (25 and 28 wt.%) corresponded to their higher CO 2 sorption values while for the MB and TF samples, lower TOCs (0.3 and 0.5 wt.%) and clay values (8.2 and 18.2 wt.%) corresponded to lower CO 2 sorption values. Here, when all 25 samples (from all lithofacies) were evaluated as a group, sorption capacities were highly correlated with TOC (Spearman's rho values +0.84 to +0.85), with rock density (rho values -0.75 to -0.76), and clay content (rho values +0.66). However, correlations with porosity and thermal maturity (Tmax) were weak. Measured sorption capacities and rock porosities showed that sorption accounted for about 10-times as much CO 2 storage as the available pore space for the UBS and LBS lithofacies. In contrast, the pore space dominated the storage resource of the TF and MB lithofacies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-temperature stability and thermal expansion behavior of equi-atomic refractory multi-principal element alloys based on MoNbTi system for Gen IV reactor applications

The present study investigates the thermal stability and thermal expansion behavior of seven equi-atomic refractory multi-principal element alloys (MPEAs) based on the MoNbTi ternary system composed of low neutron absorption cross section elements. Through an integrated approach utilizing in-situ high-temperature X-ray diffraction (HT-XRD) in conjunction with differential scanning calorimetry (DSC), dilatometry and ageing heat treatment, the thermal stability of the MPEAs was comprehensively analyzed. In-situ HT-XRD experiment confirmed the stability of the room temperature phases up to 1000 °C with no peaks observed corresponding to additional phases in the HT-XRD patterns at 500, 800 and 1100 °C. DSC thermograms showed the absence of peaks up to 1000 °C, while peaks and valleys corresponding to exothermic and endothermic events were observed above 1000 °C. Coefficient of thermal expansion (CTE) derived from second order polynomial fitting of linear thermal expansion data from the dilatometry experiment showed linear increment up to 1000 °C for all the alloys except those containing Zr. The Cr containing alloys exhibited notably higher CTE values, particularly the Al containing alloy exhibited the highest value. Ageing heat treatment at 800 and 1000 °C for 96 h and subsequent microstructural analysis revealed significant precipitation of secondary phases in MoNbTiZr, MoNbTiZrV and MoNbTiCrAl. In conclusion, a substantial increase in hardness was observed in MoNbTiZr and MoNbTiCrAl due to secondary phase precipitation, while the other alloys maintained hardness values comparable to their as-cast and homogenized states.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

DOME: Directional medical embedding vectors from Electronic Health Records

Motivation: The increasing availability of Electronic Health Record (EHR) systems has created enormous potential for translational research. Recent developments in representation learning techniques have led to effective large-scale representations of EHR concepts along with knowledge graphs that empower downstream EHR studies. However, most existing methods require training with patient-level data, limiting their abilities to expand the training with multi-institutional EHR data. On the other hand, scalable approaches that only require summary-level data do not incorporate temporal dependencies between concepts. Methods: We introduce a DirectiOnal Medical Embedding (DOME) algorithm to encode temporally directional relationships between medical concepts, using summary-level EHR data. Specifically, DOME first aggregates patient-level EHR data into an asymmetric co-occurrence matrix. Then it computes two Positive Pointwise Mutual Information (PPMI) matrices to correspondingly encode the pairwise prior and posterior dependencies between medical concepts. Following that, a joint matrix factorization is performed on the two PPMI matrices, which results in three vectors for each concept: a semantic embedding and two directional context embeddings. They collectively provide a comprehensive depiction of the temporal relationship between EHR concepts. Results: We highlight the advantages and translational potential of DOME through three sets of validation studies. First, DOME consistently improves existing direction-agnostic embedding vectors for disease risk prediction in several diseases, for example achieving a relative gain of 5.5% in the area under the receiver operating characteristic (AUROC) for lung cancer. Second, DOME excels in directional drug-disease relationship inference by successfully differentiating between drug side effects and indications, correspondingly achieving relative AUROC gain over the state-of-the-art methods by 10.8% and 6.6%. Finally, DOME effectively constructs directional knowledge graphs, which distinguish disease risk factors from comorbidities, thereby revealing disease progression trajectories. The source codes are provided at https://github.com/celehs/Directional-EHRembedding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Vibrational sum frequency generation spectroscopy reveals the inertness of chromium oxide (001) surfaces

Nanoengineered metal oxides such as Cr(III)-oxide (chromia) films have diverse potential applications in corrosion inhibition, remediation, energy generation, catalysis, data storage, and biological and environmental systems. Concerns about material degradation or oxidation to toxic chromate necessitate an understanding of chromia/aqueous interfaces, beginning with their hydroxylation and hydration behavior. Vibrational sum-frequency generation spectroscopy (vSFG) provides specific molecular-level information about water at the oxide/aqueous junction with high surface selectivity. To overcome the strong absorber problem typical of certain metal oxides in the UV–visible range during nonlinear optical studies, we employed molecular beam epitaxy to deposit transparent, epitaxial nanofilms of Cr 2 O 3 with (001) crystalline orientation on sapphire (Al 2 O 3 (001)) substrates, as confirmed by atomic force microscopy and X-ray diffraction. vSFG spectra of the air and water interfaces of the Cr 2 O 3 (001) films reveal hydroxyl features corresponding to both dissociated and molecular water on the surface. In contrast to the dangling hydroxyls found on the bare Al 2 O 3 (001) substrate, the hydroxyl groups on the deposited Cr 2 O 3 (001) nanofilm do not readily undergo isotopic H/D exchange when exposed to varying forms of D 2 O under ambient conditions. When considering the chemistries of the corresponding trivalent cations in aqueous solution, the finding is at variance with their similar acidities and proton exchange dynamics but consistent with markedly slower inner-sphere water ligand exchange of hexaquo Cr 3+ . This finding challenges the idea that proton and water exchange at the water/Cr₂O₃(001) interface are solely correlated, driven by the strength of metal-bridging oxygen bonds and surface hydroxyl distribution, without direct causation.

Boamah-Agyemang, Mavis D. [Pacific Northwest Natio↗
DOE OSTI2023

Influence of urban extent discrepancy on the estimation of surface urban heat island intensity: A global-scale assessment in 892 cities

The estimation of surface urban heat island intensity (SUHII) is crucial for studying the urban thermal environment, which is influenced not only by the commonly known definition of rural reference but also by the delineation of urban extent. Existing studies relies on various urban extent products defined in different ways, and the influence of urban extent discrepancy (UED) on SUHII estimates still remains unclear. In this study, we collected five open-source global urban extent products (GUEPs) for the year 2015 and corresponding daily land surface temperature (LST) observations (MYD11A1). Based on these products, we quantified the UED-induced uncertainty in SUHII estimates by comparing absolute difference (ΔSUHII AD ) and relative difference (ΔSUHII RD ) in SUHII among GUEPs across 892 global cities. Additionally, we introduced an ISF-constrained (ISF–C) method to reduce SUHII differences among GUEPs by constraining the impervious surface fraction (ISF) within urban and rural extents. The results show that urban extents delineated by different GUEPs are not consistent, leading to their difference in ISF as well as LST, which in turn causes uncertainties in the estimated SUHII. On average for global cities, the annual daytime and nighttime ΔSUHII AD are 0.46 ± 0.02 °C (mean ± 95% confidence interval) and 0.24 ± 0.01 °C, with corresponding ΔSUHII RD of 42.0 ± 2.7% and 35.2 ± 2.3%, respectively. The UED-induced uncertainty in SUHII estimates varies among climate zones, and the annual daytime ΔSUHII RD averaged for cities located in the arid zone reaches up to 60.8 ± 6.6%, which is nearly twice as high as that in other climate zones. More importantly, both ΔSUHII AD and ΔSUHII RD show lower values when using the ISF-C method, implying the effectiveness of this method. This study highlights the non-negligible impact of UED on the estimation of SUHII, which requires more attention due to the inconsistency of urban extents among current products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

DPM: A deep learning PDE augmentation method with application to large-eddy simulation

A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Accelerating eigenvalue computation for nuclear structure calculations via perturbative corrections

Subspace projection methods utilizing perturbative corrections have been proposed for computing the lowest few eigenvalues and corresponding eigenvectors of large Hamiltonian matrices. In this paper, we build upon these methods and introduce the term Subspace Projection with Perturbative Corrections (SPPC) method to refer to this approach. We tailor the SPPC for nuclear many-body Hamiltonians represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). We use the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hamiltonian represented in a small configuration space, whereas the second is viewed as the perturbation to the first matrix. Eigenvalues and eigenvectors of the first matrix can be computed efficiently. Because of the split, perturbative corrections to the eigenvectors of the first matrix can be obtained efficiently from the solutions of a sequence of linear systems of equations defined in the small configuration space. These correction vectors can be combined with the approximate eigenvectors of the first matrix to construct a subspace from which more accurate approximations of the desired eigenpairs can be obtained. We show by numerical examples that the SPPC method can be more efficient than conventional iterative methods for solving large-scale eigenvalue problems such as the Lanczos, block Lanczos and the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The method can also be combined with other methods to avoid convergence stagnation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗
DOE OSTI2020

Direct processing of PbZr 0.53 Ti 0.47 O 3 films on glass and polymeric substrates

This work reports on direct crystallization of PbZr 0.53 Ti 0.47 O 3 (PZT) thin films on glass and polymeric substrates, using pulsed thermal processing (PTP). Specifically, xenon flash lamps deliver pulses of high intensity, short duration, broadband light to the surface of a chemical solution deposited thin film, resulting in the crystallization of the film. Structural analysis by X-ray diffraction (XRD) and transmission electron microscopy show the existence of perovskite structure in nano-sized grains (≤5 nm). Local functional analysis by band excitation piezoelectric spectroscopy and electrostatic force microscopy confirm the presence of a ferroelectric phase and retention of voltage-written polarization for multiple days. Based on structural and functional analyses, strategies are discussed for optimization of pulse voltage and duration for the realization of crystalline ferroelectric thin films. For ~200 nm-thick PZT films on glass substrates, 500 μs-long pulses were required for crystallization, starting with 100 pulses at 350 V, 10 or 25 pulses at 400 V and in general lower number of pulses at higher voltages (resulting in higher radiant energy). Overall power densities of >6.4 kW/cm 2 were needed for appearance of peaks corresponding to the perovskite phase in the XRD. Films on glass processed at 350–400 V had a higher degree of 111-oriented perovskite grains. Higher applied radiant energy (through increased pulse voltage or count) resulted in more random and/or partially 001-oriented films. For ~1 μm-thick PZT films on polymeric substrates, 10 to 25 250 μs-long pulses at voltages ranging between 200 to 250 V, corresponding to power densities of ~2.8 kW/cm 2 , were optimal for maximized perovskite phase crystallization, while avoiding substrate damage.

36 MATERIALS SCIENCE↗