Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic coupling”

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 19 records

Coupling computational thermodynamics with density-function-theory based calculations to design L12 precipitates in Fe Ni based alloys

Achieving a high-volume fraction of thermodynamically stable L12-type precipitates that are resistant to coarsening is of great importance for the development of low-cost FeNi based austenitic steels. With the aid of computational thermodynamics, this work designed two model alloys: Fe-37.4Ni-6.1Al-2.9Ti (FNAT) and Fe-45.2Ni-5.9Al-8.5Si (FNAS). Both alloys were designed to contain a similar amount of L1 2 precipitate in Fe-Ni matrix without forming other precipitates. Density-Function-Theory (DFT) calculation was coupled with computational thermodynamics to predict the critical radius at which the precipitates change shape from spherical to cuboidal. The calculation results suggest that critical radius for the FNAT alloy is about two orders of magnitude larger than that for the FeNiAlSi alloy. Phase stability and morphology of the L1 2 precipitates in these two alloys were experimentally investigated through X-ray diffraction, atom probe tomography, and scanning and transmission electron microscopy. The L1 2 precipitates in the FeNiAlSi system were found to be cuboidal and rod shaped, with much larger size than the spherical ones in the FeNiAlTi system, agreeing with the calculation results. This work suggested that coupling computational thermodynamics with DFT calculations can be reliably used to design L1 2 precipitates in FeNi based alloys.

36 MATERIALS SCIENCE↗

Methodology for Thermodynamic Analysis Coupled with Computational Fluid Dynamics Modeling for Casting a Novel Aluminum–Cerium Alloy

In this paper, we present the results of computational fluid dynamics (CFD) analysis to assess castability, porosity, molten metal fluidity and other technological properties of a new Al-Ce alloy. A thermodynamic analysis of a new Al-Ce alloy is done to obtain the physical properties as a function of temperature across both solid and liquid phases. These properties are then used to build the CFD of model the full casting process from initial pouring through to final solidification of the part, here an example of a heavy duty mor mount is used. While several excellent commercial CFD codes exists this process shows that CFD can be used to assess the capability of the alloy to properly fill the mold, as well as give predictions where scattered porosity or large-scale defects may occur in the casting. Further, unlike the commercially available software (e.g., ProCAST, SOLIDCast , MAGMASOFT®) the complex 3D-analysis of the stress / strain fields in cast parts is not performed at this time. However, the availability of free software for assessing the required thermodynamic / thermo-physical properties of new alloys (OpenCALPHAD) and CFD codes such as OpenFOAM® makes the developed option attractive and economical, especially for the analysis of new Al-Ce alloys, for which the available data does not exist.

36 MATERIALS SCIENCE↗

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis

We report that replica-exchange molecular dynamics (REMD) has been proven to efficiently improve the convergence of free energy perturbation (FEP) calculations involving considerable reorganization of their surrounding. We previously introduced the FEP/(λ,H)-REMD algorithm for ligand binding, in which replicas along the alchemical thermodynamic coupling axis λ were expanded as a series of Hamiltonian boosted replicas along a second axis to form a two-dimensional (2D) replica-exchange exchange map [Jiang, W.; Roux, B., J. Chem. Theory Comput. 2010, 6 (9), 2559-2565]. Aiming to achieve a similar performance at a lower computational cost, we propose here a modified version of this algorithm in which only the end-states along the alchemical axis are augmented by boosted replicas. The reduced FEP/(λ,H)-REMD method with one-dimensional (1D) unbiased alchemical thermodynamic coupling axis λ is implemented on the basis of generic multiple copy algorithm (MCA) module of the biomolecular simulation program NAMD. The flexible MCA framework of NAMD enables a user to design customized replica-exchange patterns through Tcl scripting in the context of a highly parallelized simulation program without touching the source code. Two Hamiltonian tempering boosting scheme were examined with the new algorithm: a first one based on potential energy rescaling of a pre-identified “solute”, and a second one via the introduction of flattening torsional free energy barriers. As two illustrative examples with reliable experiment data, the absolute binding free energies of pxylene and n-butylbenzene to the nonpolar cavity of the L99A mutant of T4 lysozyme were calculated. Lastly, the tests demonstrate that the new protocol efficiently enhances the sampling of torsional motions for backbone and side chains around the binding pocket and accelerates the convergence of the free energy computations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel Carbon (C)-Boron (B)-Nitrogen (N)-Containing H2 Storage Materials

