Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations

Specific line shape of the lowest frequency Raman scattering modes of triethylene glycol

For both dielectric spectroscopy and light scattering spectra, the relaxation modes in the microwave region have been characterized by the Debye relaxation model, which is determined by the peak frequency, or by an empirically extended model (e.g., Cole-Davidson and Kohlrausch-Williams-Watts), which has the appropriate line shape. For light scattering from glass-forming liquids, the general line shape is a broader high frequency side in comparison with Debye relaxation. However, for triethylene glycol (TEG) in liquid form at room temperature, the lowest frequency Raman scattering (LFR) mode shows a peak at about 3.0 GHz, which is narrower than that expected for the Debye relaxation. With increasing temperature, this peak exhibits a significant blueshift and begins to resemble the Debye relaxation shape, indicating that the LFR mode of TEG is also a relaxation mode. The narrowing of the LFR mode of TEG is suggested to be caused from the increased non-whiteness of the fluctuation correlations due to increased hydrogen bonding. This is a consequence of breaking the Debye relaxation model's approximation of the overdamping and narrowing limits in the GHz region, which was found in this study by analyzing the relaxation modes of Raman scattering using the multiple random telegraph model for evaluating thermal bath correlation. The analysis results show that the LFR relaxation times of TEG and the main dielectric relaxation overlap only by 333 K. However, the second LFR mode and β-relaxation at higher frequencies coincide over a wide temperature range, suggesting that they are corresponding modes.

Surface Properties of N,N-Dimethylformamide-Water Mixtures, As Seen from Computer Simulations

The liquid-vapor interface of N,N-dimethylformamide (DMF)-water mixtures, spanning the entire composition range, is investigated in detail at 298 K by molecular dynamics simulation and intrinsic surface analysis. DMF molecules are found to adsorb strongly at the liquid surface, but this adsorption extends only to the first molecular layer. Water and DMF molecules mix with each other on the molecular scale even in the surface layer; thus, no marked self-association of any of the components is seen at the liquid surface. The major surface component prefers such orientation in which the molecular dipole vector lays parallel with the macroscopic plane of the surface. On the other hand, the preferred orientation of the minor component is determined, at both ends of the composition range, by the possibility of H-bond formation with the major component. The lack of H-donating ability of DMF leads to a rapid breakup of the percolating H-bond network at the surface; due to the strong adsorption of DMF, this breakup occurs below the bulk phase DMF mole fraction of 0.03. The disruption of the surface H-bond network also accelerates the exchange of both species between the liquid surface and bulk liquid phase, although, for water, this effect becomes apparent only above a bulk phase DMF mole fraction of 0.4. H-bonds formed by a DMF and a water molecule live, on average, 25-60% longer than those formed by two water molecules at the liquid surface. A similar, but smaller (i.e., about 10-20%) difference is seen in the bulk liquid phase. The enhanced surface mobility of the molecules results in 2-6 times larger diffusion coefficient and 2-5 times shorter H-bond lifetime values at the liquid surface than in the bulk liquid phase. The diffusion of both molecules is slowed down in the presence of the other species; in the case of DMF, this effect is caused by the formation of water-DMF H-bonds, whereas for water, steric hindrances imposed by the bulky DMF neighbors are responsible for this slowing down.

Collective Spectroscopy of Solvation Phenomena: Conflicts, Challenges, and Opportunities

Different spectroscopy types reveal different aspects of molecular processes in soft matter. In particular, collective observables can provide insights into intermolecular correlations invisible to the more popular single-particle methods. In this perspective we feature the dielectric relaxation spectroscopy (DRS) with an emphasis on the proper interpretation of this complex observable aided by computational spectroscopy. While we focus on the history and recent advances of DRS in the fields of biomolecular hydration and nanoconfinement, the discussion transcends this particular field and provides a guide for how collective spectroscopy types supported by computational decomposition can be employed to further our understanding of soft matter phenomena.

