\section{Introduction}
\label{introduction}

In the last decades light-matter interaction has been studied in a large
variety of experiments not only in many body systems such as atomic
ensembles, where coherent laser pulses are transferred and stored to a
collective excitation~\cite{lukin2003,hammerer2010}, but also on the
level of single matter particles.
Efficient coupling of single atoms to single photons has been
demonstrated by locating an atom inside an optical
cavity~\cite{raimond2001, walther2006, kimble1998, rempe1993}.
Such efficient coupling is mandatory in many proposals for
quantum information protocols which rely on quantum memories and
processors, i.e., on the efficient interaction
of a quantum state of light with the internal quantum states of matter.
Light is a fast and robust way to transmit
quantum information, but the transfer over long distances requires
quantum repeaters based on such quantum
memories~\cite{cirac1997,briegel1998,duer1999}.
The rate of information transfer
through the quantum channel depends on the process efficiency of both
storage and read-out.
Such efficient coupling can be achieved using a cavity based scheme and
accepting the corresponding modification of the mode spectrum from
continuous to discrete.

An alternative -- and quite natural -- approach is the coupling in
\emph{free space}.
In this approach, the \emph{direct} coupling of light to matter has
to be maximized instead of the coupling to a cavity.
Besides their potential use in various applications, efficient free
space coupling schemes offer large design freedom.
They enable the investigation of light--matter interaction
without modifying the density of states of the electromagnetic field
continuum.
Alternatively one can study the effect of, e.g., a modified quantum
vacuum~\cite{gardiner1986}.

Recently, one could witness increased efforts devoted to the
topic of efficient free space coupling.
When discussing the corresponding advances one can clearly distinguish
between two different processes:
One is the elastic scattering of a laser beam by an interacting
particle for which no substantial excitation of the matter system
occurs, i.e., the probability to find the matter system in its excited
state is low throughout the interaction process.
The efficient elastic scattering manifests itself in the attenuation
of a weak coherent laser beam resonant with a transition of the
interacting particle, as has been shown for single molecules
\cite{wrigge2008}, single quantum dots \cite{vamivakas2007}, single
atoms \cite{tey2008-np} or single ions \cite{slodicka2010}.
Another hallmark of efficient elastic scattering is the phase shift of
a weak laser beam detuned from resonance, which was demonstrated
recently for the free space scenario
\cite{aljunid2009,pototschnig2011}.
Corresponding theoretical contributions can be found in
Refs. \cite{kochan1994,vanenk2001,vanenk2004,zumofen2008,tey2009}.

The other process treated in recent literature is the process of
efficient absorption of a photon by a single atom, i.e., the matter
system is to be brought to its excited state with high probability.
Note that this should not be confused with a $\pi$-pulse.
While a $\pi$-pulse -- typically containing many photons -- transfers
the atom efficiently to the excited state the pulse cannot be fully
absorbed by a single atom.
Furthermore, the rotation angle of the Bloch vector (e.g. $\pi$) is given
by time integrating the electric field envelope, while a single photon
Fock state has a zero electric field expectation value.
Theoretical and conceptual contributions devoted to the absorption of
single photons can be found in
Refs. \cite{quabis2000,lindlein2007,sondermann2007,pinotsi2008,stobinska2009,wang2011}.
A recent experimental demonstration of the absorption of a single photon by
a single atom was given by Piro et al~\cite{piro2011}, whereas an
early demonstration of reduced transmission of a bright laser beam
by a single ion was performed by Wineland et al~\cite{wineland1987}.

It is common to the two processes of elastic scattering and absorption
of a photon that their efficiencies are proportional to the extent to
which the incident light mode resembles an electric dipole wave
\cite{quabis2000,vanenk2004,lindlein2007,pinotsi2008,zumofen2008}.
More precisely, the spatial properties of the incident light mode should
match the ones of the wave that would emerge from the dipole moment of the
interacting particle.
Here, we show how to experimentally generate such a dipole wave
envisioning close to optimum coupling efficiencies.
As outlined below, we do this by preparing a plane light wave that, after
reflection off a parabolic mirror, is transformed into a
dipole wave \cite{lindlein2007,sondermann2007}.
An essential property of this setup is that the parabolic mirror
covers nearly the complete solid angle surrounding the dipolar absorber.

For the specific aim of perfectly exciting a single two-level atom
with a single photon one may take guidance from the spontaneous emission
process:
The incident photon must resemble the time reversed version of a
spontaneously emitted photon
\cite{quabis2000,lindlein2007,sondermann2007}.
This includes the dipole properties discussed above.
Furthermore, time reversal arguments demand that the temporal envelope
of the incident photon is exponentially increasing with a time
constant equal to the lifetime of the atom's excited state
\cite{sondermann2007}.
This heuristic approach has been confirmed by detailed calculations
\cite{stobinska2009}.

There exist several
methods for temporal shaping of single photons in general.
One approach is to shape the photon during its creation process, for
instance by temporal modulation of the intensity of the control laser
in schemes based on electromagnetically induced transparency~\cite{eisaman2005}
or by modulation of the pump beams driving the process of four wave
mixing \cite{chen2010} or parametric
down-conversion (PDC)~\cite{kalachev2010}.
Also the shaping of single photons by spectrally filtering one of
two PDC photons has been demonstrated~\cite{haase2009}.
In the case of single photon generation with a single emitter located
inside an optical resonator~\cite{kuhn2002, keller2004, mckeever2004,
  bochmann2008} the temporal modulation of the pump laser intensity
changes the emitted single photon waveform.
Furthermore a heralded single photon can also be modulated directly
using electro-optic modulation~\cite{kolchin2008,olislager2010}.

However, with the mentioned techniques it is challenging to generate a
single photon Fock state with an exponentially increasing temporal envelop,
especially at wavelengths in the ultraviolet spectral range (which are of
interest in this paper, see below).
Instead we use coherent laser pulses which are highly
attenuated to less than one photon per pulse.
This procedure is often used as a substitute for generating single
photon Fock states~\cite{gisin2002}.
For the sake of simplicity, we will use the term
'single photon' also for strongly attenuated coherent states,
keeping in mind that all excitation probabilities mentioned
in this paper are strictly valid only for true single photon wave
packets.
Nevertheless, if the amplitude of the coherent state is low enough,
contributions from Fock states with photon numbers larger than one are
negligible.
It is then a good approximation to apply the excitation probabilities
for true single photon states to that instances of the coherent state
wave packets that would be projected onto a single photon state upon
measuring the photon number.

In what follows, we present and characterize experimental techniques
to produce the optimal optical mode for perfectly coupling light with
an atomic dipole, especially for the case of exciting a single atom
with a single photon.
The paper is organized as follows:
In the next section we motivate a measure for the coupling strength
between a single atom (or any other matter system) and light and relate
this measure to the absorption efficiency.
In Sec. \ref{Generation_optimum_doughnut} the experimental generation
of a radially polarized Laguerre-Gauss mode is shown.
We present the characterization of this beam and calculate the overlap
of the experimental data with the optimal dipole wave.
Another property of dipole waves is their uniform phase.
Therefore, it is essential to compensate for phase aberrations of the
device that focuses the incident light towards the matter particle.
This topic is treated in Sec. \ref{Aberration_correction}.
Section \ref{Temporal_pulse_shaping} is devoted to the shaping of the
temporal properties of a weak coherent state with the aim of
maximizing absorption.
This aspect as well as the spatial mode shaping is exemplified for
two specific linear dipole transitions of singly and doubly ionized
Ytterbium, but the corresponding schemes can be adapted to any dipole
transition.
The properties of our exemplary transitions are listed in
table~\ref{table_transitions}.
Finally, in the last section we estimate the coupling efficiencies and
absorption probabilities that would be achieved when interfacing the
generated light waves with the corresponding single ions.

\begin{table}
\begin{tabular}{c|c|c|c|c}
	Ion & Transition $\quad$  & Wavelength ($\lambda$)
        $\quad$ & Lifetime ($\tau$) & Label \\ \hline \hline
	$^{174}$Yb$^{+}$ & $^{2}$S$_{1/2}\rightarrow\ ^{2}$P$_{1/2}$
        & 369.5~nm & 8.1~ns & T1 \\
	$^{174}$Yb$^{2+}$  & $^{1}$S$_{0}\rightarrow\ ^{3}$P$_{1}$
        & 251.8~nm & 230~ns & T2 \\ \hline
\end{tabular}
\caption{\label{table_transitions}
Overview of the dipole transitions treated exemplarily in this paper.
Lifetime data are taken from Refs.~\cite{berends1993,zhang2001}.
The labels T1 and T2 will be used in the text for the sake of
simplicity.}
\end{table}

\section{Definition of the free space coupling strength}
\label{coupling_strength}

The coupling constant of the interaction of light and matter in the
dipole approximation is given by the Rabi frequency
$g\sim\vec{\mu}\cdot\vec{E}/\hbar$ in free space as well as in an optical
resonator, where $\vec{\mu}=\mu\cdot\vec{e}_\mu$ is the atomic dipole
moment and $\vec{E}$ the electric field.
Depending on whether the treatment is semi-classical or fully
quantized, $\vec{\mu}$ alone or both $\vec{\mu}$ and $\vec{E}$ are
operators, respectively.
In either case a large coupling constant requires a large portion of the
incident electric field to be parallel to the atomic dipole.

In cavity quantum electro-dynamics the coupling strength is
proportional to $g^2$.
Here, we adopt this measure (cf. Ref. \cite{sondermann2007}) and
define the free space coupling strength $G$ as the electric energy
density in the focus that corresponds to the field component parallel
to the atomic dipole and normalize it to the maximum possible one for
a given input power:
\begin{equation}
G= \frac{ |\vec{E}\cdot\vec{e}_\mu|^2}{|E_\textrm{max}|^2}
 \quad .
\end{equation}

In the focus, the electric field component along a specific dipole
moment is proportional to the square root of the solid angle
covered by the mediating focusing optics $\Omega_\mu$ and to the
overlap $\eta$ of the incident field with the dipole field
distribution $\vec{E}_\mu$ \cite{sondermann2008}:
$\vec{E}\cdot\vec{e}_\mu\sim \sqrt{\Omega_\mu}\cdot\eta$.
It should be noted that $\Omega_\mu$ is the solid angle weighted
with the angular dipole radiation pattern $D(\vartheta)$:
$\Omega_\mu=\int D(\vartheta)\sin\vartheta d\vartheta d\varphi$.
The overlap
\begin{equation}
\label{eq:eta_spherical}
\eta=\frac{
\int \vec{E}^\star\cdot\vec{E}_\mu \sin\vartheta\ d\vartheta d\varphi
}
{
\sqrt{\int |\vec{E}|^2 \sin\vartheta\ d\vartheta d\varphi
\cdot \int |\vec{E}_\mu|^2 \sin\vartheta\ d\vartheta d\varphi
}
}
\end{equation}
is calculated integrating over the solid angle covered by
the focusing optics \cite{sondermann2008}.
The maximum values of these quantities are
$\Omega_\mu^\textrm{max}=8\pi/3$ and $\eta^\textrm{max}=1$.
Therefore, the free space coupling strength is written as
\begin{equation}\label{eq:G}
G= \frac{3}{8\pi}\cdot \Omega_\mu \cdot \eta^2 \quad .
\end{equation}
It reaches unity only if the full solid angle is covered by the
focusing optics and the mode matching is perfect.
Likewise the maximum coupling strength for half solid angle
illumination is $G=0.5$.
These values correspond to the 'scattering ratios' of four and two,
respectively, as used by Tey et al \cite{tey2009} and Zumofen et al
\cite{zumofen2008}.

