Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

An understanding of chemoselective hydrogenation on crotonaldehyde over Pt(111) in the free energy landscape: The microkinetics study based on first-principles calculations

Microkinetic model is developed in the free energy landscape based on density functional theory (DFT) to quantitatively investigate the reaction mechanism of chemoselective partial hydrogenation of crotonaldehyde to crotyl alcohol over Pt(1 1 1) at the temperature of 353 K. Three different methods (mobile, immobile and collision theory models) were carried out to obtain free energy barrier of adsorption/desorption processes. The results from mobile and collision theory models are similar. The calculated TOFs from both models are close to the experiment value. However, for the immobile model, in which the free energy barrier of desorption approaches the energy barrier, the calculated TOF is 2 orders of magnitude lower than the other models. The difficulty of adsorption/ desorption may be overestimated in the immobile model. In addition, detailed analyses show that for the surface hydrogenation elementary steps, the entropy and internal energy effects are small under the reaction condition, while the zero-point-energy (ZPE) correction is significant, especially for the multi-step hydrogenation reaction. The total energy with the ZPE correction approaches to the full free energy calculation for the surface reaction under the reaction condition. (c) 2011 Elsevier B.V. All rights reserved.

Bio-port Free Energy City

In an age of depleting oil reserves and increasing energy demand, humanity faces a stalemate between environmentalism and politics, where crude oil is traded at record highs yet the spotlight on being ‘green’ and sustainable is stronger than ever. A key theme on today’s political agenda is energy independence from foreign nations, and the United Kingdom is bracing itself for nuclear renaissance which is hoped will feed the rapacious centralised system that the UK is structured upon. But what if this centralised system was dissembled, and in its place stood dozens of cities which grow and monopolise from their own energy? Rather than one dominant network, would a series of autonomous city-based energy systems not offer a mutually profitable alternative? Bio-Port is a utopian vision of a ‘Free Energy City’ set in Liverpool, where the old dockyards, redundant space, and the Mersey Estuary have been transformed into bio-productive algae farms. Bio-Port Free Energy City is a utopian ideal, where energy is superfluous; in fact so abundant that meters are obsolete. The city functions as an energy generator and thrives from its own product with minimal impact upon the planet it inhabits. Algaculture is the fundamental energy source, where a matrix of algae reactors swamp the abandoned dockyards; which themselves have been further expanded and reclaimed from the River Mersey. Each year, the algae farm is capable of producing over 200 million gallons of bio-fuel, which in-turn can produce enough electricity to power almost 2 million homes. The metabolism of Free-Energy City is circular and holistic, where the waste products of one process are simply the inputs of a new one. Livestock farming – once traditionally a high-carbon countryside exercise has become urbanised. Cattle are located alongside the algae matrix, and waste gases emitted by farmyards and livestock are largely sequestered by algal blooms or anaerobically converted to natural gas. Bio-Port Free Energy City mitigates the imbalances between ecology and urbanity, and exemplifies an environment where nature and the human machine can function productively and in harmony with one another. According to James Lovelock, our population has grown in number to the point where our presence is perceptibly disabling the planet, but in order to reverse the effects of our humanist flaws, it is vital that new eco-urban utopias are realised.

Exploring the free energy surfaces of clusters using reconnaissance metadynamics

A new approach is proposed for exploring the low-energy structures of small to medium-sized aggregates of atoms and molecules. This approach uses the recently proposed reconnaissance metadynamics method [G. A. Tribello, M. Ceriotti, and M. Parrinello. Proc. Natl. Acad. Sci. U.S.A. 107(41), 17509 (2010)] in tandem with collective variables that describe the average structure of the coordination sphere around the atoms/molecules. We demonstrate this method on both Lennard-Jones and water clusters and show how it is able to quickly find the global minimum in the potential energy surface, while exploring the finite temperature free energy surface. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3628676]

Simplifying the representation of complex free-energy landscapes using sketch-map

