abstract: We study the dynamics of the interaction between the solar wind ions and a partially ionized atmosphere around a comet, at a distance of 2.88 AU from the sun during a period of low nucleus activity. Comparing particle data and magnetic field data for a case study, we highlight the prime role of the solar wind electric field in the cometary ion dynamics. Cometary ion and solar wind proton flow directions evolve in a correlated manner, as expected from the theory of mass loading. We find that the main component of the accelerated cometary ion flow direction is along the anti-sunward direction, and not along the convective electric field direction. This is interpreted as the effect of an anti-sunward polarisation electric field adding up to the solar wind convective electric field.
title: Mass loading at 67P/Churyumov-Gerasimenko: a case study.

```latex
\begin{article}

%% ------------------------------------------------------------------------ %%
%
%  TEXT
%
%% ------------------------------------------------------------------------ %%

\section{Introduction}
The phenomenon of mass loading is common in space plasmas. Newly charged material added to the fast solar wind flow is accelerated by the Lorentz force. The newly added material gains energy and momentum from the solar wind. Solar wind ions experience an equal but opposite net force, thus balancing the total momentum of the system. The thin atmosphere permeated by the solar wind around a comet is one of the most evident cases where mass loading is expected to control the dynamics of the plasma environment (\citet{szego2000ssr}, section 4.1). {\it In situ} investigations of the solar wind interaction with a cometary plasma were made possible by different missions prior to Rosetta: ICE at P/Giacobini--Zinner in 1985, Giotto, Vega-1 and -2, Suisei, Sakigake at P/Halley in 1986, and Giotto at P/Grigg--Skjellerup in 1992, are some examples. However, all those measurements were performed during single flybys, at about 1 A.U. away from the Sun, and bow shocks were observed at each flyby, indicating a high nucleus activity ({\it cf.} \citet{neugebauer1990rg} and \citet{coates1997asr}).

The Rosetta mission \citep{glassmeier2007ssr} has provided a unique opportunity to continuously observe mass loading in a cometary environment over longer time scales and during varying nucleus activity. The Rosetta spacecraft reached comet 67P/Churyumov--Gerasimenko (67P/CG) in early August 2014. First results from the plasma measurements made at 67P/CG describe how the cometary environment evolves from a thin coma where only low fluxes of low energy ions are observed, to the point when the effect on the solar wind flow becomes significant \citep{nilsson2015aa}, \citep{nilsson2015science}, \citep{goldstein2015grl}. During these early observations no plasma boundaries had yet formed between the solar wind and the comet atmosphere, {\it i.e.} there was no bow shock or ionopause. The scale size of the interaction observed at comet 67P/CG was small, initial observations of water ions were made only when the spacecraft got closer to the nucleus than 100 km distance.
The solar wind -atmosphere interaction at a low activity comet may thus have similarities to artificial comets formed through barium and lithium ion releases from the AMPTE spacecraft \citep{haerendel1986nature, rodgers1986nature, coates1986jgr, coates2015jpcs}.

For a low activity comet the solar wind is undisturbed before permeating the coma, and no other acceleration process has to be taken into account to study the solar wind - atmosphere interaction. We expect to observe the simplest mass loading phenomenon: newly charged mass is simply added to the undisturbed solar wind flow. This interaction has been previously addressed by \citet{broiles2015aa}, using data from another particle instrument within the Rosetta Plasma Consortium (RPC), RPC-IES \citep{burch2007ssr}. They reported that the solar wind near the comet was deflected by a Lorenz Force opposite to that experienced by cometary pickup ions. They also found that this deflection was not well ordered by the spacecraft position relative to the comet, and was well correlated with large changes in the observed magnetic field.

We present a case study that provides new details about the dynamics of this interaction between the solar wind and the coma, based on ion and magnetic field data from the 28th of November 2014.

%------------------------------------------------

\section{Instrument description}

The Ion Composition Analyzer, part of the Rosetta Plasma Consortium (RPC-ICA), is an ion spectrometer aimed to study the interaction between the solar wind and positive cometary ions at comet 67P/CG \citep{nilsson2007ssr}. The instrument resolves energy and mass per charge of the incoming ions. The energy spans from 10 eV up to 40 keV. The instrument field of view is $360^{\circ} \times 90^{\circ}$ ($azimuth\ \times\ elevation$, illustrated in Figure \ref{fov}), with a resolution of $22.5^{\circ} \times 5.0^{\circ}$. For this data set, a full angular scan was produced every 192 s. The elevation angle of the incoming positive ions is determined by an electrostatic acceptance angle filter at the entrance of the instrument, and the azimuth angle is measured by the means of 16 anodes, part of the detection system marking the end of an ion path in the instrument. Thus the two angles are subject to different constraints, limits, and resolutions.

The magnetometer (RPC-MAG, \cite{glassmeier2007ssr-MAG}) measures the three components of the magnetic field vector in the range from DC up to 10 Hz. The measurement range is $\pm$ 16384 nT with a resolution of 31 pT. RPC-MAG is mounted on a 1.5 m long boom in order to minimise the impact of the spacecraft generated disturbance fields. The magnetometer is affected by a systematic bias field from the spacecraft, which cannot  be fully characterized as the spacecraft fields are changing related to the operation status. In this work, we take this into account by considering and propagating a $\pm$ 3 nT uncertainty on each magnetic field component in the spacecraft reference frame.

%------------------------------------------------

\section{The case}

In order to diagnose the dynamics of the interaction between the solar wind and the coma, we considered data sets with a clear and constant water ion signal, simultaneously with a clear solar wind signal. The chosen data set turned out to be the clearest in terms of dynamics, mainly because of large variations in the upstream magnetic field direction.

The following case study is based on particle and magnetic field data collected at 2.88 AU on the 28th of November 2014. The spacecraft was then flying a terminator orbit 30 km away from the center of the nucleus, with a mean speed over the day of 0.15 m/s relative to the comet. This very low speed allows us to neglect any aberration angle concerning the cometary ions.

Moreover, an aberration of 1.01$^{\circ}$ is obtained for a 400 km/s fast solar wind. Our angular resolution is larger than this angle, so we neglect this aberration angle as well.

To study the dynamics of the interaction, we focus on two species: solar wind protons $H^+$ and cometary water ions $H_2O^+$, the most abundant species in the solar wind and the ionized atmosphere respectively \citep{nilsson2015science}. A large flux of accelerated water ions is observed during the chosen day, and we consider the energy range [70 eV, 330 eV], with a peak value at 100 eV. The upper bound is set higher than the most energetic water ions observed. The lower bound isolates these accelerated water ions from the cold water ion population reported by \citet{nilsson2015aa}: this population is affected by the spacecraft potential in terms of direction, and therefore is not physically relevant for our study. Protons are observed in the energy range [350 eV, 1200 eV].

   The average magnetic field magnitude is 14 nT $\pm$ 5.2 nT over the 14 hours of the data set. The speed of the observed protons is stable, with a value of 380 km/s for a peak at 750 eV and a standard deviation $\sigma$ = 20 km/s.

 As a first approximation, we assume that the solar wind electric field is given by ${\bf E} = -{\bf v}_{H^+} \times {\bf B}$. Integrating the movement of a test particle with a velocity $v_{test}$ ,  only subject to the Lorentz force ${\bf F} = q({\bf E} + {\bf v_{test}} \times {\bf B}) $, we estimate that the most energetic cometary ions observed ($\sim$ 300 eV) were accelerated during less than 3 s, over a distance of about 40 km. With this approximation, the minimum gyro-radius is 2600 km, which means that the $H_2O^+$ we observe are on the very early phase of the gyro-motion, {\it i.e.} flowing along the local electric field. For this reason, we expect that the observed accelerated water ion flow gives us the direction of the local electric field. These numbers are summarized in Table \ref{gyro}.

 To complete the description of the environment during the measurement, we estimate the density profile of the coma using the Haser model (\citet{haser1957}). In this model, the neutral density falls off following $\sim 1/r^2$, with $r$ the distance to the nucleus. At 100 km away from the nucleus, the densities would be an order of magnitude lower than the one met along the 30 km terminator orbite.

\begin{table}
\caption{Conditions and estimations of $H_2O^+$ gyro-radius \label{gyro}}
\centering
\begin{tabular}{r | c | c}

   & Data & Undisturbed  \\
   && solar wind  \\
\hline
\hline
  $B$ & [12, 20] nT & 1 nT \\
  $v_{H^+}$ & 380 km/s & 400\ km/s \\
  $ E_{estim}$ & [4.0,\ 4.6] V/km & 3.7 V/km \\
\hline
  $R_{gyro}$ & [2600, 5900] km & -
\end{tabular}

\end{table}

%------------------------------------------------

\section{Method}

To study the dynamics of the interaction between the solar wind and the comet atmosphere, we aim to express the flow directions of $H^+$ and $H_2O^+$ in the body-Centered Solar EQuatorial (CSEQ) frame. To achieve that, we first collect the observed counts for the two ion populations, which are well separated in both energy and mass. We then compute full angular distributions for each species every 192 s.

To clearly visualize the two flows and their dynamics, we produce a sequence of pictures as seen with the instrument Field Of View (FOV). Three of the pictures composing the sequence are given in Figure \ref{fov}; the full sequence can be seen at http://irf.se/$\sim$etienne/mediaFOV.html . To that extent, the instrument is used as a camera. The produced sequence also helps us assess that the restricted field-of-view of the instrument does not impact our results.

\begin{figure*}
   \begin{center}
      \includegraphics[width=\textwidth]{FOV.pdf}
      \caption{Three examples of $H_2O^+$ and $H^+$ flows pictured in the instrument FOV at precise times during the 28th of November. Depending on the time, the two flows are coming from completely different directions (the spacecraft attitude is unchanged during this period), but always with a $180 ^{\circ}$ difference in clock angle as defined in the lower right illustration. It represents the spherical coordinate system used for this study. All vectors are normalized, ${\bf v}_{norm.} = \frac{\bf v}{|{\bf v}|}$ \label{fov}}
   \end{center}
\end{figure*}

We calculate the direction of the bulk velocity for both species, then express the two flow directions (unit vectors) in the CSEQ reference frame. The $x_{CSEQ}$ axis is along the sun-comet line, pointing to the sun. The $z_{CSEQ}$ axis is parallel to and oriented by the Sun's north pole direction, orthogonal to the +X axis. The $y_{CSEQ}$ completes the right-handed reference frame. We compute the angle of each flow direction from the sun-comet line, and refer to this angle as the cone angle. We complement it with the clock angle expressed in the $(y_{CSEQ}, \ z_{CSEQ})$ plane. It is the angle of the projected flow direction in this plane, from the $y_{CSEQ}$ axis and positive towards $z_{CSEQ}$, as shown in Figure \ref{fov}, lower right illustration. Both angles form a spherical coordinate system, and all vectors are normalized in this study, ${\bf v}_{norm.} = \frac{\bf v}{|{\bf v}|}$. The magnetic field is also expressed in this coordinate system.\\

%------------------------------------------------

\section{Results and discussion}

   If mass loading is the only mechanism in the interaction between the solar wind and the coma, we expect to see cometary ions accelerated along the local electric field direction and solar wind protons deflected with an opposite clock angle, as a result of momentum conservation. The dynamics takes place in the plane that includes both flow directions and the comet-sun line. As the direction of the solar wind electric field varies, both flows are expected to follow the field rotation, $180 ^{\circ}$ away from each other in terms of clock angle. We thus verify with only one coordinate if the two ion populations flow in the same plane or not. The cone angle quantifies the anti-sunward component in this plane. It gives the amount of solar wind proton deflection, and the anti-sunward component of the cometary ion velocity.

   \subsection{Clock angle, momentum conservation}

   In Figure \ref{comp}, clock angles for $H^+$ (red dots),  $H_2O^+$ (blue dots)  and ${\bf B}$ (black line) are given in the top panel.  All three clock angles show variations over a span of about  $200 ^{\circ}$. Despite these large variations, the difference between proton and cometary ion flow directions remain around  $180 ^{\circ}$, and is given in the center panel.
   The mean value of this angular difference is $187 ^{\circ}$ and its standard deviation is $18 ^{\circ}$. The proton flow clock angle is more correlated to the magnetic field clock angle than the cometary ion flow is. The variations thus mainly originate from the cometary ion flow direction, which is most of the time less regular, or beam-like ({\it cf.} Figure \ref{fov}, lower left FOV picture: there are two peak values instead of one). An energy dependence in the water ion flow direction could explain these variations, and will be the topic of a more advanced study focused on the cometary ion population.

   The two observed flows have very different directions. Thus the two populations, $H^+$ and $H_2O^+$, have been flowing along different paths in the coma. Observing this $180 ^{\circ}$ difference ensures us that the magnetic and electric fields have respectively very similar directions along these two different paths. As the upstream electric field direction varies, the two paths  also move inside the coma, but the angular difference remains the same. This also confirms that the flow directions of both populations remain in the same plane, containing the comet-sun line.

The direction of the deflection doesn't seem to be influenced by the position of the nucleus, the flow may be deflected away just as well as towards the nucleus. There is no indication of the ion population flowing around the nucleus. This is because the size of the the nucleus is much smaller than the gyroradii of solar wind ions. To be able to flow around an obstacle, the particles of the flow must have sufficient time  and space to interact several times inside a boundary layer. The interactions inside the boundary layer mediate the effect which "pushes away" the flow in a conventional situation. There is no such layer in the case of the new born cometary magnetosphere. The flow only experiences electromagnetic forces due to the mass loading, it feels the much larger inner coma that way, but it simply cannot feel the nucleus. Solar wind ions are bombarding the surface of the nucleus instead of flowing around it. In this regard this early stage of the interaction is very different from later stages.

   The dynamics is clearly driven by the solar wind electric and magnetic fields. The nearly constant difference in clock angle between $H_2O^+$ and $H^+$ and the correlation with magnetic field clock angles emphasize the role of mass loading as the main mechanism controlling the dynamics of this plasma environment.

   \subsection{Cone angle}

   The cone angles of $H^+$, $H_2O^+$ and ${\bf B}$ are given in the third panel of Figure \ref{comp}. An anti-correlation between the $H_2O^+$ flow direction and the magnetic field can be seen: when the magnetic field cone angle decreases (i.e. the magnetic field is more parallel to the sun-comet line), the $H_2O^+$ cone angle increases, and vice-versa.

      Concerning the proton flow, the cone angle seems to be correlated with the magnetic field cone angle in the first half of the data set. Around 9:40, a sudden and large decrease is observed in both angles. This correlation has an immediate interpretation. $B_\perp$ is the component of the magnetic field orthogonal to the proton flow. The larger $B_{\perp}$ is (for the same $B$ amplitude), the larger ${\bf E} = -{\bf v}_{H^+} \times {\bf B} $ is. The acceleration of new born ions is then larger and as a direct consequence the deflection of the protons is also larger. The correlation is not clearly observed after 12:00.

     On average, $H_2O^+$ flow direction is $20 ^{\circ}$ away from the sun-comet line, and a $40 ^{\circ}$ solar wind proton deflection is observed, as illustrated in Figure \ref{result}.\\

      \subsection{Electric field}

     The generalized Ohm's Law can be written as following:

      \[ {\bf E} = {\bf E}_{motional} + {\bf E}_{other}\]

      As previously introduced, we approximate the direction of the motional field by ${\bf E}_{motional} = -{\bf v}_{H^+} \times {\bf B}$, computing together particle data and magnetic field data. For two earlier cases (data from 2014-09-21 --{\it cf.} \citet{nilsson2015science}-- and 2014-11-16), the $H_2O^+$ ions are flowing along $-{\bf v}_{H^+} \times {\bf B}$ . ${\bf E} = {\bf E}_{motional} = -{\bf v}_{H^+} \times {\bf B}$ was then a good description of the local electric field. These data sets correspond to a lower activity and larger distances to the nucleus, resulting in a lighter mass loading of the solar wind. In the extreme case where only few cometary ions are picked up, the solar wind is almost undisturbed by the interaction, and those few pick up ions are first accelerated along the convective electric field given by ${\bf E} = -{\bf v}_{H^+} \times {\bf B}$.

       The case studied here corresponds to a more significant mass loading of the solar wind and reveals a new configuration of the electric and magnetic fields inside the coma. The result for the comparison between the $-{\bf v}_{H^+} \times {\bf B} $ direction and the $H_2O^+$ flow direction is given in Figure \ref{comp}, using the same spherical coordinate system previously introduced. In terms of clock angle (first panel), the two series are overlapping nicely. Their angular difference, given in the second panel, is centered around the solid black line at $0^{\circ}$. We note that between 09:45 and 11:00, after the sudden decrease in the magnetic field cone angle, this difference is centered around $ \sim 45 ^{\circ}$. This may indicate that a further specific study of this quasi-parallel regime is needed.

     The cone angles give us a more complex picture of the dynamics. $-{\bf v}_{H^+} \times {\bf B} $ and $H_2O^+$ flow cone angles seem to correlate with each other, but with an angular difference of about $30 ^{\circ}$. So it appears that with the activity increasing, ${\bf E} = -{\bf v}_{H^+} \times {\bf B}$ doesn't hold anymore.  Two possibilities can explain the discrepancy:
      \begin{itemize}
      \item In the description of the motional electric field ${\bf E}_{motional} = -{\bf v} \times {\bf B}$ , ${\bf v_{H^+}}$ is not a good proxy anymore. Electron and proton flows do not have the same direction anymore.
      \item ${\bf E}_{other}$ became comparable to ${\bf E}_{motional}$ .
      \end{itemize}

      Both possibilities can be directly put in relation with the increasing cometary ion density, and it seems reasonable to say that both are playing a role. The magnetic field is enhanced (with an average value of 14 nT, compared to the $\sim 1$ nT expected in the undisturbed solar wind at this heliocentric distance), while the solar wind protons are not significantly slowed down (the observed proton flow speed is 380 km/s). This implies that the electron fluid is significantly slowed down. There is no obvious reason why this should lead to an anti-sunward electric field. Therefore, we don't explore this first possibility, and confine the discussion to the possibility of another electric field developing in the coma when densities are higher, resulting in an anti-sunward acceleration of picked up cometary ions.\\

\begin{figure*}
   \begin{center}
      \includegraphics[width=\textwidth]{1128compNew.pdf}
      \caption{Clock (top panel) and cone (bottom panel) angle time series, for the $H^+$(red) and $H_2O^+$ (blue) flow directions, and for the magnetic field direction (solid black line). The difference between the two flow clock angles is given (center panel), with the two solid lines indicating $y=180 ^{\circ}$ and $y=0 ^{\circ}$. Finally, the $-{\bf v}_{H^+} \times {\bf B}$ direction is given in green, and its angular difference with cometary ion flow is given in the center panel in terms of clock angle. Between 11:00 and 13:00, no particle data are available. \label{comp}}
   \end{center}
\end{figure*}

\begin{figure}
   \begin{center}
      \includegraphics[width=.5\textwidth]{resultCase.pdf}
      \caption{Illustration of the observation: the solar wind protons (red line) are gradually deflected by the coma, reaching a $40 ^{\circ}$ angle from the comet-sun line in average, measured in the terminator plane. The main component of the cometary ion flow is along the comet-sun line. The dynamics takes place in a plane, containing the comet-sun line and rotating around it. \label{result} }
   \end{center}
\end{figure}

    The cometary atmosphere provides a constant and distributed source of ions and electrons. The electrons are picked up and swept downstream (${\bf E} \times {\bf B}$ drift), whereas the new-born ions are moving along the local electric field direction. Charge separation occurs because of the scale size of the interaction region being much smaller than the gyro radius of a new born water ions. Whatever electric fields arise from this, must in the end have the net effect of reducing this separation. The new net polarisation field is in the plane $ ({\bf v}_{H_2O^+}, {\bf v}={\bf E} \times {\bf B})$ and therefore cannot have any influence on the clock angle of the $H_2O^+$  and $H^+$ ion flow directions. It has however a large influence on their cone angle, as seen in the observation presented in the previous sub-section.

    If we now consider the solar wind, we expect that ions and electrons will also react differently to the obstacle on some scale. As previously mentioned, observations show no significant slowing down of the solar wind while we see a significant increase in the magnetic field as compared to the undisturbed solar wind. That implies decoupling between electrons and protons, which could lead to a Hall current or even charge separation.

    To get the full picture of this interaction and its dynamics, one would need to add the two components discussed here-above. The situation becomes very complex, and we believe that numerical simulations are needed to describe the interaction. We observe a mainly anti-sunward cometary ion flow, implying a mainly anti-sunward $E_{other}$. This is consistant with the mechanism discussed in the first paragraph.

     \cite{haerendel1986nature} depict such an anti-sunward electric field in the context of the AMPTE artificial comet, introducing it with the same argument of charge separation. This work was based on magnetic field data and particle measurements made by a similar 3D plasma instrument ({\it cf.} \citet{rodgers1986nature} and \citet{coates1986jgr}). This topic is also widely reviewed by \citet{szego2000ssr}.
      However, the AMPTE releases concerned transient phenomena, [3] whereas Rosetta investigate the evolution of a continuous and time varying source of neutral particles.

%------------------------------------------------

\section{Conclusion}

   Despite a limited angular resolution and uncertainties on the magnetic field components, the obvious correlations and anti-correlations presented in this study provide a detailed description of the interaction between the solar wind and the cometary atmosphere, during a period of low nucleus activity.
   The plasma environment dynamics pictured in this work does not depend on the nucleus position: when the solar wind magnetic field direction changes, the two flow directions evolve accordingly, regardless of the nucleus position, {\it i.e.} the ions don't flow around the nucleus. This is in agreement with earlier results, {\it cf.} \citet{broiles2015aa}. Both ion populations -- from the solar wind and from the coma -- flow with clock angles $180 ^{\circ}$ apart from each other: the dynamics of the interaction takes place in a single plane, that contains and rotates around the comet-sun line.

   Measured at 30 km away from the nucleus and at a distance of 2.88 AU from the Sun, the solar wind protons reach deflections higher than $50 ^{\circ}$. Even in the absence of plasma boundaries, the solar wind is already substantially disturbed, and the magnetic field is not coupled to the solar wind ions anymore, contrary to what is observed at greater distances to the Sun.

The comet ions are accelerated with a dominating anti-sunward component. We suggest that the difference of motion between the cometary ions and electrons, together with the limited scale size of the coma, smaller than the cometary ion gyro-radius, result in a mainly anti-sunward new contribution to the electric field.

%%% End of body of article:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Optional Appendix goes here
%
% \appendix resets counters and redefines section heads
% but doesn't print anything.
% After typing \appendix
%
%\section{Here Is Appendix Title}
% will show
% Appendix A: Here Is Appendix Title
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Optional Glossary or Notation section, goes here
%
%%%%%%%%%%%%%%
% Glossary is only allowed in Reviews of Geophysics
% \section*{Glossary}
% \paragraph{Term}
% Term Definition here
%
%%%%%%%%%%%%%%
% Notation -- End each entry with a period.
% \begin{notation}
% Term & definition.\\
% Second term & second definition.\\
% \end{notation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  ACKNOWLEDGMENTS

\begin{acknowledgments}
The work on RPC-ICA, as well as this PhD project, is founded by the Swedish National Space Board.
The work on RPC-MAG was financially supported by the German Ministerium f\"ur Wirtschaft und Energie and the Deutsches Zentrum f\"ur Luft- und Raumfahrt under contract 50QP 1401 .
Without the tremendous work of the Rosetta Science Ground Segment (RSGS), Rosetta Mission Operation Control (RMOC) and all instrument team planners, this study would be impossible.
Sharing data within RPC is made possible by the web-based interface AMDA, developed and made available for RPC use by Centre de Donn\'ees de la Physique des Plasmas (CDPP). This efficient interface has been of a great use for this work.
The SPICE toolkit from NASA's Navigation and Ancillary Information Facility (NAIF) is heavily used in this study, and made accessible through python by the SpiceyPy wrapper.
\end{acknowledgments}

%% ------------------------------------------------------------------------ %%
%%  REFERENCE LIST AND TEXT CITATIONS
%
% Either type in your references using
% \begin{thebibliography}{}
% \bibitem{}
% Text
% \end{thebibliography}
%
% Or,
%
% If you use BiBTeX for your references, please use the agufull08.bst file (available at % ftp://ftp.agu.org/journals/latex/journals/Manuscript-Preparation/) to produce your .bbl
% file and copy the contents into your paper here.
%
% Follow these steps:
% 1. Run LaTeX on your LaTeX file.
%
% 2. Make sure the bibliography style appears as \bibliographystyle{agufull08}. Run BiBTeX on your LaTeX
% file.
%
% 3. Open the new .bbl file containing the reference list and
%   copy all the contents into your LaTeX file here.
%
% 4. Comment out the old \bibliographystyle and \bibliography commands.
%
% 5. Run LaTeX on your new file before submitting.
%
% AGU does not want a .bib or a .bbl file. Please copy in the contents of your .bbl file here.

\begin{thebibliography}{}

\providecommand{\natexlab}[1]{#1}
\ifx\csname urlstyle\endcsname\relax
  \providecommand{\doi}[1]{doi:\discretionary{}{}{}#1}\else
  \providecommand{\doi}{doi:\discretionary{}{}{}\begingroup
  \urlstyle{rm}\Url}\fi

\bibitem[{\textit{Acton}(1996)}]{acton1996pss}
C.H. Acton
\newblock Ancillary data services of NASA's navigation and ancillary information facility.
\newblock {\em Planet.SpaceSci}, 44(65--70), 1996.

\bibitem[{\textit{Broiles et al.}(2015)}]{broiles2015aa}
T.W. Broiles, J.L. Burch, G.B. Clark, C. Koenders, E. Behar, R. Goldstein, S.A. Fuselier, K.E. Mandt, P. Mokashi, and M. Samara
\newblock Rosetta observations of solar wind interaction with the comet 67P/Churyumov--Gerasimenko.
\newblock {\em Astronomy \& Astrophysics}, 2015.

\bibitem[{\textit{Burch et al.}(2007)}]{burch2007ssr}
J.L. {Burch}, R.~{Goldstein}, T.E. {Cravens}, W.C. {Gibson}, R.N. {Lundin},
  C.J. {Pollock}, {Winningham} J.D., and D.T. {Young}.
\newblock Rpc-ies: the ion and electron sensor of the rosetta plasma
  consortium.
\newblock {\em Space Science Reviews}, 128:697--712, 2007.

\bibitem[{\textit{Coates et al.}(1986)}]{coates1986jgr}
A.~J. Coates, A.~D. Johnstone, M.~P. Smith, and D.~J. Rodgers.
\newblock Ampte/uks ion experiment observations of lithium releases in the
  solar wind.
\newblock {\em Journal of Geophysical Research}, 91(A2):1311--1319, February
  1986.

\bibitem[{\textit{Coates}(1997)}]{coates1997asr}
A.J. Coates
\newblock Ionosphere and magnetospheres of comets.
\newblock {\em Adv. Space Res.}, 20, 255-266, 1997.

\bibitem[{\textit{Coates}(2015)}]{coates2015jpcs}
A.J., Coates, J.L. Burch, R. Goldstein, H. Nilsson, G. Stenberg Wieser, E. Behar and the RPC team
\newblock Ion pickup observed at comet 67P with the Rosetta Plasma Consortium (RPC) particle sensors: similarities with previous observations and AMPTE releases, and effects of increasing activity
\newblock Proc. 14th Annual International Astrophysics Conference on Linear and Nonlinear Particle Energization throughout the Heliosphere and Beyond, Tampa (20-24 April 2015), Open Access Journal of Physics: Conference Series (JPCS), in press, Jun 2015.

\bibitem[{\textit{Glassmeier et al.}(2007a)}]{glassmeier2007ssr}
K.~H. Glassmeier, H. Boehnhardt, D. Koschny, E. K{\"u}hrt,
  and I. Richter.
\newblock The rosetta mission: Flying towards the origin of the solar system.
\newblock {\em Space Science Reviews}, 128(1-4):1--21, 2007.

\bibitem[{\textit{Glassmeier et al.}(2007b)}]{glassmeier2007ssr-MAG}
K.~H. Glassmeier, I.~Richter, A.~Diedrich, G.~Musmann, U.~Auster, and U.~et~al
  Motschmann.
\newblock Rpc-mag: The fluxgate magnetometer in the rosetta plasma consortium.
\newblock {\em Space Science Reviews}, 128:649--670, 2007.

\bibitem[{\textit{Goldstein et al.}(2015)}]{goldstein2015grl}
R.~{Goldstein}, J.L. {Burch}, P.~{Mokashi}, T.~{Broiles}, K.~{Mandt},
  T.~{Hanley}, T.~{Cravens}, A.~{Rahmati}, M.~{Samara}, G.~{Clark},
  M.~{H{\"a}ssig}, and J.M. {Webster}.
\newblock The rosetta ion and electron sensor (ies) measurement of the
  development of pickup ions from comet 67p/churyumov-gerasimenko.
\newblock {\em Geophysical Research Letters}, 42(9):3093--3099, May 2015.

\bibitem[{\textit{Haerendel et al.}(1986)}]{haerendel1986nature}
G.~{Haerendel}, G.~{Paschmann}, W.~{Baumjohann}, and C.W. {Carlson}.
\newblock Dynamics of the AMPTE artificial comet.
\newblock {\em Nature}, 320:720--723, April 1986.

\bibitem[{\textit{Hansen et al.}(2007)}]{hansen2007ssr}
K.~C. Hansen, T.~Bagdonat, U.~Motschmann, C.~Alexander, M.~R. Combi, T.~E.
  Cravens, T.~I. Gombosi, Y.-D. Jia, and I.~P. Robertson.
\newblock The plasma environment of comet 67p/churyumov-gerasimenko throughout
  the rosetta main mission.
\newblock {\em Space Science Reviews}, 2007.

\bibitem[{\textit{Haser}(1957)}]{haser1957}
L. Haser
\newblock Distribution d'intensit\'e dans la t\^ete d'une com\`ete.
\newblock {\em Acad\'emie Royale de Belgique - Bulletin de la classe des sciences}, 5\`eme s\'erie, Tome XLIII, 1957.

\bibitem{luhmann1985asr}
J.~G. Luhmann, C.~T. Russell, J.~R. Spreiter, and S.~S. Stahara.
\newblock Evidencer for mass-loading of the Venus magnetosheath.
\newblock {\em Adv. Space Res.}, 5(4):307--311, 1985.

\bibitem[{\textit{Luhman et al.}(1997)}]{luhman1997}
 J. G. Luhmann, S. H. Brecht,  J. R. Spreiter, S. S. Stahara, R. S. Steinolfson,  and A. F. Nagy.
\newblock In S. W. Bougher, D. M. Hunten, and R. J. Phillips (eds), 'Global Models of the Solar Wind Interaction with Venus', in Venus II -- Geology, Geophysics, Atmosphere, and Solar Wind Environment, University of Arizona Press, Tucson, p. 33. 1997.

\bibitem[{\textit{Neugebauer}(1990)}]{neugebauer1990rg}
M. Neugebauer.
\newblock Spacecraft Observations of the Interaction of Active Comets with the Solar Wind.
\newblock {\em Rev. Geophys.},  28, 231, 1990.

\bibitem[{\textit{Nilsson et al.}(2007)}]{nilsson2007ssr}
H.~Nilsson, R.~Lundin, K.~Lundin, S.~Barabash, H.~Borg, O.~Norberg, A.~Fedorov,
  J.-A Sauvaud, H.~Koskinen, E.~Kallio, P.~Riihel{\"a}, and J.L. Burch.
\newblock Rpc-ica: The ion composition analyzer of the rosetta plasma consortium.
\newblock {\em Space Science Reviews}, 128(1-4):671--695, 2007.

\bibitem[{\textit{Nilsson et al.}(2015a)}]{nilsson2015science}
H. Nilsson, G. Stenberg~Wieser, E. Behar, C.~S. Wedlund,
  H. Gunell, M. Yamauchi, R. Lundin, S. Barabash, M.
  Wieser, C. Carr, E. Cupido, J.~L. Burch, A. Fedorov,
  J.~A. Sauvaud, H. Koskinen, E. Kallio, J.~P. Lebreton,
  A. Eriksson, N. Edberg, R. Goldstein, P. Henri, C.
  Koenders, P. Mokashi, Z. Nemeth, I. Richter, K. Szego, M.
  Volwerk, C. Vallat, and M. Rubin.
\newblock Birth of a comet magnetosphere: A spring of water ions.
\newblock {\em Science}, 347(6220), 2015.

\bibitem[{\textit{Nilsson et al.}(2015b)}]{nilsson2015aa}
H.~Nilsson, G.~Stenberg-Wieser, and E.~Behar.
\newblock Evolution of the ion environment of comet 67p/churyumov-gerasimenko -
  observations between 3.6 and 2.2~AU.
\newblock {\em Astronomy \& Astrophysics}, 2015.

\bibitem[{\textit{Rodgers et al.}(1986)}]{rodgers1986nature}
D.J. {Rodgers}, A.J. {Coates}, A.D. {Johnstone}, M.P. {Smith}, D.A. {Bryant},
  D.S. {Hall}, and C.P. {Chaloner}.
\newblock Uks plasma measurements near the ampte artificial comet.
\newblock {\em Nature}, 320:712--716, April 1986.

\bibitem[{\textit{Szeg{\"o} et al.}(2000)}]{szego2000ssr}
K.~{Szeg{\"o}}, K.-H. {Glassmeier}, R.~{Bingham}, A.~{Bogdanov}, C.~{Fischer},
  G.~{Haerendel}, A.~{Brinca}, T.~{Cravens}, E.~{Dubinin}, K.~{Sauer},
  L.~{Fisk}, T.~{Gombosi}, N.~{Schwadron}, P.~{Isenberg}, M.~{Lee},
  C.~{Mazelle}, E.~{M{\"o}bius}, U.~{Motschmann}, V.~D. {Shapiro},
  B.~{Tsurutani}, and G.~{Zank}.
\newblock {Physics of Mass Loaded Plasmas}.
\newblock {\em Space Sci. Rev.}, 94:429--671, December 2000.

\end{thebibliography}

%Reference citation examples:

%...as shown by \textit{Kilby} [2008].
%...as shown by {\textit  {Lewin}} [1976], {\textit  {Carson}} [1986], {\textit  {Bartholdy and Billi}} [2002], and {\textit  {Rinaldi}} [2003].
%...has been shown [\textit{Kilby et al.}, 2008].
%...has been shown [{\textit  {Lewin}}, 1976; {\textit  {Carson}}, 1986; {\textit  {Bartholdy and Billi}}, 2002; {\textit  {Rinaldi}}, 2003].
%...has been shown [e.g., {\textit  {Lewin}}, 1976; {\textit  {Carson}}, 1986; {\textit  {Bartholdy and Billi}}, 2002; {\textit  {Rinaldi}}, 2003].

%...as shown by \citet{jskilby}.
%...as shown by \citet{lewin76}, \citet{carson86}, \citet{bartoldy02}, and \citet{rinaldi03}.
%...has been shown \citep{jskilbye}.
%...has been shown \citep{lewin76,carson86,bartoldy02,rinaldi03}.
%...has been shown \citep [e.g.,][]{lewin76,carson86,bartoldy02,rinaldi03}.
%
% Please use ONLY \citet and \citep for reference citations.
% DO NOT use other cite commands (e.g., \cite, \citeyear, \nocite, \citealp, etc.).

%% ------------------------------------------------------------------------ %%
%
%  END ARTICLE
%
%% ------------------------------------------------------------------------ %%
\end{article}
```


