Abstract: Strong nonlinear interactions between single photons have important applications in optical quantum information processing. Demonstrations of these interactions in cold atomic ensembles have largely been limited to exploiting slow light generated using electromagnetically induced transparency (EIT). However, these EIT implementations have limited achievable phase shifts due to spontaneous emission. Here, we demonstrate and characterize a scheme free from these limitations using gradient echo memory with inferred single photon phase shifts of $0.07\pm0.02$ $\mu\text{rad}$. Excellent agreement with theoretical modelling was observed. Degradation of memory efficiency was observed for large phase shifts but strategies to overcome that are presented.

\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*{\refname}}
\bibliographystyle{iopart-num}

\begin{document}
\title{Observation of cross phase modulation in cold atom gradient echo memory}

\author{Anthony C. Leung}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}
\author{K. S. Ida Melody}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}
\author{Aaron D. Tranter}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}
\author{Karun V. Paul}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}
\author{Geoff T. Campbell}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}
\author{Ping Koy Lam}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}
\author{Ben C. Buchler}
    \email{Corresponding author; ben.buchler@anu.edu.au}
    \affiliation{ARC Centre of Excellence for Quantum Computation and Communication Technology}
    \affiliation{Department of Quantum Science and Technology, Research School of Physics, The Australian National University, Acton, ACT 2601, Australia}

\date{\today}

\begin{abstract}
Strong nonlinear interactions between single photons have important applications in optical quantum information processing.  Demonstrations of these interactions in cold atomic ensembles have largely been limited to exploiting slow light generated using electromagnetically induced transparency (EIT).  However, these EIT implementations have limited achievable phase shifts due to spontaneous emission.  Here, we demonstrate and characterize a scheme free from these limitations using gradient echo memory with inferred single photon phase shifts of $0.07\pm0.02$ $\mu$rad.  Excellent agreement with theoretical modelling was observed.  Degradation of memory efficiency was observed for large phase shifts but strategies to overcome that are presented.
\end{abstract}

\maketitle
\section{Introduction}
Facilitating strong interactions between photons is an important step towards realizing optical quantum computing \cite{Milburn_1989}.  This can be achieved via materials possessing a very strong nonlinearity causing a photon to impart a phase shift onto another through cross phase modulation (XPM).  This process can be used for many interesting applications such as universal quantum gates \cite{Milburn_1989}, quantum non-demolition measurements of photon numbers \cite{Imoto_Haus_Yamamoto_1985} and cluster state generation \cite{Kim_Lee_Ji_Nha_Anisimov_Dowling_2015}.

Various proposals have demonstrated examples of strong nonlinear interactions using different platforms.  Early promising systems took advantage of cavity quantum electrodynamics (cQED) using single atoms \cite{Turchette_Hood_Lange_Mabuchi_Kimble_1995} and quantum dots \cite{Fushman_Englund_Faraon_Stoltz_Petroff_Vuckovic_2008}.  Over the last few years, successful demonstrations of large phase shifts were observed using cQED applied to an atomic ensemble to get phase shifts of up to $\pi/3$ per postselected single photon \cite{Beck_Hosseini_Duan_Vuletic_2016}, Rydberg  blockades \cite{Tiarks_Schmidt_Rempe_Durr_2016} for $\pi$ per photon and hollow core fibres filled with atomic gas \cite{Venkataraman_Saha_Gaeta_2013} for $0.3$mrad per photon.

The main approach for generating XPM in atomic ensemble systems has been through optical dispersion engineering \cite{Chen_Wang_Wang_Yu_2006, Lo_Su_Chen_2010,Lo_Chen_Su_Chen_Chen_Chen_Yu_Chen_2011}.  The slow-light effect in EIT can be used to strengthen the nonlinear interaction in the medium as a.c.-Stark shift imparts a phase shift and creates XPM.  Phase shifts of up to 18 $\mu$rad per post-selected single photon have been observed \cite{Feizpour_Hallaji_Dmochowski_Steinberg_2015}.  It has been shown, however, that the magnitude of phase shift obtainable between two travelling photons via EIT is limited \cite{Gea-Banacloche_2010}.  This limit could only be overcome by operating the medium as a quantum memory to avoid simultaneously propagating optical fields or in the coherent photon conversion regime (a four-wave mixing effect that causes atomic population oscillations) which was recently demonstrated \cite{sagona-stophelConditionalPPhaseShift2020}.

In this paper, we measured the nonlinearity achieved by storing one light field using gradient echo memory (GEM) \cite{hosseiniHighEfficiencyCoherent2011, choHighlyEfficientOptical2016} and propagating the other field through the medium during the storage phase.  We infer 0.07 $\mu$rad of phase shift per single photon.  While this phase shift is fairly modest, we present prospects for increasing the size of the phase shift through changes to experimental design.

