Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordination number”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Correlating 3D Surface Atomic Structure and Catalytic Activities of Pt Nanocrystals

Active sites and catalytic activity of heterogeneous catalysts is determined by their surface atomic structures. However, probing the surface structure at an atomic resolution is difficult, especially for solution ensembles of catalytic nanocrystals, which consist of heterogeneous particles with irregular shapes and surfaces. In this study, we constructed 3D maps of the coordination number (CN) and generalized CN ($\bar{C}\bar{N}$) for individual surface atoms of sub-3 nm Pt nanocrystals. Our results reveal that the synthesized Pt nanocrystals are enclosed by islands of atoms with nonuniform shapes that lead to complex surface structures, including a high ratio of low-coordination surface atoms, reduced domain size of low-index facets, and various types of exposed high-index facets. 3D maps of $\bar{C}\bar{N}$ are directly correlated to catalytic activities assigned to individual surface atoms with distinct local coordination structures, which explains the origin of high catalytic performance of small Pt nanocrystals in important reactions such as oxygen reduction reactions and CO electro-oxidation.

3D atomic structure↗
DOE OSTI2024

Highly efficient semi-hydrogenation in strained ultrathin PdCu shell and the atomic deciphering for the unlocking of activity-selectivity

Excellent ethylene selectivity in acetylene semi-hydrogenation is often obtained at the expense of activity. To break the activity-selectivity trade-off, precise control and in-depth understanding of the three-dimensional atomic structure of surfacial active sites are crucial. Here, we designed a novel Au@PdCu core–shell nanocatalyst featuring diluted and stretched Pd sites on the ultrathin shell (1.6 nm), which showed excellent reactivity and selectivity, with 100% acetylene conversion and 92.4% ethylene selectivity at 122 °C, and the corresponding activity was 3.3 times higher than that of the PdCu alloy. The atomic three-dimensional decoding for the activity-selectivity balance was revealed by combining pair distribution function (PDF) and reverse Monte Carlo simulation (RMC). The results demonstrate that a large number of active sites with a low coordination number of Pd–Pd pairs and an average 3.25% tensile strain are distributed on the surface of the nanocatalyst, which perform a pivotal function in the simultaneous improvement of hydrogenation activity and ethylene selectivity. Our work not only develops a novel strategy for unlocking the linear scaling relation in heterogeneous catalysis but also provides a paradigm for atomic 3D understanding of lattice strain in core–shell nanocatalysts.

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

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Initial hydration process of calcium silicates in Portland cement: A comprehensive comparison from molecular dynamics simulations

As the main components of Portland cement, calcium silicates show substantial differences in their hydration reactivity which have not been fully explained. A comprehensive comparison of the initial hydration process of calcium silicates, namely dicalcium silicate (C{sub 2}S) and tricalcium silicate (C{sub 3}S), was conducted using molecular dynamics simulations. The initial hydration process was divided into three stages using cut-off times of 0.001 ns and 3 ns. The hydration of M{sub 3}-C{sub 3}S (010) was more evident than that of β-C{sub 2}S (100), supported by the hydroxylation degree, radial distribution function, atomic density profile, etc. The coordination number of the surface Ca atoms might be the underlying reason for such a difference. Interactions between cement surfaces and water molecules were mainly characterised by solid OH bonding and Ca-water O bonding. Dissolution of Ca atoms was observed, although quite scarce, while no dissolution of Si atoms was observed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tuning the Interaction between Platinum Single Atoms and Ceria by Zirconia Doping for Efficient Catalytic Ammonia Oxidation

