Molecular Modeling of Phase Behavior and Microstructure of Acetone−Chloroform−Methanol Binary Mixtures

Synergistic Chloroform-Methanol Binary Solvent Mixture Is Inherently Spatially and Dynamically Heterogeneous

Nonideality in a binary solvent mixture is manifested through anomalies in various physical properties like viscosity, dielectric constant, polarity, freezing point, boiling point, and so forth. Sometimes, such anomalies become much more prominent, leading to a synergistic behavior, where the physical property of the mixture is way different from its bulk counterparts. Various alcohols/chlorinated methane binary solvent mixtures show such a synergistic behavior. The reason is attributed to the unique but diverse interactions present in the system. We speculated that these diverse interactions must manifest heterogeneity in such a binary solvent mixture. Using the improved methodology developed by our group, we investigate the presence of dynamic and spatial heterogeneity in the chloroform/methanol synergistic binary solvent mixture. To our delight, we found that our projection is accurate, and indeed, the chloroform/methanol binary solvent mixtures are heterogeneous. Two maxima for the synergistic behavior have been observed for the chloroform/methanol binary solvent mixture (at ∼0.45 and 0.75 mole fractions of methanol in chloroform) in the literature, where the extent of heterogeneity was also found to be the highest. The present study portrays the intriguing complexity of simple binary solvent mixtures, and the findings may provide valuable insights into solvent engineering for diverse applications like extraction/purification media, reaction media, polymer processing, nanomaterial synthesis, pollutant extraction, active ingredient delivery, biofuel production, and battery technology.

Integrating molecular dynamics simulations and experimental data for azeotrope predictions in binary mixtures

An azeotrope is a constant boiling point mixture, and its behavior is important for fluid separation processes. Predicting azeotropes from atomistic simulations is difficult due to the complexities and convergence problems of Monte Carlo and free-energy perturbation techniques. Here, we present a methodology for predicting the azeotropes of binary mixtures, which computes the compositional dependence of chemical potentials from molecular dynamics simulations using the S0 method and employs experimental boiling point and vaporization enthalpy data. Using this methodology, we reproduce the azeotropes, or lack thereof, in five case studies, including ethanol/water, ethanol/isooctane, methanol/water, hydrazine/water, and acetone/chloroform mixtures. We find that it is crucial to use the experimental boiling point and vaporization enthalpy for reliable azeotrope predictions, as empirical force fields are not accurate enough for these quantities. Finally, we use regular solution models to rationalize the azeotropes and reveal that they tend to form when the mixture components have similar boiling points and strong interactions.

Elucidating the structure of donor-acceptor conjugated polymer aggregates in liquid solution

High-spin donor-acceptor conjugated polymers are extensively studied for their potential applications in magnetic and spintronic devices. Inter-chain charge transfer among these high-spin polymers mainly depends on the nature of the local structure of the thin film and π-stacking between the polymer chains. However, the microscopic structural details of high-spin polymeric materials are rarely studied with an atomistic force field, and the molecular-level local structure in the liquid phase remains ambiguous. Here, we have examined the effects of oligomer chain length, side chain, and processing temperature on the organization of the high-spin cyclopentadithiophene-benzobisthiadiazole donor-acceptor conjugated polymer in chloroform solvent. We find that the oligomers display ordered aggregates whose structure depends on their chain length, with an average π-stacking distance of 3.38 ± 0.03 Å (at T = 298 K) in good agreement with the experiment. Also, the oligomers with longer alkyl side chains show better solvation and a shorter π-stacking distance. Furthermore, the clusters grow faster at higher temperature with more ordered aggregation between the oligomer chains.

Identical Diffusion Distributions and Co-Cluster Formation Dictate Azeotrope Formation: Microscopic Evidences and Experimental Signatures

