# Publications

2016
Personick ML, Montemore MM, Kaxiras E, Madix RJ, Biener J, Friend CM. Catalyst design for enhanced sustainability through fundamental surface chemistry. PHILOSOPHICAL TRANSACTIONS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES. 2016;374 (2061).Abstract
Decreasing energy consumption in the production of platform chemicals is necessary to improve the sustainability of the chemical industry, which is the largest consumer of delivered energy. The majority of industrial chemical transformations rely on catalysts, and therefore designing new materials that catalyse the production of important chemicals via more selective and energy-efficient processes is a promising pathway to reducing energy use by the chemical industry. Efficiently designing new catalysts benefits from an integrated approach involving fundamental experimental studies and theoretical modelling in addition to evaluation of materials under working catalytic conditions. In this review, we outline this approach in the context of a particular catalyst-nanoporous gold (npAu)-which is an unsupported, dilute AgAu alloy catalyst that is highly active for the selective oxidative transformation of alcohols. Fundamental surface science studies on Au single crystals and AgAu thin-film alloys in combination with theoretical modelling were used to identify the principles which define the reactivity of npAu and subsequently enabled prediction of new reactive pathways on this material. Specifically, weak van der Waals interactions are key to the selectivity of Au materials, including npAu. We also briefly describe other systems in which this integrated approach was applied.
Montemore MM, Cubuk ED, Klobas JE, Schmid M, Madix RJ, Friend CM, Kaxiras E. Controlling O coverage and stability by alloying Au and Ag. PHYSICAL CHEMISTRY CHEMICAL PHYSICS. 2016;18 (38) :26844-26853.Abstract
One of the most critical factors in oxidation catalysis is controlling the state of oxygen on the surface. Au and Ag are both effective selective oxidation catalysts for various reactions, and their interactions with oxygen are critical for determining their catalytic performance. Here, we show that the state of oxygen on a catalytic surface can be controlled by alloying Au and Ag. Using temperature programmed desorption, density functional theory (DFT), and Monte Carlo simulations, we examine how alloying Au into an Ag(110) surface affects O-2 dissociation, O coverage, and O stability. DFT calculations indicate that Au resides in the second layer, in agreement with previous experimental findings. The minimum ensemble size for O-2 dissociation is 2 Ag atoms in adjacent rows of the second layer. Surprisingly, adsorbed O-2 and the dissociation transition state interact directly with metal atoms in the adjacent trough, such that Au in this position inhibits O-2 dissociation by direct repulsion with oxygen electronic states. Using Monte Carlo simulations based on DFT energetics, we create models of the surface that agree closely with our experimental results. Both show that the O-2 uptake decreases nearly linearly as the Au concentration increases, and no O-2 uptake occurs for Au concentrations above 50%. For Au concentrations greater than 18%, increasing the Au concentration also decreases the stability of the adsorbed O. Based on these results, the O coverage and O stability can be tuned, in some cases independently. We also study how the reactivity of the surface is affected by these factors using CO2 oxidation as a simple test reaction.
Zhou Y, Chen W, Cui P, Zeng J, Lin Z, Kaxiras E, Zhang Z. Enhancing the Hydrogen Activation Reactivity of Nonprecious Metal Substrates via Confined Catalysis Underneath Graphene. NANO LETTERS. 2016;16 (10) :6058-6063.Abstract
In the hydrogen evolution reaction (HER), the reactivity as a function of the hydrogen adsorption energy on different metal substrates follows a well-known;volcano curve, peaked a the precious inetal Pt. The goal. of turning nonprecious metals into efficient catalysts for HER and other important chemical reactions is a :fundamental challenge; it is-also of technological significance. Here, we present results toward, achieving Ihis goal by exploiting the synergistic power of marginal catalysis and confined catalysis. Using density functional theory calculations, we first show that the volcano curve stays qualitatively intact when van der Waals attractions between a hydrogen adatom and different metal (111) surfaces are included. We further show that the hydrogen adsorption:energy on the metal: surfaces is weakened by 0.11-0.23 eV when hydrogen is confinedbetween graphene and the metal surfaces, with Ni exhibiting the largest change. Inparticular, we find that the graphene-modified volcano curve peaks around Ni, whose bare surface already possesses moderate (or :marginal) reactivity,i and the corresponding HER rate of grapllene-covered comparable to that of bare Pt. A hydrogen adatom has high mobility within the confined geometry. These findings demonstrate that graphene-coyeted Ni is an appealing effective stable, and economical catalytic platform for HER.
Islam MM, Kolesov G, Verstraelen T, Kaxiras E, van Duin ACT. eReaxFF: A Pseudoclassical Treatment of Explicit Electrons within Reactive Force Field Simulations. JOURNAL OF CHEMICAL THEORY AND COMPUTATION. 2016;12 (8) :3463-3472.Abstract
We present a computational tool, eReaxFF, for simulating explicit electrons within the framework of the standard ReaxFF reactive force field method. We treat electrons explicitly in a pseudoclassical manner that enables simulation several orders of magnitude faster than quantum chemistry (QC) methods, while retaining the ReaxFF transferability. We delineate here the fundamental concepts of the eReaxFF method and the integration of the Atom condensed Kohn-Sham DFT approximated to second order (ACKS2) charge calculation scheme into the eReaxFF. We trained our force field to capture electron affinities (EA) of various species. As a proof-of-principle, we performed a set of molecular dynamics (MD) simulations with an explicit electron model for representative hydrocarbon radicals. We establish a good qualitative agreement of EAs of various species with experimental data, and MD simulations with eReaxFF agree well with the corresponding Ehrenfest dynamics simulations. The standard ReaxFF parameters available in the literature are transferrable to the eReaxFF method. The computationally economic eReaxFF method will be a useful tool for studying large-scale chemical and physical systems with explicit electrons as an alternative to computationally demanding QC methods.
Montemore MM, Madix RJ, Kaxiras E. How Does Nanoporous Gold Dissociate Molecular Oxygen?. JOURNAL OF PHYSICAL CHEMISTRY C. 2016;120 (30) :16636-16640.Abstract
Nanoporous Au and other dilute AgAu alloys are highly active and selective oxidation catalysts. Their ability to dissociate O-2 is to a large extent unexplained, given that unsupported Au cannot generally dissociate O-2 while large ensembles of Ag atoms (>4) are generally necessary to lower the O-2 dissociation barrier significantly. Here, we identify a site on the surface of dilute AgAu alloys that is stable under reaction conditions and has a low O-2 dissociation barrier, in agreement with experimental measurements. Although Ag generally prefers to disperse throughout Au, the presence of adsorbed O near surface steps creates sites of high local Ag concentration, where the Ag, atoms sit in the rows next to the step Au atoms. O-2 adsorbs on the Au step atoms, but the transition state involves significant Ag-O interaction, resulting in a barrier lower than expected from the adsorption energies of either the initial or final state.
Cao Y, Luo JY, Fatemi V, Fang S, Sanchez-Yamagishi JD, Watanabe K, Taniguchi T, Kaxiras E, Jarillo-Herrero P. Superlattice-Induced Insulating States and Valley-Protected Orbits in Twisted Bilayer Graphene. PHYSICAL REVIEW LETTERS. 2016;117 (11).Abstract
Twisted bilayer graphene (TBLG) is one of the simplest van der Waals heterostructures, yet it yields a complex electronic system with intricate interplay between moire physics and interlayer hybridization effects. We report on electronic transport measurements of high mobility small angle TBLG devices showing clear evidence for insulating states at the superlattice band edges, with thermal activation gaps several times larger than theoretically predicted. Moreover, Shubnikov-de Haas oscillations and tight binding calculations reveal that the band structure consists of two intersecting Fermi contours whose crossing points are effectively unhybridized. We attribute this to exponentially suppressed interlayer hopping amplitudes for momentum transfers larger than the moire wave vector.
Huang D, Webb TA, Fang S, Song C-L, Chang C-Z, Moodera JS, Kaxiras E, Hoffman JE. Bounds on nanoscale nematicity in single-layer FeSe/SrTiO3. PHYSICAL REVIEW B. 2016;93 (12).Abstract
We use scanning tunneling microscopy (STM) and quasiparticle interference (QPI) imaging to investigate the low-energy orbital texture of single-layer FeSe/SrTiO3. We develop a T-matrix model of multiorbital QPI to disentangle scattering intensities from Fe 3d(xz) and 3d(yz) bands, enabling the use of STM as a nanoscale detection tool of nematicity. By sampling multiple spatial regions of a single-layer FeSe/SrTiO3 film, we quantitatively exclude static xz/yz orbital ordering with domain size larger than delta r(2) = 20 nm x 20 nm, xz/yz Fermi wave vector difference larger than delta k = 0.014 pi, and energy splitting larger than delta E = 3.5 meV. The lack of detectable ordering pinned around defects places qualitative constraints on models of fluctuating nematicity.
Chen W, Zeng C, Kaxiras E, Zhang Z. Dual role of Fe dopants in enhancing stability and charge transfer in (Li0.8Fe0.2)OHFeSe superconductors. PHYSICAL REVIEW B. 2016;93 (6).Abstract
The recently discovered (Li0.8Fe0.2)OHFeSe superconductor provides a new platform for exploiting the microscopic mechanisms of high-T-c superconductivity in FeSe-derived systems. Using density functional theory calculations, we first show that substitution of Li by Fe not only significantly strengthens the attraction between the (Li0.8Fe0.2)OH spacing layers and the FeSe superconducting layers along the c axis but also minimizes the lattice mismatch between the two in the ab plane, both favorable for stabilizing the overall structure. Next, we explore the electron injection into FeSe from the spacing layers and unambiguously identify the Fe-0.2 components to be the origin of the dramatically enhanced interlayer charge transfer. We further reveal that the system strongly favors collinear antiferromagnetic ordering in the FeSe layers, but the spacing layers can be either antiferromagnetic or ferromagnetic depending on the Fe-0.2 spatial distribution. Based on these insights, we predict (Li0.8Co0.2)OHFeSe to be structurally stable with even larger electron injection and potentially higher T-c.
Kaplan D, Gong Y, Mills K, Swaminathan V, Ajayan PM, Shirodkar S, Kaxiras E. Excitation intensity dependence of photoluminescence from monolayers of MoS2 and WS2/MoS2 heterostructures. 2D MATERIALS. 2016;3 (1).Abstract
A detailed study of the excitation dependence of the photoluminescence (PL) from monolayers of MoS2 and WS2/MoS2 heterostructures grown by chemical vapor deposition on Si substrates has revealed that the luminescence from band edge excitons from MoS2 monolayers shows a linear dependence on excitation intensity for both above band gap and resonant excitation conditions. In particular, a band separated by similar to 55 meV from the A exciton, referred to as the C band, shows the same linear dependence on excitation intensity as the band edge excitons. A band similar to the C band has been previously ascribed to a trion, a charged, three-particle exciton. However, in our study the C band does not show the 3/2 power dependence on excitation intensity as would be expected for a three-particle exciton. Further, the PL from the MoS2 monolayer in a bilayer WS2/MoS2 heterostructure, under resonant excitation conditions where only the MoS2 absorbs the laser energy, also revealed a linear dependence on excitation intensity for the C band, confirming that its origin is not due to a trion but instead a bound exciton, presumably of an unintentional impurity or a native point defect such as a sulfur vacancy. The PL from the WS2/MoS2 heterostructure, under resonant excitation conditions also showed additional features which are suggested to arise from the interface states at the heteroboundary. Further studies are required to clearly identify the origin of these features.
Perez D, Cubuk ED, Waterland A, Kaxiras E, Voter AF. Long-Time Dynamics through Parallel Trajectory Splicing. JOURNAL OF CHEMICAL THEORY AND COMPUTATION. 2016;12 (1) :18-28.Abstract
Simulating the atomistic evolution of materials over long time scales is a longstanding challenge, especially for complex systems where the distribution of barrier heights is very heterogeneous. Such systems are difficult to investigate using conventional long-time scale techniques, and the fact that they tend to remain trapped in small regions of configuration space for extended periods of time strongly limits the physical insights gained from short simulations. We introduce a novel simulation technique, Parallel Trajectory Splicing (Par Splice), that aims at addressing this problem through the timewise parallelization of long trajectories. The computational efficiency of Par Splice stems from a speculation strategy whereby predictions of the future evolution of the system are leveraged to increase the amount of work that can be concurrently performed at any one time, hence improving the scalability of the method. ParSplice is also able to accurately account for, and potentially reuse, a substantial fraction of the computational work invested in the simulation. We validate the method on a simple Ag surface system and demonstrate substantial increases in efficiency compared to previous methods. We then demonstrate the power of ParSplice through the study of topology changes in Ag42Cu13 core shell nanoparticles.
Heller EJ, Yang Y, Kocia L, Chen W, Fang S, Borunda M, Kaxiras E. Theory of Graphene Raman Scattering. ACS NANO. 2016;10 (2) :2803-2818.Abstract
Raman scattering plays a key role in unraveling the quantum dynamics of graphene, perhaps the most promising material of recent times. It is crucial to correctly interpret the meaning of the spectra. It is therefore very surprising that the widely accepted understanding of Raman scattering, i.e., Kramers Heisenberg Dirac theory, has never been applied to graphene. Doing so here, a remarkable mechanism we term''transition sliding'' is uncovered, explaining the uncommon brightness of overtones in graphene. Graphene's dispersive and fixed Raman bands, missing bands, defect density and laser frequency dependence of band intensities, widths of overtone bands, Stokes, anti -Stokes anomalies, and other known properties emerge simply and directly.
Hiebel F, Montemore MM, Kaxiras E, Friend CM. Direct visualization of quasi-ordered oxygen chain structures on Au(110)-(1 x 2). SURFACE SCIENCE. 2016;650 (SI) :5-10.Abstract
The Au(110) surface offers unique advantages for atomically-resolved model studies of catalytic oxidation processes on gold. We investigate the adsorption of oxygen on Au(110) using a combination of scanning tunneling microscopy (STM) and density functional theory (DFT) methods. We identify the typical (empty-states) STM contrast resulting from adsorbed oxygen as atomic-sized dark features of electronic origin. DFT-based image simulations confirm that chemisorbed oxygen is generally detected indirectly, from the binding -induced electronic structure modification of gold. STM images show that adsorption occurs without affecting the general structure of the pristine Au(110) missing-row reconstruction. The tendency to form one-dimensional structures is observed already at low coverage (<0.05 ML), with oxygen adsorbing on alternate sides of the reconstruction ridges. Consistently, calculations yield preferred adsorption on the (111) facets of the reconstruction, on a 3-fold coordination site, with increased stability when adsorbed in chains. Gold atoms with two oxygen neighbors exhibit enhanced electronic hybridization with the O states. Finally, the species observed are reactive to CO oxidation at 200 K and desorption of CO2 leaves a clean and ordered gold surface. (C) 2015 Elsevier B.V. All rights reserved,
Huang D, Webb TA, Song C-L, Chang C-Z, Moodera JS, Kaxiras E, Hoffinan JE. Dumbbell Defects in FeSe Films: A Scanning Tunneling Microscopy and First-Principles Investigation. NANO LETTERS. 2016;16 (7) :4224-4229.Abstract
The properties of iron-based superconductors (Fe-SCs) can be varied dramatically with the introduction of dopants and atomic defects. As a pressing example, FeSe, parent phase of the highest-T-c Fe-SC, exhibits prevalent defects with atomic-scale dumbbell'' signatures as imaged by scanning tunneling microscopy (STM). These defects spoil superconductivity when their concentration exceeds 2.5%. Resolving their chemical identity is a prerequisite to applications such as nanoscale patterning of superconducting/nonsuperconducting regions in FeSe as well as fundamental questions such as the mechanism of superconductivity and the path by which the defects destroy it. We use STM and density functional theory to characterize and identify the dumbbell defects. In contrast to previous speculations about Se adsorbates or substitutions, we find that an Fe-site vacancy is the most energetically favorable defect in Se-rich conditions and reproduces our observed STM signature. Our calculations shed light more generally on the nature of Se capping, the removal of Fe vacancies via annealing, and their ordering into a root 5 x root 5 superstructure in FeSe and related alkali-doped compounds.
Fang S, Kaxiras E. Electronic structure theory of weakly interacting bilayers. PHYSICAL REVIEW B. 2016;93 (23).Abstract
We derive electronic structure models for weakly interacting bilayers such as graphene-graphene and graphene-hexagonal boron nitride, based on density functional theory calculations followed by Wannier transformation of electronic states. These transferable interlayer coupling models can be applied to investigate the physics of bilayers with arbitrary translations and twists. The functional form, in addition to the dependence on the distance, includes the angular dependence that results from higher angular momentum components in the Wannier p(z) orbitals. We demonstrate the capabilities of the method by applying it to a rotated graphene bilayer, which produces the analytically predicted renormalization of the Fermi velocity, Van Hove singularities in the density of states, and moire pattern of the electronic localization at small twist angles. We further extend the theory to obtain the effective couplings by integrating out neighboring layers. This approach is instrumental for the design of van der Walls heterostructures with desirable electronic features and transport properties and for the derivation of low-energy theories for graphene stacks, including proximity effects from other layers.
Qi S, Qiao Z, Deng X, Cubuk ED, Chen H, Zhu W, Kaxiras E, Zhang SB, Xu X, Zhang Z. High-Temperature Quantum Anomalous Hall Effect in n-p Codoped Topological Insulators. PHYSICAL REVIEW LETTERS. 2016;117 (5).Abstract
The quantum anomalous Hall effect (QAHE) is a fundamental quantum transport phenomenon that manifests as a quantized transverse conductance in response to a longitudinally applied electric field in the absence of an external magnetic field, and it promises to have immense application potential in future dissipationless quantum electronics. Here, we present a novel kinetic pathway to realize the QAHE at high temperatures by n-p codoping of three-dimensional topological insulators. We provide a proof-of-principle numerical demonstration of this approach using vanadium-iodine (V-I) codoped Sb2Te3 and demonstrate that, strikingly, even at low concentrations of similar to 2% V and similar to 1% I, the system exhibits a quantized Hall conductance, the telltale hallmark of QAHE, at temperatures of at least similar to 50 K, which is 3 orders of magnitude higher than the typical temperatures at which it has been realized to date. The underlying physical factor enabling this dramatic improvement is tied to the largely preserved intrinsic band gap of the host system upon compensated n-p codoping. The proposed approach is conceptually general and may shed new light in experimental realization of high-temperature QAHE.
Shirodkar SN, Kaxiras E. Li intercalation at graphene/hexagonal boron nitride interfaces. PHYSICAL REVIEW B. 2016;93 (24).Abstract
Intercalation of Li in graphite and other layered structures is of interest for highly efficient energy storage devices. In this paper, we determine the extent to which Li intercalates at the different interfaces formed between graphene (G) and hexagonal boron nitride (hBN) heterostructures. We use ab initio calculations to explore in detail the position of the dispersed Li atoms, changes in the structure at the interfaces, energetic stability of the configurations, and the corresponding electronic structure with varying concentrations of the intercalant. We trace the origin of the energetic stability and maximum concentration of Li that intercalates into various layered structures to the ability of the interface to accept electrons. Our calculations indicate that Li intercalates easiest at G/G interfaces, followed by interfaces between G/hBN, whereas Li cannot intercalate in hBN/hBN interfaces. Our results provide a framework for the design of experimental setups with optimal Li intercalation and reveal the implications of intercalation on the dielectric properties of these materials and their possible application in plasmonics.
Falcucci G, Succi S, Montessori A, Melchionna S, Prestininzi P, Barroo C, Bell DC, Biener MM, Biener J, Zugic B, et al. Mapping reactive flow patterns in monolithic nanoporous catalysts. MICROFLUIDICS AND NANOFLUIDICS. 2016;20 (7).Abstract
The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nanoscale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domains. To address this problem, we developed a multiscale mass transport computational framework based on the lattice Boltzmann method that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.
Tritsaris GA, Shirodkar SN, Kaxiras E, Cazeaux P, Luskin M, Plechac P, Cances E. Perturbation theory for weakly coupled two-dimensional layers. JOURNAL OF MATERIALS RESEARCH. 2016;31 (7) :959-966.Abstract
A key issue in two-dimensional structures composed of atom-thick sheets of electronic materials is the dependence of the properties of the combined system on the features of its parts. Here, we introduce a simple framework for the study of the electronic structure of layered assemblies based on perturbation theory. Within this framework, we calculate the band structure of commensurate and twisted bilayers of graphene (Gr) and hexagonal boron nitride (h-BN), and of a Gr/h-BN heterostructure, which we compare with reference full-scale density functional theory calculations. This study presents a general methodology for computationally efficient calculations of two-dimensional materials and also demonstrates that for relatively large twist in the graphene bilayer, the perturbation of electronic states near the Fermi level is negligible.
S. S. Schoenholz, Cubuk ED, Sussman DM, Kaxiras E, Liu AJ. A structural approach to relaxation in glassy liquids. NATURE PHYSICS. 2016;12 (5) :469+.Abstract
In contrast with crystallization, there is no noticeable structural change at the glass transition. Characteristic features of glassy dynamics that appear below an onset temperature, T-0 (refs 1-3), are qualitatively captured by mean field theory(4-6), which assumes uniform local structure. Studies of more realistic systems have found only weak correlations between structure and dynamics(7-11). This raises the question: is structure important to glassy dynamics in three dimensions? We answer this question affirmatively, using machine learning to identify a new field, softness' which characterizes local structure and is strongly correlated with dynamics. We find that the onset of glassy dynamics at T-0 corresponds to the onset of correlations between softness (that is, structure) and dynamics. Moreover, we construct a simple model of relaxation that agrees well with our simulation results, showing that a theory of the evolution of softness in time would constitute a theory of glassy dynamics.
Cubuk ED, Schoenholz SS, Kaxiras E, Liu AJ. Structural Properties of Defects in Glassy Liquids. JOURNAL OF PHYSICAL CHEMISTRY B. 2016;120 (26) :6139-6146.Abstract
At zero temperature a disordered solid corresponds to a local minimum in the energy landscape. As the temperature is raised or the system is driven with a mechanical load, the system explores different minima via dynamical events in which particles rearrange their relative positions. We have shown recently that the dynamics of particle rearrangements are strongly correlated with a structural quantity associated with each particle, `softness'', which we can identify using supervised machine learning. Particles of a given softness have a well-defined energy scale that governs local rearrangements; because of this property, softness greatly simplifies our understanding of glassy dynamics. Here we investigate the correlation of softness with other commonly used structural quantities, such as coordination number and local potential energy. We show that although softness strongly correlates with these properties, its predictive power for rearrangement dynamics is much higher. We introduce a useful metric for quantifying the quality of structural quantities as predictors of dynamics. We hope that, in the future, authors introducing new structural measures of dynamics will compare their proposals quantitatively to softness using this metric. We also show how softness correlations give insight into rearrangements. Finally, we explore the physical meaning of softness using unsupervised dimensionality reduction and reduced curve-fitting, models, and show that softness can be recast in a form that is amenable to analytical treatment.