A Critical Appraisal of Data Augmentation Methods for Imaging-Based Medical Diagnosis Applications

Abstract: Current data augmentation techniques and transformations are well suited for improving the size and quality of natural image datasets but are not yet optimized for medical imaging. We hypothesize that sub-optimal data augmentations can easily distort or occlude medical images, leading to false positives or negatives during patient diagnosis, prediction, or therapy/surgery evaluation. In our experimental results, we found that utilizing commonly used intensity-based data augmentation distorts the MRI scans and leads to texture information loss, thus negatively affecting the overall performance of classification. Additionally, we observed that commonly used data augmentation methods cannot be used with a plug-and-play approach in medical imaging, and requires manual tuning and adjustment.

\section{Introduction}
This document is a model and instructions for \LaTeX.
Please observe the conference page limits.

\section{Ease of Use}

\subsection{Maintaining the Integrity of the Specifications}

The IEEEtran class file is used to format your paper and style the text. All margins,
column widths, line spaces, and text fonts are prescribed; please do not
alter them. You may note peculiarities. For example, the head margin
measures proportionately more than is customary. This measurement
and others are deliberate, using specifications that anticipate your paper
as one part of the entire proceedings, and not as an independent document.
Please do not revise any of the current designations.

\section{Prepare Your Paper Before Styling}
Before you begin to format your paper, first write and save the content as a
separate text file. Complete all content and organizational editing before
formatting. Please note sections \ref{AA}--\ref{SCM} below for more information on
proofreading, spelling and grammar.

Keep your text and graphic files separate until after the text has been
formatted and styled. Do not number text heads---{\LaTeX} will do that
for you.

\subsection{Abbreviations and Acronyms}\label{AA}
Define abbreviations and acronyms the first time they are used in the text,
even after they have been defined in the abstract. Abbreviations such as
IEEE, SI, MKS, CGS, ac, dc, and rms do not have to be defined. Do not use
abbreviations in the title or heads unless they are unavoidable.

\subsection{Units}
\begin{itemize}
\item Use either SI (MKS) or CGS as primary units. (SI units are encouraged.) English units may be used as secondary units (in parentheses). An exception would be the use of English units as identifiers in trade, such as ``3.5-inch disk drive''.
\item Avoid combining SI and CGS units, such as current in amperes and magnetic field in oersteds. This often leads to confusion because equations do not balance dimensionally. If you must use mixed units, clearly state the units for each quantity that you use in an equation.
\item Do not mix complete spellings and abbreviations of units: ``Wb/m\textsuperscript{2}'' or ``webers per square meter'', not ``webers/m\textsuperscript{2}''. Spell out units when they appear in text: ``. . . a few henries'', not ``. . . a few H''.
\item Use a zero before decimal points: ``0.25'', not ``.25''. Use ``cm\textsuperscript{3}'', not ``cc''.)
\end{itemize}

\subsection{Equations}
Number equations consecutively. To make your
equations more compact, you may use the solidus (~/~), the exp function, or
appropriate exponents. Italicize Roman symbols for quantities and variables,
but not Greek symbols. Use a long dash rather than a hyphen for a minus
sign. Punctuate equations with commas or periods when they are part of a
sentence, as in:
\begin{equation}
a+b=\gamma\label{eq}
\end{equation}

Be sure that the
symbols in your equation have been defined before or immediately following
the equation. Use ``\eqref{eq}'', not ``Eq.~\eqref{eq}'' or ``equation \eqref{eq}'', except at
the beginning of a sentence: ``Equation \eqref{eq} is . . .''

\subsection{\LaTeX-Specific Advice}

Please use ``soft'' (e.g., \verb|\eqref{Eq}|) cross references instead
of ``hard'' references (e.g., \verb|(1)|). That will make it possible
to combine sections, add equations, or change the order of figures or
citations without having to go through the file line by line.

Please don't use the \verb|{eqnarray}| equation environment. Use
\verb|{align}| or \verb|{IEEEeqnarray}| instead. The \verb|{eqnarray}|
environment leaves unsightly spaces around relation symbols.

Please note that the \verb|{subequations}| environment in {\LaTeX}
will increment the main equation counter even when there are no
equation numbers displayed. If you forget that, you might write an
article in which the equation numbers skip from (17) to (20), causing
the copy editors to wonder if you've discovered a new method of
counting.

{\BibTeX} does not work by magic. It doesn't get the bibliographic
data from thin air but from .bib files. If you use {\BibTeX} to produce a
bibliography you must send the .bib files.

{\LaTeX} can't read your mind. If you assign the same label to a
subsubsection and a table, you might find that Table I has been cross
referenced as Table IV-B3.

{\LaTeX} does not have precognitive abilities. If you put a
\verb|\label| command before the command that updates the counter it's
supposed to be using, the label will pick up the last counter to be
cross referenced instead. In particular, a \verb|\label| command
should not go before the caption of a figure or a table.

Do not use \verb|\nonumber| inside the \verb|{array}| environment. It
will not stop equation numbers inside \verb|{array}| (there won't be
any anyway) and it might stop a wanted equation number in the
surrounding equation.

\subsection{Some Common Mistakes}\label{SCM}
\begin{itemize}
\item The word ``data'' is plural, not singular.
\item The subscript for the permeability of vacuum $\mu_{0}$, and other common scientific constants, is zero with subscript formatting, not a lowercase letter ``o''.
\item In American English, commas, semicolons, periods, question and exclamation marks are located within quotation marks only when a complete thought or name is cited, such as a title or full quotation. When quotation marks are used, instead of a bold or italic typeface, to highlight a word or phrase, punctuation should appear outside of the quotation marks. A parenthetical phrase or statement at the end of a sentence is punctuated outside of the closing parenthesis (like this). (A parenthetical sentence is punctuated within the parentheses.)
\item A graph within a graph is an ``inset'', not an ``insert''. The word alternatively is preferred to the word ``alternately'' (unless you really mean something that alternates).
\item Do not use the word ``essentially'' to mean ``approximately'' or ``effectively''.
\item In your paper title, if the words ``that uses'' can accurately replace the word ``using'', capitalize the ``u''; if not, keep using lower-cased.
\item Be aware of the different meanings of the homophones ``affect'' and ``effect'', ``complement'' and ``compliment'', ``discreet'' and ``discrete'', ``principal'' and ``principle''.
\item Do not confuse ``imply'' and ``infer''.
\item The prefix ``non'' is not a word; it should be joined to the word it modifies, usually without a hyphen.
\item There is no period after the ``et'' in the Latin abbreviation ``et al.''.
\item The abbreviation ``i.e.'' means ``that is'', and the abbreviation ``e.g.'' means ``for example''.
\end{itemize}
An excellent style manual for science writers is \cite{b7}.

\subsection{Authors and Affiliations}
\textbf{The class file is designed for, but not limited to, six authors.} A
minimum of one author is required for all conference articles. Author names
should be listed starting from left to right and then moving down to the
next line. This is the author sequence that will be used in future citations
and by indexing services. Names should not be listed in columns nor group by
affiliation. Please keep your affiliations as succinct as possible (for
example, do not differentiate among departments of the same organization).

\subsection{Identify the Headings}
Headings, or heads, are organizational devices that guide the reader through
your paper. There are two types: component heads and text heads.

Component heads identify the different components of your paper and are not
topically subordinate to each other. Examples include Acknowledgments and
References and, for these, the correct style to use is ``Heading 5''. Use
``figure caption'' for your Figure captions, and ``table head'' for your
table title. Run-in heads, such as ``Abstract'', will require you to apply a
style (in this case, italic) in addition to the style provided by the drop
down menu to differentiate the head from the text.

