Nuclear bars and blue nuclei within barred spiral galaxies

Abstract: Multicolour near IR photometry for a sample of 32 large barred spiral galaxies is presented. By applying ellipse fitting techniques, we identify significant isophote twists with respect to the primary bar axis in the nuclear regions of $\sim$70 \%\ of the sample. These twists are identified in galaxies as late as SBbc and are clearly distinguishable from spiral arm morphology. At most seven of the galaxies with isophote twists are inferred to possess secondary (nuclear) bars, the axis ratios of which appear to correlate with morphological type. The remainder may result from triaxial bulges, or from oblate bulges misaligned with the primary bar. The near IR colour distributions in these data show evidence for (red) circumnuclear star forming rings in 4 galaxies. The majority of the sample (19) also possess striking blue nuclear regions, bluer than typical old stellar populations by $\sim$0.3 mag. in (J--H) and $\sim$0.23 mag. in (H--K). Such blue colours do not appear to correlate with the presence of nuclear rings or pseudo--rings, nor with the activity of the host galaxy (as determined from emission--line spectroscopic characteristics). Several mechanisms to explain this blue colour are considered.

\section{Introduction}

The presence of a bar is likely to be an important factor in the fuelling
of nuclear (starburst/AGN) activity in many galaxies ({\sl e.g.}\ Heckman 1980,
Hawarden {\sl et al.}\ 1986, Arsenault 1989). An increasing proportion
of barred galaxies show evidence of starbursts in the form of active star
formation within circumnuclear rings -- CNR's. Good examples of this
phenomenon are NGC1097, 4321 and 5728 (Shaw {\sl et al.}\ 1993 -- henceforth Paper I).
This is readily understood
as a direct result of gas motions within the barred potential:
the bar provides an efficient mechanism to
transport gas from the spiral arms to the inner Lindblad resonance (ILR) --
or the outer ILR if two exist -- where the resulting reservoir of gas provides
a
source for star formation activity.

Recent theoretical studies ({\sl e.g.}\ Shlosman {\sl et al.}\ 1989, Pfenniger \&\ Norman
1990) have also highlighted a mechanism whereby an independent,
rapidly rotating nuclear bar can provide an efficient means of transporting
material from the CNR onto the nucleus, thereby possibly fuelling LINER/Seyfert
activity. Some support for this hypothesis may exist in those barred
Seyfert galaxies
which also possess CNR's ({\sl e.g.}\ NGC1097: Hummel {\sl et al.}\ 1987b; NGC3783: Winge
{\sl et al.}\ 1992;
NGC5728: Schommer {\sl et al.}\ 1988, Wilson {\sl et al.}\ 1993; NGC6951: Boer
\&\ Schultz 1993). Shaw {\sl et al.}\ (Paper I) stress the
significant r\^{o}le any such mass transfer of gas may have on the
stellar dynamics in the nuclei of barred galaxies.

Inevitably, any study of the fuelling of nuclear activity by a bar
requires detailed knowledge of the mass distribution in the central
region. In general
such information is unavailable as existing observational constraints on
the stellar dynamics in the nuclear regions of barred potentials are
lacking. The influence of obscuration by dust in the nuclei of barred
galaxies ({\sl e.g.}\ Baumgart \&\ Peterson 1986) is also a problem, as this places
a fundamental limit on photometric studies of the inner regions. The
problem is particularly acute given the large number of barred galaxies
anticipated to possess nuclear isophote twists suggestive of nuclear bars
or coupled star/gas components (Paper I).  An important benefit of using
near IR photometry is that the effects of dust obscuration are greatly
reduced, leading to a more reliable determination of the morphology
of the underlying old stellar populations, and ultimately the gravitational
potential, in the nuclear regions. No near IR survey of the nuclear regions
of barred spiral galaxies has been undertaken.
It is thus unclear what proportion of barred galaxies possess nuclear bars,
and whether there is indeed the strong correlation between such bars and the
presence of CNR's implied in Paper I. The degree to which nuclear bars can
be distinguished from possible triaxial bulges
({\sl e.g.}\ Kormendy 1979, 1982) is also unknown.

In this paper, we present the results of a near IR photometric study of
32 nearby barred spirals which seeks to address such issues. For the reasons
given above, we are able to extend studies of the nuclear properties of
barred galaxies to significantly later morphological types than has been
possible to date.
The specific objectives in
our study are to determine the proportion of barred galaxies which possess
nuclear bars, CNR's and/or isophote twists in a statistically significant
sample. In searching for CNR's, we make use of the sensitivity of the K--band
(2.2 $\mu$m) to emission from M giants/supergiants residing in regions of
ongoing star formation activity. Analyses of near IR colour distributions are,
therefore, an important aspect of this study. Given the limitations inherent
in the use of rotation curves to define the pattern speed of a bar, the
identification and properties of CNR's take on added importance if, as is
suggested by N--body simulations ({\sl e.g.}\ Combes \&\ Gerin 1985),
they identify the location of one
of the principal resonances in a barred potential.

\section{The sample, data acquisition and reduction}

\subsection{Sample selection}

The data in this paper comprise near IR images of 32 large spiral
galaxies classified as SAB or SB by de Vaucouleurs {\sl et al.}\ (1976). Table 1
collates general properties of these objects.

The present dataset constitutes a subset of our complete sample of large
(D$_{25}$ \gtae\ 2 arcmin), bright (m$_B$ \ltae\ 12 mag.), non--interacting nor
irregular barred galaxies. This sample was constructed primarily from Hummel
{\sl et al.}\ (1987b), Arsenault (1989 -- both control and AGN/STB samples),
Kormendy (1979), Pogge (1989), Buta \&\ Crocker (1993) and Friedli (private
communication). It therefore includes all nearby
barred galaxies which possess nuclear rings,
pseudo--rings, lenses, that show excessive nuclear H$\alpha$ fluxes or
evidence for nuclear activity (starburst, LINER or Seyfert). The size criterion
is chosen so as to improve spatial resolution in the nuclear regions.

The sample presented in this paper
is not complete, as it constitutes the first observations acquired from our
sample. More general discussions, {\sl e.g.}\ relating to the statistical properties
of the sample, will be presented in future papers as and when the observational
data become available. Future data will improve the available spatial
resolution by utilising alternative pixel scales. However,
the number of objects studied in the present paper is sufficient to
constitute a statistically significant subset and thus allows
an exploratory investigation of the frequency of nuclear bars, isophote
twists and
investigating the possible links with CNR's. The subsample does contain an
inherent bias insofar as there was a preferential selection of the
larger (though not necessarily more face--on) barred galaxies. However,
we stress that an important aspect
of this subsample is the broad range of morphological types represented.

\subsection{Data acquisition and reduction}

Simultaneous J, H, and K (1.2 -- 2.2 $\mu$m) photometry of the present
sample was undertaken using the (256x256 element PtSi array) multichannel
IR camera (Ellis {\sl et al.}\ 1992) on
the KPNO 1.3m between 5 and 14 March 1992. The pixel scale in this
configuration is 1.35 arcsec, or 0.10--0.45 kpc/pixel over the distance
range of our sample. Integration times of 3x180 sec
on each object were interspersed by equal length sky frames, typically offset
from the galaxy by $\sim$10 arcmin. Each object frame was also offset
from the other by $\sim$10 arcsec to minimise contamination by defective
pixels. Dark current subtraction was facilitated using equal length dark
frames through the run.

Any contaminating objects were removed from the sky offset images, and the
resulting frames median filtered to create high S/N ``flat--field'' frames
for each galaxy in each colour. Such images were normalised to a mean of
unity and were divided into the corresponding object frame for that
particular filter to remove pixel--to--pixel response variations. Such
flat--fielded frames were then spatially aligned
and summed. Residual edge--effects were removed
by subsetting the images and any remaining defective pixels removed by
interpolation.

The object frames were then calibrated using observations of 10 faint (J
= 6.6--8.4 mag) near IR standards derived from Elias {\sl et al.}\ (1982).
Typically, between 4 and 8
stars were observed each night. Four of the five nights proved to be highly
photometric -- the derived zeropoints being consistent throughout
the night to within $\pm$0.01--0.02 mag. arcsec$^{-2}$ in all passbands.
Data from the single non--photometric night, which affects observations of
NGC4754 alone, were calibrated using observations of the standard star HD129653
observed immediately prior to the object frames. Differential airmass
corrections were applied to the derived zeropoints using mean extinction
coefficients appropriate for the site.
The sky background estimates used in each case were median values derived
from the object frames directly (in the majority of cases where sky
dominates the frames), or from the median filtered sky offset frames (where
the objects fill the frame). In the latter, as no significant sky background
variations were noted between each sky frame acquired, the mean flux
from each of the three frames provides a good estimate of the sky for the
objects
concerned.
Since the J, H
and K images were acquired simultaneously, no correction for relative
atmospheric extinction differences was necessary when deriving the colour
distributions of each object.

The seeing, measured directly from the object frames throughout the run,
was 2.0 ($\pm$0.3) arcsec FWHM at J,
1.9 ($\pm$0.1) arcsec at H and 1.9 ($\pm$0.2) arcsec at K. In consequence,
ellipse--fitting results derived from the innermost $\sim$2 arcsec must be
considered highly unreliable.

\section{Results}

\subsection{Primary bar position angles}

An important first step in searching for evidence of isophote twists in the
nuclei of barred galaxies is an assessment of the orientation of the primary
bar.

In each object, visual estimates derived from
contour plots and greyscale images of each J, H, K image were compared to
those derived from the ellipse fitting
procedures discussed below. This is necessary because the latter technique
can be unduely influenced by spiral arm morphology (section 3.2). The
respective
measurements are listed in Table 2, together with their corresponding
uncertainties. [Absolute N and E were defined from the known N--S
and E--W positional offsets of the standard star observations.]
The quoted errors reflect the uncertainty in a particular
position angle estimate, and the spread in position angles derived from
each passband. Also listed in
this table are angles measured from blue and/or red sky survey plates in those
objects for which the principal bars can be clearly defined.

