Graeme Henkelman

Find an error

Name: Henkelman, Graeme
Organization: The University of Texas at Austin , USA
Department: Department of Chemistry and the Institute for Computational Engineering and Sciences
Title: (PhD)
Co-reporter:Jie Lin, Jin-Myoung Lim, Duck Hyun Youn, Kenta Kawashima, Jun-Hyuk Kim, Yang Liu, Hang Guo, Graeme Henkelman, Adam Heller, and Charles Buddie Mullins
ACS Nano October 24, 2017 Volume 11(Issue 10) pp:10347-10347
Publication Date(Web):September 12, 2017
DOI:10.1021/acsnano.7b05294
Through a gelation–solvothermal method without heteroadditives, Cu–Sn–S composites self-assemble to form nanotubes, sub-nanotubes, and nanoparticles. The nanotubes with a Cu3–4SnS4 core and Cu2SnS3 shell can tolerate long cycles of expansion/contraction upon lithiation/delithiation, retaining a charge capacity of 774 mAh g–1 after 200 cycles with a high initial Coulombic efficiency of 82.5%. The importance of the Cu component for mitigation of the volume expansion and structural evolution upon lithiation is informed by density functional theory calculations. The self-generated template and calculated results can inspire the design of analogous Cu–M–S (M = metal) nanotubes for lithium batteries or other energy storage systems.Keywords: copper tin sulfide; core−shell; density functional theory; gelation−solvothermal; lithium battery; nanotube;
Co-reporter:Sarah Seraj, Pranaw Kunal, Hao Li, Graeme Henkelman, Simon M. Humphrey, and Charles J. Werth
ACS Catalysis May 5, 2017 Volume 7(Issue 5) pp:3268-3268
Publication Date(Web):March 24, 2017
DOI:10.1021/acscatal.6b03647
Well-defined palladium–gold nanoparticles (PdAuNPs) with randomly alloyed structures and broadly tunable compositions were studied in catalytic nitrite (NO2–) reduction. The catalysts were synthesized using a microwave-assisted polyol coreduction method. PdxAu100–xNPs with systematically varied compositions (x = 18–83) were supported on amorphous silica (SiO2) and studied as model catalysts for aqueous NO2– reduction in a batch reactor, using H2 as the electron donor. The reactions followed pseudo-first-order kinetics for ≥80% NO2– conversion. The PdxAu100–xNP-SiO2 catalysts showed a volcano-like correlation between NO2– reduction activity and x; the highest activity was observed for Pd53Au47, with an associated first-order rate constant of 5.12 L min–1 gmetal–1. Alloy NPs with greater proportions of Au were found to reduce the loss in catalytic activity due to sulfide fouling. Density functional theory calculations indicate that this is because Au weakens sulfur binding at PdAuNP surfaces due to atomic ensemble, electronic, and strain effects and thus reduces sulfur poisoning. The environmental relevance of the most active supported catalyst was evaluated by subjecting it to five cycles of catalytic NO2– reduction. The catalytic activity decreased over multiple cycles, but analysis of the postreaction PdxAu100–xNP-SiO2 materials using complementary techniques indicated that there were no significant structural changes. Most importantly, we show that PdxAu100–xNP-SiO2 alloys are significantly more active NO2– reduction catalysts in comparison to pure Pd catalysts.Keywords: density functional theory (DFT); heterogeneous catalysis; microwave synthesis; nitrite hydrogenation; palladium−gold alloys; sulfide poisoning; water treatment;
Co-reporter:Penghao Xiao;Renlong Liu;Xin Lian
Chemistry of Materials February 28, 2017 Volume 29(Issue 4) pp:1472-1484
Publication Date(Web):January 27, 2017
DOI:10.1021/acs.chemmater.6b02722
Co-reporter:Hao Li and Graeme Henkelman
The Journal of Physical Chemistry C December 14, 2017 Volume 121(Issue 49) pp:27504-27504
Publication Date(Web):November 20, 2017
DOI:10.1021/acs.jpcc.7b09953
Ethanol (EtOH) decomposition has been widely studied in recent years. However, the initial dehydrogenation selectivity on catalytic surfaces, which plays a crucial role in EtOH partial oxidation and steam reforming, is not well understood. Here, density functional theory (DFT) was used to calculate the initial dehydrogenation selectivities of EtOH on monometallic and X/Au (X = Pd and Rh) close-packed surfaces. The energy for the initial bond scissions of O–H and α- and β-C–H were calculated on each surface. The binding energy of EtOH is found to be a good reactivity descriptor for the scission of O–H and β-C–H bonds, while the binding energy of CH3CHOH is a good reaction descriptor for α-C–H bond scission. The scaling relationships between the activation energy barriers and binding energies on Pd/Au and Rh/Au surface alloys are significantly different from those of monometallic surfaces. Additionally, the specific atomic ensembles on the Pd/Au and Rh/Au surfaces have different initial dehydrogenation selectivities of EtOH. Our calculated scaling relationships were used to construct contour plots that provide predictive trends for the selectivity of the initial EtOH dehydrogenation. We conclude that the presence of specific atomic ensembles on the surface of alloy catalysts can efficiently control the reaction products of EtOH dehydrogenation.
Co-reporter:Long Luo, Zhiyao Duan, Hao Li, Joohoon Kim, Graeme Henkelman, and Richard M. Crooks
Journal of the American Chemical Society April 19, 2017 Volume 139(Issue 15) pp:5538-5538
Publication Date(Web):April 7, 2017
DOI:10.1021/jacs.7b01653
In this paper, we show that PtAu and PdAu random alloy dendrimer-encapsulated nanoparticles with an average size of ∼1.6 nm have different catalytic activity trends for allyl alcohol hydrogenation. Specifically, PtAu nanoparticles exhibit a linear increase in activity with increasing Pt content, whereas PdAu dendrimer-encapsulated nanoparticles show a maximum activity at a Pd content of ∼60%. Both experimental and theoretical results suggest that this contrasting behavior is caused by differences in the strength of H binding on the PtAu and PdAu alloy surfaces. The results have significant implications for predicting the catalytic performance of bimetallic nanoparticles on the basis of density functional theory calculations.
Co-reporter:Qing Chen;Penghao Xiao;Yi Pei;Yan Song;Cheng-Yan Xu;Liang Zhen
Journal of Materials Chemistry A 2017 vol. 5(Issue 32) pp:16722-16731
Publication Date(Web):2017/08/15
DOI:10.1039/C7TA03049E
Li2MnSiO4 is a promising high capacity cathode material due to the potential to extract two Li ions per formula unit. In practice, however, the use of Li2MnSiO4 is restricted by a low discharge capacity, which has been attributed to an irreversible structural change in the first charge cycle. In this work, we use density functional theory calculations to explore the details of this structural change, and our results reveal that the structural change during delithiation has two components. First, we find that the material undergoes a structural collapse upon partial delithiation, which is characterized by distortion of the MnO4 tetrahedron. Remarkably, while this transformation results in a disordered structure, our calculations show that it is reversible upon relithiation and that the transformation does not strongly impede Li de/intercalation. The calculated reversibility of the phase change is consistent with recent experimental X-ray diffraction measurements showing that peaks associated with the crystalline MnO4 order, which disappear upon delithiation, are restored upon lithiation. Additional experiments are conducted showing the reversibility of the material during cycling as a function of charging cutoff voltages. Second, we argue that, the irreversible structural degradation is primarily caused by oxygen evolution in the highly delithiated state; the oxygen deficient structure can only reincorporate half of the total Li when discharged to 1.5 V. Experimentally observed voltage profile shifts of Li2MnSiO4 during the first few cycles as well as the different electrochemical behavior exhibited by Li2FeSiO4 can be explained by this two-component structural change model.
Co-reporter:Zhiyao Duan and Graeme Henkelman  
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 7) pp:5486-5490
Publication Date(Web):01 Feb 2016
DOI:10.1039/C5CP05558J
Density functional theory calculations reveal that the work function of Au supported on MgO(001) is substantially reduced because of an interfacial dipole moment formed at the Au/MgO interface. Consequently, the Au/MgO interface plays an active role in the activation of O2 molecules by promoting charge transfer to the O2 2π* orbital. The presence of F-centers in the MgO substrate can further promote the charge transfer and bonding of O2 at the interface boundary. However, O2 dissociation is kinetically hindered. The system is then able to catalyze CO oxidation at low temperature as adsorbed CO and O2 readily react to form CO2 with a low energy barrier.
Co-reporter:Bing Liu
The Journal of Physical Chemistry C 2016 Volume 120(Issue 10) pp:5557-5564
Publication Date(Web):February 24, 2016
DOI:10.1021/acs.jpcc.6b00253
Engineering a bimetallic system with complementary chemical properties can be an effective way of tuning catalytic activity. In this work, CO oxidation on CeO2(111)-supported Pd-based bimetallic nanorods was investigated using density functional theory calculations corrected by on-site Coulomb interactions. We studied a series of CeO2(111)-supported Pd-based bimetallic nanorods (Pd–X, where X = Ag, Au, Cu, Pt, Rh, Ru) and found that Pd–Ag/CeO2 and Pd–Cu/CeO2 are the two systems where the binding sites of CO and O2 are distinct; that is, in these two systems, CO and O2 do not compete for the same binding sites. An analysis of the CO oxidation mechanisms suggests that the Pd–Ag/CeO2 system is more effective for catalyzing CO oxidation as compared to Pd–Cu/CeO2 because both CeO2 lattice oxygen atoms and adsorbed oxygen molecules at Ag sites can oxidize CO with low energy barriers. Both the Pd–Ag and Pd–CeO2 interfaces in Pd–Ag/CeO2 were found to play important roles in CO oxidation. The Pd–Ag interface, which combines the different chemical nature of the two metals, not only separates the binding sites of CO and O2 but also opens up active reaction pathways for CO oxidation. The strong metal–support interaction at the Pd–CeO2 interface facilitates CO oxidation by the Mars–van Krevelen mechanism. Our study provides theoretical guidance for designing highly active metal/oxide catalysts for CO oxidation.
Co-reporter:Long Luo, Liang Zhang, Zhiyao Duan, Aliya S. Lapp, Graeme Henkelman, and Richard M. Crooks
ACS Nano 2016 Volume 10(Issue 9) pp:8760
Publication Date(Web):September 1, 2016
DOI:10.1021/acsnano.6b04448
In this paper, we show that the onset potential for CO oxidation electrocatalyzed by ∼2 nm dendrimer-encapsulated Pt nanoparticles (Pt DENs) is shifted negative by ∼300 mV in the presence of a small percentage (<2%) of Cu surface atoms. Theory and experiments suggest that the catalytic enhancement arises from a cocatalytic Langmuir–Hinshelwood mechanism in which the small number of Cu atoms selectively adsorb OH, thereby facilitating reaction with CO adsorbed to the dominant Pt surface. Theory suggests that these Cu atoms are present primarily on the (100) facets of the Pt DENs.Keywords: bimetallic nanoparticle; CO electro-oxidation; density functional theory; isolated Cu
Co-reporter:Zhiyao Duan and Graeme Henkelman
ACS Catalysis 2015 Volume 5(Issue 3) pp:1589
Publication Date(Web):January 29, 2015
DOI:10.1021/cs501610a
Density functional theory is used to determine the reaction mechanisms of CO oxidation and the active oxygen species on a Au/TiO2 model catalyst. The model consists of a Au rod supported along the TiO2 [11̅0] direction of the TiO2(110) surface. An interfacial Au/Ti5c site at the interface boundary is identified to be particularly active toward O2 adsorption and dissociation. At this site, O2 dissociation has an energy barrier of 0.5 eV, which is facile at room temperature. The resulting adsorbed Au/O/Ti5c oxygen species are shown to be stable and active for CO oxidation under relevant reaction conditions with an activation energy of 0.24 eV. Furthermore, the adsorbed Au/O/Ti5c oxygen species functions as an electron reservoir, and it lowers the oxygen vacancy formation energy of a surface lattice oxygen (Obri), as well as the Ti interstitial formation energy, due to electron transfer from high-energy defect states to low-energy p-states of the adsorbed Au/O/Ti5c oxygen species. Hence, the Obri species is activated at the oxidized Au/TiO2 interface boundary and the energy barrier of CO oxidation with Obri is calculated to be 0.55 eV. Thus, the CO oxidation reaction can proceed at room temperature either via a Langmuir–Hinshelwood mechanism with an adsorbed Au/O/Ti5c oxygen species or via a Au-assisted Mars–van Krevelen mechanism with Obri.Keywords: Au/TiO2; CO oxidation; density functional theory; gold catalysis; heterogeneous catalysis
Co-reporter:Liang Zhang and Graeme Henkelman
ACS Catalysis 2015 Volume 5(Issue 2) pp:655
Publication Date(Web):December 12, 2014
DOI:10.1021/cs501176b
The alloy-core@shell nanoparticle structure combines the advantages of a robust noble-metal shell and a tunable alloy-core composition. In this study we demonstrate a set of linear correlations between the binding of adsorbates to the shell and the alloy-core composition, which are general across a range of nanoparticle compositions, size, and adsorbate molecules. This systematic tunability allows for a simple approach to the design of such catalysts. Calculations of candidate structures for the hydrogen evolution reaction predict a high activity for the PtRu@Pd structure, in good agreement with what has been reported previously. Calculations of alloy-core@Pt 140-atom nanoparticles reveal new candidate structures for CO oxidation at high temperature, including Au0.65Pd0.35@Pt and Au0.73Pt0.27@Pt, which are predicted to have reaction rates 200 times higher than that of Pt(111).Keywords: alloy-core/shell; binding energy; catalyst design; CO oxidation; hydrogen evolution; nanoparticles
Co-reporter:Penghao Xiao, Jie Song, Long Wang, John B. Goodenough, and Graeme Henkelman
Chemistry of Materials 2015 Volume 27(Issue 10) pp:3763
Publication Date(Web):May 12, 2015
DOI:10.1021/acs.chemmater.5b01132
The Prussian Blue analog, NaxFeMn(CN)6, is a potential new cathode material for Na-ion batteries. During Na intercalation, the dehydrated material exhibits a monoclinic to rhombohedral phase transition, while the hydrated material remains in the monoclinic phase. With density functional theory calculations, the phase transition is explained in terms of a competition between Coulomb attraction, Pauli repulsion, and d–π covalent bonding. The interstitial Na cations have a strong Coulomb attraction to the N anions in the host material, which tend to bend the Mn–N bonds and reduce the volume of the structure. The presence of lattice H2O enhances the Pauli repulsion so that the volume reduction is suppressed. The calculated volume change, as it depends upon the presence of lattice H2O, is consistent with experimental measurements. Additionally, a new LiFeMn(CN)6 phase is predicted where MnN6 octahedra decompose into LiN4 and MnN4 edge-sharing tetrahedra.
Co-reporter:Samuel T. Chill, Rachel M. Anderson, David F. Yancey, Anatoly I. Frenkel, Richard M. Crooks, and Graeme Henkelman
ACS Nano 2015 Volume 9(Issue 4) pp:4036
Publication Date(Web):April 8, 2015
DOI:10.1021/acsnano.5b00090
We present a method for quantifying the accuracy of extended X-ray absorption fine structure (EXAFS) fitting models. As a test system, we consider the structure of bare Au147 nanoparticles as well as particles bound with thiol ligands, which are used to systematically vary disorder in the atomic structure of the nanoparticles. The accuracy of the fitting model is determined by comparing two distributions of bond lengths: (1) a direct average over a molecular dynamics (MD) trajectory using forces and energies from density functional theory (DFT) and (2) a fit to the theoretical EXAFS spectra generated from that same trajectory. Both harmonic and quasi-harmonic EXAFS fitting models are used to characterize the first-shell Au–Au bond length distribution. The harmonic model is found to significantly underestimate the coordination number, disorder, and bond length. The quasi-harmonic model, which includes the third cumulant of the first-shell bond length distribution, yields accurate bond lengths, but incorrectly predicts a decrease in particle size and little change in the disorder with increasing thiol ligands. A direct analysis of the MD data shows that the particle surfaces become much more disordered with ligand binding, and the high disorder is incorrectly interpreted by the EXAFS fitting models. Our DFT calculations compare well with experimental EXAFS measurements of Au nanoparticles, synthesized using a dendrimer encapsulation technique, showing that systematic errors in EXAFS fitting models apply to nanoparticles 1–2 nm in size. Finally we show that a combination of experimental EXAFS analysis with candidate models from DFT is a promising strategy for a more accurate determination of nanoparticle structures.Keywords: Au nanoparticles; density functional theory; extended X-ray absorption fine structure; structural disorder;
Co-reporter:Samuel T. Chill, Jacob Stevenson, Victor Ruehle, Cheng Shang, Penghao Xiao, James D. Farrell, David J. Wales, and Graeme Henkelman
Journal of Chemical Theory and Computation 2014 Volume 10(Issue 12) pp:5476-5482
Publication Date(Web):November 11, 2014
DOI:10.1021/ct5008718
A set of benchmark systems is defined to compare different computational approaches for characterizing local minima, transition states, and pathways in atomic, molecular, and condensed matter systems. Comparisons between several commonly used methods are presented. The strengths and weaknesses are discussed, as well as implementation details that are important for achieving good performance. All of the benchmarks and methods are provided in an online database to make the implementation details available and the results reproducible. While this paper provides a snapshot of the benchmark results, the online framework is structured to be dynamic and incorporate new methods and codes as they are developed.
Co-reporter:Zhiyao Duan and Graeme Henkelman
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 7) pp:NaN5490-5490
Publication Date(Web):2016/02/01
DOI:10.1039/C5CP05558J
Density functional theory calculations reveal that the work function of Au supported on MgO(001) is substantially reduced because of an interfacial dipole moment formed at the Au/MgO interface. Consequently, the Au/MgO interface plays an active role in the activation of O2 molecules by promoting charge transfer to the O2 2π* orbital. The presence of F-centers in the MgO substrate can further promote the charge transfer and bonding of O2 at the interface boundary. However, O2 dissociation is kinetically hindered. The system is then able to catalyze CO oxidation at low temperature as adsorbed CO and O2 readily react to form CO2 with a low energy barrier.
Co-reporter:Zhiyao Duan, Yuanyuan Li, Janis Timoshenko, Samuel T. Chill, Rachel M. Anderson, David F. Yancey, Anatoly I. Frenkel, Richard M. Crooks and Graeme Henkelman
Catalysis Science & Technology (2011-Present) 2016 - vol. 6(Issue 18) pp:NaN6885-6885
Publication Date(Web):2016/05/26
DOI:10.1039/C6CY00559D
In this study, we present a framework for characterizing the structural and thermal properties of small nanoparticle catalysts by combining precise synthesis, extended X-ray absorption fine structure (EXAFS) spectroscopy, and density functional theory (DFT) calculations. We demonstrate the capability of this approach by characterizing the atomic structure and vibrational dynamics of Au147. With the combination of EXAFS spectroscopy and DFT, the synthesized Au147 nanoparticles are determined to have an icosahedral structure. A decrease in the Einstein temperature of the Au147 particles compared to their bulk value was observed and interpreted in terms of softer vibration modes of surface bonds.
Ferrate(3-),hexakis(cyano-kC)-,(OC-6-11)- (9CI)
Oxygen, isotope of mass18, at.
Hydroxyl