What selects azeotropic pairs and governs the azeotropic conditions (composition and temperature) is an open and intriguing question. A combined simulation and experimental work presented here investigates this by considering ethanol-water mixtures. We find identical distributions of center-of-mass diffusion coefficients for ethanol and water molecules under the azeotropic condition (95.5 wt % ethanol +4.5 wt % water, Tazeo = 351.1K). Moreover, the particle displacements show strong interspecies correlations at Tazeo. Interestingly, simulated reorientation time distributions become identical at Tazeo but at a composition different from that at which the translational diffusion distributions overlapped. Cluster analyses indicate that solutions at Tazeo with xwater ≤ 15 wt % are more microheterogeneous than those with higher water content, although no anomaly in the composition-dependent solution structural properties was detected. Ethanol-water and ethanol-ethanol interaction energies show pronounced nonideal composition dependence, but the size of the relative fluctuations in them remained small (∼0.5kBT). Rare water-water H-bonding, predominant water-ethanol H-bonding, and a sizable population of "free" water molecules characterize the azeotropic solutions. The red edge excitation spectroscopic (REES) measurements with a dissolved anionic fluorescent dye, coumarin343 (C343), support the predicted solution microheterogeneity by showing a nonmonotonic composition dependence of the excitation energy-induced changes in the fluorescence emission spectral frequencies and bandwidths, the largest changes being under the azeotropic condition. Subsequent dynamic anisotropy measurements reveal a nonmonotonic composition dependence of C343 rotation times with a peak under the azeotropic condition. In summary, equalization of the component translational diffusion coefficients and solution microheterogeneity with regular composition dependence of the solution structure appear to characterize the ethanol-water azeotrope.

Curvature-Selective Nanocrystal Surface Ligation Using Sterically-Encumbered Metal-Coordinating Ligands

Organic ligands are critical in determining the physiochemical properties of inorganic nanocrystals. However, precise nanocrystal surface modification is extremely difficult to achieve. Most research focuses on finding ligands that fully passivate the nanocrystal surface, with an emphasis on the supramolecular structure generated by the ligand shell. Inspired by molecular metal-coordination complexes, we devised an approach based on ligand anchoring groups that are flanked by encumbering organic substituents and are chemoselective for binding to nanocrystal corner, edge, and facet sites. Through experiment and theory, we affirmed that the surface-ligand steric pressures generated by these organic substituents are significant enough to impede binding to regions of low nanocurvature, such as nanocrystal facets, and to promote binding to regions of high curvature such as nanocrystal edges.

Improving Force Field Accuracy by Training against Condensed-Phase Mixture Properties

Developing a sufficiently accurate classical force field representation of molecules is key to realizing the full potential of molecular simulations as a route to gaining a fundamental insight into a broad spectrum of chemical and biological phenomena. This is only possible, however, if the many complex interactions between molecules of different species in the system are accurately captured by the model. Historically, the intermolecular van der Waals (vdW) interactions have primarily been trained against densities and enthalpies of vaporization of pure (single-component) systems, with occasional usage of hydration free energies. In this study, we demonstrate how including physical property data of binary mixtures can better inform these parameters, encoding more information about the underlying physics of the system in complex chemical mixtures. To demonstrate this, we retrain a select number of Lennard-Jones parameters describing the vdW interactions of the OpenFF 1.0.0 (Parsley) fixed charge force field against training sets composed of densities and enthalpies of mixing for binary liquid mixtures as well as densities and enthalpies of vaporization of pure liquid systems and assess the performance of each of these combinations. We show that retraining against the mixture data improves the force field's ability to reproduce mixture properties, including solvation free energies, correcting some systematic errors that exist when training vdW interactions against properties of pure systems only.

Nanoformulation-by-design: an experimental and molecular dynamics study for polymer coated drug nanoparticles