Of the 32 objects listed, 3 are too edge--on to allow determination of the
primary bar position angle. Such objects will not be discussed further in the
context of identifying possible isophote twists.
Measures for the primary
bar position angles for the remainder
are clearly equivalent between the present data and sky survey plate
estimates in all cases other than NGC3953
and NGC4536. Both galaxies are highly inclined systems (i $>$ 60$^{\circ}$,
90$^{\circ}$ corresponding to edge--on), compromising the evaluation
of the position angle from the sky survey plates. In these cases, the measures
derived from the present data have been adopted.

\subsection{Ellipse fitting}

Initial inspection of the data revealed clear evidence for isophote twists
(with respect to the primary bar position angles) in 16 objects. With a view
to quantifying this, and measuring the
degree of isophote misalignment more explicitly, ellipse fitting techniques
({\sl e.g.}\ Jedrzejewski 1987)
were applied to each (J, H, K) image independently. As a consequence, 21
objects (72\%\ of the non edge--on sample) were found to possess measureable
isophote twists, the magnitudes of which are collated in Table 2. [This
frequency does not include NGC4321 which possesses elongated
inner isophotes and will be discussed further below.] Such a
high incidence of isophote misalignment confirms the expectations in
Paper I: this is a  common property of the central near IR light
distributions of barred galaxies of all
morphological types. Figure 1 illustrates this with the 4 most
striking examples of isophote twists (NGC3941, 3953, 4613 \&\
4754), whilst the results of the ellipse fitting to these particular
objects are shown in Figure 2.

On the assumption that spiral arms are trailing, we have made an assessment
of whether the measured isophote twists trail or lead the principal bar in
a similar fashion to that undertaken by Buta \&\ Crocker (1993). Our
assessment is illustrated by the sign of the twists in Table 2. The magnitude
of the twists, and spatial regions over which they occur, are broadly
consistent between each passband (the differences being reflected in the
associated errors in this table). Only in the outer regions of some objects
are discordant ellipse fitting results derived. In part, these are due to
intrinsic colour variations within the disc components (manifest in the case
of NGC4321 -- Figure 4), but are also a reflection of the impact of the
declining quantum efficiency variations of the IR array with wavelength
(falling from 6.6\%\ at J to 3.4\%\ at K). The limited spatial coverage of the
disc component may be an additional factor for the larger galaxies.

Undoubtedly,
NGC4321 and NGC4274 present the clearest evidence of a nuclear bar within the
present
sample as Figures 3 and 4 clearly show. The nuclear bar in NGC4321 extends to a
radius $\sim$9 arcsec (1.3 kpc) -- witness the trough in ellipticity, and peak
in $\cos 4\theta$, profiles at this point (Figure 4).
This bar possesses an axis
ratio similar to that of the main galaxy at large radius and is aligned with
the
inferred primary bar major axis to within 7$^{\circ}$
(adopting a position angle from Arsenault
{\sl et al.}\ 1988 and references therein). It is also immediately interior
to the CNR evident in
our data (Section 3.3) and in emission--line images (Arsenault {\sl et al.}\ 1988).
Thus, in some respects the properties of the central $\sim$1.5 kpc of
NGC4321 bear a striking similarity to equivalent regions in NGC1097 and
NGC5728 (Paper I). Unlike those objects, however, the inner bars in NGC4321
and NGC 4274 are closely aligned with the primary bar.

A correlation exists between the scale of the measured
isophote twists ($\delta$R$_{max}$) and bar length -- larger bars
displaying isophote twists over larger radii (the derived product--moment
correlation
coefficient of 0.68 is significant at $>$ 99\% significance level). An
equally significant correlation between $\delta$R$_{max}$ and morphological
type results from
the fact that bars in the later type galaxies in this sample tend to be
larger. [Note that the anti--correlation between bar and bulge dimensions
found by Athanassoula \&\ Martinet (1980) was restricted to objects earlier
than Sa.]

Interestingly,
the nuclear bar axis ratios in NGC1097, NGC4274, NGC4321 and NGC5728
are markedly smaller than
those identified in early type galaxies ({\sl e.g.}\ Jarvis {\sl et al.}\ 1988),
suggesting a possible dependence on
morphological type or, equivalently, the nature of the primary bar.

\subsection{Colour distributions}

A striking feature of the derived colour maps is that 19 of the 32 galaxies
show particularly blue nuclei, {\sl i.e.}\ regions of the colour maps which
are considerably bluer than those areas
of each galaxy dominated by the bulge component. We refer to these objects
as comprising the ``blue nuclei'' sample. Examples of these structures
are given in Figure 5 (b--e). We quantify the
magnitude of these colours in Table 3. [The errors quoted throughout
this section correspond to standard errors on the mean values for a particular
dataset.]

As is evident in this table,
a continuum of colours exist between these ``blue nuclei'' sample
and the remainder
(henceforth referred to as the ``control'' sample). The
distinguishing property of the former sample is the
identification of substructure within the colour maps.
Typical nuclear colours within a 0.38 ($\pm$0.02) kpc radius aperture are
0.45 ($\pm$0.06) mag. in (J--H) and 0.02 ($\pm$0.04) mag. in (H--K) for the
``blue nuclei'' sample.
The colours of the ``control'' galaxies
within equivalent apertures are 0.59 ($\pm$0.06) mag.
and 0.19 ($\pm$0.04) mag. respectively. Interestingly, even the ``control''
sample
colours are somewhat bluer than those of typical spiral galaxies:
the median colours ({\sl i.e.}\ averaged
over all morphological types) from Griersmith {\sl et al.}\ (1982) are 0.71--0.76
mag. in (J--H) and 0.25--0.30 mag. in (H--K). Similarly, within 6 arcsec
diameter apertures, data from Forbes {\sl et al.}\ (1992) yield typical colours of
(J--H) = 0.83 ($\pm$0.14) mag. and (H--K) = 0.44 ($\pm$0.16) mag. for a sample
of 15 galaxies, only two of which are barred. An IR colour--colour diagram
for the galaxies listed in Table 3 is shown in Figure 6, where
each galaxy has been identified according to its spectroscopic (nuclear
activity)
characteristics as presented in Table 1.

The colour maps ({\sl e.g.}\ Figure 5) also identify the presence
of the CNR's in
NGC4303, 4314 and 4321. These rings are generally redder than all other regions
in the galaxy, having colours consistent with those of typical old stellar
populations in ellipticals and spiral bulges (Table 4). However, given the
spatial distribution of the near IR emission in the rings, a significant
contribution from M giants/supergiants in regions of active star formation
is likely. This is consistent with the situation in NGC1097 and 5728, where
an excellent spatial correspondence is seen between the 2.2$\mu$m luminosity
peaks and H$\alpha$ emission within each CNR (Paper I).
The inferred radii of the rings identified in the present data (Table 4) are in
good agreement with those measures tabulated by Pogge (1989) and Arsenault
{\sl et al.}\ (1988). The ring colours we derive for NGC4314 agree with those
measured by Benedict {\sl et al.}\ (1992).

Despite its inclusion in a list of galaxies possessing H$\alpha$ nuclear rings
(Pogge 1989), NGC3351 only shows evidence for a CNR within our (J--K)
colour map --
Figure 5 (a). The brightest emission regions within this ring have (H--K)
$\sim$0.37 mag. and (J--K) $\sim$0.99 mag., the latter value being equivalent
to that of the nucleus itself. Moreover, the major axis of this ring
is aligned with the twisted isophotes identified in Table 2. Both the ring
orientation and the location of the maxima in (J--K) coincide exactly with the
structure seen in CO and discussed by Kenney {\sl et al.}\ (1992).

We have investigated a possible link between CNR's and nuclear colours
by comparing the objects in Table 3 with the catalogue of nuclear
rings and pseudo--rings in Buta \&\ Crocker (1993). Although,
three of the galaxies
in our ``blue nuclei''
sample also have nuclear rings/pseudo--rings,
two of the galaxies in the control sample also show evidence of such
rings, implying little evidence of a correlation. Of course, such statistics
are only suggestive given the (distance dependent) selection effects inherent
in our data and those of Buta \&\ Crocker. Indeed, one must also consider
the possibility that the nuclear colours in the majority of galaxies in Table 3
are contaminated by
CNR's which remain unresolved in our data. This is particularly
important given that ILR's -- supposedly identifying the presence of CNR's
(Combes \&\ Gerin 1985) --
show a wide variation in size within the N--body simulations
conducted to date ({\sl e.g.}\ a factor of two in Paper I). There also exists the
possibility of temporal variations in the size of CNR's (Combes {\sl et al.}\ 1992).

Rings in
our data would remain unresolved if their radii were $\sim$ 3 pixels or less,
corresponding to $\sim$0.5 kpc since the distance to those blue nuclei
galaxies in Table 3 not identified as possessing CNR's is typically
28.0 ($\pm$1.5) Mpc. The possibility therefore remains that the majority
of colours in Table 3 are influenced by the presence of unresolved CNR's.
Clearly,
since the CNR's unambiguously identified in the present dataset are very red
(Table 4), their presence in the remaining data would imply that even bluer
underlying
colours exist within these nuclear regions.
Possible sources of such blue colours are
discussed in Section 4.2 below.

\section{Discussion}

\subsection{Nuclear bars and triaxial bulges}

A primary aim of the present study was to quantify the occurence of
isophote twists in the nuclear regions of barred galaxies. As Table 2
graphically illustrates, such twists are common in the near IR. Fully
$\sim$70 \%\ of the non edge--on galaxies in the present sample display
measurable twists in galaxies ranging from SB0 to SBbc.
It is therefore important to identify the mechanism(s)
by which such twists arise. Since contamination by spiral arms is minimal
in the regions concerned, nuclear bars or misaligned/triaxial bulges present
the most obvious cause.

The nature and properties of nuclear bars have been investigated in several
theoretical studies ({\sl e.g.}\ Shlosman {\sl et al.}\ 1989, Pfenniger \&\ Norman 1990).
Only two
investigations have used constraints imposed by observations of real galaxies.
In Paper I, Shaw {\sl et al.}\ argued that twists result from
the dynamical influence of the gas component (on the stars) as the gas follows
the stable x2 orbits immediately interior to the ILR. Conversely, Friedli \&\
Martinet (1993) view these systems as nested double barred galaxies. A
principal
difference between the simulations conducted in these studies is the r\^{o}le
of
dissipation. In Paper I viscosity is the mechanism by which the gas leads the
stars -- giving rise to the observed twists. Dissipation is
less important in the picture described by Friedli \&\ Martinet (1993). In
their scheme,
the twists are a direct manifestation of two misaligned bars possessing grossly
different pattern speeds.