Text heads organize the topics on a relational, hierarchical basis. For
example, the paper title is the primary text head because all subsequent
material relates and elaborates on this one topic. If there are two or more
sub-topics, the next level head (uppercase Roman numerals) should be used
and, conversely, if there are not at least two sub-topics, then no subheads
should be introduced.

\subsection{Figures and Tables}
\paragraph{Positioning Figures and Tables} Place figures and tables at the top and
bottom of columns. Avoid placing them in the middle of columns. Large
figures and tables may span across both columns. Figure captions should be
below the figures; table heads should appear above the tables. Insert
figures and tables after they are cited in the text. Use the abbreviation
``Fig.~\ref{fig}'', even at the beginning of a sentence.

\begin{table}[htbp]
\caption{Table Type Styles}
\begin{center}
\begin{tabular}{|c|c|c|c|}
\hline
\textbf{Table}&\multicolumn{3}{|c|}{\textbf{Table Column Head}} \\
\cline{2-4}
\textbf{Head} & \textbf{\textit{Table column subhead}}& \textbf{\textit{Subhead}}& \textbf{\textit{Subhead}} \\
\hline
copy& More table copy$^{\mathrm{a}}$& &  \\
\hline
\multicolumn{4}{l}{$^{\mathrm{a}}$Sample of a Table footnote.}
\end{tabular}
\label{tab1}
\end{center}
\end{table}

\begin{figure}[htbp]
\centerline{\includegraphics{fig1.png}}
\caption{Example of a figure caption.}
\label{fig}
\end{figure}

Figure Labels: Use 8 point Times New Roman for Figure labels. Use words
rather than symbols or abbreviations when writing Figure axis labels to
avoid confusing the reader. As an example, write the quantity
``Magnetization'', or ``Magnetization, M'', not just ``M''. If including
units in the label, present them within parentheses. Do not label axes only
with units. In the example, write ``Magnetization (A/m)'' or ``Magnetization
\{A[m(1)]\}'', not just ``A/m''. Do not label axes with a ratio of
quantities and units. For example, write ``Temperature (K)'', not
``Temperature/K''.

\section*{Acknowledgment}

The preferred spelling of the word ``acknowledgment'' in America is without
an ``e'' after the ``g''. Avoid the stilted expression ``one of us (R. B.
G.) thanks $\ldots$''. Instead, try ``R. B. G. thanks$\ldots$''. Put sponsor
acknowledgments in the unnumbered footnote on the first page.

\section*{References}