\section{Method}
\subsection{Preparation of ensemble}
A cloud of \textsuperscript{87}Rb atoms is captured and cooled using a magneto-optical trap (MOT).  A 2D quadrupole field is used to create an elongated geometry to maximize optical depth along the propagation path of the probe (optical depth of 450 on probe transition of $D_1$ $\ket{5S_{1/2},F{=}2} \rightarrow \ket{5P_{1/2}, F'{=}1}$).  The optical table is surrounded by three pairs of orthogonal Helmholtz coils to cancel out static background fields.  The experimental sequence was synchronized to the 50Hz mains power to ensure environmental consistency from run to run.   Details on the ensemble preparation can be found in \cite{choHighlyEfficientOptical2016, tranter2018}.

\subsection{Measuring XPM using GEM}
The experimental setup is shown in figure \ref{fig:setup}.  A weak probe was focused into the MOT along the elongated axis with a strong control field intersecting at an angle of $0.45\pm0.05^{\circ}$.  Three level GEM was implemented using these beams as illustrated in the inset of figure \ref{fig:setup}.  The probe and control fields were red detuned by $\Delta = 208$ MHz from the $\ket{5S_{1/2},F{=}2} \rightarrow \ket{5P_{1/2}, F'{=}1}$ and $\ket{5S_{1/2},F{=}1} \rightarrow \ket{5P_{1/2}, F'{=}1}$ of the \textsuperscript{87}Rb $D_1$ line respectively.  A pair of coils generated the magnetic field gradient required to inhomogeneously broaden the atomic transition to 200 kHz in the MOT to operate GEM.

The phase of the recall was measured using heterodyne detection.  An illustration of the pulse sequence used to perform this measurement is shown in figure \ref{fig:pulse}.  Before storage, a probe pulse was sent through the MOT with the control field switched off. Due to the detuning of the probe, this light passed through the atomic ensemble without interacting and was used as a phase reference.  After switching the control field on, a probe pulse was stored for 13 $\mu$s and recalled by switching the current in the gradient coils.  The XPM was generated by the signal pulse of duration $\tau$ and addressed the  $F{=}1 \rightarrow F'{=}2$ \textsuperscript{87}Rb $D_1$ transition with a detuning $\delta$.  This field was mode-matched to the probe and both coupled into the same fibre before being sent to the MOT to maximise its impact. The amount of XPM was measured by comparing recall from the memory with and without the signal pulse.

The a.c.-Stark shift can be integrated over the signal pulse duration for a simple analytical expression of the expected cross phase modulation \cite{Chen_Wang_Wang_Yu_2006}:

\begin{equation}
    \phi_{\textnormal{XPM}} = \frac{\Omega_s^2\delta\tau}{2\left(\gamma^2+\delta^2\right)}
    \label{eq:xpm}
\end{equation}
\noindent where $\Omega_s$, $\tau$ and $\gamma$ are the signal Rabi frequency, signal pulse duration and transition linewidth respectively.

\begin{figure*}
\centering
\includegraphics[scale=0.8]{setup2.pdf}
\caption{Experimental setup.  A Ti:Sapphire laser generated both the control and signal fields using acousto-optic modulators (AOM).  The probe field separated at 6.8 GHz was generated through modulation using a fibre coupled electro-optic modulator (FC-EOM) with a cavity to select the frequency shifted sideband.  The probe and signal were combined on a beamsplitter and sent through a quarter wave plate (QWP) to become circularly polarized before the probe is stored in a cold atomic ensemble of $^{87}$Rb atoms.  Gradient coils control the magnetic field gradient across the ensemble to store and recall the probe from the $\ket{5S_{1/2},F{=}1}$ state via the control beam.  The recalled probe was combined with a local oscillator (LO) beam for phase measurements via heterodyne detection on a photodetector (PD).  The inset shows the three level GEM memory scheme formed by the probe and control fields along with the signal field.}
\label{fig:setup}
\end{figure*}

\begin{figure*}
\centering
\includegraphics[scale=0.65]{Pulse_sequence_colors3.pdf}
\caption{Heterodyne measurement data for 3.7 pJ of signal energy at a detuning of $\delta = -8.7$ MHz.  a) Normalized data that has been filtered down to only the heterodyne beat frequency.  A bright pulse during the phase reference section is used to account for phase drifts between the data sets with and without signal.  A probe pulse is stored and recalled during \textquote{Memory write} and \textquote{Memory read} between which a signal pulse is sent through the memory. b) Raw data from the photodetector; the local oscillator was switched off during storage to measure the signal pulse.  The phase shift due to the signal pulse is shown with c) and d).  c) shows the phase reference pulse being in phase for both traces with and without the signal pulse and d) shows the phase of the recall being shifted due to the signal pulse.}
\label{fig:pulse}
\end{figure*}

\section{Cross phase modulation measurement}
Figure \ref{fig:evp}a shows the measured phase shift for a series of increasing signal pulse energies, along with a plot of what is expected using (1), at $\delta = -8.7$ MHz.  The energy of the pulse was varied by changing the pulse duration and measured by integrating over the photodetector signal over that duration.  The beam size in this experiment was measured to be 150$\pm$40 $\mu$m.  Good agreement was observed between the data and the model up to about 5 pJ signal energy (when fitted with a beam size of 190 $\mu$m which is within error).  Above this energy, there is significant discrepancy between the simple model and experimental data.  It appears as though the achievable phase shift saturates in a way that is not captured by (1).

