Engineering Papers⌕ Search

Engineering topics

Wang, Bo

Publications and source records attributed to Wang, Bo.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unlocking Electrostrain in Plastically Deformed Barium Titanate

Achieving substantial electrostrain alongside a large effective piezoelectric strain coefficient (d 33 *) in piezoelectric materials remains a formidable challenge for advanced actuator applications. Here, in this work, a straightforward approach to enhance these properties by strategically designing the domain structure and controlling the domain switching through the introduction of arrays of ordered {100}<100> dislocations is proposed. This dislocation engineering yields an intrinsic lock-in steady–state electrostrain of 0.69% at a low field of 10 kV cm -1 without external stress and an output strain energy density of 5.24 J cm -3 in single-crystal BaTiO 3 , outperforming the benchmark piezoceramics and relaxor ferroelectric single-crystals. Additionally, applying a compression stress of 6 MPa fully unlocks electrostrains exceeding 1%, yielding a remarkable d 33 * value over 10 000 pm V -1 and achieving a record-high strain energy density of 11.67 J cm -3 . Optical and transmission electron microscopy, paired with laboratory and synchrotron X-ray diffraction, is employed to rationalize the observed electrostrain. Phase-field simulations further elucidate the impact of charged dislocations on domain nucleation and domain switching. These findings present an effective and sustainable strategy for developing high-performance, lead-free piezoelectric materials without the need for additional chemical elements, offering immense potential for actuator technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A brief review on strain engineering of ferroelectric K x Na 1- x NbO 3 epitaxial thin films: Insights from phase-field simulations

Strains play a pivotal role in determining the phase equilibrium, domain configuration, and functional properties of the low-dimensional ferroelectrics. There is growing interest in the strain engineering of ferroelectric K x Na 1- x NbO 3 (KNN) epitaxial thin films, which exhibit excellent physical properties and promise as eco-friendly alternatives to lead-based ferroelectrics for microdevice applications. Further, advances have been made in understanding the phase equilibria and transitions, domains and domain walls, and their relations to the physical properties of KNN epitaxial thin films using a combination of experiments and theoretical modeling, particularly phase-field simulations. Here, we review recent progress in these aspects and showcase the phase-field method for establishing strain phase diagrams, elucidating the domain and domain wall structures at equilibrium, and predicting the structure–property relationships in ferroelectric KNN thin films. We also discuss challenges and opportunities to further advance our understanding of KNN thin films and potentially unlock new functionalities by leveraging phase-field simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Strain phase equilibria and phase‐field method of ferroelectric polydomain: A case study of monoclinic K x Na 1 − x NbO 3 thin films

Abstract Knowledge of the thermodynamic equilibria and domain structures of ferroelectrics is critical to establishing their structure–property relationships that underpin their applications from piezoelectric devices to nonlinear optics. Here, we establish the strain condition for strain phase separation and polydomain formation and analytically predict the corresponding domain volume fractions and wall orientations of, relatively low symmetry and theoretically more challenging, monoclinic ferroelectric thin films by integrating thermodynamics of ferroelectrics, strain phase equilibria theory, microelasticity, and phase‐field method. Using monoclinic K x Na 1 − x NbO 3 (0.5 < x < 1.0) thin films as a model system, we establish the polydomain strain–strain phase diagrams, from which we identify two types of monoclinic polydomain structures. The analytically predicted strain conditions of formation, domain volume fractions, and domain wall orientations for the two polydomain structures are consistent with phase‐field simulations and in good agreement with experimental results in the literature. The present study demonstrates a general, powerful analytical theoretical framework to predict the strain phase equilibria and domain wall orientations of polydomain structures applicable to both high‐ and low‐symmetry ferroelectrics and provide fundamental insights into the equilibrium domain structures of ferroelectric K x Na 1 − x NbO 3 thin films that are of technology relevance for lead‐free dielectric and piezoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theory of nonlinear terahertz susceptibility in ferroelectrics

An analytical theory is developed for predicting the nonlinear susceptibility of ionic polarization to continuous electromagnetic waves in both bulk and strained thin film ferroelectrics. Using a perturbation method for solving the nonlinear equation of motion for ionic polarization within the framework of Landau-Ginzburg-Devonshire theory, the full second-order nonlinear susceptibility tensor is derived as a function of frequency, temperature, and strain. Here, the theory predicts the coexistence of a significantly enhanced second-order dielectric susceptibility and a relatively low dielectric loss in BaTiO 3 films with a strain-stabilized monoclinic ferroelectric phase and in a strained SrTiO 3 film near its temperature-driven second-order ferroelectric-to-paraelectric phase transition. In this paper, we establish a theoretical framework for predicting and exploiting nonlinear interactions between terahertz waves and ferroelectric materials and, more generally, suggest exciting opportunities to strain-engineer nonlinear dynamical properties of ferroelectrics beyond the static and quasistatic limits.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrical de-poling and re-poling of relaxor-PbTiO 3 piezoelectric single crystals without heat treatment