The formulation of drug compounds into nanoparticles has many potential advantages in enhancing bioavailability and improving therapeutic efficacy. However, few drug molecules will assemble into stable, well-defined nanoparticulate structures. Amphiphilic polymer coatings are able to stabilise nanoparticles, imparting defined surface properties for many possible drug delivery applications. In the present article we explore, both experimentally and in silico, a potential methodology to coat drug nanoparticles with an amphiphilic co-polymer. Monomethoxy polyethylene glycol-polycaprolactone (mPEG-b-PCL) diblock copolymers with different mPEG lengths (M w 350, 550, 750 and 2000), designed to give different levels of colloidal stability, were used to coat the surface of indomethacin nanoparticles. Polymer coating was achieved by a flow nanoprecipitation method that demonstrated excellent batch-to-batch reproducibility and resulted in nanoparticles with high drug loadings (up to 78%). At the same time, in order to understand this modified nanoprecipitation method at an atomistic level, large-scale all-atom molecular dynamics simulations were performed in parallel using the GROMOS53a6 forcefield parameters. It was observed that the mPEG-b-PCL chains act synergistically with the acetone molecules to dissolve the indomethacin nanoparticle while after the removal of the acetone molecules (mimicking the evaporation of the organic solvent) a polymer-drug nanoparticle was formed (yield 99%). This work could facilitate the development of more efficient methodologies for producing nanoparticles of hydrophobic drugs coated with amphiphilic polymers. The atomistic insight from the MD simulations in tandem with the data from the drug encapsulation experiments thus leads the way to a nanoformulation-by-design approach for therapeutic nanoparticles.

Molecular Mechanism for Azeotrope Formation in Ethanol/Benzene Binary Mixtures through Gibbs Ensemble Monte Carlo Simulation

Azeotropes have been studied for decades due to the challenges they impose on separation processes but fundamental understanding at the molecular level remains limited. Although molecular simulation has demonstrated its capability of predicting mixture vapor-liquid equilibrium (VLE) behaviors, including azeotropes, its potential for mechanistic investigation has not been fully exploited. In this study, we use the united atom transferable potentials for phase equilibria (TraPPE-UA) force field to model the ethanol/benzene mixture, which displays a positive azeotrope. Gibbs ensemble Monte Carlo (GEMC) simulation is performed to predict the VLE phase diagram, including an azeotrope point. The results accurately agree with experimental measurements. We argue that the molecular mechanism of azeotrope formation cannot be fully understood by studying the mixture liquid-state stability at the azeotrope point alone. Rather, azeotrope occurrence is only a reflection of the changing relative volatility between the two components over a much wider composition range. A thermodynamic criterion is thus proposed on the basis of the comparison of partial excess Gibbs energy between the components. In the ethanol/benzene system, molecular energetics shows that with increasing ethanol mole fraction, its volatility initially decreases but later plateaus, while benzene volatility is initially nearly constant and only starts to decrease when its mole fraction is low. Analysis of the mixture liquid structure, including a detailed investigation of ethanol hydrogen-bonding configurations at different composition levels, reveals the underlying molecular mechanism for the changing volatilities responsible for the azeotrope.

Expanding the Solubility Parameter Method MOSCED to Pyridinium-, Quinolinium-, Pyrrolidinium-, Piperidinium-, Bicyclic-, Morpholinium-, Ammonium-, Phosphonium-, and Sulfonium-Based Ionic Liquids

MOSCED (modified separation of cohesive energy density) is a solubility parameter method that offers an improved treatment of association interactions. Solubility parameter methods are well known for their ability to both make quantitative predictions and offer a qualitative description of the underlying molecular-level driving forces, lending themselves to intuitive solvent selection and design. Currently, MOSCED parameters are available for 130 organic solvents, water, and 33 imidazolium-based room temperature ionic liquids (ILs). In this work, we expand MOSCED to cover 66 additional ILs containing the pyridinium, quinolinium, pyrrolidinium, piperidinium, bicyclic, morpholinium, ammonium, phosphonium, and sulfonium cations using 10,052 experimental limiting activity coefficients. The resulting parameters may readily be used to predict the phase behavior in mixtures involving ILs.

Studying polymer solutions with particle-based models linked to classical density functionals: co-non-solvency

