The effects of microstructure, Nb content and secondary Ruddlesden-Popper phase on thermoelectric properties in perovskite CaMn1-x Nb x O3 (x = 0-0.10) thin films.

CaMn1-x Nb x O3 (x = 0, 0.5, 0.6, 0.7 and 0.10) thin films have been grown by a two-step sputtering/annealing method. First, rock-salt-structured (Ca,Mn1-x ,Nb x )O thin films were deposited on 11̄00 sapphire using reactive RF magnetron co-sputtering from elemental targets of Ca, Mn and Nb. The CaMn1-x Nb x O3 films were then obtained by thermally induced phase transformation from rock-salt-structured (Ca,Mn1-x Nb x )O to orthorhombic during post-deposition annealing at 700 °C for 3 h in oxygen flow. The X-ray diffraction patterns of pure CaMnO3 showed mixed orientation, while Nb-containing films were epitaxially grown in [101] out of-plane-direction. Scanning transmission electron microscopy showed a Ruddlesden-Popper (R-P) secondary phase in the films, which results in reduction of the electrical and thermal conductivity of CaMn1-x Nb x O3. The electrical resistivity and Seebeck coefficient of the pure CaMnO3 film were measured to 2.7 Ω cm and -270 µV K-1 at room temperature, respectively. The electrical resistivity and Seebeck coefficient were reduced by alloying with Nb and was measured to 0.09 Ω cm and -145 µV K-1 for x = 0.05. Yielding a power factor of 21.5 µW K-2 m-1 near room temperature, nearly eight times higher than for pure CaMnO3 (2.8 µW K-2 m-1). The power factors for alloyed samples are low compared to other studies on phase-pure material. This is due to high electrical resistivity originating from the secondary R-P phase. The thermal conductivity of the CaMn1-x Nb x O3 films is low for all samples and is the lowest for x = 0.07 and 0.10, determined to 1.6 W m-1 K-1. The low thermal conductivity is attributed to grain boundary scattering and the secondary R-P phase.

Superioniclike Diffusion in an Elemental Crystal: bcc Titanium.

Recent theoretical investigations [A. B. Belonoshko et al. Nat. Geosci. 10, 312 (2017)1752-089410.1038/ngeo2892] revealed the occurrence of the concerted migration of several atoms in bcc Fe at inner-core temperatures and pressures. Here, we combine first-principles and semiempirical atomistic simulations to show that a diffusion mechanism analogous to the one predicted for bcc iron at extreme conditions is also operative and of relevance for the high-temperature bcc phase of pure Ti at ambient pressure. The mechanism entails a rapid collective movement of numerous (from two to dozens) neighbors along tangled closed-loop paths in defect-free crystal regions. We argue that this phenomenon closely resembles the diffusion behavior of superionics and liquid metals. Furthermore, we suggest that concerted migration is the atomistic manifestation of vanishingly small ω-mode phonon frequencies previously detected via neutron scattering and the mechanism underlying anomalously large and markedly non-Arrhenius self-diffusivities characteristic of bcc Ti.

Metastable silica high pressure polymorphs as structural proxies of deep Earth silicate melts.

Modelling of processes involving deep Earth liquids requires information on their structures and compression mechanisms. However, knowledge of the local structures of silicates and silica (SiO2) melts at deep mantle conditions and of their densification mechanisms is still limited. Here we report the synthesis and characterization of metastable high-pressure silica phases, coesite-IV and coesite-V, using in situ single-crystal X-ray diffraction and ab initio simulations. Their crystal structures are drastically different from any previously considered models, but explain well features of pair-distribution functions of highly densified silica glass and molten basalt at high pressure. Built of four, five-, and six-coordinated silicon, coesite-IV and coesite-V contain SiO6 octahedra, which, at odds with 3rd Pauling's rule, are connected through common faces. Our results suggest that possible silicate liquids in Earth's lower mantle may have complex structures making them more compressible than previously supposed.

Effects of configurational disorder on the elastic properties of icosahedral boron-rich alloys based on B6O, B13C2, and B4C, and their mixing thermodynamics.

The elastic properties of alloys between boron suboxide (B6O) and boron carbide (B13C2), denoted by (B6O)(1-x)(B13C2)(x), as well as boron carbide with variable carbon content, ranging from B13C2 to B4C are calculated from first-principles. Furthermore, the mixing thermodynamics of (B6O)(1-x)(B13C2)(x) is studied. A superatom-special quasirandom structure approach is used for modeling different atomic configurations, in which effects of configurational disorder between the carbide and suboxide structural units, as well as between boron and carbon atoms within the units, are taken into account. Elastic properties calculations demonstrate that configurational disorder in B13C2, where a part of the C atoms in the CBC chains substitute for B atoms in the B12 icosahedra, drastically increase the Young's and shear modulus, as compared to an atomically ordered state, B12(CBC). These calculated elastic moduli of the disordered state are in excellent agreement with experiments. Configurational disorder between boron and carbon can also explain the experimentally observed almost constant elastic moduli of boron carbide as the carbon content is changed from B4C to B13C2. The elastic moduli of the (B6O)(1-x)(B13C2)(x) system are also practically unchanged with composition if boron-carbon disorder is taken into account. By investigating the mixing thermodynamics of the alloys, in which the Gibbs free energy is determined within the mean-field approximation for the configurational entropy, we outline the pseudo-binary phase diagram of (B6O)(1-x)(B13C2)(x). The phase diagram reveals the existence of a miscibility gap at all temperatures up to the melting point. Also, the coexistence of B6O-rich as well as ordered or disordered B13C2-rich domains in the material prepared through equilibrium routes is predicted.