abstract: We experimentally show octave-spanning supercontinuum generation in a non-stoichiometric silicon-rich nitride waveguide when pumped by femtosecond pulses from an erbium fiber laser. The pulse energy and bandwidth are comparable to results achieved in stoichiometric silicon nitride waveguides, but our material platform is simpler to manufacture. We also observe wave-breaking supercontinuum generation by using orthogonal pumping in the same waveguide. Additional analysis reveals that the waveguide height is a powerful tuning parameter for generating mid-infrared dispersive waves while keeping the pump in the telecom band.
author: Xing Liu$^1$, Minhao Pu$^1$, Binbin Zhou$^1$, Clemens J. Krückel$^2$, Attila Fülöp$^2$, Victor Torres-Company$^{2,*}$, and Morten Bache$^{1,\dagger}$,
title: Octave-spanning supercontinuum generation in a silicon-rich nitride waveguide

Supercontinuum generation in photonic crystal fibers pumped with a mode-locked laser oscillator allowed generating octave-broadened coherent spectra without an amplifying stage . This remarkable achievement paved way for self-referencing in laser frequency combs and key advances in biophotonics .

There is a growing interest in supercontinuum generation in photonic integrated waveguides. This offers the prospect of supercontinuum generation in different regions of the electromagnetic spectrum with extremely low pump pulse energies delivered by near-IR ultrafast femtosecond fiber lasers. One of the most promising waveguide platforms that has demonstrated near-IR-pumped octave-spanning supercontinua is silicon nitride , in which waveguides can be manufactured directly on an oxidized silicon wafer, and which, unlike silicon, has negligible two-photon absorption in the near IR. The nonlinear Kerr coefficient is comparatively lower than other integrated waveguide technologies, but this is offset by lower linear propagation loss. As a result, large nonlinear phase shifts can be obtained at modest pump powers. The broadest supercontinuum generated on chip has been measured on a silicon nitride waveguide pumped with an Yb mode-locked fiber laser . Carefully engineered silicon nitride waveguides provided highly coherent supercontinua , suitable for detecting the carrier envelope offset frequency of frequency combs .

