Monitoring data in R with the lumberjack package

Monitoring data while it is processed and transformed can yield detailed insight into the dynamics of a (running) production system. The lumberjack package is a lightweight package allowing users to follow how an R object is transformed as it is manipulated by R code. The package abstracts all logging code from the user, who only needs to specify which objects are logged and what information should be logged. A few default loggers are included with the package but the package is extensible through user-defined logger objects.Comment: Accepted for publication in the Journal of Statistical Softwar

Studies of a Terawatt X-Ray Free-Electron Laser

The possibility of constructing terawatt (TW) x-ray free-electron lasers (FELs) has been discussed using novel superconducting helical undulators [5]. In this paper, we consider the conditions necessary for achieving powers in excess of 1 TW in a 1.5 {\AA} FEL using simulations with the MINERVA simulation code [7]. Steady-state simulations have been conducted using a variety of undulator and focusing configurations. In particular, strong focusing using FODO lattices is compared with the natural, weak focusing inherent in helical undulators. It is found that the most important requirement to reach TW powers is extreme transverse compression of the electron beam in a strong FODO lattice. The importance of extreme focusing of the electron beam in the production of TW power levels means that the undulator is not the prime driver for a TW FEL, and simulations are also described using planar undulators that reach near-TW power levels. In addition, TW power levels can be reached using pure self-amplified spontaneous emission (SASE) or with novel self-seeding configurations when such extreme focusing of the electron beam is applied.Comment: 10 pages, 12 figure

An N-body/SPH Study of Isolated Galaxy Mass Density Profiles

We investigate the evolution of mass density profiles in secular disk galaxy models, paying special attention to the development of a two-component profile from a single initial exponential disk free of cosmological evolution (i.e., no accretion or interactions). As the source of density profile variations, we examine the parameter space of the spin parameter, halo concentration, virial mass, disk mass and bulge mass, for a total of 162 simulations in the context of a plausible model of star formation and feedback (GADGET-2). The evolution of the galaxy mass density profile, including the development of a two-component profile with an inner and outer segment, is controlled by the ratio of the disk mass fraction, $m_{d}$, to the halo spin parameter, $\lambda$. The location of the break between the two components and speed at which it develops is directly proportional to $m_{d}/\lambda$; the amplitude of the transition between the inner and outer regions is however controlled by the ratio of halo concentration to virial velocity. The location of the divide between the inner and outer profile does not change with time. (Abridged)Comment: 27 pages, 31 figures. Accepted for publication at MNRAS. A high-resolution version of the paper with figures can be found here http://www.mpia-hd.mpg.de/~foyle/papers/MN-07-1491-MJ.R1.pd

Spatially Resolved Stellar Kinematics of Field Early-Type Galaxies at z=1: Evolution of the Rotation Rate

We use the spatial information of our previously published VLT/FORS2 absorption line spectroscopy to measure mean stellar velocity and velocity dispersion profiles of 25 field early-type galaxies at a median redshift z=0.97 (full range 0.6<z<1.2). This provides the first detailed study of early-type galaxy rotation at these redshifts. From surface brightness profiles from HST imaging we calculate two-integral oblate axisymmetric Jeans equation models for the observed kinematics. Fits to the data yield for each galaxy the degree of rotational support and the mass-to-light ratio M/L_Jeans. S0 and Sa galaxies are generally rotationally supported, whereas elliptical galaxies rotate less rapidly or not at all. Down to M(B)=-19.5 (corrected for luminosity evolution), we find no evidence for evolution in the fraction of rotating early-type (E+S0) galaxies between z=1 (63+/-11%) and the present (61+/-5%). We interpret this as evidence for little or no change in the field S0 fraction with redshift. We compare M/L_Jeans with M/L_vir inferred from the virial theorem and globally averaged quantities and assuming homologous evolution. There is good agreement for non-rotating (mostly E) galaxies. However, for rotationally supported galaxies (mostly S0) M/L_Jeans is on average ~40% higher than M/L_vir. We discuss possible explanations and the implications for the evolution of M/L between z=1 and the present and its dependence on mass.Comment: To appear in ApJ 683 (9 pages, 7 figures). Minor changes included to match published versio