Hyperfine splitting and room-temperature ferromagnetism of Ni at multimegabar pressure.

Magnetic and elastic properties of Ni metal have been studied up to 260 GPa by nuclear forward scattering of synchrotron radiation with the 67.4 keV Mössbauer transition of 61Ni. The observed magnetic hyperfine splitting confirms the ferromagnetic state of Ni up to 260 GPa, the highest pressure where magnetism in any material has been observed so far. Ab initio calculations reveal that the pressure evolution of the hyperfine field, which features a maximum in the range of 100 to 225 GPa, is a relativistic effect. The Debye energy obtained from the Lamb-Mössbauer factor increases from 33 meV at ambient pressure to 60 meV at 100 GPa. The change of this energy over volume compression is well described by a Grüneisen parameter of 2.09.

Charge redistribution mechanisms of ceria reduction.

Charge redistribution at low oxygen vacancy concentrations in ceria have been studied in the framework of the density functional theory. We propose a model to approach the dilute limit using the results of supercell calculations. It allows one to reproduce the characteristic experimentally observed behavior of composition versus oxygen pressure dependency. We show that in the dilute limit the charge redistribution is likely to be driven by a mechanism different from the one involving electron localization on cerium atoms. We demonstrate that it can involve charge localization on light element impurities.

High-pressure--high-temperature polymorphism in ta: resolving an ongoing experimental controversy.

Phase diagrams of refractory metals remain essentially unknown. Moreover, there is an ongoing controversy over the high-pressure melting temperatures of these metals: results of diamond anvil cell (DAC) and shock wave experiments differ by at least a factor of 2. From an extensive ab initio study on tantalum we discovered that the body-centered cubic phase, its physical phase at ambient conditions, transforms to another solid phase, possibly hexagonal omega phase, at high temperature. Hence the sample motion observed in DAC experiments is very likely not due to melting but internal stresses accompanying a solid-solid transformation, and thermal stresses associated with laser heating.

Superhard semiconducting optically transparent high pressure phase of boron.

An orthorhombic (space group Pnnm) boron phase was synthesized at pressures above 9 GPa and high temperature, and it was demonstrated to be stable at least up to 30 GPa. The structure, determined by single-crystal x-ray diffraction, consists of B12 icosahedra and B2 dumbbells. The charge density distribution obtained from experimental data and ab initio calculations suggests covalent chemical bonding in this phase. Strong covalent interatomic interactions explain the low compressibility value (bulk modulus is K300=227 GPa) and high hardness of high-pressure boron (Vickers hardness HV=58 GPa), after diamond the second hardest elemental material.

Molybdenum at high pressure and temperature: melting from another solid phase.

The Gibbs free energies of bcc and fcc Mo are calculated from first principles in the quasiharmonic approximation in the pressure range from 350 to 850 GPa at room temperatures up to 7500 K. It is found that Mo, stable in the bcc phase at low temperatures, has lower free energy in the fcc structure than in the bcc phase at elevated temperatures. Our density-functional-theory-based molecular dynamics simulations demonstrate that fcc melts at higher than bcc temperatures above 1.5 Mbar. Our calculated melting temperatures and bcc-fcc boundary are consistent with the Mo Hugoniot sound speed measurements. We find that melting occurs at temperatures significantly above the bcc-fcc boundary. This suggests an explanation of the recent diamond anvil cell experiments, which find a phase boundary in the vicinity of our extrapolated bcc-fcc boundary.

Pure iron compressed and heated to extreme conditions.

The results of a first-principles study supported by the temperature-quenched laser-heated diamond anvil-cell experiments on the high-pressure high-temperature structural behavior of pure iron are reported. We show that in contrast to the widely accepted picture, the face-centered cubic (fcc) phase becomes as stable as the hexagonal-close-packed (hcp) phase at pressures around 300-360 GPa and temperatures around 5000-6000 K. Our temperature-quenched experiments indicate that the fcc phase of iron can exist in the pressure-temperature region above 160 GPa and 3700 K, respectively. This, in particular, means that the actual structure of the Earth's core may be a complex phase with a large number of stacking faults.

Body-centered cubic iron-nickel alloy in Earth's core.