Thick waveguides are needed to achieve suitable dispersion engineering and high mode confinement. However, silicon nitride films tend to crack above 400 nm thickness. Advanced fabrication methods utilize crack barriers and in-trench growth to overcome film cracking which recently enabled the supercontinua reported in . An octave-spanning spectrum in non-stoichiometric silicon nitride waveguides was recently reported but it required a complex multi-cladding structure for dispersion engineering . A simpler approach to circumvent cracking is to reduce film stress by slightly increasing the silicon content during deposition. This allows thick crack-free film growth in a single deposition step with minimized process complexity. The resulting so-called silicon-rich nitride shows a higher refractive index and an increased material absorption loss than stoichiometric silicon nitride, but this is compensated for by a higher nonlinear Kerr coefficient .

Here we report octave-spanning supercontinuum generation in a 10 mm long silicon-rich nitride waveguide pumped in the quasi-TE mode with a mode-locked femtosecond Er-fiber laser in the telecom C band. The key waveguide aspects are that an increased nonlinearity compensates for the higher propagation loss, and that thick waveguides can be made with a simple, single-cladding, CMOS-compatible structure suitable for dispersion engineering at telecom wavelengths. This gives a supercontinuum performance similar to stoichiometric silicon nitride waveguides made with advanced fabrication techniques . When pumping in the quasi-TM mode, spectral broadening from optical wave breaking is observed for the first time in silicon nitride waveguides, which we show is caused by the TM pump mode experiencing normal dispersion.

