Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cost function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High sensitivity pressure and temperature quantum sensing in pentacene-doped p-terphenyl single crystals

Quantum sensors’ responsiveness to their physical environment enables detection of variables such as temperature (T), pressure (P), and strain. We present a molecular platform for PT sensing using para-terphenyl crystals doped with pentacene (PDP), leveraging optically detected magnetic resonance (ODMR) of photoexcited triplet electron spins. We observe maximal frequency variations of df/dP=1.8 MHz/bar from 0-8 bar and df/dT=247 kHz/K from 79–330 K, over 1200 times and threefold greater, respectively, than those seen with nitrogen-vacancy centers in diamond and > 85-fold greater pressure sensitivity over the previous record. Density functional theory calculations indicate picometer-level PT-induced molecular orbital shifts are measurable via ODMR. PDP offers additional advantages including high sensor doping levels, narrow ODMR linewidths, high contrast, and low-cost single crystal growth. Overall, this work reports low-cost, optically-interrogated PT sensors and lays the foundation for increased versatility of quantum sensors through synthetic molecular design.

Magnetic properties and materials↗
DOE OSTI2020

The joint automated repository for various integrated simulations (JARVIS) for data-driven materials design

The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and development time of materials discovery, optimization, and deployment. The major features of JARVIS are: JARVIS-DFT, JARVIS-FF, JARVIS-ML, and JARVIS-tools. To date, JARVIS consists of ≈40,000 materials and ≈1 million calculated properties in JARVIS-DFT, ≈500 materials and ≈110 force-fields in JARVIS-FF, and ≈25 ML models for material-property predictions in JARVIS-ML, all of which are continuously expanding. JARVIS-tools provides scripts and workflows for running and analyzing various simulations. We compare our computational data to experiments or high-fidelity computational methods wherever applicable to evaluate error/uncertainty in predictions. In addition to the existing workflows, the infrastructure can support a wide variety of other technologically important applications as part of the data-driven materials design paradigm. The JARVIS datasets and tools are publicly available at the website: https://jarvis.nist.gov .

36 MATERIALS SCIENCE↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗
DOE OSTI2023

Importance of nonuniform Brillouin zone sampling for ab initio Bethe-Salpeter equation calculations of exciton binding energies in crystalline solids

Excitons are prevalent in semiconductors and insulators, and their binding energies are critical for optoelectronic applications. The state-of-the-art method for first-principles calculations of excitons in extended systems is the ab initio GW-Bethe-Salpeter equation (BSE) approach, which can require a fine sampling of reciprocal space to accurately resolve solid-state exciton properties. Here, in this study, we show, for a range of semiconductors and insulators, that the commonly employed approach of uniformly sampling the Brillouin zone can lead to underconverged exciton binding energies, as impractical grid sizes are required to achieve adequate convergence. We further show that nonuniform sampling of the Brillouin zone, focused on the region of reciprocal space where the exciton wave function resides, enables efficient rapid numerical convergence of exciton binding energies at a given level of theory. We propose a well-defined convergence procedure, which can be carried out at relatively low computational cost and which in some cases leads to a correction of previous best theoretical estimates by almost a factor of 2, qualitatively changing the predicted exciton physics. These results call for the adoption of nonuniform sampling methods for ab initio GW-BSE calculations and for revisiting previously computed values for exciton binding energies of many systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

ArborX

