Computational Investigation of Carbon Based Anode Materials for Li- and Post-Li- Ion Batteries

Due to its negligible capacity with respect to sodium intercalation, graphite is not suited as anode material for sodium ion batteries. Hard carbon materials, on the other hand, provide reasonably high capacities at low insertion potential, making them a promising anode materials for sodium (and potassium) ion batteries. The particular nanostructure of these functionalized carbon-based materials has been found to be crucially linked to the material performance. However, there is still a lack of understanding with respect to the functional role of structural units, such as defects, for intercalation and storage. To overcome these problems, the intercalation of Li, Na, and K in graphitic model structures with distinct defect configurations has been investigated by density functional theory. The calculations confirm that defects are able to stabilize intercalation of larger alkali metal contents. At the same time, it is shown that a combination of phonon and band structure calculations are able to explain characteristic Raman features typically observed for alkali metal intercalation in hard carbon, furthermore allowing for the quantification of the alkali metal intercalation inbetween the layers of hard carbon anodes.

The relevance of structural variability in the time-domain for computational reflection anisotropy spectroscopy at solid-liquid interfaces

In electrochemistry, reactions and charge-transfer are to a large extent determined by the atomistic structure of the solid-liquid interface. Yet due to the presence of the liquid electrolyte, many surface-science methods cannot be applied here. Hence, the exact microscopic structure that is present under operating conditions often remains unknown. Reflection anisotropy spectroscopy (RAS) is one of the few techniques that allow for anin operandoinvestigation of the structure of solid-liquid interfaces. However, an interpretation of RAS data on the atomistic scale can only be obtained by comparison to computational spectroscopy. While the number of computational RAS studies related to electrochemical systems is currently still limited, those studies so far have not taken into account the dynamic nature of the solid-liquid interface. In this work, we investigate the temporal evolution of the spectroscopic response of the Au(110) missing row reconstruction in contact with water by combiningab initiomolecular dynamics with computational spectroscopy. Our results show significant changes in the time evolution of the RA spectra, in particular providing an explanation for the typically observed differences in intensity when comparing theory and experiment. Moreover, these findings point to the importance of structural surface/interface variability while at the same time emphasising the potential of RAS for probing these dynamic interfaces.

Zn Metal Atom Doping on the Surface Plane of One-Dimesional NiMoO4 Nanorods with Improved Redox Chemistry

The effect of zinc (Zn) doping and defect formation on the surface of nickel molybdate (NiMoO₄) structures with varying Zn content has been studied to produce one-dimensional electrodes and catalysts for electrochemical energy storage and ethanol oxidation, respectively. Zn-doped nickel molybdate (Ni_1-xZn_xMoO₄, where x = 0.1, 0.2, 0.4, and 0.6) nanorods were synthesized by a simple wet chemical route. The optimal amount of Zn is found to be around 0.25 above which the NiMoO₄ becomes unstable, resulting in poor electrochemical activity. This result agrees with our density functional theory calculations in which the thermodynamic stability reveals that Ni_1-xZn_xMoO₄ crystallized in the β-NiMoO₄ phase and is found to be stable for x≤0.25. Analytical techniques show direct evidence of the presence of Zn in the NiMoO₄ nanorods, which subtly alter the electrocatalytic activity. Compared with pristine NiMoO₄, Zn-doped NiMoO₄ with the optimized Zn content was tested as an electrode for an asymmetric supercapacitor and demonstrated an enhanced specific capacitance of 122 F g^-1 with a high specific energy density of 43 W h kg^-1 at a high power density of 384 W kg^-1. Our calculations suggest that the good conductivity from Zn doping is attributed to the formation of excess oxygen vacancies and dopants play an important role in enhancing the charge transfer between the surface and OH^- ions from the electrolyte. We report electrochemical testing, material characterization, and computational insights and demonstrate that the appropriate amount of Zn in NiMoO₄ can improve the storage capacity (∼15%) due to oxygen vacancy interactions.

Operando pH Measurements Decipher H+/Zn2+ Intercalation Chemistry in High-Performance Aqueous Zn/δ-V2O5 Batteries