Our silicon-rich nitride film is fabricated in a low-pressure chemical vapor deposition (CVD) process . The waveguide patterns are transferred by standard deep UV contact lithography followed by dry etching. The silica top cladding is done in a plasma-enhanced CVD process. The waveguides feature a Kerr parameter of $\gamma=5.7
\pm 0.5~{\rm  (Wm)}^{-1}$ around $1.55 ~\mu {\rm m}$, which is a slightly modified value to that of owing to an improved linear propagation loss estimate ($1.35\pm0.3$ dB/cm). Importantly, two-photon absorption is absent at $1.55 ~\mu {\rm m}$. The waveguide dimensions are displayed in Fig. (a). The waveguide height was around 700 nm, with about $5\%$ variation across the wafer. This uncertainty translates into a range of group-velocity dispersion (GVD) profiles bounded by the shadow areas in Fig. (a). The dispersion is calculated using finite-element simulations (COMSOL) and includes both the dispersion of the material in bulk and the waveguide's geometry. The dispersion of the film is measured up to $1.7 ~\mu {\rm m}$ and extrapolated up to $2.4 ~\mu {\rm m}$ using empirical relations from the bulk material dispersion. The slight tilt in the waveguide's wall is a consequence of the etching process. The COMSOL calculations include this tilt, and showed that tilt variations have a negligible effect in the GVD compared to the uncertainty in the waveguide height or material dispersion. For the TE mode, the resulting GVD has two zero-dispersion wavelengths (ZDWs) giving a broad region (spanning 700 nm) with flat and moderate anomalous dispersion, centered around 1550 nm, as desired for efficient supercontinuum generation.

The experimental setup is shown in Fig. (b). (Please see acronyms in the caption). An Er femtosecond fiber laser centered at 1555 nm produces a 90 MHz pulse train; each pulse has a 33 nm bandwidth (FWHM) supporting a transform-limited 105 fs (FWHM) Gaussian pulse. The pulses were coupled by an objective lens into a DCF to compensate for the accumulated chirp before the waveguide. NDFs were used to control the pulse energy, and an FPC controlled the polarization, followed by a piece of single-mode fiber spliced to a tapered lensed fiber (OZ Optics, based on SMF-28), which focused to a $2.0~\mu {\rm m}$ FWHM spot size. We used an intensity autocorrelator to estimate the pulse duration to 130 fs FWHM (assuming a Gaussian shape) at the waveguide entrance, revealing some remaining chirp in the pulse. The maximum total power before the waveguide was 42 mW. The output from the waveguide was collected with a tapered lensed fiber. To record the supercontinuum, we used two OSAs with spectral ranges 600-1700 nm and 1200-2400 nm, whose spectra were overlapped to get the final spectrum. The coupling loss is estimated at $6.5\pm1.0$ dB/facet (TE) and $5.3\pm1.0$ dB/facet (TM), providing a maximum input pulse energy of 105 pJ (TE) and 140 pJ (TM) coupled into the waveguide.

The measured power spectral density (PSD) in the TE case is shown in Fig. (a) for a range of pulse energies. The spectra for low pulse energies are clearly dominated by self-phase modulation (SPM), i.e., early stage broadening before soliton formation occurs. At 78 pJ, the soliton has formed, accompanied by two soliton-induced dispersive waves , one on each side of the two ZDWs. The strongest dispersive wave is found at low wavelengths, peaked around 820 nm, while the mid-IR dispersive wave, peaked around 2250 nm, is less powerful and changes significantly when increasing the pulse energy further to the maximum value of 105 pJ (9 mW average power).

The results are verified with numerical simulations using the so-called nonlinear analytic envelope equation (NAEE) , which resolves subcycle carrier-wave dynamics and includes a full expansion of the cubic nonlinearity. We only model the waveguide fundamental mode (FM). The mode effective index and effective area from the COMSOL calculations were extrapolated directly to the NAEE grid without polynomial expansions and were carefully extended beyond the COMSOL modeling domain to avoid unphysical dispersion scenarios. The frequency dependence of the effective mode area was modeled as shown in . A constant propagation loss was used, and noise was included as a half photon per temporal mode in the input, modeling quantum vacuum fluctuations in the Wigner representation . An extended NAEE model was used to include the Raman effect , which empirically included the Si-N asymmetric stretching mode centered at 410 $\rm cm^{-1}$ , having a broadband 70 $\rm cm^{-1}$ linewidth and a relative Raman strength $f_R=0.2$, both typical values for amorphous materials. Generally, the Raman effect gave minor contributions, in line with previous studies of silicon nitride waveguides. This might be due to the moderate anomalous dispersion range of the waveguide in which the Raman effect can influence the soliton by redshifting it.

Figure (b) shows corresponding numerical simulations where the major features of the experimental spectra are reproduced in the power sweep, except for the precise location of the dispersive waves. This is likely due to the uncertainty of the mid-IR dispersion and to the precise soliton wavelength (i.e., how much it is frequency-shifted by Raman or self-steepening effects). The dispersive waves also seem to be narrower in the experimental case. The simulations only model the FM, but from the higher-order modes (HOMs) given by the COMSOL simulations we estimate that around 55% of the input energy will be coupled to the FM, and the rest goes into the HOMs. The HOMs will not show any spectral broadening due to their lower peak powers, higher GVD, and larger mode areas. In the numerical spectra shown in Fig. (b) we have therefore added the equivalent energy of the HOMs at the end of the simulation as 45% of the total input spectrum. Additionally, combining the 33 nm bandwidth of the input pulse with the measured value of 130 fs FWHM corresponds to a group-delay dispersion of ${\rm GDD}=\pm 3500~{\rm fs}^2$. The simulations used positive chirp but similar results were found with negative chirp. We note that the simulations used lower energies than the experiments, which can be attributed to uncertainties in coupling efficiencies, waveguide nonlinearity, dispersion, input pulse pre-chirp etc.

Figure shows that the waveguide height is a powerful dispersion tuning parameter. At the pump wavelength, the TE mode has anomalous dispersion for all considered heights, while the GVD for the TM mode changes sign at $h=675$ nm. The simulations show the best match to the experimental results for $h=660$ nm. With this parameter, the long- and short-wavelength dispersive waves (DW 1 and 2) agree quite well with the experiments. This height is within the 5% variation margin across the wafer as mentioned above. Similar results were observed for $h=$`<!-- -->`{=html}650$-$`<!-- -->`{=html}675 nm, while, beyond this range, the simulations did not match the experiments very well.

Figure shows a representative simulation for the supercontinuum propagation dynamics. The pulse energy is 40 pJ, corresponding to the green curve in Fig. (b). In Fig. (a) the soliton forms as a single-cycle spike (5 fs FWHM), and, after this self-compression point ("max"), the temporal trace shows interference oscillations. These are due to the two dispersive waves phase-matched to the soliton, all temporally overlapping, but located at different wavelengths. In the spectral evolution (e) the dispersive waves are seen to be present at the soliton formation point ("max"), and, after this, they grow significantly; see also (d). In (c), the shot-to-shot coherence is shown, as calculated from the first order coherence function $g_{12}^{(1)}(\lambda)$ . This excellent coherence pertains for lower powers, while we found it degrades for higher levels (which is quite typical).

