A set of molecular models for carbon monoxide and halogenated hydrocarbons

Ionic Liquids for the Separation of Fluorocarbon Refrigerant Mixtures

This review discusses the research being performed on ionic liquids for the separation of fluorocarbon refrigerant mixtures. Fluorocarbon refrigerants, invented in 1928 by Thomas Midgley Jr., are a unique class of working fluids that are used in a variety of applications including refrigeration. Fluorocarbon refrigerants can be categorized into four generations: chlorofluorocarbons, hydrochlorofluorocarbons, hydrofluorocarbons, and hydrofluoroolefins. Each generation of refrigerants solved a key problem from the previous generation; however, each new generation has relied on more complex mixtures that are often zeotropic, near azeotropic, or azeotropic. The complexity of the refrigerants used and the fact that many refrigerants form azeotropes when mixed makes handling the refrigerants at end of life extremely difficult. Today, less than 3% of refrigerants that enter the market are recycled. This is due to a lack of technology in the refrigerant reclaim market that would allow for these complex, azeotropic refrigerant mixtures to be separated into their components in order to be effectively reused, recycled, and if needed repurposed. As the market for recovering and reclaiming refrigerants continues to grow, there is a strong need for separation technology. Ionic liquids show promise for separating azeotropic refrigerant mixtures as an entrainer in extractive distillation process. Ionic liquids have been investigated with refrigerants for this application since the early 2000s. This review will provide a comprehensive summary of the physical property measurements, equations of state modeling, molecular simulations, separation techniques, and unique materials unitizing ionic liquids for the development of an ionic-liquid-based separation process for azeotropic refrigerant mixtures.

Collisions between CO, CO[Formula: see text], H[Formula: see text]O and Ar ice nanoparticles compared by molecular dynamics simulation

Molecular dynamics simulations are used to study collisions between amorphous ice nanoparticles consisting of CO, CO[Formula: see text], Ar and H[Formula: see text]O. The collisions are always sticking for the nanoparticle size (radius of 20 nm) considered. At higher collision velocities, the merged clusters show strong plastic deformation and material mixing in the collision zone. Collision-induced heating influences the collision outcome. Partial melting of the merged cluster in the collision zone contributes to energy dissipation and deformation. Considerable differences exist-even at comparable collision conditions-between the ices studied here. The number of ejecta emitted during the collision follows the trend in triple-point temperatures and increases exponentially with the NP temperature.

Bayesian-Inference-Driven Model Parametrization and Model Selection for 2CLJQ Fluid Models

A high level of physical detail in a molecular model improves its ability to perform high accuracy simulations but can also significantly affect its complexity and computational cost. In some situations, it is worthwhile to add complexity to a model to capture properties of interest; in others, additional complexity is unnecessary and can make simulations computationally infeasible. In this work, we demonstrate the use of Bayesian inference for molecular model selection, using Monte Carlo sampling techniques accelerated with surrogate modeling to evaluate the Bayes factor evidence for different levels of complexity in the two-centered Lennard-Jones + quadrupole (2CLJQ) fluid model. Examining three nested levels of model complexity, we demonstrate that the use of variable quadrupole and bond length parameters in this model framework is justified only for some chemistries. Through this process, we also get detailed information about the distributions and correlation of parameter values, enabling improved parametrization and parameter analysis. We also show how the choice of parameter priors, which encode previous model knowledge, can have substantial effects on the selection of models, penalizing careless introduction of additional complexity. We detail the computational techniques used in this analysis, providing a roadmap for future applications of molecular model selection via Bayesian inference and surrogate modeling.

Machine Learning Directed Optimization of Classical Molecular Modeling Force Fields

Accurate force fields are necessary for predictive molecular simulations. However, developing force fields that accurately reproduce experimental properties is challenging. Here, we present a machine learning directed, multiobjective optimization workflow for force field parametrization that evaluates millions of prospective force field parameter sets while requiring only a small fraction of them to be tested with molecular simulations. We demonstrate the generality of the approach and identify multiple low-error parameter sets for two distinct test cases: simulations of hydrofluorocarbon (HFC) vapor-liquid equilibrium (VLE) and an ammonium perchlorate (AP) crystal phase. We discuss the challenges and implications of our force field optimization workflow.

