\section{Introduction}

Comet 17P/Holmes underwent the largest cometary outburst in recorded history on UT 2007 Oct.\ 23.7 $\pm$ 0.2 \citep{2007IAUC.8886....1B}, brightening from $\sim$ 17$^{th}$ magnitude to $\sim$ 2$^{nd}$ in 42 hours \citep{2008ApJ...677L..63M,2011ApJ...728...31L}.  The Jupiter family comet is known to have experienced at least two prior outbursts: the first led to its discovery on UT 1892 Nov.\ 6 \citep{1892Obs....15..441H} while the second occurred in January 1893 \citep{1893Obs....16..142D}.  \cite{1984Icar...60..522W} speculated that the unusual outburst activity might have been triggered by a grazing encounter with a natural satellite.  In this scenario, the second outburst would have been caused by the final impact of the satellite with the nucleus.  A third major outburst practically eliminates the likelihood of this scenario, suggesting that an alternative mechanism causes these unusual periods of activity in 17P/Holmes.  The two largest outbursts - those that occurred in 1892 Nov.\ and 2007 Oct.\ - both took place approximately five months after 17P/Holmes passed through perihelion at 2.1 AU, suggesting that the activity is thermally-triggered, perhaps with a lag due to the time needed for heat to reach subsurface ice deposits \citep{2008ICQ....30....3S}.

The mass of material ejected has been estimated to be between 10$^{10}$ kg and 10$^{12}$ kg \citep{2008ICQ....30....3S,2010Icar..208..276R,2011ApJ...728...31L}, corresponding to $\sim$ 0.05\% to 5\% of the nucleus mass.  This ejected mass was dominated by the contribution from the dust coma that was observed to expand around the nucleus.  The peak dust production rate approached 10$^{6}$ kg s$^{-1}$ on UT 2007 Oct.\ 24.5 $\pm$ 0.01 \citep{2011ApJ...728...31L}, about 10$^{3}$ times the dust production rates typical of Jupiter Family comets when at similar heliocentric distances.

In addition to massive quantities of dust, 17P/Holmes released significant quantities of volatiles with H$_{2}$O production rates exceeding 4.5 $\times$ 10$^{29}$ s$^{-1}$ \citep{2008ApJ...680..793D,2009AJ....138.1062S}.  Results from large-aperture (radii $>$ 10,000 km) narrowband photometry in the weeks following the outburst led \cite{2009AJ....138.1062S} to conclude that enhanced water production rates were primarily due to prolonged nucleus activity, with a possible contribution from sublimating icy grains in the inner coma.

In this work, we use a sequence of high-resolution wide-field continuum images to examine the evolution of 17P/Holmes during the four months following the outburst.  By utilizing a single instrument and observing strategy we minimize the effect of systematic errors and are able to use a uniform procedure for obtaining, reducing, and analyzing the data.  This results in a consistent, high-quality data set.  In contrast to an earlier study of the integrated light from the comet (Li et al.~2011), we here study material near the nucleus, which is the most-recently ejected and the most pristine of the observable ejecta.

\section{Observations and Data Reduction}

We used the 3.6 m Canada-France-Hawaii Telescope (CFHT) atop Mauna Kea with the wide-field imager MegaCam \citep{2003SPIE.4841...72B} to obtain high-quality, time-series imaging data.  MegaCam has a 0.96${^\circ}$ $\times$ 0.96${^\circ}$ field of view provided by an array of 36 charge-coupled devices (CCDs).  Each CCD has 2048 $\times$ 4068 pixels, with a pixel scale of $0^{\prime\prime}.185$ pixel$^{-1}$.  We obtained at least one set of five dithered images in the SDSS r$^{\prime}$ filter ($\lambda_{c} =$ 6250 \AA, bandwidth = 1210 \AA) on each of 15 nights between UT 2007 Nov.\ 6 and UT 2008 Feb.\ 12 UT.  Throughout 2007 Nov.\ and on UT 2007 Dec.\ 11 we obtained two sets of exposures with differing exposure times to allow us to probe the bright central region while still being sensitive to the faint outer coma.  In this work we use the set of exposures with the longest exposure time but no saturation of the central region.  We have analyzed all sets of the exposures and find the agreement between sets obtained on the same night to be very good, with differences in calculated magnitudes to be less than 0.02 mag.

The observations were taken in queue-scheduled mode.  In this mode, non-sidereal tracking was unavailable and the comet trailed across each image.  Exposure times varied between 5 s and 200 s, and were chosen to balance the need for sensitivity while minimizing trailing.   Table~(\ref{table:obs}) provides a journal of observations.

All images were automatically pre-processed using the Elixir pipeline \citep{2004PASP..116..449M}, which uses bias frames and twilight flatfields to remove the instrumental signature, and calculates instrumental zero-points using Landolt standard star fields \citep{1992AJ....104..340L} observed on each night.  We used SWarp, software released by the Terapix data center at the Institut d'Astrophysique de Paris, to resample the images and attain $0^{\prime\prime}.2$ astrometric accuracy.  The dithered images from each night were median-combined to cover the 80$^{\prime\prime}$ gaps between CCDs and remove cosmic rays (Figure~\ref{fig:fullFoV}).  Each combined image was convolved with a two-dimensional Gaussian function with a full-width half-maximum of 1$^{\prime\prime}$ to minimize the effect of variations in seeing from night to night.  A panel showing the central square arcminute of each image is shown in Figure~(\ref{fig:nucpan}).  The same central regions are shown with isophotal contours in Figure~(\ref{fig:nuccon}) to better highlight the evolution of the coma morphology.  The inner coma initially appears elongated in the anti-solar direction but becomes more symmetrical at later dates.

The sky transparency was monitored using the CFHT Skyprobe and  we found that only observations made on UT 2008 Feb.\ 12 experienced significant variations in attenuation ($>$ 0.1 mag).  For all other nights the weather was photometric.  Systematic errors due to flatfielding or scattered light are well-characterized for MegaCam and do not exceed 5\% \citep{2004PASP..116..449M}.  The sky background was determined in a region as far from the observable coma as possible within each stacked image, using the median signal within a $\sim$ 4$^{\prime} \times$ 4$^{\prime}$ box.  Although we sought regions on the image devoid of any noticeable brightness gradient from which to estimate the background level, it is likely that some coma contamination persisted in our sky background estimations.  Specifically, the gas coma was observed to expand at a velocity of 1 km s$^{-1}$ \citep{2008LPICo1405.8340D}, meaning that the gas coma reached the projected edges of the MegaCam field of view by 2007 mid-November.  Fine dust grains dynamically well-coupled to the gas and ejected at high velocities could contaminate the outer regions within a month of the outburst, although our observations are less sensitive to the smallest grains since their scattering efficiency is small.  In Figures~(\ref{fig:fadingtot}) and~(\ref{fig:fadingnov}) we show the effects of a 10\% error on the sky background combined with the known sky transparency and systematic errors.  Plausible sky subtraction uncertainties are negligible in most of our data, given the high surface brightness of the coma.  Only in observations from January and February 2008 do sky subtraction errors rival systematic uncertainies due to atmosphere transparency and calibration uncertainties.

\section{Photometry of the Inner Coma}

The coma was initially observed as an approximately circular dust cloud with a bright condensation south-west of the nucleus.  As it expanded, the coma became more asymmetric with a parabolic nose building in the solar direction and a diffuse tail in the anti-solar direction.  The nucleus region was always identifiable in the inner coma and stars could be seen even through the innermost coma, implying a low optical depth.  We used aperture photometry to study the material closest to the nucleus with the aim of monitoring the comet's continued activity after its massive outburst.

The position of the nucleus was determined to sub-pixel accuracy using a two-dimensional derivative search and used as the center for all aperture photometry.  Since the geocentric distance increased from 1.62 AU on UT 2007 Nov.\ 6 to 2.62 AU on UT 2008 Feb.\ 12 (Table~\ref{table:obs}), an aperture of fixed angular radius would cover an increasing physical distance on the sky, and, correspondingly, more of the coma.  For this reason, we adjusted the angular radii of the photometry apertures used for each night to compensate for the increase in geocentric distance to 17P/Holmes, and instead used apertures of fixed physical radius at the comet.  Circular apertures of projected radii between 250 km and 25000 km in increments of 250 km were used.

We chose not to subtract the flux of the nucleus from these measurements.  The expected flux contribution from the nucleus can be calculated using the observing geometry, a nucleus radius of 1.62 km, and a geometric albedo of 0.04 \citep{2006MNRAS.373.1590S,2004come.book..223L}.  The calculated nucleus contribution was typically less than 1\% of the flux measured within a 2500 km radius aperture.  In the worst case, the nucleus contribution amounted to 3\% of the measured flux but in all cases it was significantly smaller than other sources of error.

Trailing losses are negligible for observations between UT 2007 Nov.\ 6 and UT Dec.\ 11 since the motion of the comet during each exposure was less than the full-width-half-maximum of the image.  Later images were trailed by up to 7.7 pixels (1$^{\prime\prime}.4$).  To investigate the effect of trailing on our photometry we took an image with minimal trailing of 0.2 pixels (0$^{\prime\prime}$.4) and shifted and added duplicate images to mimic trailing.  We found that the loss of light from a 1.3$^{\prime\prime}$  radius aperture (2000 km for the image most badly affected by trailing) was on the order of 2\% and fell to less than 0.1\% for a 6$^{\prime\prime}$.3 (10$^{4}$ km) radius aperture.  We conclude that trailing has a negligible effect on the results from aperture photometry in our data.

We converted fluxes to apparent magnitudes, $m$, using the instrumental zero-points calculated by the Elixir pipeline, correcting for both airmass and the color of the filter.  To correct for changes in observing geometry we used Equation~(\ref{absmageq}) to determine absolute magnitudes, $H$ -
the magnitudes objects would have if observed at heliocentric and geocentric distances of 1 AU and a phase angle of 0$^{\circ}$.

\begin{equation}
H = m - 5\ log_{10} (r_{H} \Delta) - 2.5\ log_{10} \Phi_{\alpha}
\label{absmageq}
\end{equation}

\noindent where $r_{H}$ and $\Delta$ are the heliocentric and geocentric distances in AU respectively and $\Phi_{\alpha}$ is a phase function that corrects for the dependency of the brightness on the phase angle ($\alpha$).  The phase function of the observed dust coma is difficult to measure given the need for accurate photometry of the dust over a wide range of phase angles when the coma is bright enough to detect.  We use an approximation of the phase function for 1P/Halley based on a fit to data presented in \cite{1998Icar..132..397S} \citep{2011ApJ...728...31L}:

\begin{equation}
2.5\ log_{10} \Phi_{\alpha} = 0.045 \alpha - 0.0004 \alpha^{2}
\label{phieq}
\end{equation}

The phase corrections that result from Equation~(\ref{phieq}) match well with other commonly-used phase functions since the change in phase angle was relatively small ($<$ 8$^{\circ}$).  We acknowledge that this is a slight misuse of the concept of absolute magnitude given that Equation~(\ref{absmageq}) refers to unresolved point sources rather than distributed sources such as cometary comae.  An angular aperture will subtend a different distance at 1 AU than at other heliocentric distancs, resulting in an altered amount of material within the aperture.  However, since we used apertures of fixed physical distances from the nucleus we expect this effect to have little impact on our results.  Table~(\ref{table:char}) shows the results from aperture photometry for each night during our observing campaign.

The surface brightness was calculated from the total flux within concentric annuli centered on the nucleus and divided by the area measured in square kilometers.  The flux per unit area was converted to absolute magnitudes per million square kilometers at the distance of the comet, which accounted for changes in the observing geometry during the course of our observing campaign.  Figure~(\ref{fig:SBall}) shows the surface brightness profiles after the images were convolved with a 2-dimensional $1^{\prime\prime}$ Gaussian function.

\section{Results and Discussion}

