Identification of a Durability Descriptor for Molecular Oxygen Reduction Reaction Catalysts.

The development of durable platinum-group-metal-free oxygen reduction reaction (ORR) catalysts is a key research direction for enabling the wide use of fuel cells. Here, we use a combination of experimental measurements and density functional theory calculations to study the activity and durability of seven iron-based metallophthalocyanine (MPc) ORR catalysts that differ only in the identity of the substituent groups on the MPcs. While the MPcs show similar ORR activity, their durabilities as measured by the current decay half-life differ greatly. We find that the energy difference between the hydrogenated intermediate structure and the final demetalated structure (ΔEdemetalation) of the MPcs is linearly related to the degradation reaction barrier energy. Comparison to the degradation data for the previously studied metallocorrole systems suggested that ΔEdemetalation also serves as a descriptor for the corrole systems and that the high availability of protons at the active site due to the COOH group of the o-corrole decreases the durability.

Identification of High-Reliability Regions of Machine Learning Predictions Based on Materials Chemistry.

Progress in the application of machine learning (ML) methods to materials design is hindered by the lack of understanding of the reliability of ML predictions, in particular, for the application of ML to small data sets often found in materials science. Using ML prediction for transparent conductor oxide formation energy and band gap, dilute solute diffusion, and perovskite formation energy, band gap, and lattice parameter as examples, we demonstrate that (1) construction of a convex hull in feature space that encloses accurately predicted systems can be used to identify regions in feature space for which ML predictions are highly reliable; (2) analysis of the systems enclosed by the convex hull can be used to extract physical understanding; and (3) materials that satisfy all well-known chemical and physical principles that make a material physically reasonable are likely to be similar and show strong relationships between the properties of interest and the standard features used in ML. We also show that similar to the composition-structure-property relationships, inclusion in the ML training data set of materials from classes with different chemical properties will not be beneficial for the accuracy of ML prediction and that reliable results likely will be obtained by ML model for narrow classes of similar materials even in the case where the ML model will show large errors on the data set consisting of several classes of materials.

Origin of the Rarely Reported High Performance of Mn-doped Carbon-based Oxygen Reduction Catalysts.

Recent efforts to develop durable high-performance platinum-group metal (PGM)-free oxygen reduction reaction (ORR) electrocatalysts have focused on Fe- and Co-based molecular and pyrolyzed catalysts. While Mn-based catalysts have advantages of lower toxicity and higher durability, their activity has been generally poor. Nevertheless, several examples of high-performance Mn-based catalysts have been reported. Thus, it is necessary to understand why Mn-based materials much more rarely show high catalytic ORR performance and to determine the factors that can lead to the achievement of such high performance in these rare cases. We have studied the effects of the changes in the macrocycle structure, axial ligand, distance between the active sites, interactions with the dopant N atoms and the presence of an extended carbon network on the ORR catalysis of various Mn-, Fe-, and Co-based systems through the comparison of the adsorption energies of the ORR intermediates. We find that the sensitivity to the local environment changes is the largest for Mn and is the smallest for Co, with Fe between Mn and Co. Our results showed that the strong binding of OH by Mn and the strong sensitivity of the Mn to the modification of its environment necessitate a precise combination of local environment changes to achieve a high onset potential (Vonset ) in Mn-based catalysts. By contrast, the weaker binding of OH by Fe and Co and their weaker sensitivity to local environment changes lead to a wide variety of local environments with favorable catalytic activity (Vonset >0.7 V) for Co- and Fe-based systems. This explains the scarcity of reported Mn-based pyrolyzed catalysts and suggests that precise material synthesis and engineering of the active site can achieve high-performance Mn-based ORR electrocatalysts with high activity and durability.

Application of Poly-L-Lysine for Tailoring Graphene Oxide Mediated Contact Formation between Lithium Titanium Oxide LTO Surfaces for Batteries.

