Modelling atomic and nanoscale structure in the silicon-oxygen system through active machine learning.

Silicon-oxygen compounds are among the most important ones in the natural sciences, occurring as building blocks in minerals and being used in semiconductors and catalysis. Beyond the well-known silicon dioxide, there are phases with different stoichiometric composition and nanostructured composites. One of the key challenges in understanding the Si-O system is therefore to accurately account for its nanoscale heterogeneity beyond the length scale of individual atoms. Here we show that a unified computational description of the full Si-O system is indeed possible, based on atomistic machine learning coupled to an active-learning workflow. We showcase applications to very-high-pressure silica, to surfaces and aerogels, and to the structure of amorphous silicon monoxide. In a wider context, our work illustrates how structural complexity in functional materials beyond the atomic and few-nanometre length scales can be captured with active machine learning.

Tailoring high-energy storage NaNbO3-based materials from antiferroelectric to relaxor states.

Reversible field-induced phase transitions define antiferroelectric perovskite oxides and lay the foundation for high-energy storage density materials, required for future green technologies. However, promising new antiferroelectrics are hampered by transition´s irreversibility and low electrical resistivity. Here, we demonstrate an approach to overcome these problems by adjusting the local structure and defect chemistry, delivering NaNbO3-based antiferroelectrics with well-defined double polarization loops. The attending reversible phase transition and structural changes at different length scales are probed by in situ high-energy X-ray diffraction, total scattering, transmission electron microcopy, and nuclear magnetic resonance spectroscopy. We show that the energy-storage density of the antiferroelectric compositions can be increased by an order of magnitude, while increasing the chemical disorder transforms the material to a relaxor state with a high energy efficiency of 90%. The results provide guidelines for efficient design of (anti-)ferroelectrics and open the way for the development of new material systems for a sustainable future.

SiCO Ceramics as Storage Materials for Alkali Metals/Ions: Insights on Structure Moieties from Solid-State NMR and DFT Calculations.

Polymer-derived silicon oxycarbide ceramics (SiCO) have been considered as potential anode materials for lithium- and sodium-ion batteries. To understand their electrochemical storage behavior, detailed insights into structural sites present in SiCO are required. In this work, the study of local structures in SiCO ceramics containing different amounts of carbon is presented. 13 C and 29 Si solid-state MAS NMR spectroscopy combined with DFT calculations, atomistic modeling, and EPR investigations, suggest significant changes in the local structures of SiCO ceramics even by small changes in the material composition. The provided findings on SiCO structures will contribute to the research field of polymer-derived ceramics, especially to understand electrochemical storage processes of alkali metal/ions such as Na/Na+ inside such networks in the future.

The origin of jerky dislocation motion in high-entropy alloys.

Dislocations in single-phase concentrated random alloys, including high-entropy alloys (HEAs), repeatedly encounter pinning during glide, resulting in jerky dislocation motion. While solute-dislocation interaction is well understood in conventional alloys, the origin of individual pinning points in concentrated random alloys is a matter of debate. In this work, we investigate the origin of dislocation pinning in the CoCrFeMnNi HEA. In-situ transmission electron microscopy studies reveal wavy dislocation lines and a jagged glide motion under external loading, even though no segregation or clustering is found around Shockley partial dislocations. Atomistic simulations reproduce the jerky dislocation motion and link the repeated pinning to local fluctuations in the Peierls friction. We demonstrate that the density of high local Peierls friction is proportional to the critical stress required for dislocation glide and the dislocation mobility.

Dislocation-toughened ceramics.

Functional and structural ceramics have become irreplaceable in countless high-tech applications. However, their inherent brittleness tremendously limits the application range and, despite extensive research efforts, particularly short cracks are hard to combat. While local plasticity carried by mobile dislocations allows desirable toughness in metals, high bond strength is widely believed to hinder dislocation-based toughening of ceramics. Here, we demonstrate the possibility to induce and engineer a dislocation microstructure in ceramics that improves the crack tip toughness even though such toughening does not occur naturally after conventional processing. With modern microscopy and simulation techniques, we reveal key ingredients for successful engineering of dislocation-based toughness at ambient temperature. For many ceramics a dislocation-based plastic zone is not impossible due to some intrinsic property (e.g. bond strength) but limited by an engineerable quantity, i.e. the dislocation density. The impact of dislocation density is demonstrated in a surface near region and suggested to be transferrable to bulk ceramics. Unexpected potential in improving mechanical performance of ceramics could be realized with novel synthesis strategies.