The following summarizes the research conducted for DOE project DE-EE0005658 “Novel Carbon(C)-Boron(B)-Nitrogen(N)-Containing H2 Storage Materials”. This work focused in part on the continued study of two materials identified from the preceding project DE-FG360GO18143 (“Hydrogen Storage by Novel CBN Heterocycle Materials”) as lead candidates to meet the DOE technical targets for either vehicular or non-automotive hydrogen storage applications. Specifically, a room-temperature liquid, 3-methyl-1,2-cyclopentane (B), and a high H2 capacity solid, 1,2-BN-cyclohexane (J), were selected for further characterization and performance optimization. In addition to these compounds, the current project also aimed to prepare several new materials predicted to be disposed towards direct reversibility of H2 release and uptake, a feature deemed critical to achieving efficient recycling of spent fuel end products. To assist in the rational design of these and other next-generation materials, this project undertook to investigate the mechanism of hydrogen release from established compounds (mainly B and J) using a combined experimental/computational approach. Among this project’s signature accomplishments, the preliminary synthetic route to B was optimized for production on decagram scale. With such quantities of material available, its performance in powering an actual 30 W proton exchange membrane (PEM) fuel cell stack was tested and found to be identical to that of facility H2. Despite this positive proof-of-concept achievement, however, further consideration of neat B as a potential hydrogen storage material was abandoned due to evidence of thermal instability. Specifically, mass spectrometry-coupled thermogravimetric analysis (TGA-MS) revealed significant H2 release from B to initiate at 50 °C, well below the 60 °C minimum threshold set by the DOE. This result prompted a more extensive investigation in the decomposition mechanism of B vis-à-vis that of J, which exhibited in neat form a substantially higher onset temperature for spontaneous H2 release (70 °C). Solution-phase kinetic experiments using ReactIR established a second-order dependence for the initial loss of H2 from both B and J; Arrhenius analysis, however, revealed the activation barrier for this reaction was lower for B than for J, which presumably contributes to the diminished thermal stability of the former. On the basis of these and other experimental results, extensive computational efforts yielded a reasonable mechanistic model for the dehydrogenation of 1,2-BN-cycloalkane materials. While the prospect of neat B as a suitable hydrogen storage material was discarded, it was proposed that the combination of B with more thermally stable amine-borane-based materials might afford mixtures with improved properties. Indeed, when B was combined with ammonia borane (AB) in a 2:1 molar ratio, the two materials formed a liquid. More significantly, this mixture remained liquid even after complete dehydrogenation, thus establishing the potential for a single-phase fuel cycle. (In contrast, the dehydrogenation product of neat B is a low melting solid (mp = 28-30 °C).) Another advantage conferred by the blend formulation was a dramatic reduction in the amount of borazine produced by AB. Borazine is a well-known contaminant of H2 produced by the thermal decomposition of neat AB, and exerts deleterious effects on fuel cell performance. Residual gas analysis (RGA) of the gas stream generated from the B-AB blend, however, detected just 0.01% borazine content when a Pt-Ni nanoparticle dehydrogenation catalyst was used. In all ii then, the 2:1 B-AB blend marks a major achievement in the effort to develop a suitable liquid amine-borane hydrogen storage material, and merits further investigation into the optimization for practical adoption. Similar realization of the potential of J as a high % wt. H2 material required a method to dehydrogenate the carbonaceous components of the molecule without the use of a sacrificial hydrogen acceptor, as had been reported in the previous project. Ultimately, this reaction was achieved for a B,N-disubstituted BN-cyclohexene model substrate using a gas flow system with a fixed Pd/C catalyst bed. Considerable work remains, however, to translate these initial results into a general protocol for complete dehydrogenation of fully saturated BN-cycloalkane materials such as J. With concrete confirmation of the possibility to perform both BN and CC dehydrogenation on a single theoretical substrate, COMSOL modeling was used to evaluate the effects of thermodynamically coupling the two reactions. It was hypothesized that the heat generated from exothermic BN dehydrogenation would partially drive the endothermic CC dehydrogenation reaction; this additional heat consumption was expected to in turn confer the benefit of lowering the maximum reactor temperature. A two-dimensional model of an axisymmetric reactor including experimental kinetic and calculated thermodynamic parameters for both reactions did indeed predict these outcomes. The extent to which the effects of thermodynamic coupling actually manifested, however, were also revealed to depend strongly on the relative rates of the two reactions, as well as the magnitude of the equilibrium constant governing the progress of the endothermic process. Given the evident complexity of attaining high effective % wt. H2 capacity with J, alternative systems were investigated for greater facility of extensive H2 release. Among those studied, 1,2,4,5-bis-BN-cyclohexane (H) demonstrated the most favorable properties, particularly with respect to thermal stability: rather than decompose, a neat sample instead sublimed when heated above 150 °C. Nevertheless, two commercially available catalytic systems were identified to effect release of two H2 equivalents from H. Release of further equivalents were apparently impeded by the formation of either polymeric material or one of two dimeric cage compounds depending on the catalyst used. Notably, a method to regenerate H from these product mixtures remains to be developed. Thus, while H may prove useful for certain long-term energy storage needs, it is currently less suited applications involving frequent fuel consumption. Similar difficulties were also encountered in attempts to realize the complete fuel cycle of 1,3-BN-cyclohexane (E) and B,N-substituted derivatives thereof. It had been initially proposed that E would provide for readily reversible BN dehydrogenation through a measure of frustrated Lewis pair-type character. Indeed, computations predicted this reaction would be essentially thermoneutral in solution. In the course of attempts to fully hydrogenate the spent fuel, however, dimeric species formed and proved resistant to further BN reduction. While a number of monomeric cyclic compounds were also successfully synthesized as formal boron-nitrogen frustrated Lewis pairs, none demonstrated any capacity to split H2 across the BN unit. The challenge of developing a practical amine-borane-based material for readily reversible hydrogen storage thus remains unresolved at this time. As such, it deserves consideration as a major objective of any future work.