Vanadium oxides have been recognized to be among the most promising positive electrode materials for aqueous zinc metal batteries (AZMBs). However, their underlying intercalation mechanisms are still vigorously debated. To shed light on the intercalation mechanisms, high-performance δ-V₂O₅ is investigated as a model compound. Its structural and electrochemical behaviors in the designed cells with three different electrolytes, i.e., 3 m Zn(CF₃SO₃)₂/water, 0.01 M H₂SO₄/water, and 1 M Zn(CF₃SO₃)₂/acetonitrile, demonstrate that the conventional structural and elemental characterization methods cannot adequately clarify the separate roles of H⁺ and Zn²⁺ intercalations in the Zn(CF₃SO₃)₂/water electrolyte. Thus, an operando pH determination method is developed and used toward Zn/δ-V₂O₅ AZMBs. This method indicates the intercalation of both H⁺ and Zn²⁺ into δ-V₂O₅ and uncovers an unusual H⁺/Zn²⁺-exchange intercalation-deintercalation mechanism. Density functional theory calculations further reveal that the H⁺/Zn²⁺ intercalation chemistry is a consequence of the variation of the electrochemical potential of Zn²⁺ and H⁺ during the electrochemical intercalation/release.

Design and Tuning of the Electrochemical Properties of Vanadium-Based Cation-Disordered Rock-Salt Oxide Positive Electrode Material for Lithium-Ion Batteries

Disordered rock-salt compounds are becoming increasingly important due to their potential as high-capacity positive electrode materials for lithium-ion batteries. Thereby, a significant number of studies have focused on increasing the accessible Li capacity, but studies to manipulate the electrochemical potential are limited. This work explores the effect of transition-metal substitution on the electrochemistry of ternary disordered rock-salt-type compounds with LiM²⁺_0.5V_0.5⁴⁺O₂ stoichiometry (M = Mn, Fe, Co) directly synthesized through mechanochemistry. Rietveld refinements of synchrotron X-ray diffraction patterns confirm the disordered rock-salt structures. First-principles density functional theory study is used to predict the impact of the cation substitution on the expected average voltage and the electronic structures of these materials are used to analyze the underlying redox processes. For LiM²⁺_0.5V⁴⁺_0.5O₂ (M = Mn, Fe, Co), discharge voltages increase in the order of Mn < Fe < Co with 2.28, 2.41, and 2.51 V, exhibiting discharge capacities of 219, 207, and 234 mAh g^-1, respectively. In comparison, for the disordered rock-salt Li₂VO_3, an average discharge voltage of ∼2.2 V with V^5+/4+ redox couple has been reported. However, detrimental electrode-electrolyte interactions manifested as transition-metal dissolution has been found to result in severe capacity fading. Thereto, the use of a concentrated 5.5 M LiFSI increased the cycling stability significantly, effectively reducing transition-metal dissolution. The underlying reasons for the capacity fading of disordered rock salts are yet unclear. We stress the importance of cathode-electrolyte interactions, thus opening new directions for the improvement of cation-disordered materials.

Design of Nickel-Based Cation-Disordered Rock-Salt Oxides: The Effect of Transition Metal (M = V, Ti, Zr) Substitution in LiNi0.5M0.5O2 Binary Systems

Cation-disordered oxides have been ignored as positive electrode material for a long time due to structurally limited lithium insertion/extraction capabilities. In this work, a case study is carried out on nickel-based cation-disordered Fm3 ̅m LiNi_0.5M_0.5O₂ positive electrode materials. The present investigation targets tailoring the electrochemical properties for nickel-based cation-disordered rock-salt by electronic considerations. The compositional space for binary LiM⁺³O₂ with metals active for +3/+4 redox couples is extended to ternary oxides with LiA_0.5B_0.5O₂ with A = Ni²⁺ and B = Ti⁴⁺, Zr⁴⁺, and V⁺⁴ to assess the impact of the different transition metals in the isostructural oxides. The direct synthesis of various new unknown ternary nickel-based Fm3̅ m cation-disordered rock-salt positive electrode materials is presented with a particular focus on the LiNi_0.5V_0.5O₂ system. This positive electrode material for Li-ion batteries displays an average voltage of ∼2.55 V and a high discharge capacity of 264 mAhg^-1 corresponding to 0.94 Li. For appropriate cutoff voltages, a long cycle life is achieved. The charge compensation mechanism is probed by XANES, confirming the reversible oxidation and reduction of V⁴⁺/V⁵⁺. The enhancement in the electrochemical performances within the presented compounds stresses the importance of mixed cation-disordered transition metal oxides with different electronic configuration.

