The unprecedented optical outburst of the quasar 3C 454.3. The WEBT campaign of 2004-2005

Abstract: The radio quasar 3C 454.3 underwent an exceptional optical outburst lasting more than 1 year and culminating in spring 2005. The maximum brightness detected was R = 12.0, which represents the most luminous quasar state thus far observed (M_B ~ -31.4). In order to follow the emission behaviour of the source in detail, a large multiwavelength campaign was organized by the Whole Earth Blazar Telescope (WEBT). Continuous optical, near-IR and radio monitoring was performed in several bands. ToO pointings by the Chandra and INTEGRAL satellites provided additional information at high energies in May 2005. The historical radio and optical light curves show different behaviours. Until about 2001.0 only moderate variability was present in the optical regime, while prominent and long-lasting radio outbursts were visible at the various radio frequencies, with higher-frequency variations preceding the lower-frequency ones. After that date, the optical activity increased and the radio flux is less variable. This suggests that the optical and radio emissions come from two separate and misaligned jet regions, with the inner optical one acquiring a smaller viewing angle during the 2004-2005 outburst. Moreover, the colour-index behaviour (generally redder-when-brighter) during the outburst suggests the presence of a luminous accretion disc. A huge mm outburst followed the optical one, peaking in June-July 2005. The high-frequency (37-43 GHz) radio flux started to increase in early 2005 and reached a maximum at the end of our observing period (end of September 2005). VLBA observations at 43 GHz during the summer confirm the

\section{Introduction}
The flat-spectrum and highly polarized radio quasar \object{3C 454.3}
at $z=0.859$ belongs to the blazar class
of active galactic nuclei. The dominant non-thermal emission of blazars is ascribed
to a plasma jet oriented at a small angle to the line of sight, so that the radiation
is relativistically beamed.
The broad-band spectral energy distribution (SED) of blazars shows two components:
the lower-energy one is attributed to synchrotron radiation by relativistic electrons, while the
higher-energy one is commonly thought to be due to inverse-Compton scattering
of soft photons by the same relativistic electrons.
   \begin{figure*}
   \centering
   \includegraphics{4817fig1.eps}
      \caption{Optical (top) and radio (bottom) light curves from 1966 to the end of September 2005;
      the yellow strip indicates the period of the 2004--2005 WEBT campaign, while grey
      vertical lines show the times of the satellite pointings (see Fig.\ \ref{sed}).
      Bottom panel: black points indicate the 8 GHz light curve (1489 data points);
      red, green, and blue curves represent the cubic spline interpolations through the
      30-day binned light curves at 5, 14.5, and 37 GHz, respectively;
      blue points are sparse 37 GHz data.
      Top panel: besides the $R$-band light curve (2709 data points, 1502 of which are from
	the WEBT campaign, blue points), the radio hardness ratios
      $H_{37/8}$ (blue line) and $H_{14.5/8}$ (green line) are plotted as $18-H$; the light blue
      shading highlights where $H_{37/8}$ is harder than average.}
         \label{radop}
   \end{figure*}

In addition to these two components,
quasars can also show a ``big blue bump" in the UV band, which is interpreted
as the signature of the thermal disc feeding the supermassive black hole.
This component is mostly visible when the synchrotron emission is faint \citep[see e.g.][]{hag94,pia99}.

One of the main features of blazars is their strong emission variability at all wavelengths,
from the radio band to $\gamma$-ray energies, on time scales ranging
from a few minutes to several years.

