Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain energy density”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High Efficiency Multijunction Devices: Solar Cells, Thermophotovoltaics, LEDs

High-efficiency III-V multijunction solar cells require multiple alloys with optimal bandgaps and excellent material quality. However, many useful alloys are lattice-mismatched to available substrates, and require strain-engineering tricks to enable efficient use. To this extent, we have developed low-dislocation density metamorphic material and optically-thick stress-balanced superlattices, enabling access to a wide variety of III-V alloys. Combined with transparent tunnel junctions, these components allow optimizing multijunction devices for a variety of applications. We have designed multijunction cells for the G173G-terrestrial and AM0-space spectra with record 39.5% and 34.2% efficiency, respectively, and also thermophotovoltaic cells for blackbody spectra that reach over 40% TPV-efficiency. Finally, using these PV components, we demonstrate a multijunction LED with high quantum efficiency.

high-efficiency↗
DOE OSTI2020

Beyond Volume Variation: Anisotropic and Protrusive Lithiation in Bismuth Nanowire

Materials storing energy via an alloying reaction are promising anode candidates in rechargeable lithium-ion batteries (LIBs) due to their much higher energy density than the current graphite anode. Until now, the volumetric expansion of such electrode particles during lithiation has been considered as solely responsible for cycling-induced structural failure. Here, we report different structural failure mechanisms using single-crystalline bismuth nanowires as the alloying-based anode. The Li–Bi alloying process exhibits a two-step transition, that is, Bi–Li 1 Bi and Li 1 Bi–Li 3 Bi. Interestingly, the Bi–Li 1 Bi phase transition occurs not only in the bulk Bi nanowire but also on the particle surface showing its characteristic behavior. The bulk alloying kinetics favors a Bi-(012)-facilitated anisotropic lithiation, whose mechanism and energetics are further studied using the density functional theory calculations. More importantly, the protrusion of Li 1 Bi nanograins as a result of anisotropic Li–Bi alloying is found to dominate the surface morphology of Bi particles. The growth kinetics of Li 1 Bi protrusions is understood atomically with the identification of two different controlling mechanisms, that is, the dislocation-assisted strain relaxation at the Bi/Li 1 Bi interface and the short-range migration of Bi supporting the off-Bi growth of Li 1 Bi. As loosely rooted to the bulk substrate and easily peeled off and detached into the electrolyte, these nanoscale protrusions developed during battery cycling are believed to be an important factor responsible for the capacity decay of such alloying-based anodes at the electrode level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Epitaxial growth of oriented CoO films by radio-frequency sputtering deposition

Rock-salt CoO is a p-type semiconductor and its Neel temperature is close to room temperature. CoO-based compounds are known as promising systems for renewable energy harvest with high efficiency. CoO with catalytic and exchange bias properties can be widely used for industrial applications. In this paper, we report high-quality stable CoO(111) and (100) films epitaxially grown on c-cut (0001) and r-cut ($10\bar1$2) alpha-Al 2 O 3 substrates, respectively, using radio-frequency sputtering deposition. X-ray diffraction (XRD) measurements revealed that the CoO films had a rock-salt structure (Fm3m) with lattice constants of 4.2477 Å and 4.2617 Å for film grown on (0001) and ($10\bar1$2) α-Al 2 O 3 substrates, respectively. CoO films with the best crystal quality were grown at a substrate temperature of similar to ~700°C. XRD measurements of CoO(111) films indicated a lack of structural residual strain, whereas CoO(100) films had substantial amounts of structural strain. X-ray reflectivity (XRR) and transmission electron microscopy measurements showed neither oxygen vacancy nor defects in both CoO(111) and (100) films. Further, XRR revealed that the mean electron density of the CoO films was nearly identical to a pure CoO and that the films were considerably stable under the atmosphere.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Acoustic Radiation Pressure

