@article { ISI:000446141900028,
	title = {Anhydrous Methanol and Ethanol Dehydrogenation at Cu(111) Step Edges},
	journal = {JOURNAL OF PHYSICAL CHEMISTRY C},
	volume = {122},
	number = {38},
	year = {2018},
	month = {SEP 27},
	pages = {21952-21962},
	abstract = {Oxidative methanol dehydrogenation is a major industrial reaction with global formaldehyde production exceeding 30 million tonnes per year. Unfortunately, oxidative dehydrogenation produces water aldehyde mixtures that require subsequent distillation. Anhydrous alcohol dehydrogenation is a promising alternative that produces H-2 instead of water. Pursuant to recent experimental work showing that highly stepped Cu(111) surfaces exhibit anhydrous dehydrogenation activity, we present first-principles density functional theory calculations for methanol and ethanol dehydrogenation at Cu(111) step edges to provide an atomistic understanding of the catalytic mechanism; these sites stabilize all intermediates while reducing activation energies. We find that van der Waals contributions to the energy account for more than 50\% of adsorption energies, and essential in achieving good agreement with experimental desorption temperatures. Furthermore, vibrational zero-point energy corrections significantly reduce the activation energy for all reaction steps considered here. Hydrogen bonding among ethanol intermediates at step edges is weakened by geometric frustration. These insights lead us to propose several suggestions for further research on undercoordinated Cu sites as anhydrous alcohol dehydrogenation catalysts.},
	issn = {1932-7447},
	doi = {10.1021/acs.jpcc.8b06730},
	author = {Hoyt, Robert A. and Montemore, Matthew M. and Sykes, E. Charles H. and Efthimios Kaxiras}
}
@article { ISI:000445597200001,
	title = {Berry curvature dipole current in the transition metal dichalcogenides family},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {12},
	year = {2018},
	month = {SEP 25},
	abstract = {We study the quantum nonlinear Hall effect in two-dimensional (2D) materials with time-reversal symmetry. When only one mirror line exists, a transverse charge current occurs in the second-order response to an external electric field, as a result of the Berry curvature dipole in momentum space. Candidate 2D materials to observe this effect are two-dimensional transition metal dichalcogenides (TMDCs). First, we use an ab initio based tight-binding approach to demonstrate that monolayer T-d-structure TMDCs exhibit a finite Berry curvature dipole. In the 1H and 1T{\textquoteright} phase of TMDCs, we show the emergence of a finite Berry curvature dipole with the application of strain and an electrical displacement field, respectively.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.121109},
	author = {You, Jhih-Shih and Fang, Shiang and Xu, Su-Yang and Efthimios Kaxiras and Low, Tony}
}
@article { ISI:000430255900029,
	title = {A Comparative Ab Initio Study of Anhydrous Dehydrogenation of Linear-Chain Alcohols on Cu(110)},
	journal = {JOURNAL OF PHYSICAL CHEMISTRY C},
	volume = {122},
	number = {14},
	year = {2018},
	month = {APR 12},
	pages = {7806-7815},
	abstract = {The catalytic behavior of Cu surfaces in the anhydrous production of aldehydes from alcohols, a process of industrial significance, is puzzling: the two simplest alcohols (methanol and ethanol) show dramatically different decomposition behavior on Cu. Here, we study the thermodynamic and kinetic processes involved in the anhydrous dehydrogenation of linear-chain alcohols including methanol, ethanol, 1-propanol, and 1-butanol on the Cu(110) surface using multiscale approaches. First, we obtain the adsorption structures and energies of the reaction intermediates, in which van der Waals (vdW) interactions play a crucial role. Then, we determine the kinetic barriers for the two dehydrogenation steps, namely, the O-H and the subsequent C-H bond-breaking on Cu. The reaction of methoxy-to-formaldehyde has a rather high-energy transition state, in contrast to that of alkoxide-to-aldehyde in the longer-chain systems. This difference qualitatively explains the lower production efficiency of formaldehyde on Cu. Finally, we simulate the production rates of aldehydes based on which we optimize reaction conditions and propose possible avenues for enhancing the production of anhydrous formaldehyde using Cu-based catalysts.},
	issn = {1932-7447},
	doi = {10.1021/acs.jpcc.8b01698},
	author = {Chen, Wei and Ekin Dogus Cubuk and Montemore, Matthew M. and Reece, Christian and Madix, Robert J. and Cynthia M. Friend and Efthimios Kaxiras}
}
@article { ISI:000437380000004,
	title = {Compression of Wannier functions into Gaussian-type orbitals},
	journal = {COMPUTER PHYSICS COMMUNICATIONS},
	volume = {230},
	year = {2018},
	month = {SEP},
	pages = {27-37},
	abstract = {We propose a greedy algorithm for the compression of Wannier functions into Gaussian-polynomials orbitals. The so-obtained compressed Wannier functions can be stored in a very compact form, and can be used to efficiently parameterize effective tight-binding Hamiltonians for multilayer 2D materials for instance. The compression method preserves the symmetries (if any) of the original Wannier function. We provide algorithmic details, and illustrate the performance of our implementation on several examples, including graphene, hexagonal boron-nitride, single-layer FeSe, and bulk silicon in the diamond cubic structure. (C) 2018 Elsevier B.V. All rights reserved.},
	issn = {0010-4655},
	doi = {10.1016/j.cpc.2018.04.011},
	author = {Bakhta, Athmane and Cances, Eric and Cazeaux, Paul and Fang, Shiang and Efthimios Kaxiras}
}
@article { ISI:000429103300038,
	title = {Correlated insulator behaviour at half-filling in magic-angle graphene superlattices},
	journal = {NATURE},
	volume = {556},
	number = {7699},
	year = {2018},
	month = {APR 5},
	pages = {80+},
	abstract = {A van der Waals heterostructure is a type of metamaterial that consists of vertically stacked two-dimensional building blocks held together by the van der Waals forces between the layers. This design means that the properties of van der Waals heterostructures can be engineered precisely, even more so than those of two-dimensional materials(1). One such property is the {\textquoteleft}twist{\textquoteright} angle between different layers in the heterostructure. This angle has a crucial role in the electronic properties of van der Waals heterostructures, but does not have a direct analogue in other types of heterostructure, such as semiconductors grown using molecular beam epitaxy. For small twist angles, the moire pattern that is produced by the lattice misorientation between the two-dimensional layers creates long-range modulation of the stacking order. So far, studies of the effects of the twist angle in van der Waals heterostructures have concentrated mostly on heterostructures consisting of monolayer graphene on top of hexagonal boron nitride, which exhibit relatively weak interlayer interaction owing to the large bandgap in hexagonal boron nitride(2-5). Here we study a heterostructure consisting of bilayer graphene, in which the two graphene layers are twisted relative to each other by a certain angle. We show experimentally that, as predicted theoretically(6), when this angle is close to the {\textquoteleft}magic{\textquoteright} angle the electronic band structure near zero Fermi energy becomes flat, owing to strong interlayer coupling. These flat bands exhibit insulating states at half-filling, which are not expected in the absence of correlations between electrons. We show that these correlated states at half-filling are consistent with Mott-like insulator states, which can arise from electrons being localized in the superlattice that is induced by the moire pattern. These properties of magic-angle-twisted bilayer graphene heterostructures suggest that these materials could be used to study other exotic many-body quantum phases in two dimensions in the absence of a magnetic field. The accessibility of the flat bands through electrical tunability and the bandwidth tunability through the twist angle could pave the way towards more exotic correlated systems, such as unconventional superconductors and quantum spin liquids.},
	issn = {0028-0836},
	doi = {10.1038/nature26154},
	author = {Cao, Yuan and Fatemi, Valla and Demir, Ahmet and Fang, Shiang and Tomarken, Spencer L. and Luo, Jason Y. and Sanchez-Yamagishi, Javier D. and Kenji Watanabe and Takashi Taniguchi and Efthimios Kaxiras and Ashoori, Ray C. and Pablo Jarillo-Herrero}
}
@article { ISI:000449909600001,
	title = {Density functional theory beyond the Born-Oppenheimer approximation: Accurate treatment of the ionic zero-point motion},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {19},
	year = {2018},
	month = {NOV 12},
	abstract = {We introduce a method to carry out zero-temperature calculations within density functional theory (DFT) but without relying on the Born-Oppenheimer (BO) approximation for the ionic motion. Our approach is based on the finite-temperature many-body path-integral formulation of quantum mechanics by taking the zero-temperature limit and treating the imaginary-time propagation of the electronic variables in the context of DFT. This goes beyond the familiar BO approximation and is limited from being an exact treatment of both electrons and ions only by the approximations involved in the DFT component. We test our method in two simple molecules, H-2 and benzene. We demonstrate that the method produces a difference from the results of the BO approximation which is significant for many physical systems, especially those containing light atoms such as hydrogen; in these cases, we find that the fluctuations of the distance from its equilibrium position, due to the zero-point motion, is comparable to the interatomic distances. The method is suitable for use with conventional condensed-matter approaches and currently is implemented on top of the periodic pseudopotential code SIESTA.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.195112},
	author = {Grigory Kolesov and Efthimios Kaxiras and Manousakis, Efstratios}
}
@article { ISI:000452007300003,
	title = {Dihedral-angle-corrected registry-dependent interlayer potential for multilayer graphene structures},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {23},
	year = {2018},
	month = {DEC 4},
	abstract = {The structural relaxation of multilayer graphene is essential in describing the interesting electronic properties induced by intentional misalignment of successive layers, including the recently reported superconductivity in twisted bilayer graphene. This is difficult to accomplish without an accurate interatomic potential. Here, we present a new, registry-dependent Kolmogorov-Crespi-type interatomic potential to model interlayer interactions in multilayer graphene structures. It consists of two parts, representing attractive interaction due to dispersion and repulsive interaction due to anisotropic overlap of electronic orbitals. An important new feature is a dihedral-angle-dependent term that is added to the repulsive part to describe correctly several distinct stacking states that the original Kolmogorov-Crespi potential cannot distinguish. We refer to the new model as the dihedral-angle-corrected registry-dependent interlayer potential (DRIP). Computations for several test problems show that DRIP correctly reproduces the binding, sliding, and twisting energies and forces obtained from ab initio total-energy calculations based on density-functional theory. We use the new potential to study the structural properties of a twisted graphene bilayer and the exfoliation of graphene from graphite. Our potential is available through the OpenKIM interatomic potential repository at https://openkim org.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.235404},
	author = {Wen, Mingjian and Carr, Stephen and Fang, Shiang and Efthimios Kaxiras and Tadmor, Ellad B.}
}
@article { ISI:000430866200005,
	title = {Dry Dehydrogenation of Ethanol on Pt-Cu Single Atom Alloys},
	journal = {TOPICS IN CATALYSIS},
	volume = {61},
	number = {5-6, SI},
	year = {2018},
	month = {MAY},
	pages = {328-335},
	abstract = {The non-oxidative dehydrogenation of ethanol to acetaldehyde and hydrogen is an industrially relevant chemical conversion. Although Cu-based catalysts show high reactivity toward oxidative ethanol dehydrogenation, the flat Cu(111) surface is rather inactive for ethanol dehydrogenation in the absence of water, surface oxygen or defects. Herein we show, using experimental and theoretical studies of model systems, that adding 1\% Pt into the surface of Cu(111) to form dilute Pt-Cu single atom alloys (SAAs) increases the activity of Cu(111) for ethanol dehydrogenation sixfold. The mechanism of ethanol dehydrogenation was investigated at the molecular level using scanning tunneling microscopy, temperature programmed experiments and density functional theory calculations. Our results demonstrate that Pt-Cu SAAs are much more active than Cu(111) for converting ethanol to acetaldehyde and hydrogen in the absence of surface oxygen and water. Specifically, the O-H bond of ethanol is activated at Pt sites below 160 K, followed by ethoxy spillover to Cu sites which results in a significant increase of the ethoxy intermediate yield. The C-H bond of ethoxy is then activated at 310 K, and the final product, acetaldehyde, desorbs from Cu(111) in a reaction rate limited process. Finally, we show that the Cu model surfaces exhibit stability with respect to poisoning as well as 100\% selectivity in the alcohol dehydrogenation to acetaldehyde and hydrogen.},
	issn = {1022-5528},
	doi = {10.1007/s11244-017-0856-3},
	author = {Wang, Zhi-Tao and Hoyt, Robert A. and El-Soda, Mostafa and Madix, Robert J. and Efthimios Kaxiras and Sykes, E. Charles H.}
}
@article { ISI:000440820800001,
	title = {Electronic structure theory of strained two-dimensional materials with hexagonal symmetry},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {7},
	year = {2018},
	month = {AUG 6},
	abstract = {We derive electronic tight-binding Hamiltonians for strained graphene, hexagonal boron nitride, and transition-metal dichalcogenides based on Wannier transformation of ab initio density functional theory calculations. Our microscopic models include strain effects to leading order that respect the hexagonal crystal symmetry and local crystal configuration and are beyond the central force approximation which assumes only pairwise distance dependence. Based on these models, we also derive and analyze the effective low-energy Hamiltonians. Our ab initio approaches complement the symmetry group representation construction for such effective low-energy Hamiltonians and provide the values of the coefficients for each symmetry-allowed term. These models are relevant for the design of electronic device applications since they provide the framework for describing the coupling of electrons to other degrees of freedom including phonons, spin, and the electromagnetic field. The models can also serve as the basis for exploring the physics of many-body systems of interesting quantum phases.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.075106},
	author = {Fang, Shiang and Carr, Stephen and Cazalilla, Miguel A. and Efthimios Kaxiras}
}
@article { ISI:000435680100001,
	title = {Emergence and dynamical properties of stochastic branching in the electronic flows of disordered Dirac solids},
	journal = {EPL},
	volume = {122},
	number = {2},
	year = {2018},
	month = {APR},
	abstract = {Graphene as well as more generally Dirac solids constitute two-dimensional materials where the electronic flow is ultra-relativistic. When a Dirac solid is deposited on a different substrate surface with roughness, a local random potential develops through an inhomogeneous charge impurity distribution. This external potential affects profoundly the charge flow and induces a chaotic pattern of current branches that develops through focusing and defocusing effects produced by the randomness of the surface. An additional bias voltage may be used to tune the branching pattern of the charge carrier currents. We employ analytical and numerical techniques in order to investigate the onset and the statistical properties of carrier branches in Dirac solids. We find a specific scaling-type relationship that connects the physical scale for the occurrence of branches with the characteristic medium properties, such as disorder and bias field. We use numerics to test and verify the theoretical prediction as well as a perturbative approach that gives a clear indication of the regime of validity of the approach. This work is relevant to device applications and may be tested experimentally.},
	issn = {0295-5075},
	doi = {10.1209/0295-5075/122/27003},
	author = {Marios Mattheakis and G. P. Tsironis and Efthimios Kaxiras}
}
@article { ISI:000443671700003,
	title = {Energetics and kinetics of vacancy defects in 4H-SiC},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {10},
	year = {2018},
	month = {SEP 4},
	abstract = {Defect engineering in wide-gap semiconductors is important in controlling the performance of single-photon emitter devices. The effective incorporation of defects depends strongly on the ability to control their formation and location, as well as to mitigate attendant damage to the material. In this study, we combine density functional theory, molecular dyamics (MD), and kinetic Monte Carlo (KMC) simulations to study the energetics and kinetics of the silicon monovacancy V-Si and related defects in 4H-SiC. We obtain the defect formation energy for V-Si in various charge states and use MD simulations to model the ion implantation process for creating defects. We also study the effects of high-temperature annealing on defect position and stability using KMC and analytical models. Using a larger (480-atom) supercell than previous studies, we obtain the temperature-dependent diffusivity of V-Si in various charge states and find significantly lower barriers to diffusion than previous estimates. In addition, we examine the recombination with interstitial Si and conversion of V-Si into CSiVC during annealing and propose methods for using strain to reduce changes in defect concentrations. Our results provide guidance for experimental efforts to control the position and density of V-Si defects within devices, helping to realize their potential as solid-state qubits.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.104103},
	author = {Defo, Rodrick Kuate and Zhang, Xingyu and Bracher, David and Gunn Kim and Hu, Evelyn and Efthimios Kaxiras}
}
@article { ISI:000439299200003,
	title = {Enhanced superconductivity upon weakening of charge density wave transport in 2H-TaS2 in the two-dimensional limit},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {3},
	year = {2018},
	month = {JUL 20},
	abstract = {Layered transition-metal dichalcogenides that host coexisting charge-density wave (CDW) and superconducting orders provide ideal systems for exploring the effects of dimensionality on correlated electronic phases. Dimensionality has a profound effect on both superconductivity and CDW instabilities. Here we report a substantial enhancement of the superconducting T-c to 3.4 K for 2H-TaS2 in the monolayer limit, compared to 0.8 K in the bulk. In addition, the transport signature of a CDW phase transition vanishes in the two-dimensional limit. In our analysis of electronic and vibrational properties of this material, we show that a reduction of the CDW amplitude results in a substantial increase of the density of states at the Fermi energy, which can boost T, by an amount similar to that seen in experiment. Our results indicate competition between CDW order and superconductivity in ultrathin 2H-TaS2 down to the monolayer limit, providing insight toward understanding correlated electronic phases in reduced dimensions.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.035203},
	author = {Yang, Yafang and Fang, Shiang and Fatemi, Valla and Ruhman, Jonathan and Navarro-Moratalla, Efren and Kenji Watanabe and Takashi Taniguchi and Efthimios Kaxiras and Pablo Jarillo-Herrero}
}
@article { ISI:000423140600019,
	title = {General Effect of van der Waals Interactions on the Stability of Alkoxy Intermediates on Metal Surfaces},
	journal = {JOURNAL OF PHYSICAL CHEMISTRY B},
	volume = {122},
	number = {2, SI},
	year = {2018},
	month = {JAN 18},
	pages = {555-560},
	abstract = {The critical role of noncovalent van der Waals (vdW) interactions in determining the relative thermodynamic stability of alkoxy intermediates has been demonstrated for the Cu(110) surface using a combination of experiment and theory. The results may be significant for the selectivity control of copper-based reactions of alcohols. Previous examination of this effect on Au(110) was also extended to include higher molecular weight alcohols; on Cu(110) and Au(110) the hierarchy for the strength of binding of the alkoxys was found to be the same within experimental accuracy, with alkoxy species of greater chain length being more stable. The equilibrium constants governing the competition of alcohol pairs for binding sites of the alkoxys are also similar on the two surfaces. These results reveal the generality of such vdW effects. This work expands the understanding of the role of vdW interactions on the binding efficacy of key reactive intermediates on metal surfaces, a key factor in the rational design of complex and selective catalytic processes.},
	issn = {1520-6106},
	doi = {10.1021/acs.jpcb.7b04952},
	author = {Yunfei Xu and Chen, Wei and Efthimios Kaxiras and Cynthia M. Friend and Madix, Robert J.}
}
@article { ISI:000436037800046,
	title = {Heterointerface effects in the electrointercalation of van der Waals heterostructures},
	journal = {NATURE},
	volume = {558},
	number = {7710},
	year = {2018},
	month = {JUN 21},
	pages = {425+},
	abstract = {Molecular-scale manipulation of electronic and ionic charge accumulation in materials is the backbone of electrochemical energy storagel(1-4). Layered van der Waals (vdW) crystals are a diverse family of materials into which mobile ions can electrochemically intercalate into the interlamellar gaps of the host atomic lattice(5,6). The structural diversity of such materials enables the interfacial properties of composites to be optimized to improve ion intercalation for energy storage and electronic devices(7-12). However, the ability of heterolayers to modify intercalation reactions, and their role at the atomic level, are yet to be elucidated. Here we demonstrate the electrointercalation of lithium at the level of individual atomic interfaces of dissimilar vdW layers. Electrochemical devices based on vdW heterostructures(13) of stacked hexagonal boron nitride, graphene and molybdenum dichalcogenide (MoX2;X=S, Se) layers are constructed. We use transmission electron microscopy, in situ magnetoresistance and optical spectroscopy techniques, as well as low-temperature quantum magneto-oscillation measurements and ab initio calculations, to resolve the intermediate stages of lithium intercalation at heterointerfaces. The formation of vdW heterointerfaces between graphene and MoX2 results in a more than tenfold greater accumulation of charge in MoX2 when compared to MoX2/MoX2 homointerfaces, while enforcing a more negative intercalation potential than that of bulk MoX2 by at least 0.5 V. Beyond energy storage, our combined experimental and computational methodology for manipulating and characterizing the electrochemical behaviour of layered systems opens new pathways to control the charge density in two-dimensional electronic and optoelectronic devices.},
	issn = {0028-0836},
	doi = {10.1038/s41586-018-0205-0},
	author = {Bediako, D. Kwabena and Rezaee, Mehdi and Yoo, Hyobin and Larson, Daniel T. and Zhao, S. Y. Frank and Takashi Taniguchi and Kenji Watanabe and Brower-Thomas, Tina L. and Efthimios Kaxiras and Philip Kim}
}
@article { ISI:000431389300017,
	title = {Identifying key descriptors in surface binding: interplay of surface anchoring and intermolecular interactions for carboxylates on Au(110)},
	journal = {CHEMICAL SCIENCE},
	volume = {9},
	number = {15},
	year = {2018},
	month = {APR 21},
	pages = {3759-3766},
	abstract = {The relative stability of carboxylates on Au(110) was investigated as part of a comprehensive study of adsorbate binding on Group IB metals that can be used to predict and understand how to control reactivity in heterogeneous catalysis. The binding efficacy of carboxylates is only weakly dependent on alkyl chain length for relatively short-chain molecules, as demonstrated using quantitative temperature-programmed reaction spectroscopy. Corresponding density functional theory (DFT) calculations demonstrated that the bidentate anchoring geometry is rigid and restricts the amount of additional stabilization through adsorbate-surface van der Waals (vdW) interactions which control stability for alkoxides. A combination of scanning tunneling microscopy (STM) and low-energy electron diffraction (LEED) shows that carboxylates form dense local islands on Au(110). Complementary DFT calculations demonstrate that adsorbate-adsorbate interactions provide additional stabilization that increases as a function of alkyl chain length for C-2 and C-3 carboxylates. Hence, overall stability is generally a function of the anchoring group to the surface and the inter-adsorbate interaction. This study demonstrates the importance of these two important factors in describing binding of key catalytic intermediates.},
	issn = {2041-6520},
	doi = {10.1039/c7sc05313d},
	author = {O{\textquoteright}Connor, Christopher R. and Fanny Hiebel and Chen, Wei and Efthimios Kaxiras and Madix, Robert J. and Cynthia M. Friend}
}
@article { ISI:000427815800001,
	title = {Implanted neural network potentials: Application to Li-Si alloys},
	journal = {PHYSICAL REVIEW B},
	volume = {97},
	number = {9},
	year = {2018},
	month = {MAR 20},
	abstract = {Modeling the behavior of materials composed of elements with different bonding and electronic structure character for large spatial and temporal scales and over a large compositional range is a challenging problem. Cases in point are amorphous alloys of Si, a prototypical covalent material, and Li, a prototypical metal, which are being considered as anodes for high-energy-density batteries. To address this challenge, we develop a methodology based on neural networks that extends the conventional training approach to incorporate pre-trained parts that capture the character of different components, into the overall network; we refer to this model as the {\textquoteleft}{\textquoteleft}implanted neural network{\textquoteright}{\textquoteright} method. We show that this approach works well for the Si-Li amorphous alloys for a wide range of compositions, giving good results for key quantities like the diffusion coefficients. The method is readily generalizable to more complicated situations that involve two or more different elements.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.97.094106},
	author = {Onat, Berk and Ekin D. Cubuk and Malone, Brad D. and Efthimios Kaxiras}
}
@article { ISI:000449308400005,
	title = {Lithium Intercalation in Graphene-MoS2 Heterostructures},
	journal = {JOURNAL OF PHYSICAL CHEMISTRY C},
	volume = {122},
	number = {43},
	year = {2018},
	month = {NOV 1},
	pages = {24535-24541},
	abstract = {Two-dimensional (2D) heterostructures are interesting candidates for efficient energy storage devices due to their high carrier capacity by reversible intercalation. We employ here density functional theory calculations to investigate the structural and electronic properties of lithium intercalated graphene/molybdenum disulfide (Gr/MoS2) heterostructures. We explore the extent to which Li intercalates at the interface formed between graphene (Gr) and molybdenum disulfide (MoS2) layers by considering the adsorption and diffusion of Li atoms, the energetic stability, and the changes in the structural morphology of MoS2. We investigate the corresponding electronic structure and charge distribution within the heterostructure at varying concentrations of Li. Our results indicate that the maximum energetically allowed ratio of Li to Mo (Li to C) is 1:1 (1:3) for both the 2H and 1T{\textquoteright} phases of MoS2. This is double the Li concentration allowed in graphene bilayers. We find that there is 60\% more charge transfer to MoS2 than to Gr in the bilayer heterostructure, which results in a maximum doping of Gr and MoS2 of n(C) = 3.6 x 10(14) cm(-2) and n(MoS2) = 6.0 x 10(14) cm(-2), respectively.},
	issn = {1932-7447},
	doi = {10.1021/acs.jpcc.8b07548},
	author = {Larson, Daniel T. and Ioanna Fampiou and Gunn Kim and Efthimios Kaxiras}
}
@conference { ISI:000451429100013,
	title = {Modeling Electronic Properties of Twisted 2D Atomic Heterostructures},
	booktitle = {COUPLED MATHEMATICAL MODELS FOR PHYSICAL AND BIOLOGICAL NANOSCALE SYSTEMS AND THEIR APPLICATIONS},
	series = {Springer Proceedings in Mathematics \& Statistics},
	volume = {232},
	year = {2018},
	note = {Banff International Research Station (BIRS) Workshop on Coupled Mathematical Models for Physical and Biological Nanoscale Systems and Their Applications, Banff, CANADA, AUG 28-SEP 02, 2016},
	pages = {245-265},
	publisher = {Banff Int Res Stat},
	organization = {Banff Int Res Stat},
	abstract = {We present a general method for the electronic characterization of aperiodic 2D materials using ab-initio tight binding models. Specifically studied is the subclass of twisted, stacked heterostructures, but the formalism provided can be implemented for any 2D system without long-range interactions. This new method provides a multi-scale approach for dealing with the ab-initio calculation of electronic transport properties in stacked nanomaterials, allowing for fast and efficient simulation of multi-layered stacks in the presence of twist angles, magnetic field, and defects. We calculate the electronic density of states in twisted bilayer systems of graphene and MX2 transition metal dichalcogenides (TMDCs). We comment on the interesting features of their density of states as a function of twist-angle and local configuration and how these features are experimentally observable. These results support the bilayer twist-angle as a new variable for controlling electronic properties in artificial nanomaterials ({\textquoteright}{\textquoteright}Twistronics{\textquoteright}{\textquoteright}).},
	isbn = {978-3-319-76599-0; 978-3-319-76598-3},
	issn = {2194-1009},
	doi = {10.1007/978-3-319-76599-0\_13},
	author = {Carr, Stephen and Massatt, Daniel and Fang, Shiang and Cazeaux, Paul and Luskin, Mitchell and Efthimios Kaxiras},
	editor = {Bonilla, LL and E. Kaxiras and Melnik, R}
}
@article { ISI:000444793500098,
	title = {New Pathway for Hot Electron Relaxation in Two-Dimensional Heterostructures},
	journal = {NANO LETTERS},
	volume = {18},
	number = {9},
	year = {2018},
	month = {SEP},
	pages = {6057-6063},
	abstract = {Two-dimensional (2D) heterostructures composed of transition-metal dichalcogenide atomic layers are the new frontier for novel optoelectronic and photovoltaic device applications. Some key properties that make these materials appealing, yet are not well understood, are ultrafast hole/electron dynamics, interlayer energy transfer and the formation of interlayer hot excitons. Here, we study photoexcited electron/hole dynamics in a representative heterostructure, the MoS2/WSe2 interface, which exhibits type II band alignment. Employing time-dependent density functional theory in the time domain, we observe ultrafast charge dynamics with lifetimes of tens to hundreds of femtoseconds. Most importantly, we report the discovery of an interfacial pathway in 2D heterostructures for the relaxation of photoexcited hot electrons through interlayer hopping, which is significantly faster than intralayer relaxation. This finding is of particular importance for understanding many experimentally observed photoinduced processes, including charge and energy transfer at an ultrafast time scale (,
	issn = {1530-6984},
	doi = {10.1021/acs.nanolett.8b03005},
	author = {Zhang, Jin and Hong, Hao and Zhang, Jia and Fu, Huixia and You, Peiwei and Lischner, Johannes and Liu, Kaihui and Efthimios Kaxiras and Sheng Meng}
}
@article { ISI:000445713200021,
	title = {Nonadiabatic Hydrogen Dissociation on Copper Nanoclusters},
	journal = {JOURNAL OF PHYSICAL CHEMISTRY LETTERS},
	volume = {9},
	number = {18},
	year = {2018},
	month = {SEP 20},
	pages = {5339+},
	abstract = {Copper surfaces exhibit high catalytic selectivity but have poor hydrogen dissociation kinetics; therefore, we consider icosahedral Cu-13 nanoclusters to understand how nanoscale structure might improve catalytic prospects. We find that the spin state is a surprisingly important design consideration. Cu-13 clusters have large magnetic moments due to finite size and symmetry effects and exhibit magnetization-dependent catalytic behavior. The most favorable transition state for hydrogen dissociation has a lower activation energy than that on single-crystal copper surfaces but requires a magnetization switch from 5 to 3 mu(B). Without this switch, the activation energy is higher than that on single-crystal surfaces. Weak spin-orbit coupling hinders this switch, decreasing the kinetic rate of hydrogen dissociation by a factor of 16. We consider strategies to facilitate magnetization switches through optical excitations, substitution, charge states, and co-catalysts; these considerations demonstrate how control of magnetic properties could improve catalytic performance.},
	issn = {1948-7185},
	doi = {10.1021/acs.jpclett.8b02133},
	author = {Hoyt, Robert A. and Montemore, Matthew M. and Efthimios Kaxiras}
}
@article { ISI:000442875900001,
	title = {Pressure dependence of the magic twist angle in graphene superlattices},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {8},
	year = {2018},
	month = {AUG 28},
	abstract = {The recently demonstrated unconventional superconductivity [Cao et al., Nature (London) 556, 43 (2018)] in twisted bilayer graphene (tBLG) opens the possibility for interesting applications of two-dimensional layers that involve correlated electron states. Here we explore the possibility of modifying electronic correlations by the application of uniaxial pressure on the weakly interacting layers, which results in increased interlayer coupling and a modification of the magic angle value and associated density of states. Our findings are based on first-principles calculations that accurately describe the height-dependent interlayer coupling through the combined use of density functional theory and maximally localized Wannier functions. We obtain the relationship between twist angle and external pressure for the magic angle flat bands of tBLG. This may provide a convenient method to tune electron correlations by controlling the length scale of the superlattice.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.085144},
	author = {Carr, Stephen and Fang, Shiang and Pablo Jarillo-Herrero and Efthimios Kaxiras}
}
@article { ISI:000433012900007,
	title = {Quantum plasmons with optical-range frequencies in doped few-layer graphene},
	journal = {PHYSICAL REVIEW B},
	volume = {97},
	number = {19},
	year = {2018},
	month = {MAY 21},
	abstract = {Although plasmon modes exist in doped graphene, the limited range of doping achieved by gating restricts the plasmon frequencies to a range that does not include the visible and infrared. Here we show, through the use of first-principles calculations, that the high levels of doping achieved by lithium intercalation in bilayer and trilayer graphene shift the plasmon frequencies into the visible range. To obtain physically meaningful results, we introduce a correction of the effect of plasmon interaction across the vacuum separating periodic images of the doped graphene layers, consisting of transparent boundary conditions in the direction perpendicular to the layers; this represents a significant improvement over the exact Coulomb cutoff technique employed in earlier works. The resulting plasmon modes are due to local field effects and the nonlocal response of the material to external electromagnetic fields, requiring a fully quantum mechanical treatment. We describe the features of these quantum plasmons, including the dispersion relation, losses, and field localization. Our findings point to a strategy for fine-tuning the plasmon frequencies in graphene and other two-dimensional materials.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.97.195435},
	author = {Shirodkar, Sharmila N. and Marios Mattheakis and Cazeaux, Paul and Prineha Narang and Marin Soljacic and Efthimios Kaxiras}
}
@article { ISI:000440717600005,
	title = {Raman spectrum of CrI3: An ab initio study},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {8},
	year = {2018},
	month = {AUG 3},
	abstract = {We study the Raman spectrum of CrI3, a material that exhibits magnetism in a single layer. We employ first-principles calculations within density functional theory to determine the effects of polarization, strain, and incident angle on the phonon spectra of the three-dimensional bulk and the single-layer two-dimensional structure, for both the high- and low-temperature crystal structures. Our results are in good agreement with existing experimental measurements and serve as a guide for additional investigations to elucidate the physics of this interesting material.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.085406},
	author = {Larson, Daniel T. and Efthimios Kaxiras}
}
@article { ISI:000449723900049,
	title = {Reaction-Induced Excitations and Their Effect on Surface Chemistry},
	journal = {ACS CATALYSIS},
	volume = {8},
	number = {11},
	year = {2018},
	month = {NOV},
	pages = {10358-10363},
	abstract = {Despite intensive study of reactions on metals, it is unclear whether electronic excitations play an important role. Here, we show that nonadiabatic effects do indeed play a significant role in N-2 and H-2 dissociation on Ru nanoparticles. We employ nonadiabatic dynamical calculations based on realtime, time-dependent density functional theory to study energy dissipation during these exothermic reaction steps. We find that dissipation of the excess energy into excitation of electrons exceeds thermal dissipation into phonons. For isolated dissociation events, electronic friction can increase reaction barriers; furthermore, the excitations induced by a dissociation event can affect other reacting molecules. Our studies suggest that, for exothermic reactions, metal catalysts in reaction conditions may be constantly experiencing electronic excitations, and these excitations can significantly affect surface chemistry.},
	issn = {2155-5435},
	doi = {10.1021/acscatal.8b03266},
	author = {Montemore, Matthew M. and Robert Hoyt and Grigory Kolesov and Efthimios Kaxiras}
}
@article { ISI:000452322000001,
	title = {Relaxation and domain formation in incommensurate two-dimensional heterostructures},
	journal = {PHYSICAL REVIEW B},
	volume = {98},
	number = {22},
	year = {2018},
	month = {DEC 5},
	abstract = {We introduce configuration space as a natural representation for calculating the mechanical relaxation patterns of incommensurate two-dimensional (2D) bilayers. The approach can be applied to a wide variety of 2D materials through the use of a continuum model in combination with a generalized stacking fault energy for interlayer interactions. We present computational results for small-angle twisted bilayer graphene and molybdenum disulfide (MoS2), a representative material of the transition-metal dichalcogenide family of 2D semiconductors. We calculate accurate relaxations for MoS2 even at small twist-angle values, enabled by the fact that our approach does not rely on empirical atomistic potentials for interlayer coupling. The results demonstrate the efficiency of the configuration space method by computing relaxations with minimal computational cost. We also outline a general explanation of domain formation in 2D bilayers with nearly aligned lattices, taking advantage of the relationship between real space and configuration space. The configuration space approach also enables calculation of relaxations in incommensurate multilayer systems.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.98.224102},
	author = {Carr, Stephen and Massatt, Daniel and Torrisi, Steven B. and Cazeaux, Paul and Luskin, Mitchell and Efthimios Kaxiras}
}
@article { ISI:000429103300032,
	title = {Unconventional superconductivity in magic-angle graphene superlattices},
	journal = {NATURE},
	volume = {556},
	number = {7699},
	year = {2018},
	month = {APR 5},
	pages = {43+},
	abstract = {The behaviour of strongly correlated materials, and in particular unconventional superconductors, has been studied extensively for decades, but is still not well understood. This lack of theoretical understanding has motivated the development of experimental techniques for studying such behaviour, such as using ultracold atom lattices to simulate quantum materials. Here we report the realization of intrinsic unconventional superconductivity-which cannot be explained by weak electron-phonon interactions-in a two-dimensional superlattice created by stacking two sheets of graphene that are twisted relative to each other by a small angle. For twist angles of about 1.1 degrees-the first {\textquoteleft}magic{\textquoteright} angle-the electronic band structure of this {\textquoteleft}twisted bilayer graphene{\textquoteright} exhibits flat bands near zero Fermi energy, resulting in correlated insulating states at half-filling. Upon electrostatic doping of the material away from these correlated insulating states, we observe tunable zero-resistance states with a critical temperature of up to 1.7 kelvin. The temperature-carrier-density phase diagram of twisted bilayer graphene is similar to that of copper oxides (or cuprates), and includes dome-shaped regions that correspond to superconductivity. Moreover, quantum oscillations in the longitudinal resistance of the material indicate the presence of small Fermi surfaces near the correlated insulating states, in analogy with underdoped cuprates. The relatively high superconducting critical temperature of twisted bilayer graphene, given such a small Fermi surface (which corresponds to a carrier density of about 1011 per square centimetre), puts it among the superconductors with the strongest pairing strength between electrons. Twisted bilayer graphene is a precisely tunable, purely carbon-based, two-dimensional superconductor. It is therefore an ideal material for investigations of strongly correlated phenomena, which could lead to insights into the physics of high-critical-temperature superconductors and quantum spin liquids.},
	issn = {0028-0836},
	doi = {10.1038/nature26160},
	author = {Cao, Yuan and Fatemi, Valla and Fang, Shiang and Kenji Watanabe and Takashi Taniguchi and Efthimios Kaxiras and Pablo Jarillo-Herrero}
}
@article { ISI:000423427400002,
	title = {Universal behavior of a dispersive Dirac cone in gradient-index plasmonic metamaterials},
	journal = {PHYSICAL REVIEW B},
	volume = {97},
	number = {3},
	year = {2018},
	month = {JAN 25},
	abstract = {We demonstrate analytically and numerically that the dispersive Dirac cone emulating an epsilon-near-zero (ENZ) behavior is a universal property within a family of plasmonic crystals consisting of two-dimensional (2D) metals. Our starting point is a periodic array of 2D metallic sheets embedded in an inhomogeneous and anisotropic dielectric host that allows for propagation of transverse-magnetic (TM) polarized waves. By invoking a systematic bifurcation argument for arbitrary dielectric profiles in one spatial dimension, we show how TM Bloch waves experience an effective dielectric function that averages out microscopic details of the host medium. The corresponding effective dispersion relation reduces to a Dirac cone when the conductivity of the metallic sheet and the period of the array satisfy a critical condition for ENZ behavior. Our analytical findings are in excellent agreement with numerical simulations.},
	issn = {2469-9950},
	doi = {10.1103/PhysRevB.97.035307},
	author = {Matthias Maier and Marios Mattheakis and Efthimios Kaxiras and Luskin, Mitchell and Dionisios Margetis}
}
@article { ISI:000423505500007,
	title = {Water facilitates oxygen migration on gold surfaces},
	journal = {PHYSICAL CHEMISTRY CHEMICAL PHYSICS},
	volume = {20},
	number = {4},
	year = {2018},
	month = {JAN 28},
	pages = {2196-2204},
	abstract = {The water-oxygen-gold interface is important in many surface processes and has potential influence on heterogeneous catalysis. Herein, it is shown that water facilitates the migration of atomic oxygen on Au(110), demonstrating the dynamic nature of surface adsorption. We demonstrate this effect for the first time, using in situ scanning tunnelling microscopy (STM), temperature-programmed reaction spectroscopy (TPRS) and first-principles theoretical calculations. The dynamic interaction of water with adsorbed O maintains a high dispersion of O on the surface, potentially creating reactive transient species. At low temperature and pressure, isotopic experiments show that adsorbed oxygen on the Au(110) surface exchanges with oxygen in (H2O)-O-18. The presence of water modulates local electronic properties and facilitates oxygen exchange. Combining experimental results and theory, we propose that hydroxyl is transiently formed via proton transfer from the water to adsorbed oxygen. Hydroxyl groups easily recombine to regenerate water and adsorbed oxygen atoms, the net result of which is migration of the adsorbed oxygen without significant change in its overall distribution on the surface. The presence of water creates a dynamic surface where mobile surface oxygen atoms and hydroxyls are present, which can lead to a better performance of gold catalysis in oxidation reactions.},
	issn = {1463-9076},
	doi = {10.1039/c7cp06451a},
	author = {Fang Xu and Ioanna Fampiou and O{\textquoteright}Connor, Christopher R. and Stavros Karakalos and Fanny Hiebel and Efthimios Kaxiras and Madix, Robert J. and Cynthia M. Friend}
}
