Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulation”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-salinity brine desalination with amine-based temperature swing solvent extraction: A molecular dynamics study

Solvent extraction desalination (SED) is a promising technique for efficiently treating high salinity brines. However, only a few studies in the literature have investigated the influence of the molecular characteristics of the organic solvents on the emergent performance of the separation process. Thus, in this work we investigate the influence of the molecular structure of four different secondary amines on the desalination process of a concentrated aqueous solution of NaCl. Molecular dynamics simulations and free energy calculations are used to obtain the thermodynamic and structural properties of the pure solvent bulk phase, as well as a homogeneous solvent-water bulk mixture. Furthermore, solvent-water and solvent-brine interfacial systems are used to explicitly analyze the water partitioning behavior in the relevant temperature range for SED. The water solubilization in these secondary amines is found to be an exothermic process. Additionally, although the hydrogen-bonding interaction between the water and solvent molecules are beneficial, this interaction is not necessarily the dominant factor for ranking the water solubility performance of the different solvents. Finally, the excess volume of mixing of the bulk solvent-water mixture, as well as the electrostatic surface features of the solvents are also found to provide important clues for screening solvent behavior.

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

Robust and tunable oxide nanoscrolls for solar-driven H 2 generation and storage

Hydrogen gas is a promising alternative to fossil fuels due to its high energy output and environmentally safe byproducts. Various morphologies of photocatalytic materials have been explored for high-efficiency H 2 production, for instance, quasi-1D nanoscroll structures that provide a larger surface-to-volume ratio. Recently, we predicted layer-by-layer formation of stable oxide nanoscrolls directly from dichalcogenide precursors, eliminating the need for costly formation of two-dimensional oxides for a roll-up synthesis of nanoscrolls. Here, in this study, we evaluate the suitability of those oxide nanoscroll materials—MoO 3 , WO 3 , PdO 2 , HfO 2 , and GeO 2 —for solar-driven photocatalytic H 2 production and storage. Using excited state theory coupled with Bethe–Salpeter equation simulations, we discern their electronic and optical properties as a function of interlayer scroll spacing and find them to be highly conducive for solar-driven photocatalysis. Additionally, using ab initio molecular dynamics simulations, we show that they are also suitable for H 2 storage as the nanoscrolls exhibit an effective trapping of hydrogen, even in the presence of defects and vacancies in the oxides. This work thus demonstrates the discovery of robust and tunable oxide nanoscrolls as materials for advancing solar-driven hydrogen technologies.

Chemical compounds↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2020

Addition of Short Polymer Chains Mechanically Reinforces Glassy Poly(2-vinylpyridine)–Silica Nanoparticle Nanocomposites

The addition of hard fillers to a polymer matrix is a well-known process for achieving mechanical reinforcement. With a decrease in the size of the fillers, the contribution from polymer–particle nanometer-sized interfaces become significant, and these interfaces affect the mechanical performance of polymer nanocomposites (PNCs) beyond the limits established for conventional composites. However, the molecular mechanisms underlying the improvement in the mechanical performance of glassy PNCs remain unresolved, necessitating a deeper understanding of the structure–property relationships in these intrinsically heterogeneous systems. In this effort, by using Brillouin light scattering (BLS) and dynamic mechanical analysis (DMA), we demonstrated that adding shorter chains to a PNC prepared with high molecular weight polymers significantly improved the mechanical properties of the PNC in the glassy state. The strongest enhancement of mechanical properties occurred at an optimum concentration of short chains. This is in contrast to behavior of the glass transition temperature of PNCs which shows a monotonic decrease with an increase in the concentration of shorter chains. Using experimental data and coarse-grained molecular dynamics (MD) simulations, we have identified the molecular mechanism leading to the observed non-monotonic changes in mechanical reinforcement. Here, this mechanism includes changes in the nanoscale organization at the interface, combined with chain stretching amplified by the addition of the short chains. Overall, our approach paves a simple, cost-effective pathway to fabricating glassy PNCs with significantly improved mechanical properties that will fill various practical needs.

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