\subsection{The Fading of Comet 17P/Holmes}
\label{sec:fading}
The absolute magnitude within a 2500 km radius aperture increases from 8.1 mag on UT 2007 Nov.\ 6 to 12.4 mag on UT 2008 Feb.\ 12, indicating a $\sim$ 50-fold decrease in brightness (Figure~\ref{fig:fadingtot}).  Both the initial fading observed throughout 2007 November and the fading during the entire observing campaign can be reasonably well fitted with an exponential function of the form
$H = a + b e^{c t}$
where $H$ is absolute magnitude, $t$ is time since the first observation in days, and $a$, $b$, and $c$ are constants.  The best fits, as determined via a non-linear least-squares fit, are shown in Figures (\ref{fig:fadingtot}) and~(\ref{fig:fadingnov}).  We note that there is no clear physical reason for the decrease in brightness to be fitted with an exponential function in log-linear space.  However, \cite{2009AJ....138.1062S} found that similar functions matched the changing production rates for OH, CN, NH, C$_{2}$, C$_{3}$, and the continuum over a period of several months.

The expansion of the dust coma was studied by numerous observers and the expansion velocity was found to be $\sim$ 550 m s$^{-1}$ \citep{2007CBET.1123....1G,2008ICQ....30....3S,2008LPICo1405.8278S,2010MNRAS.407.1784H}, although \cite{2008A&A...479L..45M} reported a lower value of 200 m s$^{-1}$.  For this paper we use the relationship between grain velocity, $v$ (m s$^{-1}$), and radius, $a_{\mu}$ ($\mu$m), derived by \cite{1951ApJ...113..464W} and assume that a 1 $\mu$m radius dust grain has a velocity of 550 m s$^{-1}$ so that:

\begin{equation}
v = \frac{550}{\sqrt{a_{\mu}}}.
\label{veleq}
\end{equation}

The high measured velocities imply that the majority of dust particles ejected would exit even a 25000 km radius aperture in $\sim$ 1 day.  Even 1 cm sized particles ($a_{\mu}$ = 10$^4$) would leave a 25000 km radius aperture in about 50 days by Equation~(\ref{veleq}).  Given the brightness observed for several months following the outburst, the nucleus must have continued to supply the inner coma with fresh material after the 2007 Oct.\ outburst.  This result initially seems to be in conflict with the production curve calculated by \cite{2011ApJ...728...31L}, who found that significant activity persisted for only a few days, after which relatively little new material was injected into the coma.  However, their analysis refers to the entire coma of 17P/Holmes, which spanned several hundred thousand km, whereas our work focuses on a very small fraction of the coma close to the nucleus, and is consequently much more sensitive to small, localized changes.

The apparent magnitude, $m$, can be used to estimate the total cross-section of the dust within an aperture using:

\begin{equation}
\sigma = \frac{2.25\times10^{22}}{g_{\lambda}  \Phi_{\alpha}} ~ \pi~ r_{H}^{2} ~\Delta^{2} ~ 10^{-0.4 (m - m_{\Sun})}
\label{crosssectioneq}
\end{equation}

where $\sigma$ is the total geometric cross-section of the material (m$^{2}$), $g_{\lambda}$ is the geometric albedo, and $m_{\sun}$ is the apparent magnitude of the Sun.  All other variables have the same definitions as in previous equations.  The variables $r_{H}$ and $\Delta$ are given in AU.  The albedo of the ejected material was measured to be between 0.03 $\pm$ 0.01 and 0.12 $\pm$ 0.04  at a phase angle of 16$^{\circ}$ \citep{2010ApJ...714.1324I}.  The corresponding geometric albedos - as would be observed at a phase angle of 0$^{\circ}$ - are estimated to be 0.05 and 0.20 after using Equation~(\ref{phieq}) to calculate a phase correction.

Given the uncertainties, we adopt a nominal value of 0.1 for the geometric albedo of dust grains and accept that the actual value may differ from this by a factor of two.  The apparent magnitude of the Sun in the SDSS r$^{\prime}$ filter is -26.95 mag \citep{2001AJ....122.2749I}.  We then convert the cross-section into a crude mass, assuming that spherical dust grains with radii of 1 $\mu$m and a bulk density of 1000 kg m$^{3}$ are scattering the light.  Figure~(\ref{fig:sigmass}) shows the change in cross-section and mass of the material within a 2500 km radius aperture between UT 2007 Nov.\ 6 and UT 2008 Feb.\ 12.  The cross-section of material contained within the aperture decreases from 6.6 $\times$ 10$^{9}$ m$^{2}$ to 1.3 $\times$ 10$^{8}$ m$^{2}$.  These cross-sections correspond to effective circular radii of 46 km and 6 km and masses of 8.9 $\times$ 10$^{6}$ kg and 1.7 $\times$ 10$^{5}$ kg, respectively, which was only a tiny fraction of the total mass ejected (2 - 90 $\times$ 10$^{10}$ kg; \citealt{2011ApJ...728...31L}).  The mass estimates are strongly dependent on the assumed size distribution of the particles, and are used here for order of magnitude estimates only.

For small apertures we can derive mass loss rates using the aperture crossing times of the dust grains.  Large apertures will contain a mix of old, slow-moving large grains and fresh, high-velocity small grains, while the reflected light contained in small apertures is dominated by new material that may be only hours old.
We continue to use our earlier assumptions that micron-sized particles traveling at 550 m s$^{-1}$ dominate the cross-section of observed material.  We infer a mass loss rate of 4900 kg s$^{-1}$ at the beginning of our observations, which had decreased to 90 kg s$^{-1}$ by UT 2008 Feb.\ 12.  Peak mass loss rates of 3 $\times$ 10$^{5}$ kg s$^{-1}$ and 6-7 $\times$ 10$^{5}$ kg s$^{-1}$ were estimated by \cite{2011ApJ...728...31L} and \cite{2008ICQ....30....3S}, respectively.  Our results are lower by several orders of magnitude, indicating a rapid decrease following the initial outburst.

\subsection{A Second Outburst}

The fading inner coma of 17P/Holmes brightened unexpectedly on UT 2007 Nov.\ 12 by $\sim$ 25\% (Figure~\ref{fig:fadingnov}).  The expected absolute magnitude, as determined from a fading function fit to the November data of the exponential form discussed previously, is 9.3 mag, whereas the absolute magnitude calculated from the observed brightness was 9.0 mag.  Given the associated photometric uncertainties, this is a 3.0 $\sigma$ (99.8\% probability) result.  We note that this significance increases to 3.7 $\sigma$ (99.97\% probability) if the UT 2007 Nov.\ 12 data point is excluded from the fit.
A background object bright enough to cause the increase would need to have an apparent magnitude of 14.6 mag; an inspection of the STScI Digitized Sky Survey showed that there was no such object near 17P/Holmes.  We attribute the brightening to a second, albeit much smaller, outburst on the surface of the nucleus.  The subsequent observation of 17P/Holmes on 2007 Nov.\ 13 also lies above the fit shown in Figure~(\ref{fig:fadingnov}), suggesting that some of the material released during the second outburst remained in the aperture at this time.

The surface brightness profiles of UT 2007 Nov.\ 11, 12 and 13 provide independent evidence for a second outburst (Figure~\ref{fig:SBall}).  The coma within $\sim$ 5000 km of the nucleus is brighter than expected on UT 2007 Nov.\ 12 when compared to the surface brightness profiles of 17P/Holmes on the previous and following nights.  We can estimate the onset time of the second outburst if we assume that material ejected was dominated by dust grains traveling at $\sim$ 550 m s$^{-1}$.  We find that this material could have traveled 5000 km in 2.5 hours, leading us to conclude that the outburst started at approximately 2007 Nov.\ 12.4.  It seems likely that the outburst shut down within hours, given that subsequent observations revealed no unusual activity.

The difference between the expected magnitude and the observed magnitude, as shown in Figure~(\ref{fig:fadingnov}) using a 2500 km radius aperture, amounts to 0.3 mag.  Using Equation~(\ref{crosssectioneq}) and assuming that the material consists of dust grains with radii of 1 $\mu$m, albedos of 0.1, and densities of 1000 kg m$^{-3}$, the extra material has a scattering cross-section of 6 $\times$ 10$^{8}$ m$^{2}$, and a mass of 8 $\times$ 10$^{5}$ kg, corresponding to a sphere of material with a radius of 6 m.

The short duration of the Nov.\ 12 event suggests that the loss of material did not cause significant exposure of volatile-rich patches on the surface of the nucleus.  Ejection of a surface layer by gas pressure from within the nucleus would likely have left an icy crater exposed that would have remained active for days or weeks as volatile deposits diminished.
The event could have involved the dislodging of a fragment that was only weakly bound to the surface.  The fragment would then have broken apart into a collection of particles with a much greater total cross-section.  Rotational spin-up upon ejection, or the sublimation of volatiles that were holding the fragment together could have caused the break-up.
The nucleus may have remained unstable for a significantly longer period of time given that \cite{2009SASS...28...51M} detected a 0.85 mag brightening on UT 2009 Jan.\ 4.  It is likely that 17P/Holmes will continue to undergo minor outbursts for the next few orbital periods as unstable regolith settles and freshly-exposed volatiles sublimate.

\subsection{Surface Brightness Profile}

The distribution of material in the inner coma can be studied using the surface brightness profile.  We fitted the gradient of each surface brightness profile between 10000 km and 25000 km using a power law.  The gradient, $m_{p}$, is given by

\begin{equation}
m_{p} = \frac{d\ ln\ B(R)}{d\ ln\ R}
\label{meq}
\end{equation}

where $B(R)$ is the flux per 10$^{6}$ km$^{2}$, and $R$ (km) is the distance from the nucleus.  A dust coma generated by an isotropic, steady source would have a slope $m_{p}$ = -1 \citep{1991ASSL..167...19J}.  Instead, the measured slope varied between -0.2 and -0.3 with no significant temporal trend (Table~\ref{table:mtable}).
We can be confident that the shallow slopes are not a result of inaccurate sky subtraction caused by coma-contamination.  If the estimated sky value actually included a contribution from the dust coma, the images would be over-subtracted, resulting in an artificially steep slope.  The shallow slope suggests that either the dust production rate varied in such a way as to temporarily produce a smaller scattering cross-section close to the nucleus, or that the cross-section increased with increasing nucleo-centric distances, possibly through a process such as dust fragmentation.  The likelihood of the activity on the nucleus varying with time in such a way as to consistently produce a shallow surface brightness profile for several months is extremely low. Another process is needed to account for the shallow surface brightness gradients.

Cometary ejecta is likely to include agglomerate grains, which are clusters of smaller particles that are loosely held together by electrostatic forces or volatile material \citep{2004JGRE..10912S03C}.  Heating of these tenuous aggregates by incident solar radiation could cause them to break apart, greatly increasing their scattering cross-section.  For example, the brightness of a 1 cm radius cluster of micron-sized dust grains could increase by a factor of $\sim$ 10$^{4}$ if it broke apart into its constituent particles.   Thus, fragmenting or disaggregating particles could cause a shallow surface brightness profile where the scattering cross-section of material increases with distance from the nucleus.  \cite{2009AJ....137.4538Y} identified populations of both hot and cold grains in the inner coma of 17P/Holmes using infrared spectra.  They postulated that the cold ice grains and hot dust grains could have originally been contained within a porous aggregate that broke apart when exposed to sunlight.  We conjecture that fragmenting aggregates were responsible for the shallow profile at small nucleo-centric distances ($<$ 25000 km), \cite{2010MNRAS.407.1784H} suggested that this process continued out to distances of greater than 10$^{5}$ km.

\subsection{Ice Grain Halo}

The surface brightness profiles consistently show a bump in the innermost region before the profiles straighten out to a decreasing slope (Figure~\ref{fig:SBall}).  To investigate the nature of this bump we created a model surface brightness profile for each night with a slope matching that measured at large angular radii in the CFHT data, plus an unresolved central source.  The model was forced to match those brightness levels observed at 250 km and 10000 km, and then convolved with a 1$^{\prime\prime}$ Gaussian function to simulate the effects of seeing.  Figure~(\ref{fig:SBcomp}) compares the results of our model (dashed line) to the data (solid line).  The comparison shows that the brightness of the central region cannot simply be attributed to a smooth coma with a contribution from the nucleus.  One possible interpretation is that a halo of icy sublimating grains persisted around the nucleus.  Ice grain halos in comets were suggested decades ago \citep{1970P&SS...18..717D} and evidence for halos has been reported in several comets since then, including  C/1996 B2 Hyakutake \citep{1997Sci...277..676H} and 103P/Hartley 2 \citep{2011ApJ...734L...1M}.  \cite{2009AJ....138.1062S} found that derived water production rates for 17P/Holmes during the month following the outburst were too large to be explained by sublimation of water ice from the surface of the nucleus.  He concluded that a population of icy grains must have been sublimating near the nucleus, which fits well with the ice grain halo scenario.

