1 |
|
\documentclass[11pt]{article} |
2 |
|
\usepackage{amsmath} |
3 |
|
\usepackage{amssymb} |
4 |
+ |
\usepackage{times} |
5 |
+ |
\usepackage{mathptm} |
6 |
|
\usepackage{setspace} |
7 |
< |
\usepackage{endfloat} |
7 |
> |
\usepackage{float} |
8 |
|
\usepackage{caption} |
9 |
|
|
10 |
|
%\usepackage{tabularx} |
13 |
|
%\usepackage{booktabs} |
14 |
|
%\usepackage{bibentry} |
15 |
|
%\usepackage{mathrsfs} |
14 |
– |
%\usepackage[ref]{overcite} |
16 |
|
\usepackage[square, comma, sort&compress]{natbib} |
17 |
|
\usepackage{url} |
18 |
|
\pagestyle{plain} \pagenumbering{arabic} \oddsidemargin 0.0cm |
20 |
|
9.0in \textwidth 6.5in \brokenpenalty=10000 |
21 |
|
|
22 |
|
% double space list of tables and figures |
23 |
< |
\AtBeginDelayedFloats{\renewcommand{\baselinestretch}{1.66}} |
23 |
> |
%\AtBeginDelayedFloats{\renewcomand{\baselinestretch}{1.66}} |
24 |
|
\setlength{\abovecaptionskip}{20 pt} |
25 |
|
\setlength{\belowcaptionskip}{30 pt} |
26 |
|
|
27 |
< |
%\renewcommand\citemid{\ } % no comma in optional reference note |
27 |
< |
\bibpunct{[}{]}{,}{n}{}{;} |
27 |
> |
\bibpunct{}{}{,}{s}{}{;} |
28 |
|
\bibliographystyle{achemso} |
29 |
|
|
30 |
|
\begin{document} |
49 |
|
%% |
50 |
|
|
51 |
|
%Title |
52 |
< |
\title{Investigation of the Pt and Au 557 Surface Reconstructions |
53 |
< |
under a CO Atmosphere} |
52 |
> |
\title{Molecular Dynamics simulations of the surface reconstructions |
53 |
> |
of Pt(557) and Au(557) under exposure to CO} |
54 |
> |
|
55 |
|
\author{Joseph R. Michalka, Patrick W. McIntyre and J. Daniel |
56 |
|
Gezelter\footnote{Corresponding author. \ Electronic mail: gezelter@nd.edu} \\ |
57 |
|
Department of Chemistry and Biochemistry,\\ |
58 |
|
University of Notre Dame\\ |
59 |
|
Notre Dame, Indiana 46556} |
60 |
+ |
|
61 |
|
%Date |
62 |
< |
\date{Dec 15, 2012} |
62 |
> |
\date{Dec 15, 2012} |
63 |
> |
|
64 |
|
%authors |
65 |
|
|
66 |
|
% make the title |
120 |
|
|
121 |
|
|
122 |
|
|
120 |
– |
|
121 |
– |
|
122 |
– |
|
123 |
|
\section{Simulation Methods} |
124 |
|
The challenge in modeling any solid/gas interface problem is the |
125 |
|
development of a sufficiently general yet computationally tractable |
248 |
|
%where did you actually get the functionals for citation? |
249 |
|
%scf calculations, so initial relaxation was of the four layers, but two layers weren't kept fixed, I don't think |
250 |
|
%same cutoff for slab and slab + CO ? seems low, although feibelmen had values around there... |
251 |
< |
The Au-C and Au-O interaction parameters were also fit to a Lennard-Jones |
252 |
< |
and Morse potential respectively, to reproduce Au-CO binding energies. |
253 |
< |
These energies were obtained from quantum calculations carried out using |
254 |
< |
the PBE GGA exchange-correlation functionals\cite{Perdew_GGA} for gold, carbon, and oxygen |
255 |
< |
constructed by Rappe, Rabe, Kaxiras, and Joannopoulos. \cite{RRKJ_PP}. |
256 |
< |
All calculations were run using the {\sc Quantum ESPRESSO} package. \cite{QE-2009} |
257 |
< |
First, a four layer slab of gold comprised of 32 atoms displaying a (111) surface was |
258 |
< |
converged using a 4X4X4 grid of Monkhorst-Pack \emph{k}-points.\cite{Monkhorst:1976} |
259 |
< |
The kinetic energy of the wavefunctions were truncated at 20 Ry while the |
260 |
< |
cutoff for the charge density and potential was set at 80 Ry. This relaxed |
261 |
< |
gold slab was then used in numerous single point calculations with CO at various heights |
262 |
< |
to create a potential energy surface for the Au-CO interaction. |
251 |
> |
The Au-C and Au-O cross-interactions were fit using Lennard-Jones and |
252 |
> |
Morse potentials, respectively, to reproduce Au-CO binding energies. |
253 |
|
|
254 |
+ |
The fits were refined against gas-surface calculations using DFT with |
255 |
+ |
a periodic supercell plane-wave basis approach, as implemented in the |
256 |
+ |
{\sc Quantum ESPRESSO} package.\cite{QE-2009} Electron cores are |
257 |
+ |
described with the projector augmented-wave (PAW) |
258 |
+ |
method,\cite{PhysRevB.50.17953,PhysRevB.59.1758} with plane waves |
259 |
+ |
included to an energy cutoff of 20 Ry. Electronic energies are |
260 |
+ |
computed with the PBE implementation of the generalized gradient |
261 |
+ |
approximation (GGA) for gold, carbon, and oxygen that was constructed |
262 |
+ |
by Rappe, Rabe, Kaxiras, and Joannopoulos.\cite{Perdew_GGA,RRKJ_PP} |
263 |
+ |
Ionic relaxations were performed until the energy difference between |
264 |
+ |
subsequent steps was less than 0.0001 eV. In testing the CO-Au |
265 |
+ |
interaction, Au(111) supercells were constructed of four layers of 4 |
266 |
+ |
Au x 2 Au surface planes and separated from vertical images by six |
267 |
+ |
layers of vacuum space. The surface atoms were all allowed to relax. |
268 |
+ |
Supercell calculations were performed nonspin-polarized, and energies |
269 |
+ |
were converged to within 0.03 meV per Au atom with a 4 x 4 x 4 |
270 |
+ |
Monkhorst-Pack\cite{Monkhorst:1976,PhysRevB.13.5188} {\bf k}-point |
271 |
+ |
sampling of the first Brillouin zone. The relaxed gold slab was then |
272 |
+ |
used in numerous single point calculations with CO at various heights |
273 |
+ |
(and angles relative to the surface) to allow fitting of the empirical |
274 |
+ |
force field. |
275 |
+ |
|
276 |
|
%Hint at future work |
277 |
|
The fit parameter sets employed in this work are shown in Table 2 and their |
278 |
|
reproduction of the binding energies are displayed in Table 3. Currently, |