Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Deterrent Concentration Measurement with FTIR and Subsequent Ballistic Performance in Medium Caliber Propellant

To meet ballistic requirements, medium and small caliber propellants use deterrent coatings to obtain burn rate progressivity. The required amount and distribution of deterrent varies between gun systems, propellant types, and often between lots. Micro Fourier Transform Infrared (FTIR) spectroscopy was used to measure deterrent gradients in RP36 propellants coated with methyl centralite (MC) at different deterrent levels and different processing conditions. The aromatic C-C bonds at 1496 cm(exp -1) wavenumber were used to monitor the deterrent profiles through the grain. Deterrent gradients measured with FTIR spectroscopy were then used to estimate burn rate gradients in the deterred grains. Burn rates were calculated from literature models and from closed bomb data of RP36 containing uniform deterrent concentration. Finally, the burn rate gradients were input into an IBHFG2 model of a 200 cc-closed bomb. The early flame spreading portion of the closed bomb ballistic cycle (0 to 0.2 P/Pmax) was roughly modeled by dividing the charge up into five propellant decks and igniting them at different times in the ballistic cycle. Pressure traces and vivacity curves from closed bomb shots were compared to predictions. In addition to the burn rate gradient, the closed bomb pressure trace was heavily dependent on ignition and flame spread. These two phenomena were not readily distinguishable from one another in deterred grains. The same RP-36 propellant was shot in a 25 mm M793TP round which was again modeled with IBHVG2. Peak pressure and muzzle velocity were accurately modeled when erosive burning effects were empirically factored into the model.

Furrow, Keith W.↗
DOE OSTI2014

Plasma response measurements of non-axisymmetric magnetic perturbations on DIII-D via soft x-ray imaging

Recent observations on DIII-D have advanced the understanding of plasma response to applied resonant magnetic perturbations (RMPs) in both H-mode and L-mode plasmas. Three distinct 3D features localized in minor radius are imaged via filtered soft x-ray emission: (i) the formation of lobes extending from the unperturbed separatrix in the X-point region at the plasma boundary, (ii) helical kink-like perturbations in the steep-gradient region inside the separatrix, and (iii) amplified islands in the core of a low-rotation L-mode plasma. In this study, these measurements are used to test and to validate plasma response models, which are crucial for providing predictive capability of edge-localized mode control. In particular, vacuum and two-fluid resistive magnetohydrodynamic(MHD) responses are tested in the regions of these measurements. At the plasma boundary in H-mode discharges with n = 3 RMPs applied, measurements compare well to vacuum-field calculations that predict lobe structures. Yet in the steep-gradient region, measurements agree better with calculations from the linear resistive two-fluid MHD code, M3D-C1. Relative to the vacuum fields, the resistive two-fluid MHD calculations show a reduction in the pitch-resonant components of the normal magnetic field (screening), and amplification of non-resonant components associated with ideal kink modes. However, the calculations still over-predict the amplitude of the measuredperturbation by a factor of 4. In a slowly rotating L-mode plasma with n = 1 RMPs, core islands are observed amplified from vacuum predictions. Finally, these results indicate that while the vacuum approach describes measurements in the edge region well, it is important to include effects of extended MHD in the pedestal and deeper in the plasma core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Minimizing inner product data dependencies in conjugate gradient iteration

The amount of concurrency available in conjugate gradient iteration is limited by the summations required in the inner product computations. The inner product of two vectors of length N requires time c log(N), if N or more processors are available. This paper describes an algebraic restructuring of the conjugate gradient algorithm which minimizes data dependencies due to inner product calculations. After an initial start up, the new algorithm can perform a conjugate gradient iteration in time c*log(log(N)).

Vanrosendale, J.↗
NASA NTRS1985

Far-infrared line intensities of H2O and CO from warm molecular clouds

