Theoretical and Computational Analysis of Static and Dynamic Anomalies in Water-DMSO Binary Mixture at Low DMSO Concentrations

Experiments have repeatedly observed both thermodynamic and dynamic anomalies in aqueous binary mixtures, surprisingly at low solute concentration. Examples of such binary mixtures include water-DMSO, water-ethanol, water-tertiary butyl alcohol (TBA), and water-dioxane, to name a few. The anomalies have often been attributed to the onset of a structural transition, whose nature, however, has been left rather unclear. Here we study the origin of such anomalies using large scale computer simulations and theoretical analysis in water-DMSO binary mixture. At very low DMSO concentration (below 10%), small aggregates of DMSO are solvated by water through the formation of DMSO-(H2O)(2) moieties. As the concentration is increased beyond 10-12% of DMSO, spanning clusters comprising the same moieties appear in the system. Those clusters are formed and stabilized not only through H-bonding but also through the association of CH3 groups of DMSO. We attribute the experimentally observed anomalies to a continuum percolation-like transition at DMSO concentration X-DMSO approximate to 12-15%. The largest cluster size of CH3-CH3 aggregation clearly indicates the formation of such percolating clusters. As a result, a significant slowing down is observed in the decay of associated rotational auto time correlation functions (of the S = O bond vector of DMSO and O-H bond vector of water). Markedly unusual behavior in the mean square fluctuation of total dipole moment again suggests a structural transition around the same concentration range. Furthermore, we map our findings to an interacting lattice model which substantiates the continuum percolation model as the reason for low concentration anomalies in binary mixtures where the solutes involved have both hydrophilic and hydrophobic moieties.

Applicability of Butler's equation in interpreting the thermodynamic behavior of surfaces and adsorption in Fe-S-O melts

Expressions for various second-order derivatives of surface tension with respect to composition at infinite dilution in terms of the interaction parameters of the surface and those of the bulk phases of dilute ternary melts have been presented. A method of deducing the parameters, which consists of repeated differentiation of Butler's equations with subsequent application of the appropriate boundary conditions, has been developed. The present investigation calculates the surface tension and adsorption functions of the Fe-S-O melts at 1873 and 1923 K using the modified form of Butler's equations and the derived values for the surface interaction parameters of the system. The calculated values are found to be in good agreement with those of the experimental data of the system. The present analysis indicates that the energetics of the surface phase are considerably different from those of the bulk phase. The present research investigates a critical compositional range beyond which the surface tension increases with temperature. The observed increase in adsorption of sulfur with consequent desorption of oxygen as a function of temperature above the critical compositional range has been ascribed to the increase of activity ratios of oxygen to sulfur in the surface relative to those in the bulk phase of the system.

Integral treatment for the representation of thermodynamic properties in multicomponent systems using interaction parameters

The partial thermodynamic functions of the solvent component of a ternary system have been deduced in terms of the interaction parameters by integration of several series which emerge from the Maclaurin infinite series based on the integral property of the system and subjected to appropriate boundary conditions. The series integration shows that the resulting partial functions are suitable for interpreting the thermodynamic properties of the system and are independent of compositional paths. In the present analysis, the higher order terms of these series are found to make insignificant contributions.

Interaction of sheep liver apo-serine hydroxymethyltransferase with pyridoxal-5'-phosphate: A physicochemical, kinetic, and thermodynamic study

Sheep liver serine hydroxymethyltransferase (EC 2.1.2.1) is a homotetramer of M(r) 213,000 requiring pyridoxal-5'-phosphate (PLP) as cofactor, Removal of PLP from the holoenzyme converted the enzyme to the apo form which, in addition to being inactive, was devoid of the characteristic absorption spectrum. Upon the addition of PLP to the apoenzyme, complete activity was restored and the visible absorption spectrum with a maximum at 425 nm was regained. The interaction of PLP with the apoenzyme revealed two phases of reaction with pseudo-first-order rate constants of 20 +/- 5 s(-1) and 12.2 +/- 2.0 x 10(-3) s(-1), respectively. However, addition of PLP to the apoenzyme did not cause gross conformational changes as evidenced by circular dichroic and fluorescence spectroscopy. Although conformationally apoenzyme and holoenzyme were indistinguishable, they had distinct apparent melting temperatures of 51 +/- 2 and 58 +/- 2 degrees C, respectively, and the reconstituted holoenzyme was thermally as stable as the native holoenzyme. These results suggested that there was no apparent difference in the secondary structure of holoenzyme, apoenzyme, and reconstituted holoenzyme, However, sedimentation analysis of the apoenzyme revealed the presence of two peaks of S-20,S-w values of 8.7 +/- 0.5 and 5.7 +/- 0.3 S, respectively. A similar pattern was observed when the apoenzyme was chromatographed on a calibrated Sephadex G-150 column. The first peak corresponded to the tetrameric form (M(r) 200,000 +/- 15,000) while the second peak had a M(r) of 130,000 +/- 10,000. Reconstitution experiments revealed that only the tetrameric form of the apoenzyme could be converted into an active holoenzyme while the dimeric form could not be reconstituted into an active enzyme. These results demonstrate that PLP plays an important role in maintaining the structural integrity of the enzyme by preventing the dissociation of the enzyme into subunits, in addition to its function in catalysis. (C) 1996 Academic Press, Inc.

