Atomistic molecular dynamics simulation and COSMO-SAC approach for enhanced 1,3-propanediol extraction with imidazolium-based ionic liquids

CONTEXT

1,3-Propanediol (1,3-PDO) is a key chemical in various industries, including pharmaceuticals and material sciences, and is projected to see significant market growth. However, the current challenges in its downstream processing, particularly in terms of cost and efficiency, highlight the need for innovative solutions. Our study delves into using ionic liquids (ILs) as a potential alternative, aiming to address these critical separation challenges more sustainably and efficiently. In this study, we utilized molecular dynamics (MD) simulations and the COSMO-SAC to examine 1,3-propanediol (1,3-PDO) extraction using four imidazolium-based ionic liquids with 1-butyl-3-methylimidazolium [Bmim] cation and with different anions bis(pentafluoroethanesulfonyl)imide [NPF2]-, bis(trifluoromethylsulfonyl)imide [NTF2]-, thiocyanate [SCN]-, and trifluoromethanesulfonate [TFO]-. Molecular dynamics simulations, incorporating analysis of radial distribution functions (RDF) and spatial distribution functions (SDF), revealed that [Bmim][SCN] and [Bmim][TFO] exhibit enhanced interactions with 1,3-PDO. Notably, [Bmim][SCN] formed the most hydrogen bonds, averaging 1.639 per molecule, due to its coordinating [SCN]- anion. This was in contrast to the fewer hydrogen bonds formed by non-coordinating anions in [Bmim][NPF2] and [Bmim][NTF2]. In ternary systems, [Bmim][SCN] and [Bmim][TFO] demonstrated superior selectivity for 1,3-PDO extraction compared to the other ionic liquids, with selectivity values around 29. These findings, supported by COSMO-SAC predictive modeling, highlight the potential of [Bmim][SCN] as a promising candidate for 1,3-PDO extraction, emphasizing the importance of anion selection in optimizing ionic liquid properties for this application.

METHODS

In our study, we employed MD simulations, incorporating the OPLS-AA force field, and COSMO-SAC to investigate the extraction of 1,3-PDO using imidazolium-based ionic liquids: [Bmim][NTF2], [Bmim][NPF2], [Bmim][SCN], and [Bmim][TFO]. The MD simulations were conducted using LAMMPS software, focusing on elucidating the RDF, SDF, and hydrogen bonding. Analysis of the distribution coefficient (β) and selectivity (S) for the ternary mixture was also conducted. These aspects of the simulation were analyzed using TRAVIS and VMD software. Additionally, the COSMO-SAC model was employed to determine the activity coefficients of 1,3-PDO in the ionic liquids, with molecular optimization conducted using Gaussian16 and sigma profile calculations performed using COSMO-SAC.

Theoretical study of the saturation and nature of the hydrogen bonds to gold

Traditional hydrogen bonds are well-known to exhibit directionality and saturation. By contrast, gold involved hydrogen bonds (GHBs) have been extensively studied but remain lack of in-depth understanding towards the intrinsic nature and saturation property. This work exemplifies three series of complexes: [L-Au-L]-⋯(HF)n (L = H, CH3, (CH3)3; n = 1-8) containing GHBs to dig into the intrinsic nature with the aid of multiple theoretical analysis methods, finding that the formation of GHB is highly subject to orbital interactions along with steric hindrance. Moreover, the saturation level of GHBs largely depends on the ligand attached to the gold center, since different ligands typically possess varying electron-giving ability and steric volume. This work confirms the coexistence of as many as 6 GHBs for one Au atom and thoroughly studies the saturation level of GHBs, which will provide new insights into GHBs and facilitate future synthesis of more complicated gold complexes.

Mechanism, Kinetics, and Potential of Mean Force of Evaporation of Water from Aqueous Sodium Chloride Solutions of Varying Concentrations

The mechanism, kinetics, and potential of mean force of evaporation of water from aqueous NaCl solutions are investigated through both unbiased molecular dynamics simulations and also biased simulations using the umbrella sampling method. The results are obtained for aqueous solutions of three different NaCl concentrations ranging from 0.6 to 6.0 m and also for pure water. The rate of evaporation is found to decrease in the presence of ions. It is found that the process of evaporation of a surface water molecule from ionic solutions can be triggered through its collision with another water or chloride ion. Such collisions provide the additional kinetic energy that is required for evaporation. However, when the collision takes place with a Cl^- ion, the evaporation of the escaping water also involves a collision with water in the vicinity of the ion at the same time along with the ion-water collision. These two collisions together provide the required kinetic energy for escape of the evaporating water molecule. Thus, the mechanism of evaporation process of ionic solutions can be more complex than that of pure water. The potential of mean force (PMF) of evaporation is found to be positive and it increases with increasing ion concentration. Also, no barrier in the PMF is found to be present for the condensation of water from vapor phase to the surfaces of the solutions. A detailed analysis of the unsuccessful evaporation attempts by surface water molecules is also made in the current study.