Local Structure of DMF-Water Mixtures, as Seen from Computer Simulations and Voronoi Analysis

Molecular dynamics simulations of mixtures of N,N-dimethylformamide (DMF) with water of various compositions, covering the entire composition range, are performed on the canonical (N,V,T) ensemble. The local structure of the mixtures is analyzed in terms of radial distribution functions and the contributions of the first five neighbors to them, various order parameters of the water molecules around each other, and properties of the Voronoi polyhedra of the molecules. The analyses lead to the following main conclusions. The two molecules are mixing with each other even on the molecular scale; however, small self-aggregates of both components persist even at their small mole fraction values. In particular, water-water H-bonds exist in the entire composition range, while water clusters larger than 3 and 2 molecules disappear above the DMF mole fraction values of about 0.7 and 0.9, respectively. The O atoms of the DMF molecules can well replace water O atoms in the hydrogen-bonding network. Further, the H-bonding structure is enhanced by the presence of the hydrophobic CH₃ groups of the DMF molecules. On the other hand, the H-bonding network of the molecules gradually breaks down upon the addition of DMF to the system due to the lack of H-donating groups of the DMF molecules. Finally, in neat DMF, the molecules form weak, CH-donated H-bonds with each other; however, these H-bonds disappear upon the addition of water due to the increasing competition with the considerably stronger OH-donated H-bonds DMF can form with the water molecules.

Computer simulation investigation of the adsorption of acetamide on low density amorphous ice. An astrochemical perspective

The adsorption of acetamide on low density amorphous (LDA) ice is investigated by grand canonical Monte Carlo computer simulations at the temperatures 50, 100, and 200 K, characteristic of certain domains of the interstellar medium (ISM). We found that the relative importance of the acetamide-acetamide H-bonds with respect to the acetamide-water ones increases with decreasing temperature. Thus, with decreasing temperature, the existence of the stable monolayer, characterizing the adsorption at 200 K, is gradually replaced by the occurrence of marked multilayer adsorption, preceding even the saturation of the first layer at 50 K. While isolated acetamide molecules prefer to lay parallel to the ice surface to maximize their H-bonding with the surface water molecules, this orientational preference undergoes a marked change upon saturation of the first layer due to increasing competition of the adsorbed molecules for H-bonds with water and to the possibility of their H-bond formation with each other. As a result, molecules stay preferentially perpendicular to the ice surface in the saturated monolayer. The chemical potential value corresponding to the point of condensation is found to decrease linearly with increasing temperature. We provide, in analogy with the Clausius-Clapeyron equation, a thermodynamic explanation of this behavior and estimate the molar entropy of condensed phase acetamide to be 34.0 J/mol K. For the surface concentration of the saturated monolayer, we obtain the value 9.1 ± 0.8 µmol/m², while the heat of adsorption at infinitely low surface coverage is estimated to be -67.8 ± 3.0 kJ/mol. Our results indicate that the interstellar formation of peptide chains through acetamide molecules, occurring at the surface of LDA ice, might well be a plausible process in the cold (i.e., below 50 K) domains of the ISM; however, it is a rather unlikely scenario in its higher temperature (i.e., 100-200 K) domains.

Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study

Picosecond fast motions and their involvement in the biochemical processes such as protein-ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein-ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme-triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues' spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein-ligand binding chemistry and dynamics.

Adsorption of Formamide at the Surface of Amorphous and Crystalline Ices under Interstellar and Tropospheric Conditions. A Grand Canonical Monte Carlo Simulation Study

