Search

A non-polarizable force field for lithium (Li+) and bis(trifluoromethanesulfonyl)imide (TFSI−) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Overall, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level. [ABSTRACT FROM AUTHOR]

A trace amount of interfacial water is required to initiate hydrosilation reactions of trifunctional organosilanes to form surface assemblies. In recent studies, we have learned that water also has a critical role in directing molecular placement on surfaces because water can react with silicon to provide oxygenated sites for surface binding. Consequently, the wettability nature of substrates influences the placement and density of organosilane films formed by vapor-phase reactions. Nanopatterning protocols were designed using vapor-phase organosilanes and colloidal lithography to compare the wettability differences of hydrophilic mica(0001) compared to relatively hydrophobic Si(100) as a strategy for tracking the location of water on surfaces. The competition between hydrophobic and hydrophilic domains for the adsorption and coalescence of water condensed from vapor can be mapped indirectly by mapping the organosilanes, which bind to water at the solid interface, using atomic force microscopy. Trifunctional octadecyltrichlorosilane (OTS) was used as a marker molecule to map out the areas of the surface where water was deposited. The effect of systematic changes in film thickness and surface coverage of OTS was evaluated at the vapor/solid interface by adding an incremental amount of water to sealed reaction vessels to wet the surface and assessing the outcome after reaction with vapor-phase trichlorosilane. Reactive molecular dynamics simulations of the silicon–water vapor interface combined with electronic structure calculations of oxygenated silicon clusters with methyltrichlorosilane provided insight of the mechanism for surface binding, toward understanding the nature of the interface and wettability factors, which influence the association and placement of silane molecules on surfaces. [ABSTRACT FROM AUTHOR]

Sodium-based rechargeable battery technologies are being pursued as an alternative to lithium, in part due to the relative abundance of sodium compared to lithium. Despite their low dielectric constant, glyme-based electrolytes are particularly attractive for these sodium-based batteries due to their ability to chelate with the sodium ion and their high electrochemical stability. While the glyme chain length is a parameter that can be tuned to modify solvation properties, charge transport behavior, reactivity, and ultimately battery performance, anion identity provides another tunable variable. Trifluoromethanesulfonate (triflate/OTf) and bis(trifluoromethane)sulfonamide (TFSI) are chemically similar anions, which are often used in battery electrolytes for lithium-based batteries. In this paper, molecular simulations are used to examine the differences in ion association and charge transport between sodium salts of these two anions at different salt concentrations in glymes with the increasing chain length. The use of the modified force field developed for NaOTf in glymes for the NaTFSI electrolytes was validated by comparing the TFSI–sodium ion radial distribution functions to the results from ab initio molecular dynamics simulations on 1.5 M NaTFSI in diglyme. While the ion association behavior as a function of salt concentration showed similar trends for both NaOTf and NaTFSI in tetraglyme and triglyme electrolytes, the dominant solvation structures for the two sets of electrolytes are distinctly different in the monoglyme and diglyme cases. The conductivity is impacted by both the ion association behavior in these electrolytes and the non-vehicular or hopping transport of the anions in these systems. [ABSTRACT FROM AUTHOR]

We explore the connection between the solvation dynamics of a chromophore upon photon excitation and equilibrium fluctuations of the solvent. Using molecular dynamics simulations, fluorescence Stokes shift for the tryptophan in Staphylococcus nuclease was examined using both nonequilibrium calculations and linear response theory. When the perturbed and unperturbed surfaces exhibit different solvent equilibrium fluctuations, the linear response approach on the former surface shows agreement with the nonequilibrium process. This agreement is excellent when the perturbed surface exhibits Gaussian statistics and qualitative in the case of an isomerization induced non-Gaussian statistics. However, the linear response theory on the unperturbed surface breaks down even in the presence of Gaussian fluctuations. Experiments also provide evidence of the connection between the excited state solvent fluctuations and the total fluorescence shift. These observations indicate that the equilibrium statistics on the excited state surface characterize the relaxation dynamics of the fluorescence Stokes shift. Our studies specifically analyze the Gaussian fluctuations of the solvent in the complex protein environment and further confirm the role of solvent fluctuations on the excited state surface. The results are consistent with previous investigations, found in the literature, of solutes dissolved in liquids. [ABSTRACT FROM AUTHOR]

