Local Density Augmentation in Supercritical Solvents: Electronic Shifts of Anthracene Derivatives

Temperature-Dependent Dielectric Relaxation Measurements of (Betaine + Urea + Water) Deep Eutectic Solvent in Hz-GHz Frequency Window: Microscopic Insights into Constituent Contributions and Relaxation Mechanisms

A combined experimental and simulation study of dielectric relaxation (DR) of a deep eutectic solvent (DES) composed of betaine, urea, and water with the composition [Betaine:Urea:Water = 11.7:12:1 (weight ratio) and 9:18:5 (molar ratio)] was performed to explore and understand the interaction and dynamics of this system. Temperature-dependent (303 ≤ T/K ≤ 343) measurements were performed over 9 decades of frequency, combining three different measurement setups. Measured DR, comprising four distinct steps with relaxation times spreading over a few picoseconds to several nanoseconds, was found to agree well with simulations. The simulated total DR spectra, upon dissection into three self (intraspecies) and three cross (interspecies) interaction contributions, revealed that the betaine-betaine self-term dominated (∼65%) the relaxation, while the urea-urea and water-water interactions contributed only ∼7% and ∼1%, respectively. The cross-terms (betaine-urea, betaine-water, and urea-water) together accounted for <30% of the total DR. The slowest DR component with a time constant of ∼1-10 ns derived dominant contribution from betaine-betaine interactions, where betaine-water and urea-water interactions also contributed. The subnanosecond (0.1-0.6 ns) time scale originated from all interactions except betaine-water interaction. An extensive interaction of water with betaine and urea severely reduced the average number of water-water H-bonds (∼0.7) and heavily decreased the static dielectric constant of water in this DES (εs ∼ 2). Furthermore, simulated first rank collective single particle reorientational relaxations (C1(t)) and the structural H-bond fluctuation dynamics (CHB (t)) exhibited multiexponential kinetics with time scales that corresponded well with those found both in the simulated and measured DR.

Interactions and Dynamics in Aqueous Solutions of pH-Responsive Polymers: A Combined Fluorescence and Dielectric Relaxation Study

Interaction and dynamics of aqueous solutions of pH-responsive smart polymers are investigated via steady-state, time-resolved fluorescence emission spectroscopy with the help of external local reporter coumarin 153 (C153), while MHz to GHz dielectric relaxation spectroscopic (DRS) measurement reports the intrinsic medium relaxation features. A series of pH-responsive random copolymers (DPL-DP60) comprising of a pH-responsive moiety 2-((leucinyl)oxy)ethyl methacrylate (l-Leu-HEMA) and hydrophobic methyl methacrylate (MMA) are synthesized and characterized. A balance between the pH-responsive (l-Leu-HEMA) and the hydrophobic (MMA) content dictates the phase transition pH, which is found to be ∼5-7 for these aqueous copolymer solutions (1 mg/mL). Dynamic light scattering measurements in aqueous solutions of these polymers reflect a small particle size (∼2-8 nm) at solution pH below their individual phase transition pH, while a large particle size (∼140-340 nm) forms beyond their phase transition pH. No signature of a phase transition pH-driven abrupt change in static and dynamic properties of aqueous polymer solutions has been registered from pH-dependent dielectric relaxation as well as solute (C153)-centric fluorescence measurements. A significant impact of varying the l-Leu-HEMA/MMA segment ratio on steady-state fluorescence emission and rotational anisotropy decay of the fluorophore solute (C153) has been observed. MHz to GHz DRS in aqueous solutions of these pH-responsive polymers reflects bulk water-like dielectric features.

Water in biodegradable glucose–water–urea deep eutectic solvent: modifications of structure and dynamics in a crowded environment

Molecular dynamics simulations have been performed on a highly viscous (η ∼ 255 cP) naturally abundant deep eutectic solvent (NADES) composed of glucose, urea and water in a weight ratio of 6 : 4 : 1 at 328 K.