Aiming at the development of an efficient NH 3 oxidation catalyst to eliminate the harmful NH 3 slip from the stationary flue gas denitrification system and diesel exhaust aftertreatment system, a facile ZrO 2 doping strategy was proposed to construct Pt 1 /Ce x Zr 1–x O 2 catalysts with a tunable Pt-CeO 2 interaction strength and Pt–O–Ce coordination environment. According to the results of systematic characterizations, Pt species supported on Ce x Zr 1–x O 2 were mainly in the form of single atoms when x ≥ 0.7, and the strength of the Pt-CeO 2 interaction and the coordination number of Pt–O–Ce bond (CN Pt–O–Ce ) on Pt 1 /Ce x Zr 1–x O 2 showed a volcanic change as a function of the ZrO 2 doping amount. Here, it was proposed that the balance between the reasonable concentration of oxygen defects and limited surface Zr–O x species well accounted for the strongest Pt-CeO 2 interaction and the highest CN Pt–O–Ce on Pt/Ce 0.9 Zr 0.1 O 2 . It was observed that the Pt/Ce 0.9 Zr 0.1 O 2 catalyst exhibited much higher NH 3 oxidation activity than other Pt/Ce x Zr 1–x O 2 catalysts. The mechanism study revealed that the Pt 1 species with the stronger Pt-CeO 2 interaction and higher CN Pt–O–Ce within Pt/Ce 0.9 Zr 0.1 O 2 could better activate NH 3 adsorbed on Lewis acid sites to react with O 2 thus resulting in superior NH 3 oxidation activity. This work provides a new approach for designing highly efficient Pt/CeO 2 based catalysts for low-temperature NH 3 oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Moving beyond the constraints of chemistry via crystal structure discovery with isotropic multiwell pair potentials

The rigid constraints of chemistry-dictated by quantum mechanics and the discrete nature of the atom-limit the set of observable atomic crystal structures. What structures are possible in the absence of these constraints? Here, we systematically crystallize one-component systems of particles interacting with isotropic multiwell pair potentials. We investigate two tunable families of pairwise interaction potentials. Our simulations self-assemble a multitude of crystal structures ranging from basic lattices to complex networks. Sixteen of the structures have natural analogs spanning all coordination numbers found in inorganic chemistry. Fifteen more are hitherto unknown and occupy the space between covalent and metallic coordination environments. The discovered crystal structures constitute targets for self-assembly and expand our understanding of what a crystal structure can look like.

crystal structures↗
DOE OSTI2022

Large-scale frictionless jamming with power-law particle size distributions

Due to significant computational expense, discrete element method simulations of jammed packings of size-dispersed spheres with size ratios greater than 1:10 have remained elusive, limiting the correspondence between simulations and real-world granular materials with large size dispersity. Here, invoking a recently developed neighbor binning algorithm, we generate mechanically stable jammed packings of frictionless spheres with power-law size distributions containing up to nearly 4 000 000 particles with size ratios up to 1:100. By systematically varying the width and exponent of the underlying power laws, we analyze the role of particle size distributions on the structure of jammed packings. The densest packings are obtained for size distributions that balance the relative abundance of large-large and small-small particle contacts. Although the proportion of rattler particles and mean coordination number strongly depend on the size distribution, the mean coordination of nonrattler particles attains the frictionless isostatic value of six in all cases. The size distribution of nonrattler particles that participate in the load-bearing network exhibits no dependence on the width of the total particle size distribution beyond a critical particle size for low-magnitude exponent power laws. This signifies that only particles with sizes greater than the critical particle size contribute to the mechanical stability. However, for high-magnitude exponent power laws, all particle sizes participate in the mechanical stability of the packing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Understanding the Instability of the Halide Perovskite CsPbI 3 through Temperature-Dependent Structural Analysis

Despite the tremendous interest in halide perovskite solar cells, the structural reasons that cause the all-inorganic perovskite CsPbI 3 to be unstable at room temperature remain mysterious, especially since many tolerance-factor-based approaches predict CsPbI 3 should be stable as a perovskite. Here single-crystal X-ray diffraction and X-ray pair distribution function (PDF) measurements characterize bulk perovskite CsPbI3 from 100 to 295 K to elucidate its thermodynamic instability. While Cs occupies a single site from 100 to 150 K, it splits between two sites from 175 to 295 K with the second site having a lower effective coordination number, which, along with other structural parameters, suggests that Cs rattles in its coordination polyhedron. PDF measurements reveal that on the length scale of the unit cell, the Pb-I octahedra concurrently become greatly distorted, with one of the I-Pb-I angles approaching 82° compared to the ideal 90°. The rattling of Cs, low number of Cs-I contacts, and high degree of octahedral distortion cause the instability of perovskite-phase CsPbI 3 . These results reveal the limitations of tolerance factors in predicting perovskite stability and provide detailed structural information that suggests methods to engineer stable CsPbI 3 -based solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of nanoconfinement and surface charge on iron adsorption on mesoporous silica