When producing stable electrodes, polymeric binders are highly functional materials that are effective in dispersing lithium-based oxides such as Li4Ti5O12 (LTO) and carbon-based materials and establishing the conductivity of the multiphase composites. Nowadays, binders such as polyvinylidene fluoride (PVDF) are used, requiring dedicated recycling strategies due to their low biodegradability and use of toxic solvents to dissolve it. Better structuring of the carbon layers and a low amount of binder could reduce the number of inactive materials in the electrode. In this study, we use computational and experimental methods to explore the use of the poly amino acid poly-L-lysine (PLL) as a novel biodegradable binder that is placed directly between nanostructured LTO and reduced graphene oxide. Density functional theory (DFT) calculations allowed us to determine that the (111) surface is the most stable LTO surface exposed to lysine. We performed Kubo-Greenwood electrical conductivity (KGEC) calculations to determine the electrical conductivity values for the hybrid LTO-lysine-rGO system. We found that the presence of the lysine-based binder at the interface increased the conductivity of the interface by four-fold relative to LTO-rGO in a lysine monolayer configuration, while two-stack lysine molecules resulted in 0.3-fold (in the plane orientation) and 0.26-fold (out of plane orientation) increases. These outcomes suggest that monolayers of lysine would specifically favor the conductivity. Experimentally, the assembly of graphene oxide on poly-L-lysine-TiO2 with sputter-deposited titania as a smooth and hydrophilic model substrate was investigated using a layer-by-layer (LBL) approach to realize the required composite morphology. Characterization techniques such as X-ray photoelectron spectroscopy (XPS), atomic force microscopy (AFM), Kelvin probe force microscopy (KPFM), scanning electron microscopy (SEM) were used to characterize the formed layers. Our experimental results show that thin layers of rGO were assembled on the TiO2 using PLL. Furthermore, the PLL adsorbates decrease the work function difference between the rGO- and the non-rGO-coated surface and increased the specific discharge capacity of the LTO-rGO composite material. Further experimental studies are necessary to determine the influence of the PLL for aspects such as the solid electrolyte interface, dendrite formation, and crack formation.

Understanding hydrazine oxidation electrocatalysis on undoped carbon.

Carbons are ubiquitous electrocatalytic supports for various energy-related transformations, especially in fuel cells. Doped carbons such as Fe-N-C materials are particularly active towards the oxidation of hydrazine, an alternative fuel and hydrogen carrier. However, there is little discussion of the electrocatalytic role of the most abundant component - the carbon matrix - towards the hydrazine oxidation reaction (HzOR). We present a systematic investigation of undoped graphitic carbons towards the HzOR in alkaline electrolyte. Using highly oriented pyrolytic graphite electrodes, as well as graphite powders enriched in either basal planes or edge defects, we demonstrate that edge defects are the most active catalytic sites during hydrazine oxidation electrocatalysis. Theoretical DFT calculations support and explain the mechanism of HzOR on carbon edges, identifying unsaturated graphene armchair defects as the most likely active sites. Finally, these findings explain the 'double peak' voltammetric feature observed on many doped carbons during the HzOR.

CVD-Assisted Synthesis of 2D Layered MoSe2 on Mo Foil and Low Frequency Raman Scattering of Its Exfoliated Few-Layer Nanosheets on CaF2 Substrates.

Transition-metal dichalcogenides (TMDCs) are unique layered materials with exotic properties. So, examining their structures holds tremendous importance. 2H-MoSe2 (analogous to MoS2; Gr. 6 TMDC) is a crucial optoelectronic material studied extensively using Raman spectroscopy. In this regard, low-frequency Raman (LFR) spectroscopy can probe this material's structure as it reveals distinct vibration modes. Here, we focus on understanding the microstructural evolution of different 2H-MoSe2 morphologies and their layers using LFR scattering. We grew phase-pure 2H-MoSe2 (with variable microstructures) directly on a Mo foil using a two-furnace ambient-pressure chemical vapor deposition (CVD) system by carefully controlling the process parameters. We analyzed the layers of exfoliated flakes after ultrasonication and drop-cast 2H-MoSe2 of different layer thicknesses by choosing different concentrations of 2H-MoSe2 solutions. Further detailed analyses of the respective LFR regions confirm the presence of newly identified Raman signals for the 2H-MoSe2 nanosheets drop-cast on Raman-grade CaF2. Our results show that CaF2 is an appropriate Raman-enhancing substrate compared to Si/SiO2 as it presents new LFR modes of 2H-MoSe2. Therefore, CaF2 substrates are a promising medium to characterize in detail other TMDCs using LFR spectroscopy.