We explore the connection between the solvation dynamics of a chromophore upon photon excitation and equilibrium fluctuations of the solvent. Using molecular dynamics simulations, fluorescence Stokes shift for the tryptophan in Staphylococcus nuclease was examined using both nonequilibrium calculations and linear response theory. When the perturbed and unperturbed surfaces exhibit different solvent equilibrium fluctuations, the linear response approach on the former surface shows agreement with the nonequilibrium process. This agreement is excellent when the perturbed surface exhibits Gaussian statistics and qualitative in the case of an isomerization induced non-Gaussian statistics. However, the linear response theory on the unperturbed surface breaks down even in the presence of Gaussian fluctuations. Experiments also provide evidence of the connection between the excited state solvent fluctuations and the total fluorescence shift. These observations indicate that the equilibrium statistics on the excited state surface characterize the relaxation dynamics of the fluorescence Stokes shift. Our studies specifically analyze the Gaussian fluctuations of the solvent in the complex protein environment and further confirm the role of solvent fluctuations on the excited state surface. The results are consistent with previous investigations, found in the literature, of solutes dissolved in liquids. [ABSTRACT FROM AUTHOR]

The role of the solvated excess proton and hydroxide ions in interfacial properties is an interesting scientific question with applications in a variety of aqueous behaviors. The role that charge transfer (CT) plays in interfacial behavior is also an unsettled question. Quantum calculations are carried out on clusters of water with an excess proton or a missing proton (hydroxide) to determine their CT. The quantum results are applied to analysis of multi-state empirical valence bond trajectories. The polyatomic nature of the solvated excess proton and hydroxide ion results in directionally dependent CT, depending on whether a water molecule is a hydrogen bond donor or acceptor in relation to the ion. With polyatomic molecules, CT also depends on the intramolecular bond distances in addition to intermolecular distances. The hydrated proton and hydroxide affect water's liquid/vapor interface in a manner similar to monatomic ions, in that they induce a hydrogen-bonding imbalance at the surface, which results in charged surface waters. This hydrogen bond imbalance, and thus the charged waters at the surface, persists until the ion is at least 10 Å away from the interface. [ABSTRACT FROM AUTHOR]

Zwitterions are naturally occurring molecules that have a positive and a negative charge group in its structure and are of great importance in many areas of science. Here, the vibrational and hydration dynamics of the zwitterionic system betaine (N,N,N-trimethylglycine) is reported. The linear infrared spectrum of aqueous betaine exhibits an asymmetric band in the 1550-1700 cm-1 region of the spectrum. This band is attributed to the carboxylate asymmetric stretch of betaine. The potential of mean force computed from ab initio molecular dynamic simulations confirms that the two observed transitions of the linear spectrum are related to two different betaine conformers present in solution. A model of the experimental data using non-linear response theory agrees very well with a vibrational model comprising of two vibrational transitions. In addition, our modeling shows that spectral parameters such as the slope of the zeroth contour plot and central line slope are both sensitive to the presence of overlapping transitions. The vibrational dynamics of the system reveals an ultrafast decay of the vibrational population relaxation as well as the correlation of frequency-frequency correlation function (FFCF). A decay of ~0.5 ps is observed for the FFCF correlation time and is attributed to the frequency fluctuations caused by the motions of water molecules in the solvation shell. The comparison of the experimental observations with simulations of the FFCF from ab initio molecular dynamics and a density functional theory frequency map shows a very good agreement corroborating the correct characterization and assignment of the derived parameters. [ABSTRACT FROM AUTHOR]