$G$ can be interpreted as a geometric quantity that is determined by
the focusing setup and the spatial properties of the incident light mode.
For a classical field or a coherent state wave packet $\sqrt{G}$
scales the field amplitude at the focus (or its expectation value,
respectively).
If the incident light field is a Fock state, $\sqrt{G}$ scales the
magnitude of the uncertainty of the incident field at the focus.

While $G=1$ is necessary for reaching the maximum effect
in all sorts of free space interaction (phase shifts, scattered
power, absorption, etc.) it is not sufficient for reaching unit probability
during the absorption of a single photon by a single atom.
As mentioned above, the temporal profile of the single photon wave
packet plays a decisive role.
The optimum effect, i.e., 100\% excitation of a two-level atom can
only be obtained if the temporal profile matches an increasing
exponential with proper time
constant \cite{sondermann2007,stobinska2009,wang2011}.
To quantify deviations from this profile we introduce the temporal
overlap $\eta_t$ of the incident field amplitude $E_\textrm{inc}(t)$
with the amplitude
$E_\textrm{ideal}(t)=e^{\frac{\Gamma}{2}t}\cdot\theta(-t)$ of the
ideal pulse shape:
\begin{equation}
\label{eq:eta_t}
\eta_t=\frac{\int_{-\infty}^\infty E_\textrm{inc}(t) \cdot
   E_\textrm{ideal}(t) dt}
{\sqrt{\int_{-\infty}^\infty |E_\textrm{inc}(t)|^2 dt / \Gamma}}
\quad .
\end{equation}
$\theta(t)$ is the Heaviside function and $\Gamma$ the
spontaneous emission rate of the two-level system to be excited.
This overlap can be considered as a generalization of Eq. (8) of
Ref. \cite{stobinska2009} which treats the case of
an exponentially increasing pulse amplitude.
It is assumed that the incident signal is shifted in time by an amount
that maximizes $\eta_t$.
Since the respective equation in Ref. \cite{stobinska2009} determines
the probability amplitude to find an atom in the excited state,
we take $\eta_t^2$ as a multiplicative factor determining the absorption
probability
\begin{equation}
\label{eq:Pabs}
P_a= G \cdot \eta_t^2 = \frac{3}{8\pi}\cdot \Omega_\mu \cdot \eta^2
\cdot \eta_t^2 \quad .
\end{equation}

The prospective coupling efficiencies $G$ and absorptions
probabilities $P_a$ based on the experimentally achieved overlaps will
be given in Sec. \ref{discussion}.

\section{Generation of the optimum spatial mode}
\label{Generation_optimum_doughnut}

A dipole wave incident from nearly the full solid angle is
generated by reflection of a collimated laser beam from a deep
parabolic mirror.
The ideal beam shape incident on the parabolic mirror is the time
reversed version of a dipole wave that emerges from the focus of the
mirror and is reflected at the parabolic surface.
After reflection off the parabolic mirror the field vector
$\vec{E}_\mu$ of light emitted by a linear dipole
located at the focal point and with the quantization axis along the
mirror axis is described by~\cite{lindlein2007, sondermann2008}
\begin{equation}
\label{eq:dipole}
\vec{E}_\mu\left( \dfrac{r}{f} \right) =
E_{0}\cdot
\frac{
  \dfrac{r}{f}}{ \left[ \left( \dfrac{r}{2f} \right) ^2 +1 \right] ^2}
\cdot\vec{e}_{\text{r}}
\end{equation}
where $f$ is the focal length of the mirror, $\vec{e}_{\text{r}}$ the
unit vector in radial direction, $E_0$ is an amplitude factor and a
factor $\exp(i\vec{k}\vec{r})$ has been omitted.
This field distribution has to be generated by an optical element from
an incoming plane wave.
The amplitude distribution of this ideal mode can in
principle be created using diffractive optical elements. However,
there are several drawbacks to this method~\cite{lindlein2007}.
Instead we generate a radially polarized doughnut beam which can have
a large overlap with the ideal mode and can be created without using
diffractive optics.

The field vector $\vec{E}_{\text{rpd}}$
of a radially polarized doughnut mode can be expressed as
\begin{equation} \label{eq_rad_pol_doughnut}
\vec{E}_{\text{rpd}}\left( \dfrac{r}{f}\right)=E_0\cdot
\dfrac{r}{f}
\cdot e^{-\frac{(r/f)^2}{(w/f)^2}} \cdot \vec{e}_{\text{r}} \quad.
\end{equation}
The overlap of the radially polarized doughnut mode with the ideal
distribution is optimized by adjusting the beam radius $w$ with the
maximum achievable overlap depending on the magnitude of the weighted
solid angle $\Omega_\mu$ (i.e., the depth of the mirror with respect
to its focal length or the half opening angle) \cite{sondermann2008}.
Our parabolic mirrors have a focal length of $f=$2.1~mm and a front
aperture radius of 10~mm.
This corresponds to a length of $5.67\cdot f$ (aperture radius
$r_\textrm{max}/f=4.76$, half opening angle of 134$^{\circ}$) and a
weighted solid angle for a linear dipole of
$\Omega_\mu$=0.94$\cdot\frac{8\pi}{3}$.
The maximum possible overlap for this geometric parameters is
$\eta$=0.982 for a beam radius $w=2.26\cdot f$ \cite{sondermann2008}.
This limits the experimentally achievable coupling strength $G$ and
absorption probability $P_a$ to 0.91.
A bore of 0.75~mm radius concentric to the optical axis of the
mirror (cf. Sec. \ref{Aberration_correction}) induces only negligible
changes to $\Omega_\mu$ and $\eta$.

\begin{figure}
\begin{center}
\includegraphics{Golla_fig1.pdf}
\caption{\label{fig1}
  Experimental setup.
  The temporal profile is formed  using an arbitrary waveform
  generator (AWG) which determines the temporal envelope of the light
  power scattered into the first diffraction order of an acousto
  optical modulator (AOM).
  The output of the AOM is filtered by a single mode fiber (SMF).
  Linearly polarized light is sent through a segmented half wave plate
  (HWP).
  Subsequently the produced doughnut mode is filtered by a pinhole
  aperture and expanded to appropriate dimensions.
  The aberrations of the parabolic mirror (PM) are compensated using
  a phase plate.
}
\end{center}
\end{figure}

There are different methods for producing a radially polarized doughnut
beam. Examples are the use of sub-wavelength dielectric
gratings~\cite{bomzon2002,ghadyani2011}, nematic liquid
crystals~\cite{stalder1996} and the coherent addition of the two modes
TEM$_{01}$ and TEM$_{10}$ inside laser cavities~\cite{oron2000} or
interferometers~\cite{tidwell1990, tidwell1993}.
We use a segmented half-wave plate similar to the one
described in Refs.~\cite{dorn2003,quabis2005} to generate the radially
polarized doughnut beam.
It consists of $n=8$ segments  of a low order
half-wave plate suitable for both wavelengths (370~nm and 252~nm).
The experimental setup to generate the doughnut mode is shown in
Fig.~\ref{fig1}.
A continuous wave laser beam at the wavelength of interest first
passes the setup for temporal pulse shaping
(cf. Sec. \ref{Temporal_pulse_shaping}).
It is then sent through a polarization maintaining single mode fiber
(370~nm beam) or a filtering telescope with a pinhole aperture (252~nm
beam), respectively.
The resulting linearly polarized Gaussian beam is centered onto the
segmented half-wave plate where its polarization is rotated in every
segment by an angle of $2\pi/n$, obtaining a full $2\pi$ rotation of
the electrical field vector over the whole beam cross section.
After propagation through the segmented half-wave plate the beam contains
the radially polarized doughnut mode and modes of higher order.
The latter ones are filtered out using a pinhole aperture in the
focal plane of a Keplerian telescope.
With this setup the doughnut mode is generated with a power
efficiency of 70\%.

\begin{figure*}
\begin{center}
\includegraphics{Golla_fig2.pdf}
\caption{ \label{fig2}
  Characterization of the radially polarized doughnut beam at
  370~nm (transition T1):
  (a) Stokes parameter $S_0$ (total intensity).
  (b) Orientation angle $\psi$.
  (c) Ellipticity angle $\chi$.
  The shown data points cover the aperture of the parabolic mirror.
}
\end{center}
\end{figure*}

\begin{figure*}
\begin{center}
\includegraphics{Golla_fig3.pdf}
\caption{ \label{fig3}
  Characterization of the radially polarized doughnut beam at
  252~nm (transition T2):
  (a) Stokes parameter $S_0$ (total intensity).
  (b) Orientation angle $\psi$.
  (c) Ellipticity angle $\chi$.
  The shown data points cover the aperture of the parabolic mirror.
}
\end{center}
\end{figure*}

The polarization properties of the generated doughnut beam are
determined performing spatially resolved  measurements of the Stokes
parameters based on the method described in Ref. \cite{schaefer2007}.
The profile of a beam sent through a rotatable quarter-wave plate and a
non rotatable linear polarizer is measured with a CMOS camera.
Due to the technical limitations of the CMOS camera the Stokes
measurements are not performed with single photon pulses but
with continuous wave laser light.

Once the Stokes parameters are measured the total intensity $S_0$ and
the characteristic angles of the polarization ellipse are
calculated for every camera pixel.
Since the CMOS camera chip used in the Stokes measurements is much
smaller than the front aperture of the parabolic mirror, the
measurements have been performed before the final magnification step.
Nevertheless, the data are presented in a coordinate system scaled
such that the beam radius of the measured
doughnut mode matches the optimum beam radius predicted by theory.
The results are shown in Figs.~\ref{fig2}
and~\ref{fig3} for 370~nm and 252~nm, respectively.

The orientation angle $\psi \in \left[0 ,\pi\right]$ describes the
rotation of the main axis of the polarization ellipse with respect to
the x-axis.
As shown in Figs.~\ref{fig2}b) and
\ref{fig3}b), at both wavelengths $\psi$ covers the full range
of $0$ to $\pi$ in the lower and the upper half plane of the beam
cross section.
The value of $\psi$ changes continuously with rotation around the beam
center.
The relative phase of the electric field is the same in both half
planes, otherwise a defect line would be observed in $S_0$ at the
$y=0$ line where $\psi$ jumps from zero to $\pi$.
Hence, the orientation of the polarization vector is radial on the
majority of the beam cross section.
Deviations from this behavior occur mainly for radii $\gtrsim 3.5\cdot
f$, i.e., in regions of low intensity.
We attribute these deviations to measurement noise originating from the
discretization of the camera pixels' intensity values.
Furthermore, measurement noise is particularly cumbersome in the
vicinity of the $x$-axis, since there it leads to fluctuations between
values of $\psi=0$ and $\psi=\pi$.
When calculating the overlap with an ideal, noise-free state of
polarization (performed below) these fluctuations lead to pixels on
which the scalar product of the field vectors becomes negative.
For pixels with a large value of $S_0$ this results in a
measurement-induced reduction of the overlap integral.