The theoretical foundation of acoustic radiation pressure in plane wave beams is reexamined. It is shown from finite deformation theory and the Boltzmann-Ehrenfest Adiabatic Principle that the Brillouin stress tensor (BST) is the radiation stress in Lagrangian coordinates (not Eulerian coordinates) and that the terms in the BST are not the momentum flux density and mean excess Eulerian stress but are simply contributions to the variation in the wave oscillation period resulting from changes in path length and true wave velocity, respectively, from virtual variations in the strain. It is shown that the radiation stress in Eulerian coordinates is the mean Cauchy stress (not the momentum flux density, as commonly assumed) and that Langevin's second relation does not yield an assessment of the mean Eulerian pressure, since the enthalpy used in the traditional derivations is a function of the thermodynamic tensions - not the Eulerian pressure. It is shown that the transformation between Lagrangian and Eulerian quantities cannot be obtained from the commonly-used expansion of one of the quantities in terms of the particle displacement, since the expansion provides only the difference between the value of the quantity at two different points in Cartesian space separated by the displacement. The proper transformation is obtained only by employing the transformation coefficients of finite deformation theory, which are defined in terms of the displacement gradients. Finite deformation theory leads to the result that for laterally unconfined, plane waves the Lagrangian and Eulerian radiation pressures are equal with the value (1/4)(2K) along the direction of wave propagation, where (K) is the mean kinetic energy density, and zero in directions normal to the propagation direction. This is contrary to the Langevin result that the Lagrangian radiation pressure in the propagation direction is equal to (2K) and the BST result that the Eulerian radiation pressure in that direction is the momentum flux density.

Finite deformation theory↗
DOE OSTI2025

Alloying multiple halide perovskites on the same sublattice in search of stability and target band gaps

Single-component halide perovskites (HPs) rarely satisfy all the necessary criteria for optoelectronic applications, such as achieving an optimal band gap while maintaining high chemical and structural stability. Alloying halide perovskites has emerged as a promising strategy, not only to enhance stability but also to fine-tune their electronic and optical properties. In this work, we explore multiple degrees of freedom in alloy design, considering different substitution sublattices sites (A, B, or X in ABX3 perovskites), various chemical species (isovalent and hetero-valent elements), and multi-component compositions on a given sublattice. Using first-principles calculations based on density functional theory (DFT), we investigate how compositional variations influence the electronic (band gap) and structural properties (mixing enthalpy) of HP alloys. Our approach employs the polymorphous cell model, allowing full local relaxation which breaks local symmetry while preserving global cubic symmetry—an essential framework for accurately modeling HPs. Our results reveal that X-site mixing (halogen substitution) primarily affects the valence band maximum, allowing target band gap engineering. Additionally, variations in halogen radii introduce internal strain through octahedral distortions, influencing the mixing enthalpy. A-site substitution, while not directly contributing to the band edge states, modifies structural stability via volume effects, indirectly impacting the band gap. B-site alloying plays a dominant role in band gap modulation, leading to either positive or negative band gap bowing. Specifically, isovalent B-site mixing (Sn–Pb) induces strong positive bowing, where the alloy band gap is smaller than the average gap of parent compounds, whereas hetero-valent mixing (Cd–Pb) results in pronounced negative bowing. As an aside, we investigate the competition between the excess energy of disordered alloys vs. that of long-range ordered double perovskites of the same compositions, seeking examples of ordered phases emerging from disordered alloys. Furthermore, our findings provide fundamental insights into the electronic and structural behavior of HP alloys, offering valuable design principles for the development of stable and efficient materials for next-generation photovoltaic and optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2020

Ferroelectricity in Si-Doped Hafnia: Probing Challenges in Absence of Screening Charges