We demonstrate the potential of hybrid particle-based models, where interactions are introduced through functionals of local order parameters, in describing multicomponent polymer solutions. The link to a free-energy-like functional is advantageous for controlling the thermodynamics of the model. We focus on co-non-solvency - the collapse of polymer chains in dilute mixtures with two miscible good solvents, having different affinities towards the polymer. We employ a simple model where polymers and solvents are represented, respectively, by worm-like chains and single particles. Non-bonded interactions are captured by a polynomial which is third order in local densities and can, therefore, describe liquid-vapour coexistence. The parameterisation of the functional benefits from an elementary mean-field approximation to the statistical mechanics of the model. The model provides a framework for Monte Carlo simulations using a particle-to-mesh algorithm. Studies with conventional generic bead-spring and all-atom models have demonstrated that co-non-solvency is caused by preferential binding of the better solvent (termed cosolvent) with polymer. Hence, segmental loops bridged by cosolvent molecules are formed, initiating polymer collapse. The mesoscopic hybrid model differs conceptually from the conventional microscopic descriptions. Yet, it reproduces the same co-non-solvency mechanism supporting its universality. Films of adsorbed ternary solutions, showing co-non-solvency in the dilute regime, are considered at high concentrations. In this case, chains do not collapse. The properties of loops and tails of the adsorbed polymer agree with early theoretical predictions obtained for concentrated binary solutions.

Microstructures of negative and positive azeotropes

Azeotropes famously impose fundamental restrictions on distillation processes, yet their special thermodynamic properties make them highly desirable for a diverse range of industrial and technological applications. Using neutron diffraction, this study provides first insights into the microstructures of azeotropes.

Molecular Solution Behaviour of an Intermediate Biofuel Feedstock: Acetone-Butanol-Ethanol (ABE)

Mixtures of acetone, butanol, and ethanol (ABE) are common intermediate products in the production of biofuels via biomass fermentation. Their separation to yield, for example, bio-butanol, is still difficult due to the lack of a fundamental understanding of these mixtures at the molecular level. In order to bridge this gap, a detailed analysis of characteristic features of the vibrational spectrum is carried out. A systematic study of the binary solutions of acetone with ethanol and butanol does not only reveal a universal behaviour at the molecular level when acetone is mixed with short-chain alcohols, it also shows that the phenomena at a length scale between the molecules and in the macroscopic solution need to be taken into account to understand the structure-property relationships. The size of self-associated molecule clusters seems to determine whether or not a system exhibits an azeotrope. When a second alcohol is added to an acetone/alcohol solution, no additional non-idealities are induced, which is advantageous for modelling ternary ABE mixtures and for improving their processing in the production of biofuels.

Properties of the liquid–vapor interface of acetone–methanol mixtures, as seen from computer simulation and ITIM surface analysis

The intrinsic surface of acetone–methanol mixtures is studied by computer simulation and ITIM analysis.

Direct observation of coherent oscillations in solution due to microheterogeneous environment

We report, for the first time, direct observation of coherent oscillations in the ground-state of IR775 dye due to microheterogeneous environment. Using ultrafast near-infrared degenerate pump-probe technique centered at 800 nm, we present the dynamics of IR775 in a binary mixture of methanol and chloroform at ultra-short time resolution of 30 fs. The dynamics of the dye in binary mixtures, in a time-scale of a few fs to ~740 ps, strongly varies as a function of solvent composition (volume fraction). Multi-oscillation behavior of the coherent vibration was observed, which increased with decreasing percentage of methanol in the dye mixture. Maximum number of damped oscillations were observed in 20% methanol. The observed vibrational wavepacket motion in the ground-state is periodic in nature. We needed two cosine functions to fit the coherent oscillation data as two different solvents were used. Dynamics of the dye molecule in binary mixtures can be explained by wavepacket motion in the ground potential energy surface. More is the confinement of the dye molecule in binary mixtures, more is the number of damped oscillations. The vibrational cooling time, τ₂, increases with increase in the confinement of the system. The observed wavepacket oscillations in ground-state dynamics continued until 1.6 ps.

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

