Aggregation of FeCl2 clusters in supercritical water investigated by molecular dynamics simulations

We have carried out molecular dynamics (MD) simulations of the aggregation of FeCl 2 clusters in supercritical water. The particle formation in systems of 2048 water molecules (rigid SPC/E-model) and 120 Fe (2+) ions and 240 Cl (-) ions has been investigated for 250 ps at five different state points at temperatures from 798 to 873 K and system densities from 0.18 g/cm (3) to 0.13 g/cm (3). We describe the particle growth by means of properties of the largest cluster in a system as well as cluster size averaged and time averaged observables. From preexisting or immediately formed units of Fe (2+)-Cl (-), Fe (2+)-Cl (-) 2, Fe (2+)-Cl (-) 3 etc., the further growth of clusters is dominated by aggregation of such small building blocks. Clusters up to 10 ions in size with large charge imbalances are found during the growth process while a balanced positive to negative charge ratio is found on the average with time and cluster size development. Water molecules are found within the FeCl 2 clusters during the whole time interval covered by the simulations, which is in agreement with the existence of crystal water in FeCl 2 crystals grown from aqueous solutions. The radial distribution functions obtained from the simulation data are in good agreement with experimental results of slightly distorted FeCl 2.4H 2O crystals.     

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.     

Numerical simulation of the solvate structures of acetylsalicylic acid in supercritical carbon dioxide containing polar co-solvents

Hydrogen-bonded complexes of acetylsalicylic acid with polar co-solvents in supercritical carbon dioxide, modified by methanol, ethanol, and acetone of 0.03 mole fraction concentration, are studied by numerical methods of classical molecular dynamics simulation and quantum chemical calculations. The structure, energy of formation, and lifetime of hydrogen-bonded complexes are determined, along with their temperature dependences (from 318 to 388 K at constant density of 0.7 g cm^?3). It is shown that the hydrogen bonds between acetylsalicylic acid and methanol are most stable at 318 K and are characterized by the highest value of absolute energy. At higher supercritical temperatures, however, the longest lifetime is observed for acetylsalicylic acid–ethanol complexes. These results correlate with the known literature experimental data showing that the maximum solubility of acetylsalicylic acid at density values close to those considered in this work and at temperatures of 318 and 328 K is achieved when using methanol and ethanol as co-solvents, respectively.     

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.     

Solubilities of lactic acid and 2-hydroxyhexanoic acid in supercritical CO2

Solubilities of lactic acid and 2-hydroxyhexanoic acid in supercritical CO₂ have been measured at T=(311 or 313, 318, 328 K) in the pressure range from 50 to 200 bar. The measurements have been performed using a flow-type apparatus. The solute solubility in compressed carbon dioxide increased with pressure at all investigated temperatures. At pressures below 130 bar, a solubility decrease on temperature increase was observed. An accurate correlation method for the solubility of low volatile substances in supercritical CO₂ has been applied for the interpolation of the experimental results.     

Hydration of hydrogentungstate anions at different ph conditions: a car-parrinello molecular dynamics study

Standard density functional theory calculations with a continuous model of solvation as well as Car-Parrinello molecular dynamics simulations with explicit solvent molecules are carried out to analyze the effect of the pH of the solution on the coordination sphere of the W (VI) ion. Both methodologies agree in predicting an expansion of the coordination sphere of the W (VI) ion upon a decrease in the pH. Continuous solvation models, however, are unable to predict as stable some structural isomers of a hydrated hydrogentungstate anion and tungstic acid.     

Size and structure of molecular clusters in supercritical water

The type and topology of hydrogen-bonded molecular clusters of water are investigated by the molecular dynamics method for five models of water in supercritical conditions. Small clusters (of the order of 10 molecules) are present in all models, even at densities of less than 0.2g/cm³. When the density increases, a phase transition occurs from vapor-like to fluid-like state. Among small clusters, linear structures are predominant.     

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.     

Investigations of the hydrogen bond in the crystals of tropolone and thiotropolone via car-parrinello and path integral molecular dynamics

Car-Parrinello and path integrals molecular dynamics (CPMD and PIMD) simulations were carried out for the 10π-electron aromatic systems: 2-hydroxy-2,4,6-cycloheptatrien-1-one, commonly known as Tropolone (I) and 2-hydroxy-2,4,6-cycloheptatriene-1-thione, called Thiotropolone (II) in vacuo and in the solid state. The extremely fast proton transfer (FPT) and “prototropy” tautomerism in the keto-enol (thione-enethiol) systems have been analyzed on the basis of CPMD and PIMD methods level. Comparisons of two-dimensional (2D) free-energy landscapes of reaction coordinate δ-parameter and R_O…O or R_O…S distances shows that the OH… tautomer to be more favorable in the Thiotropolone. The hydrogen between the oxygen and the sulfur atoms adopts a starkly asymmetrical position in the double potential well. The values of the energy barriers for the FPT were calculated and suggested a strong hydrogen bond with low barrier for FPT mechanism. These studies and the 2D average index of π-delocalization 〈λ〉 landscape of time evolutions of R_O1…O2 and R_C7O2 or R_C7S1 distances for the both crystals indicate that hydrogen bonds in the crystals of Tropolone (I) and Thiotropolone (II) have characteristic properties for the type of bonding model resonance-assisted hydrogen bonds and also low-barrier hydrogen bonds. In the crystal of the Thiotropolone (II), we found the hydrogen bond O H…S existing without the equilibrium of the two tautomers whereas in the crystal of the Tropolone (I) has been confirmed the hydrogen bond O H…O existing with the equilibrium of the two tautomers. It was also found the significant differences in frequency, speed, and the image of the FPT in the studied crystals. © 2018 Wiley Periodicals, Inc.     

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.     