Direct measurement of individual phonon lifetimes in the clathrate compound Ba7.81Ge40.67Au5.33

Engineering lattice thermal conductivity requires to control the heat carried by atomic vibration waves, the phonons. The key parameter for quantifying it is the phonon lifetime, limiting the travelling distance, whose determination is however at the limits of instrumental capabilities. Here, we show the achievement of a direct quantitative measurement of phonon lifetimes in a single crystal of the clathrate Ba_7.81Ge_40.67Au_5.33, renowned for its puzzling ‘glass-like’ thermal conductivity. Surprisingly, thermal transport is dominated by acoustic phonons with long lifetimes, travelling over distances of 10 to 100 nm as their wave-vector goes from 0.3 to 0.1 Å⁻¹. Considering only low-energy acoustic phonons, and their observed lifetime, leads to a calculated thermal conductivity very close to the experimental one. Our results challenge the current picture of thermal transport in clathrates, underlining the inability of state-of-the-art simulations to reproduce the experimental data, thus representing a crucial experimental input for theoretical developments.

Phonon lifetime is a fundamental parameter of thermal transport however its determination is challenging. Using inelastic neutron scattering and the neutron resonant spin-echo technique, Lory et al. determine the acoustic phonon lifetime in a single crystal of clathrate Ba7.81Ge40.67Au5.33.

Atomic structure and phason modes of the Sc-Zn icosahedral quasicrystal

The detailed atomic structure of the binary icosahedral (i) ScZn7.33 quasicrystal has been investigated by means of high-resolution synchrotron single-crystal X-ray diffraction and absolute scale measurements of diffuse scattering. The average atomic structure has been solved using the measured Bragg intensity data based on a six-dimensional model that is isostructural to the i-YbCd5.7 one. The structure is described with a quasiperiodic packing of large Tsai-type rhombic triacontahedron clusters and double Friauf polyhedra (DFP), both resulting from a close-packing of a large (Sc) and a small (Zn) atom. The difference in chemical composition between i-ScZn7.33 and i-YbCd5.7 was found to lie in the icosahedron shell and the DFP where in i-ScZn7.33 chemical disorder occurs on the large atom sites, which induces a significant distortion to the structure units. The intensity in reciprocal space displays a substantial amount of diffuse scattering with anisotropic distribution, located around the strong Bragg peaks, that can be fully interpreted as resulting from phason fluctuations, with a ratio of the phason elastic constants K 2/K 1 = -0.53, i.e. close to a threefold instability limit. This induces a relatively large perpendicular (or phason) Debye-Waller factor, which explains the vanishing of 'high-Q perp' reflections.

Vacancy-dependent stability of cubic and wurtzite Ti1-x Al x N

While it is well-known that supersaturated cubic-structured Ti_1-x Al _x N can be prepared by physical vapor deposition, the impact of point defects on formation process and cubic to wurtzite transition is largely unexplored. Irrespective of point defects, ab initio calculations correctly predict the Al concentration of the cubic to wurtzite transition. By means of density functional theory we show that vacancies on metal and/or non-metal sites only slightly affect the cubic to wurtzite transition region, whereas they clearly affect the physical properties.

Designing thin film materials - Ternary borides from first principles

Exploiting the mechanisms responsible for the exceptional properties of aluminum based nitride coatings, we apply ab initio calculations to develop a recipe for designing functional thin film materials based on ternary diborides. The combination of binary diborides, preferring different structure types, results in supersaturated metastable ternary systems with potential for phase transformation induced effects. For the exemplary cases of M _x W_1 - x B₂ (with M = Al, Ti, V) we show by detailed ab initio calculations that the respective ternary solid solutions are likely to be experimentally accessible by modern depositions techniques.

Circular test structure for the determination of piezoelectric constants of Sc x Al1-x N thin films applying Laser Doppler Vibrometry and FEM simulations