The extent of an ice grain halo is dependent on the sublimation rate of the ice grains, which in turn depends on several factors including the size, albedo, temperature, and density of the grains.  To explore these factors, we calculate the sublimation rate assuming thermal equilibrium.  The sublimation rate is expressed in meters per second, representing the rate at which an icy surface would retreat under these conditions.  We assume two cases: isothermal grains where the heat absorbed is radiated isotropically back into space (``cold grains''), and grains at the temperature of the sub-solar point (``hot grains'').  The sublimation rate is strongly dependent on the Bond albedo, which is the fraction of total power incident on a body that is reflected from the surface.  High-albedo pure ice grains have low sublimation rates since they absorb little heat and remain cold.  However, adding even a small amount of refractory material to the ice lowers the albedo significantly, resulting in greater absorption of heat and an increase in the temperature of the grain \citep{1982Icar...49..244C}.

The radius of the ice grain halo, $L$, is related to the grain velocity, $v(a)$ (m s$^{-1}$), radius, $a$ (m), and sublimation rate, $\dot{a}$ (m s$^{-1}$), by:

\begin{equation}
L = v(a) \frac{a}{\dot{a}}
\label{extenteq}
\end{equation}

We utilized the velocity-radius relationship given in Equation~(\ref{veleq}) and measured $L$ by comparing the surface brightness profile on each night to the model profile.  We defined the edge of the halo as where the brightness fell to within 10\% of the model.  We acknowledge that this approach is model-dependent but use it to explore the concept of an ice grain halo in this context.  The radius determined varied between $\sim$ 2900 km and 4900 km (Table~\ref{table:mtable}).  The median radius was 3450 km, consistent with results from \cite{2009AJ....137.4538Y} who found that ice grains persisted up to 3500 km from the nucleus in 2007 Oct.\  We plot the radius-albedo relations for hot and cold ice grains that can survive across a 3450 km radius aperture in Figure~(\ref{fig:avspv}).  Cold ice grains with a radius of 1 $\mu$m (c.f. \citealt{2009AJ....137.4538Y}) must have a Bond albedo of $\sim$ 0.35 in order to survive, while micron-sized grains in the high temperature limit must have a higher Bond albedo, of $\sim$ 0.85.  Of these two approximations, the cold grain case is probably more relevant to ice grains in 17P/Holmes, given their small size and the likelihood that they were rotating rapidly.

The results of \cite{2010ApJ...714.1324I} suggest that the geometric albedo of the coma was between $\sim$ 0.05 and 0.20.  We note that the results from \cite{2010ApJ...714.1324I} are based on measurement of the flux reflected from all material in a 5300 km radius aperture, which would have included a significant contribution from lower-albedo dust grains as well as high-albedo ice grains.  In this sense, our result that the Bond albedo of ice grains in the halo must have been $\sim$ 0.35 appears to be consistent with their findings.

Comparatively little contaminating material is required to significantly darken icy material.  \cite{1982Icar...49..244C} found that contaminating 1 mm snow grains with 0.1 $\mu$m soot particles that amounted to only 10$^{-5}$ of the sample weight lowered the albedo from 0.94 to 0.40.  Results from \cite{2009AJ....137.4538Y}  suggested that only ice grains with $\lesssim$ 10\% impurities could survive long enough to travel 3500 km from the nucleus before sublimating.  \cite{2006Icar..180..473B} found that micron-sized pure ice grains beyond $\sim$ 2 AU should survive for weeks or even months, while dirty ice grains (10\% dirt by mass) should survive for only hours, which agrees well with our results.

The ice grain halo radius varied during our observing campaign, although it is difficult to discern any obvious trends (Figure~\ref{fig:extentfit}).  We note that the halo was larger after periods of known activity, such as on 2007 Nov.\ 6 - our observation closest to the original outburst - and on 2007 Nov.\ 12 after the much smaller second outburst.  The radius appeared to decrease after these periods of activity until the heliocentric distance of 17P/Holmes began to change appreciably.  As the heliocentric distance increased from 2.63 AU on 2007 Dec.\ 11 to 2.90 AU on 2008 Feb.\ 12, the temperature of the ice grains must have dropped, causing a decrease in sublimation rates and thus, an increase in the survival time of the grains.  We attempted to fit the increase of the ice grain halo radius using the thermal equilibrium model described previously.  Figure~(\ref{fig:extentfit}) shows the data with results from models of the halo radii expected during our observations, assuming the halo consists of isothermal grains.  The increase in heliocentric distance by itself cannot account for the halo size variations.  We suggest that changes in the level of activity, size and albedo of the ice grains ejected, or the grain velocities may also play important roles in the extent of the ice grain halo.

Finally, we compare the expected mass loss rate from a sublimating ice grain halo to the observed production rates of OH molecules using aperture photometry from this work and production rates taken from \cite{2009AJ....138.1062S}.  The production rates were derived from narrowband photometry of 17P/Holmes' inner coma between UT 2007 Nov.\ 1 and 2008 Mar.\ 5.  Since the published production rates were corrected for 10$^{4}$ km radius apertures, we also use 10$^{4}$ km radius apertures.  We convert the apparent magnitudes to geometric scattering cross-sections using Equation~\ref{crosssectioneq} and assumptions detailed in Section~\ref{sec:fading}.  To estimate the mass loss rate expected from a sublimating ice grain halo, we multiply the cross-section of material (presumed to be 100\% ice) by the density of the ice (assumed to be 1000 kg m$^{-3}$) and the sublimation rate (m s$^{-1}$) for ice grains at the corresponding heliocentric distances of 17P/Holmes.  Figure~\ref{fig:schcomp} compares the expected mass loss rates for ice grain halos composed of grains with albedos of 0 and 0.3 to the derived production rates from \cite{2009AJ....138.1062S}.  We find that the mass loss rate for an ice halo composed of grains with albedo = 0.3 drops off more quickly than the derived water production rate.  The changing mass loss rate for an ice grain halo with grains of albedo = 0 matches the trend in the water production rate well but the albedo conflicts with our previous determination that the extent of the halo is best modeled by sublimating ice grains with albedos of $\sim$ 0.2 - 0.4.  We acknowledge that the model is underconstrained given that the ice grains may not be isothermal, would have a range of sizes and albedos, and those that determine the observed radius of the ice grain halo may not be the grains that supply the sublimated gas.  We conclude by simply saying that the decrease in water production rates agrees loosely with the expected decrease in mass loss rates for an ice grain halo consisting of low-albedo ice grains.

\section{Summary}

We have used aperture photometry to examine the inner coma of 17P/Holmes between UT 2007 Nov.\ 6 and 2008 Feb.\ 12 following its massive outburst in 2007 Oct.\  Our results include:

\begin{enumerate}

\item{The nucleus of 17P/Holmes remained active for at least several months following its spectacular outburst in 2007 Oct., although at a level very low compared to the peak outburst rates.  The inner coma decreased in brightness by a factor of $\sim$ 50 during our observations, but an examination of aperture-crossing times shows that fresh material must have been continually supplied from active regions on the nucleus.}

\item{A second much smaller outburst occurred on UT 2007 Nov.\ 12, ejecting approximately 8 $\times$ 10$^{5}$ kg of material that had traveled up to 5000 km by the time of our observations.  We propose that a fragment broke away from the surface and rapidly broke apart, releasing a significant amount of dust into the coma.}

\item{The surface brightness profile of 17P/Holmes had a measured slope between -0.2 and -0.3.  This is less steep than expected for a steady-state coma consistent with the fragmentation of icy grain aggregates.}

\item{A halo consisting of freshly released ice grains persisted around the nucleus, extending to radii of several thousand km.  This halo was likely responsible for the anomalously high water production rates observed by \cite{2009AJ....138.1062S} several months after the outburst, although the bulk of gas and dust production could originate from distinctly different populations of grains.}

\item{The radius of the ice grain halo is compatible with an origin in isothermal micron-sized grains having a Bond albedo of $\sim$ 0.35.  The ice grains may have been contaminated with small amounts of refractory material.}

\end{enumerate}

Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l'Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii.  This work is based in part on data products produced at the TERAPIX data center located at the Institut d'Astrophysique de Paris.  We appreciate grant support from the NASA Outer Planets Research Program to DJ.  RS acknowledges support from the NASA Postdoctoral Fellowship Program.  Part of the research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.

\clearpage


\section{{Introduction and Motivation}}
\begin{prob}
$7$ divides $\binom{n}{7}-\lfloor\frac{n}{7}\rfloor$, $\forall n \in \mathbb{N}$.
\end{prob}

The above appeared as a problem in the Regional Mathematical Olympiad, India in 2003. Later in 2007, a similar type of problem was set in the undergraduate admission test of Chennai Mathematical Institute, a premier research institute of India where $7$ was replaced by $3$.

This became the basis of the following

\begin{thm}[\cite{mps}, Saikia-Vogrinc]\label{mps-jv}
A natural number $p>1$ is a prime if and only if $\binom{n}{p}-\lfloor\frac{n}{p}\rfloor$ is divisible by $p$ for every non-negative $n$, where $n>p+1$ and the symbols have their usual meanings.
\end{thm}

\section{{Proof of Theorem \ref{mps-jv}}}

In \cite{mps}, the above theorem is proved. The authors give three different proofs, however the third proof is incomplete. We present below a completed version of that proof.

\begin{proof}
 First we assume that $p$ is prime. Now we consider $n$ as $n=ap+b$ where $a$ is a non-negative integer and $b$ an integer $0\leq b<p$. Obviously,
\begin{equation}
\left\lfloor \frac{n}{p}\right\rfloor=\lfloor \frac{ap+b}{p}\rfloor\equiv a~(mod~p).
\end{equation}
Now let us calculate $\binom{n}{p}~(mod~p)$.
\begin{eqnarray*}
\binom{n}{p}&=&\binom{ap+b}{p}\\
&=&\frac{(ap+b)\cdot(ap+b-1)\cdots(ap+1)\cdot ap\cdot(ap-1)\cdots(ap+b-p+1)}{p\cdot(p-1)\cdots 2\cdot1}\\
&=&\frac{a\cdot(ap+b)\cdot(ap+b-1)\cdots(ap+1)\cdot(ap-1)\cdots(ap+b-p+1)}{(p-1)\cdot(p-2)\cdots 2\cdot 1}\\
&=&\frac{a X}{(p-1)!}
\end{eqnarray*}
where $X=(ap+b)\cdot(ap+b-1)\cdots(ap+1)\cdot(ap-1)\cdots(ap+b-p+1)$.

We observe that there are $(p-1)$ terms in $X$ and each of them has one of the following forms,\\ (a) $ap+r_1$, or\\ (b) $ap-r_2$\\ where $1\leq r_1\leq b$ and $1\leq r_2\leq (p-1-b)$.\\ Thus any two terms from either $(a)$ or $(b)$ differs by a number strictly less than $p$ and hence not congruent modulo $p$. Similarly, if we take two numbers - one from $(a)$ and the other from $(b)$, it is easily seen that the difference between the two would be $r_1+r_2$ which is at most $(p-1)$ (by the bounds for $r_1$ and $r_2$); thus in this case too we find that the two numbers are not congruent modulo $p$. Thus the terms in $X$ forms a reduced residue system modulo $p$ and so, we have,
\begin{eqnarray}
X\equiv(p-1)!\;(mod\;p)
\end{eqnarray}
Thus using $(2.2)$ we obtain,
\begin{eqnarray}
\binom{n}{p}=a\frac{X}{(p-1)!}\equiv a\;(mod\;p)
\end{eqnarray}
So, $(2.1)$ and $(2.3)$ combined gives
\begin{eqnarray}
\left\lfloor \frac{n}{p}\right\rfloor\equiv \binom{n}{p}\;(mod\;p).
\end{eqnarray}
So, forward implication is proved.