Figure (a) shows the TM case. While less broadening was observed, the spectrum broadens sufficiently to enter the important 1600-1850 nm wavelength range for 3-photon absorption microscopy . The spectral shapes are indicative of optical wave breaking , which happens when the pump dispersion is normal, so no soliton can form. The spectral broadening is typically much weaker than in the soliton case because no dispersive waves are formed and the pulse does not self compress. Instead, the pulse becomes highly chirped and eventually develops steep temporal shock fronts. In frequency domain, this is accompanied by highly coherent SPM-induced spectral broadening. To understand when TM has normal GVD at 1555 nm, we refer to Fig. , which shows that the TM mode has normal GVD at 1555 nm when the waveguide height is taken well below 700 nm. Indeed, the numerical simulations shown in Fig. (c) agree well with the experimental results using the exact same parameters as in the TE case concerning waveguide size and pulse pre-chirp. In time domain (not shown), the pulses were chirped, but the powers were too low to observe severe shock front formation. The experimental data show extremely symmetrical SPM broadening, which could only be observed in the simulations when using positively chirped input pulse, while a negative chirp tended to give more asymmetric spectra. Finally, Fig. (b) shows the coherence of the simulated TM continua, demonstrating perfect coherence across the generated bandwidth, as is typical for optical wave breaking.

Our simulations show that for a transform-limited pump, the results in Fig. (a) could be obtained using substantially lower pump energies. Additionally, the long pump wavelength of Er-fiber lasers allows generating a mid-IR dispersive wave in the TE case for thick waveguides (Fig. ). Combining these efforts, Fig. shows the promising potential of these waveguides: with a $h=750$ nm thick waveguide and a TE transform-limited pump with only 20 pJ energy, over 2 octaves of coherent supercontinua are generated, extending well into the mid-IR.

Summarizing, we studied supercontinuum generation in a silicon-rich nitride waveguide pumped with low-energy femtosecond pulses (80-140 pJ) from an Er fiber laser. In the TE case, a soliton and two dispersive waves were excited to give a 1.5 octave supercontinuum (820-2250 nm at $-30$ dB). In the TM case, a continuum was generated by optical wave breaking as the TM mode had normal dispersion at the pump wavelength. Numerical results indicate an excellent coherence of the supercontinua and that they could extend into the mid-IR with a slightly thicker waveguide. These results are promising for short-range near-IR and mid-IR coherent supercontinuum generation in an easily manufacturable CMOS-compatible waveguide.

**Funding.** Vetenskapsrådet; Teknologi og Produktion, Det Frie Forskningsråd (FTP, DFF) (11-106702); European Research Council (ERC) (ERC-AdvG 291618).


abstract: In our paper, we construct a real-algebraic function on some closed manifold whose *Reeb* (*Kronrod-Reeb*) *graph* is a graph respecting some *algebraic domain*: a graph for this is called a *Poincaré-Reeb* graph.
  .
  The *Reeb graph* of a smooth function is defined as a natural graph which is the quotient space of the manifold of the domain under a natural equivalence relation for some wide and nice class of smooth functions. The vertex set is defined as the set of all connected components containing some singular points of the function: a *singular* point of a smooth function is a point where the differential vanishes. Morse-Bott functions give very specific cases. The relation is to contract each connected component of each preimage to a point.
  .
  Sharko has asked a natural and important problem: can we construct a nice smooth function whose Reeb graph is a given graph? Explicit answers have been given first by Masumoto-Saeki in a generalized manner for closed surfaces. After that various answers have been presented by various researchers and most of them are essentially for functions on closed surfaces and Morse functions such that connected components of preimages containing no singular points are spheres. Recently the author has also considered questions and solved in cases the preimages are general manifolds.
address: Institute of Mathematics for Industry, Kyushu University, 744 Motooka, Nishi-ku Fukuoka 819-0395, Japan
  TEL (Office): +81-92-802-4402
  FAX (Office): +81-92-802-4405
author: Naoki kitazawa
title: Real algebraic functions on closed manifolds whose Reeb spaces are given graphs

# Introduction

The *Reeb graph* (*Kronrod-Reeb graph*) of a differentiable function $c:X \rightarrow Y$ is a graph whose underlying topology is the quotient space of the manifold of the domain defined in the following way: two points $x_1, x_2 \in X$ are equivalent if and only if they are in a same connected component of some preimage $c^{-1}(y)$.

Let $W_c$ denote this. We define the quotient map $q_c:X \rightarrow W_c$ and we can define a map $\bar{c}$ enjoying the relation $c=\bar{c} \circ q_c$ uniquely.

This is regarded as a graph for some wide class of smooth functions. We explain about the structure of our Reeb graph for $W_c$.

Such topological and combinatorial objects have been fundamental and strong tools in understanding the manifolds. is a related pioneering paper.

We define fundamental terminologies, notions, and notation.

For a topological space $X$ having the structure of some cell complex whose dimension is finite, its dimension is uniquely defined as the dimension of the cell complex. Let $\dim X$ denote this.

A polyhedron and a CW complex is of course of such a class and a topological manifold is of such a class, being known to have the structure of a CW complex of a finite dimension.

A graph is a $1$-dimensional CW complex where the vertex set is the set of all $0$-dimensional cells and the edge set is the set of all $1$-dimensional cells. If the closure of an edge is homeomorphic to a circle, then it is called a *loop*. Hereafter, we do not consider such graphs and a graph is always a $1$-dimensional polyhedron. An *isomorphism* from a graph $K_1$ onto $K_2$ means a piecewise smooth homeomorphism mapping the edge set and the vertex set of $K_1$ onto those of $K_2$.

For a differentiable map $c:X \rightarrow Y$, a *singular* point $x \in X$ is a point where the differential is smaller than both $\dim X$ and $\dim Y$. The value at some singular point of $c$ is a *singular value* of $c$. If a smooth function on a closed manifold has finitely many singular values, then the Reeb space of it is a graph where the vertex set consists of all points $p$ whose preimages ${q_c}^{-1}(p)$ contain some singular points of $c$. Morse-Bott functions and smooth functions of some considerably wide classes satisfy this. This is due to .

has asked this first. Smooth functions on closed surfaces have been explicitly constructed for graphs satisfying some nice conditions there. generalizes this for arbitrary graphs. Later, have given answers. These study cases for smooth functions on closed surfaces and Morse functions such that each connected component of each preimage containing no singular points is always a sphere for example. In , the author has studied cases where such connected components are general manifolds with mild conditions on singularities of the functions for example. appears as a paper motivated by and through related informal discussions by us.

The category we discuss has been always the smooth one. In this paper, we try Problem and obtain related answers under some specific or explicit situations.

To present our main result, we need some terminologies, notions and notation from the theory of smooth manifolds and (real) algebraic manifolds for example.

${\mathbb{R}}^k$ is the $k$-dimensional Euclidean space and this is a simplest smooth manifold of dimension $k$ for an arbitrary integer $k>0$. This is also a simplest real algebraic manifold, which is also called the ($k$-dimensional) *real affine space*. It is also a Riemannian manifold equipped with the so-called Euclidean metric. For a point $x \in {\mathbb{R}}^k$, $||x|| \geq 0$ is the distance between $x$ and the original $0 \in {\mathbb{R}}^k$. $\mathbb{R}:={\mathbb{R}}^1$. $S^k:=\{x \in {\mathbb{R}}^{k+1} \mid ||x||=1 \}$ is the $k$-dimensional unit sphere. This is a $k$-dimensional algebraic submanifold of ${\mathbb{R}}^{k+1}$ which is compact and has no boundary. It is also a smooth submanifold. It is connected for $k \geq 1$ and it is a discrete set with exactly two points for $k=0$.

An *algebraic domain* $D$ of ${\mathbb{R}}^k$ is some open subset there such that the boundary of the closure $\overline{D}$ consists of finitely many smooth algebraic submanifolds of dimension $k-1$ or smooth algebraic hypersurfaces with no boundaries.

To simplify our arguments, let us assume the following where $l \geq 0$ is an integer.

-   For each hypersurface $S_j$ in the family $\{S_j\}_{j=1}^l$, a real polynomial $f_{{\rm P}, S_j}$ is given so that the zero set and $S_j$ coincide and that the polynomial function $f_{{\rm P}, S_j}:S_j \rightarrow \mathbb{R}$ defined canonically has no singular points on $S_j$.

-   $D$ is assumed to be the intersection ${\bigcap}_{j=1}^l  \{x \in {\mathbb{R}}^k \mid  f_{{\rm P}, S_j}(x)>0 \}$.

For example, the interior ${\rm Int}\ D^k$ of $D^k:=\{x \in {\mathbb{R}}^{k} \mid ||x|| \leq 1 \}$ is a simplest example and $D^k$ is the $k$-dimensional unit disk. This is also a $k$-dimensional smooth, compact and connected submanifold. Note that $||x||={\Sigma}_{j=1}^k {x_j}^2$ where $x:=(x_1,\cdots,x_k)$.

A *Poincaré-Reeb graph* is defined for a pair of an algebraic domain $D$ of the real affine space of dimension $k>1$ and a canonical projection ${\pi}_{k,1}$ mapping $(x_1,x_2) \in {\mathbb{R}}^{k}$ to $x_1 \in {\mathbb{R}}$. This can be presented in a more general manner. Hereafter, we mainly respect the preprint and there such cases are discussed. Note that terminologies and situations are different in considerable cases and that here we can argue in a self-contained way.

See also for related theory for example. We present our main result. In the following section we prove this and present related comments as our main content.

The author is a member of the project JSPS KAKENHI Grant Number JP22K18267 "Visualizing twists in data through monodromy" (Principal Investigator: Osamu Saeki). Our present study is supported by the project.

**Data availability.**
Data essentially supporting our present study are all contained in our present paper.

# On Main Theorem .

*A proof of Main Theorem .* Let $k_0$ be an arbitrary integer satisfying $k_0>k$ as presented.

Let us use $x:=(x_1,\cdots,x_k)$ for (local) coordinates for ${\mathbb{R}}^k$.

Let us use $y:=(y_1,\cdots,y_{k^{\prime}})$ for (local) coordinates for ${\mathbb{R}}^{k^{\prime}}$ where $k^{\prime}:=k_0-k$.

We take two steps to complete the proof.

STEP 1 Defining a set in $M_D \subset {\mathbb{R}}^{k_0}$, which is a real algebraic hypersurface and a smooth regular compact submanifold of dimension $k_0-1$ with no boundary.
First we define $M_{D_0}:= \{(x,y) \in {\mathbb{R}}^k \times {\mathbb{R}}^{k^{\prime}}={\mathbb{R}}^{k_0} \mid  {\prod}_{j=1}^l (f_{{\rm P}, S_j}(x))-||y||^2=0\}$.
For the notation, remember the rule $||y||={\Sigma}_{j=1}^l {y_j}^2$.

We show this is also a smooth regular submanifold in ${\mathbb{R}}^{k_0}$. We consider the partial derivative of the function ${\prod}_{j=1}^l  (f_{{\rm P}, S_j}(x))-{\Sigma}_{j=1}^l {y_j}^2$ for variants $x_j$ and $y_j$. First we take a point $(x_0,y_0) \in M_{D_0}$ such that ${\prod}_{j=1}^l (f_{{\rm P}, S_j}(x_0))>0$. We use $x_0:=(x_{0,1},\cdots,x_{0,k})$ and $y_0:=(y_{0,1},\cdots,y_{0,k})$. Here we consider the partial derivative of the function for some $y_j$ and we have the value $2{y_j}=2y_{0,j} \neq 0$. The differential of the restriction of the function ${\prod}_{j=1}^l  (f_{{\rm P}, S_j}(x))-{\Sigma}_{j=1}^l {y_j}^2$ at $(x_0,y_0) \in M_{D_0}$ is not of rank $0$ and this is not a singular point of the function.

