\section{Introduction}

Given a sequence $\{z_n\}\subset{\mathbb D}$ subject to the Blaschke
condition
\begin{equation}\label{blacon}
\alpha:=\sum_{n=1}^\infty (1-|z_n|)<\infty,
\end{equation}
let
\begin{equation*}
B(z)=\prod_{n=1}^\infty b_n(z), \qquad b_n(z)=\frac{\overline
z_n}{|z_n|}\frac{z_n-z}{1-\overline z_n z}
\end{equation*}
be a Blaschke product with the zero set $Z(B)=\{z_n\}$. With no loss
of generality we will assume that $B(0)\not=0$.

One of the central problems with Blaschke products is that of the
membership of their derivatives in classical function spaces, $B'\in
X$. There is a vast literature on the problem, starting from
investigations of P. Ahern and his collaborators \cite{A79, A83,
AC74, AC76} and D. Protas \cite{P73} in 1970s, up to quite recent
results of the Spanish school \cite{GP06, GPV07, GPV08}, see also
\cite{MS09, P10}. The above mentioned spaces $X$ are primarily the
Hardy spaces $H^p$, the Bergman spaces $A^p$, the Banach envelopes
of the Hardy spaces $B^p$ etc. Recall the definition of the Bergman
spaces $A^p$, $p>0$:
$$ A^p=\{f\in{\mathcal{A}}({\mathbb D}): \ \ \int\limits_{\mathbb D} |f(z)|^p\,dxdy<\infty\},
\quad z=x+iy.
$$

In this paper we will add to the list some new spaces
$X={\mathcal{A}}(E,\rho)$ of analytic functions in the unit disk, introduced
recently in \cite{FG}. Given a closed set $E=\overline E\subset{\mathbb T}$ and
$\rho>0$, we say that an analytic function $f$ belongs to
${\mathcal{A}}(E,\rho)$ if
\begin{equation*}
|f(z)|\le C_1\exp\left(\frac{C_2}{d^\rho(z,E)}\right)\,, \quad
d(z,E)={\rm dist}(z,E)
\end{equation*}
is the distance from $z\in{\mathbb D}$ to $E$, $C_{1,2}$ are positive
constants.

The simplest and most general result drops out immediately from the
Schwarz--Pick lemma for functions $g$ from the unit ball of
$H^\infty$:
$$ |g'(z)|\le\frac{1-|g(z)|^2}{1-|z|^2}\,, $$
and states, that $g'\in A^p$ for all $0<p<1$. The result is sharp:
there exists a Blaschke product $B$ such that $B'\notin A^1$ (W.
Rudin).

To proceed further, one should impose some additional restrictions
either on absolute values $|z_n|$, stronger than \eqref{blacon}, or
on location, distribution of arguments of zeros $z_n$ etc. For
instance, a typical result of the first type is due to Protas
\cite{P73}:
$$ \sum_{n=1}^\infty (1-|z_n|)^r<\infty, \quad 0<r<\frac12 \ \
\Rightarrow \ \ B'\in H^{1-r}, $$ and $\frac12$ is sharp.

We are more interested in the second direction, related to the
location of zeros. A typical assumption here is that $Z(B)$ belongs
to certain regions inside the unit disk.

Let $t\in{\mathbb T}$, $\gamma\ge1$. Following \cite{C62, AC74, GPV08} we
introduce regions
\begin{equation}\label{genstol}
R(t,\gamma,K):=\{\lambda\in{\mathbb D}: |t-\lambda|^{\gamma}\le K(1-|\lambda|)\}, \quad K\ge1.
\end{equation}
For $\gamma=1$, $K>1$ this is the standard Stolz angle. When $\gamma>1$ the
region touches the circle ${\mathbb T}$ at the vertex $t$ with the power
degree of tangency. The following result claims that $B'$ belongs to
$H^p$ or $A^p$ as soon as $Z(B)\subset R(t,1,K)$.

{\bf Theorem A}. Let $Z(B)\subset R(t,1,K)$. Then
\begin{enumerate}
\item $B'\in H^p$, $p<\frac12$, and $\frac12$ is sharp;
\smallskip
\item $B'\in A^p$, $p<\frac32$, and $\frac32$ is sharp.
\end{enumerate}

\noindent The first statement is proved in \cite[Theorem
2.3]{GPV07}, for the second one see \cite{GPV07, GP06}. For related
results in the case $Z(B)\subset R(t,\gamma,K)$ with $\gamma>1$, see
\cite[Section 3]{GPV08}.

We study the same problem for more general Stolz--type regions.

A function $\phi$ on the right half-line will be called a {\it model
function}, if it is nonnegative, continuous and increasing, and
\begin{equation}\label{model}
\phi(x)\le Cx, \qquad x\ge0, \quad C=C(\phi)>0.
\end{equation}
We define a {\it Stolz angle associated with a model function
$\phi$} with the vertex at $t\in{\mathbb T}$ as
\begin{equation}\label{expstol}
{\mathcal{S}}_\phi(t,K)={\mathcal{S}}(t,K):=\{\lambda\in{\mathbb D}: \phi(|t-\lambda|)\le K(1-|\lambda|)\},
\quad K>0.
\end{equation}
Since $|t-\lambda|\le2$ for $t,\lambda\in\overline{\mathbb D}$, it is clear that regions
\eqref{genstol} are of the form \eqref{expstol} for an appropriate
$\phi$. Precisely, one can put $\phi(x)=x^\gamma$ for $0\le x\le 2$,
and $\phi(x)=2^{\gamma-1}x$ for $x\ge2$.
Next, given a closed set $E=\overline E\subset{\mathbb T}$ we define a
{\it Stolz region, associated with a model function $\phi$ and the set
$E$}, as
\begin{equation*}
{\mathcal{S}}(E,K):=\{\lambda\in{\mathbb D}: \phi(d(\lambda,E))\le
K(1-|\lambda|)\}=\bigcup\limits_{t\in E} {\mathcal{S}}(t,K).
\end{equation*}

Here is our main result.
\begin{theorem}\label{mainres}
Let $B$ be a Blaschke product such that $Z(B)\subset{\mathcal{S}}(E,K)$. Then
\begin{equation}\label{bla4}
|B'(z)|\le
2(2C+K)^2
\sum_{n=1}^\infty(1-|z_n|)\,\phi^{-2}\left(\frac{d(z,E)}6\right)\,.
\end{equation}
\end{theorem}

For the standard Stolz angle and $E=\{t\}$ we take $\phi(x)=x$, so
$$ |B'(z)|\le \frac{C_3}{|t-z|^2}, $$
and part (1) in Theorem A follows. Similarly, for the region $R(t,\gamma,K)$, $\gamma>1$, \eqref{bla4} implies
$$ |B'(z)|\le \frac{C_4}{|t-z|^{2\gamma}}, $$
and we come to the following result (cf. \cite[Remark 1]{GPV08}).
\begin{corollary}
If $Z(B)\subset R(t,\gamma,K)$, $\gamma>1$, then $B'\in H^p$ for all $p<1/2\gamma$.
\end{corollary}

\medskip

We are particularly interested in the model function
$\phi(x)=\exp\{-x^{-\rho}\}$, $\rho>0$. In this case \eqref{bla4}
says that $B'\in{\mathcal{A}}(E,\rho)$.