Takahashi et al. (1983) have studied the role of H2O in the overall energetics of warm molecular clouds with large optical depths. It was found that in molecular clouds associated with nearby or embedded luminosity sources, the dust is likely to be warmer than the molecular gas. Thus, the warm dust represents the primary source of heat for the gas. In some cases, H2O was found to cool the gas. Takahashi et al. (1983) have presented the formalism for treating H2O in an optically thick molecular cloud with embedded sources of luminosity. The total H2O heating and cooling rates were also calculated. In the present investigation, the emergent intensities of far-infrared lines of H2O and excited, rotational CO transitions from relatively warm molecular clouds based on the large-velocity gradient model of the previous study are calculated. It is found that the dominant submillimeter CO transitions are related to emission, while the H2O transitions can be either absorption or emission lines.

Takahashi, T.↗
DOE OSTI2022

Techno-Economic Performance of Eavor Loop 2.0

This project evaluated techno-economic performance for a sample Eavor-Loop 2.0 design for electricity production and direct-use heating. The Eavor-Loop 2.0 design investigated is a 7.5-km deep closed-loop geothermal system consisting of 12 laterals for a total of more than 90 km of downhole well and lateral length. Both a high geothermal gradient scenario of 60 degrees C/km and a low geothermal gradient scenario of 30 degrees C/km were considered. With pure water injected at 60 degrees C and 80 kg/s, reservoir simulations with the Slender-Body Theory simulator indicate average production temperatures over a 30-year lifetime of ~125 degrees C and ~210 degrees C for the low and high geothermal gradient scenario, respectively. These correspond to heat production of ~22 M Wth and ~51 M Wth, respectively. Using IPSEpro simulations, we find average power production of ~2.2 M We and ~8.6 M We, respectively, for a subcritical organic Rankine cycle power plant with air-cooled condensers. Cost estimates indicate the overall capital and levelized costs are dominated by the lateral drilling cost. Obtaining a levelized cost of electricity below $70/M Wh requires a geothermal gradient of 60 degrees C/km, a discount rate below 9%, and lateral drilling cost below $400/m. A well cost model indicates that ~$400/m for the Eavor-Loop 2.0 design investigated can be obtained for a drilling rate of penetration about 40 ft/hr (with bit life of 50 hours), and omitting casing and cement. Traditional (geothermal) well drilling has achieved these drilling rate conditions, including the Utah FORGE project where the rate of penetration has exceeded 50 ft/hr in granite. However, it is unclear if these conditions are still valid for drilling the Eavor-Loop 2.0 laterals (i.e., ~82 km of laterals at 4 to 7.5-km vertical depth with rock temperatures up to 460 degrees C), as such downhole completion has never been developed before. Competitive levelized cost of heat values ($1.2-$8.2/GJ) are calculated, even for the low geothermal gradient scenario (30 degrees C/km) and lateral drilling cost of $600/m.

advanced geothermal system↗
DOE OSTI2021

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4a Raw Data

Calculated formation energy of (0001 ̅) c-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (0001 ̅) c-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗
DOE OSTI2021

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4b Raw Data

Calculated formation energy of (1010 ̅) m-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (1010 ̅) m-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
DOE OSTI2025

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗
DOE OSTI2016

Spin and charge thermopower effects in the ferromagnetic graphene junction

Using wave function matching approach and employing the Landauer-Buttiker formula, a ferromagnetic graphene junction with temperature gradient across the system is studied. We calculate the thermally induced charge and spin current as well as the thermoelectric voltage (Seebeck effect) in the linear and nonlinear regimes. Our calculation revealed that due to the electron-hole symmetry, the charge Seebeck coefficient is, for an undoped magnetic graphene, an odd function of chemical potential while the spin Seebeck coefficient is an even function regardless of the temperature gradient and junction length. We have also found with an accurate tuning external parameter, namely, the exchange filed and gate voltage, the temperature gradient across the junction drives a pure spin current without accompanying the charge current. Another important characteristic of thermoelectric transport, thermally induced current in the nonlinear regime, is examined. It would be our main finding that with increasing thermal gradient applied to the junction the spin and charge thermovoltages decrease and even become zero for non zero temperature bias.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

TU-AB-BRC-06: Dose Calculation in Curved Space

