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.

Measuring the Characteristic Function of the Work Distribution

We propose an interferometric setting for the ancilla-assisted measurement of the characteristic function of the work distribution following a time-dependent process experienced by a quantum system. We identify how the configuration of the effective interferometer is linked to the symmetries enjoyed by the Hamiltonian ruling the process and provide the explicit form of the operations to implement in order to accomplish our task. We finally discuss two physical settings, based on hybrid optomechanical-electromechanical devices, where the theoretical proposals discussed in our work could find an experimental demonstration.

Detecting the work statistics through Ramsey-like interferometry

Out-of-equilibrium statistical mechanics is attracting considerable interest due to the recent advances in the control and manipulations of systems at the quantum level. Recently, an interferometric scheme for the detection of the characteristic function of the work distribution following a time-dependent process has been proposed [L. Mazzola et al., Phys. Rev. Lett. 110 (2013) 230602]. There, it was demonstrated that the work statistics of a quantum system undergoing a process can be reconstructed by effectively mapping the characteristic function of work on the state of an ancillary qubit. Here, we expand that work in two important directions. We first apply the protocol to an interesting specific physical example consisting of a superconducting qubit dispersively coupled to the field of a microwave resonator, thus enlarging the class of situations for which our scheme would be key in the task highlighted above. We then account for the interaction of the system with an additional one (which might embody an environment), and generalize the protocol accordingly.

Nonequilibrium Quantum Landauer Principle

Using the operational framework of completely positive, trace preserving operations and thermodynamic fluctuation relations, we derive a lower bound for the heat exchange in a Landauer erasure process on a quantum system. Our bound comes from a nonphenomenological derivation of the Landauer principle which holds for generic nonequilibrium dynamics. Furthermore, the bound depends on the nonunitality of dynamics, giving it a physical significance that differs from other derivations. We apply our framework to the model of a spin-1/2 system coupled to an interacting spin chain at finite temperature.

A simple model of atomic interactions in noble metals based explicitly on electronic structure

A total energy tight-binding model with a basis of just one s state per atom is introduced. It is argued that this simplest of all tight-binding models provides a surprisingly good description of the structural stability and elastic constants of noble metals. By assuming inverse power scaling laws for the hopping integrals and the repulsive pair potential, it is shown that the density matrix in a perfect primitive crystal is independent of volume, and structural energy differences and equations of state are then derived analytically. The model is most likely to be of use when one wishes to consider explicitly and self-consistently the electronic and atomic structures of a generic metallic system, with the minium of computation expense. The relationship to the free-electron jellium model is described. The applicability of the model to other metals is also considered briefly.

Binding of Calcium and Carbonate to Polyacrylates

Polyacrylate molecules can be used to slow the growth of calcium carbonate. However, little is known about the mechanism by which the molecules impede the growth rate. A recent computational study (Bulo et al. Macromolecules 2007, 40, 3437) used metadynamics to investigate the binding of calcium to polyacrylate chains and has thrown some light on the coiling and precipitation of these polymers. We extend these simulations to examine the binding of calcium and carbonate to polyacrylate chains. We show that calcium complexed with both carbonate and polyacrylate is a very stable species. The free energies of calcium-carbonate-polyacrylate complexes, with different polymer configurations, are calculated, and differences in the free energy of the binding of carbonate are shown to be due to differences in the amount of steric hindrance about the calcium, which prevents the approach of the carbonate ion.

Protection of DNA against low-energy electrons by amino acids: A first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, we have investigated the notion that amino acids can play a protective role when DNA is exposed to excess electrons produced by ionizing radiation. In this study we focus on the interaction of glycine with the DNA nucleobase thymine. We studied thymine-glycine dimers and a condensed phase model consisting of one thymine molecule solvated in amorphous glycine. Our results show that the amino acid acts as a protective agent for the nucleobase in two ways. If the excess electron is initially captured by the thymine, then a proton is transferred in a barrier-less way from a neighboring hydrogen-bonded glycine. This stabilizes the excess electron by reducing the net partial charge on the thymine. In the second mechanism the excess electron is captured by a glycine, which acts as a electron scavenger that prevents electron localization in DNA. Both these mechanisms introduce obstacles to further reactions of the excess electron within a DNA strand, e.g. by raising the free energy barrier associated with strand breaks.

Understanding the Interaction between Low-Energy Electrons and DNA Nucleotides in Aqueous Solution

Reactions that can damage DNA have been simulated using a combination of molecular dynamics and density functional theory. In particular, the damage caused by the attachment of a low energy electron to the nucleobase. Simulations of anionic single nucleotides of DNA in an aqueous environment that was modeled explicitly have been performed. This has allowed us to examine the role played by the water molecules that surround the DNA in radiation damage mechanisms. Our simulations show that hydrogen bonding and protonation of the nucleotide by the water can have a significant effect on the barriers to strand breaking reactions. Furthermore, these effects are not the same for all four of the bases.

A Phase-field Model for Computer Simulation of Lamellar Structure Formation in y-TiAI

In the present paper, a phase-field model is developed to simulate the formation and evolution of lamellar microstructure in γ-TiAl alloys. The mechanism of formation of TiAl lamellae proposed by Denquin and Naka is incorporated into the model. The model describes the formation and evolution of the face-centered cubic (fcc) stacking lamellar zone followed by the subsequent appearance and growth of the γ-phase, involving both the chemical composition change by atom transfer and the ordering of the fcc lattice. The thermodynamics of the model system and the interaction between the displacive and diffusional transformations are described by a non-equilibrium free energy formulated as a function of concentration and structural order parameter fields. The long-range elastic interactions, arising from the lattice misfit between the α, fcc (A1) and the various orientation variants of the γ-phase are taken into account by incorporating of the elastic strain energy into the total free energy. Simulation studies based on the model successfully predicted some essential features of the lamellar structure. It is found that the formation and evolution of the lamellar structure are predominantly controlled by the minimization of the elastic energy of the interfaces between the different fcc stacking groups, low-symmetry product phase γ and the high-symmetry α-phase, as well as between the various orientation variants of the product phase.

Stability of Sr adatom model structures for SrTiO3(001) surface reconstructions

We report results of first-principles calculations on the thermodynamic stability of different Sr adatom structures that have been proposed to explain some of the observed reconstructions of the (001) surface of strontium titanate (Kubo and Nozoye 2003 Surf Sci. 542 177). From surface free energy calculations, a phase diagram is constructed indicating the range of conditions over which each structure is most stable. These results are compared with Kubo and Nozoye's experimental observations. It is concluded that low Sr adatom coverage structures can only be explained if the surface is far from equilibrium. Intermediate coverage structures are stable only if the surface is in or very nearly in equilibrium with the strontium oxide.