To fully understand this phase saturation effect, the full model Maxwell-Bloch equations of this four level system were numerically solved.  These equations and more details on how they were used can be found in the Appendix (A.1).  It was observed that this saturation effect only occurs when there is some impurity in the initial state of the system; the storage state ($\ket{5S_{1/2},F{=}1}$) had some population before the memory operation.  Due to the relatively small control detuning of this experiment, the control field was capable of incoherently scattering a very small level of light into the probe mode.  As the signal energy is increased, the size of the recall is reduced as seen in figure \ref{fig:evp}b but the level of light incoherently scattered by the control field remains constant.  The phase saturation occurs when the control scattered light started to dominate over the reduced recall signal.  The initial population in the storage state was calculated by measuring the signal absorption at various signal frequencies (figure \ref{fig:sigabs} in the appendix) and was found to be 2\% of the total population (the simulation ignored other $m_F$ states in the same manifold so the 2\% population could be interpreted as there being 2 atoms in the storage state for every 98 atoms in the initial ground state).  As plotted in the figure \ref{fig:evp}a, the full numerical model using the above initial atomic population state ratio predicts phase shift saturation at $\sim$1.5 rad which agrees well with the experimental data.  The precise parameters of this model were optimized to fit experimental results using the same machine learner as in \cite{tranter2018} within realistic bounds set by the experiment.

The effect of changing the signal detuning was also investigated at 2.3 pJ signal energies as shown in figure \ref{fig:fvp} and good agreement between theory and experiment was also observed.  The error in the phase shift in both figures \ref{fig:evp}a and \ref{fig:fvp} were attributed to phase fluctuations at the heterodyne detector; the probe and local oscillator had to travel through separate optical fibres so their relative phase suffered from significant fluctuations, which can reach more than $2\pi$ rad from data set to data set.  The phase shift uncertainty in each data point on both plots were propagated from the standard deviation over 20 data-sets of the phase difference between the reference pulse and the recall (for with and without signal pulse).

\begin{figure}[h]
\centering
\includegraphics[scale=0.8]{energyvsphaseeff2.pdf}
\caption{a) Measured phase shift due to cross phase modulation at a series of signal pulse energies along with the theoretically predicted phase shifts using (1) at $\delta = -8.7$ MHz and the phase shift expected using a full model that accounts for an impure initial state using Maxwell-Bloch equations. b) Normalized recall efficiency for increasing signal pulse energies overlaid with the efficiency expected by the full model Maxwell-Bloch equations.}
\label{fig:evp}
\end{figure}

\begin{figure}[h]
\centering
\includegraphics[scale=0.7]{freqvsphase_fixed.pdf}
\caption{XPM phase shifts for different detunings, $\delta$, from the Stark shift transition at 2.3 pJ of signal energy overlaid with the theoretically predicted phase shifts using (1). }
\label{fig:fvp}
\end{figure}

\section{Discussion}
Extrapolating from figure \ref{fig:evp}a, a signal pulse containing a single photon will cause 0.07$\pm$0.02 $\mu$rad of phase shift in the stored probe.  This level of phase shift is still not large enough to perform the weak Kerr nonlinearities gate protocol Munro et al. proposed in \cite{Munro_Nemoto_Spiller_2005} that requires $10^{-5}$ rad.  Moreover our experiment shows absorption of the probe and signal fields that would further disadvantage a quantum gate operation.  We will now consider how these limitations may be addressed.

As seen in figure \ref{fig:evp}b, memory efficiency drops as the signal energy increases.  This is  because the atoms in the storage state ($F{=}1$) are pumped by the signal field leading to incoherent absorption of both probe and signal fields.  This effect was also observed in the numerical simulations using the Maxwell-Bloch equations.  The effect of the signal energy on the memory efficiency was particularly sensitive to memory bandwidth.  This was because the signal beam was focused through the ensemble and therefore had a spatially varying intensity across the ensemble; this imparted a spatially varying phase shift.  This hindered the rephasing of the atoms as the stored state was recalled.  Due to the nature of GEM storing frequency components of the input pulse across a spatial Fourier transform, smaller memory bandwidths meant the written atomic coherence was more spread out through the ensemble.  This caused the signal pulse to degrade the memory efficiency more for smaller memory bandwidths.

This probe and signal absorption, however, is not an inherent property of the a.c.-Stark shift scheme presented here.  The signal pulse reducing memory efficiency could be mitigated by either using larger memory bandwidths to localize the stored pulse in a smaller physical volume around the signal beam focus or use a less tightly focused signal beam.  As mentioned previously, it was determined that the atomic ensemble was initialized with some residual population in the $F{=}1$ storage state of the memory.  On the basis of our numerical model, this state impurity is responsible for most of the signal absorption (figure \ref{fig:sigabs}), as well as the saturation of the phase shift (seen in figure~\ref{fig:evp}a).  The signal absorption could, therefore, be strongly mitigated if the ensemble were prepared purely in the $F{=}2$ state.  An effectively 100\% pure initial state could be obtained via RF evaporation which is commonly used to create Bose-Einstein condensates \cite{andersonObservationBoseEinsteinCondensation1995}.  Crucially, the atoms transferred to the storage state via the memory operation will have a negligible effect on the signal absorption in a quantum gate application where the probe is a weak excitation. Signal absorption can be further suppressed if larger signal detunings were used at the expense of some phase shift magnitude.

