N-body simulations with two-orders-of-magnitude higher performance using wavelets

Noise is a problem of major concern for N-body simulations of structure formation in the early Universe, of galaxies and plasmas. Here for the first time we use wavelets to remove noise from N-body simulations of disc galaxies, and show that they become equivalent to simulations with two orders of magnitude more particles. We expect a comparable improvement in performance for cosmological and plasma simulations. Our wavelet code will be described in a following paper, and will then be available on request.Comment: Mon. Not. R. Astron. Soc., in press. The interested reader is strongly recommended to ignore the low-resolution Fig. 3 (and Fig. 4), and to download the full-resolution paper (700 kb) from http://www.oso.chalmers.se/~romeo/Paper_VI.ps.g

Pressure-energy correlations and thermodynamic scaling in viscous Lennard-Jones liquids

We use molecular dynamics simulation results on viscous binary Lennard-Jones mixtures to examine the correlation between the potential energy and the virial. In accord with a recent proposal [U. R. Pedersen et. al. Phys. Rev. Lett. 100, 015701 (2008)], the fluctuations in the two quantities are found to be strongly correlated, exhibiting a proportionality constant, Gamma, numerically equal to one-third the slope of an inverse power law approximation to the intermolecular potential function. The correlation is stronger at higher densities, where interatomic separations are in the range where the inverse power law approximation is more accurate. These same liquids conform to thermodynamic scaling of their dynamics, with the scaling exponent equal to Gamma. Thus, the properties of strong correlation between energy and pressure and thermodynamic scaling both reflect the ability of an inverse power law representation of the potential to capture interesting features of the dynamics of dense, highly viscous liquids.Comment: 5 pages, 4 figures; published version, one figure remove

Use of Iverson's language APL for statistical computing

APL computer language for statistical computatio

A wavelet add-on code for new-generation N-body simulations and data de-noising (JOFILUREN)

Wavelets are a new and powerful mathematical tool, whose most celebrated applications are data compression and de-noising. In Paper I (Romeo, Horellou & Bergh 2003, astro-ph/0302343), we have shown that wavelets can be used for removing noise efficiently from cosmological, galaxy and plasma N-body simulations. The expected two-orders-of-magnitude higher performance means, in terms of the well-known Moore's law, an advance of more than one decade in the future. In this paper, we describe a wavelet add-on code designed for such an application. Our code can be included in common grid-based N-body codes, is written in Fortran, is portable and available on request from the first author. The code can also be applied for removing noise from standard data, such as signals and images.Comment: Mon. Not. R. Astron. Soc., in press. The interested reader is strongly recommended to ignore the low-resolution Figs 10 and 11, and to download the full-resolution paper (800 kb) from http://www.oso.chalmers.se/~romeo/Paper_VII.ps.g

Events, processes, and the time of a killing

The paper proposes a novel solution to the problem of the time of a killing (ToK), which persistently besets theories of act-individuation. The solution proposed claims to expose a crucial wrong-headed assumption in the debate, according to which ToK is essentially a problem of locating some event that corresponds to the killing. The alternative proposal put forward here turns on recognizing a separate category of dynamic occurents, viz. processes. The paper does not aim to mount a comprehensive defense of process ontology, relying instead on extant defenses. The primary aim is rather to put process ontology to work in diagnosing the current state of play over ToK, and indeed in solving it

Seriously misleading results using inverse of Freeman-Tukey double arcsine transformation in meta-analysis of single proportions.

Standard generic inverse variance methods for the combination of single proportions are based on transformed proportions using the logit, arcsine, and Freeman-Tukey double arcsine transformations. Generalized linear mixed models are another more elaborate approach. Irrespective of the approach, meta-analysis results are typically back-transformed to the original scale in order to ease interpretation. Whereas the back-transformation of meta-analysis results is straightforward for most transformations, this is not the case for the Freeman-Tukey double arcsine transformation, albeit possible. In this case study with five studies, we demonstrate how seriously misleading the back-transformation of the Freeman-Tukey double arcsine transformation can be. We conclude that this transformation should only be used with special caution for the meta-analysis of single proportions due to potential problems with the back-transformation. Generalized linear mixed models seem to be a promising alternative

Source detection using a 3D sparse representation: application to the Fermi gamma-ray space telescope

The multiscale variance stabilization Transform (MSVST) has recently been proposed for Poisson data denoising. This procedure, which is nonparametric, is based on thresholding wavelet coefficients. We present in this paper an extension of the MSVST to 3D data (in fact 2D-1D data) when the third dimension is not a spatial dimension, but the wavelength, the energy, or the time. We show that the MSVST can be used for detecting and characterizing astrophysical sources of high-energy gamma rays, using realistic simulated observations with the Large Area Telescope (LAT). The LAT was launched in June 2008 on the Fermi Gamma-ray Space Telescope mission. The MSVST algorithm is very fast relative to traditional likelihood model fitting, and permits efficient detection across the time dimension and immediate estimation of spectral properties. Astrophysical sources of gamma rays, especially active galaxies, are typically quite variable, and our current work may lead to a reliable method to quickly characterize the flaring properties of newly-detected sources.Comment: Accepted. Full paper will figures available at http://jstarck.free.fr/aa08_msvst.pd

Economics, Agency, and Causal Explanation

The paper considers three questions. First, what is the connection between economics and agency? It is argued that causation and explanation in economics fundamentally depend on agency. So a philosophical understanding of economic explanation must be sensitive to an understanding of agency. Second, what is the connection between agency and causation? A causal view of agency-involving explanation is defended against a number of arguments from the resurgent noncausalist tradition in the literature on agency and action-explanation. If agency is fundamental to economic explanation, it is argued, then so is causation. Third, what is the connection between causal explanation and the natural sciences? It is argued that, though the explanations given in economics and other social sciences are causal explanations, they are different in kind from the causal explanations of the natural sciences. On the one hand, then, the causal explanations of the social sciences are irreducible to those found in the natural sciences. On the other hand, the causal relations described by the social sciences are not completely autonomous; they do not float free of, or operate independently from, the causal relations charted by the natural sciences