Please number citations consecutively within brackets {\cite{b1}. The
sentence punctuation follows the bracket \cite{b2}. Refer simply to the reference
number, as in \cite{b3}---do not use ``Ref. \cite{b3}'' or ``reference \cite{b3}'' except at
the beginning of a sentence: ``Reference \cite{b3} was the first $\ldots$''

Number footnotes separately in superscripts. Place the actual footnote at
the bottom of the column in which it was cited. Do not put footnotes in the
abstract or reference list. Use letters for table footnotes.

Unless there are six authors or more give all authors' names; do not use
``et al.''. Papers that have not been published, even if they have been
submitted for publication, should be cited as ``unpublished'' \cite{b4}. Papers
that have been accepted for publication should be cited as ``in press'' \cite{b5}.
Capitalize only the first word in a paper title, except for proper nouns and
element symbols.

For papers published in translation journals, please give the English
citation first, followed by the original foreign-language citation \cite{b6}.

\section{Introduction}

In this study, we investigate how commonly used data augmentation methods affect medical imaging based diagnosis tasks, the following methods are used: RandAugment~\cite{randaug}, AutoAugment~\cite{aa}, Fast AutoAugment~\cite{faa}, Trivial Augment~\cite{trivialaug} and AugMix~\cite{augmix}.

\begin{figure}[!h]
\centerline{\includegraphics[width=0.8\linewidth]{figs/arch1.jpg}}
\vspace{-3mm}
\caption{Overall architecture for ResNeST~\cite{zhang2020resnest} training with data augmentation.
\label{fig:archicture}
\vspace{-6mm}
\end{figure}

 \section{Method}
 Figure~\ref{fig:archicture} shows the block diagram of the proposed method.
In this study, we tested commonly used data augmentation methods RandAugment \cite{randaug}, AutoAugment~\cite{aa}, Fast AutoAugment~\cite{faa}, Trivial Augment~\cite{trivialaug} and AugMix~\cite{augmix}, and their impact on the MRI based IPMN classification problem. The classification problem was modelled via a deep learning architecture, called ResNeST~\cite{zhang2020resnest}, where we classified the images into normal, low grade and high grade categories.

\begin{table}
\caption{MRI based IPMN classification results with respect to data augmentation methods.}
\vspace{-7mm}
\begin{center}
\begin{tabular}{l|c}
\hline
\textbf{Method} & \textbf{Accuracy} \\
\hline
Baseline & $61.70\% \pm 9.42$ \\
\hline
RandAug & $55.30\% \pm 9.29$ \\
\hline
RandAug - geometric & $61.67\% \pm 8.37$ \\
\hline
Auto Augment &  $51.03\% \pm 10.31$ \\
\hline
Fast Auto Augment &  $55.34\% \pm 12.78$ \\
\hline
Trivial Augment &  $51.75\% \pm 9.30$ \\
\hline
AugMix & $54.56\% \pm 8.19$ \\
\hline
\end{tabular}
\label{table:results}
\end{center}
\vspace{-7mm}
\end{table}

\subsection{Results}

Table~\ref{table:results} exhibits the performance of each augmentation method and the baseline. It can be observed that the baseline achieves an accuracy of $61.70\% \pm 9.42$. However, other augmentation methods such as RandAug, Auto Augment, Fast Auto Augment, Trivial Augment and AugMix are causing significant performance drop (refer Table~\ref{table:results}). Table shows that when the intensity based augmentation is removed performance stays close to the baseline in the ``RandAug - geometric" experiments where it achieves $61.67\% \pm 8.37$.

\section{Conclusion}
In this study, we raised critical appraisals for the role of data augmentation for medical imaging tasks. We analyzed five commonly used data augmentation approaches and their effect on the performance of the MRI based IPMN classification problem. Our study in the controlled experiments showed that the commonly used data augmentation methods are designed  specifically for natural images and they can have adverse effects in medical diagnosis tasks if used without modification.

 \subsubsection*{Acknowledgement}
 This project is supported by the NIH funding: R01-CA246704 and R01-CA240639.
\vspace{-2mm}
\bibliographystyle{IEEEtran}

Image and Video Processing (eess.IV); Computer Vision and Pattern Recognition (cs.CV)


Doping of graphene adsorbed on the a-SiO$_2$ surface

Abstract: We have performed an {\it ab initio} theoretical investigation of graphene sheet adsorbed on amorphous SiO$_2$ surface (G/a-SiO$_2$). We find that graphene adsorbs on the a-SiO$_2$ surface through van der Waals interactions. The inhomogeneous topology of the a-SiO$_2$ clean surface promotes a total charge density displacement on the adsorbed graphene sheet, giving rise to electron-rich as well as hole-rich regions on the graphene. Such anisotropic distribution of the charge density may contribute to the reduction of the electronic mobility in G/a-SiO$_2$ systems. Furthermore, the adsorbed graphene sheet exhibits a net total charge density gain. In this case, the graphene sheet becomes n-type doped, however, with no formation of chemical bonds at the graphene--SiO$_2$ interface. The electronic charge transfer from a-SiO$_2$ to the graphene sheet occurs upon the formation of a partially occupied level lying above the Dirac point. We find that such partially occupied level comes from the three-fold coordinated oxygen atoms in the a-SiO$_2$ substrate.

\section{Introduction}

Considerable  progress has  been  made addressing  the electronic  and
structural  properties  of   graphene,  however,  the  development  of
nanodevices based on  graphene is still in its  infancy. Graphene is a
semiconductor material  with zero energy bandgap. The  valence and the
conduction  bands  are  connected  through  a  linear  energy-momentum
relation  at  the K  points  lying at  the  corners  of the  hexagonal
Brillouin zone (Dirac  points)~\cite{castronetoRMP2009}. The very high
electronic  mobility, up  to  200,000 cm$^2$/(V.s),  makes graphene  a
quite  suitable   material  to  the  development   of  new  electronic
(nano)devices.  Indeed, there are  experimental realizations  of field
effect       transistors      (FETs)      based       on      graphene
sheets~\cite{lemmeIEEE2007,liScience2008,wangPRL2008,oostingaNatMat2008,romeroACSNano2008}.
In  those  systems, SiO$_2$  has  been  used  as the  gate  dielectric
material. Lemme {\it et al.}~\cite{lemmeIEEE2007} verified a reduction
of  the  electron mobility  upon  the  graphene  interaction with  the
SiO$_2$  gate,  while  Romero  {\it  et  al.}~\cite{romeroACSNano2008}
obtained a n-type doping for  graphene in contact with a SiO$_2$ gate.
In this case, the authors attributed the presence of surface states to
the electronic charge transfer from  SiO$_2$ to the graphene sheet. In
contrast, there  are experimental  results indicating a  p-type doping
for                 graphene                 on                SiO$_2$
gate~\cite{berciaudNanoLett2009,songNanotech2010}.  Very  recently  n-
and p-type  doping of graphene,  lying on SiO$_2$ substrate,  has been
tuned    through     a    suitable    graphene--SiO$_2$    ``interface
engineering''~\cite{wangACSNano2011},  where  foreign  molecules  have
been  placed  between the  graphene  sheet  and  the SiO$_2$  surface.
Indeed,  there  are experimental  and  theoretical studies  addressing
``interface  engineering''  applied   to  modify  the  electronic  and
structural    properties    of    carbon   nanotubes    adsorbed    on
SiO$_2$~\cite{soaresNanoLett2010}   as   well    as   on   the   other
semiconductor   surfaces   like   Si(001)~\cite{OrellanaPRL2003}   and
InAs(111)~\cite{kimPRL2004}.

The  most of  the  experimental  results indicate  that  there are  no
covalent bonds  between the graphene  sheet and the  amorphous SiO$_2$
substrate   (a-SiO$_2$).  It  has   been  suggested   a  physisorption
process~\cite{stolyarovaPNAS2007,ishigamiNanoLett2007,sinitskiiACSNano2010}
where  van  der  Waals  (vdW)  interactions play  an  important  role.
Graphene sheet feels the corrugations of the a-SiO$_2$ surface, giving
rise to  electron-rich as  well as hole-rich  regions on  graphene, so
called ``electron-hole  puddles''~\cite{martinNatPhys2008}. Based upon
the aforementioned  findings, the understanding of  the electronic and
the  energetic properties,  within  an atomic  scale,  of graphene  on
a-SiO$_2$ surface is  a important issue to the  examined by using {\it
  ab initio} total energy calculations.

\begin{figure}[h]
\includegraphics[width= 7cm]{structure4bb.eps}
\caption{Atomic structure of graphene adsorbed a-SiO$_2$ surface,
  G/a-SiO$_2$ (a) side-view, (b) top-view, details of the equilibrium
  geometry of O$_{\rm three-fold}$ (c).  In (a) the bottommost O
  dangling bond was saturated by a H atom and the dashed region
  indicate a periodic unit cell. In (c) the atomic distances are in
  \AA.  Filled (Empty) circles represent the O (Si) atoms. (d) Total
  charge density fluctuation ($\Delta\rho$) within a region of
  $\pm$0.5~\AA\ with respect to the (equilibrium) vertical position of
  the n-type doped graphene sheet. Within the parentheses we present
  $\Delta\rho$ for the neutral G/a-SiO$_2$ graphene sheet. $\Delta\rho$
  is in unit of $10^{14} e/{\rm cm}^2$.}
\label{structure}
\end{figure}

In this  letter we  present {\it ab  initio} calculations  of graphene
adsorbed on the a-SiO$_2$  surface, G/a-SiO$_2$. The calculations were
performed  within the  density functional  theory (DFT)  including the
long-range vdW interactions.  Our total energy results reveal that vdW
interactions   rule   the  graphene   adsorption   on  the   a-SiO$_2$
surface.  The  inhomogeneous  topology  of  a-SiO$_2$  gives  rise  to
electron-rich and hole-rich regions on the adsorbed graphene sheet, in
accordance with the  experimentally observed ``electron-hole puddles''
for  graphene lying  on  a-SiO$_2$.  However, there  is  a net  charge
density  gain on the  graphene sheet.  The electronic  charge transfer
from a-SiO$_2$ to graphene occurs  due to the formation of a partially
occupied  state above  the Dirac  point. We  find that  such partially
occupied  state   is  attributed   to  the  formation   of  three-fold
coordinated O atoms within a-SiO$_2$.

The  amorphous  structure  was   generated  through  {\it  ab  initio}
molecular  dynamics (MD)  simulations  based on  the  DFT approach  as
implemented                 in                 the                VASP
code~\cite{kressePRB1993,kresseJPhys1994,kressePRB1996}.  Initially we
obtained the  a-SiO$_2$ bulk phase. We  have used a  supercell with 96
atoms and  we made a suitable  choice of the lattice  constants of the
supercell,  12.47  and  12.96~\AA\  along  the  x  and  y  directions,
respectively,  in  order  to  match with  the  calculated  equilibrium
lattice     constant     of     an     isolated     graphene     sheet
(Fig.~\ref{structure}). While the supercell volume was minimized along
the z direction, thus, avoiding any artificial strain in the a-SiO$_2$
bulk.  We  next include a  vacuum region of $\sim$12~\AA\  parallel to
the z  direction, within the  slab method, and  then we repeat  the MD
simulations  in   order  to  get  an  a-SiO$_2$   clean  surface.   In
Ref.~\cite{scopelPRB2008}  we   present  details  on   the  generation
procedure of amorphous SiO$_2$ bulk structure. Once we obtained a well
described a-SiO$_2$  surface, we start to  investigate the equilibrium
geometry,  energetic stability,  and  the electronic  properties of  a
graphene sheet adsorbed onto  a-SiO$_2$ surface. Here the calculations
were performed by using the  DFT approach as implemented in the SIESTA
code~\cite{siesta}.  The  generalized  gradient approximation  due  to
Perdew, Burke and Ernzerhof (PBE)~\cite{PBE} was used, and the van der
Waals  interaction  was  described  within a  semiempirical  approach,
following    the   Grimme    formula~\cite{grimmeJCompChem2006}.   The
electron-ion  interactions were  calculated  by using  norm-conserving
pseudopotentials~\cite{pseudoTM}.   All the  atomic  positions of  the
G/a-SiO$_2$  system  were relaxed  by  using  the conjugated  gradient
scheme, within a force convergence criterion of 20~meV/\AA.

Figure~\ref{structure}(a)  presents  the   structural  model  for  the
physisorbed  G/a-SiO$_2$   system,  where  we   have  considered  nine
G/a-SiO$_2$   configurations  on   two  different   a-SiO$_2$  surface
structures.     At    the    equilibrium    geometry   we    find    a
graphene~--~a-SiO$_2$    equilibrium   vertical    distance   ($d_{\rm
  G-SiO_2}$) of 3.64~(3.30)~\AA\  for an a-SiO$_2$ surface corrugation
of  1.7   (1.5)~\AA~\cite{corrug}.   Although  the   relatively  short
distance between the graphene sheet and the topmost O atoms, we do not
find the  formation of C--O  chemical bonds. The strength  of graphene
adsorption  on  a-SiO$_2$  can  be  measured by  comparing  the  total
energies  of the  isolated  systems, {\it  viz.}:  graphene sheet  and
a-SiO$_2$  surface, with  the total  energy of  the  graphene adsorbed
system, G/a-SiO$_2$.  We find that  the formation of G/a-SiO$_2$ is an
exothermic   process  with   an  adsorption   energy   ($E^{ads}$)  of
6.3$\pm$0.4~{\rm  meV/\AA$^2$} Our  adsorption energy  and equilibrium
geometry  results   are  in  good  agreement   with  the  experimental
estimative of Ishigami {\it  et al.}, $E^{ads}=6~{\rm meV/\AA^2}$, and
$d_{\rm   G-SiO_2}$  =   4.2~\AA~\cite{ishigamiNanoLett2007}.   Recent
theoretical  calculations,   within  the  DFT-LDA   approach  with  no
inclusion of vdW interations, obtained $E^{ads}=1~{\rm meV/\AA^2}$ and
$d_{\rm      G-SiO_2}       =3.6~$~\AA\      for      graphene      on
a-SiO$_2$~\cite{romeroACSNano2008}.   In contrast,  other experimental
measurements obtained  $d_{\rm G-SiO_2}$ of  5.8~\AA\ for G/a-SiO$_2$,
which may indicate the presence  of foreign elements at the G--SiO$_2$
interface~\cite{songNanotech2010}.   We  next  calculate $E^{ads}$  by
turning off  the vdW contribution  from our total energy  results.  In
this case, the formation  of G/a-SiO$_2$ becomes no longer exothermic,
$E^{ads}=-0.4$~meV/\AA$^2$. Thus, showing  that graphene sheet adsorbs
on a-SiO$_2$ is mediated by vdW interactions.

Recent experimental findings indicate a corrugation of $\sim$5~\AA\ on
the     a-SiO$_2$    surface,     for    a     lateral     scale    of
$\sim$100~\AA~\cite{stolyarovaPNAS2007}.          Meanwhile         in
Ref.~\cite{ishigamiNanoLett2007}   the   authors   suggest  that   the
corrugation of the adsorbed graphene sheet should correspond to around
60\% of  the one  of the a-SiO$_2$  surface. They measured  a graphene
corrugation  of   1.9~\AA\  and  estimate  a   deformation  energy  of
1~meV/\AA$^2$.  We  obtained a very small corrugation  (up to 0.3~\AA)
and  deformation energy  (0.7~meV/\AA$^2$) for  the  adsorbed graphene
sheet, whereas our simulated  a-SiO$_2$ surface exhibits a corrugation
of $\sim$1.6~\AA. Here we would be able to describe the corrugation of
G/a-SiO$_2$ in  a suitable way by  increasing the size  of the surface
unit cell,  and performing an additional total  energy minimization of
the  surface  area.  In  this  case,  we may  find  higher  values  of
$E^{ads}$,  since the  graphene--SiO$_2$ contact  area  will increase,
however, such  increase of  $E^{ads}$ will be  limited by  the induced
strain energy upon the deformation of the graphene sheet.

It is  worth to  note that graphene  adsorbed on a-SiO$_2$  exhibits a
quite   different    picture   when   compared    with   graphene   on
$\alpha$-quartz.   In  the  latter,  due  to the  presence  of  oxygen
dangling bonds,  there is the  formation of C--O chemical  bonds (with
$E^{ads}=300~{\rm      meV/\AA^2}$)      on      the      O-terminated
surface~\cite{kangPRB2008}.  In  addition,  it  has been  suggested  a
p-type doping  of the adsorbed graphene sheet.   Whereas in a-SiO$_2$,
which also exhibits an O-rich surface [Fig.~\ref{structure}(a)], we do
not find the formation C--O chemical bonds. In our simulated annealing
process  all   the  surface   O  atoms  become   twofold  coordinated,
suppressing  the (energetically unfavorable)  O dangling  bonds.  Very
recent  DFT-LDA calculations  indicate that  the  graphene interaction
with  the (0001)  surface of  $\alpha$-quartz  is weak,  with no  C--O
chemical bonds~\cite{nguyenPRL2011}.

Figure~\ref{structure}(d) depicts the  net charge density displacement
($\Delta\rho$) on  the adsorbed  graphene sheet.  $\Delta\rho$  can be
written as,
$$
\Delta\rho = \rho[{\rm G_{SiO_2}}] - \rho[{\rm G}].
$$ Where $\rho[{\rm G_{SiO_2}}]$ and $\rho[{\rm G}]$ represent the planar
averaged  total  charge  density  of the  adsorbed  (G/a-SiO$_2$)  and
isolated  graphene  sheets,  respectively.   The  planar  average  was
performed  by considering  a region  of 0.5~\AA\  above and  below the
graphene  sheet.   The  last  term, $\rho[{\rm  G}]$,  was  calculated
keeping the same equilibrium geometry  as that of G/a-SiO$_2$. We find
graphene regions  with $\Delta\rho > 0$  as well as  $\Delta\rho < 0$,
electron-     and    hole-rich    regions,     respectively.     Thus,
Fig.~\ref{structure}(d) indicates that there is a total charge density
fluctuation on the graphene sheet, being in accordance with the recent
experimental verification of ``electron--hole puddles'' in G/a-SiO$_2$
systems~\cite{stampferAPL2007,martinNatPhys2008}.   Here we  can infer
that  the inhomogeneous  a-SiO$_2$  surface topology  rules the  total
charge  density  fluctuation of  the  adsorbed  graphene sheet.   Such
fluctuation  on  the electronic  charge  density  distribution on  the
graphene  sheet should play  an important  role on  the experimentally
verified reduction  of the electronic  mobility for graphene  lying on
a-SiO$_2$ surface~\cite{lemmeIEEE2007}.

\begin{figure}[h]
\includegraphics[width=6.5cm]{pdos3.eps}
\caption{Electronic band structure (a) and the projected density of
  states (b) of G/SiO$_2$ system. Dark shaded regions indicate the
  total density of states, and the light shaded region indicate the
  projected density of states of a-SiO$_2$ surface, and the solid line
  indicates the projected density of state of (adsorbed) graphene
  sheet. The inset in (a) present the energy dispersion of the
  partially occupied level $c1$.}
\label{pdos}
\end{figure}

There is a net electronic charge density gain of about 3~$\times
10^{13}~e/{\rm cm^2}$, Fig.~\ref{structure}(d), thus, indicating a
n-type doping on the adsorbed graphene sheet.  Our electronic band
structure calculation, Fig.~\ref{pdos}(a), confirms statement above.
The shape of the Dirac cone has been maintained, since the graphene
sheet weakly interacts with the a-SiO$_2$ substrate.  Due to the
electronic charge transfer from a-SiO$_2$ to the graphene sheet, the
Dirac point lies below the Fermi energy ($\rm E_F$). The n-type doping
occurs due to the presence of a partially occupied state ($c1$) below
the conduction band minimum (CBM) pinning the Fermi level above the
Dirac point, see inset of Fig.~\ref{pdos}(a).  Figure~\ref{pdos}(b)
presents the density of states (DOS) of G/a-SiO$_2$ (gray region), and
separately the contribution from the adsorbed graphene sheet (solid
line), and the a-SiO$_2$ surface (light gray region). We find, in
contrast with previous investigations~\cite{romeroACSNano2008}, that
$c1$ comes from the formation of three-fold coordinated O atoms
(O$_{\rm three-fold}$) in a-SiO$_2$. Those O$_{\rm three-fold}$ sites
may occur when we have a slightly hight local concentration of Si
atoms or nearby oxygen vacancies in
a-SiO$_2$~\cite{luPRL2002,sokolovNonCrystSol1997,winerPRB1988,feustonJChemPhys1989}.
The Si--O$_{\rm three-fold}$ bond lengths (1.82--1.86~\AA) are
stretched in comparison with the ones of the two fold coordinated O
atoms, $\sim$1.66~\AA, Fig.~\ref{structure}(c).  Those larger bond
lengths reduce the binding energy of the electronic states along
Si--O$_{\rm three-fold}$, giving rise to the partially occupied state,
$c1$, just below the CBM and pinning the Fermi level. We find that the
electronic states of $c1$ are localized on the Si atoms nearby O$_{\rm
  three-fold}$.  In order to verify the adequacy of the proposed
origin of $c1$, we calculate the electronic structure of G/a-SiO$_2$
by suppressing the presence of O$_{\rm three-fold}$, {\it i.  e.} we
have only two fold coordinated O atoms within the a-SiO$_2$ slab.  In
this case, we do not find the formation of $c1$, and the Dirac point
lies on the Fermi level.  That is, the adsorbed graphene sheet becomes
neutral.  Those results support the proposed role of O$_{\rm
  three-fold}$ on the non-covalent n-type doping of graphene on
a-SiO$_2$ substrate.  It is worth to note that, for the neutral
graphene sheet, we find the same picture for the electron- and
hole-rich regions, as depicted in Fig.~\ref{structure}(d), however
with different values of $\Delta\rho$.  The total charge density
fluctuation on the neutral graphene is about 25\% smaller in
comparison with the (n-type) doped ones.

 Finally, we have considered  the presence of oxygen vacancies
  (O$_{\sf V}$) in a-SiO$_2$. We have examined three plausible O$_{\sf V}$
configurations, on the topmost surface site, and on the subsurface
sites, indicated as {\sf V1}--{\sf V3} in Fig.~\ref{structure}(a),
respectively. We verify that those vacancies do not give rise to a
donor level on the a-SiO$_2$ surface. For {\sf V1} and {\sf V2} we
find the formation of Si-dimers on the surface (Si--Si bond length of
2.24~\AA) inducing the formation of occupied (empty) electronic states
near the valence (conduction) band maximum (minimum). Whereas {\sf V3}
gives rise to a deep occupied state within the a-SiO$_2$ bandgap.

In conclusion, we find that vdW interactions rules the graphene
adsorption on the a-SiO$_2$ surface.  Through a charge density map on
the adsorbed graphene sheet we verify that the inhomogeneous a-SiO$_2$
surface topology promotes an electronic charge density displacement on
the graphene, giving rise to the experimentally verified
``electron-hole puddles''~\cite{stampferAPL2007,martinNatPhys2008}.
Such electron-rich and hole-rich regions should contribute to the
experimentally verified reduction on the electronic mobility of
graphene sheets lying on a-SiO$_2$ surfaces.  We find a net charge
density gain on the adsorbed graphene sheet, thus, characterizing a
n-type non-covalent doping on the graphene.  The electronic charge
transfer from a-SiO$_2$ to the graphene sheet comes from the formation
of a partially occupied state above the Dirac point, which is
attributed to the formation of three-fold coordinated oxygen atoms in
a-SiO$_2$.  Those three-fold coordinated oxygen atoms may occurs upon
slightly hight local concentration of Si atoms or nearby oxygen
vacancies.

\acknowledgments

The authors thanks Dr. Renato Pontes for fruitful discussions, the
financial support from the Brazilian agencies CNPq, CAPES, FAPEMIG,
FAPESP, CNPq/INCT and the computational support from CENAPAD/SP.


The Spectroscopic Age of 47 Tuc

Abstract: High signal-to-noise integrated spectra of the metal-rich globular cluster 47 Tuc, spanning the H-gamma(HR) and Fe4668 line indices, have been obtained. The combination of these indices has been suggested (Jones & Worthey 1995, ApJ, 446, L31) as the best available mechanism for cleanly separating the age-metallicity degeneracy which hampers the dating of distant, unresolved, elliptical galaxies. For the first time, we apply this technique to a nearby spheroidal system, 47 Tuc, for which independent ages, based upon more established methods, exist. Such an independent test of the technique's suitability has not been attempted before, but is an essential one before its application to more distant, unresolved, stellar populations can be considered valid. Because of its weak series of Balmer lines, relative to model spectra, our results imply a spectroscopic ``age'' for 47 Tuc well in excess of 20 Gyr, at odds with the colour-magnitude diagram age of 14+/-1 Gyr. The derived metal abundance, however, is consistent with the known value. Emission ``fill-in'' of the H-gamma line as the source of the discrepancy cannot be entirely excluded by existing data, although the observational constraints are restrictive.

\section{Introduction}
\label{introduction}

The star formation history of elliptical galaxies still remains, at many
levels, a mystery.  The fundamental uncertainty afflicting
attempts to determine the ages of such (presumably) old stellar populations
lies in the unfortunate fact-of-life that age and metallicity have almost
identical effects upon broad-band colours and on most line indices
(O'Connell 1976).  This age-metallicity \it degeneracy \rm (sometimes referred
to as Worthey's
(1994) ``3/2-rule''\footnote{The 3/2-rule simply states that if two
stellar populations have $\Delta\log({\rm Age})/\Delta\log({\rm
Metallicity})\equiv -3/2$, their broad-band colours and (most)
line indices will be nearly identical.})
has led many ({\rm e.g.,}\ O'Connell 1976; Worthey 1994; Rose 1994) to
search for individual absorption features which may be
sensitive to either age \it or \rm metallicity, but not both.

Unfortunately, line indices have not proven to be the panacea once hoped for.
As Worthey (1994) has shown (see his Table 6), the vast majority of
low-resolution lines defined by the Lick IDS survey ({\rm e.g.,}\ Worthey {\rm et~al.}\
1994) lie uncomfortably close to the 3/2-rule degeneracy ``zone''.  An
encouraging result to come out of the Worthey (1994) IDS study, however,
was the
identification of one particularly useful, and strong,
indicator, apparently free of any strong age influence.  Fe4668, whose bandpass
spans $\sim 4635\rightarrow 4722$\,\AA, tracks not only Fe, but contains a
strong
contribution from C$_2$, Ti, Cr, and Mg, and as Worthey, Trager \& Faber (1996)
note, is the best current candidate tracer for the mean metallicity ``Z'',
almost exclusively sensitive to this alone.

On the other hand, identifying a matching ``pure'' age discriminator has proven
more difficult.  H$_\alpha$ and H$_\beta$ have long been popular indices, in
this regard ({\rm e.g.,}\ O'Connell 1976), but as stressed by Gonz\'alez
(1993), both lines can be easily distorted (even filled in) by emission from
ionized gas, \it severely \rm
compromising their usage.  Higher-order Balmer emission
lines ({\rm e.g.,}\ H$_\delta$ and H$_\gamma$)
are admittedly weaker, yet free of such ``emission'' complications, and appear
to provide a very sensitive age diagnostic.  Rose (1994) recognized this latter
fact, and designed a set of higher-resolution line indices
in order to take advantage of these age discriminating
``powers''.  The Rose system employs a 3.74\AA\ window over which the line
indices are measured; in contrast, the Lick IDS system (typically) uses an
$\sim$40\AA\ window.  One index in particular,
H$_{\gamma_{\rm HR}}$ (centred at 4340.5\,\AA), was strongly sensitive to age
alone.\footnote{The ``HR'' subscript stands for ``high resolution'', in order
to distinguish the Rose-style equivalent widths, with 3.74\,\AA\, bandpasses,
from the Lick-style indices which are an order of magnitude broader.}

The near-orthogonal behaviour, in the age-metallicity plane,
of the H$_{\gamma_{\rm HR}}$ and Fe4668 line
indices, was demonstrated by Jones \& Worthey (1995; Figure 2), and suggests
a viable mechanism for breaking the aforementioned age-metallicity
degeneracy.  An extension of this Lick+Rose index comparison was
undertaken by Jones (1995,1996ab), but now concentrating solely upon pairs
of high resolution Rose indices -- {\rm e.g.,}\
Jones examined the H$_{\gamma_{\rm HR}}$--Ca\,I$_{\rm HR}$ age-metallicity
diagnostic diagram.

Jones \& Worthey (1995) stress several caveats regarding their
proposed optimal age-dating technique.  First, while the parameter-space
coverage ({\rm e.g.,}\ T$_{\rm eff}$, $\log g$, and [Fe/H]) of their stellar library,
upon which the model sequences are ultimately based, is good, it is
not complete.  Specifically, there is an admitted lack of metal-poor stars,
especially hot horizontal branch stars, which could contribute to the Balmer
line absorption.  Because of this,
Jones \& Worthey emphasize that their sequences are not directly
applicable to systems with metallicities outside $-1.0\le{\rm [Fe/H]}\le +0.6$.
Second, in order to apply this technique accurately
to integrated spectra of spheroidal systems, very high signal-to-noise ratios
($\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$} 100$), at moderately high spectral resolutions ($\sim 1$\,\AA), are a
necessity.

To this end, Jones \& Worthey (1995) present such a spectrum for M32,
and find an age of $\sim 6\rightarrow 7$ Gyr, in agreement with
O'Connell (1980) and Rose (1994).  While this result is a necessary proof
of the technique, it is not a sufficient one, as it effectively compares
two population \it
synthesis \rm methods, both of which are \it indirectly \rm determining
an age estimate.  Optimally, any testing of a new, and
potentially exciting, technique should be completely ``independent'', at some
level.  Unfortunately, there are simply no independently age-dated ellipticals
known.
On the other hand, there does exist another family of systems which are
ideally suited to such a study -- the metal-rich globular clusters of the Milky
Way. \it
Metal-rich globular clusters are the only practical systems with which to
independently assess
the reliability of the H$_{\gamma_{\rm HR}}$--Fe4668 age-dating technique. \rm
This comes about because
accurate ages have already been determined for at least four
clusters in the metal-rich regime to which the Jones \& Worthey technique is
applicable ({\rm e.g.,}\ Richer {\rm et~al.}\ 1996).
With its high metallicity ({\rm i.e.,}\
[Fe/H]$\approx -0.7$) and colour-magnitude diagram (CMD)-derived
age of $14\pm 1$ Gyr
(Richer {\rm et~al.}\ 1996), the southern globular cluster, 47~Tuc, is an ideal
laboratory within which to test this new age-dating technique.

Existing high-resolution integrated spectra of 47~Tuc in the literature
include those of Smith (1979) and Rose (1994).
In both cases though, the observations were
restricted to $\lambda<4400$\,\AA, meaning the metallicity-discriminant
Fe4668 line index could not be measured.

In Section
\ref{observations}, we briefly describe our new observations of 47~Tuc,
which now provide simultaneous coverage of both the
age-sensitive feature
H$_{\gamma_{\rm HR}}$, and the
metallicity-sensitive features Ca\,I$_{\rm HR}$ and Fe4668.
Such data will allow us to test how
accurately the spectroscopic age of this simple stellar population agrees
with its accurately determined CMD age.
The results of our analysis are
discussed in Section \ref{discussion}, and summarised in Section \ref{summary}.

\section{Observations and Analysis}
\label{observations}

Two high signal-to-noise ($S/N\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$} 100$) integrated spectra of
47~Tuc were acquired during the nights of October 20 and 21,
1996, using
the Coud\'e spectrograph, with the 32$^{\prime\prime}$ Schmidt camera
(Grating E; 1$^{\rm st}$ order),
on the 74$^{\prime\prime}$ telescope at Mount Stromlo Observatory,
was employed.
The $1^\prime.1$ long slit was
scanned back and forth (E-W) across the slit, to limits of $\pm 0^\prime.55$
from the cluster center (total coverage
in excess of 2.5 core radii, ensuring stochastic
fluctuations due to finite sampling were minimized), for
$\sim 30$ minutes per exposure.
The 24 $\mu$m pixels of the
Tek 2k$\times$2k CCD translated to 0.48\,\AA\, per pixel -- {\rm i.e.,}\ 983\,\AA\, of
spectral coverage, which for a central wavelength setting of $\sim 4500$\,\AA,
allowed simultaneous determination of
Ca I$_{\rm HR}$, H$_{\gamma_{\rm HR}}$, and Fe4668 indices.  Companion
blank-sky exposures (again, of $\sim 30$ minutes duration,
taken $\sim 20^\prime$ to the north of
the cluster core) were acquired subsequent to each 47~Tuc exposure. The
relative strength of the sky versus cluster signal was negligible.
For calibration purposes, 13 stars (spanning
a range in metallicity and spectral type) in common with
Jones' Coud\'e Feed Spectral
Library (see Section 2.4 of Leitherer {\rm et~al.}\ 1996 and Jones 1996b)
and the Lick/IDS Library
(Worthey {\rm et~al.}\ 1994) were observed, as were two
spectrophotometric standards.  Standard IRAF tasks were employed for
the flat-fielding, flux, and wavelength calibration.  Observations taken
throughout the run of bright, metal-poor halo giants imply that sky
subtraction has not introduced any unforeseen spectral artifacts.
Because the central one-dimensional velocity dispersion of 47~Tuc is only
$\sigma_c=11.6\pm 1.4$\,km\,s$^{-1}$ (Meylan \& Mayor 1986), no correction for
internal dispersion was applied to the spectra; as Trager {\rm et~al.}\ (1998; Figure
3) demonstrate, such a correction, for globular clusters, is
negligible.

Before measuring the H$_{\gamma_{\rm HR}}$ and Ca\,I$_{\rm HR}$ line indices,
our observed standards
and target 47~Tuc spectra were first smoothed
with a $\sigma=103$\,km\,s$^{-1}$
Gaussian, and re-binned to 0.623\AA/pixel, in both cases
to parallel Jones (1996b).
The H$_{\gamma_{\rm HR}}$ and Ca\,I$_{\rm HR}$ equivalent widths are measured
as detailed in Section 3 of Jones \& Worthey (1995).  Flanking blue and red
pseudo-continuum peaks, usually associated with the immediately-adjacent blue
and red shoulders of the relevant absorption line, were identified; the
equivalent widths were then measured relative to the continuum formed by
connecting the two pseudo-continuum peaks.  For H$_{\gamma_{\rm HR}}$, the blue
pseudo-continuum was associated with the local maximum between 4329\,\AA\,
and 4336\,\AA, while the red pseudo-continuum was associated with the maximum
in the 4346$-$4350\,\AA\,
window; for Ca\,I$_{\rm HR}$, the correspondings windows
were 4218$-$4222\,\AA\, and, either 4230$-$4232\,\AA\, (early spectral types)
or 4242$-$4246\,\AA\, (later spectral types).
Figure \ref{fig:fig1} shows the pair of integrated 47~Tuc spectra in the
region of both the H$_{\gamma_{\rm HR}}$ and Ca\,I$_{\rm HR}$ indices; the
associated flanking pseudo-continua are similarly noted.
The equivalent width was then derived via a simple Simpson's Rule algorithm,
summing the area associated with the seven 0.623\,\AA\, pixels
centred on the minimum of the high resolution absorption
feature\footnote{A quirk of the line-measuring code employed by Jones
(1995,1996ab), itself an offspring of that used by Rose (1994, and earlier
papers), is that
the equivalent width derivation requires three pixels on either side of
absorption feature minimum (itself one pixel) -- {\rm i.e.,}\ while the pixel
centre-to-centre window is the 3.74\,\AA\ referred to throughout this
paper (and earlier Rose and Jones papers), the effective area over which the
integration is acting corresponds to seven pixels, or 4.36\,\AA.}
-- Gaussian-fitting to the line profiles was not employed in our analysis,
again, in order to maintain consistency with Jones (1996b).
The derived
equivalent widths of Figure \ref{fig:fig1} are
consistent, from night-to-night, at the $2\rightarrow 3$\% level.

\placefigure{fig:fig1}

For the Fe4668 measurement, the 13 standards and 47~Tuc spectra were smoothed
with a $\sigma=250$\,km\,s$^{-1}$ Gaussian, to match the Lick IDS resolution.
Adjacent blue and red continua were determined from the mean flux in the
4611.50$-$4630.25\,\AA\, and 4742.75$-$4756.500\,\AA\, windows, respectively.
Using the continuum defined by the line connecting the mean of these two
windows, the
equivalent width was then derived by integrating the absorption
feature over the 4634.00$-$4720.25\,\AA\, window.  The adopted windows are
marginally different from those defined by Worthey {\rm et~al.}\ (1994; Table 1), but
reflect those adopted in the 1997 edition of the Lick/IDS stellar catalog
(Trager {\rm et~al.}\ 1998).

Comparison of the equivalent widths derived from our 13 standards with those
tabulated in
the Coud\'e Feed Spectral Library (H$_{\gamma_{\rm HR}}$ and Ca\,I$_{\rm HR}$)
and Lick/IDS Library (Fe4668) revealed the two systems were on the
scale, and no further corrections were
needed to bring the systems into agreement -- {\rm i.e.,}\ both groups were measuring
identical equivalent widths for the same standards -- and therefore, our
measured 47~Tuc equivalent widths could be used \it as is \rm in comparing
against the model age-metallicity diagnostic diagrams of Jones \& Worthey
(1995) and Jones (1996a).  Evidence to this effect is presented in
Figure \ref{fig:fig2}, in which our standard star Ca\,I$_{\rm HR}$
equivalent widths (x-axis) are plotted against those from Jones' (1996b)
Coud\'e Feed Spectral Library (y-axis) - the mean offset (in the sense of
Jones$-$Present Study) is only $+0.010\pm 0.006$\AA, a factor of
$\sim$3$\times$ smaller than the typical uncertainty associated with
pseudo-continua placement and night-to-night variations (recall
Figure \ref{fig:fig1}).  Identical conclusions hold for the other line indices
in our study.
Our measured (and adopted)
47~Tuc equivalent widths are therefore
$1.032\pm 0.040$\,\AA\, (H$_{\gamma_{\rm HR}}$),
$1.042\pm 0.030$\,\AA\, (Ca\,I$_{\rm HR}$), and
$1.95\pm 0.45$\,\AA\, (Fe4668), where the uncertainties reflect those due to
pseudo-continuum placement, and, to a lesser degree, the scatter derived from
repeat measurements of the set of standards.

\placefigure{fig:fig2}

\section{Discussion}
\label{discussion}

Following the prescription outlined by Jones \& Worthey (1995) and
Jones (1996b), the Worthey (1994)
population synthesis models were used to predict H$_{\gamma_{\rm HR}}$,
Ca\,I$_{\rm HR}$, and Fe4668 line index strengths, as a function of age and
metallicity [Fe/H], for a grid of single burst stellar populations.
Twenty-four realizations, spanning ages $2\rightarrow
17$\,Gyr and metallicities [Fe/H]=$-0.50\rightarrow +0.50$, were made; the full
set of 24 models were then incorporated into the H$_{\gamma_{\rm
HR}}$--Ca\,I$_{\rm HR}$ age-metallicity diagnostic grid shown in Figure
\ref{fig:fig3}.  While not perfectly orthogonal, this pair of high resolution
Rose-style indices do provide adequate separation of age and
metallicity, with a minimum of degeneracy.  The four filled symbols represent
the same set of low luminosity ellipticals shown in the earlier analysis of
Jones (1996a), but now smoothed to the same $\sigma=103$\,km\,s$^{-1}$
resolution adopted by Jones (1996b) and the present study.  As discussed
previously (Jones 1996a), these four galaxies represent a (luminosity-weighted)
solar metallicity age sequence spanning the range $2\rightarrow 12$\,Gyr.

\placefigure{fig:fig3}

The open starred
symbol represents our new 47~Tuc data point; unlike the case for the
galaxies shown in Figure \ref{fig:fig3}, 47~Tuc has an independent
age\footnote{Salaris \& Weiss (1998), using $\alpha$-element enhanced stellar
models, have recently revised the age of 47~Tuc from the canonical
(scaled-solar abundance) Richer {\rm et~al.}\
(1996) value of $14\pm 1$\,Gyr, to the lower value of $9\pm 1$\,Gyr.  On the
surface this only worsens the spectroscopic versus CMD age discrepancy of
Figure \ref{fig:fig3}, although the comparison is perhaps not a fair one, since
Worthey's (1994) synthesis models are (necessarily) tied to models and stellar
data which are likewise of scaled-solar abundances.  Similarly, the
post-Hipparcos globular cluster distance
scale appears to favor a 2-3\,Gyr reduction
in the mean age of the oldest globulars (Chanoyer {\rm et~al.}\ 1998); again, this
only strengthens our conclusion regarding the spectroscopic age of 47~Tuc,
although to be fair, as before, we retain the Richer {\rm et~al.}\ (1996) result
as the most representative one against which to compare our results.}
($14\pm 1$\,Gyr) and metallicity ([Fe/H]$\approx -0.7$) determination, derived
from isochrone fitting to the CMD (Richer {\rm et~al.}\ 1996; Salaris \& Weiss 1998).
It should be immediately clear that the spectroscopic age of 47~Tuc, as
determined by the age-sensitive H$_{\gamma_{\rm HR}}$ line index is
inconsistent with its CMD age, at the $\sim$4$\sigma$
level, being $\sim$0.2\,\AA\, lower than expected for this 14\,Gyr old simple
stellar population.  This $\sim$0.2\,\AA\, H$_\gamma$ discrepancy was
previously noted by Rose (1994; Section 5.1), although not in the context of
age-dating its integrated spectrum.  On the other hand,
it is somewhat re-assuring that the
metallicity inferred from Figure \ref{fig:fig3} ({\rm i.e.,}\ [Fe/H]$<$-0.5) is
consistent with that derived from the CMD (Richer {\rm et~al.}\ 1996).

Figure \ref{fig:fig4} shows the
H$_{\gamma_{\rm HR}}$--Fe4668 age-metallicity diagnostic diagram as
presented by Jones \& Worthey (1995; Figure 2);  the same nineteen models
used in their analysis are incorporated in our grid.  The only (subtle)
difference between Figure \ref{fig:fig4}, and that of Jones \& Worthey,
is that, in the latter, the spectra were smoothed to
$\sigma=80$\,km\,s$^{-1}$ (to match the resolution of their M32 spectrum),
while in the present study, recall that
 a $\sigma=103$\,km\,s$^{-1}$ Gaussian was
employed (to match the spectra and models of Jones 1996b).  Again,
the model sequences of Figure \ref{fig:fig4} clearly demonstrate the
near perfect orthogonal nature of the H$_{\gamma_{\rm HR}}$ and Fe4668
line indices in the age-metallicity plane.

\placefigure{fig:fig4}

As in Figure \ref{fig:fig3}, the filled square represents the Jones \&
Worthey (1995) observation of M32, while the open starred symbol shows our
new 47~Tuc data point.  As was the case for the pair of Rose-style
high resolution line indices (Figure \ref{fig:fig3}), in the
H$_{\gamma_{\rm HR}}$ versus Fe4668, the anomalously low
H$_{\gamma_{\rm HR}}$ in the integrated spectrum of 47~Tuc results in a
spectroscopic age ($\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$}$20\,Gyr) which is clearly discrepant with
CMD age ($14\pm 1$\,Gyr).

