Two-color coherent control in photoemission from gold needle tips

Abstract: We demonstrate coherent control of photoemission from a gold needle tip using a two-color laser field. The relative phase between a fundamental field and its second harmonic imprints a strong modulation on the emitted photocurrent with up to 96.5 % contrast. The contrast as a function of the second harmonic intensity can be described by three interfering quantum pathways. Increasing the bias voltage applied to the tip reduces the maximum achievable contrast and modifies the weights of the involved pathways. Simulations based on the time-dependent Schrödinger equation reproduce the characteristic cooperative signal and its dependence on the second harmonic intensity, which further confirms the involvement of three emission pathways.

\section{Introduction: file preparation and submission}

The \verb"iopart" \LaTeXe\ article class file is provided to help authors prepare articles for submission to IOP Publishing journals.
  This document gives advice on preparing your submission, and specific instructions on how to use \verb"iopart.cls" to follow this advice.  You
  do not have to use \verb"iopart.cls"; articles prepared using any other common class and style files can also be submitted.
    It is not necessary to mimic the appearance of a published article.

The advice
on \LaTeX\ file preparation in this document applies to
the journals listed in table~\ref{jlab1}.  If your journal is not listed please go to the journal website via \verb"http://iopscience.iop.org/journals" for specific
submission instructions.

\begin{table}
\caption{\label{jlab1}Journals to which this document applies, and macros for the abbreviated journal names in {\tt iopart.cls}. Macros for other journal titles are listed in appendix\,A.}
\footnotesize
\begin{tabular}{@{}llll}
\br
Short form of journal title&Macro name&Short form of journal title&Macro name\\
\mr
2D Mater.&\verb"\TDM"&Mater. Res. Express&\verb"\MRE"\\
Biofabrication&\verb"\BF"&Meas. Sci. Technol.$^c$&\verb"\MST"\\
Bioinspir. Biomim.&\verb"\BB"&Methods Appl. Fluoresc.&\verb"\MAF"\\
Biomed. Mater.&\verb"\BMM"&Modelling Simul. Mater. Sci. Eng.&\verb"\MSMSE"\\
Class. Quantum Grav.&\verb"\CQG"&Nucl. Fusion&\verb"\NF"\\
Comput. Sci. Disc.&\verb"\CSD"&New J. Phys.&\verb"\NJP"\\
Environ. Res. Lett.&\verb"\ERL"&Nonlinearity$^{a,b}$&\verb"\NL"\\
Eur. J. Phys.&\verb"\EJP"&Nanotechnology&\verb"\NT"\\
Inverse Problems&\verb"\IP"&Phys. Biol.$^c$&\verb"\PB"\\
J. Breath Res.&\verb"\JBR"&Phys. Educ.$^a$&\verb"\PED"\\
J. Geophys. Eng.$^d$&\verb"\JGE"&Physiol. Meas.$^{c,d,e}$&\verb"\PM"\\
J. Micromech. Microeng.&\verb"\JMM"&Phys. Med. Biol.$^{c,d,e}$&\verb"\PMB"\\
J. Neural Eng.$^c$&\verb"\JNE"&Plasma Phys. Control. Fusion&\verb"\PPCF"\\
J. Opt.&\verb"\JOPT"&Phys. Scr.&\verb"\PS"\\
J. Phys. A: Math. Theor.&\verb"\jpa"&Plasma Sources Sci. Technol.&\verb"\PSST"\\
J. Phys. B: At. Mol. Opt. Phys.&\verb"\jpb"&Rep. Prog. Phys.$^{e}$&\verb"\RPP"\\
J. Phys: Condens. Matter&\verb"\JPCM"&Semicond. Sci. Technol.&\verb"\SST"\\
J. Phys. D: Appl. Phys.&\verb"\JPD"&Smart Mater. Struct.&\verb"\SMS"\\
J. Phys. G: Nucl. Part. Phys.&\verb"\jpg"&Supercond. Sci. Technol.&\verb"\SUST"\\
J. Radiol. Prot.$^a$&\verb"\JRP"&Surf. Topogr.: Metrol. Prop.&\verb"\STMP"\\
Metrologia&\verb"\MET"&Transl. Mater. Res.&\verb"\TMR"\\
\br
\end{tabular}\\
$^{a}$UK spelling is required; $^{b}$MSC classification numbers are required; $^{c}$titles of articles are required in journal references; $^{d}$Harvard-style references must be used (see section \ref{except}); $^{e}$final page numbers of articles are required in journal references.

\end{table}
\normalsize

Any special submission requirements for the journals are indicated with footnotes in table~\ref{jlab1}.
Journals which require references in a particular format will need special care if you are using BibTeX, and you might need to use a \verb".bst" file
that gives slightly non-standard output in order to supply any extra information required.  It is not
necessary to give references in the exact style of references used in published articles, as long as all of
the required information is present.

Also note that there is an incompatibility
between \verb"amsmath.sty" and \verb"iopart.cls" which cannot be completely worked around.  If your article relies
on commands in \verb"amsmath.sty" that are not available in \verb"iopart.cls", you may wish to consider using a different
class file.

Whatever journal you are submitting to, please look at recent published articles (preferably
articles in your subject area) to familiarize yourself with the features of the journal.  We do not demand
that your \LaTeX\ file closely resembles a published article---a generic `preprint' appearance of the sort
commonly seen on \verb"arXiv.org" is fine---but your submission should be presented
in a way that makes it easy for the referees to form an opinion of whether it is suitable for the journal.
The generic advice in this document---on what to include in an abstract, how best to present complicated
mathematical expressions, and so on---applies whatever class file you are using.

\subsection{What you will need to supply}
Submissions to our journals are handled via the ScholarOne web-based submission system.  When you submit
a new article to us you need only submit a PDF of your article.  When you submit a revised version,
we ask you to submit the source files as well.  Upon acceptance for publication we will use the source files to produce a proof of your article in the journal style.

\subsubsection{Text.}When you send us the source files for a revised version of your submission,
you should send us the \LaTeX\ source code of your paper with all figures read in by
the source code (see section \ref{figinc}).  Articles can be prepared using almost any version of \TeX\ or \LaTeX{},
not just \LaTeX\ with the class file \verb"iopart.cls".  You may split your \LaTeX\ file into several parts, but please show
which is the `master' \LaTeX\ file that reads in all of the other ones by naming it appropriately.  The `master'
\LaTeX\ file must read in all other \LaTeX\ and figure files from the current directory.  {\it Do not read in files from a different directory, e.g. \verb"\includegraphics{/figures/figure1.eps}" or
\verb"\include{../usr/home/smith/myfiles/macros.tex}"---we store submitted files
all together in a single directory with no subdirectories}.
\begin{itemize}
\item {\bf Using \LaTeX\ packages.} Most \LaTeXe\ packages can be used if they are
available in common distributions of \LaTeXe; however, if it is essential to use
a non-standard package then any extra files needed to process the article must
also be supplied.  Try to avoid using any packages that manipulate or change the standard
\LaTeX\ fonts: published articles use fonts in the Times family, but we prefer that you
use \LaTeX\ default Computer Modern fonts in your submission.  The use of \LaTeX\ 2.09, and of plain
\TeX\ and variants such as AMSTeX is acceptable, but a complete PDF of your submission should be supplied in these cases.
\end{itemize}
\subsubsection{Figures.} Figures should ideally be included in an article as encapsulated PostScript files
(see section \ref{figinc}) or created using standard \LaTeX\ drawing commands.
 Please name all figure files using the guidelines in section \ref{fname}.
We accept submissions that use pdf\TeX\ to include
PDF or bitmap figures, but please ensure that you send us a PDF that uses PDF version 1.4 or lower
(to avoid problems in the ScholarOne system).
You can do this by putting \verb"\pdfminorversion=4" at the very start of your TeX file.

\label{fig1}All figures should be included within the body of the text
at an appropriate point or grouped together with their captions at the end of the article. A standard graphics inclusion package such as \verb"graphicx" should be used for figure inclusion, and the package should be declared in the usual
way, for example with \verb"\usepackage{graphicx}", after the \verb"\documentclass" command.
Authors should avoid using special effects generated by including verbatim
PostScript code in the submitted \LaTeX\ file. Wherever possible, please try to use standard \LaTeX\ tools
and packages.

\subsubsection{References.\label{bibby}}
You can produce your bibliography in the standard \LaTeX\ way using the \verb"\bibitem" command. Alternatively
you can use BibTeX: our preferred  \verb".bst" styles are:

\begin{itemize}
\item For the numerical (Vancouver) reference style we recommend that authors use
 \verb"unsrt.bst"; this does not quite follow the style of published articles in our
 journals but this is not a problem.  Alternatively \verb"iopart-num.bst" created by Mark A Caprio
 produces a reference style that closely matches that in published articles.  The file is available from