Denote $Z(B')=\{z_n'\}$ the zero set of $B'$. Each result of the
form $B'\in X$ provides some information about the critical points
of $B$ (zeros of $B'$), as long as the information about zero sets
of functions from $X$ is available. The most general condition
applied to an arbitrary Blaschke product arises from the fact that
$B'\in A^p$, $p<1$, so (cf. \cite[Theorem 4.7]{HKZ})
$$ \sum_{n=1}^\infty
\frac{1-|z_n'|}{\left(\log\frac1{1-|z_n'|}\right)^{1+\varepsilon}}<\infty,
\qquad \forall\varepsilon>0. $$
 On the other hand there are Blaschke products $B$ such that $\sum
1-|z_n'|=\infty$ (see, e.g., \cite{Pom}).

A Blaschke--type condition for zeros of functions from ${\mathcal{A}}(E,\rho)$
is given in a recent paper \cite{FG}. To present its main result we
define, following P. Ahern and D. Clark \cite[p.113]{AC76}, the type
$\beta(E)$ of a closed subset $E$ of the unit circle as
$$ \beta(E):=\sup\{\beta\in{\mathbb R}:\ |E_x|=O(x^\beta), \ x\to 0\}, $$
where
$$ E_x:=\{t\in{\mathbb T}:\ d(t,E)<x\}, \quad x>0, $$
is an $x$-neighborhood of $E$, $|E_x|$ its normalized Lebesgue
measure. For the equivalent definition and properties of the type
see also \cite{FG, FG10}.

\begin{theorem}\label{mainres1}
Given a closed set $E\subset{\mathbb T}$ and a Blaschke product $B$, assume
that $Z(B)\subset{\mathcal{S}}_\phi(E,K)$, $\phi(x)=\exp\{-x^{-\rho}\}$,
$\rho>0$. Then
\begin{equation*}
\sum_{n=1}^\infty (1-|z_n'|)d^{(\rho-\beta(E)+\varepsilon)_+}(z_n',E)<\infty,
\qquad \forall\varepsilon>0,
\end{equation*}
$\beta(E)$ is the type of $E$, $(a)_+=\max(a,0)$.
\end{theorem}
It is clear that $\beta(E)=1$ for each finite set $E$.
\begin{corollary}
Let $Z(B)\subset{\mathcal{S}}_\phi(t,K)$ with the same $\phi$. Then
\begin{equation*}
\sum_{n=1}^\infty (1-|z_n'|)|t-z_n'|^{(\rho-1+\varepsilon)_+}<\infty, \qquad
\forall\varepsilon>0.
\end{equation*}
\end{corollary}

The authors thank the referee for a number of comments that improved the paper write-up.

\section{Main results}

A model function $\phi$ is nonnegative
and increasing, so for all $x,y,u\ge 0$
\begin{equation}\label{phi}
\phi\left(\frac{x+y+u}3\right)\le \phi(x)+\phi(y)+\phi(u).
\end{equation}

We begin with the following result, which is similar to Vinogradov's
lemma from \cite{V94}.
\begin{lemma}\label{mainlem}
Let $z\in\overline{\mathbb D}$, $t\in{\mathbb T}$ and $\lambda\in {\mathcal{S}}(t,K)$. Then
\begin{equation}\label{mainbound}
\frac1{|1-\overline\lambda z|}\, \phi\left(\frac{\left
|t-z|\lambda|\right|}3\right)\le 2C+K.
\end{equation}\end{lemma}
\begin{proof}
With no loss of generality we assume that $t=1$. Since
\begin{equation*}
\begin{split}
\left|1-z|\lambda|\right|&=\left|1-\overline\lambda z+z(\overline\lambda-|\lambda|)\right|\le
|1-\overline\lambda z|+|\overline\lambda-|\lambda|| \\ &\le |1-\overline\lambda z|+(1-|\lambda|)+|1-\lambda|,
\end{split}
\end{equation*}
then by \eqref{phi}
\begin{equation*}
\phi\left(\frac{\left |1-z|\lambda|\right|}3\right)\le
\phi(\left|1-\overline\lambda z \right|)+\phi(1-|\lambda|)+\phi(|1-\lambda|),
\end{equation*} so
$$ \frac1{|1-\overline\lambda z|}\,
\phi\left(\frac{\left |1-z|\lambda|\right|}3\right)\le A_1+A_2+A_3. $$
By \eqref{model}, for the first two terms we have
$$ A_1=\frac{\phi(\left
|1-\overline\lambda z\right|)}{|1-\overline\lambda z|}\le C, \quad
A_2\le \frac{\phi(1-|\lambda|)}{1-|\lambda|}\le C. $$
As for the third one,
\begin{equation*}
A_3
\le\frac{\phi(|1-\lambda|)}{1-|\lambda|}\le K
\end{equation*}
according to the assumption $\lambda\in {\mathcal{S}}(1,K)$. The proof is complete.
\end{proof}

Our main result gives a bound for the derivative $B'$ in the case
when the zero set $Z=Z(B)\subset{\mathcal{S}}(E,K)$.

\noindent {\it Proof of Theorem \ref{mainres}}.
Denote $Z(t):=Z\bigcap {\mathcal{S}}(t,K)$, $t\in E$. Then there is an
at most countable set $\{t_k\}_{k=1}^\omega$, $\omega\le\infty$, $t_k\in E$,
so that $z_k\in Z(t_k)$, and $Z=\bigcup_k Z(t_k)$. It is clear that
there is a disjoint decomposition
$$ Z=\bigcup_k Z_k, \quad Z_k\not=\emptyset, \qquad Z_k\subset Z(t_k),
\qquad Z_j\bigcap Z_k=\emptyset, \quad j\not=k. $$ Let us label the
set $Z$ in such a way that
$$ Z=\{z_{kj}\}, \quad k=1,2,\ldots,\omega, \quad j=1,2,\ldots\omega_k, \quad
\{z_{kj}\}_{j=1}^{\omega_k}\subset Z_k. $$

We proceed with the expression
$$ B'(z)=\sum_{j=1}^\infty b_n'(z)B_n(z), \qquad
B_n(z)=\frac{B(z)}{b_n(z)}\,, $$ so
\begin{equation*}
|B'(z)|\le \sum_{n=1}^\infty \frac{1-|z_n|^2}{|1-\overline z_n
z|^2}\,|B_n(z)|\le \sum_{n=1}^\infty \frac{1-|z_n|^2}{|1-\overline z_n
z|^2}=\sum_{k=1}^\omega\,\sum_{j=1}^{\omega_k}\frac{1-|z_{kj}|^2}{|1-\overline z_{kj}
z|^2}\,.
\end{equation*}
Note that $|t-z|\le 2\left|t-z|\lambda|\right|$ for all $z,\lambda\in{\mathbb D}$ and
$t\in{\mathbb T}$. Indeed,
$$ \left|t-z|\lambda|\right|\ge 1-|z\lambda|\ge 1-|\lambda| $$
and
$$ |t-z|\le \left|t-z|\lambda|\right|+|z|(1-|\lambda|)\le 2\left|t-z|\lambda|\right|, $$
as claimed. Hence
$$
\phi\left(\frac{|t-z|}6\right)\le
\phi\left(\frac{|t-z|\lambda||}3\right)\,,
$$ and by \eqref{mainbound}
\begin{equation*}
\begin{split}
\phi^2\left(\frac{d(z,E)}6\right)\,|B'(z)| &\le
\sum_{k=1}^\omega\,\sum_{j=1}^{\omega_k}\frac{1-|z_{kj}|^2}{|1-\overline z_{kj}
z|^2}\,\phi^2\left(\frac{|t_k-z|}6\right)
\\
&\le \sum_{k=1}^\omega\,\sum_{j=1}^{\omega_k}\frac{1-|z_{kj}|^2}{|1-\overline
z_{kj} z|^2}\,\phi^2\left(\frac{|t_k-z|z_{kj}||}3\right)
\\ &\le 2(2C+K)^2\sum_{n=1}^\infty (1-|z_n|),
\end{split}\end{equation*}
which is \eqref{bla4}. The proof is complete.