Is there some way in which to reconcile the spectroscopic versus CMD
age discrepancy?  As discussed by Rose (1994), virtually every line index
in the integrated spectrum of 47~Tuc appears perfectly well-behaved,
except the anomalously strong CN bands and anomalously weak H$_\gamma$.
Jones (1996b) speculated that perhaps the H$_\gamma$ ``problem'' was
being caused by CH contamination in the H$_\gamma$ line; because the
central regions of
47~Tuc ({\rm i.e.,}\ the region over which our integrated spectrum
was obtained) is replete with CN-strong giants (Norris \& Freeman 1979;
Da Costa 1997),
and there exists the well-known CN-CH anticorrelation for CN-strong
giants (Norris, Freeman \& Da Costa 1984), a scenario was envisaged whereby
H$_\gamma$ might be weaker in CN-strong (CH-weak) stars.
Jones (1996b), though,
was able to show, using a sample of 11 CN-strong and CN-weak
giant and horizontal branch stars from M71 (a cluster similar to 47~Tuc),
that this initial hypothesis was not supported by the data.

The most obvious speculative source of the discrepancy is that some
unidentified source of
H$_\gamma$ \it emission \rm is acting to ``fill in'' the
H$_\gamma$ \it absorption \rm feature.  Hesser \& Shawl (1977) provide an upper
limit to the emission flux at H$_\alpha$ of F(H$_\alpha$)$\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$} 1.3\times
10^{-11}$\,erg\,cm$^{-2}$\,s$^{-1}$.  Here we have used their most restrictive
dataset - {\rm i.e.,}\ $C_1=26.3$\,s$^{-1}$ (from their 22/08/75 observations), coupled
with the definition F(H$_\alpha$)$\approx 5\times
10^{-13}C_1$\,erg\,cm$^{-2}$\,s$^{-1}$, from Smith, Hesser \& Shawl (1976).
Now, the core of 47~Tuc, with M$_V\approx +6$ (and ignoring reddening),
corresponds to a flux of
$\sim 1.6\times 10^{-11}$\,erg\,cm$^{-2}$\,s$^{-1}$\,\AA$^{-1}$.  So, the
Hesser \& Shawl upper limit of $1.3\times 10^{-11}$\,erg\,cm$^{-2}$\,s$^{-1}$
corresponds to $1.3/1.6=0.8$\,\AA\ of equivalent width at H$_\alpha$.
Several simplifying assumptions must be made to go from H$_\alpha$ to
H$_\gamma$ - specifically,
under case B conditions and a temperature of $10000$\,K,
radiative recombination leads to an intrinsic ratio
L(H$_\alpha$)/L(H$_\gamma$)$\approx 6.5$, which gets reduced by a factor 1.7,
reflecting the ratio of continuum flux at H$_\alpha$ to that at H$_\gamma$
(employing the Worthey 1994 models).  In other words, H$_\gamma$ is a factor of
$3.8\times$ weaker than H$_\alpha$; recalling the constraint at H$_\alpha$ of
$\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$}
0.8$\,\AA, this corresponds to $\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$}
0.21$\,\AA\ of equivalent width at H$_\gamma$.