Here, we present a combined molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) spectroscopic investigation of aqueous iron adsorption on nanoconfined amorphous silica surfaces. The simulation models examine the effects of pore size, pH (surface charge), iron valency, and counter-ion (chloride or hydroxide). The simulation methods were validated by comparing the coordination environment of adsorbed iron with coordination numbers and bond lengths derived from XAFS. In the MD models, nanoconfinement effects on local iron coordination were investigated by comparing results for unconfined silica surfaces and in confined domains within 2 nm, 4 nm, and 8 nm pores. Experimentally, coordination environments of iron adsorbed onto mesoporous silica with 4 nm and 8 nm pores at pH 7.5 were investigated. The effect of pH in the MD models was included by simulating Fe(II) adsorption onto negatively charged SiO 2 surfaces and Fe(III) adsorption on neutral surfaces. The simulation results show that iron adsorption depends significantly on silica surface charge, as expected based on electrostatic interactions. Adsorption on a negatively charged surface is an order of magnitude greater than on the neutral surface, and simulated surface coverages are consistent with experimental results. Pore size effects from the MD simulations were most notable in the adsorption of Fe(II) at deprotonated surface sites (SiO - ), but adsorption trends varied with concentration and aqueous Fe speciation. The coordination environment of adsorbed iron varied significantly with the type of anion. Considerable ion pairing with hydroxide anions led to the formation of oligomeric surface complexes and aqueous species, resulting in larger iron hydroxide clusters at higher surface loadings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Structural and physical properties of Ni 1- x V x alloys around and away from quantum critical point

We investigate the room temperature structure (global and local), temperature dependent magnetic and transport behaviour of Ni 1-x V x (0$⩽$x$ ⩽$ $0.13$) alloys. Our Energy Dispersive Analysis of x-rays results show that the prepared compositions are stoichiometric. With increase in V doping, the compounds exhibit a quantum phase transition around x c = 0.12, where the ferromagnetic phase is suppressed. Further, our results show that all the compounds stabilize in face centred cubic structure at RT and the lattice parameter shows unusual behaviour close to x c . The magnetic and heat capacity studies show signature of Griffiths phase on either side of x c . From 25 K to the lowest collected temperature, we observe a linear T dependence of resistivity at x = 0.1 and around x c , which is separated by a Fermi-liquid region around x = 0.106. This suggests that the origin of the transport behaviour is different around the quantum critical point and away from it. Our Ni K-edge x-ray Absorption Spectroscopy results show that there is a significant reduction in the first coordination number around Ni central atom on doping. Further, with doping, there is distortion in the first coordination shell around Ni. This suggests, with V doping, the local structure around Ni is different from the global structure as obtained from the x-ray Diffraction results. Interestingly, with doping, we observe a direct connection between the extent of distortion at RT and the magnetic disorder obtained at 2 K. We believe our results will motivate the scientific community to further study the interplay between the structural disorder and quantum fluctuations with temperature at the local level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Understanding Metal–Organic Framework Nucleation from a Solution with Evolving Graphs

A mechanistic understanding of metal–organic framework (MOF) synthesis and scale-up remains underexplored due to the complex nature of the interactions of their building blocks. In this work, we investigate the collective assembly of building units at the early stages of MOF nucleation, using MIL-101(Cr) as a prototypical example. Using large-scale molecular dynamics simulations, we observe that the choice of solvent (water and N,N-dimethylformamide), the introduction of ions (Na+ and F–) and the relative populations of MIL-101(Cr) half-secondary building unit (half-SBU) isomers have a strong influence on the cluster formation process. Additionally, the shape, size, nucleation and growth rates, crystallinity, and short and long-range order largely vary depending on the synthesis conditions. We evaluate these properties as they naturally emerge when interpreting the self-assembly of MOF nuclei as the time evolution of an undirected graph. Solution-induced conformational complexity and ionic concentration have a dramatic effect on the morphology of clusters emerging during assembly. While pure solvents lead to the rapid formation of a small number of large clusters, the presence of ions in aqueous solutions results in smaller clusters and slower nucleation. This diversity is captured by the key features of the graph representation. Principle component analysis on graph properties reveals that only a small number of molecular descriptors is needed to deconvolute MOF self-assembly. Furthermore, descriptors such as the average coordination number between half-SBUs and fractal dimension are of particular interest as they can be can be followed experimentally by techniques like by time-resolved spectroscopy. Ultimately, graph theory emerges as an approach that can be used to understand complex processes revealing molecular descriptors accessible by both simulation and experiment.

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

