ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/COonPt/firstTry.tex
(Generate patch)

Comparing trunk/COonPt/firstTry.tex (file contents):
Revision 3872 by jmichalk, Mon Mar 11 22:37:32 2013 UTC vs.
Revision 3875 by gezelter, Wed Mar 13 19:18:03 2013 UTC

# Line 1 | Line 1
1 < \documentclass[11pt]{article}
1 > \documentclass[journal = jpccck, manuscript = article]{achemso}
2 > \setkeys{acs}{usetitle = true}
3 > \usepackage{achemso}
4 > \usepackage{caption}
5 > \usepackage{float}
6 > \usepackage{geometry}
7 > \usepackage{natbib}
8 > \usepackage{setspace}
9 > \usepackage{xkeyval}
10 > %%%%%%%%%%%%%%%%%%%%%%%
11   \usepackage{amsmath}
12   \usepackage{amssymb}
13   \usepackage{times}
# Line 6 | Line 15
15   \usepackage{setspace}
16   \usepackage{endfloat}
17   \usepackage{caption}
18 < %\usepackage{tabularx}
18 > \usepackage{tabularx}
19 > \usepackage{longtable}
20   \usepackage{graphicx}
21   \usepackage{multirow}
22 < %\usepackage{booktabs}
23 < %\usepackage{bibentry}
24 < %\usepackage{mathrsfs}
25 < \usepackage[square, comma, sort&compress]{natbib}
22 > \usepackage{multicol}
23 >
24 > \usepackage[version=3]{mhchem}  % this is a great package for formatting chemical reactions
25 > % \usepackage[square, comma, sort&compress]{natbib}
26   \usepackage{url}
27   \pagestyle{plain} \pagenumbering{arabic} \oddsidemargin 0.0cm
28   \evensidemargin 0.0cm \topmargin -21pt \headsep 10pt \textheight
# Line 22 | Line 32
32   %\AtBeginDelayedFloats{\renewcomand{\baselinestretch}{1.66}}
33   \setlength{\abovecaptionskip}{20 pt}
34   \setlength{\belowcaptionskip}{30 pt}
35 + % \bibpunct{}{}{,}{s}{}{;}
36  
37 < \bibpunct{}{}{,}{s}{}{;}
38 < \bibliographystyle{achemso}
37 > %\citestyle{nature}
38 > % \bibliographystyle{achemso}
39  
40 < \begin{document}
40 > \title{Molecular Dynamics simulations of the surface reconstructions
41 >  of Pt(557) and Au(557) under exposure to CO}
42  
43 + \author{Joseph R. Michalka}
44 + \author{Patrick W. McIntyre}
45 + \author{J. Daniel Gezelter}
46 + \email{gezelter@nd.edu}
47 + \affiliation[University of Notre Dame]{251 Nieuwland Science Hall\\
48 +  Department of Chemistry and Biochemistry\\ University of Notre
49 +  Dame\\ Notre Dame, Indiana 46556}
50  
51 + \keywords{}
52 +
53 + \begin{document}
54 +
55 +
56   %%
57   %Introduction
58   %       Experimental observations
# Line 47 | Line 71
71   %Summary
72   %%
73  
50 %Title
51 \title{Molecular Dynamics simulations of the surface reconstructions
52  of Pt(557) and Au(557) under exposure to CO}
74  
54 \author{Joseph R. Michalka, Patrick W. McIntyre and J. Daniel
55 Gezelter\footnote{Corresponding author. \ Electronic mail: gezelter@nd.edu} \\
56 Department of Chemistry and Biochemistry,\\
57 University of Notre Dame\\
58 Notre Dame, Indiana 46556}
59
60 %Date
61 \date{Mar 5, 2013}
62
63 %authors
64
65 % make the title
66 \maketitle
67
68 \begin{doublespace}
69
75   \begin{abstract}
76   We examine surface reconstructions of Pt and Au(557) under
77   various CO coverages using molecular dynamics in order to
# Line 300 | Line 305 | zone.\cite{Monkhorst:1976,PhysRevB.13.5188} The relaxe
305   performed until the energy difference between subsequent steps
306   was less than $10^{-8}$ Ry.   Nonspin-polarized supercell calculations
307   were performed with a 4~x~4~x~4 Monkhorst-Pack {\bf k}-point sampling of the first Brillouin
308 < zone.\cite{Monkhorst:1976,PhysRevB.13.5188} The relaxed gold slab was
308 > zone.\cite{Monkhorst:1976} The relaxed gold slab was
309   then used in numerous single point calculations with CO at various
310   heights (and angles relative to the surface) to allow fitting of the
311   empirical force field.
# Line 346 | Line 351 | future work.
351    \multirow{2}{*}{\textbf{Pt-CO}} & \multirow{2}{*}{-1.9} & -1.4 \bibpunct{}{}{,}{n}{}{,}
352    (Ref. \protect\cite{Kelemen:1979}) \\
353   & &  -1.9 \bibpunct{}{}{,}{n}{}{,} (Ref. \protect\cite{Yeo}) \\ \hline
354 <  \textbf{Au-CO} & -0.39 & -0.40 \bibpunct{}{}{,}{n}{}{,}  (Ref. \protect\cite{TPD_Gold}) \\
354 >  \textbf{Au-CO} & -0.39 & -0.40 \bibpunct{}{}{,}{n}{}{,}  (Ref. \protect\cite{TPDGold}) \\
355    \hline
356   \end{tabular}
357   \label{tab:co_energies}
# Line 364 | Line 369 | The different bulk melting temperatures (1337~K for Au
369   1200~K were performed to confirm the relative
370   stability of the surfaces without a CO overlayer.  
371  
372 < The different bulk melting temperatures (1337~K for Au
373 < and 2045~K for Pt) suggest that any possible reconstruction should happen at
372 > The different bulk melting temperatures (1337~K for Au\cite{Au:melting}
373 > and 2045~K for Pt\cite{Pt:melting}) suggest that any possible reconstruction should happen at
374   different temperatures for the two metals.  The bare Au and Pt surfaces were
375   initially run in the canonical (NVT) ensemble at 800~K and 1000~K
376   respectively for 100 ps. The two surfaces were relatively stable at these
# Line 378 | Line 383 | phase.  These systems were allowed to reach thermal eq
383   Because of the difference in binding energies, nearly all of the CO was bound to the Pt surface, while
384   the Au surfaces often had a significant CO population in the gas
385   phase.  These systems were allowed to reach thermal equilibrium (over
386 < 5 ns) before being run in the microcanonical (NVE) ensemble for
387 < data collection. All of the systems examined had at least 40 ns in the
386 > 5~ns) before being run in the microcanonical (NVE) ensemble for
387 > data collection. All of the systems examined had at least 40~ns in the
388   data collection stage, although simulation times for some Pt of the
389   systems exceeded 200~ns.  Simulations were carried out using the open
390   source molecular dynamics package, OpenMD.\cite{Ewald,OOPSE}
# Line 393 | Line 398 | metal system. The surfaces to which no CO was exposed
398   \subsection{Structural remodeling}
399   The surfaces of both systems, upon dosage of CO, began
400   to undergo remodeling that was not observed in the bare
401 < metal system. The surfaces to which no CO was exposed
402 < did experience minor roughening of the step-edge, but the
401 > metal system. The surfaces which were not exposed to CO
402 > did experience minor roughening of the step-edge because
403 > of the elevated temperatures, but the
404   (557) lattice was well-maintained throughout the simulation
405   time. The Au systems were limited to greater amounts of
406   roughening, i.e. breakup of the step-edge, and some step
407   wandering. The lower coverage Pt systems experienced
408   similar restructuring but to a greater extent when
409   compared to the Au systems. The 50\% coverage
410 < Pt system formed double layers at numerous spots upon its surface.
410 > Pt system was unique among our simulations in that it
411 > formed numerous double layers through step coalescence,
412 > similar to results reported by Tao et al.\cite{Tao:2010}
413  
414  
415   \subsubsection{Step wandering}
416 < The 0\% coverage surfaces for both metals showed
417 < minimal movement at their respective run temperatures.
418 < As the coverage increased, the mobility of the surface
419 < also increased. Additionally, at the higher coverages
420 < on both metals, there was a large increase in the amount
421 < of observed step-wandering. Previous work by
422 < Williams\cite{Williams:1993} highlighted the entropic
423 < contribution to the repulsion felt between step-edges,
424 < and situations were that repulsion could be negated, or
425 < overcome, to allow for step coalescence or facet formation.
416 > The 0\% coverage surfaces for both metals showed minimal
417 > movement at their respective run temperatures. As the CO
418 > coverage increased however, the mobility of the surface,
419 > adatoms and step-edges alike, also increased. Additionally,
420 > at the higher coverages on both metals, there was more
421 > step-wandering. Except for the 50\% Pt system, the step-edges
422 > did not coalesce in any of the other simulations, instead preferring
423 > to keep nearly the same distance between steps as in the
424 > original (557) lattice. Previous work by Williams et al.\cite{Williams:1991, Williams:1994}
425 > highlights the repulsion that exists between step-edges even
426 > when no direct interactions are present in the system. This
427 > repulsion arises because the entropy of the step-edges is constrained,
428 > since step-edge crossing is not allowed. This entropic repulsion
429 > does not completely define the interactions between steps,
430 > which is why some surfaces will undergo step coalescence,
431 > where additional attractive interactions can overcome the
432 > repulsion\cite{Williams:1991} and others will not. The presence and concentration
433 > of adsorbates, as shown in this work, can affect these step interactions, potentially
434 > leading to a new surface structure as the thermodynamic minimum.
435  
436   \subsubsection{Double layers}
437   Tao et al. have shown experimentally that the Pt(557) surface
438 < undergoes two separate reconstructions upon CO
439 < adsorption.\cite{Tao:2010} The first involves a doubling of
440 < the step height and plateau length. Similar behavior has been
441 < seen to occur on numerous surfaces at varying conditions: Ni(977), Si(111).
442 < \cite{Williams:1994,Williams:1991,Pearl} Of the two systems
443 < we examined, the Pt system showed a greater propensity for
444 < reconstruction when compared to the Au system. The amount
445 < of reconstruction is correlated to the amount of CO
438 > undergoes two separate reconstructions upon CO adsorption.\cite{Tao:2010}
439 > The first involves a doubling of the step height and plateau length.
440 > Similar behavior has been seen to occur on numerous surfaces
441 > at varying conditions: Ni(977), Si(111).\cite{Williams:1994,Williams:1991,Pearl}
442 > Of the two systems we examined, the Pt system showed a greater
443 > propensity for reconstruction when compared to the Au system
444 > because of the larger surface mobility and extent of step wandering.
445 > The amount of reconstruction is correlated to the amount of CO
446   adsorbed upon the surface.  This appears to be related to the
447 < effect that adsorbate coverage has on edge breakup and on the surface
448 < diffusion of metal adatoms. While both systems displayed step-edge
449 < wandering, only the Pt surface underwent the doubling seen by
450 < Tao et al. within the time scales studied here.  
451 < Only the 50\% coverage Pt system exhibited
452 < a complete doubling in the time scales we
453 < were able to monitor. Over longer periods (150~ns) two more double layers formed on this interface.
454 < Although double layer formation did not occur in the other Pt systems, they show
455 < more lateral movement of the step-edges
456 < compared to their Au counterparts. The 50\% Pt system is highlighted
440 < in Figure \ref{fig:reconstruct} at various times along the simulation
441 < showing the evolution of a step-edge.
447 > effect that adsorbate coverage has on edge breakup and on the
448 > surface diffusion of metal adatoms. While both systems displayed
449 > step-edge wandering, only the 50\% Pt surface underwent the
450 > doubling seen by Tao et al.\cite{Tao:2010} within the time scales studied here.
451 > Over longer periods, (150~ns) two more double layers formed
452 > on this interface. Although double layer formation did not occur
453 > in the other Pt systems, they show more step-wandering and
454 > general roughening compared to their Au counterparts. The
455 > 50\% Pt system is highlighted in Figure \ref{fig:reconstruct} at
456 > various times along the simulation showing the evolution of a step-edge.
457  
458   The second reconstruction on the Pt(557) surface observed by
459   Tao involved the formation of triangular clusters that stretched
460   across the plateau between two step-edges. Neither system, within
461 < the 40~ns time scale, experienced this reconstruction.
461 > the 40~ns time scale or the extended simulation time of 150~ns for
462 > the 50\% Pt system, experienced this reconstruction.
463  
464   \subsection{Dynamics}
465   Previous atomistic simulations of stepped surfaces dealt largely
# Line 451 | Line 467 | of ignoring the dynamics of the system. Previous work
467   \cite{Williams:1991,Williams:1994}. Consequently, the most common
468   technique utilized to date has been Monte Carlo sampling. Monte Carlo gives an efficient
469   sampling of the equilibrium thermodynamic landscape at the expense
470 < of ignoring the dynamics of the system. Previous work by Pearl and
471 < Sibener\cite{Pearl}, using STM, has been able to show the coalescing
470 > of ignoring the dynamics of the system. Previous experimental work by Pearl and
471 > Sibener\cite{Pearl}, using STM, has been able to capture the coalescing
472   of steps on Ni(977). The time scale of the image acquisition,
473 < $\sim$70 s/image provides an upper bound for the time required for
473 > $\sim$70~s/image provides an upper bound for the time required for
474   the doubling to occur. In this section we give data on dynamic and
475   transport properties, e.g. diffusion, layer formation time, etc.
476  
# Line 462 | Line 478 | arising from the individual movements, primarily throu
478   \subsubsection{Transport of surface metal atoms}
479   %forcedSystems/stepSeparation
480   The movement or wandering of a step-edge is a cooperative effect
481 < arising from the individual movements, primarily through surface
466 < diffusion, of the atoms making up the steps. An ideal metal surface
481 > arising from the individual movements of the atoms making up the steps. An ideal metal surface
482   displaying a low index facet, (111) or (100), is unlikely to experience
483   much surface diffusion because of the large energetic barrier that must
484 < be overcome to lift an atom out of the surface. The presence of step-edges
485 < on higher-index surfaces provide a source for mobile metal atoms.
484 > be overcome to lift an atom out of the surface. The presence of step-edges and other surface features
485 > on higher-index facets provide a lower energy source for mobile metal atoms.
486   Breaking away from the step-edge on a clean surface still imposes an
487   energetic penalty around $\sim$~40 kcal/mol, but this is significantly easier than lifting
488   the same metal atom vertically out of the surface,  \textgreater~60 kcal/mol.
489   The penalty lowers significantly when CO is present in sufficient quantities
490   on the surface. For certain distributions of CO, the penalty can fall as low as
491   $\sim$~20 kcal/mol. Once an adatom exists on the surface, the barrier for
492 < diffusion is negligible ( \textless~4 kcal/mol) and these adatoms are
492 > diffusion is negligible ( \textless~4 kcal/mol for a Pt adatom). These adatoms are
493   able to explore the terrace before rejoining either the original step-edge or
494   becoming a part of a different edge. It is a more difficult process for an atom
495   to traverse to a separate terrace although the presence of CO can lower the
# Line 490 | Line 505 | movement of buried atoms. Diffusion on  a surface is s
505   between saved configurations of the system (typically 10-100 ps). An atom that was
506   truly mobile would typically travel much greater distances than this, but the 2~\AA~cutoff
507   was used to prevent swamping the diffusion data with the in-place vibrational
508 < movement of buried atoms. Diffusion on  a surface is strongly affected by
508 > movement of buried atoms. Diffusion on a surface is strongly affected by
509   local structures and in this work, the presence of single and double layer
510   step-edges causes the diffusion parallel to the step-edges to be different
511   from the diffusion perpendicular to these edges. Parallel and perpendicular
512   diffusion constants are shown in Figure \ref{fig:diff}.
513  
514 + The lack of a definite trend in the Au diffusion data is likely due
515 + to the weaker bonding between Au and CO. This leads to a lower
516 + coverage ({\it x}-axis) when compared to dosage amount, which
517 + then further limits the affects of the surface diffusion. The correlation
518 + between coverage and Pt diffusion rates conversely shows a
519 + definite trend marred by the highest coverage surface. Two
520 + explanations arise for this drop. First, upon a visual inspection of
521 + the system, after a double layer has been formed, it maintains its
522 + stability strongly and is no longer a good source for adatoms. By
523 + performing the same diffusion calculation but on a shorter run time
524 + (20~ns), only including data before the formation of the double layer,
525 + provides a $\mathbf{D}_{\perp}$ diffusion constant of $1.69~\pm~0.08$
526 + and a $\mathbf{D}_{\parallel}$ diffusion constant of $6.30~\pm~0.08$.
527 + This places the parallel diffusion constant more closely in line with the
528 + expected trend, while the perpendicular diffusion constant does not
529 + drop as far. A secondary explanation arising from our analysis of the
530 + mechanism of double layer formation show the affect that CO on the
531 + surface has with respect to overcoming surface diffusion of Pt. If the
532 + coverage is too sparse, the Pt engages in minimal interactions and
533 + thus minimal diffusion. As coverage increases, there are more favorable
534 + arrangements of CO on the surface allowing the formation of a path,
535 + a minimum energy trajectory, for the adatom to explore the surface.
536 + As the CO is constantly moving on the surface, this path is constantly
537 + changing. If the coverage becomes too great, the paths could
538 + potentially be clogged leading to a decrease in diffusion despite
539 + their being more adatoms and step-wandering.
540 +
541   \subsubsection{Dynamics of double layer formation}
542   The increased diffusion on Pt at the higher
543   CO coverages plays a primary role in double layer formation. However, this is not
544   a complete explanation -- the 33\%~Pt system
545   has higher diffusion constants but did not show
546 < any signs of edge doubling. On the
547 < 50\%~Pt system, three separate layers were formed over
548 < 150~ns of simulation time. Previous experimental
546 > any signs of edge doubling in the observed run time. On the
547 > 50\%~Pt system, one layer formed within the first 40~ns of simulation time, while two more were formed as the system was run for an additional
548 > 110~ns (150~ns total). Previous experimental
549   work gives insight into the upper bounds of the
550   time required for step coalescence.\cite{Williams:1991,Pearl}
551   In this system, as seen in Figure \ref{fig:reconstruct}, the first
552   appearance of a double layer, appears at 19~ns
553   into the simulation. Within 12~ns of this nucleation event, nearly half of the step has
554 < formed the double layer and by 86 ns, the complete layer
554 > formed the double layer and by 86~ns, the complete layer
555   has been flattened out. The double layer could be considered
556   ``complete" by 37~ns but remains a bit rough. From the
557   appearance of the first nucleation event to the first observed double layer, the process took $\sim$20~ns. Another
# Line 526 | Line 568 | take longer at lower temperatures.
568   \begin{figure}[H]
569   \includegraphics[width=\linewidth]{ProgressionOfDoubleLayerFormation_yellowCircle.png}
570   \caption{The Pt(557) / 50\% CO system at a sequence of times after
571 <  initial exposure to the CO: (a) 258 ps, (b) 19 ns, (c) 31.2 ns, and
572 <  (d) 86.1 ns. Disruption of the (557) step-edges occurs quickly.  The
571 >  initial exposure to the CO: (a) 258~ps, (b) 19~ns, (c) 31.2~ns, and
572 >  (d) 86.1~ns. Disruption of the (557) step-edges occurs quickly.  The
573    doubling of the layers appears only after two adjacent step-edges
574    touch.  The circled spot in (b) nucleated the growth of the double
575    step observed in the later configurations.}
# Line 574 | Line 616 | From simulations which exhibit a double layer, the tim
616   This combination of growth and decay of the step-edges is in a state of
617   dynamic equilibrium. However, once two previously separated edges
618   meet as shown in Figure 1.B, this nucleates the rest of the edge to meet up, forming a double layer.
619 < From simulations which exhibit a double layer, the time delay from the initial appearance of a nucleation point to a fully formed double layer is $\sim$35 ns.
619 > From simulations which exhibit a double layer, the time delay from the initial appearance of a nucleation point to a fully formed double layer is $\sim$35~ns.
620  
621   A number of possible mechanisms exist to explain the role of adsorbed
622   CO in restructuring the Pt surface. Quadrupolar repulsion between adjacent
# Line 605 | Line 647 | are displayed in Table \ref{tab:energies}. These value
647   of Pt atoms was then examined to determine possible barriers. Because
648   the movement was forced along a pre-defined reaction coordinate that may differ
649   from the true minimum of this path, only the beginning and ending energies
650 < are displayed in Table \ref{tab:energies}. These values suggest that the presence of CO at suitable
650 > are displayed in Table \ref{tab:rxcoord} with the corresponding beginning and ending reaction coordinates in Figure \ref{fig:lambdaTable}. These values suggest that the presence of CO at suitable
651   locations can lead to lowered barriers for Pt breaking apart from the step-edge.
652   Additionally, as highlighted in Figure \ref{fig:lambda}, the presence of CO makes the
653   burrowing and lifting of adatoms favorable, whereas without CO, the process is neutral
# Line 613 | Line 655 | in terms of energetics.
655  
656   %lambda progression of Pt -> shoving its way into the step
657   \begin{figure}[H]
658 < \includegraphics[width=\linewidth]{lambdaProgression_atopCO.png}
658 > \includegraphics[width=\linewidth]{lambdaProgression_atopCO_withLambda.png}
659   \caption{A model system of the Pt(557) surface was used as the framework
660   for exploring energy barriers along a reaction coordinate. Various numbers,
661   placements, and rotations of CO were examined as they affect Pt movement.
# Line 624 | Line 666 | in terms of energetics.
666   \label{fig:lambda}
667   \end{figure}
668  
669 + \begin{figure}[H]
670 + \includegraphics[totalheight=0.9\textheight]{lambdaTable.png}
671 + \caption{}
672 + \label{fig:lambdaTable}
673 + \end{figure}
674  
675  
676 +
677 + \begin{table}[H]
678 + \caption{}
679 + \centering
680 + \begin{tabular}{| c || c | c | c | c |}
681 + \hline
682 + \textbf{System} & 0.5~\AA & 2~\AA & 4~\AA & 6~\AA \\
683 + \hline
684 + A & 6.38 & 38.34 & 44.65 & 47.60 \\
685 + B & -20.72 & 0.67 & 17.33 & 24.28 \\
686 + C & 4.92 & 27.02 & 41.05 & 47.43 \\
687 + D & -16.97 & 21.21 & 35.87 & 40.93 \\
688 + E & 5.92 & 30.96 & 43.69 & 49.23 \\
689 + F & 8.53 & 46.23 & 53.98 & 65.55 \\
690 + \hline
691 + \end{tabular}
692 + \label{tab:rxcoord}
693 + \end{table}
694 +
695 +
696   \subsection{Diffusion}
697   The diffusion parallel to the step-edge tends to be
698   much larger than that perpendicular to the step-edge. The dynamic
# Line 644 | Line 711 | more likely a growth point is to be formed.
711   %breaking of the double layer upon removal of CO
712   \begin{figure}[H]
713   \includegraphics[width=\linewidth]{doubleLayerBreaking_greenBlue_whiteLetters.png}
714 < %:
648 < \caption{(A)  0 ps, (B) 100 ps, (C) 1 ns, after the removal of CO. The presence of the CO
714 > \caption{(A)  0~ps, (B) 100~ps, (C) 1~ns, after the removal of CO. The presence of the CO
715   helped maintain the stability of the double layer and upon removal the two layers break
716   and begin separating. The separation is not a simple pulling apart however, rather
717   there is a mixing of the lower and upper atoms at the edge.}
# Line 700 | Line 766 | In this work we have shown the reconstruction of the P
766   % \end{tabular}
767   % \end{table}
768  
769 < \section{Acknowledgments}
769 > \begin{acknowledgement}
770   Support for this project was provided by the National Science
771   Foundation under grant CHE-0848243 and by the Center for Sustainable
772   Energy at Notre Dame (cSEND). Computational time was provided by the
773   Center for Research Computing (CRC) at the University of Notre Dame.
774 <
774 > \end{acknowledgement}
775   \newpage
776   \bibliography{firstTryBibliography}
777 < \end{doublespace}
777 > %\end{doublespace}
778 >
779 > \begin{tocentry}
780 > %\includegraphics[height=3.5cm]{timelapse}
781 > \end{tocentry}
782 >
783   \end{document}

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines