Thermodynamics of mixing water with dimethyl sulfoxide, as seen from computer simulations

Structural and Electronic Properties of Novel Azothiophene Dyes: A Multilevel Study Incorporating Explicit Solvation Effects

Among azobenzene derivatives, azothiophenes represent a relatively recent family of compounds that exhibit similar characteristics as dyes and photoreactive systems. Their technological applications are extensive thanks to the additional design flexibility conferred by the heteroaromatic ring. In this study, we present a comprehensive investigation of the structural and electronic properties of novel dyes derived from 3-thiophenamine, utilizing a multilevel approach. We thoroughly examined the potential energy surfaces of the E and Z isomers for three molecules, each bearing different substituents on the phenyl ring at the para position relative to the diazo group. This exploration was conducted through quantum chemistry calculations at various levels of theory, employing a continuum solvent model. Subsequently, we incorporated an explicit solvent (a dimethyl sulfoxide-water mixture) to simulate the most stable isomers using classical molecular dynamics, delivering a clear picture of the local solvation structure and intermolecular interactions. Finally, a hybrid quantum mechanics/molecular mechanics (QM/MM) approach was employed to accurately describe the evolution of the solutes' properties within their environment, accounting for finite temperature effects.

Correlating Substrate Reactivity at Electrified Interfaces with the Electrolyte Structure in Synthetically Relevant Organic Solvent/Water Mixtures

Optimizing electrosynthetic reactions requires fine tuning of a vast chemical space, including charge transfer at electrocatalyst/electrode surfaces, engineering of mass transport limitations, and complex interactions of reactants and products with their environment. Hybrid electrolytes, in which supporting salt ions and substrates are dissolved in a binary mixture of organic solvent and water, represent a new piece of this complex puzzle as they offer a unique opportunity to harness water as the oxygen or proton source in electrosynthesis. In this work, we demonstrate that modulating water-organic solvent interactions drastically impacts the solvation properties of hybrid electrolytes. Combining various spectroscopies with synchrotron small-angle X-ray scattering (SAXS) and force field-based molecular dynamics (MD) simulations, we show that the size and composition of aqueous domains forming in hybrid electrolytes can be controlled. We demonstrate that water is more reactive for the hydrogen evolution reaction (HER) in aqueous domains than when strongly interacting with solvent molecules, which originates from a change in reaction kinetics rather than a thermodynamic effect. We exemplify novel opportunities arising from this new knowledge for optimizing electrosynthetic reactions in hybrid electrolytes. For reactions proceeding first via the activation of water, fine tuning of aqueous domains impacts the kinetics and potentially the selectivity of the reaction. Instead, for organic substrates reacting prior to water, aqueous domains have no impact on the reaction kinetics, while selectivity may be affected. We believe that such a fine comprehension of solvation properties of hybrid electrolytes can be transposed to numerous electrosynthetic reactions.

Promotive effects of marine-derived dimethyl sulfoxide on the photodegradation of phenanthrene in the atmosphere

Dimethyl sulfoxide (DMSO), a versatile medium, is a particular component in the marine atmosphere that possibly causes polycyclic aromatic hydrocarbons (PAHs) to degrade differently than they do in the continental atmosphere. In this study, phenanthrene (Phe) was used as a model PAH in batch photochemical experiments to investigate the chemical actions of DMSO and the underlying mechanisms. The photodegradation of Phe in aqueous solutions with DMSO volume fractions from 0 % to 100 % was initiated by ultraviolet (UV) radiation and promoted by singlet oxygen, which was consistent with pseudo-first-order kinetics. Phe photodegraded faster in a mixture of DMSO and water than in water or DMSO alone, and the rate constant showed a unimodal distribution over the DMSO fraction range, peaking at 33 % DMSO (0.0333 ± 0.0009 min-1) and 40 % DMSO (0.0199 ± 0.0005 min-1) under 254 nm and 302 nm UV radiation, respectively. This interesting phenomenon was attributed to the competition of DMSO for UV radiation and singlet oxygen and changes in dissolved oxygen and free water contents caused by the interaction between DMSO and water molecules. In addition, 9,10-phenanthrenequinone (9,10-PhQ) with high cytotoxicity was the main photodegradation product of Phe under various conditions. The photodegradation rate of Phe in the mixtures of DMSO and water was comparable to its reaction rate with OH radicals, suggesting that 9,10-PhQ can be rapidly generated in the marine atmosphere, driven by a mechanism different from that in the continental or urban atmosphere. Under the presented experimental conditions, UV intensity and DMSO fraction were the primary factors that affected the photodegradation rate of Phe and 9,10-PhQ and altered their integrated toxicity. The findings of this study support the conclusion that the marine atmosphere is an essential field in the atmospheric transport of PAHs, in which DMSO is an important component that affects their photodegradation.