Black carbon enriches short-range-order ferrihydrite in Amazonian Dark Earth: Interplay mechanism and environmental implications

Our study underpins the mechanism of organo-mineral interaction between black carbon (BC, biochar) and associated minerals in the historical BC-rich Amazonian Dark Earth (ADE) by using synchrotron-based microscopic (TXM), microspectroscopic (μFTIR) and spectroscopic (XAS and μ-diffraction) approaches. The BC-rich ADE contained over 100% more poorly crystalline minerals than the adjacent tropical soil. Linear combination fitting of k-spacing in the X-ray Absorption Spectra (XAS) revealed that ferrihydrite contributed to 81.1% of the Fe-minerals in BC. A small but distinct peak was observed at 5.7 Å –1 in the extended X-ray absorption fine structure k oscillation of BC, revealing the presence of Fe—C (including Fe-O-C) covalent bonds. No Fe—C path was yielded by the XAS fitting when an obvious peak downshift of the first (Fe—Fe 1 ) shell was observed, suggesting that the availability of inner-sphere Fe—C complexation was limited to the BC surface and interphase region. In this study, the main minerals for organo-mineral complexation were short-range-order (SRO) ferrihydrite on BC instead of corner-sharing FeO 6 octahedra. Compared to ADE, the coordination number of the first (Fe—Fe 1 ) and second (Fe—Fe 2 ) shell was higher in BC, revealing a higher degree of order in coordination between the neighboring Fe mineral crystals. Black C limited the progressive aging of amorphous Fe phases and greatly enriched SRO ferrihydrite in the redox-fluctuating and high-leaching environment. The transformation of SRO ferrihydrite into the more crystalline Fe oxides was controlled by the local pH environment. A strong signal from the complexed phenolic group (aryl-OH, 1241 cm –1 ) and a distinct band of inner-sphere complexation (Fe-aryl C, 1380–1384 cm –1 ) were identified in the FTIR spectra. The enrichment of poorly crystalline minerals can have positive feedback on the long-term stabilization of BC. The scale-up application of biochar to agricultural and ecological systems may have a long-lasting impact on the enrichment and transformation of the SRO minerals in the soil.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Relaxation Dynamics and Ion Conduction of Poly(Ethylene Carbonate/Ethylene Oxide) Copolymer-Based Electrolytes

Poly(ethylene carbonate/ethylene oxide) P(EC/EO) copolymer was synthesized from EC monomer as a new matrix component of solid polymer electrolytes (SPEs), and the correlation between polymer dynamics and ionic conductivity in the P(EC/EO)-based SPEs at various lithium salt weight fractions (w Li ) was investigated by dielectric spectroscopy, rheology, and molecular dynamics (MD) simulations. With the addition of lithium salt, the molecular motions at various scales (i.e., from local to segmental and global motion scales) were found to change as follows: (i) the local motion of P(EC/EO) was slightly accelerated, and at w Li ≥ 0.22, it was separated into two relaxation modes, one faster and the other slower than that of P(EC/EO), (ii) the segmental motion of P(EC/EO) became significantly slower with the addition of lithium salt, and (iii) the global terminal relaxation became slightly slower with the increase of w Li . Thus, the motion of P(EC/EO) at various scales in SPEs changed with different trends with increasing lithium salt. It was also found that the temperature dependence of the ionic conductivity was described by a Vogel–Fulcher–Tammann-type function, which was correlated with the segmental motion, as is known in other SPE systems. Furthermore, from MD simulations, lithium ions are coordinated with the EC units more likely and with a closer distance than EO units in P(EC/EO), and the coordination number of oxygens around one lithium ion was estimated as 6–7 at the distance of 0.30 nm. Here, these results suggest that the coordination structure of the oxygens around the lithium ion is slightly disrupted due to the presence of EC units compared to the chelating structure of pure PEO electrolytes.

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