Second we take a point $(x_{S_a},y_{S_a}) \in M_{D_0}$ such that $f_{{\rm P}, S_a}(x_{S_a})=0$. By the assumption on the hypersurfaces $S_b$ and the polynomials $f_{{\rm P}, S_j}(x_{S_b})$, $f_{{\rm P}, S_{a^{\prime}}}(x_{S_a})> 0$ for $a^{\prime} \neq a$. The polynomial function defined canonically from the polynomial $f_{{\rm P}, S_{a}}$ is assumed to have no singular points on $S_{a}$. We use $x_{S_a}:=(x_{S_a,1},\cdots,x_{S_a,k})$ and $y_{S_a}:=(y_{S_a,1},\cdots,y_{S_a,k})$. Here we consider the partial derivative of the function for some $x_{j}$ and we have the non-zero value represented as the product of the partial derivative of the function $f_{{\rm P}, S_a}(x)$ for $x_{j}$ at $(x_{S_a},y_{S_a})$ and the product of $l-1$ numbers defined as the values of polynomials (or the canonically defined polynomial functions) in the family $\{f_{{\rm P}, S_j}\}_{j=1}^l$ at $x_{S_a}$ except the number $j \neq a$. The differential of the restriction of the function at $(x_{S_a},y_{S_a}) \in M_{D_0}$ is not of rank $0$ and this is not a singular point of the function.

We have shown that $M_{D_0}$ is a smooth regular submanifold by the implicit function theorem.

We define $M_{D}$ as the set of all points in $M_{D_0}$ such that $x \in \overline{D} \supset D$. We investigate a small neighborhood of each point in $M_{D}$.

First we consider a point $p_1 \in D$ and a point $(p_1,q_1) \in M_{D}$ and take its sufficiently small open neighborhood $U_{p_1,q_1}$ in ${\mathbb{R}}^{k_0}$. For any point in $M_{D_0} \bigcap U_{p_1,q_1}$, by the definition, it is also a point in $M_{D}$. Second we consider a point $p_2 \in \partial \overline{D}$ in the boundary $\partial \overline{D} \subset \overline{D}$ and a point $(p_2,q_2) \in M_{D}$ and take its sufficiently small open neighborhood $U_{p_2,q_2}$ in ${\mathbb{R}}^{k_0}$. Take an arbitrary point $(p^{\prime},q^{\prime})$ in $M_{D_0} \bigcap U_{p_2,q_2}$. By the definition and the assumption on the hypersurfaces $S_b$ and the polynomials $f_{{\rm P}, S_j}(x_{S_b})$, we can know $f_{{\rm P}, S_{b^{\prime}}}(p^{\prime})>0$ for $1 \leq b^{\prime} \leq l$ except for one $b^{\prime}:={b_0}^{\prime}$. $f_{{\rm P}, S_{{b_0}^{\prime}}}(p^{\prime})<0$ cannot occur by the form of the function ${\prod}_{j=1}^l  (f_{{\rm P}, S_j}(x))-{\Sigma}_{j=1}^l {y_j}^2$. $(p^{\prime},q^{\prime})$ is also a point in $M_{D}$.

We have shown that $M_{D}$ is also a smooth regular submanifold with $M_{D}=M_{D_0}$. By the form of the function ${\prod}_{j=1}^l (f_{{\rm P}, S_j}(x))-{\Sigma}_{j=1}^l {y_j}^2$ and the compactness of the closure $\overline{D}$, it is also a smooth compact manifold with no boundary.

STEP 2 The composition of the restriction of the canonical projection mapping $(x_1,x_2) \in {\mathbb{R}}^{k_0}$ to $x_1 \in {\mathbb{R}}^k$ to the submanifold $M_D$ with the restriction of the given projection ${\pi}_{k,1}$ or ${\pi}_D$ in Definition .
First restrict the canonical projection mapping $(x_1,x_2) \in {\mathbb{R}}^{k_0}$ to $x_1 \in {\mathbb{R}}^k$ to the submanifold $M_D$. We thus have a surjection onto $\overline{D}$. We restrict this to the preimage of $D$. This is, by the form of the function, regarded as a projection and a submersion. If we restrict this to the preimage of the boundary $\partial \overline{D}$, then, by the form of the function, we have a smooth and real algebraic embedding onto $\partial \overline{D}$. We compose the surjection onto $\overline{D}$ with the ${\pi}_D$ to have a new real algebraic function.

By our definitions and situations, we can see that the composition obtained before is regarded as a desired function for $k_0>k+1$ where we need to respect connectedness of the preimages.

We note some. Some may help us to understand our arguments more rigorously. First, our resulting function is a function for , having finitely many singular values. Second, our map on $M_D$ onto $\overline{D}$ is, topologically, regarded as a so-called *special generic* map. The class of special generic maps contains Morse functions with exactly two singular points on spheres, or Morse functions in the so-called Reeb's theorem, and canonical projections of unit spheres. See for fundamental theory on special generic maps and some advanced results on manifolds admitting such maps. For construction of special generic maps related to our construction of the map on $M_D$ onto $\overline{D}$, consult also the preprints of the author for example.

This completes the proof.

. ◻

We present remarks related to our result.

We present another remark which is not directly related to our study in the present situation. We also wish this is closely related to our study in the future.

# Acknowledgement

