Modeling solution vapor equilibria with solvation and solute assembly

Modeling Henry's law and phase separations of water-NaCl-organic mixtures with solvation and ion-pairing

Empirical measurements of solution vapor pressure of ternary acetonitrile (MeCN) H2O-NaCl-MeCN mixtures were recorded, with NaCl concentrations ranging from zero to the saturation limit, and MeCN concentrations ranging from zero to an absolute mole fraction of 0.64. After accounting for speciation, the variability of the Henry's law coefficient at vapor-liquid equilibrium (VLE) of MeCN ternary mixtures decreased from 107% to 5.1%. Solute speciation was modeled using a mass action solution model that incorporates solute solvation and ion-pairing phenomena. Two empirically determined equilibrium constants corresponding to solute dissociation and ion pairing were utilized for each solute. When speciation effects were considered, the solid-liquid equilibrium of H2O-NaCl-MeCN mixtures appear to be governed by a simple saturation equilibrium constant that is consistent with the binary H2O-NaCl saturation coefficient. Further, our results indicate that the precipitation of NaCl in the MeCN ternary mixtures was not governed by changes in the dielectric constant. Our model indicates that the compositions of the salt-induced liquid-liquid equilibrium (LLE) boundary of the H2O-NaCl-MeCN mixture correspond to the binary plateau activity of MeCN, a range of concentrations over which the activity remains largely invariant in the binary water-MeCN system. Broader comparisons with other ternary miscible organic solvent (MOS) mixtures suggest that salt-induced liquid-liquid equilibrium exists if: (1) the solution displays a positive deviation from the ideal limits governed by Raoult's law; and (2) the minimum of the mixing free energy profile for the binary water-MOS system is organic-rich. This work is one of the first applications of speciation-based solution models to a ternary system, and the first that includes an organic solute.

Solvent Organization and Electrostatics Tuned by Solute Electronic Structure: Amide versus Non-Amide Carbonyls

The ability to exploit carbonyl groups to measure electric fields in enzymes and other complex reactive environments by using the vibrational Stark effect has inspired growing interest in how these fields can be measured, tuned, and ultimately designed. Previous studies have concentrated on the role of the solvent in tuning the fields exerted on the solute. Here, we explore instead the role of the solute electronic structure in modifying the local solvent organization and electric field exerted on the solute. By measuring the infrared absorption spectra of amide-containing molecules, as prototypical peptides, and contrasting them with non-amide carbonyls in a wide range of solvents, we show that these solutes experience notable differences in their frequency shifts in polar solvents. Using vibrational Stark spectroscopy and molecular dynamics simulations, we demonstrate that while some of these differences can be rationalized by using the distinct intrinsic Stark tuning rates of the solutes, the larger frequency shifts for amides and dimethylurea primarily result from the larger solvent electric fields experienced by their carbonyl groups. These larger fields arise due to their stronger p-π conjugation, which results in larger C═O bond dipole moments that further induce substantial solvent organization. Using electronic structure calculations, we decompose the electric fields into contributions from solvent molecules that are in the first solvation shell and those from the bulk and show that both of these contributions are significant and become larger with enhanced conjugation in solutes. These results show that structural modifications of a solute can be used to tune both the solvent organization and electrostatic environment, indicating the importance of a solute-centric paradigm in modulating and designing the electrostatic environment in condensed-phase chemical processes.

Solvent-driven fractional crystallization for atom-efficient separation of metal salts from permanent magnet leachates

This work reports a dimethyl ether-driven fractional crystallization process for separating rare earth elements and transition metals. The process has been successfully applied in the treatment of rare earth element-bearing permanent magnet leachates as an atom-efficient, reagent-free separation method. Using ~5 bar pressure, the solvent was dissolved into the aqueous system to displace the contained metal salts as solid precipitates. Treatments at distinct temperatures ranging from 20-31 °C enable crystallization of either lanthanide-rich or transition metal-rich products, with single-stage solute recovery of up to 95.9% and a separation factor as high as 704. Separation factors increase with solution purity, suggesting feasibility for eco-friendly solution treatments in series and parallel to purify aqueous material streams. Staged treatments are demonstrated as capable of further improving the separation factor and purity of crystallized products. Upon completion of a crystallization, the solvent can be recovered with high efficiency at ambient pressure. This separation process involves low energy and reagent requirements and does not contribute to waste generation.

Mass action model of solution activity via speciation by solvation and ion pairing equilibria

Solutes and their concentrations influence many natural and anthropogenic solution processes. Electrolyte and solution models are used to quantify and predict such behavior. Here we present a mechanistic solution model based on mass action equilibria. Solvation and ion pairing are used to model speciated solute and solvent concentrations such that they correlate to a solution's vapor pressure (solvent activity) according to Raoult's law from dilute conditions to saturation. This model introduces a hydration equilibrium constant (K_ha) that is used with either an ion dissociation constant (K_id) or a hydration modifier (m) with an experimentally determined ion dissociation constant, as adjustable parameters to fit vapor-liquid equilibrium data. The modeled solvation equilibria are accompanied by molecular dynamics (MD) studies that support a decline in the observed degree of solvation with increased concentration. MD calculations indicate this finding is a combination of a solvent that solvates multiple solutes, and changes in a solute's solvation sphere, with the dominant factor changing with concentration. This speciation-based solution model is lateral to established electrostatics-based electrolyte theories. With its basis in mass action, the model can directly relate experimental data to the modeled solute and solvent speciated concentrations and structures.