Photodissociation and Infrared Spectroscopy of U + (CO 2 ) n , UO + (CO 2 ) n , and UO 2+ (CO 2 ) n Cation-Molecular Complexes

Laser vaporization of uranium in a pulsed supersonic expansion of carbon dioxide is used to produce complexes of the form U + (CO 2 ) n , UO + (CO 2 ) n , and UO 2+ (CO 2 ) n . These ions are selected in a reflectron time-offlight mass spectrometer and studied with visible laser photodissociation and tunable infrared laser photodissociation spectroscopy in the region of the CO 2 antisymmetric stretch. The dissociation patterns and spectroscopy of these ions indicate that CO 2 ligands are intact molecules. Although reaction products that form oxide-carbonyl or oxalate species are predicted to be stable, there is no direct evidence in the frequency range studied for the formation of these species. There is no clear indication for the coordination numbers for singly charged uranium and its oxide complexes with CO 2 . However, there is strong support in the vibrational patterns for an eight-coordinate complex of the doubly charged UO 2+ species, i.e., UO 2+ (CO 2 ) 8 .

Cluster chemistry↗
NASA NTRS1976

Numerical Airfoil Optimization Using a Reduced Number of Design Coordinates

A method is presented for numerical airfoil optimization whereby a reduced number of design coordinates are used to define the airfoil shape. The approach is to define the airfoil as a linear combination of shapes. These basic shapes may be analytically or numerically defined, allowing the designer to use his insight to propose candidate designs. The design problem becomes one of determining the participation of each such function in defining the optimum airfoil. Examples are presented for two-dimensional airfoil design and are compared with previous results based on a polynomial representation of the airfoil shape. Four existing NACA airfoils are used as basic shapes. Solutions equivalent to previous results are achieved with a factor of more than 3 improvements in efficiency, while superior designs are demonstrated with an efficiency greater than 2 over previous methods. With this shape definition, the optimization process is shown to exploit the simplifying assumptions in the inviscid aerodynamic analysis used here, thus demonstrating the need to use more advanced aerodynamics for airfoil optimization.

Vanderplaats, G. N.↗
DOE OSTI2023

Disordered interfaces of alkaline aluminate salt hydrates provide glimpses of Al 3+ coordination changes

Hypothesis: The precipitation and dissolution of aluminum-bearing mineral phases in aqueous systems often proceed via changes in both aluminum coordination number and connectivity, complicating molecular-scale interpretation of the transformation mechanism. Here, the thermally induced transformation of crystalline sodium aluminum salt hydrate, a phase comprised of monomeric octahedrally coordinated aluminate which is of relevance to industrial aluminum processing, has been studied. Because intermediate aluminum coordination states during melting have not previously been detected, it is hypothesized that the transition to lower coordinated aluminum ions occurs within a highly disordered quasi-two-dimensional phase at the solid-solution interface. Experiments and simulations: In this work, in situ X-ray diffraction (XRD), Raman and 27 Al nuclear magnetic resonance (NMR) spectroscopy were used to monitor the melting transition of nonasodium aluminate hydrate (NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6H 2 O). A mechanistic interpretation was developed based on complementary classical molecular dynamics (CMD) simulations including enhanced sampling. A reactive forcefield was developed to bridge speciation in the solution and in the solid phase. Findings: In contrast to classical dissolution, aluminum coordination change proceeds through a dynamically stabilized ensemble of intermediate states in a disordered layer at the solid-solution interface. In both melting and dissolution of NSA, octahedral, monomeric aluminum transition through an intermediate of pentahedral coordination. The intermediate dehydroxylates to form tetrahedral aluminate (Al(OH)$_4^–$) in the liquid phase. This coordination change is concomitant with a breaking of the ionic aluminate-sodium ion linkages. The solution phase Al(OH)$_4^–$ ions subsequently polymerize into polynuclear aluminate ions. However, there are some differences between bulk melting and interfacial dissolution, with the onset of the surface-controlled process occurring at a lower temperature (~30 °C) and the coordination change taking place more gradually as a function of temperature. This work to determine the local structure and dynamics of aluminum in the disordered layer provides a new basis to understand mechanisms controlling aluminum phase transformations in highly alkaline solutions.

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

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.

36 MATERIALS SCIENCE↗