Theorem \ref{mainres1} is a direct consequence of Theorem
\ref{mainres} and \cite[Theorem 3]{FG}.


\chapter
Dinner Speech: Follow-up in the age of surveys}

\author{\raggedright \baselineskip=10p
{\bf Paul J. Groot,$^{1}$}\\
{\small \it
(1) Department of Astrophysics/IMAPP, Radboud University Nijmegen,
  P.O. Box 9010, 6500 GL Nijmegen, The Netherlands
}
}

\AuthorContents{Paul Groot}

\AuthorIndex{Groot}{P.}

     \baselineskip=10pt
     \parindent=10pt

\section{Introduction}

\medskip
Many dinner speeches are not prepared but still given, so this is
probably one of the few that was prepared but never given. The reason
for this is the fact that the conference dinner was held on the roof
of the Kyoto Holiday Inn, where the `beer garden' provided a wonderful
ambience, but where any attempt to words of wisdom would have been immediately
lost on the wind. Committing them to paper is my effort to make them
longer lasting.

\medskip
First of all I would like to thank the organisers of this meeting, and
in particular Daisaku-san for bringing together such a diverse and
broad community on the topic of accretion physics. Accretion is one of
those `universal' processes that is still poorly understood, but at
the same time encountered in a large variety of astronomical
settings. From the compact binaries discussed here to young stellar
objects, supernovae fall back disks and gamma ray bursts and the
large, massive disks found in active galactic nuclei.

\medskip
Looking around I see that it is a very active and young community. If
not always young of age then at least of heart! This is a very good
sign, since many challenges still lie before us, and over the next few
years the possibilities in our field will grow as rapidly as the
number of known objects. The field of accreting compact objects is
reinventing itself. After a flurry of activity in the 1980s and 1990s
there appeared to be a slump in activity in the early 2000s. A change
is now occuring from an `object-based' view to a `population-based'
view. No longer are our conferences dominated by talks on single, pet
objects, but more and more the view is widened to address the question
`What does the big picture look like?'. This is a very good and
healthy progression. The field of Cataclysmic Variables (CVs) in
particular was often stigmatized as splitting itself up into yet
another subfield as soon as two new systems were found. Although
questions related to the physics and structure of the myriad of
subpopulations have not always been answered, more and more often a
helicopter view is taken that tries to encompass the populations of
compact binaries as a whole.

\medskip
How do the different populations of compact
binaries (CVs, symbiotics, novae, AM CVns, LMXBs, AMXPs) fit together?
Can we understand underlying processes such as the physics of
the common-envelope by comparing the sizes and characteristics of
different populations with each other? For that matter: do the
properties of a single population already give answers? What
is the relation of compact binaries to other fields in astronomy?
Examples that stand out in this respect are of course the progenitors
of supernovae Type Ia, but also the populations of gravitational wave
emitters and the possibility of new or surviving exoplanets around
white dwarfs. Such a branching out is essential for the continued
health of the field and the anchoring of compact binaries within the
larger community of astrophysics.

\medskip
During the conference we have already seen very nice examples of this
population view. The almost Sisyphus-like effort of Paula Szkody and
collaborators on uncovering the population of CVs in the Sloan Digital
Sky Survey (SDSS) stands out in this respect. This is a prime example of how
a population based view can answer long-standing questions. The
orbital period distribution of the Sloan CVs as compiled by Boris
G\"ansicke and collaborators clearly shows the existence of the
long-awaited and expected orbital period spike. It also shows the
extreme care one has to take in taking observational selection effects
properly into account: a point already shown quite nicely by the work
of Reta Pretorius and Christian Knigge on the pre-SDSS population of
CVs.

\section{Big Surveys: the future is near}

\medskip
The SDSS is, however, only the beginning. The near future will see
enormously large synoptic and static surveys that will allow us to
significantly increase the known populations of compact binaries. Not
only of CVs, but also of other, even more rare or relatively
less-studied populations such as AM CVns, pre-CVs, sdBs, LMXBs,
detached binary white dwarfs, AXMPs, novae, etc. Even more importantly, the
systematic design and execution of most of these surveys allows us to
construct observationally well-understood samples of compact binaries
where it is possible to properly model the selection effects and
correct for them. This is absolutely essential for a proper comparison
of observational results with those from e.g. population synthesis
modeling. Any of us designing and/or performing these large surveys
should therefore aim to execute them as systematically and
homogeneously as possible and to properly quantify the selection
effects that went into the design.

\medskip
Examples of these new and upcoming large surveys that will be very
relevant for our field are the European Galactic Plane Surveys (EGAPS:
IPHAS, UVEX and VPHAS+), Gaia, ATLAS (a Southern Sloan-like survey
with the VST), the Galactic Bulge Survey and the large
synoptic surveys: the Palomar Transient Survey, Pan-Starrs, SkyMapper,
the VVV survey, and also of course Gaia and in the longer run the
LSST. In the radio domain LOFAR will start operation in 2011,
including the Transient Key Project, and proposals for variability
surveys on the SKA-precursors ASKAP (Australia) and MeerKAT (South
Africa) have been submitted. The work presented during this conference
by Brian Warner and Patrick Warner on the Catalina Real Time Transient
(CRTS) Survey and by Martin Still on the Kepler data shows just a
glimpse of what is ahead in these synoptic surveys.

\section{Necessary preparation}

\medskip
As a community we will have to be prepared for these surveys. Not only
by designing efficient algorithms to mine these large databases for
`our' preferred populations, or by having the appropriate theoretical
models to interpret the results, but most importantly to provide
enough resources to execute the task of the photometric and
spectroscopic follow-up on these systems. Many of the targets provided
by the large surveys will require follow-up that is not a part of the
surveys themselves. To characterize systems in terms of masses,
orbital periods, mass-transfer rates, secular evolution and chemical
composition will require vast photometric and spectroscopic
follow-up. We need to be able not only to bring these resources to the
field, but should also try to organize ourselves in such a way that we
can do this follow-up in an efficient and homogeneous way.

