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.

Anharmonic Helmholtz free energy to O(lambdaâ ´) in the non- leading term approximation /|nby Leslie Wilk. -- 260 St. Catharines, Ont. : [s. n.],

Expressions for the anharmonic Helmholtz free energy contributions up to o( f ) ,valid for all temperatures, have been obtained using perturbation theory for a c r ystal in which every atom is on a site of inversion symmetry. Numerical calculations have been carried out in the high temperature limit and in the non-leading term approximation for a monatomic facecentred cubic crystal with nearest neighbour c entralforce interactions. The numbers obtained were seen to vary by a s much as 47% from thos e obtai.ned in the leading term approximati.on,indicating that the latter approximati on is not in general very good. The convergence to oct) of the perturbation series in the high temperature limit appears satisfactory.

On the Ewald method and the quartic Helmholtz free energy of an anharmonic metallic crystal

Thesis (M. Sc.) - Brock University, 1978.

An anharmonic contribution to the Helmholtz free energy O(lambda 6)

The anharmonic contributions of order A6 to the Helmholtz free energy for a crystal in which every atom is on a site of inversion symmetry, have been evaluated The cor~esponding diagrams in the various orders of the perturbation theory have been presented The validity of the expressions given is for high temperatures. Numerical calculations for the diagrams which contribute to the free energy have been worked out for a nearest-n~ighbour central-force model of a facecentered cubic lattice in the high-temperature limit and in the leading term and the Ludwig approximations. The accuracy of the Ludwig approximation in evaluating the Brillouin-zone sums has been investigated. Expansion for all diagrams in the high-temperature limit has been carried out The contribution to the specific heat involves a linear as well as cubic term~ We have applied Lennard-Jones, Morse and Exponential 6 types of potentials. A comparison between the contribution to the free energy of order A6 to that of order A4 has been made.

The Landau free energy of hard ellipses obtained from microscopic simulations

Systems of two-dimensional hard ellipses of varying aspect ratios and packing fractions are studied by Monte Carlo simulations in the generalised canonical ensemble. From this microscopic model, we extract a coarse-grained macroscopic Landau-de Gennes free energy as a function of packing fraction and orientational order parameter. We separate the free energy into the ideal orientational entropy of non-interacting two-dimensional spins and an excess free energy associated with excluded volume interactions. We further explore the isotropic-nematic phase transition using our empirical expression for the free energy and find that the nature of the phase transition is continuous for the aspect ratios we studied.

Solvatochromism in Binary Mixtures: First Report on a Solvation Free Energy Relationship between Solvent Exchange Equilibrium Constants and the Properties of the Medium

We have employed UV-vis spectroscopy in order to investigate details of the solvation of six solvatochromic indicators, hereafter designated as ""probes"", namely, 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr,, respectively. These can be divided into three pairs, each includes two probes of similar pK(a) in water and different lipophilicity. Solvation has been studied in binary mixtures, BMs, of water, W, with 12 protic organic solvents, S, including mono- and bifunctional alcohols (2-alkoxyethanoles, unsaturated and chlorinated alcohols). Each medium was treated as a mixture of S, W, and a complex solvent, S-W, formed by hydrogen bonding. Values of lambda(max) (of the probe intramolecular charge transfer) were converted into empirical polarity scales, E(T)(probe) in kcal/mol, whose values were correlated with the effective mole fraction of water in the medium, chi w(effective). This correlation furnished three equilibrium constants for the exchange of solvents in the probe solvation shell; phi(W/S) (W substitutes S): phi(S-W/W) (S-W substitutes W), and phi(S-W/S) (S-W substitutes S), respectively. The values of these constants depend on the physicochemical properties of the probe and the medium. We tested, for the first time, the applicability of a new solvation free energy relationship: phi = constant + a alpha(BM) + b beta(BM) + s(pi*(BM) + d delta) + p log P(BM), where a, b, s, and p are regression coefficients alpha(BM), beta(BM), and pi*(BM) are solvatochromic parameters of the BM, delta is a correction term for pi*, and log P is an empirical scale of lipophilicity. Correlations were carried out with two-, three-, and four-medium descriptors. In all cases, three descriptors gave satisfactory correlations; use of four parameters gave only a marginal increase of the goodness of fit. For phi(W/S), the most important descriptor was found to be the lipophilicity of the medium; for phi(S-W/W) and phi(S-W/S), solvent basicity is either statistically relevant or is the most important descriptor. These responses are different from those of E(T)(probe) of many solvatochromic indicators in pure solvents, where the importance of solvent basicity is usually marginal, and can be neglected.

Non-quantized minimum free energy in untranslated region exons

In an attempt to improve automated gene prediction in the untranslated region of a gene, we completed an in-depth analysis of the minimum free energy for 8,689 sub-genetic DNA sequences. We expanded Zhang's classification model and classified each sub-genetic sequence into one of 27 possible motifs. We calculated the minimum free energy for each motif to explore statistical features that correlate to biologically relevant sub-genetic sequences. If biologically relevant sub-genetic sequences fall into distinct free energy quanta it may be possible to characterize a motif based on its minimum free energy. Proper characterization of motifs can lead to greater understanding in automated genefinding, gene variability and the role DNA structure plays in gene network regulation.

Our analysis determined: (1) the average free energy value for exons, introns and other biologically relevant sub-genetic sequences, (2) that these subsequences do not exist in distinct energy quanta, (3) that introns exist however in a tightly coupled average minimum free energy quantum compared to all other biologically relevant sub-genetic sequence types, (4) that single exon genes demonstrate a higher stability than exons which span the entire coding sequence as part of a multi-exon gene and (5) that all motif types contain a free energy global minimum at approximately nucleotide position 1,000 before reaching a plateau. These results should be relevant to the biochemist and bioinformatician seeking to understand the relationship between sub-genetic sequences and the information behind them.

Machine to die for: the quest for free energy

This is about the search for 'perpetual motion' and 'free energy'. Conventional science claims that it is impossible, yet generations of inventors have been mesmerised by the promise of an engine that powers itself. The world's reliance on diminishing fossil fuel resources and the associated problems of pollution serve to spur them on. It showcases a number of dedicated, sometimes eccentric, and always obsessive individuals who have devoted their lives to this quest.