Far-Infrared and Sub-Millimeter Observations and Physical Models of the Reflection Nebula Ced 201

ISO [C II] 158 micron, [O I] 63 micron, and H_2 9 and 17 micron observations are presented of the reflection nebula Ced 201, which is a photon-dominated region illuminated by a B9.5 star with a color temperature of 10,000 K (a cool PDR). In combination with ground based [C I] 609 micron, CO, 13CO, CS and HCO+ data, the carbon budget and physical structure of the reflection nebula are constrained. The obtained data set is the first one to contain all important cooling lines of a cool PDR, and allows a comparison to be made with classical PDRs. To this effect one- and three-dimensional PDR models are presented which incorporate the physical characteristics of the source, and are aimed at understanding the dominant heating processes of the cloud. The contribution of very small grains to the photo-electric heating rate is estimated from these models and used to constrain the total abundance of PAHs and small grains. Observations of the pure rotational H_2 lines with ISO, in particular the S(3) line, indicate the presence of a small amount of very warm, approximately 330 K, molecular gas. This gas cannot be accommodated by the presented models.Comment: 32 pages, 7 figures, in LaTeX. To be published in Ap

Computational science and re-discovery: open-source implementations of ellipsoidal harmonics for problems in potential theory

We present two open-source (BSD) implementations of ellipsoidal harmonic expansions for solving problems of potential theory using separation of variables. Ellipsoidal harmonics are used surprisingly infrequently, considering their substantial value for problems ranging in scale from molecules to the entire solar system. In this article, we suggest two possible reasons for the paucity relative to spherical harmonics. The first is essentially historical---ellipsoidal harmonics developed during the late 19th century and early 20th, when it was found that only the lowest-order harmonics are expressible in closed form. Each higher-order term requires the solution of an eigenvalue problem, and tedious manual computation seems to have discouraged applications and theoretical studies. The second explanation is practical: even with modern computers and accurate eigenvalue algorithms, expansions in ellipsoidal harmonics are significantly more challenging to compute than those in Cartesian or spherical coordinates. The present implementations reduce the "barrier to entry" by providing an easy and free way for the community to begin using ellipsoidal harmonics in actual research. We demonstrate our implementation using the specific and physiologically crucial problem of how charged proteins interact with their environment, and ask: what other analytical tools await re-discovery in an era of inexpensive computation?Comment: 25 pages, 3 figure

The geomorphological setting of some of Scotland's east coast freshwater mills: a comment on Downward and Skinner (2005) ‘Working rivers: the geomorphological legacy...’

Many of the water mills on Scotland's east coast streams, unlike those discussed recently by Downward and Skinner (2005 Area 37 138–47), are found in predominantly bedrock reaches immediately downstream of knickpoints (i.e. bedrock steps). Bedrock knickpoints in the lower reaches of Scottish rivers are a widespread fluvial response to the glacio-isostatic rebound of northern Britain. These steps in the river profile propagate headward over time, but for intervals of a few centuries or so they are sufficiently stable to be exploited for the elevational fall necessary to power the mill wheel. Many of these mills were apparently powered by ‘run-of-the-river’, as are some today that formerly had mill dams. The typical lack of sediment storage along the erosional lower reaches of many Scottish rivers means that failure of mill structures in Scotland will probably have less dramatic geomorphological and management implications than those suggested by Downward and Skinner for southern English rivers

A background-priority discrete boundary triangulation method

Discrete boundary triangulation methods generate triangular meshes through the centers of the boundary voxels of a volumetric object. At some voxel configurations it may be arbitrary whether a part of the volume should be included in the object or could be classified as background. Consequently, important details such as concave and convex edges and corners are not consistently preserved in the describing geometry. We present a "background priority" version of an existing "object priority" algorithm [6]. We show that the ad hoc configurations of the well-known Discretized Marching Cubes algorithm [13] can be derived from our method and that a combined triangulation with "object priority" and "background priority" better would preserve object details

Fractional Statistics in Three Dimensions: Compact Maxwell-Higgs System

We show that a (3+1)-dimensional system composed of an open magnetic vortex and an electrical point charge exhibits the phenomenon of Fermi-Bose transmutation. In order to provide the physical realization of this system we focus on the lattice compact scalar electrodynamics $SQED_c$ whose topological excitations are open Nielsen-Olesen strings with magnetic monopoles attached at their ends.Comment: 8 page

Phase Transitions Driven by Vortices in 2D Superfluids and Superconductors: From Kosterlitz-Thouless to 1st Order

The Landau-Ginzburg-Wilson hamiltonian is studied for different values of the parameter $\lambda$ which multiplies the quartic term (it turns out that this is equivalent to consider different values of the coherence length $\xi$ in units of the lattice spacing $a$). It is observed that amplitude fluctuations can change dramatically the nature of the phase transition: for small values of $\lambda$ ($\xi/a > 0.7$), instead of the smooth Kosterlitz-Thouless transition there is a {\em first order} transition with a discontinuous jump in the vortex density $v$ and a larger non-universal drop in the helicity modulus. In particular, for $\lambda$ sufficiently small ($\xi/a \cong 1$), the density of bound pairs of vortex-antivortex below $T_c$ is so low that, $v$ drops to zero almost for all temperature $T<Tc$.Comment: 8 pages, 5 .eps figure

Performance of the EUDET-type beam telescopes

Test beam measurements at the test beam facilities of DESY have been conducted to characterise the performance of the EUDET-type beam telescopes originally developed within the EUDET project. The beam telescopes are equipped with six sensor planes using MIMOSA26 monolithic active pixel devices. A programmable Trigger Logic Unit provides trigger logic and time stamp information on particle passage. Both data acquisition framework and offline reconstruction software packages are available. User devices are easily integrable into the data acquisition framework via predefined interfaces. The biased residual distribution is studied as a function of the beam energy, plane spacing and sensor threshold. Its standard deviation at the two centre pixel planes using all six planes for tracking in a 6\,GeV electron/positron-beam is measured to be (2.88\,\pm\,0.08)\,\upmu\meter.Iterative track fits using the formalism of General Broken Lines are performed to estimate the intrinsic resolution of the individual pixel planes. The mean intrinsic resolution over the six sensors used is found to be (3.24\,\pm\,0.09)\,\upmu\meter.With a 5\,GeV electron/positron beam, the track resolution halfway between the two inner pixel planes using an equidistant plane spacing of 20\,mm is estimated to (1.83\,\pm\,0.03)\,\upmu\meter assuming the measured intrinsic resolution. Towards lower beam energies the track resolution deteriorates due to increasing multiple scattering. Threshold studies show an optimal working point of the MIMOSA26 sensors at a sensor threshold of between five and six times their RMS noise. Measurements at different plane spacings are used to calibrate the amount of multiple scattering in the material traversed and allow for corrections to the predicted angular scattering for electron beams