The second angle describing the polarization ellipse is the
ellipticity angle $\chi\in[-\pi/4,\pi/4]$ giving the amount of
elliptical polarization.
Figs. \ref{fig2}c) and \ref{fig3}c) show
the ellipticity angles over the cross section of the generated
doughnut beams.
On the majority of the beam cross section the ellipticity angle is
rather low, i.e., the local state of polarization can be regarded as
linear with a radial orientation.

Next we calculate the overlap of the ideal electric dipole field and
the generated radially polarized doughnut mode in the entrance plane
of the focusing parabolic mirror, assuming perfect reflection at the
mirror surface.
This gives the same value as
calculating the mode overlap on a sphere close to the focal point, as
described in Ref.~\cite{sondermann2008}.
The integration is performed over the mirror aperture, i.e., only
pixels lying within the boundary $r\le r_\textrm{max}=4.76\cdot f$ are
contributing.
The spatial field overlap $\eta$ of the measured radially
polarized doughnut beam with the ideal dipole distribution reflected
at the mirror is calculated by taking the normalized scalar product of
the electric fields
\begin{equation} \label{eq_overlap_doughnut_dipole}
\eta=\dfrac{\int E_\text{meas}^\star \cdot
  \vec{e}_\text{meas}^\star \cdot \vec{E}_\mu r d r
  d\phi }{\sqrt{\int \vert E_{\text{meas}}\vert^{2} r d r d\phi
    \int\vert \vec{E}_\mu\vert^{2} r d r d\phi}} \quad ,
\end{equation}
where $\vec{e}_{\text{meas}}$  is the normalized measured polarization
vector constructed from  the angles $\psi$ and $\chi$.
$E_{\text{meas}}$  is calculated  by taking the square root of the
measured intensity distribution (the Stokes parameter $S_0$) and
$\vec{E}_\mu$ is given by Eq. (\ref{eq:dipole}).
Equation (\ref{eq_overlap_doughnut_dipole}) is the equivalent to
Eq. (\ref{eq:eta_spherical}).

For the mode driving transition T1 we calculate the field overlap of the
generated radially polarized doughnut mode with the ideal dipole wave
to be 0.963.
If we assume perfect orientation angles for the measured beam, i.e.,
we disregard the discretization effects discussed above, we find
$\eta_\textrm{370}=0.979$.
The overlap of the generated doughnut mode for transition T2 with the
ideal dipole mode is calculated to be 0.939. Assuming ideal
orientation angles we arrive at $\eta_\textrm{252}=0.975$.

Our parabolic mirrors are made from aluminum.
One might wonder whether the radially varying amplitude of the
reflection coefficient, which is due to the  radially varying angle of
incidence, might influence the obtainable overlaps (see also
Sec. \ref{Aberration_correction}).
Accounting for the radial variation of the reflectivity in
calculating the overlap of the modes given by Eqs. (\ref{eq:dipole})
and (\ref{eq_rad_pol_doughnut}) leads to a slight increase by 0.03\%
in comparison to the maximum obtainable mode overlap for constant
reflectivity.

Therefore, we conclude that the obtainable overlaps are not
influenced by the material properties of aluminum in a significant way.

\section{Aberration correction}
\label{Aberration_correction}

An essential prerequisite for generating dipole modes is an
aberration free phase front:
The phase of the focused wave should be uniform on a
spherical surface enclosing the focal point.
Such an ideal spherical wave can only be generated if the parabolic
mirror is free of surface deviations. Since any realistic device has
aberrations (see, e.g.,
Refs.~\cite{leuchs2008,drechsler2001,stadler2008}), these need to be
characterized and corrected for.
Ideally the wavefront incident onto the parabolic mirror can be shaped
such as to create an aberration free inward moving spherical wave
after reflection off the mirror.

One method of measuring the imperfections of a parabolic mirror is to
scan its surface using a
profilometer~\cite{stadler2008,drechsler2001}. But the huge relative
depth of our parabola makes it hardly accessible for profilometer
heads.
We measure the shape deviations by
interferometry with nanometer accuracy~\cite{leuchs2008}.
Let us briefly recall the core features of the interferometer
setup, details can be found in Ref.~\cite{leuchs2008}.
We use an interferometer of the Fizeau type at the standard
helium-neon laser wavelength $\lambda_{\text{HeNe}}=632.8$~nm.
A reference wave is produced at a wedged fused silica plate with
surface flatness $\le\lambda/$100.
A small steel (or glass) sphere located concentric to the focus of the
parabolic mirror is acting as a null element.
Any light ray traveling parallel to the optical axis of the parabola
is reflected towards the sphere, reflected back onto its
path of incidence and interferometrically overlapped with the
reference beam.
Owing to this double pass geometry, we measure twice the phase
aberration introduced by the mirror.
Possible aberrations of the sphere are eliminated by averaging over
different spheres and different rotational orientations of each sphere
relative to the parabolic mirror.
The deviations from a perfect parabola obtained after removal of
misalignment aberrations are finally represented by Zernike
polynomials.

Our parabolic mirrors are made of diamond turned aluminum and
are produced at the Fraunhofer Institute for Applied Optics and
Precision Engineering at Jena (see
Sec.~\ref{Generation_optimum_doughnut} for the geometric parameters).
A central on axis bore of 1.5~mm diameter is used for locating the
sphere or any other object such as an ion~\cite{maiwald2009} in the
focal point.

The interferometric measurements are performed with a radially
polarized doughnut beam centered on the optical axis of the parabolic
mirror.
In this particular case light is polarized parallel to the
plane of incidence at every spot on the surface of the parabolic reflector.

\begin{figure*}
\centering
\includegraphics{Golla_fig4.pdf}
\caption{
\label{fig4}
(a) Measured phase aberrations in single pass reflection. The
resulting aberrations are shown averaged over two different
spheres and six rotational settings per sphere.
The PV value is 0.53 waves.
(b) Measured aberrations with a compensating phase plate (single
experimental run).
The PV value is 0.45 waves.
In both panels the data result from a Zernike fit of degree 10.
Adjacent contour lines are separated by 0.05 waves at
$\lambda_\textrm{HeNe}$.
}
\end{figure*}

Figure~\ref{fig4}a) shows the result of such a measurement
for one of our parabolic mirrors.
The peak-to-valley (PV) value of the aberrations is 0.53 waves at
$\lambda_{\textrm{HeNe}}$.
The root mean square (RMS) value is 0.09 waves at
$\lambda_{\textrm{HeNe}}$.
The data shown correspond to maximum deviations of about 130~nm from
the ideal parabolic shape~\cite{leuchs2008}.
This is of the same order of magnitude as reported for the parabolic
mirror used in Ref.~\cite{drechsler2001}.

Based on these interferometry results we have produced a
phase plate which nominally imprints the phase conjugate of the
measured aberrations onto an otherwise flat phase front.
The fused silica phase plate is produced in house by laser lithography
and (selective) reactive ion etching.
In order to check the quality of aberration compensation the phase
plate is inserted into the interferometric setup between the parabolic
mirror and the flat reference plate.
While constantly repeating the interferometric measurement, position and
orientation of the phase plate are optimized.
As a first optimization criterion we take the reduction of the
PV value of the aberrations while for the fine adjustment the RMS
value is a better measure, since this is the quantity determining the
Strehl ratio~\cite{born-wolf1991}.
The final outcome of such a measurement series is depicted in
Fig.~\ref{fig4}b).
The PV value of the aberrations is reduced to 0.45 waves at
$\lambda_{\text{HeNe}}$.
The RMS value shrinks to 0.033 waves at $\lambda_{\text{HeNe}}$.
It should be noted that the seemingly small reduction of the PV value
is due to effects occurring at the boundaries of the phase plate and
the parabolic mirror.
Neglecting the outer 0.5 mm of the aperture the PV value shrinks to
0.16 waves at $\lambda_{\text{HeNe}}$ and the RMS value is reduced to
0.023 waves.

To judge the quality of these improvements we have performed
simulations of the focal intensities based on a generalization of the
method by Richards and Wolf~\cite{richards1959}. We assume a wave
front according to the aberrations shown in Fig.~\ref{fig4}
to be incident on a perfectly reflecting parabolic mirror.
The incident light wave is a radially polarized doughnut mode (see
Eq.~(\ref{eq_rad_pol_doughnut})) with the optimum beam radius.
The resulting focal intensities are normalized to the focal intensity
obtained for an aberration free wave front.
This yields the Strehl ratios.
At the testing wavelength of our interferometer (633~nm), the Strehl
ratio is improved from 68\% to 99\% by using the correcting phase
plate.
In other words, the combined system 'parabolic-mirror-plus-phase-plate'
is practically capable of diffraction limited focusing.

At the wavelengths for experiments with Ytterbium ions the Strehl
ratios of the uncorrected mirror are 32\% for 370~nm and 9\% for
252~nm.
Close to diffraction limited performance seems feasible in
these cases if the phase plate is designed for the respective
wavelength considering the corresponding refractive index
of fused silica.

Besides that, a phase plate that perfectly compensates the aberrations
at the testing wavelength $\lambda_{\text{HeNe}}$ is also rather
suitable at the wavelengths of the transitions T1 and T2.
Assuming a phase plate working perfectly at the testing wavelength, we
simulated the performance of the compound system
'phase-plate-plus-paraboloid' at the other two wavelengths.
Since the relative refractive index changes of fused silica at the
wavelengths of interest are on
the order of 1\% -- 3\%, our simulations reveal Strehl ratios ranging
from 97\% to more than 99\%, depending on the wavelength and the
specific parabolic mirror.

So far we have assumed a phase shift upon
reflection that does not depend on the distance to the mirror's
optical axis.
But due to the complex refractive index of aluminum this assumption
is generally not justified.
The radially varying phase shift upon reflection imposes a defocus.
Nevertheless we argue that the influence of the refractive index of
aluminum can be neglected.
We calculate the phase shift upon reflection by taking the respective
refractive indices from Ref. \cite{smith1998} and determining the
argument of the complex reflection coefficient for electric fields
parallel to the plane of incidence.
For this we use the relations between angle of incidence and
radial distance to the optical axis given in
Ref.~\cite{lindlein2007}.
The obtained phase shifts are then used in simulations with an
otherwise aberration free paraboloid.
For the wavelengths considered here the simulations show an axial
shift of the intensity maximum by less than a tenth of a wavelength.
The maximum intensity is reduced by amounts on the order of 0.2\% in
the case of transition T2 and even less for the other wavelengths.
The intensity in the nominal focus is reduced by
0.5\% (1.5\%, 2.8\%) at the testing (transition T1, transition T2)
wavelength.

The phase shifts induced by the above effects can in principle be
measured by our interferometric setup.
However, a substantial part of the treatment of the interferometric
data is the removal of misalignment aberrations of the null sphere
from the exact focal position~\cite{leuchs2008}.
Since the monotonous radial variation of the phase shift upon
reflection is equivalent to a slight displacement of the null sphere,
this effect is removed in the measurement procedure.

We finish this section by noting that the orders of magnitude of the
experimental results shown are generally reproduced for other
parabolic mirrors, with the shape of the aberrations of the
uncorrected mirrors differing from device to device.

\section{Generation of the optimum temporal mode}
\label{Temporal_pulse_shaping}