A Predictive Theory for Domain Walls in Oxide Ferroelectrics Based on Interatomic Interactions and its Implications for Collective Material Properties.

Domain walls separating regions of ferroelectric material with polarization oriented in different directions are crucial for applications of ferroelectrics. Rational design of ferroelectric materials requires the development of a theory describing how compositional and environmental changes affect domain walls. To model domain wall systems, a discrete microscopic Landau-Ginzburg-Devonshire (dmLGD) approach with A- and B-site cation displacements serving as order parameters is developed. Application of dmLGD to the classic BaTiO3 , KNbO3, and PbTiO3 ferroelectrics shows that A-B cation repulsion is the key interaction that couples the polarization in neighboring unit cells of the material. dmLGD decomposition of the total energy of the system into the contributions of the individual cations and their interactions enables the prediction of different properties for a wide range of ferroelectric perovskites based on the results obtained for BaTiO3 , KNbO3, and PbTiO3 only. It is found that the information necessary to estimate the structure and energy of domain-wall "defects" can be extracted from single-domain 5-atom first-principles calculations, and that "defect-like" domain walls offer a simple model system that sheds light on the relative stabilities of the ferroelectric, antiferroelectric, and paraelectric bulk phases. The dmLGD approach provides a general theoretical framework for understanding and designing ferroelectric perovskite oxides.

A Multilevel Analytical Theory for Prediction of Ferroelectric Perovskite Oxide Properties from Composition.

Prediction of properties from composition is a fundamental goal of materials science that is particularly relevant for ferroelectric perovskite oxide solid solutions where compositional variation is a primary tool for material design. Design of ferroelectric oxide solid solutions has been guided by heuristics and first-principles and Landau-Ginzburg-Devonshire theoretical methods that become increasingly difficult to apply in ternary, quaternary, and quintary solid solutions. To address this problem, a multilevel model is developed for the prediction of the ferroelectric-to-paraelectric transition temperature (Tc ), coercive field (Ec ), and polarization (P) of PbTiO3 -derived ferroelectric solid solutions from composition. The characteristics of the materials at different length scales, starting at the level of the electronic structure and chemical bonding of the constituent ions and ending at the level of collective behavior, are analytically related by using ferroelectric domain walls and cationic off-center displacements as the key links between the different levels of the model. The obtained composition-structure-property relationships provide a unified quantitatively predictive theory for understanding PbTiO3 -derived solid solutions. Such a multilevel analytical modeling approach is likely to be generally applicable to different classes of ferroelectric perovskite oxides and to other functional properties, and to materials and properties beyond the field of ferroelectrics.

Fe-N-C electrocatalysts in the oxygen and nitrogen cycles in alkaline media: the role of iron carbide.

Fe-N-C electrocatalysts hold a great promise for Pt-free energy conversion, driving the electrocatalysis of oxygen reduction and evolution, oxidation of nitrogen fuels, and reduction of N2, CO2, and NOx. Nevertheless, the catalytic role of iron carbide, a component of nearly every pyrolytic Fe-N-C material, is at the focus of a heated controversy. We now resolve the debate by examining a broad range of Fe3C sites, spanning across many typical size distributions and carbon environments. Removing Fe3C selectively by a non-oxidizing acid reveals its inactivity towards two representative reactions in alkaline media, oxygen reduction and hydrazine oxidation. The activity is assigned to other pre-existing sites, most probably Fe-Nx. DFT calculations prove that the Fe3C surface binds O and N intermediates too strongly to be catalytic. By settling the argument on the catalytic role of Fe3C in alkaline electrocatalysis, we hope to spur innovation in this critical field.