Molecular Simulations of Adsorption and Energy Storage of R1234yf, R1234ze(z), R134a, R32, and their Mixtures in M-MOF-74 (M = Mg, Ni) Nanoparticles

The refrigerant circulation heat can be enhanced through the mutual transformation between thermal energy and surface energy during the adsorption and separation process of fluid molecules in porous materials. In this paper, the adsorption and energy storage of R1234ze(z), R1234yf, R32 and R134a, as well as their mixed refrigerants in Mg-MOF-74 and Ni-MOF-74 nanoparticles were investigated by means of molecular dynamics simulations and grand canonical Monte Carlo simulations. The results suggested that, in the case of pure refrigerant adsorption, the adsorption quantities of R32 and R134a in MOFs were higher than those of R1234yf and R1234ze(z). However, in the case of saturation adsorption, the desorption heat of R32 was lower than that of R1234yf and R1234ze(z). The addition of MOF-74 nanoparticles (NPs) could enhance the energy storage capacity of the pure refrigerant; besides, R1234yf and R1234ze(z) nanofluids had superior enhancement effect to that of R32 nanofluid. In mixed refrigerant adsorption, the adsorption quantities of R1234ze(z) and R1234yf were lower than those of R32 and R134a; with the increase in temperature, the adsorption of R1234ze(z) and R1234yf showed a gradually increasing trend, while that of R32 was gradually decreased.

Relative Populations of Some Tautomeric Forms of 2'-Deoxyguanosine-5-Fluorouridine Mismatch

The importance of the 2'-deoxyguanosine-uridine mispair as the most occurring mismatch in transcriptional studies of RNAs from DNAs is multiplied when 5-halo-substituted uridine species cause a serious increase in the probability of its occurrence. Many studies relate this higher probability to the existence of possible tautomeric and ionic forms of its constituent bases. According to these statements, relative populations of mismatches between 5-fluorouridine and both keto and enol forms of 2'-deoxyguanosine are computed by using a conformational search. In order to have a complete scan of all of the highly probable conformers in a moderate computational time, an extensive conformational search methodology is employed here, which benefits from the advantages of both the molecular dynamics simulations and quantum mechanics calculations. The population of an enolic tautomer of normal wobble orientation is about 0.057% of that of its keto tautomer, whereas the population of an enolic tautomer of reverse wobble orientation is about 0.0054% of that of its keto tautomer. Totally, the reverse wobble orientation is about six times more populated than the normal wobble orientation. Calculated populations are in good agreement with experimental populations of closely related compounds. The reliability of the applied methodology is certified, in part, by a good agreement obtained between some experimental data and corresponding Boltzmann-weighted average data of most probable conformers such as NMR parameters. The validation of this methodology is certified with high accuracy by applying it on the substituted diuridine pairs, where experimental populations are available. Not only are the calculated populations and NMR parameters of this test in very good agreement with the experimental data, but also they are free of the ambiguities mentioned by experimentalists.

Gas Separation through Bilayer Silica, the Thinnest Possible Silica Membrane

Membrane-based gas separation processes can address key challenges in energy and environment, but for many applications the permeance and selectivity of bulk membranes is insufficient for economical use. Theory and experiment indicate that permeance and selectivity can be increased by using two-dimensional materials with subnanometer pores as membranes. Motivated by experiments showing selective permeation of H₂/CO mixtures through amorphous silica bilayers, here we perform a theoretical study of gas separation through silica bilayers. Using density functional theory calculations, we obtain geometries of crystalline free-standing silica bilayers (comprised of six-membered rings), as well as the seven-, eight-, and nine-membered rings that are observed in glassy silica bilayers, which arise due to Stone-Wales defects and vacancies. We then compute the potential energy barriers for gas passage through these various pore types for He, Ne, Ar, Kr, H₂, N₂, CO, and CO₂ gases, and use the data to assess their capability for selective gas separation. Our calculations indicate that crystalline bilayer silica, which is less than a nanometer thick, can be a high-selectivity and high-permeance membrane material for ³He/⁴He, He/natural gas, and H₂/CO separations.