Representation of thermodynamic properties of ternary systems and its application to the system silver-gold-copper at 1350 K

The article presents a generalized analytical expression for description of the integral excess Gibbs free energy of mixing of a ternary system. Twelve constants of the equation are assessed by the least mean squares regressional analysis of the experimental integral excess data of the constituent binaries; three ternary parameters are evaluated by a regressional analysis based on the partial experimental data of a component of the ternary system. The assessed values of the ternary parameters describe the nature of the ternary interaction in the system. Activities and isoactivities of the components in the Ag-Au-Cu system at 1350 K are calculated and found to be in good agreement with the experimental data. This analytical treatment is particularly useful to ternary systems where the thermodynamic data are available from different sources.

Thermodynamic modelling of absorption-resorption heating cycles with some new working pairs

This paper presents the results of a thermodynamic cycle analysis of single stage resorption heat pump (RHP) and resorption heat transformer (RHT) cycles with the new working pairs R22-NMP and R22-DMA. The coefficients of performance (COP) are correlated with the low grade source temperature, temperature at which useful heat is obtained and ambient temperature. The COPs are in the range 1.20–1.60 for the RHP mode and 0.25–0.45 for the RHT mode. Absorber temperatures (useful temperatures) as high as 50°C in the RHP mode and 87°C in the RHT mode have been obtained. It is observed that absorption-resorption systems are inflexible in their range of operating temperature and necessitate a higher pump work as compared with simple single-stage absorption heating systems. However, single stage RHTs show higher temperature boosts than simple absorption heat transformers.

Phase Relations in the System La-Rh-O and Thermodynamic Properties of LaRhO3

Phase relations in the system La-Rh-O at 1223 Ii have been determined by examination of equilibrated samples by optical and scanning electron microscopy, powder X-ray diffraction (XRD), and energy-dispersive analysis of X-rays (EDAX). Only one ternary oxide, LaRhO3, with distorted orthorhombic perovskite structure (Pbnm, a = 0.5525, b = 0.5680, and c = 0.7901 nm) was identified. The alloys and intermetallics along the La-Rh binary are in equilibrium with La2O3. The thermodynamic properties of LaRhO3 were determined in the temperature range 890 to 1310 K, using a solid-state cell incorporating yttria-stabilized zirconia as the electrolyte. A new four-compartment design of the emf cell was used to enhance the accuracy of measurement. For the reaction 1/2La(2)O(3) + 1/2Rh(2)O(3) --> LaRhO3, Delta G degrees = - 70 780 + 4.89T (+/- 90) J.mol(-1) The compound decomposes on heating to a mixture of La2O3, Ph and O-2. The calculated decomposition temperatures are 1843 (+/- 5) K in pure O-2 and 1728 (+/- 5) K in air at a pressure of 1.01 x 10(5) Pa. The phase diagrams for the system La-Rh-O at different partial pressures of oxygen are calculated from the thermodynamic information.

The CaO-SrO-CuO-O-2 system: Phase equilibria and thermodynamic properties at 1123 K

Phase relationships in the CaO-SrO-CuO system in pure oxygen at 1.01 x 10(5) Pa pressure were determined by equilibrating different compositions at 1123 K for similar to 120 h and analyzing the phases present in the quenched samples using X-ray diffraction (XRD), optical and scanning electron microscopy, and energy dispersive analysis of X-rays (EDAX). Four solid solution series were observed in the system, The CawSr1-wO monoxide solid solution with rock-salt structure was found to exhibit an asymmetric miscibility gap, The mixing properties of the monoxide system were deduced using a subregular solution model, For the (CaxSr1-x)(2)CuO3 series, a complete solid solution range with orthorhombic space group Immm was obtained. Calcium substituted for strontium up to 68 at. % in SrCuO2+delta and 51.5 at. % in Sr14Cu24O41-delta. The tie lines between the solid solutions were determined accurately, The activity-composition relations in (CaxSr1-x)(2)CuO3, CaySr1-yCuO2+delta, and (Ca2Sr1-z)(14)Cu24O41-delta solid solutions were determined from experimental tie lines. Activities in the (CaxSr1-x)(2)CuO3 and CaySr1-yCuO2+delta series were close to the predictions of the Temkin model, The behavior of the (CazSr1-(z))(14)Cu24O41-delta solid solution was more complex, with the activity of SrCu(24/14)O-(41-delta/14) exhibiting both positive and negative deviations from ideality. Gibbs energy of formation of the CaCuO2+delta metastable phase at 1123 K was deduced from an analysis of the phase diagram.

