Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Configuration interaction”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bottlebrush Block Copolymer Thin Films

The self-assembly of block copolymers (BCPs) in thin films is governed by interfacial interactions (enthalpy) between the blocks and interfaces (at both the substrate and the surface). Advantage can also be taken of entropy to control the polymer orientation. Here, we synthesized a series of bottlebrush block copolymers (BBCPs), with deuterated polystyrene (DPS) as the core block and poly­(solketal acrylate) (PSA) as the corona block, where the backbone length (N BB) and grafting densities (GDs) were varied. The hydrophobic PSA block was converted to a hydrophilic poly­(glyceryl acrylate) (PGA) block by solid-state hydrolysis, bringing the BBCP from a disordered state into an ordered state with a lamellar microdomain morphology. The orientation of the morphology as a function of distance from the interfaces was systematically investigated by atomic force microscopy (AFM), interference microscopy, grazing-incidence small-angle X-ray scattering (GISAXS), and grazing-incidence small angle neutron scattering (GISANS). For N BB = 1 (diblock BCP), a mixed lamellar orientation was found. For N BB = 5 (star-like BBCP), a vertical orientation of the lamellar microdomains was found dominantly at the interfaces that propagated into the bulk of the film due to the entropic penalty associated with constraints on the junction points. Increasing the N BB to 50 (rod-like BBCP) and 100 (worm-like BBCP) yielded a random orientation, as the side chains became overcrowded when aligned in the same direction. By reducing the grafting density (GD) of a N BB = 50 BBCP from 100% to 75% and 50% by the incorporation of spacers into the backbone, the steric hindrance of side chains is reduced markedly, facilitating a vertical orientation of the lamellar microdomains. However, the vertical orientation becomes more random within the bulk of the thin film for different BBCP architectures. By modifying the substrate with a layer of PS, the preferential interactions between the core block of the BBCP and the interfaces induce a reconfiguration of the BBCP, drawing the core block out to contact substrate directly driven by enthalpy. Under these conditions, N BB = 5 showed a predominantly parallel orientation of the lamellar microdomains with the corona block minimizing contact with the substrate. Higher N BB values of 50 and 100 showed a more random orientation of the lamellar microdomains at different depths on the PS-modified substrate. Overall, thin-film morphology was regulated by tuning the configuration of the BBCP, varying the entropic contribution to the microdomain orientation.

Chen, Zhan↗
DOE OSTI2023

Structure determination of ζ-N 2 from single-crystal X-ray diffraction and theoretical suggestion for the formation of amorphous nitrogen

The allotropy of solid molecular nitrogen is the consequence of a complex interplay between fundamental intermolecular as well as intramolecular interactions. Understanding the underlying physical mechanisms hinges on knowledge of the crystal structures of these molecular phases. That is especially true for ζ-N 2 , key to shed light on nitrogen’s polymerization. Here, we perform single-crystal X-ray diffraction on laser-heated N 2 samples at 54, 63, 70 and 86 GPa and solve and refine the hitherto unknown structure of ζ-N 2 . In its monoclinic unit cell (space group C 2/ c ), 16 N 2 molecules are arranged in a configuration similar to that of ε-N 2 . The structure model provides an explanation for the previously identified Raman and infrared lattice and vibrational modes of ζ-N 2 . Density functional theory calculations give an insight into the gradual delocalization of electronic density from intramolecular bonds to intermolecular space and suggest a possible pathway towards nitrogen’s polymerization.

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

Reference lattice, sound, stiffness, and magnetic transitions of Ising monolayers

A reference lattice, away from which elastic distortions induced by the spin texturing of two-dimensional (2D) magnets take hold, is motivated from a picture of pairwise Biot-Savart interactions among identical solenoids that either elongate or compress a (“zero-current”) spring lattice. Applied to a paradigmatic CrSiTe 3 monolayer, the reference is given by the average between the atomic positions of ferromagnetic (FM) and Néel antiferromagnetic (AFM) lattices; such an atomic disposition permits understanding structural distortions and elastic energies due to magnetism readily. Furthermore, the anisotropic speed of sound in the magnetic ground state explains an observed anisotropy of vibrational frequencies on similar magnets. Elastic stiffness constants are reported, too. Magnetic energies in four Ising structural configurations were calculated and the strain needed for those 2D magnets to undergo an AFM to FM quantum phase transition was determined as well.

density functional theory↗
DOE OSTI2023

Determining the Electronic Structure and Thermoelectric Properties of MoS 2 /MoSe 2 Type‐I Heterojunction by DFT and the Landauer Approach