Cloud Point Driven Dynamics in Aqueous Solutions of Thermoresponsive Copolymers: Are They Akin to Criticality Driven Solution Dynamics?

Cloud point driven interaction and relaxation dynamics of aqueous solutions of amphiphilic thermoresponsive copolymers were explored through picosecond resolved and steady state fluorescence measurements employing hydrophilic (coumarin 343, C343) and hydrophobic (coumarin 153, C153) solute probes of comparable sizes. These thermoresponsive random copolymers, with tunable cloud point temperatures (T_cp's) between 298 and 323 K, were rationally designed first and then synthesized via reversible addition-fragmentation chain transfer (RAFT) copolymerization of methyl methacrylate (MMA) and poly(ethylene glycol) monomethyl ether methacrylate (PEGMA). Subsequently, copolymers were characterized by NMR spectroscopy and size exclusion chromatography (SEC). A balance between the hydrophilic (PEGMA) and the hydrophobic (MMA) content dictates the critical aggregation concentration (CAC), with CAC ∼ 2-14 mg/L for these copolymers in aqueous media. No abrupt changes in the steady state spectral features of both C153 and C343 in the aqueous solutions of these polymers near but below the cloud point temperatures were observed. Interestingly, spectral properties of C153 in these solutions show the impact of hydrophobic/hydrophilic interaction balance but not by those of C343. More specifically, C153 reported a blue shift (relative to that in neat water) and heterogeneity in its local environment. This suggested different locations for the hydrophilic (C343) and the hydrophobic (C153) probes. In addition, the excited state fluorescence lifetime (⟨τ_life⟩) of C153 increased with the increase of hydrophobic (MMA) content in these copolymers. However, C343 reported no such variations, although fluorescence anisotropy decays for both solutes were significantly slowed down in these aqueous solutions compared to neat water. Anisotropy decays indicated bimodal time-dependent friction for these solutes in aqueous solutions of these copolymers but monomodal in neat water. A linear dependence of the average rotational relaxation rates (⟨k_rot⟩ = ⟨τ_rot⟩^-1) of the type ⟨k_rot⟩ ∝ (|T - T_cp|/T_cp)^γ with negative values for the exponent γ was observed for both solutes. No slowing down of the solute rotation with temperature approaching the T_cp was detected; rather, rotation became faster upon increasing the solution temperature, suggesting domination of the local friction.

Interaction and Dynamics in a Fully Biodegradable Glucose-Containing Naturally Abundant Deep Eutectic Solvent: Temperature-Dependent Time-Resolved Fluorescence Measurements

A new room-temperature deep eutectic solvent (DES) composed of glucose, urea, and water has been prepared and its relaxation dynamics explored via temperature-dependent time-resolved fluorescence measurements employing hydrophilic and hydrophobic solute probes. Differential scanning calorimetry measurements indicate a glass transition temperature (T_g) of ∼236 K. Measured viscosity coefficients (η) vary from ∼600 to ∼100 cP in the temperature range 318 ≤ T/K ≤ 343 and exhibit Arrhenius-type temperature dependence with an activation energy of ∼65 kJ mol^-1. Interestingly, this DES forms a stable liquid at ∼300 K but is too viscous to be accurately measured by us below 318 K. Temperature-dependent dynamic fluorescence anisotropy measurements using hydrophobic and hydrophilic solutes of similar sizes reveal bi-exponential kinetics and Arrhenius-type temperature dependence for solute rotation times (⟨τ_r⟩) but with significantly decreased activation energies, ∼31 kJ mol^-1 (hydrophobic) and ∼21 kJ mol^-1 (hydrophilic). Deviation from hydrodynamics is further reflected in the strong fractional viscosity dependence of ⟨τ_r⟩: ⟨τ_r⟩ ∝ (η/T)^p with p ≈ 0.3-0.5, indicating pronounced temporal heterogeneity in the relaxation dynamics. Dynamic fluorescence Stokes shift measurements (temporal resolution ∼85 ps) produce dynamic shifts of ∼500-700 cm^-1, bi-exponential solvation energy relaxation with time constants in the range ∼0.2 ns and ∼4 ns, and estimated missing amplitudes of ∼65-75%. Impact of the density difference between a nonpolar solvent and this DES on the estimated missing amplitudes is explored via measuring the temperature-dependent densities and refractive indices of this DES. Lifetime measurements suggest considerable temperature dependence for the hydrophobic solute but no such dependence for the hydrophilic one. Excitation energy dependence of fluorescence emission of various solutes with widely different lifetimes indicates mild spatial heterogeneity for this DES.

Widom Delta of Supercritical Gas-Liquid Coexistence

Density fluctuations and the Widom line are of great importance in understanding the critical phenomena and the behaviors of supercritical fluids (SCFs). We report on the direct classification of liquid-like and gas-like molecules coexisting in the SCF, identified by machine learning analysis on simulation data. The deltoid coexistence region encloses the Widom line and may therefore be termed the Widom delta. Number fractions of gas-like and liquid-like particles are found to undergo continuous transition across the delta, following a simplified two-state model. These fractions are closely related to the magnitude of supercritical anomaly, which originates from the fluctuation between the two types. This suggests a microscopic view of the SCF as a mixture of liquid-like and gas-like structures, providing an integrative explanation to the anomalous behaviors near the critical point and the Widom line.

Modeling Solvation in Supercritical CO2

In recent decades, a microscopic understanding of solute-solvent intermolecular interactions has been key to advances in technologies based on supercritical carbon dioxide. In many cases, computational work has provided the impetus for new discoveries, shedding new light on important concepts such as the local structure around the solute in the supercritical medium, the influence of the peculiar properties of the latter on the molecular behavior of dissolved substances and, importantly, CO₂ -philicity. In this Review, the theoretical work that has been relevant to these developments is surveyed and, by presenting some crucial open questions, the possible routes to achieving further progress based on the interplay between theory and experiments is discussed.

How Heterogeneous Are Trehalose/Glycerol Cryoprotectant Mixtures? A Combined Time-Resolved Fluorescence and Computer Simulation Investigation

Heterogeneity and molecular motions in representative cryoprotectant mixtures made of trehalose and glycerol are investigated in the temperature range 298 ≤ T (K) ≤ 353, via time-resolved fluorescence Stokes shift and anisotropy measurements, and molecular dynamics simulations of four-point density-time correlations and H-bond relaxations. Mixtures containing 5 and 20 wt % of trehalose along with neat glycerol are studied. Viscosity coefficients for these systems lie in the range 0.30 < η (P) < 23. Measured solute (Coumarin 153) rotation and solvation times reveal a substantial departure from the hydrodynamic viscosity dependence, suggesting the strong microheterogeneous nature of these systems. Fluorescence anisotropy decays are highly nonexponential, reflecting a non-Markovian character of the medium friction. A complete missing of the Stokes shift dynamics in these systems at 298 K but partial detection of it at other higher temperatures (shift magnitude being ∼400-600 cm^-1) indicates rigid solute environments. An amorphous solid-like feature emerges in the simulated radial distribution functions at these temperatures. Analyses of mean squared displacements reveal rattling-in-a-cage motion, non-Gaussian displacement distributions, and strong dynamic heterogeneity features. Simulated dynamic structure factors and four-point correlations hint, respectively, at very long α-relaxation and correlated time scales at 298 K. This explains the long solute rotation times (∼80-200 ns) measured at 298 K. Stretched exponential decay of the simulated H-bond relaxations with long time scales further highlights the strong temporal heterogeneity and slow dynamics inherent to these systems. In summary, this work provides the first insight into the molecular motions and interspecies interaction in a representative cryoprotectant mixture, and stimulates further study to investigate the interconnection between cryoprotection and dynamic heterogeneity.

Colloidal interactions in liquid CO2--a dry-cleaning perspective

Liquid CO(2) is a viable alternative for the toxic and environmentally harmful solvents traditionally used in dry-cleaning industry. Although liquid CO(2) dry-cleaning is being applied already at a commercial scale, it is still a relatively young technique which poses many challenges. The focus of this review is on the causes of the existing problems and directions to solve them. After presenting an overview of the state-of-the-art, we analyze the detergency challenges from the fundamentals of colloid and interface science. The properties of liquid CO(2) such as dielectric constant, density, Hamaker constant, refractive index, viscosity and surface tension are presented and in the subsequent chapters their effects on CO(2) dry-cleaning operation are delineated. We show, based on theory, that the van der Waals forces between a model soil (silica) and model fabric (cellulose) through liquid CO(2) are much stronger compared to those across water or the traditional dry-cleaning solvent PERC (perchloroethylene). Prevention of soil particle redeposition in liquid CO(2) by electrostatic stabilization is challenging and the possibility of using electrolytes having large anionic parts is discussed. Furthermore, the role of different additives used in dry-cleaning, such as water, alcohol and surfactants, is reviewed. Water is not only used as an aid to remove polar soils, but also enhances adhesion between fabric and soil by forming capillary bridges. Its role as a minor component in liquid CO(2) is complex as it depends on many factors, such as the chemical nature of fabrics and soil, and also on the state of water itself, whether present as molecular solution in liquid CO(2) or phase separated droplets. The phenomena of wicking and wetting in liquid CO(2) systems are predicted from the Washburn-Lucas equation for fabrics of various surface energies and pore sizes. It is shown that nearly complete wetting is desirable for good detergency. The effect of mechanical action and fluid dynamic conditions on dry-cleaning is analyzed theoretically. From this it follows that in liquid CO(2) an order of magnitude higher Reynold's number is required to exceed the binding forces between fabric and soil as opposed to PERC or water, mainly due to the strong van der Waals forces and the low viscosity of CO(2) at dry-cleaning operational conditions.

Morphology and stability of CO2-in-water foams with nonionic hydrocarbon surfactants

The morphologies, stabilities, and viscosities of high-pressure carbon dioxide-in-water (C/W) foams (emulsions) formed with branched nonionic hydrocarbon surfactants were investigated by in situ optical microscopy and capillary rheology. Over two dozen hydrocarbon surfactants were shown to stabilize C/W foams with Sauter mean bubble diameters as low as 1 to 2 microm. Coalescence of the C/W foam bubbles was rare for bubbles larger than about 0.5 microm over a 60 h time frame, and Ostwald ripening became very slow. By better blocking of the CO(2) and water phases with branched and double-tailed surfactants, the interfacial tension decreases, the surface pressure increases, and the C/W foams become very stable. For branched surfactants with propylene oxide middle groups, the stabilities were markedly lower for air/water foams and decane-water emulsions. The greater stability of the C/W foams to coalescence may be attributed to a smaller capillary pressure, lower drainage rates, and a sufficient surface pressure and thus limiting surface elasticity, plus small film sizes, to hinder spatial and surface density fluctuations that lead to coalescence. Unexpectedly, the foams were stable even when the surfactant favored the CO(2) phase over the water phase, in violation of Bancroft's rule. This unusual behavior is influenced by the low drainage rate, which makes Marangoni stabilization of less consequence and the strong tendency of emerging holes in the lamella to close as a result of surfactant tail flocculation in CO(2). The high distribution coefficient toward CO(2) versus water is of significant practical interest for mobility control in CO(2) sequestration and enhanced oil recovery by foam formation.

Molecular dynamics simulation of benzene-propene-cumene mixtures in different phases

Molecular dynamics simulations have been performed to study the microscopic configuration and dynamic behavior of mixtures of benzene, propene, and cumene for the cumene synthesis process. The comparisons have been made for the intermolecular radial distribution functions of the binary and ternary mixtures at the conditions that are near, below, and above their respective critical points. The results have shown that in both binary and ternary mixtures propene molecules have a small tendency to cluster in the liquid state, but at supercritical conditions they tend to be uniformly distributed. Contrary to propene, cumene molecules have a tendency to cluster in ternary mixtures. A moderate local density augmentation is also found in the benzene-propene binary supercritical fluid. The excess functions for benzene-propene binary mixtures have shown that there exists an enhancement of the potential energy when benzene mixes with propene. This enhancement provides a rational explanation for the experimental critical properties, which exhibit the behavior of the nonmonotonous dependence of critical pressure on compositions.

Intramolecular charge transfer reaction, polarity, and dielectric relaxation in AOT/water/heptane reverse micelles: pool size dependence

Intramolecular charge transfer (ICT) reaction in a newly synthesized molecule, of 4-(1-morpholenyl) benzonitrile (M6C), in AOT/water/heptane reverse micelles at different pool sizes has been studied by using steady-state and time-resolved fluorescence emission spectroscopy. The pool size dependences of the reaction equilibrium constant and reaction rate have been explained in terms of the average polarity of the confined solvent pools estimated from the fluorescence emission Stokes shift of a nonreactive probe, coumarin 153, dissolved in these microemulsions. The complex permittivity measurements in the frequency range 0.01<or=nu/GHz<or=2 for these microemulsions at different pool sizes (0<or=w0<or=40) and AOT concentrations (0.1<or=c/M<or=0.5) at 298.15 K have also been performed. At sufficient water content, a large dispersion with a relaxation time of approximately 600 ps has been observed at approximately 300 MHz and attributed to the average reorientation of water molecules residing in the close vicinity of the polar interface of the AOT headgroup and n-heptane. The reorientation of these interfacial water molecules is probably responsible for the nanosecond component observed in numerous polar solvation dynamics experiments in these reverse micelles. Subsequently, the estimated polarity and the measured reorientational time scale have been used to explain the dramatic slowing down of the ICT reaction rate and its dependence on pool size in these confined environments.

Molecular simulations of outersphere reorganization energies in polar and quadrupolar solvents. The case of intramolecular electron and hole transfer

Outersphere reorganization energies (lambda) for intramolecular electron and hole transfer are studied in anion- and cation-radical forms of complex organic substrates (p-phenylphenyl-spacer-naphthyl) in polar (water, 1,2-dichloroethane, tetrahydrofuran) and quadrupolar (supercritical CO2) solvents. Structure and charge distributions of solute molecules are obtained at the HF/6-31G(d,p) level. Standard Lennard-Jones parameters for solutes and the nonpolarizable simple site-based models of solvents are used in molecular dynamics (MD) simulations. Calculation of lambda is done by means of the original procedure, which treats electrostatic polarization of a solvent in terms of a usual nonpolarizable MD scheme supplemented by scaling of reorganization energies at the final stage. This approach provides a physically relevant background for separating inertial and inertialless polarization responses by means of a single parameter epsilon(infinity), optical dielectric permittivity of the solvent. Absolute lambda values for hole transfer in 1,2-dichloroethane agree with results of previous computations in terms of the different technique (MD/FRCM, Leontyev, I. V.; et al. Chem. Phys. 2005, 319, 4). Computed lambda values for electron transfer in tetrahydrofuran are larger than the experimental values by ca. 2.5 kcal/mol; for the case of hole transfer in 1,2-dichloroethane the discrepancy is of similar magnitude provided the experimental data are properly corrected. The MD approach gives nonzero lambda values for charge-transfer reaction in supercritical CO2, being able to provide a uniform treatment of nonequilibrium solvation phenomena in both quadrupolar and polar solvents.

Local Density Inhomogeneities and Dynamics in Supercritical Water: A Molecular Dynamics Simulation Approach

Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.