Microwave Sensor for the Determination of DMSO Concentration in Water-DMSO Binary Mixture

This research aims to develop a microwave sensor to accurately measure the concentration of dimethyl sulfoxide (DMSO) in water-DMSO binary mixtures. The proposed sensor will utilize microwave frequency measurements to determine the DMSO concentration, providing a non-invasive and efficient method for analyzing DMSO solutions. The research will involve the design, fabrication, and testing of the sensor, as well as the development of an appropriate calibration model. The outcomes of this study will contribute to improved monitoring and quality control in various fields, including pharmaceuticals, chemical synthesis, and biomedical research. The binary mixtures of dimethyl sulfoxide (DMSO) and water with varying concentrations were investigated in the frequency range of 1 GHz to 5 GHz at room temperature using a microwave sensor. The proposed microwave sensor design was based on an interdigital capacitor (IDC) microstrip antenna loaded with a hexagonal complementary ring resonator (HCRR). The performance of the sensor, fabricated using the print circuit board (PCB) technique, was validated through simulations and experiments. The reflection coefficient (S₁₁) and resonance frequency (F_r) of binary mixtures of DMSO and water solutions were recorded and analyzed for DMSO concentrations ranging from 0% v/v to 75% v/v. Mathematical models were developed to analyze the data, and laboratory tests showed that the sensor can detect levels of DMSO/water binary mixtures. The sensor is capable of detecting DMSO concentrations ranging from 0% v/v to 75% v/v, with a maximum sensitivity of 0.138 dB/% for S₁₁ and ΔS₁₁ and 0.2 MHz/% for F_r and ΔF_r at a concentration of 50% v/v. The developed microwave sensor can serve as an alternative for detecting DMSO concentrations in water using a simple and cost-effective technique. This method can effectively analyze a wide range of concentrations, including highly concentrated solutions, quickly and easily.

Redefining the Molecular Interplay between Dimethyl Sulfoxide, Lipid Bilayers, and Dehydration

The potentially damaging action of dimethyl sulfoxide (DMSO) on phospholipid bilayers remains a matter of controversy. We have conducted a series of long-scale molecular dynamics simulations of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayers at various levels of hydration in the presence of variable quantities of DMSO. These simulations provide evidence for a non-destructive dehydrating mechanism of action for DMSO on DOPC bilayers across a wide concentration range and levels of hydration. Specifically, under full- and low-hydration conditions, the bilayer underwent a minor lateral contraction, coinciding with surface dehydration in the presence of dilute DMSO solutions (X_DMSO < 0.3). At higher DMSO concentrations, this bilayer structure was retained despite a progressive deterioration of the hydration structure at the interface. A similar convergence of bilayer structural properties was observed under dehydration conditions for 0.3 < X_DMSO < 0.7. Destabilization occurred for dehydrated bilayers in the presence of X_DMSO ≥ 0.7, suggesting the existence of a DMSO concentration and/or dehydration threshold. However, such DMSO concentrations far exceed those established as toxic to other cellular components. Our findings represent a computational model for DMSO-DOPC interactions that is consistent with a range of experimental characterizations, offering new molecular insights into the cryoprotective mechanisms of action of DMSO.

Calculation of the Free Energy of Mixing as a Tool for Assessing and Improving Potential Models: The Case of the N,N-Dimethylformamide-Water System

The Helmholtz free energy, energy, and entropy of mixing of N,N-dimethylformamide (DMF) and water are calculated in the entire composition range by means of Monte Carlo computer simulations and thermodynamic integration using all possible combinations of five DMF and three widely used water models. Our results reveal that the mixing of DMF and water is highly non-ideal. Thus, in their dilute solutions, both molecules induce structural ordering of the major component, as evidenced by the concomitant decrease in the entropy. Among the 15 model combinations considered, only 4 reproduce the well-known full miscibility of DMF and water, 3 of which strongly exaggerate the thermodynamic driving force of the miscibility. Thus, the combination of the CS2 model of DMF and the TIP4P/2005 water model reproduces the properties of the DMF-water mixtures far better than the other combinations tested. Our results also reveal that moving a fractional negative charge from the N atom to the O atom of the DMF molecule, leading to the increase in its dipole moment, improves the miscibility of the model with water. Starting from the CS2 model and optimizing the charge to be moved, we propose a new model of DMF that reproduces very accurately both the Helmholtz free energy of mixing of aqueous DMF solutions in the entire composition range (when used in combination with the TIP4P/2005 water model) and also the internal energy of neat DMF.

Thermodynamics of mixing methanol with supercritical CO2 as seen from computer simulations and thermodynamic integration

