Hydrogen bond, electron donor-acceptor dimer, and residence dynamics in supercritical CO(2)-ethanol mixtures and the effect of hydrogen bonding on single reorientational and translational dynamics: A molecular dynamics simulation study

Role of Terminal Groups of cis-1,4-Polyisoprene Chains in the Formation of Physical Junction Points in Natural Rubber

The terminal structures of cis-1,4-polyisoprene (PI) chains play a vital role in the excellent comprehensive performance of Hevea natural rubber (NR) with properties such as high toughness, tear-resistance, and wet skid resistance. The cis-1,4-polyisoprene chain constituting NR exhibits a distinct composition of terminal groups comprising two distinct types, namely, the ω and α terminal groups. The structures of the ω terminal [dimethyl allyl (DMA)-(trans-1,4-isoprene)2] and six kinds of α end groups of the polymer chain of NR have been explored by utilizing a newly developed 2D NMR method. In the present work, we examine different kinds of PI melt systems, and we choose various combinations of terminal groups: Hydrogen, one DMA unit with two trans isoprene units as ω end groups and ester-terminated isopentene (α1), hydroxy-terminated isopentene (α2), ester-terminated isobutane (α3), hydroxy-terminated isobutane (α4), ester-terminated 1,4-cis-isoprene (α5), and hydroxy-terminated 1,4-cis-isoprene (α6), i.e., HPIH (PI0)-pure PI (Hydrogen terminal), ωPIα1 (PII), ωPIα2 (PIII), ωPIα3 (PIIII), ωPIα4 (PIIV), ωPIα5 (PIV), and ωPIα6 (PIVI). We evaluated dynamic and static properties of PI chains such as the end-to-end vector autocorrelation function (C(t)), its average relaxation time (τ), end-to-end distance (Ree), and radius of gyration (Rg). We also estimated the diffusion coefficients of polyisoprene chains and pair correlation functions [radial distribution functions (RDFs)], potentials of mean force (PMFs) in between end residues, and survival probability (P(τ)) of end groups around the end group by analyzing the equilibrated trajectories of full-atom MD simulations. As per the examination of C(t), rotational relaxation time τ, and RDFs, we discovered that the existence of a strong hydrogen bond in α2-α2, α4-α4, and α6-α6 residues makes the dynamics of hydroxy-terminated polyisoprene chains in ωPIα2,α4,α6 melt systems slower. From the analyses of RDFs and PMFs (W(r)), the association between [α2]-[α2], [α4]-[α4], and [α6]-[α6] terminals in ωPIα2,α4,α6 melt systems is significantly stronger than in [ISO]-[ISO] [Hydrogen terminated 1,4-cis-isoprene:(ISO)] in HPIH and ω-ω, [α1]-[α1], [α3]-[α3], and [α5]-[α5] in ωPIα1,α3,α5 systems. We quantified the fraction of cluster formation of terminal groups of a given size in the seven PI melt systems by employing the criteria of PMFs. It is revealed that no stable cluster exists in the HPIH, ωPIα1, ωPIα3, and ωPIα5 melt systems. Conversely, in the ωPIα2, ωPIα4, and ωPIα6 systems, we perceived stable clusters of [(α2)p] [(α4)p] and [(α6)p] end groups where p (2 ≤ x ≤ 6). These stable clusters validate the presence of physical junction points in between hydroxy-terminated polyisoprene chains through their α2, α4, and α6 terminals. These physical junction points might be crucial for superior properties of NR such as high toughness, crack growth resistance, and strain-induced crystallization.

Understanding the performance of graphdiyne membrane for the separation of nitrate ions from aqueous solution at the atomistic scale

A molecular dynamics simulation study is conducted to investigate the capability of the pristine graphdiyne nanosheet for nitrate ion separation from water. The removal of nitrate ion contaminants from water is of critical importance as it represents an environmental hazard. The graphdiyne is a carbon-based membrane with pore density of 2.4 × 10¹⁸ pores/m² and incircle radius of 2.8 Å. We show that the efficient water flow is accurately controlled through fine regulation of the exerted hydrostatic pressure. The high water permeability of 6.19 L.Day^-1cm^-2MPa^-1 with 100% nitrate ions rejection suggests that the graphdiyne can perform as a suitable membrane for nitrate separation. The potential of mean force analysis of the single water molecule and nitrate ion indicated the free energy barriers for nitrate of about 4 times higher than that of water molecules. The results reveal the weak interaction of the water molecules and the membrane which aid to high water flux.

