MuMax: a new high-performance micromagnetic simulation tool

We present MuMax, a general-purpose micromagnetic simulation tool running on Graphical Processing Units (GPUs). MuMax is designed for high performance computations and specifically targets large simulations. In that case speedups of over a factor 100x can easily be obtained compared to the CPU-based OOMMF program developed at NIST. MuMax aims to be general and broadly applicable. It solves the classical Landau-Lifshitz equation taking into account the magnetostatic, exchange and anisotropy interactions, thermal effects and spin-transfer torque. Periodic boundary conditions can optionally be imposed. A spatial discretization using finite differences in 2 or 3 dimensions can be employed. MuMax is publicly available as open source software. It can thus be freely used and extended by community. Due to its high computational performance, MuMax should open up the possibility of running extensive simulations that would be nearly inaccessible with typical CPU-based simulators.Comment: To be published in JMM

Near Real-Time Data Labeling Using a Depth Sensor for EMG Based Prosthetic Arms

Recognizing sEMG (Surface Electromyography) signals belonging to a particular action (e.g., lateral arm raise) automatically is a challenging task as EMG signals themselves have a lot of variation even for the same action due to several factors. To overcome this issue, there should be a proper separation which indicates similar patterns repetitively for a particular action in raw signals. A repetitive pattern is not always matched because the same action can be carried out with different time duration. Thus, a depth sensor (Kinect) was used for pattern identification where three joint angles were recording continuously which is clearly separable for a particular action while recording sEMG signals. To Segment out a repetitive pattern in angle data, MDTW (Moving Dynamic Time Warping) approach is introduced. This technique is allowed to retrieve suspected motion of interest from raw signals. MDTW based on DTW algorithm, but it will be moving through the whole dataset in a pre-defined manner which is capable of picking up almost all the suspected segments inside a given dataset an optimal way. Elevated bicep curl and lateral arm raise movements are taken as motions of interest to show how the proposed technique can be employed to achieve auto identification and labelling. The full implementation is available at https://github.com/GPrathap/OpenBCIPytho

ROMA: a map-making algorithm for polarised CMB data sets

We present ROMA, a parallel code to produce joint optimal temperature and polarisation maps out of multidetector CMB observations. ROMA is a fast, accurate and robust implementation of the iterative generalised least squares approach to map-making. We benchmark ROMA on realistic simulated data from the last, polarisation sensitive, flight of BOOMERanG.Comment: Accepted for publication in Astronomy & Astrophysics. Version with higher quality figures available at http://www.fisica.uniroma2.it/~cosmo/ROM

Galaxy-galaxy(-galaxy) lensing as a sensitive probe of galaxy evolution

The gravitational lensing effect provides various ways to study the mass environment of galaxies. We investigate how galaxy-galaxy(-galaxy) lensing can be used to test models of galaxy formation and evolution. We consider two semi-analytic galaxy formation models based on the Millennium Run N-body simulation: the Durham model by Bower et al. (2006) and the Garching model by Guo et al. (2011). We generate mock lensing observations for the two models, and then employ Fast Fourier Transform methods to compute second- and third-order aperture statistics in the simulated fields for various galaxy samples. We find that both models predict qualitatively similar aperture signals, but there are large quantitative differences. The Durham model predicts larger amplitudes in general. In both models, red galaxies exhibit stronger aperture signals than blue galaxies. Using these aperture measurements and assuming a linear deterministic bias model, we measure relative bias ratios of red and blue galaxy samples. We find that a linear deterministic bias is insufficient to describe the relative clustering of model galaxies below ten arcmin angular scales. Dividing galaxies into luminosity bins, the aperture signals decrease with decreasing luminosity for brighter galaxies, but increase again for fainter galaxies. This increase is likely an artifact due to too many faint satellite galaxies in massive group and cluster halos predicted by the models. Our study shows that galaxy-galaxy(-galaxy) lensing is a sensitive probe of galaxy evolution.Comment: 11 pages, 8 figures, accepted in A&