The changes in extensive thermodynamic quantities, such as volume, energy, Helmholtz free energy and entropy, occurring upon mixing liquid methanol with supercritical CO₂, are calculated using Monte Carlo simulations and thermodynamic integration for all eight combinations of four methanol and two CO₂ potential models in the entire composition range at 313 K. The obtained results are also compared with experimental data whenever possible. The transition of the system from liquid to a supercritical state is found to occur at this temperature around a CO₂ mole fraction value of 0.95 with all model combinations considered. This liquid to supercritical transition is always accompanied by positive Helmholtz free energy of mixing values and, consequently, by the non-miscibility of the two components. Furthermore, both this non-miscibility around the liquid to supercritical transition and also the miscibility of the two components below this transition, in the liquid regime, are found to be primarily of the energetic rather than entropic origin; the entropy of mixing turns out to be very close to zero, and around the liquid to supercritical transition even its qualitative behaviour is strongly model dependent. Finally, it is found that the methanol expansion coefficient is not sensitive to the details of the potential models, and it is always in excellent agreement with the experimental data. On the other hand, both the volume and the energy of mixing depend strongly on the molar volume of neat CO₂ in the model being used, and in this respect the TraPPE model of CO₂ [J. J. Potoff and J. I. Siepmann, AIChE J., 2001, 47, 1676] performs considerably better than that of Zhang and Duan [Z. Zhang and Z. Duan, J. Chem. Phys., 2005, 122, 214507].

Molecular Mechanism of Water Reorientation Dynamics in Dimethyl Sulfoxide Aqueous Mixtures

Nonmonotonic composition dependence is often observed for numerous properties in the aqueous mixtures of small amphiphilic molecules. The molecular picture underlying this structure-activity relationship, however, remains largely elusive. We herein studied water reorientation dynamics in the aqueous mixture of dimethyl sulfoxide (DMSO), which has a significant nonmonotonic composition dependence, using molecular dynamic simulation and an extended molecular jump model. The analysis indicates that this nonideal behavior is driven by the collective frame diffusion component of water reorientation, which decelerates in the water-rich regime because of the strengthened hydrogen bonds and accelerates in the water-poor regime as the hydrogen bonding network is broken into smaller aggregates. The current work therefore connects the microheterogeneity in the solvation structure of DMSO-water with its nonmonotonic hydration dynamics and sheds new light on how microsegregation leads to the multiscale hydration nonideality in general.

All-atom molecular dynamics study on the non-solvent induced phase separation: Thermodynamics of adding water to poly(vinylidene fluoride)/N-methyl-2-pyrrolidone solution

The change in the thermodynamics when adding water in poly(vinylidene fluoride) (PVDF)/N-methyl-2-pyrrolidone (NMP) solution is studied from all atom molecular dynamics (MD) simulations. This is done by estimating the free energy of mixing of PVDF/NMP solution with increasing volume fraction of water (ϕ_w) using an appropriately chosen thermodynamic cycle and the Bennett acceptance ratio method. The MD calculations predict the thermodynamic phase separation point of water/NMP/PVDF to be at ϕ_w = 0.08, in close agreement with the experimental cloud point measurement (ϕ_w = 0.05). Examining the enthalpic and entropic components of the free energy of mixing reveals that at low concentrations of water, the enthalpy term has the most significant contribution to the miscibility of the ternary system, whereas at higher concentrations of water, the entropy term dominates. Finally, the free energy of mixing was compared with the Flory-Huggins (FH) free energy of mixing by computing the concentration-dependent interaction parameters from MD simulations. The FH model inadequately predicted the miscibility of the PVDF solution, mainly due to its negligence of the excess entropy of mixing.

Thermodynamics of Mixing Primary Alkanolamines with Water

The volume, energy, entropy, and Helmholtz free energy of mixing of the seven simplest primary alkanolamine molecules, i.e., monoethanolamine, monoisopropanolamine, 2-amino-propan-1-ol, 2-amino-butan-1-ol, 2-amino-2-methyl-propan-1-ol, 1-amino-2-methyl-propan-2-ol, and 1-amino-butan-2-ol, with water is investigated by extensive computer simulations and thermodynamic integration. To check the force field dependence of the results, all calculations are repeated with two commonly used water models, namely, SPC/E and TIP4P. The obtained results show that the thermodynamics of mixing of alkanolamines and water is largely independent from the type of the alkanolamine molecule. The Helmholtz free energy of mixing is found to be negative for all alkanolamines at every composition, in accordance with the experimentally known full miscibility of these molecules and water. This free energy decrease occurring upon mixing is found to be clearly of energetic origin, as the energy of mixing always turns out to be negative in the entire composition range, while the entropy of mixing is also negative up to high alkanolamine mole fractions. The obtained results suggest that alkanolamines form, on average, stronger hydrogen bonds with water than what is formed by two water molecules, and they induce some ordering of the hydrating water molecules both through the hydrophobic hydration of their side chains and through the strong hydrogen bonding.