Purpose: Microbeam Radiation Therapy is a preclinical method in radiation oncology that modulates radiation fields on a micrometre scale. Dose calculation is challenging due to arising dose gradients and therapeutically important dose ranges. Monte Carlo (MC) simulations, often used as gold standard, are computationally expensive and hence too slow for the optimisation of treatment parameters in future clinical applications. On the other hand, conventional kernel based dose calculation leads to inaccurate results close to material interfaces. The purpose of this work is to overcome these inaccuracies while keeping computation times low. Methods: A point kernel superposition algorithm is modified to account for tissue inhomogeneities. Instead of conventional ray tracing approaches, methods from differential geometry are applied and the space around the primary photon interaction is locally warped. The performance of this approach is compared to MC simulations and a simple convolution algorithm (CA) for two different phantoms and photon spectra. Results: While peak doses of all dose calculation methods agreed within less than 4% deviations, the proposed approach surpassed a simple convolution algorithm in accuracy by a factor of up to 3 in the scatter dose. In a treatment geometry similar to possible future clinical situations differences between Monte Carlo and the differential geometry algorithm were less than 3%. At the same time the calculation time did not exceed 15 minutes. Conclusion: With the developed method it was possible to improve the dose calculation based on the CA method with respect to accuracy especially at sharp tissue boundaries. While the calculation is more extensive than for the CA method and depends on field size, the typical calculation time for a 20×20 mm{sup 2} field on a 3.4 GHz and 8 GByte RAM processor remained below 15 minutes. Parallelisation and optimisation of the algorithm could lead to further significant calculation time reductions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Pressure induced magnetic phase transition in RhFe{sub 3}N and IrFe{sub 3}N: An ab-initio study

The structural, electronic, elastic and magnetic properties of RhFe{sub 3}N and IrFe{sub 3}N are investigated using ab-initio calculations based on density functional theory as implemented in VASP code within the gradient generalized approximation. The non-spin polarized and spin polarized calculations are performed for these nitrides at normal and high pressures. It is found that these ternary nitrides are stable in ferromagnetic state at normal pressure. The lattice constant and bulk modulus values are calculated. The electronic structure reveals that these nitrides are metallic at normal pressure. The calculated elastic constants indicate that they are mechanically stable at ambient pressure. Ferromagnetic to nonmagnetic phase transition is observed in RhFe{sub 3}N and IrFe{sub 3}N at high pressure. Ferromagnetism is quenched in these nitrides at high pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

The running coupling of 8 flavors and 3 colors

We compute the renormalized running coupling of SU(3) gauge theory coupled to N f = 8 flavors of massless fundamental Dirac fermions. The recently proposed finite volume gradient flow scheme is used. The calculations are performed at several lattice spacings allowing for a controlled continuum extrapolation. The results for the discrete β-function show that it is monotonic without any sign of a fixed point in the range of couplings we cover. Finally, as a cross check the continuum results are compared with the well-known perturbative continuum β-function for small values of the renormalized coupling and perfect agreement is found.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

First-Principle Study on the Origin of Ferromagnetism in Non-magnetic Halide Perovskite CsCdM{sub 3} (M = Cl and Br) Doped with B, C, and N

Based on the density functional theory (DFT), we elucidate the origin of magnetization in non-magnetic cubic halide perovskite CsCdM{sub 3} (M = Cl and Br) induced by substitution of non-magnetic 2p-impurities (B, C, and N) substitution. The calculations were done within the generalized gradient approximation (GGA) and GGA plus-modified Becke and Johnson (mBJ-GGA) as the exchange correlation. The results reveal that half-metallic ferromagnetism can be obtained for C- and N-dopings with the integer magnetic moment of 3.00 and 2.00 µ{sub B} per cell. However, B substitution does not induce magnetism in both CsCdCl{sub 3} and CsCdBr{sub 3} systems. The localized magnetic moments substantially come from impurity atoms. The origin of ferromagnetism can be attributed to the p - p hybridization between 2p-impurities and its surrounding halide atoms. Our results show that doped perovskites could provide a new promising class of materials for future spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

