Solvatochromic behavior of phenol blue in CO2+ethanol and CO2+n-pentane mixtures in the critical region and local composition enhancement

The UV-Vis spectra of probe phenol blue in CO(2)+ethanol and CO(2)+n-pentane binary mixtures were studied at 308.15 K and different pressures. The experiments were conducted in both supercritical region and subcritical region of the mixtures by changing the compositions of the mixed solvents. On the basis of the experimental results the local compositions of the solvents about phenol blue were estimated by neglecting the size difference of CO(2) and the cosolvents. Then the local composition data were corrected by a method proposed in this work, which is mainly based on Lennard-Jones sphere model. It was demonstrated that the local mole fraction of the cosolvents is higher than that in the bulk solution at all the experimental conditions. In the near critical region of the mixed solvents the local composition enhancement, defined as the ratio of cosolvent mole fraction about the solute to that in the bulk solution, increased significantly as pressure approached the phase boundary from high pressure. The local composition enhancement was not considerable as pressure was much higher than the critical pressure. In addition, in subcritical region the degree of composition enhancement was much smaller and was not sensitive to pressure in the entire pressure range as the concentration of the cosolvents in the mixed solvents was much higher than the concentration at the critical point of the mixtures.     

Assessment of the potential models of acetone/CO2 and ethanol/CO2 mixtures by computer simulation and thermodynamic integration in liquid and supercritical states

Binary mixtures of CO(2) with ethanol and with acetone are studied by computer simulation, including extensive free energy calculations done by the method of thermodynamic integration, at 313 K, i.e., above the critical point of CO(2) in the entire composition range. The calculations are repeated with three different models of acetone and ethanol, and two models of CO(2). Comparisons of the molar volume of the different systems as well as of the change of their molar volume accompanying the mixing of the two components with experimental data reveal that, among the model pairs tested, the best results are obtained if both components are described by the Transferable Potentials for Phase Equilibria (TraPPE) force field. Around the ethanol/acetone mole fraction of 0.05 all ethanol/CO(2) and almost all acetone/CO(2) model pairs considered predict the existence of a sharp maximum of the molar volume. Due to the lack of experimental data in this composition range, however, these predictions cannot be verified/falsified yet. Most of the model pairs considered also predict limited miscibility of these compounds, as seen from the positive values of the free energy change accompanying their mixing, and the miscibility gap is located at the same composition range as the aforementioned molar volume maximum.     

Structural features of binary mixtures of supercritical CO2 with polar entrainers by molecular dynamics simulation

Computer simulations of supercritical carbon dioxide and its mixtures with polar cosolvents: water, methanol, and ethanol (concentration, 0.125 mole fractions) at T = 318 K and ρ = 0.7 g/cm³ are performed. Atom-atom radial distribution functions are calculated by classical molecular dynamics, while the probability distributions of relative orientation of CO₂ molecules in the first and second coordination spheres describing the geometry of the nearest environment of CO₂ molecules and the trajectories of cosolvent molecules are found using Car-Parrinello molecular dynamics. Based on the latter, the conclusions regarding structure and interactions of polar entrainers in their mixtures with supercritical CO₂ are made. It is shown that the microstructure of carbon dioxide varies only slightly upon the introduction of cosolvents.     

Pressure and temperature dependence of the hydrogen bonding in supercritical ethanol: a computer simulation study