\verb"http://ctan.org/tex-archive/biblio/bibtex/contrib/iopart-num/" .
\item For alphabetical (Harvard) style references we recommend that authors use the \verb"harvard.sty"
in conjunction with the \verb"jphysicsB.bst" BibTeX style file.  These, and accompanying documentation, can be downloaded
from \penalty-10000 \verb"http://www.ctan.org/tex-archive/macros/latex/contrib/harvard/".
Note that the \verb"jphysicsB.bst" bibliography style does not include article titles
in references to journal articles.
To include the titles of journal articles you can use the style \verb"dcu.bst" which is included
in the \verb"harvard.sty" package.  The output differs a little from the final journal reference
style, but all of the necessary information is present and the reference list will be formatted
into journal house style as part of the production process if your article is accepted for publication.
\end{itemize}

\noindent Please make sure that you include your \verb".bib" bibliographic database file(s) and any
\verb".bst" style file(s) you have used.

\subsection{\label{copyright}Copyrighted material and ethical policy} If you wish to make use of previously published material for which you do not own the copyright then you must seek permission from the copyright holder, usually both the author and the publisher.  It is your responsibility to obtain copyright permissions and this should be done prior to submitting your article. If you have obtained permission, please provide full details of the permission granted---for example, copies of the text of any e-mails or a copy of any letters you may have received. Figure captions must include an acknowledgment of the original source of the material even when permission to reuse has been obtained.  Please read our ethical policy before writing your article.

\subsection{Naming your files}
\subsubsection{General.}
Please name all your files, both figures and text, as follows:
\begin{itemize}
\item Use only characters from the set a to z, A to Z, 0 to 9 and underscore (\_).
\item Do not use spaces or punctuation characters in file names.
\item Do not use any accented characters such as
\'a, \^e, \~n, \"o.
\item Include an extension to indicate the file type (e.g., \verb".tex", \verb".eps", \verb".txt", etc).
\item Use consistent upper and lower case in filenames and in your \LaTeX\ file.
If your \LaTeX\ file contains the line \verb"\includegraphics{fig1.eps}" the figure file must be called
\verb"fig1.eps" and not \verb"Fig1.eps" or \verb"fig1.EPS".  If you are on a Unix system, please ensure that
there are no pairs of figures whose names differ only in capitalization, such as \verb"fig_2a.eps" and \verb"fig_2A.eps",
as Windows systems will be unable to keep the two files in the same directory.
\end{itemize}
When you submit your article files, they are manipulated
and copied many times across multiple databases and file systems. Including non-standard
characters in your filenames will cause problems when processing your article.
\subsubsection{\label{fname}Naming your figure files.} In addition to the above points, please give each figure file a name which indicates the number of the figure it contains; for example, \verb"figure1.eps", \verb"figure2a.eps", etc. If the figure file contains a figure with multiple parts, for example figure 2(a) to 2(e), give it a name such as \verb"figure2a_2e.eps", and so forth.
\subsection{How to send your files}
Please send your submission via the ScholarOne submission system.  Go to the journal home
page, and use the `Submit an article' link on the right-hand side.

\section{Preparing your article}

\subsection{Sample coding for the start of an article}
\label{startsample}
The code for the start of a title page of a typical paper in the \verb"iopart.cls" style might read:
\small\begin{verbatim}
\documentclass[12pt]{iopart}
\begin{document}
\title[The anomalous magnetic moment of the
neutrino]{The anomalous magnetic moment of the
neutrino and its relation to the solar neutrino problem}

\author{P J Smith$^1$, T M Collins$^2$,
R J Jones$^3$\footnote{Present address:
Department of Physics, University of Bristol, Tyndalls Park Road,
Bristol BS8 1TS, UK.} and Janet Williams$^3$}

\address{$^1$ Mathematics Faculty, Open University,
Milton Keynes MK7~6AA, UK}
\address{$^2$ Department of Mathematics,
Imperial College, Prince Consort Road, London SW7~2BZ, UK}
\address{$^3$ Department of Computer Science,
University College London, Gower Street, London WC1E~6BT, UK}
\ead{williams@ucl.ac.uk}

\begin{abstract}
...
\end{abstract}
\keywords{magnetic moment, solar neutrinos, astrophysics}
\submitto{\jpg}
\maketitle
\end{verbatim}
\normalsize

At the start of the \LaTeX\ source code please include
commented material to identify the journal, author, and (if you are sending a revised
version or a resubmission) the reference number that the journal
has given to the submission. The first non-commented line should be
\verb"\documentclass[12pt]{iopart}"  to load the preprint class
file.  The normal text will be in the Computer Modern 12pt font.
It is possible to specify 10pt font size by passing the option \verb"[10pt]" to the class file.
Although it is possible to choose a font other than Computer Modern by loading external packages, this is not recommended.

The article text begins after \verb"\begin{document}".
Authors of very long articles may find it convenient to separate
their article into a series of \LaTeX\ files each containing one section, and each of which is called
in turn by the primary file.  The files for each section should be read in from the current directory;
please name the primary file clearly so that we know to run \LaTeX\ on this file.

Authors may use any common \LaTeX\ \verb".sty" files.
Authors may also define their own macros and definitions either in the main article \LaTeX\ file
or in a separate \verb".tex" or \verb".sty" file that is read in by the
main file, provided they do not overwrite existing definitions.
It is helpful to the production staff if complicated author-defined macros are explained in a \LaTeX\ comment.
The article class \verb"iopart.cls" can be used with other package files such
as those loading the AMS extension fonts
\verb"msam" and \verb"msbm", which provide the
blackboard bold alphabet and various extra maths symbols as well as symbols useful in figure
captions.  An extra style file \verb"iopams.sty" is provided to load these
packages and provide extra definitions for bold Greek letters.

\subsection{\label{dblcol}Double-column layout}
The \verb"iopart.cls" class file produces single-column output by default, but a two-column layout can be obtained by
using \verb"\documentclass[10pt]" at the start of the file and \verb"\ioptwocol" after the \verb"\maketitle" command.  Two-column output will begin
on a new page (unlike in published double-column articles, where the two-column material
starts on the same page as the abstract).

In general we prefer to receive submissions in single-column format even for journals
published in double-column style; however, the \verb"\ioptwocol" option may be useful to test figure sizes
and equation breaks for these journals.  When setting material
in two columns you can use the asterisked versions of \LaTeX\ commands such as \verb"\begin{figure*} ... \end{figure*}"
to set figures and tables across two columns.  If you have any problems or any queries about producing two-column output, please contact us at \verb"submissions@iop.org".

\section{The title and abstract page}
If you use \verb"iopart.cls", the code for setting the title page information is slightly different from
the normal default in \LaTeX.  If you are using a different class file, you do not need to mimic the appearance of
an \verb"iopart.cls" title page, but please ensure that all of the necessary information is present.

\subsection{Titles and article types}
The title is set using the command
\verb"\title{#1}", where \verb"#1" is the title of the article. The
first letter
of the title should be capitalized with the rest in lower case.
The title appears in bold case, but mathematical expressions within the title may be left in light-face type.

If the title is too long to use as a running head at the top of each page (apart from the
first) a short
form can be provided as an optional argument (in square brackets)
before the full title, i.e.\ \verb"\title[Short title]{Full title}".

For article types other than papers, \verb"iopart.cls"
has a generic heading \verb"\article[Short title]{TYPE}{Full title}"
and some specific definitions given in table~\ref{arttype}. In each case (apart from Letters
to the Editor and Fast Track Communications) an
optional argument can be used immediately after the control sequence name
to specify the short title; where no short title is given, the full title
will be used as the running head.  Not every article type has its own macro---use \verb"\article" for
any not listed.  A full list of the types of articles published by a journal is given
in the submission information available via the journal home page.
The generic heading could be used for
articles such as those presented at a conference or workshop, e.g.
\small\begin{verbatim}
\article[Short title]{Workshop on High-Energy Physics}{Title}
\end{verbatim}\normalsize
Footnotes to titles may be given by using \verb"\footnote{Text of footnote.}" immediately after the title.
Acknowledgment of funding should be included in the acknowledgments section rather than in a footnote.

\begin{table}
\caption{\label{arttype}Types of article defined in the {\tt iopart.cls}
class file.}
\footnotesize\rm
\begin{tabular*}{\textwidth}{@{}l*{15}{@{\extracolsep{0pt plus12pt}}l}}
\br
Command& Article type\\
\mr
\verb"\title{#1}"&Paper (no surtitle on first page)\\
\verb"\ftc{#1}"&Fast Track Communication\\
\verb"\review{#1}"&Review\\
\verb"\topical{#1}"&Topical Review\\
\verb"\comment{#1}"&Comment\\
\verb"\note{#1}"&Note\\
\verb"\paper{#1}"&Paper (no surtitle on first page)\\
\verb"\prelim{#1}"&Preliminary Communication\\
\verb"\rapid{#1}"&Rapid Communication\\
\verb"\letter{#1}"&Letter to the Editor\\
\verb"\article{#1}{#2}"&Other articles\\\ & (use this for any other type of article; surtitle is whatever is entered as {\tt
\#1})\\
\br
\end{tabular*}
\end{table}

\subsection{Authors' names and addresses}
For the authors' names type \verb"\author{#1}",
where \verb"#1" is the
list of all authors' names. Western-style names should be written as initials then
family name, with a comma after all but the last
two names, which are separated by `and'. Initials should {\it not} be followed by full stops. First (given) names may be used if
desired.  Names in Chinese, Japanese and Korean styles should be written as you want them to appear in the published article. Authors in all IOP Publishing journals have the option to include their names in Chinese, Japanese or Korean characters in addition to the English name: see appendix B for details.

