A portable platform for accelerated PIC codes and its application to GPUs using OpenACC

We present a portable platform, called PIC_ENGINE, for accelerating Particle-In-Cell (PIC) codes on heterogeneous many-core architectures such as Graphic Processing Units (GPUs). The aim of this development is efficient simulations on future exascale systems by allowing different parallelization strategies depending on the application problem and the specific architecture. To this end, this platform contains the basic steps of the PIC algorithm and has been designed as a test bed for different algorithmic options and data structures. Among the architectures that this engine can explore, particular attention is given here to systems equipped with GPUs. The study demonstrates that our portable PIC implementation based on the OpenACC programming model can achieve performance closely matching theoretical predictions. Using the Cray XC30 system, Piz Daint, at the Swiss National Supercomputing Centre (CSCS), we show that PIC_ENGINE running on an NVIDIA Kepler K20X GPU can outperform the one on an Intel Sandybridge 8-core CPU by a factor of 3.4

ALMA observations of the variable 12CO/13CO ratio around the asymptotic giant branch star R Sculptoris

[abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....Comment: 6 pages, 5 figures, online data available at http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=J/A+A/556/L

Higher-Order Angular Galaxy Correlations in the SDSS: Redshift and Color Dependence of non-Linear Bias

We present estimates of the N-point galaxy, area-averaged, angular correlation functions $\bar{\omega}_{N}$($\theta$) for $N$ = 2,...,7 for galaxies from the fifth data release of the Sloan Digital Sky Survey. Our parent sample is selected from galaxies with $18 \leq r < 21$, and is the largest ever used to study higher-order correlations. We subdivide this parent sample into two volume limited samples using photometric redshifts, and these two samples are further subdivided by magnitude, redshift, and color (producing early- and late-type galaxy samples) to determine the dependence of $\bar{\omega}_{N}$($\theta$) on luminosity, redshift, and galaxy-type. We measure $\bar{\omega}_{N}$($\theta$) using oversampling techniques and use them to calculate the projected, $s_{N}$. Using models derived from theoretical power-spectra and perturbation theory, we measure the bias parameters $b_1$ and $c_2$, finding that the large differences in both bias parameters ($b_1$ and $c_2$) between early- and late-type galaxies are robust against changes in redshift, luminosity, and $\sigma_8$, and that both terms are consistently smaller for late-type galaxies. By directly comparing their higher-order correlation measurements, we find large differences in the clustering of late-type galaxies at redshifts lower than 0.3 and those at redshifts higher than 0.3, both at large scales ($c_2$ is larger by $\sim0.5$ at $z > 0.3$) and small scales (large amplitudes are measured at small scales only for $z > 0.3$, suggesting much more merger driven star formation at $z > 0.3$). Finally, our measurements of $c_2$ suggest both that $\sigma_8 < 0.8$ and $c_2$ is negative.Comment: 46 pages, 19 figures, Accepted to Ap

Developing a Pilot Case and Modelling the Development of a Large European CO₂ Transport Infrastructure -The GATEWAY H2020 Project

The H2020 GATEWAY project aims to develop a comprehensive model Pilot Case which, intentionally, will pave the ground for CCS deployment in Europe. It will result from the assessment of, technical, commercial, judicial and societal issues related to a future CO2 transport infrastructure. The Pilot Case derived on this basis, will emphasize a gateway for CO2 transport in the North Sea Basin. Four potential pilot cases have been evaluated through a combination of techno-economic modelling of the individual cases and evaluation against more qualitative criteria. The chosen Pilot Case, Rotterdam Nucleus, will be refined and developed during the remaining period of the GATEWAY project. To maximise impact, the GATEWAY project adapts its work to lay the foundation for a future application to a European ‘Project of Common Interest’ (PCI). Continuous dialogue with the most relevant stakeholders is an important part of GATEWAY, as a Coordination and Support Action (CSA) H2020 project

Biases affecting injected doses of an experimental drug during clinical trials.

During clinical trials, researchers rarely question nominal doses specified on labels of investigational products, overlooking the potential for inaccuracies that may result when calculating pharmacokinetic and pharmacodynamic parameters. This study evaluated the disparity between nominal doses and the doses actually administered in two Phase I trials of a biosimilar drug. In Trial A, 12 healthy volunteers received various doses of an interferon β-1a biosimilar via either subcutaneous or intravenous injection, prepared by partially emptying 0.53 ml syringes supplied by the manufacturer. In Trial B, 12 volunteers received three different formulations of the drug via intravenous injection (biosimilar with and without albumin and a comparator), followed by multiple subcutaneous injections. In both trials, the dose administered was calculated as D = C × V - losses, where C is the drug concentration assessed using ELISA, V is the volume administered calculated using syringe weighing and losses are deduced from in-vitro experiments. Interferon binding to added albumin and infusion lines was evaluated using a (125)I-interferon tracer with gel-filtration chromatography. In Trial A, measured concentrations were close to the nominal strength indicated by the manufacturer (median bias: -6 %), whereas in Trial B they differed significantly for all three formulations (median biases: +67 %, +73 % and +31 % for the biosimilar with albumin, the biosimilar without albumin and the comparator, respectively). In Trial A, the doses actually administered showed large variability and biases, especially at the lowest doses. Indeed, actually injected volumes differed by as much as 74 % from theoretical volumes - a phenomenon mainly attributed to unnoticed fluid re-aspiration through the syringe needle. This was corrected in Trial B. Interferon was not significantly adsorbed on the infusion lines used for intravenous administration. Its binding to albumin was slow, reaching 50 % after a 16 h incubation. These examples illustrate the importance of assessing the actual doses administered in clinical trials, to ensure accuracy in the determination of clearance, distribution volume, bioavailability and dose-response relationships. Clinicaltrials.gov NCT02515695 (Trial A) and NCT02517788 (Trial B). Registered on 24 July and 5 August 2015, respectively

First Penning-trap mass measurement in the millisecond half-life range: the exotic halo nucleus 11Li

In this letter, we report a new mass for $^{11}$Li using the trapping experiment TITAN at TRIUMF's ISAC facility. This is by far the shortest-lived nuclide, $t_{1/2} = 8.8 \rm{ms}$, for which a mass measurement has ever been performed with a Penning trap. Combined with our mass measurements of $^{8,9}$Li we derive a new two-neutron separation energy of 369.15(65) keV: a factor of seven more precise than the best previous value. This new value is a critical ingredient for the determination of the halo charge radius from isotope-shift measurements. We also report results from state-of-the-art atomic-physics calculations using the new mass and extract a new charge radius for $^{11}$Li. This result is a remarkable confluence of nuclear and atomic physics.Comment: Formatted for submission to PR

The circumstellar envelope around the S-type AGB star W Aql Effects of an eccentric binary orbit

The CO(J=3-2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE. The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed. The estimated average mass-loss rate of W~Aql agrees with previous results. The size of the emitting region is consistent with photodissociation models. The CO(J=3-2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of 10" and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yrs) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate.Comment: 10 pages, 8 figure

Clarifications to the limitations of the s-α equilibrium model for gyrokinetic computations of turbulence

In the context of gyrokinetic flux-tube simulations of microturbulence in magnetized toroidal plasmas, different treatments of the magnetic equilibrium are examined. Considering the Cyclone DIII-D base case parameter set [Dimits et al., Phys. Plasmas 7, 969 (2000)], significant differences in the linear growth rates, the linear and nonlinear critical temperature gradients, and the nonlinear ion heat diffusivities are observed between results obtained using either an $s-\alpha$ or an MHD equilibrium. Similar disagreements have been reported previously [Redd et al., Phys. Plasmas 6, 1162 (1999)]. In this paper it is shown that these differences result primarily from the approximation made in the standard implementation of the $s-\alpha$ model, in which the straight field line angle is identified to the poloidal angle, leading to inconsistencies of order $\varepsilon$ ($\varepsilon=a/R$ is the inverse aspect ratio, $a$ the minor radius and $R$ the major radius). An equilibrium model with concentric, circular flux surfaces and a correct treatment of the straight field line angle gives results very close to those using a finite $\varepsilon$, low $\beta$ MHD equilibrium. Such detailed investigation of the equilibrium implementation is of particular interest when comparing flux tube and global codes. It is indeed shown here that previously reported agreements between local and global simulations in fact result from the order $\varepsilon$ inconsistencies in the $s-\alpha$ model, coincidentally compensating finite $\rho^*$ effects in the global calculations, where $\rho^*=\rho _s / a$ with $\rho_s$ the ion sound Larmor radius. True convergence between local and global simulations is finally obtained by correct treatment of the geometry in both cases, and considering the appropriate $\rho^* \rightarrow 0$ limit in the latter case