To prove the reverse implication, we adopt a contrapositive argument meaning that if $p$ were not prime (that is composite) then we must construct an $n$ such that $(4)$ does not hold. So, let $q$ be a prime factor of $p$. We write $p$ as $p=q^x k$, where $(q,k)=1$. In other words, $x$ is the largest power of $q$ such that $q^x|p$ but $q^{x+1}\!\not|\,p$ (in notation, $q^x || p$). By taking, $n=p+q=q^xk+q$, we have
\begin{eqnarray*}
\binom{p+q}{p}=\binom{p+q}{q}=\frac{(q^xk+q)(q^xk+q-1)\dots (q^xk+1)}{q!}
\end{eqnarray*}
which after simplifying the fraction equals $(q^{x-1}k+1)\frac{(q^xk+q-1)\dots (q^xk+1)}{(q-1)!}$. Clearly, $(q^xk+q-1)\dots (q^xk+1)\equiv (q-1)!\not\equiv 0\;(mod\;q^x)$. Therefore,
\begin{eqnarray*}
\frac{(q^xk+q-1)\dots (q^xk+1)}{(q-1)!}\equiv 1~(mod~q^x)
\end{eqnarray*}
 and
 \begin{eqnarray*}
\binom{p+q}{p}\equiv q^{x-1}k+1~(mod~q^x).
 \end{eqnarray*}
On the other hand obviously,
\begin{eqnarray*}
\left\lfloor\frac{p+q}{p}\right\rfloor=\left\lfloor\frac{q^xk+q}{q^xk}\right\rfloor\equiv 1~(mod~q^x).
\end{eqnarray*}
Now, since $(q,k)=1$, it follows that $q^{x-1}k+1\not\equiv 1\;(mod\;q^x)$. So we conclude,
\begin{eqnarray}
\binom{p+q}{p}\not\equiv \left\lfloor\frac{p+q}{p}\right\rfloor\;(mod\;q^x).
\end{eqnarray}
So, $p\nmid (\binom{p+q}{p}-\lfloor\frac{p+q}{p}\rfloor)$, for if $p|(\binom{p+q}{p}-\lfloor\frac{p+q}{p}\rfloor)$, then since $q^{x}|p$, we would have $q^x|(\binom{p+q}{p}-\lfloor\frac{p+q}{p}\rfloor)$, a contradiction to  $(5)$. Thus, $\binom{p+q}{p}\not\equiv\lfloor\frac{p+q}{p}\rfloor\;(mod\;p).$ Hence we are through with the reverse implication too.

This completes the proof of Theorem \ref{mps-jv}.
\end{proof}

\section{{Another simple result}}

We state and prove the following simple result which generalizes one part of Theorem \ref{mps-jv}

\begin{thm}\label{gen}
For $n=ap+b=a_{(k)}p^k+b_{(k)}$, we have  $$
{a_{(k)}p^k+b_{(k)}\choose p^k}
-\left\lfloor\frac{a_{(k)}p^k+b_{(k)}}{p^k}\right\rfloor\equiv
0~(\textup{\textup{mod}}~p)
$$
with $p$ a prime, $0\leq b_{(k)}\leq p^k-1$ and $k$ a positive integer
such that $1\leq k\leq l$, where $$n=a_0+a_1p+\ldots+a_kp^k+a_{k+1}p^{k+1}+\ldots+a_lp^l$$ and
for $k\geq 1$
$$
a_{(k)}=a_k+a_{k+1}p+\ldots+a_lp^{l-k}
$$
and
$$
b_{(k)}=a_0+a_1p+\ldots+a_{k-1}p^{k-1}.
$$.
\end{thm}

The proof of this follows from the reasoning of the proof of Theorem \ref{mps-jv} although there are some subtleties.

In particular, we have
$$
a=a_{(1)}=a_1+a_{2}p+\ldots+a_lp^{l-1}
$$
and
$$
b=b_{(0)}=a_0.
$$
For $k=0$, we set the convention that $a_{(0)}=n=a_0+a_{1}p+\ldots+a_lp^{l}
$ and $b_{(0)}=0$.
Notice that Theorem \ref{gen} is obviously true for $k=0$. But the
case $k=0$ doesn't correspond really to a power of $p$ where $p$ is a
prime.

\begin{proof}
We have
$$\binom{n}{p^k}=\binom{a_{(k)}p^k+b_{(k)}}{p^k}$$
$$=\frac{(a_{(k)}p^k+b_{(k)})\cdot(a_{(k)}p^k+b_{(k)}-1)\cdots(a_{(k)}p^k+1)\cdot
a_{(k)}p^k\cdot(a_{(k)}p^k-1)\cdots(a_{(k)}p^k+b_{(k)}-p^k+1)}
{p^k\cdot(p^k-1)\cdots 2\cdot1}$$
$$=\frac{a_{(k)}\cdot(a_{(k)}p^k+b_{(k)})\cdot(a_{(k)}p^k+b_{(k)}-1)\cdots
(a_{(k)}p^k+1)\cdot(a_{(k)}p^k-1)\cdots(a_{(k)}p^k+b_{(k)}-p^k+1)}
{(p^k-1)\cdot(p^k-2)\cdots 2\cdot 1}.$$ Thus we obtain
$$
(p^k-1)!\binom{n}{p^k}=a_{(k)}\left({\displaystyle\prod^b_{r=1}(a_{(k)}p^k+r)}
\right)\,
\left({\displaystyle\prod^{p^k-1-b}_{r=1}(a_{(k)}p^k-r)}\right).
$$ Or $a_{(k)}p^k+r\equiv r~(\textup{\textup{mod}}~p^k)$ and
$a_{(k)}p^k-r\equiv -r\equiv p^k-r~(\textup{\textup{mod}}~p^k)$ with
$0<r<p^k$. It follows
$$
\left({\displaystyle\prod^b_{r=1}(a_{(k)}p^k+r)}
\right)\,
\left({\displaystyle\prod^{p^k-1-b}_{r=1}(a_{(k)}p^k-r)}\right)
\equiv\left({\displaystyle\prod^b_{r=1}r}\right)\,
\left({\displaystyle\prod^{p^k-1-b}_{r=1}(p^k-r)}
\right)~(\textup{\textup{mod}}~p^k).
$$ Since
$$
\left({\displaystyle\prod^b_{r=1}r}\right)\,
\left({\displaystyle\prod^{p^k-1-b}_{r=1}(p^k-r)}\right)
=\left({\displaystyle\prod^b_{r=1}r}\right)\,
\left({\displaystyle\prod^{p^k-1}_{r=b+1}r}\right)
={\displaystyle\prod^{p^k-1}_{r=1}r}=(p^k-1)!
$$ we have
$$
\left({\displaystyle\prod^b_{r=1}(a_{(k)}p^k+r)}
\right)\,
\left({\displaystyle\prod^{p^k-1-b}_{r=1}(a_{(k)}p^k-r)}\right)
\equiv(p^k-1)!~(\textup{\textup{mod}}~p^k).
$$

We can notice that,
$$
(p^k-1)!=q(p-1)!\,p^{1+p+\ldots+p^{k-1}-k}
$$
with $\textup{gcd}(p,q)=1$ and because $\textup{ord}_{p}((p^k-1)!)=1+p+\ldots+p^{k-1}-k$.
Therefore we have
$$
a_{(k)}c_{(k)}p^{k(p^k-1)}+(p^k-1)!\left\{a_{(k)}-\binom{n}{p^k}\right\}=0.
$$
Equivalently
$$
a_{(k)}c_{(k)}p^{k(p-1)(1+p+\ldots+p^{k-1})}+q(p-1)!\,p^{1+p+\ldots+p^{k-1}-k}
\left\{a_{(k)}-\binom{n}{p^k}\right\}=0
$$

Dividing the above equation by $p^{1+p+\ldots+p^{k-1}-k}$ we have
$$
q(p-1)!\left\{a_{(k)}-\binom{n}{p^k}\right\}
+a_{(k)}c_{(k)}p^{k+(k(p-1)-1)(1+p+\ldots+p^{k-1})}=0.
$$
Thus
$$
q(p-1)!\left\{a_{(k)}-\binom{n}{p^k}\right\}\equiv 0~(\textup{\textup{mod}}~p^k)
$$
Since if $m\equiv
n~(\textup{\textup{mod}}~p^k)$ implies $m\equiv
n~(\textup{\textup{mod}}~p)$ (the converse is not always
true), we also have
$$
q(p-1)!\left\{a_{(k)}-\binom{n}{p^k}\right\}\equiv 0~(\textup{\textup{mod}}~p).
$$
As $q(p-1)!$ with $\textup{gcd}(p,q)=1$ and $p$ are relatively prime, we get
$$
\binom{n}{p^k}-a_{(k)}\equiv 0~(\textup{\textup{mod}}~p).
$$

We finally have
$$
\binom{n}{p^k}\equiv
\left\lfloor\frac{n}{p^k}\right\rfloor~(\textup{\textup{mod}}~p).
$$
\end{proof}

\begin{thm}\label{t3.2}
Let $p$ be a prime number, let $k$ be a natural number and
let $x$ be a positive integer such that
$$
x\equiv r\pmod {p^k}
$$
with $0\leq r<p^k$. Denoting $q=\lfloor\frac{x}{p^k}\rfloor$ the
quotient of the division of $x$ by $p^k$, if
there exists $s\in\mathbb{N}^{\star}$ for which
$$
\lfloor\frac{x}{p^{ks}}\rfloor=q^s
$$
then we have
$$
x\equiv r\pmod {p^{ks}}.
$$
\end{thm}

\begin{proof}
Given $p$ a prime number, let $x$ be a positive integer such that
$x\equiv r\pmod {p^k}$ with $0\leq r<p^k$. Denoting
$q=\lfloor\frac{x}{p^k}\rfloor$, we assume that
there exists $s\in\mathbb{N}^{\star}$ for which
$$
\lfloor\frac{x}{p^{ks}}\rfloor=q^s.
$$
If $k=0$, the result is obvious since for all integers $x,r$, we have
$x\equiv r\pmod 1$. In the following, we assume that $k\in\mathbb{N}^{\star}$.
\\[0.1in]
Then, we have
$$
x=qp^k+r
$$
and
$$
x=q^sp^{ks}+r'
$$
with $0\leq r'<p^{ks}$. It comes that
$$
qp^k+r=q^sp^{ks}+r'.
$$
So ($s\in\mathbb{N}^{\star}$)
$$
q^sp^{ks}-qp^k=r-r'\geq 0.
$$
Since $0\leq r<p^k$, we have
$$
0\leq r-r'<p^k.
$$
Moreover, rewriting the equality $q^sp^{ks}-qp^k=r-r'$ as
$$
qp^k(q^{s-1}p^{k(s-1)}-1)=r-r'
$$
we can notice that $p^k|r-r'$. Since $0\leq r-r'<p^k$, it is only possible if
$r-r'=0$ and so
$$
r=r'.
$$
From the equality $x=q^sp^{ks}+r'$, we deduce that
$$
x\equiv r\pmod {p^{ks}}.
$$
\end{proof}

A consequence of the Theorem \ref{t3.2} is that if an integer $y$ is
congruent to a positive integer $x$ modulo $p^k$ such that $x\equiv
r\pmod {p^k}$, provided the
conditions stated in the Theorem \ref{t3.2} are fulfilled, we have also
$y\equiv r\pmod {p^{ks}}$.

It can be verified easily that the product
$\left(\prod^b_{r=1}(a_{(k)}p^k+r)\right)\,
\left(\prod^{p^k-1-b}_{r=1}(a_{(k)}p^k-r)\right)$ contains the term
$(a_{(k)}p^k)^{p^k-1}=a_{(k)}^{p^k-1}p^{k(p^k-1)}$. The term
$(a_{(k)}p^k)^{p^k-1}$ is the only
term in $p^{k(p^k-1)}$ which appears in the decomposition of this product into
sum of linear combination of powers of $p$. Notice also that the number
$k(p^k-1)$ is the greatest exponent
of $p$ in this product
when we decompose this product into sum of linear combination
of powers of $p$ (like a polynomial expression in variable
$p$). Afterwards, we write
$a_{(k)}^{p^k-1}$ as $c_{(k)}$ in order to simplify the notation. Thus,
the quotient of the division of this product by
$p^{k(p^k-1)}$ is $c_{(k)}=a_{(k)}^{p^k-1}$.