A new experimental snow avalanche test site at Seehore peak in Aosta Valley (NW Italian Alps) - Part II: Engineering aspects

The estimate of the effects produced by the impact of a snow avalanche against an obstacle is of the utmost importance in designing safe mountain constructions. For this purpose, an ad-hoc instrumented obstacle was designed and built in order to measure impact forces of small and medium snow avalanches at Seehore peak (NW Italian Alps). The structural design had to consider several specific and unusual demands dictated by the difficult environment. In this article, the new test facility is described from the engineering point of view, discussing the most important aspects of the analyzed problems which were solved before and after the construction. The performance of the instrumented obstacle in the first two operating seasons, and some proposals for future upgrading are eventually illustrate

Impact of modulation on CMB B-mode polarization experiments

We investigate the impact of both slow and fast polarization modulation strategies on the science return of upcoming ground-based experiments aimed at measuring the B-mode polarization of the CMB. Using simulations of the Clover experiment, we compare the ability of modulated and un-modulated observations to recover the signature of gravitational waves in the polarized CMB sky in the presence of a number of anticipated systematic effects. The general expectations that fast modulation is helpful in mitigating low-frequency detector noise, and that the additional redundancy in the projection of the instrument's polarization sensitivity directions onto the sky when modulating reduces the impact of instrumental polarization, are borne out by our simulations. Neither low-frequency polarized atmospheric fluctuations nor systematic errors in the polarization sensitivity directions are mitigated by modulation. Additionally, we find no significant reduction in the effect of pointing errors by modulation. For a Clover-like experiment, pointing jitter should be negligible but any systematic mis-calibration of the polarization coordinate reference system results in significant E-B mixing on all angular scales and will require careful control. We also stress the importance of combining data from multiple detectors in order to remove the effects of common-mode systematics (such as 1/f atmospheric noise) on the measured polarization signal. Finally we compare the performance of our simulated experiment with the predicted performance from a Fisher analysis. We find good agreement between the Fisher predictions and the simulations except for the very largest scales where the power spectrum estimator we have used introduces additional variance to the B-mode signal recovered from our simulations.Comment: Replaced with version accepted by MNRAS. Analysis of half-wave plate systematic (differential transmittance) adde

Ray-tracing through the Millennium Simulation: Born corrections and lens-lens coupling in cosmic shear and galaxy-galaxy lensing

(abridged) We study the accuracy of various approximations to cosmic shear and weak galaxy-galaxy lensing and investigate effects of Born corrections and lens-lens coupling. We use ray-tracing through the Millennium Simulation to calculate various cosmic-shear and galaxy-galaxy-lensing statistics. We compare the results from ray-tracing to semi-analytic predictions. We find: (i) The linear approximation provides an excellent fit to cosmic-shear power spectra as long as the actual matter power spectrum is used as input. Common fitting formulae, however, strongly underestimate the cosmic-shear power spectra. Halo models provide a better fit to cosmic shear-power spectra, but there are still noticeable deviations. (ii) Cosmic-shear B-modes induced by Born corrections and lens-lens coupling are at least three orders of magnitude smaller than cosmic-shear E-modes. Semi-analytic extensions to the linear approximation predict the right order of magnitude for the B-mode. Compared to the ray-tracing results, however, the semi-analytic predictions may differ by a factor two on small scales and also show a different scale dependence. (iii) The linear approximation may under- or overestimate the galaxy-galaxy-lensing shear signal by several percent due to the neglect of magnification bias, which may lead to a correlation between the shear and the observed number density of lenses. We conclude: (i) Current semi-analytic models need to be improved in order to match the degree of statistical accuracy expected for future weak-lensing surveys. (ii) Shear B-modes induced by corrections to the linear approximation are not important for future cosmic-shear surveys. (iii) Magnification bias can be important for galaxy-galaxy-lensing surveys.Comment: version taking comments into accoun

Superpulsed low-level laser therapy protects skeletal muscle of mdx mice against damage, inflammation and morphological changes delaying dystrophy progression.