A new scheme, sketch-map, for obtaining a low-dimensional representation of the region of phase space explored during an enhanced dynamics simulation is proposed. We show evidence, from an examination of the distribution of pairwise distances between frames, that some features of the free-energy surface are inherently high-dimensional. This makes dimensionality reduction problematic because the data does not satisfy the assumptions made in conventional manifold learning algorithms We therefore propose that when dimensionality reduction is performed on trajectory data one should think of the resultant embedding as a quickly sketched set of directions rather than a road map. In other words, the embedding tells one about the connectivity between states but does not provide the vectors that correspond to the slow degrees of freedom. This realization informs the development of sketch-map, which endeavors to reproduce the proximity information from the high-dimensionality description in a space of lower dimensionality even when a faithful embedding is not possible.

Solid-liquid interfacial free energy out of equilibrium

The properties of the interface between solid and melt are key to solidification and melting, as the interfacial free energy introduces a kinetic barrier to phase transitions. This makes solidification happen below the melting temperature, in out-of-equilibrium conditions at which the interfacial free energy is ill defined. Here we draw a connection between the atomistic description of a diffuse solid-liquid interface and its thermodynamic characterization. This framework resolves the ambiguities in defining the solid-liquid interfacial free energy above and below the melting temperature. In addition, we introduce a simulation protocol that allows solid-liquid interfaces to be reversibly created and destroyed at conditions relevant for experiments. We directly evaluate the value of the interfacial free energy away from the melting point for a simple but realistic atomic potential, and find a more complex temperature dependence than the constant positive slope that has been generally assumed based on phenomenological considerations and that has been used to interpret experiments. This methodology could be easily extended to the study of other phase transitions, from condensation to precipitation. Our analysis can help reconcile the textbook picture of classical nucleation theory with the growing body of atomistic studies and mesoscale models of solidification.

Volumetric properties and enthalpies of solution of alcohols CkH2k+1OH (k = 1, 2, 6) in 1-methyl-3-alkylimidazolium bis(trifluoromethylsulfonyl)imide {[C1CnIm][NTf2] n = 2, 4, 6, 8, 10} ionic liquids

We present a study on the effect of the alkyl chain length of the imidazolium ring in 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquids, [C1CnIm][NTf2] (n = 2 to 10), on the mixing properties of (ionic liquid + alcohol) mixtures (enthalpy and volume). We have measured small excess molar volumes with highly asymmetric curves as a function of mole fraction composition (S-shape) with more negative values in the alcohol-rich regions. The excess molar volumes increase with the increase of the alkyl-chain length of the imidazolium cation of the ionic liquid. The values of the partial molar excess enthalpy and the enthalpy of mixing are positive and, for the case of methanol, do not vary monotonously with the length of the alkyl side-chain of the cation on the ionic liquid – increasing from n = 2 to 6 and then decreasing from n = 8. This non-monotonous variation is explained by a more favourable interaction of methanol with the cation head group of the ionic liquid for alkyl chains longer than eight carbon atoms. It is also observed that the mixing is less favourable for the smaller alcohols, the enthalpy of mixing decreasing to less positive values as the alkyl chain of the alcohol increases. Based on the data from this work and on the knowledge of the vapour pressure of {[C1CnIm][NTf2] + alcohol} binary mixtures at T = 298 K reported in the literature, the excess Gibbs free energy, excess enthalpy and excess entropy could be then calculated and it was observed that these mixtures behave like the ones constituted by a non-associating and a non-polar component, with its solution behaviour being determined by the enthalpy.

Understanding the Performance of Melt-Extruded Poly(ethylene oxide)–Bicalutamide Solid Dispersions: Characterisation of Microstructural Properties Using Thermal, Spectroscopic and Drug Release Methods

In this article, we have prepared hot-melt-extruded solid dispersions of bicalutamide (BL) using poly(ethylene oxide) (PEO) as a matrix platform. Prior to preparation, miscibility of PEO and BL was assessed using differential scanning calorimetry (DSC). The onset of BL melting was signi?cantly depressed in the presence of PEO, and using Flory– Huggins (FH) theory, we identi?ed a negative value of -3.4, con?rming miscibility. Additionally, using FH lattice theory, we estimated the Gibbs free energy of mixing which was shown to be negative, passing through a minimum at a polymer fraction of 0.55. Using these data, solid dispersions at drug-to-polymer ratios of 1:10, 2:10 and 3:10 were prepared via hot-melt extrusion. Using a combination of DSC, powder X-ray diffractometry and scanning electron
microscopy, amorphous dispersions of BL were con?rmed at the lower two drug loadings. At the 3:10 BL to PEO ratio, crystalline BL was detected. The percent crystallinity of PEO was reduced by approximately 10% in all formulations following extrusion. The increased amorphous content within PEO following extrusion accommodated amorphous BL at drug to polymer loadings up to 2:10; however, the increased amorphous domains with PEO following extrusion were not suf?cient to fully accommodate BL at drug-to-polymer ratios of 3:10.