Abstract The electronic structure and thermoelectric properties of MoX 2 (X = S, Se) Van der Waals heterojunctions are reported, with the intention of motivating the design of electronic devices using such materials. Calculations indicate the proposed heterojunctions are thermodynamically stable and present a band gap reduction from 1.8 eV to 0.8 eV. The latter effect is highly related to interactions between metallic d‐character orbitals and chalcogen p‐character orbitals. The theoretical approach allows to predict a transition from semiconducting to semi‐metallic behavior. The band alignment indicates a type‐I heterojunction and band offsets of 0.2 eV. Transport properties show clear n‐type nature and a high Seebeck coefficient at 300 K, along with conductivity values (σ/τ) in the order of 10 20 . Lastly, using the Landauer approach and ballistic transport, the proposed heterojunctions can be modeled as a channel material for a typical one‐gate transistor configuration predicting subthreshold values of ≈60 mV dec −1 and field–effect mobilities of ≈160 cm −2 V −1 s −1 .

López‐Galán, Oscar A.↗
DOE OSTI2024

The use of idealised experiments in testing a new convective parametrization: Performance of CoMorph-A

CoMorph is a new mass-flux convection parametrization under development at the Met Office designed for use within the Unified Model and its successor model, LFRic. Use of a three-dimensional idealised model enables controlled tests of the performance of the scheme across different regimes. This includes the interaction between the physical parametrizations and the resolved dynamics, allowing study of the emergent organisation of convection on the resolved scale. A selection of well-known cases is revisited here, with the purpose of documenting the extent to which CoMorph captures a range of important, but challenging, behaviour such as the diurnal cycle and sensitivity to tropospheric moisture. Simulations using CoMorph-A, a new physics package, that has been demonstrated to perform well at numerical weather prediction (NWP) and climate scales, are compared against the current global atmosphere configuration and high-resolution results. In addition to an entirely new convection scheme, the package of changes includes significant changes to the cloud, microphysics, and boundary-layer parametrizations. Recognising that CoMorph-A is the first version of a scheme that will continue to be substantially developed and to obtain good performance, compromises in tuning have had to be made. These idealised tests therefore show what works well in this configuration, and what areas will require further work. As such, it is quite a demanding testbed and could be viewed as some of the equipment required for a “convective playground”.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Nanometric modulations of the magnetic structure of the element Nd

The rare earth neodymium arguably exhibits the most complex magnetic ordering and series of magnetic phase transitions of the elements. Here we report the results of small-angle neutron scattering (SANS) measurements as a function of temperature and applied magnetic field to study magnetic correlations on nanometer length scales in Nd. We support these results with additional neutron diffraction, thermodynamic and synchrotron x-ray measurements. The SANS measurements reveal the presence of modulation vectors characterizing the ordered spin configuration which exhibit changes in magnitude and direction that are phase dependent. Between 5.9 and 7.6 K an additional modulation vector is observed with magnitude Q = 0.12 Å -1 . This scattering appears to originate from order of the Nd layers which contain a center of inversion. Tracking this modulation vector as a function of magnetic field indicates a phase boundary at ≈1 T. The modulation vectors observed in this study establish the presence of nanometer length scale spin textures which are likely stabilized by frustrated Ruderman–Kittel–Kasuya–Yosida (RKKY) interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗
DOE OSTI2021

Evolution of non-Kramers doublets in magnetic field in PrNi 2 Cd 20 and PrPd 2 Cd 20

Praseodymium-based 1-2-20 cage compounds PrT 2 X 20 (T is generally Ti, V, Nb, Ru, Rh, Ir; and X is either Al, Zn, or Cd) provide yet another platform to study nontrivial electronic states of matter ranging from topological and magnetic orders to unconventional multipolar orders and superconductivity. In this paper, we report measurements of the electronic heat capacity in two Pr-based 1-2-20 materials: PrNi 2 Cd 20 and PrPd 2 Cd 20 . We find that the lowest-energy multiplet of the Pr 4f 2 valence configuration is a Γ 3 non-Kramers doublet and the first excited triplet is assumed to be a magnetic Γ 5 . By analyzing the dependence of the energy splitting between the ground and first excited singlet states on an external magnetic field, we found that the maximum in the heat capacity corresponding to the Schottky anomaly in PrPd 2 Cd 20 , unlike PrNi 2 Cd 20 , shows a pronounced linear dependence on an external magnetic field at higher field values. In conclusion, this effect is associated with the exchange interactions between the field-induced magnetic dipole moments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

On Simulating the Formation of Structured, Crystalline Systems via Non-classical Pathways

