-
Type II radio bursts are evidence of shocks in the solar atmosphere and inner
heliosphere that emit radio waves ranging from sub-meter to kilometer lengths.
These shocks may be associated with CMEs and reach speeds higher than the local
magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now
possible with LOFAR, opening a new radio window in which to study coronal
shocks that leave the inner solar corona and enter the interplanetary medium
and to understand their association with CMEs. To this end, we study a coronal
shock associated with a CME and type II radio burst to determine the locations
at which the radio emission is generated, and we investigate the origin of the
band-splitting phenomenon.
-
Context. The Sun is an active star that produces large-scale energetic events
such as solar flares and coronal mass ejections and numerous smaller-scale
events such as solar jets. These events are often associated with accelerated
particles that can cause emission at radio wavelengths. The reconfiguration of
the solar magnetic field in the corona is believed to be the cause of the
majority of solar energetic events and accelerated particles. Aims. Here, we
investigate a bright J-burst that was associated with a solar jet and the
possible emission mechanism causing these two phenomena. Methods. We used data
from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio
data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph
(NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July
2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic
components and it was associated with a solar jet observed at extreme
ultraviolet wavelengths with SDO. The solar jet occurred at a time and location
coincident with the radio burst, in the northern hemisphere, and not inside a
group of complex active regions in the southern hemisphere. The jet occurred in
the negative polarity region of an area of bipolar plage. Newly emerged
positive flux in this region appeared to be the trigger of the jet.
Conclusions. Magnetic reconnection between the overlying coronal field lines
and the newly emerged positive field lines is most likely the cause of the
solar jet. Radio imaging provides a clear association between the jet and the
J-burst which shows the path of the accelerated electrons.
-
We present a study of carbon radio recombination lines towards Cassiopeia A
using LOFAR observations in the frequency range 10-33 MHz. Individual carbon
$\alpha$ lines are detected in absorption against the continuum at frequencies
as low as 16 MHz. Stacking several C$\alpha$ lines we obtain detections in the
11-16 MHz range. These are the highest signal-to-noise measurements at these
frequencies. The peak optical depth of the C$\alpha$ lines changes considerably
over the 11-33 MHz range with the peak optical depth decreasing from
4$\times10^{-3}$ at 33 MHz to 2$\times10^{-3}$ at 11 MHz, while the line width
increases from 20 km s$^{-1}$ to 150 km s$^{-1}$. The combined change in peak
optical depth and line width results in a roughly constant integrated optical
depth. We interpret this as carbon atoms close to local thermodynamic
equilibrium.
In this work we focus on how the 11-33 MHz carbon radio recombination lines
can be used to determine the gas physical conditions. We find that the ratio of
the carbon radio recombination lines to that of the 158 $\mu$m [CII]
fine-structure line is a good thermometer, while the ratio between low
frequency carbon radio recombination lines provides a good barometer. By
combining the temperature and pressure constraints with those derived from the
line width we are able to constrain the gas properties (electron temperature
and density) and radiation field intensity. Given the 1$\sigma$ uncertainties
in our measurements these are; $T_{e}\approx68$-$98$ K,
$n_{e}\approx0.02$-$0.035$ cm$^{-3}$ and $T_{r,100}\approx1500$-$1650$ K.
Despite challenging RFI and ionospheric conditions, our work demonstrates that
observations of carbon radio recombination lines in the 10-33 MHz range can
provide insight into the gas conditions.
-
We use the Low Frequency Array to perform a systematic high spectral
resolution investigation of the low-frequency 33-78 MHz spectrum along the line
of sight to Cassiopeia A. We complement this with a 304-386 MHz Westerbork
Synthesis Radio telescope observation. In this first paper we focus on the
carbon radio recombination lines.
We detect Cn$\alpha$ lines at -47 and -38 km s$^{-1}$ in absorption for
quantum numbers n=438-584 and in emission for n=257-278 with high signal to
noise. These lines are associated with cold clouds in the Perseus spiral arm
component. Hn$\alpha$ lines are detected in emission for n=257-278. In
addition, we also detect Cn$\alpha$ lines at 0 km s$^{-1}$ associated with the
Orion arm.
We analyze the optical depth of these transitions and their line width. Our
models show that the carbon line components in the Perseus arm are best fit
with an electron temperature 85 K and an electron density 0.04 cm$^{-3}$ and
can be constrained to within 15\%. The electron pressure is constrained to
within 20\%. We argue that much of these carbon radio recombination lines arise
in the CO-dark surface layers of molecular clouds where most of the carbon is
ionized but hydrogen has made the transition from atomic to molecular. The
hydrogen lines are clearly associated with the carbon line emitting clouds, but
the low-frequency upperlimits indicate that they likely do not trace the same
gas. Combining the hydrogen and carbon results we arrive at a firm lower limit
to the cosmic ray ionization rate of 2.5$\times$10$^{-18}$ s$^{-1}$, but the
actual value is likely much larger.
-
In the first paper of this series, we study the level population problem of
recombining carbon ions. We focus our study on high quantum numbers
anticipating observations of Carbon Radio Recombination Lines to be carried out
by the LOw Frequency ARray (LOFAR). We solve the level population equation
including angular momentum levels with updated collision rates up to high
principal quantum numbers. We derive departure coefficients by solving the
level population equation in the hydrogenic approximation and including low
temperature dielectronic recombination effects. Our results in the hydrogenic
approximation agree well with those of previous works. When comparing our
results including dielectronic recombination we find differences which we
ascribe to updates in the atomic physics (e.g., collision rates) and to the
approximate solution method of the statistical equilibrium equations adopted in
previous studies. A comparison with observations is discussed in an
accompanying article, as radiative transfer effects need to be considered.
-
In the second paper of the series, we have modeled low frequency carbon radio
recombination lines (CRRL) from the interstellar medium. Anticipating the LOw
Frequency ARray (LOFAR) survey of Galactic CRRLs, we focus our study on the
physical conditions of the diffuse cold neutral medium (CNM). We have used the
improved departure coefficients computed in the first paper of the series to
calculate line-to-continuum ratios. The results show that the line width and
integrated optical depths of CRRL are sensitive probes of the electron density,
gas temperature, and the emission measure of the cloud. Furthermore, the ratio
of CRRL to the [CII] at 158 $\mu$m line is a strong function of the temperature
and density of diffuse clouds. Guided by our calculations, we analyze CRRL
observations and illustrate their use with data from the literature.
-
We present the first wide area (19 deg$^2$), deep ($\approx120-150$ {\mu}Jy
beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna
image of the Bo\"otes field made at 130-169 MHz. This image is at least an
order of magnitude deeper and 3-5 times higher in angular resolution than
previously achieved for this field at low frequencies. The observations and
data reduction, which includes full direction-dependent calibration, are
described here. We present a radio source catalogue containing 6276 sources
detected over an area of $19$\,deg$^2$, with a peak flux density threshold of
$5\sigma$. As the first thorough test of the facet calibration strategy,
introduced by van Weeren et al., we investigate the flux and positional
accuracy of the catalogue. We present differential source counts that reach an
order of magnitude deeper in flux density than previously achieved at these low
frequencies, and show flattening at 150 MHz flux densities below 10 mJy
associated with the rise of the low flux density star-forming galaxies and
radio-quiet AGN.
-
Air showers induced by cosmic rays create nanosecond pulses detectable at
radio frequencies. These pulses have been measured successfully in the past few
years at the LOw Frequency ARray (LOFAR) and are used to study the properties
of cosmic rays. For a complete understanding of this phenomenon and the
underlying physical processes, an absolute calibration of the detecting antenna
system is needed. We present three approaches that were used to check and
improve the antenna model of LOFAR and to provide an absolute calibration of
the whole system for air shower measurements. Two methods are based on
calibrated reference sources and one on a calibration approach using the
diffuse radio emission of the Galaxy, optimized for short data-sets. An
accuracy of 19% in amplitude is reached. The absolute calibration is also
compared to predictions from air shower simulations. These results are used to
set an absolute energy scale for air shower measurements and can be used as a
basis for an absolute scale for the measurement of astronomical transients with
LOFAR.
-
Context. Observing Jupiter's synchrotron emission from the Earth remains
today the sole method to scrutinize the distribution and dynamical behavior of
the ultra energetic electrons magnetically trapped around the planet (because
in-situ particle data are limited in the inner magnetosphere). Aims. We perform
the first resolved and low-frequency imaging of the synchrotron emission with
LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV)
which map a broad region of Jupiter's inner magnetosphere. Methods (see article
for complete abstract) Results. The first resolved images of Jupiter's
radiation belts at 127-172 MHz are obtained along with total integrated flux
densities. They are compared with previous observations at higher frequencies
and show a larger extent of the synchrotron emission source (>=4 $R_J$). The
asymmetry and the dynamic of east-west emission peaks are measured and the
presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux
density measurements are on the low side of previous (unresolved) ones,
suggesting a low-frequency turnover and/or time variations of the emission
spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The
observations at 127 MHz depict an extended emission up to ~4-5 planetary radii.
The similarities with high frequency results reinforce the conclusion that: i)
the magnetic field morphology primarily shapes the brightness distribution of
the emission and ii) the radiating electrons are likely radially and
latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent
of the brightness combined with the overall lower flux density, yields new
information on Jupiter's electron distribution, that may shed light on the
origin and mode of transport of these particles.
-
We present the results of a four-month campaign searching for low-frequency
radio transients near the North Celestial Pole with the Low-Frequency Array
(LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data
were recorded between 2011 December and 2012 April and comprised 2149 11-minute
snapshots, each covering 175 deg^2. We have found one convincing candidate
astrophysical transient, with a duration of a few minutes and a flux density at
60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or
high-energy counterpart, as a result of which its nature is unclear. The
detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7)
x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2,
at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data
were also searched for transients at a range of other time-scales, from 0.5 to
297 min, which allowed us to place a range of limits on transient rates at 60
MHz as a function of observation duration.
-
LOFAR offers the unique capability of observing pulsars across the 10-240 MHz
frequency range with a fractional bandwidth of roughly 50%. This spectral range
is well-suited for studying the frequency evolution of pulse profile morphology
caused by both intrinsic and extrinsic effects: such as changing emission
altitude in the pulsar magnetosphere or scatter broadening by the interstellar
medium, respectively. The magnitude of most of these effects increases rapidly
towards low frequencies. LOFAR can thus address a number of open questions
about the nature of radio pulsar emission and its propagation through the
interstellar medium. We present the average pulse profiles of 100 pulsars
observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100
profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with
Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations
at higher frequencies (350 and1400 MHz) in order to study the profile
evolution. The profiles are aligned in absolute phase by folding with a new set
of timing solutions from the Lovell Telescope, which we present along with
precise dispersion measures obtained with LOFAR. We find that the profile
evolution with decreasing radio frequency does not follow a specific trend but,
depending on the geometry of the pulsar, new components can enter into, or be
hidden from, view. Nonetheless, in general our observations confirm the
widening of pulsar profiles at low frequencies, as expected from
radius-to-frequency mapping or birefringence theories. We offer this catalog of
low-frequency pulsar profiles in a user friendly way via the EPN Database of
Pulsar Profiles (http://www.epta.eu.org/epndb/).
-
Context. The Sun is an active source of radio emission that is often
associated with energetic phenomena ranging from nanoflares to coronal mass
ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond
duration radio bursts have been reported, such as radio spikes or solar S
bursts (where S stands for short). To date, these have neither been studied
extensively nor imaged because of the instrumental limitations of previous
radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were
used to study the spectral and spatial characteristics of a multitude of S
bursts, as well as their origin and possible emission mechanisms. Methods. We
used 170 simultaneous tied-array beams for spectroscopy and imaging of S
bursts. Since S bursts have short timescales and fine frequency structures,
high cadence (~50 ms) tied-array images were used instead of standard
interferometric imaging, that is currently limited to one image per second.
Results. On 9 July 2013, over 3000 S bursts were observed over a time period of
~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and
narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a
wide range of circular polarisation degrees (2-8 times more polarised than the
accompanying Type III bursts). Extrapolation of the photospheric magnetic field
using the potential field source surface (PFSS) model suggests that S bursts
are associated with a trans-equatorial loop system that connects an active
region in the southern hemisphere to a bipolar region of plage in the northern
hemisphere. Conclusions. We have identified polarised, short-lived solar radio
bursts that have never been imaged before. They are observed at a height and
frequency range where plasma emission is the dominant emission mechanism,
however they possess some of the characteristics of electron-cyclotron maser
emission.
-
The Square Kilometre Array (SKA) is called to revolutionise essentially all
areas of Astrophysics. With a collecting area of about a square kilometre, the
SKA will be a transformational instrument, and its scientific potential will go
beyond the interests of astronomers. Its technological challenges and huge cost
requires a multinational effort, and Europe has recognised this by putting the
SKA on the roadmap of the European Strategy Forum for Research Infrastructures
(ESFRI). The Spanish SKA White Book is the result of the coordinated effort of
120 astronomers from 40 different research centers. The book shows the enormous
scientific interest of the Spanish astronomical community in the SKA and
warrants an optimum scientific exploitation of the SKA by Spanish researchers,
if Spain enters the SKA project.
-
Using observations obtained with the LOw Fequency ARray (LOFAR), the
Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA)
data, we have traced the radio emission to large scales in the complex source
4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial
resolution (~4"), the source was known to have two inner radio lobes spanning
31 kpc and a diffuse, low-brightness extension running parallel to them, offset
by about 11 kpc (in projection).
At 62 MHz, we detect the radio emission of this structure extending out to
210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure
appears to have a helical morphology.
We have derived the characteristics of the radio spectral index across the
source. We show that the source morphology is most likely the result of at
least two episodes of AGN activity separated by a dormant period of around 35
Myr.
The AGN is hosted by one of the galaxies located in the cluster core of Abell
407. We propose that it is intermittently active as it moves in the dense
environment in the cluster core. Using LOFAR, we can trace the relic plasma
from that episode of activity out to greater distances from the core than ever
before.
Using the the WSRT, we detect HI in absorption against the center of the
radio source. The absorption profile is relatively broad (FWHM of 288 km/s),
similar to what is found in other clusters.
Understanding the duty cycle of the radio emission as well as the triggering
mechanism for starting (or restarting) the radio-loud activity can provide
important constraints to quantify the impact of AGN feedback on galaxy
evolution. The study of these mechanisms at low frequencies using morphological
and spectral information promises to bring new important insights in this
field.
-
We present measurements of radio emission from cosmic ray air showers that
took place during thunderstorms. The intensity and polarization patterns of
these air showers are radically different from those measured during
fair-weather conditions. With the use of a simple two-layer model for the
atmospheric electric field, these patterns can be well reproduced by
state-of-the-art simulation codes. This in turn provides a novel way to study
atmospheric electric fields.
-
The Square Kilometre Array (SKA) will transform our understanding of the role
of the cold, atomic gas in galaxy evolution. The interstellar medium (ISM) is
the repository of stellar ejecta and the birthsite of new stars and, hence, a
key factor in the evolution of galaxies over cosmic time. Cold, diffuse, atomic
clouds are a key component of the ISM, but so far this phase has been difficult
to study, because its main tracer, the HI 21 cm line, does not constrain the
basic physical information of the gas (e.g., temperature, density) well. The
SKA opens up the opportunity to study this component of the ISM through a
complementary tracer in the form of low-frequency (<350 MHz) carbon radio
recombination lines (CRRL). These CRRLs provide a sensitive probe of the
physical conditions in cold, diffuse clouds. The superb sensitivity, large
field of view, frequency resolution and coverage of the SKA allows for
efficient surveys of the sky, that will revolutionize the field of
low-frequency recombination line studies. By observing these lines with the SKA
we will be able determine the thermal balance, chemical enrichment, and
ionization rate of the cold, atomic medium from degree-scales down to scales
corresponding to individual clouds and filaments in our Galaxy, the Magellanic
Clouds and beyond. Furthermore, being sensitive only to the cold, atomic gas,
observations of low-frequency CRRLs with the SKA will aid in disentangling the
warm and cold constituents of the HI 21 cm emission.
-
The nuclear starburst in the nearby galaxy M82 provides an excellent
laboratory for understanding the physics of star formation. This galaxy has
been extensively observed in the past, revealing tens of radio-bright compact
objects embedded in a diffuse free-free absorbing medium. Our understanding of
the structure and physics of this medium in M82 can be greatly improved by
high-resolution images at low frequencies where the effects of free-free
absorption are most prominent.
The aims of this study are, firstly, to demonstrate imaging using
international baselines of the Low Frequency Array (LOFAR), and secondly, to
constrain low-frequency spectra of compact and diffuse emission in the central
starburst region of M82 via high-resolution radio imaging at low frequencies.
The international LOFAR telescope was used to observe M82 at 110-126MHz and
146-162MHz. Images were obtained using standard techniques from very long
baseline interferometry. images were obtained at each frequency range: one only
using international baselines, and one only using the longest Dutch (remote)
baselines.
The 154MHz image obtained using international baselines is a new imaging
record in terms of combined image resolution (0.3$"$) and sensitivity
($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects
at 154MHz, six of these also at 118MHz. Four weaker but resolved features are
also found: a linear (50pc) filament and three other resolved objects, of which
two show a clear shell structure. We do not detect any emission from either
supernova 2008iz or from the radio transient source 43.78+59.3. The images
obtained using remote baselines show diffuse emission, associated with the
outflow in M82, with reduced brightness in the region of the edge-on
star-forming disk.
-
Aims. An efficient means of locating calibrator sources for International
LOFAR is developed and used to determine the average density of usable
calibrator sources on the sky for subarcsecond observations at 140 MHz.
Methods. We used the multi-beaming capability of LOFAR to conduct a fast and
computationally inexpensive survey with the full International LOFAR array.
Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density
using existing catalogues. By observing 30 different sources in each of the 12
sets of pointings per hour, we were able to inspect 630 sources in two hours to
determine if they possess a sufficiently bright compact component to be usable
as LOFAR delay calibrators. Results. Over 40% of the observed sources are
detected on multiple baselines between international stations and 86 are
classified as satisfactory calibrators. We show that a flat low-frequency
spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz.
We extrapolate from our sample to show that the density of calibrators on the
sky that are sufficiently bright to calibrate dispersive and non-dispersive
delays for the International LOFAR using existing methods is 1.0 per square
degree. Conclusions. The observed density of satisfactory delay calibrator
sources means that observations with International LOFAR should be possible at
virtually any point in the sky, provided that a fast and efficient search using
the methodology described here is conducted prior to the observation to
identify the best calibrator.
-
LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer
constructed in the north of the Netherlands and across europe. Utilizing a
novel phased-array design, LOFAR covers the largely unexplored low-frequency
range from 10-240 MHz and provides a number of unique observing capabilities.
Spreading out from a core located near the village of Exloo in the northeast of
the Netherlands, a total of 40 LOFAR stations are nearing completion. A further
five stations have been deployed throughout Germany, and one station has been
built in each of France, Sweden, and the UK. Digital beam-forming techniques
make the LOFAR system agile and allow for rapid repointing of the telescope as
well as the potential for multiple simultaneous observations. With its dense
core array and long interferometric baselines, LOFAR achieves unparalleled
sensitivity and angular resolution in the low-frequency radio regime. The LOFAR
facilities are jointly operated by the International LOFAR Telescope (ILT)
foundation, as an observatory open to the global astronomical community. LOFAR
is one of the first radio observatories to feature automated processing
pipelines to deliver fully calibrated science products to its user community.
LOFAR's new capabilities, techniques and modus operandi make it an important
pathfinder for the Square Kilometre Array (SKA). We give an overview of the
LOFAR instrument, its major hardware and software components, and the core
science objectives that have driven its design. In addition, we present a
selection of new results from the commissioning phase of this new radio
observatory.
-
We report results from a study of the HI content and stellar properties of
nearby galaxies detected by the Arecibo Legacy Fast ALFA blind 21-cm line
survey and the Sloan Digital Sky Survey in two declination strips covering a
total area of 9 hr X 16 deg. Our analysis seeks to assemble a control sample of
galaxies suitable for providing absolute measures of the HI content of gaseous
objects. From a database of ~15,000 HI detections, we have assembled three
samples of gas-rich galaxies expected to show little or no evidence of
interaction with their surroundings that could provide adequate HI standards.
The most reliable results are obtained with a sample of 5647 sources found in
low density environments, as defined by a nearest neighbor approach. The other
two samples contain several hundred relatively isolated galaxies each, as
determined from standard isolation algorithms. We find that isolated objects
are not particularly gas-rich compared to their low-density-environment
counterparts, while they suffer from selection bias and span a smaller dynamic
range. All this makes them less suitable for defining a reference for HI
content. We have explored the optical morphology of gaseous galaxies in quiet
environments finding that, within the volume surveyed, the vast majority of
them display unequivocal late-type galaxy features. In contrast, bona fide
gas-rich early-type systems account only for a negligible fraction of the 21-cm
detections. We argue that HI emission provides the most reliable way to
determine the morphological population to which a galaxy belongs. We have also
observed that the color distribution of flux-limited samples of
optically-selected field HI emitters does not vary significantly with
increasing distance, while that of non-detections becomes notably redder. This
result suggests that the colors and HI masses of gas-rich galaxies cannot be
very closely related.
-
This is the second paper of two reporting results from a study of the HI
content and stellar properties of nearby galaxies detected by the Arecibo
Legacy Fast ALFA blind 21-cm line survey and the Sloan Digital Sky Survey in a
2160 deg^2 region covered by both surveys. We apply strategies of multivariate
data analysis to a complete HI flux-limited subset of 1624 objects extracted
from the control sample of HI emitters assembled by Toribio et al. (2011a) in
order to: i) investigate the correlation structure of the space defined by an
extensive set of observables describing gas-rich systems; ii) identify the
intrinsic parameters that best define their HI content; and iii) explore the
scaling relations arising from the joint distributions of the quantities most
strongly correlated with the HI mass. The principal component analysis
performed over a set of five galaxy properties reveals that they are strongly
interrelated, supporting previous claims that nearby HI emitters show a high
degree of correlation. The best predictors for the expected value of MHI are
the diameter of the stellar disk, D25r, followed by the total luminosity (both
in the r-band), and the maximum rotation speed, while morphological proxies
such as color show only a moderately strong correlation with the gaseous
content attenuated by observational error. The simplest and most accurate
prescription is log(MHI/Msun)= 8.72 + 1.25*log(D25r/kpc). We find a slope of
$-8.2 \pm 0.5$ for the relation between optical magnitude and log rotation
speed, in good agreement with Tully-Fisher studies, and a log slope of $1.55
\pm 0.06$ for the HI mass-optical galaxy size relation. Given the homogeneity
of the measurements and the completeness of our dataset, the latter outcome
suggests that the constancy of the average (hybrid) HI surface density
advocated by some authors for the spiral population is a crude approximation.