To design the temporal mode of a single photon wave packet that
would be absorbed efficiently by a two level system, we exploit the
time reversal symmetry of the compound system atom and free space
field modes~\cite{sondermann2007}. An excited two level atom
spontaneously decays into the ground state under emission of a photon
with a probability amplitude that is exponentially decreasing in time
\cite{weisskopf1930}.
If a photon wave packet with such a
decreasing envelope is sent onto a two level system the maximum
possible excitation probability is  calculated to be
54\%~\cite{stobinska2009}.

An excitation probability of 100\% is achieved
by generating a photon wave packet with the time reversed shape of the
spontaneously emitted photon: an exponentially increasing field
amplitude with a time constant which is twice the lifetime of the
excited state.
Furthermore the ideal pulse length is infinite and has a sharp edge at
the end.
The excitation probability depends not only on a good overlap of the
time constant and the sharp edge at the end of the pulse, but also on
the total pulse length. In Ref.~\cite{stobinska2009} it is shown
theoretically that for pulse lengths of five lifetimes 99\% excitation
probability can be achieved.
These numbers are valid for single photon Fock states.
However, as outlined already in Sec.~\ref{introduction}, coherent state
pulses with average photon numbers $\ll$1 are a good approximation to
single photon wave packets.

In the following, we discuss the generation of these coherent state wave
packets.
The experimental apparatus for the generation of arbitrarily shaped
pulses is shown in Fig. \ref{fig1}.
A continuous wave laser beam is cut into pulses using an acousto-optic
modulator (AOM).
The intensity diffracted into the first order of the AOM depends on
the amplitude $U_0$ of the applied radio frequency (RF) signal as
$I\left(t\right)\sim\sin^{2}(U_0(t))$.
An exponentially increasing wave packet is obtained by modulating the
amplitude of the applied RF signal $U_{\text{RF}}(t)=U_0(t)\cdot
\sin(\omega_\text{RF} t)$ as
\begin{equation}
U_0(t)=\arcsin\left(\exp\left[\dfrac{t}{2\tau}\right]\right)
\end{equation}
where $\omega_\textrm{RF}=2\pi\cdot400$~MHz is the RF driving
frequency of the AOM.
The modulated RF signal is generated with a digital arbitrary waveform
generator and amplified before being applied to the AOM.
The first diffraction order of the AOM is coupled to a polarization
maintaining single mode fiber and sent to the setup for the spatial beam
shaping.
The attenuation to a mean photon number per pulse of much less
than one is done using a set of half-wave plates and linear
polarizers.
The temporal profile of these highly
attenuated laser pulses is measured using a photomultiplier tube with
a quantum efficiency of 15\% and a frequency counter in start-stop
configuration.
Since we do not measure a photon in every pulse, the start trigger is
given by the detection of a photon and the stop trigger by the
following clock signal activating the output of the arbitrary waveform
generator.
We account for this trigger configuration by reversing the
temporal coordinate of the experimental data.
In Fig.~\ref{fig5} the measurements of the pulses designed to match
the lifetimes of transition T1 and T2 are shown.
The mean photon number per pulse is 0.09.
For simplicity reasons the 370~nm laser was used for the creation
of the pulses for both time constants.

\begin{figure}
\begin{center}
\includegraphics{Golla_fig5.pdf}
\caption{ \label{fig5}
  Generated exponentially increasing pulses with a mean photon
  number per pulse of 0.09~photons. The dashed line indicates the
  perfect pulse shapes. The time constant (life time) is 8~ns in panel
  (a) and 230~ns in panel (b).
  The total number of counts are $430\cdot10^3$ and
  $487\cdot10^3$, respectively. }
\end{center}
\end{figure}

The temporal overlap $\eta_{\text{t}}$ of the field amplitudes is
calculated according to Eq. (\ref{eq:eta_t}).
The incident field amplitude $E_\textrm{inc}(t)$
is obtained by taking the square root of the photon count statistics
and the integration is performed over the time intervals covered by
the experimental data shown in Fig. \ref{fig5}.
The results of this procedure are $\eta_t=$0.96 for the pulse
with 8~ns time constant and $\eta_t=$0.99 for the 230~ns pulse.
The lower overlap value in the case of the 8~ns pulse is attributed to
the finite build-up/decay time (5~ns) of the optical grating inside
the AOM crystal.
This results in non-negligible deviations from a sharp pulse end when
compared to the time constant of the rising exponential.

\section{Discussion}
\label{discussion}

\begin{table}
\begin{tabular}{l|c|c}
	Transition  & T1 & T2 \\ \hline \hline
	Spatial overlap $\eta$ & 0.979 & 0.975  \\ \hline
        Solid angle $\Omega_\mu/\frac{8\pi}{3}$ & 0.94 & 0.94 \\ \hline
	Strehl ratio  & 0.99 & 0.99 \\ \hline
	Temporal overlap $\eta_t$& 0.96 & 0.99 \\ \hline \hline
        Coupling strength $G$ & 0.892 & 0.885 \\ \hline
        Absorption probability $P_a$ & 0.812 & 0.867  \\ \hline
\end{tabular}
\caption{ \label{estimation_overlaps}
Overview of the achieved figures of merit.
 The Strehl ratio is taken from the proof of principle experiment at
 $\lambda_\textrm{HeNe}$.
}
\end{table}

With the acquired experimental data at hand we are now able to give an
estimate of the coupling efficiency and absorption probability which
would be achieved in experiments employing the generated field modes.
Table~\ref{estimation_overlaps} summarizes the field overlaps of the
generated light modes in time and in space domain, as well as the
Strehl ratio.
According to Eqs. (\ref{eq:G}) and (\ref{eq:Pabs}) we calculate $G$
and $P_a$ from these values and let the Strehl ratio enter as a
multiplicative factor accounting for losses caused by the
mirror imperfections remaining after compensation.
We thus arrive at prospective coupling strengths of $G=0.892$
($G=0.885$) for transition T1 (T2).
These values are only about two percent below the upper limit set by the
experimental geometry when using doughnut modes.
The corresponding absorption probabilities are $P_a=0.812$
($P_a=0.867$).

These absorption probabilities are valid under the assumption of a
closed two-level system (e.g., transition T2 of $^{174}$Yb$^{2+}$).
As an example for a more complex level structure we take the
transition T1 of $^{174}$Yb$^{+}$:
Both the excited and the ground state are split into two
Zeeman sub-levels with the magnetic quantum numbers
$m=\pm 1/2$ (we neglect the dipole allowed
$^{2}$P$_{1/2}\rightarrow\ ^{2}$D$_{3/2}$ transition due to the
small branching ratio of 0.5\% \cite{migdalek1980}).
Two linear and two circular dipole transitions are possible.
The probability for spontaneous decay via a linear dipole
transition is $1/3$ according to the corresponding Clebsch-Gordan
coefficient and $2/3$ via a circular dipole transition.
Exciting the linear dipole transition leads to a
maximum possible excitation probability of $1/3$, assuming the
excitation light mode is ideal.
In this sense, one has to understand the value $P_a=0.812$ given
for a mode designed for this transition as a relative value, i.e., the
de facto absorption probability would be $P_a/3$.
However, the coupling strength $G$ itself is not affected by such
considerations.
For example, an elastic scattering experiment on any of these
transitions should be feasible with the coupling strength derived for
an ideal two-level system because no dipole moments are induced for
transitions other than the driven one (cf. Ref. \cite{jin2011} for a
calculation applying the level structure discussed above).

We conclude this paper by emphasizing that the method presented here
to generate a dipole mode can be applied to a dipole transition at any
optical wavelength.
An application to circular dipole transitions is
straightforward requiring corresponding amplitude and polarization
distributions.
However, these polarization patterns become quite complex -- especially
when focusing from a large solid angle \cite{sondermann2008}.

In summary, we have demonstrated the generation of modes of the
electromagnetic field allowing for efficient coupling of single
photons and single two-level systems in free space using a
deep parabolic mirror.
The free space coupling method described here should compete
favorably with other coupling schemes.

\acknowledgments
The authors acknowledge fruitful discussions with
A. Villar, S. Heugel, M. Fischer and J. Schwider.
M.S. acknowledges financial support from the Deutsche
Forschungsgemeinschaft (DFG).


\section{Introduction}
\label{Introduction}
The efficient market hypothesis assumes that financial markets discount immediately all information available about the listed assets. In general the flux of information is both endogenous and exogenous and this leads to a classification of different forms of efficiency. It is natural therefore that a large empirical literature exists with the aim of measuring the role of exogenous and endogenous sources of information in explaining price dynamics, both in absolute and in relative terms. The seminal work of \cite{Cutler1989} started a stream of research trying to connect exogenous news with price movements. As detailed in the next literature review section, more recent papers have investigated stock price reaction to news, e.g. \cite{Chan2003,Vega2006}, the correlation between high/low pessimism of media and high market trading volume \citep{Tetlock2007}, the relation between the sentiment of news, earnings and return predictability \citep{Tetlock2008}, the correlation between the volume of searching of news in the Google searching engine and many financial indicators of stocks \citep{Da2011}, the role of news in the trading action of short sellers \citep{Engelberg2010}, the role of macroeconomic news in the performance of stock returns \citep{Birz2011}, and the high frequency market reaction to news, \citep{Joulin2008,Gross2011}. All these papers are concerned with the relation between news and price movements. A different and less explored stream of research, to which the present paper belongs, investigates the role of news on the trading and investment decisions of single investors. The main difficulty of this type of research is the availability of micro data about the activity of single investors. Recently, some papers have investigated how news affect the selection of stocks performed by single investors \citep{Barber2008}, the role of individual investor decisions in causing post-earning announcement drift \citep{Hirshleifer2008},  and the relation between high news attention and the level of trading of single investors \citep{Yuan2008}.

Financial markets are extremely heterogeneous systems, and investors represent an important source of heterogeneity. Investors are different in many respects, including their risk profile, the size of their investment, the regulatory constraints to which they are subject, the information they have access to, etc. One clear form of heterogeneity among investors, which incorporates most of the aforementioned differences, is the {\it category} to which an investor belongs. A large financial institution is clearly different from an household, or from a governmental institution and this difference might be reflected in the way each investor reacts to exogenous (news) or endogenous (price returns or volatility) factors. The concept of category is related to a specific classification and it is not unique. In this paper, we will use a classification given by the data we are using, which allows us to discriminate between non-financial corporations, financial and insurance corporations, general governmental organizations, non-profit institutions, households, and foreign organizations.

More specifically, in the present study, we investigate how different categories of single investors react to exogenous and endogenous factors by looking at their trading activity, both in terms of being active in the market and in terms of the decision to buy or to sell an asset.
To this aim we make use of two very detailed datasets. We consider the Nokia stock and we have access to a database containing the trading activity of all the investors whose financial ownership is recorded by the Finnish Central Securities Depository.  These data allow us to classify single investors in terms of the category mentioned above, and to characterize the buying and selling activity on Nokia with a daily resolution. For the exogenous news we consider all the Thomson Reuters news released during the time period 2003-2008 and containing information about the company Nokia.