Aim: To evaluate the effects of preventive treatment with low-level laser therapy (LLLT) on progression of dystrophy in mdx mice. Methods: Ten animals were randomly divided into 2 experimental groups treated with superpulsed LLLT (904 nm, 15 mW, 700 Hz, 1 J) or placebo-LLLT at one point overlying the tibialis anterior muscle (bilaterally) 5 times per week for 14 weeks (from 6th to 20th week of age). Morphological changes, creatine kinase (CK) activity and mRNA gene expression were assessed in animals at 20th week of age. Results: Animals treated with LLLT showed very few morphological changes in skeletal muscle, with less atrophy and fibrosis than animals treated with placebo-LLLT. CK was significantly lower (p = 0.0203) in animals treated with LLLT (864.70 U.l−1, SEM 226.10) than placebo (1708.00 U.l−1, SEM 184.60). mRNA gene expression of inflammatory markers was significantly decreased by treatment with LLLT (p<0.05): TNF-α (placebo-control = 0.51 µg/µl [SEM 0.12], - LLLT = 0.048 µg/µl [SEM 0.01]), IL-1β (placebo-control = 2.292 µg/µl [SEM 0.74], - LLLT = 0.12 µg/µl [SEM 0.03]), IL-6 (placebo-control = 3.946 µg/µl [SEM 0.98], - LLLT = 0.854 µg/µl [SEM 0.33]), IL-10 (placebo-control = 1.116 µg/µl [SEM 0.22], - LLLT = 0.352 µg/µl [SEM 0.15]), and COX-2 (placebo-control = 4.984 µg/µl [SEM 1.18], LLLT = 1.470 µg/µl [SEM 0.73]). Conclusion: Irradiation of superpulsed LLLT on successive days five times per week for 14 weeks decreased morphological changes, skeletal muscle damage and inflammation in mdx mice. This indicates that LLLT has potential to decrease progression of Duchenne muscular dystrophy

Kinematic modelling of disk galaxies I. A new method to fit tilted rings to data cubes

This is the first of a series of papers in which the kinematics of disk galaxies over a range of scales is scrutinised employing spectroscopy. A fundamental aspect of these studies is presented here: the new publicly available software tool TiRiFiC (http://www.astro.uni-bonn.de/~gjozsa/tirific.html) enables a direct fit of a ``tilted-ring model'' to spectroscopic data cubes. The algorithm generates model data cubes from the tilted-ring parametrisation of a rotating disk, which are automatically adjusted to reach an optimum fit via a chi-squared minimisation method to an observed data cube. The structure of the new software, the shortcomings of the previously available programs to produce a tilted-ring model, and the performance of TiRiFiC are discussed. Our method is less affected by the well-known problem of beam smearing that occurs when fitting to the velocity field. Since with our method we fit many data points in a data cube simultaneously, TiRiFiC is sensitive to very faint structures and can hence be used to derive tilted-ring models significantly extending in radius beyond those derived from a velocity field. The software is able to parametrise HI disks of galaxies that are intersected by the line-of-sight twice or more, i.e. if the disks are heavily warped, and/or with a significant shift of the projected centre of rotation, and/or if seen edge-on. Furthermore, our method delivers the surface-brightness profile of the examined galaxy in addition to the orientational parameters and the rotation curve. In order to derive kinematic and morphological models of disk galaxies, especially reliable rotation curves, a direct-fit method as implemented in our code should be the tool of choice.Comment: 37 pages, 24 figures. Accepted for publication in Astronomy and Astrophysic

O(N) continuous electrostatics solvation energies calculation method for biomolecules simulations

We report a development of a new fast surface-based method for numerical calculations of solvation energy of biomolecules with a large number of charged groups. The procedure scales linearly with the system size both in time and memory requirements, is only a few percent wrong for any molecular configurations of arbitrary sizes, gives explicit value for the reaction field potential at any point, provides both the solvation energy and its derivatives suitable for Molecular Dynamics simulations. The method works well both for large and small molecules and thus gives stable energy differences for quantities such as solvation energies of molecular complex formation.Comment: 6 pages, 4 figures, more results, examples and references adde