08 HYDROGEN↗

Shaping the composition profiles in heteroepitaxial quantum dots: Interplay of thermodynamic and kinetic effects

Atomistic Monte Carlo simulations, coupling thermodynamic and kinetic effects, resolve a longstanding controversy regarding the origin of composition profiles in heteroepitaxial SiGe quantum dots. It is shown that profiles with cores rich in the unstrained (Si) component derive from near-equilibrium processes and intraisland diffusion. Profiles with cores rich in the strained (Ge) component are of nonequilibrium nature, i.e., they are strain driven but kinetically limited. They are shaped by the distribution of kinetic barriers of atomic diffusion in the islands. The diffusion pathways are clearly revealed for the first time. Geometrical kinetics play a minor role.

36 MATERIALS SCIENCE↗

On computing the evolution of temperature for materials under dynamic loading

Modeling and simulation of the dynamic response of materials is important to many applications including the development of armor systems, understanding the safety of explosives, and assessing the crashworthiness of vehicles. Within such applications it is often critical to accurately compute the evolution of the temperature because it is a state variable that affects the kinetics of competing active processes within the material (e.g., dislocation motion, phase transformation, decomposition). Depending on the selection of an independent state variable, e.g. temperature or entropy, the approach for computing temperature is well understood based on the thermodynamic framework attributed to Coleman and Noll. However, different computational codes used for modeling the dynamic response of materials adopt different independent state variables. In this work, two thermodynamically consistent strategies for computing the temperature of a coupled thermodynamic state are compared and implemented into two different Lagrangian computational codes. The equivalence of these two approaches is established through the numerical solutions of several test problems. Finally, the implication of various approximations made to each of these approaches within the literature are assessed in the context of uniaxial stress conditions (split Hopkinson pressure bar experiments) and uniaxial strain conditions (plate impact experiments). Here, it is shown that the temperature rate or energy partition approaches are equivalent when implemented in their complete forms, but that several common simplifying assumptions, that are warranted in the case of uniaxial stress, lead to significant errors in the resulting Hugoniot state for plate impact.

36 MATERIALS SCIENCE↗

Electron Bifurcation: Thermodynamics and Kinetics of Two-Electron Brokering in Biological Redox Chemistry

How can proteins drive two electrons from a redox active donor onto two acceptors at very different potentials and distances? And how can this transaction be conducted without dissipating very much energy or violating the laws of thermodynamics? Nature appears to have addressed these challenges by coupling thermodynamically uphill and downhill electron transfer reactions, using two-electron donor cofactors that have very different potentials for the removal of the first and second electron. Although electron bifurcation is carried out with near perfection from the standpoint of energy conservation and electron delivery yields, it is a biological energy transduction paradigm that has only come into focus recently. This Account provides an exegesis of the biophysical principles that underpin electron bifurcation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Differential Ligation Alters Electronic State and Coupling Signals of Iron-Sulfur Clusters in Flavin-Based Electron Bifurcation

Flavin-based electron bifurcation (FBEB) is employed by microorganisms for controlling pools of redox equivalents by reversibly splitting electron pairs into high- and low-energy levels from an initial midpoint potential. Our ability to harness this phenomenon is crucial for biocatalytic design which is limited by our understanding of energy coupling in the bifurcation system. In Pyrococcus furiosus, FBEB is carried out by the NADH-dependent ferredoxin:NADP+-oxidoreductase (NfnSL), coupling the uphill reduction of ferredoxin in NfnL to the downhill reduction of NAD+ in NfnS from oxidation of NADPH. Flanking the bifurcating flavin are two site-differentiated iron-sulfur clusters; the nearest is a glutamate-ligated [4Fe-4S] cluster in NfnL. Recent biochemical experiments substituting the native glutamate with cysteine led to loss of coupling between the uphill and downhill pathways, in contrast to the tight thermodynamic coupling in the native system. To understand how this decoupling is biochemically manifested by the cysteine-substituted [4Fe-4S] in NfnL, we employed electron paramagnetic resonance (EPR) spectroscopy to identify changes in electronic architecture and square wave voltammetry (SWV) to probe thermodynamic shifts produced by the substitution. We observed notable g-value shifts in the EPR for the cysteine-substituted iron-sulfur cluster in addition to significant downward shifts in the redox potential, as well as the disappearance of several low-field signals observed in the native NfnSL complex. These results suggest the site-differentiated glutamate residue facilitates higher spin states in the [4Fesingle bond4S] cluster to bridge energetic gaps in electron transfer to the bifurcating flavin in the native complex, preventing unwanted short-circuiting seen in the cysteine-substituted complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic implications of excited high-spin states, spin–spin coupling, and differential [2Fe–2S] + cluster temperature relaxations in the electron-bifurcating NfnABC from Thermococcus sibiricus