Re-poling of unexpected partially depoled piezoelectric materials conventionally needs to be first fully depoled through annealing above their Curie temperature to revive piezoelectric performances. Here, we investigated de-poling and re-poling of Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 single crystals under electric fields at room temperature. We found that alternating current electric fields with amplitudes near the coercive field at low frequencies (<10 Hz) can be employed to successfully depolarize poled crystals at room temperature. We also demonstrated a reversible polarization switching process with a relaxor-PbTiO 3 single crystal ultrasound transducer without device performance degradations. This experimental observation is supported by phase-field simulation, showing that alternating current electric fields can readily induce de-poling at room temperature, while direct current electric fields induce a transient depoled state only within an uncontrollable short period of time. The findings suggest new strategies for unprecedented in-device tailoring of the polarization states of ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrothermal mineralization of per- and polyfluoroalkyl substances for soil remediation

Abstract Per- and polyfluoroalkyl substances (PFAS) are persistent and bioaccumulative pollutants that can easily accumulate in soil, posing a threat to environment and human health. Current PFAS degradation processes often suffer from low efficiency, high energy and water consumption, or lack of generality. Here, we develop a rapid electrothermal mineralization (REM) process to remediate PFAS-contaminated soil. With environmentally compatible biochar as the conductive additive, the soil temperature increases to >1000 °C within seconds by current pulse input, converting PFAS to calcium fluoride with inherent calcium compounds in soil. This process is applicable for remediating various PFAS contaminants in soil, with high removal efficiencies ( >99%) and mineralization ratios ( >90%). While retaining soil particle size, composition, water infiltration rate, and cation exchange capacity, REM facilitates an increase of exchangeable nutrient supply and arthropod survival in soil, rendering it superior to the time-consuming calcination approach that severely degrades soil properties. REM is scaled up to remediate soil at two kilograms per batch and promising for large-scale, on-site soil remediation. Life-cycle assessment and techno-economic analysis demonstrate REM as an environmentally friendly and economic process, with a significant reduction of energy consumption, greenhouse gas emission, water consumption, and operation cost, when compared to existing soil remediation practices.

Science & Technology - Other Topics↗
DOE OSTI2024

Theory of strain phase equilibria and diagrams

Chemical phase equilibria and phase diagrams are well established and have served as indispensable guides for designing materials through chemical tuning. However, a general thermodynamic theory for strain phase equilibria remains elusive despite extensive studies which have revealed dramatic impacts of mechanical strain on the stability of phases and domain states in solid-state materials. Here, we establish the thermodynamic theory of strain equilibria and a general framework for efficiently constructing multidomain and multiphase diagrams of arbitrarily strained solids under incoherent conditions. As examples, we obtain temperature-strain phase diagrams of ferroelectric PbTiO 3 , strongly correlated VO 2 , and unconventional ferroelectric Hf 0.5 Zr 0.5 O 2 . We reveal the analogs of the Gibbs phase rules, multiple- and multi-critical points, common-tangent construction, and level rule in strain equilibria to those in the familiar chemical phase equilibria. Finally, our strain equilibria theory offers a powerful framework for predicting the thermodynamic stability of structural phases and domains in strained solids and for guiding the strain engineering of their functional properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase-field-guided design of record-high piezoelectricity and discovery of simultaneous high light transparency and high piezoelectricity in relaxor ferroelectrics

Here, the phase-field method has been extensively applied to predicting the domain structures and their responses to external fields and understanding experimentally observed domain states under different electromechanical conditions in ferroelectric heterostructures. This article highlights the successful examples of phase-field applications in guiding the design of relaxor ferroelectric ceramics and crystals with record-high piezoelectricity and the discovery of simultaneous high light transparency and piezoelectricity of relaxor ferroelectric crystals for a wide range of biomedical, robotics, and optoelectronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Primordial black hole mass functions as a probe of cosmic origin

Here, we discuss a novel window to probe the origin of our universe via the mass functions of primordial black holes (PBHs). The mass functions of PBHs are simply estimated using the conventional Press-Schechter formalism for two paradigms of cosmic origin, including inflationary ΛCDM and bounce cosmology. The standard inflationary ΛCDM model cannot generate an appreciable number of massive PBHs; however, non-trivial inflation models with blue-tilted power spectra at small scales and matter bounce cosmology provide formation mechanisms for heavy PBHs, which in turn, may seed the observed supermassive black holes (SMBHs). By fitting the SMBH mass functions at high redshift (z ~ 6) derived from Sloan Digital Sky Survey (SDSS) and Canada-France High-z Quasar Survey (CFHQS) quasars, for two paradigms of cosmic origin, we derive constraints on the PBH density fraction $f$ PBH at z ~ 6 and the characteristic mass M $\ast$ , with the prior assumption that all SMBHs stem from PBHs. We demonstrate that this newly proposed procedure, relying on astronomical measurements that utilize deep-field surveys of SMBHs at high redshift, can be used to constrain models of cosmic origin. Additionally, although not the main focus of this paper, we evolve the mass function from z ~ 6 to z ~ 0 through an assumption of 3 × 10 8 -year Eddington’s accretion, and give a rough estimation of $f$ PBH at z ~ 0.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Polarization near dislocation cores in SrTiO 3 single crystals: The role of flexoelectricity