Piezoelectric scandium aluminium nitride (Sc _x Al_1-x N) thin films offer a large potential for the application in micro electromechanical systems, as advantageous properties of pure AlN thin films are maintained, but combined with an increased piezoelectric actuation and sensing potential. Sc _x Al_1-x N thin films with x = 27% have been prepared by DC reactive magnetron sputtering to find optimized deposition parameters to maximize the piezoelectric constants d₃₃ and d₃₁. For the accurate and simultaneous measurement of these constants Laser Doppler Vibrometry has been applied and compared to finite element (FEM) simulations. The electrode design has been optimized to rotational symmetric structures enabling a 180° phase shifted excitation, so that a straight-forward comparison of experimental displacement curves with those obtained from FEM is feasible.

Complementary ab initio and X-ray nanodiffraction studies of Ta2O5

The complex structure of Ta₂O₅ led to the development of various structural models. Among them, superstructures represent the most stable configurations. However, their formation requires kinetic activity and long-range ordering processes, which are hardly present during physical vapor deposition. Based on nano-beam X-ray diffraction and concomitant ab initio studies, a new metastable orthorhombic basic structure is introduced for Ta₂O₅ with lattice parameters a = 6.425 Å, b = 3.769 Å and c = 7.706 Å. The unit cell containing only 14 atoms, i.e. two formula unit blocks in the c direction, is characterized by periodically alternating the occupied oxygen site between two possible positions in succeeding 002-planes. This structure can be described by the space group 53 (Pncm) with four Wyckoff positions, and exhibits an energy of formation of -3.209 eV atom^-1. Among all the reported basic structures, its energy of formation is closest to those of superstructures. Furthermore, this model exhibits a 2.5 eV band gap, which is closer to experimental data than the band gap of any other basic-structure model. The sputtered Ta₂O₅ films develop only a superstructure if annealed at temperatures >800 °C in air or vacuum. Based on these results and the conveniently small unit cell size, it is proposed that the basic-structure model described here is an ideal candidate for both structure and electronic state descriptions of orthorhombic Ta₂O₅ materials.

Localization of propagative phonons in a perfectly crystalline solid

Perfectly crystalline solids are excellent heat conductors. Prominent counterexamples are intermetallic clathrates, guest-host systems with a high potential for thermoelectric applications due to their ultralow thermal conductivities. Our combined experimental and theoretical investigation of the lattice dynamics of a particularly simple binary representative, Ba(8)Si(46), identifies the mechanism responsible for the reduction of lattice thermal conductivity intrinsic to the perfect crystal structure. Above a critical wave vector, the purely harmonic guest-host interaction leads to a drastic transfer of spectral weight to the guest atoms, corresponding to a localization of the propagative phonons.

Lattice dynamics of the icosahedral quasicrystals i-ZnMgSc and i-ZnAgSc and the cubic 1/1-approximant Zn6Sc

A comparison of periodic approximants and their quasicrystalline counterparts offers the opportunity to better understand the structure, physical properties and stabilizing mechanisms of these complex phases. We present a combined experimental and computational study of the lattice dynamics of the icosahedral quasicrystals i-ZnMgSc and i-ZnAgSc and compare these to the lattice dynamics of the cubic 1/1-approximant Zn6Sc. The two phases, quasicrystal and approximant, are built up from the same atomic clusters, which are packed either quasiperiodically or on a body centered cubic lattice, respectively. Using inelastic neutron scattering and atomic scale simulations, we show that the vibrational spectra of these three systems are very similar, however, they contain a clear signature of the increasing structural complexity from approximant to quasicrystal.

Short- and long-range ordering during the phase transition of the Zn6Sc 1/1 cubic approximant