We have identified what proportion of those galaxies displaying
isophote twists (other than NGC4321)
are candidates for possessing nuclear/misaligned bars. In doing so, three
methods are available to us. In the first instance,
we have compared the spatial dimensions of the observed twists
($\delta$R) to the
lengths (L) of the principal bars (Table 2). In the mean, $\delta$R/L
extends from 0.16 ($\pm$0.02) to 0.67 ($\pm$0.05), neglecting NGC4321 -- {\sl i.e.}\
considerably in excess of the 0.1--0.2 inferred by Friedli \&\ Martinet (1993)
or the values of 0.14 and 0.09 observed in NGC1097 and 5728 respectively
(Paper I). The only likely nuclear bar candidates on these grounds would be
NGC4643 and, possibly, NGC4274.
Moreover, in the hypothesis advanced in Paper I, the observed twists
must lie immediately interior to the radius of the CNR. Again neglecting
NGC4321, the mean dimensions of the observed twists have inner and outer
radii of 0.8 ($\pm$0.1)
kpc and 2.9 ($\pm$0.3) kpc respectively, seemingly larger than the typical
radii of CNR's (Paper I and references therein).

A second discriminant comes from the misalignment angles in Table 2. In
agreement
with the (limited) sample presented by Buta \&\ Crocker (1993), the isophote
twists lead or trail the spiral arms in roughly equal proportions, even for
the most face--on galaxies in the sample. [The measures in Table
2 have not been deprojected as the measured  inclination angles
are derived from projected axis ratios and are notoriously unreliable --
depending, for example, on spiral arm morphology.] This suggests that
relatively
few of the objects in Table 2 are candidates for double barred galaxies,
although a more reliable conclusion would require accurate deprojection.

Finally, we consider the orbital families in the barred potential. Stable
periodic orbits cannot cross (Sparke \&\ Sellwood 1987, Athanassoula
1992). Therefore, the minor axis dimension of the primary bar (identifying the
x1 orbital family) corresponds to the maximum dimension of any nuclear bar and
thus the radius of any CNR. We have determined the axis ratios of the
principal bars from the J--band images of each galaxy with measurable isophote
twists in Table 2. No deprojection was undertaken -- we assume the nuclear
and primary bars are coplanar and that projection foreshortens the dimensions
of each equally. Comparison of the minor axis dimensions to the radial scales
of the twists listed in this table suggests that nuclear bars are present in
NGC4262, 4274, 4314 and 4643, with NGC4371 and 4754 being additional, though
less
well--defined, candidates.

In conclusion, at most 7 of the 22 objects displaying near IR isophote twists
are likely to possess nuclear/misaligned bars of the type envisaged in
Paper I or Friedli \&\ Martinet (1993). Using the present data, a
discrimination
between these two studies would require a detailed investigation of the
r\^{o}le
of projection. As one must take into account the intrinsic figure of the
bar itself (Paper I), this analysis is beyond the scope of the present paper.
Observationally, a more definitive conclusion would result from detailed 2D
stellar and gas dynamics of galaxies possessing twisted isophotes.

It is possible that the observed twists in the remaining galaxies
result from the presence of triaxial bulges ({\sl e.g.}\ Kormendy 1979, 1982).
However, the
complex transition between bulge and primary bar luminosity components renders
unreliable any results derived from photometry alone. In fact, the observed
twists could readily result from a misalignment between a conventional (oblate)
bulge and the principal bar. The stability of such a configuration remains to
be determined, although the frequency of twists in the present sample implies
such a configuration must be stable. A definitive identification of triaxiality
in these galaxies could only be undertaken by combining full 2D luminosity
decomposition with extensive stellar and gas kinematics.

\subsection{The significance of blue nuclear colours}

Many galaxies in the present sample display blue nuclear colours. Indeed,
even in those galaxies displaying CNR's,
blue nuclei are evident on scales considerably smaller than the inferred
dimensions of rings (Figure 5).
We consider two mechanisms by which this blue colour
could arise.

\subsubsection{Extinction variations}

The first results from a {\it reduction}
in extinction within the
nuclear regions. The reddening vector in Figure 6 suggests a change of
A$_v$ $\sim$1 mag. would be sufficient to bring the IR colours of the majority
of the sample into agreement with those of typical old stellar populations.
Reduced extinction may result from the highly collimated nature of the gas
distribution in a bar. The gas is principally restricted to the leading edges
of the bar (as in, for example, NGC613: Hummel {\sl et al.}\ 1987a, and NGC1097:
Gerin {\sl et al.}\ 1988). The region interior to the corotation radius, but outside
the bar itself, is thus largely devoid of gas in many barred galaxies.
An aperture $\sim$0.4 kpc in
radius (Section 3.3) would encompass such regions, possibly leading to a
reduction in the near IR colour indices when compared to non--barred galaxies.
Optical colours would not support this assertion, since the nuclear colours
derived are invariably very red. However, large scale surveys of near IR
galaxian colour distributions remain to be undertaken.

\subsubsection{Nuclear activity}

The blue colour could be an observational consequence of nuclear
activity if it implies the existence of starbursts with suitable star formation
properties. This hypothesis is favoured by the fact that, of the 13 galaxies in
the sample whose nuclei have the spectral characteristics of starbursts/LINERS
(Arsenault 1989 and references therein), 10 have unusually blue nuclei.
Conversely, only 9 galaxies not showing such activity also possess blue
colours. Formally, sampling theory suggests that these ratios differ at the
93.6\%\ significance level using the ``two--tailed'' test. In
reality, however, the results are only suggestive given the limited
samples involved and the large apertures used in the spectral
classification (4--8 arcsec diameter -- Keel (1983) -- such that much of the
emission from the starburst candidates may come from CNR's rather than the
nuclei directly).

A link between blue colours and nuclear activity could be understood as
resulting
from the featureless blue continuum (FBC) detected in many Seyferts
and LINERS ({\sl e.g.}\ Yee 1983). This FBC is bluer (at optical wavelengths)
than the colours of spiral galaxy bulges, and may correspond to the
continuum in a young starburst (Terlevich \&\ Melnick 1985). A difficulty with
this
hypothesis results from the fact that
the blue colours in Table 3 are often distributed in a collimated structure
(Figure 5), and thus are more extended than the (unresolved) FBC in most
active galaxies. Moreover, the contribution of the FBC is likely to be small
in the near IR. The blue colours in our sample are observed even in galaxies
which are not classified as active whilst, in the scheme of Terlevich \&\
Melnick (1985), ``blue LINERS'' are galaxies which have already evolved beyond
a Seyfert 2 phase of activity.

