Monte Carlo simulations of vapour–liquid phase equilibrium and microstructure for the system containing azeotropes

Abstract

Phase equilibrium data of the mixtures including alcohols, esters and organic acids are of first interest particularly to design and optimise biodiesel production and reactive distillation processes. In this work, vapour–liquid phase equilibrium of these systems was simulated at low pressure using Gibbs ensemble Monte Carlo method. All Lennard–Jones parameters of pseudo-atoms involved in the systems were derived from previous parametrisations of TraPPE-UA force field. The Fourier coefficients of dihedrals encountered in ethyl acetate molecule have been obtained from the quantum calculations. Using this force field, temperature-composition diagrams are well reproduced for ethyl acetate + ethanol, ethyl acetate + methanol at 70.00 kPa and ethyl acetate + acetic acid mixtures at 77.33 kPa. The transferability of this force field to mixtures in these systems is noticeable. Analysis of the microstructure for the ethyl acetate + ethanol and ethyl acetate + acetic acid mixtures was presented. We found that the hydrogen bond networks consist of autoassociation and cross-association and autoassociation occupies the main position as compared with cross association in the ethyl acetate + ethanol mixture. OC_HAc–H_HAc and OC_EtOAc–H_HAc hydrogen bond interactions play a significant role in the phase behaviours or structures of ethyl acetate + acetic acid mixture.     

Preparation of sulfoacetate derivatives of cellulose by direct esterification

A new and efficient method for the preparation in one step of water-soluble cellulose acetate sulfate derivatives (CAS) is reported. Acetylation and sulfation were carried out simultaneously, using a mixture of acetic anhydride and sulfuric acid in glacial acetic acid. The reaction time and the amount of acetic anhydride were optimized and the method provided water-soluble esters, with a degree of acetylation in the range 1.6 and 2.4 and a degree of sulfation of 0.3. This method has been successfully applied to pure cellulose and to cellulose-enriched materials obtained from agricultural by-products. The product exhibited a high viscosity in aqueous solution suggesting interesting rheological properties.     

Parameterization and optimization of the menthol force field for molecular dynamics simulations

Menthol’s various biological properties render it a useful component for medical and cosmetological applications, while its three centers of asymmetry mean that it can be used in a range of organic reactions. Menthol-substituted ionic liquids (ILs) have been found to exhibit promising antimicrobial and antielectrostatic properties, as well as being useful in organic catalysis and biochemical studies. However, so far, a force field designed and validated specifically for the menthol molecule has not been constructed. In the present work, the validation and optimization of force field parameters with regard to the ability to reproduce the macroscopic properties of menthol is presented. The set of optimized potentials for liquid simulations all atom (OPLS-AA) compatible parameters was tested and carefully tuned. The refinement of parameters included fitting of partial atomic charges, optimization of Lennard-Jones parameters, and recalculation of the dihedral angle parameters needed to reproduce quantum energy profiles. To validate the force field, a variety of physicochemical properties were calculated for liquid menthol. Both thermodynamic and kinetic properties were taken into account, including density, surface tension, enthalpy of vaporization, and shear viscosity. The obtained force field was proven to accurately reproduce the properties of the investigated compound while being fully compatible with the OPLS-AA force field.

     

MolMod – an open access database of force fields for molecular simulations of fluids

The MolMod database is presented, which is openly accessible at http://molmod.boltzmann-zuse.de and contains intermolecular force fields for over 150 pure fluids at present. It was developed and is maintained by the Boltzmann-Zuse Society for Computational Molecular Engineering (BZS). The set of molecular models in the MolMod database provides a coherent framework for molecular simulations of fluids. The molecular models in the MolMod database consist of Lennard-Jones interaction sites, point charges, and point dipoles and quadrupoles, which can be equivalently represented by multiple point charges. The force fields can be exported as input files for the simulation programmes ms2 and ls1 mardyn, GROMACS, and LAMMPS. To characterise the semantics associated with the numerical database content, a force field nomenclature is introduced that can also be used in other contexts in materials modelling at the atomistic and mesoscopic levels. The models of the pure substances that are included in the database were generally optimised such as to yield good representations of experimental data of the vapour–liquid equilibrium with a focus on the vapour pressure and the saturated liquid density. In many cases, the models also yield good predictions of caloric, transport, and interfacial properties of the pure fluids. For all models, references to the original works in which they were developed are provided. The models can be used straightforwardly for predictions of properties of fluid mixtures using established combination rules. Input errors are a major source of errors in simulations. The MolMod database contributes to reducing such errors.     