Since there are clear ways to mitigate signal absorption, absorption of the probe light will be the limit on fidelity when using this system as a quantum gate. State impurity in this case will not be a factor since the scheme uses Fock states with very few photons in the signal field, as seen in figure~\ref{fig:evp}b.  The GEM protocol has to date shown a best efficiency of 87\% \cite{choHighlyEfficientOptical2016}.  This will be a limit if we are seeking to build unconditional quantum gates, where the loss of a photon means a loss of fidelity. In principle, GEM can attain efficiencies approaching 100\%.  The current operational limit, we believe, is due to inhomogeneity in the control field.  In future work we will address this issue as we aim for higher GEM efficiency.  Machine learning \cite{tranter2018} will also be explored as way to optimise the pulse timing and shaping.

In terms of the achievable nonlinearity, a simple path forward towards larger phase shift is to focus the signal field more tightly at the MOT to achieve higher signal intensities.  The relatively large beam size in this experiment leaves much room for improvement.  If the beam waist was tightened to $1 \mu$m, up to $2.5$ mrad of single photon phase shift could be observed.  Further improvements may also be observed if slow light or stationary light effects are applied to the signal field.

\section{Conclusion}
In conclusion, we have characterized the phase shift achievable via cross phase modulation facilitated by a cold-atom quantum memory.  0.07$\pm$0.02 $\mu$rad of phase shift per single photon was inferred by our measurement with excellent agreement with our modeling.  This measurement is a promising first step towards using quantum memories as deterministic quantum gates which is free from the limitations detailed in \cite{Gea-Banacloche_2010} for travelling wave cross phase modulation implementations.  Further improvements in phase shift of orders of magnitude is possible with this scheme through tighter beam focusing and application of slow or stationary light.

This research was conducted by the Australian Research Council Centres of Excellence Centre for Quantum Computation and Communication Technology (Grant No. CE170100012).

\section{Introduction}
Facilitating strong interactions between photons is an important step towards realizing optical quantum computing \cite{Milburn_1989}.  This can be achieved via materials possessing a very strong nonlinearity causing a photon to impart a phase shift onto another through cross phase modulation (XPM).  This process can be used for many interesting applications such as universal quantum gates \cite{Milburn_1989}, quantum non-demolition measurements of photon numbers \cite{Imoto_Haus_Yamamoto_1985} and cluster state generation \cite{Kim_Lee_Ji_Nha_Anisimov_Dowling_2015}.

Various proposals have demonstrated examples of strong nonlinear interactions using different platforms.  Early promising systems took advantage of cavity quantum electrodynamic (cQED) using single atoms \cite{Turchette_Hood_Lange_Mabuchi_Kimble_1995} and quantum dots \cite{Fushman_Englund_Faraon_Stoltz_Petroff_Vuckovic_2008}.  Over the last few years, successful demonstrations of large phase shifts were observed using cQED applied to an atomic ensemble to get phase shifts of up to $\pi/3$ per postselected single photon \cite{Beck_Hosseini_Duan_Vuletic_2016}, Rydberg  blockades \cite{Tiarks_Schmidt_Rempe_Durr_2016} for $\pi$ per photon and hollow core fibres filled with atomic gas \cite{Venkataraman_Saha_Gaeta_2013} for $0.3$mrad per photon.

The main approach for generating XPM in atomic ensemble systems has been through electromagnetically induced transparency (EIT) \cite{Chen_Wang_Wang_Yu_2006, Lo_Su_Chen_2010,Lo_Chen_Su_Chen_Chen_Chen_Yu_Chen_2011}.  The slow-light effect in EIT can be used to strengthen the nonlinear interaction in the medium as ac-stark shift imparts a phase shift and creates XPM.  Phase shifts of up to $18\mu$rad per post-selected single photon have been observed \cite{Feizpour_Hallaji_Dmochowski_Steinberg_2015}.  It has been shown, however, that the magnitude of phase shift obtainable between two travelling photons via EIT is limited \cite{Gea-Banacloche_2010}.  This limit could only be mitigated by operating the medium as a quantum memory to avoid simultaneously propagating optical fields or in the coherent photon conversion regime (a four-wave mixing effect that causes atomic population oscillations) which was recently demonstrated \cite{sagona-stophelConditionalPPhaseShift2020}.

In this paper, we measured the nonlinearity achieved by storing one light field using gradient echo memory (GEM) and propagating the other field through the medium during the storage phase.  We infer $0.07\mu$rad of phase shift per single photon.  While this phase shift is fairly modest, we present prospects for increasing the size of the phase shift through changes to experimental design.