As a step toward deeper insight on the "hydrogen bonding" in supercritical ethanol (scEtOH), we carried out NVT molecular dynamics simulations of the fluid over a wide range of temperatures and pressures. The fluid was studied at SC conditions for which thermodynamic and spectroscopic (NMR, infrared, Raman, dielectric) data are available. The various site-site pair distribution functions (pdf's) were calculated, and their temperature and pressure dependence was obtained. It was found that over the thermodynamic conditions investigated here, scEtOH remains highly structured. Moreover, the characteristic behavior of the first peaks in H-H, O-O, and H-O pdf's reveals that hydrogen bonds still exist in scEtOH. The analysis focuses also on the reorientational dynamics of the bond unit vectors O-H, C-O, and of the permanent dipole moment of the molecules as well as the total dipole moment of the sample. The corresponding Legendre time correlation functions were discussed in connection to the "hydrogen bonding" in the fluid and in the context of experimental results. Specifically, the behavior of the O-H dynamics exhibits the well-known associative nature of the molecules in the system. A further analysis of the hydrogen bonds was carried out, and the degree of aggregation (average number of H-bonds per molecule) was obtained and compared with results from NMR chemical shift studies. Also the estimated monomer and free O-H groups in the fluid were compared with results from IR and Raman vibrational spectroscopy. The percentage analysis fi of the liquid and scEtOH molecules, with i = 0, 1, 2, 3, ... hydrogen bonds per molecule, has been obtained. The results show the existence of small, linear-chain oligomers formed mainly by two molecules, whereas the number of the three body oligomers, and specifically that of four body oligomers in the sample, is relatively small.     

Molecular dynamics simulation of a reverse micelle self assembly in supercritical CO2

In this communication we report on molecular dynamics computer simulations of self-assembly of reverse micelles in supercritical carbon dioxide. The reverse micelles contain perfluoropolyether ammonium carboxylate surfactants and an aqueous core. We observed a quick self-assembly of these micelles over time periods of approximately 5 ns, irrespective of initial conditions. In most cases, the self-assembled perfluorinated reverse micelles have a nice spherical shape and properties consistent with experiments. When the fluorinated surfactant is replaced by its hydrogenated analogue, the assembled aggregate contains a region of direct contact between water and carbon dioxide, indicating that hydrogenated surfactant is not a good agent for creation of microemulsions in water/carbon dioxide mixtures.     

Partial enthalpies and related quantities in mixtures from computer simulation

We report a method of calculating partial molar quantities in mixtures by computer simulation. The method is based on an extension of Widom's potential distribution theorem and provides an alternative way of computing partial enthalpies and volumes.     

Excited-state proton transfer in gas-expanded liquids: the roles of pressure and composition in supercritical CO2/methanol mixtures

Excited-state proton transfer from 5,8-dicyano-2-naphthol to methanol takes place in CO2/methanol mixtures, in the pressure and temperature ranges of supercritical CO2. The efficiency of the proton-transfer step decreases with the pressure. This is assigned to the perturbation of the methanol clusters solvating the naphthol.     

Solvation structure and dynamics for passivated Au nanoparticle in supercritical CO2: a molecular dynamic simulation

All-atomic molecular dynamics simulations have been performed to study the interfacial structural and dynamical properties of passivated gold nanoparticles in supercritical carbon dioxide (scCO(2)). Simulations were conducted for a 55-atom gold nanocore with thiolated perfluoropolyether as the packing ligands. The effect of solvent density and surface coverage on the structural and dynamical properties of the self-assembly monolayer (SAM) has been discussed. The simulation results demonstrate that the interface between nanoparticle and scCO(2) solvent shows a depletion region due to the preclusion of SAM. The presence of scCO(2) solvent around the passivated Au nanoparticle can lead to an enhanced extension of the surface SAM. Under full coverage, the structure and conformation of SAM are insensitive to the density change of scCO(2) fluid. This simulation results clarify the microscopic solvation mechanism of passivated nanoparticles in supercritical fluid medium and is expected to be helpful in understanding the scCO(2)-based nanoparticle dispersion behavior.     

Adsorption and self-assembly of surfactant/supercritical CO2 systems in confined pores: a molecular dynamics simulation

A coarse-grained molecular dynamics simulation has been carried out to study the adsorption and self-organization for a model surfactant/supercritical CO2 system confined in the slit-shape nanopores with amorphous silica-like surfaces. The solid surfaces were designed to be CO2-philic and CO2-phobic, respectively. For the CO2-philic surface, obviously surface adsorption is observed for the surfactant molecules. The various energy profiles were used to monitor the lengthy dynamics process of the adsorption and self-assembly for surfactant micelles or monomers in the confined spaces. The equilibrium properties, including the morphologies and micelle-size distributions of absorbed surfactants, were evaluated based on the equilibrium trajectory data. The interaction between the surfactant and the surface produces an obvious effect on the dynamics rate of surfactant adsorption and aggregation, as well as the final self-assembly equilibrium structures of the adsorbed surfactants. However, for the CO2-phobic surfaces, there are scarcely adsorption layers of surfactant molecules, meaning that the CO2-phobic surface repels the surfactant molecules. It seems to conclude that the CO2 solvent depletion near the interfaces determines the surface repellence to the surfactant molecules. The effect of the CO2-phobic surface confinement on the surfactant micelle structure in the supercritical CO2 has also been discussed. In summary, this study on the microscopic behaviors of surfactant/Sc-CO2 in confined pores will help to shed light on the surfactant self-assembly from the Sc-CO2 fluid phase onto solid surfaces and nanoporous media.     

Tautomeric equilibrium of ethyl acetoacetate in compressed CO2 + ethanol and CO2 + methanol mixtures

Tautomerism equilibrium of ethyl acetoacetate (EAA) in compressed CO2 + methanol and CO2 + ethanol mixtures was studied by UV-Vis spectroscopy at 308.15 K and different pressures. The volume expansion coefficient (alpha) of the solvents at different pressures was also determined. The relative permittivity (epsilon) of CO2 + methanol and CO2 + ethanol mixtures at different conditions was calculated using the Kc and Onsager solvent parameter. The equilibrium constant (Kc) of EAA in the binary mixtures increases considerably with increasing pressure or volume expansion coefficient. The relative permittivity or the polarity of the binary mixtures decreases sharply with increasing volume expansion coefficient in the range of 0 < alpha < 1.5. However, as the volume expansion coefficient exceeds 1.5, the relative permittivity decreases slowly. In other words, the dissolution of CO2 in the polar solvents can reduce the polarity of the solvents significantly in the low volume expansion coefficient range, and the polarity of the solution is not sensitive to the volume expansion coefficient as its value is large enough. The difference in polarity of the two solvents reduces with increasing pressure and becomes negligible after volume expansion coefficient exceeds about 2.5.     

Novel numerical simulation of drug solubility in supercritical CO2 using machine learning technique: Lenalidomide case study

In pharmaceutical industry, finding promising ways to enhance the solubility of disparate types of drugs is an important challenge for the orally administered drug delivery system. Disparate techniques based on drug characteristics, nature of dosage form and properties of excipients have recently been under extensive evaluation all over the world to improve the solubility of poorly water-soluble drugs. Among them, supercritical fluid carbon dioxide (SC-CO₂) has received paramount attentions due to having considerable advantages like cost-effectiveness and low flammability. Lenalidomide belongs is an orally administered anti-cancer agent, which has recently received indication for the treatment of adult patients with different bone marrow-related malignancies such as multiple myeloma, mantle cell lymphoma and follicular lymphoma. Predicting the optimized value of Lenalidomide inside the SC-CO₂ in a wide range of pressure and temperature via developing mathematical models based on artificial intelligence (AI) is the main objective of this paper. In this study, three different machine learning based models are selected to predict and optimized the drug solubility. The available data includes 28 rows of data with two inputs including temperature and pressure and two outputs including density and solubility. Selected models are Kernel Ridge Regression (KRR), least angle regression (LAR), and Multilayer Perceptron (MLP). After optimizing models and comparing the results, the MLP was selected as the primary model of this research. The models illustrated R-squared scores of 0.999 and 0.994 for density and solubility. The maximum errors are also 2.92 and 6.44 × 10^-2 for these outputs, which shows the accuracy and significant generality of the model.     

Density dependence of the entropy and the solvation shell structure in supercritical water via molecular dynamics simulation

We perform molecular dynamics simulations of supercritical water (SCW) with a wide range of densities along a near critical isotherm using the simple point charge extended (SPC/E) pair potential in order to study the entropy and the solvation shell structure around a central water molecule. It is shown that both the translational and orientational two-particle correlation entropy terms can serve as the metrics of the translational-orientational structural orders in water and it is revealed that the translational structural order is very sensitive to the density variation in the gas-like and liquid-like region, while the orientational structural order is much more dependent upon compression in the medium-density SCW region. The comparison of the magnitudes of the full thermodynamic excess entropy and two-particle correlation entropy confirms the recent findings that the many-body terms other than two-body ones also present significant and non-neglectable contributions to the full excess entropy for the highly anomalous fluids like water. The analysis of entropy terms as a function of intermolecular distance and the orientational distribution functions as well as the three-dimensional spatial distribution functions indicate that the structural order occurs only in a much more diffused first solvation shell due to the elongated hydrogen bonds under supercritical conditions. It is revealed that no obvious second or higher neighbor shells occur in SCW, in contrast with the feature of normal liquid water that the anomalous decrease of translational order upon compression occurs mainly in the second shell.     

Synthesis of TiO2 photocatalysts in supercritical CO2 via a non-hydrolytic route

Nanoscaled TiO2 powders with narrow size dispersion were prepared in supercritical carbon dioxide via non-hydrolytic acylation/deacylation of titanium alkoxide precursors with or without tris-fluorination. The microstructures of these powders were characterized by spectroscopic (FTIR, TGA, and XRD), microscopic (SEM or TEM), and surface area (BET) measurements. Photocatalytic oxidation of 1-octanol on these calcined TiO2 powders and on commercial T805 TiO2 suspended in aerated supercritical carbon dioxide revealed relative reactivity controlled by the powder microstructures. Calcined TiO2 prepared from titanium(IV) isopropoxide and trifluoroacetic anhydride was effectively dispersed in aerated supercritical carbon dioxide under stirring and exhibited high photocatalytic oxidation activity.     

Direct calculation of the crystal-melt interfacial free energy via molecular dynamics computer simulation

We review our recent work on the direct calculation of the interfacial free energy, gamma, of the crystal-melt interface via molecular dynamics computer simulation for a number of model systems. The value of gamma as a function of crystal orientation is determined using a thermodynamic integration technique employing moving cleaving walls [Phys. Rev. Lett. 2000, 85, 4751]. The calculation is sufficiently precise to resolve the small anisotropy in gamma, which is crucial in determining the kinetics and morphology of dendritic growth. We report values of gamma for the hard-sphere and Lennard-Jones systems, as well as recent results on the series of inverse-power potentials. For the inverse sixth-, seventh-, and eighth-power systems, we determine gamma for both fcc and bcc crystal structures. For these systems, the bcc-melt gamma is lower than that for fcc crystals by about 25%, consistent with recent experiments and computer simulations on fcc-forming systems that show preferential formation of bcc nuclei in the initial stages of crystallization. In addition, we show that our results give a molecular interpretation of Turnbull's rule, which is the empirical relationship between gamma and the enthalpy of fusion.     

Molecular simulation of interaction between passivated gold nanoparticles in supercritical CO2

Molecular dynamics simulations have been performed to study the potential of mean force (PMF) between passivated gold nanoparticles (NPs) in supercritical CO(2) (scCO(2)). The nanoparticle model consists of a 140 atom gold nanocore and a surface self-assembled monolayer, in which two kinds of fluorinated alkanethiols were considered. The molecular origin of the thermodynamics interaction and the solvation effect has been comprehensively studied. The simulation results demonstrate that increasing the solvent density and ligand length can enhance the repulsive feature of the free energy between the passivated Au nanoparticles in scCO(2), which is in good agreement with previous experimental results. The interaction forces between the two passivated NPs have been decomposed to reveal various contributions to the free energy. It was revealed that the interaction between capping ligands and the interaction between the capping ligands and scCO(2) solvent molecules cooperatively determine the total PMF. A thermodynamic entropy-energy analysis for each PMF contribution was used to explain the density dependence of PMF in scCO(2) fluid. Our simulation study is expected to provide a novel microscopic understanding of the effect of scCO(2) solvent on the interaction between passivated Au nanoparticles, which is helpful to the dispersion and preparation of functional metal nanoparticles in supercritical fluids.     

Hydrogen-bonded clusters and solvate structures in the supercritical CO2-water-o-hydroxybenzoic acid system: the car-parrinello molecular dynamics

Hydrogen-bonded clusters and solvate structures formed by o-hydroxybenzoic acid (o-HBA) and water in supercritical CO₂ were studied (T = 318 K, 348 K, ρ = 0.7 g/cm³). The atom-atom radial distribution functions, coordination numbers, average numbers of hydrogen bonds for individual atomic groups, and power spectrum were calculated by the Car-Parrinello molecular dynamics. Despite the high polarity of the cosolvent, the hydroxyl group of o-HBA predominantly forms intramolecular hydrogen bond, while hydrogen bonds with water involve only the atoms of carboxyl groups. The temperature effect on the stability of these bonds showed itself in different ways. The intermolecular interactions of o-HBA with carbon dioxide were found to be weaker than those with water. It was established that the Lewis acid-Lewis base interactions between CO₂ and the hydroxyl group of the solute increase with increasing temperature. Instantaneous configurations illustrating the temperature effects on the molecular structures were obtained.     

Relationship between the density of supercritical CO_2 + ethanol binary system and its critical properties

The solvent strength and selectivity of supercritical fluids (SCF) can be greatly enhanced by addition of one or two entrainers into the system. The amount of entrainer added is usually less than 5% (mole fraction). However, even with such slight amount, solubility of organic solutes has been observed to increase by several orders magnitude[1]. Therefore, critical pressure and tem-perature data of these supercritical fluid + cosolvent systems are imperative for the reasonable design of effici…     

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.     

The effects of solute-solvent electrostatic interactions on solvation dynamics in supercritical CO2

We present here the results of molecular-dynamics simulation of solvation dynamics in supercritical CO(2) at a temperature of about 1.05T(c), where T(c) is the critical temperature, and at a series of densities ranging from 0.4 to 2.0 of the critical density rho(c). We focus on electrostatic solvation dynamics, representing the electronic excitation of the chromophore as a change in its charge distribution from a quadrupolar-symmetry ground state to a dipolar excited state. Two perturbations are considered, corresponding to different magnitudes of solute excited-state dipoles, denoted as d5 and d8. The d8 solute is more attractive, leading to a larger enhancement in CO(2) clustering upon solute electronic excitation. This has a large impact on solvation dynamics, especially at densities below rho(c). At these densities, solvation dynamics is much slower for the d8 than for the d5 solute. For both solutes, solvation dynamics becomes faster at densities above rho(c) at which solvent clustering diminishes. We show that the slowest solvation time scale is associated with solvent clustering and we relate it to solute-solvent mutual translational diffusion and the extent of change in effective local density resulting from solute electronic excitation.