Spontaneous Formation of Hydrogen Peroxide in Water Microdroplets

There is accumulating evidence that many chemical reactions are accelerated by several orders of magnitude in micrometer-sized aqueous or organic liquid droplets compared to their corresponding bulk liquid phase. However, the molecular origin of the enhanced rates remains unclear as in the case of spontaneous appearance of 1 μM hydrogen peroxide in water microdroplets. In this Letter, we consider the range of ionization energies and whether interfacial electric fields of a microdroplet can feasibly overcome the high energy step from hydroxide ions (OH^-) to hydroxyl radicals (OH^•) in a primary H₂O₂ mechanism. We find that the vertical ionization energies (VIEs) of partially solvated OH^- ions are greatly lowered relative to the average VIE in the bulk liquid, unlike the case of the Cl^- anion which shows no reduction in the VIEs regardless of solvation environment. Overall reduced hydrogen-bonding and undercoordination of OH^- are structural features that are more readily present at the air-water interface, where the energy scale for ionization can be matched by statistically probable electric field values.

D-Mannosamine hydrochloride (2-amino-2-deoxy-D-mannose hydrochloride): ionic hydrogen bonding in saccharides involving chloride and aminium ions

D-Mannosamine hydrochloride (2-amino-2-deoxy-D-mannose hydrochloride), C6H14NO5 +·Cl−, (I), crystallized from a methanol/ethyl acetate/n-hexane solvent mixture at room temperature in a 4 C 1 chair conformation that is slightly distorted towards the C3,O5 B form. A comparison of the structural parameters of (I) with the corresponding parameters in α-D-glucosamine hydrochloride, (II), and β-D-galactosamine hydrochloride, (III)/(III′), was undertaken to evaluate the effects of ionic hydrogen bonding on structural properties. Three types of ionic hydrogen bonds are present in the crystals of (I)–(III)/(III′), i.e. N+—H...O, N+—H...Cl−, and O—H...Cl−. The exocyclic structural parameters in (I), (II), and (III)/(III′) appear to be most influenced by this bonding, especially the exocyclic hydroxy groups, which adopt eclipsed conformations enabled by ionic hydrogen bonding to the chloride anion. Anomeric disorder was observed in crystals of (I), with an α:β ratio of 37:63. However, anomeric configuration appears to exert minimal structural effects; that is, bond lengths, bond angles, and torsion angles are essentially identical in both anomers. The observed disorder at the anomeric C atom of (I) appears to be caused by the presence of the chloride anion and atom O3 or O4 in proximal voids, which provide opportunities for hydrogen bonding to atom O1 in both axial and equatorial orientations.

Layer-by-Layer Cell Encapsulation for Drug Delivery: The History, Technique Basis, and Applications

The encapsulation of cells with various polyelectrolytes through layer-by-layer (LbL) has become a popular strategy in cellular function engineering. The technique sprang up in 1990s and obtained tremendous advances in multi-functionalized encapsulation of cells in recent years. This review comprehensively summarized the basis and applications in drug delivery by means of LbL cell encapsulation. To begin with, the concept and brief history of LbL and LbL cell encapsulation were introduced. Next, diverse types of materials, including naturally extracted and chemically synthesized, were exhibited, followed by a complicated basis of LbL assembly, such as interactions within multilayers, charge distribution, and films morphology. Furthermore, the review focused on the protective effects against adverse factors, and bioactive payloads incorporation could be realized via LbL cell encapsulation. Additionally, the payload delivery from cell encapsulation system could be adjusted by environment, redox, biological processes, and functional linkers to release payloads in controlled manners. In short, drug delivery via LbL cell encapsulation, which takes advantage of both cell grafts and drug activities, will be of great importance in basic research of cell science and biotherapy for various diseases.

Insights on choline chloride-based deep eutectic solvent (reline) + primary alcohol mixtures: a molecular dynamics simulation study