Low pressure carbon dioxide solubility in pure electrolyte solvents for lithium-ion batteries as a function of temperature. Measurement and prediction

Experimental values for the carbon dioxide solubility in eight pure electrolyte solvents for lithium ion batteries – such as ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), ethyl methyl carbonate (EMC), diethyl carbonate (DEC), ?-butyrolactone (?BL), ethyl acetate (EA) and methyl propionate (MP) – are reported as a function of temperature from (283 to 353) K and atmospheric pressure. Based on experimental solubility data, the Henry’s law constant of the carbon dioxide in these solvents was then deduced and compared with reported values from the literature, as well as with those predicted by using COSMO-RS methodology within COSMOthermX software and those calculated by the Peng–Robinson equation of state implemented into Aspen plus. From this work, it appears that the CO2 solubility is higher in linear carbonates (such as DMC, EMC, DEC) than in cyclic ones (EC, PC, ?BL). Furthermore, the highest CO2 solubility was obtained in MP and EA solvents, which are comparable to the solubility values reported in classical ionicliquids. The precision and accuracy of the experimental values, considered as the per cent of the relative average absolute deviations of the Henry’s law constants from appropriate smoothing equations and from literature values, are close to (1% and 15%), respectively. From the variation of the Henry’s law constants with temperature, the partial molar thermodynamic functions of dissolution such as the standard Gibbs free energy, the enthalpy, and the entropy are calculated, as well as the mixing enthalpy of the solvent with CO2 in its hypothetical liquid state.

Molecular dynamics simulation of nanoindentation of Fe3C and tetrahedral Fe4C

Study of nanomechanical response of iron carbides is important because presence of iron carbides greatly influences the performance and longevity of steel components. This work contributes to the literature by exploring nanoindentation of nanocrystalline Fe3C and tetrahedral-Fe4C using molecular dynamics simulation. The chemical interactions of iron and carbon were described through an analytical bond order inter-atomic potential (ABOP) energy function. The indentations were performed at an indentation speed of 50 m/sec and a repeat trial was performed at 5 m/sec. Load-displacement (P-h) curve for both these carbides showed residual indentation depth and maximum indentation depth (hf/hmax) ratio to be higher than 0.7 i.e. a circumstance where Oliver and Pharr method was not appropriate to be applied to evaluate the material properties. Alternate evaluation revealed Fe3C to be much harder than Fe4C. Gibbs free energy of formation and radial distribution function, coupled with state of the average local temperature and von Mises stresses indicate the formation of a new phase of iron-carbide. Formation of this newer phase was found to be due to deviatoric strain rather than the high temperature induced in the substrate during nanoindentation

Low pressure methane solubility in lithium-ion batteries based solvents and electrolytes as a function of temperature. Measurement and prediction