The author would like to thank Osamu Saeki and Shuntaro Sakurai for interesting discussions on and construction of explicit smooth maps which are also real algebraic. This with an interesting talk in the conference "Singularity theory and its applications (RIMS-Sing 1)"
(http://www.math.kobe-u.ac.jp/HOME/saji/math/conf2022/spmon/index.html)
by M. S. Sorea on and her interesting answers to questions on construction of nice domains and extensions of the presented study and results to graphs of wider classes by the author has motivated us to start studies on our present paper.


= 20pt

**The Status of String Theory** [^1]

John H. Schwarz[^2]
*California Institute of Technology, Pasadena, CA 91125 USA*

**Abstract**

> There have been many remarkable developments in our understanding of superstring theory in the past few years, a period that has been described as "the second superstring revolution." Several of them are discussed here. The presentation is intended primarily for the benefit of nonexperts.

# Introduction

This manuscript presents a brief overview of some of the advances in understanding superstring theory that have been achieved in the last few years. It is aimed at physicists who are not experts in string theory, but who are interested in hearing about recent developments. Where possible, the references cite review papers rather than original sources.

It is now clear that what had been regarded as five distinct superstring theories in ten dimensions are better viewed as five special points in the moduli space of consistent vacua of a single theory. Morover, another special limit corresponds to a vacuum with Lorentz invariance in eleven dimensions. Some of the evidence that supports this picture is reviewed in Sect. 2. The "second superstring revolution" is characterized by the discovery of various non-perturbative properties of superstring theory. An important aspect of this is the occurrence of $p$-dimensional excitations, called $p$-branes. Their properties are under good mathematical control when they preserve some of the underlying supersymmetry. The maximally supersymmetric $p$-branes that occur in 10 or 11 dimensions are surveyed in Sect. 3. For detailed reviews of the material in Sect. 2 and Sect. 3, see – .

Sect. 4 describes how suitably constructed brane configurations can be used to derive, and make more geometrical, some of the non-perturbative properties of supersymmetric gauge theories that have emerged in recent years. Sect. 5 presents evidence for the existence of new non-gravitational quantum theories in six dimensions. In particular, there are pairs of theories with (2,0) and (1,1) supersymmetry that are related by T duality. Finally, in Sect. 6, the Matrix Theory proposal, which is a candidate for a non-perturbative description of M theory in a certain class of backgrounds, is sketched. This subject has been reviewed recently in . This is a rapidly developing subject, which appears likely to be a major focus of research in the next couple of years.

There have been other interesting developments in the past few years, which are omitted from this survey. The most remarkable, perhaps, is the application of D-brane technology to the study of black hole physics. This has led to a microscopic explanation of the origin of black hole thermodynamics in wide classes of examples. (For reviews see and .) Other omitted topics include applications to particle physics phenomenology and to cosmology. For two other surveys of recent developments in string theory, see and .

# M Theory

A schematic representation of the relationship between the five superstring vacua in 10d and the 11d vacuum, characterized by 11d supergravity at low energy, is given in Fig. 1. The idea is that there is some large moduli space of consistent vacua of a single underlying theory – denoted by M here. The six limiting points, represented as circles, are special in the sense that they are the ones with (super) Poincaré invariance in ten or eleven dimensions. The letters on the edges refer to the type of transformation relating a pair of limiting points. The numbers 16 or 32 refer to the number of unbroken supersymmetries. In 10d the minimal spinor is Majorana–Weyl and has 16 real components, so the conserved supercharges correspond to just one MW spinor in three cases (type I, HE, and HO). Type II superstrings have two MW supercharges, with opposite chirality in the IIA case and the same chirality in the IIB case. In 11d the minimal spinor is Majorana with 32 real components.

The 11d vacuum, including 11d supergravity, is characterized by a single scale – the 11d Planck scale $m_p$. It is proportional to $G_N^{-1/9}$, where $G_N$ is the 11d Newton constant. The connection to type IIA theory is obtained by taking one of the ten spatial dimensions to be a circle ($S^1$ in the diagram) of radius $R$. Type IIA string theory in 10d has a dimensionless coupling constant $g_s$, which is given by the vev of $e^\phi$, where $\phi$ is the dilaton field – a massless scalar field belonging to the IIA supergravity multiplet. In addition, the IIA theory has a mass scale, $m_s$, whose square gives the tension of the fundamental IIA string. The relationship between the parameters of the 11d and IIA descriptions is given by $$\label{ems}
m_s^2 = Rm_p^3$$ $$\label{gs}
g_s = Rm_s.$$ Numerical factors (such as $2\pi$) are not important for present purposes and have been dropped. The significance of these equations will emerge later. However, one point can be made immediately. The conventional perturbative analysis of the IIA theory is an expansion in powers of $g_s$ with $m_s$ fixed. The second relation implies that this is an expansion about $R = 0$, which accounts for the fact that the 11d interpretation was not evident in studies of perturbative string theory. The radius $R$ is a modulus – the vev of a massless scalar field with a flat potential. One gets from the IIA point to the 11d point by continuing this vev from zero to infinity. This is the meaning of the edge of Fig. 1 labeled $S^1$.

The relationship between the $E_8 \times E_8$ heterotic string vacuum (denoted HE) and 11d is very similar. The difference is that the compact spatial dimension is a line interval (denoted I in the diagram) instead of a circle. The same relations in eqs. () and () apply in this case. This compactification leads to an 11d space-time that is a slab with two parallel 10d faces. One set of $E_8$ gauge fields is confined to each face, whereas the gravitational fields reside in the bulk. There is a nice generalization of the 10d anomaly cancellation mechanism to this 11d setting . It only works for $E_8$ gauge groups.

The two edges of Fig. 1 labeled T connect vacua related by T duality. For example, if the IIA theory is compactified on (another) circle of radius $R_A$ leaving nine noncompact dimensions, this is equivalent to compactifying the IIB theory on a circle of radius $$R_B = (m_s^2 R_A)^{-1}.$$ Thus, continuing the modulus $R_A$ from infinity to zero (or $R_B$ from zero to infinity) gives an interpolation between the IIA and IIB theories. The T duality relating the two heterotic theories (HE and HO) is essentially the same, though there are additional technical details in this case.

The edge connecting the HO vacuum and the type I vacuum is labeled by S in the diagram, since these two vacua are related by S duality. Specifically, denoting the two string coupling constants by $g_s^{(HO)}$ and $g_s^{(I)}$, the relation is $$g_s^{(I)} g_s^{(HO)} = 1.$$ In other words, the two dilatons satisfy $\phi^{(I)} + \phi^{(HO)} = 0$, and the edge connecting the HO and I points in Fig. 1 represents a continuation from weak coupling ($\phi = - \infty$) to strong coupling ($\phi = + \infty$). It has been known for a long time that the two vacua have the same gauge symmetry ($SO(32)$) and the same supersymmetry, but it was unclear how they could be equivalent because type I strings and heterotic strings are very different. The explanation is that heterotic strings appear as nonperturbative excitations in the type I description. The converse is not quite true, because type I strings disintegrate at strong coupling.

The final link, labeled $\Omega$ in Fig. 1, connects the type IIB and type I vacua. $\Omega$ represents an "orientifold projection," which involves modding out by a particular $Z_2$ discrete symmetry. Starting from the IIB picture, the $Z_2$ in question is an orientation reversal of the IIB string ($\sigma \rightarrow - \sigma$) . This results in unoriented closed strings ("untwisted sector") and unoriented open strings carrying $SO(32)$ gauge symmetry ("twisted sector"). In the modern viewpoint, the open strings can be regarded as ending on 32 superimposed D9-branes. We will say more about D-branes later.

# $p$-branes

Supersymmetry algebras with central charges admit "short representations", the existence of which is crucial for testing conjectured non-perturbative properties of theories that previously were only defined perturbatively. Schematically, when a state carries a central charge $Q$, the supersymmetry algebra implies that its mass is bounded below $(M\geq |Q|)$. Moreover, when the state is "BPS saturated," *i.e.*, $M = |Q|$, the representation theory changes, and a state can belong to a short representation of the algebra. This phenomenon is already familiar for the case of Poincaré symmetry in 4d, which allows a massless photon to have just two helicity states (a short representation), whereas a massive vector boson must have three helicity states.

This BPS saturation property arises not only for point particles, characterized by a mass $M$, but for extended objects with $p$ spatial dimensions, called $p$-branes. In this case the central charge is a rank $p$ tensor. At first sight, this might seem to be in conflict with the Coleman–Mandula theorem, which forbids finite tensorial central charges. However, the $p$-branes carry a finite charge per unit volume, so that the total charge is infinite for a BPS $p$-brane that is an infinite hyperplane, and there is no contradiction. The BPS saturation condition in this case implies that the tension (or mass per unit volume) of the $p$-brane equals the charge density. Another way of viewing BPS $p$-branes is as solitons that preserve some of the supersymmetry of the underlying theory.

The theories in question (I will focus on the ones with 32 supercharges) are approximated at low energy by supergravity theories that contain various antisymmetric tensor gauge fields. They are conveniently represented by differential forms $$A_n \equiv A_{\mu_{1}\mu_{2} \ldots \mu_{n}} dx^{\mu_{1}} \wedge dx^{\mu_{2}}
\wedge \ldots \wedge dx^{\mu_{n}}.$$ In this notation, the corresponding gauge-invariant field strength is given by an $(n + 1)$-form $F_{n+1} = dA_n$ plus possible additional terms. A type II or 11d supergravity theory with such a gauge field has two kinds of BPS $p$-brane solutions, which preserve one-half of the supersymmetry. One, which can be called "electric," has $p = n - 1$. The other, called "magnetic," has $p = D - n - 3$, where $D$ is the space-time dimension (ten or eleven for the cases considered here).

A hyperplane with $p$ spatial dimensions in a space-time with $D-1$ spatial dimensions can be surrounded by a sphere $S^{D-p-2}$. If $A$ is a $(p + 1)$-form potential for which a $p$-brane is the source, the electric charge $Q_E$ of the $p$-brane is given by a straightforward generalization of Gauss's law: $$Q_E \sim \int_{S^{D-p-2}} * F,$$ where $S^{D-p-2}$ is a sphere surrounding the $p$-brane and $*F$ is the Hodge dual of the $(p + 2)$-form field strength $F$. Similarly, a dual $(D - p - 4)$-brane has magnetic charge given by $$Q_M \sim \int_{S^{p+2}} F.$$ The Dirac quantization condition, for electric and magnetic 0-branes in $D =
4$, has a straightforward generalization to a $p$-brane and a dual $(D-p-4)$-brane in $D$ dimensions $${1\over 2\pi} Q_E Q_M \in {\bf Z}.$$

An approximate description of the classical dynamics of a "thin" $p$-brane is given by a generalized Nambu–Goto formula $$S_p = T_p \int \left(\sqrt{-\det G_{\mu\nu}} + \ldots\right) d^{p+1} \sigma,$$ where $$G_{\mu\nu} = g_{MN} (X) \partial_\mu X^M \partial_\nu X^N.$$ Here $G_{\mu\nu}(\sigma)$ is a metric on the $(p+1)$-dimension world-volume of the $p$-brane obtained as a pullback of the D-dimensional space-time metric $g_{MN}(X)$. The functions $X^M (\sigma)$ describe the embedding of the $p$-brane in space-time. The coefficient $T_p$ is the $p$-brane tension – its universal mass per unit volume. Note that (for $\hbar = c = 1$) $T_p \sim
({\rm mass})^{p+1}$. This integral is just the volume of the embedded $p$-brane, generalizing the invariant length of the world-line of a point particle or the area of the world-sheet of a string. The dots represent terms involving other world-volume degrees of freedom required by supersymmetry.

Superstring theories in 10d have three distinct classes of $p$-branes. These are distinguished by how the tension $T_p$ depends on the string coupling constant $g_s$. A "fundamental" $p$-brane has $T_p \sim (m_s)^{p+1}$, with no dependence on $g_s$. Such $p$-branes only occur for $p = 1$ – the fundamental strings. Since these are the only objects that survive at $g_s = 0$, they are the only ones that can be used as the fundamental degrees of freedom in a perturbative description. A second class of $p$-branes, called "solitonic," have $T_p \sim (m_s)^{p+1}/g_s^2$. These only occur for $p = 5$, the five-branes that are the magnetic duals of the fundamental strings. This dependence on the coupling constant is familiar from field theory. A good example is the mass of an 't Hooft–Polyakov monopole in gauge theory. The third class of $p$-branes, called "Dirichlet" (or D$p$-branes), have $T_p \sim
(m_s)^{p+1}/g_s$. This behavior, intermediate between "fundamental" and "solitonic," was not previously known in field theory. In 10d type II theories D-branes occur for all $p\leq 9$ – even values in the IIA case and odd ones in the IIB case. They are all interrelated by T dualities; moreover, the magnetic dual of a D$p$-brane is a D$p'$-brane with $p' = 6 - p$. D-branes are very important, and so we will have more to say about them later.

Eleven-dimensional supergravity contains a three-form potential. Therefore, the 11d vacuum admits two basic kinds of $p$-branes – the M2-brane (also known as the supermembrane) and the M5-brane. These are EM duals of one another. Since the only parameter of the 11d vacuum is the Planck mass $m_p$, their tensions are necessarily $T_{M2} = (m_p)^3$ and $T_{M5} = (m_p)^6$, up to numerical coefficients.

We can use the relation between the 11d theory compactified on a circle of radius $R$ and the IIA theory in 10d to deduce the tensions of certain IIA $p$-branes. Starting with the M2-brane we can either allow one of its dimensions to wrap the circular dimension, leaving a string in the remaining dimensions, or we can simply embed it in the non-compact dimensions, where it is then still viewed as a 2-brane. In the latter case, the tension remains $m_p^3$. Using eqs.  () and (), we can recast this as $T = (m_s)^3/g_s$, which we recognize as the tension of the D2-brane of IIA theory. On the other hand, the wrapped M2-brane leaves a string with tension $T = m_p^3 R = m_s^2$. Thus we see that eq. () reflects the fact that a fundamental IIA string is actually a wrapped M2 brane. Starting with the M5-brane, we can carry out analogous calculations. If it is not wrapped we obtain a IIA 5-brane with tension $T = m_p^6 =
m_s^6/g_s^2$, which is the correct relation for the solitonic 5-brane (usually called the NS5-brane). If it is wrapped on the circle, one is left with a IIA 4-brane with tension $T = m_p^6 R = m_s^5/g_s$. This has the correct tension to be identified as a D4-brane. In other words, the D4-brane is actually a wrapped M5-brane.

There are a couple basic facts about D-branes in type II superstring theories that should be pointed out. First of all, they can be understood in the weak coupling limit (which makes them heavy) as surfaces on which fundamental type II strings can end. This is where the Dirichlet boundary conditions come in. This has a number of implications. One is that the dynamics of D-branes at weak coupling can be deduced from that of fundamental strings using perturbative methods. Another is that since a type II string carries a conserved charge that couples to a two-form potential, the end of a string must carry a point charge, which gives rise to electric flux of a Maxwell field. This implies that the world-volume theory of a D-brane contains a $U(1)$ gauge field. In fact, for strong fields that vary slowly it is actually a non-linear theory of the Born–Infeld type. The $U(1)$ gauge field can be regarded as arising as the lowest excitation of an open string with both ends attached to the D-brane.

Consider now $k$ parallel D$p$-branes, which are $(p+1)$-dimensional hyperplanes in ${\bf R}^{10}$. In this case, open strings can end on two different branes. The lowest mode of a string connecting the $i$th and $j$th D-brane is a gauge field that carries $i$ and $j$ type electric charges at its two ends. Altogether one has a $U(k)$ gauge theory in $p + 1$ dimensions. Classically, this can be constructed as the dimensional reduction of $U(k)$ super Yang–Mills theory in 10d. The separations of the D-branes are given by the vevs of scalar fields, which break the gauge group to a subgroup. For $p \leq 3$, these gauge theories have a straightforward quantum interpretation, but for $p > 3$ the gauge theories are non-renormalizable. I will return to this issue in Sect. 5.

# Brane-Configuration Constructions of SUSY Gauge Theories

In the last section we saw that a collection of $k$ parallel D-branes gives a supersymmetric $U(k)$ gauge theory. The unbroken supersymmetry in this case is maximal (16 conserved supercharges). In this section we describe more complicated brane configurations, which break additional supersymmetries, and give susy gauge theories in 4d with a richer structure. This is an active subject, which can be approached in several different ways. Here we will settle for two examples in one particular approach. (For a different approach see .)

The first example is a configuration of NS5-branes and D4-branes in type IIA theory depicted in Fig. 2. This configuration gives rise to an $SU(N_C)$ gauge theory in 4d with $N = 2$ supersymmetry (8 conserved supercharges). To explain why, one must first describe the geometry. All of the branes are embedded in 10d so as to completely fill the dimensions that will be identified as the 4d space-time with coordinates $x^0, x^1, x^2, x^3$. In addition, the NS5-branes also fill the $x^4$ and $x^5$ dimensions, which are represented by the vertical direction in the figure, and they have fixed values of $x^6, x^7, x^8,
x^9$. The D4-branes, on the other hand, have a specified extension in the $x^6$ direction, depicted horizontally in the figure, and they have fixed values of $x^4, x^5, x^7,
x^8, x^9$. The idea is that the gauge theory lives on the $N_C$ D4-branes, which are suspended between the NS5-branes. The $x^6$ extension becomes negligible for energies $E \ll 1/L$, where $L$ is the separation between the NS5-branes. In this limit the 5d theory on the D4-branes is effectively four dimensional. In addition there are $N_F$ semi-infinite D4-branes, which result in $N_F$ hypermultiplet flavors in the fundamental representation of the gauge group. These states arise as the lowest modes of open strings connecting the two types of D4-branes. The presence of the NS5-branes is responsible for breaking the supersymmetry from $N = 4$ to $N = 2$.

This picture is valid at weak coupling, because the gauge coupling constant $g_{YM}$ is given by $g_{YM}^2 =
g_s/(Lm_s)$, and the IIA picture is valid for small $g_s$. Substituting eq. (), we see that $g_{YM}^2 = R/L$, where $R$ is the radius of a circular eleventh dimension. So far, the description of the geometry omits consideration of this eleventh dimension, but by taking it into account we can see what happens to the gauge theory when $g_{YM}^2$ is not small and quantum effects become important. The key step is to recall that a D4-brane is actually an M5-brane wrapped around the circular eleventh dimension. Thus, reinterpreted as a brane configuration embedded in 11d, the entire brane configuration corresponds to a single smooth M5-brane. The junctions are now smoothed out in a way that can be made quite explicit. The correct configuration is one that is a stable static solution of the M5-brane equation of motion, which degenerates to the IIA configuration we have described in the limit $R \rightarrow 0$. There is a simple method, based on complex analysis, for finding such solutions. If space is described as a complex manifold, with a specific choice of complex structure, then the brane configuration is a stable static solution if its spatial dimensions are embedded holomorphically. In the example at hand, the relevant dimensions are two dimensions of the M5-brane, which are embedded in the four dimensions denoted $x^4, x^5, x^6, x^{10}$, where $x^{10}$ is the circular eleventh dimension. A complex structure is specified by choosing as holomorphic coordinates $v = x^4
+ i x^5$ and $t = \exp [(x^6 + i x^{10})/R]$, which is single-valued. Then a holomorphically embedded submanifold is specified by a holomorphic equation of the form $F(t,v) = 0$. The appropriate choice of $F$ is a polynomial in $t$ and $v$ with coefficients that correspond in a simple way to the positions of the NS5-branes and D4-branes. (For further details see Ref. .) This 2d surface is precisely the Seiberg–Witten Riemann surface (or "curve") that characterizes the exact non-perturbative low-energy effective action of the gauge theory. When first discovered, this curve was introduced as an auxiliary mathematical construct with no evident geometric significance. We now see that the Seiberg–Witten solution is given by an M5-brane with four of its six dimensions giving the space time and the other two giving the Seiberg–Witten curve! This simple picture makes the exact non-perturbative low energy physics of a wide class of $N = 2$ gauge theories almost trivial to work out.

Let me briefly mention how the brane configuration described above can be modified to describe certain $N = 1$ susy gauge theories. One way to achieve this is to rotate one of the two NS5-branes so that it fills the dimensions $x^8, x^9$ and has fixed $x^5, x^6$ coordinates. When this is done the $N_C$ D4-branes running between the NS5-branes are forced to be coincident. The rotation breaks the supersymmetry to $N = 1$. One of the remarkable discoveries of Seiberg is that an $N =
1$ susy gauge theory with gauge group $SU(N_C)$ and $N_F\geq N_C$ flavors is equivalent in the infrared to an $SU(N_F - N_C)$ gauge theory with a certain matter content. This duality can be realized geometrically in the brane configuration picture by smoothly deforming the picture so as to move one NS 5-brane to the other side of the other one . Such a move certainly changes the exact quantum vacuum described by the configuration. However, the parameters involved are irrelevant in the infrared limit, so one achieves a simple understanding of Seiberg duality.

# New Non-gravitational 6d Quantum Theories

We have seen that it is interesting and worthwhile to consider the world volume theory of a collection of coincident or nearly coincident branes. For such a theory to be regarded in isolation in a consistent way, it is necessary to define a limit in which the brane degrees of freedom decouple from those of the surrounding space-time "bulk." Such a limit was implicitly involved in the discussion of the preceding section. (This involves some subtleties, which we did not address.) In this section we wish to consider the 6d world-volume theory that lives on a set of (nearly) coincident 5-branes. If one can define a limit in which the degrees of freedom of the world-volume theory decouple from those of the bulk, but still remain self-interacting, then we will have defined a consistent non-trivial 6d quantum theory . (The only assumption that underlies this is that M theory/superstring theory is a well-defined quantum theory.) The 6d quantum theories that are obtained this way do not contain gravity. The existence of consistent quantum theories without gravity in dimensions greater than four came as quite a surprise to many people.

As a first example consider $k$ parallel M5-branes embedded in flat 11d space-time. This neglects their effect on the geometry, which is consistent in the limit that will be considered. The only parameters are the 11d Planck mass $m_p$ and the brane separations $L_{ij}$. In 11d an M2-brane is allowed to terminate on an M5-brane. Therefore, a pair of M5 branes can have an M2-brane connect them. When the separation $L_{ij}$ becomes small, this M2-brane is well approximated by a string of tension $T_{ij} = L_{ij} m_p^3$. The limit that gives decoupling of the bulk degrees of freedom is $m_p \rightarrow \infty$. By letting the separations approach zero at the same time, this limit can be carried out holding the string tensions $T_{ij}$ fixed. In the limit one obtains a chiral 6d quantum theory with $(2,0)$ supersymmetry containing $k$ massless tensor supermultiplets and a spectrum of strings with tensions $T_{ij}$. There are five massless scalars associated to each brane (parametrizing their transverse excitations). They are coordinates for the moduli space of the resulting theory, which is $({\bf R}^5)^k/S_k$. The permutation group $S_k$ is due to quantum statistics for identical branes. String tensions depend on position in moduli space, and specific ones approach zero at its singularities.

A closely related construction is to consider $k$ parallel NS 5-branes in the IIA theory. The difference in this case is that one of the transverse directions (parametrized by one of the five scalars) is the circular eleventh dimension. In carrying out the decoupling limit one can send the radius $R$ to zero at the same time, holding the fundamental type IIA string tension $T =
m_s^2 = m_p^3 R$ fixed. The resulting decoupled 6d theory contains this string in addition to the ones described above. It becomes bound to the NS5-branes in the limit, as the amplitude to come free vanishes in the limit $g_s \rightarrow 0$. The resulting theory has the moduli space $({\bf R}^4 \times S^1)^k/S_k$. This theory contains fundamental strings and has a chiral extended supersymmetry, features that are analogous to type IIB superstring theory in 10d. However, it is actually a class of non-gravitational theories (labeled by $k$) in 6d. Because of the analogy some authors refer to this class of theories as $iib$ string theories. Six-dimensional non-gravitational analogs of type IIA string theory, denoted $iia$ string theories, are obtained by means of a similar decoupling limit applied to a set of parallel NS5-branes in IIB theory. These $iia$ and $iib$ string theories are related by T duality. Explicitly, compactifying one spatial dimension on a circle of radius $R_a$ or $R_b$, the theories (with given $k$) become equivalent for the identification $m_s^2 R_a R_b = 1$. This feature is directly inherited from the corresponding property of the IIA and IIB theories.

There are various generalizations of these theories that will not be described here. There are also 6d non-gravitational counterparts of the two 10d heterotic theories. These have chiral ($1,0$) supersymmetry. In the notation of Fig. 1, they could be referred to as $he$ and $ho$ theories. They, too, are related by T duality. Although the constructions make us confident about the existence and certain general properties of these theories, they are not very well understood. The 10d string theories have been studied for many years, whereas these 6d string theories are only beginning to be analyzed. Like their 10d counterparts, the fact that they have T dualities implies that they are not conventional quantum field theories.

# The Matrix Theory Proposal

The discovery of string dualities and the connection to 11d has taught us a great deal about non-perturbative properties of superstring theories, but it does not constitute a complete non-perturbative formulation of the theory. In October 1996, Banks, Fischler, Shenker, and Susskind made a specific conjecture for a complete nonperturbative definition of the theory in eleven uncompactified dimensions called 'Matrix Theory' . In this approach, as we will see, other compactification geometries require additional inputs. It is far from obvious that the BFSS proposal is well-defined and consistent with everything we already know. However, it seems to me that there is enough that is right about it to warrant the intense scrutiny that it has received and is continuing to receive. At the time of this writing, the subject is in a state of turmoil. On the one hand, there is a new claim that the BFSS prescription (as well as a variant due to Susskind ) can be derived from previous knowledge . On the other, some people are (cautiously) claiming to have found specific settings in which it gives wrong answers! In the following, we do not comment further upon these claims. Instead, we describe the basic ideas of Matrix Theory, as well as some of its successes and limitations.

One of the $p$-branes that has not been discussed yet is the D0-brane of type IIA theory in 10d. Being a D-brane, its mass is $M = m_s/g_s$. Using eq. (), one sees that $M = 1/R$, which means that it can be understood as the first Kaluza–Klein excitation of the 11d supergravity multiplet on the circular eleventh dimension. In fact, this is a good way of understanding (and remembering) eq. (). Like all the type II D-branes it is a BPS state that preserves half of the supersymmetry, so one has good mathematical control. From the 11d viewpoint it can be viewed as a wave going around the eleventh dimension with a single quantum of momentum. Higher Kaluza–Klein excitations with $M = N/R$ are also BPS states. From the IIA viewpoint these are bound states of $N$ D0-branes with zero binding energy. The existence of a bound state at a threshold is a very subtle dynamical question, which must be true in this case. This has in fact been proved for $N=2$ in ref. and for all prime values of $N$ in .

By the prescription given in Sect. 4, the dynamics of $N$ D0-branes is described by the dimensional reduction of $U(N)$ super Yang–Mills theory in 10d to one time dimension only. When this is done, the spatial coordinates of the $N$ D0-branes are represented by $N\times N$ Hermitian matrices! This theory has higher order corrections, in general. However, one can speculate that these effects are suppressed by viewing the $N$ D0-branes in the infinite momentum frame (IMF). This entails letting $p_{11} = {N/ R}$ approach infinity at the same time as $R \rightarrow \infty$. The techniques involved here are reminiscent of those developed in connection with the parton model of hadrons in the late 1960's. The BFSS conjecture is that this IMF frame $N \rightarrow
\infty$ limit of the D0-brane system constitutes an exact non-perturbative description of the 11d quantum theory. The $N \rightarrow \infty$ limit is awkward, to say the least, for testing this conjecture. A stronger version of the conjecture, due to Susskind, is applicable to finite $N$. It asserts that the IMF D0-brane system, with fixed $N$, provides an exact non-perturbative description of the 11d theory compactified on a light-like circle with $N$ units of (null) momentum along the circle.

One of the first issues to be addressed was how this conjecture should be generalized when additional dimensions are compact, specifically if they form an $n$-torus $T^n$. The reason this is a non-trivial problem is that open strings connecting pairs of D0-branes can lie along many topologically distinct geodesics. It turns out that all these modes can be taken into account very elegantly by replacing the 1-dimensional quantum theory of the D0-branes by an $(n + 1)$-dimensional quantum theory, where the $n$ spatial dimensions lie on the dual torus $\tilde T^n$. The extra dimensions precisely account for all the possible stretched open strings. This picture had some immediate successes. For example, it nicely accounted for all the duality symmetries for various values of $n$. However, $(n + 1)$-dimensional super Yang–Mills theory is non-renormalizable for $n>3$, so this description of the theory is certainly incomplete in those cases. The new theories described in Sect. 5 provide natural candidates when $n = 4$ or $5$, but when $n>5$ there are no theories of this type, and so we seem to be stuck.

In conclusion, Matrix Theory is a very interesting proposal for defining M theory non-perturbatively. Whether it is correct, or needs to be modified, is very much up in the air at the present time. However, even if it is right, it is unclear how to define vacua with more than five compact dimensions. This fact is very intriguing, since this is precisely what is required to describe the world that we observe.

[^1]: Lecture presented at the 31st International Ahrenshoop Symposium in Buckow, Germany.

[^2]: Work supported in part by the U.S. Dept. of Energy under Grant No. DE-FG03-92-ER40701.