Deep eutectic solvents (DESs) emerged as green solvents for new generation technologies owing to their high chemical and thermal stability. Addition of restricted amount of organic solvents into the DESs plays a significant role in the improvement of thermodynamic and the transport properties to work as a potential solvent in process industries. In this paper, molecular dynamics (MD) simulations were performed to understand the thermophysical and transport properties of choline chloride-based DES (reline) and primary alcohol (methanol and ethanol) mixture in relation to microscopic structure. Density, radial distribution function, coordination number, average number of H-bond, diffusion coefficient and spatial distribution function was calculated in order to understand the structure and involvement of H-bond network at an atomic level. H-bond and spatial distribution function analyses revealed that the chloride ion prefers to be spatially distributed around hydroxyl group of alcohol and found to be more pronounced upon increase in alcohol concentration. As a consequence, it was observed that the H-bonds between Cl^- and urea decreases overall with the loading of alcohol and effect is more pronounced beyond a concentration of 0.4. Self-diffusion values for choline, Cl^- and urea were found to be increased significantly upon increase in concentration of alcohol from 0.6 to 0.8. Overall, our simulation points to the interplay and interactions between the chloride ions and the solvents in determining the structural and transport properties of choline chloride-based DES.

In Silico Elucidation of Molecular Picture of Water-Choline Chloride Mixture

Choline chloride (ChCl) is a component of several deep eutectic solvents (DESs) having numerous applications. Recent studies have reported manifold promising use of aqueous choline chloride solution as an alternative to DES, where water plays the role of the hydrogen-bond donor. The characteristic physical properties of the DESs and aqueous DES originate from the "inter-" and intraspecies hydrogen-bond network formed by the constituents. However, a detailed molecular-level picture of choline chloride and water mixture is largely lacking in the literature. This motivates us to carry out extensive all-atom molecular dynamics simulations of the ChCl-water mixture of varying compositions. Our analyses clearly show an overall increase in the interspecies association with an increase in ChCl concentration. At higher concentrations, the trimethylammonium groups of choline are stabilized by a nonpolar interaction, whereas the hydroxyl groups preferentially interact with water. Chloride ions are found to be involved in two types of interactions: one where chloride ions intercalate two or more choline cations, and the other one where they are surrounded by five to six water molecules forming solvated chloride ions. However, the relative fractions of these two types of associations depend on the concentration of ChCl in the mixture. Another important structural aspect is the disruption of the hydrogen-bonded water network due to the presence of both choline cations and chloride ions. However, chloride ions participate to partially restore the tetrahedral arrangement of partners around water molecules.

Stimuli-Responsive, Hydrolyzable Poly(Vinyl Laurate-co-vinyl Acetate) Nanoparticle Platform for In Situ Release of Surfactants

A stimuli-responsive, sub-100 nm nanoparticle (NP) platform with a hydrolyzable ester side chain for in situ generation of surfactants is demonstrated. The NPs were synthesized via copolymerization of vinyl-laurate and vinyl-acetate [p-(VL-co-VA), 3:1 molar ratio] and stabilized with a protective poly(ethylene-glycol) shell. The NPs are ∼55 nm in diameter with a zeta potential of -54 mV. Hydrolysis kinetics in an accelerated, base-catalyzed reaction show release of about 11 and 30% of the available surfactant at 25 and 80 °C, respectively. The corresponding values in seawater are 22 and 76%. The efficiency of the released surfactant in reducing the interfacial tension, altering wettability, and stabilizing oil-water emulsion was investigated through contact angle measurements and laser confocal scanning microscopy and benchmarked to sodium laurate, a commercially available surfactant. All these measurements demonstrate both the efficacy of the NP system for surfactant delivery and the ability of the released surfactant to alter wettability and stabilize an oil-water emulsion.

Revealing the Frank-Evans "Iceberg" Structures within the Solvation Layer around Hydrophobic Solutes

Using computer simulations, the structural properties of solvation water of three model hydrophobic molecules, methane and two fullerenes (C60 and C80), were studied. Systems were simulated at temperatures in the range of 250-298 K. By analyzing both the local ordering of the molecules of water in the solvation layers and the structure of hydrogen bond network, it is shown that in the solvation layer of hydrophobic molecules, ordered aggregates consisting of water molecules are formed. Even though it is difficult to define the exact structure of these aggregates, their existence alone is clearly noticeable. Moreover, these aggregates become more pronounced with the decrease of temperature. The existence of the ordered aggregates around the hydrophobic solutes complies with the concept of "icebergs" proposed by Frank and Evans.

Dielectric relaxation of deep eutectic solvent + water mixtures: structural implications and application to microwave heating

All dipolar species at their full individual strengths but synchronized in motion: structural implications of cooperative dynamics in glyceline/water and reline/water mixtures.