From Figures \ref{fig:fig3} and \ref{fig:fig4}, we can see that the
H$_\gamma$ ``discrepancy'' amounts to $\sim 0.17$\,\AA, in that we \it a priori
\rm expected an equivalent width of $\sim 1.20$\,\AA, to match the known CMD
age of 47~Tuc, as opposed to the $\sim 1.03$\,\AA\ observed.  Recalling that
the inferred upper limit to emission fill-in at H$_\gamma$ (from the previous
paragraph) was $\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$} 0.21$\,\AA, we can see that the requisite $\sim
0.17$\,\AA\ of emission needed to reconcile the CMD and spectroscopic ages is
not entirely excluded by the Hesser \& Shawl (1977) observations, although
little room is left to maneuver should new H$_\alpha$ emission observations
push this limit even lower.

The situation regarding the anomalously low H$_\gamma$ index in 47~Tuc
remains as it was at the time of the Rose (1994) and Jones (1996b) studies -
{\rm i.e.,}\ 47~Tuc appears normal in virtually every line index, except H$_\gamma$,
and to date no satisfactory explanation has been put forth.  Emission fill-in
has not been entirely excluded, but the constraints are tight, as shown above.
Should such emission prove to be the ultimate resolution of the discrepancy
though, it will beg the question as to exactly which mechanism is responsible,
and why it differs from the local calibrating stars used to derive the model
sequences of Figures \ref{fig:fig3} and \ref{fig:fig4}.