ArborX library tackles a problem of efficiently finding geometric objects that are close in space. Variations of this problem, such as finding the nearest neighbors of a point, or finding all objects within a certain distance, are inherent components of applications in many fields. The data may be large so that solving the problem efficiently may require significant computational resources, such as multiple processors or accelerators such as general purpose GPUs. ArborX' main advantage in its ability to solve large problems efficiently utilizing a combination of distributed and on-node parallelism. ArborX can be run efficiently on a wide variety of hardware, including GPUs from different vendors, which distinguishes it from other available libraries which typically choose only few of these. The other advantage is that it supports both types of user problems: spatial problems (useful for intersections and finding objects within certain distance), and nearest neighbor problems. ArborX also supports flexible interface in its interaction with a user. Particularly, it allows a user to call user's own function on a positive match, a functionality not rarely available in other libraries. ArborX implements construction and traversal algorithms using efficient tree structures, such as bounding volume hierarchy (BVH). At its core, it uses linear BVH for its low construction cost and sufficient quality. ArborX is written using C++, and is parallelized using the message passing interface (MPI) for the distributed communication, and the Kokkos library for on-node parallelism. This approach allows ArborX to be run on a wide variety of hardware, from common laptops and desktops to supercomputers while using the same codebase. ArborX also implements several advanced algorithms using geometric search, such as density-based clustering algorithm DBSCAN.

ECP↗
DOE OSTI2025

ArborX 2.0

ArborX library tackles a problem of efficiently finding geometric objects that are close in space. Variations of this problem, such as finding the nearest neighbors of a point, or finding all objects within a certain distance, are inherent components of applications in many fields. The data may be large so that solving the problem efficiently may require significant computational resources, such as multiple processors or accelerators such as general purpose GPUs. ArborX' main advantage in its ability to solve large problems efficiently utilizing a combination of distributed and on-node parallelism. ArborX can be run efficiently on a wide variety of hardware, including GPUs from different vendors, which distinguishes it from other available libraries which typically choose only few of these. The other advantage is that it supports both types of user problems: spatial problems (useful for intersections and finding objects within certain distance), and nearest neighbor problems. ArborX also supports flexible interface in its interaction with a user. Particularly, it allows a user to call user's own function on a positive match, a functionality not rarely available in other libraries. ArborX implements construction and traversal algorithms using efficient tree structures, such as bounding volume hierarchy (BVH). At its core, ArborX uses linear BVH for its low construction cost and sufficient quality. ArborX implements both spatial and nearest-neighbor traversal algorithms. ArborX also provides several clustering algorithms (minimum spanning tree, DBSCAN, HDBSCAN*), interpolation using minimum least squares and ray tracing. ArborX is written using C++, and is parallelized using the message passing interface (MPI) for the distributed communication, and the Kokkos library for on-node parallelism. This approach allows ArborX to be run on a wide variety of hardware, from common laptops and desktops to supercomputers while using the same codebase.