The adsorption of formamide is studied both at the surface of crystalline (I_h) ice at 200 K and at the surface of low density amorphous (LDA) ice in the temperature range of 50-200 K by grand canonical Monte Carlo (GCMC) simulation. These systems are characteristic of the upper troposphere and of the interstellar medium (ISM), respectively. Our results reveal that while no considerable amount of formamide is dissolved in the bulk ice phase in any case, the adsorption of formamide at the ice surface under these conditions is a very strongly preferred process, which has to be taken into account when studying the chemical reactivity in these environments. The adsorption is found to lead to the formation of multimolecular adsorption layer, the occurrence of which somewhat precedes the saturation of the first molecular layer. Due to the strong lateral interaction acting between the adsorbed formamide molecules, the adsorption isotherm does not follow the Langmuir shape. Adsorption is found to be slightly stronger on LDA than I_h ice under identical thermodynamic conditions, due to the larger surface area exposed to the adsorption. Indeed, the monomolecular adsorption capacity of the LDA and I_h ice surfaces is found to be 10.5 ± 0.7 μmol/m² and 9.4 μmol/m², respectively. The first layer formamide molecules are very strongly bound to the ice surface, forming typically four hydrogen bonds with each other and the surface water molecules. The heat of adsorption at infinitely low surface coverage is found to be -105.6 kJ/mol on I_h ice at 200 K.

Water-like Behavior of Formamide: Jump Reorientation Probed by Extended Depolarized Light Scattering

Water is a strong self-associated liquid with peculiar properties that crucially depend on H-bonding. As regards its molecular dynamics, only recently has water reorientation been successfully described based on a jump mechanism, which is responsible for the overall H-bonding exchange. Here, using high-resolution broad-band depolarized light scattering, we have investigated the reorientational dynamics of formamide (FA) as a function of concentration from the neat liquid to diluted aqueous solutions. Our main findings indicate that in the diluted regime the water rearrangement can trigger the motion of FA solute molecules, which are forced to reorient at the same rate as water. This highlights an exceptional behavior of FA, which perfectly substitutes water within its network. Besides other fundamental implications connected with the relevance of FA, its water-like behavior provides rare experimental evidence of a solute whose dynamics is completely slaved to the solvent.

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.

Aqueous solvation of amphiphilic molecules by extended depolarized light scattering: the case of trimethylamine-N-oxide

Hydrophilic and hydrophobic interactions strongly affect the solvation dynamics of biomolecules. To understand their role, small model systems are generally employed to simplify the investigations. In this study the amphiphile trimethylamine N-oxide (TMAO) is chosen as an exemplar, and studied by means of extended frequency range depolarized light scattering (EDLS) experiments as a function of solute concentration. This technique proves to be a suitable tool for investigating different aspects of aqueous solvation, being able at the same time to provide information about relaxation processes and vibrational modes of solvent and solute. In the case study of TMAO, we find that the relaxation dynamics of hydration water is moderately retarded compared to the bulk, and the perturbation induced by the solute on surrounding water is confined to the first hydration shell. The results highlight the hydrophobic character of TMAO in its interaction with water. The number of molecules taking part in the solvation process decreases as the solute concentration increases, following a trend consistent with the hydration water-sharing model, and suggesting that aggregation between solute molecules is negligible. Finally, the analysis of the resonant modes in the THz region and the comparison with the corresponding results obtained for the isosteric molecule tert-butyl alcohol (TBA) allow us to provide new insights into the different solvating properties of these two biologically relevant molecules.

Hydration shells of proteins probed by depolarized light scattering and dielectric spectroscopy: orientational structure is significant, positional structure is not

Water interfacing hydrated proteins carry properties distinct from those of the bulk and is often described as a separate entity, a "biological water." We address here the question of which dynamical and structural properties of hydration water deserve this distinction. The study focuses on different aspects of the density and orientational fluctuations of hydration water and the ability to separate them experimentally by combining depolarized light scattering with dielectric spectroscopy. We show that the dynamics of the density fluctuations of the hydration shells reflect the coupled dynamics of the solute and solvent and do not require a special distinction as "biological water." The orientations of shell water molecules carry dramatically different physics and do require a separation into a sub-ensemble. Depending on the property considered, the perturbation of water's orientational structure induced by the protein propagates 3-5 hydration shells into the bulk at normal temperature.