Pitfalls in the Evaluation of the Thermodynamic Consistency of Experimental VLE Data Sets

The thermodynamic consistency of almost 90 VLE data series, including isothermal and isobaric conditions for systems of both total and partial miscibility in the liquid phase, has been examined by means of the area and point-to-point tests. In addition, the Gibbs energy of mixing function calculated from these experimental data has been inspected, with some rather surprising results: certain data sets exhibiting high dispersion or leading to Gibbs energy of mixing curves inconsistent with the total or partial miscibility of the liquid phase, surprisingly, pass the tests. Several possible inconsistencies in the tests themselves or in their application are discussed. Related to this is a very interesting and ambitious initiative that arose within the NIST organization: the development of an algorithm to assess the quality of experimental VLE data. The present paper questions the applicability of two of the five tests that are combined in the algorithm. It further shows that the deviation of the experimental VLE data from the correlation obtained by a given model, the basis of some point-to-point tests, should not be used to evaluate the quality of these data.

Improving some thermodynamic consistency tests for VLE data evaluation

Thermodynamics Conference 2013 (Statistical Mechanics and Thermodynamics Group of the Royal Society of Chemistry), The University of Manchester, 3-6 September 2013.

Thermodynamic concepts in adaptive resolution simulations

This thesis work is devoted to the conceptual and technical development of the Adaptive Resolution Scheme (AdResS), a molecular dynamics method that allows the simulation of a system with different levels of resolution simultaneously. The simulation domain is divided into high and low resolution zones and a transition region that links them, through which molecules can freely diffuse.rnThe first issue of this work regards the thermodynamic consistency of the method, which is tested and verified in a model liquid of tetrahedral molecules. The results allow the introduction of the concept of the Thermodynamic Force, an external field able to correct spurious density fluctuations present in the transition region in usual AdResS simulations.rnThe AdResS is also applied to a system where two different representations with the same degree of resolution are confronted. This simple test extends the method from an Adaptive Resolution Scheme to an Adaptive Representation Scheme, providing a way of coupling different force fields based on thermodynamic consistency arguments. The Thermodynamic Force is successfully applied to the example described in this work as well.rnAn alternative approach of deducing the Thermodynamic Force from pressure consistency considerations allows the interpretation of AdResS as a first step towards a molecular dynamics simulation in the Grand Canonical ensemble. Additionally, such a definition leads to a practical way of determining the Thermodynamic Force, tested in the well studied tetrahedral liquid. The effects of AdResS and this correction on the atomistic domain are analyzed by inspecting the local distribution of velocities, radial distribution functions, pressure and particle number fluctuation. Their comparison with analogous results coming from purely atomistic simulations shows good agreement, which is greatly improved under the effect of the external field.rnA further step in the development of AdResS, necessary for several applications in biophysics and material science, consists of its application to multicomponent systems. To this aim, the high-resolution representation of a model binary mixture is confronted with its coarse-grained representation systematically parametrized. The Thermodynamic Force, whose development requires a more delicate treatment, also gives satisfactory results.rnFinally, AdResS is tested in systems including two-body bonded forces, through the simulation of a model polymer allowed to adaptively change its representation. It is shown that the distribution functions that characterize the polymer structure are in practice not affected by the change of resolution.rnThe technical details of the implementation of AdResS in the ESPResSo package conclude this thesis work.

Estudo do equilíbrio líquido-vapor do sistema água+etanol+líquido iônico visando a separação do álcool anidro

Anhydrous ethanol is used in chemical, pharmaceutical and fuel industries. However, current processes for obtaining it involve high cost, high energy demand and use of toxic and pollutant solvents. This problem occurs due to the formation of an azeotropic mixture of ethanol + water, which does not allow the complete separation by conventional methods such as simple distillation. As an alternative to currently used processes, this study proposes the use of ionic liquids as solvents in extractive distillation. These are organic salts which are liquids at low temperatures (under 373,15 K). They exhibit characteristics such as low volatility (almost zero/ low vapor ), thermal stability and low corrosiveness, which make them interesting for applications such as catalysts and as entrainers. In this work, experimental data for the vapor pressure of pure ethanol and water in the pressure range of 20 to 101 kPa were obtained as well as for vapor-liquid equilibrium (VLE) of the system ethanol + water at atmospheric pressure; and equilibrium data of ethanol + water + 2-HDEAA (2- hydroxydiethanolamine acetate) at strategic points in the diagram. The device used for these experiments was the Fischer ebulliometer, together with density measurements to determine phase compositions. The experimental data were consistent with literature data and presented thermodynamic consistency, thus the methodology was properly validated. The results were favorable, with the increase of ethanol concentration in the vapor phase, but the increase was not shown to be pronounced. The predictive model COSMO-SAC (COnductor-like Screening MOdels Segment Activity Coefficient) proposed by Lin & Sandler (2002) was studied for calculations to predict vapor-liquid equilibrium of systems ethanol + water + ionic liquids at atmospheric pressure. This is an alternative for predicting phase equilibrium, especially for substances of recent interest, such as ionic liquids. This is so because no experimental data nor any parameters of functional groups (as in the UNIFAC method) are needed