\section{Method}
\subsection{Preparation of ensemble}
A cloud of \textsuperscript{87}Rb atoms is capture and cooled using a magneto-optical trap (MOT).  A 2D quadrupole field is used to create an elongated geometry to maximize optical depth along the propagation path of the probe (optical depth of 450 on probe transition of $D_1$ $\ket{5S_{1/2},F = 2} \rightarrow \ket{5P_{1/2}, F' = 1}$).  The optical table is surrounded by three pairs of orthogonal Helmholtz coils to cancel out static background fields.  The experimental sequence was synchronized to the 50Hz mains power to ensure environmental consistency from run to run.   Details on the ensemble preparation can be found in Ref. \cite{Sparkes_Bernu_Hosseini_Geng_Glorieux_Altin_Lam_Robins_Buchler_2013}.

\subsection{Measuring XPM using GEM}
The experimental setup is shown in Figure \ref{fig:setup}.  A weak probe was focused into the MOT along the elongated axis with a strong control field intersecting at an angle of $0.45\pm0.5^{\circ}$.  Three level GEM was implemented using these beams as illustrated in the inset of Figure \ref{fig:setup}.  The probe and control fields were red detuned by $\Delta = 208$MHz from the $\ket{5S_{1/2},F = 2} \rightarrow \ket{5P_{1/2}, F' = 1}$ and $\ket{5S_{1/2},F = 1} \rightarrow \ket{5P_{1/2}, F' = 1}$ of the \textsuperscript{87}Rb $D_1$ line respectively.  A pair of coils generated the magnetic field gradient required to inhomogeneously broaden the atomic transition to 200 kHz in the MOT to operate GEM.

The phase of the recall was measured using heterodyne detection.  An illustration of the pulse sequence used to perform this measurement is shown in Figure \ref{fig:pulse}.  Before storage, a probe pulse was sent through the MOT with the control field switched off. Due to the detuning of the probe, this light passed through the atomic ensemble without interacting and was used as a phase reference.  After switching the control field on, a probe pulse was stored for 13$ \mu$s and recalled by switching the current in the gradient coils.  The XPM was generated by the signal pulse of duration $\tau$ and addressed the  $F = 1 \rightarrow F' = 2$ \textsuperscript{87}Rb $D_1$ transition with a detuning $\delta$.  This field was mode-matched to the probe and both coupled into the same fibre before being sent to the MOT to maximise its impact. The amount of XPM was measured by comparing recall from the memory with and without the signal pulse.

To quantify the strength of the non-linearity in the memory, the phase shift was measured as a function of the signal pulse energy and frequency detuning.  The expected phase shift is modelled by:
\begin{equation}\label{starkeqn}
    \phi_{\textnormal{XPM}} = \frac{\Omega_s^2\delta\tau}{2\left(\gamma^2+\delta^2\right)}
\end{equation}

\noindent where $\Omega_s$, $\tau$ and $\gamma$ are the signal Rabi frequency, signal pulse duration and transition linewidth respectively.  This was obtained by integrating over the A.C. Stark shift over the signal duration \cite{Chen_Wang_Wang_Yu_2006}.

\begin{figure*}
\centering
\includegraphics[scale=0.8]{setup.pdf}
\caption{Experimental setup.  FC-EOM: Fibre coupled electro-optic modulator, AOM: acoustic-optic modulator, QWP: quarter waveplate, LO: local oscillator, PD: photodetector.  The probe, stark and control fields were controlled using AOMs and were all dervied from a Ti:sapphire laser running at 795nm.  The inset shows the $\Lambda$-GEM memory level scheme formed by the probe and control fields along with the signal field.}
\label{fig:setup}
\end{figure*}

\begin{figure*}
\centering
\includegraphics[scale=0.6]{Pulse_sequence_colors.pdf}
\caption{Heterodyne measurement data for 3.7pJ of signal energy at a detuning of $\delta = -8.7$MHz.  a) is data that has been normalized and filtered down to only the heterodyne beat frequency.  b) is raw data from the photodetector; the local oscillator was switched off during storage to measure the signal pulse.  The phase shift due to the signal pulse is shown with c) and d).  c) shows the phase reference pulse being in phase for both traces with and without the signal pulse and d) shows the phase of the recall being shifted due to the signal pulse.}
\label{fig:pulse}
\end{figure*}

\section{Cross phase modulation measurement}
Figure \ref{fig:evp} shows the measured phase shift for a series of increasing signal pulse energies along with a plot of what is expected by the model at $\delta = -8.7\textrm{MHz}$.  The energy of the pulse was varied by changing the pulse duration and measured by integrating over the photodetector signal over that duration.  The beam size in this experiment was measured to be $150\pm40\mu m$.  Good agreement was observed between the data and the model up to about 5pJ signal energy (when fitted with a beam size of $190\mu m$ which is within error).  The amount of phase shift achieved with higher signal energies above that seems to saturate.  To fully understand this phase saturation effect, the complete Maxwell-Bloch equations of this four level system was numerically solved.  It was observed that this saturation effect only occurs when there is some impurity in the initial state of the system; the storage state ($\ket{5S_{1/2},F = 1}$) had some population before the memory operation.  Due to the relatively small control detuning of this experiment, the control field was capable of incoherently scattering a very small level of light into the probe mode.  As the signal energy is increased, the recall efficiency is reduced as seen in Figure \ref{fig:eve}.  The phase saturation occurs when the control scattered light started to dominate over the reduced recall signal.  The initial population in the storage state was calculated by measuring the signal absorption at various signal frequencies and was found to be 2\% of the total population (the simulation ignored other $m_F$ states in the same manifold so the 2\% population could be interpreted as there being 2 atoms in the storage state for every 98 atoms in the initial ground state).  Using this as the initial state in the simulation gave the calculated full model phase shift in Figure \ref{fig:evp}.

The effect of changing the signal detuning was also investigated at 2.3pJ signal energies as shown in Figure \ref{fig:fvp} and good agreement between theory and experiment was also observed.

\begin{figure}
\centering
\includegraphics[scale=0.6]{energyvsphasefull.pdf}
\caption{Measured phase shift due to cross phase modulation at a series of signal pulse energies along with the theoretically predicted phase shifts using Equation \ref{starkeqn} at $\delta = -8.7$MHz and the phase shift expected from an impure initial state as modelled using Maxwell-Bloch equations.}
\label{fig:evp}
\end{figure}