A hybridized discontinuous Galerkin method on mapped deforming domains

Here in this paper we present a hybridized discontinuous Galerkin (HDG) discretization for unsteady simulations of convection-dominated flows on mapped deforming domains. Mesh deformation is achieved through an arbitrary Lagrangian-Eulerian transformation with an analytical mapping. We present details of this transformation applied to the HDG system of equations, with focus on the auxiliary gradient equation, viscous stabilization, and output calculation. We discuss conditions under which optimal unsteady output convergence rates can be attained, and we show that both HDG and discontinuous Galerkin (DG) achieve these rates in advection-dominated flows. Results for scalar advection-diffusion and the Euler equations verify the implementation of the mesh motion mapping for both discretizations, show that HDG and DG yield similar results on a given mesh, and demonstrate differences in output convergence rates depending on the choice of HDG viscous stabilization. We note that such similar results bode well for HDG, which has fewer globally-coupled degrees of freedom compared to DG. A simulation of the unsteady compressible Navier-Stokes equations demonstrates again very similar results for HDG and DG and illustrates a pitfall of using steady-state adapted meshes for accurate unsteady simulations.

42 ENGINEERING↗
DOE OSTI2018

Analysis of the structural, electronic, elastic and thermodynamic properties of CuAl2X4 (X = O, S) spinel structure

Highlights: • Fundamental properties of CuAl{sub 2}X{sub 4} (X = O, S) aluminate spinels have been investigated. • The pressure dependence of the elastic constants C{sub ij} can be successfully fit by a straight line. • CuAl{sub 2}X{sub 4} (X = O, S) compounds are mechanically stable at ambient conditions. • Temperature and pressure dependences of some macroscopic parameters are achieved. - Abstract: The structural, electronic, elastic and thermodynamic properties of oxo- and thio-spinels have been predicted based on the ultrasoft pseudo-potential scheme as implemented in the CASTEP code. The exchange-correlation potential was treated within the generalized gradient approximation (GGA-PBE). The calculated structural parameters such as lattice constants and internal parameters are in good agreement with the available experimental data. The ground state electronic band structure revealed that CuAl{sub 2}X{sub 4} (XO, S) spinels are direct band gap semiconductors with a transition along Γ-Γ symmetry points. The calculated values of elastic constants satisfy the Born criteria at ambient conditions and confirm that the spinel structures of CuAl{sub 2}O{sub 4} and CuAl{sub 2}S{sub 4} are mechanically stable. The thiospinel structure may have a phase transition under pressure, due to the negative value of C{sub 44}. The influence of temperature and pressure on macroscopic characteristics of CuAl{sub 2}O{sub 4} and CuAl{sub 2}S{sub 4} compounds was estimated using the quasi-harmonic Debye model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Thermodynamic Stability of Low- and High-Index Spinel LiMn 2 O 4 Surface Terminations

Density functional theory calculations are performed within the generalized gradient approximation (GGA+U) to determine stable terminations of both low- and high-index spinel LiMn2O4 (LMO) surfaces. A grand canonical thermodynamic approach is employed, permitting a direct comparison of offstoichiometric surfaces with previously reported stoichiometric surface terminations at various environmental conditions. Within this formalism, we have identified trends in the structure of the low-index surfaces as a function of the Li and O chemical potentials. The results suggest that, under a range of chemical potentials for which bulk LMO is stable, Li/O and Li-rich (111) surface terminations are favored, neither of which adopts an inverse spinel structure in the subsurface region. This thermodynamic analysis is extended to identify stable structures for certain high-index surfaces, including (311), (331), (511), and (531), which constitute simple models for steps or defects that may be present on real LMO particles. The low- and high-index results are combined to determine the relative stability of each surface facet under a range of environmental conditions. The relative surface energies are further employed to predict LMO particle shapes through a Wulff construction approach, which suggests that LMO particles will adopt either an octahedron or a truncated octahedron shape at conditions in which LMO is thermodynamically stable. These results are in agreement with the experimental observations of LMO particle shapes.

DFT↗