\section{Summary}
\label{summary}

From the outset, our goal was a simple one - to provide an independent
test of the Jones \& Worthey (1995) H$_{\gamma_{\rm HR}}$--
Fe4668 diagnostic (and the Jones 1996a
H$_{\gamma_{\rm HR}}$--Ca\,I$_{\gamma_{\rm HR}}$ diagnostic) as a tool for
breaking the long-standing age-metallicity degeneracy plaguing integrated
spectra of stellar populations.  Such a test can only be done
(properly) upon systems with independently-derived ages and metallicities;
local metal-rich globular clusters with accurate, deep, CMDs,
provide an ideal test-bed.
Our new high $S/N$, high resolution, integrated spectrum of 47~Tuc shows
anomalously low H$_\gamma$ absorption, which when tied to the diagnostic
models of Jones \& Worthey (1995), shows, for the first time, that
the inferred spectroscopic age is well in excess
of 20\,Gyr. This is
clearly at odds with the CMD age of $14\pm 1$\,Gyr.  No suitable
explanation as to the source of this discrepancy has been found, despite
the best efforts of Rose (1994) and Jones (1996b).  Emission, amounting to
$\sim 0.17$\,\AA of equivalent width at H$_\gamma$, is required to reconcile
the discrepancy.  Such emission is not excluded by the Hesser \& Shawl (1977)
observations of 47~Tuc
(which set an upper limit of $\sim 0.21$\,\AA\ at H$_\gamma$), but
if it is the responsible mechanism, must be lurking just marginally
below the limit set by these early Fabry-Perot observations.  Clearly, a deeper
re-visit is called for.

It is true that the Jones \& Worthey (1995)
H$_{\gamma_{\rm HR}}$--Fe4668 age-metallicity diagnostic will remain an
invaluable tool for age-dating integrated spectra of early-type systems, but
perhaps the most important conclusion to take from our simple test of the
technique is that it is not infallible.
Clearly, it will be imperative to extend our analysis to a larger sample of
metal-rich globulars, Galactic open clusters, and star clusters in the
Magellanic Clouds -- such an analysis will
demonstrate whether 47~Tuc is a pathological case or
the inconsistent spectroscopic age derived here is indicative of
an endemic, and currently underappreciated, flaw in the technique.

\acknowledgments

We wish to thank Guy Worthey, Jim Rose, and Ken Freeman for their helpful
advice.  A special word of thanks is aimed at the referee, who went above and
beyond the call of duty in aiding us in the interpretation of the Hesser \&
Shawl (1977) observations.
DSM acknowledges the financial support provided
through the ANU/MSSSO Summer Research Fellowship Program.

The Spectroscopic Age of 47 Tuc