Theoretical Analysis of the Electromotive Force of a Cell Incorporating a Composition Gradient Solid Electrolyte

New composition gradient solid electrolytes have been designed for application in high temperature solid-state galvanic sensors and in thermodynamic measurements. The functionally gradient electrolyte consists of a solid solution between two or more ionic conductors with a common ion and gradual variation in composition of the other ionic species. Unequal rates of migration of the ions, caused by the presence of the concentration gradient, may result in the development of space charge, manifesting as diffusion potential. Presented is a theoretical analysis of the EMF of cells incorporating gradient solid electrolytes. An analytical expression is derived for diffusion potential, using the thermodynamics of irreversible processes, for different types of concentration gradients and boundary conditions at the electrode/electrolyte interfaces. The diffusion potential of an isothermal cell incorporating these gradient electrolytes becomes negligible if there is only one mobile ion and the transport numbers of the relatively immobile polyionic species and electrons approach zero. The analysis of the EMF of a nonisothermal cell incorporating a composition gradient solid electrolyte indicates that the cell EMF can be expressed in terms of the thermodynamic parameters at the electrodes and the Seebeck coefficient of the gradient electrolyte under standard conditions when the transport number of one of the ions approaches unity.

The effect of temperature on the surface tension and adsorption functions of the Fe-S-O melts - an analysis based on interaction parameters

The present investigation analyses the thermodynamic behaviour of the surfaces and adsorption as a function of temperature and composition in the Fe-S-O melts based on the Butler's equations. The calculated-values of the surface tensions exhibit an elevation or depression depending on the type of the added solute at a concentration which coincides with that already present in the system. Generally, the desorption of the solutes as a function of temperature results in an initial increase followed by a decrease in the values of the surface tension. The observations are analyzed based on the surface interaction parameters which are derived in the present research.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

Phase diagram for the system Nd-Mn-O and thermodynamic properties of NdMnO3 and NdMn2O5

Studies on the phase relations in the system Nd-Mn-O at 1223 K showed two stable ternary compounds, NdMnO3 and NdMn2O5. An isothermal section of the ternary phase diagram for the system Nd-Mn-O was constructed based on phase analysis of samples quenched after equilibration using XRPD and EDS. An advanced version of the solid-state cell incorporating a buffer electrode was used to determine the Gibbs energies of decomposition of NdMnO3 and NdMn2O5 in the temperature range from 925 to 1400 K. Pure oxygen gas at 0.1 MPa was used as the reference electrode, and yttria-stabilized zirconia as the solid electrolyte. The buffer electrode was designed to prevent polarization of the three-phase electrode and ensure accurate data. The measured oxygen potential corresponding to the reaction,2 Nd2O3 + 4 MnO + O-2 --> 4 NdMnO3 can be represented by the equation,Amu(o2) / J.mol(-1) (+/-580) = -523 960 + 170.96 (T/K)Similarly, for the formation of NdMn2O5 according to the reaction,3 NdMnO3 + Mn3O4 + O-2 --> 3 NdMn2O5 Amu(o2) / J.mol(-1) (+/-660) = - 269 390 + 181.74 (T/K) (C) 2002 Elsevier Science Ltd. All rights reserved.

System Dy-Fe-O: thermodynamic properties of ternary oxides using Calvet calorimetry and solid-state electrochemical cell