Phase behavior in model homopolymer/CO2 and surfactant/CO2 systems: discontinuous molecular dynamics simulations

Discontinuous molecular dynamics simulations are performed on surfactant (HmTn)/solvent systems modeled as a mixture of single-sphere solvent molecules and freely jointed surfactant chains composed of m slightly solvent-philic head spheres (H) and n solvent-philic tail spheres (T), all of the same size. We use a square-well potential to account for the head-head, head-solvent, tail-tail, and tail-solvent interactions and a hard-sphere potential for the head-tail and solvent-solvent interactions. We first simulate homopolymer/supercritical CO2 (scCO2) systems to establish the appropriate interaction parameters for a surfactant/scCO2 system. Next, we simulate surfactant/scCO2 systems and explore the effect of the surfactant volume fraction, packing fraction, and temperature on the phase behavior. The transition from the two-phase region to the one-phase region is located by monitoring the contrast structure factor of the equilibrated surfactant/scCO2 system, and the micelle to unimer transition is located by monitoring the aggregate size distribution of the equilibrated surfactant/scCO2 system. We find a two-phase region, a micelle phase, and a unimer phase with increasing packing fraction at fixed temperature or with increasing temperature at fixed packing fraction. The phase diagram for the surfactant/scCO2 system in the surfactant volume fraction-packing fraction plane and the density dependence of the critical micelle concentration are in qualitative agreement with experimental observations. The phase behavior of a surfactant/scCO2 system can be directly related to the solubilities of the corresponding homopolymers that serve as the head and tail blocks for the surfactant. The influence of surfactant structure (head and tail lengths) on the phase transitions is explored.     

Investigation of the local composition enhancement and related dynamics in supercritical CO2-cosolvent mixtures via computer simulation: the case of ethanol in CO2

The supercritical mixture ethanol-carbon dioxide (EtOH-CO2) with mole fraction of ethanol X(EtOH) congruent with 0.1 was investigated at 348 K, by employing the molecular dynamics simulation technique in the canonical ensemble. The local intermolecular structure of the fluid was studied in terms of the calculated appropriate pair radial distribution functions. The estimated average local coordination numbers and mole fractions around the species in the mixture reveal the existence of local composition enhancement of ethanol around the ethanol molecules. This finding indicates the nonideal mixing behavior of the mixture due to the existence of aggregation between the ethanol molecules. Furthermore, the local environment redistribution dynamics have been explored by analyzing the time correlation functions (TCFs) of the total local coordination number (solvent, cosolvent) around the cosolvent molecules in appropriate parts. The analysis of these total TCFs in the auto-(solvent-solvent, cosolvent-cosolvent) and cross-(solvent-cosolvent, cosolvent-solvent) TCFs has shown that the time dependent redistribution process of the first solvation shell of ethanol is mainly determined by the redistribution of the CO2 solvent molecules. These results might be explained on the basis of the CO2-CO2 and EtOH-CO2 intermolecular forces, which are sufficiently weaker in comparison to the EtOH-EtOH hydrogen bonding interactions, creating in this way a significantly faster redistribution of the CO2 molecules in comparison with EtOH. Finally, the self-diffusion coefficients and the single reorientational dynamics of both the cosolvent and solvent species in the mixture have been predicted and discussed in relationship with the local environment around the species, which in the case of the EtOH molecules seem to be strongly affected.     

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.     

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.     

Infrared and molecular-dynamics studies of the rotational dynamics of water highly diluted in supercritical CO2