\medskip
Luckily we do, collectively, have the means to do this. We are a world
community. I count participants to this conference from 22 different
countries on 5 different continents. The only continent that is
missing is Australia, but luckily we have Martin Still from outer
space to make up for this. This large and widely-spread community
offers us an enormous opportunity to tackle the challenge before
us. The world will see only 2 or 3 Extremely Large Telescopes, but
there are still many smaller telescopes in the 1-4m class that are
perfectly suited for a large part of the required follow-up. In the
western world the 2-4m class telescopes are under severe funding
pressure, but I am delighted to see that in many other countries new
telescopes in this category are coming into operation. Here we should
tackle the challenge of systematic follow-up by trying to
standardize our follow-up observations as much as possible. We should
be like Japanese vending machines: available on every corner with a
standard output, but they work like a charm! We should try to set-up a
cheap, efficient network of `follow-up vending machines' that deliver
a homogeneous set of data that can easily be intercompared.

\section{Multi-band photometry}

\medskip
My proposal is {\sl not} to build a new network of telescopes. These
we already have to a large extent. The proposal is to standardize our
follow-up methods and instruments. It strikes me again and again how
often we `waste' astronomical information because our instruments are
not designed to take full advantage of what modern technology can
offer. Also during this conference we have seen beautiful photometric follow-up
studies where on a particular object, for instance, 3000 frames have been
taken, of which 2800 in, say, $R$, 200 in $V$ and 2 in $I$. This is
certainly not unique to our field. I see the same happening in the
gamma-ray burst and supernovae follow-up, although the last is
slightly better since the MCSL method for calibrating Supernovae Type
Ia requires multi-band data.

\medskip
Modern technology, however, should allow us to build a simple
multi-band photometer that covers, e.g. the SDSS bands simultaneously
and allows for a 1 second read-out of the CCDs. This is the type of
instrument that we should have as a standard 'photometer' on as many
of the 1-4m class telescopes in the world. Only when we these
instruments will it be possible to come to a homogeneous approach to
the photometric follow-up of the thousands of targets that will be
provided by the next generation of large scale surveys.

\section{Closing}

\medskip
In closing I would like to say again how much I am encouraged by the
diversity of researchers present here. It has been a week of high
spirits, and enthousiasm, and in this, I would like to say that we, as
a community, should also take our Japanese hosts as an example. I have
rarely been in a country where the cultural and language barriers
appeared to be so large, but where the crossing of these barriers is
undertaken with such an enthousiasm and such friendliness.

\bigskip
{\sl Domo arigato} and {\sl itadakimasu}

\end{document}


\section{Introduction}

For the last three decades,
the first-order phase transition remains one of the most important topics
in the cosmology, particle physics, and condensed matter.
In particular, this topic is widely related to many phenomena about the early universe,
e.g., the electroweak phase transition \cite{EWbaryogenesis},
the quark confinement~\cite{QGP},
and inflationary models \cite{inflation}.
Since the phase transition is the phenomena closely related to bosonic degree,
it is an inevitably important element for supersymmetric models.
Recently, this topic is discussed
in relation to charge and/or color breaking minima
and supersymmetry breaking metastable minima,
which is intensively studied after the recent work by Intriligator et al.~\cite{Intriligator:2006dd}.

For first-order phase transitions,
the decay rate of a false vacuum is the key ingredient.
The estimate of the rate has been studied in many works,
which is usually calculated using the Euclidean action of the bounce solution~\cite{Coleman:1977py,Coleman:1977th}.
The accuracy of the estimate has also been improved
with considering the effects of dissipation and finite temperature~\cite{dissipation,Linde:1981zj}.

In this Letter, we find that a novel correction to the effective action
in view of quantum field theory,
which is induced by a field coupling to the scalar field
related to a first-order phase transition.
Hence, this effect is considered to be induced by the environment.
This correction cannot be treated in the method of the effective potential
and acts as dissipative and fluctuation effects on tunneling phenomena.
For a phase transition induced by classical thermal fluctuation,
some effects of environment have been studied so far.
For the quantum tunneling,
the effect of the environment has been studied in viewpoints of quantum mechanics.
We study this effect in a framework of quantum field theory
and find a novel correction.
This effect changes drastically the decay rate of the false vacuum,
since it is the correction to the action,
where the decay rate depends on the action exponentially.
It is also shown that this effect drastically makes the lifetime of the false vacuum
either longer or shorter depending on the interaction which couples the scalar field.

\section{Tunneling without Environment}

\begin{figure}[t]
\scalebox{0.7}{\includegraphics{V.eps}}
\caption{\small Potential describing the tunneling of the $\phi$ field.}
\label{fig: Potential}
\end{figure}

We briefly summarize how we calculate the probability of the tunneling without the environment using the bounce method~\cite{Coleman:1977py} and clarify our notation. With the potential shown in Fig.~\ref{fig: Potential}, the Lagrangian which describes a first-order phase transition of the $\phi$ field is given by
\begin{eqnarray}
{\cal L}_\phi
=
\frac{1}{2}
\left(\partial_\mu \phi \right)
\left(\partial^\mu \phi \right)
-
V(\phi).
\label{Lagrangian without environment}
\end{eqnarray}
The tunneling probability from the false vacuum to the true one, per unit time per unit volume, in the semi-classical approximation is given by the following formula,
\begin{eqnarray}
\Gamma/V \simeq K \exp\left(-{\cal S}_{\rm E}[\phi_{\rm B}]\right),
\end{eqnarray}
where the configuration $\phi_{\rm B}$ is called the bounce solution, and it satisfies the equation of motion of the Euclidean action ${\cal S}_{\rm E}[\phi]$ deduced from the Lagrangian (\ref{Lagrangian without environment}),
\begin{eqnarray}
-\left(\frac{d^2}{dr^2} + \frac{3}{r}\frac{d}{dr}\right)\phi_{\rm B}
+ V^\prime (\phi_{\rm B}) = 0.
\end{eqnarray}
Here, the prime denotes the derivative by the argument. The O(4) symmetry of the solution, $\phi_{\rm B}(x_{\rm E}) = \phi_{\rm B}(r = |x_{\rm E}|)$, is assumed in the above equation, where $x_{\rm E}$ is the coordinates of the Euclidean space-time~\cite{Coleman:1977th}. The solution also satisfies the boundary conditions, $\phi_{\rm B}(\infty) = 0$ and $\dot{\phi}_{\rm B}(0) = 0$, where the dot indicates the derivative by $r$.

On the other hand, the coefficient $K$ has the form of
\begin{eqnarray}
K = \frac{{\cal S}_{\rm E}[\phi_{\rm B}]^2}{4\pi^2}
\sqrt{
\left|
\frac{\det[-\partial_\mu\partial_\mu + V^{\prime\prime}(0)]}
{\det^\prime[-\partial_\mu\partial_\mu + V^{\prime\prime}(\phi_{\rm B})]}
\right|},
\label{Coefficient K}
\end{eqnarray}
where the primed determinant means that the zero modes of the operator, $-\partial_\mu\partial_\mu + V^{\prime\prime}$, which correspond to the translational invariance about the bounce's location, is subtracted. Instead of the subtraction, the bounce action ${\cal S}_{\rm E}[\phi_{\rm B}]$ appears in the right-hand side of Eq.~(\ref{Coefficient K}).

\section{Tunneling with Environment}