\begin{figure}
\centering
\includegraphics[scale=0.6]{energyvseff.pdf}
\caption{Normalized recall efficiency for increasing signal pulse energies}
\label{fig:eve}
\end{figure}

\begin{figure}
\centering
\includegraphics[scale=0.6]{freqvsphase_fixed.pdf}
\caption{XPM phase shifts for different detunings, $\delta$, from the stark shift transition at 2.3pJ of signal energy.}
\label{fig:fvp}
\end{figure}

Extrapolating from Figure \ref{fig:evp}, a signal pulse containing a single photon will cause $0.07\pm0.02 \textrm{$\mu$rad}$ of phase shift in the stored probe.  A simple path forward towards larger phase shift is to focus the signal field more tightly at the MOT to achieve higher signal intensities.  The relatively large beam size in this experiment leaves much room for improvement.  If the beam waist was tightened to $1 \mu$m, up to $2.5$ mrad of single photon phase shift could be observed.  Further improvements may also be observed if slow light or stationary light effects are applied to the signal field.

\section{Conclusion}
In conclusion, we have characterized the phase shift achievable via cross phase modulation facilitated by a cold-atom quantum memory.  $0.07\pm0.02\textrm{$\mu$rad}$ of phase shift per single photon was inferred by our measurement with excellent agreement with our modeling.  This measurement is a promising first step towards using quantum memories as deterministic quantum gates which is free from the limitations detailed in \cite{Gea-Banacloche_2010} for travelling wave cross phase modulation implementations.  Further improvements in phase shift of orders of magnitude is possible with this scheme through tighter beam focusing and application of slow or stationary light.

\section*{References}

\section{Introduction}
Facilitating strong interactions between photons is an important step towards realizing optical quantum computing \cite{Milburn_1989}.  This can be achieved via materials possessing a very strong nonlinearity causing a photon to impart a phase shift onto another through cross phase modulation (XPM).  This process can be used for many interesting applications such as universal quantum gates \cite{Milburn_1989}, quantum non-demolition measurements of photon numbers \cite{Imoto_Haus_Yamamoto_1985} and cluster state generation \cite{Kim_Lee_Ji_Nha_Anisimov_Dowling_2015}.

Various proposals have demonstrated examples of strong nonlinear interactions using different platforms.  Early promising systems took advantage of cavity quantum electrodynamics (cQED) using single atoms \cite{Turchette_Hood_Lange_Mabuchi_Kimble_1995} and quantum dots \cite{Fushman_Englund_Faraon_Stoltz_Petroff_Vuckovic_2008}.  Over the last few years, successful demonstrations of large phase shifts were observed using cQED applied to an atomic ensemble to get phase shifts of up to $\pi/3$ per postselected single photon \cite{Beck_Hosseini_Duan_Vuletic_2016}, Rydberg  blockades \cite{Tiarks_Schmidt_Rempe_Durr_2016} for $\pi$ per photon and hollow core fibres filled with atomic gas \cite{Venkataraman_Saha_Gaeta_2013} for $0.3$mrad per photon.

The main approach for generating XPM in atomic ensemble systems has been through electromagnetically induced transparency (EIT) \cite{Chen_Wang_Wang_Yu_2006, Lo_Su_Chen_2010,Lo_Chen_Su_Chen_Chen_Chen_Yu_Chen_2011}.  The slow-light effect in EIT can be used to strengthen the nonlinear interaction in the medium as ac-stark shift imparts a phase shift and creates XPM.  Phase shifts of up to 18 $\mu$rad per post-selected single photon have been observed \cite{Feizpour_Hallaji_Dmochowski_Steinberg_2015}.  It has been shown, however, that the magnitude of phase shift obtainable between two travelling photons via EIT is limited \cite{Gea-Banacloche_2010}.  This limit could only be mitigated by operating the medium as a quantum memory to avoid simultaneously propagating optical fields or in the coherent photon conversion regime (a four-wave mixing effect that causes atomic population oscillations) which was recently demonstrated \cite{sagona-stophelConditionalPPhaseShift2020}.

In this paper, we measured the nonlinearity achieved by storing one light field using gradient echo memory (GEM) and propagating the other field through the medium during the storage phase.  We infer 0.07 $\mu$rad of phase shift per single photon.  While this phase shift is fairly modest, we present prospects for increasing the size of the phase shift through changes to experimental design.