The enthalpy increments and the standard molar Gibbs energies of formation-of DyFeO3(s) and Dy3Fe5O12(s) have been measured using a Calvet micro-calorimeter and a solid oxide galvanic cell, respectively. A co-operative phase transition, related to anti-ferromagnetic to paramagnetic transformation, is apparent. from the heat capacity data for DyFeO3 at similar to 648 K. A similar type of phase transition has been observed for Dy3Fe5O12 at similar to 560 K which is related to ferrimagnetic to paramagnetic transformation. Enthalpy increment data for DyFeO3(s) and Dy3Fe5O12(s), except in the vicinity of the second-order transition, can be represented by the following polynomial expressions:{H(0)m(T) - H(0)m(298.15 K)) (Jmol(-1)) (+/-1.1%) = -52754 + 142.9 x (T (K)) + 2.48 x 10(-3) x (T (K))(2) + 2.951 x 10(6) x (T (K))(-1); (298.15 less than or equal to T (K) less than or equal to 1000) for DyFeO3(s), and {H(0)m(T) - H(0)m(298.15 K)} (Jmol(-1)) (+/-1.2%) = -191048 + 545.0 x (T - (K)) + 2.0 x 10(-5) x (T (K))(2) + 8.513 x 10(6) x (T (K))(-1); (208.15 less than or equal to T (K) less than or equal to 1000)for Dy3Fe5O12(s). The reversible emfs of the solid-state electrochemical cells: (-)Pt/{DyFeO3(s) + Dy2O3(s) + Fe(s)}/YDT/CSZ//{Fe(s) + Fe0.95O(s)}/Pt(+) and (-)Pt/{Fe(s) + Fe0.95O(s)}//CSZ//{DyFeO3(s) + Dy3Fe5O12(s) + Fe3O4(s)}/Pt(+), were measured in the temperature range from 1021 to 1250 K and 1035 to 1250 K, respectively. The standard Gibbs energies of formation of solid DyFeO3 and Dy3Fe5O12 calculated by the least squares regression analysis of the data obtained in the present study, and data for Fe0.95O and Dy2O3 from the literature, are given by Delta(f)G(0)m(DyFeO3,s)(kJmol(-1))(+/-3.2)= -1339.9 + 0.2473 x (T(K)); (1021 less than or equal to T (K) less than or equal to 1548)and D(f)G(0)m(Dy3Fe5O12,s) (kJmol(-1)) (+/-3.5) = -4850.4 + 0.9846 x (T (K)); (1035 less than or equal to T (K) less than or equal to 1250) The uncertainty estimates for Delta(f)G(0)m include the standard deviation in the emf and uncertainty in the data taken from the literature. Based on the thermodynamic information, oxygen potential diagram and chemical potential diagrams for the system Dy-Fe-O were developed at 1250 K. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Thermodynamic studies on NdFeO3(s)

The enthalpy increments and the standard molar Gibbs energy of formation of NdFeO3(s) have been measured using a hightemperature Calvet microcalorimeter and a solid oxide galvanic cell, respectively. A lambda-type transition, related to magnetic order-disorder transformation (antiferromagnetic to paramagnetic), is apparent from the heat capacity data at similar to 687 K. Enthalpy increments, except in the vicinity of transition, can be represented by a polynomial expression: {Hdegrees(m)(T)-Hdegrees(m) (298.15 K)} /J(.)mol(-1) (+/- 0.7%)=-53625.6+146.0(T/K) +1.150 X 10(-4)(T/K)(2) +3.007 x 10(6)(T/K)(-1); (298.15 less than or equal to T/K less than or equal to 1000). The heat capacity, the first differential of {Hdegrees(m)(T)-Hdegrees(m)(298.15 K)}with respect to temperature, is given by Cdegrees(pm)/J(.)K(-1.)mol(-1)=146.0+ 2.30x10(-4) (T/K) - 3.007 X 10(6)(T/K)(-2). The reversible emf's of the cell, (-) Pt/{NdFeO3(s) +Nd2O3(s)+Fe(s)}//YDT/CSZ// Fe(s)+'FeO'(s)}/Pt(+), were measured in the temperature range from 1004 to 1208 K. It can be represented within experimental error by a linear equation: E/V=(0.1418 +/- 0.0003)-(3.890 +/- 0.023) x 10(-5)(T/K). The Gibbs energy of formation of solid NdFeO, calculated by the least-squares regression analysis of the data obtained in the present study, and data for Fe0.95O and Nd2O3 from the literature, is given by Delta(f)Gdegrees(m)(NdFeO3 s)/kJ (.) mol(-1)( +/- 2.0)=1345.9+0.2542(T/K); (1000 less than or equal to T/K less than or equal to 1650). The error in Delta(f)Gdegrees(m)(NdFeO3, s, T) includes the standard deviation in emf and the uncertainty in the data taken from the literature. Values of Delta(f)Hdegrees(m)(NdFeO3, s, 298.15 K) and Sdegrees(m) (NdFeO3 s, 298.15 K) calculated by the second law method are - 1362.5 (+/-6) kJ (.) mol(-1) and 123.9 (+/-2.5) J (.) K-1 (.) mol(-1), respectively. Based on the thermodynamic information, an oxygen potential diagram for the system Nd-Fe-O was developed at 1350 K. (C) 2002 Elsevier Science (USA).