The ability to develop ferroelectric materials using binary oxides is critical to enable novel low-power, high-density non-volatile memory and fast switching logic. The discovery of ferroelectricity in hafnia-based thin films, has focused the hopes of the community on this class of materials to overcome the existing problems of perovskite-based integrated ferroelectrics. However, both the control of ferroelectricity in doped-HfO2 and the direct characterization at the nanoscale of ferroelectric phenomena, are increasingly difficult to achieve. The main limitations are imposed by the inherent intertwining of ferroelectric and dielectric properties, the role of strain, interfaces and electric field-mediated phase, and polarization changes. In this work, using Si-doped HfO2 as a material system, we performed a correlative study with four scanning probe techniques for the local sensing of intrinsic ferroelectricity on the oxide surface. Putting each technique in perspective, we demonstrated that different origins of spatially resolved contrast can be obtained, thus highlighting possible crosstalk not originated by a genuine ferroelectric response. By leveraging the strength of each method, we showed how intrinsic processes in ultrathin dielectrics, i.e., electronic leakage, existence and generation of energy states, charge trapping (de-trapping) phenomena, and electrochemical effects, can influence the sensed response. We then proceeded to initiate hysteresis loops by means of tip-induced spectroscopic cycling (i.e., “wake-up”), thus observing the onset of oxide degradation processes associated with this step. Finally, direct piezoelectric effects were studied using the high pressure resulting from the probe’s confinement, noticing the absence of a net time-invariant piezo-generated charge. Our results are critical in providing a general framework of interpretation for multiple nanoscale processes impacting ferroelectricity in doped-hafnia and strategies for sensing it.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemical Oxidation in Garnet-Type Solid Electrolyte by Formation of Point Defects

All-solid-state batteries hold greater promise for improving safety and energy density over conventional battery technology employing organic liquid electrolytes. One of the required features of a Li + conducting solid electrolyte is electrochemical stability, attained thermodynamically or kinetically, within the targeted operating voltage and temperature ranges. Therefore, understanding of the oxidative or reductive degradation mechanism is important to allow the design of stable solid electrolyte materials. This work contributes to building an understanding of the oxidative degradation mechanism in lithium solid electrolytes at cell operating conditions. Here, we have focused on resolving the oxidative decomposition mechanism of Al-doped lithium garnet Li 6.28 Al 0.24 La 3 Zr 2 O 12 (LLZO) as a state-of-the-art inorganic ceramic electrolyte. By combining experimental and computational analyses, we show that oxidation of LLZO occurs by simultaneous loss of oxygen and lithium from the structure, resulting in substoichiometric LLZO, at a moderate temperature (80 °C) and a high electrode potential (4.3 V vs Li/Li + ). Based on X-ray absorption and diffraction analyses, we find that the zirconium coordination shells in LLZO contract while the crystal structure experiences positive chemical strain upon electrochemical oxidation. The results from ex situ structural characterization of both the local structure and crystal symmetry are supported by a substoichiometric LLZO with lithium and oxygen vacancies, modeled by density functional theory (DFT) calculations. These chemical and structural changes in LLZO suppress effective lithium-ion conductivity by an order of magnitude. Formation of lithium and oxygen vacancies in LLZO upon electrochemical oxidation is different from prior thermodynamic predictions of phase decomposition of LLZO. The difference here is that the experiments were conducted at near-room temperature, which can hinder the kinetics of phase separation, and thus, the resultant LLZO solid electrolyte is still single-phase but substoichiometric in Li and O. In conclusion, these findings contribute an important degradation mechanism of the electrolyte, relevant for practical operational conditions of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Strain-driven surface reconstruction and cation segregation in layered Li(Ni 1– x – y Mn x Co y )O 2 (NMC) cathode materials

The composition, structure and phase transformations occurring on cathode surfaces greatly affect the performance of Li-ion batteries. Li-Ion diffusion and surface–electrolyte interaction are two major phenomena that impact the capacity and cell impedance. The effects of surface reconstruction (SR) of cathode materials on the performance of Li-ion batteries are of current interest. However, the origin and evolution of the SR are still not well understood. In this work, density functional theory (DFT) calculations are used to investigate the processes taking place during surface segregation and reconstruction. Facet dependent segregation was found in Li(Ni 1– x – y Mn x Co y )O 2 (NMC) cathodes. Specifically, Co tends to segregate to the (104) surface of the primary particles within the transition metal layer, while Ni ions tend to segregate to the (012) surface in the Li layer, forming a SR. Experimental evidence shows the SR to be epitaxial with the bulk of the as-synthesized material, and the new SR phase is pinned to the NMC unit cell leading to a strained SR. Here, we demonstrate that strain can stabilize a spinel structure of the SR layers. Understanding the effects of surface strain opens a new avenue for the design of cathode materials with enhanced surface properties.