Miscibility and Thermodynamics of Mixing of Different Models of Formamide and Water in Computer Simulation

The thermodynamic changes that occur upon mixing five models of formamide and three models of water, including the miscibility of these model combinations itself, is studied by performing Monte Carlo computer simulations using an appropriately chosen thermodynamic cycle and the method of thermodynamic integration. The results show that the mixing of these two components is close to the ideal mixing, as both the energy and entropy of mixing turn out to be rather close to the ideal term in the entire composition range. Concerning the energy of mixing, the OPLS/AA_mod model of formamide behaves in a qualitatively different way than the other models considered. Thus, this model results in negative, while the other ones in positive energy of mixing values in combination with all three water models considered. Experimental data supports this latter behavior. Although the Helmholtz free energy of mixing always turns out to be negative in the entire composition range, the majority of the model combinations tested either show limited miscibility, or, at least, approach the miscibility limit very closely in certain compositions. Concerning both the miscibility and the energy of mixing of these model combinations, we recommend the use of the combination of the CHARMM formamide and TIP4P water models in simulations of water-formamide mixtures.

Accuracy Test of the OPLS-AA Force Field for Calculating Free Energies of Mixing and Comparison with PAC-MAC

We have calculated the excess free energy of mixing of 1053 binary mixtures with the OPLS-AA force field using two different methods: thermodynamic integration (TI) of molecular dynamics simulations and the Pair Configuration to Molecular Activity Coefficient (PAC-MAC) method. PAC-MAC is a force field based quasi-chemical method for predicting miscibility properties of various binary mixtures. The TI calculations yield a root mean squared error (RMSE) compared to experimental data of 0.132 k_BT (0.37 kJ/mol). PAC-MAC shows a RMSE of 0.151 k_BT with a calculation speed being potentially 1.0 × 10⁴ times greater than TI. OPLS-AA force field parameters are optimized using PAC-MAC based on vapor–liquid equilibrium data, instead of enthalpies of vaporization or densities. The RMSE of PAC-MAC is reduced to 0.099 k_BT by optimizing 50 force field parameters. The resulting OPLS-PM force field has a comparable accuracy as the OPLS-AA force field in the calculation of mixing free energies using TI.

Local lateral environment of the molecules at the surface of DMSO-water mixtures

Molecular dynamics simulations of the liquid-vapour interface of dimethyl sulphoxide (DMSO)-water mixtures of 11 different compositions, including two neat systems are performed on the canonical (N, V, T) ensemble at 298 K. The molecules constituting the surface layer of these systems are selected by means of the identification of the truly interfacial molecules (ITIM) method, and their local lateral environment at the liquid surface is investigated by performing Voronoi analysis. The obtained results reveal that both molecules prefer to be in a mixed local environment, consisting of both kinds of molecules, at the liquid surface, and this preference is even stronger here than in the bulk liquid phase. Neat-like patches, in which a molecule is surrounded by like neighbours, are not found. However, vacancies that are surrounded solely by water molecules are observed at the liquid surface. Our results show that strongly hydrogen bonded DMSO·H2O complexes, known to exist in the bulk phase of these mixtures, are absent from the liquid surface.

Force-Field Based Quasi-Chemical Method for Rapid Evaluation of Binary Phase Diagrams

We present the Pair Configurations to Molecular Activity Coefficients (PAC-MAC) method. The method is based on the pair sampling technique of Blanco (Fan, C. F.; Olafson, B. D.; Blanco, M.; Hsu, S. L. Application of Molecular Simulation to Derive Phase Diagrams of Binary Mixtures. Macromolecules 1992, 25, 3667-3676) with an extension that takes the packing of the molecules into account by a free energy model. The intermolecular energy is calculated using classical force fields. PAC-MAC is able to predict activity coefficients and corresponding vapor-liquid equilibrium diagrams at least 4 orders of magnitude faster than molecular simulations. The accuracy of the PAC-MAC method is tested by comparing the results with experimental data and with the results of the COSMO-SAC model (Lin, S.-T.; Sandler, S. I. A Priori Phase Equilibrium Prediction from a Segment Contribution Solvation Model. Ind. Eng. Chem. Res. 2002, 41, 899-913). PAC-MAC (using the OPLS-aa force field) is shown to be comparable in accuracy to COSMO-SAC, at the considerable advantage that PAC-MAC in principle does not require quantum calculation, provided proper force fields to be available.

The local environment of the molecules in water–DMSO mixtures, as seen from computer simulations and Voronoi polyhedra analysis

The local structure of DMSO–water mixtures is studied by computer simulation and Voronoi analysis.