Effects of acetylation in vapor phase and mercerization on the properties of sugarcane fibers

Chemical modification of sugarcane bagasse fiber was achieved by mercerization reaction and esterification reaction with anhydride acetic vapor. This is a new acetylation procedure. The results show that the fiber length and diameter are reduced after the reactions. Fourier transform infrared spectroscopy (FT-IR) studies produced clear evidence of the partial acetylation reaction. Optical microscopy revealed fibrillation in the acetylated fiber attributed to hemicellulose dissolution. The thermal stability measured by thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) increased after acetylation and decreased after mercerization. The higher thermal stability of the acetylated fiber as compared with modified fibers in liquid medium was attributed to the small quantity of water and acetic acid present for the reaction in vapor phase. The lesser tensile strength of the acetylated fiber was due to fibrillation. The porous structure obtained favors migration of the polymer chains into the fiber acetylated, and thus it should enhance the polymer–fiber adhesion in polymer composites.     

Accurate representation of B-DNA double helical structure with implicit solvent and counterions

High-resolution nuclear magnetic resonance (NMR) and crystallographic data have been taken to refine the force field used in the torsion angle space nucleic acids molecular mechanics program DUPLEX. The population balance deduced from NMR studies of two carcinogen-modified DNA conformers in equilibrium was used to fine tune a sigmoidal, distance-dependent dielectric function so that reasonable relative energies could be obtained. In addition, the base-pair and backbone geometry from high-resolution crystal structures of the Dickerson-Drew dodecamer was used to re-evaluate the deoxyribose pseudorotation profile and the Lennard-Jones nonbonded energy terms. With a modified dielectric function that assumes a very steep distance-dependent form, a deoxyribose pseudorotation profile with reduced energy barriers between C2'- and C3'-endo minima, and a shift of the Lennard-Jones potential energy minimum to a distance approximately 0.4 A greater than the sum of the van der Waals' radii, the sequence-dependent conformational features of the Dickerson-Drew dodecamer in both the solid state and the aqueous liquid crystalline phase are well reproduced. The robust performance of the revised force field, in conjunction with its efficiency through implicit treatment of solvent and counterions, provides a valuable tool for elucidating conformations and structure-function relationships of DNA, including those of molecules modified by carcinogens and other ligands.     

Comparison of simple perturbation-theory estimates for the liquid–solid and the liquid–vapor interfacial free energies of Lennard-Jones systems

The most naive perturbation method to estimate interfacial free energies is based on the assumption that the interface between coexisting phases is infinitely sharp. Although this approximation does not yield particularly accurate estimates for the liquid–vapor surface tension, we find that it works surprisingly well for the interface between a dense liquid and a solid. As an illustration we estimate the liquid–solid interfacial free energy of a Lennard-Jones system with truncated and shifted interactions and compare the results with numerical data that have been reported in the literature. We find that the agreement between theory and simulation is excellent. In contrast, if we apply the same procedure to estimate the variation of the liquid–vapor surface tension, for different variants of the Lennard-Jones potential (truncated/shifted/force-shifted), we find that the agreement with the available simulation data is, at best, fair. The present method makes it possible to obtain quick and easy estimate of the effect on the surface free energy of different potential-truncation schemes used in computer simulations.     

The investigation of effect of organic carbon sources addition in anaerobic–aerobic (low dissolved oxygen) sequencing batch reactor for nutrients removal from wastewaters

The effect of addition of organic carbon sources (acetic acid and waste activated sludge alkaline fermentation liquid) on anaerobic–aerobic (low dissolved oxygen, 0.15–0.45 mg/L) biological municipal wastewater treatment was investigated. The results showed that carbon source addition affected not only the transformations of polyhydroxyalkanoates (PHA), glycogen, nitrogen and phosphorus, but the net removal of nitrogen and phosphorus. The removal efficiencies of TN and TP were, respectively, 61% and 61% without organic carbon source addition, 81% and 95% with acetic acid addition, and 83% and 97% with waste activated sludge alkaline fermentation liquid addition. It seems that the alkaline fermentation liquid of waste biosolids generated in biological wastewater treatment plant can be used to replace acetic acid as an additional carbon source to improve the anaerobic–aerobic (low dissolved oxygen) municipal wastewater nutrients removal although its use was observed to cause a slight increase of effluent BOD and COD concentrations.     

Process integration involving lipase-catalyzed ester synthesis reactions