Vapor Liquid Equilibria of Hydrofluorocarbons Using Dispersion-Corrected and Nonlocal Density Functionals

Recent developments in dispersion corrected and nonlocal density functionals are aimed at accurately capturing dispersion interactions, a key shortcoming of local and semilocal approximations of density functional theory. These functionals have shown significant promise for dimers and small clusters of molecules as well as crystalline materials. However, their efficacy for predicting vapor liquid equilibria is largely unexplored. In this work, we examine the accuracy of dispersion-corrected and nonlocal van der Waals functionals by computing the vapor liquid coexistence curves (VLCCs) of hydrofluoromethanes. Our results indicate that the PBE-D3 functional performs significantly better in predicting saturated liquid densities than the rVV10 functional. With the PBE-D3 functional, we also find that as the number of fluorine atoms increase in the molecule, the accuracy of saturated liquid density prediction improves as well. All the functionals significantly underpredict the saturated vapor densities, which also result in an underprediction of saturated vapor pressure of all compounds. Despite the differences in the bulk liquid densities, the local microstructures of the liquid CFH3 and CF2H2 are relatively insensitive to the density functional employed. For CF3H, however, rVV10 predicts slightly more structured liquid than the PBE-D3 functional.

Force-Field Parameters from the SAFT-γ Equation of State for Use in Coarse-Grained Molecular Simulations

A description of fluid systems with molecular-based algebraic equations of state (EoSs) and by direct molecular simulation is common practice in chemical engineering and the physical sciences, but the two approaches are rarely closely coupled. The key for an integrated representation is through a well-defined force field and Hamiltonian at the molecular level. In developing coarse-grained intermolecular potential functions for the fluid state, one typically starts with a detailed, bottom-up quantum-mechanical or atomic-level description and then integrates out the unwanted degrees of freedom using a variety of techniques; an iterative heuristic simulation procedure is then used to refine the parameters of the model. By contrast, with a top-down technique, one can use an accurate EoS to link the macroscopic properties of the fluid and the force-field parameters. We discuss the latest developments in a top-down representation of fluids, with a particular focus on a group-contribution formulation of the statistical associating fluid theory (SAFT-γ). The accurate SAFT-γ EoS is used to estimate the parameters of the Mie force field, which can then be used with confidence in direct molecular simulations to obtain thermodynamic, structural, interfacial, and dynamical properties that are otherwise inaccessible from the EoS. This is exemplified for several prototypical fluids and mixtures, including carbon dioxide, hydrocarbons, perfluorohydrocarbons, and aqueous surfactants.

Microscopic structure of liquid 1-1-1-2-tetrafluoroethane (R134a) from Monte Carlo simulation

1-1-1-2-tetrafluoroethane (R134a) is one of the most commonly used refrigerants. Its thermophysical properties are important for evaluating the performance of refrigeration cycles. These can be obtained via computer simulation, with an insight into the microscopic structure of the liquid, which is not accessible to experiment. In this paper, vapour-liquid equilibrium properties of R134a and its liquid microscopic structure are investigated using coupled-decoupled configurational-bias Monte Carlo simulation in the Gibbs ensemble, with a recent potential [J. Phys. Chem. B 2009, 113, 178]. We find that the simulations agree well with the experimental data, except at the vicinity of the critical region. Liquid R134a packs like liquid argon, with a coordination number in the first solvation shell of 12 at 260 K. The nearest neighbours prefer to be localized in three different spaces around the central molecule, in such a manner that the dipole moments are in a parallel alignment. Analysis of the pair interaction energy shows clear association of R134a molecules, but no evidence for C-HF type hydrogen bonding is found. The above findings should be of relevance to a broad range of fluoroalkanes.