If the authors are at different addresses a superscripted number, e.g. $^1$, \verb"$^1$", should be used after each
name to reference the author to his/her address.
If an author has additional information to appear as a footnote, such as
a permanent address, a normal \LaTeX\ footnote command
should be given after the family name and address marker
with this extra information.

The authors' affiliations follow the list of authors.
Each address is set by using
\verb"\address{#1}" with the address as the single parameter in braces.
If there is more
than one address then the appropriate superscripted number, followed by a space, should come at the start of
the address.

E-mail addresses are added by inserting the
command \verb"\ead{#1}" after the postal address(es) where \verb"#1" is the e-mail address.
See section~\ref{startsample} for sample coding. For more than one e-mail address, please use the command
\verb"\eads{\mailto{#1}, \mailto{#2}}" with \verb"\mailto" surrounding each e-mail address.  Please ensure
that, at the very least, you state the e-mail address of the corresponding author.

\subsection{The abstract}
The abstract follows the addresses and
should give readers concise information about the content
of the article and indicate the main results obtained and conclusions
drawn. It should be self-contained---there should be no references to
figures, tables, equations, bibliographic references etc.  It should be enclosed between \verb"\begin{abstract}"
and \verb"\end{abstract}" commands.  The abstract should normally be restricted
to a single paragraph of around 200 words.

\subsection{Subject classification numbers}
We no longer ask authors to supply Physics and Astronomy Classification System (PACS)
classification numbers.  For submissions to {\it Nonlinearity}\/ we ask that you should
supply Mathematics Subject Classification (MSC) codes.  MSC numbers are included after the abstract
using \verb"\ams{#1}".

The command
\verb"\submitto{#1}" can be inserted, where \verb"#1" is the journal name written in full or the appropriate control sequence as
given in table~\ref{jlab1}. This command is not essential to the running of the file and can be omitted.

\subsection{Keywords}
Keywords are required for all submissions. Authors should supply a minimum of three (maximum seven) keywords appropriate to their article as a new paragraph starting \verb"\noindent{\it Keywords\/}:" after the end of the abstract.

\subsection{Making a separate title page}
To keep the header material on a separate page from the
body of the text insert \verb"\maketitle" (or \verb"\newpage") before the start of the text.
If \verb"\maketitle" is not included the text of the
article will start immediately after the abstract.

\section{The text}
\subsection{Sections, subsections and subsubsections}
The text of articles may be divided into sections, subsections and, where necessary,
subsubsections. To start a new section, end the previous paragraph and
then include \verb"\section" followed by the section heading within braces.
Numbering of sections is done {\it automatically} in the headings:
sections will be numbered 1, 2, 3, etc, subsections will be numbered
2.1, 2.2,  3.1, etc, and subsubsections will be numbered 2.3.1, 2.3.2,
etc.  Cross references to other sections in the text should, where
possible, be made using
labels (see section~\ref{xrefs}) but can also
be made manually. See section~\ref{eqnum} for information on the numbering of displayed equations. Subsections and subsubsections are
similar to sections but
the commands are \verb"\subsection" and \verb"\subsubsection" respectively.
Sections have a bold heading, subsections an italic heading and
subsubsections an italic heading with the text following on directly.
\small\begin{verbatim}
\section{This is the section title}
\subsection{This is the subsection title}
\end{verbatim}\normalsize

The first section is normally an introduction,  which should state clearly
the object of the work, its scope and the main advances reported, with
brief references to relevant results by other workers. In long papers it is
helpful to indicate the way in which the paper is arranged and the results
presented.

Footnotes should be avoided whenever possible and can often be included in the text as phrases or sentences in parentheses. If required, they should be used only for brief notes that do not fit conveniently into the text. The use of
displayed mathematics in footnotes should be avoided wherever possible and no equations within a footnote should be numbered.
The standard \LaTeX\ macro \verb"\footnote" should be used.  Note that in \verb"iopart.cls" the \verb"\footnote" command
produces footnotes indexed by a variety of different symbols,
whereas in published articles we use numbered footnotes.  This
is not a problem: we will convert symbol-indexed footnotes to numbered ones during the production process.