Modelling of high-temperature order-disorder phase transitions of non-stoichiometric Mo2C and Ti2C from first principles.

High-temperature order-disorder phase transitions play an important role in determining the structure and physical and chemical properties of non-stoichiometric transition metal carbides. Due to the large number of possible carbon vacancy arrangements, it is difficult to study these systems with first-principles calculations. Here, we construct a simple atomistic potential capable of accurately reproducing the energetics of the carbon vacancy arrangements in cubic Mo2C and Ti2C obtained from density functional theory calculations. We show that this potential can be applied to correctly predict the transition temperatures between the ordered and disordered states in Monte Carlo simulations on large supercells and reveal the extend of local order in the disordered phases of Mo2C and Ti2C that show interesting physical and chemical properties. We find that even the high-temperature disordered phase exhibit a relatively high degree of local order as indicated by the relatively small change in the root mean square number of C atom neighbours of Mo/Ti compared to the ordered phase (from 3.0 to 3.1-3.2). This atomistic potential enables the study of how the structure of these carbides can be tuned through the synthesis temperature to control the properties of carbide materials that are related to the degree of disorder in the system such as catalytic activity and electrical conductivity and play an important role in applications of these carbides. Fundamentally, the successful modelling of these carbides suggests that despite the presence of metallic, covalent and ionic interactions, bonding in carbides can be modelled by simple and physically intuitive interatomic potentials.

Sustainable existence of solid mercury (Hg) nanoparticles at room temperature and their applications.

Although liquid mercury (Hg) has been known since antiquity, the formation of stable solid nano forms of Hg at room temperature has not been reported so far. Here, for the first time, we report a simple sonochemical route to obtain solid mercury nanoparticles, stabilized by reduced graphene oxide at ambient conditions. The as-formed solid Hg nanoparticles were found to exhibit remarkable rhombohedral morphology and crystallinity at room temperature. Extensive characterization using various physicochemical techniques revealed the unique properties of the solid nanoparticles of Hg compared to its bulk liquid metal phase. Furthermore, the solid nature of the Hg nanoparticles was studied electrochemically, revealing distinctive properties. We believe that solid Hg nanoparticles have the potential for important applications in the fields of electroanalytical chemistry and electrocatalysis.

Control of Molecular Catalysts for Oxygen Reduction by Variation of pH and Functional Groups.

In the search for replacement of the platinum-based catalysts for fuel cells, MN4 molecular catalysts based on abundant transition metals play a crucial role in modeling and investigation of the influence of the environment near the active site in platinum-group metal-free (PGM-free) oxygen reduction reaction (ORR) catalysts. To understand how the ORR activity of molecular catalysts can be controlled by the active site structure through modification by the pH and substituent functional groups, the change of the ORR onset potential and the electron number in a broad pH range was examined for three different metallocorroles. Experiments revealed a switch between two different ORR mechanisms and a change from 2e- to 4e- pathway in the pH range of 3.5-6. This phenomenon was shown by density functional theory (DFT) calculations to be related to the protonation of the nitrogen atoms and carboxylic acid groups on the corroles indicated by the pKa values of the protonation sites in the vicinity of the ORR active sites. Control of the electron-withdrawing nature of these groups characterized by the pKa values could switch the ORR from the H+ to e- rate-determining step mechanisms and from 2e- to 4e- ORR pathways and also controlled the durability of the corrole catalysts. The results suggest that protonation of the nitrogen atoms plays a vital role in both the ORR activity and durability for these materials and that pKa of the N atoms at the active sites can be used as a descriptor for the design of high-performance, durable PGM-free catalysts.