Diffusion of the carbon dioxide-ethanol mixture in the extended critical region

The effect of traces of ethanol in supercritical carbon dioxide on the mixture's thermodynamic properties is studied by molecular simulations and Taylor dispersion measurements. This mixture is investigated along the isobar p = 10 MPa in the temperature range between T = 304 and 343 K. Along this path, the mixture undergoes two transitions: First, the Widom line is crossed, marking the transition from liquid-like to gas-like conditions. A second transition occurs from the supercritical gas-like domain to a subcritical gas. The Widom line crossover entails inflection points for most of the studied properties, i.e. density, enthalpy, shear viscosity, Maxwell-Stefan and intradiffusion coefficients. On the other hand, the transition between the super- and subcritical regions is found to be generally smooth, an observation that is qualitatively confirmed by experimental Taylor dispersion measurements. Dedicated atomistic simulations show the presence of microheterogeneities due to ethanol self-association along the investigated path, which lead to the mixture's anomalous behavior in its extended critical region.

Computational development and validation of a representative MDI-BDO-based polyurethane hard segment model

Segmented polyurethanes show extraordinary physicochemical properties, mainly owing to the nature and the chemistry of the hard segment domains. There are yet many inexplicable physiochemical properties of MDI-BDO-based hard polyurethane segments such as the geometry, cis-trans isomerism, electronic structure, chemical reactivity, the inter-hard-segment interactions, and the photo-response. In the present study, it was attempted to develop and validate a model system that would facilitate further research on the structural and chemical properties of the MDI-BDO hard segments. It was found that the trans isomer of urethane bond is more stable than the cis isomer, and it is argued here that thermal transformation from trans to cis not possible due to the high rotational energy barrier. The differences between the calculated IR spectra of the cis and trans isomers are proposed as a powerful differentiation tool. The calculated Fukui indices show that cis and trans isomers are different in their chemical reactivity. The findings of the present study suggest intermolecular and intramolecular pi-stacking and highly plausible two significant types of hydrogen bond types between hard segments. In the present study, a model system for MDI-BDO hard segment was developed and successfully validated via computational experiments. Further calculations done with the new model provided an indispensable understanding of the structure, cis-trans isomerism, reactivity, and intermolecular interactions of the MDI-BDO hard segments. The proposed model can be further improved in the future by incorporating suitable soft segments. In summary, the model system developed and validated in the present study has provided new opportunities to understand and further study the structural and chemical features of the hard segments of the MDI-BDO-based polyurethane.

Growth Pattern, Stability, and Properties of Complexes of C2H5OH and nCO2 (n = 1-5) Molecules: A Theoretical Study

This work is dedicated to theoretically investigate the formation process of C2H5OH···nCO2 (n = 1-5) complexes and to shed light on the nature of interactions formed under the variation of CO2 concentration. It is found that CO2 molecules tend to locate around the polarized -OH group to interact with the lone pairs of the O atom. The interaction of ethanol with three CO2 molecules (C2H5OH···3CO2) induces the most stable structure in the sequence considered. The atoms in molecules (AIM), NCIplot, and natural bond orbital (NBO) analyses point out that the Oethanol···CCO2 tetrel bond overcomes hydrogen, chalcogen, and CO2···CO2 tetrel-bonded interactions and mainly contributes to the strength of C2H5OH···nCO2 (n = 1-5) complexes. All intermolecular interactions in the examined complexes are weakly noncovalent, and their positive cooperativity is evaluated to be slightly weaker than that of CO2 pure systems. SAPT2+ and molecular electrostatic potential (MEP) calculations indicate that the electrostatic force is the main factor underlying the attractive interplay in the complexes of C2H5OH and CO2.

MD simulation study of the diffusion and local structure of n-alkanes in liquid and supercritical methanol at infinite dilution