Influence of Br-/S2- site-exchange on Li diffusion mechanism in Li6PS5Br: a computational study.

We investigate how low degrees of [Formula: see text] site-exchange influence the [Formula: see text] diffusion in the argyrodite-type solid electrolyte [Formula: see text] by ab initio molecular dynamics simulations. Based on the atomic trajectories of the defect-free material, a new mechanism for the internal [Formula: see text] reorganization within the [Formula: see text] cages around the [Formula: see text] sites is identified. This reorganization mechanism is highly concerted and cannot be described by just one rotation axis. Simulations with [Formula: see text] defects reveal that [Formula: see text] interstitials ([Formula: see text]) are the dominant mobile charge carriers and originate from Frenkel pairs. These are formed because [Formula: see text] defects on the [Formula: see text] sites donate one or even two [Formula: see text] to the neighbouring cages. The [Formula: see text] then carry out intercage jumps via interstitial and interstitialcy mechanisms. With that, one single [Formula: see text] defect enables [Formula: see text] diffusion over an extended spatial area explaining why low degrees of site-exchange are sufficient to trigger superionic conduction. The vacant sites of the Frenkel pairs, namely [Formula: see text], are mostly immobile and bound to the [Formula: see text] defect. Because [Formula: see text] defects on [Formula: see text] sites act as sinks for [Formula: see text] they seem to be beneficial only for the local [Formula: see text] transport. In their vicinity T4 tetrahedral sites start to get occupied. Because the [Formula: see text] transport was found to be rather confined if [Formula: see text] and [Formula: see text] defects are direct neighbours, their relative arrangement seems to be crucial for effective long-range transport. This article is part of the Theo Murphy meeting issue 'Understanding fast-ion conduction in solid electrolytes'.

On the Origin of Zero Interface Resistance in the Li6.25Al0.25La3Zr2O12|Li0 System: An Atomistic Investigation.

Understanding the nature of ion transfer at the interface between Li metal and solid electrolytes (SE) is essential for further optimization of all-solid-state Li-ion batteries. Thus, the Li transfer across the SE|Li metal interface is investigated by means of ab initio calculations based on density functional theory in this work. The aluminum-doped garnet Li6.25Al0.25La3Zr2O12 (LLZO) is considered as a model SE due to its practical stability against Li metal. A low-energy interface model in bicrystal geometry is constructed and investigated by nudged elastic band calculations as well as ab initio molecular dynamics (AIMD) simulations. In order to distinguish between interface and bulk transport in the AIMD simulations, a post-processing protocol is developed. We find that the activation energies and diffusivities of Li are comparable in bulk LLZO and across the interface, substantiating that the interface kinetics are not rate-limiting. Moreover, electronic structure analysis indicates that charge transfer occurs gradually. Finally, Al3+ loss of LLZO at the interface rationalizes the experimentally observed phase transition from cubic to tetragonal observed close to Li metal contacts.

Experimental and theoretical investigation of the chemical exfoliation of Cr-based MAX phase particles.

Two-dimensional carbides/nitrides, typically called MXenes, are an emerging member of the ever-growing family of two-dimensional materials. The prediction of a ferromagnetic groundstate in chromium-containing MXenes has triggered growing interest in their chemical exfoliation from Cr-based MAX phases. However, the exfoliation poses serious difficulties using standard etching agents such as hydrofluoric acid (HF). Here, we investigate the exfoliability of Cr2GaC particles by chemical etching with aqueous HF both experimentally and theoretically. Structural and microstructural analyses show that the Cr2GaC particles decompose into chromium carbide and oxide without the formation of a Cr-based MXene. A thermodynamic analysis based on ab initio electronic structure calculations reveals that the exfoliation of Cr-based MXene from Cr2GaC by HF-etching is inhibited by more favorable competing reactions. This result confirms the experimental finding and suggests that HF is an unsuitable etching agent for a successful exfoliation of Cr2GaC.