Tuning of ORR activity through the stabilization of the adsorbates by hydrogen bonding with substituent groups.

Metallocorroles and metalloporphyrins (M-N-C) are some of the best alternative molecular catalysts for the replacement of the expensive platinum-group metals (PGM) in oxygen reduction reaction (ORR) catalysis in polymer electrolyte membrane (PEM) fuel cells. To date, Co-based corroles have shown the best performance, but still suffer from the poor stability and the toxicity of the Co metal. Mn-N-C are more stable than Co-N-C, and are also less reactive towards peroxide formation. In this work, using first-principles density functional theory calculations, we study the improvement of the Mn-based corrole ORR activity by exploiting hydrogen bonding with substituent groups to modify the adsorption energies of the ORR intermediates and obtain higher onset potential (Vonset) values. We found that by using phenyl acetic acid as a substituent, Vonset increased from 0.54 V for the unsubstituted corrole to ∼0.9 V which is competitive with the Vonset of the Co-based corroles. Our results suggest that hydrogen bonding with substituent groups should be considered in the analysis and design of the reactivity of active sites in non-PGM ORR catalysts.

Resonant domain-wall-enhanced tunable microwave ferroelectrics.

Ordering of ferroelectric polarization1 and its trajectory in response to an electric field2 are essential for the operation of non-volatile memories3, transducers4 and electro-optic devices5. However, for voltage control of capacitance and frequency agility in telecommunication devices, domain walls have long been thought to be a hindrance because they lead to high dielectric loss and hysteresis in the device response to an applied electric field6. To avoid these effects, tunable dielectrics are often operated under piezoelectric resonance conditions, relying on operation well above the ferroelectric Curie temperature7, where tunability is compromised. Therefore, there is an unavoidable trade-off between the requirements of high tunability and low loss in tunable dielectric devices, which leads to severe limitations on their figure of merit. Here we show that domain structure can in fact be exploited to obtain ultralow loss and exceptional frequency selectivity without piezoelectric resonance. We use intrinsically tunable materials with properties that are defined not only by their chemical composition, but also by the proximity and accessibility of thermodynamically predicted strain-induced, ferroelectric domain-wall variants8. The resulting gigahertz microwave tunability and dielectric loss are better than those of the best film devices by one to two orders of magnitude and comparable to those of bulk single crystals. The measured quality factors exceed the theoretically predicted zero-field intrinsic limit owing to domain-wall fluctuations, rather than field-induced piezoelectric oscillations, which are usually associated with resonance. Resonant frequency tuning across the entire L, S and C microwave bands (1-8 gigahertz) is achieved in an individual device-a range about 100 times larger than that of the best intrinsically tunable material. These results point to a rich phase space of possible nanometre-scale domain structures that can be used to surmount current limitations, and demonstrate a promising strategy for obtaining ultrahigh frequency agility and low-loss microwave devices.

First-Principles Investigation of the Formation of Pt Nanorafts on a Mo2C Support and Their Catalytic Activity for Oxygen Reduction Reaction.

We use first-principles calculations to study the formation of Pt nanorafts and their oxygen reduction reaction (ORR) catalytic activity on Mo2C. Due to the high Pt binding energy on C atoms, Pt forms sheet-like structures on the Mo2C surface instead of agglomerating into particles. We find that the disordered Mo2C surface carbon arrangement limits the Pt sheet growth, leading to the formation of 4-6 atom Pt nanorafts. The O-O repulsion between the O atoms on the Mo2C and O adsorbate enhances the ORR activity by weakening the O adsorption energy. We find a significant change from the usual scaling of the energies of the intermediates in the ORR pathway and a strong interaction between the nanoraft and water that lead to a high activity of the Pt nanorafts. Fundamentally, our work demonstrates that the activity of metal catalysts can be strongly affected by manipulation of the atomic arrangement of the supporting carbide surface.

First-Principles Investigation of Black Phosphorus Synthesis.