\subsection{Acknowledgments}
Authors wishing to acknowledge assistance or encouragement from
colleagues, special work by technical staff or financial support from
organizations should do so in an unnumbered `Acknowledgments' section
immediately following the last numbered section of the paper. In \verb"iopart.cls" the
command \verb"\ack" sets the acknowledgments heading as an unnumbered
section.

Please ensure that you include all of the sources of funding and the funding contract reference numbers that you are contractually obliged to acknowledge. We often receive requests to add such information very late in the production process, or even after the article is published, and we cannot always do this. Please collect all of the necessary information from your co-authors and sponsors as early as possible.

\subsection{Appendices}
Technical detail that it is necessary to include, but that interrupts
the flow of the article, may be consigned to an appendix.
Any appendices should be included at the end of the main text of the paper, after the acknowledgments section (if any) but before the reference list.
If there are
two or more appendices they should be called Appendix A, Appendix B, etc.
Numbered equations will be in the form (A.1), (A.2), etc,
figures will appear as figure A1, figure B1, etc and tables as table A1,
table B1, etc.

The command \verb"
\section{Introduction}
Two-color laser fields formed by a strong fundamental laser pulse and its weak, phase-locked second harmonic became an important tool to manipulate and probe the electron emission dynamics in atomic and molecular gases. Complex and highly asymmetric waveforms can be generated by tuning the relative phase between both fields, changing their individual intensities or polarization angles. These sculpt fields can manipulate and probe the yield \cite{schumacher1994,muller1990}, angular structure \cite{yin1992,yin1995two,skruszewicz2015}, interference in momentum spectrum \cite{xie2012} and trajectories \cite{shafir2012} of emitted electrons with applications in the generation of high-harmonic \cite{mashiko2008} and terahertz \cite{dai2009} radiation. Metallic nanostructures, on the other hand, play an important technological role as DC and laser-triggered electron sources \cite{kruger2018}. Their usage ranges from high-resolution electron microscopes \cite{spence2013} and compact x-ray sources \cite{swanwick2014,graves2012} to recently developed on-chip, light-driven electronics \cite{krausz2014,rybka2016,bionta2020}. Thus, it is highly desirable to combine metallic solid-state systems with tailored electron emission in two-color laser fields.

So far, it has been shown that sweeping the relative phase between a fundamental field and its second harmonic modulates the current emitted from a tungsten needle tip with an unprecedented contrast of up to 97.5\,\% \cite{foerster2016,paschen2017,li2021}. The high contrast in this so called coherent control scheme applied to a tungsten needle tip can be well explained by the interference of two quantum pathways. Tungsten has many advantageous material properties such as the highest melting point among all metals, but it is a non-plasmonic material in contrast to, for example, gold. The creation of propagating surface plasmons and localized surface plasmons in the case of antenna structures, however, is often desired in light-matter interaction. The plasmon resonance helps to additionally increase optical near-fields at nanostructures \cite{dombi2020,thomas2013,thomas2015} and therefore reduces the required incident field strengths for non-linear photoemission. Propagating surface plasmons, on the other hand, can spatially separate the optical excitation from the electron emission as used in nanofocused electron point sources \cite{vogelsang2015,mueller2016}.

Here, we show that the coherent control scheme can be applied to the plasmonic material gold, too. A visibility of up to $96.5\%$ is obtained in the emitted photocurrent from a gold needle tip. In order to explain the scaling of the coherent signal with the second harmonic intensity, an additional, third pathway needs to be included. Simulations based on the time-dependent Schrödinger equation (TDSE) \cite{seiffert2018} reproduce the observed scaling and further validate the three-pathway model. Furthermore, we find that the achievable visibility can be suppressed by increasing the applied bias voltage \cite{paschen2017}, which is accommodated by changing contributions of the individual quantum pathways to the coherent signal.

\section{Experimental setup}
\begin{figure}[h!]
	\centering
	\includegraphics[width=0.75\linewidth]{Figures/Setup1}
	\caption{Experimental setup. Multi-cycle laser pulses centered at 1560 nm ($\omega$) and their second harmonic ($2\omega$) are collinearly focused onto the apex of a gold needle tip by an off-axis parabolic mirror within an UHV vessel. The emitted electrons are detected by a multi-channel plate (MCP) and counted as function of the optical delay $\tau$ between both laser fields. A bias voltage $\textrm{U}_\textrm{DC}$ is applied to the gold tip. Inset: Scanning electron micrograph of the gold tip with deduced apex radius of 25 nm and half-opening angle of $12^\circ$.}
	\label{fig:setup}
\end{figure}

The two-color field used in the experiment is composed of fundamental pulses centered around 1560\,nm with 74\,fs pulse duration and their second harmonics centered at $780\,$nm with $64\,$fs pulse duration. The fundamental pulses are delivered by an Erbium-doped fiber laser at a repetition rate of $100\,$MHz and the second harmonic is generated in a $100\,\mu$m thick $\beta$-barium borate crystal. Due to the parametric process the second harmonic has a fixed phase relation with respect to the fundamental field even if the carrier-envelope phase is not stabilized. A dichroic Mach-Zehnder type interferometer provides a tuneable optical delay between both fields, allows to adjust the individual intensities and to match the polarizations of both colors with the symmetry axis of the tip. A more detailed description of the optical setup is given in \cite{foerster2016}. The two-color field is tightly focused onto the apex of a gold needle tip inside a UHV vessel with a base pressure of $\approx 1 \times 10^{-10}\,$hPa as indicated in figure \ref{fig:setup}. The beam waists of the fundamental and second harmonic fields are $w_{\omega} = 4.0 \,\mu$m and $w_{2\omega} = 3.1 \,\mu$m ($1/e^2$ intensity radii). Electrons emitted from the tip are counted by a multi-channel plate (MCP) detector and the tip can be biased with a static voltage $\textrm{U}_\textrm{DC}$.

The gold tip was etched from a 0.1 mm thick poly-crystalline gold wire using a 90\% saturated aqueous potassium chloride solution \cite{eisele2011}. A tip apex radius of 25 nm and half-opening angle of $12^\circ$ can be estimated from a scanning electron micrograph (inset in figure \ref{fig:setup}). Nearfield enhancement factors of $FE_{\omega} = 6.5$ and $FE_{2\omega} = 5$ for the fundamental and second harmonic field are extracted from finite-difference time-domain simulations \cite{thomas2015} for the measured tip geometry. The maximum incident powers in this study are $P_\omega = 48\,$mW and $P_{2\omega} = 2.0\,$mW, which correspond to maximum near-field intensities of $I_{\omega} \approx 9 \times 10^{11} \textrm{ W}/\textrm{cm}^2$ and $I_{2\omega} \approx 4 \times 10^{10} \textrm{ W}/\textrm{cm}^2$. The work function for a clean gold surface $W = 5.4\,$eV \cite{kahn2016} results in minimal Keldysh parameters of $\gamma_\omega = 3.6$ and $\gamma_{2\omega} = 33$ placing our experiment clearly in the perturbative photoemission regime.

\begin{figure*}[h!]
	\centering
	\includegraphics[width=1\linewidth]{Figures/Experiment_power_sweep2a_2f}
	\caption{Coherent control of photocurrents from a gold needle tip. (a) Photocurrent as function of the delay between the $\omega$ and $2\omega$ field at intensities $I_{\omega} = 900 \textrm{ GW}/\textrm{cm}^2$ and $I_{2\omega}= 32 \, \textrm{GW}/\textrm{cm}^2$ with bias voltage $\textrm{U}_\textrm{DC} = -150 \, $V. (b) Count rate modulation with up to $96.5\% \pm 0.9 \%$ visibility in the center of the temporal overlap. The black line indicates a smoothing spline for the visibility estimation. (c) Amplitude of Fourier transform of the time domain data in (a). The cooperative signal is formed by Fourier components within regions of interest (ROIs) around the DC, $2\omega/(2\pi)$ and $4\omega/(2\pi)$ frequencies labeled $\textrm{ROI}_\textrm{DC}$, $\textrm{ROI}_{2\omega}$ and $\textrm{ROI}_{4\omega}$. (d) Inverse Fourier transform with additional Hilbert transform of $\textrm{ROI}_{2\omega}$ and $\textrm{ROI}_{4\omega}$. Gaussian fits according to equation (\ref{eq:Gauss}) are indicated by black dashed lines. (e) Peak heights $B_\textrm{DC}$, $B_{2\omega}$ and $B_{4\omega}$ of the Gaussian fits together with extracted visibility as function of the second harmonic intensity and admixture. Dashed line is a square root fit and solid lines are fits according to the three-pathway model. (f) Double-logarithmic representation of peak heights from (e) with power law exponents obtained by linear fits (solid lines).}
	\label{fig:power}
\end{figure*}

\section{Coherent control at gold needle tips}
Changing the delay between the fundamental and second harmonic field near the center of the temporal overlap causes strong oscillations in the electron count rate as shown in figure \ref{fig:power}(a).

The two-color field increases the count rate by more than five times or nearly fully suppresses the emission compared to separated pulses. By zooming into the center of the temporal overlap (figure \ref{fig:power}(b)) we observe an oscillation period of $T = 2.6$ fs  matching the period of the second harmonic field and determine a contrast of $96.5\% \pm 0.9 \%$. We use four local maxima and minima of a smoothing spline through the measured data to calculate the contrast or visibility $V$ and its standard deviation at the center of the time trace.

A Fourier analysis (figure \ref{fig:power}(c)) reveals that the time domain signal is composed of three characteristic frequency components centered around DC, $2\omega/(2\pi)$ and $4\omega/(2\pi)$ with $\omega= 192 \, \textnormal{THz} \times 2 \pi$ being the angular frequency of the fundamental field. To further identify the impact of the components on the time domain signal, we define regions of interest ($\textrm{ROI}_\textrm{DC}$, $\textrm{ROI}_{2\omega}$ and $\textrm{ROI}_{4\omega}$) around these Fourier components and transform them back individually. An additional Hilbert transformation provides their envelopes (figure \ref{fig:power}(d)), which can be well approximated by Gaussian functions
\begin{equation}
	G_i = A_i + B_i \exp(-(\tau-\tau_{c,i})^2/D_i^2) \label{eq:Gauss}
\end{equation}
with offsets $A_i$, peak heights $B_i$, delay shifts $\tau_{c,i}$ and widths $D_i$ for $i = \left\{ \textnormal{DC},2\omega,4\omega \right\}$. As the peak visibility is obtained in the center of the trace, it only depends on the peak heights $B_i$ and the offset $A_\textnormal{DC}$.

 In figure \ref{fig:power}(e) we observe a steep rise in visibility as the second harmonic intensity $I_{2\omega}$ is increased. $B_\textrm{DC}$ and $B_{4\omega}$ depend linearly on $I_{2\omega}$. $B_{2\omega}$ deviates from an exact square-root scaling. Extracting the corresponding power law exponents in figure \ref{fig:power}(f) confirms this observation.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.75\linewidth]{Figures/Pathway_model3}
	\caption{Quantum-pathway model. The pathway II (III) is created by exchanging two (four) photons of the fundamental field (red, weight $\alpha$) with one (two) photon of the second harmonic field (blue, weight $\beta$). The arrows depict partial amplitudes as indicated on the right. $\zeta^2 \neq \alpha^2 \beta$ is introduced as a new weight for pathway II labeled as pathway II'. }
	\label{fig:pathways}
\end{figure}

A model based on the interference between two quantum pathways is used to explain the experimental observations in the case of tungsten \cite{foerster2016}. The model considers the pathways I and II, illustrated in figure \ref{fig:pathways}, where two $(\omega)$-photons (red) are replaced by one $(2\omega)$-photon (blue) in pathway II. The weight $\alpha$ is related to the n-photon emission law $\alpha^n I_{\omega}^n$ and the weight $\beta$ is defined analogously for the second harmonic. The resulting (oscillating) transition probability is calculated by multiplying all partial amplitudes represented by the arrows in figure \ref{fig:pathways} within each path, adding the considered pathways and calculating the squared modulus to shift from amplitudes to rates.

Thus, the two-channel model predicts an oscillating term $(\alpha^{n-4} I_{\omega}^{n-4})$  $\left( 2 \alpha^3 I_{\omega}^3 \sqrt{\beta} \sqrt{I_{2\omega}} \right) \cos(2\omega\tau)$, whose prefactor is identified with $B_{2\omega}$. The resulting term $\alpha^n I_{\omega}^n$ corresponds to $A_{\textrm{DC}}$, where the small rate added by the second harmonic field is neglected in this model. The third term $(\alpha^{n-4} I_{\omega}^{n-4}) \times  \alpha^2 I_{\omega}^2 \beta I_{2\omega}$ is the incoherent background $B_{\textrm{DC}}$, which scales linear with $I_{2\omega}$. The two-pathway model cannot explain a frequency component at $4\omega/(2\pi)$ and expects $B_{2\omega}$ to follow an exact square-root dependence.

High second harmonic admixtures make higher order substitutions likely. We now include pathway III, where four $(\omega)$-photons are replaced by two $(2\omega)$-photons. By evaluating all interference terms with the amplitudes given in figure \ref{fig:pathways} we obtain a component oscillating at $4\omega/(2\pi)$, which we identify as  $B_{4\omega}$, and the scaling laws
\begin{eqnarray}
B_\textrm{DC} &=  (\alpha^{n-4} I_{\omega}^{n-4}) \times \alpha^2 I_{\omega}^2 \beta I_{2\omega} = F_\textrm{DC} I_{2\omega} \label{eq:F_DC} \\
B_{2\omega} &= (\alpha^{n-4} I_{\omega}^{n-4})  \nonumber \\
&\times \left( 2 \alpha^3 I_{\omega}^3 \sqrt{\beta} \sqrt{I_{2\omega}} + 2 \alpha I_{\omega} \sqrt{\beta}^3 \sqrt{I_{2\omega}}^3 \right) \\
&= F_{2\omega,1} \sqrt{I_{2\omega}} + F_{2\omega,2} \sqrt{I_{2\omega}}^3   \\
B_{4\omega} &= (\alpha^{n-4} I_{\omega}^{n-4}) \times 2 \alpha^2 I_{\omega}^2 \beta I_{2\omega} = F_{4\omega} I_{2\omega}. \label{eq:F_4}
\end{eqnarray}

The three-channel model can explain the linear scalings of $B_\textrm{DC}$ and $B_{4\omega}$ as well as an increased order for $B_{2\omega}$. The corresponding fits match the data in figure \ref{fig:power}(e) and evaluating the fit coefficients $F_i$ results in

\begin{eqnarray*}
 F_\textrm{DC} &= (1.70 \pm 0.02) \times 10^{-6} \,\textnormal{Hz} \,\textnormal{cm}^2 \,\textnormal{W}^{-1} \\
 F_{2\omega,1} &= (4.47 \pm 0.06) \times 10^{-1} \,\textnormal{Hz} \,\textnormal{cm} \,\textnormal{W}^{-1/2} \\
 F_{2\omega,2} &= (2.7 \pm  0.2)\times 10^{-12} \,\textnormal{Hz} \,\textnormal{cm}^3\, \textnormal{W}^{-3/2} \\
 F_{4\omega} &=  (3.40 \pm 0.07) \times 10^{-7} \,\textnormal{Hz} \,\textnormal{cm}^2 \,\textnormal{W}^{-1}.
\end{eqnarray*}

These four coefficients, however, cannot be retrieved by the simple relations for $\alpha$ and $\beta$. This can be seen by comparing $F_\textrm{DC}$ with $F_{4\omega}$, which do not fulfill $2F_\textrm{DC} = F_{4\omega}$ as expected from equations (\ref{eq:F_DC}) and (\ref{eq:F_4}). Allowing a new weight $\zeta$ for the exchange of one $(2\omega)$-photon with two $(\omega)$-photons (pathway II') we get

\begin{eqnarray}
B_\textrm{DC} &=  (\alpha^{n-4} I_{\omega}^{n-4}) \times \zeta^2 I_{\omega}^2 I_{2\omega} \\
B_{2\omega} &= (\alpha^{n-4} I_{\omega}^{n-4}) \nonumber  \\
&\times \left( 2 \zeta I_{\omega} \alpha^2 I_{\omega}^2  \sqrt{I_{2\omega}} + 2 \zeta I_{\omega} \beta \sqrt{I_{2\omega}}^3 \right) \label{eq:comp}\\
B_{4\omega} &= (\alpha^{n-4} I_{\omega}^{n-4}) \times 2 \alpha^2 I_{\omega}^2 \beta I_{2\omega}.
\end{eqnarray}
Using $F_\textrm{DC}$, $F_{2\omega,1}$ and $F_{2\omega,2}$ we can evaluate

\begin{figure*}[h]
	\centering
	\includegraphics[width=1\linewidth]{Figures/TDSE4a_4d}
	\caption{TDSE simulation of the cooperative signal. (a) Time domain representation of simulated photocurrent for pulse durations of $\tau_\omega = 25 \textrm{ fs}$ and $\tau_{2\omega} = 22 \textrm{ fs}$ and intensities of $I_{\omega} = 900 \,\textrm{GW}/\textrm{cm}^2$ and $I_{2\omega}= 4.5 \,\textrm{GW}/\textrm{cm}^2$.
		Left inset: Fourier amplitude and regions of interest $\textrm{ROI}_\textrm{DC}$, $\textrm{ROI}_{2\omega}$ and $\textrm{ROI}_{4\omega}$. Right inset: peak heights of Gaussian fits for the corresponding back transformed regions as function of the second harmonic intensity. The double-logarithmic representation clearly shows a change in the power law exponent for the peak height in $\textrm{ROI}_{2\omega}$. (b) Visibility together with peak heights $B_\textrm{DC}$, $B_{2\omega}$ and $B_{4\omega}$ (same as right inset of (a)) in linear scale as function of the second harmonic intensity and intensity admixture. Here, all data points are used to estimate the average power law exponent for $B_{2\omega}$. (c) Retrieved time domain signal (light gray) and original time domain signal form (a) by stepwise adding the Fourier components to the offset $A_\textrm{DC}$. For clarity only the envelopes are indicated. Small deviations are visible as the values for $B_\textrm{DC}$, $B_{2\omega}$ and $B_{4\omega}$ are taken from the fit curves in (b). (d) Fourier components with peak heights $\tilde{B}_\textrm{DC}$, $\tilde{B}_{2\omega}$ and $\tilde{B}_{4\omega}$ showing the corresponding frequency domain structure. The blue curve indicates the simulation.}
	\label{fig:tdse}
\end{figure*}

\begin{eqnarray*}
\alpha  &=  (8.70 \pm 0.05) \times 10^{-12} \,\textnormal{cm}^2 \,\textnormal{W}^{-1} \,\textnormal{Hz}^{-5} \\
\beta  &=  (  3.7 \pm  0.3) \times 10^{-10}  \,\textnormal{cm}^2 \,\textnormal{W}^{-1} \,\textnormal{Hz}^{-3}  \\
\zeta &=  (5.18 \pm    0.04) \times 10^{-16} \,\textnormal{cm}^3 \,\textnormal{W}^{-3/2} \,\textnormal{Hz}^{-13/2}.
\end{eqnarray*}

$F_\textrm{DC}$, $F_{2\omega}$ and $F_{4\omega}$ are by definition exactly reproduced, but also the independently retrieved value $F_{4\omega,r} = 2 \alpha^2 I_{\omega}^2 \beta =  (3.5 \pm  0.3) \times 10^{-7} \,\textnormal{Hz} \,\textnormal{cm}^2 \textnormal{W}^{-1}$ matches. The retrieved offset $A_\textrm{DC,r} \approx 2.9 \times 10^4$ Hz is also close to $A_\textrm{DC} \approx 3.5 \times 10^4$ Hz in the experiment. The maximum visibility in figure \ref{fig:power}(e) can be estimated using
\begin{eqnarray}
 V = \frac{B_{2\omega}}{A_\textrm{DC}+B_\textrm{DC}+B_{4\omega}}
\end{eqnarray}
and reproduces the experiment well, which is depicted by the solid black curve in figure \ref{fig:power}(e). We conclude that the third quantum pathway is needed in the presence of a frequency component around $(4\omega)/(2\pi)$ and allows us to accurately describe the experimental data.

\section{TDSE-simulations}
In this section we compare our experimental findings with model simulations based on the 1-dimensional time-dependent Schrödinger equation (TDSE), which is explained in detail in \cite{seiffert2018}. We use the same near-field intensities as in the experiment and assume a clean gold surface with work function $W=5.4$ eV and Fermi energy $E_\textnormal{F} = 5.2$ eV \cite{wu2020} without an applied bias voltage. In the experiment the effective barrier height is lower, as Schottky-lowering and possible residuals at the gold surface can reduce the work function \cite{kahn2016}. Further, we had to decrease the pulse durations to $\tau_\omega = 25 \textrm{ fs}$ and $\tau_{2\omega} = 22 \textrm{ fs}$ to keep the numerical effort manageable for a delay step size of 0.26 fs. Nevertheless, both pulses are clearly within the multi-cycle regime.

A typical simulated time trace is shown in figure \ref{fig:tdse}(a) and resembles the experiment closely. For a qualitative comparison we chose a trace with similar visibility and scaled the peak count rate to match the one from figure \ref{fig:power}(a). The Fourier amplitudes in the left inset of figure \ref{fig:tdse}(a) contain again the three characteristic components and show no additional features. The peak heights show a clear trend in the power law exponent of $B_{2\omega}$ exceeding the square root scaling for high second harmonic intensities as expected from the three-pathway model. $B_{2\omega}$ and $B_\textnormal{DC}$ maintain a linear scaling. Figure \ref{fig:tdse}(b) shows the same trends as in the experiment, even the average order of $B_{2\omega}$ coincides. The peak visibility in the simulation is reached for lower second harmonic admixtures, which is caused by the different effective barrier heights and uncertainties in the exact near field strengths in the experiment. The fit values

\begin{eqnarray*}
	F_\textrm{DC} &= (  1.18 \pm 0.01) \times 10^{-5} \,\textnormal{Hz} \,\textnormal{cm}^2 \,\textnormal{W}^{-1} \\
	F_{2\omega,1} &= (1.264 \pm 0.002) \,\textnormal{Hz} \,\textnormal{cm} \,\textnormal{W}^{-1/2} \\
	F_{2\omega,2} &= (8.08 \pm  0.03)\times 10^{-11} \,\textnormal{Hz} \,\textnormal{cm}^3 \,\textnormal{W}^{-3/2} \\
	F_{4\omega} &=  (4.72 \pm  0.03) \times 10^{-6} \,\textnormal{Hz} \,\textnormal{cm}^2 \,\textnormal{W}^{-1}
\end{eqnarray*}
allow again to estimate

\begin{eqnarray*}
	\alpha &= (  4.931 \pm 0.007) \times 10^{-12}  \,\textnormal{cm}^2 \,\textnormal{W}^{-1}\, \textnormal{Hz}^{-7} \\
	\beta &= (1.26  \pm  0.02) \times 10^{-9} \,\textnormal{cm}^2\, \textnormal{W}^{-1} \,\textnormal{Hz}^{-4} \\
	\zeta &= ( 4.08 \pm  0.02)\times 10^{-16} \,\textnormal{cm}^3 \,\textnormal{W}^{-3/2} \,\textnormal{Hz}^{-9}.
\end{eqnarray*}

The obtained offset $A_{\textrm{DC}} = 4.25 \times 10^{4} \,\textnormal{Hz}$ slightly deviates from the retrieved value of $A_{\textrm{DC},r} =  3.4 \times 10^{4} \,\textnormal{Hz}$ as well as $F_{4\omega,r} =  (4.33 \pm  0.06) \times 10^{-6} \,\textnormal{Hz} \,\textnormal{cm}^2 \,\textnormal{W}^{-1}$. Furthermore, we can compare the values of $\alpha$ and $\beta$ against those estimated from an intensity-dependent yield-sweep providing $\alpha_p =   4.687  \times 10^{-12} \, \textnormal{cm}^2 \,\textnormal{W}^{-1} \,\textnormal{Hz}^{-7}$ and $\beta_p =  3.31 \times 10^{-11} \,\textnormal{cm}^2 \,\textnormal{W}^{-1} \,\textnormal{Hz}^{-4}$. The strong deviation in $\beta$ suggests that although all scalings are correctly predicted by the quantum-path model, the prefactors within the temporal overlap cannot be simply constructed from those of the individual fields.

In figure \ref{fig:tdse}(c) we investigate the formation of the time domain signal by adding up the corresponding regions of interest in the frequency domain in figure \ref{fig:tdse}(d) including their estimated widths. The main part of the signal is constructed by the offset $A_{\textrm{DC}}$ and $B_\textnormal{DC}$ defining the base line around which the (symmetric) oscillations with amplitude $B_{2\omega}$ appear. $B_{4\omega}$ flattens the trace for count rates close to zero and has to be included to prevent negative count rates. Adding up all these contributions gives nearly full agreement between retrieved and analyzed trace, even though $B_\textnormal{DC}$, $B_{2\omega}$ and $B_{4\omega}$ are obtained from the fit values for $\alpha$, $\beta$ and $\zeta$. In the Fourier domain, see figure \ref{fig:tdse}(d), this agreement is intuitively explained as the Fourier components are simply replaced by their fits and no complicated Fourier phase is present.

\section{Influence of the tip bias voltage}

\begin{figure*} [h]
	\centering
	\includegraphics[width=1\linewidth]{Figures/Experiment_bias_sweep5a_5b}
	\caption{ Influence of the tip bias voltage. (a) Visibility as function of second harmonic intensity for bias voltages $\textrm{U}_\textrm{DC}$ from -100 V to -1300 V in -200 V steps at fundamental intensity $I_{\omega} = 540 \, \textrm{ GW}/\textrm{cm}^2$. Inset: visibility at the highest second harmonic intensity ($\approx$ maximal visibility). (b) Power law exponent for the DC- and $2\omega$ component as function of the bias voltage.}

	\label{fig:bias}
\end{figure*}

Gold tips with moderate apex radii ($>20$ nm) have high field emission thresholds ($<-1300$ V) while preserving sufficient photocurrents due to strong optical near-fields. Therefore, the big range of possible bias voltages is another important tuning knob in the coherent control scheme for gold. In figure \ref{fig:bias}(a) we investigate the influence of the bias voltage on the maximum achievable visibility and its dependence on the second harmonic intensity. Although all voltages are below the field emission threshold, a clear reduction of the maximum visibility is observed and the high visibility domain requires higher second harmonic intensities.
The reduced visibility is accompanied by a shift in the power law exponent for the $2\omega$ component as shown in figure \ref{fig:bias}(b). The effective power law exponent is connected to the ratio of the components $2 \zeta I_{\omega} \alpha^2 I_{\omega}^2$ and $ 2 \zeta I_{\omega} \beta I_{2\omega}$, which scale differently with $I_{2\omega}$ (equation (\ref{eq:comp})). Their ratio is given by $(\alpha^2 I_{\omega}^2)/ \beta I_{2\omega}$ and shows that higher negative voltages further increase $\alpha I_\omega$ with respect to $\beta I_{2\omega}$. Thus, also the modifications of the visibility are most likely caused by changes in the weights of the respective quantum pathways.

 The bias voltage modifies the shape and height of the potential barrier simultaneously and therefore affects not only the values of $\alpha$,  $\beta$ and $\zeta$, but also their corresponding multi-photon orders. This interplay makes a model for the cooperative signal including bias voltages a challenging task, which is beyond the scope of this work.

\section{Conclusion and outlook}
In this manuscript we have shown that the coherent control scheme can be applied to gold needle tips and, as in the case of tungsten, a near 100$\%$ modulation of the photocurrent it possible. The quantum pathway model was extended by a third emission channel to describe the experiment and is supported by TDSE simulations. Applying bias voltages can modify the weights of the individual emission channels and leads to a reduced maximum visibility.

The requirement of independent weights for the three channels and the discrepancy in the extracted prefactors from those of independent pulses suggest modifications to the simple quantum pathway model. Analyzing the quantum pathway model using analytical solutions of the TDSE including bias voltages \cite{luo2019,luo2021} are therefore of highest interest for coherent control schemes in general and photoemission from metals in two-color fields in particular.

Using shorter and stronger fundamental pulses should in addition allow to observe trajectory modifications \cite{seiffert2018} driven by the two-color field and the spatial emission profile \cite{yanagisawa2009} could be modified by polarization shaped pulses. Thus, the combination of gold nanostructures with tailored emission in two-color fields may transfer electron wavepacket shaping capabilities previously reserved for atomic gases to nano-plasmonics, on-chip lightwave electronics and electron sources.

\section*{Acknowledgments}
This project was funded in part by the ERC grant ``Near Field Atto''and DFG SPP 1840 ``QUTIF''.

\section*{References}

\section{Introduction}

The \texttt{iopart-num} Bib\TeX{} style is intended for use in
preparing manuscripts for Institute of Physics Publishing journals,
including Journal of Physics.  It provides numeric citation with
Harvard-like formatting, based upon the specification in ``How to
prepare and submit an article for publication in an IOP journal using
\LaTeXe'' by Graham Douglas (2005).

The \texttt{iopart-num} package is available on the Comprehensive
\TeX{} Archive Network (CTAN) as \texttt{/biblio/bibtex/contrib/iopart-num}.

\section{General instructions}

To use the \texttt{iopart-num} style, include the command
\verb+\bibliographystyle{iopart-num}+ in the
document preamble.  The reference section is then inserted into the
document with the command \verb+


Effect of water-wall interaction potential on the properties of nanoconfined water

Abstract: Much of the understanding of bulk liquids has progressed through study of the limiting case in which molecules interact via purely repulsive forces, such as a hard-core potential. In the same spirit, we report progress on the understanding of confined water by examining the behavior of water-like molecules interacting with planar walls via purely repulsive forces and compare our results with those obtained for Lennard-Jones (LJ) interactions between the molecules and the walls. Specifically, we perform molecular dynamics simulations of 512 water-like molecules which are confined between two smooth planar walls that are separated by 1.1 nm. At this separation, there are either two or three molecular layers of water, depending on density. We study two different forms of repulsive confinements, when the interaction potential between water-wall is (i) $1/r^9$ and (ii) WCA-like repulsive potential. We find that the thermodynamic, dynamic and structural properties of the liquid in purely repulsive confinements qualitatively match those for a system with a pure LJ attraction to the wall. In previous studies that include attractions, freezing into monolayer or trilayer ice was seen for this wall separation. Using the same separation as these previous studies, we find that the crystal state is not stable with $1/r^9$ repulsive walls but is stable with WCA-like repulsive confinement. However, by carefully adjusting the separation of the plates with $1/r^9$ repulsive interactions so that the effective space available to the molecules is the same as that for LJ confinement, we find that the same crystal phases are stable. This result emphasizes the importance of comparing systems only using the same effective confinement, which may differ from the geometric separation of the confining surfaces.

\section{Introduction}

Confinement of water in nanopores affects many properties of water,
such as freezing temperature, crystal
structure~\cite{bellisent1,koga1,koga2,koga2005,zangi-rev,zangi1}, the
glass transition temperature, and the position of the hypothesized
liquid-liquid (LL) critical
point~\cite{zanotti1999,ChenPRVT,bellisent,xu2005pnas,ku2005sub}.
Indeed, water confined in nanoscale geometries has received much recent
attention, in part because of its importance in biology, engineering,
geophysics and atmospheric sciences.  The effects of different kinds
of confinement have been studied, both using experiments and
simulations~\cite{antognozzi,gilijamse,seemasingh,koga1,koga2,zangi1,zangi2,kumar2005,nicolas,matri,gallo2002}.

Bulk supercooled water -- water cooled below the equilibrium freezing
temperature -- shows many anomalous
properties~\cite{Debenedetti03a,Debenedetti03b,angell-review,bellisent1}.
Experiments find that at low temperatures, various response functions,
such as isothermal compressibility and specific heat, increase
sharply. There has been comparatively less research on confined water.
Using the ST2 potential to model water confined between smooth plates
\cite{meyer}, a LL phase transition has been proposed.  A
liquid-to-amorphous transition is seen in simulations of water using
the TIP4P potential \cite{tip4p} confined in carbon nanotubes
\cite{koga2}.  Recent theoretical work \cite{truskett} suggests that
hydrophobic Lennard-Jones (LJ) confinement shifts the LL transition to
lower temperature and lower pressure compared to bulk water, a feature
also found in simulations of water confined between hydrophobic
plates~\cite{kumar2005}.

Confinement is known to enhance solidification of molecules that are
more or less spherical~\cite{rhykerd,brushan,persson}. However,
careful experiments on thin films of water show that water performs
extremely well as a lubricant, suggesting that confined water may be
more fluid than bulk water~\cite{fluidity}.  Recent experiments show
that water in hydrophilic confinement, when cooled to very low $T$,
does not freeze~\cite{bellisent} -- a phenomenon also supported by
simulation studies~\cite{Gallo00,Gallo00b}.  In contrast,
simulations~\cite{koga1,koga2,zangi1,zangi2,kumar2005,nicolas} show
that hydrophobically confined water does freeze into different
crystalline structures, which do not have counterparts in bulk water.
Indeed monolayer, bilayer and trilayer ice have all been found in
simulations~\cite{koga1,koga2,zangi1,zangi2,kumar2005,nicolas}.  Thus
hydrophobic confinement seems to facilitate the freezing of
water. However the reason for this facilitation is not yet fully
understood. The hydrogen-bond interaction between water molecules is
an order of magnitude stronger than Van der Waals attraction with the
hydrophobic walls. Thus one may hypothesize that freezing in
hydrophobic confinement depends critically on the separation between
confining walls that may distort or facilitate a particular
crystalline structure, rather than on the weak details of the
water-wall interaction potential. To test this hypothesis we perform
molecular dynamics (MD) simulations of water in two different forms of
repusive confinement. Specifically, we study:

\begin{itemize}

\item The $1/r^9$ repulsive part of the LJ potential studied in
Ref.~\cite{kumar2005}.

\item The same potential used in Ref.~\cite{kumar2005} but truncated
and shifted such that there is no attractive part in the potential,
analogous to the Weeks-Chandler-Anderson (WCA)
potential~\cite{weeks}. This potential allows us to examine the role,
if any, the attractive part of the water-wall LJ potential plays in
determining the thermodynamics and structure of confined water.

\end{itemize}

We compare the case when the water-wall interactions are purely
repulsive (``repulsive confinement'') with the studied case of
pure LJ confinement~\cite{kumar2005}. We also compare the freezing in
repulsive confinements with the freezing found when the water-wall
interactions are represented by an LJ interaction (``LJ hydrophobic
confinement'')~\cite{kumar2005,koga1}.

This paper is organized as follows: In Sec.~II, we provide details of
our simulations and analysis methods. Simulation results for the
liquid state are provided in Sec.~III. In Sec.~IV we discuss the
freezing properties of our system.

\section{Simulation and Analysis Methods}

We perform MD simulations of a system composed of
water-like molecules confined between two smooth walls.  The molecules
interact via the TIP5P pair potential \cite{jorgensen1} which, like
the ST2~\cite{st2} potential, treats each water molecule as a
tetrahedral, rigid, and non-polarizable unit consisting of five point
sites~\cite{note1}. The TIP5P potential predicts many of the anomalies
of bulk water \cite{masako}. For example, TIP5P reproduces the density
anomaly at $T=277$~K and $P=1$~atm and its structural properties
compare well with experiments
\cite{masako,jorgensen1,jorgensen2,sorenson,kuOdessa}.  TIP5P is known
to crystallize at high pressures~\cite{masako} within accessible
computer simulation time scales, and shows a ``nose-shaped'' curve of
temperature versus crystallization time~\cite{masako}, a feature found
in experimental data on water solutions \cite{baez}.

In our simulations, $N=512$ water molecules are confined between two
 smooth planar walls, as shown schematically in
 Fig.~\ref{fig:schematic}.  The walls are located at $z_w=\pm0.55$~nm
 (wall-wall separation of $1.1$~nm), which results in $\approx 2-3$
 layers of water molecules.  Periodic boundary conditions are used in
 the $x$ and $y$ directions, parallel to the walls.

We study two different forms of purely repulsive water-wall
interaction. The first uses only the $r^{-9}$ repulsive core, which we
call the $1/r^9$ repulsive potential,
\begin{equation}
U\left(z-z_{W}\right) =
4\epsilon_{\rm OW}\left[\left(\frac{\sigma_{\rm OW}}{|z-z_W|}\right)^9\right].
\label{Udz}
\end{equation}
Here $|z-z_W|$ is the distance from the oxygen atom of a water
molecule to the wall, while $\epsilon_{\rm OW} = 0.25$ kJ/mol and
$\sigma_{\rm OW}= 0.25$ nm are potential parameters
(Fig.~\ref{fig:potentials}). Similar but different parameter values
were used in previous confined water simulations using the TIP5P
interaction potential~\cite{kumar2005}. Specifically, in Ref.~\cite
{kumar2005}, the water-wall interaction was modelled using a 9-3
LJ potential with $\epsilon_{\rm OW}=1.25$~kJ/mol and
$\sigma_{OW}=0.25 $~nm. We choose a different $\epsilon_{\rm OW}$ in
the case of repulsive confinement so that the repulsion between the
water and wall decays to almost zero where the 9-3 LJ -potential has a
minimum.

The second purely repulsive potential uses both attractive and
repulsive terms of the 9-3 LJ potential, but truncates and shifts the
potential at the position of the minimum to create a repulsive
potential that exactly mimics the repulsion of ref.~\cite{kumar2005},
in analogy to the WCA potential,

\begin{displaymath}
U(z-z_W) = \left\{ \begin{array}{ll}
    4\epsilon_{\rm OW} [ (\frac{\sigma_{\rm OW}}{z-z_W})^9-(\frac{\sigma_{\rm OW}}{z-z_W})^3 ] + \frac{8\epsilon_{\rm OW}}{3^{3/2}} & \textrm{if $|z-z_W| < 3^{1/6}\sigma_{\rm OW}$}\\
 0 & \textrm{if $|z-z_W| > 3^{1/6}\sigma_{\rm OW}$},\\
\end{array} \right.
\end{displaymath}

where $\epsilon_{\rm OW}=1.25$~kJ/mol and $\sigma_{OW}=0.25 $~nm. For
the $1/r^9$ repulsive potential, we perform simulations for $56$ state
points, corresponding to seven temperatures $T=220$~K, $230$~K,
$240$~K, $250$~K, $260$~K, $280$~K, and $300$~K, and eight ``geometric
densities'' $\rho_g=0.60$~g/cm$^3$, $0.655$~g/cm$^3$,
$0.709$~g/cm$^3$, $0.764$~g/cm$^3$, $0.818$~g/cm$^3$,
$0.873$~g/cm$^3$, $0.927$~g/cm$^3$, and $0.981$~g/cm$^3$ -- the same
as studied in Ref.~\cite{kumar2005}.  The geometric values of density
do not take into account the fact that the repulsive interactions of
molecules with the walls increases the overall amount of available
space, since the $\epsilon_{OW}$ parameter of the $1/r^9$ repulsive
potential is smaller than the $\epsilon_{OW}$ used for the LJ confined
system~\ref{fig:potentials}. For systems confined by LJ interactions,
there is a well-defined preferred distance from the wall, making it
relatively straightforward to evaluate the ``effective'' density of
molecules confined by the attractive wall.  In our system with only
repulsive interactions, there is no such preferred distance, as
emphasized by Fig.~\ref{fig:potentials}.

We can approximate the effective density by examining the local
density $\rho(z)$ (Fig.~\ref{fig:rhoz}). We utilize the fact that
$\rho(z)$ has an inflection, and estimate the effective $L_{z}$ by the
location where the second derivative of $\rho(z) = 0$, or when first
derivative of $\rho(z)$ has a maximum.  We must also add to this value
of $L_{z}$ the molecular diameter of water(0.278~nm) to calculate the
real space available along $z$-direction.  The resulting ``effective
densities'' are $\rho=0.715$~g/cm$^3$, $0.777$~g/cm$^3$,
$0.829$~g/cm$^3$, $0.890$~g/cm$^3$, $0.949$~g/cm$^3$,
$1.000$~g/cm$^3$, $1.060$~g/cm$^3$, and $1.115$~g/cm$^3$.  We will use
these effective densities throughout the paper, since they will be
most comparable to the effective densities with LJ confinement.

For the WCA potential, we perform simulations for 32 state points,
corresponding to 8 different temperatures and 4 different ``geometric
densities'', $0.60$~g/cm$^3$,
$0.655$~g/cm$^3$,$0.709$~g/cm$^3$,$0.764$~g/cm$^3$ respectively. These
geometric densities correspond to ``effective densities''
$0.80$~g/cm$^3$, $0.88$g/cm$^3$, $0.95$g/cm$^3$, $1.02$g/cm$^3$
respectively. Note that these effective densities were calculated
using the method described in ~\cite{meyer,kumar2005}.

We control the temperature using the Berendsen thermostat with a time
constant of $5$ ps \cite{berend} and use a simulation time step of
$1$~fs, just as in the bulk system \cite{masako}. Water-water
interactions are truncated at a distance~$0.9$~nm as discussed in
Ref.~\cite{jorgensen1}.

\section{Thermodynamics and Structure}

One of the defining characteristics of water is the existence of a
temperature of maximum density (TMD). Relative to bulk water -- LJ
confinement shifts the locus of the TMD to lower $T$ by $\approx
40$~K~\cite{kumar2005}. Additionally, the sharpness of the density
maximum is markedly decreased in comparison to the bulk.
Fig.~\ref{fig:pxy-comp} shows isochores of $P$ for LJ confinement,
$1/r^9$ repulsive confinement and WCA confinement, for similar
densities. A TMD in this plot is coincident with the minimum in the
isochore.  For $1/r^9$ repulsive confinement, the minimum is very
weak, but the location of the flatness in the isochore is near to that
of the system with LJ confinement.  This result suggests the the
$1/r^9$ repulsive confinement further suppresses the structural
ordering of the molecules that is known to be responsible to the
presence of a density maximum. The TMD for the case of WCA confinement
again appears at the same $T$ as the $1/r^9$ repulsive confinement and
LJ confinement cases but the isochore in the TMD region is flatter
than for case the of LJ confinement confinemen. Hence both kinds of
the repulsive confinement supprsess the structural ordering in lateral
directions compared to the case of bulk and LJ confinement. We further
notice that the value of the lateral pressure $P_{\|}$~in the case of
the LJ confinement approaches the value of $P_{\|}$ in case of WCA
confinement at high temperatures. This behavior of $P_{\|}$ for
LJ confinement should be expected since at very high temperatures the
molecules will not feel the potential minimum of the water-wall
interaction.

In order to compare the structural properties of repulsive
confinement with those of LJ confinement, we calculate the lateral
oxygen-oxygen radial distribution function (RDF) defined by
\begin{equation}
g_{\|}(r) \equiv \frac{1}{\rho^2 V}\sum_{i \ne j} \delta(r-r_{ij})
\left[\theta \left(|z_i-z_j|+\frac{\delta z}{2}\right) -\theta
\left(|z_i-z_j|-\frac{\delta z}{2}\right)\right].
\end{equation}
Here $V$ is the volume, $r_{ij}$ is the distance parallel to the walls
between molecules $i$ and $j$, $z_i$ is the $z$-coordinate of the oxygen
atom of molecule $i$, and $\delta(x)$ is the Dirac delta function. The
Heaviside functions, $\theta(x)$, restrict the sum to a pair of oxygen
atoms of molecules located in the same slab of thickness $\delta z =
0.1$~nm.  The physical interpretation of $g_{\|}(r)$ is that $g_{\|}(r)
2 \pi r dr \delta z$ is proportional to the probability of finding an
oxygen atom in a slab of thickness $\delta z$ at a distance $r$ (parallel
to the walls) from a randomly chosen oxygen atom. In a bulk liquid, this
would be identical to $g(r)$, the standard RDF.

Figure~\ref{grxy} shows the temperature and density dependence of the
lateral oxygen-oxygen pair correlation function for both $1/r^9$
repulsive (Fig.~\ref{grxy}a, Fig.~\ref{grxy}b) and WCA
(Fig.~\ref{grxy}c, Fig.~\ref{grxy}d) confinements. For both repulsive
confinements, the qualitative behavior of the dependence of $g_\|$ is
the same. At low temperature and low density, the first two peaks in
$g_{\|}$ appear at $r=2.78$~\AA and $r=4.5$~\AA, but at high densities
the second peak moves to a larger distance.  This behavior is nearly
identical to that observed for water confined between LJ surfaces, and
is discussed in detail in Ref.~\cite{kumar2005}.

We also confirm the structural similarity with LJ confinement by
calculating the lateral static structure factor $S_{\|}(q)$, defined
as the Fourier transform of the lateral RDF
$g_{\|}(r)$,
\begin{equation}
S_{\|}\left(q\right) \equiv \frac{1}{N}\sum_{j,k}\left<e^{i\vec{\bf
q}.(\vec{\bf r_j} - {\bf \vec{r_k}})}\right>.
\end{equation}
Here the $q$-vector is the corresponding wave vector in the $xy$ plane
and $r$ is the projection of the position vector on the $xy$ plane. In
Fig.~\ref{fig:sqxy}, we show the temperature and pressure dependence
of lateral structure factors for both repulsive confinements. For both
forms of repulsive confinement, the temperature and density dependence
of $S_{\|}$ is similar. We find that confined water has a weaker
pre-peak at $\approx 18$nm$^{-1}$ compared to bulk water
(Fig.~\ref{fig:sqxy-comp}), consistent with the possibility that the
local tetrahedrality is weakened by repulsive confinements. Of the
three forms of confinement, the $S_{\|}$ for LJ confinement is most
like bulk water (Fig.~\ref{fig:sqxy-comp}). Local tetrahedrality
becomes weaker in case of repulisive confinements compared to LJ
confinement. Further, we see that the water in $1/r^9$ repulsive
confinement is less structured in lateral directions compared to water
in WCA confinement, indicated by a less sharp and broad pre-peak at
$\approx 18$~nm$^{-1}$~ in $S_{\|}(q)$~(Fig.~\ref{fig:sqxy-comp}).

\section{Freezing of TIP5P water}
\label{sec:freezing}

Bulk TIP5P water crystallizes within the simulation time for $\rho
\gtrsim 1.15$~g/cm$^3$ at low temperatures~\cite{masako}. Crystallization
of confined water is seen in some simulations
\cite{zangi1,koga1,kumar2005}.  A similar crystallization appears in
simulations when an electric field is applied in lateral directions to
a system of water confined between silica walls~\cite{zangi2}.

At plate separation of 1.1~nm with hydrophobic LJ confinement, water
crystallizes to trilayer ice~\cite{kumar2005}.  From our simulations
of TIP5P water in repulsive confinements with the same plate
separation of 1.1~nm, we find that the system does not freeze within
accessible simulation time scales for $1/r^9$ repulsive confinement;
however the system freezes for WCA confinement.  As a more stringent
confirmation of this fact, we also use a starting ice configuration
obtained from simulations with LJ confinement for the same thickness,
and confirm that the ice melts to a liquid with $1/r^9$ repulsive
confinement. In Fig.~\ref{fig:melting-repulsive}, we show the
evolution of potential energy and lateral structure factor with time,
when the crystal formed in LJ confinement~\cite{kumar2005} is kept
between the $1/r^9$ repulsive walls. The potential energy first
increases and then reaches its equilibrium value of the liquid
potential energy accompanied by a structural change from a crystal
(presence of sharp Bragg peaks) to a liquid (absence of Bragg peaks).

Based on this observation, it is tempting to claim that repulsion
inhibits crystallization, and that a preferable distance from the wall
determined by the attractive portion of the LJ potential is necessary
to induce crystallization.  However, as discussed above for the same
plate separation $1.1$~nm, $1/r^9$ repulsive confinement with the
chosen parameters increases the available space for molecules relative
to LJ confinement.  Hence to properly compare the crystallization
behavior, we must adjust the separation of the wall so that the
available space for the water molecules is the same in both systems.
We can make the available space the same by tuning the separation of
the plates or by tuning the potential. By tuning the parameters (see
Fig.~\ref{fig:potentials}, where the values of parameter
$\epsilon_{\rm OW}$ and $\sigma_{\rm OW}$ for the modified
$1/r^9$-repulsive potential are $1.25$~kJ/mol and $0.23$~nm
respectively) of the $1/r^9$repulsive potential in the repulsive
system such that the density profile along the $z$-axis becomes
similar, we have identical values of the available space between the
plates (see Fig.~\ref{fig:rhoz-comp}), and we find that an initial
crystal configuration does {\it not} melt, emphasizing that the
presence of the crystal is very sensitive to density {\it and} to
plate separation -- since the separation determines the accessible
packing arrangements between the plates.  Similar sensitivity to plate
separation for monolayer ice was seen in ref.~\cite{zangi2}.

In addition to examining the stability of initially crystalline
structures, we also consider whether freezing from the liquid state
occurs when we have the same effective plate separation.  We find that
for the repulsively confined systems, the crystal will also
spontaneously form if the available space between the plates is the
same as that for which initially crystalline configurations are stable.
Hence plate separation appears to be the dominant cause in
determining whether or not a crystal will form.

\section{Acknowledgments}

We thank F. Sciortino for helpful discussions, and the NSF for
support under grants CHE~0096892, CHE~0404699, and DMR-0427239. We
also thank the Boston University computational facility for
computational time and the Office of Academic Affairs at Yeshiva
University for support.