We discuss how the effect of the environment, namely, interactions between $\phi$ and other fields alter the tunneling probability. For quantitative discussion, we consider two specific cases described by the following Lagrangians,
\begin{eqnarray}
{\cal L}_{\rm tot}
&=&
{\cal L}_\phi
+
{\cal L}^{(I)}_{\rm env}
+
{\cal L}_c,
\label{Lagrangian with environment}
\\
{\cal L}_{\rm env}^{(S)}
&=&
-S^*(\Box + m_S^2 + y_S M \phi)S,
\label{Lagrangian S}
\\
{\cal L}_{\rm env}^{(F)}
&=&
+\bar{F}\left(i\Slash{\partial} - m_F - y_F \phi \right) F,
\label{Lagrangian F}
\end{eqnarray}
where $(I)$ is $(S)$ or $(F)$, and the field $S$ ($F$) is a complex scalar boson (a Dirac fermion) with $m_S$ ($m_F$) being its mass, and $M$ is the curvature of the potential minimum at $\phi = 0$, namely, the mass of $\phi$ in the false vacuum. The strength of the interaction between $\phi$ and $S$ ($F$) is characterized by the coupling constant $y_S$ ($y_F$). The counter terms ${\cal L}_c$ are introduced in the Lagrangians, which regularize the environmental effect of the interactions. The on-shell renormalization condition is adopted for ${\cal L}_c$.

The probability of the tunneling with the environment is obtained from the effective Euclidean action. This action is derived by integrating the environment field ($S$ or $F$) out from the total Euclidean action deduced from the Lagrangian (\ref{Lagrangian with environment}). This action has the form of
\begin{align}
{\cal S}_{\rm E}^{(I)}[\phi]
=&
\int d^4x_{\rm E} \,
\left[
\frac{1}{2} \dot{\phi}^2
+
\frac{1}{2} \left(\nabla \phi\right)^2
+
V(\phi)
\right]
\label{Effective Euclidean action}
\\
&
+\frac{1}{2}
\int d^4x_{\rm E} \int d^4x^\prime_{\rm E} \,
\phi(x_{\rm E})
f_0^{(I)}(x_{\rm E} - x^\prime_{\rm E})
\phi(x^\prime_{\rm E}),
\nonumber
\end{align}
where we implicitly assume that coupling constants ($y_S$ and $y_F$) are small enough and hence neglect higher order terms of $\phi$ than $\phi^2$. The effect of the environment is summarized in the function $f_0^{(I)}(x_{\rm E})$, and, for later use, it is convenient to express it in the Fourier transform,
\begin{eqnarray}
f_0^{(I)}(x_{\rm E})
=
\int \frac{d^4p}{(2\pi)^4} \tilde{f}^{(I)}(p) e^{-i p x_{\rm E}}
+
\tilde{f}_0^{(I)}(0) \, \delta^{(4)}(x_{\rm E}),
\nonumber
\end{eqnarray}
where $\tilde{f}^{(I)}(p) = \tilde{f}_0^{(I)}(p) - \tilde{f}_0^{(I)}(0)$. The term proportional to a delta-function gives a finite correction to the potential $V(\phi)$, which is renormalized by the redefinition of the mass term in the potential. On the other hand, the term involving $\tilde{f}^{(I)}(p)$ cannot be a correction to $V(\phi)$, because this term vanishes when $\phi$ is independent of $x_{\rm E}$. In fact, when we expand $\phi(x_{\rm E}) \int d^4x^\prime_{\rm E} f^{(I)}(x_{\rm E} - x^\prime_{\rm E}) \phi(x^\prime_{\rm E})$ in terms of local operators including two $\phi$'s, these operators have to involve the derivative of $\phi$ such as $\phi \Box \phi$. This effect of the environment, therefore, is never taken into account in the method of the effective potential. In this Letter, we focus mainly on how this effect alters the tunneling probability.

Explicit expression of the function $\tilde{f}^{(I)}(p)$ obtained from the interactions (\ref{Lagrangian S}) and (\ref{Lagrangian F}) are given by
\begin{eqnarray}
\tilde{f}^{(I)}(p)
&=&
\int^\Lambda_{2m_I} d\omega
\frac{p^2 \, r^{(I)}(\omega)}{\omega^2(\omega^2 + p^2)}
-
p^2 \Pi^{(I)\prime}(M^2),
\label{Function f}
\\
r^{(S)}(\omega)
&=&
\frac{\alpha_S}{2\pi} M^2 \left(\omega^2 - 4m_S^2\right)^{1/2},
\label{r of S}
\\
r^{(F)}(\omega)
&=&
\frac{\alpha_F}{\pi} \left(\omega^2 - 4m_F^2\right)^{3/2},
\label{r of F}
\end{eqnarray}
where $\alpha_I = y_I^2/(4\pi)$ and $\Lambda$ is the cutoff parameter, which is eventually taken to be infinity. The function $\Pi^{(I)}(M^2)$ comes from the counter terms ${\cal L}_c$, which is defined by
\begin{eqnarray}
\Pi^{(I)}(M^2)
=
{\cal P} \int^\Lambda_{2m_I} d\omega \, \frac{r^{(I)}(\omega)}{\omega^2 - M^2},
\end{eqnarray}
where ${\cal P}$ denotes a principal valued integral.