The diffusion coefficients of 14 n-alkanes (ranging from methane to n-tetradecane) in liquid and supercritical methanol at infinite dilution (at a pressure of 10.5 MPa and at temperatures of 299 K and 515 K) were deduced via molecular dynamics simulations. Values for the radial distribution function, coordination number, and number of hydrogen bonds were then calculated to explore the local structure of each fluid. The flexibility of the n-alkane (as characterized by the computed dihedral distribution, end-to-end distance, and radius of gyration) was found to be a major influence and hydrogen bonding to be a minor influence on the local structure. Hydrogen bonding reduces the flexibility of the n-alkane, whereas increasing the temperature enhances its flexibility, with temperature having a greater effect than hydrogen bonding on flexibility. Graphical abstract The flexibility of the alkane is a major influence and the hydrogen bonding is a minor influence on the first solvation shell; the coordination numbers of long-chain n-alkanes in the first solvation shell are rather low.

Structure-Property Relationships in CO2-philic (Co)polymers: Phase Behavior, Self-Assembly, and Stabilization of Water/CO2 Emulsions

This Review provides comprehensive guidelines for the design of CO2-philic copolymers through an exhaustive and precise coverage of factors governing the solubility of different classes of polymers. Starting from computational calculations describing the interactions of CO2 with various functionalities, we describe the phase behavior in sc-CO2 of the main families of polymers reported in literature. The self-assembly of amphiphilic copolymers of controlled architecture in supercritical carbon dioxide and their use as stabilizers for water/carbon dioxide emulsions then are covered. The relationships between the structure of such materials and their behavior in solutions and at interfaces are systematically underlined throughout these sections.

Dynamic and thermodynamic characteristics associated with the glass transition of amorphous trehalose-water mixtures

The glass transition temperature Tg of biopreservative formulations is important for predicting the long-term storage of biological specimens. As a complementary tool to thermal analysis techniques, which are the mainstay for determining Tg, molecular dynamics simulations have been successfully applied to predict the Tg of several protectants and their mixtures with water. These molecular analyses, however, rarely focused on the glass transition behavior of aqueous trehalose solutions, a subject that has attracted wide scientific attention via experimental approaches. Important behavior, such as hydrogen-bonding dynamics and self-aggregation has yet to be explored in detail, particularly below, or in the vicinity of, Tg. Using molecular dynamics simulations of several dynamic and thermodynamic properties, this study reproduced the supplemented phase diagram of trehalose-water mixtures (i.e., Tg as a function of the solution composition) based on experimental data. The structure and dynamics of the hydrogen-bonding network in the trehalose-water systems were also analyzed. The hydrogen-bonding lifetime was determined to be an order of magnitude higher in the glassy state than in the liquid state, while the constitution of the hydrogen-bonding network exhibited no noticeable change through the glass transition. It was also found that trehalose molecules preferred to form small, scattered clusters above Tg, but self-aggregation was substantially increased below Tg. The average cluster size in the glassy state was observed to be dependent on the trehalose concentration. Our findings provided insights into the glass transition characteristics of aqueous trehalose solutions as they relate to biopreservation.

Distinctly Different Glass Transition Behaviors of Trehalose Mixed with Na2HPO 4 or NaH 2PO 4: Evidence for its Molecular Origin

PURPOSE

The present study is aimed at understanding how the interactions between sugar molecules and phosphate ions affect the glass transition temperature of their mixtures, and the implications for pharmaceutical formulations.

METHODS

The glass transition temperature (Tg) and the α-relaxation temperature (Tα) of dehydrated trehalose/sodium phosphate mixtures (monobasic or dibasic) were determined by differential scanning calorimetry and dynamic mechanical analysis, respectively. Molecular dynamics simulations were also conducted to investigate the microscopic interactions between sugar molecules and phosphate ions. The hydrogen-bonding characteristics and the self-aggregation features of these mixtures were quantified and compared.

RESULTS