Using in situ x-ray scattering and synchrotron radiation, we have experimentally elucidated the mechanism of the cubic to monoclinic phase transition in the Zn6Sc 1/1 approximant to an icosahedral quasicrystal. The high-temperature cubic phase is described as a bcc packing of a large Tsai-type icosahedral cluster whose center is occupied by an orientationally disordered Zn4 tetrahedron. A clear monoclinic distortion has been found to take place within 2 K around Tc = 157 K, in excellent agreement with the observed anomalies in the electrical resistivity and heat capacity. Also, a rapid variation of the super-structure reflection intensity is observed. The low-temperature monoclinic phase, as determined by single-crystal x-ray diffraction at 40 K, has been confirmed to consist of ordered Zn4 tetrahedra, oriented in an anti-parallel way along the [[Formula: see text]] direction. Above Tc, a diffuse scattering signal is observed at the position of the super-structure reflections, which evidences that a short-range ordering of the Zn4 tetrahedra takes place. In a way similar to a second-order phase transition, the correlation length describing this short-range ordering increases rapidly when the temperature diminishes and almost diverges when the temperature is close to Tc, going from 200 Å at 220 K to reach the very large value of 1200 Å at 161 K. Finally, using single-crystal x-ray diffraction, the atomic structure of the low-temperature monoclinic super-structure (space group C2/c) could be solved. The ordering of the Zn4 tetrahedra is accompanied by a strong distortion of the surrounding shells.

Tetrahedron dynamics in the icosahedral quasicrystals i-ZnMgSc and i-ZnAgSc and the cubic 1/1-approximant Zn6Sc

A comparison of periodic approximants and their quasicrystalline counterparts offers the opportunity to better understand the structure, physical properties and stabilizing mechanisms of these complex phases. We present a combined experimental and molecular dynamics study of the lattice dynamics of the icosahedral quasicrystals i-ZnMgSc and i-ZnAgSc and compare it to recently published results obtained for the cubic 1/1-approximant Zn(6)Sc. Both phases, quasicrystal and approximant, are built up from large atomic clusters which contain a tetrahedral shell at the cluster centre and are packed either quasiperiodically or on a bcc lattice. Using quasielastic neutron scattering and atomic scale simulations, we show that in the quasicrystal the tetrahedra display a dynamics similar to that observed in the 1/1-approximant: the tetrahedra behave as a 'single molecule' and reorient dynamically on a timescale of the order of a few ps. The tetrahedra reorientation is accompanied by a large distortion of the surrounding cluster shells which provide a unique dynamical flexibility to the quasicrystal. However, whereas in the 1/1-approximant the tetrahedron reorientation is observed down to T(c) = 160 K, where a phase transition takes place, in the quasicrystal the tetrahedron dynamics is gradually freezing from 550 to 300 K, similarly to a glassy system.

Ordering and dynamics of the central tetrahedron in the 1/1 Zn6Sc periodic approximant to quasicrystal

Periodic approximants to quasicrystals offer a unique opportunity to better understand the structure, physical properties and stabilizing mechanisms of their quasicrystal counterparts. We present a detailed study of the order-disorder phase transition occurring at about 160 K in the Zn(6)Sc cubic approximant to the icosahedral quasicrystal i-MgZnSc. This transition goes along with an anti-parallel ordering of the tetrahedra located at the centres of large atomic clusters, which are packed on a bcc lattice. Single crystal x-ray diffuse scattering shows that the tetrahedra display pre-transitional short range ordering above T(c) (Yamada et al 2012 in preparation). Using quasielastic neutron scattering (QENS) we clearly evidence this short range order to be dynamical in nature above T(c). The QENS data are consistent with a model of tetrahedra 'jumping' between almost equivalent positions, which is supported by molecular dynamics simulations. This demonstrates a unique dynamical flexibility of the Zn(6)Sc structure even at room temperature.

Vibrational properties of MgZn2

We present here simulation results on the dynamical structure factor of the C14 Laves Phase of MgZn2, the simplest of the Mg–(Al,Zn) Frank-Kasper alloy phases. The dynamical structure factor was determined in two ways. Firstly, the dynamical matrix was obtained in harmonic approximation from ab-initio forces. The dynamical structure factor can then be computed from the eigenvalues of the dynamical matrix. Alternatively, Molecular Dynamics simulations of a larger sample were used to measure the correlation function corresponding to the dynamical structure factor. Both results are compared to data from neutron scattering experiments. This comparison also includes the intensity distribution, which is a very sensitive test. We find that the dynamical structure factor determined with either method agrees reasonably well with the experiment. In particular, the intensity transfer from acoustic to optic phonon modes can be reproduced correctly. This shows that simulation studies can complement phonon dispersion measurements.