Scanamorphos: a map-making software for Herschel and similar scanning bolometer arrays

Scanamorphos is one of the public softwares available to post-process scan observations performed with the Herschel photometer arrays. This post-processing mainly consists in subtracting the total low-frequency noise (both its thermal and non-thermal components), masking high-frequency artefacts such as cosmic ray hits, and projecting the data onto a map. Although it was developed for Herschel, it is also applicable with minimal adjustment to scan observations made with some other imaging arrays subjected to low-frequency noise, provided they entail sufficient redundancy; it was successfully applied to P-Artemis, an instrument operating on the APEX telescope. Contrary to matrix-inversion softwares and high-pass filters, Scanamorphos does not assume any particular noise model, and does not apply any Fourier-space filtering to the data, but is an empirical tool using purely the redundancy built in the observations -- taking advantage of the fact that each portion of the sky is sampled at multiple times by multiple bolometers. It is an interactive software in the sense that the user is allowed to optionally visualize and control results at each intermediate step, but the processing is fully automated. This paper describes the principles and algorithm of Scanamorphos and presents several examples of application.Comment: This is the final version as accepted by PASP (on July 27, 2013). A copy with much better-quality figures is available on http://www2.iap.fr/users/roussel/herschel

Geometrodynamics in a spherically symmetric, static crossflow of null dust

The spherically symmetric, static spacetime generated by a crossflow of non-interacting radiation streams, treated in the geometrical optics limit (null dust) is equivalent to an anisotropic fluid forming a radiation atmosphere of a star. This reference fluid provides a preferred / internal time, which is employed as a canonical coordinate. Among the advantages we encounter a new Hamiltonian constraint, which becomes linear in the momentum conjugate to the internal time (therefore yielding a functional Schr\"{o}dinger equation after quantization), and a strongly commuting algebra of the new constraints.Comment: Section on boundary behavior and fall-off conditions of canonical variables added. New references, 1 new figure, 12 pages. Version accepted in Phys.Rev.

Dome-building volcanic activity in the Oas-Gutai Neogene Volcanic Area, Eastern Carpathians, Romania

A complex dome-building volcanic activity developed during a 5 Myr time interval (13.2- 8.0 Ma) in Oaş-Gutâi Mts., associated to the intermediate volcanism of the Oaş-Gutâi Neogene volcanic area (OG). Numerous domes were built up in the entire volcanic region also triggering both non-explosive and explosive fragmentation volcanic processes. The volcanic forms consist of extrusive domes, lava domes and dome- flows/coulées and cryptodomes predominantly as solitary domes, or compound domes and dome complexes. The domes are comprised of andesites, dacites and rhyolites (acid andesites and dacites are prevalent). The volcanic rocks show a calc- alkaline and medium to high-K character and typical subduction-zone geochemical signatures. Overall, either subaerial or subaqueous, the dome growth and collapse associated with fragmental explosive or non-explosive processes, was dominantly responsible for most of the volcanic products. Dome emplacement in submarine setting is commonly associated with marginal auto- brecciation, much subordinated explosive events and subsequent resedimentation. Overall, the dome-building volcanic activity in OG is recorded to a monogenetic-type of volcanism. The series of dome-building events which were triggered and controlled by magma-mixing and -mingling processes developed from time to time in different locations of O

Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise II. Fully discrete schemes

We present an abstract framework for analyzing the weak error of fully discrete approximation schemes for linear evolution equations driven by additive Gaussian noise. First, an abstract representation formula is derived for sufficiently smooth test functions. The formula is then applied to the wave equation, where the spatial approximation is done via the standard continuous finite element method and the time discretization via an I-stable rational approximation to the exponential function. It is found that the rate of weak convergence is twice that of strong convergence. Furthermore, in contrast to the parabolic case, higher order schemes in time, such as the Crank-Nicolson scheme, are worthwhile to use if the solution is not very regular. Finally we apply the theory to parabolic equations and detail a weak error estimate for the linearized Cahn-Hilliard-Cook equation as well as comment on the stochastic heat equation

Far-Infrared Properties of Spitzer-selected Luminous Starbursts

We present SHARC-2 350 micron data on 20 luminous z ~ 2 starbursts with S(1.2mm) > 2 mJy from the Spitzer-selected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S(350um) > 25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass--temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6um--2mm measurements of local starbursts. We find characteristic single-component temperatures T1 ~ 35.5+-2.2 K and integrated infrared (IR) luminosities around 10^(12.9+-0.1) Lsun for the SWIRE-selected sources. Molecular gas masses are estimated at 4 x 10^(10) Msun, assuming kappa(850um)=0.15 m^2/kg and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply >~2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 um flux densities, provides an effective means for the study of SMGs at z ~ 1.5--2.5.Comment: 13 pages, 4 figures, edited to match published version in ApJ 717, 29-39 (2010

Vortex waistlines and long range fluctuations

We examine the manner in which a linear potential results from fluctuations due to vortices linked with the Wilson loop. Our discussion is based on exact relations and inequalities between the Wilson loop and the vortex and electric flux order parameters. We show that, contrary to the customary naive picture, only vortex fluctuations of thickness of the order of the spatial linear size of the loop are capable of producing a strictly linear potential. An effective theory of these long range fluctuations emerges naturally in the form of a strongly coupled Z(N) lattice gauge theory. We also point out that dynamical fermions introduced in this medium undergo chiral symmetry breaking.Comment: 17 pages, LaTex file with 7 eps figures, revised references, minor comments adde

Center Vortices, Nexuses, and Fractional Topological Charge

It has been remarked in several previous works that the combination of center vortices and nexuses (a nexus is a monopole-like soliton whose world line mediates certain allowed changes of field strengths on vortex surfaces) carry topological charge quantized in units of 1/N for gauge group SU(N). These fractional charges arise from the interpretation of the standard topological charge integral as a sum of (integral) intersection numbers weighted by certain (fractional) traces. We show that without nexuses the sum of intersection numbers gives vanishing topological charge (since vortex surfaces are closed and compact). With nexuses living as world lines on vortices, the contributions to the total intersection number are weighted by different trace factors, and yield a picture of the total topological charge as a linking of a closed nexus world line with a vortex surface; this linking gives rise to a non-vanishing but integral topological charge. This reflects the standard 2\pi periodicity of the theta angle. We argue that the Witten-Veneziano relation, naively violating 2\pi periodicity, scales properly with N at large N without requiring 2\pi N periodicity. This reflects the underlying composition of localized fractional topological charge, which are in general widely separated. Some simple models are given of this behavior. Nexuses lead to non-standard vortex surfaces for all SU(N) and to surfaces which are not manifolds for N>2. We generalize previously-introduced nexuses to all SU(N) in terms of a set of fundamental nexuses, which can be distorted into a configuration resembling the 't Hooft-Polyakov monopole with no strings. The existence of localized but widely-separated fractional topological charges, adding to integers only on long distance scales, has implications for chiral symmetry breakdown.Comment: 15 pages, revtex, 6 .eps figure