\section{The WEBT campaign}
The historical optical light curve of 3C 454.3 starts from $\sim 1900$.
In the radio bands the source has been monitored since the 1960s.
Figure \ref{radop} shows the optical (top) and radio (bottom) light curves over the last $\sim 40$ years.
The historical radio data are from the University of Michigan Radio Astronomy Observatory (UMRAO) and
the Mets\"ahovi Radio Observatory (\citealt{ter05}, references therein, and unpublished data),
while data from the Crimean (RT-22), MDSCC (PARTNeR), Medicina, Noto, and
SAO RAS (RATAN-600) Observatories as well as data from the VLA/VLBA Polarization Calibration
Database\footnote{{\tt http://www.vla.nrao.edu/astro/calib/polar/}.}
contribute to the most recent part of the light curves.
In the figure the radio data are represented as black points (8 GHz), red line (cubic spline through the
30-day binned 5 GHz data), green line (the same for 14.5 GHz), while blue points/spline indicate the 37
GHz fluxes (including also 36 GHz data).
The radio light curves show prominent, long-lasting outbursts, with
higher-frequency variations leading the lower-frequency ones, as typically found in many blazars.

The optical light curve shows two very different phases.
Until $\sim 2001.0$ only moderate variability is observed in the range $R \sim 15$--17
(an exception being the higher activity shown in $\sim 1940$--1955, when $R$ reached $\sim 14$;
see \citealt{ang68}).
From 2001 the amplitude of the variations starts to increase.

On May 9.36, 2005 (at the start of the observing season) the source was observed at $R=12.0$\footnote{
This outburst peak likely represents the most luminous quasar state ever observed: $M_{B} \sim -31.4$
($H_0=70$ km $\rm s^{-1}$ $\rm Mpc^{-1}$, $\Omega_M=0.3$, $\Omega_\Lambda=0.7$), comparable only
to the 3C 279 outburst of 1937 ($M_{B} \sim -31.2$).
For comparison, the bright but nearby quasar 3C 273
has ``only" $M_{B} \sim -26.4$, while the most luminous quasars of the SDSS and of the $z>4$ sample
of \citet{vig03} reach $M_{B} \sim -29.1$ and $M_{B} \sim -30.2$, respectively.},
thus triggering a  multifrequency campaign by the Whole Earth Blazar Telescope
(WEBT)\footnote{{\tt http://www.to.astro.it/blazars/webt/}; see e.g.\ \citet{vil04a,rai05}.}.
In order to have information on both the rising and decreasing phases of the outburst,
data were collected from June 2004 to the end of September 2005.
In total, 5584 $UBVRI$ observations from 18 telescopes\footnote{Kyoto,
Osaka Kyoiku, Lulin, Mt.\ Maidanak, Abastumani, Crimean, MonteBoo,
Catania, Armenzano, Michael Adrian, Torino, Sabadell, Teide (BRT),
Roque de los Muchachos (KVA and Liverpool), Foggy Bottom, Kitt Peak (MDM), and
San Pedro Martir Observatories.}
were performed in this period;
moreover, $JHK$ data were taken at Campo Imperatore, Calar Alto, and Roque de los Muchachos (NOT).
Radio data from 1 to 43 GHz were acquired at several telescopes, mentioned above.
In the top panel of Fig.\ \ref{radop} the $R$-band light curve from the WEBT campaign
is displayed as blue dots (see also Fig.\ \ref{colori}): it is composed of 1502 points. In particular,
1146 data points were taken in 143.59 days,
from May 9.36 to September 29.95, 2005, during the outburst decreasing phase,
with a mean time separation of 3.0 hours, and
only 4 gaps longer than 36 hours (2--4 days).
Previous data are shown as red points; they come from the literature (including the recent
paper by \citealt{fuh06}) as well as from some WEBT observatories.
Here we present some results of this campaign;
a more detailed study is deferred to a forthcoming paper.
   \begin{figure*}
   \centering
   \includegraphics{4817fig2.eps}
      \caption{Time evolution of the $B-R$ colour index (top) during the WEBT campaign
      compared with the $R$-band light curve (bottom); orange points in the
      light curve highlight the $R$ data used to derive the $B-R$ indices;
      the red curve is a cubic spline interpolation through the 43 GHz light curve,
      arbitrarily scaled to fit the figure; the light blue strips indicate the epochs
      of the INTEGRAL (left) and Chandra (right) pointings, while the light green and red ones
      refer to spectra shown in Fig.\ \ref{sed}; the red triangles mark the VLBA epochs
      of Fig.\ \ref{maps}.}
         \label{colori}
   \end{figure*}

By looking at Fig.\ \ref{radop} no evident correlation appears between the
optical and radio events. In particular, a low radio state characterizes
the 2004--2005 optical outburst as well as
the previous optical activity, while the past strong radio
outbursts do not show any optical counterpart. Even the hardness ratio
between the 37 GHz and 8 GHz splines ($H_{37/8}=F_{37}/F_8$, plotted in the top panel as $18-H_{37/8}$,
blue line) does not show correlation with the optical, as is found by \citet{vil04b} for BL Lac.
Only during the WEBT campaign can one see a relation:
a fast rise of $H_{37/8}$ to its maximum value (1.43),
delayed by several months with respect to the rise of the optical outburst.
This growth of $H_{37/8}$ is essentially due to the increase of the 37 GHz flux
visible in the bottom panel. The 43 GHz light curve,
not plotted there (but in Fig.\ \ref{colori}, see below), confirms this trend.
In the top panel $H_{14.5/8}$ is also plotted (green line) for comparison.
In general, its trend mimics fairly well that of $H_{37/8}$, with some delay and
reduced variability amplitude. In particular, just a slight increase can be seen in 2005.
The 22 GHz flux and the corresponding $H_{22/8}$ display an intermediate trend, as expected.
   \begin{figure}
   \centering
   \includegraphics[height=13cm]{4817fig3.eps}
      \caption{VLBA maps at 43 GHz; the total intensity is shown by contours from
0.15\% up to 76.8\% of the August 2005 peak value, increasing by a factor of 2;
the straight lines show the electric vector and their length is proportional to the local
polarized intensity.}
         \label{maps}
   \end{figure}

The general lack of correlation between the radio and optical behaviour suggests that
the optically emitting region is not transparent to the radio frequencies, so that these must come from
another, outer part of the jet, possibly misaligned with respect to the former. It seems that before 2001
the radio part was more aligned with the line of sight and the radio variability was thus enhanced
by Doppler effects. On the contrary,
in the last 5 years the optical radiation would dominate the scenario due to a
smaller viewing angle of the corresponding emitting region.
Thus the flux increase we are seeing at the high radio frequencies should be produced
close to the optical region. The outburst will probably
propagate further out, towards lower-frequency emitting regions.
However, we cannot foresee whether it will be
visible, since this depends on the jet curvature.

A confirmation of the above general picture comes from the source behaviour in the mm bands, where
a strong outburst was observed to peak in June--July 2005, reaching $\sim 45$ Jy at 1 mm and $\sim
26$ Jy at 3 mm
(data from the IRAM 30-m telescope), while previous observations in 1985--2004 did not show any comparable
activity, the flux density ranging between $\sim 2$ and $\sim 15$ Jy in both bands
\citep[data from][ and unpublished data from IRAM and SEST]{ter92,ter98,ter04,ste93,tor96,reu97}.

The time evolution of the $B-R$ colour index during the optical outburst is shown in Fig.\
\ref{colori} (top panel), compared with the $R$-band light curve (bottom panel).
Orange dots in the latter panel highlight the $R$-band points used
to derive the $B-R$ colour indices.
In general, a ``redder when brighter" trend can be recognized.
Since 3C 454.3 is a quasar, this can be interpreted
as due to the contribution given by the thermal emission from the accretion disc,
which mainly affects the bluer region of the optical spectrum when the jet emission is faint.
In other words, there would be two optical components: a variable one from the jet and the other
from the disc, not necessarily variable. When the jet is brighter, its redder spectrum dominates,
and vice versa.
However, there seems to be a kind of ``saturation" effect in the $B-R$ trend:
in the brightest part of the outburst ($R \la 14$) the trend is not visible any longer.
Most likely, this is due to the dominance of the non-thermal jet emission,
which would display its usual ``bluer when brighter" behaviour,
thus balancing the opposite trend.

In the bottom panel of Fig.\ \ref{colori} the 43 GHz flux-density spline is also shown
(arbitrarily scaled to fit the figure, the range is $\sim 4$--13 Jy) for comparison with
the optical light curve. It peaks in September 2005, i.e.\ 2--3 months
after the mm outburst, which seems to have a similar delay with respect to the optical event.

\section{VLBA observations}
Figure \ref{maps} shows three recent VLBA maps at 43 GHz compared with an older map
dated April 2001\footnote{Multiepoch and multifrequency VLBA monitoring of 3C 454.3
(with polarization) was performed between May and December 2005 (Savolainen et al., in preparation).}.
The plot is constructed with both total and polarized intensity, normalized
to the peak values achieved at the epoch of August 2005 ($S_{\rm tot}=9.62$ Jy/beam,
$S_{\rm pol}=0.162$ Jy/beam).
The total intensity is shown by contours from
0.15\% up to 76.8\% of the August 2005 peak value, increasing by a factor of 2.
The straight lines show the
electric vector and their length is proportional to the local polarized intensity.

There is a significant increase in both total and polarized intensity of the VLBI core
compared with the values in April 2001 (see also Table \ref{vlba}), thus confirming the expectation
that the 43 GHz rising flux in 2005 seen in Fig.\ \ref{colori} comes from the VLBI core.
The three VLBA epochs are indicated in Fig.\ \ref{colori} as red triangles: there is a strict agreement
between the trend of the 43 GHz light curve and the total intensity values reported in Table \ref{vlba}.
Since the VLBI core of 3C 454.3 at 43 GHz used to have very low polarization
compared to the western feature at $\sim 0.7$ mas from the core, a strong increase
of the core polarization along with the total flux suggests
a new component emerging from the core. However, as discussed in the previous section,
the emerging feature will be more or less visible depending on the jet bending.
\begin{table}
\caption{Results of VLBA observations of 3C 454.3 at 43 GHz.}
\label{vlba}
\centering
\begin{tabular}{l c c }
\hline\hline
Epoch    &  Total Int.\ (mJy)  &  Polarized Int.\ (mJy)\\
\hline
April 14, 2001  &     1710    &       15  \\
June 25, 2005   &     6560    &       74  \\
July 16, 2005   &     6290    &      103  \\
August 28, 2005 &    10200    &      200  \\
\hline
\end{tabular}
\end{table}
   \begin{figure*}
   \centering
   \includegraphics[width=13cm]{4817fig4.eps}
      \caption{Spectral energy distribution of 3C 454.3 showing contemporaneous radio, near-IR, optical,
      and X-ray (Chandra and INTEGRAL) data during May 15--20, 2005. Previous data are also
      plotted for comparison, together with two other spectra from the WEBT campaign (see text for
      further details).}
         \label{sed}
   \end{figure*}

\section{Chandra observation}
3C 454.3 was observed by Chandra on May 19.7--21.0, 2005 with the HRC-LETG
instrumental configuration. The Chandra pointing was performed as
part of a ToO program on blazars in outburst, and lasted 112 ks.
We used version 3.3 of the CIAO software to reduce and analyse the data.
We extracted source and background spectra from the standard
HRC-LETG {\em bow-tie} regions, and then grouped the source spectrum to contain
at least 20 counts per channel. We used the fitting package {\tt Sherpa}
to fit the background-subtracted source spectrum in the range $\sim 0.2$--8 keV.
Here we report the results of the continuum spectral analysis, while
the study of narrow absorption features due to possible intervening Warm-Hot
Intergalactic Medium filaments is deferred to a forthcoming paper (Nicastro et al.,
in preparation).

The 0.2--8 keV spectrum is well described ($\chi^2/{\rm d.o.f.} = 1817/2368$)
by a rather flat power law, with photon index $\Gamma = 1.477 \pm 0.017$, absorbed by a large
column density $N_{\rm H} = (13.40 \pm 0.05) \times 10^{21}$ $\rm cm^{-2}$,
exceeding by more than a factor of 2 the Galactic value.
We measured de-absorbed fluxes
of $F_{0.2-2} = (5.5 \pm 0.2) \times 10^{-11}$ erg $\rm cm^{-2}$ $\rm s^{-1}$ and
$F_{2-8} = (8.4 \pm 0.2) \times 10^{-11}$ erg $\rm cm^{-2}$ $\rm s^{-1}$.

\section{Spectral energy distribution and conclusions}
In Fig.\ \ref{sed} we show the broad-band SED of the source.
The small black crosses in the radio--optical range indicate archival data taken from NED.
In the X-ray band we plotted data from old
observations by ROSAT in November 1991 \citep{sam97} and in May 1992 \citep{pri96},
as well as data from a BeppoSAX observation in June 2000\footnote{{\tt http://www.asdc.asi.it/bepposax/}.},
and data from a pile-up affected Chandra observation in November 2002
\citep{mar05}.

During the brightest phases of the optical outburst, in May 2005,
there were other observations by high-energy satellites (RXTE, Swift, INTEGRAL)
besides the Chandra one presented in the previous section.
The ToO INTEGRAL observation was performed on May 15.8--18.4.
In Fig.\ \ref{sed} we show the corresponding 3--200 keV spectrum \citep{pia06}.
The Chandra spectrum of May 19.7--21.0 is also plotted, together with radio, near-IR and
optical data from the WEBT in the period of the two pointings (blue symbols).
Since in this period the optical brightness varied by $\Delta R = 0.9$,
we plot two $UBVRIJHK$ spectra (each composed of simultaneous data) corresponding to a high and a low state.
No significant variation
in the spectral slope is seen between the two states, despite their different brightness levels, in
agreement with the rather stable colour index found in this period (see Fig.\ \ref{colori}, light blue strips).
The contemporaneous radio data follow a power law from 1 to 22 GHz, which then breaks due to the
37 and 43 GHz fluxes that begin to rise (see Figs.\ \ref{radop} and \ref{colori}). The grey
vertical strips in Fig.\ \ref{sed} show the historical variation range of the radio, mm and
$R$-band data collected for this paper (and partly displayed in Figs.\ \ref{radop} and \ref{colori}).

One month after the May satellite pointings, around June 21.0, the reddest observed optical state is
achieved (light green line in Fig.\ \ref{colori}). The corresponding near-IR--optical spectrum is
displayed in Fig.\ \ref{sed} with green symbols: it is indeed very steep, crossing the lower May spectrum
around the $H$ band, and most likely connecting with the highest parts of the mm strips, since
we should be close to the mm outburst peak. The 37--43 GHz flux has also increased
(see Figs.\ \ref{radop} and \ref{colori}), while no appreciable variation is seen at lower frequencies
(where different-colour symbols often overlap in the figure, so that only the blue circles are visible).

In late September, at the foot of the outburst, the bluest optical state of 2005 is found
(light red strip in Fig.\ \ref{colori}). The corresponding (red) spectrum in Fig.\ \ref{sed} is the average
over 5 days (September 23.0--28.0) of intensive $UBVRIJHK$ monitoring: the thermal component seems to strongly
affect the spectrum. At this time the 22 GHz flux starts to increase and the radio spectral break shifts
downwards,
while a peak of the outburst is seen at 37--43 GHz (see also Figs.\ \ref{radop} and \ref{colori}).

The main feature of the plotted SEDs is the strong variability of the optical--IR, mm
and X-ray fluxes.
In Fig.\ \ref{radop} vertical grey lines indicate the times of the X-ray pointings.
The X-ray flux seems to correlate neither with radio flux (always rather low),
nor with the optical level
(in 2002 it was brighter than during the previous pointings, but the X-ray flux was lower),
nor with the radio hardness (higher during the ROSAT pointings and intermediate during the others).
It seems instead to correlate with the mm flux, since a strong mm outburst was close to its peak
during the May 2005 X-ray observations. This is expected
due to the probable inverse-Compton nature of the hard X-ray spectrum.

Thus, this unprecedented quasar optical outburst preceded an equally unprecedented mm outburst (accompanied
by an exceptionally bright X-ray state) by 2--3 months, and, 2--3 months later, it had fully propagated to the
high radio frequencies. This occurred in the VLBI radio core. We cannot foresee whether and when
further propagation will be visible at lower frequencies and outside the radio core,
due to probable jet bending. Further multifrequency monitoring
is needed to follow the next phases of this exceptional event.

\begin{acknowledgements}
The mm data from the IRAM 30-m telescope were provided by H.\ Ungerechts
based on preliminary results from regular flux monitoring observations
by the IRAM Granada staff.
This work is partly based on observations made with the Nordic Optical Telescope, operated
on the island of La Palma jointly by Denmark, Finland, Iceland,
Norway, and Sweden, in the Spanish Observatorio del Roque de los
Muchachos of the Instituto de Astrof\'{\i}sica de Canarias. It is partly based also
on observations collected at the German-Spanish Astronomical
Center (DSAZ), Calar Alto, operated by the Max-Planck-Institut
f\"ur Astronomie Heidelberg jointly with the Spanish National
Commission for Astronomy.
We thank Calar Alto for allocation director's discretionary time
to this programme.
This research has made use of data from the University of Michigan Radio Astronomy Observatory,
which is supported by the National Science Foundation and by funds from the University of Michigan,
and of the NASA/IPAC Extragalactic Database (NED),
which is operated by the Jet Propulsion Laboratory, California Institute of Technology,
under contract with the National Aeronautics and Space Administration.
This work was partly supported by the European Community's Human Potential Programme
under contract HPRN-CT-2002-00321 (ENIGMA). The St.\ Petersburg team
acknowledges support from Russian Federal Program for Basic Research under grant 05-02-17562.
RATAN-600 observations were partly supported by the
Russian Foundation for Basic Research grant 05-02-17377.
C.\ D\'{\i}az acknowledges support from the Spanish Programa Nacional de
Astronom\'{\i}a y Astrof\'{\i}sica under grant AYA2003-1676.
\end{acknowledgements}


Where the granular flows bend

Abstract: Based on IMaX/Sunrise data, we report on a previously undetected phenomenon in solar granulation. We show that in a very narrow region separating granules and intergranular lanes the spectral line width of the Fe I 5250.2 A line becomes extremely small. We offer an explanation of this observation with the help of magneto-convection simulations. These regions with extremely small line widths correspond to the places where the granular flows bend from mainly upflow in granules to downflow in intergranular lanes. We show that the resolution and image stability achieved by IMaX/Sunrise are important requisites to detect this interesting phenomenon.

\section{Introduction}

Solar granulation is one of the most studied phenomena in the
lower solar atmosphere and its properties are thought to be well
understood \citep{Bray1984}. Spectral and imaging observational
data have allowed solar granulation to be studied in detail.
Commonly considered parameters include sizes and lifetimes of
granules, velocity and intensity fluctuations together with their
correlations, proper motions of individual features, etc.
\citep[e.g.,][]{Title+others1989, Karpinsky1990, Nesis+etal1993,
Espagnet+Muller+Roudier+Mein+Mein1995,
Hirzberger+Vazquez+others1997, Hirzberger1999, Kostik+etal2009}.
At the same time, state of the art simulations of solar convection
and magneto-convection reproduce very satisfactorily many
observational granulation properties, which has made it possible
to explain self-consistently the mean shapes of spectral line
profiles \citep{Nordlund+etal2009}. Simulations have shown that
parameters as micro- and macroturbulent velocities (used widely in
the past to match the spectral lines shapes to observations) are
not necessary any more as the spectral lines are broadened by the
velocity field of the convection motions.

In this Letter, we report on spatially resolved observations of
the line broadening by the granular velocity field (or,
equivalently, the spatially resolved action of solar
microturbulence). The data used in the Letter were obtained with
the Imaging Magnetograph eXperiment \citep[IMaX;
see][]{MartinezPillet+etal2010} that flew on board the {\sc
Sunrise} balloon telescope in June, 2009 \citep{Barthol+etal2010,
Solanki+etal2010}. We study the spatial variations of the full
width at half maximum (FWHM) of the \FeI\ 5250.2 \AA\ line. We
found that the FWHM is larger in intergranular lanes and smaller
in granules. But the most interesting finding is that the FWHM is
particularly small in narrow lanes located at
granular/intergranular borders.
We found only few works that discuss the relation between the line
FWHM and granulation \citep[\eg ][]{Nesis+etal1992,
Hanslmeier+etal2008}, pointing out that the FWHM is
anti-correlated with brightness. The observations of
\citet{Nesis+etal1992} were interpreted in terms of supersonic
flows, which was questioned by
\citet{Solanki+Rueedi+Bianda+Steffen1996}. We did not find any
previous mention of the detection of the narrow regions with
extremely small FWHM adjacent to intergranular borders.
The aim of this Letter is to describe their properties and to
offer an explanation for their origin.

\begin{figure*}
\center
\includegraphics[width=16cm]{fig1a.eps}
\includegraphics[height=8cm]{fig1b.eps}
\caption{(a) Image of the full width at half maximum (FWHM) of the
\FeI\ 5250.2 \AA\ line in the IMaX field of view; the grey scaling
ranges from 90 to 150 m\AA; the mean FWHM is equal to 115 m\AA;
(b) same image after application of the unsharp masking filter.
Panels (c) and (d) give an enlarged view of the area of 6 Mm
$\times$ 6 Mm marked by a white square in panels (a) and (b). Blue
contours in panel (d) correspond to normalized continuum intensity
$I_c /\bar{I_c}=1$. Red dashed contours mark locations with
zero line of sight velocity. }\label{fig:fwhm} \vspace{-0.5cm}
\end{figure*}

\begin{figure*}
\center
\includegraphics[width=8cm]{fig2a.eps}
\includegraphics[width=8cm]{fig2b.eps}
\caption{{\it Left panel:} temporal variations of the FWHM in the
IMaX data along a virtual slit placed at one selected horizontal
location; the grey scaling ranges from 90 to 150 m\AA; the
mean FWHM is equal to 115 m\AA.  {\it Middle panel} gives the
same image after application of the unsharp masking filter. {\it
Right panel:} histograms of the LOS velocity in the IMaX data at
the locations where the continuum intensity is above unity, i.e.
in granules (red solid line); at the locations where the continuum
intensity is below unity, i.e. in intergranular lanes (red dashed
line); at the locations with enhanced FWHM marked white in
Fig.~\ref{fig:fwhm}b (black dashed line); and at locations with
extremely small FWHM marked black in Fig.~\ref{fig:fwhm}b (black
solid line).}\label{fig:time} \vspace{-0.5cm}
\end{figure*}

\section{Observed line width}

The IMaX instrument and the details of the data acquisition and
reduction are described in \citet{MartinezPillet+etal2010}. Here
we use a time series of images of about 23 min duration taken near
the solar disc center (UT 00:37). The observational sequence
consists of five polarized filtergrams taken at $\pm40$, $\pm80$
and $+227$ m\AA\ from the \FeI\ 5250.2 \AA\ line center. We use
Level 2 data (reconstructed by deconvolution). The spatial
resolution of the reconstructed data is estimated to be 0\farcs15.
The field of view is about 33\arcsec $\times$ 33\arcsec\ (after
removing the apodized areas), with a pixel size of 0\farcs055. We
use the following parameters extracted from the Stokes $I$ data:
continuum intensity, line of sight velocity, and FWHM. The latter
two parameters were calculated from Gaussian fits to the five
intensity wavelength points. The statistical errors of the
fitting procedure were within the range of 20--40 \ms\ for the
velocities and a few m\AA\ to 10 m\AA\ (at most) for the FWHM.

Figure \ref{fig:fwhm} shows an image of the FWHM at one moment of
time. This image demonstrates clearly that the width of the \FeI\
5250.2 \AA~intensity profile follows the granulation structure.
The FWHM increases in intergranular lanes (appearing as bright
regions) and becomes smaller in granules. A more careful
inspection of the image reveals that adjacent to almost each
intergranular lane there are narrow black regions (that also have
elongated lane-like shapes) of smaller widths.
The average FWHM of the \FeI\ 5250.2 line in the granular
regions is of 120 m\AA, in intergranular regions is of 140 m\AA,
but only 100 m\AA\ at the borders between granules and
intergranular lanes.
The unsharp mask filter\footnote{The algorithm that works by
enhancing the contrast between neighboring pixels in an image}
applied to this image enhances the gradients and makes the
locations with extremely small FWHM more evident
(Fig.~\ref{fig:fwhm}b).
Panels c and d with enlarged view of the region of 6 Mm $\times$ 6
Mm size (marked by the white squares in Figs.~\ref{fig:fwhm}a and
b) give some very clear examples of the dark lanes surrounding
bright intergranular spaces. These dark features sometimes appear
only on one side of the bright intergranular lane, and sometimes
on both sides. In many cases they appear in regions where the
normalized continuum intensity $I/\bar{I_c} \approx 1$, i.e. just
at the borders between granules and intergranular lanes. The width
of these dark features is about $3-4$ IMaX pixels, i.e.
$\sim$0\farcs2. Their lifetime is the same as the lifetime of the
intergranular lane they are attached to, see Fig.~\ref{fig:time}
(left and middle panels).
These dark lanes in the FWHM image appear also in the
non-reconstructed data (with less contrast), proving that they are
not artifacts of the phase diversity reconstruction procedure.

\begin{figure*}
\center
\includegraphics[width=17cm]{fig3.eps}
\caption{(a) Image of the FWHM of the \FeI\ 5250.2 \AA\ line
obtained from simulations degraded to IMaX resolution, the grey
scaling ranges from 90 to 150 m\AA; the mean FWHM is equal to 112
m\AA; (b) same image after application of the unsharp masking
filter; blue contours correspond to $I_c /\bar{I_c}=1$; red
dashed contours correspond to zero line of sight velocity; green
lines indicate locations considered in Fig. 5; (c) image of the
FWHM of the \FeI\ 6301.5 \AA\ line obtained from simulations
degraded to Hinode resolution; the grey scaling ranges from
123 to 160 m\AA; the mean FWHM is equal to 132 m\AA.
}\label{fig:fwhm_sims} \vspace{-0.5cm}
\end{figure*}

\begin{figure}
\center
\includegraphics[width=11cm]{fig4.eps}
\caption{Left panel: image of the FWHM obtained from simulations
when the spectral synthesis is done after setting horizontal
variations of thermodynamics variables to zero; the grey scaling
ranges from 90 to 150 m\AA; the average value is 112 m\AA. Right
panel: the same for the case that the synthesis is performed after
setting velocities to zero in the simulation snapshot; the grey
scaling ranges from 98 to 105 m\AA, the average value is 100 m\AA.
}\label{fig:std} \vspace{-0.9cm}
\end{figure}

\section{Origin of the dark lanes}

To find the origin of these dark lanes in FWHM images, we have
performed spectral line synthesis of the \FeI\ 5250.2 \AA\ line in
a snapshot of magneto-convection simulations \citep{Vogler2005}.
We considered a simulation run with a bi-polar magnetic field
structure and an average unsigned magnetic field in the box of 30
G, corresponding to a quiet solar region \citep[for more details
on this simulation snapshot see][]{Khomenko+etal2005}. The grid
size of the simulations is about 20 km. The synthesis was done in
LTE using the SIR code \citep{RuizCobo+delToroIniesta1992}. No
additional line broadening mechanisms, such as macro- or
microturbulence were introduced.

We found that the areas with extremely narrow intensity profiles
adjacent to intergranular lanes also appear in the simulations.
For an adequate comparison with the observations we have degraded
the simulation spectra in the following way. Firstly, we performed
a convolution of each monochromatic image with an Airy function
imitating the action of a 100 cm telescope. Secondly, we convolved
the spectra in wavelength with the transmission profiles of the
IMaX etalon. Finally, we fit a Gaussian function to the five
wavelength points and calculated the width of the Gaussian
profiles as a proxy to the intensity profile width, as well
as velocity, similar to the observations. The resulting image of
the simulated FWHM is shown in Figure \ref{fig:fwhm_sims}a.

The FWHM image shows structures similar to the observed ones, i.e.
larger FWHM in intergranular lanes and very low FWHM in the narrow
regions adjacent to lanes. Note that the average value and the
range of variations agree very well in the simulations and
observations.
Similar to the observations, the dark lanes in the FWHM image
correspond to the normalized continuum intensity around unity,
i.e. they mark the transition between granules and intergranular
lanes (Fig.~\ref{fig:fwhm_sims}b). Their width is a few pixels
($\sim$0\farcs2).

\begin{figure*}
\center
\includegraphics[width=17cm]{fig5.eps}
\caption{Examples of the velocity field in the simulations at
locations marked (A), (B) and (C) in Fig.~\ref{fig:fwhm_sims}b.
The color coding indicates vertical velocity, with red being a
downflow and blue being an upflow, both saturated at 3 \kms.
The white and black dotted lines indicate the effective formation
height of the continuum and line core of the \FeI\ 5250.2 \AA,
respectively, estimated from the calculation of the Response
Functions. The white solid line shows horizontal variations of
FWHM in arbitrary units. }\label{fig:flows}
\end{figure*}

Due to the small width of these dark lanes, about 0\farcs2, high
resolution and image stability are required for their detection.
Figure \ref{fig:fwhm_sims}c shows the FWHM after the simulated
spectra were degraded to Hinode resolution and pixel size (note
that the \FeI\ 6301.5 \AA\ is used in this case). This image
preserves brightenings in intergranular lanes corresponding to
large FWHM, but the adjacent dark lanes do not show up at this
resolution. On the one hand, this comparison demonstrates that
IMaX/{\sc Sunrise} data quality is necessary to detect these
features. On the other hand, it proves that the instrument reached
indeed very high spatial resolution \citep{Berkefeld+etal2010}.

The IMaX spectral line of  \FeI\ 5250.2 at \AA\ is sensitive to
the velocity broadening, but also to magnetic fields via the
Zeeman effect, and to temperature. To find out which effect is
primarily responsible for the spatial variations of the FWHM, we
performed additional calculations. We repeated the spectral
synthesis for three cases: (1) magnetic field was set to zero,
preserving the variations in the rest of the variables; (2)
horizontal variations of the temperature and gas pressure were set
to zero using for all the points the mean stratification, but
preserving variations of the velocity and magnetic field;
(3) velocity was set to zero at all points, while preserving
variations in the rest of the variables.

We found that the magnetic field's influence on the horizontal
variations of the FWHM is negligible except for a few isolated
locations in intergranular lanes where the field reaches kG
values. There the FWHM is additionally enhanced by the Zeeman
broadening. Figure \ref{fig:std} gives the variations of the FWHM
in the cases (2) and (3). Case (2), where only velocity variations
are preserved, looks almost identical to the original case in
Fig.~\ref{fig:fwhm_sims}a. While case (3), where only variations
of the thermodynamic variables are preserved, looks very
different, showing almost an opposite behaviour to the case (2).
Note that, curiously, in the case (3) the broadening due to
magnetic field becomes evident in some intergranular lanes where
the field is strongest (locations at $X =1-2$ Mm; $Y=1-2$ Mm and
$X =4$ Mm; $Y=3$ Mm in the right panel of Fig.\ref{fig:std}). The
broadening due to magnetic field is much weaker compared to the
velocity broadening, but is comparable in magnitude to the thermal
broadening.
The average value of the FWHM and its standard deviation is much
lower for the case (3) as compared to the case (2) and to the
observations. Thus, we conclude that velocity broadening is the
prime responsible for the horizontal variations of the FWHM.

Figure~\ref{fig:flows} shows several examples of the velocity
field at the locations marked by letters in
Fig.~\ref{fig:fwhm_sims}b. At these locations the red line crosses
an intergranular lane with enhanced FWHM as well as the transition
between granule and intergranule with very low FWHM. As follows
from Fig.~\ref{fig:flows}, all these locations have features in
common. At the center of an intergranular lane, in the deep
photospheric layers, the vertical velocity is mainly a downflow
with a magnitude reaching about $3-4$ \kms. The intergranular
lanes are also characterized by a strong vertical velocity
gradient. In a few hundred kilometers, the velocity decreases from
$3-4$ \kms\ to zero, following the change of its sign at the top
of the atmosphere. These strong velocity gradients have immediate
implications on the width of the \FeI\ 5250.2 \AA\ line profiles,
whose core forms around 400 km \citep{Gurtovenko+Kostik1989}. The
action of these vertical velocity gradients over the height range
of formation of the line is similar to the action of
microturbulence velocity, producing additional broadening of the
profiles. We conclude that the large FWHM in intergranular lanes
is produced by a microturbulence-like effect due to strong
vertical velocity gradients, unlike the mechanism based on
super-sonic flows proposed by \citet{Nesis+etal1992}. Similar
conclusion, as for the enhancement of the FWHM in intergranular
lanes, was previously reached in the work of
\citet{Gadun+etal1997} from studying the line-parameter variations
in the two-dimensional hydrodynamic simulations.

Figure~\ref{fig:flows} shows that the regions with extremely low
FWHM correspond to the locations where the granular flow bends.
These locations are the transitions from granules with upflowing
velocities to intergranular lanes with downflows. As the velocity
vector bends there, the vertical velocity component is close to
zero and the vertical velocity gradient is very small as well. It
makes the velocity broadening at these transition locations
extremely low, thus producing narrow profiles, as detected in the
IMaX observations. In Fig.~\ref{fig:fwhm_sims}b the red
dashed contours mark the location with zero vertical velocity in
simulations. It shows that almost always these locations coincide
with dark lanes in FWHM image.

The IMaX data only have five points in wavelength, making it
impossible to check the conclusions about the velocity gradients,
obtained from simulations. However, the absolute values of the LOS
velocities are reliably obtained from these data. According to
simulations, the regions with extremely low FWHM correspond to
locations with almost zero vertical velocities. In
Fig.~\ref{fig:time} (right panel), we plot the LOS velocity
histograms obtained from the IMaX data over the whole time series.
This figure shows the existence of the three distinct families of
features. Granules, selected as locations with $I_c /\bar{I_c}>1$,
have the velocity histogram peak at $-0.6$ \kms. Intergranular
lanes, selected as locations with $I_c /\bar{I_c}<1$, show
downflowing velocities with the histogram peak at 0.4 \kms. This
histogram almost coincides with the one obtained by selecting
locations with enhanced FWHM given by white regions in
Fig.~\ref{fig:fwhm}b. Finally, there is a third distribution
corresponding to the borders between granules and intergranular
lanes. It is obtained by selecting locations with low FWHM given
by the black regions in Fig.~\ref{fig:fwhm}b. Its maximum appears
at zero LOS velocity. This observational confirmation supports the
conclusion based on simulations that the regions with extremely
small FWHM correspond to locations where the velocity vector
bends, producing a zero LOS velocity component.

\section{Conclusions}

In this Letter we have investigated the spatial distribution of
the full width at half maximum of the \FeI\ 5250.2 \AA\ intensity
profiles in the IMaX/{\sc Sunrise} data. We have shown that the
FWHM becomes large in intergranular lanes and is very small in
narrow 0\farcs2 wide strips adjacent to the lanes. The lifetime of
these narrow regions is the same as the intergranular lanes they
are attached to and they evolve in parallel. Based on the
similarity between IMaX observations and simulations of
magneto-convection, we have offered an explanation for the
formation of these profiles with small FWHM. We suggest that these
profiles appear in regions where the granular flow bends. There,
the velocity broadening is extremely small, as the LOS velocity
component is close to zero and has very small gradients. This
conclusion is supported by the analysis of the LOS velocity
histograms in the IMaX data. We claim that the {\sc Sunrise}/IMaX
resolution and image stability is the necessary condition to
detect these narrow transition zones between granules and
intergranular lanes.

\acknowledgments The German contribution to {\sc Sunrise} is
funded by the Bundesministerium f\"{u}r Wirtschaft und Technologie
through Deutsches Zentrum f\"{u}r Luft- und Raumfahrt e.V. (DLR),
Grant No. 50~OU~0401, and by the Innovationsfond of the President
of the Max Planck Society (MPG). The Spanish contribution has been
funded by the Spanish MICINN under projects ESP2006-13030-C06 and
AYA2009-14105-C06 (including European FEDER funds). The HAO
contribution was partly funded through NASA grant number
NNX08AH38G. This work has been partially supported by WCU grant
No. R31-10016 funded by the Korean Ministry of Education, Science,
and Technology.

Solar and Stellar Astrophysics (astro-ph.SR)

Where the granular flows bend


Cognitive Reflection Test and the Polarizing Force-Identification Questions in the FCI

Abstract: The set of polarizing force-identification (PFI) questions in the FCI consists of six items all basically asking only one question: the set of forces acting on a given body. Although it may sound trivial, these questions are among the most challenging in the FCI. In this work involving 163 students, we investigate the correlation between student performance on the set of PFI questions and the Cognitive Reflection Test. We find that for both scores in the FCI as a whole and in the PFI questions, the range of values of the Pearson coefficient at 95\% confidence interval, is suggestive that cognitive reflection may be one of the contributing factors in the student performance in the FCI. This is consistent with the idea that high level of cognitive reflection may help in eliminating seemingly valid choices (misconceptions) in the FCI that are intuitive from everyday experience or "common sense" but otherwise misleading. The ability to activate System 2 in Dual Process Theory, whether from System 1 or right after reading a physics problem, may contribute in narrowing down the set of prospective valid answers in a given physics problem. Complementary to cognitive reflection are other factors associated with deep understanding of physics whose effects are expected to become more evident with the level of difficulty of a set of physics problems. Given two students with the same level of cognitive reflection, the one with deeper understanding of physics is more likely to get the correct answer. In our analysis, the range of correlation coefficient for the set of PFI questions is downshifted with respect to that for the FCI as a whole. This may be attributed to the more challenging nature of the latter compared to a significant fraction of the remaining questions in the former.

\section{Introduction}
In spite of the enormity and complexity of possible information configurations that it can process, it appears that the way human mind thinks or decides on a given problem or situation, can often be simplified into merely two systems. Dual process theory (DPT) in psychology distinguishes these two systems as \textit{intuitive} (`System 1'), one that is considered fast and autonomous, and \textit{analytic} (`System 2'), the other one that is considered as deliberative and slow; see e.g., Refs. \cite{Neys, Stanovich, Kahneman}. System 1 allows us to make quick decisions with minimal strain on our mental resources (e.g., recognizing a friend in a classroom). On the other hand, System 2 enables us to solve problems where higher-order thinking is required (e.g., solving a $10 \times 10$ \textit{sudoku}).

In dealing with different situations, it is helpful to efficiently choose the appropriate process to use. More properly, with System 1 often proceeding ``unconsciously,'' it is important to be cognizant of the need to activate System 2 on top of System 1. The Cognitive Reflection Test (CRT) \cite{Frederick} is a widely used measure for the propensity to override System 1 by System 2. It is composed of three-item set of questions designed to initially draw the problem solver into activating System 1. Considering the bat-and-ball problem in the CRT, for instance, using System 1, the tendency is to answer 10 cents. A quick reflection however, indicates that this cannot be correct; for then, the total cost would be $\$0.10 + \$1.10 = \$1.20$! The correct answer using short algebraic manipulation or trial-and-error is 5 cents.

Scores in the CRT are found to have moderate correlation with performance on Wonderlic Personnel Test \cite{Frederick}, heuristics-and-biases tasks \cite{Toplak}, incidences of conjunction fallacy \cite{Oechssler,Liberali}, and susceptibility to some behavioral biases \cite{Hoppe}, amongst others. From a simplified perspective, these studies suggest that people who tend to think deeply, as may be measured by CRT, are able to make better decisions and solutions to problems. In the field of Physics Education, to which this paper belongs, such a perspective although simplified, offers an attractive avenue to look into the possible correlation between CRT and student performance on some standardized tests. The idea is that DPT, as it relates to CRT, may be able to explain the way students approach and solve problems in Physics. This in turn, may give us an insight about the Physics learning process and as educators, the ways by which we may be able to improve it.

Following similar line of thinking, Wood, Galloway, and Hardy \cite{Wood} (see also Refs. \cite{Kryjevskaia,Gette}) investigated the question ``Can dual processing theory explain physics students' performance on the Force Concept Inventory? [FCI]'' They examined the relationship between student performance on CRT and FCI \cite{FCI,Savi} and found a moderately positive linear correlation between the two for both pretest and post test. The ``findings indicate that students who are more likely to override the system 1 intuitive response and to engage in the more demanding cognitive reflection needed to answer the CRT question correctly are also more likely to score highly on the FCI, implying that similar cognitive processes account for at least some of the cognitive abilities needed for each test.'' \cite{Wood}

This study intends to further probe the possible relationship between CRT, DPT, and the FCI. In particular, while not necessarily neglecting FCI as a whole, we wish to investigate student performance in the set of six Polarizing Force-Identification (PFI) questions (see Refs. \cite{Alinea, Alinea2}) in the FCI, namely, questions  5, 11, 13, 18, 29 and 30, with their scores in the CRT. This subset of the FCI effectively asks only one basic question: identification of the force(s) acting on a given body. Surprisingly, out of the five choices for each question, the majority of the students tend to only select two choices (polarizing choices). One of these choices is the correct answer containing the right set of forces acting on a given body while the other one contains the right set of forces plus at least one erroneous force; see \textbf{Fig.} \ref{fig_1} for an illustration\footnote{This question is not part of the actual FCI used in this study. It is presented here to illustrate the underlying idea (within the context of this work) about PFI questions. If it may be of help to limit the proliferation of FCI questions for easy student access, we do not quote here actual FCI questions.}. The inclusion of this erroneous force confuses students causing the ``polarization''\footnote{The word polarization is inspired by charge polarization in Physics. A charged body brought near a conducting rod will make one end of this rod positively charged while the other end negatively charged. In a sense, the set of PFI questions ``induces'' this polarization in student response: positive for right answer and negative for wrong answer in the pair of ``polarizing'' choices.} of student response.
\end{multicols}

\begin{figure}[htb!]
    \centering
    \includegraphics[scale=1.0]{sample_pol_quest.pdf}
    \caption{\it In this sample question, one is asked to identify the force(s) acting on the sphere. Whereas most students may be able to identify forces i and ii correctly, the addition of erroneous force iii may confused them, causing ``polarization'' of student answers between choices (D) and (E).}
    \label{fig_1}
\end{figure}

\begin{multicols}{2}
Dual process theory is especially appealing because of the wide gamut of its applicability from decision-making process in dealing with usual day-to-day encounters to the way we address challenging problems of academic nature. To the extent that CRT could be a link between dual process theory as a possible good measure of the inclination to shift from System 1 to System 2, and the process of learning Physics remains to be further investigated. The pioneering study for the investigation of this link for DPT, CRT, and FCI in Ref. \cite{Wood} is just a beginning. While this work is not a replica of the mentioned study, as part of this paper, we wish to further examine the interconnections among DPT, CRT, and FCI as a whole. In addition to this, we want to probe the set of PFI questions in the FCI, first identified in Ref. \cite{Alinea}, as it relates to CRT and DPT.

This paper is organized as follows. In the following section (Sec. \ref{secmethod}), we summarize the methodology on data gathering and simple statistical analysis. Following this is the results and discussion in Sec. \ref{secresults}. This section is divided into two subsections for separate discussions of the results for the FCI as a whole and CRT, and then for the PFI questions and CRT. In the last section, we summarize our findings and suggest future prospects in line with this study.

\section{Methodology}
\label{secmethod}
The subjects of our study were 163 students under engineering, chemistry, and physics degree programs offered by the University of the Philippines Los Ba\~nos. All of these students have already passed the prescribed calculus-based mechanics course and were taking other fundamental physics course under the author of this work,  upon sitting for the FCI and CRT.\footnote{Because our subjects have already passed a mechanics course, we only needed to administer one FCI. This effectively allowed us to avoid some challenges in explaining the possible shift of correlation between scores in CRT and FCI from pretest to post test, and in analyzing the correlation between normalized gain \cite{Hake} and CRT scores. While these concerns are interesting in their own right, we leave it for further studies.} The two tests were administered as pen-and-paper test and students were given incentives to take the tests seriously. Data gathering took place at the University of the Philippines Los Ba\~nos from Academic Years 2017-18 to 2018-2019.

To minimize possible bias with the CRT scores, we opted to remove results for students who did not take the test for the first time.\footnote{There was no deduction in the incentives provided for students with discarded results.} The remaining 163 test results\footnote{Interested readers may see Ref. \cite{Bialek} for a study on the robustness of the correlation of CRT scores with different variables with multiple exposure to the test and Ref. \cite{Garza} for the effect of \textit{visibility} on the performance in CRT} for the FCI and CRT were subjected to item analysis with the help of ZipGrade \cite{zipgrade} and a spreadsheet application. Common statistical parameters such as mean, standard error, and Pearson coefficients were calculated to find possible meaningful relationships between scores in FCI as a whole, set of PFI questions, and CRT.

\section{Results and Discussion}
\label{secresults}

\subsection{FCI and CRT}
\label{subsecfcicrt}
Table \ref{table_crt} shows the distribution of students with respect to CRT scores ranging from zero out of three (all wrong answers) to three out of three (perfect score). About half (47\%) of the students got a perfect score in the test while only one tenth (10\%) got the lowest score. With 43\% scoring one or two, the skewed frequency distribution drives the average to 2.1/3 (closed to that of Ref. \cite{Wood}, 2.3, with similar cohorts). This is about one and a half times higher than that in Ref. \cite{Garza} covering more than 40,000 students. The difference may be due to stronger inclination of our subjects to mathematics as suggested by their chosen degree programs in engineering, chemistry, and physics, compared to the population with more varied interests investigated in Ref. \cite{Garza}. In Ref. \cite{Sinayev}, the authors found a moderately positive correlation between cognitive reflection and numeric ability.

Majority of the students who got the correct answers perform some sort of short  mathematical calculations on paper. This constitute the majority of our test subjects based on Table \ref{table_crt2}. For the bat-and-ball problem, for instance, majority of the nearly 80\% who got the correct answer, approached the problem by solving some form of algebraic equation(s); some did try trial-and-error to fit the conditions of the problem. For the lilypad problem, geometric sequence either from day 1 to day 48 or backwards from day 48, can be found in the student solutions with correct answers. For the machine problem, short solutions involving the use of ratio and proportion can be found on student papers. Admittedly, few students who got the wrong answers made some marks on paper in addition to the final answer. However, the ``seriousness'' of these scribbles are not at par with those of students who got the correct answer(s). In the lilypad problem for instance, one can find `48/2' yielding 24 days---the intuitive but wrong answer.
\end{multicols}

\begin{table}[t]
\centering
\begin{tabular}{lccccl}
\hline\hline
\small CRT Score                          & \small 0/3         & \small 1/3         & \small 2/3         & \small 3/3         & \small Mean
\\
\hline
\small Percent Students (No. of Students) & \small 10\% (16) & \small 22\% (36) & \small 21\% (34) & \small 47\% (77) & \multicolumn{1}{c}{\small 2.1} \\
\hline\hline
\end{tabular}
\caption{Percent (and number) of students who got CRT scores from 0/3 (all wrong) to 3/3 (perfect score).}
\label{table_crt}
\end{table}

\bigskip

\begin{table}[t]
\begin{center}
\begin{tabular}{lccccl}
\hline\hline
\small CRT Problem                          & \small bat and ball         & \small machine         & \small lilypad
\\
\hline
\small Percent Students (No. of Students) & \small 79\% (128) & \small 59\% (96) & \small 68\% (111) \\
\hline\hline
\end{tabular}
\end{center}
\caption{Percent (and number) of students who answered each CRT \cite{Frederick} question correctly. Problems 1, 2, and 3 are labeled \textit{bat and ball}, \textit{machine}, and \textit{lilypad}, respectively.}
\label{table_crt2}
\end{table}

\begin{multicols}{2}
It is worth noting that students were not required to present a mathematical solution nor any other form of solution to the CRT; only the final answers were required in the closely supervised test. This has been made clear in the instructions explained in class before students could start the tests. Any mathematical solution (be it serious or simply light scribbles) in the form of algebraic manipulation, ratio and proportion, or geometric sequence, is out of their own volition. We find that the majority of students who got the wrong answer in any of the CRT problems simply gave the wrong intuitive answers (i.e., 10 cents, 100 minutes, and 24 days for the bat-and-ball, machine, and lilypad problems respectively) without any supplemental solution at all, or tried to perform some light mathematical scribbles on paper only to end up with or confirm the same wrong intuitive answers; e.g., writing 1:1:1 for the machine problem then yielding the intuitive but wrong answer of 100 minutes. The intuitive answers account for 86\%, 69\%, and 54\% of all the wrong answers in the bat-and-ball, machine, and lilypad problems respectively; this is the majority of wrong answers as in Ref. \cite{Wood}. On the other hand, among the correct student responses are clear cases with erasure marks covering the wrong intuitive answers right beside the correct one. This indicates the transition from System 1 to System 2, consistent with the aim of CRT.

Having said this, there is a possibility that one could have actually activated System 2, as may be evident in their scratch works, for a given CRT question, without finding the correct answer nor the intuitive but wrong answer.\setcounter{footnote}{0}\footnote{This could be avoided by using an online or computer-based administration of CRT. However, some students who are accustomed to doing some scratch works on paper may find some level of discomfort in getting the correct answers even on activation of System 2.} This may then weigh against the accuracy of CRT score as a measure of the tendency to activate System 2. Assuming that students who got the correct answers in the CRT activated System 2 and those who answered the intuitive but wrong answers (10 cents, 100 minutes, 24 hours) activated only System 1, this could possibly account for a maximum of 14\%, 31\%, and 46\% of all the wrong answers in the bat-and-ball, machine, and lilypad problems, respectively. However, consistent with the immediately preceding two paragraphs, we find this group of students to be somewhat of a polar opposite to the group of students who got the correct answer(s)---most students who got the wrong but non-intuitive answers did not write any form of solution at all. Some scribbled one-liner numeric calculation (e.g., `100/5' for the machine problem and `48/4' for the lilypad problem) which may hardly be considered as good evidence of activating System 2. For the three CRT problems we find only three cases (two for the bat-and-ball problem, one for the machine problem, and zero for the lilypad problem) with some tractable short algebraic calculation corresponding to wrong and non-intuitive answers. This gives us confidence to hold on to the raw CRT score as our measure of the level of cognitive reflection at least as far the scope of this study is concerned.

After having discussed the CRT results, let us now look into the student performance in the FCI with respect to CRT. Figure \ref{fig_2} shows the distribution of student scores for the two tests. For CRT scores ranging from 1 to 3, the FCI mean score exhibits an upward trend suggestive of a good linear correlation with the former for the mentioned range. However, the FCI mean score for the CRT score of zero, in addition to the errors associated with the mean, somewhat smears this prospect for a possible good positive linear correlation. The Pearson coefficient for the CRT-FCI scores turns out to be 0.24 in the approximate range [0.09, 0.38]\footnote{It is desirable to have this range of \textit{r} as narrow as possible. We approximated that given the same base value of the correlation coefficient, we may need a sample size of $\gtrsim 5000$.} (computed through the Fisher transformation) at 95\% confidence interval. Assuming the same normality condition, this is consistent with the results in Ref. \cite{Wood} for the pretest, $r = 0.38$ in the approximate range [0.23, 0.51] and post test, $r = 0.32$ in the approximate range [0.17, 0.46]\footnote{In the calculation of the approximate interval for the correlation coefficient in Ref. \cite{Wood}, we made used of the based values of this coefficient and the sample size 148 students provided in the reference.}, at the same confidence interval; see Ref. \cite{pearsoncomp} about comparing correlation coefficients.

We provide a two-fold explanation for this finding. Firstly, if CRT is a good measure of the tendency to shift from System 1 to System 2, within the context of DPT, it means that we cannot rule out the possibility that one of the contributing factors in gaining a high score in the FCI is the high level of cognitive reflection. The questions in the FCI each have four out of five incorrect answers in the choices, many of which invoke student misconception or misleading preconception \cite{halloun,finegold} in fundamental classical mechanics. Many of these misconceptions or preconceptions in turn, are from ``common-sense'' everyday experiences that constitute their intuition. For instance, students may have a predisposition to decide right away from ``common sense'' that heavier objects tend to fall faster than a lighter object of the same size and shape. A low level of cognitive reflection can drive a student to select intuitive answers that are wrong. This may then lead, as one of the contributing factors, to low FCI score.

Secondly, while a shift from System 1 to System 2 may contribute in the more ``serious'' drive to find the correct answer in one FCI question, this may not be enough. After System 2 is activated, some incorrect choices may be discarded after a short but relatively deeper reflection compared to System 1. However, there could remain further hurdles---the other incorrect answers---before a student could pinpoint the correct answer. This means that even if all students start out with activating System 2 without necessarily passing System 1 (similar to that when encountering a math problem to find $\sqrt {365.25}$), they could still end up with the wrong answer, and the student scores could still significantly vary. With the recognition that cognitive reflection as one of the contributing factors in narrowing down the set of correct answers in the FCI, we believe that other factors (e.g., industry measured by study hours, attention span, natural talent, etc.) linked to sufficiently deep understanding of physics also play a significant role in pinpointing the correct answers in the FCI or any other physics tests for that matter.
\end{multicols}

\begin{figure}[htb!]
    \centering

    \includegraphics[scale=1.0]{fci_n_crt.pdf}
    \caption{\it Distribution of FCI scores with respect to CRT scores. The error bars are set one standard error of the mean above and below the mean.}
    \label{fig_2}
\end{figure}

\begin{multicols}{2}
\subsection{PFI Questions and CRT}
\label{polcrt}
We identified and elaborated the six PFI questions in the FCI in our earlier study presented in Ref. \cite{Alinea}. Back then our sample size was only about 50 (international) students. And although relatively small, the pattern for the PFI questions was sharp enough worthy of publication. Figure \ref{fig_3} confirms the existence of these PFI questions, now with a much larger sample size of 163 students---triple that of the former study. As can be seen in the figure, majority of the students, for the six force identification questions, effectively answer only two (polarizing choices) out of the five choices. Except possibly for the identified PFI question 3, where the number of students who chose letter B is similar\footnote{Further study with larger sample size should resolve the difference between the number of students who elected choices B and D.} to that for letter D (the correct answer), the six bar graphs are consistent with the result in Ref. \cite{Alinea}.
\end{multicols}

\begin{figure}[htb!]
    \centering

    \includegraphics[scale=1.0]{pol_graphs_v2p20.pdf}
    \caption{\it Distribution of scores for the six PFI questions (Questions 5, 11, 13, 18, 29 and 30 in the FCI) with respect to CRT scores. Of the five choices, with X meaning no answer, majority of the students chose the two polarizing choices---one that contains the correct answer and the other one that is effectively a superset of the correct answer.}
    \label{fig_3}
\end{figure}

\begin{multicols}{2}
With the PFI questions at hand, let us look into its possible relationship with the CRT. Figure \ref{fig_4} shows the distribution of student scores for the set of PFI questions with respect to the CRT score. The graph shown is effectively a subset of that shown in Fig. \ref{fig_2} involving all the questions in the FCI. For the PFI questions with respect to CRT, the error bars are wider. There seems to be a rising trend between PFI mean score and CRT score from CRT score of 0 up to 2. However, the downshift at CRT score equal to 3 seems to have spoiled this possible trend. All in all, we get a correlation coefficient of $r = 0.096$ in the approximate range [-0.06, 0.25] at 95\% confidence interval.

The base value of the Pearson coefficient is small but its range tells us that we cannot simply set aside cognitive reflection in view of student performance in the PFI questions. The identification of forces acting on a given body, as simple as it may sound, is one of the most basic skills necessary in the study of dynamics. Yet even the acquisition of this very basic skill is plagued by misconceptions or misleading preconceptions directly or indirectly from everyday experiences that form our intuition. Instances of these concepts include the requirement for a force to sustain the motion of a body (in a vacuum) and the existence of centrifugal force (in an inertial frame). When confronted with questions asking for a set of forces acting on a given body, low level of cognitive reflection may lead to the inclusion of intuitive but misleading forces. Deep thinkers on the other hand, may rule out these wrong forces, effectively narrowing down the set of prospective correct answers.
\end{multicols}

\begin{figure}[htb!]
    \centering

    \includegraphics[scale=1.0]{pol_fci_n_crt_v1p05.pdf}
    \caption{\it Distribution of scores for the set of PFI questions in the FCI with respect to CRT scores. The error bars are set one standard error of the mean above and below the mean.}
    \label{fig_4}
\end{figure}

\begin{multicols}{2}
Having said this, our range of the Pearson coefficient is evidently on the lower half of its absolute spectrum from 0 to 1. Similar to that for the FCI as a whole, we contend that on top of the ability to shift from System 1 to System 2, is the need for deep understanding of Physics concepts or ideas to pinpoint the right set of forces acting on a given body. In other words, we see that high level of cognitive reflection is some sort of initial ``push'' needed to submerge one into a sea of thought, and complemented by will, natural talent, and/or industry, may lead them to acquire the right understanding or decision in solving physics problems be it as basic as force identification or as complex flying a real rocket.

Before we leave this sub-section, we take cognizance of the downwshifted range for the Pearson coefficient for PFI score vs. CRT score, with respect to FCI score (as a whole) vs. CRT score; that is, from [0.09, 0.38] to [-0.06, 0.25]. Although the intervals are still overlapping, we take our freedom to account for the possible difference (hopefully to be resolved in future studies). Figure \ref{fig_5} shows the distribution of the number of students who got the correct answer for each FCI question. Based on the figure, our students found FCI item numbers 5, 11, 13, 14, 17, 18, 21, 25, 26, and 30, to be the top 10 most challenging FCI questions. Of the six-item set of PFI questions, five belongs to these top 10 questions. We may see then that on average, the subjects of our study found the set of PFI questions to be more challenging compared to the rest of the FCI questions and this possibly caused the downward shift in the correlation coefficient. We are inclined to think that cognitive reflection stands as an important contributing factor in the student performance in the set of PFI questions and in the FCI as a whole. When students activate System 2, this is where the other factors associated with deep understanding of physics come into play. The gravity of these other factors become more evident with the level of difficulty faced by students in answering a given set of physics problems.
\end{multicols}

\begin{figure}[htb!]
    \centering

    \includegraphics[scale=0.9]{distri_correct_ans_v1p05.pdf}
    \caption{\it Distribution of the number of students who got the correct answer for each question in the FCI. The top 10 (out of 30) most challenging problems based on student performance are identified in graph as the `lowest 10'. Of these 10 questions, five belong to the set of PFI questions.}
    \label{fig_5}
\end{figure}

\begin{multicols}{2}
\section{Concluding Remarks and Future Prospects}
\label{seconclusion}
The operation of the human mind is still one of most complex processes far from our complete comprehension. From the perspective of a Physics Educator, there is a need to understand it so as to optimize student learning experience. But as the goal post of complete understanding is still far beyond the horizon, we are delighted of large strides leading to this end. Dual Process Theory may be seen as one these strides telling us of a significant simplification employing two systems of thought processes: one that is intuitive (System 1) and the other one being analytic (System 2). The Cognitive Reflection Test is a good measuring instrument of cognitive reflection indicating the tendency to shift from System 1 to System 2. The use of CRT offers an attractive avenue to look into the relationship between cognitive reflection and student performance in Physics tests such as the FCI.

In this study we look into the idea that students who can easily transition from System 1 to System 2 or start right away with System 2 in DPT, may be able to perform better in the set of polarizing force-identification questions in the FCI and in the FCI as a whole. The result of our analysis of tests involving 163 students is suggestive that cognitive reflection is one of the contributing factors involved in student performance in the FCI and the set of PFI questions. Our insight is that high level of cognitive reflection may enable students to cross out seemingly correct choices in the test that are normally part of our intuition from everyday experience or ``common sense'' but otherwise misleading. Complementary to cognitive reflection are other factors associated with deep understanding in physics. These become more evident with the level of difficulty of a given set of physics problems. We find that the range of correlation coefficient between CRT score and PFI score is downshifted with respect to that for CRT score and FCI score. The possible difference may be attributable to the set of PFI questions being more challenging on average compared the rest of FCI questions.

Looking ahead, with all its efforts, this study has only covered a small (but significant) portion of student learning and performance in physics in relation to cognitive reflection and DPT. From here, we foresee future studies involving further tests hopefully with larger number of participants. Other Physics inventories or tests related to scientific reasoning ability (see Refs. \cite{coletta,ates}) may be explored to find correlations with the level of cognitive reflection. Considering CRT, a higher resolution test (see the proposed expanded CRT in Ref. \cite{ToplakME})  in tandem with an expanded PFI questionnaire, may be used to better resolve differences in the level of cognitive reflection. Regarding the administration of the CRT, the use of ``mild subterfuge'' as in Ref. \cite{Wood} may be employed and studied to see its effect on the correlation between cognitive reflection and FCI or any other standardized Physics Test.


Effect of Cosmic UV Background on Star Formation in High Redshift Galaxies

Abstract: The effect of the cosmic UV background on star formation in high redshift galaxies is explored by means of high resolutions cosmological simulations. The simulations include star formation, 3D radiative transfer, and a highly detailed ISM model, and reach spatial resolution sufficient to resolve formation sites for molecular clouds. In the simulations the local radiation field in the Lyman-Werner band around star-forming molecular clouds dominates over the cosmic UV background by a factor of 100, similarly to the interstellar radiation field in the Milky Way and in a few high redshift galaxies for which measurements exist. The cosmic UV background, therefore, is essentially irrelevant for star formation in normal galaxies.

\section{Introduction}
\label{sec:intro}

Since the end of cosmic reionization, the intergalactic space has been
filled with the accumulated ultra-violet radiation from the previous
generations of massive stars and quasars, the so-called cosmic UV
background. The UV background completely controls the ionization
state of the intergalactic medium (IGM) that manifests itself in the
numerous Ly-$\alpha$ absorption lines in the spectra of distant
quasars, the Ly-$\alpha$ forest \citep[see][for a recent review]{igm:m09}.

In the interstellar medium of local galaxies (ISM) the situation
appears to be reversed: the Haardt-Madau model for the UV background
\citep{jnu:hm01} predicts the $z=0$ radiation field at $1000\AA$ of
about
$2\times10^3\dim{photons}/\dim{cm}^2/\dim{s}/\dim{ster}/\dim{eV}$,
some 500 smaller than the interstellar radiation field in the Milky
Way at the solar circle \citep{h2:d78,h2:mmp83}.

At higher redshift the situation is, however, less clear. For example,
the UV background at $z\sim3$ is expected to be 50 to 100 times higher
than at $z=0$ \citep{jnu:hm01,jnu:flzh09}. On the other hand,
estimates of the interstellar radiation field in high redshift
gamma-ray burst hosts also give some 100 times higher values
than the Milky Way field \citep{hizgal:cppp09}, but the observational
constraints remain sparse.

In the cosmological simulation community, the effect of the cosmic UV
background is often included in the cooling rates
\citep[e.g.][etc]{sims:co92,sims:kwh96,sims:ns97,sims:k03,sims:ck09,sims:svbw09,sims:honj09}. However,
the effect of the local interstellar radiation field (that dominates
over the cosmic UV background by a larger factor at least at $z=0$,
and may be at high redshift too) has not yet been included in those
simulations. In this paper the effect of the cosmic UV background is
critically reassessed with numerical simulations that both include the
full 3D, time-dependent and spatially variable treatment of radiative
transfer and have high enough spatial resolution to resolve the sites
of molecular clouds and associated star formation.

\section{Simulations and star formation model}
\label{sec:sims}

The physical ingredients and computational setup for the simulations
used in this paper have been recently described in great detail
elsewhere \citep{ng:gk10a,ng:gk10b}. As a brief reminder, the
simulations have been performed with the Adaptive Refinement Tree (ART) code
\citep{misc:k99,misc:kkh02,sims:rzk08} that uses adaptive mesh
refinement in both the gas dynamics and gravity calculations to
achieve high dynamic range in spatial scale.

The simulations include star formation and supernova enrichment and
thermal energy feedback, as well as a highly detailed ISM model. The
3D radiative transfer of UV radiation from individual stellar
particles is followed self-consistently with the OTVET approximation
\citep{ng:ga01}. The simulations incorporate non-equilibrium chemical
network of hydrogen and helium and non-equilibrium,
metallicity-dependent cooling and heating rates, and a
phenomenological model of molecular hydrogen formation on and
shielding by cosmic dust, as well as ${{\rm H}_2}$ self-shielding
\citep{ng:gk10a,ng:gk10b}.

Particular simulations used in this paper model a small region
including a Milky-Way progenitor galaxy and a number of smaller
galaxies with the mass resolution of $1.3\times10^6\, {\rm M}_{\odot}$ in dark
matter, $2.2\times 10^5\, {\rm M}_{\odot}$ in baryons, and with the spatial resolution of
$65\dim{pc}\times[4/(1+z)]$ (in physical units) within the fully
refined region.

Star formation in the simulation is occurring in the molecular gas
only, using the prescriptions of \citet{sfr:km05} and
\citet{sfr:kt07}. The exact formulation of the star formation recipe
is shown in Equation (2) of \citet{ng:gk10b}.

For the purpose of this paper, two simulations are considered that
differ only by the inclusion of the cosmic UV background from
\citet{jnu:hm01}. In the simulation without the cosmic UV background,
the local radiation field produced by nearby massive stars is still
included in exactly the same manner as in the simulation with the UV
background. Thus, these two simulations can be used to evaluate the
particular effect of the cosmic UV background on the properties of
model galaxies. Because of computational expense, the simulations are
not continued beyond $z=2$.

The simulation with the cosmic UV background is the same one as
described in \citet{ng:gk10a}. That reference also shows good agreement of
that simulation with several observational constraints on the
properties of high redshift galaxies.

\section{Results}
\label{sec:res}

\begin{figure}[t]
\includegraphics[scale=0.43]{\figname{sfrmsmt}{f1.eps}}
\caption{The stellar mass (top) and the star formation rate (bottom)
  vs the total mass for all fully resolved model galaxies at
  $z=3$. Open squares show the run without the cosmic UV background,
  while filled squares show the fully self-consistent run, with the UV
  background.  The star formation rates are averaged over
  $20\dim{Myr}$.\label{fig:gals}}
\end{figure}

A direct comparison of stellar masses and star formation rates between
the galaxies in the two simulations is shown in Figure
\ref{fig:gals}. In order to minimize the effect of finite numerical
resolution, only highly resolved galaxies - i.e.\ galaxies that reach
the 8th level of mesh refinement ($130\dim{pc}$ spatial resolution at $z=3$)
within their gaseous disks - are shown in this and all subsequent
figures. At this spatial resolution the sub-cell model for ${{\rm H}_2}$
formation performs reliably, as is demonstrated in Figure 13 of
\citet{ng:gk10b}. Galaxies with masses below $M_{\rm tot}=10^{10}\, {\rm M}_{\odot}$ are not
sufficiently resolved in these simulations - none of such galaxies
achieves 8 levels of mesh refinement within their gaseous disks.

While the two simulations do not produce identical
results, the difference between the two runs is not dramatic and is
fully consistent with the timing differences in two simulations that
have slightly different time-steps and, hence, times of intermediate
outputs. Thus, no significant effect of the cosmic UV background on
the global properties of simulated galaxies is observed in these
simulations.

\begin{figure}[t]
\includegraphics[scale=0.43]{\figname{pzcomp}{f2.eps}}
\caption{Comparison between the proximity zones (filled diamonds) and
  radii containing 90\% of the neutral (top) and molecular (bottom)
  gas in model galaxies at $z=3$ (open squares). Vertical error-bars
  give 10\% - 90\% range for the proximity zone sizes as viewed in
  different directions from the center of each galaxy.
  \label{fig:pzcomp}}
\end{figure}

The origin for this conclusion becomes apparent from Figure
\ref{fig:pzcomp}, where the sizes of the proximity zones -
regions around model galaxies that are dominated by the local
radiation from the galaxies themselves and not by the cosmic UV
background - are shown together with radii containing most of atomic
and molecular gas.

A size of a proximity zone depends, of course, on the wavelength of
radiation, because, in general, spectral shapes of the local
interstellar radiation field and the cosmic UV background are
different. Figure \ref{fig:pzcomp}, therefore, shows two important
wavelengths: the Lyman limit and $\lambda=1000\AA$, in the middle of
the Lyman-Werner band. Radiation above the Lyman limit ionizes neutral
hydrogen, and so is important for determining which part of the
galactic gas can cool rapidly. Radiation in the Lyman-Werner band
destroys molecular hydrogen; only when gas is shielded from that
radiation by cosmic dust and molecular hydrogen self-shielding can it
become fully molecular (and, hence form stars).

In a given frequency band, absorption of radiation emitted by stars in
a galaxy will, generally, be different in different
directions. Therefore, the proximity zone of a given galaxy is not
necessarily spherical, but can have a complex shape.  In order to
quantify variation in the shapes of proximity zones, a HEALPix
tessellation \citep{misc:ghbw05} of the celestial sphere (as seen from
the center of a particular galaxy) is constructed and the proximity
zone size is measured separately for each pixel on the sky. Because
HEALPix provides uniform sampling of all possible directions, it
offers a convenient way to quantify variations in sizes of proximity
zones in different wavebands. Error-bars in Figure \ref{fig:pzcomp}
show the 10\% - 90\% range in the size of the proximity zone as a
function of the angle on the sky.

To evaluate the effect of the local radiation field on the gas
physics, Figure \ref{fig:pzcomp} also displays radii that contain 90\%
of atomic (top) and molecular (bottom) gas in the model galaxies. As
can be seen, the proximity zone in the Lyman-Werner band extends over
10 times beyond the edge of the molecular gas. That explains why the
star formation in the simulations is insensitive to the presence of
the cosmic UV background - the local radiation field in the sites of
star formation (i.e.\ in molecular clouds) always dominates over the
cosmic background by a factor of about 100.

The situation appears to be reversed for atomic hydrogen. The escape
fraction for ionizing radiation is small in these simulations,
consistent with observational estimates \citep{ng:gkc08}. It is
particularly small along the disk, so that little local ionizing
radiation shines on the outer parts of ${\rm H\,I}$ disks; instead, the edges
of ${\rm H\,I}$ disks are determined by the cosmic background.

\begin{figure}[t]
\includegraphics[scale=0.43]{\figname{pzevol}{f3.eps}}
\caption{Comparison between the proximity zones and radii containing
  90\% of the neutral and molecular gas for the fully self-consistent
  run, with the UV background, at three different redshifts. Symbol
  and panels are similar to Figure \ref{fig:pzcomp}. Three shades of
  gray show three different redshifts, as marked in the figure.
\label{fig:pzevol}}
\end{figure}

Finally, to verify the robustness of these results,  Figure
\ref{fig:pzevol} presents proximity zones and 90\%-mass radii for the full
self-consistent simulation (with the cosmic UV background) at 3 different
redshifts. While there is substantial scatter between individual
points, the trends observed at $z=3$ are reproduced at other redshifts
as well.

\section{Conclusions}
\label{sec:discussion}

The results of this paper can be summarized in just a few lines: the
local radiation field in the Lyman-Werner band around star-forming
molecular clouds dominates over the cosmic UV background by a factor
of 100, similarly to the interstellar radiation field in the Milky Way
\citep{h2:d78,h2:mmp83} and in a few high redshift galaxies for which
measurements exist \citep{hizgal:cppp09}. The cosmic UV background,
therefore, is essentially irrelevant for star formation in normal
galaxies. Only in the lowest mass dwarfs, where the radiation can
actually photo-evaporate the ISM, may the cosmic background affect
star formation. Such galaxies, however, contribute very little to the
total star formation history of the universe.

The situation is more complex for ionizing radiation - in that band
the cosmic UV background may play an important role in determining the
exact locations of the edges of extended ${\rm H\,I}$ disks, in agreement
with the earlier results of \citet{igm:s04a}. Those outside regions of
galactic disks, however, contain little molecular gas and, therefore,
are inert to star formation.

The simulations presented here do not continue beyond $z=2$, and so
can not be used for testing the results of \citet{sims:honj09}, who
found little effect of the cosmic UV background on star formation in
normal galaxies at $z>2$, but a larger effect at lower
redshifts. Simulations of \citet{sims:honj09}, however, did not
include radiative transfer and, therefore, failed to account for the
local radiation and its dominance over the cosmic UV background within
the galactic proximity zones. For example, if the Milky Way galaxy is
representative of low redshift ($z<2$) normal galaxies, then one may
expect that even at low redshifts local radiation dominates over the
cosmic background within the galactic ISM, rendering the background
irrelevant to star formation.

\acknowledgements

I am grateful to Andrey Kravtsov for comprehensive comments and
enlightening discussion. This work was supported in part by the DOE at
Fermilab, by the NSF grant AST-0908063, and by the NASA grant
NNX-09AJ54G. The simulations used in this work have been performed on
the Joint Fermilab - KICP Supercomputing Cluster, supported by grants
from Fermilab, Kavli Institute for Cosmological Physics, and the
University of Chicago.  This work made extensive use of the HEALPix
spherical tessellation package
and the NASA Astrophysics Data System and
{\tt arXiv.org} preprint server.

\bibliographystyle{apj}

Cosmology and Nongalactic Astrophysics (astro-ph.CO)