Secondary-Phase-Assisted Grain Boundary Migration in CuInSe_{2}.

Significant structural evolution occurs during the deposition of CuInSe_{2} solar materials when the Cu content increases. We use in situ heating in a scanning transmission electron microscope to directly observe how grain boundaries migrate during heating, causing nondefected grains to consume highly defected grains. Cu substitutes for In in the near grain boundary regions, turning them into a Cu-Se phase topotactic with the CuInSe_{2} grain interiors. Together with density functional theory and molecular dynamics calculations, we reveal how this Cu-Se phase makes the grain boundaries highly mobile.

Reactivity of Isocyanate-Functionalized Lignins: A Key Factor for the Preparation of Lignin-Based Polyurethanes.

Using isocyanate-functionalized Kraft lignin as a reactive macromonomer for the preparation of polyurethane foams by a prepolymer technique is a well-known strategy to incorporate the biomacromolecule into a higher value polymer material. However, as of today the mechanical properties of the resulting materials are still insufficient for a number of possible applications. One reason for this limitation is that the reaction pathway and the morphological arrangement of such foams is of uttermost complexity and depends on a large number of influencing material-intrinsic factors. One important parameter is the reactivity of the functionalized lignin, which has a great impact on the interphase reaction kinetics and thus, on the geometry and mechanical properties of the resulting polyurethane foams. The reactivity is implied, amongst others, by the electron affinity of the isocyanate moiety. Herein, we investigate the reactivity of Kraft lignin modified with different commercially used isocyanates in the reaction with conventional polyols. Therefore, differently reactive prepolymers were synthesized, characterized and polyurethane foams were prepared thereof by using these compounds and the foam formation kinetics, morphological as well as mechanical properties were investigated. Finally, the results were supported by quantum mechanical calculations of the electron affinities of representative model compounds for the lignin-based prepolymers. This work gives rise to a better understanding of the effect of the reactivity and isocyanate structure linked to Kraft lignin on the polyurethane formation and enables rational choice of the isocyanate for pre-functionalization of lignin to prepare materials with better mechanical performance.

Origins of the Inverse Electrocaloric Effect.

The occurrence of the inverse (or negative) electrocaloric effect, where the isothermal application of an electric field leads to an increase in entropy and the removal of the field decreases the entropy of the system under consideration, is discussed and analyzed. Inverse electrocaloric effects have been reported to occur in several cases, for example, at transitions between ferroelectric phases with different polarization directions, in materials with certain polar defect configurations, and in antiferroelectrics. This counterintuitive relationship between entropy and applied field is intriguing and thus of general scientific interest. The combined application of normal and inverse effects has also been suggested as a means to achieve larger temperature differences between hot and cold reservoirs in future cooling devices. A good general understanding and the possibility to engineer inverse caloric effects in terms of temperature spans, required fields, and operating temperatures are thus of fundamental as well as technological importance. Here, the known cases of inverse electrocaloric effects are reviewed, their physical origins are discussed, and the different cases are compared to identify common aspects as well as potential differences. In all cases the inverse electrocaloric effect is related to the presence of competing phases or states that are close in energy and can easily be transformed with the applied field.

Highly Porous Silicon Embedded in a Ceramic Matrix: A Stable High-Capacity Electrode for Li-Ion Batteries.