Experimental investigation of abnormal transverse flow enhancement of α particles in heavy-ion collisions

The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Complex Structure of Molten FLiBe (2 Li F – Be F 2 ) Examined by Experimental Neutron Scattering, X-Ray Scattering, and Deep-Neural-Network Based Molecular Dynamics

The use of molten salts as coolants, fuels, and tritium breeding blankets in the next generation of fission and fusion nuclear reactors benefits from furthering the characterization of the molecular structure of molten halide salts, paving the way to predictive capability of the chemical and thermophysical properties of molten salts. Due to its neutronic, chemical, and thermochemical properties, 2 Li F - Be F 2 is a candidate molten salt for several fusion- and fission-reactor designs. We performed neutron and x-ray total-scattering measurements to determine the atomic structure of liquid 2 Li F - Be F 2 . We also performed and neural-network molecular-dynamics simulations to predict the structure obtained by neutron- and x-ray-diffraction experiments. The use of machine learning provides improvements to the efficiency in predicting the structure at a longer length scales than is achievable with simulations at significantly lower computational expense while retaining near accuracy. We found that the NNMD simulations accurately predicted the Be F 4 2 − oligomer formations seen in the experimental first-structure-factor peak. Our combination of high-resolution measurements with large-scale molecular dynamics provided an avenue to explore and experimentally verify the intermediate-range ordering beyond the first-nearest neighbor that has posed too many experimental and computational challenges in previous works. With a deeper understanding of the salt structure and ion ordering, the evolution of salt chemistry over the lifetime of a reactor can be better predicted, which is crucial to the licensing and operation of advanced fission and fusion reactors that employ molten salts. To this end, this work will serve as a reference for future studies of salt structure and macroscopic properties with and without the addition of solutes. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

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

Computational Prediction of Coiled–Coil Protein Gelation Dynamics and Structure

Protein hydrogels represent an important and growing biomaterial for a multitude of applications, including diagnostics and drug delivery. We have previously explored the ability to engineer the thermoresponsive supramolecular assembly of coiled–coil proteins into hydrogels with varying gelation properties, where we have defined important parameters in the coiled–coil hydrogel design. Using Rosetta energy scores and Poisson–Boltzmann electrostatic energies, we iterate a computational design strategy to predict the gelation of coiled–coil proteins while simultaneously exploring five new coiled–coil protein hydrogel sequences. Provided this library, we explore the impact of in silico energies on structure and gelation kinetics, where we also reveal a range of blue autofluorescence that enables hydrogel disassembly and recovery. As a result of this library, we identify the new coiled–coil hydrogel sequence, Q5, capable of gelation within 24 h at 4 °C, a more than 2-fold increase over that of our previous iteration Q2. The fast gelation time of Q5 enables the assessment of structural transition in real time using small-angle X-ray scattering (SAXS) that is correlated to coarse-grained and atomistic molecular dynamics simulations revealing the supramolecular assembling behavior of coiled–coils toward nanofiber assembly and gelation. This work represents the first system of hydrogels with predictable self-assembly, autofluorescent capability, and a molecular model of coiled–coil fiber formation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗
DOE OSTI2026

Solvent effects on triplet yields in BODIPY-based photosensitizers

Here, we employ molecular dynamics simulations and quantum rate theories to elucidate the complex condensed-phase dynamics underpinning triplet-state formation in organic photosensitizers. Using models informed by first-principles calculations complete with a molecular representation of solvents of different polarities, we elucidate the interplay of the internal and environmental interactions underlying triplet yield. We find that triplet yields depend sensitively on the dielectric stabilization of the charge transfer intermediate that facilitates a transition into the triplet manifold. Our results illustrate the importance of molecularly detailed models in understanding the excited-state internal charge-transfer dynamics of photochemically relevant organic molecules.