\section{Method}
\subsection{Preparation of ensemble}
A cloud of \textsuperscript{87}Rb atoms is capture and cooled using a magneto-optical trap (MOT).  A 2D quadrupole field is used to create an elongated geometry to maximize optical depth along the propagation path of the probe (optical depth of 450 on probe transition of $D_1$ $\ket{5S_{1/2},F{=}2} \rightarrow \ket{5P_{1/2}, F'{=}1}$).  The optical table is surrounded by three pairs of orthogonal Helmholtz coils to cancel out static background fields.  The experimental sequence was synchronized to the 50Hz mains power to ensure environmental consistency from run to run.   Details on the ensemble preparation can be found in \cite{Sparkes_Bernu_Hosseini_Geng_Glorieux_Altin_Lam_Robins_Buchler_2013}.

\subsection{Measuring XPM using GEM}
The experimental setup is shown in figure \ref{fig:setup}.  A weak probe was focused into the MOT along the elongated axis with a strong control field intersecting at an angle of $0.45\pm0.5^{\circ}$.  Three level GEM was implemented using these beams as illustrated in the inset of figure \ref{fig:setup}.  The probe and control fields were red detuned by $\Delta = 208$MHz from the $\ket{5S_{1/2},F{=}2} \rightarrow \ket{5P_{1/2}, F'{=}1}$ and $\ket{5S_{1/2},F{=}1} \rightarrow \ket{5P_{1/2}, F'{=}1}$ of the \textsuperscript{87}Rb $D_1$ line respectively.  A pair of coils generated the magnetic field gradient required to inhomogeneously broaden the atomic transition to 200 kHz in the MOT to operate GEM.

The phase of the recall was measured using heterodyne detection.  An illustration of the pulse sequence used to perform this measurement is shown in figure \ref{fig:pulse}.  Before storage, a probe pulse was sent through the MOT with the control field switched off. Due to the detuning of the probe, this light passed through the atomic ensemble without interacting and was used as a phase reference.  After switching the control field on, a probe pulse was stored for 13 $\mu$s and recalled by switching the current in the gradient coils.  The XPM was generated by the signal pulse of duration $\tau$ and addressed the  $F{=}1 \rightarrow F'{=}2$ \textsuperscript{87}Rb $D_1$ transition with a detuning $\delta$.  This field was mode-matched to the probe and both coupled into the same fibre before being sent to the MOT to maximise its impact. The amount of XPM was measured by comparing recall from the memory with and without the signal pulse.

The a.c. stark shift can be integrated over the signal pulse duration for a simple analytical expression of the expected cross phase modulation \cite{Chen_Wang_Wang_Yu_2006}:

\begin{equation}\label{starkeqn}
    \phi_{\textnormal{XPM}} = \frac{\Omega_s^2\delta\tau}{2\left(\gamma^2+\delta^2\right)}
    \label{eq:xpm}
\end{equation}
\noindent where $\Omega_s$, $\tau$ and $\gamma$ are the signal Rabi frequency, signal pulse duration and transition linewidth respectively.

\begin{figure*}[h]
\centering
\includegraphics[scale=0.8]{setup2.pdf}
\caption{Experimental setup.  FC-EOM: Fibre coupled electro-optic modulator, AOM: acoustic-optic modulator, QWP: quarter waveplate, LO: local oscillator, PD: photodetector.  The probe, stark and control fields were controlled using AOMs and were all dervied from a Ti:sapphire laser running at 795nm.  The inset shows the three level GEM memory scheme formed by the probe and control fields along with the signal field.}
\label{fig:setup}
\end{figure*}

\begin{figure*}
\centering
\includegraphics[scale=0.6]{Pulse_sequence_colors2.pdf}
\caption{Heterodyne measurement data for 3.7pJ of signal energy at a detuning of $\delta = -8.7$MHz.  a) is data that has been normalized and filtered down to only the heterodyne beat frequency.  b) is raw data from the photodetector; the local oscillator was switched off during storage to measure the signal pulse.  The phase shift due to the signal pulse is shown with c) and d).  c) shows the phase reference pulse being in phase for both traces with and without the signal pulse and d) shows the phase of the recall being shifted due to the signal pulse.}
\label{fig:pulse}
\end{figure*}

\section{Cross phase modulation measurement}
Figure \ref{fig:evp} shows the measured phase shift for a series of increasing signal pulse energies, along with a plot of what is expected using (1), at $\delta = -8.7\textrm{MHz}$.  The energy of the pulse was varied by changing the pulse duration and measured by integrating over the photodetector signal over that duration.  The beam size in this experiment was measured to be 150$\pm$40 $\mu$m.  Good agreement was observed between the data and the model up to about 5 pJ signal energy (when fitted with a beam size of 190 $\mu$m which is within error).  Above this energy, there is significant discrepancy between the simple model and experimental data.  It appears as though the achievable phase shift saturates in a way that is not captured by (1).

To fully understand this phase saturation effect, the complete Maxwell-Bloch equations of this four level system was numerically solved.  These equations and more details on how they were used can be found in the Appendix (A.1).  It was observed that this saturation effect only occurs when there is some impurity in the initial state of the system; the storage state ($\ket{5S_{1/2},F{=}1}$) had some population before the memory operation.  Due to the relatively small control detuning of this experiment, the control field was capable of incoherently scattering a very small level of light into the probe mode.  As the signal energy is increased, the size of the recall is reduced as seen in figure \ref{fig:eve} but the level of light incoherently scattered by the control field remains constant.  The phase saturation occurs when the control scattered light started to dominate over the reduced recall signal.  The initial population in the storage state was calculated by measuring the signal absorption at various signal frequencies (figure \ref{fig:sigabs} in the appendix) and was found to be 2\% of the total population (the simulation ignored other $m_F$ states in the same manifold so the 2\% population could be interpreted as there being 2 atoms in the storage state for every 98 atoms in the initial ground state).  Using this as the initial state in the simulation gave the calculated full model phase shift in figure \ref{fig:evp}.

The effect of changing the signal detuning was also investigated at 2.3 pJ signal energies as shown in figure \ref{fig:fvp} and good agreement between theory and experiment was also observed.