The Helmholtz free energy, energy, and entropy of mixing of eight different models of dimethyl sulfoxide (DMSO) with four widely used water models are calculated at 298 K over the entire composition range by means of thermodynamic integration along a suitably chosen thermodynamic path, and compared with experimental data. All 32 model combinations considered are able to reproduce the experimental values rather well, within RT (free energy and energy) and R (entropy) at any composition, and quite often the deviation from the experimental data is even smaller, being in the order of the uncertainty of the calculated free energy or energy, and entropy values of 0.1 kJ/mol and 0.1 J/(mol K), respectively. On the other hand, none of the model combinations considered can accurately reproduce all three experimental functions simultaneously. Furthermore, the fact that the entropy of mixing changes sign with increasing DMSO mole fraction is only reproduced by a handful of model pairs. Model combinations that (i) give the best reproduction of the experimental free energy, while still reasonably well reproducing the experimental energy and entropy of mixing, and (ii) that give the best reproduction of the experimental energy and entropy, while still reasonably well reproducing the experimental free energy of mixing, are identified.

Unveiling the mechanism of selective gate-driven diffusion of CO2 over N2 in MFU-4 metal–organic framework

Schematic view of the pore system of MFU-4. Two small pores, formed by chlorine ligands in cubic shape, are highlighted. Left: CO₂ molecule in a small pore. Right: N₂ molecule in a small pore.

SINGLE-PARTICLE AND PAIR DYNAMICAL PROPERTIES OF ACETONE–METHANOL MIXTURES CONTAINING CHARGED AND NEUTRAL SOLUTES: A MOLECULAR DYNAMICS STUDY

The dynamical properties of acetone–methanol mixtures containing either an ionic or a neutral hydrophobic solute are investigated by means of a series of molecular dynamics simulations. The primary goal has been to study how the solute and solvent dynamical properties change with variation of composition of the mixture ranging from pure acetone to pure methanol. The variations of structure and energetics of the mixture with composition are also calculated. The diffusion coefficients of both ionic and neutral solutes are found to show nonlinear variation with composition of the mixture, although the extent of nonlinearity in the diffusion of the neutral solute is much weaker. Calculations of appropriate solute-solvent distribution functions reveal the extent and nature of selective solvation of these solute species which play a role in determining the nonideal dynamical characteristics of these solutes. The free energies of solvation of the ionic solutes are also calculated and the results are discussed in the context of their dynamical behavior. The hydrogen bond statistics and dynamics of these mixtures are also calculated over their entire composition range. The energies and lifetimes of hydrogen bonds between an acetone and a methanol molecule or between two methanol molecules are found to increase with increase of acetone mole fraction of the mixture. Residence times of methanol molecules in solvation shells of acetone and methanol are also found to follow the same trend as relaxation times. However, these pair dynamical properties show essentially linear dependence on composition, thus behave almost ideally with respect to changes in composition of the mixture.

Computer simulation of volatile organic compound adsorption in atomistic models of molecularly imprinted polymers

Molecularly imprinted polymers (MIPs) offer a unique opportunity to significantly advance volatile organic compound (VOC) sensing technologies and a number of other applications. However, the development of these applications using MIPs has been hindered by poor understanding of the microstructure of MIPs, geometry of binding sites, and the details of molecular recognition processes in these materials. This is further complicated by the vast number of optimization parameters such as building components and processing conditions. Computer simulations and molecular modeling can help us understand adsorption and binding phenomena in MIPs on the molecular level and thus provide a route to more efficient MIP design strategies. So far, molecular models have been either oversimplified or severely limited in length scale, essentially focusing on a single binding site. Here, we propose a more general, atomistically detailed model that describes the microstructure of MIPs. We apply this model to investigate adsorption of pyridine, benzene, and toluene in MIPs and demonstrate that it is able to capture a number of essential experimental features. Therefore, this model can serve as a starting point in computational design and optimization of MIPs.