We demonstrate a cost-effective synthesis route that provides Si-based anode materials with capacities between 2000 and 3000 mAh·gSi-1 (400 and 600 mAh·gcomposite-1), Coulombic efficiencies above 99.5%, and almost 100% capacity retention over more than 100 cycles. The Si-based composite is prepared from highly porous silicon (obtained by reduction of silica) by encapsulation in an organic carbon and polymer-derived silicon oxycarbide (C/SiOC) matrix. Molecular dynamics simulations show that the highly porous silicon morphology delivers free volume for the accommodation of strain leading to no macroscopic changes during initial Li-Si alloying. In addition, a carbon layer provides an electrical contact, whereas the SiOC matrix significantly diminishes the interface between the electrolyte and the electrode material and thus suppresses the formation of a solid-electrolyte interphase on Si. Electrochemical tests of the micrometer-sized, glass-fiber-derived silicon demonstrate the up-scaling potential of the presented approach.

Reconciling Local Structure Disorder and the Relaxor State in (Bi1/2Na1/2)TiO3-BaTiO3.

Lead-based relaxor ferroelectrics are key functional materials indispensable for the production of multilayer ceramic capacitors and piezoelectric transducers. Currently there are strong efforts to develop novel environmentally benign lead-free relaxor materials. The structural origins of the relaxor state and the role of composition modifications in these lead-free materials are still not well understood. In the present contribution, the solid-solution (100-x)(Bi1/2Na1/2)TiO3-xBaTiO3 (BNT-xBT), a prototypic lead-free relaxor is studied by the combination of solid-state nuclear magnetic resonance (NMR) spectroscopy, dielectric measurements and ab-initio density functional theory (DFT). For the first time it is shown that the peculiar composition dependence of the EFG distribution width (ΔQISwidth) correlates strongly to the dispersion in dielectric permittivity, a fingerprint of the relaxor state. Significant disorder is found in the local structure of BNT-xBT, as indicated by the analysis of the electric field gradient (EFG) in (23)Na 3QMAS NMR spectra. Aided by DFT calculations, this disorder is attributed to a continuous unimodal distribution of octahedral tilting. These results contrast strongly to the previously proposed coexistence of two octahedral tilt systems in BNT-xBT. Based on these results, we propose that considerable octahedral tilt disorder may be a general feature of these oxides and essential for their relaxor properties.

Influence of grain size and composition, topology and excess free volume on the deformation behavior of Cu-Zr nanoglasses.

The influence of grain size and composition on the mechanical properties of Cu-Zr nanoglasses (NGs) is investigated by molecular dynamics simulations using two model glasses of different alloy composition, namely Cu64Zr36 (Cu-rich) and Cu36Zr64 (Zr-rich). When the grain size is increased, or the fraction of interfaces in these NGs is decreased, we find a transition from a homogeneous to an inhomogeneous plastic deformation, because the softer interfaces are promoting the formation shear transformation zones. In case of the Cu-rich system, shear localization at the interfaces is most pronounced, since both the topological order and free volume content of the interfaces are very different from the bulk phase. After thermal treatment the redistribution of free volume leads to a more homogenous deformation behavior. The deformation behavior of the softer Zr-rich nanoglass, in contrast, is only weakly affected by the presence of glass-glass interfaces, since the interfaces don't show topological disorder. Our results provide clear evidence that the mechanical properties of metallic NGs can be systematically tuned by controlling the size and the chemical composition of the glassy nanograins.

Low temperature heat capacity of a severely deformed metallic glass.

The low temperature heat capacity of amorphous materials reveals a low-frequency enhancement (boson peak) of the vibrational density of states, as compared with the Debye law. By measuring the low-temperature heat capacity of a Zr-based bulk metallic glass relative to a crystalline reference state, we show that the heat capacity of the glass is strongly enhanced after severe plastic deformation by high-pressure torsion, while subsequent thermal annealing at elevated temperatures leads to a significant reduction. The detailed analysis of corresponding molecular dynamics simulations of an amorphous Zr-Cu glass shows that the change in heat capacity is primarily due to enhanced low-frequency modes within the shear band region.

Plasticity of nanocrystalline alloys with chemical order: on the strength and ductility of nanocrystalline Ni-Fe.