25 ENERGY STORAGE↗
DOE OSTI2023

Density functional study of atomic arrangements in CrMnFeCoNi high-entropy alloy and their impact on vacancy formation energy and segregation

Here, using the density functional theory-coupled Monte Carlo approach, we explored the chemical short-range order (SRO) and element segregation in equimolar CrMnFeCoNi alloy. We found that state-of-the-art approximation of random element distribution is only applicable at > 1100 K close to the melting temperature, while the Cr-Cr repulsion driving the system stabilization and accompanying the formation of cubic Cr sublattice, and mild Ni-Ni attraction are the most prominent pair interactions at lower temperatures. Chemical potential and vacancy formation energy calculations indicate that Cr is most sensitive to the local chemical environment, making Cr atoms most stabilized when the preferred SRO is introduced. While the vacancy formation is predicted equally probable among five constituting elements in the random solid solution, Cr and Ni atoms show the lowest vacancy formation energies in the structure with SRO. Furthermore, distinct element segregation was predicted in the vicinity of planar defects, including symmetric tilt grain boundary and stacking fault, which we correlated to the site- and chemistry-dependent atomic volume and bond lengths. It suggests that the local mechanical strain and bond energy induce the SRO development and element segregation: Namely, the system takes advantage of segregation of Ni atoms having large atomic volume or Cr-Cr pairs having elongated bond lengths to fill in the excess volume at defects that relaxes the mechanical strain field and optimizes bond energy distribution. The correlation between the SRO and properties of CrMnFeCoNi alloy needs further investigations, which is expected to greatly help understand and control the properties of high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Catch bond kinetics are instrumental to cohesion of fire ant rafts under load

Dynamic networks composed of constituents that break and reform bonds reversibly are ubiquitous in nature owing to their modular architectures that enable functions like energy dissipation, self-healing, and even activity. While bond breaking depends only on the current configuration of attachment in these networks, reattachment depends also on the proximity of constituents. Therefore, dynamic networks composed of macroscale constituents (not benefited by the secondary interactions cohering analogous networks composed of molecular-scale constituents) must rely on primary bonds for cohesion and self-repair. Toward understanding how such macroscale networks might adaptively achieve this, we explore the uniaxial tensile response of 2D rafts composed of interlinked fire ants ( S. invicta ). Through experiments and discrete numerical modeling, we find that ant rafts adaptively stabilize their bonded ant-to-ant interactions in response to tensile strains, indicating catch bond dynamics. Consequently, low-strain rates that should theoretically induce creep mechanics of these rafts instead induce elastic-like response. Our results suggest that this force-stabilization delays dissolution of the rafts and improves toughness. Nevertheless, above 35 % strain low cohesion and stress localization cause nucleation and growth of voids whose coalescence patterns result from force-stabilization. These voids mitigate structural repair until initial raft densities are restored and ants can reconnect across defects. However mechanical recovery of ant rafts during cyclic loading suggests that—even upon reinstatement of initial densities—ants exhibit slower repair kinetics if they were recently loaded at faster strain rates. These results exemplify fire ants’ status as active agents capable of memory-driven, stimuli-response for potential inspiration of adaptive structural materials.

Science & Technology - Other Topics↗
DOE OSTI2025

Analysis of soot formation from aviation fuels in laminar counterflow flames

Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.

Counterflow↗
DOE OSTI2022

Single-molecule magnet Mn 12 on GaAs-supported graphene: Gate field effects from first principles