Equilíbrio líquido-líquido para o sistema biodiesel de óleo de semente de melão (Cucumis melo l.) + metanol + glicerina

Rio Grande do Norte, northeast state from Brazil, it is the greatest producer and exporter of yellow melon, well known as Spanish melon. Despite the consumption of this fruit to be mainly its pulp, melon seeds are an important source of lipids considered an industrial residue it has been discharge product. The use of oilseeds in order to produce biodiesel establishes an important raw material and the increase of its production promotes the national development of the agriculture. In this background, the aim of this work has been to use oil from seeds of yellow melon to produce biodiesel and to accomplish a study of the phase equilibrium of the system evolving biodiesel, methanol and glycerin. The biodiesel was obtained by oil transesterification through methylic route with molar ratio 1:9.7 (oil:alcohol) and with a mass of NaOH of 0.5% from the oil mass; the reaction time was 73 minutes at 55 °C. A yield of 84.94% in biodiesel was achieved. The equilibria data present a well-characterized behavior with a great region of two phases. The tie lines indicate that methanol has a best solubility in the phase that is rich in glycerin. Consistency of the experimental data was made based on Othmer-Tobias and Hand correlations which values above 0.99 were found to correlation coefficients, this fact confers a good thermodynamic consistency to the experimental data. NRTL and UNIQUAC models were employed to predict liquid-liquid equilibrium of this system. It was observed a better concordance of the results when NRTL was applied (standard deviation 1.25%) although the UNIQUAC model has presented a quite satisfactory result either (standard deviation 2.70%). The NRTL and UNIQUAC models were also used to evaluate the effect of temperature in the range of 328 K to 358 K, in which a little change in solubility with respect to the data obtained at 298 K was observed, thus being considered negligible the effect of temperature

Determinação de dados de equilíbrio de fases para sistemas aquosos com eletrólitos

In the present work are established initially the fundamental relationships of thermodynamics that govern the equilibrium between phases, the models used for the description of the behavior non ideal of the liquid and vapor phases in conditions of low pressures. This work seeks the determination of vapor-liquid equilibrium (VLE) data for a series of multicomponents mixtures of saturated aliphatic hydrocarbons, prepared synthetically starting from substances with analytical degree and the development of a new dynamic cell with circulation of the vapor phase. The apparatus and experimental procedures developed are described and applied for the determination of VLE data. VLE isobarics data were obtained through a Fischer s ebulliometer of circulation of both phases, for the systems pentane + dodecane, heptane + dodecane and decane + dodecane. Using the two new dynamic cells especially projected, of easy operation and low cost, with circulation of the vapor phase, data for the systems heptane + decane + dodecane, acetone + water, tween 20 + dodecane, phenol + water and distillation curves of a gasoline without addictive were measured. Compositions of the equilibrium phases were found by densimetry, chromatography, and total organic carbon analyzer. Calibration curves of density versus composition were prepared from synthetic mixtures and the behavior excess volumes were evaluated. The VLE data obtained experimentally for the hydrocarbon and aqueous systems were submitted to the test of thermodynamic consistency, as well as the obtained from the literature data for another binary systems, mainly in the bank DDB (Dortmund Data Bank), where the Gibbs-Duhem equation is used obtaining a satisfactory data base. The results of the thermodynamic consistency tests for the binary and ternary systems were evaluated in terms of deviations for applications such as model development. Later, those groups of data (tested and approved) were used in the KijPoly program for the determination of the binary kij parameters of the cubic equations of state original Peng-Robinson and with the expanded alpha function. These obtained parameters can be applied for simulation of the reservoirs petroleum conditions and of the several distillation processes found in the petrochemistry industry, through simulators. The two designed dynamic cells used equipments of national technology for the determination of VLE data were well succeed, demonstrating efficiency and low cost. Multicomponents systems, mixtures of components of different molecular weights and also diluted solutions may be studied in these developed VLE cells

Determinação de Dados de Equilíbrio Líquido-Vapor para Sistemas Hidrocarbonetos e Desenvolvimento de uma Nova Célula Dinâmica

