\section{Introduction}
\label{Motiv_SDI}
Since the discovery of the first extrasolar planet around 51 Pegasi by \citeauthor{Mayor95} in 1995, more than 480 extrasolar planets (exoplanets) in more than 410 planetary systems have been identified as of August 2010\,\footnote{http://exoplanet.eu/catalog.php}. However, most of these detections were made using indirect detection methods like radial velocity (RV) surveys, the transit method or microlensing events. Further, the RV technique can in general only provide minimum masses (M\,sin\,$i$) and thus limited knowledge of the astrophysical properties of the companion.
From the beginning, these exoplanets revealed a remarkable diversity in their physical and orbital characteristics like separation, eccentricities and masses and thus challenged the planet formation scenario. In particular, the period-mass distribution shows tendencies which are hard to explain within the framework of current models (e.\,g.\,massive "hot Jupiter" were found at very small orbits, even down to a semi-major axis of only 0.02\,AU). A revision of the standard scenario seemed to be required.
However, the discovery space of RV searches does not yet cover planets at distances larger than about 5\,AU (for statistics see \citealt{Marcy03}). The question is, if the observed period-mass distribution is, at least partly, a selection effect (higher sensitivity at shorter orbital periods) and thus might provide only an incomplete picture of the overall planet population, or if giant planets are really rare at larger separations. In addition, the indirect RV and microlensing techniques do not enable any photometric or spectroscopic information of the exoplanet companions and their physical properties such as brightness, color, effective temperature or composition.

Therefore, a lot of attention has been directed towards direct imaging of planetary mass objects (PMO) in the last years. A direct detection would enable extensive follow-up studies to characterize the planets, but the main challenge of this method is to deal with the huge brightness contrast and the small angular separation between the faint planet companion and its host star. With the possibility of high-resolution space-based observations and the improvement of adaptive optics (AO) systems on large ground-based telescopes, the necessary sub-arcsecond spatial resolution became achievable. In addition, planetary mass objects and brown dwarfs are hotter at younger ages and hence the brightness contrast between the primary and a possible companion will be much smaller. Thus, in order to push the direct detection threshold down to the planetary mass regime, systematic searches focused especially on young stars in nearby associations.

First extensive \emph{Hubble Space Telescope} (HST) programs to search for substellar companions to post-TTauri stars  \citep{Brand00}, as well as to very low mass (VLM) stars and brown dwarfs in the Pleiades open cluster \citep{Mart00} did not reveal any companion at intermediate orbits. Since then, several other extended HST surveys of free-floating ultra-cool dwarfs and brown dwarfs showed that roughly 20$\%$ of them have substellar companions. Almost all of them constitute close to equal-mass systems (e.\,g.\,\citealt{Bouy03, Reid01, Burgasser03}) and hence, rather resemble more massive, stellar binaries than planetary systems. Consequently, some interesting questions arose: \\[2mm]
\hspace*{2cm}\parbox{14cm}{\textsl{Do "single" brown dwarfs evolve similar to single stars?}}\\[1mm]
\hspace*{2cm}\parbox{14cm} {\textsl{Are they capable of forming planetary systems?}}\\[2mm]
The detection of circumstellar disks around young brown dwarfs revealed that they at least posses the raw material to form planets (e.\,g.\,\citealt{Apai02, Pascu03, Liu03, RayJay03}; later confirmed by results of e.\,g.\,\citealt{Apai05, Luhman06, Luhman08_02, Scholz08}, see also reviews in \citealt{Luhman07} and \citealt{Henning08}).\\
On that account, we started in 2004 the first spectral differential imaging (SDI) survey with HST to directly detect planetary mass objects in close orbits around nearby free-floating brown dwarfs.  A direct image detection would have important implications on our understanding of planetary mass objects, their origin, the formation of substellar bodies in general and binarity in the brown dwarf domain in particular.

Brown dwarf multiplicity properties exhibit significant differences in frequency, separation and mass ratio distribution compared to Sun-like stars. While the binary fraction for Sun-like stars is $\approx$\,57\% \citep{Duquennoy1991}, this fraction steadily decreases to $\approx$\,26\,-\,42\,\% for early- to mid-type M dwarfs (e.g. \citealp{Reid1997, Bergfors10}) down to $\approx$\,10\,-\,30\,\% for the very low mass stars and brown dwarfs (e.g. \citealp{Bouy03, Burgasser07_2, Joergens08, Goldman08}). The same steady decline is true for the semi-major axis distribution which is spread over a wide range of separations for G dwarfs with a peak around 30\,AU \citep{Duquennoy1991}, smaller separations for M dwarfs where the distribution peaks between 4 and 30 AU \citep{Fischer1992}, and a very narrow separation distribution for the brown dwarfs which peaks at 3\,-\,10\,AU (e.g.\,\citealt{Allen07}). The mass ratio distribution for very low mass stars and brown dwarfs reveals a clear tendency to equal mass system components with a distinct peak around \emph{q}\,$\sim$\,1. In contrast, low \emph{q} values are much more common at all separations for binary systems of solar-type stars (see review by \citealt{Burgasser07_3} and references therein).

Recent analysis of radial velocity surveys confirm a 6\,-\,8\,\% frequency of planetary mass companions to Sun-like stars (e.g. \citealt{Grether06} and \citealt{Mordasini09} with references therein) at relatively small separations which peak around 1 AU \citep{Udry07}. While commonly observed mass ratios between solar type stars and substellar (planetary) mass objects range between 1:1000\,-\,1:100 (cf.\,Extrasolar Planets Encyclopedia\,\footnote{http://exoplanet.eu/catalog.php})  one would expect a much higher mass ratio of 1:20\,-\,1:5 for a system consisting of a planetary mass object orbiting a brown dwarf due to the lower initial mass of the primary. \citet{Apai08} point out that depending on the planet formation process, the frequency for lower-mass planets is expected to be higher than for solar type stars. For the semi-major axis distribution of such a system one would also expect relatively small separations of a few AU due to the much smaller sizes of the circumstellar disks around young brown dwarfs compared to solar type stars.

\section{Observing strategy}
\label{Obs_strat_SDI}
Because of the higher achievable contrast and the lack of suitable AO reference stars, the observations were conducted from space with HST as explained in the following.
The observing strategy that suited best to reach our goal of detecting planetary mass objects around brown dwarfs and to overcome the remaining brightness contrast, was the spectral differential imaging (SDI) technique \citep{Smith87}. This technique takes advantage of the fact that a cool ($\le$\,1300\,K) giant planet or low-mass brown dwarf (i.\,e.\,T dwarf) is much fainter in certain molecular bands due to strong absorption than in the neighboring continuum.
Already \citet{Allard95} stated that the spectral energy distribution (SED) of brown dwarfs is very peculiar. The molecular opacities, which globally define the continuum of brown dwarfs, cause the SED to peak around 1.1\,$\mu$m for solar metallicities, almost independent of their effective temperature. Further, the observed \emph{J}--\emph{K} colors for L dwarfs are red (about +\,2 mag), while they are blue (\,-\,0.1 to +\,0.4 mag) for the T dwarfs. This is caused by dust grain formation \citep{Allard01}. Therefore, the brightness difference of any L dwarf\, -- T dwarf pair is smallest around 1.1\,$\mu$m, making the \emph{J}\,-\,band the best spectral regime to detect T\,-\,type or even cooler companions to L dwarfs.

Figure~\ref{SDI_bandpasses} shows the spectral features of L and T dwarfs in this wavelength regime in detail. It reveals that L dwarfs have a comparative flat SED between 1.06\,$\mu$m and 1.15\,$\mu$m, while the flux of a T dwarf drops significantly at the beginning of the molecular water absorption band at $\sim$ 1.10\,$\mu$m. The T dwarf plot in this Figure~\ref{SDI_bandpasses} additionally includes a synthetic spectrum from \citet{Burrows03_2} for a 5\,$M_\mathrm{Jup}$, 1 Gyr old planetary mass object with T$_{\mathrm{eff}}$\,=\,312\,K and shows, that this effect is still present at these very low effective temperatures. Thus, a subtraction of two images obtained simultaneously, one in and one off, yet near this molecular band, will cancel out most of the equally bright PSF structure of the primary (L dwarf), revealing the much fainter signal from a cool, possibly planetary mass companion (see \citealt{Rosenthal96}; \citealt{Racine99}).
The narrow-band filters F108N and F113N of the NICMOS1 camera (NIC1) on board HST match particularly well the relevant bands (as shown in Figure~\ref{SDI_bandpasses}) and are ideally suited for the differential imaging technique.

\section{Sample selection}
Given the angular separation and luminosity ratio problem described above, the list of targets was assembled with the focus on brown dwarfs in proximity to the Sun and with a relatively young age. A closer distance to the Sun translates in a larger angular separation for a given physical separation between the brown dwarf and an eventual planetary-mass object, thus companions on smaller orbits can be detected. In addition, the detection threshold caused by sky background decreases for a fainter planetary mass object, since its apparent brightness becomes larger, the closer it is to the Sun.

Concerning the age, priority was given to targets with estimated ages younger than 1 Gyr. Giant planets, as well as young brown dwarfs are more luminous when they are young due to the remaining gravitational contraction energy from their formation process. Left without a resource of producing new energy they cool down rapidly and become dimmer with time. Based on evolutionary models and synthetic spectra, a young (100 Myr) giant planet or brown dwarf is $\sim$\,100\,-\,200 times more self-luminous than a 5 Gyr old one (e.\,g.\,\citealt{Burrows97}; \citealt{Baraffe03}). Hence any planetary mass companion to an old brown dwarf would be relatively cool and too faint for detection.

The age estimation of the targets in this survey was not straightforward since all of them are field brown dwarfs without any hints of being a member of an association or moving group. In general, \citet{Allen05} showed that proceeding down the spectral sequence from L to T dwarfs, the average age decreases and the relative proportion of young dwarfs ($\le$\,2 Gyr) increases if compared to M dwarfs. This is mainly caused by the lack of stable hydrogen burning in the brown dwarf core and therefore their rapid cooling through these temperature regimes with time \citep{Burrows01}. Our targets were selected according to spectral age indicators, such as the presence of lithium in the atmosphere. The detection and strength of lithium absorption is a clear confirmation of the substellar nature of these objects, as no lithium is expected for masses $\ge$ 65 $M_\mathrm{Jup}$ due to lithium destruction by proton capture in the deeper layers of the convection zone (e.g. \citealp{Bodenheimer1965, Nelson1993}). Combined with the still relatively high effective temperature representative for the spectral type L this is a clear sign of youth $\le$\,1 Gyr (\citealt{Basri98_2}; \citealt{Magazzu93}). Moreover, \citet{Gizis00} point out that chromospheric activity is primarily observed in older and hence stellar mass L dwarfs. They also find a clear anti-correlation between the presence of lithium and H$_\alpha$ emission. Thus, early L dwarfs with a very weak or non-existent H$_\alpha$ emission line are more likely young brown dwarfs rather than stars.

The final sample was compiled out of  $\sim$\,250 L dwarfs cataloged in D.\,Kirkpatrick's online archive\footnote{http://www.DwarfArchives.org} as of 2004 (start of this project) and consists of 12 young L dwarfs which fulfilled the following requirements:
\begin{itemize}
\item the object is within 30\,pc to the Sun
\item the \ion{Li}{1} absorption at 6708 \r{A} is present and therefore a clear confirmation of their substellar nature, and indication of youth.
\item lack of strong chromospheric activity (H$_\alpha$ emission) which is preferentially observed in stellar mass L dwarfs.
\item the objects are isolated in the sense of no known close companion so far.
\end{itemize}

These age and proximity criteria helped to enlarge the range of mass and separation over which the survey would be sensitive and to maximize the discovery probability. The properties of all selected targets are summarized in Table~\ref{Diff_imag_targets}. Although there are no Lithium measurements available for the two mid\,-\,L dwarfs 2MASSW\,004521+1634 and 2MASSI\,0835--0819, they were included since they belonged to the 10 closest L dwarfs known and because of their brightness compared to other L dwarfs of the same spectral type. A similar decision was made for Kelu-1 which had only a weak detection of the Lithium absorption line but was apparently overluminos \citep{Mart99_2, Leggett01, Goli04_2} and showed photometric and spectroscopic variability \citep{Clarke02, Clarke03}. Furthermore, the nearby T dwarf 2MASSW\,055919--1404 was added to calibrate the spectral differential imaging method for HST/NICMOS.

\section{Observations and SDI data reduction}
The observations of the twelve isolated L dwarfs and the calibration T dwarf were scheduled as HST program GO 10208 (PI: W. Brandner) from September 2004 to July 2005. The sources were observed with the NICMOS1 (NIC1) camera, providing a high-resolution pixel scale of 0\,\farcs0432 and a field of view (FoV) of 11\arcsec \,x 11\arcsec, in the two narrowband filters F108N (1.08 $\mu$m) and F113N (1.13 $\mu$m). The observation sequences\- consisted of 4 exposures at two different detector positions in each filter. This two\,-\,point dither pattern facilitates an optimized PSF sampling, effectively rejects bad pixels and gains redundancy against cosmic ray events. Altogether, two orbits per target were necessary to achieve the required S/N ratio of 6\,-\,8 for a limiting magnitude of \emph{J}\,=\,20 mag in the narrowband filters, and a brightness difference of $\approx$ 6 mag at a separation of 0.3$\arcsec$. All data were acquired in MULTIACCUM mode and the total integration times per filter were 2560\,s (F108N) and 2816\,s (F113N), respectively. In order to limit possible HST breathing variations in the PSF, differential imaging observations in the F108N and F113N filters were obtained in each of the orbits.

For the analysis of the data sets, the HST pipeline reduced images were used, followed by an additional bad pixel masking in order to compensate for the non-optimal bad pixel mask available in the pipeline. The SDI reduction was then accomplished with a custom made IDL program: First, a square aperture around the brown dwarf is extracted from all four images in each filter. To make use of the full integration time, the extracted frames were then co-added after they were aligned with respect to the first image of the related filter. This alignment was done with a 2D cross-correlation and a FFT-shift based combination, which provides a sub-pixel accuracy. In the next step, a sub-pixel re-sampling was performed in order to transform the observation in the F108N filter to the same $\lambda$/D scale as in the F113N filter (where $\lambda$ is the respective observing wavelength and \emph{D} is the HST primary mirror diameter). The final difference between the two filters was then calculated after additional alignment of the two master images with the shift algorithm, followed by a flux calibration to account for any flux losses during the previous re-scale process. In cases where a clear positive residual remained in the reduced image, the difference was also obtained separately for the individual dither positions (hence only two images per filter). This provides two independent detections of a real signal, or the exclusion, if it only exists at one position and is most likely caused by a cosmic ray remnant or hot pixel buried in the central core of the PSF.

\section{Results}
The visual examination of the reduced images reveals that for nine of our targets the residuals are at the level expected from photon noise and no significant positive signal can be detected.
 The final images of these data sets are shown in Figures \ref{Diff_imaging_result_1} and \ref{Diff_imaging_result_2}. The remaining single positive signals in some of the targets in Figure \ref{Diff_imaging_result_1}  just close to the original PSF center, are examples for residuals caused by bad pixels which could not be corrected by the additional bad pixel masking, since they were buried in the bright core of the PSF. They are not visible in both of the ``single" position SDI images simultaneously and thus do not correspond to a real signal of a possible faint companion. For the four brightest and closest sample targets 2MASSW\,0045+1634, LSR\,0602+3910, 2MASSI\,0652+4710 and 2MASSI\,0835-0819 (Figure \ref{Diff_imaging_result_2}) the final SDI results show systematic mirrored dark and bright features in the $\sim$ 6\,-\,7 central pixels of the images. These are common residuals in contrast-limited cases and are most probably caused by the systematic variations in the PSF structure due to differential aberrations of the filters. Since the other five targets are much fainter and their residuals are mostly dominated by photon noise, this pattern is no longer visible. If they were brighter their residuals would look the same.

To asses the overall sensitivity for the survey, the standard deviation of the background residuals after PSF subtraction was calculated for each data set as a function of radial separation from the host brown dwarf. The limits are then based on the median of the statistical errors at the various radii and the primary F108N magnitude. Figure~\ref{SDI_sensitivity_mag} shows the derived 3$\sigma$ narrowband detection limits. The bold solid line represents the mean of the data sets, while the shaded area corresponds to the standard deviation. At a separation of 0.05$\arcsec$ from the primary, a mean contrast of $\Delta$\,F108N\,=\,3.5 mag can already be achieved. For the brighter brown dwarfs in the survey we reach a contrast as large as $\sim$\,6.7 mag at separations $\ge$\,0.4$\arcsec$ from the primary. For the fainter brown dwarfs in the sample we still reach a contrast of $\sim$\,5.4 mag for separations $\ge$\,0.4$\arcsec$.
In order to determine the lower mass limit sensitivity of the survey, one can translate the contrast (in combination with the brightness of the primary) into a mass using the theoretical  COND evolutionary models from \citet{Baraffe03}. First of all, we converted the known \emph{J}-\,band magnitudes of the primary L dwarfs into \emph{Y}-\,band magnitudes. \citet{Hillenbrand02} and \citet{Hewett06} showed that the \emph{Y}\,--\,\emph{J} colors remain relatively constant at 1.0\,$\pm$\,0.15 mag from early-L through late-T type dwarfs. Thereafter, we applied the mass-luminosity relationships of the \emph{Y}-\,band (I.\,Baraffe, private communication) which are closest to the F108N bandpass. Figure~\ref{SDI_sensitivity_mass} illustrates the resulting mass limits for 1 Gyr (left image) and 0.5 Gyr (right image). Similar to Figure~\ref{SDI_sensitivity_mag}, the bold solid line corresponds to the mean value of the individual observing runs and the shaded areas represent the standard deviation.

Altogether, the smooth results after the SDI reduction indicate the absence of a planetary mass companion with a minimum mass of 6\,-\,11\,$M_\mathrm{Jup}$ at an age of 1 Gyr and between 0.07$\arcsec$ and $\sim$\,1.4$\arcsec$ around the brown dwarfs ($\sim$ 1 - 40 AU, depending on their distance to the Sun). This is an area usually most strongly affected by PSF residuals for non-differential imaging observations. For primaries with an age of 0.5 Gyr our observations even exclude any planetary mass object down to 5\,-\,7\,$M_\mathrm{Jup}$. This result validates our observing and data analysis strategy. For the three remaining L dwarfs, the observations revealed very interesting results which will be discussed in the following sections.

\subsection{Two newly resolved brown dwarf binaries}
\label{Diff_binaries}
Even though we tried to avoid known binaries during the sample selection process, the observations revealed that two of the targets are close brown dwarf binary systems. Figure \ref{HST_binaries} shows the NICMOS detection images where both systems are clearly resolved in two components.

\subsubsection{Kelu-1\,AB}
Kelu-1 was discovered to be a binary system by \citet{Liu} and also independently resolved by \citet {Gelino}, both with AO from the ground and just four months before our own HST observation in July 2005. This solved the by then unexplained overluminosity of the ``single" object. The HST observation presented here, is now the first high-resolution confirmation of the physical association of both components, by proofing that they are a common proper motion pair. The system showed clear evidence for orbital motion with a significant increase in separation by 15\,mas in four months out to 299.8 $\pm$ 0.2\,mas and a position angle (PA) of 221.33$\degr$ $\pm$ 0.04$\degr$, hinting at  a short orbit binary. This makes Kelu-1\,AB a good target for a monitoring program. The extensive follow up observations derived with the \emph{Very Large Telescope} (VLT) and all results concerning this system are discussed in Stumpf et\,al.\,(\citeyear{Stumpf09}; 2010b, submitted).

\subsubsection{2MASSW\,0310+1648\,AB}
The second discovered very low mass binary is 2MASSW\,0310+1648\,AB. With an original spectral type of L9, this brown dwarf belongs to the still peculiar L/T transition objects which exhibit several unusual characteristics. Resolved for the first time  as a binary system with almost equally bright components, it adds another system to the already higher binary fraction among the brown dwarfs that span the transition between the L and T spectral class (e.g.\,\citealp{Burgasser06_1, Burgasser10_1}). The first astrometric measurements yield a separation of 204.3 $\pm$ 0.4\,mas and a PA of 206.4$\degr$ $\pm$ 0.1$\degr$ on 2004 September 24. Including follow up observations, this target is discussed in \citet{Stumpf10_1}.

\subsection{A planetary companion candidate}
\label{Diff_companion}
The HST observations of the L4.5 brown dwarf 2MASSW\,033703-175807 (hereafter 2M0337-1758) revealed a clear residual signal in the final SDI reduced image, close to the center of the brown dwarf. Unlike for all other L dwarfs in this survey, this signal is also visible in the two individually reduced detector position images. Therefore, it cannot be caused by a single residual bad pixel. Further, with a S/N ratio of 7 at its peak in the final image, it is well above the noise level. Figure~\ref{SDI_candidate_image} displays the resulting images with the single detector position results in the middle and the final result, including all observation exposures, at the bottom.

\subsubsection{Evidence for binarity}
The F108N image of 2M0337-1758 and its final SDI reduced image were used to calculate the magnitude difference between the brown dwarf and the companion candidate, as well as their separation and position angle. This was achieved by applying aperture photometry with the IRAF \emph{phot} routine in the \emph{apphot} package.
From the identified centroid positions of the photometric apertures a separation of 0.087$\arcsec$\,$\pm$\,0.015$\arcsec$ (corresponding to 2.52\,$\pm$\,0.44\,AU at a distance of 29\,pc) and a position angle of 195.6$\degr$ $\pm$ 5.5$\degr$ was calculated.

For the photometry, the preliminary assumption that the flux of the companion only contributes to the F108N filter and not the F113N filter was necessary, since the result of the applied SDI technique does not provide any information about the actual flux distribution of the companion candidate in each of the filters. This approximation is roughly valid for very late T dwarfs and even later spectral type objects (see Figure \ref{SDI_bandpasses}), but implies that the derived magnitude in F108N for the companion is only a lower limit.
For further comparison with other observations and theoretical models, the instrumental count rates were converted into the Vega magnitude scale after aperture correction and using the most recent photometric keyword-value as provided by the STScI webpage\footnote{http://www.stsci.edu/hst/nicmos/performance/photometry/postncs\_keywords.html}.  Using this method, we derived a F108N magnitude of 16.54\,$\pm$\,0.02 mag for 2M\,0337-1758\,A, a brightness difference of 3.88\,$\pm$\,0.03 mag between the brown dwarf and its companion candidate and hence a F108N magnitude of 20.42\,$\pm$\,0.04 mag for the candidate. With a system distance of 29\,pc, this corresponds to absolute magnitudes of M$_{F108N}$\,=\,14.23\,$\pm$\,0.02 mag and 18.11\,$\pm$\,0.04 mag for A and B, respectively.

Together with the age estimation of the primary brown dwarf, the comparison of the derived values with the COND evolutionary models for brown dwarfs and EGPs (extrasolar giant planets) by \citet{Baraffe03} provides a first approximation of the mass and T$_{\mathrm{eff}}$ for the companion candidate. Since the F108N wavelength range lies much closer to the \emph{Y}-\,band ($\lambda_{\mathrm{C}}$\,=\,1.03\,$\mu$m) than to the \emph{J}\,-\,band, it was assumed that the objects are equally bright in F108N and \emph{Y}. For the following comparison with the isochrones, an extended version of the above mentioned evolutionary models including the \emph{Y}-\,band (I.\,Baraffe, private communication) was used. A comparison of the derived mass and effective temperature for the primary 2M0337-1758\,A with those one would derive with its 2MASS \emph{J} magnitude, proves the assumption of M$_{F108N}$ = M$_{Y}$ to be acceptable. The result is given in Table \ref{Mass_est_HST_comp}. For the companion candidate we derive a minimum mass of 10\,-\,15 $M_\mathrm{Jup}$ with a T$_{\mathrm{eff}}$ between $\sim$ 600\,K and $\sim$ 630\,K  for an age of 0.5 and 1 Gyr, respectively. To investigate how these results would change if the companion magnitudes in F108N and \emph{Y} are in fact not equally bright, we can assume two different scenarios as one can see in the right image of Figure~\ref{SDI_bandpasses}. In the two relevant bandpasses the flux of the T dwarfs peaks right in the F108N filter and these objects will therefore appear brighter in F108N than in the \emph{Y}-\,band. Assuming an offset of $\sim$\,0.5 mag for late T dwarfs, similar to the offset \citet{Biller06} derived between the \emph{H}-\,band and their SDI narrow band, the minimum mass of the companion would decrease down to 9\,-\,12 $M_\mathrm{Jup}$. In contrast, for a planetary mass object with 5 $M_\mathrm{Jup}$ at 1 Gyr the water absorption band extends into the F108N filter and this object would appear dimmer than in the \emph{Y}-\,band. Here, the offset might not be as large but to derive a valuable upper error limit, we use the same value. This results in an increase of the minimum mass to 11\,-\,17 $M_\mathrm{Jup}$ for ages of 0.5 and 1 Gyr, respectively. The results for the companion candidate are also summarized in Table \ref{Mass_est_HST_comp}, including the larger uncertainties due to the different assumptions.

As of fall 2009, there are only nine T dwarfs known, which are classified as T8 or later\footnote{http://www.DwarfArchives.org} and only five of them have a T$_{\mathrm{eff}}$ $\sim$ 600\,K (\citealt{Warren07}; \citealt{Delorme08}; \citealt{Burningham08}; \citealt{Burgasser08_2}; \citealt{Burningham09}). Hence, the companion adds up to a still very small number of just recently detected very cool brown dwarfs and the mass estimation places it right at the theoretical mass limit of $\sim$\,13 $M_\mathrm{Jup}$ for deuterium burning, which is mostly used as the boundary to distinguish between brown dwarfs and giant planets.
Nevertheless, the dominant uncertainty in translating the derived magnitude difference to brown dwarf and even planetary masses, is the uncertainty of the evolutionary models, which are not yet fully calibrated by observations, especially the mass-luminosity relation as a function of age (see further discussions below).

\subsubsection{A possible background object?}
To clarify if the detection might be caused by a background object, the position of 2M0337-1758 on previous HST/WFPC2 images was analyzed. These images were taken on 2000 December 25  in the broadband filter F814W and the medium band filter F1024M. No source could be detected in the WFPC2 frames above the noise level (corresponding to m$_{F814W}$\,=\,24\,mag) at the 2004 NICMOS position of the brown dwarf, but the comparison of the two positions allowed the first proper motion determination for 2M\,0337-1758 with $\mu_{\alpha}$\,cos\,$\delta$ = 0.247\,$\pm$\,0.02$\arcsec$/yr and $\mu_{\delta}$ = -\,0.024\,$\pm$\,0.028$\arcsec$/yr. This result will also be crucial to confirm the common proper motion nature of the binary in follow-up observations.

Analyzing the probability of a background source from a more physical side, we can already exclude a faint, stellar background source (like a M dwarf). These objects have, like L dwarfs, a relatively flat continuum between 1.08 -\,1.13\,$\mu$m and would be therefore not detected by our SDI technique. On the contrary, this method does not exclude a distant galaxy or a quasar. Based on the fact that we have a point-like remaining signal with an approximate magnitude of \emph{J}\,$\sim$\,19.5 mag, the candidate is unlikely a ``normal" galaxy and even too bright to be an active galactic nuclei (AGN). However, we cannot exclude a quasar or other peculiar background object (e.g.\,Wolf-Rayet galaxy). Therefore, the next step will be to establish common proper motion of the two components and to exclude any background object hypothesis.

\section{Discussion}
The direct detection of planetary mass objects still remains a challenging endeavor. In the past few years many systematic high-contrast direct imaging surveys for planetary mass companions around a large number of young, nearby stars have been conducted with HST (\citealt{Luhman05}; \citealt{Lowrance05}) and with AO at large ground\,-\,based telescopes (e.\,g.\,\citealt{Masciadri05}; \citealt{Metchev06}; \citealt{Lafreniere07}; \citealt{Biller07}; \citealt{Kasper07_2}, the latter two employing SDI). However, only very few planetary mass or very low mass brown dwarf companions have been discovered in these surveys. The first directly imaged planetary mass companion was detected around the young brown dwarf \object[2MASS J12073347-3932540]{2MASS\,1207-3932} (hereafter 2M1207) by \citet{Chauvin04, Chauvin05_01} with an estimated mass of 5\,-\,8 $M_\mathrm{Jup}$ and a separation of $\sim$ 55 AU. Shortly after, several more planetary mass object detections were reported to very young primaries with an age $\le$ 60 Myr: GQ Lupi\,b (\citealt{Neuh05}), AB Pic\,B (\citealt{Chauvin05_2}), DH Tau\,B (\citealt{Itoh05}; \citealt{Luhman06_2}), Oph\,1622-2405\,B (\citealt{RayJay06}; \citealt{Close07}), CHXR\,73\,B (\citealt{Luhman06_2}) and 1RXS\,J1609-2105\,B (\citealt{Lafreniere08}). With mass estimates between $\sim$\,7\,-\,20 $M_\mathrm{Jup}$ (thus at or just above the brown dwarf/planet boundary) and projected separations of 55\,-\,300 AU from their primary, their nature is still extensively debated and the exact number of planetary mass companions remains controversial. Such large separations suggest a formation process more similar to brown dwarfs (via fragmentation) rather than giant planets (via core accretion, out to a maximum distance of $\sim$ 10\,-\,15 AU; or via gravitational disk instability which demands a very high massive initial disk). Alternatively, planets could form close-in to the star and then migrate outward to a larger radii. \citet{MartinR07} show that such an outward migration could in principle be possible.
Recently, \citet{Marois08} announced the first directly imaged multiple planetary system (HR\,8799\,b,c,d) with masses between $\sim$ 7\,-\,10 $M_\mathrm{Jup}$ at separations between 24 and 68 AU. At almost the same time, \citet{Kalas08} detected the so far lowest mass companion Fomalhaut\,b ($\sim$ 3 $M_\mathrm{Jup}$ ) and \citet{Lagrange09} claimed the discovery of the closest companion ever directly imaged to a star ($\beta$\,Pictoris\,b with $\sim$ 8 AU). While Fomalhaut\,b is again at a relatively large separation, HR\,8799\,b,c,d as well as $\beta$\,Pictoris\,b, are the first directly imaged companions around a star, whose formation could be explained by the common core accretion or disk instability scenarios.

But not only the formation process hinders the exact assignment of all these companions to the planetary regime. So far, the masses of these planet candidates are entirely based on evolutionary models. These models have been developed in significant detail over the past years (\citealt{Burrows97, Burrows05}; \citealt{Baraffe03, Baraffe08}) but remain purely theoretical, especially for planets and very young ages due to the lack of direct detections of well\,-\,characterized giant exoplanets. For example, the original mass estimate for GQ Lup\,b ranged from 1 to 42 $M_\mathrm{Jup}$ (\citealt{Neuh05}). This wide range in mass estimate is primarily caused by the use of different theoretical evolutionary models (with different initial conditions, see \citealt{Burrows97, Baraffe02, Wuchterl03}) to determine the mass as a function of luminosity, temperature and estimated age. The resulting huge discrepancy shows the difficulty of determining an absolute mass for such directly imaged low mass objects and especially at such young ages.

As mentioned in \S\,\ref{Motiv_SDI}, the first surveys for brown dwarf binaries already revealed that the resolved systems showed a tendency toward a high mass ratio distribution (M$_{2}$/M$_{1}$\,=\,\emph{q} $\ge$ 0.8) with a clear peak at \emph{q}\,$\sim$\,1. Further high-resolution surveys on ultracool dwarfs searching for binary systems strengthened this trend (see Figure\,\ref{mass_ratio_dist} on the left). Although an underestimation of the number of low mass ratio systems cannot be excluded due to the lack of sensitivity for \emph{q} $\le$ 0.5 in most of these surveys, the clear peak to equal mass binary systems is not an observational bias effect (for a detailed review see \citealt{Burgasser07_3} and for an extensive Bayesian study \citealt{Allen07}). The same surveys unfolded a semi-major axis distribution which peaks around 3\,-\,5\,AU with only very few separations beyond $\sim$ 15 AU (e.\,g.\,\citealt{Bouy03}; Table\,1 in \citealt{Burgasser07_3}; supported by numerical simulations of \citealt{Umbreit05} and \citealt{Allen07}).

In the following, the results of our HST survey will be compared with respect to the findings discussed above. The resolved and confirmed brown dwarf binary Kelu-1\,AB with its mass ratio \emph{q} = 0.82 and a separation of 5.59\,AU at the time of discovery, as well as 2MASSW\,0310+1648\,AB with its mass ratio \emph{q} $\sim$ 1 and a separation of 5.17 AU, are in full agreement with the previous findings and are an important contribution to the brown dwarf binary statistics. The non-detection of any companion in this survey at separations larger than 10 AU is also consistent with the separation statistics. The general observed binary fraction of very low mass stars (VLMS) and brown dwarfs, as found in direct imaging surveys is $\sim$\,15\,-\,20\%, and might be as high as $\sim$\,25\,-\,30\% (for extended discussions see e.\,g. \citealt{Burgasser07_3, Allen07}). With a total sample of 13 brown dwarfs, the two confirmed binary systems yield a binary fraction of $15^{+18}_{-10}$\,\% (or $23^{+19}_{-13}$\,\% if 2M0337-1758\,AB is confirmed as a third binary system\,\footnote{The uncertainty of the measured binary fraction was estimated by searching for the fraction of binaries at which the binomial distribution P$_B$(x,n,p) for x positive events in n trials falls off to 1/$e$ of its maximum value (cf. \citealt{Basri06}).}) for our survey. This is consistent with the results from previous VLMS/BD surveys.

The discussion about the companion candidate to 2M0337-1758 can only be preliminary until second epoch observations confirm the companionship. Yet, it shows some very interesting properties. With a preliminary mass ratio of 0.18\,-\,0.22 (depending on its estimated age), this system would add to the very small number of resolved brown dwarf binaries with \emph{q} $\le$ 0.5. One of them, which has a resolved and confirmed  brown dwarf primary and a similar mass ratio is 2M1207\,AB (\emph{q} $\approx$ 0.3 if the most actual mass estimation of \citet{Mohanty07} is taken into account). But even if its secondary has a mass just below the deuterium burning limit, and therefore appears to be similar to the possible 2M0337-1758\,AB system, the \object[NAME 2M1207b]{2M1207\,AB} system has an unusual large separation and is much younger ($\sim$\,8 Myr). At these very young ages the surveys suggest a flatter mass ratio distribution than that for field or older cluster binaries (see shaded bins in Figure \ref{mass_ratio_dist} on the left), thus the low \emph{q} value might not be too unusual in the case of 2M1207\,AB. Another interesting system in this account is \object[SCR J1845-6357B]{SCR\,1845\,AB}. This M8.5\,-\,T6 field binary (where the secondary was discovered during a SDI survey as well) reveals also a  relatively low mass ratio with \emph{q} $\approx$ 0.41 (\citealt{Kasper07}). Additionally, the \object[SCR J1845-6357]{SCR\,1845\,AB} system is, with a separation of $\sim$\,4.5 AU, similar tightly bound as 2M0337-1758\,AB. Hence, both systems might be potential candidates for a dynamical mass determination on a reasonable time scale. This will further help to calibrate evolutionary models under a different aspect. The mass ratio distribution in binary systems exhibits a clear correlation with the primary mass (see Figure \ref{mass_ratio_dist} on the right). While G dwarf primaries show a relatively flat mass ratio distribution, the ratio for M dwarfs, VLMS and brown dwarfs peaks more towards unity (see also \citealt{Bergfors10} and references therein). Although a \emph{q} value of 0.18\,-\,0.22 is, independent of age and separation, still much larger than the mass ratios known from general star\,-\,giant planet systems (determined with RV, transit etc.\,methods), it does not completely exclude a brown dwarf\,-\,PMO system. Due to the lower initial mass of the primary the mass ratio of such a system would be a priori higher than for a main sequence star\,-\,giant planet system. Thus, such a \emph{q} value for our system probes an interesting ratio regime which thus far has not been explored for brown dwarfs.

The minimum mass estimate of 10\,-\,15 $M_\mathrm{Jup}$ for the companion candidate of 2M0337-1758 does not a priori exclude a classification as a planet. Nevertheless, a comparison of four of the $\sim$\,600\,K  T dwarfs (all field dwarfs) in \citet{Leggett09} shows, that at least three of them (so far single objects) have similar mass estimates, just right below or above 13 $M_\mathrm{Jup}$. Therefore, 2M\,0337-1758 is more likely a brown dwarf -- brown dwarf binary system with a so far unusual very low mass secondary which possibly straddles the T/Y transition, rather than a system including an exoplanet.
Follow-up observations of the system are proposed to achieve a definite classification. Even more important, a larger number of further detections of objects in this mass and temperature regime are necessary to reliably calibrate the theoretical models.

\section{Conclusions}
While previous direct imaging surveys focused especially on searches for planetary mass objects around young stars in nearby associations, we obtained high-resolution observations with HST/NICMOS of 12 L dwarfs within 30 pc of the sun, using the spectral differential imaging technique in the two narrow-band filters F108N and F113N. To this date, this is the only SDI survey for planetary mass companions around field brown dwarfs (not related to any young association or moving group) conducted from space or ground-based telescopes. Further, it is the only one which uses the spectroscopic feature of water absorption instead of the methane feature.

The survey resolved two brown dwarf binaries (Kelu-1\,AB and 2MASS\,0310+1648AB) for which common proper motion has already been confirmed in follow-up observations (\citealp{Stumpf09, Stumpf10_1}). The overall binary fraction of $15^{+18}_{-10}$\,\%, as well as separations of the binaries smaller than 6 AU are consistent with previous brown dwarf binary statistics. The mass ratios of \emph{q} $\ge$\,0.8 confirm the preference for equal mass systems as found in previous direct imaging surveys.

Furthermore, tentative evidence was found for a very low mass companion close to the planetary mass regime around 2MASS\,0337-1758. It has a very low T$_{\mathrm{eff}}$ $\approx$ 600\,-\,630 K and an unusual mass ratio (\emph{q} $\approx$ 0.2) compared to the vast majority of previously found brown dwarf binaries. Follow-up observations are proposed to confirm common proper motion of the system and to derive more physical properties like colors for spectral type classification. In addition, these new observations will provide second epoch astrometry for the determination of orbital parameters and hence an initial dynamical mass estimation. If gravitationally bound and depending on the age, the 2MASSW\,0337-1758 system will be one of the closest resolved brown dwarf binaries (0.087$\arcsec \pm$\,0.015$\arcsec$, corresponding to 2.52\,$\pm$\,0.44\,AU at a distance of 29\,pc) with the so far coldest companion ever directly imaged and the least massive companion to any L or T dwarf.
Therefore, including its multiplicity status, it will then be an important testbed in the newly explored $\le$ 700 K temperature regime and might imply new constraints on the existing formation scenarios.

\acknowledgments

M.B.Stumpf and W. Brandner acknowledge support by the \emph{DLR Verbundforschung} project numbers 50 OR 0401 and 50 OR 0902.
We are grateful to Tricia Royle at STScI for her great and efficient support before and during observations. We would like to thank the anonymous referee for very helpful comments to improve the paper.
This research has benefitted from the M, L and T dwarf compendium housed at DwarfArchives.org and maintained by Chris Gelino, Davy Kirkpatrick and Adam Burgasser. This publication made also use of the Very-Low-Mass Binaries Archive housed at http://www.vlmbinaries.org and maintained by Nick Siegler, Chris Gelino and Adam Burgasser; and it made use of the SIMBAD database, operated at CDS, Strasbourg, France.

{\it Facilities:}  \facility{HST (NICMOS)}.


\section{Introduction

While magnetic fields are known to play a prominent role for the
dynamics and in the energy budget of astrophysical systems on
galactic and smaller scales, their role on larger scales is still
elusive~\citep{r02,r08}.
Extragalactic magnetic fields (EGMF) are notoriously difficult to
measure and the data are incomplete. So far only in a few galaxy
clusters observational constraints have been obtained, either by
observing their synchrotron radiation halos or by performing Faraday
rotation measurements (RMs). Within galaxy clusters the inferred magnetic
fields are between 0.1 and 1.0\,$\mu$G on scales as large as 1 Mpc,
and can be as strong as 30\,$\mu$G localized inside cluster cool cores.
However, as both observational methods need a prerequisite to
measure magnetic fields (high thermal density for RMs and presence of
relativistic particles for radio emission), they have been successfully
applied only to high density regions of collapsed objects as galaxies and
galaxy clusters. Fields significantly below $\mu$G level are barely
detectable with these methods. Also other constraints, for instance the
absence of distortions  in the spectrum and the polarization properties
of the cosmic microwave background  radiation implies only a fairly
large, global upper limit on the EGMF at the level of  $10^{-9}$\,G.

The observed magnetic fields in galaxies and galaxy clusters are
assumed to result from the amplification of much weaker seed
fields. Such seeds could be created in the early universe,
e.g.\ during phase transitions, and then amplified by plasma processes.
Alternatively, an early population of starburst galaxies or AGN could
have generated the seeds of the EGMFs at redshift between five
and six, before galaxy clusters formed as gravitationally
bound systems.
In both cases, a large fraction of the material collapsing to today's
visible structures will be seeded by such fields, but whether a significant
fraction of the volume of the universe can be filled by such EGMFs depends
on the efficiency of the transport processes involved.
A quite different possibility is that the ejecta of AGN magnetized the
intracluster medium only at low redshifts, and that thus the EGMF are confined
within galaxy clusters and groups.
The detection of EGMFs outside clusters is therefore crucial in
discriminating different models for the origin of their seed fields.

An alternative approach to obtain information about the EGMFs is
to use its effect on the radiation from TeV gamma-ray sources.
The multi-TeV $\gamma$-ray\ flux from distant blazars is strongly attenuated
by pair production on the infrared/optical extragalactic
background light (EBL),  initiating electromagnetic cascades in the
intergalactic space. The charged component of these cascades is
deflected by the EGMF. Potentially observable effects of such
electromagnetic cascades in the EGMF include the delayed ``echoes" of
multi-TeV $\gamma$-ray\  flares or gamma-ray bursts~\citep{plaga,jap} and the
appearance of extended
emission around initially point-like $\gamma$-ray\ sources
\citep{coppi,neronov07,halo1,el09,jet}.

An additional way to derive lower limits on the EGMF has been
pointed out recently by \citet{ne10} and \citet{ta10}: Since the
deflection of the cascade flux into an extended halo weakens the
point-like image, the non-observation of TeV blazars in the GeV range
by Fermi-LAT can been used to derive a lower limit on the
EGMF. Particular suitable candidates are blazars with a very
hard TeV spectrum like 1ES 0229+200\ that show a low
intrinsic GeV emission. In this way, \citet{ne10} and \citet{ta10}
derived the lower bound $B\raise0.3ex\hbox{$\;>$\kern-0.75em\raise-1.1ex\hbox{$\sim\;$}} 5\times 10^{-15}$\,G on the EGMF.
Note that \citet{ak} found evidence for gamma-ray halos in the
stacked images of the 170 brightest active galactic nuclei
observed with Fermi-LAT. The size of these halos are consistent with
$B\sim 10^{-15}$\,G, however it is unclear if the excess found is real
or due to the imperfect knowledge of the Fermi-LAT point spread function
(PSF) as argued by \citet{Neronov:2010bi}.
All these analyses assumed a stationary source and did not discuss
the effect of time delays induced by EGMFs.

In this work, we will concentrate on 1ES 0229+200\ and follow
essentially the same assumptions about the source and the sensitivity
of Fermi-LAT as \citet{ta10}. We will improve on previous analyses in
two respects: First, we use a Monte Carlo simulation for the
development of electromagnetic cascades in the EBL that includes the
effects of magnetic fields like synchrotron radiation and deflections
of electrons. Moreover, we calculate the time delay that is
induced by the EGMF deflections.
Second, we examine the influence of a more realistic,
structured magnetic field on the EGMF limit.

\section{Simulation procedure

We inject at the redshift $z=0.14$ photons with an energy
distribution given by ${\cal F}\propto E^{-2/3}$ and a maximal
energy $E_{\max}=20$\,TeV. Such a hard injection spectrum gives a
reasonable fit to the HESS data and is consistent with recent
SWIFT~\citep{swift} observations as argued by \citet{ta09}.
We assume also  a relatively low Lorentz factor, $\Gamma=10$,
i.e.\  $\Theta_{\rm jet}=6^\circ$, in order to derive a conservative
lower bound on the EGMF.

We follow electromagnetic cascades initiated by these high energy photons with
the Monte Carlo code introduced by \citet{CenA}. Since Fermi-LAT provides
only an upper limit on the photon flux from 1ES 0229+200\, it is sufficient to
assume a radial-symmetric halo and to calculate the deflections
and time delays
as described in~\citep{halo1}. We use as EBL background the best-fit
model from the calculations of \citet{kneiske10}.

We divide the photons arriving at the observer into a point-like
flux and a halo outside a circle with angle $\vartheta_{\ast}$ around the
source. For Fermi-LAT, we use for $\vartheta_{\ast}$ an analytical
approximation of the angle $\vartheta_{95}$ containing 95\% of the
flux emitted by a point-source
given by
$ \vartheta_{95}\approx 1.68^\circ (E/{\rm GeV})^{-0.77} + 0.2^\circ
 \exp(-10\,{\rm GeV}/E)$.
Above 300 GeV, $\vartheta_{95}$ is set equal to
$0.11^\circ$, a value typical for the resolution of atmospheric
Cherenkov telescopes.

\section{Results for different EGMF models

In Figs.~\ref{fig:uni} and \ref{fig:fill} we show our results for the
fluence contained inside the 95\% confidence contour of the PSF of
Fermi-LAT. Additionally, these figures contain the HESS observations
\citep{HESS} as black dots with error bars and the Fermi-LAT upper
limits derived by \citet{ta10}. The fluences have been
normalized fitting them to the HESS data.

\begin{figure}
\vskip-2.2cm
\includegraphics[width=\linewidth,angle=0]{fig1}
\caption{Fluence contained inside the 95\% confidence contour of the
PSF of  Fermi-LAT as function of energy together with Fermi-LAT upper
limits and HESS observations for a uniform EGMF with strengths
varying from $B=10^{-16}$\,G  to $B=10^{-13}$\,G with
 $E_{\rm max}=20$\,TeV (solid) and 100 TeV (dashed). The direct
  component for $B=10^{-14}$\,G is also shown.}
\label{fig:uni}
\end{figure}

\begin{figure}
\vskip-2.2cm
\includegraphics[width=\linewidth,angle=0]{fig2}
\caption{
Fluence contained inside the 95\% confidence contour of the
PSF of  Fermi-LAT as function of energy for a EGMF with top-hat
profile and filling factor $f$ varying from $f=0.1$ to $f=0.9$
with $E_{\rm max}=20$\,TeV (solid) and 100 TeV (dashed).}
\label{fig:fill}
\end{figure}

In Fig.~\ref{fig:uni}, we have used a uniform magnetic field so that
our results can be directly compared with the analytical estimates of
\citet{ta10}.  Note that a turbulent field with correlation length
$L_{\rm cr}$ much larger than the mean free path $l_{\rm IC}$ of
electrons in the Thomson regime, $L_{\rm cr}\gg l_{\rm IC}\sim
1$\,kpc, is well approximated by a
uniform field. For smaller correlation lengths, $L_{\rm cr}\ll l_{\rm
  IC}$, the electron diffuses in the small-angle deflection regime,
requiring larger magnetic fields for the same deflection angle.
Demanding that the cascade flux is below the upper limits of Fermi-LAT
leads to a lower limit on the magnetic field strength of
$\sim 10^{-14}$\,G. For this case the direct component, i.e.\
photons arriving at the detector without cascading, is also
  shown.
Note that for small $E_{\max}$ the transition from the direct to the
cascade contribution leads  to a break at $\sim$\,TeV in the
spectrum, as suggested by the HESS data.

While our limit agrees
reasonably well with the analytical estimate of \citet{ta09}, the
shape of the cascade flux obtained differs. There are several reasons
responsible for these differences: First, \citet{ta09} assume that the
spectral shape of the cascade flux below the threshold energy $\sim
10^{11}$\,eV is given by ${\cal F}\propto E^{-0.5}$ for negligible magnetic
fields. Such a slope typical for the regime of Thomson cooling is
however restricted to energies $E\raise0.3ex\hbox{$\;<$\kern-0.75em\raise-1.1ex\hbox{$\sim\;$}} 10^8$\,eV, while at higher
energies a plateau ${\cal F}\propto E^\alpha$ with $\alpha\sim -0.9$ is
expected~\citep{casc}. Second, deflections in the EGMF lead even for
isotropic emission to extended images of point-like sources.  This
effect has been neglected in \citet{ta09}.  Third, using full
probability distributions for the interactions there is a
non-negligible probability for photons not to
interact, especially towards the low-energy end of the
  injected energy range.  Finally, the energy dependent PSF of Fermi
introduces an artificial energy dependence of the point-like flux
${\cal F}\propto \vartheta_{95}^2$.
Note that an increase in $E_{\rm max}$ from 20 to 100 TeV {\em reduces\/}
the limit on the magnetic field strength to $\sim 5\times 10^{-15}$\,G,
see Fig.~1, while a further increase of $E_{\rm max}$ strengthens the limit
again. The counter-intuitive behavior between 20 to 100\,TeV is caused
by the dominance of direct photons at the high-energy part of the
spectrum for small $E_{\max}$.

We discuss next the consequence of time delays induced by the
  EGMF.  We stress that the non-observation of 1ES 0229+200\ by Fermi-LAT is a
  signature for non-zero EGMFs, even if the delay of the Fermi signal
  would exceed the life time $\tau$ of 1ES 0229+200\ as TeV source. However,
  the numerical value of the limit deduced using deflections would be
  modified. For $B=5\times 10^{-15}$\,G, we find as average time delay
  $1\times 10^6$, $1\times 10^5$, and 5000\,yr for the three energy
  bins (0.1--10), (1--10), and (10--30)\,GeV. Thus the high-energy bin
  of the Fermi data can be used to derive the constraint as above, if
  1ES 0229+200\ is a relatively stable source over a time scale longer than
  $10^4$\,yr.  Even if $\tau$ is much smaller, say around 100\,yr, the
  lower limit on $B$ weakens only by a factor 10 to $B={\cal
    O}(10^{-16}\,{\rm G})$.

Since the EGMF is strongly structured, one may wonder how  a non-uniform
field modifies this limit. In particular, we want to
address the question whether the presence of relatively strong fields
concentrated inside cosmic structures like filaments could mimic the
effect of an EGMF present also in voids. As simplest possible test,
we use first a top-hat profile for the structure of the EGMF:
We set the field strength to zero in a fraction $1-f$ of space and use a value
which in general is assumed to be representative for filaments,
$B=10^{-10}$\,G, in the remaining part. For the
separation of the peaks we use $D=10$\,Mpc motivated by the typical
distances between cosmological structures,
although the exact value of $D$ plays
no role as long as $(1-f)D\ll l_\gamma$, with  $l_\gamma$ as
mean free path of photons.
The dependence of the fluence contained inside the PSF of Fermi-LAT
on the filling factor $f$ is shown in Fig.~\ref{fig:fill}. To be
consistent with the Fermi upper limits, sufficiently strong magnetic
fields should fill $\raise0.3ex\hbox{$\;>$\kern-0.75em\raise-1.1ex\hbox{$\sim\;$}} 80\%$ of space. The derived limit on the
filling factor is practically independent of the source life time $\tau$ and
$B$, as long as the
field is stronger than $\gtrsim 5\times 10^{-15}$\,G.
As in the previous case, by assuming a higher injected $E_{\rm max}$ the
required filling factor is slightly reduced to 60\% .

The failure of strong fields filling only a small fraction of the universe
to suppress sufficiently the point-like cascade flux can be understood
as follows: The HESS observations of 1ES 0229+200\ cover the energy range
0.5--11\,TeV. In the same energy range, the mean free path $l_\gamma$
of VHE $\gamma$-ray s through the EBL varies between 1000 and 50\,Mpc
and is thus always much larger than the typical extension of
regions with large fields, $(1-f)D$. For the energies considered,
the cascade consists typically of only three steps,
$\gamma\to e^\pm\to\gamma$. Since the mean free path $l_{\rm IC}$
of electrons in the Thomson regime is very small, $l_{\rm IC}\sim 1$\,kpc,
all cascades with electrons created outside the strong-field regions are
undeflected. Thus it is not possible to trade smaller values of
$f$ against larger values of $B$: Increasing the field strength
beyond $\sim 10^{-13}$\,G leads only to an increase of the deflection,
while the fraction of cascades deflected outside the Fermi PSF
remains constant.

To mimic the even more complex situation in presence of cosmological
structures, we have used magnetic fields derived from several
cosmological MHD simulations using constrained initial conditions
(Dolag et al. 2004, Donnert et al. 2009; for alternative simulations
see e.g.\ Ryu et al. 2008). Therefore they reproduce fairly
well the true large scale matter distribution, out to a distance of
$\approx 114$\,Mpc, cf.\  for more details Dolag et al. (2005).
We extracted the EGMF
as predicted by these simulations along the line of sight towards the
position of 1ES 0229+200 for several models for the magnetic seed field.
A comparison
of the EGMF component perpendicular to the line-of-sight
towards 1ES 0229+200\ obtained in the models is shown
in Fig.~\ref{fig:Bprofiles}.

\begin{figure}
\vskip-2.2cm
\includegraphics[width=\linewidth,angle=0]{fig3}
\caption{The EGMF component perpendicular to the line-of-sight
towards 1ES 0229+200\ from four different MHD simulations.}
\label{fig:Bprofiles}
\end{figure}

\begin{figure}
\vskip-2.2cm
\includegraphics[width=0.9\linewidth,angle=0]{fig4bis}
\caption{Cumulative volume filling factor $C(B)$
for the four different EGMF models found in MHD simulations.
\label{fig:Bcumulative}}
\end{figure}

\begin{figure}
\vskip-2.2cm
\includegraphics[width=\linewidth,angle=0]{fig5}
\caption{Fluence contained inside the 95\% confidence contour of the
PSF of  Fermi-LAT as function of energy for EGMFs from
four different MHD simulations  with
  $E_{\rm max}=20$\,TeV (solid) and 100 TeV (dashed).}
\label{fig:Bresult}
\end{figure}

Most magnetic seed models used in cosmological simulations
assumed a space filling (e.g.\ primordial) magnetic seed field of order
$10^{-12}$G, which inside galaxy clusters is amplified by the collapse
and additionally by shear turbulence up to the observed $\mu$G level. For the
purpose of this work
we have selected the available simulation with the lowest primordial seed
field of $10^{-13}$\,G (MHDx), which inside voids is even diluted to values
of order of $10^{-14}$\,G due to cosmic expansion. To construct a second model
with magnetic fields as low as $10^{-17}$\,G inside voids,
we scaled down the magnetic field of the original simulation by a factor
proportional to the local density to the power 2/3 (Model 3x) with the
fix point at the core density of galaxy clusters. Thereby, the Faraday
RMs of the newly constructed model within galaxy
clusters does not significantly change, whereas the magnetic field in
filaments and voids drop by nearly three orders of magnitude compared
to the original simulation.

Models which intrinsically produce lower magnetic fields (and filling
factors) in filaments and voids are those where the seed
field originates from the galaxies formed. In
such simulations, galaxies are identified at early times, and their
halos are magnetized according to the magnetic field observed in
nearby galaxies.  So far such simulations do not include any transport
processes for the magnetic field (like galactic outflows or AGN
activity) beyond the large scale potential well in which galaxies
live.
Nevertheless, these models reproduce well the
  observed magnetic field properties within galaxy clusters, since
  the gas component of the galactic halos is stripped
  during the formation process of the cluster and the magnetic field
  gets well mixed and subsequently amplified by the turbulence within
  the intercluster medium.
The result of such a simulation is the
model (MHD Gal3) from Donnert et al.\ 2009. As
galaxies in filaments form later than their counterparts in high
density environment, we have also extracted the magnetic field from a
model which subsequently put magnetic fields into newly formed halos
also at lower redshift (MHD Gal5). The main difference to the previous
model is the presence of more magnetized structures within filaments
whereas the magnetic field properties in galaxy clusters do not change
significantly.

The simulation used covers only up to distances of
114\,Mpc, and therefore we are forced to extend the magnetic field
profile beyond that point. For this purpose we assume
that the extended field has the same statistical properties
as the simulated patch.
In practice, we simply mirror and repeat the field beyond the
first 114\,Mpc.
In Fig.~\ref{fig:Bcumulative}  we show
  the cumulative volume filling factor, i.e.\ the fraction $C(B)$ of
  volume occupied by regions with field strength $>B$. The vertical
  and horizontal bands stand for the minimum magnetic field and
  filling factor, respectively, required to avoid the Fermi limit, as
  obtained in Figs.\,\ref{fig:uni} and \ref{fig:fill}.

The fluence contained inside the PSF of Fermi-LAT for the EGMF models
MHDx, MHD Gal3 and MHD Gal5 is shown in Fig.~\ref{fig:Bresult}.  From
the results for the top-hat model, we anticipate that the EGMF in the
models MHD3x and MHD Gal3 is too weak and does not lead to a
significant modification of the fluence, see
  Fig.~\ref{fig:Bcumulative}:
Since the model
  MHD Gal5 is not far from the intersection of the two bands,
  we expect it not to lie far from the Fermi limit. In
  Fig.~\ref{fig:Bresult} it is seen how in the particular case of
  $E_{\rm max}=100$\,TeV this model is marginally consistent with the Fermi-LAT
  observations.

Finally, the fluence is several orders of magnitude below the
observational limits for fields as strong as in the model MHDx.
In such a case, extended images in the GeV range as predicted e.g.\
by \citet{jet} will be impossible to observe with Fermi-LAT.

We would like to point out, that on the other hand, all these models
are in agreement with observations of magnetic fields in galaxy clusters --
which
so far is the only place where the strength and structure of magnetic fields
has been directly detected -- and therefore demonstrate the constraining
power of using high energy photons to infer the origin and the processes
leading tho EGMFs.

\section{Summary

We have calculated the fluence of 1ES 0229+200\ as seen by Fermi-LAT using a
Monte Carlo simulation for the cascade development.  We have discussed
the effect of different EGMF profiles on the resulting suppression of
the point-like flux seen by Fermi-LAT.  Since the electron cooling
length is much smaller than the mean free path of the TeV photons, a
sufficient suppression of the point-like flux requires that the EGMF
fills a large fraction along the line-of-sight towards 1ES 0229+200,
$f\raise0.3ex\hbox{$\;>$\kern-0.75em\raise-1.1ex\hbox{$\sim\;$}} 0.6$. The lower limit on the magnetic field
strength in this volume is $B\sim {\cal O}(10^{-15})$\,G,
assuming 1ES 0229+200\ is stable at least for $10^4$\,yr, weakening by a factor of
10 for $\tau=10^2$\,yr.
These limits put very stringent constraints on the
origin of EGMFs: Either the seeds for EGMFs have to be produced by a
volume filling process (e.g.\ primordial) or very efficient transport
processes have to be present which redistribute magnetic fields that
were generated locally (e.g.\ in galaxies) into filaments and voids
with a significant volume filling factor.

\vskip0.2cm \textit{Acknowledgments.}  We are grateful to L. Costamante and
D. Paneque for discussions and to the referee
for his/her useful remarks. K.D.\ is supported partly by
the DFG SFB 1177 and a Cluster of Excellence, while S.O.\ acknowledges
a fellowship from the program Romforskning of Norsk
Forsknigsradet.

\vskip0.2cm
\textit{Note added.}
After the completion of this work the preprint of Tavecchio et al.\ (2010b)
appeared in which these authors study the four BL Lac objects (RGB J0152+017,
1ES 0229+200, 1ES 0347-121 and PKS 0548-322) and derive a bound of
$B\sim 5\times 10^{-15}$\,G on a space filling magnetic field.