Dynamical Monte Carlo Study of Equilibrium Polymers (II): The Role of Rings

We investigate by means of a number of different dynamical Monte Carlo simulation methods the self-assembly of equilibrium polymers in dilute, semidilute and concentrated solutions under good-solvent conditions. In our simulations, both linear chains and closed loops compete for the monomers, expanding on earlier work in which loop formation was disallowed. Our findings show that the conformational properties of the linear chains, as well as the shape of their size distribution function, are not altered by the formation of rings. Rings only seem to deplete material from the solution available to the linear chains. In agreement with scaling theory, the rings obey an algebraic size distribution, whereas the linear chains conform to a Schultz--Zimm type of distribution in dilute solution, and to an exponentional distribution in semidilute and concentrated solution. A diagram presenting different states of aggregation, including monomer-, ring- and chain-dominated regimes, is given

The excitation of near-infrared H2 emission in NGC 253

Because of its large angular size and proximity to the Milky Way, NGC 253, an archetypal starburst galaxy, provides an excellent laboratory to study the intricacies of this intense episode of star formation. We aim to characterize the excitation mechanisms driving the emission in NGC 253. Specifically we aim to distinguish between shock excitation and UV excitation as the dominant driving mechanism, using Br\gamma, H_2 and [FeII] as diagnostic emission line tracers. Using SINFONI observations, we create linemaps of Br\gamma, [FeII]_{1.64}, and all detected H_2 transitions. By using symmetry arguments of the gas and stellar gas velocity field, we find a kinematic center in agreement with previous determinations. The ratio of the 2-1 S(1) to 1-0 S(1) H_2 transitions can be used as a diagnostic to discriminate between shock and fluorescent excitation. Using the 1-0 S(1)/2-1 S(1) line ratio as well as several other H_2 line ratios and the morphological comparison between H_2 and Br\gamma and [FeII], we find that excitation from UV photons is the dominant excitation mechanisms throughout NGC 253. We employ a diagnostic energy level diagram to quantitatively differentiate between mechanisms. We compare the observed energy level diagrams to PDR and shock models and find that in most regions and over the galaxy as a whole, fluorescent excitation is the dominant mechanism exciting the H_2 gas. We also place an upper limit of the percentage of shock excited H_2 at 29%. We find that UV radiation is the dominant excitation mechanism for the H_2 emission. The H_2 emission does not correlate well with Br\gamma but closely traces the PAH emission, showing that not only is H_2 fluorescently excited, but it is predominately excited by slightly lower mass stars than O stars which excite Br\gamma, such as B stars

Potential environmental impact of tidal energy extraction in the Pentland Firth at large spatial scales : results of a biogeochemical model

A model study was carried out of the potential large-scale (> 100 km) effects of marine renewable tidal energy generation in the Pentland Firth, using the 3-D hydrodynamics–biogeochemistry model GETM-ERSEM-BFM. A realistic 800 MW scenario and a high-impact scenario with massive expansion of tidal energy extraction to 8 GW scenario were considered. The realistic 800 MW scenario suggested minor effects on the tides, and undetectable effects on the biogeochemistry. The massive-expansion 8 GW scenario suggested effects would be observed over hundreds of kilometres away with changes of up to 10 % in tidal and ecosystem variables, in particular in a broad area in the vicinity of the Wash. There, waters became less turbid, and primary production increased with associated increases in faunal ecosystem variables. Moreover, a one-off increase in carbon storage in the sea bed was detected. Although these first results suggest positive environmental effects, further investigation is recommended of (i) the residual circulation in the vicinity of the Pentland Firth and effects on larval dispersal using a higher-resolution model and (ii) ecosystem effects with (future) state-of-the-art models if energy extraction substantially beyond 1 GW is planned