Summary The technical feasibility of integrating bioconversion (using a lipase from Mucor miehei immobilized by ion exchange onto a macroporous resin) and product recovery (via mass transfer towards liquid/vapor equilibrium) was experimentally verified for the production of ethyl acetate from acetic acid and ethanol at 30 and 40 °C. The configuration selected consisted of a packed bed thermostatted reactor inserted between the condenser and the liquid reflux to a distillation column. The concentrations of the aforementioned components and water were measured periodically for a 9 to 10-h period. The extents of conversion of acetic acid ranged from 1 to 32% at 30 °C, and from 13 to 25% at 40 °C, and they were maximum at initial molar fractions of acid of ca. 0.3 and 0.4, respectively. The ratio of the reaction product to the equilibrium reaction product increased from 2 to 32 at 30 °C, and from 5 to 27 at 40 °C when the initial molar fraction of acid in the feed increased from 0.2 to 0.7.     

Non Equilibrium Molecular Dynamics Simulation of Coupled Heat- and Mass Transport Across a Liquid/Vapor Interface

Abstract

The transport properties of bulk liquid, gas and at the gas/liquid interface were studied for two binary Lennard-Jones/spline mixtures by use of nonequilibrium molecular dynamics. One of the mixtures was an ideal isotope mixture, the other a non-ideal mixture. The simulations gave the thermal conductivity, mutual diffusion coefficient, heat flux, mass flux, and the changes in these quantities across the interface. The local entropy production was expressed in terms of fluxes and thermodynamic forces, and numercial estimates are given. It was shown that the largest contribution to the total entropy production occurs in the vapor phase under the chosen conditions. We expect, however that if the mass flux were larger, the major contribution to the entropy production would come from the liquid phase.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Water activity coefficients in aqueous amino acid solutions by molecular dynamics simulation: 1. Force field development

New force fields for molecular dynamics (MD) simulation of aqueous zwitterionic amino acid simulations were developed. These were especially designed to calculate activity coefficient of water in amino acid solutions with high accuracy. For example, aqueous solutions of the following amino acids were considered: glycine, alanine, α-aminobutyric acid, α-aminovalerianic acid, valine and leucine. The force fields were obtained by quantum chemical calculations using B3LYP/6-31G and MP2/6-311(d,p) model theories in combination with the Merz–Kollmann–Singh scheme. To further increase the accuracy of the force field, a polarised continuum was considered in all quantum chemical calculations. Water activity coefficients obtained from MD using different all-purpose literature force fields, namely, OPLS, AMBER ff03 and GROMOS 53A6 as well as experimental data are compared with the results utilising the new force field. The new force field is shown to give better results compared with experimental data than existing force fields.     

Calculation of the surface tension of planar interfaces by molecular simulations: from Lennard-Jones fluids to binary mixtures

The relative longevity of the research in the field of the molecular simulations of the liquid–vapour interfaces of Lennard-Jones (LJ) particles can be explained by the dependence of the surface tension on many methodological factors. After a few illustrations on the parameters that can impact the results of surface tension on the LJ interfaces, we establish the ability of the current methodologies to quantitatively predict the surface tension of various liquid–vapour interfaces of pure components at different temperatures. We also show that the methods perform very well for the reproduction of the interfacial tension of binary mixtures in a wide range of pressures.     

Gibbs Ensemble Calculations with an Equation of State: An Application to Vapor-Liquid Equilibria

Abstract

A new modification of the Gibbs ensemble Monte Carlo computer simulation method for fluid phase equilibria is described. The modification is based on a thermodynamic model for the vapor phase, and uses an equation of state to account for the weak interactions between the vapor phase molecules. Reductions in the computational time by 30–40% as compared to the original Gibbs ensemble method are obtained. The algorithm is applied to Lennard-Jones - (12,6) fluids and their mixtures and the results are in good agreement with results obtained from simulations using the full Gibbs ensemble method.     

Analysis of solvent-free ethyl oleate enzymatic synthesis at equilibrium conditions

This work reports experimental equilibrium data for the esterification of pure oleic acid and a fatty acid mixture with ethanol, using an immobilized Candida antarctica B lipase as catalyst. Reactions are performed in a solvent-free system, containing a mixture of substrates and different amounts of distilled water. According to the initial amount of water and the extent of the reaction, one or two liquid phases are present. Therefore, when the equilibrium is achieved, the liquid–liquid and chemical reaction equilibria have to be simultaneously satisfied.