Cosmochemical, geochemical, and geophysical studies provide evidence that Earth's core contains iron with substantial (5 to 15%) amounts of nickel. The iron-nickel alloy Fe(0.9)Ni(0.1) has been studied in situ by means of angle-dispersive x-ray diffraction in internally heated diamond anvil cells (DACs), and its resistance has been measured as a function of pressure and temperature. At pressures above 225 gigapascals and temperatures over 3400 kelvin, Fe(0.9)Ni(0.1) adopts a body-centered cubic structure. Our experimental and theoretical results not only support the interpretation of shockwave data on pure iron as showing a solid-solid phase transition above about 200 gigapascals, but also suggest that iron alloys with geochemically reasonable compositions (that is, with substantial nickel, sulfur, or silicon content) adopt the bcc structure in Earth's inner core.

Noblest of all metals is structurally unstable at high pressure.

In a series of experiments in externally electrically heated diamond anvil cells we demonstrate that at pressures above approximately 240 GPa gold adopts a hexagonal-close-packed structure. Ab initio calculations predict that at pressures about 250 GPa different stacking sequences of close-packed atomic layers in gold become virtually degenerate in energy, strongly supporting the experimental observations.

Electron transport in stretched monoatomic gold wires.

The conductance of monoatomic gold wires containing 3-7 gold atoms has been obtained from ab initio calculations. The transmission is found to vary significantly depending on the wire stretching and the number of incorporated atoms. Such oscillations are determined by the electronic structure of the one-dimensional (1D) part of the wire between the contacts. Our results indicate that the conductivity of 1D wires can be suppressed without breaking the contact.

First-principles theory of intermediate-valence f-electron systems.

We propose a first-principles based method for calculating the electronic structure and total energy of solids in an intermediate-valence configuration. The method takes into account correlation effects (d-f Coulomb interaction) and many-body renormalization of the effective hybridization parameter of the f system. As an example, the formation of a pressure-induced intermediate-valence state in Yb is considered and its electronic structure and equation of state are calculated and compared to experimental data. The agreement is found to be excellent for both properties, and we argue that the developed method, which applies to any element or compound, provides for the first time a quantitative theoretical treatment of intermediate-valence materials.

Fluctuating lattice constants of indium under high pressure.

Recent high-pressure investigations of elemental In have yielded controversial results. We show that the observed high-pressure face-centered orthorhombic (fco) structure can be explained as an intermediate state between two body-centered tetragonal (bct) structures with different c/a ratios (c/a < square root [2] and c/a > square root [2], respectively). In a pressure range from about 50 to 200 GPa these two bct structures correspond to local minima of the total energy with respect to orthorhombic distortion of the ground-state bct In structure. The fco saddle point represents a tiny barrier and even at low temperatures rapid structural fluctuations should occur. Such a situation has not been identified in any other elemental metal.

High-pressure melting of molybdenum.

The melting curve of the body-centered cubic (bcc) phase of Mo has been determined for a wide pressure range using both direct ab initio molecular dynamics simulations of melting as well as a phenomenological theory of melting. These two methods show very good agreement. The simulations are based on density functional theory within the generalized gradient approximation. Our calculated equation of state of bcc Mo is in excellent agreement with experimental data. However, our melting curve is substantially higher than the one determined in diamond anvil cell experiments up to a pressure of 100 GPa. An explanation is suggested for this discrepancy.

Origin of magnetic anisotropy of Gd metal.

Using first-principles theory, we have calculated the energy of Gd as a function of spin direction, theta, between the c and a axes and found good agreement with experiment for both the total magnetic anisotropy energy and its angular dependence. The calculated low temperature direction of the magnetic moment lies at an angle of 20 degrees to the c axis. The calculated magnetic anisotropy energy of Gd metal is due to a unique mechanism involving a contribution of 7.5 microeV from the classical dipole-dipole interaction between spins plus a contribution of 16 microeV due to the spin-orbit interaction of the conduction electrons. The 4f spin polarizes the conduction electrons via exchange interaction, which transfers the magnetic anisotropy of the conduction electrons to the 4f spin.

Van der Waals density functional for layered structures.

To understand sparse systems, we must account for both strong local atom bonds and weak nonlocal van der Waals forces between atoms separated by empty space. A fully nonlocal functional form [Phys. Rev. B 62, 6997 (2000)]] of density-functional theory (DFT) is applied here to the layered systems graphite, boron nitride, and molybdenum sulfide to compute bond lengths, binding energies, and compressibilities. These key examples show that the DFT with the generalized-gradient approximation does not apply for calculating properties of sparse matter, while use of the fully nonlocal version appears to be one way to proceed.

Metal-nonmetal transition in the boron group elements.

Structural competition in boron group elements has been studied by means of ab initio calculations. For boron we predict a structural change alpha-B-->alpha-Ga accompanied by a nonmetal-metal transition at a pressure of about 74 GPa. For Al and Ga we find an icosahedron based elemental modification (alpha-B) 0.22 and 0.05 eV/atom, respectively, higher in energy than the corresponding metallic ground state structures. In particular, the low energy difference for Ga raises expectations into the experimental feasibility of new modifications for these elements, especially in nanosized systems.