Observations in crystal growth and assembly from recent in situ methods suggest alternative, non-classical crystallization pathways play an important role in the determination of the micro- and meso- structures in crystalline systems. These processes display parallels that cross-cut multiple disciplines investigating crystallization across four orders of magnitude in size scales and widely differing environments, hinting that alternative crystal growth pathways may be a fundamental scheme in natural crystal formation. Using a system of short-range attractive microbeads, we demonstrate that the addition of a small concentration of sub-species incommensurate with the lattice spacing of the dominant species results in a stark change in crystal size and morphology. These changes are attributed to the presence of fleeting, amorphous-like configurations of beads that ultimately change the melting and growth dynamics in preferred directions. From these real-time observations, we hypothesize the amorphous mineral precursors present in biological mineralized tissues undergo similar non-classical crystallization processes resulting in the complex structures found in biomineralization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structure and Desorption Kinetics of Acetonitrile Thin Films on Pt(111) and on Graphene on Pt(111)

Acetonitrile thin films were prepared on Pt(111) and on graphene on Pt(111) using molecular beam techniques. Temperature programed desorption (TPD) experiments of acetonitrile on Pt(111) displayed first-order kinetics for monolayer desorption and zero-order kinetics for multilayer desorption. We used reflection adsorption infrared spectroscopy (RAIRS) to characterize the orientation of the monolayer on Pt(111). Acetonitrile molecules have been shown to adopt a bridging structure where the C=N bond breaks to form two single bonds to underlying platinum atoms. When Tdeposition = 120 K our RAIRS spectra are consistent with this bridging configuration. However, when Tdeposition was lowered to 60 K molecules retained the C=N bond and adopted a vertical orientation with the nitrogen closest to the platinum. On graphene, the TPD experiments show that both the monolayer and multilayer desorb with zero-order kinetics. Our RAIRS spectra suggest that the monolayer consists of acetonitrile molecules oriented with their molecular axes parallel to the graphene substrate. This orientation did not change when Tdeposition was varied from 60 to 120 K. Our results provide a more complete picture of how acetonitrile molecules interact with Pt(111) and with graphene surfaces.

desorption tin films, molecules, thin films, desor↗
DOE OSTI2023

Role of Tropical Cyclones in Determining ENSO Characteristics

El Niño-Southern Oscillation (ENSO) can effectively modulate global tropical cyclone (TC) activity, but the role TCs may play in determining ENSO characteristics remains unclear. Here we investigate the impact of TC winds on ENSO using a suite of Earth system model experiments where we insert TC winds, extracted from a TC-permitting high-resolution simulation, into a low-resolution model configuration with nearly no intrinsic TCs. The presence of TC winds in the model increases ENSO power and shifts ENSO frequency closer to what we observe. TCs lead to an increase of strong to extreme El Niño events seen in observations and not simulated in the low-resolution model without intrinsic TCs, mainly through enhanced zonal advection feedback and thermocline feedback. Our results indicate that TCs play a fundamental role in producing the ENSO characteristics we experience today in the climate system and point to a two-way climatological interaction between TCs and ENSO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assembly of polyelectrolyte star block copolymers at the oil–water interface

To understand and resolve adsorption, reconfiguration, and equilibrium conformations of charged star copolymers, we carried out an integrated experimental and coarse-grained molecular dynamics simulation study of the assembly process at the oil–water interface. This is important to guide development of novel surfactants or amphiphiles for chemical transformations and separations. The star block copolymer consisted of arms that are comprised of hydrophilic–hydrophobic block copolymers that are covalently tethered via the hydrophobic blocks to one point. The hydrophobic core represents polystyrene (PS) chains, while the hydrophilic corona represents quaternized poly(2-vinylpyridine) (P2VP) chains. The P2VP is modeled to become protonated when in contact with an acidic aqueous phase, thereby massively increasing the hydrophilicity of this block, and changing the nature of the star at the oil–water interface. This results in a configurational change whereby the chains comprising the hydrophilic corona are significantly stretched into the aqueous phase, while the hydrophobic core remains solubilized in the oil phase. In the simulations, we followed the kinetics of the anchoring and assembly of the star block copolymer at the interface, monitoring the lateral assembly, and the subsequent reconfiguration of the star via changes in the interfacial tension that varies as the degree-of-protonation increases. At low fractions of protonation, the arm cannot fully partition into the aqueous side of the interface and instead interacts with other arms in the oil phase forming a network near the interface. These insights were used to interpret the non-monotonic dependence of pH with the asymptotic interfacial tension from pendant drop tensiometry experiments and spectral signatures of aromatic stretches seen in vibrational sum frequency generation (SFG) spectroscopy. We describe the relationship of interfacial tension to the star assembly via the Frumkin isotherm, which phenomenologically describes anti-cooperativity in adsorbing stars to the interface due to crowding. Although our model explicitly considers long-range electrostatics, the contribution of electrostatics to interfacial tension is small and brought about by strong counterion condensation at the interface. Finally, these results provide key insights into resolving the adsorption, reconfiguration, and equilibrium conformations of charged star block copolymers as surfactants.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