Electron bifurcation (EB) is a mechanism of biological energy transduction in which multiple oxidation–reduction (redox) reactions are thermodynamically coupled within a single enzyme, enabling the enzyme to harness the excess free energy from an exergonic process to drive an endergonic process. Because of this unprecedented chemistry, there is interest to translate EB principles to artificial and bioengineered systems, but a hurdle is that knowledge pertaining to the fundamental design principles of EB enzymes remains scarce. Here, we investigated the fundamental physical and electronic properties of electron transfer sites in a spectroscopically uncharacterized member of the BfuABC family of EB enzymes, the NADH-dependent reduced-ferredoxin:NADP + oxidoreductase from Thermococcus sibiricus (Tsi NfnABC). Cryo-EM structures of Tsi NfnABC previously demonstrated that it contains twelve redox cofactors: two flavins (one FAD and one FMN), eight [4Fe–4S] clusters, and two [2Fe–2S] clusters. The FMN, one [4Fe–4S] cluster, and one [2Fe–2S] cluster comprise the bifurcating active site termed the electron-bifurcating flavobicluster (BF-FBC), which is found in all BfuABC family members. By using electron paramagnetic resonance spectroscopy, we identified spectral signatures originating from interactions between the FMN radical and [4Fe–4S] + cluster in the BF-FBC and observed temperature dependent behavior of the BF-FBC's [2Fe–2S] + cluster indicative of moderately slow spin–lattice relaxation. Additionally, we uncovered numerous spectral features corresponding to half-integer, S > ½ spin states of [4Fe–4S] + clusters, including one attributable to the consequences of lysine-ligation of a [4Fe–4S] cluster unique to NfnABC. We contextualize these findings to electron transfer theory and NfnABC's structure. Our insights further the understanding of how enzymes are designed to exert control over electron transfer to conduct thermodynamically challenging reactions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Pi-CO₂ aqueous post-combustion CO₂ capture: Proof of concept through thermodynamic, hydrodynamic, and gas-lift pump modeling

Partnering in Innovation, Inc. (Pi-Innovation) introduces an aqueous post-combustion carbon dioxide (CO₂) capture system (Pi-CO₂) that offers high market value by directly addressing the primary constraints limiting beneficial re-use markets (lowering parasitic energy costs, reducing delivered cost of capture, eliminating the need for special solvents, etc.). A highly experienced team has completed initial design, modeling, manufacturing verification, and financial analysis for commercial market entry. Coupled thermodynamic and thermal-hydraulic mass transfer modeling results fully support proof of concept. Pi-CO₂ has the potential to lower total cost and risk to levels sufficient to stimulate global demand for CO₂ from local industrial sources.

54 ENVIRONMENTAL SCIENCES↗

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulation of Enhanced Geothermal Systems: A Benchmarking and Code Intercomparison Study

Numerical simulation codes have become critical tools for understanding complex geologic processes, as applied to technology assessment, system design, monitoring, and operational guidance. Recently the need for quantitatively evaluating coupled Thermodynamic, Hydrologic, geoMechanical, and geoChemical (THMC) processes has grown, driven by new applications such as geologic sequestration of greenhouse gases and development of unconventional energy sources. Here we focus on Enhanced Geothermal Systems (EGS), which are man-made geothermal reservoirs created where hot rock exists but there is insufficient natural permeability and/or pore fluids to allow efficient energy extraction. In an EGS, carefully controlled subsurface fluid injection is performed to enhance the permeability of pre-existing fractures, which facilitates fluid circulation and heat transport. EGS technologies are relatively new, and pose significant simulation challenges. To become a trusted analytical tool for EGS, numerical simulation codes must be tested to demonstrate that they adequately represent the coupled THMC processes of concern. This presentation describes the approach and status of a benchmarking and code intercomparison effort currently underway, supported by the U. S. Department of Energy’s Geothermal Technologies Program. This study is being closely coordinated with a parallel international effort sponsored by the International Partnership for Geothermal Technology (IPGT). We have defined an extensive suite of benchmark problems, test cases, and challenge problems, ranging in complexity and difficulty, and a number of modeling teams are applying various simulation tools to these problems. The descriptions of the problems and modeling results are being compiled using the Velo framework, a scientific workflow and data management environment accessible through a simple web-based interface.