So, from the Theorem \ref{t3.2}, we can now write
$$
\left({\displaystyle\prod^b_{r=1}(a_{(k)}p^k+r)}
\right)\,
\left({\displaystyle\prod^{p^k-1-b}_{r=1}(a_{(k)}p^k-r)}\right)
=c_{(k)}p^{k(p^k-1)}+(p^k-1)!.
$$

\section{{Acknowledgements}}
The authors are grateful to Professor Nayandeep Deka Baruah for going through an earlier version of this work and offering various helpful suggestions which
helped to make the presentation much clearer. The authors are also thankful to Ankush Goswami for cleaning up the exposition in the proof of Theorem \ref{mps-jv}. The second named author would also like to thank Bishal Deb for pointing out to him a factual error in \cite{mps} and also to Parama Dutta for a careful reading of the first two sections of this paper.


\section{Vague thoughts about varieties}

\begin{bibdiv}
\begin{biblist}

\bib{aw3}{article}{
    author = {Antieau, Benjamin},
    author = {Williams, Ben},
    title = {The topological period-index problem for $6$-complexes},
    journal = {J. Top.},
    eprint = {http://arxiv.org/abs/1208.4430},
    note = {To appear.}
}

\bib{aw4}{article}{
    author = {Antieau, Benjamin},
    author = {Williams, Ben},
    title = {Unramified division algebras do not always contain Azumaya maximal orders},
    journal = {Inv. Math.},
    eprint = {http://dx.doi.org/doi:10.1007/s00222-013-0479-7},
    year = {2013},
}

\bib{bott}{article}{
    author={Bott, Raoul},
    title={The space of loops on a Lie group},
    journal={Michigan Math. J.},
    volume={5},
    date={1958},
    pages={35--61},
    issn={0026-2285},
}

\bib{cadek-vanzura}{article}{
    author={\v{C}adek, Martin},
    author={Van{\v{z}}ura, Ji{\v{r}}{\'{\i}}},
    title={On the classification of oriented vector bundles over $5$-complexes},
    journal={Czechoslovak Math. J.},
    volume={43(118)},
    date={1993},
    number={4},
    pages={753--764},
    issn={0011-4642},
    eprint={http://www.math.muni.cz/~cadek/list.html},
}

\bib{dold-whitney}{article}{
    author={Dold, A.},
    author={Whitney, H.},
    title={Classification of oriented sphere bundles over a $4$-complex},
    journal={Ann. Math.},
    date={1959},
    volume={69},
    pages={667--677},
}

\bib{serre}{article}{
    author={Serre, Jean-Pierre},
    title={Cohomologie modulo $2$ des complexes d'Eilenberg-MacLane},
    journal={Comment. Math. Helv.},
    volume={27},
    date={1953},
    pages={198--232},
    issn={0010-2571},
}

\bib{woodward}{article}{
    author={Woodward, L. M.},
    title={The classification of principal ${\rm PU}_{n}$-bundles over a $4$-complex},
    journal={J. London Math. Soc. (2)},
    volume={25},
    date={1982},
    number={3},
    pages={513--524},
    issn={0024-6107},
}

\bib{woodward-orientable}{article}{
    author={Woodward, L. M.},
    title={The classification of orientable vector bundles over CW-complexes
    of small dimension},
    journal={Proc. Roy. Soc. Edinburgh Sect. A},
    volume={92},
    date={1982},
    number={3-4},
    pages={175--179},
    issn={0308-2105},
}

\end{biblist}
\end{bibdiv}

\end{document}


\section{Introduction}

The ability of scanning probe microscopy (SPM) techniques to probe and manipulate electronic states on a scale smaller than typical coherence lengths has allowed to directly visualize electron interference, the main fingerprint of quantum phenomena.  The results can be spectacular, as witnessed by the coherent electron flow pictures observed in earlier scanning gate microscopy (SGM) measurements by the Westervelt group in Harvard.\cite{Topinka2000,Topinka2001,Leroy2005}  In these experiments, the negatively biased tip of an atomic force microscope (AFM) was exploited to backscatter electrons transmitted across a quantum point contact (QPC) in a two-dimensional electron gas (2DEG). The interference among different backscattering paths produces a modulation of the transmission probability which is visible in the SGM images as fringe structures with spacing $\lambda_F/2\approx 20$~nm, where $\lambda_F$ is the Fermi wavelength. The high resolution required to resolve such structures  is \textit{not limited by the width of the tip-induced potential}: it just depends on the accuracy of tip positioning. This idea is at the basis of most experiments exploiting the SGM technique.

In this paper we apply the same technique to study the backscattering mechanisms in wide constrictions defined in quantum Hall (QH) systems. When the bulk filling factor $\nu_b$ is integer, the current is carried by chiral edge channels.\cite{Halperin1982} Our constriction is sufficiently wide to allow full transmission when the SGM tip is far away from the constriction axis. When the negatively biased  tip is moved towards the constriction, backscattering is induced between counter-propagating edge-channels. The resulting SGM images, obtained by mapping the transmitted linear conductance  $G_T$ as a function of tip position, show that $G_T$  monotonically decreases when the tip approaches the constriction, with plateaus in correspondence to multiples of $G_0= e^2/h$. This behaviour is well known for short QPCs.\cite{Paradiso2010} However, our SGM scans over wide constrictions show a surprising fine structure consisting in small spatial oscillations of $G_T$ that lead to arc features in the SGM maps. Based on a detailed comparison with a numerical model, the emergence of such arc features can be attributed to backscattering through the constriction mediated by naturally-occurring antidot (AD) structures originating from potential fluctuations present in the constriction.

Arc structures have been observed before in SGM experiments performed on a variety of different systems, with or without magnetic field (e.g.~carbon nanotubes,\cite{Woodside2002} quantum rings,\cite{Hackens2006,Hackens2010} InAs nanowires,\cite{Bleszynski2007} and quantum dots\cite{Fallahi2005,Gildemeister2007}) and can be caused by a tip-induced modulation of either quantum interference \cite{Hackens2006} or Coulomb blockade (CB).\cite{Woodside2002,Bleszynski2007,Fallahi2005,Gildemeister2007}  Indeed, both mechanisms can be expected to play a role in artificial AD devices, as highlighted by a number of investigations in past years reviewed in Ref.~\onlinecite{Ford} and by recent experiments specifically designed to address this issue.\cite{Kou,Markus2} In our experiment, the SGM tip modulates the transport through random localized states, inducing a suitable coupling between them and the extended edge modes. For specific tip configurations which will be discussed in the paper, we induce in a highly controlled manner AD backscattering paths which would appear as conductance fluctuations in transport experiments upon pinching-off the constriction. This mechanism is absent or negligible in other tip modulation schemes such as, for instance, remote gating in quantum dots.\cite{Woodside2002,Bleszynski2007,Fallahi2005,Gildemeister2007}  It should be noted that similar arc structures have also been observed by the scanning charge accumulation technique in the QH regime.\cite{Finkelstein2000,Ashoori1}  Even if these results are linked to edge-edge scattering mediated by localized states, the mechanism behind the arc formation is different: in Refs.~\onlinecite{Finkelstein2000,Ashoori1} arcs have an angular extension which reflects the incompressible stripe modulation given by the combined action of the tip electric field and the local average field at the active edge-edge scattering center; in our work the arc angular extension rather reflects the position of the antidot with respect to the constriction borders, as discussed in the following sections.

In Ref.~\onlinecite{Paradiso2011} we studied the nature and the impact of impurity potentials on the co-propagating edge channel mixing. In our constrictions we found typically a few strong scattering centers per square micron. In order to study the behaviour of these centers as ADs linking counter-propagating edge channels crossing the constriction, we designed a constriction which contains several scattering centers between the channels. For this reason, the constriction area ($\approx$ 3~$\mu$m$^2$) in the present experiment is much larger than in our earlier work on short QPCs,\cite{Paradiso2010} which aimed at revealing the inner edge structure.

\section{Experimental details}

The samples for this study were fabricated starting from an Al$_{0.33}$Ga$_{0.67}$As/GaAs heterostructure with a two-dimensional
electron gas (2DEG), which is confined 55~nm underneath the surface. The electron sheet density and mobility at low temperature are $n = 3.12\times 10^{15}$~m$^{-2}$ and $\mu = 4.2\times 10^{2}$~m$^2$/Vs, respectively, as determined by Shubnikov-de Haas measurements. The Hall bar was patterned via optical lithography and wet etching. Source and drain contacts were fabricated at the ends of the Hall bar by evaporation and thermal annealing of a standard Ni/AuGe/Ni/Au multilayer (10/200/10/100~nm). In order to define a constriction in the 2DEG, we fabricated two split-gates via electron beam lithography. They consist of a Ti/Au bilayer (10/20~nm). The nominal gap between the gates is 1.2~$\mu$m wide and 2.5~$\mu$m long.

Our measurements were performed with the 2DEG at bulk filling factor $\nu_b=2$ ($B= 6.450$~T). The split--gate bias was set to $-0.350$~V, so that the filling factor under the gates is $g=0$, and both edge channels of the $\nu_b=2$ bulk phase are
sent into the constriction. The width of the constriction has been chosen large enough to allow full transmission of both edge channels, so that the source-drain transmitted  conductance, in the absence of the tip, is $G_T=\nu_b G_0 = 2e^2/h$.

 The SGM was operated in a $^3$He cryostat (base temperature 300~mK). Sample temperature was 400~mK, as calibrated with a Coulomb blockade thermometer. The cryostat is equipped with a superconducting coil which provides magnetic fields of up to 9~T.
Details of our SGM setup are reported in Ref.~\onlinecite{Paradiso2010}. SGM maps are obtained by scanning a negatively biased AFM tip over the constriction area and acquiring for each position $(x_t,y_t)$ the corresponding source-drain transmitted  conductance $G_T(x_t,y_t)$.

\section{Experimental results}

Figure~\ref{fig:general}(a) shows the result of a scan over the whole constriction area, performed at $B=6.450$~T and $V_{tip}=-6.0$~V. When the SGM tip is far from the constriction, there is no backscattering between counter-propagating edge states, hence the source-drain conductance transmitted across the constriction is $G_T = 2 G_0 = 2e^2/h$.  When the tip approaches the constriction axis, tunneling between the counter-propagating edges is induced, and $G_T$ is reduced, similarly to what was observed in previous SGM experiments on short QPCs.\cite{Aoki2005,Paradiso2010} The detailed shape of contour lines at equal conductance is different from device to device (for six devices measured in our experiments) and depends on the actual disorder potential inside the constriction. An extended area with constant conductance belongs to the $G_T=G_0$ plateau and is indicated in Fig.~\ref{fig:general}(a) with a dashed white line.

\begin{figure}[tb]
\includegraphics[width=\columnwidth]{Fig1.eps}
\caption{(a) SGM scan over a 1.2~$\mu$m wide and 2.5~$\mu$m long constriction defined in a $\nu_b=2$ quantum Hall system ($B=6.450$~T, $V_{tip}=-6.0$~V). The map displays the transmitted conductance $G_T(x_t,y_t)$ plotted as a function of the tip position. As a reference, we indicated the middle of the $G_T=G_0$ plateau with a dashed white line. (b) Gradient map $|\mathbf{\nabla}G_T(x_t,y_t)|$ which emphasizes short scale modulations of the $G_T$ signal. A fine pattern of arc structures is observed, indicating the occurrence of peculiar backscattering processes.}
\label{fig:general}
\end{figure}

The present QPC design allows to further investigate the backscattering mechanism which governs transport in the constriction. A number of interesting experimental features can be highlighted by plotting the $|\mathbf{\nabla}G_T(x_t,y_t)|$ map (shown in Fig.~\ref{fig:general}(b)), which allows to emphasize short-range variations of $G_T$. Clear resonance patterns can be observed along a set of arc-shaped regions. The arcs are regular, and their centers typically lie in the $G_T=0$ region, i.e.~close to the center of the constriction. Interestingly, the arc structures are completely suppressed within the $G_T=G_0$ plateau.

Transport data in one of the resonance regions is analyzed in finer details in Fig.~\ref{fig:zoom}(a), reporting a scan ($B=6.510$~T and $V_{tip}=-6.0$~V) over the region corresponding to the dashed blue rectangle in Fig.~\ref{fig:general}(a).
We study the evolution of the resonances as a function of applied magnetic field $B$ and tip voltage $V_{tip}$. A careful analysis of these images shows that both an increasing magnetic field and an increasingly more depleting (i.e.~more negative) scanning tip induce a shift of the resonances towards the outer side of the arcs. The effect can be appreciated in Figs.~\ref{fig:zoom}(b) and (c), where the position of the maxima of the conductance gradient
modulus ($|\mathbf{\nabla} G_T(x,y)|$) along the section indicated in Fig.~\ref{fig:zoom}(a) is reported. The distance between two adjacent maxima in $|\mathbf{\nabla} G_T(x,y)|$ corresponds to \textit{half} an oscillation period in $G_T(x,y)$. We can estimate from the roughly linear curves in Fig.~\ref{fig:zoom}(b,c) the typical periodicity of conductance oscillations in the radial space direction ($\lambda_R\simeq 80\pm20$~nm, roughly independent of $B$ and $V_{tip}$), as well as in the magnetic field ($\lambda_B\simeq 73\pm18$~mT) and tip bias ($\lambda_V\simeq 0.64\pm0.16$ V) for fixed tip position. As we shall argue in  Sec.\ref{sec:origin}, oscillations over these scales provide evidence for the presence of backscattering events mediated by localized AD channel states.

\begin{figure}[tb]
\includegraphics[width=0.8\columnwidth]{Fig2a.eps}
\includegraphics[width=0.8\columnwidth]{Fig2b.eps}
\includegraphics[width=0.8\columnwidth]{Fig2c.eps}
\caption{(a) High resolution scan over the region indicated by the dashed blue box in Fig.~\ref{fig:general}(a). $B=6.510$~T and $V_{tip}=-6.0$~V. Measuring a series of such images allows to follow the position of five arcs, labeled  A--E, when either tip voltage $V_{tip}$ or magnetic field $B$ are varied.
(b) Displacements of the A--E peaks as a function of $V_{tip}$ at constant $B= 6.450$~T. Displacements are linear in $V_{tip}$, with approximately the same slope for all five structures, so that their spacing is nearly constant.
(c) The corresponding arc displacements when $B$ is varied, at constant $V_{tip}=-6.0$~V. Also in this case, an approximately  linear dependence is observed.}
\label{fig:zoom}
\end{figure}

\section{Discussion}

\begin{figure}[!ht]
\includegraphics[width=\columnwidth]{Fig3.eps}
\caption{(a) Artificial constriction landscape with disorder and (b) the corresponding edge structure. (c,d) Modified constriction after the addition of a tip depletion: the tip displaces the edge trajectory and activates two backscattering channels through AD2 and AD3. (e) Direct and antidot-mediated regions in the scan and (f) heuristic calculation of the $G_T(x_t,y_t)$ map.}
\label{fig:model1}
\end{figure}

In order to analyze the mechanism behind the tip-induced backscattering, we have performed a numerical calculation of $G_T$ as a function of the tip position $(x_t,y_t)$ on the basis of a realistic potential landscape $U(x,y)$ (defined on a lattice with $10\,{\rm~nm}$ stepsize) which models a single-edge constriction and the disorder present in the sample (we neglect the spin degree of freedom). While a specific potential landscape is considered in the following, qualitatively similar results could be obtained regardless of the precise form of $U(x,y)$. The constriction is obtained by imposing a potential step of $25\,{\rm meV}$, sufficient to strongly deplete the 2DES. The electron liquid is characterized by a Landau level gap $\Delta\approx 10\,{\rm meV}$, of the order of the cyclotron gap in our experiment $\hbar\omega_c=11.1\,{\rm meV}$.
The 25~meV potential step is imposed under the gate position plus a smooth depletion corona of $\approx 200$~nm, while impurities are introduced by adding a set of random Gaussian bumps. A specific realization $U(x,y)$ is reported in Fig.~\ref{fig:model1}(a), with three prominent potential hills (light blue) on the center-left side of the constriction which is defined by the top and bottom gate barriers (orange-red). The guiding center (trajectory) of the first QH edge channel in the constriction is obtained as the equipotential line satisfying $U(x,y)=\Delta/2$, where we take the constriction saddle point as the energy origin. Such an edge-extraction procedure is expected to hold exactly in the $B\to\infty$ limit and to be reasonable as long as the relevant edge features are not much smaller than the magnetic length $\ell=\sqrt{\hbar/eB}$ ($\approx10\,{\rm~nm}$ in our experimental configuration). The edge channel trajectory for the potential of Fig.~\ref{fig:model1}(a) is presented in Fig.~\ref{fig:model1}(b) and shows the formation of three localized edge states around the locally-depleted regions denoted by AD1, AD2, and AD3.

The influence of the SGM tip is introduced in the calculations as an additional moveable depletion spot located at the tip position $(x_t,y_t)$ and described by the potential
\begin{equation}
U_{tip}(x,y)=A_t\exp\left[-\Delta r^2/2R_t^2\right] \; ,
\label{eq:potential}
\end{equation}
where $\Delta r^2=(x-x_t)^2+(y-y_t)^2$, $R_t=200\,{\rm~nm}$, and $A_t=$25~meV, i.e.~inducing complete depletion under the tip. It is worth to remark that in the experiment we rather observe the effect of the tail of the Gaussian potential of Eq.~(\ref{eq:potential}). Therefore the energy dispersion and the electron density slope at the edge of the constriction are considerably smoother when the tip is present. The resulting potential landscape for $x_t=-650\,{\rm~nm}$ and $y_t=-300\,{\rm~nm}$ is shown in Fig.~\ref{fig:model1}(c) and the corresponding edge channel trajectories in Fig.~\ref{fig:model1}(d). As one could naively expect, the tip is able to deviate the edge trajectory and induce new backscattering channels between the top edge (TE) and bottom edge (BE) of the constriction. The determination of the backscattering amplitude for a given edge configuration is highly non-trivial and will depend on a host of (often unknown) experimental parameters, including edge chirality, drift velocity, dephasing times, edge potentials and energy distributions, etc. Despite this, the key features of the conductance scan $G_T(x,y)$ can be captured by taking into account the {\it direct} backscattering BE $\to$ TE as well as the {\it indirect} backscattering involving a {\it single} antidot, for instance BE $\to$ AD2 $\to$ TE. For sake of simplicity and because of their expected smaller and less likely contribution, we instead neglected backscattering channels involving conduction through multiple AD in series.

As the simplest possible model we assume that an appreciable backscattering takes only place if the minimum distance between two edge trajectories is less than a cut-off interaction length of 200~nm. Above this limit the two edges will be considered as non--interacting. In Fig.~\ref{fig:model1}(e) we show how this criterion, combined with the calculated guiding centers of the edge channels described above, allows us to determine the regions where $G_T=0$ (in blue, when the two edges are fully backscattered), those where $G_T=\nu_be^2/h$ (in white, when no interaction is present) and the regions where a \textit{direct} (red) or \textit{AD-mediated} (green shade) backscattering is taking place. The presence of arcs in the conductance map can now be accounted for by assuming that each AD-mediated backscattering path gives rise to conductance oscillations arising from either coherent interference between multiple trajectories carrying different Aharonov-Bohm (AB) phases or charging effects. Both mechanisms can give rise to oscillations, as discussed in Sec.~\ref{sec:origin}.  The resulting plot is shown in   Fig.~\ref{fig:model1}(f) -- see Appendix for details. By comparing it with Fig.~\ref{fig:general}(b) one realizes that our basic model reproduces the general features, including the arc structures, of the SGM scan. In this respect it is important to note that while the chosen phenomenological cut-off length of 200~nm can be considered to be arbitrary, the general features of the simulations are largely independent of the precise value of this parameter.

In particular, the simulation shows that the centers of the arcs observed in the simulated SGM maps are approximately located on an AD, as indicated by the overlays in Fig.~\ref{fig:model1}(d) and in the shaded regions in Fig.~\ref{fig:model1}(e). This proves the origin of the resonances and allows us to map the location of the most relevant scattering centers in our constriction potential landscape. As clarified by Fig.~\ref{fig:model1}(d), when moving along these arcs the tip mediates a suitable interaction between the TE and the BE channels through the AD, thus activating the backscattering mechanism. Different tip positions along the arc correspond to similar values of AD-mediated edge channels coupling. On the contrary, moving the tip towards or away from the AD will modify the coupling and also the area of the AD because of the long range depletion due to the tail of the tip potential. This can also be seen from a comparison between Fig.~\ref{fig:model1}(b) and Fig.~\ref{fig:model1}(d): when the tip is moved towards an AD, for example AD3, the AD area increases.

\section{Interpretation of the conductance oscillations}\label{sec:origin}

In this section we focus on the origin of the observed resonances. Backscattering paths mediated by a localized electronic state are generally expected to give rise to periodic oscillations as a function of the magnetic field as well as of tip position and bias. As argued in the following, this kind of transport features can originate from two different phenomena: coherent interference between different paths encircling, for instance, an isolated antidot (AB effect); and/or CB oscillations due to charging effects. The exact prevalence of one or the other regime can be subtle, and different working conditions have been identified and studied in recent experiments depending on sample details.\cite{Kou,Markus2} Our experimental data do not allow to conclusively rule out one of the two possibilities, therefore in the following we propose a basic model which considers both effects.

The AB effect can be accounted for by assuming that an AD is induced by a relatively smooth potential fluctuation, and that the tunneling region between the AD and an edge channel could be locally approximated by a parabolic saddle-point potential.\cite{TunnelingJain} In this case the conductance $G_T$ due to backscattering through the AD can be written as:\cite{TunnelingKivelson}
\begin{equation}
\Delta G_T=\frac{e^2}{h}{\left| \frac{t_{1}t_{2}}{1-r_{1}r_{2}\exp({i\varphi})} \right|}^2 \; ,
\label{eq:DGt}
\end{equation}
\noindent where $t_1$ and $t_2$ are the transmission amplitudes between BE and AD, and between AD and TE, respectively, while
(under the assumption of a symmetric tunneling region) $r_j=i\sqrt{1-|t_j|^2}$ ($j=1,2$) are the relative reflection amplitudes.\cite{Zeilinger} Here $\varphi=2\pi\Phi/\Phi_0$ is the AB phase corresponding to the magnetic flux $\Phi=B\Omega$ enclosed within the effective area $\Omega$ of the AD, $\Phi_0=h/e\approx4.13\,{\rm mT \cdot\mu m^2}$ being the flux quantum, giving rise to oscillations with a period $\Delta B = \Phi_0/\Omega$.

Coulomb interaction can influence this simple picture in a number of different ways, in particular due to the emergence of compressible stripes and to the importance of CB phenomena in the transport. Charging effects in the localized backscattering center can be taken into account (see Refs.~\onlinecite{Sim,Kataoka,Ford,Ihn,HoutenBenakkerStaring}) resulting in a periodicity of the Aharonov-Bohm oscillations which is rescaled with a factor $\alpha$, i.e.~$\Delta{B}\rightarrow \Delta{B}/\alpha$. Taking $\alpha$ as the number of edge states around the AD (see Refs.~\onlinecite{FordSimp,Goldman,KataokaFord,Kou}), in our case we have $\alpha=2$. The observed periodicity of $73\pm18$ mT would thus correspond, in the case of pure AB oscillations, to an area of $0.057\pm0.014\,{\rm \mu m^2}$ or a circle with radius $R\approx 140$~nm. Considering also charging effects, we obtain an area of $0.113\pm0.028\,{\rm \mu m^2}$ or a circle with radius $R\approx 190$~nm.

Important information about the backscattering center can also be gained by comparing the field, tip position, and tip bias dependence of the oscillations. For instance, looking at data reported in Fig.~\ref{fig:zoom}, one can deduce that an increase of $B$ by $\lambda_B$ can be compensated by moving the tip away from the AD by a radial distance $\lambda_R$. This allows to determine how the tip radial distance $x$ influences the effective area of the scatterer ($d\Omega/dx = -\Phi_0/B\lambda_R<0$) and to conclude that the backscattering center becomes larger when the tip depletion spot is moved towards it, i.e.~that it behaves consistently with our AD interpretation.

\begin{figure}[tb]
\includegraphics[width=\columnwidth]{Fig4.eps}
\caption{(Color online) Conductance oscillations (solid line) and tunneling amplitudes $|t_1 t_2|^2$ (dashed line) as a function of the tip position $x$, assuming $\eta\simeq 10^{-5}$~nm$^{-2}$ and $|t_2|^2=0.75$.}
\label{fig:oscill}
\end{figure}

Figure~\ref{fig:oscill} shows a simple prediction for the AD conductance as a function of the tip position $x$ assuming \cite{TunnelingJain} a dependence of the tunneling parameter $t_1 = e^{-\eta x^2}$ with a heuristic $\eta=10^{-5}\,{\rm nm^{-2}}$ (on the side of the edge deflected by the tip potential) and a constant $|t_2|^2=0.75$ (describing the coupling to the opposite edge which is to first order unperturbed by the tip).
It should be noted that we have not taken into account any effect of temperature, dephasing, or self-averaging mechanisms, which are expected to further broaden the shape of the resonances as well as to reduce the visibility of the oscillations. This can explain the fact that only few small-amplitude oscillations are detectable in our experiment.

\section{Conclusions}

We performed a SGM mapping of the transmission through a wide constriction at filling factor $\nu_b=2$. Beyond observing the usual selective backscattering of the two spin-split edges, we consistently observe arc-shaped resonances with a typical limited extension in both the radial and angular direction. We demonstrate based on a numerical model that such transport features are consistent with AD-mediated backscattering events modulated by the tip potential tails and activated by the deflection of the edge trajectory caused by the core of the tip potential. The evolution of the resonances as a function of the SGM position, the tip potential, and the magnetic field was analyzed also taking into account effects due to the AB phase and charging effects. Our imaging technique allows to reconstruct prominent features of the confinement potential including the effective profile of the gate and the approximate position of strong scattering centers which are nucleating AD structures in the measurement.

\begin{acknowledgments}
We acknowledge financial support from the Italian Ministry of Research (MIUR-FIRB projects RBIN045MNB and RBID08B3FM). We gratefully acknowledge discussions with Rosario Fazio.
\end{acknowledgments}


\section{Introduction}

Asymmetries in spectral profiles of Stokes $V$ are commonly measured on the solar surface, from the penumbrae of sunspots to the quiet Sun; they have been known for decades \citep[e.g.][]{sami_84, grossmann_96, valentin_97, sigwarth_99, sigwarth_01, rezaei_steiner_07, marian_08, bartolomeo_11}. The generally accepted model to produce asymmetric Stokes $V$ profiles is by having anticorrelated velocity and magnetic field gradients along the line of sight \citep{illing_75, jorge_89, arturo_02}. Therefore, circular polarization asymmetries encode the structuring of magnetic fields along the line of sight. However, the limited spatial resolution of most observations mixes the contributions of different regions to the asymmetries, hence, allowing only indirect inferences of gradients along the line of sight.

Outside of sunspots the classical model to produce asymmetries is at the canopies of the magnetic feaures expanding with height \citep{grossmann_88, sami_89}. This model has been able to reproduce the simultaneous
presence of Stokes $V$ area asymmetry, the absence of a Stokes $V$ zero-crossing shift \citep{sami_86, valentin_97}, and the change in sign of the area asymmetry near the solar limb \citep{bunte_93} in obervations of the quiet Sun and active region plage taken at rather low spatial resolution \citep{stenflo_solanki_87, valentin_97}. Alternatives have also been suggested, such as the presence of bulk velocities in the flux tube \citep{luis_97, luis_00,sigwarth_01}, or the MISMA and MESMA models \citep{jorge_96, thorsten_markus_07}.

This letter is devoted to employing Stokes $V$ asymmetry in order to diagnose the internal structure of network magnetic features, i.e., the stongest flux accumulations in the quiet Sun. It is known, since the work of \cite{stenflo_73}, that network patches are composed of magnetic elements with kG fields. Later, two very different models have been proposed to explain the observed asymmetries. First, the thin flux tube model that assumed a rather vertical magnetic tube expanding with height. Second, the MISMA hypothesis that assumes an infinite number of magnetic atmospheres in the resolution element. In both cases, the cause for the asymmetries is the correlation between the distribution of magnetic fields and velocities. However, the MISMA model is a one-dimensional model, having no predictive power for the spatial variation of asymmetries in magnetic structures. In order to understand how magnetic fields are really organized and to discern between these two models (or create new ones), we must observe at higher spatial resolutions, which will allow us to directly map the substructure of magnetic features and to trace gradients along the optical depth.

\section{Data and analysis}

We map the Stokes $V$ profile asymmetries in the quiet Sun using data obtained by IMaX \citep{valentin_11} on Sunrise \citep{sami_10, barthol_11}. The data analyzed here consist of Stokes $I$ and $V$ profiles of the Fe\,{\sc i} 5250.2 \AA\ line (Land\'e factor of 3) sampled at 11 wavelegth positions that were recorded on June 10th at the disk center. The effective field of view (FOV) is 46.8$''\times$46.8$''$  and the spatial resolution is about 0.2$''$. The noise level in Stokes $V$ is $10^{-3}$ in units of the continuum intensity, $I_c$, prior to the application
of the phase diversity reconstruction technique. After this procedure, the noise of the reconstructed data is spatially and temporaly correlated, but still follows a Gaussian distribution, with a standard deviation of $2.5 \times 10^{-3}$ $I_c$.

For the analysis presented in this letter, we select a small region of the FOV containing the largest and the most intense circular polarization patch (see Fig. \ref{mapas_region}). It has a magnetic flux of $\sim 2\times 10^{19}$ Mx contained in an area of $\sim 5.5$ Mm$^2$. The pattern of the granulation is distorted within the magnetic patch, which is typical of kG field accumulations. In order to minimize the effect of the noise in our analysis, we select only those profiles with amplitudes larger than 5 times the noise level. 46.6\% of the non-reconstructed, and 29.3\% of the reconstructed profiles belonging to the small selected area meet this criterion. The selected Stokes $V$ profiles in this region are mostly regular profiles. More precisely, 93.6\% and 88.9\% of the profiles have two lobes in the non-reconstructed and in the phase diversity reconstructed data, respectively. We define the area ($\delta A$) and amplitude ($\delta a$) asymmetries of regular profiles as:

\begin{figure}[!t]
\centering{
\includegraphics[width=0.25\textwidth, bb=173 434 500 686]{cont.ps}
\hspace{-1.5cm}
\includegraphics[width=0.25\textwidth, bb=173 434 500 686]{cont_rec.ps}\\
\includegraphics[width=0.25\textwidth, bb=173 434 500 686]{ampl_v.ps}
\hspace{-1.5cm}
\includegraphics[width=0.25\textwidth, bb=173 434 500 686]{ampl_v_rec.ps}\\
\includegraphics[width=0.25\textwidth, bb=173 434 500 686]{anchura_int.ps}
\hspace{-1.5cm}
\includegraphics[width=0.25\textwidth, bb=173 434 500 686]{anchura_int_rec.ps}}
\caption{From top to bottom, continuum intensity, circular polarization amplitude and full width at half maximum of the intensity profile for the stongest flux region in the IMaX field of view. Left panels display the observables of the original data, while the right ones show the same scenes after phase-diversity reconstruction. The white contours represent the isocontours of log$[A_V/I_c]=-2.3,-1.9$ for the non recostructed and the reconstructed data, respectively. The symbol $A_V$ represents the amplitude of the Stokes $V$ profile.}
\label{mapas_region}
\end{figure}

\begin{eqnarray}
\delta A &=&\sigma \frac{\int_{\lambda_i}^{\lambda_f} V(\lambda)\mathrm{d}\lambda}{\int_{\lambda_i}^{\lambda_f} \mid V(\lambda)\mid \mathrm{d}\lambda} \\
\delta a &=& \frac{a_b-a_r}{a_b+a_r}
\label{eq}
\end{eqnarray}
The symbols $\lambda_i=5250.0164$ \AA\ and $\lambda_f=5250.4014$ \AA\ are the
initial and final wavelengths of the integral over the observed profile. The
symbol $\sigma$ is the sign of the Stokes $V$ blue lobe. The $a_b$ and $a_r$ symbols are
the unsigned blue and red lobe amplitudes, respectively. The non-linearity of
the equations to compute asymmetries make it difficult to evaluate the
relevance of instrument errors. The first error source to take into account is
the photon noise. We have verified (simulating Stokes $V$ profiles as derivatives of a gaussian that match IMaX intensity data and including gaussian noise)
that the photon noise does not introduce a
bias in the calculation, only spreads the results, the relative errors in the
area and amplitude asymmetries being at most $\pm 15$\%.

Calibrations performed at the launch station in Kiruna, proved that there
exist considerable amounts of crosstalk between linear and circular polarization
introduced by the ISLiD \citep{gandorfer_11}. The relative amounts of Stokes $Q$ and Stokes $U$
change with time due to image rotation and cannot be predicted for a specific data set.
When IMaX observes in a vectorial mode, all four Stokes parameters
are measured and multiplication by the inverse of the Mueller matrix takes
into account all crosstalk sources (leaving as unknown only the zero of
the azimuth angle). In the longitudinal observing mode used here \citep[L12-2,
refer to][]{valentin_11} only Stokes $I$ and Stokes $V$ are measured making a
proper correction of crosstalk impossible. However, we expect the impact of this
crosstalk to be negligible for the investigation performed in this work that
concentrates on a relatively large network patch. First of all, the fields
in these regions are known to be predominantly vertical and all horizontal
fields observed with IMaX in other full vector maps show that the
linear polarization signals mainly corresponded to internetwork elements
\citep{danilovic_10}. Second,
inspection of similar network patches observed in the IMaX vector mode V5-6
\citep[bottom part of Figure 3 in][]{sami_10}, show that they never display
linear polarization signals above 3 sigma levels at the network concentration, neither
inside it nor at the boundaries, where a possible canopy could
generate some linear polarization. We thus conclude
that in the region under study, the known crosstalk in the longitudinal mode of
IMaX can be ignored even though it cannot be accounted for as is done
in the IMaX vector modes.

Finally, the effect of the finite spectral (and spatial) resolution, as well as
the effect of the phase-diversity reconstruction, is studied in
\cite{andres_11} for the special case of MHD simulations. These authors
conclude that the spatial degradation tends to understimate both amplitude and
area asymmetries. Degrading the spectral resolution, however, tends to
oveerstimate the area asymmetry while it does not affect strongly the amplitude
asymmetries, confirming earlier results of \cite{sami_stenflo_86}. After
applying the reconstruction algorithm, the computed asymmetries resemble the
original ones, i.e., for MHD simulations at full spectral and spatial
resolution. However, we must note that these conclusions depend on the actual
shape of the profiles. The work of \cite{sami_86} predicts (for the IMaX case)
an overstimate of the area asymmetries by a factor $\sim 1.6$ and an
understimate of the amplitude asymmetries by a factor $\sim 0.8$. Consequently,
strongly asymmetric/antisymmetric profiles observed by IMaX do correspond to
profiles with large/small asymmetries at full spectral resolution, although the
exact values may differ.

\begin{figure*}[!t]
\centering{
\includegraphics[width=0.33\textwidth, bb=127 406 472 700]{mapa_asim_area_region_kg.ps}
\hspace{-1cm}
\includegraphics[width=0.33\textwidth, bb=127 406 472 700]{mapa_asim_area_region_kg_rec.ps}\\
\includegraphics[width=0.33\textwidth, bb=127 406 472 700]{mapa_asim_ampl_region_kg.ps}
\hspace{-1cm}
\includegraphics[width=0.33\textwidth, bb=127 406 472 700]{mapa_asim_ampl_region_kg_rec.ps}}
\caption{Top panels represent the area asymmetry of the circular polarization profiles of the selected strong flux patch. Bottom panels display the amplitude asymmetry. Right (left) panels present the results showing the reconstructed (non-reconstructed) data. The areas in white represent those pixels where the signal to noise ratio of the Stokes $V$ profiles is lower than 5. Grey arrows are the locations of the profiles displayed in Fig. \ref{perfiles}.}
\label{mapas_asim}
\end{figure*}

\section{Circular polarization asymmetries and the magnetic field geometry}

Figure \ref{mapas_asim} shows maps of the area and amplitude asymmetries of the observed Stokes $V$ profiles in the selected patch. The left panels correspond to the non-reconstructed data and the right  ones to the phase-diversity reconstructed ones. All four panels display the same pattern: the smaller asymmetries tend to be located in the inner parts of the structure, while the larger ones lie at the borders, close to where the signal drops below the noise level. The larger asymmetries situated at the borders of the structure have dominantly positive values, but negative values are also found (e.g. along the top edge of the feature). The spatial pattern for the two types of asymmetries is very similar, the correlation coefficient between amplitude and area asymmetries being 0.9.

Figure \ref{perfiles} displays three sets of Stokes $I$ and $V$ profiles, extracted from the positions marked with small grey arrows in Fig. \ref{mapas_asim}. The Stokes $V$ profile displayed in the left panels ($\delta A=0.5$\% and $\delta a=4$\%) is representative of the nearly antisymmetric profiles found in the central part of the patch. The profile in the middle panels has larger area and amplitude asymmetries (38\% and 33\%, respectively), as found closer to the boundary of the magnetic feature. Finally, the profile in the rightmost panels is a one-lobed profile, with extreme asymmetries: 81\% area asymmetry, and 80\% amplitude asymmetry. Such profiles are found almost exclusively at the very edge of the feature.

In order to check if the expanding flux-tube model applies to our observations,
we use the SIRJUMP code, which is a modification of the SIRGAUS code
described in \cite{sirgaus}. The two codes follow the inversion
strategy of SIR \citep{basilio_92}. SIRJUMP has recently been
employed to interpret single-lobed Stokes $V$ profiles in the quiet Sun by \cite{alberto_12}.
With SIRJUMP we retrieve the parameters of a model
atmosphere showing  discontinuous stratifications of some of the parameters
along the line of sight. The position of the jump (i.e., the discontinuity) is a free
parameter in the inversion scheme. The code can put it below optical
depth unity, which would indicate that no discontinuity is needed to explain the
observed profiles.  Above and below the jump, the magnetic field and the line-of-sight
velocity are assumed to be constant with optical depth.
The temperature stratification is modified with respect to an initial guess by the inversion algorithm in three nodes and
interpolated between them. The microturbulent velocity and
the inclination of the magnetic field vector to the line of sight
are height independent. Since we do not have the information of
the linear polarization, we cannot retrieve the  magnetic field azimuth, and can barely infer the inclination.
Therefore, we perform two different inversions, one forcing the field to be aligned with the line-of-sight
and another one leaving the inclination as a free parameter. The macroturbulent
velocity is fixed to 2.06 km s$^{-1}$, which matches the broadening
induced by the Point Spread Function of the IMaX instrument.

Figure \ref{mapas_inv} displays the results of the two different inversions. Top left panels show the optical depth at which the discontinuity occurs for the inversion with the inclination as a free parameter (panel 4a) and the inversion forcing the inclination to 0$^\circ$ (panel 4b). In both cases, the results are very similar, i.e., they do not depend much on the assumed geometry of the magnetic field. For log$\tau_{5000}  \gtrsim 0.5$ the contribution to the spectral line is small and the discontinuity is not well constrained. In spite of the noise some trends can be seen. In the central part of the structure we retrieve a magnetic field with no clear discontinuity. As we move from the center towards the edges of the structure, the discontinuity appears at increasingly higher layers. Therefore, the inversions made pixel by pixel have naturally retrieved the expansion of a magnetic flux concetration with height in the solar atmosphere. At all heights (below and above the discontinuity), the pattern of the velocity is qualitatively similar, following the local distorted granulation.

In the lower panels of the same figure we display the relative differences of the longitudinal components of the magnetic field (panels e and f) and the velocity (panels g and h) above and below the discontinuity. Note that, since we do not have linear polarization, the only trustable quantity is the longitudinal component of the magnetic field. We have used the following definitions:

\begin{eqnarray}
\delta B_{long} &=& \frac{|B_{long}^{a}| - |B_{long}^{b}|}{|B_{long}^{a}| + |B_{long}^{b}|} \nonumber \\
\delta v_{long} &=& \frac{|v_{long}^{a}| - |v_{long}^{b}|}{|v_{long}^{a}| + |v_{long}^{b}|},
 \end{eqnarray}
where the superindex $a$ and $b$ refer to the quantities above and below the discontinuity, respectively.
From these figures, we can see that both the longitudinal magnetic field and the line-of-sight velocity
are almost constant with optical depth at the core of the structure, coinciding with strong kG fields. As we move towards the borders of the structure (where the longitudinal component of the magnetic field is of the order of 100 G), the relative variation of both quantities increases. Note that the gradient of the longitudinal magnetic field changes sign at the regions where the area and amplitude asymmetries were negative.

\begin{figure}[!t]
\includegraphics[width=0.5\textwidth]{perfiles_canopy.ps}
\caption{Examples of Stokes $V$ profiles. The left panels display a respresentative of the relatively antisymmetric profiles found in the center of the structure, middle panels the more asymmetric ones found near the borders of the patch, while the right panels depict a one-lobed profile at the edge of the magnetic patch. The exact locations are marked in Fig. \ref{mapas_asim} by small grey arrows. The observed profiles are represented by the rhombs and the best fit using an inversion algorithm based on a discontinous model along the optical depth (see the text) is displayed with a blue line. The inversion with the inclination forced to zero is represented by the green lines and the inversion with the inclination as a free parameter is displayed with the blue lines.}
\label{perfiles}
\end{figure}

\begin{figure*}[!t]
\centering{
\includegraphics[width=0.28\textwidth, bb=80 50 260 340]{salto_canopy_cociente_campo.eps}
\hspace{-1.65cm}
\includegraphics[width=0.28\textwidth, bb=80 50 260 340]{salto_canopy_cociente_campo_gamma0.eps}
\hspace{0.5cm}
\includegraphics[width=0.28\textwidth, bb=80 50 260 340]{temp_logtau-1.5_cociente_vel.eps}
\hspace{-1.65cm}
\includegraphics[width=0.28\textwidth, bb=80 50 260 340]{temp_logtau-1.5_cociente_vel_gamma0.eps}}
\caption{Results of the two inversions of the original data based on a model that has a discontinuity along the optical axis. Panels a, c, e, and g display the inferred quantities for the inversion in which the inclination of the magnetic field is a free parameter. Panels b, d, f, and h represent the retrieved parameters from the inversion in which we have fixed the inclination of the magnetic field with respect to the line-of-sight to 0$^\circ$. Similar results can be derived from the phase diversity reconstructed data. Panels a-d show the optical depth of the discontinuity (left) and the temperature at log${\tau_{5000}}=-1.5$ (right). Panels e-h display the relative differences of the line-of-sight magnetic field (left) and velocity (right) above and below the discontinuity (see Eq. 3)}
\label{mapas_inv}
\end{figure*}

\section{Conclusions}

We have computed the area and amplitude asymmetries of the Stokes $V$ profiles of the Sunrise/IMaX dataset with the best spectral sampling (35 m\AA). We selected the strongest magnetic flux accumulation, which belongs to the network pattern, in the IMaX FOV. The core of the structure is, on average, characterized by an area asymmetry of 0.2\% and an amplitude asymmetry of 4.4\% (with an upper limit for the error of $\pm 1.3$\% for both averaged area and amplitude asymmetries). In the outermost parts of the network patch, the Stokes $V$ profiles have the largest area and amplitude asymmetries, even aproaching 100\% (one-lobed profiles). The fact that we are seeing a spatial variation
of Stokes $V$ asymmetries means that, for the first time, we are resolving the internal structure of a magnetic patch belonging to the photospheric network. The non negligible amplitude asymmetry inside the magnetic feature suggests that unresolved motions are still present (definitely gradients along the line of sight and probably horizontal gradients).

We inverted our data using the SIRJUMP code, which models the gradients of magnetic field and line-of-sight velocity with an abrupt discontinuity at a certain optical depth. From two different assumptions on the inclination of the field, we inferred that the observed network patch is made of resolved kG fields which is consistent with the work of \cite{lagg_10}. Moreover, the inversions naturally retrieve a geometry corresponding to a magnetic field expanding with height.

The spatial distribution of the Stokes $V$ asymmetries and the inversions applied to these data suggest a geometry similar to the classical model of magnetic features expanding with height \citep{spruit_76, sami_93} that takes pressure balance and magnetic flux conservation with height into account and is consistent with three-dimensional magneto-hydrodynamical simulations \citep{lotfi_09}. In such a model, the central core of the structure traces the body of the flux tube, containing the vertical kG fields, while the external parts correspond to the canopy, i.e. expanding, inclined field overlying the quiet photosphere (this expansion is caused by the need to conserve magnetic flux for decreasing gas pressure with height).

The standard model for the creation of asymmetries in the presence of such a canopy is based on the presence
of a discontinuity in both the magnetic field strength and the line-of-sight velocity at the base of the canopy (i.e. abrupt gradients in such quantities along the optical depth). In the past, it was assumed that there is no velocity flow above the canopy's base but there is a downflow below it, with the magnetic field being present only above the canopy base. This geometry produced a positive area (and amplitude) asymmetry \citep{grossmann_88, sami_89}, with nearly one-lobed profiles in the outer parts of the flux tubes \citep{sami_89} in agreement with the observations.

Our results for the asymmetries (see Fig. \ref{mapas_asim}) are in
agreement with this general picture. However, there are some locations along the edge of our feature where the area (and amplitude) asymmetry is negative, indicating that either the magnetic field or the velocity field gradient must change sign there. We know that the kG network fields are embedded in the more disorganised weakly magnetized quiet photosphere. Hence, we would expect a smooth transition between the organized kG fields to the
tangled weak photospheric fields (note that we have selected the profiles above the noise level, i.e., the borders of the structure depend on the amplitude recquired for the profiles). If those regions are mostly tracing the quiet photosphere (note that the height of the discontinuity is located at lower layers and that the temperature is colder than that of the central part of the structure), it is not surprinsing to find gradients of magnetic field and velocity that change sign as compared to the classical canopy overlying a non-magnetized photosphere.

The observations presented in this letter support the interpretation of the network as kG fields expanding with height. In contrast to what was assumed in the past, we find velocity patterns above and below the canopy that follow the local distorted granulation. To strengthen the model for strong flux concentrations as the one studied in this letter, in the future we need to include in our analysis the constraints introduced by the linear polarization and, more importantly, the simultaneous information of many spectral lines to trace the canopy at different heights.

\begin{acknowledgements}

We specially thank Manuel Collados Vera for very helpful discussions and for carefully reading the manuscript. The Spanish contribution has been funded by the Spanish Ministry of Science and Innovation under projects ESP2006-13030-C06, AYA2009-14105-C06 (including European FEDER funds), and  AYA2010-18029 (Solar Magnetism and Astrophysical Spectropolarimetry). The German contribution to Sunrise is funded by the Bundesministerium f\"ur Wirtschaft und Technologie through Deutsches Zentrum f\"ur Luft- und Raumfahrt e.V. (DLR), Grant No. 50 OU 0401, and by the Innovation fond of the President of the Max Planck Society (MPG). This work has been partially supported by the WCU grant (No. R31-10016) funded by the Korean Ministry of Education, Science and Technology.

\end{acknowledgements}