There also exists the intriguing
possibility of a relationship between the blue nuclear
structures (Figure 5) and the anisotropic radiation field in active galaxies.
Specifically, these elongated blue structures may correspond to the blue,
collimated and extended emission--line regions (EELR's) commonly
observed in Seyferts ({\sl e.g.}\ Haniff {\sl et al.}\ 1988, Wilson {\sl et al.}\ 1988). The
suggestion of a strong causal relationship between nuclear activity and gas
motions within the innermost regions of a barred potential (Paper I and
references therein) would imply a correspondence between the orientations
of the EELR and any nuclear bar.

To investigate this possibility further, we have measured the alignments
of the (7) extended blue regions in our colour maps. They are compared
to the orientations of the nuclear isophotes in Table 5. There appears to be
a particularly close correspondence between these respective P.A.'s in NGC3412,
3945,
4536 and 4643. Even in the remaining cases, it is clear that there
exists a considerable misalignment between the extended blue nuclear regions
and the orientation of the primary bars. More extensive observational data
are required to further investigate this possible correspondence.

Near IR data for a small sample of active galaxies are presented by Hunt \&\
Giovanardi (1992). More detailed studies of the 2D near IR
luminosity distributions in LINERS/starbursts (Forbes {\sl et al.}\ 1992) suggest
the nuclei of such galaxies to be redder than typical old stellar populations
-- particularly in (H--K). The colours of the LINERS/starbursts in Forbes
{\sl et al.}\
are also consistent with those of Seyfert I nuclei (Kotilainen {\sl et al.}\ 1992b,
Kotilainen \&\
Ward 1994), although the emission mechanisms are likely
to be quite different.
These red colours have most recently
been interpreted as implying re--radiation
from hot dust at 2.2$\mu$m, influenced by
the proximity to the active nucleus (Kotilainen \&\ Ward 1994). The
contribution
from M giants/supergiants is, however, likely to be an important factor in
starbursts.

Unfortunately, there appear considerable uncertainties regarding the colours
derived for these LINER/starburst and Seyfert I samples.
Estimates of the sky background are unreliable in data of such limited
spatial coverage. Furthermore,
the removal of a nuclear (non--stellar) contribution
to each passband is highly uncertain given the assumptions of azimuthal
symmetry imposed, difficulties in the assignment of a suitable PSF and the
unreliable nature of complex luminosity profile deconvolution given the
limited data available (Kotilainen {\sl et al.}\ 1992a). The
apertures used by Kotilainen \&\ Ward (1994) are not sufficiently
small to eliminate significant contributions from circumnuclear star
formation in some objects. Interestingly, of the 4 Seyfert galaxies in
our sample (Table 1), only NGC4593
has an unusually red (H--K) colour. Our values of (J--H)=0.67
mag. and (H--K)=0.44 mag. for this object are in reasonable agreement with
those
of Kotilainen \&\ Ward (1994) over similar apertures.

Clearly, detailed studies of the near IR luminosity distributions of
LINER/starburst and Seyfert galaxies must await observations with high
spatial resolution
and large areal coverage.

\section{Conclusions}

In this paper, we present near IR photometry of 32 large, non--interacting
barred galaxies. Fully 21 of the 29 non edge--on objects in this sample
display measurable inner isophote twists with respect to the primary bars,
thereby confirming the frequency of this phenomenon inferred by Shaw {\sl et al.}\
(1993). Such twists are identified in galaxies as late as SBbc.

Comparisons have been undertaken between the scale of such isophote twists and
the typical dimensions of nuclear bars, circumnuclear rings (CNR's)
and the inferred
maximum dimensions of the x2 orbital family within each primary bar. We
conclude
that 7 of these galaxies may possess a
secondary (nuclear) bar misaligned from the primary.
NGC4274, 4321 and 4643 present
the clearest evidence of this morphology. There appears to be a correlation
between the axis ratios of the inner and primary bars, as the nuclear bars are
much
thinner and more distinct in later type galaxies.
However, we find no evidence of a
correlation between twisted near IR isophotes and the presence of CNR's.
Twists in the
remaining galaxies may reflect the presence of triaxial bulges or of
oblate bulges misaligned with the principal bar.

The near IR colour distributions in our sample yield evidence for very blue
structures within the
nuclei of 19 galaxies -- these regions being
bluer than typical old stellar populations by
$\sim$0.30 mag. in (J--H) and $\sim$0.23 mag. in (H--K). These blue colours
do not appear to correlate with the presence of CNR's.

It is unlikely that these blue colours
result from a contribution from young giants/supergiants within unresolved
circumnuclear star forming rings as the CNR's observed in the present data
are all red and probably dominated by emission from M stars.
If such blue regions reflect
reddening variations within the nuclei for objects with typical
elliptical/bulge
colours, a {\it reduction} in extinction of A$_{v}$ $\sim$1.0 mag. is implied.
This may result from the geometry of the gas distribution within the corotation
radius
of the barred potential.

Interestingly, 3 of the 4 galaxies
showing clear evidence for resolved CNR's in our data also possess blue
nuclei -- one of these being a Seyfert, one a LINER and one a starburst. This
may be suggestive of a link with nuclear activity, although
the extended
nature of the blue nuclear colours, even in galaxies not classified as active,
argues against their resulting from the featureless blue continuum observed
in many Seyferts/LINERs. However, a general equivalence between the
orientations of these blue regions and the innermost nuclear isophotes
may suggest that the former are a manifestation of the extended, blue,
emission--line regions commonly identified in some Seyfert galaxies.

\section*{Acknowledgements}

We thank Dr. F. Combes and D. Friedli for valuable input. MAS was supported
by PPARC throughout the course of this work. We acknowledge support, in
the form of data analysis software, provided by the Starlink Project.

\newpage
\parindent 0.0cm
\section*{References}

Arsenault R., 1989. A \&\ A, 217, 66

Arsenault R., Boulesteix J., Georgelin Y., Roy J.--R., 1988. A \&\ A,
200, 29

Athanassoula E., 1992. MNRAS, 259, 328

Athanassoula E., Martinet L., 1980. A \&\ A, 87, L10

Baumgart C., Peterson B., 1986. PASP, 98, 56

Benedict G., Higdon J., Tollestrup E., Hahn J., Harvey P., 1992. AJ,
103, 757

Boer B., Schultz H., 1993. A \&\ A, 277, 397

Bottinelli L., Gouguenheim L., Paturel G., \&\ de Vaucouleurs G.,
1984. A \&\ A, 56, 381

Buta R., Crocker D., 1993. AJ, 105, 1344

Combes F., Gerin M., 1985. A \&\ A, 150, 327

Combes F., Gerin M., Nakai N., Kawabe R., Shaw M., 1992. A \&\ A, 259, L27

Devereux N., Kenney J., Young J., 1992. AJ, 103, 784

Elias J., {\sl et al.}, 1982. AJ, 87, 1029

Ellis T., {\sl et al.}, 1992. Proc. SPIE, 1765, 94

Forbes D., Ward M., DePoy D., Boisson C., Smith M., 1992. MNRAS, 254,
509

Friedli D., Martinet L., 1993. A \&\ A, 277, 27

Garcia--Barreto J., Downes D., Combes F., Gerin M., Magri C., Carrasco L.,
Cruz--Gonzalez I., 1991. A \&\ A, 244, 257

Gerin M., Nakai N., Combes F., 1988. A \&\ A, 203, 44

Griersmith D., Hyland A., Jones T., 1982. AJ, 87, 1106

Haniff C., Wilson A., Ward M., 1988. ApJ, 334, 104

Hawarden T., Mountain C., Leggett S., Puxley P., 1986. MNRAS, 221, 41p

Heckman T., 1980. A \&\ A, 88, 365

Hummel E., J\"{o}rs\"{a}ter S., Lindblad P., Sandqvist A., 1987a.
A \&\ A, 172, 51

Hummel E., van der Hulst J., Keel W., 1987b. A \&\ A, 172, 32

Hunt L. K., Giovanardi C., 1992. AJ, 104, 1018

Jarvis B., Dubath P., Martinet L., Bacon R., 1988. A \&\ A supp.,
74, 513

Jedrzejewski R., 1987. MNRAS, 226, 747

Keel W., 1983. ApJ suppl., 52, 229

Kenney J., Wilson C., Scoville N., Devereux N., Young J., 1992. ApJ,
395, L79

Kormendy J., 1979. ApJ, 227, 714

Kormendy J., 1982. ApJ, 257, 75

Kotilainen J., Ward M., Boisson C., DePoy D., Bryant L., Smith
M., 1992a. MNRAS, 256, 125

Kotilainen J., Ward M., Boisson C., DePoy D., Smith
M., 1992b. MNRAS, 256, 149

Kotilainen J., Ward M., 1994. MNRAS, 266, 953

Palumbo G., Tanzella--Nitti G., Vettolani G., 1980. Catalogue of
Radial Velocities, Gordon \&\ Breach, London

Pfenniger D., Norman C., 1990. ApJ, 363, 391

Pogge R., 1989. ApJ supp., 71, 433

Schommmer R., Caldwell N., Wilson A., Baldwin J., Phillips M,
Williams T., Turtle A., 1988. ApJ, 324, 154

Shaw M., Combes F., Axon D., Wright G., 1993. A \&\ A, 273, 31 (Paper I)

Shlosman I., Frank J., Begelman M., 1989. Nature, 338, 45

Sparke L., Sellwood J., 1987. MNRAS, 225, 653

Terlevich R., Melnick J., 1985. MNRAS, 213, 841

de Vaucouleurs G., de Vaucouleurs A., Corwin H., 1976. Second
Reference Catalogue of Bright Galaxies, University of Texas press

de Vaucouleurs G., Buta R., 1980. AJ, 85, 637

Veron--Citty M., Veron P., 1989. Catalogue of Active Galaxies and
Quasars, ESO scientific report no. 13

Wilson A., Ward M., Haniff C., 1988. ApJ, 334, 121

Wilson A., Braatz J., Heckman T., Krolik J., Miley G., 1993. ApJ,
419, L61

Winge C., Pastoriza M., Storchi--Bergmann T., Lipari S., 1992. ApJ, 393,
98

Yee H., 1983. ApJ, 272, 473

\newpage
\parindent 1.0cm
\normalsize
\section*{Figure captions}

\begin{description}
\item [Fig. 1 :] Contour plots from subregions of the J--band images of
NGC3941, 3953, 4643 and 4754. In all cases, N is down, E to the left. The
axes are marked in pixels (1 pixel = 1.4 arcsec) and the data are
contoured from 14.0 to 18.4 mag. arcsec$^{-2}$ in 0.4 mag. arcsec$^{-2}$
intervals. The primary bar position angles, derived as discussed in the text,
are shown by the dashed lines.

\item [Fig. 2 :] Plots of the radial variation of position angle,
ellipticity, $\cos 4\theta$ and azimuthally averaged surface brightness
($\mu$)
derived from ellipse fitting to the galaxies displayed in Figure 1. The
profiles
for each passband are distinguished by the symbols shown in the legend
accompanying the lower right hand panel for each galaxy. The horizontal
dashed lines in the position angle profiles denote the P.A. of the primary
bars,
whilst the arrow identifies the inferred minor axis dimensions of these bars.

\item [Fig. 3 :] As Figure 1, but for NGC4274 (a) and NGC4321 (b). An expanded
plot of the central region of (b) is given in (c). Data in (a)
are contoured using the levels adopted in Figure 1. In (b)
contours are plotted from 14.0
to 17.6 (J) mag. arcsec$^{-2}$ in 0.4 mag. arcsec$^{-2}$ intervals, whilst in
(c)
contours are from 14.0 to 16.2 mag. arcsec$^{-2}$ in 0.2 mag. arcsec$^{-2}$
intervals.

\item [Fig. 4 :] As Figure 2, but for NGC4321.

\item [Fig. 5 :] Greyscale plots of the 2D colour distributions within:
(a) NGC3351 (J--K); (b) NGC3945 (H--K); (c) NGC4303 (J--H); (d) NGC4321 (J--H)
and (e) NGC4754 (J--H).
Display levels are: (a) 0.5 (black) to 1.0 (white) mag.; (b) 0.0 to 0.35 mag.;
(c) 0.35 to 0.7 mag.; (d) 0.25 to 1.2 mag. and (e) 0.3 to 0.65 mag. The
spatial coverage is: (a) 91x92 arcsec; (b) 117x109 arcsec; (c) 127x139 arcsec;
(d) 206x186 arcsec and (e) 104x135 arcsec.

\item [Fig. 6 :] Near IR colour-colour diagram of the present sample.
Also displayed is the region occupied by typical old stellar populations
(circle),
and the track delineated by a $\sim$30 \%\ contribution from a
600 K black body. The influence of a 10\%\ contribution from hot (T=200K)
dust at 2.2$\mu$m and a reddening vector of A$_v$ = 1 mag. are also shown. The
galaxies in the present sample are plotted according to their classifications
inferred from their emission--line spectra in Keel (1983). Typical
uncertainties
on both colours are plotted in the lower right corner. Galaxies identified as
Seyferts in Veron--Cetty \&\ Veron (1989) are marked by filled symbols.

\end{description}

\section*{Table captions}

\begin{description}
\item [Table 1 :] Properties of the present barred galaxy sample.

\item [Table 2 :] Summary of results derived from ellipse fitting to
the sample galaxies.

\item [Table 3 :] Results derived from analyses of the colour
distributions in the sample galaxies.

\item [Table 4 :] Colour indices for the nuclear rings in NGC3351, 4303,
4314 and 4321. Ring dimensions are quoted as major x minor axes, derived
from the (J--H) colour maps or, in the case of NGC3351, (J--K).

\item [Table 5 :] Comparison of the orientations of the extended blue regions
in the colour maps and the nuclear isophote orientations listed in Table 2.

\end{description}

\newpage
\section*{Notes to Table 1}

\begin{description}
\item [source of type :] RC2 -- de Vaucouleurs, de Vaucouleurs \&\ Corwin
(1976);
dVB -- de Vaucouleurs \&\ Buta (1980).

\item [dist. :] Object distance, taken from recessional velocities in
Palumbo, Tanzella--Nitti \&\ Vettolani (1980), using the distance moduli in
Bottinelli {\it et al.} (1984), or -- for NGC4314 -- the distance estimated
by Garcia--Barreto et al. (1991). An
H$_{\circ}$ of 50 km s$^{-1}$ Mpc$^{-1}$ is assumed throughout.

\item [classification :] Object classification from inclusion within the
starburst/A.G.N. or control samples of Arsenault (1989).

\item [comments :] nuc. ring -- nuclear ring; opt. twists -- twists in optical
isophotes with radius.

\item [source of comments :] BC93 -- Buta \&\ Crocker (1993); VV89 --
Veron--Citty \&\ Veron (1989); D92 -- Devereux, Kenney \&\ Young (1992);
K79 -- Kormendy (1979); J88 -- Jarvis {\it et al.} (1988); A88 -- Arsenault
{\it et al.} (1988).

\end{description}

\newpage
\section*{Notes to Table 2}

\begin{description}
\item [col. 2 :] A = anticlockwise, C = clockwise, measured from the POSS
or the present dataset assuming spiral arms are trailing.

\item [cols. 3--6 :] POSS corresponds to position angles defined
from the POSS survey plates, whilst ``ellipse fit'' lists those derived
from ellipse fitting to the present dataset (averages of all 3 near IR
frames with the exception of NGC3945
which is measured from the J--band image alone).

\item [col. 7 : ] Galaxy inclination angles derived from log R$_{25}$
measures in de Vaucouleurs, de Vaucouleurs \&\ Corwin (1976) assuming
an intrinsic disc axial ratio of 0.2.

\item [cols. 8, 9 :] Offset between inner near--IR isophotes and those
of the principal bar component. Measures derived from ellipse fitting
results except for NGC3945, NGC4340 and NGC4371 where $\theta$ is measured
directly from the images. For NGC3351, values refer to
the J and H images only.
Sign of $\theta$ identifies whether the inner
misaligned isophotes lead (+) or trail (--) the primary bar. In galaxies
for which the sense of rotation is undefined, the isophote twists are
denoted by ``$\pm$''.

\item [cols. 10, 11 :] Spatial scale over which the isophote twists are
observed.

\item [col. 12 :] Radius of the primary bar,
measured from the J--band images except for
NGC4321 which comes from the estimate of Arsenault et al. (1988).

\item [col. 13 :] The spatial scale of the misaligned isophotes ($\delta$R)
expressed as a fraction of the radius (L) of the primary bar.

\end{description}

\newpage
\section*{Notes to Table 3}

\begin{description}
\item [col. 2 :] Aperture radius (kpc) over which mean colour
indices in cols. 3,4 correspond. For the ``control'' sample, an
aperture of 2 pixels was adopted throughout.

\item [cols. 3, 4 :] Colour indices with errors corresponding to standard
deviations about these mean values within the specified apertures.

\item [cols. 5, 6 :] Magnitudes of colour differences between
nuclear structures in colour maps and the colours of the surrounding
regions.

\item [col. 7 :] Comments on structures evident in colour maps.

\end{description}

\newpage
\section*{Notes to Table 5}

\begin{description}
\item [col. 2 :] Colour map used in evaluation of entries in cols. 3 and 4.
All corresponds to J--H, H--K and J--K colour maps.

\item [col. 3 :] Position angle of the elongated blue structures evident
in the colour maps cited in col. 2. Where several colour maps have been
used, values quoted are weighted means.

\item [col. 4 :] Differences between PA$_{col}$ in col. 3
and the position angle of the primary bar (as defined in Table 2). Measured
in the sense PA$_{col}$ -- PA (primary bar).

\item [col. 5 :] $\theta$ measures in Table 2, expressed in the sense
PA (inner isophotes) -- PA (primary bar). The sign of $\theta$ may differ
from that in Table 2 as the latter measures identify whether the inner
isophotes trail or lead the primary bars.

\end{description}
\end{document}


Magnetic field intensification: comparison of 3D MHD simulations with Hinode/SP results

Abstract: Recent spectro-polarimetric observations have provided detailed measurements of magnetic field, velocity and intensity during events of magnetic field intensification in the solar photosphere. We consider the temporal evolution of the relevant physical quantities for three cases of magnetic field intensification in a numerical simulation. We determine the evolution of the intensity, magnetic flux density and zero-crossing velocity derived from the synthetic Stokes parameters by taking into account the spectral and spatial resolution of the spectropolarimeter (SP) on board Hinode. The three events considered show a similar evolution: advection of magnetic flux to a granular vertex, development of a strong downflow, evacuation of the magnetic feature, increase of the field strength and the appearance of the bright point. We find that synthetic and real observations are qualitatively consistent and, for one of the cases considered, agree very well also quantitatively. The effect of finite resolution (spatial smearing) is most pronounced in the case of small features, for which the synthetic Hinode/SP observations miss the bright point formation and also the high-velocity downflows during the formation of the smaller magnetic features.

\section{Introduction}

Magnetic field is ubiquitously present in the solar photosphere
\citep{deWijn:etal:2008}. On granular scales, it undergoes
continual deformation and displacement. It is swept by the
horizontal flows and concentrated in the intergranular lanes.
Flows are able to compress the field so that the magnetic energy
density $B^{2}/8\pi$ approaches the kinetic energy density $\rho
v^{2}/2$ of the flow \citep{Parker:1963,Weiss:1966}. This results
in a magnetic field strength of a few hundred Gauss at the solar
surface. Further intensification to kG strength is driven by the
mechanism referred to as: \textit{superadiabatic effect}
\citep{Parker:1978}, \textit{convective collapse}
\citep{Webb:Roberts:1978,Spruit:Zweibel:1979} or
\textit{convective intensification}
\citep{GrossmannDoerth:etal:1998}.

The first two concepts are a theoretical idealization of the
process. The superadiabatic effect contains the basic idea.
\citet{Parker:1978} pointed out that a thermally isolated
dowflowing gas within the flux tube in a superadiabatically
stratified environment will be accelerated, which would lead to
evacuation of the flux tube. Because of the resulting pressure
deficit, the gas inside the flux tube will then be pressed
together (together with the frozen-in magnetic field) by the
surrounding gas, causing the magnetic pressure to increase until a
balance of total pressure (magnetic + gas) is reached.

The convective collapse extents the concept to the convective
instability. It starts with a flux tube in thermal and mechanical
equilibrium with the surrounding hydrostatically superadiabaticly
stratified plasma. Since external stratification is convectively
unstable, any vertical motion within the flux tube can be
amplified. Downward flow will grow in amplitude and drain the
material from the flux tube. The process continues until a new
equilibrium with a strong field is reached. Different aspects of
the concept have been the subject of extensive research
\citep[see][for reviews]{Schuessler:1990,Steiner:1999}.
\begin{figure*}
    \centering
    \includegraphics[width=0.95\hsize,]{fig0.ps}
    \caption{Maps of the whole simulation domain at $t=140$~s.
    Normalized continuum intensity at $630$~nm(left), vertical component of magnetic field (middle) and velocity (right) at a geometrical height roughly corresponding
     to the level of $\langle\tau_{500}\rangle = 1$ ($\approx 930$~km above the bottom of the computational box).}
    \label{fig:snap}
\end{figure*}
The term convective intensification is used for magnetic field
intensification in realistic MHD simulations, where the process
occurs in its full complexity. It is driven by the thermal effect
in the surface layer of the magnetic concentration. There, due to
the presence of magnetic field, heat transport by convection is
reduced. The material inside the concentration radiates more that
it receives. This leads to cooling of material which thus starts
to sink and partial evacuation of the concentration occurs.
Contraction of the magnetic concentration by the surroundings
(result of the pressure imbalance) leads to an increase in
magnetic field strength. Thus, the simulations
\citep{Nordlund:1983,GrossmannDoerth:etal:1998,Gadun:etal:2001,Voegler:etal:2005,Cheung:etal:2008}
bear out the basic properties described by idealized concepts.
That is the downflow, the evacuation of the magnetic structure,
the field increase and, in some cases, establishment of a new
equilibrium. The 3D MHD simulations show that the strong magnetic
concentrations form as the horizontal flows in the intergranular
lanes advect weak, nearly vertical field and concentrate it at the
vertices of granular and mesogranular downflow lanes
\citep{Stein:Nordlund:1998,Stein:Nordlund:2006}. Larger magnetic
structures form at sites where a granule submerges and the
surrounding field is pushed into the resulting dark region.
Whether the formed concentration appears dark or bright in the
continuum intensity depends on whether the vertical cooling is
compensated or not by the lateral heating due to horizontal energy
exchange \citep{Bercik:etal:2003,Voegler:etal:2005}. This
formation scenario is consistent with the observations described
by \citet{Muller:1983} and \citet{Muller:Roudier:1992}. Their
observations show that network bright points form in intergranular
spaces, at the junction of converging granules as the magnetic
field gets compressed by the converging granular flow.

The 2D simulations by \citet{GrossmannDoerth:etal:1998} revealed
that magnetic flux concentrations formed by convective
intensification can evolve in different ways. They present two
possible outcomes. Depending on the initial magnetic flux, a
magnetic concentration can reach a stable state after the process,
or can be dispersed due to an upflow that develops as high speed
downflowing material rebounces from the dense bottom of the tube.
Similar results were presented by \citet{Takeuchi:1999} and
\citet{Sheminova:Gadun:2000}.

Observational evidence was found for both cases.
\citet{BelloGonzalez:etal:2008} reported on the formation of a
magnetic feature at the junction of intergranular lanes, without
any significant upflow observed. \citet{BellotRubio:etal:2001}, on
the other hand, detected a strongly blueshifted Stokes V profile
originating in a upward propagating shock, $13$ minutes after the
amplification of magnetic field.
\citet{SocasNavarro:MansoSainz:2005} found that supersonic upflows
are actually quite common.

Events that are interpreted as convective collapse were detected
also with the spectropolarimeter (SP) \citep{Lites:etal:2001} of
the Solar Optical telescope \citep{Tsuneta:etal:2008} on board
Hinode \citep{Kosugi:etal:2007}. Both, \citet{Shimizu:etal:2008}
and \citet{Nagata:etal:2008} show cases of high speed downflows
followed by magnetic field intensification and bright point
appearance. The event described by \citet{Nagata:etal:2008} shows
stronger field strength and upflow in the final phase of
evolution.

In this paper we give three examples of magnetic field
intensification from MURaM simulations and make detailed
comparison with the results of \citet{Nagata:etal:2008} and
\citet{Shimizu:etal:2008}.

\begin{SCfigure*}
    \centering
    \includegraphics[width=0.75\textwidth]{fig1_g25l.ps}
    \caption{Evolution of the continuum intensity at $630$~nm and the magnetic field in a $4\arcsec\times
    3\arcsec$ sized region. \textit{Left-hand side} (double column): original spatial resolution; red and blue contours outline the
    downflow of 6~km/s at 80~km and upflow of 4~km/s at 400~km above $\tau = 1$, respectively;
    horizontal lines mark the positions of vertical cuts shown in Figs.~\ref{fig:cut} and ~\ref{fig:cut_f23}; boxes corresponding to cases I, II and III at
    coordinates $[2\arcsec,1.5\arcsec]$, $[0.5\arcsec,1\arcsec]$ and $[3.5\arcsec,1\arcsec]$, respectively, enlarged in Fig.~\ref{fig:evol_zoom}.
    \textit{Right-hand side}: synthetic Hinode/SP observations; left column shows continuum intensity while
    the right column shows the apparent field strength (see the text); red contour marks the region
    with $0.01$~pm of signal excess (see text); yellow crosses mark the positions of pixels that we studied in detail in
    Fig.~\ref{fig:evol_red}. Grey scales cover the range of 0.6-1.6 and 0.8-1.2 for normalized intensity
    at original and reduced resolution, respectively, $\pm 1000$~G for the vertical component of magnetic field (original resolution) and
    $\pm200$~Mx/cm$^{2}$  for the apparent longitudinal magnetic flux density (Hinode resolution).}
    \label{fig:evol}
\end{SCfigure*}

\section{Simulation data and spectral synthesis}
\label{sec:data}

We use 3D radiative MHD simulations of a thin layer containing the
solar surface carried out with the MURaM code \citep{Voegler:2003,
Voegler:etal:2005} in a $6\times6\times1.68$~Mm domain with
non-grey radiative transfer included. Vertical and horizontal
spatial resolution is 10 km and 14 km respectively. The bottom and
top boundaries are open, permitting free in and outflow of matter.
The initial magnetic field of $\langle|B|\rangle = 200$~G, is
introduced in a checkerboard-like $2\times2$ pattern, with
opposite polarity in adjacent parts. As the field is redistributed
by convective motions, opposite polarities are pushed together and
dissipated, so that the unsigned magnetic flux decreases with
time. In this way, the run simulates the decay of the magnetic
field in a mixed polarity region. Local dynamo action
\citep{Voegler:2007} does not occur since the magnetic Reynolds
number is below the threshold for dynamo action. We do not expect
that the intensification process would be significantly different
at higher Reynolds numbers (smaller grid cells). The simulated
magnetic flux concentrations considered here are well resolved.
The main effect of an increased resolution would be a decrease in
the width of the boundary layers between the flux concentration
and the surrounding downflows.

We examine a 30 min sequence of simulation snapshots with a
cadence of approximately 90 s. The snapshots shown in this paper
have $\langle|B|\rangle~\approx~150$~G at $\tau =1$. In
Fig.~\ref{fig:snap}, we show the continuum image and maps of the
vertical components of the magnetic field and velocity for a
snapshot at time $t=140$~s ($t=0$~s corresponds to the first
snapshot considered in this paper).

In order to synthesize the Stokes profiles, the physical
parameters from the simulation are used as an input for the 1D LTE
radiative transfer code, SPINOR \citep{Frutiger:etal:2000}. A
spectral range that contains Fe I lines 630.15 and 630.25 nm is
sampled in steps of $7.5$ m\AA. The Fe abundance used for the
synthesis has been taken from \citet{Thevenin:1989} and the values
of the oscillator strengths from the VALD database
\citep{Piskunov:etal:1995}. Before comparing with Hinode/SP
\citep{Lites:etal:2001} observations, the synthetic line profiles
have been treated to bring them to the same resolution as Hinode
data. Firstly, a realistic point spread function PSF
\citep{Danilovic:etal:2008} has been applied to the synthesized
Stokes profiles. The PSF takes into account the basic optical
properties of the Hinode SOT/SP system and a slight defocus which
brings the continuum contrast of the simulation to the observed
value of 7.5\%. Secondly, to take into account spectral resolution
of the spectropolarimeter, the profiles are convolved with a
Gaussian function of $25$~m\AA~FWHM and resampled to a wavelength
spacing of $21.5$ m\AA. Thirdly, noise of $10^{-3} I_{c}$ is
added. Velocities are determined from the shift between the Stokes
$V$ zero-crossing wavelength and the line core position of the
Stokes $I$ profile averaged over pixels with polarization signal
amplitudes less than $10^{-3} I_{c}$. We correct these velocities
for convective blueshift by subtracting $150$~m~s$^{-1}$. Finally,
the procedure by \citet{Lites:etal:2008} is used to calculate the
longitudinal apparent magnetic flux density (shown in
Fig.~\ref{fig:evol}).

\begin{figure}
    \centering
    \includegraphics[width=0.95\hsize,angle=0]{hor_flow.ps}
    \caption{Enlargement of the regions, having a size of $0.75\arcsec\times 0.75\arcsec$, outlined by yellow squares in Fig.~\ref{fig:evol}. From left to
right: cases I, II and III. Continuum images (left column for each
case) and maps of the vertical components of the magnetic field
(right)  with overplotted horizontal velocities at approximately
80 km above $\tau = 1$ (white arrows). Red contours outline
locations of dowflows of 5 km/s or larger and yellow contours the
vertical component of vorticity exceeding $\pm0.2s^{-1}$ (cases I
and II) or $\pm0.1s^{-1}$ (case III).}
    \label{fig:evol_zoom}
\end{figure}

\section{Results}

Fig.~\ref{fig:evol} shows the evolution of the continuum intensity
and the magnetic field in a $2.9\times2.2$~Mm subdomain of the
simulation, over approximately $7$ minutes. Shown are maps at the
original resolution of the simulation (left-hand side) and at the
spatial resolution of Hinode (right-hand side). For the original
resolution, we show the vertical component of the magnetic field
from the simulations (right column on the left-hand side). The
longitudinal apparent flux density \citep{Lites:etal:2008}
retrieved from synthetic Stokes profiles is shown at the Hinode
resolution (right column on the right-hand side). During the
period shown, three bright points appear near the coordinates
$[2\arcsec,1.5\arcsec]$, $[0.5\arcsec,1\arcsec]$ and
$[3.5\arcsec,1\arcsec]$, in the regions outlined by yellow
squares. We refer to them here as cases I, II and III,
respectively. They are identified in the top left frame.

\subsection{Horizontal flows}
\label{sec:hor_flows}

In all three cases, magnetic field is advected by the flow to the
junction of multiple granules, where it is confined and
concentrated. The evolution of magnetic concentrations is examined
more closely in Fig.~\ref{fig:evol_zoom}, where the regions
outlined by yellow squares in Fig.~\ref{fig:evol}, are enlarged.
Horizontal velocities, at approximately $80$~km above
$\langle\tau\rangle=1$, are represented by arrows in continuum
maps and magnetograms. Regions with a strong vertical component of
vorticity are also outlined.

Vortex flows around strong downflows are quite common in
simulations of granulation \citep{Nordlund:1986}. They are formed
at the vertices between multiple granules, where flows converge
and angular velocities with respect to the center of downflow
increase due to angular momentum conservation. The lifetime of
vortex flows depends on the dynamical behavior of the neighboring
granules. In a $30$~min run that we examined, they last from less
then $90$~s to approximately $10$~min. The ones shown in
Fig.~\ref{fig:evol_zoom} are short-lived. In case II, the swirling
motion persists for at least 150 s, until the shape of the
neighboring granules is changed. The formed flux concentration is
then squeezed between two granules and stretched into a
flux-sheet-like feature. Later on, as the granules evolve, the
magnetic feature is advected to the left and trapped again in the
junction between newly formed granules. In case III, a strong
preexisting flux sheet is carried by the flow toward the junction
of the granules. There, the field is caught in a vortex flow that
lasts for at least $100$~s. The flow becomes disturbed and then
starts to swirl in the opposite direction. The vortex axes are
inclined, in both cases II and III, with respect to the vertical
direction. Also in both cases, there is a spatial and temporal
coincidence between the existence of vortex flows and the
formation of high speed downflows near the surface layer.

The vortices shown in Fig.~\ref{fig:evol_zoom} are less then
$0.25\arcsec$ in diameter, much smaller then the ones found in
observations \citep{Bonet:etal:2008}. The size of the observed
vortex flows is of order $1$\arcsec and the average lifetime is
around 5 min. They seem to outline supergranulation and
mesogranulation cells. As the authors of the observational study
suggest, most vortices could have been missed because of the
limited spatial resolution (Swedish Solar Telescope) and the
method used (tracking the motion of bright points in G band). A
short-lived vortex flow preceding formation of a network bright
point has been observed by \citet{Roudier:etal:1997}.

No vortex is detected during the formation of flux concentration
in the case I. In this case, a strong downflow, and later upflow
develops inside the flux concentration, whose shape and contrast
in the intensity map is greatly affected by the convective
motions. Spatial and temporal fluctuations across the magnetic
feature make it difficult to define its center. The bright points,
visible in continuum maps, form where fragments of flux
concentrations become narrow so that lateral heating becomes
important and in the regions with higher field strength, i.e.
higher evacuation.

During their further evolution after the intensification phase ($t
> 420$~s), the flux concentrations in cases I, II and III have
different fates. In case II, it meets a magnetic feature of
opposite polarity and vanishes in $\sim 3$~minutes. The feature
formed in case I is gradually stretched by the flow and fragmented
after approximately $9$~minutes. A few minutes after that, the
larger fragment encounters the opposite-polarity flux
concentration from case III. The cancellation continues for $\sim
12$~minutes until both features disappear. The smaller fragment of
the feature from case I is caught in a granular vortex flow
approximately $5$ minutes later which leads to further
intensification.

\begin{SCfigure*}
    \centering
    \includegraphics[width=0.75\textwidth]{hor_cut_rad_cc.ps}
    \caption{Time sequence of vertical cuts through the magnetic concentration corresponding to case I
    (positions of the cuts are marked by yellow horizontal lines in right-hand images in
    Fig.~\ref{fig:evol}). The plotted cuts correspond to the same instances of time as sampled in
    Fig.~\ref{fig:evol}. From left to right: density excess with respect to the the mean density at every
    height, vertical component of magnetic field, vertical component of velocity (overplotted arrows show the
    components of the velocity field in the x-z plane), and radiative heating rate.  Horizontal solid lines follow the levels of $\log~\tau = 0, -1$ and $-2$.
    Vertical dashed lines mark the position of the height profiles given
    in Fig.~\ref{fig:1dh}.}
    \label{fig:cut}
\end{SCfigure*}

\begin{figure}
   \centering
\includegraphics[width=0.85\hsize,angle=90]{1d_1p_cc.ps}
   \caption{Height profiles of various quantities in the flux concentration corresponding to the case I (magnetic field, vertical component of
   velocity, density and temperature) at the positions marked by dashed lines in
   Fig.~\ref{fig:cut} at different time instances. Thick vertical lines mark the positions of the optical depth $\tau_{500} = 1$ for each atmosphere.
    Thin vertical lines in the top two panels mark the positions of the height of formation of Fe~I~630.25~nm line (see text).}
   \label{fig:1dh}
\end{figure}

\subsection{Field intensification}

The evolution of the magnetic structures formed during the studied
time interval can be followed in more detail by taking vertical
slices through the regions marked by horizontal lines in
Fig.~\ref{fig:evol_zoom}. We will examine three cases.

\subsubsection{Case I}

Vertical cuts through the regions corresponding to this case are
given in Fig.~\ref{fig:cut}. Shown are the density excess
(relative to the mean density at each geometrical height), the
vertical components of magnetic field and velocity, and the
radiative heating rate. Cuts at the instant $t=0$ show the
fluctuation in density outlining the granulation structure, with
less dense, hotter granules and a cooler intergranular lane, where
magnetic field of a few hundred G is accumulated and a downflow of
a few km/s is present. The fluctuation in radiative heating rate
is highest in the thin layer near optical depth unity, where
upflowing granular material cools most strongly. The regions along
the intergranular lane seems to undergo cooling in the layers
above $\tau=1$. In the next instant, $140$~s later, a strong
magnetic concentration has already formed. Variations in magnetic
field strength at a given geometrical height are accompanied by
variations in the gas pressure. Regions with increased field
strength have lower density, which, in turn, shifts the level of
optical depth unity downwards relative to neighboring regions with
weaker field. The depression of the visible surface inside the
magnetic structure gives rise to radiative heating through the
sidewalls from the hot neighboring material. The downflow velocity
is significantly increased from $t=0$ to $t=140$~s, becoming
supersonic in the layer between 100 km above and 350 km below the
surface. By $t=280$~s it is greatly diminished again inside the
flux concentration. The diameter and the shape of the flux
concentration continuously evolves and does not reach a stationary
state. At $t=280$~s strong radiative heating is present inside the
concentration. This is also the moment when continuum intensity is
the highest as shown later (Fig.~\ref{fig:evol_red}).

In the last two sets of snapshots, shown in Fig.~\ref{fig:cut}, an
upflow is visible inside the flux concentration, extending almost
from the bottom of the simulation domain. It reaches supersonic
velocities in the upper layers of the photosphere. Its direction
is inclined with respect to the vertical. The upflowing material
inside the concentration exhibits strong cooling at the surface.
As the material refills the region inside the flux concentration,
the level of optical depth unity shifts upwards and the magnetic
field strength decreases. However, the upflow does not lead to a
complete dispersal of the field, as described in the last
paragraph of Sec.~\ref{sec:hor_flows}.

The evolution of this flux concentration is further illustrated by
the vertical profiles given in Fig.~\ref{fig:1dh}. Three phases
before and after the field intensification are shown. The
horizontal positions of the profiles are marked by vertical dashed
lines in Fig.~\ref{fig:cut}. The profiles of magnetic field
strength, vertical component of velocity, density and temperature
for each phase are shown, together with the corresponding level of
optical depth unity, marked by thick vertical lines. The level
$h=0$ corresponds to the level of the mean optical depth unity
$\langle\tau_{500}\rangle=1$ for the whole snapshot. Thin vertical
lines in the upper plots mark the position of the average height
of formation of the Fe I 630.25 nm line, which is defined by
calculating the centroid of the contribution function (CF) for
line intensity depression \citep{Solanki:Bruls:1994}, at the
wavelength of the line core. The definition gives only a rough
estimate of the heights sampled, since the lines are formed over a
large portion of the photosphere and the CFs are asymmetric in
wavelength due to the presence of high velocity gradients.

The figure shows that, in the interval from $t=0$~s to $t=140$~s,
the downflow extends to deeper layers and increases in amplitude,
reaching 10 km/s at 200 km below the surface. The magnetic field
strength increases from a few hundred to more than 2000 G at the
level of $\tau_{500}=1$. The significant reduction in density
results in a shift of 200 km in the optical depth unity level. As
a result of radiative heating the temperature gradient is flatter
and the temperature is around 400 K higher at $h=0$. The
temperature difference at equal optical depth is even larger,
exceeding 1000 K at $\tau_{500}=1$. The estimated height of
formation of Fe I 6302 shifts from roughly 300 km to 150 km above
the surface when magnetic field concentration forms, which is in
agreement with \citet{Khomenko:Collados:2007}. This effect
contributes to the increase of observed magnetic field strength in
the first two phases.

In the third phase, the upflow leads to a density enhancement and
a weakening of the magnetic field by a few hundred Gauss in the
lower photosphere. The level of $\tau_{500}=1$ shifts upward by
100 km. No discontinuity is visible in the vertical component of
velocity, in contrast to the event studied by
\citet{GrossmannDoerth:etal:1998}.

\begin{SCfigure*}
    \centering
    \includegraphics[width=0.75\textwidth]{hor_cut_f23_rad_cc.ps}
    \caption{The same as Fig.~\ref{fig:cut} but for cases II and
    III. Each frame is subdivided into a left part displaying case
    II and a right part showing case III.
    Vertical lines mark the positions of the height profiles given
    in Fig.~\ref{fig:1dh_f2} and Fig.~\ref{fig:1dh_f3}.}
    \label{fig:cut_f23}
\end{SCfigure*}

\begin{figure}
   \centering
  \includegraphics[width=0.85\hsize,angle=90]{1d_2p_cc.ps}
   \caption{The same as Fig.~\ref{fig:1dh}, but for case II.}
   \label{fig:1dh_f2}
\end{figure}

\begin{figure}
   \centering
\includegraphics[width=0.85\hsize,angle=90]{1d_3p_cc.ps}
   \caption{The same as Fig.~\ref{fig:1dh}, but for case III.}
   \label{fig:1dh_f3}
\end{figure}

\subsubsection{Cases II and III}

Vertical cuts through the regions corresponding to these cases are
shown in Fig.~\ref{fig:cut_f23}. The positions of the cuts are
marked by horizontal lines in Fig.~\ref{fig:evol_zoom}. The
parameters shown are the same as in Fig.~\ref{fig:cut}.

At $t=0$, weak magnetic structures in the intergranular lanes are
already present in both cases. Strong downflows along the magnetic
field lines persist over the entire time interval shown. They
reach supersonic velocities in some instances, but only in the
layers below $\tau =1$. In case II, the downflow is strongest at
$t=140$~s, just before the vortex flow becomes deformed owing to
evolution of the neighboring granule. Then it gets somewhat
reduced at t=280s, as the structure is squeezed between the
granules. Later, at $t=380$~s, a strong downflow is present again
in the newly formed granular junction, to which the flux
concentration is carried and where it is confined again. The
rightmost columns show strong cooling of the layers above the
level of $\tau_{500}=1$ in both cases. At $t=420$~s, the interior
of the flux concentrations reaches radiative equilibrium. At that
moment, the flux concentrations in both cases are almost vertical,
with diameters of approximately $0.1\arcsec$ at $\tau_{500}=1$.

Figures~\ref{fig:1dh_f2} and \ref{fig:1dh_f3} show height profiles
along the positions marked by the dashed lines in
Fig.~\ref{fig:cut_f23}. Since the flux concentrations are
inclined, the choice of curved lines, along the axes of magnetic
structures, would be more logical here. However, we are interested
in what would be observed at the disc center and the observables
(i.e. the Stokes profiles) depend on the vertical profiles of the
physical parameters which are shown here. In the case II, we
choose to show profiles at instants $t=0$, $280$ and $t=420$~s
(solid, dotted, dashed respectively). For the case III, profiles
at the last three instants from Fig.~\ref{fig:cut_f23} are shown.
Both cases exhibit a similar behavior: there is a significant
evacuation of the magnetic concentration due to the strong
downflow and an increase of the magnetic field strength from a few
hundred G to kG values in the process. There is a persistent
dowflow that reaches 5 km/s (in the vertical direction) at
$\tau_{500}=1$ and a decrease in density with a corresponding
shift in the height of $\tau_{500}=1$. In both cases, the shift in
$\tau_{500}=1$ reaches 150 km and the magnetic field strength
increases to 2000 G at that depth. The height of formation of the
Fe I 630.25 nm line shifts to deeper layers as the magnetic field
strength increases. In the last instant, the Fe I 630.25 nm line
probes layers with kG fields and vertical velocity of 5 km/s, in
both cases.

\begin{figure*}
    \centering
    \includegraphics[width=0.6\hsize,angle=90]{evolution_all_g25.ps}
    \caption{Temporal change of normalized continuum intensity (left), vertical component of magnetic field (middle) and vertical component of velocity
    (right), for cases I (top row), II (middle row) and III (bottom row) at original (stars/solid line) and Hinode resolution (crosses/dashed line).
    Vertical component of magnetic field and velocity at log~$\tau_{500} = -2$ (for original resolution) are each plotted together with
    the longitudinal apparent magnetic flux density (note the separate scale employed for B$^{L}_{app}$ in the middle column) and velocity retrieved from Stokes V profiles (Hinode resolution), respectively. }
    \label{fig:evol_red}
\end{figure*}

\begin{figure}
    \centering
    \includegraphics[width=0.4\hsize,angle=90]{stokes_i_or.ps}
    \includegraphics[width=0.4\hsize,angle=90]{stokes_v_or.ps}
    \includegraphics[width=0.4\hsize,angle=90]{stokes_i_g25.ps}
    \includegraphics[width=0.4\hsize,angle=90]{stokes_v_g25.ps}
    \caption{Synthetic Stokes I (left) and V (right) profiles at original (upper row) and Hinode resolution
    (lower row) for the times $t=0,140$ and $420$ s marked by dotted, solid and dashed,
    respectively. The intensity is normalized to the mean continuum intensity in
    the full snapshot domain.}
    \label{fig:stokes}
\end{figure}

\subsection{Comparison with synthetic Hinode observations}

The intensity contrast is considerably reduced at Hinode spatial
resolution, as can be seen in Fig.~\ref{fig:evol} (note that the
color scales are adjusted). The maps of longitudinal apparent
magnetic flux density are shown for comparison with the
magnetograms at the original resolution. Fine structure in both
maps, continuum and magnetograms, is lost due to smearing. The
bright points that correspond to cases II and III are barely
discernible.

\subsubsection{Case I}

In order to compare case I with the results of
\citet{Shimizu:etal:2008}, we calculate the parameter they refer
to as \textit{signal excess} (SE). It is defined as the Stokes $V$
profile integrated over the spectral range of $250-400$~m\AA~
redward from the nominal line center. This wavelength range
corresponds to dowflows of $7-14$ km/s, so that high values of SE
indicate the presence of hight-speed downflows. The red contour in
the right-hand side of Fig.~\ref{fig:evol} outlines a region with
SE of $0.01$~pm, which is a factor of 10 higher than in the
surroundings. It fits the location and time of the strong
downflows visible at the original resolution (left hand-side of
Fig.~\ref{fig:evol}). Both, synthesized (here) and observed
\citep{Shimizu:etal:2008} SE occur simultaneously with the
appearance of the bright point and the intensification of the
magnetic field. However, the event described by
\citet{Shimizu:etal:2008} (Fig. 8) lasts at least $6$ minutes,
while in our simulations it is present for less than $4$ minutes.

Figure~\ref{fig:evol_red} compares the observable parameters with
the corresponding values at original resolution. The locations
chosen for the plots are marked by yellow crosses in
Fig.~\ref{fig:evol}. The pixels are selected such that regions of
downflow and upflow are covered, as well as the evolution of the
bright point. The plots show the temporal change of normalized
intensity, longitudinal apparent flux density
\citep{Lites:etal:2008}, and zero-crossing velocity retrieved from
Stokes $V$ profiles of the Fe I 630.25 nm line, so that a direct
comparison with the results of \citet{Nagata:etal:2008} can be
made. Overplotted (stars/solid lines) are the normalized intensity
at the original resolution, the vertical component of magnetic
field and velocity at $\log\tau_{500}=-2$.

The comparison shows that the correlation between the parameters
at original and reduced resolution is not so good. This is a
consequence of different effects. Firstly, since the formation
region of the Fe I 630 nm lines extends over a large part of the
atmosphere, there is smearing along the line of sight and the
resulting profiles carry information from different layers with
different values of the physical parameters. Secondly, the
smearing due to spatial resolution affects every pixel differently
depending on its neighboring pixels. The profiles could be more
asymmetric and broader due to averaging over inhomogeneous regions
and the amplitude of Stokes $V$ can be affected.
Figure~\ref{fig:stokes} shows examples of Stokes $I$ and $V$
profile at original and reduced resolution at  $t=0,~140$ and
$420$~s. The profiles show stronger asymmetries and higher Doppler
shifts at original resolution. The amplitudes of Stokes $V$
profiles at Hinode resolution are reduced by a factor of 10. They
are of the same order of magnitude as the observed ones (Fig. 2 of
\citet{Nagata:etal:2008}). Both, synthesized and observed Stokes
$I$ profiles show an extended red wing at the phase with strong
downflow.

Figure~\ref{fig:evol_red} reveals that intensity contrast is
significantly reduced by spatial smearing. Thus at $t=230$~s the
contrast of the bright point drops from 100\% at original
resolution to below 30\% at the reduced resolution. The velocity
determined from Stokes V profiles is at most times somewhat lower
than the velocities at $\log\tau_{500}=-2$. However, the general
trends of the parameters are preserved. For both, simulation and
observation \citep{Nagata:etal:2008}, the intensity reaches its
maximum approximately 100 s later than velocity. This is
consistent with the idea that the enhanced brightness results from
the evacuation due to the downflow. When the velocity changes
sign, the brightness starts to fade. The maximum brightness of the
simulated and observed bright points are similar after degrading
the simulated intensity to Hinode resolution, while the velocities
derived from the simulations are consistent, but somewhat smaller.
They reach approximately $5$~km/s (downflow) and $0.5$~km/s
(upflow), while values up to $6$~km/s and $2$~km/s, respectively,
are observed. The magnetic field strength obtained by
\citet{Nagata:etal:2008}, by inversions, increases from a few
hundred to $2000$~G in approximately $200$~s, which is in
agreement with simulations.

\subsubsection{Cases II and III}

The middle and bottom rows in Figure~\ref{fig:evol_red} show the
temporal evolution for cases II and III, respectively. We find a
simultaneous increase in all parameters (continuum intensity,
apparent flux density and zero-crossing velocity). In both cases,
the downflow is suppressed at $t=280$~s possible because of the
evolution of the surrounding granules. The velocity at
$\log\tau_{500}=-2$ reaches 6 and 4 km/s for case II and III,
respectively. As the magnetic field strength increases, the
intensity follows and, in both cases, the intensity contrast
reaches 50\% with respect to the mean at original resolution.

The influence of spatial smearing is on the whole similar,
although a few differences relative to case I are noteworthy.
Firstly, due to the small size of the features, the correlation of
the magnetic and brightness signals is largely destroyed. Although
the magnetic signals are comparable to the values in case I, the
brightness is lower, so that the features are inconspicuous in
continuum images. The increase in magnetic flux density is
followed by just a small fluctuation in intensity. Secondly,
although the magnetic field strengths and diameters of the
features are comparable for cases II and III at original
resolution, the case II shows lower magnetic flux density after
convolution. Thirdly, the signature of the strong dowflow in case
II is lost due to spatial smearing.

\section{Summary}

Our case study of three examples of magnetic field intensification
has shown that in all three cases, the field is advected to the
junction of several granules. There, it is confined by converging
granular flows, which, in two cases, form a vortex. Owing to the
presence of the magnetic field, the thermal effect (radiative
cooling) induces evacuation of the flux concentration. The
evacuation leads to a downward shift of the optical depth scale
within the flux concentrations. The shift is smaller for the
smaller features due to the lateral radiative heating, which
inhibits further evacuation \citep{Venkatakrishnan:1986}. As a
result, the magnetic field at $\tau=1$, in the smaller features,
is weaker than in the case of the feature with more flux. This is
in accordance with recent numerical \citep{Cheung:etal:2007} and
observational \citep{Rueedi:etal:1992,Solanki:etal:1996} results.

During the evacuation, the dowflow velocities reach maximum values
of 5-10 km/s at $\tau=1$. In the case of the largest feature, the
downflow extends from the upper boundary of the simulation domain
and becomes supersonic in the lower photosphere. The magnetic
features formed have diameters of 0.1-0.2~\arcsec. In the case of
the biggest feature, a supersonic upflow develops approximately
$200$~s after the formation of the flux concentration. The upflow
does not lead to a complete dispersal of the field, but the
feature persists until 9 minutes later when it undergoes a
fragmentation. The disappearance of the features in all three
cases occurs when they meet opposite polarity features, between 3
and 20 minutes after their formation.

We also show what happens with the observables when the effects of
smearing to observational spatial resolution is taken into
account. An important result is that, in the case of small
features, Hinode/SP would miss the bright point formation and, in
some cases, also the high velocity downflows that develop in the
process. On the other hand, the signatures of the evolution of
large features are detectable even after the spatial smearing. We
show that this case can be quantitatively compared with Hinode/SP
observation \citep{Nagata:etal:2008,Shimizu:etal:2008} and
exhibits a very similar evolution. This suggests that the magnetic
field intensification process in the MURaM simulations is a
faithful description of the process taking place on the Sun.
Furthermore, our study indicates that the analysis and
interpretation of the observations in terms of the convective
intensification process is well-founded.

\begin{acknowledgements}
Hinode is a Japanese mission developed and launched by ISAS/JAXA,
with NAOJ as domestic partner and NASA and STFC (UK) as
international partners. It is operated by these agencies in
co-operation with ESA and NSC (Norway). We thank R. Cameron for
valuable suggestions. This work was partially supported by WCU
grant No. R:31-0016 funded by the Korean Ministry of Education,
Science and Technology.
\end{acknowledgements}

\bibliographystyle{aa}

Solar and Stellar Astrophysics (astro-ph.SR)