Hydrogen Bonding of Methanol in Supercritical CO2 Studied by 13C Nuclear Spin−Lattice Relaxation

Existence of a liquid-liquid phase transition in methanol

A simple model is constructed to study the phase diagram and thermodynamic properties of methanol, which is described as a dimer of an apolar sphere mimicking the methyl group and a sphere with core-softened potential as the hydroxyl group. Performing classical Monte Carlo simulations, we obtained the phase diagram, showing a second critical point between two different liquid phases. Evaluating systems with a different number of particles, we extrapolate to infinite size in accordance with Ising universality class to obtain bulk values for critical temperature, pressure, and density. Strong evidence that the structure of the liquid changes upon transition from high- to low-density phase was provided. From the experimentally determined hydrogen bond strength and length in methanol and water, we propose where the second critical point of methanol should be.

Hydrogen bond, electron donor-acceptor dimer, and residence dynamics in supercritical CO(2)-ethanol mixtures and the effect of hydrogen bonding on single reorientational and translational dynamics: A molecular dynamics simulation study

The hydrogen bonding and dynamics in a supercritical mixture of carbon dioxide with ethanol as a cosolvent (X(ethanol) approximately 0.1) were investigated using molecular dynamics simulation techniques. The results obtained reveal that the hydrogen bonds formed between ethanol molecules are significantly more in comparison with those between ethanol-CO(2) molecules and also exhibit much larger lifetimes. Furthermore, the residence dynamics in the solvation shells of ethanol and CO(2) have been calculated, revealing much larger residence times for ethanol molecules in the ethanol solvation shell. These results support strongly the ethanol aggregation effects and the slow local environment reorganization inside the ethanol solvation shell, reported in a previous publication of the authors [Skarmoutsos et al., J. Chem. Phys. 126, 224503 (2007)]. The formation of electron donor-acceptor dimers between the ethanol and CO(2) molecules has been also investigated and the calculated lifetimes of these complexes have been found to be similar to those corresponding to ethanol-CO(2) hydrogen bonds, exhibiting a slightly higher intermittent lifetime. However, the average number of these dimers is larger than the number of ethanol-CO(2) hydrogen bonds in the system. Finally, the effect of the hydrogen bonds formed between the individual ethanol molecules on their reorientational and translational dynamics has been carefully explored showing that the characteristic hydrogen bonding microstructure obtained exhibits sufficiently strong influence upon the behavior of them.

Novel 1,3-dipolar cycloadditions of dinitraminic acid: implications for the chemical stability of ammonium dinitramide

Density functional theory at the B3LYP/6-31+G(d,p) level and ab initio calculations at the CBS-QB3 level have been used to analyze 1,3 dipolar cycloaddition reactions of dinitraminic acid (HDN) and its proton transfer isomer (HO(O)NNNO2). It is shown that the nitro group of HDN and the -N-N=O functionality of the isomer react readily with carbon-carbon double bonds. Cycloadditions of HDN are compared with the corresponding reactions with azides and nitrile oxides as 1,3 dipoles. It is shown that the reactivities of HDN and its proton transfer isomer decrease with increasing electron withdrawing power of the substituents adjacent to the carbon-carbon double bond. In contrast, for azides and nitrile oxides, the highest reactivity is obtained with dipolarophiles with strongly electron withdrawing substituents. The observed reactivity trends allow for the design of unsaturated compounds that are highly reactive toward azides and chemically inert toward dinitramides. This may be of relevance for the development of binder materials for ammonium dinitramide based propellants.

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.

Molecular Interactions and CO2-Philicity in Supercritical CO2. A High-Pressure NMR and Molecular Modeling Study of a Perfluorinated Polymer in scCO2