In the present work are established initially the fundamental relationships of thermodynamics that govern the equilibrium between phases, the models used for the description of the behavior non ideal of the liquid and vapor phases in conditions of low pressures. This work seeks the determination of vapor-liquid equilibrium (VLE) data for a series of multicomponents mixtures of saturated aliphatic hydrocarbons, prepared synthetically starting from substances with analytical degree and the development of a new dynamic cell with circulation of the vapor phase. The apparatus and experimental procedures developed are described and applied for the determination of VLE data. VLE isobarics data were obtained through a Fischer's ebulliometer of circulation of both phases, for the systems pentane + dodecane, heptane + dodecane and decane + dodecane. Using the two new dynamic cells especially projected, of easy operation and low cost, with circulation of the vapor phase, data for the systems heptane + decane + dodecane, acetone + water, tween 20 + dodecane, phenol + water and distillation curves of a gasoline without addictive were measured. Compositions of the equilibrium phases were found by densimetry, chromatography, and total organic carbon analyzer. Calibration curves of density versus composition were prepared from synthetic mixtures and the behavior excess volumes were evaluated. The VLE data obtained experimentally for the hydrocarbon and aqueous systems were submitted to the test of thermodynamic consistency, as well as the obtained from the literature data for another binary systems, mainly in the bank DDB (Dortmund Data Bank), where the Gibbs-Duhem equation is used obtaining a satisfactory data base. The results of the thermodynamic consistency tests for the binary and ternary systems were evaluated in terms of deviations for applications such as model development. Later, those groups of data (tested and approved) were used in the KijPoly program for the determination of the binary kij parameters of the cubic equations of state original Peng-Robinson and with the expanded alpha function. These obtained parameters can be applied for simulation of the reservoirs petroleum conditions and of the several distillation processes found in the petrochemistry industry, through simulators. The two designed dynamic cells used equipments of national technology for the determination Humberto Neves Maia de Oliveira Tese de Doutorado PPGEQ/PRH-ANP 14/UFRN of VLE data were well succeed, demonstrating efficiency and low cost. Multicomponents systems, mixtures of components of different molecular weights and also diluted solutions may be studied in these developed VLE cells

Liquid–Liquid, Vapor–Liquid, and Vapor–Liquid–Liquid Equilibrium Data for the Water–n-Butanol–Cyclohexane System at Atmospheric Pressure: Experimental Determination and Correlation

The temperature and the composition of the vapor–liquid–liquid equilibrium (VLLE) and the vapor–liquid equilibrium (VLE) of a ternary mixture of water–n-butanol–cyclohexane were measured at atmospheric pressure (101.32 kPa) in a modified dynamic recirculating still. As found in the literature, the experimental data obtained reveal a ternary azeotrope at 341.86 K with a mole fraction composition of 0.281, 0.034, and 0.685 water, n-butanol, and cyclohexane, respectively. The liquid–liquid equilibrium (LLE) compositions were measured at a constant temperature of 313.15 K and compared with data in the literature collected at other temperatures. Thermodynamic consistency of all the experimental data was demonstrated. The universal quasichemical (UNIQUAC) and the nonrandom two-liquid (NRTL) thermodynamic models were used to correlate the VLE and LLE data, while the original universal functional (UNIFAC) model was used to compare the predicted data.

Phase equilibria of the water + 1-butanol + toluene ternary system at 101.3 kPa

Isobaric vapour–liquid and vapour–liquid–liquid equilibrium data for the water + 1-butanol + toluene ternary system were measured at 101.3 kPa with a modified VLE 602 Fischer apparatus. In addition, the liquid–liquid equilibrium data at 313.15 K were measured and compared with data from other authors at different temperatures. The system exhibits a ternary heterogeneous azeotrope whose temperature and composition have been determined by interpolation. The thermodynamic consistency of the experimental vapour–liquid and vapour–liquid–liquid data was checked by means of the Wisniak’s Li/Wi consistency test. Moreover, the vapour–liquid and the liquid–liquid equilibrium correlation for the ternary system with NRTL and UNIQUAC models, together with the prediction made with the UNIFAC model, were studied and discussed.

Vapour-liquid equilibrium studies at atmospheric to moderate pressures

A study of vapour-liquid equilibria is presented together with current developments. The theory of vapour-liquid equilibria is discussed. Both experimental and prediction methods for obtaining vapour-liquid equilibria data are critically reviewed. The development of a new family of equilibrium stills to measure experimental VLE data from sub-atmosphere to 35 bar pressure is described. Existing experimental techniques are reviewed, to highlight the needs for these new apparati and their major attributes. Details are provided of how apparatus may be further improved and how computer control may be implemented. To provide a rigorous test of the apparatus the stills have been commissioned using acetic acid-water mixture at one atmosphere pressure. A Barker-type consistency test computer program, which allows for association in both phases has been applied to the data generated and clearly shows that the stills produce data of a very high quality. Two high quality data sets, for the mixture acetone-chloroform, have been generated at one atmosphere and 64.3oC. These data are used to investigate the ability of the new novel technique, based on molecular parameters, to predict VLE data for highly polar mixtures. Eight, vapour-liquid equilibrium data sets have been produced for the cyclohexane-ethanol mixture at one atmosphere, 2, 4, 6, 8 and 11 bar, 90.9oC and 132.8oC. These data sets have been tested for thermodynamic consistency using a Barker-type fitting package and shown to be of high quality. The data have been used to investigate the dependence of UNIQUAC parameters with temperature. The data have in addition been used to compare directly the performance of the predictive methods - Original UNIFAC, a modified version of UNIFAC, and the new novel technique, based on molecular parameters developed from generalised London's potential (GLP) theory.