In this work, we study gate field effects on the Mn 12 O 12 (COOH) 16 (H 2 O) 4 | graphene | GaAs heterostructure via first-principles calculations. We find that under moderate doping levels electrons can be added to but not taken from the single-molecule magnet Mn 12 O 12 (COOH) 16 (H 2 O) 4 (Mn 12 ). The magnetic anisotropy energy (MAE) of Mn 12 decreases as the electron doping level increases, due to electron transfer from graphene to Mn 12 and change in the band alignment between Mn 12 and graphene. At an electron doping level of -5.00×10 13 cm -2 , the MAE decreases by about 18% compared with zero doping. The band alignment between graphene and GaAs is more sensitive to electron doping than to hole doping, since the valence band of GaAs is close to the Fermi level. The GaAs substrate induces a small band gap in the supported graphene under zero gate field and a nearly strain-free configuration. Finally, we propose a vertical tunnel junction for probing the gate dependence of MAE via electron transport measurements.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Accounting for Sustainability and Performance: Contrasting Effects of PET-Derived Modifiers on Asphalt Durability

The valorization of waste plastics in roadway infrastructure requires careful balance between sustainability and performance, as ineffective modifiers may incur both environmental and financial costs. This study shows that two modifiers derived from the same waste plastic source (PET)-DOTP (N1,N4-(dioctyl)terephthalamide) and HETP (N1,N4-Bis(2-hydroxyethyl)terephthalamide)-exhibit opposite effects on asphalt properties, despite having nearly the same carbon footprint. DOTP improves elastic recovery in multiple creep recovery tests (MCRT), while HETP increases strain accumulation and reduces elasticity. Density functional theory (DFT) shows that DOTP's molecular structure can disperse well in bitumen mixtures by binding to both polar and non-polar components. In contrast, HETP binds only to the binder's polar sites, lacking adequate dispersion. This increases the risk of aggregation of polar compounds in the bitumen mixture, resulting in a more deformation-prone binder. A techno-economic analysis highlights trade-offs between sustainability and performance in industrial-scale synthesis. This study advocates for applying first-principles thinking in material design, coupled with experimental validation, to engineer plastic waste-derived additives that achieve both performance and sustainability goals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nitrogen-Doped PtNi Catalysts on Polybenzimidazole-Functionalized Carbon Support for the Oxygen Reduction Reaction in Polymer Electrolyte Membrane Fuel Cells

PtM (M = 3d transition metals) alloys are known as the promising oxygen reduction reaction catalysts and have been considered as the replacement of pure Pt catalysts for the commercialization of proton exchange membrane fuel cells. Although great progress has been made in the past three decades, the performance and durability of PtM catalysts still face stringent challenges from practical applications. Functionalization of a catalyst carbon support with nitrogen-contained groups can add charges onto its surface, which can be utilized to build a more complete ionomer/catalyst interface, to reduce the catalyst particle size, and to improve particle size distribution. Nitriding of PtNi catalysts can effectively improve the catalyst activity and stability by the modification of lattice strain. Hereby, we propose a synergistic approach of combining polybenzimidazole-grafted Vulcan XC72 carbon as the catalyst carbon support and the nitriding of PtNi to develop PtNiN/XC72-polybenzimidazole catalysts. Such PtNiN/XC72-PBI catalysts exhibit the excellent performance of fuel cell membrane electrode assembly (i.e., mass activity, 440 mA mgPt –1 ; electrochemical surface area, 51 m 2 gPt –1 ; and rated power density, 836 mW cm –2 ) as well as promising catalyst stability. The developed PtNiN/XC72-PBI meets the US DOE 2020 targets of mass activity for the fuel cell catalysts. This work provides a novel approach and a promising pathway on the development of the catalyst using such a synergistic approach-modification of the catalyst structure by nitrogen doping and functionalization of carbon support by polybenzimidazole for both high performance and high durability.

25 ENERGY STORAGE↗
DOE OSTI2021

Nature of Oxygen Adsorption on Defective Carbonaceous Materials