benchmark↗

Layer thickness and period as design parameters to tailor pyroelectric properties in ferroelectric superlattices

We theoretically examine the pyroelectric properties of ferroelectric-paraelectric superlattices as a function of layer thickness and configuration using non-linear thermodynamics coupled with electrostatic and electromechanical interactions between layers. We specifically study PbZr{sub 0.3}Ti{sub 0.7}O{sub 3}/SrTiO{sub 3} superlattices. The pyroelectric properties of such constructs consisting of relatively thin repeating units are shown to exceed the pyroelectric response of monolithic PbZr{sub 0.3}Ti{sub 0.7}O{sub 3} films. This is related to periodic internal electric fields generated due to the polarization mismatch between layers that allows tailoring of the shift in the transition temperature. Our results indicate that higher and electric field sensitive pyroresponse can be achieved from layer-by-layer engineered ferroelectric heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A direct numerical simulation-based investigation and modeling of pressure Hessian effects on compressible velocity gradient dynamics

The pressure Hessian tensor plays a key role in shaping the behavior of the velocity gradient tensor, and in turn, that of many incumbent non-linear processes in a turbulent flow field. In compressible flows, the role of pressure Hessian is even more important because it represents the level of fluid-thermodynamic coupling existing in the flow field. In this work, we first perform a direct numerical simulation-based study to clearly identify, isolate, and understand various important inviscid mechanisms that govern the evolution of the pressure Hessian tensor in compressible turbulence. The ensuing understanding is then employed to introduce major improvements to the existing Lagrangian model of the pressure Hessian tensor (the enhanced Homogenized Euler equation or EHEE) in terms of (i) non-symmetric, non-isentropic effects and (ii) improved representation of the anisotropic portion of the pressure Hessian tensor. Finally, we evaluate the new model extensively by comparing the new model results against known turbulence behavior over a range of Reynolds and Mach numbers. Indeed, the new model shows much improved performance as compared to the EHEE model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Design for Highly Piezoelectric and Visible/Near-Infrared Photoresponsive Perovskite Oxides

Defect–engineered perovskite oxides that exhibit ferroelectric and photovoltaic properties are promising multifunctional materials. Though introducing gap states by transition metal doping on the perovskite B–site can obtain low bandgap (i.e., 1.1–3.8 eV), the electrically leaky perovskite oxides generally lose piezoelectricity mainly due to oxygen vacancies. Therefore, the development of highly piezoelectric ferroelectric semiconductor remains challenging. Here, inspired by point–defect–mediated large piezoelectricity in ferroelectrics especially at the morphotropic phase boundary (MPB) region, an efficient strategy is proposed by judiciously introducing the gap states at the MPB where defect–induced local polar heterogeneities are thermodynamically coupled with the host polarization to simultaneously achieve high piezoelectricity and low bandgap. A concrete example, Ni 2+ –mediated (1– x )Na 0.5 Bi 0.5 TiO 3– x Ba(Ti 0.5 Ni 0.5 )O 3–δ ( x = 0.02–0.08) composition is presented, which can show excellent piezoelectricity and unprecedented visible/near–infrared light absorption with a lowest ever bandgap ≈0.9 eV at room temperature. In particular, the MPB composition x = 0.05 shows the best ferroelectricity/piezoelectricity ( d 33 = 151 pC N –1 , P r = 31.2 μC cm –2 ) and a largely enhanced photocurrent density approximately two orders of magnitude higher compared with classic ferroelectric (Pb,La)(Zr,Ti)O 3 . Furthermore, this research provides a new paradigm for designing highly piezoelectric and visible/near–infrared photoresponsive perovskite oxides for solar energy conversion, near–infrared detection, and other multifunctional applications.

36 MATERIALS SCIENCE↗