Vapor–liquid equilibrium of binary mixtures containing isopropyl acetate and alkanols at 101.32 kPa

Isobaric vapor-liquid equilibria of binary mixtures of isopropyl acetate plus an alkanol (1-propanol, 2-propanol, 1-butanol, or 2-butanol) were measured at 101.32 kPa, using a dynamic recirculating still. An azeotropic behavior was observed only in the mixtures of isopropyl acetate + 2-propanol and isopropyl acetate + 1-propanol. The application of four thermodynamic consistency tests (the Herington test, the Van Ness test, the infinite dilution test, and the pure component test) showed the high quality of the experimental data. Finally, both NRTL and UNIQUAC activity coefficient models were successfully applied in the correlation of the measured data, with the average absolute deviations in vapor phase composition and temperature of 0.01 and 0.16 K, respectively.

Thermodynamic self-consistency issues related to the Cluster Variation Method: The case of the BCC Cr-Fe (Chromium-Iron) system

The Cluster Variation Method (CVM), introduced over 50 years ago by Prof. Dr. Ryoichi Kikuchi, is applied to the thermodynamic modeling of the BCC Cr-Fe system in the irregular tetrahedron approximation, using experimental thermochemical data as initial input for accessing the model parameters. The results are checked against independent data on the low-temperature miscibility gap, using increasingly accurate thermodynamic models, first by the inclusion of the magnetic degrees of freedom of iron and then also by the inclusion of the magnetic degrees of freedom of chromium. It is shown that a reasonably accurate description of the phase diagram at the iron-rich side (i.e. the miscibility gap borders and the Curie line) is obtained, but only at expense of the agreement with the above mentioned thermochemical data. Reasons for these inconsistencies are discussed, especially with regard to the need of introducing vibrational degrees of freedom in the CVM model. (C) 2008 Elsevier Ltd. All rights reserved.

Thermodynamic Study Of The Micellization Of Zwitterionic Surfactants And Their Interaction With Polymers In Water By Isothermal Titration Calorimetry.

The micellization of a homologous series of zwitterionic surfactants, a group of sulfobetaines, was studied using isothermal titration calorimetry (ITC) in the temperature range from 15 to 65 °C. The increase in both temperature and the alkyl chain length leads to more negative values of ΔGmic(0) , favoring the micellization. The entropic term (ΔSmic(0)) is predominant at lower temperatures, and above ca. 55-65 °C, the enthalpic term (ΔHmic(0)) becomes prevalent, figuring a jointly driven process as the temperature increases. The interaction of these sulfobetaines with different polymers was also studied by ITC. Among the polymers studied, only two induced the formation of micellar aggregates at lower surfactant concentration: poly(acrylic acid), PAA, probably due to the formation of hydrogen bonds between the carboxylic group of the polymer and the sulfonate group of the surfactant, and poly(sodium 4-styrenesulfonate), PSS, probably due to the incorporation of the hydrophobic styrene group into the micelles. The prevalence of the hydrophobic and not the electrostatic contributions to the interaction between sulfobetaine and PSS was confirmed by an increased interaction enthalpy in the presence of electrolytes (NaCl) and by the observation of a significant temperature dependence, the latter consistent with the proposed removal of hydrophobic groups from water.

On The Consistency Of Monte Carlo Track Structure Dna Damage Simulations.

Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV-220 keV were superimposed on a geometric DNA model composed of 2.7 × 10(6) nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy ETH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when ETH ≈ 10.79 eV, but deviate significantly for higher ETH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors' show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.

Anomalous thermodynamic properties in ferropericlase throughout its spin crossover transition

The thermodynamics properties of ferropericlase (Mg(1-x)Fe(x)O where x=0.1875) (Fp) throughout its spin crossover were investigated by first principles. Fp was treated as an ideal solid solution of pure high-spin and low-spin states. The Gibbs free energies of the pure states were addressed using the LDA+U method. A vibrational virtual-crystal model was developed to address the vibrational properties of the pure spin cases and used in conjunction with quasiharmonic theory to compute their vibrational free energies. The thermodynamics properties of Fp display significant anomalies that should be typical of spin crossover systems in general. In Fp, in particular, they are fundamental for understanding the state of earth's interior, where the pressure and temperature conditions of the crossover are realized.