Prokopenko, Andrey [Oak Ridge National Laboratory ↗
DOE OSTI2024

Development of Human and Technology Integration Guidance for Work Optimization and Effective Use of Information

Existing nuclear power contributes to roughly 20% of the total electricity generation, and consistently provides the highest capacity factor of any other electricity generating resource in the United States. Despite these advantages, the existing nuclear power plant fleet in the United States has been enduring significant challenges in providing electricity in a cost-competitive manner, which has ultimately threatened the long-term economic viability of these plants. A major contributor to these increased operating and maintenance costs has been the continued reliance of a large workforce who perform their work under an operating model that has largely remained unchanged since the commissioning of these nuclear power plants. Unfortunately, while this operating model has provided safe and reliable electricity, other industries have already began transforming their workforce through the use of advanced digital technology and automation that has reduced their operating and maintenance cost significantly. In order for the United States nuclear industry to remain economically viable, a similar transformation must be considered in which digital technologies and automation capabilities are brought in to support key plant functions across all work functions across plant operation, maintenance, and support. To effectively integrate digital technology and automation in the existing nuclear power plant operating model, a multidisciplinary approach is needed that addresses technological and sociotechnical (i.e., human and technology integration) considerations. This report describes an extension to the human and technology integration methodology, herein referred to as Human Integration and Technology Task Force for Work Management Optimization (HITT), to support the safe, reliable, and efficient use of proposed innovations with the intended users in their intended environment to perform their intended tasks. By effectively incorporating human and technology integration into a plant modernization effort by using HITT, we believe that a utility can significantly improve work performance and overall workforce quality of life. For instance, if HITT is performed to optimize work management, scheduling performance and scope stability can be improved by up to 10%. HITT enables these benefits through developing a rich understanding of the work being performed, the utility’s vision, and the opportunities that provide greatest value to optimize performance through a 10-step process illustrated below as a quick guide. The quick guide contains links in the righthand side of each step that allows for navigation to the detailed section (Section 4) of HITT in this technical report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Non-reactive facet specific adsorption as a route to remediation of chlorinated organic contaminants

The present work quantifies metal-contaminant interactions between palladium substrates and three salient chlorinated organic contaminants, namely trichloroethylene 1,3,5-trichlorobenzene (TCB), and 3,3′,4,4′-tetrachlorobiphenyl (PCB77). Given that Pd is one of the conventional catalytically active materials known for contaminant removal, maximizing catalytic efficiency through optimal adsorption dynamics reduces the cost of remediation of contaminants that are persistent water pollutants chronically affecting public health. Adsorption efficiency analyses from all-atom molecular dynamics (MD) simulations advance the understanding of reaction mechanisms available from density functional theory (DFT) calculations to an extractable feature scale that can fit the parametric design of supported metal catalytic systems and feed into high throughput catalyst selection. Data on residence time, site-specific adsorption, binding energies, packing geometries, orientation profiles, and the effect of adsorbate size show the anomalous behaviour of organic contaminant adsorption on the undercoordinated {110} surface as compared to the {111} and {100} surfaces. The intermolecular interaction within contaminants from molecular dynamics simulation exhibits refreshing results than ordinary single molecule density functional theory calculation. Since complete adsorption and dechlorination is an essential step for chlorinated organic contaminant remediation pathways, the presented profiles provide essential information for designing efficient remediation systems through facet-controlled palladium nanoparticles.

Guo, Hao↗
DOE OSTI2024

Semi-Lagrangian nodal discontinuous Galerkin method for the BGK model

In this work, we propose a semi-Lagrangian (SL) nodal discontinuous Galerkin (DG) solver for the BGK equation. The BGK model was introduced by Bhatnagar, Gross, and Krook [1] as a relaxation model for the fundamental Boltzmann equation [5], which describes the kinetic dynamic of rarefied gases with a probability distribution function. The challenges of designing efficient numerical schemes for the Boltzmann equation mainly come from its high dimensionality and complicated nonlinear collision operator. The BGK model gains interests since it has much lower computational cost, due to the relatively simple structure of the relaxation operator in replacement of the collision operator, while simultaneously preserving several important physical properties, such as macroscopic quantities and dissipation of entropy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗
DOE OSTI2023

Hybrid renewable energy systems: the value of storage as a function of PV-wind variability

As shares of variable renewable energy (VRE) on the electric grid increase, sources of grid flexibility will become increasingly important for maintaining the reliability and affordability of electricity supply. Lithium-ion battery energy storage has been identified as an important and cost-effective source of flexibility, both by itself and when coupled with VRE technologies like solar photovoltaics (PV) and wind. In this study, we explored the current and future value of utility-scale hybrid energy systems comprising PV, wind, and lithium-ion battery technologies (PV-wind-battery systems). Using a price-taker model with simulated hourly energy and capacity prices, we simulated the revenue-maximizing dispatch of a range of PV-wind-battery configurations across Texas, from the present through 2050. Holding PV capacity and point-of-interconnection capacity constant, we modeled configurations with varying wind-to-PV capacity ratios and battery-to-PV capacity ratios. We found that coupling PV, wind, and battery technologies allows for more effective utilization of interconnection capacity by increasing capacity factors to 60%–80%+ and capacity credits to close to 100%, depending on battery capacity. We also compared the energy and capacity values of PV-wind and PV-wind-battery systems to the corresponding stability coefficient metric, which describes the location-and configuration-specific complementarity of PV and wind resources. Our results show that the stability coefficient effectively predicts the configuration-location combinations in which a smaller battery component can provide comparable economic performance in a PV-wind-battery system (compared to a PV-battery system). These PV-wind-battery hybrids can help integrate more VRE by providing smoother, more predictable generation and greater flexibility.

14 SOLAR ENERGY↗
DOE OSTI2023

Integrated Sustainable Product Design With Warranty and End-of-Use Considerations

The concept of integrated sustainable product design has recently emerged, aiming to incorporate downstream life cycle performance into the initial product design to enhance sustainability. Various sustainable product design tools based on life cycle assessment or quality function deployment have been established while the impact of reliability on circular practices has received limited attention. Recognizing the critical role of product reliability in post-design performance, this paper develops a product design optimization model that considers the warranty performance and the effect of end-of-use options. The model takes into account the effect of uncertain operating conditions on product reliability. Two optimization goals including the minimization of expected unit life cycle cost and environmental impact are achieved by the model. To demonstrate the benefits of the integrated approach, the model is applied to an electric motor design problem. The results highlight that integrating end-of-use options in the early design phase leads to adjustments in component selection and reliability design. Moreover, the circular utilization of used products enables cost savings throughout the product’s life cycle and contributes to environmental impact reduction. Finally, the study analyzes the effects of operating conditions, warranty policies, and take-back prices for used products on design decisions, providing valuable insights for product designers.

Integrated design↗
DOE OSTI2020

Advanced sorbents for modular oxygen production for REMS gasifiers

Under a DOE funded effort, Thermosolv LLC has been developing a sorbent-based oxygen production technology for gasification and oxy-combustion applications. The sorbents utilize oxygen-storage properties of certain perovskites to (1) selectively adsorb oxygen at moderate temperature from compressed air and (2) release the adsorbed oxygen into a vacuum or a sweep gas such as CO2 and/or steam. Through cyclic operations of multiple sorbent beds such a process can be made continuous. Pressure drop across the sorbent bed and other similar design considerations dictate that the sorbent be in the form of a high surface-area-to-volume pellet and yet still possess sufficient crush strength and integrity. Process efficiency is determined by the sorption and desorption kinetics and sorbent capacity. Typical process involves cycle times in the order of 100 seconds, a time too short to fully utilize the full volume of the sorbent and the sorption capacity of the sorbent pellets. As a part of this project, Thermosolv LLC undertook the development of oxygen sorbents as a high-surface area supported sorbent on a light-weight inert support to (1) reduce the cost, (2) increase the productivity, and (3) reduce the overall weight of the reactor. Using a previously developed perovskite (LSCF 1991) with its well-characterized performance, composite sorbent pellets each consisting of an inert core coated by a thin layer of the functional perovskite material were produced. Several low-density, inexpensive inert solids supports in the 1/8”-1/4” size were selected to keep the pressure drop across the sorbent bed in the acceptable range. A number of commercial technologies including spray coating, dip coating, incipient impregnation of porous substrate followed by thermal annealing at various temperatures and duration were employed and tested to produce robust composite supported sorbent pellets. A sense of optimum coating thickness was developed by testing sorbent pellets of various increasing diameter pellets. LSCF-1991extrudates ranging in size from 1/32” to 3/16” were tested in a TGA and in a fixed-bed pressure swing test set-up for sorption/desorption cycles of interest. The data show that for the operational conditions of interest the coating thickness for a supported sorbent pellet needs to be approximately 1/32” (about 0.8 mm). Subsequent work thereby concentrated on developing composite pellets of various substrates, shapes and sizes with coating layer of about this size range. Candidate support materials were chosen based on cost, inertness, mechanical strength, thermal expansion and chemical stability with sorbent material, and in a size range to give an acceptable pressure drop in fixed-bad reactor configurations. Coating application methods used for application of the sorbent onto the support included precipitation, spray coating and dip coating from sorbent slurries and sol gels. For all coated supports where we could successfully apply a uniform coating of desired thickness with an acceptable handling performance in terms of exfoliation, TGA-based cyclic sorption/desorption testing was performed to determine cycling capacity of oxygen. Among nearly twenty different support materials tested, best adhesion performance was obtained from stainless saddle supports. In the bench-scale fixed-bad tests stainless steel saddle support composite pellets with a coating thickness in the 0.6 mm or so range, the composite sorbent pellet performance approached up to 95% of that of the 2 mm parent material pellets. In a parallel approach to reducing the cost of the sorbent, alternate sorbent formulations replacing/reducing the amount of cobalt in the LSCF family were also investigated. Successful formulations that could match the performance of LSCF 1991 were identified based on TGA cyclic tests as LSCF 1919 and LSF 1910. In the range of operational envelop of cycle times and other relevant process operating conditions, the reduced Co formulations showed comparable performance in the bench-scale fixed-bad cyclic operations. As a part of this project, we also attempted substitution of La and Sr with Ca and Ba, and Co with Mn, Ni and Cu with little success, but the overall project goal of reducing sorbent cost in terms of raw material and manufacturing expenses was successful.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2020

Perspectives on three-dimensional printing of self-assembling materials and structures

Spiders, silks, and webs are marvels of nature that have inspired high-performance material design, resilient structural design, and innovative manufacturing in the materials science and engineering fields. In this article, we review the roles and mechanisms of self-assembly and spinning processes for making hierarchical structures of synthetic and natural fibers with multifunctionality. We discuss the ability of spider webs to maintain their functionality under construction while withstanding predator, human, and environmental threats, which comes from their self-assembled structures from nanoscales to macroscales. We conclude by exploring the future of novel mobile, autonomous, time- and cost-effective spider-inspired three-dimensional printing technology for making resilient biomaterials, structures, and synthetic silk materials.

42 ENGINEERING↗
DOE OSTI2020

Polycaprolactone: A Promising Addition to the Sequential Infiltration Synthesis Polymer Family Identified through In Situ Infrared Spectroscopy

Infiltration of inorganic oxides inside polymers using sequential infiltration synthesis (SIS) is an effective method for creating materials for a broad range of applications. The reactions between various polymer functional groups and organometallic/inorganic precursors are unique, which makes it essential to understand the specific interactions for a range of precursors and polymers to enable predictive process design and to extend the utility of SIS to applications. In this paper, in situ Fourier transform infrared spectroscopy (FTIR) measurements have been performed during Al 2 O 3 and TiO 2 SIS in three different homopolymers: poly(methyl methacrylate) (PMMA), poly(epsilon-caprolactone) (PCL), and poly(2-vinylpyridine) (P2VP). From the FTIR intensity changes after precursor exposure and during the subsequent purge times, it is demonstrated quantitatively that the interaction dynamics and the stability of the intermediate complexes of these polymers with the metal precursors are substantially different. A key finding from this comparative study is that PCL interacts far more strongly with metal precursors, even though its carbonyl (C=O) and ester (C-O-R) functional groups are similar to those of the more weakly interacting PMMA. This behavior suggests that additional factors beyond the identity of the functional groups dictate how polymers interact with metal compounds in SIS. PCL, which has not previously been reported in SIS processes, may be an attractive polymeric template for implementing SIS with improved uniformity and cost-effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Voltage Phosphate Cathodes for Rechargeable Ca-Ion Batteries

Calcium-ion batteries (CIBs) are under investigation as next-generation energy storage devices due to their theoretically high operating potentials and lower costs tied to the high natural abundance of calcium. However, the development of CIBs has been limited by the lack of available positive electrode materials. In this study, for the first time, we report two functional polyanionic phosphate materials as high-voltage cathodes for CIBs at room temperature. NaV 2 (PO 4 ) 3 electrodes were found to reversibly intercalate 0.6 mol of Ca 2+ (81 mA h g –1 ) near 3.2 V (vs Ca 2+ /Ca) with stable cycling performance at a current density of 3.5 mA g –1 . The olivine framework material FePO 4 reversibly intercalates 0.2 mol of Ca 2+ (72 mA h g –1 ) near 2.9 V (vs Ca 2+ /Ca) at a current density of 7.5 mA g –1 in the first cycle. Structural, electronic, and compositional changes are consistent with reversible Ca 2+ intercalation into these two materials.

25 ENERGY STORAGE↗