The role of higher CO-multipole moments in understanding the dynamics of photodissociated carbonmonoxide in myoglobin

The influence of electrostatic multipole moments up to hexadecapole on the dynamics of photodissociated carbon monoxide (CO) in myoglobin is investigated. The CO electrostatic potential is expressed as an expansion into atomic multipole moments of increasing order up to octopole which are obtained from a distributed multipole analysis. Three models with increasingly accurate molecular multipoles (accurate quadrupole, octopole, and hexadecapole moments, respectively) are developed and used in molecular dynamics simulations. All models with a fluctuating quadrupole moment correctly describe the location of the B-state whereas the sign of the octopole moment differentiates between the Fe...CO and Fe...OC orientation. For the infrared spectrum of photodissociated CO, considerable differences between the three electrostatic models are found. The most detailed electrostatic model correctly reproduces the splitting, shift, and width of the CO spectrum in the B-state. From an analysis of the trajectories, the spectroscopic B(1) and B(2) states are assigned to the Fe...CO and Fe...OC substates, respectively.

Microscopic and thermodynamic properties of the HFA134a-water interface: atomistic computer simulations and tensiometry under pressure

A combined computational and experimental approach is used to determine the interfacial thermodynamic and structural properties of the liquid 1,1,1,2-tetrafluoroethane (HFA134a)-vapor and liquid HFA134a-water (HFA134a|W) interfaces at 298 K and saturation pressure. Molecular dynamics (MD) computer simulations reveal a stable interface between HFA134a and water. The "10-90" interfacial thickness is comparable with those typically reported for organic-water systems. The interfacial tension of the HFA134a|W interface obtained from the pressure tensor analysis of the MD trajectory is in good agreement with the experimental value determined using in situ high-pressure tensiometry. These results indicate that the potential models utilized are capable of describing the intermolecular interactions between these two fluids. The tension of the HFA134a|W interface is significantly lower than those typically observed for conventional oil-water interfaces and similar to that of the compressed CO(2)-water interface, observed at moderate CO(2) pressures. The MD and tensiometric results are also compared and contrasted with the HFA134a|W and chlorofluorocarbon-water tension values estimated from a parametric relationship. This represents the first report of the interfacial and microscopic properties of the (propellant) hydrofluoroalkanes (HFA)|W interface. The results presented here are of relevance in the design of surfactants capable of forming and stabilizing water-in-HFA microemulsions. Reverse aqueous microemulsions in HFA-based pressurized metered-dose inhalers are candidate formulations for the systemic delivery of biomolecules to and through the lungs.

Photoinduced Isomerization Kinetics of Diiodomethane in Supercritical Fluid Solution: Local Density Effects

The density dependence of diiodomethane photoinduced isomerization in supercritical (sc) CO2, CHF3, and C2H6 was investigated by transient absorption spectroscopy, covering a fluid density range from 0.7 to 2.5 (in reduced units). The solvent-caged photoproduct iso-diiodomethane is formed even at the lowest density, and its yield increases about 4-fold over the whole range. At the same time, isomer formation rate constants increase by roughly an order of magnitude and show little variation between CO2, C2H6, and CHF3. Furthermore, the formation rate constant decreases significantly with increasing excitation energy. We propose an isomer formation mechanism involving a rapidly established preequilibrium between a solvent-caged iodine atom-methyliodide radical pair and a loosely bound iodine-methyliodide radical complex, from which the reaction subsequently proceeds to the isomer. The latter step seems to be controlled by collisional stabilization of the initially hot radical moiety, as the formation rate constant increases linearly with sc solvent viscosity. The model predicts a quadratic dependence of relative isomer yield on fluid density. A corresponding correlation is found with the local fluid density, calculated via solute-solvent radial distribution functions obtained from molecular dynamics (MD) simulations.