The equation of motion for the bounce solution with the environment is obtained from the effective Euclidean action (\ref{Effective Euclidean action}). Imposing the O(4) symmetry on the solution, the equation turns out to be
\begin{eqnarray}
&&
-\left(\frac{d^2}{dr^2} + \frac{3}{r}\frac{d}{dr}\right)\phi^{(I)}_{\rm B}
+
V^\prime (\phi^{(I)}_{\rm B})
\label{EOM with environment}
\\
&&
+ \int dr' \int dp \,
\frac{p \, r^{\prime 2} J_1(pr) J_1(pr')}{r}
\tilde{f}^{(I)}(p) \phi^{(I)}_{\rm B}(r') = 0,
\nonumber
\end{eqnarray}
with the boundary conditions, $\phi^{(I)}_{\rm B}(\infty) = 0$ and $\dot{\phi}^{(I)}_{\rm B}(0) = 0$. Here, $J_1(x)$ is the Bessel function of the first kind. Once the bounce solution $\phi^{(I)}_{\rm B}$ is computed by solving the above integro-differential equation, the probability of the tunneling with the environment is obtained as
\begin{eqnarray}
\Gamma(\alpha_I)/V
\simeq
K^{(I)} \exp \left(-{\cal S}^{(I)}_{\rm E}[\phi^{(I)}_{\rm B}]\right).
\label{Probability with environment}
\end{eqnarray}
Since we are interested in cases where interactions between $\phi$ and environment fields are weak enough, the effect of the environment in the coefficient $K^{(I)}$ is weak, at most, ${\cal O}$(10)\% level. We therefore neglect such a correction in the following discussion. On the other hand, the effect in ${\cal S}^{(I)}_{\rm E}[\phi^{(I)}_{\rm B}]$ can be very significant because of its exponential sensitivity as seen in the probability (\ref{Probability with environment}).

\begin{figure}[t]
\hspace{-8mm}
\scalebox{0.7}{\includegraphics{S.eps}}
\scalebox{0.7}{\includegraphics{F.eps}}
\qquad\qquad
\caption{\small Ratio of the tunneling probability between with and without the environment.}
\label{fig: Tunneling ratio}
\end{figure}

In Fig.~\ref{fig: Tunneling ratio}, the ratio of the tunneling probability between with and without the environment is shown as a function of $\alpha_I$ with several choices of $m_I$. A cubic potential,
\begin{eqnarray}
V(\phi)
=
\frac{1}{2} M^2 \phi^2 \left(1 - \frac{\phi}{\phi_0}\right),
\label{Cubic potential}
\end{eqnarray}
is used to calculate the ratio. The metastable well and the potential barrier are assumed to be well approximated by the above potential. The parameter $\phi_0/M$ is fixed so that the bounce action without the environment takes the value of 400, which is required for the stability of the false vacuum. As can be seen in the figure, the effect of the environment becomes significant when $m_I$ is small, especially, in the case $(S)$.

\section{Effect of Environment}

We are now in position to discuss the effect of the environment shown in Fig.~\ref{fig: Tunneling ratio} in detail. For the intuitive understanding of the effect, it is convenient to rewrite the integro-differential equation (\ref{EOM with environment}) to a simple integral equation with the help of the following transformations,
\begin{eqnarray}
\phi^{(I)}_{\rm B}(r)
&=&
\int^\infty_0 dp \, \frac{p}{r} J_1(pr) \, \tilde{\phi}^{(I)}_{\rm B}(p),
\\
\tilde{\phi}^{(I)}_{\rm B}(p)
&=&
\int^\infty_0 dr \, r^2 J_1(pr) \, \phi^{(I)}_{\rm B}(r).
\end{eqnarray}
It is easy to confirm that the above transformations satisfy the boundary conditions for $\phi^{(I)}_{\rm B}$. It then turns out that Eq.~(\ref{EOM with environment}) is translated into
\begin{eqnarray}
&&
\frac{\pi \phi_0}{3M^2}
\left[p^2 + M^2 + \tilde{f}^{(I)}(p)\right]
\tilde{\phi}^{(I)}_{\rm B}(p)
\label{Integral equation}
\\
&&
=
\int^\infty_0 dp' \int^{p + p'}_{|p - p'|} dp'' \,
\tilde{\phi}^{(I)}_{\rm B}(p') \, \tilde{\phi}^{(I)}_{\rm B}(p'')
\frac{\Delta(p, p', p'')}{p},
\nonumber
\end{eqnarray}
where the function $\Delta(p, p', p'')$ gives the area of the triangle with the sides $p$, $p'$, and $p''$. The cubic potential (\ref{Cubic potential}) is assumed in the above equation, while it is also possible to apply the translation to more general cases. Using the function $\tilde{\phi}^{(I)}_{\rm B}(p)$, the bounce action is given by
\begin{eqnarray}
{\cal S}_{\rm E}^{(I)}[\phi^{(I)}_{\rm B}]
&=&
\frac{\pi M^2}{\phi_0} \int dp \, dp' \, dp''
\label{Bounce action}
\\
&\times&
\tilde{\phi}^{(I)}_{\rm B}(p) \,
\tilde{\phi}^{(I)}_{\rm B}(p') \,
\tilde{\phi}^{(I)}_{\rm B}(p'') \,
\Delta(p, p', p'').
\nonumber
\end{eqnarray}

We first notice that $\tilde{\phi}^{(I)}_{\rm B}(p)$ is rapidly decreased when $p \gg M$, because Eq.~(\ref{Integral equation}) gives the relation $(\phi_0/M^2) \, p^2 \, \tilde{\phi}^{(I)}_{\rm B} \sim p^3 \, (\tilde{\phi}^{(I)}_{\rm B})^2$, which is obtained from a simple dimensional analysis. The function $\tilde{\phi}^{(I)}_{\rm B}(p)$ has, therefore, a non-negligible value in the region $p \lesssim M$. On the other hand, Eq.~(\ref{Integral equation}) also indicates, when the effect of the environment, $\tilde{f}^{(I)}(p)$, is larger, $\tilde{\phi}^{(I)}_{\rm B}(p)$ is larger in order to keep a balance between the right- and left-hand sides of Eq.~(\ref{Integral equation}). It leads to larger action (\ref{Bounce action}), namely, smaller tunneling probability.

\begin{figure}[t]
\hspace{-8mm}
\scalebox{0.7}{\includegraphics{gS.eps}}
\scalebox{0.7}{\includegraphics{gF.eps}}
\qquad\qquad
\caption{\small Effect of the environment, $\tilde{f}^{(I)}(p)$,
in the unit of $\alpha_I M^2$.}
\label{fig: Environment effect}
\end{figure}

In Fig.~\ref{fig: Environment effect}, the effect of the environment, $\tilde{f}^{(I)}(p)$, is shown as a function of $p$ with several choices of $m_I$. It can be seen that the effect is larger when $m_I$ is smaller, while no enhancement can be found in the case of $m_I = M$. This fact can be intuitively understood, because the $\phi$ field cannot decay into environment fields when $M \leq 2m_I$ and the dissipative effect on the tunneling disappears. When $M > 2m_I$, this effect is more significant for smaller $m_I$ as expected. The difference between the cases $(S)$ and $(F)$ on this effect comes from the fact that the $\phi$ field decays to $SS^*$ through the S-wave process, while it decays to $F\bar{F}$ through the P-wave one. The function $\tilde{f}^{(I)}(p)$, therefore, grows quickly at $p \lesssim M$ in the case $(S)$, which leads to a strong dissipative effect on the tunneling as shown in Fig.~\ref{fig: Tunneling ratio}. On the other hand, in the case $(F)$, $\tilde{f}^{(I)}(p)$ grows slowly, which result in a small dissipative effect as also shown in Fig.~\ref{fig: Tunneling ratio}.

There is another important effect of the environment, which can be especially seen in the case $(F)$. The function $\tilde{f}^{(I)}(p)$ can be negative, leading to the suppression of the bounce action, namely, the enhancement of the tunneling probability. Since the effect does not vanish even in the case of $M = m_I$, it cannot be a dissipative effect. Instead, it should be regarded as a effect of fluctuation due to the interaction between $\phi$ and an environment field. This fact can be understood by considering the case $(F)$ at $m_F = 0$ limit, where $\tilde{f}^{(I)}(p)$ is analytically given by $\tilde{f}^{(F)}(p) = (\alpha_F/2\pi) [p^2 - p^2 \ln(p^2/M^2)]$. The logarithmic term is the origin of the suppression. The combination of $p^2 + \tilde{f}^{(F)}(p)$ in Eq.~(\ref{Integral equation}) is rewritten as $\sim M^2(p^2/M^2)^{1 - \alpha_F/(2\pi)}$, which is nothing but the anomalous dimension from the Yukawa interaction $y_F\phi\bar{F}F$. Since the existence of this effect and the weak dissipative effect in the case $(F)$, the tunneling probability is enhanced compared to that without the environment as shown in Fig.~\ref{fig: Tunneling ratio}.

\section{Summary}

We have investigated the effect of the environment on the tunneling probability which cannot be taken into account in the method of the effective potential. We have formulated how the effect is involved in the calculation of the tunneling probability within the framework of the bounce method and shown that the effect indeed can be very significant. We have also clarified what kind of interaction enhance or suppress the tunneling probability. It is, therefore, important to consider this effect for the estimation of the probability, especially in models built on a false vacuum. It may also be interesting to consider the scenario on a false vacuum which is stabilized using the environment discussed in this Letter.

\section*{Acknowledgements}

The work of S. M. was supported in part by the Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) (Nos.\ 21740174 and 22244021).
The work of K. I. N. was supported in part
by Grant-in-Aid for Nagoya University Global COE Program from the MEXT
and the Grant-in-Aid for JSPS Fellows.


\section{Impact of ELSA}

First, what does it take to establish and maintain a leadership position in the
field of space astrometry? I believe the following ingredients are important.
The basic condition is of course to have a wide scientific community interested
in the data that can be obtained from dedicated astrometric surveys. The
interest should not just be in the classical applications of astrometry (i.e.,
catalogues of star positions) but in pushing the limits of achievable
astrometric, photometric and radial velocity accuracies in order to expand our
horizons and learn more about our universe. This strong scientific interest
should be complemented by an expert community of astronomers capable of
designing and delivering missions such as Gaia. Because of the nature of space
astrometric missions, in terms of hardware (spacecraft and payload design) and
software (data processing), these experts should welcome contributions from the
fields of engineering, mathematics, and nowadays, computer science. The third
important ingredient is the existence of competent industrial partners capable
of delivering the very complex satellites required for space astrometry.
Finally, it is crucial that there are good contacts between these three
communities. The scientific users of the space astrometric survey have to
understand the limitations of what can be achieved, while the space astrometry
experts should not work in isolation but be motivated by the science that will
be done with the data they produce. The latter will also very much facilitate
the transmission of scientific requirements to ESA and its industrial partners.

The existence of a scientific community in Europe interested in space astrometry
is aptly demonstrated by the numerous contributions to this volume, covering a
broad range of science topics, that show the eager anticipation of the Gaia
data. The long astrometry tradition in Europe combined with the efforts to
realize the Hipparcos mission have led to a world leading community of space
astrometry experts. The Hipparcos experience was also incorporated within ESA
and its industrial partners. How has ELSA contributed to maintaining and
extending the four above mentioned components of Europe's leading role in space
astrometry?

\paragraph{Maximizing the scientific return from Gaia data}

If the Gaia catalogue with its 1 billion sources --- including high accuracy
astrometry, photometry, and radial velocities for over 100 million sources ---
would land on one's desk today it is not obvious how to most efficiently exploit
this huge data set. Creative thinking and advance preparation is required in
order to be ready to apply novel and efficient data analysis methods to the Gaia
results. Moreover, the unprecedented levels of accuracy make it essential to
think about possible complications in the interpretation of the data. Eight of
the ELSA research projects were directly relevant to these areas. The work
described by V\'aradi makes a strong case for publishing per-CCD photometry for
the (short-period) variable stars observed by Gaia, and illustrates how the
survey requirements evolve in response to preparatory research. Re Fiorentin
provides an excellent example of how one can do current science as well as
prepare for the Gaia data by analyzing the combination of existing astrometric
(GSC-II), photometric (SDSS), and spectroscopic surveys (RAVE). The works of
Santoro and Saguner illustrate how Gaia has stimulated new theoretical research
into stellar structure and atmospheres as well as preparatory ground-based
observational campaigns, which can already be exploited scientifically. Further
preparations for the scientific analysis of Gaia data are described in the
contributions by Czekai and Belcheva which focus on the non-trivial question of
how to use the Gaia data to study our galaxy or the Magellanic clouds as a
whole. These works also strengthened the Gaia universe model (see contribution
by Luri) leading to more reliable simulations of the sky for mission and
scientific preparations. The investigations by Pasquato into the effects of star
spots on the astrometry from Gaia highlight the importance of re-examining the
established methods of interpreting observational data when the accuracies are
pushed to the limit.  Finally the research into inversion methods by Oskiewicz
has stimulated much discussion about the way in which the results from Gaia
should be transmitted (more on this topic below).

\paragraph{The next generation of space astrometry experts}

One of the central aims of ELSA was of course to transfer the expertise from the
generation involved in the Hipparcos project to a new generation who will be
involved in the preparation and running of the Gaia mission. The network has
been very successful in this respect, not only by directly training a new
generation of space astrometry experts but also by serving as a focal point for
research in support of the Gaia mission preparations. This is most clearly the
case in the areas of the astrometric solution and the radiation damage
mitigation. The AGISLab package developed by Holl is an essential tool in the
development of enhancements to the astrometric iterative solution and will be
used to predict the complicated behaviour of the astrometric errors and the
correlations between them. The conjugate gradient solver for the iterative
astrometric solution was first implemented and tried in AGISLab before it was
transferred to AGIS. This led to a large improvement in the speed of the
astrometric solution as described by Bombrun which in turn facilitates research
into the properties of the solution itself. The efforts by Risquez resulted in
an accurate dynamical model of the Gaia spacecraft attitude which can be used to
furnish realistically simulated attitude data for mission preparations. During
the mission such a tool can aid the understanding of effects seen in the
attitude reconstructed from the observations (see also the contribution by Van
Leeuwen). Finally, in the case of Gaia there is no practical possibility to have
two independent consortia carry out the data processing, yet it remains
important to verify of the reliability of the Gaia catalogue through independent
checks. This task is made possible by the work described by Abbas which will
lead to an independent astrometric data reduction for a subset of the Gaia
sources.

An important problem for the Gaia mission that was identified already at the
time the mission was conceived is the issue of radiation damage to Gaia's CCDs
caused by Solar wind protons. This leads to increased charge transfer
inefficiency in Gaia's detectors and will affect all the data collected by Gaia
by causing systematic errors that have to be accounted for. Again, ELSA served
as a focus for efforts to mitigate the effects of radiation damage. The work by
Prod'homme will result in the most sophisticated model to date of the effects of
radiation damage on the charge collection and read-out process in CCDs. The
contribution by Weiler describes how this work has stimulated the development of
fast analytical models to describe radiation damage effects. Such models will
play a central role in the data processing for Gaia.

\paragraph{Strengthening the partnership with industry}

The demands of the Gaia mission are pushing the limits of what can be achieved
with conventional CCD technology and input from industrial partners has proved
essential in the studies of radiation damage mitigation. The supplier of Gaia's
CCDs, e2v, supported with their expertise the development of microscopic models
of the electron distribution within a Gaia CCD pixel (work by G.\ Seabroke at
Open University) which led to much improved modelling at the pixel level (see
contribution by Prod'homme). The modelling work and the development of the
radiation damage mitigation scheme for Gaia also relied heavily on the
laboratory experiments that EADS-Astrium carried out with irradiated Gaia CCDs
(see contribution by Pasquier). In turn the industrial partners were pushed by
the constant questioning from the scientist to enhance their own understanding
of radiation damaged CCDs. The new knowledge on CCDs was absorbed by ESA as the
facilitating partner and this will certainly be of great benefit to future
European space missions.

\paragraph{Networking and bringing in new expertise}

The opening and closing conferences of the ELSA network were aimed at bringing
together the scientist that will use the Gaia data and the community that will
produce the data. Where the opening conference served to provide the ELSA
fellows with the scientific motivation for their work, the conference of which
this volume forms the proceedings served to show the European astronomical
community how the Gaia data will be produced. I was pleased to hear an
astronomer not involved in Gaia state; ``I'm convinced I came to the right
conference!''. These important contacts between `user' and `producer' were also
fostered within the ELSA network, for example through the use of AGISLab in the
study of Gaia's capabilities with respect to variable stars.

Lastly, the ELSA network also successfully brought much needed outside expertise
into the astronomical community. Examples are the engineering expertise of
Prod'homme and the industrial partners, the mathematical expertise of Bombrun
and the knowledge of high performance computing that is brought in through the
work described by Fries in this volume.

\section{Future proofing the Gaia data}

Gaia will provide an unprecedented stereoscopic map of our Milky Way and the
nearby universe. The catalogue will contain over 1 billion stars, $\sim300\,000$
solar system objects, millions of galaxies, $\sim500\,000$ quasars and thousands
of exoplanets. For all these objects accurate astrometry, photometry, and (for a
subset) spectroscopy will be available as `basic' data. In addition the
classification, variability characterization, and astrophysical parameters of
each object will be provided.  When this catalogue is `finished' around 2020 and
combined with other large sky surveys it will become {\em the} astronomical data
resource for decades thereafter, representing a tremendous discovery potential.

However, I strongly believe that the true potential of the Gaia data can only be
unlocked if we take an ambitious and innovative approach to data publication and
access, including the provision of advanced data analysis tools. This means that
we should not plan the catalogue publication based on what we can image is
possible with current resources but rather base ourselves on what is possible
after 2020. I therefore advocate the following guidelines:

\paragraph{Publish early and publish often} The experience from the current
large sky surveys, notably SDSS and RAVE, has shown that early and regular
releases of the data are a very successful approach to survey publication. As
pointed out in the contribution by Juric, the astronomical community very much
appreciates and makes heavy use of the early releases (even if they are not as
polished as one would like); the survey producers themselves benefit from
immediate user feedback allowing them to correct important errors early on; it
enables the rapid reaction by the community to new discoveries, and will
facilitate synergies between Gaia and other projects. In this context it is
worth keeping mind the enormous synergy possible between Gaia, LSST and
Pan-STARRS, where especially in the case of LSST there is a smooth connection to
the Gaia survey in terms of the astrometric and photometric accuracies achieved.
The goal should thus be to set an ambitious publication schedule with a first
Gaia data release foreseen as soon as the sky has been surveyed once and an
all-sky catalogue of positions, and broad-band photometry can be released ($G$,
$G_\mathrm{BP}$, $G_\mathrm{RP}$). Such a trivial sounding data release (no
parallaxes, motions, or detailed astrophysical characterization yet) will in
fact constitute the highest spatial resolution all-sky map ever produced and I
am convinced that it will lead to exciting discoveries. Such a release would in
addition greatly enhance target selection for surveys that specifically aim at
complementing Gaia with, for example, high resolution spectroscopy for chemical
abundance determinations. This release should be followed with regular releases
containing ever more, more complex, and better data.

\paragraph{Keep raw data, calibration data, and processing software available}
The contribution by Van Leeuwen shows how better insights into the attitude
modelling for Hipparcos combined with present-day computing power enabled a
higher quality re-processing of the entire Hipparcos data set. This resulting
new version of the Hipparcos catalogue features very much reduced error
correlations and improved astrometric accuracies (by up to a factor of 4) for
the bright stars. This is the best illustration of the fact that the raw Gaia
data, all the calibration data, and the processing software should be stored
such that they are permanently accessible and readable, just as the catalogue
itself will be. It is in fact a prerequisite for the next guideline.

\paragraph{Facilitate (re-) processing of the (raw) data} Already in the case of
Hipparcos there are numerous examples of the re-processing of the data, notably
to improve the astrometry of binaries and very red giant stars. Other examples
include the re-processing of intermediate data for groups of stars in order to
derive a common radial velocity or parallax, the re-processing of data for
objects that are discovered or confirmed to be binaries following a data
release, or the re-determination of astrophysical parameters for stars following
improvements in atmosphere modelling. In principle also for Gaia the
re-processing of {\em all} the raw data might be warranted at some point in the
future. In addition to the re-processing of the data the Gaia archive should
also facilitate very complex operations on large chunks of the catalogue (say an
all-sky search for stellar streams). Both these aims may be best served by
implementing the idea of `bringing the processing to the data' by offering users
a virtual machine at the data centre hosting the Gaia archive. On this machine
one could code whatever analysis or processing algorithm is called for and run
it in a way specified by the user.

\paragraph{Make the archive `live'} A concept closely related to the previous
item is that of making the Gaia data archive a `living entity'. By this I mean
that it should be possible to incorporate new information into the catalogue.
Examples are complementary ground-based spectroscopy, updated classifications or
parameterizations of stars based on independent information, better distance
estimates for faint stars, etc. In addition the Gaia archive should seamlessly
integrate with other large sky surveys including ones not foreseen at the time
of the Gaia data publication. As an example, it should be possible to query the
catalogue for sources brighter and fainter than the $G=20$ survey limit of Gaia,
where behind the scenes the work is done to combine Gaia and, for example, LSST
data.

\paragraph{Don't publish a catalogue} This guideline seems in complete
contradiction to the discussion above but is meant to capture in a few words the
ideas presented in the contribution by Hogg \& Lang in this volume. I think
serious efforts should be invested in their proposal to consider a `catalogue'
as only our best current model that explains all the raw data and to publish the
data such that users can test their hypotheses (almost) against the raw image
pixels. I do not expect that this approach will be possible anytime soon but it
should not stop us from being ambitious and trying.

\section{What about ELSA-II?}

Is there a need for a follow-up to the ELSA research activities, i.e.\
investigations at the interface of astronomy, engineering, computer science and
mathematics? I think the answer is yes and that there are in fact at least two
different topics an `ELSA-II' could focus on. The Gaia data will surely
stimulate demands for future space astrometry missions that can open up new
wavelength domains to precision astrometry, push the astrometric accuracy
limits, or extend existing accuracies to much fainter stars. All these would
require technological innovations to break the current scaling laws that would
prevent us from simply building a Gaia$++$. A second and possibly more important
avenue would be to pursue the ideas outlined in the previous section. Questions
to investigate would be: how does one take care of the data curation for Gaia,
including software and calibration data, such that the stored data are easily
retrievable in the future? How do we best go about bringing the processing to
the data or combining Gaia data with other archives in a transparent manner? How
does one implement the testing of hypotheses against the raw data? Is this at
all practically possible? Note that Pfenniger argues in his contribution that
modelling the Galaxy may be best achieved through a very large particle number
$N$-body simulation, possibly containing a particle for every star in the
Galaxy. How does one decide for such models which is best? Do we convert these
models to predictions of the Gaia pixels?

Although these research topics are oriented toward the longer term future they
will generate many spin-offs that can be applied immediately to existing
astronomical instruments, astronomical data, or in industry, thus maintaining
and extending European Leadership in Space Astrometry.

\begin{acknowledgements}
  I thank William O'Mullane and David Hogg for discussions that shaped the ideas
  presented in this contribution. Special thanks go to Stefan Jordan for
  providing the nice portraits of the ELSA fellows which livened up what would
  otherwise have been a rather dull set of closing remarks. Finally, leadership
  requires funding and I thank again the European Community for funding the ELSA
  Marie Curie research training network through contract MRTN-CT-2006-033481.
\end{acknowledgements}

\end{document}