The number of daily news about Nokia gives us a signal about the intensity of exogenous information without interpreting the content of the news.
In order to have a semantic interpretation of the each news and to classify it in terms of good or bad news we use the General Inquirer (http://www.wjh.harvard.edu/~inquirer/), a well-known content analysis program which is using the General Inquire categories from the Harvard psychosocial dictionary. We construct a simple proxy of the sentiment of news arriving into the market by applying the General Inquirer and measuring the absolute or relative difference between positive and negative words in the headline.

The main analysis of the paper is a linear regression and partial correlation analysis which allows us to assess for each category of investors the absolute and relative role of endogenous (price return and volatility) and exogenous (number of news and sentiment indicator) factors in explaining the decision to trade and, when they trade, the decision to buy or to sell. As detailed below we find a different behavior among different categories of investors. Compared to other categories, governmental and  non profit institutions are in absolute terms less affected by endogenous and exogenous factors. In relative terms, they are more affected by news than by price dynamics. On the contrary, trading action of households and non financial companies is significantly explained by the regression, but returns and volatility are more important than exogenous news. Finally, we show that financial institutions and foreign institutions display an intermediate behavior, but endogenous factors are more important than exogenous ones.

The paper is organized as follows:
In Section II we review the empirical literature on news, price dynamics, and trading activity. In Section III we describe the databases used in our study. In Section IV we introduce the variables used to characterize the trading activity of the investors and the proxies used for the flux of news, volatility and sentiment indicators. Section V presents and discuss the results obtained from the regression analysis, and Section VI concludes.

\section{Literature review}
 \label{Lr}

There is a vast literature about the role of news in financial markets. We can divide the literature in two streams of research. The first and older one considered the problem of how news affect asset price. The second stream of research, more closely related to the present paper, considered the problem on how single investors react to news announcements. This type of analysis has been possible only starting recently because of the availability of large datasets with records of the trading history of single investors.
In this section we review these two streams of literature.

Many papers have investigated price reaction to news since the pioneering work of \cite{Cutler1989} that estimated for the first time the fraction of the variation that can be attributed to economic news in aggregated stock returns. A few years later, \cite{Ederington1993} studied the impact of scheduled macroeconomic news announcements on interest rate and foreign exchange futures markets. \cite{Engle1993} investigated how new information is incorporated into volatility estimates in the presence of asymmetry in the impact of news. \cite{Mitchell1994} studied the relation between the daily number of Dow Jones news and aggregate measures of market activity such as trading volume and market returns.

Starting from 2003, new studies using comprehensive databases of news appeared in the literature. \cite{Chan2003} showed that stocks experiencing negative returns concurrent with the arrival of a news story continued to underperform their peers. The same public news database was successfully used by \cite{Vega2006}, together with the estimation of the probability of private information-based trading, to empirically measure the effect of private and public information on post-announcement drift. The sentiment carried by news impacting the market was first investigated by using the daily content from a popular {\it Wall Street Journal} column \citep{Tetlock2007}. In his study, Paul Tetlock found that news with negative sentiment predicts downward pressure on market price followed by a reversion to fundamentals. In a successive study, \cite{Tetlock2008} examined whether quantitative measures of sentiment of the text of news can be used to predict individual firms' accounting earnings and stock returns. In their study they concluded that linguistic content of news captures otherwise hard to quantify aspects of firms' fundamentals that are quickly reflected into stock prices. Another study investigated the role of dissemination of information on security pricing \citep{Fang2009}, showing that stocks with no media coverage earn higher returns than stocks with high media coverage. The role of investors' attention was also considered from the different perspective of information demand in a study of the Google Search Volume Index \citep{Da2011}. In this study, authors related the Google search Volume Index to a sample of Russel 3000 stocks showing that an increase in the Search Volume Index predicts higher stock prices in the next two weeks and an eventual price reversal.

The role of news in short sales was investigated by \cite{Engelberg2010}, who found that a negative relation between short sales and future returns is order twice larger on news days than on days without a significant flux of news. The relationship is of the order of four times on days with negative news. The analysis of a large electronic database of news allowed to investigate the role of news in high frequency trading on both volatility \citep{Joulin2008} and price formation and book dynamics \citep{Gross2011}. \cite{Joulin2008} found that volatility patterns around market endogenous jumps and around exogenous news are quite different with endogenous jumps followed by increased volatility and news triggering periods of lower than average volatility. The role of high frequency sentiment indicators on future price trends and bid-ask spreads has been studied by considering the high frequency price evolution of twenty stocks traded at the London Stock Exchange  \citep{Gross2011} and the Reuters NewScope Sentiment Engine, which is a pre-processed set of news data and electronic tools analyzing textual information using linguistic pattern recognition algorithms. Recent studies have also considered the role of specific categories of news, such as macroeconomic news \citep{Birz2011}, and the role of investor sentiment on the market's mean-variance trade off \citep{Yu2011}.

The literature on the role of news on trading decision and activity of single investors is more limited. \cite{Barber2008} tested and confirmed the hypothesis that individual investors are net buyers of stocks frequently discussed in the news. They proposed a model of decision making in which individual investors consider primarily those stocks having attention-grabbing qualities and preferential selection among them is exercised only after attention has limited the choice set. The response of individual investors to news was investigated by Yuan in a study that considered the trading and position information of 78,000 households investing in US markets from January 1991 to December 1996 \citep{Yuan2008}. In his study, Yuan showed that the impact of attention is pervasive across market. High attention causes individual investors to reduce stock positions in good times and moderately increase stock positions in bad times. His results also indicate that attention is one source of the cost of monitoring portfolios and that investors sensitive to news are more subject to the disposition effect. Another investigation on the behavior of individual investors in the presence of public news studied the role of individual investors in causing post-earnings announcement drift \citep{Hirshleifer2008}. Authors found that individuals are net buyers after both negative and positive bold earning announcements.

\section{Data}
\label{Data}
In this paper we investigate the database maintained by the Euroclear Finland (previously Nordic Central Securities Depository Finland). The database is the central register of shareholdings for Finnish stocks and financial assets in the Finnish Central Securities Depository. Practically all major publicly traded Finnish companies have joined the register. The register reports the shareholdings of all Finnish investors and of all foreign investors asking to exercise their vote right. Both retail and institutional investors are included. The database records official ownership of companies and financial assets and the trading records are updated on a daily basis according to the Finnish Book Entry System. The records include all the transactions, executed in worldwide stock exchanges and in other venues, which change the ownership of the assets.

The database classifies investors into six main categories: non-financial corporations, financial and insurance corporations, general governmental organizations, non-profit institutions, households, and foreign organizations.
The database is collected since January 1, 1995.  In the present study we investigate the market activity of investors trading the Nokia stock, which was, across the years under investigation, either the most capitalized stock or one of the most capitalized stocks in the Finnish stock market.

While the database contains very detailed information about the Finnish domestic investors, foreign investors can choose to use nominee registration. In this case, the investor's book entry account provider aggregates all the transactions from all of its accounts, and a single nominee register coded identity contains the holdings of many foreign investors.
This means that our results describe in a detailed way the actions of all the Finnish domestic investors and those foreign investors who do not use nominee registration, while a very small fraction of the coded identities correspond to a large aggregated ownership.

For this reason, in the present study, we consider only the set of single investors trading the Nokia stock during the period of time from January 2, 2003  to December 30, 2008 (a set of $1,510$ daily records) and we investigate all the market transactions performed by them. Single investor means here a retail or an institutional investor that do not use nominee registration (essentially all the Finnish investors). The total number of investors is $141,190$ and the total number of transactions is $7,494,104$. Table \ref{TableDescr1} reports the number of investors, the number of transactions, and the traded volume for the six categories.

\begin{table*}
\caption{\label{TableDescr1} Summary of the number of investors (\# ids), the number of transactions ($N$), and the exchanged volume ($V$) for the Nokia single investors in the period Jan. 2, 2003 - Dec. 30, 2008.  Volume is given in millions of shares. The investors are divided in the six categories. Nominee registered investors are not considered. Note that for transactions between two single investors the volume is counted twice, once for the buyer and once for the seller.}
\begin{center}
\begin{tabular}{|l|rrr|}
  \hline
  Category				& \# ids &$N$~~~ 	&$V$~~ 	\\
  \hline
Companies			 	&8,396	&	 1,009,226 	& 4,825		\\
 Financial 			 		&392		&	 4,079,174 & 21,402	 	\\
 Governamental			&124		&	 39,278	& 1,985		\\
 Non profit					&922		&	 21,778 & 248		\\
 Households				&129,952	&	 1,555,096	& 1,993		\\
 Foreign 					&1,405	&	789,552	 	& 7,685		\\
 \hline
Total						&141,190	&	 7,494,104 & 38,138		\\
\hline
\end{tabular}
\end{center}
\end{table*}

In this paper we investigate the relation between the trading  of Nokia investors, the Nokia price dynamics, and the flux of news about Nokia. As a source of news reaching financial markets worldwide we use the Headlines of the NewsScope archive of news released in English by Thomson Reuters during the investigated time period. Specifically, from the complete NewsScope archive we have extracted all headlines in English language labeled with at least one Nokia Reuters Instrument Code \footnote{The RICs used to extract the headlines  are NOK.W, NOK1V.HE, NOK1V.AS, NOKN.MX, NOKA.BA, NOKy.BE, NOK.MW, NOKy.F, NOK1VEUR.VIp, NOK1VEUR.Ip, NOK1VEUR.STp, NOKS.HA, NOKS.H, NOKS.DE, NOKy.D, NOKAc.BA, NOKy.MU, NOKy.DE, NOK1VM0110.HE, NOK1VEUR.PZ, NOK, NOK1VEUR.DEp, NOKI.ST, NOKS.BE, NOKS.F, NOK.DF, NOK.N, NOKAd.BA, NOK.P, NOK.C, NOK1VEUR.MIp, 0HAF.L, NOK1VEUR.PAp, NOK1V.MI ,NOKS.D, NOKS.MU}. The set comprises 11,484 unique headlines. Each headline is associated with one or more release time (multiple releases of the same headline are frequent). In case of multiple releases of the same headline we use as time of the headline the time of the first release.
\begin{figure}
\begin{center}
\includegraphics[scale=0.45,angle=0]{Fig1news2.pdf}
\caption{Average daily pattern of the arrival rate of news on the Nokia company. The rate is measured in number of headlines per minute. The vertical lines indicate the time of opening and closing of European and New York Stock Exchange market. Data are adjusted for the daylight saving time.}
\label{DFhead}
\end{center}
\end{figure}
In Fig. \ref{DFhead} we show the average daily pattern of the arrival rate of  Nokia news per minute as a function of the time of the day. Time is computed in coordinated universal time (UTC) and is corrected for the setting of daylight saving time in UK and for the difference between UK and US daylight saving time. The figure shows that news start to arrive at an high rate around 5.00 am (UTC time) and the distribution is roughly flat until 4.30 pm, which is the time of market closing in Europe. Spikes of the probability density function are observed around the time of market opening (8.00 am) and closing (4.30 pm) in Europe and opening (2.30 pm) and closing (9.00 pm) of US NYSE and NASDAQ markets. In the figure, opening and closing of markets are indicated with vertical lines.

We assume that the largest majority of the Finnish investors are trading in European markets and for this reason, in the present study, we consider only the headlines reaching financial markets during European trading hours (from 8.00 am to 4.30 pm UTC time). Note that a clear peak of news arrival is observed around the closing time of European markets (4.30 pm). Since most of these news are probably telling how the market is closing, we repeated the analysis by removing the last 10 minutes of the trading day, but the results presented below remain essentially unchanged.

\section{Variables characterizing trading activity, price dynamics, and the flux of news}\label{categorical}
\label{Variables}
In our analysis we consider three sets of variables, one containing variables describing the trading action of the investors, one describing the price dynamics, and one describing the news feed. In this section we define the variables and describe some of their statistical properties.

\subsection{Definition of the variables}

The first set of variables characterize the trading activity of single investors belonging to different categories.
The high degree of heterogeneity of investors in the frequency and volume of trading makes difficult to compare, for example, the activity of an household trading small volumes once every three months with the one of a financial institution that trades every day large volumes. Since we are primarily interested in comparing the impact of news on the daily trading of single investors, we use categorical variables that describes their trading activity. Similarly to what \cite{TumminelloNJP} did, we use the daily categorical variables of Selling investors (S), Buying investors (B) and Buying and Selling investors (BS).

The classification is obtained as follows: for each investor $i$ and each trading day $t$, we consider the Nokia volume sold $V_s(i,t)$ and the Nokia volume purchased $V_b(i,t)$ by the investor $i$ in that day. This information is then converted into a categorical variable with 3 states: primarily buying \emph{B}, primarily selling \emph{S}, buying and selling approximately closing the position \emph{BS}.  The conversion is done by using the ratio
\begin{equation}
q(i,t)=\frac{V_b(i,t)-V_s(i,t)}{V_b(i,t)+V_s(i,t)}.
\end{equation}
We assign an investor a primarily buying state \emph{B} when $q(i,t)>\theta$,  a primarily selling state \emph{S} when $q(i,t)<-\theta$, and a buying and selling state \emph{BS} when $ -\theta \le q(i,t) \le \theta$  with $V_b(i,t)>0$ and $V_s(i,t)>0$. When $V_b(i,t)=V_s(i,t)=0$ we consider the investor not active on day $t$. According to \cite{TumminelloNJP}, in the present study, we set $\theta=0.01$. Roughly speaking, investors in a buy (sell) state can be seen as acting as net buyers (sellers), while investors in a buy and sell state can be thought as intermediaries or day traders.

We use the categorical variables associated to each investor $i$ for each trading day $t$ to compute the time evolution of the number of investors of a given category $K$ performing a specific trading action (buying, selling or buying and selling). Specifically, $N_B^{K}(t)$ is the number of investors of category $K$ classified as buyers at day $t$, $N_S^{K}(t)$ is the number of investors of category $K$ classified as sellers at day $t$, and $N_{BS}^{K}(t)$ is the number of investors of category $K$ classified as buying and selling at day $t$.

From these variables we obtain the derived variables
\begin{eqnarray}
N^{K}(t)=N_B^{K}(t)+N_S^{K}(t)+N_{BS}^{K}(t) \\
\Delta N^K_{A} (t)=N_B^{K}(t)-N_S^{K}(t)\\
\Delta N^K_R(t)=\frac{N_B^{K}(t)-N_S^{K}}{N^{K}(t)}
\end{eqnarray}
The variable $N^K(t)$ quantifies the number of trading investors of category $K$ without discriminating the nature of the trading action. The variables $\Delta N_{A}(t)$ and $\Delta N_R (t)$ quantify the polarization of the trading choices of investors of category $K$ towards a buying decision (positive values) or a selling decision (negative values) in absolute and relative terms, respectively. Note that $N^K(t)$ includes all the investors, including those in \emph{BS} state, whereas $\Delta N_{A}(t)$ and the numerator of $\Delta N_R (t)$ are calculated by using only those in \emph{B} and in \emph{S} state. However the denominator of $\Delta N_R (t)$ contains all the active investors.

Market price dynamics is quantified by considering daily return of Nokia stock traded at the Nordic Stock Exchange, i.e.
\begin{equation}
Ret(t)=\log P(t)-\log P(t-1)
\end{equation}
where $P(t)$ is the closing price at day $t$. We also consider a proxy of the daily volatility of Nokia at the Nordic Stock exchange defined as
\begin{equation}
Vol(t)=2 \frac{P_{max}(t)-P_{min}(t)}{P_{max}(t)+P_{min}(t)}
\end{equation}
where $P_{max}(t)$ and $P_{min}(t)$ are the highest and lowest price of Nokia at day $t$, respectively.

We finally consider two variables quantifying the flux of news about Nokia arriving at day $t$. First of all, we consider the number $H(t)$ of Nokia headlines released by Thomson Reuters between 8.00 am and 4.30 pm during trading day $t$. This number is a measure of the ``intensity" of news reaching the market in a given day. However headlines news and the associated stories may bear positive, negative or doubtful information. In analogy to \cite{Tetlock2008}, we quantify the sentiment carried by the headlines by constructing a sentiment proxy using the number of positive and negative words present in them. Positive and negative words are detected by using the General Inquirer (http://www.wjh.harvard.edu/~inquirer/), a well-known content analysis program which is using the General Inquire categories from the Harvard psychosocial dictionary. Once the number of positive ($G(t)$) and negative ($B(t)$) words contained in all the headline news at day $t$ have been computed, we determine the variables
\begin{equation}
S_{A} (t)= G(t)-B(t)~~~~~~~~S_R (t)=\frac{G(t)-B(t)}{G(t)+B(t)}
\end{equation}
giving the absolute and relative, respectively, sentiment of the news in a given day. To avoid spurious discretization effects in the calculation of $S_R (t) \in[-1,1]$ we require $G(t)+B(t)\geqslant5$ to compute the sentiment indicator. When $G(t)+B(t)<5$ we set $S_R (t)=0$. The time evolution of $S_R (t)$ is therefore different from zero only when a significant number of positive and negative words are detected in the headlines.

 \subsection{Descriptive analysis}
 \begin{sidewaystable}
\caption{Summary statistics of the investigated variables.}
\begin{tabular}{l|l||rrrrrrr}
\hline
Variable& Category& Min&	5\% quant.&	95\% quant.&	 Max&	Mean&	Median&	Std Dev\\
\hline
$N^K$    &Companies	  &0	&25	&171.55	&630	&73.96	&59	&57.3\\
            &Financial	  &4	&14	&31	&69	&20.29	&19	&5.798\\
            &Governmental &0 	&1	&12	&32	&4.713	&4	&3.769\\
            &Non profit	   &0	&1	&14	&43	&5.278	&4	&4.883\\
           &Households	   &0	&202	&1234.65	&4619	&538	&424	&423\\
           &Foreign	            &0	&3	&16	&47	&8.356	&7	&4.443\\
\hline
$\Delta N_A^K$ &Companies	   &-332	&-64.55	&83	&564	&4.826	&1	&49.8 \\
            &Financial	   & -25	&-11	         &10	& 45          &0.1576	&0	&6.835\\
            &Governmental &-23	&-6	         &7	& 32	         & 0.1781	&0	&4.127\\
            &Non profit	   &-35	&-8.55	&6	&25	         &-0.9609	&-1	&4.648\\
            &Households	   &-2143	&-494.05	&662.55	&4147	&26.5	&-15	&406.6\\
            &Foreign	            &-43	&-8	         &6	&28	          & -0.7152 &-1	&4.629\\
\hline
$\Delta N_R^K$& Companies &-0.8977	&-0.6475	&0.6917	&0.9224	&0.03161	&0.03252	&0.423\\
                         &Financial	    &-0.8125	&-0.4783	&0.5	&0.8947	&0.009647	&0	&0.3012\\
                      &Governmental &-1	&-1	&1	&1	&-0.003318	&0	&0.6574\\
                      &Non profit	      &-1	&-1	&1	&1	&-0.1573	&-0.2	&0.6798\\
                     &Households	      &	-1	&-0.7428	&0.7164	&0.913	&-0.02797	&-0.05142	&0.4637\\
                     &Foreign	      &-1	&-0.7778	&0.6667	&1	&-0.08613	&-0.1181	&0.4468\\
\hline
\hline
$Ret$	&& -0.1843	&-0.0365	&0.0347	&0.1361	&-0.000206	&0	&0.0234\\
$Vol$  &&	0.0050	&0.01007	&0.0609	&0.2296	&0.02637	&0.0207	&0.02041\\
\hline
\hline
$H$   &&	0	&0	&11	&46	&3.742	&3	&4.449\\
$S_A$&&	-12	&-2.55	&4	&23	&0.492	&0	&2.512\\
$S_R$&&	-1	&-0.2	&0.4286	&1	&0.01795	&0	&0.2172\\
\hline
\end{tabular}
\label{summarystat}
 \end{sidewaystable}

 \begin{figure}
\begin{center}
\includegraphics[scale=0.5,angle=0]{Fig2news.pdf}
\caption{From top to bottom the figure shows the time series of the number of Nokia headlines $H(t)$, the daily volatility $Vol(t)$ of Nokia stock, and the time series of $N^K(t)$ for the category of Financial investors and for the category of Households investors.}
\label{ExNa}
\end{center}
\end{figure}

Table \ref{summarystat} shows the summary statistics of the investigated variables.  We note that for the number of active investors $N^K$ the mean is always larger than the median, indicating a positive skew of the distributions. Moreover the standard deviations are typically quite large. The Table shows that for the governmental,  non profit, and foreign categories, the relative imbalance $\Delta N^K_R$ reaches the minimum and maximum value of $-1$ and $+1$, respectively.
The presence of these values indicates that, in some days, all the active investors of the considered category took the same market position. However, it is worth noting that the three categories presenting this behavior are the ones having an average number of active investors lower than ten (4.713, 5.278 and 8.356 for governmental, non profit and foreign category respectively). In other words the complete market polarization in most cases involves a limited number of single investors. The same is valid for the occurrence of the minimum value $-1$ observed for households category. This occurrence happened only once during a market session that involved only four households investors.

Figure \ref{ExNa} shows the time series of the number of Nokia headlines $H(t)$, the daily volatility $Vol(t)$ of Nokia stock, and the time series of $N^K(t)$ for the category of Financial investors and for the category of Households investors.  We note that the time series of $H(t)$ presents a background around 4 headlines per day and a series of spikes jumping to 20 headlines or more per day. The flux of news is not clustered in time, in fact the autocorrelation function is significant only at one lag.  The time evolution of the volatility proxy $Vol(t)$ shows that also this quantity is characterized by a typical value (of the order of 2 percent) and by days of huge swings with values of $Vol(t)$ of the order of $20\%$. The spikes of news and volatility are clearly correlated with the spikes of $N^K(t)$ for all the categories of investors. In the figure we show the time evolution of the number of financial and households investors. The overall behavior of the other categories is similar. It is worth noticing that the autocorrelation of $N^K$ is quite persistent and it is statistically significant at $2\sigma$ for more than 30 days.

\begin{figure}
\begin{center}
\includegraphics[scale=0.5,angle=0]{Fig3news.pdf}
\caption{From top to bottom the figure shows the time series of the difference $S_{A} (t)$ between the number of positive and negative words in Nokia headlines, of the daily return $Ret(t)$ of Nokia stock, and of  $\Delta N^K_{A} (t)$ for the category of Financial investors and of Households investors.}
\label{ExNbNs}
\end{center}
\end{figure}

In Fig. \ref{ExNbNs} we show the time series of the difference $S_{A} (t)$ between the number of positive and negative words in Nokia headlines, of the daily return $Ret(t)$ of Nokia stock, and of  $\Delta N^K_{A} (t)$ for the category of Financial investors and of Households investors.  The time evolution of $S_{A}(t)$ fluctuates around zero but also presents a series of positive and negative spikes jumping to the level of order 20 for positive or 10 for negative words per day.  The time evolution of Nokia return is characterized by a non Gaussian profile of the return probability density function and volatility clustering.
Spikes of $\Delta N^K_{A} (t)$ are also detected but they are in general less pronounced than in the case of $N^K(t)$ (see the bottom two panels of Fig. \ref{ExNa}) suggesting that the interpretation of news and/or endogenously extracted market information is usually different among investors of the same category. However, some pronounced spikes are still observed showing that in some occasion investors' categories take the same kind of trading action. Moreover autocorrelation analysis shows that for financial institutions $\Delta N^K_A$ is persistent only at one lag, while for households the autocorrelation function is statistical significant up to 30 days. This fact indicates the presence of persistent "moods" the household investors are following in their trading actions.

\begin{figure}
\begin{center}
\includegraphics[scale=0.5,angle=0]{Fig4news.pdf}
\caption{From top to bottom the figure shows the time series of the relative sentiment indicator  $S_R (t)$, the daily return $Ret(t)$ of Nokia stock, and the time series of $\Delta N^K_R (t)$ for the category of Financial investors and of Households investors. }
\label{ExDN}
\end{center}
\end{figure}

In order to investigate the difference between absolute and relative variables we show in Fig. \ref{ExDN} the time series of the relative sentiment indicator  $S_R (t)$, the daily return $Ret(t)$ of Nokia stock, and the time series of $\Delta N^K_R (t)$ for the category of Financial investors and of Households investors.   The time series of $S_R (t)$ presents a series of positive and negative spikes often clustered in time and value.  Quite interestingly, the time series of $\Delta N^K_R (t)$ does not present spikes but rather a noisy oscillation around zero. As in the case of $\Delta N_A^K$, the time evolution of households shows a slow dynamics lasting up to several trading months.

\section{Regression results}
\label{Results}
In this Section we use regression and partial correlation analysis to assess the role of endogenous or exogenous factors in determining the trading behavior of single investors. More specifically, we consider first of all the decision to trade (irrespectively of the specific position taken) and we regress it against an endogenous factor, namely volatility, and against an exogenous factor, namely the number of news. Similarly, in a second step we consider the imbalance (absolute or relative) between buyers and sellers and we regress it against contemporaneous return and the sentiment indicator. Also in this case, the first regressor can be considered endogenous and the second one as exogenous. Before presenting the results of the regression, two comments are in order. First of all, the distinction between exogenous and endogenous is not clear cut and in fact, as we will see, in both cases the regressors are not independent. Second, we shall consider contemporaneous variables and therefore it should be clear that no causality can be attached to our results. To be more explicit, the fact that we find a significant relation in the regression between number of investors and volatility does not necessarily imply that volatility triggers people to trade, but it can be the other way around, i.e. a large number of investors increases the volatility through their trading. Another possibility is that the two variables are influenced by a third unobserved factor.

\subsection{Volatility and  number of news}

Let us consider first how the decision to trade (irrespectively on being a buyer or a seller) is related with volatility $Vol$ and the flux of news $H$ for different categories of investors. Volatility and news are different sources of information, one primarily endogenous and one primarily exogenous to the market, but they are not mutually independent. In fact the empirical correlation between the two variables is Corr$[H,Vol]=0.501$. To interpret the role of the two variables in trading, we consider the linear model
\begin{equation}
\label{linearNa}
\widehat{N}^K (t)=\alpha_{H} \widehat{H} (t)+ \alpha_{Vol} \widehat{Vol} (t)+\epsilon(t)
\end{equation}
where $\widehat{N}^K$, $\widehat{H}$, and $\widehat{Vol}$ are standardized versions with zero mean and unit variance of $N^K$, $H$, and $Vol$, respectively.

 \begin{sidewaystable}
\caption{Summary of the results of the linear regression of Eq. \ref{linearNa} of the number $N^K$  of trading investors versus the news intensity signal $H$ and the volatility proxy $Vol$. The number in parentheses are the 5\%-95\% confidence intervals under Gaussian hypothesis and by using bootstrap analysis. The last two columns show the results of the partial correlation analysis.}
\begin{tabular}{|l|c|c|c|c|c|}
\hline
Investor& $\alpha_{H}$&$\alpha_{Vol}$&\% variance& $\rho(N^K,H|Vol)$ &$\rho(N^K,Vol|H)$  \\
category&~&~&of residual of $N^K$&&\\
\hline
Companies&0.271 (0.229,0.313)&0.517 (0.475,0.559)&51.8 \%&0.309&0.534\\
bootstrap&~~~~~~~ (0.205,0.335)&~~~~~~~ (0.437,0.597)&~~~~&&\\
\hline
Financial&0.195 (0.149,0.242)&0.479 (0.433,0.526)&63.8 \%&0.207&0.461\\
bootstrap&~~~~~~~ (0.125,0.264)&~~~~~~~ (0.407,0.558)&~~~~&&\\
\hline
Governmental&0.238 (0.183,0.292)&0.192 (0.138,0.246)&86.0 \%&0.215&0.180\\
bootstrap&~~~~~~~ (0.164,0.303)&~~~~~~~ (0.119,0.262)&~~~~&&\\
\hline
Non profit&0.319 (0.269,0.369)&0.270 (0.220,0.320)&73.9 \%&0.305&0.264\\
bootstrap&~~~~~~~ (0.249,0.394)&~~~~~~~ (0.199,0.344)&~~~~&&\\
\hline
Households&0.226 (0.188,0.263)&0.627 (0.589,0.664)&41.4 \%&0.289&0.651\\
bootstrap&~~~~~~~ (0.165,0.285)&~~~~~~~ (0.554,0.697)&~~~~&&\\
\hline
Foreign org.&0.158 (0.109,0.207)&0.442 (0.393,0.492)&70.9 \%&0.160&0.416\\
bootstrap&~~~~~~~ (0.094,0.224)&~~~~~~~ (0.374,0.517)&~~~~&&\\
\hline
\hline
\end{tabular}
\label{NaCtg}
 \end{sidewaystable}

In Table \ref{NaCtg} we show the values of the $\alpha_{H}$ and $\alpha_{Vol}$ coefficients together with the variance of the residual obtained by ordinary least squares. In the Table we also report the 5\%-95\% confidence interval of each coefficient under two null hypotheses. The first is the customary assumption of Gaussian errors, while the second ones are obtained by bootstrapping the data and therefore taking into account the distributional properties of the data. We show the results for each category of investors separately. All the regression coefficients are statistically not consistent with zero (at the given confidence level).
The $\alpha_{H}$ coefficient is ranging from a minimum value of $0.158$, which is observed for the foreign organizations, to a maximum value of $0.319$, which is observed for the non profit organizations. The $\alpha_{Vol}$ coefficient is ranging from a minimum value of $0.192$, which is observed for the governmental organizations to a maximum value of $0.627$, which is observed for households. The variance of the residuals of $\widehat{N}^K$ is ranging from 41.4\% (households) to 86\% (governmental organizations) indicating that in most cases the explanatory value of the two variables is quite significant.

For four categories of investors (companies, financial institutions, households and foreign organizations), the $\alpha_{Vol}$ coefficient is higher than $\alpha_{H}$ indicating that for these investors the market endogenous information has on average an higher explanatory role than market exogenous information conveyed by news on their decision to trade. The case of governmental organizations and non profit organizations is different. For these categories $\alpha_{H}$ is greater than $\alpha_{Vol}$ but the two values are within the 5\%-95\% intervals of the estimated coefficients. Moreover, for these categories the variance of the residual assumes the maximal values which are observed being 86.0 \% for governmental organizations and 73.9 \% for non-profit organizations.

We complement the regression analysis by computing the partial correlation coefficients of the three variables $N^K$, $Vol$, and $H$. The partial correlation coefficient $\rho(x,y|z)$  between variables $x$ and $y$ conditioning on the variable $z$ is the Pearson correlation coefficient between the residuals of $x$ and $y$ that are uncorrelated with $z$. Partial correlation is clearly related to linear regression. However, the information one obtains from the analysis of regression coefficients is not identical to the one obtained by considering partial correlation. In fact, for normalized variables, the best linear fit $z=\alpha_x x+\alpha_y y+ \epsilon$ gives
\begin{equation}
\label{lineartopartial}
\frac{\alpha_x}{\alpha_y}=\frac{\rho(z,x|y)}{\rho(z,y|x)}\sqrt{\frac{1-\rho^2_{y,z}}{1-\rho^2_{x,z}}}
\end{equation}
where $\rho_{x,z}$ and $\rho_{y,z}$ are the correlation coefficients between $x$ and $z$ and between $y$ and $z$, respectively (please refer to Appendix A fort more details). Only when  $\rho_{x,z}=\rho_{y,z}$ the ratio between the regression coefficients is equal to the ratio between the partial correlations.

Table \ref{NaCtg} shows the values of $\rho(N^K,H|Vol)$ and $\rho(N^K,Vol|H)$ for the six categories of investors. Consistently  with the previous results, we notice that $\rho(N^K,H|Vol)>\rho(N^K,Vol|H)$ for Governmental and Non profit organization, while the reverse is true for the other categories of investors.

In conclusion, volatility and the flux of news are correlated with the decision to trade of single investors and their relative role, when properly disentangled, is different for different categories of investors. Governmental and non profit institutions are the categories for which news and volatility give less explanatory power of their presence in the market. Moreover, these investors are more sensitive to news than to volatility. Households and companies are much more sensitive to volatility than to news and the variance of their activity explained by these factors is quite high. Also for financial institutions, volatility is more important than news, but the variance is relatively smaller. Finally, foreign organizations are more affected by volatility than by news, but the variance of the regression is quite small.

\subsection{Returns and sentiment indicator}

By having verified that news play an important role in the decisions of single investors to trade, we now focus our attention on the impact of the sentiment carried by news on Nokia return and on the trading behavior of the single investor to buy or sell a certain amount of Nokia stock. As sentiment indicators we consider both the absolute sentiment indicator $S_{A}(t)$ and the relative sentiment indicator $S_R(t)$. Both daily sentiment indicators are correlated with the daily Nokia return of the Nordic Stock Exchange. The values of the correlation are Corr$[S_{A},Ret]=0.155$ and  Corr$[S_R,Ret]=0.118$. These values are small but statistically significant. In fact the average correlation observed for shuffled time series of $S_{A} (t)$ and $S_R (t)$ with Nokia return is smaller than $0.002$ with a standard deviation of $0.026$.

 \begin{sidewaystable}
\caption{Summary of the results of the linear regression of Eq. \ref{linearDNa} of the difference $\Delta N^K_{A}$   between buying and selling investors versus the absolute sentiment indicator $S_A$ and the stock return $Ret$. The number in parentheses are the 5\%-95\% confidence intervals under Gaussian hypothesis and by using bootstrap analysis. The last two columns show the results of the partial correlation analysis.}
\begin{tabular}{|l|c|c|c|c|c|}
\hline
Investor& $\alpha_{S_{A}}$&$\alpha_{Ret}$&\% variance of & $\rho(N^K_A,S_A|Ret)$ &$\rho(N^K_A,Ret|S_A)$\\
category&~&~&residual of $\Delta N^K_{A}$&&\\
\hline
Companies&0.041 (0.002,0.080)&-0.653 (-0.692,-0.614)&58.0 \%&0.0528&-0.6463\\
bootstrap&~~~~~~~ (-0.031,0.102)&~~~~~~~ (-0.749,-0.572)&~~~~&&\\
\hline
Financial&0.008 (-0.034,0.050)&-0.576 (-0.618,-0.534)&66.9 \%&0.0095&-0.5709\\
bootstrap&~~~~~~~ (-0.055,0.069)&~~~~~~~ (-0.651,-0.516)&~~~~&&\\
\hline
Governmental&-0.019 (-0.070,0.031)&-0.196 (-0.246,-0.146)&95.9 \%&-0.0196&-0.1940\\
bootstrap&~~~~~~~ (-0.083,0.047)&~~~~~~~ (-0.270,-0.125)&~~~~&&\\
\hline
Non profit&-0.026 (-0.076,0.024)&-0.220 (-0.270,-0.170)&94.8 \%&-0.0263&-0.2178\\
bootstrap&~~~~~~~ (-0.116,0.056)&~~~~~~~ (-0.303,-0.136)&~~~~&&\\
\hline
Households&0.060 (0.021,0.099)&-0.655 (-0.694,-0.616)&57.9 \%&0.0781&-0.6474\\
bootstrap&~~~~~~~ (-0.012,0.131)&~~~~~~~ (-0.745,-0.568)&~~~~&&\\
\hline
Foreign org.&-0.006 (-0.050,0.038)&-0.504 (-0.548,-0.460)&74.4 \%&-0.0072&-0.4999\\
bootstrap&~~~~~~~ (-0.072,0.058)&~~~~~~~ (-0.571,-0.438)&~~~~&&\\
\hline
\hline
\end{tabular}
\label{DNAbCtg}
 \end{sidewaystable}

\begin{sidewaystable}
\caption{Summary of the results of the linear regression of Eq. \ref{linearDNR} of the relative difference $\Delta N^K_{R}$   between buying and selling investors versus the relative sentiment indicator $S_R$ and the stock return $Ret$. The number in parentheses are the 5\%-95\% confidence intervals under Gaussian hypothesis and by using bootstrap analysis. The last two columns show the results of the partial correlation analysis.}
\begin{tabular}{|l|c|c|c|c|c|}
\hline
Investor& $\alpha_{S_R}$&$\alpha_{Ret}$&\%  variance of & $\rho(N^K_R,S_R|Ret)$ &$\rho(N^K_R,Ret|S_R)$\\
category&~&~&residual of $\Delta N^K_R$&&\\
\hline
Companies&0.055 (0.014,0.095)&-0.610 (-0.650,-0.569)&63.3 \%&0.0685&-0.6056\\
bootstrap&~~~~~~~ (0.015,0.100)&~~~~~~~ (-0.685,-0.548)&~~~~&&\\
\hline
Financial&0.018 (-0.025,0.062)&-0.520 (-0.564,-0.477)&73.1 \%&0.0212&-0.5170\\
bootstrap&~~~~~~~ (-0.030,0.064)&~~~~~~~ (-0.587,-0.463)&~~~~&&\\
\hline
Governmental&0.021 (-0.029,0.071)&-0.179 (-0.230,-0.129)&96.8 \%&0.0215&-0.1782\\
bootstrap&~~~~~~~ (-0.027,0.075)&~~~~~~~ (-0.225,-0.136)&~~~~&&\\
\hline
Non profit&0.025 (-0.025,0.075)&-0.175 (-0.225,-0.125)&96.9 \%&0.0256&-0.1738\\
bootstrap&~~~~~~~ (-0.028,0.079)&~~~~~~~ (-0.227,-0.130)&~~~~&&\\
\hline
Households&0.068 (0.026,0.110)&-0.565 (-0.608,-0.523)&68.4 \%&0.0811&-0.5615\\
bootstrap&~~~~~~~ (0.025,0.111)&~~~~~~~ (-0.629,-0.512)&~~~~&&\\
\hline
Foreign org.&0.030 (-0.017,0.077)&-0.400 (-0.446,-0.353)&84.2 \%&0.0323&-0.3970\\
bootstrap&~~~~~~~ (-0.015,0.076)&~~~~~~~ (-0.449,-0.354)&~~~~&&\\
\hline
\hline
\end{tabular}
\label{DNRCtg}
\end{sidewaystable}

We analyze the explanatory role of $S_{A}$ and of the Nokia return $Ret$ by considering the linear model
\begin{equation}
\label{linearDNa}
\widehat{\Delta}N^K_{A} (t)=\alpha_{S_{A}} \widehat{S}_{A} (t)+ \alpha_{Ret} \widehat{R}et (t)+\epsilon(t)
\end{equation}
where $\widehat{\Delta}N^K_{A}$, $\widehat{S}_{A}$, and $\widehat{R}et$ are standardized versions with zero mean and unitary variance of $\Delta N^K_{A}$, $S_{A}$, and $Ret$, respectively.

In Table \ref{DNAbCtg} we show the values of the $\alpha_{S_{A}}$ and $\alpha_{Ret}$ coefficients together with the variance of the residuals obtained by ordinary least squares. As in the previous case, we also report the 5\%- 95\% confidence interval of each coefficient under Gaussian hypothesis and bootstrap analysis and the results are shown for each category of investors separately.
Most of the $\alpha_{S_{A}}$ coefficients are consistent with zero within the 5\%-95\% confidence interval. In fact, for bootstrap confidence intervals all the $\alpha_{S_{A}}$ are consistent with zero, while under normality hypothesis values (slightly) different from zero are observed in the cases of companies and households investors.
The $\alpha_{Ret}$ coefficient is always statistically significant and it is ranging from a minimum value of $-0.653$ observed for households to a maximum value of $-0.196$ observed for governmental organizations. The variance of the residual of $\widehat{\Delta}N^K_{A}$ is ranging from 57.9\% (households) to 95.9\% (governmental organizations) indicating that still for some of the categories (companies, financial institutions, households and foreign organizations) the explanatory value of the two variables, and especially of the return, is significant.

The $\alpha_{Ret}$ coefficients are in several cases (companies, financial institutions, households and foreign organizations) negative large values indicating that, for these categories, the market polarization of trading actions is strongly anticorrelated with the Nokia return. The majority of single investors of these categories are therefore buying when the Nokia price goes down and selling when the price goes up. This is reminiscent of a contrarian behavior, even if in this case we are considering contemporaneous rather than lagged correlation between return and decision to buy or sell.
A similar behavior is also  observed, although to a less pronounced level, for governmental and non-profit organizations investors.
The partial correlations showed in Table  \ref{DNAbCtg} confirm this results (see below for an extended  comment) .

These results are essentially confirmed by the analysis concerning the explanatory role of the relative sentiment indicator $S_{R}$ and of the Nokia return $Ret$ for the relative variable $\Delta N^K_{R}$ investigated by considering the linear model
\begin{equation}
\label{linearDNR}
\widehat{\Delta}N^K_{R} (t)=\alpha_{S_{R}} \widehat{S}_{R} (t)+ \alpha_{Ret} \widehat{R}et (t)+\epsilon(t)
\end{equation}
where $\widehat{\Delta}N^K_{R}$, $\widehat{S}_{R}$, and $\widehat{R}et$ are standardized versions with zero mean and unitary variance of $\Delta N^K_{R}$, $S_{R}$, and $Ret$, respectively. By comparing the results of Table  \ref{DNRCtg} with the results reported in Table \ref{DNAbCtg}  we conclude that the observations done for $\widehat{\Delta}N^K_{A}$, $\widehat{R}et$ and $\widehat{S}_{A}$ are very similar to the ones obtained for the relative variables $\widehat{\Delta}N^K_{R}$, $\widehat{R}et$ and $\widehat{S}_{R}$, and therefore do not depend significantly on the specific definition of the sentiment indicator. The only difference we detect is observed for companies and households. It concerns the values of the coefficient  $\alpha_{S_{R}} $ which are not consistent with zero within the 5\%-95\% confidence interval.

The partial correlations shown in Table \ref{DNAbCtg} and \ref{DNRCtg} indicate that $|\rho(N^K_A,S_A|Ret)|$ $\ll$ $|\rho(N^K_A,Ret|S_A)|$ and $|\rho(N^K_R,S_R|Ret)|$ $\ll$ $|\rho(N^K_R,Ret|S_R)|$ for all the categories of investors. The highest values of partial correlations, $\rho(N^K_A,S_A|Ret)$ and $\rho(N^K_R,S_R|Ret)$, are observed for companies and households. However, in general, these correlations are quite small in absolute value (between $5\%$ and $8\%$), and very close to the noise level. For the other categories, these partial correlations are even smaller, and statistically not significant.

The joint analysis of absolute and relative imbalance between buyers and sellers leads us to the following conclusions. The activity of governmental and non profit organizations is very poorly explained by return and news sentiment. Of the two factors, return plays clearly a major role. Households and companies are those for which sentiment and returns have the best explanatory power of their trading action. Return is clearly more important, but sentiment has also some explanatory power, especially when one consider the relative imbalance between buyers and sellers. For financial and foreign organizations the variance explained by the regressions is somewhat intermediate between the two pairs of categories above, but in general returns have a much higher explanatory power and sentiment plays a negligible role.

In conclusion, the regression analysis shows that the total flux of news is significantly correlated with the decision to trade.
This seems to indicate that, on a daily time scale, news move investors to trade.  However, we find that most of the times the sentiment indicator is not significantly correlated with the imbalance between buyers and sellers. One possible explanation is that the sentiment indicator does not discriminate accurately a good news  from bad news. As mentioned above, sentiment indicator has a correlation with returns of value around $0.15$ and therefore it seems possible that a more accurate semantic indicator could better capture how good or bad is the news. Another possibility is that most of the times investors do not agree on the positivity or negativity of the news and/or that they react in an heterogeneous way to the same news.

\section{Conclusions}
\label{Conclusions}
By using a linear regression model to describe the number of investors trading the Nokia stock as a linear combination of a proxy of news and a proxy of stock volatility, we have shown that  the trading activity of all the investigated categories of investors is significantly correlated with both the flux of news and the daily volatility. The relative role of exogenous (news) and endogenous (volatility) factors is, in general, difficult to disentangle because the two variables are highly correlated. By assuming a linear model, we have shown that the relative relevance of the two variables changes across the different categories of investors. The dependency from volatility turned out to be more pronounced than the one from news, for companies, financial institutions, households, and foreign organizations.

The second conclusion concerns the relationship between the sentiment of news, Nokia return, and the market polarization towards buying or selling for the different investors' categories. Our results show that the sentiment time series is correlated with Nokia return. We have also shown that both the Nokia return and the sentiment of news explain part of the trading polarization dynamics of some categories, according to a linear model. Specifically, we observed that all categories have a contrarian like response to return, and that companies and households have a positive, although small sized, reaction to positive sentiment of news.\\

{\bf Acknowledgments} Authors acknowledge Thomson Reuters for providing the NewScope archive data and Euroclear for granting access to the ownership data of single investors. FL acknowledges partial support by the grant SNS11LILLB ``Price formation, agent's heterogeneity, and market efficiency". JP acknowledges Magnus Ehrnrooth Foundation, and Jenny and Antti Wihuri Foundation for financial support.