\begin{figure}[h]
\centering
\includegraphics[scale=0.8]{energyvsphasefull.pdf}
\caption{Measured phase shift due to cross phase modulation at a series of signal pulse energies along with the theoretically predicted phase shifts using (1) at $\delta = -8.7$MHz and the phase shift expected from an impure initial state as modelled using Maxwell-Bloch equations.}
\label{fig:evp}
\end{figure}

\begin{figure}[h]
\centering
\includegraphics[scale=0.8]{energyvseff3.pdf}
\caption{Normalized recall efficiency for increasing signal pulse energies overlaid with the efficiency simulated using the Maxwell-Bloch equations.}
\label{fig:eve}
\end{figure}

\begin{figure}[h]
\centering
\includegraphics[scale=0.8]{freqvsphase_fixed.pdf}
\caption{XPM phase shifts for different detunings, $\delta$, from the stark shift transition at 2.3pJ of signal energy.}
\label{fig:fvp}
\end{figure}

\section{Discussion}
Extrapolating from figure \ref{fig:evp}, a signal pulse containing a single photon will cause 0.07$\pm$0.02 $\mu$rad of phase shift in the stored probe.  This level of phase shift is still not large enough to perform the weak Kerr nonlinearities gate protocol Munro et al. proposed in \cite{Munro_Nemoto_Spiller_2005} that requires $10^{-5}$ rad.  Moreover our experiment shows absorption of the probe and signal fields that would further disadvantage a quantum gate operation.  We will now consider how these limitations may be addressed.

As seen in figure \ref{fig:eve}, memory efficiency drops as the signal energy increases.  This is  because the atoms in the storage state ($F{=}1$) are pumped by the signal field leading to incoherent absorption of both probe and signal fields.  This effect was also observed in the numerical simulations using the Maxwell-Bloch equations.  The effect of the signal energy on the memory efficiency was particularly sensitive to memory bandwidth.  This was because the signal beam was focused through the ensemble and therefore had a spatially varying intensity across the ensemble; this imparted a spatially varying phase shift.  This hindered the rephasing of the atoms as the stored state was recalled.  Due to the nature of GEM storing frequency components of the input pulse across a spatial fourier transform, smaller memory bandwidths meant the written atomic coherence was more spread out through the ensemble.  This caused the signal pulse to degrade the memory efficiency more for smaller memory bandwidths.

This probe and signal absorption, however, is not an inherent property of the ac stark shift scheme presented here.  The signal pulse reducing memory efficiency could be mitigated by either using larger memory bandwidths to localize the stored pulse in a smaller physical volume around the signal beam focus or use a less tightly focused signal beam.  As mentioned previously, it was determined that the atomic ensemble was initialized with some residual population in the $F{=}1$ storage state of the memory.  On the basis of our numerical model, this state impurity is responsible for most of the signal absorption (figure \ref{fig:sigabs}), as well as the saturation of the phase shift (seen in figure~\ref{fig:evp}).  The signal absorption could, therefore, be strongly mitigated if the ensemble were prepared purely in the $F{=}2$ state.  An effectively 100\% pure initial state could be obtained via RF evaporation which is commonly used to create Bose-Einstein condensates \cite{andersonObservationBoseEinsteinCondensation1995}.  Crucially, the atoms transferred to the storage state via the memory operation will have a negligible effect on the signal absorption in a quantum gate application where the probe is a weak excitation. Signal absorption can be further suppressed if larger signal detunings were used at the expense of some phase shift magnitude.

Since there are clear ways to mitigate signal absorption, absorption of the probe light will be the limit on fidelity when using this system as a quantum gate. State impurity in this case will not be a factor since the scheme uses Fock states with very few photons in the signal field, as seen in figure~\ref{fig:eve}.  The GEM protocol has to date shown a best efficiency of 87\% \cite{Cho_Campbell_Everett_Bernu_Higginbottom_Cao_Geng_Robins_Lam_Buchler_2016}.  This will be a limit if we are seeking to build unconditional quantum gates, where the loss of a photon means a loss of fidelity. In principle, GEM can attain efficiencies approaching 100\%.  The current operational limit, we believe, is due to inhomogeneity in the control field.  In future work we will address this issue as we aim for higher GEM efficiency.  Machine learning \cite{tranter2018} will also be explored as way to optimise the pulse timing and shaping.

In terms of the achievable nonlinearity, a simple path forward towards larger phase shift is to focus the signal field more tightly at the MOT to achieve higher signal intensities.  The relatively large beam size in this experiment leaves much room for improvement.  If the beam waist was tightened to $1 \mu$m, up to $2.5$ mrad of single photon phase shift could be observed.  Further improvements may also be observed if slow light or stationary light effects are applied to the signal field.

\section{Conclusion}
In conclusion, we have characterized the phase shift achievable via cross phase modulation facilitated by a cold-atom quantum memory.  0.07$\pm$0.02 $\mu$rad of phase shift per single photon was inferred by our measurement with excellent agreement with our modeling.  This measurement is a promising first step towards using quantum memories as deterministic quantum gates which is free from the limitations detailed in \cite{Gea-Banacloche_2010} for travelling wave cross phase modulation implementations.  Further improvements in phase shift of orders of magnitude is possible with this scheme through tighter beam focusing and application of slow or stationary light.
\ack
This research was conducted by the Australian Research Council Centres of Excellence Centre for Quantum Computation and Communication Technology (Grant No. CE170100012).