The methane solubility in five pure electrolyte solvents and one binary solvent mixture for lithium ion batteries – such as ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), ethyl methyl carbonate (EMC), diethyl carbonate (DEC) and the (50:50 wt%) mixture of EC:DMC was studied experimentally at pressures close to atmospheric and as a function of temperature between (280 and 343) K by using an isochoric saturation technique. The effect of the selected anions of a lithium salt LiX (X = hexafluorophosphate,

Levoglucosan formation mechanisms during cellulose pyrolysis

Levoglucosan is one important primary product during cellulose pyrolysis either as an intermediate or as a product. Three available mechanisms for levoglucosan formation have been studied theoretically in this paper, which are free-radical mechanism; glucose intermediate mechanism; and levoglucosan chain-end mechanism. All the elementary reactions included in the pathway of every mechanism were investigated; thermal properties including activation energy, Gibbs free energy, and enthalpy for every pathway were also calculated. It was concluded that free-radical mechanism has the highest energy barrier during the three levoglucosan formation mechanisms, glucose intermediate mechanism has lower energy barrier than free-radical mechanism, and levoglucosan chain-end mechanism is the most reasonable pathway because of the lowest energy barrier. By comparing with the activation energy obtained from the experimental results, it was also concluded that levoglucosan chain-end mechanism fits better with the experimental data for the formation of levoglucosan. © 2013 Elsevier B.V. All rights reserved.

Probing The effects of Experimental Conditions on the Character of Drug-Polymer Phase Diagrams Constructed Using Flory-Huggins Theory

Purpose: Amorphous drug-polymer solid dispersions have been found to result in improved drug dissolution rates when compared to their crystalline counterparts. However, when the drug exists in the amorphous form it will possess a higher Gibb’s free energy than its associated crystalline state and can recrystallize. Drug-polymer phase diagrams constructed through the application of the Flory Huggins (F-H) theory contain a wealth of information regarding thermodynamic and kinetic stability of the amorphous drug-polymer system. This study was aimed to evaluate the effects of various experimental conditions on the solubility and miscibility detections of drug-polymer binary system. Methods: Felodipine (FD)-Polyvinylpyrrolidone (PVP) K15 (PVPK15) and FD-Polyvinylpyrrolidone/vinyl acetate (PVP/VA64) were the selected systems for this research. Physical mixtures with different drug loadings were mixed and ball milled. These samples were then processed using Differential Scanning Calorimetry (DSC) and measurements of melting point (Tend) and glass transition (Tg) were detected using heating rates of 0.5, 1.0 and 5.0°C/min. Results: The melting point depression data was then used to calculate the F-H interaction parameter (χ) and extrapolated to lower temperatures to complete the liquid–solid transition curves. The theoretical binodal and spinodal curves were also constructed which were used to identify regions within the phase diagram. The effects of polymer selection, DSC heating rate, time above parent polymer Tg and polymer molecular weight were investigated by identifying amorphous drug miscibility limits at pharmaceutically relevant temperatures. Conclusion: The potential implications of these findings when applied to a non-ambient processing method such as Hot Melt Extrusion (HME) are also discussed.

Density, conductivity, viscosity, and excess properties of (pyrrolidinium nitrate-based Protic Ionic Liquid + propylene carbonate) binary mixture

Density, ?, viscosity, ?, and conductivity, s, measurements of binary mixtures containing the pyrrolidinium nitrate Protic Ionic Liquid (PIL) and propylene carbonate (PC), are determined at the atmospheric pressure as a function of the temperature from (283.15 to 353.15) K and within the whole composition range. The temperature dependence of both the viscosity and conductivity of each mixture exhibits a non-Arrhenius behaviour, but is correctly fitted by using the Vogel–Tamman–Fulcher (VTF) equation. In each case, the best-fit parameters, such as the pseudo activation energy, View the MathML source and ideal glass transition temperature, T0 are then extracted. The excess molar volumes VE, and viscosity deviations from the ideality, ??, of each investigated mixture were then deduced from the experimental results, as well as, their apparent molar volumes, V?, thermal expansion coefficients ap, and excess Gibbs free energies (?G*E) of activation of viscous flow. The VE, apE, ?? values are negative over the whole composition range for each studied temperature therein. According to the Walden rule, the ionicity of each mixture was then evaluated as a function of the temperature from (283.15 to 353.15) K and of the composition. Results have been then discussed in terms of molecular interactions and molecular structures in this binary mixture.

The phase diagram of water at negative pressures: Virtual ices

The phase diagram of water at negative pressures as obtained from computer simulations for two models of water, TIP4P/2005 and TIP5P is presented. Several solid structures with lower densities than ice Ih, so-called virtual ices, were considered as possible candidates to occupy the negative pressure region of the phase diagram of water. In particular the empty hydrate structures sI, sII, and sH and another, recently proposed, low-density ice structure. The relative stabilities of these structures at 0 K was determined using empirical water potentials and density functional theory calculations. By performing free energy calculations and Gibbs-Duhem integration the phase diagram of TIP4P/2005 was determined at negative pressures. The empty hydrates sII and sH appear to be the stable solid phases of water at negative pressures. The phase boundary between ice Ih and sII clathrate occurs at moderate negative pressures, while at large negative pressures sH becomes the most stable phase. This behavior is in reasonable agreement with what is observed in density functional theory calculations.