Plastic deformation and alloying of nanocrystalline Ni-Fe is studied by means of atomic scale computer simulations. By using a combination of Monte-Carlo and molecular dynamics methods we find that solutes have an ordering tendency even if grain sizes are in the nanometer regime, where the phase field of the ordered state is widened as compared to larger grain sizes. Tensile testing of disordered structures with various elemental distributions and the simultaneous analysis of intragranular defects reveal that solid solution strengthening is absent for the studied grain sizes. The composition and relaxation state of the grain boundary control the strength of the material, which is also found for ordered structures (L12), where dislocation activity is suppressed.

Deformation-induced grain growth and twinning in nanocrystalline palladium thin films.

The microstructure and mechanical properties of nanocrystalline Pd films prepared by magnetron sputtering have been investigated as a function of strain. The films were deposited onto polyimide substrates and tested in tensile mode. In order to follow the deformation processes in the material, several samples were strained to defined straining states, up to a maximum engineering strain of 10%, and prepared for post-mortem analysis. The nanocrystalline structure was investigated by quantitative automated crystal orientation mapping (ACOM) in a transmission electron microscope (TEM), identifying grain growth and twinning/detwinning resulting from dislocation activity as two of the mechanisms contributing to the macroscopic deformation. Depending on the initial twin density, the samples behaved differently. For low initial twin densities, an increasing twin density was found during straining. On the other hand, starting from a higher twin density, the twins were depleted with increasing strain. The findings from ACOM-TEM were confirmed by results from molecular dynamics (MD) simulations and from conventional and in-situ synchrotron X-ray diffraction (CXRD, SXRD) experiments.

Plasticity of Cu nanoparticles: Dislocation-dendrite-induced strain hardening and a limit for displacive plasticity.

The plastic behaviour of individual Cu crystallites under nanoextrusion is studied by molecular dynamics simulations. Single-crystal Cu fcc nanoparticles are embedded in a spherical force field mimicking the effect of a contracting carbon shell, inducing pressure on the system in the range of gigapascals. The material is extruded from a hole of 1.1-1.6 nm radius under athermal conditions. Simultaneous nucleation of partial dislocations at the extrusion orifice leads to the formation of dislocation dendrites in the particle causing strain hardening and high flow stress of the material. As the extrusion orifice radius is reduced below 1.3 Å we observe a transition from displacive plasticity to solid-state amorphisation.

Size-dependent lattice expansion in nanoparticles: reality or anomaly?

Size-dependent lattice expansion of nanoparticles is observed for many ionic compounds, including metal oxides, while lattice contraction prevails for pure metals. However, the physical origin of this effect, which is of importance for the thermodynamic, chemical and electronic properties of nanoparticles, is discussed controversially. After a survey of the experimental literature, revealing a wide variety of materials with size-dependent lattice expansion, we show that the negative surface stress is the key reason for lattice expansion, while the excess of lattice sums or point defects of various charge states can be excluded as general explanations. Ab initio calculations of surface stresses for various surface structures of metal oxides confirm the model of a surface-induced lattice expansion.

Size-dependent phase diagrams of metallic alloys: A Monte Carlo simulation study on order-disorder transitions in Pt-Rh nanoparticles.

Nanoparticles of Pt-Rh were studied by means of lattice-based Monte Carlo simulations with respect to the stability of ordered D0(22)- and 40-phases as a function of particle size and composition. By thermodynamic integration in the semi-grand canonical ensemble, phase diagrams for particles with a diameter of 7.8 nm, 4.3 nm and 3.1 nm were obtained. Size-dependent trends such as the lowering of the critical ordering temperature, the broadening of the compositional stability range of the ordered phases, and the narrowing of the two-phase regions were observed and discussed in the context of complete size-dependent nanoparticle phase diagrams. In addition, an ordered surface phase emerges at low temperatures and low platinum concentration. A decrease of platinum surface segregation with increasing global platinum concentration was observed, when a second, ordered phase is formed inside the core of the particle. The order-disorder transitions were analyzed in terms of the Warren-Cowley short-range order parameters. Concentration-averaged short-range order parameters were used to remove the surface segregation bias of the conventional short-range order parameters. Using this procedure, it was shown that the short-range order in the particles at high temperatures is bulk-like.