The carbon and fluorine chemical shifts of mixtures of carbon dioxide and Krytox, a carboxylic acid end-capped perfluorinated polyether used as stabilizer for the dispersion polymerization of methyl methacrylate, have been studied using high-pressure, high-resolution nuclear magnetic resonance. 13C and 19F spectra were measured in the density region between 0.54 and 0.73 g.cm(-3) at 334 K for different solutions of Krytox in scCO2 (0.22, 1.13 and 1.72 w/w %). An in-house developed high-pressure apparatus with the capability to change in situ the sample composition was used for this purpose using a 10 mm polyether ketone NMR tube. The nature of CO2-Krytox interaction was assessed both by comparing the CO2 deltaC variation of neat CO2 with that of mixtures with increasing surfactant composition and by the analysis of Krytox 19F corrected chemical shifts in terms of medium magnetic susceptibility. Ab initio calculations, at the second-order Møller-Plesset level of theory to include the effects of electron correlation, were performed to access and compare the nature of the interactions between CO2 and perfluorinated and nonfluorinated analogue model molecules. Both experimental 13C and 19F HP-NMR results and molecular modeling studies support a F...CO2 site-specific Lewis acid-Lewis base interaction model. A positive entropic variation for the formation of CO2-fluorinated solute complex is advanced as an explanation for the higher solubility of perfluorinated molecules when compared to the nonfluorinated analogues.

Elucidating the Vibrational Spectra of Hydrogen-Bonded Aggregates in Solution: Electronic Structure Calculations with Implicit Solvent and First-Principles Molecular Dynamics Simulations with Explicit Solvent for 1-Hexanol in n-Hexane

Fourier transform infrared spectroscopy is a popular method for the experimental investigation of hydrogen-bonded aggregates, but linking spectral information to microscopic information on aggregate size distribution and aggregate architecture is an arduous task. Static electronic structure calculations with an implicit solvent model, Car-Parrinello molecular dynamics (CPMD) using the Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and classical molecular dynamics simulations for the all-atom version of the optimized parameters for liquid simulations (OPLS-AA) force field were carried out for an ensemble of 1-hexanol aggregates solvated in n-hexane. The initial configurations for these calculations were size-selected from a distribution of aggregates obtained from a large-scale Monte Carlo simulation. The vibrational spectra computed from the static electronic structure calculations for monomers and dimers and from the CPMD simulations for aggregates up to pentamers demonstrate the extent of the contribution of dangling or nondonating hydroxyl groups found in linear and branched aggregates to the "monomeric" peak. Furthermore, the computed spectra show that there is no simple relationship between peak shift and aggregate size nor architecture, but the effect of hydrogen-bond cooperativity is shown to differentiate polymer-like (cooperative) and dimer-like (noncooperative) hydrogen bonds in the vibrational spectrum. In contrast to the static electronic structure calculations and the CPMD simulations, the classical molecular dynamics simulations greatly underestimate the vibrational peak shift due to hydrogen bonding.

Binary phase behavior and aggregation of dilute methanol in supercritical carbon dioxide: A Monte Carlo simulation study

Configurational-bias Monte Carlo simulations in the Gibbs and isobaric-isothermal ensembles using the transferable potentials for phase equilibria force field were carried out to investigate the thermophysical properties of mixtures containing supercritical carbon dioxide and methanol. The binary vapor-liquid coexistence curves were calculated at 333.15 and 353.15 K and are in excellent agreement with experimental measurements. The self-association of methanol in supercritical carbon dioxide was investigated over a range of temperatures and pressures near the mixture critical point. The temperature dependence of the equilibrium constants for the formation of hydrogen-bonded aggregates (from dimer to heptamer) allowed for the determination of the enthalpy of hydrogen bonding, DeltaHHB, in supercritical carbon dioxide with values for DeltaHHB of about 15 kJ mol(-1) falling within the range of previously proposed values. No strong pressure dependence was observed for the formation of aggregates. Apparently the decrease of the entropic penalty and of the enthalpic benefit upon increasing pressure or solvent density mostly cancel each other's effect on aggregate formation.