\section{Introduction} \label{sec:intro}
{For the past} decade, exoplanetary systems are being discovered at a spectacular rate \citep[e.g.][]{mayor04,borucki11}. Driven by these discoveries, there is an increasing interest in the global properties of planetary systems, from the epoch of their formation in protoplanetary discs ({PPDs}) to their long-term stability. While there is a natural focus on {internal processes} that govern the {evolution} of such systems {\citep[e.g.][]{lee03,dullemond05,gorti09,morbidelli09,brasser09,blum10,williams11}}, it is clear that {not all planetary systems form in isolation and} environmental effects should be considered as well. Particularly, theoretical studies show that {external} photoevaporation {\citep{scally01,adams04, adams06, fatuzzo08}} and dynamical interactions {\citep{bonnell01,pfalzner05,olczak06,olczak10,spurzem09,lestrade11,dukes12,parker12,bate12}} can lead to the truncation of {PPDs} and planetary systems.

{While the external photoevaporation of PPDs has been studied observationally \citep{odell93,robberto08,rigliaco09,mann10}, there is no conclusive evidence of dynamical effects \citep{eisner06,olczak08,reche09}.} In part, this {is likely} due to the relatively short lifetimes of {PPDs} \citep[up to $\sim8$~Myr, e.g.][]{haisch01,hernandez08,ercolano11,smith12} compared to the time it takes stellar encounters to have an observable effect on the disc {population \citep[$\sim0.1$--$1$~Gyr, e.g.][]{adams10}. Previous observational studies on the Orion Nebula Cluster {(ONC)} \citep[e.g.][]{eisner06,olczak08} did aim to find traces of dynamical interactions in the population of PPDs, but lacked a sufficient number of sources and/or suffered from uncertainties on the disc mass measurements.}

We address the problem statistically by considering the sizes of {PPDs} as a function of {their} ambient stellar density, using samples of {PPDs} and young stellar objects {(YSOs)} from the {latest infrared surveys}. If stellar encounters truncate {PPDs} by tidally stripping the outskirts of the discs {\citep[e.g.][]{clarke93,heller95,hall96}}, {this should be observable} above some characteristic {ambient stellar} density {\citep[$\sim10^3~{\rm pc}^{-3}$][]{adams10}}, because the {encounter rate} increases with density {\citep[e.g.][Eq.~7-61]{binney87}}. In this Letter, we find model-independent evidence of {a change in the {PPD} size distribution for ambient stellar surface densities $\Sigma\mathrel{\copy\simgreatbox} 10^{3.5}~{\rm pc}^{-2}$} at the $>97\%$ confidence level.

\section{{Protoplanetary} discs and their environment} \label{sec:method}

\subsection{Data selection} \label{sec:data}
To verify whether a relationship exists between the sizes of {PPD}s and their {ambient stellar surface density, $\Sigma$, we} combine existing catalogues of {PPD}s and young stellar objects {in star-forming regions (SFRs) of} the solar neighbourhood.