Here, plane-wave density functional theory has been used to study oxygen adsorption on graphene, graphite, and (12,0) zigzag single-walled carbon nanotubes with and without Stone–Wales (SW) and single-vacancy (SV) defects to understand the role of defects on carbonaceous material reactivity. Atomic oxygen adsorption leads to the formation of an epoxide on defect-free graphene and graphite and an ether on the exterior wall of carbon nanotubes and SW-defected materials. O 2 chemisorption is endothermic on defect-free graphene and graphite and slightly exothermic on defect-free nanotubes. O 2 chemisorption energies are predicted to be –1.1 to –1.4 eV on an SW defect and –6.0 to –8.0 eV on an SV defect. An SW defect lowers the energy barriers by 0.90 and 0.50 eV for O 2 chemisorption on graphene and nanotubes, respectively. The formation of a C–O–O–C group is important for O 2 dissociation on defect-free and SW-defected materials. The energy barrier is less than 0.30 eV on an SV defect. The more reactive SW defect toward O adsorption on graphene is mostly due to the strained defective carbon atoms being able to donate more electrons to an O to form an ether. The larger 2s character in the hybrid orbitals in an ether than in an epoxide makes the ether C–O bond stronger. Stronger C–O binding on an SW-defective carbon nanotube than on a defect-free nanotube is in part due to more flexibility of the defect to release the epoxide ring strain to form an ether.

25 ENERGY STORAGE↗
DOE OSTI2023

A VOI Web Application for Distinct Geothermal Domains: Statistical Evaluation of Different Data Types within the Great Basin

The Great Basin region contains different domains that have different structural and hydrothermal flow patterns. Depending on the characteristics of these patterns, certain data types may be more successful at detecting hidden geothermal resources. In this paper, we quantitatively evaluate if certain data types are more successful in certain domains. Given different aquifer, strain and structural conditions, we explore which data types statistically reveal positively labeled geothermal sites. We utilize value of information (VOI) metrics to help quantify the reliability of data types to discriminate against "positive" and "negative" labeled geothermal sites. We also evaluate how kernel density estimation can help generalize the statistics that inform VOI, which is necessary given the limited data in geothermal exploration. Except for the Carbonate Aquifer, the highest ranking of the Vimperfect is the Local Structural Setting. Next, the slip and dilation tendency is first for Carbonate Aquifer and second for Central Nevada Seismic Belt and Western Great Basin. For the Carbonate Aquifer, heat flow is has the lowest Vimperfect value compared to the other three domains, which is consistent with the understanding of how heat flow measurements are masked by regional groundwater flow.

Bayesian analysis↗
DOE OSTI2020

First-Principles Modeling of Sodium Ion and Water Intercalation into Titanium Disulfide Interlayers for Water Desalination

Recent experiments revealed the possibility of using titanium disulfide (TiS 2 ) as the cathode material in capacitive deionization (CDI) devices for water desalination. Although it performed stably up to 70 cycles with a salt removal capacity of 14 mg/g (corresponding to a sodium ion removal capacity of 35.8 mg/g) at high molar concentration (600 mM NaCl), the maximum capacity of TiS 2 as a CDI electrode was much lower (~70 mAh/g) than as a supercapacitor (239 mAh/g). Understanding why the ion capacities of these two configurations of the same material differ will entail elucidating detailed charge/discharge mechanisms at the atomic scale. Here, we present density functional theory simulations of sodium intercalation into TiS 2 interlayers in order to gain such understanding. We systematically investigated TiS 2 stacking patterns and ion intercalation sites, energetics of the intercalated compounds, and phase transformation during sodium intercalation. The calculated structural evolution and capacitance–voltage curve agree quite well with previous measurements. We conclude that the different maximum capacities of TiS 2 measured in aqueous and dry environments originate from weaker interlayer interactions with respect to shear strain after ~33% intercalation of Na + –H 2 O pairs, which is detrimental to the mechanical stability of TiS 2 . In conclusion, this study sheds light on the underlying mechanisms of ion intercalation into layered materials and contributes to understanding requirements for future design and optimization of CDI electrode materials for water desalination.

36 MATERIALS SCIENCE↗