Thermal analysis measurements demonstrated that the addition of NaH2PO4 decreased both the glass transition temperature and the α-relaxation temperature of the dehydrated trehalose/NaH2PO4 mixture compared to trehalose alone while both Tg and Tα were increased by adding Na2HPO4 to pure trehalose. The hydrogen-bonding interactions between trehalose and HPO4(2-) were found to be stronger than both the trehalose-trehalose hydrogen bonds and those formed between trehalose and H2PO4(-). The HPO4(2-) ions also aggregated into smaller clusters than H2PO4(-) ions.

CONCLUSIONS

The trehalose/Na2HPO4 mixture yielded a higher T g than pure trehalose because marginally self-aggregated HPO4(2-) ions established a strengthened hydrogen-bonding network with trehalose molecules. In contrast H2PO4(-) ions served only as plasticizers, resulting in a lower Tg of the mixtures than trehalose alone, creating large-sized ionic pockets, weakening interactions, and disrupting the original hydrogen-bonding network amongst trehalose molecules.

Polymerization effect of electrolytes on hydrogen-bonding cryoprotectants: ion-dipole interactions between metal ions and glycerol

Protectants which are cell membrane permeable, such as glycerol, have been used effectively in the cryopreservation field for a number of decades, for both slow cooling and vitrification applications. In the latter case, the glass transition temperature (T_g) of the vitrification composition is key to its application, dictating the ultimate storage conditions. It has been observed that the addition of some electrolytes to glycerol, such as MgCl₂, could elevate the T_g of the mixture, thus potentially providing more storage condition flexibility. The microscopic mechanisms that give rise to the T_g-enhancing behavior of these electrolytes are not yet well understood. The current study focuses on molecular dynamics simulation of glycerol mixed with a variety of metal chlorides (i.e., NaCl, KCl, MgCl₂, and CaCl₂), covering a temperature range that spans both the liquid and glassy states. The characteristics of the ion–dipole interactions between metal cations and hydroxyl groups of glycerol were analyzed. The interruption of the original hydrogen-bonding network among glycerol molecules by the addition of ions was also investigated in the context of hydrogen-bonding quantity and lifetime. Divalent metal cations were found to significantly increase the T_g by strengthening the interacting network in the electrolyte/glycerol mixture via strong cation–dipole attractions. In contrast, monovalent cations increased the T_g insignificantly, as the cation–dipole attraction was only slightly stronger than the original hydrogen-bonding network among glycerol molecules. The precursor of crystallization of NaCl and KCl was also observed in these compositions, potentially contributing to weak T_g-enhancing ability. The T_g-enhancing mechanisms elucidated in this study suggest a structure-enhancing role for divalent ions that could be of benefit in the design of protective formulations for biopreservation purposes.

Assessing the non-ideality of the CO2-CS2 system at molecular level: a Raman scattering study

The dense phase of CO2-CS2 mixtures has been analysed by Raman spectroscopy as a function of the CO2 concentration (0.02-0.95 mole fractions) by varying the pressure (0.5 MPa up to 7.7 MPa) at constant temperature (313 K). The polarised and depolarised spectra of the induced (ν2, ν3) modes of CS2 and of the ν1-2ν2 Fermi resonance dyad of both CO2 and CS2 have been measured. Upon dilution with CO2, the evolution of the spectroscopic observables of all these modes displays a "plateau-like" region in the CO2 mole fraction 0.3-0.7 never previously observed in CO2-organic liquids mixtures. The bandshape and intensity of the induced modes of CS2 are similar to those of pure CS2 up to equimolar concentration, after which variations occur. The preservation of the local ordering from pure CS2 to equimolar concentration together with the non-linear evolution of the spectroscopic observables allows inferring that two solvation regimes exist with a transition occurring in the plateau domain. In the first regime, corresponding to CS2 concentrated mixtures, the liquid phase is segregated with dominant CS2 clusters, whereas, in the second one, CO2 monomers and dimers and CO2-CS2 hetero-dimers coexist dynamically on a picosecond time-scale. It is demonstrated that the subtle interplay between attractive and repulsive interactions which provides a molecular interpretation of the non-ideality of the CO2-CS2 mixture allows rationalizing the volume expansion and the existence of the plateau-like region observed in the pressure-composition diagram previously ascribed to the proximity of an upper critical solution temperature at lower temperatures.