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.

Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

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.

A stabilization mechanism of zirconia based on oxygen vacancies only

The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO2) is analyzed by means of a self-consistent tight-binding model. Using this model, energies and structures of zirconia containing different vacancy concentrations are calculated, equivalent in concentration to the charge compensating vacancies associated with dissolved yttria (Y2O3) in the tetragonal and cubic phase fields (3.2 and 14.4% mol, respectively). The model is shown to predict the large relaxations around an oxygen vacancy, and the clustering of vacancies along the 111 directions, in good agreement with experiments and first principles calculations. The vacancies alone are shown to explain the stabilization of cubic zirconia, and the mechanism is analyzed. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Electron energy-loss near-edge shape as a probe to investigate the stabilization of yttria-stabilized zirconia

The electron energy-loss near-edge structure (ELNES) at the O K edge has been studied in yttria-stabilized zirconia (YSZ). The electronic structure of YSZ for compositions between 3 and 15 mol % Y2O3 has been computed using a pseudopotential-based technique to calculate the local relaxations near the O vacancies. The results showed phase transition from the tetragonal to cubic YSZ at 10 mol % of Y2O3, reproducing experimental observations. Using the relaxed defect geometry, calculation of the ELNES was carried out using the full-potential linear muffin-tin orbital method. The results show very good agreement with the experimental O K-edge signal, demonstrating the power of using ELNES to probe the stabilization mechanism in doped metal oxides.

Effect of relaxation on the oxygen K-edge electron energy-loss near-edge structure in yttria-stabilized zirconia

The electron energy-loss near-edge structure (ELNES) at the oxygen K-edge has been investigated in a range of yttria-stabilized zirconia (YSZ) materials. The electronic structure of the three polymorphs of pure ZrO2 and of the doped YSZ structure close to the 33 mol %Y2O3 composition have been calculated using a full-potential linear muffin-tin orbital method (NFP-LMTO) as well as a pseudopotential based technique. Calculations of the ELNES dipole transition matrix elements in the framework of the NFP-LMTO scheme and inclusion of core hole screening within Slater's transition state theory enable the ELNES to be computed. Good agreement between the experimental and calculated ELNES is obtained for pure monoclinic ZrO2. The agreement is less good with the ideal tetragonal and cubic structures. This is because the inclusion of defects is essential in the calculation of the YSZ ELNES. If the model used contains ordered defects such as vacancies and metal Y planes, agreement between the calculated and experimental O K-edges is significantly improved. The calculations show how the five different O environments of Zr,Y,O, are connected with the features observed in the experimental spectra and demonstrate clearly the power of using ELNES to probe the stabilization mechanism in doped metal oxides.

A reactive force field simulation of liquid-liquid phase transitions in phosphorus

A force field model of phosphorus has been developed based on density functional (DF) computations and experimental results, covering low energy forms of local tetrahedral symmetry and more compact (simple cubic) structures that arise with increasing pressure. Rules tailored to DF data for the addition, deletion, and exchange of covalent bonds allow the system to adapt the bonding configuration to the thermodynamic state. Monte Carlo simulations in the N-P-T ensemble show that the molecular (P-4) liquid phase, stable at low pressure P and relatively low temperature T, transforms to a polymeric (gel) state on increasing either P or T. These phase changes are observed in recent experiments at similar thermodynamic conditions, as shown by the close agreement of computed and measured structure factors in the molecular and polymer phases. The polymeric phase obtained by increasing pressure has a dominant simple cubic character, while the polymer obtained by raising T at moderate pressure is tetrahedral. Comparison with DF results suggests that the latter is a semiconductor, while the cubic form is metallic. The simulations show that the T-induced polymerization is due to the entropy of the configuration of covalent bonds, as in the polymerization transition in sulfur. The transition observed with increasing P is the continuation at high T of the black P to arsenic (A17) structure observed in the solid state, and also corresponds to a semiconductor to metal transition. (C) 2004 American Institute of Physics.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Phase transitions in epitaxial Ba0.5Sr0.5TiO3 thin films

Ba0.5Sr0.5TiO3 (BST) thin-film capacitor structures with various thicknesses, (50-1200 nm) and different strain conditions (on lanthanum strontium cobalt oxide La0.5Sr0.5CoO3 and strontium ruthenate SrRuO3 buffer layers) were made using pulsed laser deposition, and characterized by x-ray diffraction. The out-of-plane lattice parameter was followed as a function of temperature within the 100-300 K temperature interval. The phase sequence (cubic-tetragonal-orthorhombic-rhombohedral) known to exist in the bulk analog is shown to be strongly affected by both the stress conditions imposed by the buffer layer and the thickness of the BST film itself. Thus, no phase transition was found for the in-plane compressed BST films. On the stress-free BST films, on the contrary, more phase transitions were observed. It appeared that the complexity of structural phase transitions increased as the film thickness in this system was reduced.

Hydrogenation of S to H2S on Pt(111): A first-principles study

Density-functional theory has been used to investigate the chemisorption of S, SH, and H2S as well as the coadsorption of S and H and SH and H on Pt(111). In addition reaction pathways and energy profiles for the conversion of adsorbed S and H into gas-phase H2S have been determined. It has been found that S, SH, and H2S bind preferentially at face-centered-cubic (fcc), bridge, and top sites, respectively. Both the S+H and SH+H reactions have high barriers (similar to1 eV) and high exothermicities (similar to1 eV). This reveals that adsorbed H2S and SH are highly unstable adsorbates on Pt(111) and that adsorbed S (and H) is the most stable SHX (X=0,1,2) intermediate on Pt(111) (C) 2001 American Institute of Physics.

Origin of Ferroelastic domains in Free-Standing Single Crystal Ferroelectric Films

The origin of the unusual 90 degrees ferroelectric/ferroelastic domains, consistently observed in recent studies on mesoscale and nanoscale free-standing single crystals of BaTiO3 [Schilling , Phys. Rev. B 74, 024115 (2006); Schilling , Nano Lett. 7, 3787 (2007)], has been considered. A model has been developed which postulates that the domains form as a response to elastic stress induced by a surface layer which does not undergo the paraelectric-ferroelectric cubic-tetragonal phase transition. This model was found to accurately account for the changes in domain periodicity as a function of size that had been observed experimentally. The physical origin of the surface layer might readily be associated with patterning damage, seen in experiment; however, when all evidence of physical damage is removed from the BaTiO3 surfaces by thermal annealing, the domain configuration remains practically unchanged. This suggests a more intrinsic origin, such as the increased importance of surface tension at small dimensions. The effect of surface tension is also shown to be proportional to the difference in hardness between the surface and the interior of the ferroelectric. The present model for surface-tension induced twinning should also be relevant for finely grained or core-shell structured ceramics.