Coello Escalante, Leonardo [University of Californ↗
DOE OSTI2021

Cs retention and diffusion in C-S-H at different Ca/Si ratio

Cement and concrete are commonly used in the construction of repository sites for radioactive wastes. The correct isolation of those contaminants requires good adsorption and low diffusion rates. Both parameters are highly affected by many factors, such as pH, temperature or composition. The large variability of experimental conditions and formulations makes it extraordinarily difficult to tackle the influence of each of them independently in experimental samples. To this effect, molecular dynamics simulations have been employed in this study to investigate the role of the composition in the capacity to retain Cs and diffusivity of these ions in calcium silicate hydrate (C-S-H) gel pores. The results indicate that the adsorption of Cs ions is worsened at high Ca/Si ratios due to a lower interaction of the cations with the C-S-H surface, while hydrophilicity of the C-S-H nanopore rises, resulting in higher long-range ordering and lower diffusion coefficients of Cs ions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular dynamics study of carbon dioxide and nitrogen selectivity through poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) membrane

Here, a molecular dynamics simulation model was developed to comprehensively understand carbon dioxide over nitrogen (CO 2 /N 2 ) selectivity through polyphosphazene-based membrane comprising of poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) selective layer at the molecular level. The effects of temperature, pressure, and initial feed gas composition on the CO 2 transport on a polymer membrane were studied. The computed free energy and density profile of the permeating gas mixture exhibited that CO 2 molecules express higher interactions with the membrane than N 2 molecules, resulting in higher CO 2 /N 2 selectivity. Statistical analysis of gas molecules (CO 2 , water (H 2 O), and N 2 ) transportation suggested that hydro- and CO 2 -philic functional groups in the membrane significantly impact CO 2 permeability and CO 2 /N 2 selectivity. This study suggested that Lewis acid–base and hydrogen bonding combinations contribute to CO 2 permeation and CO 2 /N 2 selectivity. An equal CO 2 /N 2 selectivity was observed with and without water vapor in the feed gas suggesting that water does not hinder CO 2 transport through the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

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

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

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

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation–polyanion complexation and polyion pair–pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

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

Solvation Dynamics of Wet Ethaline: Water is the Magic Component

The past two decades witnessed the development of a new type of solvent system, named deep eutectic solvents, which have become increasingly investigated because they offer new and potentially favorable properties, such as wide tunability in electrochemical, mechanical, and transport properties. Deep eutectic solvent (DES) systems are composed of at least one main solvent and an additional component that is meant to interrupt the original solvent/solvent interactions, thereby introducing lower melting points relative to each individual component. Ethaline (a 1:2 mole% mixture of choline chloride and ethylene glycol) is one of the most promising DES systems. However, it is also known to be very hygroscopic, which is a constant concern because water absorption during the use of ethaline alters its properties. Within this work, we demonstrate that modest amounts of water addition (1-10%) to ethaline are of little concern for practical use and can even lead to performance improvements, such as accelerated relaxation and solvation. In contrast, very small amounts of <1% of water lead to additional slowing of the solvent response. Thus, we suggest that the attempt to dry ethaline below 1% moisture is rather counter-productive if one attempts to achieve effective solvation and charge transport properties from DESs. This study investigates the effect of water content on the diffusional relaxation dynamics of ethaline. A set of independent spectroscopic experiments and computational simulations are aimed to provide insight into the solvent response of the DES system using femtosecond time-resolved absorption spectroscopy (fs-TA), broadband dielectric spectroscopy (BDS), nuclear magnetic resonance (NMR) diffusometry and broadband relaxometry, and molecular dynamics simulations (MDS) on ethaline with 0, 0.1, 1, 10, and 28.5wt% added water. For dry ethaline, we identify choline chloride as the rate-limiting solvation component in ethaline. However, the role of the solvent components changes gradually as water is added. Further, we provide quantitative solvent relaxation rates using the different presented time-resolved spectroscopic techniques and find remarkable agreement between them. Based on the solvent relaxation rates and combined with MDS, we develop a molecular understanding of the individual solvent components and their interactions in dry and wet ethaline with varying amounts of water content.

25 ENERGY STORAGE↗