Several reports dealing with enzymatic reactions performed in two-phase systems have found that the value of the reaction equilibrium constant calculated from overall experimental concentrations varies not only with temperature but also with substrate ratio and water content. Although this approach is a valuable way to explore equilibrium shifts in biphasic systems, it is limited to ideal systems with constant partition coefficients. The aim of this work is to consider the biphasic nature of the reactive mixture through a computational procedure that simultaneously takes into account liquid–liquid and reaction equilibria. This approach enables the determination of a classical temperature-dependent thermodynamic equilibrium constant, which accurately fits experimental equilibrium conversions over a wide range of operating conditions.     

Direct Sampling of Local Density Fluctuations in Monte Carlo Simulations

A novel method for accelerating Monte Carlo simulations of fluids based on a direct sampling of local density fluctuations by a multiparticle move is proposed. The method is expected to be particularly efficient for inhomogeneous pure fluids consisting of spherical or moderately nonspherical molecules which is confirmed by a sample simulation of a Lennard-Jones fluid in a slit pore. An analogous method for a mixture, a direct sampling of local concentration fluctuations by swapping particles of different species, is successfully tested on a liquid mixture of argon and nitrogen.     

Enzymatic synthesis of unsaturated fatty acid glucoside esters for dermo-cosmetic applications.

Unsaturated fatty acid alpha-butylglucoside esters were prepared by enzymatic esterification of alpha-butylglucoside in nonaqueous media. Conditions were firstly optimized using oleic acid as acyl group. Synthesis was possible in several solvents but the presence of water co-product in the medium limited the reaction to a thermodynamic equilibrium corresponding to a maximal conversion yield of 62%. In pure molten substrates, the removal of water under reduced pressure enabled yields superior to 95% to be obtained. Product profiles depended on enzyme origin : whatever the support, immobilized lipase B from Candida antarctica proved to be far more regioselective for the primary hydroxyl group of glucose than immobilized lipase from Rhizomucor miehei. Results obtained could be easily transposed to the acylation of alpha-butylglucoside with a commercial mixture of unsaturated fatty acids containing more than 60% of linoleic acid. The biocatalyst could be recycled more than ten times without any significant activity loss.     

4-O-α-d-galactopyranosyl-d-galactose: Efficient synthetic routes from “polygalacturonic acid”

Enzymic hydrolysis of “polygalacturonic acid” gave a mixture of oligomers which was fractionated by ion-exchange chromatography. The resulting di- and tri-saccharides were treated, respectively, with methanol and ethylene oxide, and the resulting esters were reduced with sodium borohydride. Treatment of the products with acetic anhydride and sulfuric acid, followed by deacetylation, produced the title compound.     

The Extrapolation of Vapour–liquid Equilibrium Curves of Pure Fluids in Alternative Gibbs Ensemble Monte Carlo Implementations

Extrapolation schemes based on Taylor series expansion to determine the vapour–liquid equilibrium (VLE) curves of pure molecular fluids are presented for the NpH and μVL versions of the Gibbs ensemble Monte Carlo (GEMC) simulations. The coexistence curves of the various configurational quantities can be expressed as Taylor series around the simulated equilibrium point as a function of pressure in the NpH version and chemical potential in the μVL version. The coefficients of the Taylor series are calculated from single GEMC simulations using Clapeyron-like equations and fluctuation formulas. A Padè approximant is used to widen the range where the extrapolation is accurate. These methods are demonstrated on atomic Lennard-Jones fluid. The procedure is found to be an accurate and useful tool to calculate wide sections of the VLE curves. With this procedure the saturation heat capacity can be directly determined using the calculated derivatives.     

用生料制备功能性双歧醋生产工艺的实验室研究

目的将双歧杆菌、醋酸菌、酵母菌和粉碎的制醋原料及麸曲共同发酵,通过生料制醋的方法来制备功能性双歧醋。方法将粉碎的玉米与麸曲、酵母液、麸皮和水搅拌均匀,使其经过液态糖化和酒精发酵后,接入醋酸菌和双歧杆菌(二者比例为1∶1),同时加入辅料,进行醋酸发酵,当检测到醋酸酸度为5.0%~7.5%时,加入食盐终止发酵,经过过滤,除菌澄清得到功能性双歧醋。结果双歧醋的最终醋酸度为3.2%,外观红棕色,光泽度好,清澈透明,无沉淀和悬浮物。总菌数:醋酸菌为3.3×1011/m l,双歧杆菌为1.9×107/m l;活菌数:醋酸菌为1.7×1011/m l,双歧杆菌为6.8×106/m l;大肠菌群数3个/100 m l;致病菌:不得检出。结论双歧杆菌及其代谢物可以在双歧醋中存活,生料固态发酵制备双歧醋的方法可行。