E3SM‐Arctic: Regionally Refined Coupled Model for Advanced Understanding of Arctic Systems Interactions

Earth system models are essential tools for climate projections, but coarse resolutions limit regional accuracy, especially in the Arctic. Regionally refined meshes (RRMs) enhance resolution in key areas while maintaining computational efficiency. This paper provides an overview of the United States (U.S.) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2.1 with an Arctic RRM, hereafter referred to as E3SMv2.1-Arctic, for the atmosphere (25 km), land (25 km), and ocean/ice (10 km) components. We evaluate the atmospheric component and its interactions with land, ocean, and cryosphere by comparing the RRM (E3SM2.1-Arctic) historical simulations (1950–2014) with the uniform low-resolution (LR) counterpart, reanalysis products, and observational data sets. The RRM generally reduces biases in the LR model, improving simulations of Arctic large-scale mean fields, such as precipitation, atmospheric circulation, clouds, atmospheric river frequency, and sea ice thickness. However, it introduces a seasonally dependent surface air temperature bias, reducing the LR cold bias in summer but enhancing the LR warm bias in winter, which contributes to the underestimated winter sea ice area and volume. Radiative feedback analysis shows similar climate feedback strengths in both model configurations, with the RRM exhibiting a more positive surface albedo feedback and contributing to a stronger surface warming than LR. These findings underscore the importance of high-resolution modeling for advancing our understanding of Arctic climate changes and their broader global impacts, although some persistent biases appear to be independent of model resolution at 10–100 km scales.

Energy Exascale Earth System Model (E3SM)↗
DOE OSTI2022

Cross-beam energy transfer saturation: ion heating and pump depletion

Cross-beam energy transfer (CBET) was measured in two regimes where the energy transfer saturation mechanism was determined by the plasma and laser beam conditions. Here, linear kinetic CBET theory agrees well with the measured energy transfer in all experimental configurations and at all probe beam intensities when accounting for pump depletion and the plasma conditions measured using Thomson-scattering. Simultaneous CBET and Thomson-scattering measurements enabled uncertainties in the plasma conditions to be isolated from CBET theory, which allowed the saturation mechanisms to be identified. In the perpendicular-beam configuration the saturation mode was through ion heating, which resulted from ion trapping in the driven waves and subsequent ion-ion collisional heating. In the co-propagating beam configuration there was minimal ion heating and the saturation mode was through pump depletion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Specificity determinants revealed by the structure of glycosyltransferase Campylobacter concisus PglA

Abstract In selected Campylobacter species, the biosynthesis of N‐linked glycoconjugates via the pgl pathway is essential for pathogenicity and survival. However, most of the membrane‐associated GT‐B fold glycosyltransferases responsible for diversifying glycans in this pathway have not been structurally characterized which hinders the understanding of the structural factors that govern substrate specificity and prediction of resulting glycan composition. Herein, we report the 1.8 Å resolution structure of Campylobacter concisus PglA, the glycosyltransferase responsible for the transfer of N ‐acetylgalatosamine (GalNAc) from uridine 5′‐diphospho‐ N ‐acetylgalactosamine (UDP‐GalNAc) to undecaprenyl‐diphospho‐ N , N ′‐diacetylbacillosamine (UndPP‐diNAcBac) in complex with the sugar donor GalNAc. This study identifies distinguishing characteristics that set PglA apart within the GT4 enzyme family. Computational docking of the structure in the membrane in comparison to homologs points to differences in interactions with the membrane‐embedded acceptor and the structural analysis of the complex together with bioinformatics and site‐directed mutagenesis identifies donor sugar binding motifs. Notably, E113, conserved solely among PglA enzymes, forms a hydrogen bond with the GalNAc C6″‐OH. Mutagenesis of E113 reveals activity consistent with this role in substrate binding, rather than stabilization of the oxocarbenium ion transition state, a function sometimes ascribed to the corresponding residue in GT4 homologs. The bioinformatic analyses reveal a substrate‐specificity motif, showing that Pro281 in a substrate binding loop of PglA directs configurational preference for GalNAc over GlcNAc. This proline is replaced by a conformationally flexible glycine, even in distant homologs, which favor substrates with the same stereochemistry at C4, such as glucose. The signature loop is conserved across all Campylobacter PglA enzymes, emphasizing its importance in substrate specificity.

Vuksanovic, Nemanja↗
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↗