The data for the {PPD}s is taken from {circumstellardisks.org} {(Karl Stapelfeldt, NASA/JPL)}. This catalogue gathers resolved {PPD}s that have been confirmed and described in the literature. {If a {PPD} is resolved} in different wavelengths {\citep[probing different parts of the disc, see e.g.][]{lada06}}, the catalogue lists the largest {measured} diameter, implying that the disc radii used in this work are lower limits. {About $75\%$ of the {PPD} radii in our final sample {(see below)} are measured at wavelengths around $1~\mu{\rm m}$, with only $\sim25\%$ of the {PPD}s (all in low-density regions) being observed at mm wavelengths {(see Appendix~\ref{sec:app0})}. The {catalogue contains} an estimate for how well-resolved each disc is by listing the number of diffraction-limited beams that fit within its diameter. We only consider the discs for which this value is greater than unity. {This provides us with $133$ {PPD} sources from which we exclude those {whose host star is} classified as a main sequence star, weak-line T-Tauri star, {or} Class~0 YSO. Sources at distances $>500$~pc {(which covers all our YSOs) are also excluded}. These criteria reduce the sample to a total of $101$ sources.

To estimate the local ambient surface density of each {PPD} source, we use publicly available near-infrared data of nearby SFRs (see Table~\ref{tab:regions} in Appendix~\ref{sec:app0}). The ambient surface density of stars around each {PPD} is estimated as in \citet{casertano85} {-- see Appendix~\ref{sec:app0} for a detailed explanation}. The thus-obtained angular {ambient} surface densities are converted to {physical ambient stellar} surface densities, $\Sigma$, using the distances listed in Table~\ref{tab:regions}. In cases where the listed {PPD} distance differs from the distance to the region that it is a member of, we adjust the {PPD} distance and size. Since a low surface density can be due to incomplete YSO coverage, any discs with $\Sigma<0.1~{\rm pc}^{-2}$ are omitted from our analysis. {Moreover, the minimum {PPD} radius that can be resolved increases with distance and hence introduces a distance-dependent selection bias. To avoid this, we exclude any {PPD}s smaller than the smallest radius ($\sim50$~AU) that is resolved in the most distant region of our sample, which is the Orion Nebula Cluster (ONC). The final sample thus contains $67$ sources (see Table~\ref{tab:sources} in Appendix~\ref{sec:app0})}.}{ Completeness does not affect the densities because the surveys of Table~\ref{tab:regions} are complete down to the hydrogen-burning limit.}

\subsection{A simple theoretical estimate for the truncation of discs} \label{sec:theory}
To interpret the data, we include a rough theoretical estimate for the expected truncation radii of {PPD}s as a function {of $\Sigma$. This is obtained by combining the truncation induced by each particular encounter with the stellar encounter rate}. The derivation is presented in detail in {Appendix~\ref{sec:app1}.

We use the numerical simulations of disc perturbations in clustered environments by \citet{olczak06} to obtain the disc radius as a function of the encounter parameters. {We convert their expressions for the disc mass loss to a radial truncation assuming that it occurs by stripping the outer disc layers to the Lagrangian point between both stars, and adopting a PPD surface density profile $\Sigma_{\rm d}\propto r^{-1}$ \citep{olczak06}}. Under these assumptions, {we write for the upper limit to the disc radius}
\begin{equation}
\label{eq:rd}
r_{\rm d}(r_{\rm p},m_1,m_2)=\frac{r_{\rm p}}{\sqrt{m_2/m_1}+1} ,
\end{equation}
where $r_{\rm p}$ is the pericentre radius at which the perturber passes, $m_1$ the mass of the perturbed system, and $m_2$ the mass of the {perturber.} This approximation is validated in Appendix~\ref{sec:app1}.

The encounter radius $r_{\rm p}$ is obtained from the impact parameter $b$, encounter velocity $v$, and masses $m_1$ and $m_2$ by accounting for gravitational focussing {(see Appendix~\ref{sec:app1})}. The masses are assumed to follow a \citet{salpeter55} type initial mass function in the range 0.1--100~\mbox{M$_\odot$}, and the velocity distribution is taken to be Maxwellian {with a} velocity dispersion of $\sigma=2$~km~s$^{-1}$, as is typical for SFRs \citep{hillenbrand98,covey06}. This enables the derivation of the encounter rate as a function of $b$ and $v$ \citep{binney87}, which for a given age provides the total number of encounters $n$. Because encounters with pericentres at inclination angles $\theta\simgreat45^\circ$ with respect to the disc plane affect the disc only mildly, about $30\%$ of the $n$ encounters lead to the disc truncation described by Eq.~\ref{eq:rd} \citep{pfalzner05}. We use the probability distribution functions (PDFs) for $b$, $v$, and $m_2$ to calculate the PDF of the `most disruptive' encounter, i.e.~the parameter set that gives the smallest disc size, according {to \citet[Eq.~A5]{maschberger08}. We} then integrate the product of the disc radius $r_{\rm d}$, the PDF of the `most disruptive encounter', and the mass PDF of the perturbed object $m_1$ to obtain the expected disc truncation radius $r_{\rm tr}$. {To} compare the theoretical estimate to the observations, we relate the stellar volume density $\nu$ to the surface density $\Sigma$ as $\nu=\Sigma/2R$, where $R\approx 2$~pc is a typical radius for the SFRs in our sample \citep{hillenbrand98,evans09}.

\section{Results}\label{sec:result}
\subsection{Evidence for environmental effects} \label{sec:theplot}
\begin{figure}
\center\resizebox{9cm}{!}{\includegraphics{theplot.ps}}\\
\caption[]{\label{fig:sizes}
      {{\it Top panel}: Protoplanetary disc (PPD) radius distribution versus ambient stellar density. Squares, triangles, stars, and crosses correspond to ONC sources, T-Tauri stars, Herbig stars, and YSOs, respectively. Dotted lines represent the theoretical truncation for ages of $0.3$ (top) and $1$~Myr (bottom). {\it Bottom panel}: KS test $p_{\rm KS}$-values for the PPD radii at that density to be consistent with the sample at lower densities, using $\{6,7,8\}$ (triangles, squares, diamonds) objects per density bin (see text). The grey thick line shows the $p_{\rm KS}$-value for ONC PPDs only when dividing the subsample in two at the indicated density. The horizontal dashed line and grey area indicate the median $p_{\rm KS}$ and its $1\sigma$ dispersion retrieved from 30,000 Monte Carlo experiments to verify the significance of the results (see text). Dotted lines mark $95\%$ (top) and $99\%$ (bottom) confidence levels.
                 }}
\end{figure}
The upper panel of Fig.~\ref{fig:sizes} {shows the observed {PPD}} radii  {versus $\Sigma$. The} distribution is} relatively insensitive to the ambient density {until $\Sigma\sim10^{3.5}~{\rm pc}^{-2}$, where it} {appears truncated at large radii. This is consistent with the simple theoretical approximation from Sect.~\ref{sec:theory}, which predicts a truncation at these densities for ages between $0.3$ and $1$ Myr} ({see} Fig.~\ref{fig:sizes}). The affected {PPD}s are all {in the ONC. If the truncation} is interpreted as being due to dynamical interactions, the theoretical curves {suggest that youngest sources in the ONC have ages} of $\sim0.6$~Myr, which is reasonably consistent with observations \citep{palla99,dario10,jeffries11}. {However, this is not a unique explanation, since the truncation might also be due to external photoevaporation {by nearby massive stars} \citep[e.g.][]{clarke07}. It is also important to keep in mind that some fraction of the {PPD}s at this projected surface density will actually reside in a region of lower volume density, either behind or in front of the high-density core of the ONC.} The distribution also shows some evidence of a density-independent {upper limit to the {PPD} radius} of $\sim10^3$~AU, which could be related to binarity \citep{artymowicz94,kraus12} or {be intrinsic} \citep{basu98}.

To test the statistical significance of the {change {in} the radius distribution,} we perform a Kolmogrov-Smirnov (KS) test, {which }is theory-independent {and} thus insensitive to model assumptions. {It gives the {probability $p_{\rm KS}$ that two samples were drawn from the same parent distribution}. Starting at the high-density end, we first divide the sample in density bins {of} $\{6,7,8\}$ objects per bin, which {represents a balance between good statistics} and enough bins to resolve the regime where the radius distribution changes}. The KS test is then carried out {comparing radii} in a bin at density $\Sigma_{\rm bin}$ with {those at lower densities (i.e.~$10^{-1}~{\rm pc}^{-2}<\Sigma<\Sigma_{\rm bin}$)}. We do not include bins at $\Sigma_{\rm bin}<10^1~{\rm pc}^{-2}$ to avoid low-number statistics in the reference sample. {The bottom panel of Fig.~\ref{fig:sizes} shows the results of the test}. At intermediate densities, {all} KS tests give high {$p_{\rm KS}$-values, but for densities higher than $\Sigma\simgreat10^{3.5}~{\rm pc}^{-2}$ there is} a pronounced drop. As such, the KS test yields {a detection} of a change in the {PPD} radius distribution {at the $>97\%$ confidence level}. {Note that this does not change when excluding the largest {PPD} in the sample (at $\Sigma\sim200~{\rm pc}^{-2}$). {The result} also holds within the ONC only {(grey} line in Fig.~\ref{fig:sizes}), when} dividing the ONC subsample in two at each density and running a KS test for the radii at both sides of the separation. {To check the result, we performed 30,000 Monte Carlo experiments in which the KS test was applied in the same way to distributions of randomly paired radii and surface densities, i.e. erasing any possible correlation. Figure~\ref{fig:sizes} indicates the resulting median $p_{\rm KS}$ and its $1\sigma$ dispersion, showing that the results obtained for the original sample are unlikely to be due to the adopted statistical method.}

\subsection{Implications for habitable zone occupancy lifetimes}\label{sec:hz}
\begin{figure}
\center\resizebox{8.5cm}{!}{\includegraphics{hz.ps}}\\
\caption[]{\label{fig:hz}
      {Habitable zone lifetime} as a function of stellar mass {for} different ambient stellar densities. The dotted line indicates the total lifetime of the host star, and the diamond indicates Earth.
}
\end{figure}
If the drop of the $p_{\rm KS}$-value at high densities is indeed caused by dynamical truncation, then the model can be used to give the maximum {time} during which the {habitable zone (HZ)} can host a planet or a {PPD} {(`HZ lifetime')}. {We calculate} $r_{\rm tr}$ as a function of $\Sigma$ and $\tau$ as in Fig.~\ref{fig:sizes}, but without averaging over the host stellar mass to retain the mass dependence. We then compare it to the inner radius of the HZ, $r_{\rm HZ}^{\rm in}$, and determine until which age they overlap, {as a function of mass and ambient density}. The radius $r_{\rm HZ}^{\rm in}$ {depends on} the required radiative equilibrium temperature $T_{\rm eq}$, on the {stellar} luminosity $L$ and on the {properties} of the planet \citep[i.e.~$f$, a proxy for atmospheric thermal circulation, and the Bond albedo $A$, see][]{kasting93,tarter07}:
\begin{equation}
\label{eq:rhzmin}
r_{\rm{HZ}}^{\rm in}=\frac{1}{2} \sqrt{\frac{L f (1-A)}{4 \pi \sigma_{\rm{T}} T_{\rm{eq}}^4}} ,
\end{equation}
where $\sigma_{\rm{T}}$ is the Stefan-Boltzman constant. {Following \cite{borucki11}, we calculate $r_{\rm HZ}^{\rm in}$ for an Earth-like planet using $f=1$, $A=0.3$, and $T_{\rm eq}=307~K$, which is the maximum temperature that allows for the presence of liquid water when accounting for the greenhouse effect.} The luminosity is taken from stellar evolution models at solar metallicity \citep{marigo08}. {The results are} shown in Figure~\ref{fig:hz} for ambient densities of $\Sigma=\{10^2,10^3,10^4,10^5\}~{\rm pc}^{-2}$. The HZ lifetime decreases with density and stellar mass, due to the enhanced encounter rate and the large $r_{\rm HZ}^{\rm in}$, respectively. {These estimates for the HZ lifetime hold both for {PPD}s and planetary systems, {since} dynamical interactions {would have comparable effects in both cases} {\citep[e.g.][]{olczak06,parker12}}.}

Figure~\ref{fig:hz} shows that, based on the Earth's existence alone, the solar system cannot have formed in a dense ($\Sigma>10^3~{\rm pc}^{-2}$) {environment}, unless the ambient density decreased on a short ($\tau\mathrel{\copy\simlessbox} 100$~Myr) timescale. Conversely, meteoritic evidence indicates that the young solar system must have endured nearby supernovae \citep[e.g.][]{cameron77}, which provides a lower limit to the product {$\Sigma\tau\simgreat10^{3.8}~{\rm Myr}~{\rm pc}^{-2}$} \citep[see the review by][]{adams10}. {A plausible scenario is thus} that the solar system formed in a massive {($\sim10^4~\mbox{M$_\odot$}$, $\Sigma\sim10^3~{\rm pc}^{-2}$)}, but unbound association, which dispersed on a short ($\tau\sim10$~Myr) timescale \citep[see e.g.][although they refer to such a system as a `cluster']{dukes12}. Our results seem to disagree with \citet{eisner06} who derive the disc fraction {in the ONC} and find no evidence of disc truncations. {However, their conclusion may result from low-number statistics, and PPD mass estimates are more uncertain than radius measurements. On the other hand, our results agree with the studies of \citet{bonnell01} and \citet{spurzem09}, and would explain} why no planets have been found in the globular clusters {47~Tuc and} NGC~6397 \citep{gilliland00,nascimbeni12}, where $\tau\Sigma\sim10^8~{\rm Myr}~{\rm pc}^{-2}$ {within the half-mass radius}. This implies such a high number of encounters that it is {improbable} that {\textit{any}} bound planets survived, {most of them likely to have escaped the cluster due to two-body relaxation \citep[e.g.][]{kruijssen09c}}.

{We} note that our theoretical estimates are conservative and provide upper limits to the disc sizes, because (1) we do not account for potentially higher ambient densities in the past \citep[e.g.][]{bastian08b}, (2) {we neglect the presence of massive stars at earlier ages of the SFRs}, (3) we only consider the most disruptive encounters and {ignore} the cumulative effect of weak perturbations. Figure~\ref{fig:hz} thus provides strict upper limits to the HZ lifetimes.

\section{Further observational avenues} \label{sec:outlook}
We {present evidence} for a change in the {PPD} radius distribution at ambient stellar densities of $\Sigma>10^{3.5}~{\rm pc}^{-2}$ {at the $>97\%$ confidence level}, in line with the expected range due to close encounters with other stars {on a $\sim1$~Myr timescale}. These densities are only reached in the  densest parts of the ONC, which {is not only consistent with the detection of reduced {PPD} masses in the centre of the region \citep{mann10}, but also with} {studies concluding} that encounters are not important in the ONC as a whole \citep[e.g.][]{scally01}. Our results demonstrate that {the} stellar environment {can be} an important factor in setting the habitability of planetary systems. {For instance, the existence of unbound, free-floating planets \citep[see e.g.][]{bihain09,sumi11,strigari12} is a natural outcome of our results}. However, a ubiquity of Earth-like planets in the HZ {of stars remains likely} because a large fraction {($\sim90\%$)} of stars forms in unbound associations {(see \citealt{kruijssen12d} for a recent review, and observational references therein)}, of which the density quickly decreases after their formation.

\begin{figure}
\center\resizebox{8.5cm}{!}{\includegraphics{mdf2.ps}}\\
\caption[]{\label{fig:mdf}
      {{Evolution of the disc mass function (DMF) due to dynamical encounters at an age of 3 Myr for densities $\{2.25\times10^4, 2.5\times10^3, 5\times10^2\}~{\rm pc}^{-2}$ \textbf{(solid, dashed, dash-dotted)}, sampling NGC~3603 at \textbf{radius intervals of $\{0"-5",16.5"-43.5",46.5"-73.5"\}$. The dotted curve} shows the initial DMF, and the vertical dotted line indicates the $4\sigma$ sensitivity limit of ALMA Cycle~1 (30~min at 1.3~mm) at the distance of NGC~3603.}
                 }}
\end{figure}
To verify our results, more observations of {PPD}s in clustered environments are desirable. {A} fruitful approach would be to probe {PPD} truncations using disc mass measurements, {which would provide much larger samples of {PPDs. Figure~\ref{fig:mdf} shows the effect of (only) dynamical encounters on the disc mass function (DMF) using the disc mass loss description from \citet{olczak06} (see Appendix~\ref{sec:app2})}. {This} will be easily observable in dense and young stellar clusters with ALMA. While the full ALMA array will be able to directly measure disc sizes, we predict that, already from Early Science Cycle 1, the sensitivity will be sufficient to detect the variation of the DMF caused by the truncation (see Fig.~\ref{fig:mdf}), which would verify the result of Fig.~\ref{fig:sizes}.}

\begin{acknowledgements}
{We are grateful to the anonymous referee for a thoughtful and constructive report.} We thank Cathie Clarke, Michael Meyer, and Christoph Olczak for insightful comments on the manuscript, and David Jewitt, Thomas Maschberger, and Niels Oppermann for helpful discussions. NB is supported by the DFG cluster of excellence `Origin and Structure of the Universe' and HC by the Millennium Science Initiative, Chilean Ministry of Economy, Nucleus P10-022-F.
\end{acknowledgements}

\bibliographystyle{aa}


\section{Introduction}

High-field  magnetic white  dwarfs  have field  strengths larger  than
10$^6$~G  and  up  to  10$^9$~G \citep{Schmidtetal03}.   The  galactic
population   of   these   white   dwarfs   exhibits   two   remarkable
peculiarities.  The  first one is  that very few  of them belong  to a
non-interacting  binary system  \citep{Kawkaetal07}, while  the second
one    is    that    they    are    more    massive    than    average
\citep{Silvestrietal07}.  Generally  speaking there are  two competing
scenarios to explain the formation of these stars.  One possibility is
that these  white dwarfs  descend from single  stars, so  the magnetic
field is a fossil of previous evolution \citep{Angeletal81}.  However,
this scenario cannot explain  why they are preferentially massive, and
why they  are not found in non-interacting  binary systems.  Recently,
it has  been suggested \citep{Toutetal08,  Nordhausetal11} that strong
magnetic fields can be produced  during a common envelope episode in a
close  binary system  in which  one of  the components  is degenerate.
During this phase, spiral-in of the secondary in the extended envelope
results in differential rotation  in the convective envelope, which in
turn  produces a  stellar  dynamo  which could  give  raise to  strong
magnetic fields.  However,  it has been shown that  the magnetic field
produced  in this  way  does not  penetrate  in the  white dwarf,  but
instead   decays  rapidly   when  the   common  envelope   is  ejected
\citep{PotteryTout10}.

The coalescence of  double degenerate objects has been  the subject of
numerous   theoretical   studies  during   the   last  years.    Large
computational efforts  have been  devoted to study  this intrinsically
three-dimensional problem because it is  thought to be at the heart of
several  interesting astrophysical phenomena,  among which  we mention
Type   Ia  supernova   \citep{Webbink84,   IbenyTutukov84},  magnetars
\citep{Kingetal01},  and   hydrogen-deficient  carbon  stars   and  of
R~Corona Borealis stars \citep{Longlandetal11, Longlandetal12}.  Also,
the  large  metal abundances  found  around  some hydrogen-rich  white
dwarfs with  dusty disks  around them could  also be explained  by the
merger    of   a    carbon-oxygen   and    a   helium    white   dwarf
\citep{Garciaberroetal07}.  Finally, during  the phase previous to the
coalescence of a double  white-dwarf close binary system gravitational
waves are  emitted, and it  has been shown  that LISA will be  able to
detect them  \citep{Lorenaguilaretal05}.  Here we  quantitatevely show
that the merger of two  degenerate cores can also explain the presence
of  very high  magnetic fields  in some  white dwarfs  ---  a scenario
initially suggested by \cite{WickramasingheyFerrario00}.

\section{An $\alpha\omega$ dynamo in the merger remnant}

Full three-dimensional  simulations of the merger of  two white dwarfs
\citep{Guerreroetal04, Lorenaguilaretal09}  have shown that  the final
remnant of the coalescence is a central white dwarf containing all the
mass of the undisrupted primary. On top of this white dwarf there is a
hot  corona  made  of  about   half  of  the  mass  of  the  disrupted
secondary.  Finally,  surrounding this  structure  a rapidly  rotating
Keplerian disk can be found. This disk contains nearly all the mass of
the secondary which has not been incorporated to the hot corona, since
little  mass ($\sim 10^{-3}\,  M_{\sun}$) is  ejected from  the system
during the  merger process.  The structure of  the system is  shown in
Fig.~1.

\begin{figure}[t]
\vspace{9.5cm}
\begin{center}
\special{psfile=fig1.eps hscale=62 vscale=62 voffset=-105 hoffset=10}
\caption{Temperature (solid line,  left scale) and rotational velocity
  (dashed line, right scale) stratifications in the final remnant of a
  double white dwarf merger in  which the binary system is composed of
  two stars  of 0.6 and  $0.8\, M_{\sun}$. The central  spinning white
  dwarf rotates as a rigid body with a rotational velocity $\omega\sim
  0.26$~s$^{-1}$, and  the corona rotates differentially,  with a peak
  rotation  rate $\omega\sim  0.33$~s$^{-1}$.  These  velocities arise
  from energy and angular  momentum conservation, since little mass is
  ejected  from the  system, so  the orbital  angular momentum  of the
  binary  system is  invested in  spinning up  the remnant,  while the
  available energy  is primarily invested in heating  the corona.  The
  location of the convective region is displayed by the shaded area.}
\end{center}
\label{corona}
\end{figure}

The temperature  gradient in the  hot corona is  high, and thus  it is
convective. Using the Schwarzschild  criterion we found that the inner
and  outer edges  of the  convective region  are located  at $R\approx
0.012\, R_{\sun}$, and $R \approx 0.026\, R_{\sun}$, respectively, and
that the  total mass  inside this region  is $\sim 0.24  \, M_{\sun}$.
Assuming  energy equipartition,  the  resulting $\alpha\omega$  dynamo
produces  a magnetic field  $B^2/8\pi\approx \rho(\omega  R)^2/2$. For
the  typical values  found in  our simulations  the magnetic  field is
$B\sim 3.2\times  10^{10}$~G. For this  mechanism to be  efficient the
dynamo must work for several convective turnovers before the energy of
the hot corona is radiated away. The temperature of the corona is very
high,  and hence it  is preferentially  cooled by  neutrinos.  Typical
luminosities  are  $L_\nu\sim 4.0\times  10^2\,  L_{\sun}$, while  the
total thermal  energy of the corona is  $U\sim 8.8\times 10^{48}$~erg.
Consequently,  the  convective  shell  lasts for  $\tau_{\rm  hot}\sim
1.8\times 10^5$~yr.   The convective turnover  timescale is $\tau_{\rm
conv}\approx H_P/v_{\rm  conv}$, where $H_P\approx  2.7\times 10^8$~cm
is the pressure scale  height and $v_{\rm conv}\approx 8.0\times 10^7$
cm/s is  the convective velocity.  Thus,  $\tau_{\rm conv}\sim 3.3$~s,
and during  the lifetime  of the hot  corona the number  of convective
cycles is  sufficiently large.  Thus, the  $\alpha\omega$ mechanism is
able to produce strongs magnetic fields.

Now that we have proved  that large magnetic fields can be established
in the  aftermath of  a double degenerate  merger, it is  necessary to
assess  if  during the  lifetime  of  the  resulting white  dwarf  the
magnetic field diffuses outwards, to the surrounding disk, or inwards,
to the degenerate primary. Thus,  we solved the diffusion equation for
both the disk and the degenerate core, and we found that the timescale
for  diffusion of  the magnetic  field across  the disk  is $\tau_{\rm
disk}\sim 2.0\times 10^{11}$~yr, while for the diffusion timescale for
the central white  dwarf turns out to be  $\tau_{\rm WD}\sim 4.3\times
10^9$~yr.  Thus, it  can safely stated that as  the white dwarf cools,
the magnetic field is not allowed to diffuse across the keplerian disk
nor  to penetrate  in the  degenerate core,  and  consequently remains
confined to the surface layers.

Our model predicts that the masses of high-field magnetic white dwarfs
should be larger than the average  of field white dwarfs and that they
should be  observed as single white dwarfs,  as observationally found.
However, high-field  magnetic white dwarfs  are generally found  to be
slow  rotators \citep{WickramasingheyFerrario00}.   Nevertheless, this
apparent  shortcoming  can be  easily  solved.   If  the rotation  and
magnetic axes  are misaligned, magneto-dipole  radiation rapidly spins
down  the  white dwarf.   The  evolution  of  the rotational  velocity
\citep{Benacquistaetal03} is given by $\dot{\omega} = -2\mu^2 \omega^3
\sin^2\alpha /(3Ic^3)  $, where  $I$ is the  moment of inertia  of the
white dwarf, $\alpha$  is the angle between the  magnetic and rotation
axes and $\mu=BR_{\rm WD}^3$.   Adopting typical values resulting from
our  SPH  simulations  we  obtain  a  spin-down  timescale  $\tau_{\rm
  MDR}\sim  2.4\times10^{8}/\sin^2\alpha$~yr,  when  a field  strength
$B=10^{7}$~G is  adopted.  Hence, if  both axes are  perfectly aligned
the remnant of the coalescence  will be a high-field, rapidly rotating
magnetic  white dwarf.   On  the  contrary, if  both  axes are  nearly
perpendicular magneto-dipole radiation efficiently brakes the remnant.
Consequently,   very  young,   hot,   ultramassive,  slowly   rotating
high-field magnetic  white dwarfs can  also be easily  accommodated in
our model.

On the other hand, spectro-polarimetric observations show that in most
cases high-field magnetic white  dwarfs have fields with both toroidal
and  poloidal components.   Our  scenario can  also  account for  this
observational fact, since  in the $\alpha\omega$ mechanism, convection
is responsible for the generation of poloidal fields, whereas rotation
is responsible  for the generation of toroidal  fields. In particular,
the energy available to generate the poloidal field component is $\rho
v_{\rm conv}^2/2\sim  4.0\times 10^{20}$~erg, which is  $\sim 10\%$ of
the  energy available  to build  the toroidal  component, $\rho(\omega
R)^2/2\sim 5.5\times 10^{21}$~erg.  Thus,  we expect that the magnetic
field geometry of the remnant of the merger will be complex.

\section{The number of high-field magnetic white dwarfs in the solar
neighborhood}

To  assess  the  number of  mergers  that  could  occur in  the  solar
neighborhood    we   expanded   an    existing   Monte    Carlo   code
\citep{Garciaberroetal99, Torresetal02, Garciaberroetal04} designed to
study the galactic population of single white dwarfs to deal with that
of double  degenerates. The population synthesis code  is described in
full  depth  in  \cite{Garciaberroetal12},  where  all  the  necessary
ingredients are detailed.  Here, for  the sake of conciseness, we only
present the most significant results of our calculations.

Our  population synthesis  calculations predict  that the  fraction of
merged  double degenerate  cores in  the solar  neighborhood  is $\sim
2.9\%$ of the  total population.  This number includes  not only white
dwarf  mergers ($\sim  0.3\%$), but  also the  coalescence of  a white
dwarf and a giant star with  a degenerate core ($\sim 1.1\%$), and the
merger of two  giants with degenerate cores ($\sim  1.5\%$).  In these
two last cases the mergers  obviously occur during the common envelope
phase,  while  for  the  cases  in which  two  white  dwarfs  coalesce
gravitational  wave  radiation is  the  final  driver  of the  merger.
Finally, it is  important to mention as well  that the distribution of
remnant  masses  is  nearly  flat,  in  agreement  with  the  observed
distribution     of     masses     of    magnetic     white     dwarfs
\citep{NaleztyyMadej04}.

Within   20~pc    of   the   Sun   there   are    122   white   dwarfs
\citep{Holbergetal08},     and    contains    14     magnetic    stars
\citep{Kawkaetal07}.  Although  scarce, this sample  is 80\% complete,
and allows reliably determine the true incidence of magnetism in white
dwarfs.  Mass determinations are available for 121 of them.  Of the 14
magnetic white dwarfs  in local sample, 8 have  magnetic fields larger
than $10^7$~G,  and 3 have masses  larger than $0.8\,  M_{\sun}$ --- a
value which  is $\sim 2.5\sigma$ away  from the average  mass of field
white dwarfs. This  mass cut is chosen to cull  only white dwarfs that
are  expected to  be the  result  of stellar  mergers.  Actually,  our
calculations  predict that  $\sim 4$  white dwarfs  are the  result of
double  degenerate  mergers,  and   have  masses  larger  than  $0.8\,
M_{\sun}$,   in  agreement   with  observations.    Additionally,  our
simulations  predict that the  fraction of  white dwarfs  more massive
than $\sim 0.8\, M_{\sun}$  resulting from single stellar evolution is
$\sim  10\%$.   Consequently, the  expected  number  of massive  white
dwarfs in  the local sample should  be $\sim 12$.   Instead, the local
sample contains 20, pointing  towards a considerable excess of massive
white dwarfs, which could be the  progeny of mergers.  The rest of the
population of  magnetic white dwarfs  ($\sim 5$) --- those  with small
magnetic field strengths --- would be well the result of the evolution
of single stars \citep{AznarCuadradoetal04}.

\section{Conclusions}

It has  been shown that  the hot, convective,  differentially rotating
corona   predicted  by   detailed   Smoothed  Particle   Hydrodynamics
simulations  of the coalescence  of two  degenerate stellar  cores can
produce  very high  magnetic fields.   We have  also shown  that these
magnetic fields are confined to the outer layers of the remnant of the
coalescence, and do not propagate neither to the interior of the white
dwarf or  to the debris region.   Our model has  two clear advantages,
which meet the requirements imposed  by observations. The first one is
that  it  naturally predicts  that  high-field  magnetic white  dwarfs
should preferentially  have high masses. The second  advantage is that
it  explains why  high-field magnetic  white  dwarfs are  found to  be
single stars.   Moreover, our scenario  predicts that if  the rotation
and magnetic  axes are  not aligned the  magnetic white  dwarf rapidly
spins down due  to the emission of magnetodipole  radiation.  Thus, we
expect  that high-field  magnetic white  dwarfs should  have different
rotation periods,  depending on their  evolutionary status and  on the
geometry of  the magnetic field.  Moreover,  in the case  in which the
masses of the merging objects white dwarfs are not equal the keplerian
disk  can  eventually   form  a  second-generation  planetary  system.
Disruption of small bodies  in this planetary system could contaminate
the atmospheric layers  of the magnetic white dwarf.   This finding is
important,  as it can  explain the  recently discovered  population of
metallic magnetic  white dwarfs.  If,  on the contrary, the  masses of
the  merging  white  dwarfs  are  similar the  remnant  has  spherical
symmetry  and rotates  very rapidly,  as observed  in  some high-field
magnetic  white dwarfs.  Also, the  geometry of  the  surface magnetic
fields  can be well  explained by  our model.   Finally, we  have also
shown that the expected number of double degenerate mergers is roughly
consistent with the number of  high-field magnetic white dwarfs in the
local sample.   In summary, our  calculations indicate that  a sizable
fraction of all  high-field magnetic white dwarfs could  be the result
of double  degenerate mergers, a hypothesis  previously anticipated by
\cite{WickramasingheyFerrario00},  but  not  hitherto demonstrated  by
a quantitative analysis.

\acknowledgements
This research  was supported by AGAUR, by  MCINN grants AYA2011--23102
and  AYA08-1839/ESP, by  the European  Union FEDER  funds, by  the ESF
EUROGENESIS  project (grant  EUI2009-04167), by  AGENCIA:  Programa de
Modernizaci\'on  Tecnol\'ogica BID 1728/OC-AR,  and by  PIP 2008-00940
from CONICET.


\section{Introduction}
This article briefly summarizes works published by \cite{mesz11-aa} and \cite{mesz11-ap}.
It can be shown that for the observed peak-flux or fluence $P(z)$ of gamma-ray bursts (GRBs) it holds:
$
P(z) = \frac{(1+z)^N \tilde{L}(z)}{4\pi d_\mathrm{l}^2(z)},
$
$z$ is the redshift; $d_\mathrm{l}(z)$ is the luminosity distance;
$\tilde{L}(z)$ is isotropic peak-luminosity or emitted energy.
$N=0;\, 1;\, 2$ depending on the units of the observables:
$N=0$ if $P(z)$ is peak-flux in units erg/(cm$^2$s) and then $\tilde{L}(z)$ in units erg/s;
$N=1$ if $P(z)$ is fluence in units erg/cm$^2$ and then $\tilde{L}(z)$ in ergs or
$N=1$ if $P(z)$ is peak-flux in units ph/(cm$^2$s) and then $\tilde{L}(z)$ in units ph/s;
$N=2$ if $P(z)$ is fluence in units ph/cm$^2$ and then $\tilde{L}(z)$ is in photons.
$P(z)$ is given by photons with  energies $E_1 \leq E \leq E_2$, where
$E_{1,2}$ is the detector range. $\tilde{L}_\mathrm{ph}(z)$ is from
energy range $E_1(1+z) \leq E \leq E_2(1+z)$.

It is a standard cosmology that for small redshifts ($z \ll 0.1$) $d_\mathrm{l}(z)\propto z$, and for high redshifts
$\lim_{z \rightarrow \infty} \frac{d_\mathrm{l}(z)}{1+z} =$ finite and positive number
for any $H_o$, $\Omega_M$, $\Omega_{\Lambda}$.
If we assume that $\tilde{L}(z) \propto (1+z)^{q}; N+q>2$ for $z \rightarrow \infty$
then $\frac{dP(z)}{dz} > 0$, for $z \rightarrow \infty$.
If $\tilde{L}(z)$ increases with $z$ faster than $\propto (1+z)^{2-N}$ for high redshifts,
then inverse behaviour can happen, and the apparently brighter GRBs can be at higher redshifts; $dP(z)/dz > 0$
can occur.

For small redshifts $z$, $P(z)$ decreases as $z^{-2}$, but if
$ \tilde{L}(z) \propto (1+z)^{q}$ and $N+q > 2$
then $P(z)$ starts to increase as $z^{N+q-2}$ for high $z$.
To answer the question where the $z_\mathrm{turn}$ is, i.e. where $dP(z)/dz = 0$, one can
simply search for the minimum of $Q(z) = (1+z)^{N+q}/d_\mathrm{l}(z)^2$, i.e., when $dQ(z)/dz = 0$.

\section{Samples of long GRBs}

We studied 134 long Swift GRBs with know redshifts. The distribution of observed flunces/peak fluxes together with
the curves of constant $\tilde{E}_\mathrm{iso}$ and $\tilde{L}_\mathrm{iso}$ are shown in Fig.\,\ref{fig1}.
We also studied 6 Fermi GRBs, 9 BATSE GRBs with measured redshifts and 13 BATSE GRBs with calculated pseudo-redshifts.
These samples also demonstrate that fainter bursts can well have smaller redshifts (for details see \cite{mesz11-aa}).

\begin{figure}
\begin{center}
\includegraphics[width= {0.95\columnwidth}]{ripa_fig1.eps}
\caption{Distributions of fluences (left panel) and peak-fluxes (right panel) Swift long GRBs.
On the left panel the curves denote the values of fluences for $\tilde{E}_\mathrm{iso} = \tilde{E_o}
(1+z)^q $ (constants $\tilde{E_o}$ are in $10^{51}$\,erg: I. 0.1; II. 1.0;
III.  10.0). On the right panel the curves denote the values of peak-fluxes for $\tilde{L}_\mathrm{iso} =
\tilde{L_o} (1+z)^q$ (constants $\tilde{L_o}$ are in $10^{58}$\,ph/s: I.
0.01; II. 0.1; III. 1.0). Here $N=1$ thus value $q=1$ is the limiting case.}
\label{fig1}
\end{center}
\end{figure}

\section{Conclusion}

The theoretical study of the $z$-dependence of the observed fluences and peak-fluxes of
GRBs have shown that fainter bursts could well have smaller redshifts.
This is really fulfilled for the four different samples of long GRBs.
These results do not depend on the cosmological parameters and on GRB models.

\acknowledgements
This study was supported by the OTKA grant K77795, by the Grant Agency of the Czech Republic grants No. P209/10/0734, by the Research Program MSM0021620860 of the Ministry of Education of the Czech Republic, and by Creative Research Initiatives (RCMST) of MEST/NRF and the World Class University grant no R32-2008-000-101300.