Far-infrared (FIR) and mid-infrared (MIR) profiles of D2O infinitely dilute in supercritical CO2 have been studied using molecular-dynamics simulations. For this purpose, we have proposed an intermolecular potential model taking implicitly into account electron donor-acceptor (EDA) interactions between water and CO2 evaluated from ab initio calculations of the intermolecular potential-energy surface (IPS). Interaction-induced dipole mechanisms have been also taken into account in addition to the water permanent dipole to evaluate the simulated FIR profiles of water and CO2 polarizable molecules. They were found to play a minor role in the genesis of the FIR profiles of water/CO2 under supercritical conditions. The analysis of the reorientational dynamics of D2O shows that the rotational dynamics of water is weakly anisotropic due to the EDA interactions which affect more specifically the reorientational motions of the C2 symmetry axis of solute. These results have been used to assess the contribution of the vibrational relaxation in the experimental mid-infrared profiles associated with the nu1 symmetric and nu3 antisymmetric stretching and nu2 bending modes of D2O. It was found that the rotational dynamics mainly contribute to the broadening of the infrared (IR) profiles. Nevertheless, the vibrational processes play a role in the frequency shifts of the band centers and the relative intensity enhancements of the nu1 and nu3 modes of D2O. In particular, the EDA interactions between water and CO2 lead to the appearance of a well-defined IR band of the nu1 mode of D2O. Finally, a comparison with another model taking only into account dipole-quadrupole electrostatic interactions between water and CO2 molecules clearly reveals that EDA interactions have to be considered to reproduce both MIR and FIR measurements. From this point of view CO2 can be classified on a hydrophilic solvent scale based upon the solubility criterion as an intermediate solvent between "inert" xenon and carbon tetrachloride.     

Nucleation of NaCl nanoparticles in supercritical water: molecular dynamics simulations

Formation of NaCl nanoparticles in supercritical water is studied using molecular dynamics simulation method. We have simulated particle nucleation and growth in NaCl-H2O fluids, with salt concentration of 5.1 wt %, in the temperature and density range of 673-1073 K and 0.17-0.34 g/cm(3), respectively. The cluster size distributions, the size of critical nuclei and cluster lifetimes are reported. The size distribution of emerging clusters shows a very strong dependence on the system's density, with larger clusters forming at lower densities. Clusters consisting of approximately 14-24 ions appear critical for the thermodynamic states examined. The local structures of critical clusters are found to be amorphous. The lifetime values for clusters containing more than 20 ions are in the range of 10-50 ps. We have calculated the NaCl nucleation rates, which appear to be on the order of 10(28) cm(-3) s(-1).     

Fluorinated surfactants in supercritical CO2

Liquid or supercritical carbon dioxide has important environmental and economic advantages over petrochemical solvents currently used for industrial processes. However, low solubility in CO₂, particularly of polar compounds, is a hurdle to its implementation as an acceptable alternative. These solubility problems have been overcome by employing specialised fluorinated surfactants to stabilise water nano-droplets as water-in-supercritical/liquid CO₂ microemulsions. Such novel microemulsions can now facilitate innovative ‘green-and-clean’ applications of carbon dioxide technology.     

Extraction of ferulic acid from Angelica sinensis with supercritical CO2

Extraction of pharmacologically active ingredient of ferulic acid from the root of Angelica sinensis with supercritical CO2 was investigated. The experimental results show that the extract yields were 0.87-4.06% at temperatures from 45 to 65 degrees Celsius and pressures from 30 to 50 MPa, and the maximum content of ferulic aicd in the extracts was about 0.35-0.37%, which is lower than that of 0.61-0.85% by conventional percolation methods. Ethanol was used as co-solvent in different ratios to raw materials in order to increase the content of ferulic acid in the extracts. The experimental results show that both the extract yields and the content of ferulic acid in the extracts increase greatly compared with pure CO2 extraction. When the ratio of ethanol to the raw material was 1.6, the content of ferulic acid in the extracts was 0.91-1.27%, indicating that supercritical fluid extraction (SFE) with CO2 in the presence of suitable co-solvent is superior to percolation in extracting polar ferulic acid from Angelica sinensis.     

The dynamics of phase changes in molecular clusters

The significant advantages offered by systems of molecular clusters in the study of homogeneous nucleation are discussed. Determinations of nucleation rates in clusters can be followed experimentally in supersonic jets or computationally in molecular dynamics simulations. Extraordinarily high rates may be encountered, both in freezing and in solid-state transitions. From such information can be inferred the interfacial free energies, σ_sl and σ_SS, mechanisms of solid-state transitions. and an explanation of why certain crystalline phases not found in bulk systems can be seen in large molecular clusters.     

Heterometallic heterochalcogenmethylarsenide clusters: Synthesis,molecular structures,and thermolysis

Reactions of the arsinechalcogenide complexes [Fe₃(μ₃-X)(μ₃-AsCH₃)(CO)₉] (X = Se (Ia) or Te (Ib)) with (PPh₃)₂Pt(PhC≡CPh) (transmetalation reaction) and Cp₂Cr₂(SCME₃)₂S (Cp = π-C₅H₅) (photochemical reaction) gave the heterometallic (heterochalcogen)(methylarsine) clusters [(PPh₃)₂Pt(μ₃-X)(μ₃-AsCH₃)Fe₂(CO)₆] (II and III, respectively), as well as Fe₃(μ₃-X)(μ₃-AsCH₃)(CO)₈(C₅H₅)₂Cr₂(μ₃-S)(μ₂-S ^t Bu)₂ (IV and V, respectively). The structures of complexes II, IV, and V were determined by X-ray diffraction analysis. Thermolysis of all the complexes yielded no metal carbides or oxides.