Black phosphorus (BP) is a layered semiconductor with outstanding properties, making it a promising candidate for optoelectronic and other applications. BP synthesis is an intriguing task largely due to the insufficient understanding of the synthesis mechanism. In this work, we use density functional theory calculations to examine BP and its precursor red phosphorus as they are formed from P4 building blocks. Our results suggest that, without external effects such as pressure or addition of a catalyst, the precursor is energetically favored in the initial steps of the synthesis, even though BP is the more stable allotrope. The higher energy of BP is dictated by its 2D geometry that gives rise to the higher number of high-energy strained bonds at the edge compared to the 1D geometry of red phosphorus. The elucidated BP formation pathway provides a natural explanation for the effectiveness of the recently discovered Sn/I catalyst used in BP synthesis.

Slush-like polar structures in single-crystal relaxors.

Despite more than 50 years of investigation, it is still unclear how the underlying structure of relaxor ferroelectrics gives rise to their defining properties, such as ultrahigh piezoelectric coefficients, high permittivity over a broad temperature range, diffuse phase transitions, strong frequency dependence in dielectric response, and phonon anomalies. The model of polar nanoregions inside a non-polar matrix has been widely used to describe the structure of relaxor ferroelectrics. However, the lack of precise knowledge about the shapes, growth and dipole patterns of polar nanoregions has led to the characterization of relaxors as "hopeless messes", and no predictive model for relaxor behaviour is currently available. Here we use molecular dynamics simulations of the prototypical Pb(Mg1/3,Nb2/3)O3-PbTiO3 relaxor material to examine its structure and the spatial and temporal polarization correlations. Our simulations show that the unusual properties of relaxors stem from the presence of a multi-domain state with extremely small domain sizes (2-10 nanometres), and no non-polar matrix, owing to the local dynamics. We find that polar structures in the multi-domain state in relaxors are analogous to those of the slush state of water. The multi-domain structure of relaxors that is revealed by our molecular dynamics simulations is consistent with recent experimental diffuse scattering results and indicates that relaxors have a high density of low-angle domain walls. This insight explains the recently discovered classes of relaxors that cannot be described by the polar nanoregion model, and provides guidance for the design and synthesis of new relaxor materials.

Intrinsic ferroelectric switching from first principles.

The existence of domain walls, which separate regions of different polarization, can influence the dielectric, piezoelectric, pyroelectric and electronic properties of ferroelectric materials. In particular, domain-wall motion is crucial for polarization switching, which is characterized by the hysteresis loop that is a signature feature of ferroelectric materials. Experimentally, the observed dynamics of polarization switching and domain-wall motion are usually explained as the behaviour of an elastic interface pinned by a random potential that is generated by defects, which appear to be strongly sample-dependent and affected by various elastic, microstructural and other extrinsic effects. Theoretically, connecting the zero-kelvin, first-principles-based, microscopic quantities of a sample with finite-temperature, macroscopic properties such as the coercive field is critical for material design and device performance; and the lack of such a connection has prevented the use of techniques based on ab initio calculations for high-throughput computational materials discovery. Here we use molecular dynamics simulations of 90° domain walls (separating domains with orthogonal polarization directions) in the ferroelectric material PbTiO3 to provide microscopic insights that enable the construction of a simple, universal, nucleation-and-growth-based analytical model that quantifies the dynamics of many types of domain walls in various ferroelectrics. We then predict the temperature and frequency dependence of hysteresis loops and coercive fields at finite temperatures from first principles. We find that, even in the absence of defects, the intrinsic temperature and field dependence of the domain-wall velocity can be described with a nonlinear creep-like region and a depinning-like region. Our model enables quantitative estimation of coercive fields, which agree well with experimental results for ceramics and thin films. This agreement between model and experiment suggests that, despite the complexity of ferroelectric materials, typical ferroelectric switching is largely governed by a simple, universal mechanism of intrinsic domain-wall motion, providing an efficient framework for predicting and optimizing the properties of ferroelectric materials.