Surface phase diagrams from nested sampling.

Studies in atomic-scale modeling of surface phase equilibria often focus on temperatures near zero Kelvin due to the challenges in calculating the free energy of surfaces at finite temperatures. The Bayesian-inference-based nested sampling (NS) algorithm allows for modeling phase equilibria at arbitrary temperatures by directly and efficiently calculating the partition function, whose relationship with free energy is well known. This work extends NS to calculate adsorbate phase diagrams, incorporating all relevant configurational contributions to the free energy. We apply NS to the adsorption of Lennard-Jones (LJ) gas particles on low-index and vicinal LJ solid surfaces and construct the canonical partition function from these recorded energies to calculate ensemble averages of thermodynamic properties, such as the constant-volume heat capacity and order parameters that characterize the structure of adsorbate phases. Key results include determining the nature of phase transitions of adsorbed LJ particles on flat and stepped LJ surfaces, which typically feature an enthalpy-driven condensation at higher temperatures and an entropy-driven reordering process at lower temperatures, and the effect of surface geometry on the presence of triple points in the phase diagrams. Overall, we demonstrate the ability and potential of NS for surface modeling.

Factors Governing Oxygen Vacancy Formation in Oxide Perovskites.

The control of oxygen vacancy (VO) formation is critical to advancing multiple metal-oxide-perovskite-based technologies. We report the construction of a compact linear model for the neutral VO formation energy in ABO3 perovskites that reproduces, with reasonable fidelity, Hubbard-U-corrected density functional theory calculations based on the state-of-the-art, strongly constrained and appropriately normed exchange-correlation functional. We obtain a mean absolute error of 0.45 eV for perovskites stable at 298 K, an accuracy that holds across a large, electronically diverse set of ABO3 perovskites. Our model considers perovskites containing alkaline-earth metals (Ca, Sr, and Ba) and lanthanides (La and Ce) on the A-site and 3d transition metals (Ti, V, Cr, Mn, Fe, Co, and Ni) on the B-site in six different crystal systems (cubic, tetragonal, orthorhombic, hexagonal, rhombohedral, and monoclinic) common to perovskites. Physically intuitive metrics easily extracted from existing experimental thermochemical data or via inexpensive quantum mechanical calculations, including crystal bond dissociation energies and (solid phase) reduction potentials, are key components of the model. Beyond validation of the model against known experimental trends in materials used in solid oxide fuel cells, the model yields new candidate perovskites not contained in our training data set, such as (Bi,Y)(Fe,Co)O3, which we predict may have favorable thermochemical water-splitting properties. The confluence of sufficient accuracy, efficiency, and interpretability afforded by our model not only facilitates high-throughput computational screening for any application that requires the precise control of VO concentrations but also provides a clear picture of the dominant physics governing VO formation in metal-oxide perovskites.

Chemical Pressure-Driven Enhancement of the Hydrogen Evolving Activity of Ni2P from Nonmetal Surface Doping Interpreted via Machine Learning.

The activity of Ni2P catalysts for the hydrogen evolution reaction (HER) is currently limited by strong H adsorption at the Ni3-hollow site. We investigate the effect of surface nonmetal doping on the HER activity of the Ni3P2 termination of Ni2P(0001), which is stable at modest electrochemical conditions. Using density functional theory (DFT) calculations, we find that both 2 p nonmetals and heavier chalcogens provide nearly thermoneutral H adsorption at moderate surface doping concentrations. We also find, however, that only chalcogen substitution for surface P is exergonic. For intermediate surface concentrations of S, the free energy of H adsorption at the Ni3-hollow site is -0.11 eV, which is significantly more thermoneutral than the undoped surface (-0.45 eV). We use the regularized random forest machine learning algorithm to discover the relative importance of structure and charge descriptors, extracted from the DFT calculations, in determining the HER activity of Ni2P(0001) under different doping concentrations. We discover that the Ni-Ni bond length is the most important descriptor of HER activity, which suggests that the nonmetal dopants induce a chemical pressure-like effect on the Ni3-hollow site, changing its reactivity through compression and expansion.

Role of proton hopping in surface charge transport on tin dioxide as revealed by the thermal dependence of conductance.

The presence of water on an oxide surface can dramatically alter its electrical properties with important consequences for electrical measurements by scanning probe microscopy, and for the use of semiconducting oxides in sensing applications. Here, the thermal dependence of the conductance of tin dioxide is interpreted by combining semiconductor equilibrium carrier statistics with a proton hopping mechanism. First, the functional form of this charge transport model is fit to experimental conductance data for tin dioxide. Next, the important energy parameters in the model are computed with density functional theory. Comparing the values of the energy parameters obtained by fitting, to the values for the same parameters obtained from electronic structure calculations, yields new insight into the surface charge transport in tin dioxide. In particular, it is found that mobile protons, freed from the dissociative adsorption of water on the [110] surface, are an essential component of the observed thermal dependence of conductance in tin dioxide.

Engineering Relaxor Behavior in (BaTiO3 )n /(SrTiO3 )n Superlattices.

Complex-oxide superlattices provide a pathway to numerous emergent phenomena because of the juxtaposition of disparate properties and the strong interfacial interactions in these unit-cell-precise structures. This is particularly true in superlattices of ferroelectric and dielectric materials, wherein new forms of ferroelectricity, exotic dipolar textures, and distinctive domain structures can be produced. Here, relaxor-like behavior, typically associated with the chemical inhomogeneity and complexity of solid solutions, is observed in (BaTiO3 )n /(SrTiO3 )n (n = 4-20 unit cells) superlattices. Dielectric studies and subsequent Vogel-Fulcher analysis show significant frequency dispersion of the dielectric maximum across a range of periodicities, with enhanced dielectric constant and more robust relaxor behavior for smaller period n. Bond-valence molecular-dynamics simulations predict the relaxor-like behavior observed experimentally, and interpretations of the polar patterns via 2D discrete-wavelet transforms in shorter-period superlattices suggest that the relaxor behavior arises from shape variations of the dipolar configurations, in contrast to frozen antipolar stripe domains in longer-period superlattices (n = 16). Moreover, the size and shape of the dipolar configurations are tuned by superlattice periodicity, thus providing a definitive design strategy to use superlattice layering to create relaxor-like behavior which may expand the ability to control desired properties in these complex systems.

Large-area synthesis of high-quality monolayer 1T'-WTe2 flakes.

Large-area growth of monolayer films of the transition metal dichalcogenides is of the utmost importance in this rapidly advancing research area. The mechanical exfoliation method offers high quality monolayer material but it is a problematic approach when applied to materials that are not air stable. One important example is 1T'-WTe2, which in multilayer form is reported to possess a large non saturating magnetoresistance, pressure induced superconductivity, and a weak antilocalization effect, but electrical data for the monolayer is yet to be reported due to its rapid degradation in air. Here we report a reliable and reproducible large-area growth process for obtaining many monolayer 1T'-WTe2 flakes. We confirmed the composition and structure of monolayer 1T'-WTe2 flakes using x-ray photoelectron spectroscopy, energy-dispersive x-ray spectroscopy, atomic force microscopy, Raman spectroscopy and aberration corrected transmission electron microscopy. We studied the time dependent degradation of monolayer 1T'-WTe2 under ambient conditions, and we used first-principles calculations to identify reaction with oxygen as the degradation mechanism. Finally we investigated the electrical properties of monolayer 1T'-WTe2 and found metallic conduction at low temperature along with a weak antilocalization effect that is evidence for strong spin-orbit coupling.