Spontaneous polarization as large as ~28 μC/cm 2 was recently observed around the dislocation cores in non-polar SrTiO 3 bulk crystals, and its origin was attributed to the flexoelectric effect, i.e., polarization induced by strain gradients. However, the roles of flexoelectricity, relative to other electromechanical contributions, and the nature of dislocations, i.e., edge vs screw dislocations in the induced polarization, are not well understood. In this work, we study the role of flexoelectricity in inducing polarization around three types of dislocation cores in SrTiO 3 : b= a(100) edge dislocation, b= a(110) edge dislocation, and b= a(010) screw dislocation, where b is the Burgers vector. For the edge dislocations, polarization can be induced by electrostriction alone, while flexoelectricity is essential for stabilizing the symmetric polarization pattern. The shear component of the flexoelectric tensor has a dominant effect on the magnitude and spatial distribution of the flexoelectric polarization. In contrast, no polarization is induced around the b= a(010) screw dislocation through either electrostriction or flexoelectricity. Finally, our findings provide an in-depth understanding of the role of flexoelectricity in inducing polarization around dislocation cores and offer insights into the defect engineering of dielectric/ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, we theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Domain-dependent strain and stacking in two-dimensional van der Waals ferroelectrics

Van der Waals (vdW) ferroelectrics have attracted significant attention for their potential in next-generation nano-electronics. Two-dimensional (2D) group-IV monochalcogenides have emerged as a promising candidate due to their strong room temperature in-plane polarization down to a monolayer limit. However, their polarization is strongly coupled with the lattice strain and stacking orders, which impact their electronic properties. Here, we utilize four-dimensional scanning transmission electron microscopy (4D-STEM) to simultaneously probe the in-plane strain and out-of-plane stacking in vdW SnSe. Specifically, we observe large lattice strain up to 4% with a gradient across ~50 nm to compensate lattice mismatch at domain walls, mitigating defects initiation. Additionally, we discover the unusual ferroelectric-to-antiferroelectric domain walls stabilized by vdW force and may lead to anisotropic nonlinear optical responses. Our findings provide a comprehensive understanding of in-plane and out-of-plane structures affecting domain properties in vdW SnSe, laying the foundation for domain wall engineering in vdW ferroelectrics.

42 ENGINEERING↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Generalized analytical and numerical modeling of optical second harmonic generation in anisotropic crystals and complex heterostructures using #SHAARP package

Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. The closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, here we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP) for single interface (si) and in multilayers (ml) for homogeneous crystals. The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials. SHAARP.si and SHAARP.ml are available through GitHub https://github.com/Rui-Zu/SHAARP and https://github.com/bzw133/SHAARP.ml, respectively.

complex systems↗
DOE OSTI2023

Metabolic engineering of Synechococcus elongatus 7942 for enhanced sucrose biosynthesis

The capability of cyanobacteria to produce sucrose from CO 2 and light has a remarkable societal and biotechnological impact since sucrose can serve as a carbon and energy source for a variety of heterotrophic organisms and can be converted into value-added products. However, most metabolic engineering efforts have focused on understanding local pathway alterations that drive sucrose biosynthesis and secretion in cyanobacteria rather than analyzing the global flux re-routing that occurs following induction of sucrose production by salt stress. Here, we investigated global metabolic flux alterations in a sucrose-secreting (cscB-overexpressing) strain relative to its wild-type Synechococcus elongatus 7942 parental strain. Here we used targeted metabolomics, 13 C metabolic flux analysis (MFA), and genome-scale modeling (GSM) as complementary approaches to elucidate differences in cellular resource allocation by quantifying metabolic profiles of three cyanobacterial cultures - wild-type S. elongatus 7942 without salt stress (WT), wild-type with salt stress (WT/NaCl), and the cscB-overexpressing strain with salt stress (cscB/NaCl) - all under photoautotrophic conditions. We quantified the substantial rewiring of metabolic fluxes in WT/NaCl and cscB/NaCl cultures relative to WT and identified a metabolic bottleneck limiting carbon fixation and sucrose biosynthesis. This bottleneck was subsequently mitigated through heterologous overexpression of glyceraldehyde-3-phosphate dehydrogenase in an engineered sucrose-secreting strain. Our study also demonstrates that combining 13 C-MFA and GSM is a useful strategy to both extend the coverage of MFA beyond central metabolism and to improve the accuracy of flux predictions provided by GSM.

09 BIOMASS FUELS↗
DOE OSTI2023

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