Detecting multivariate interactions in spatial point patterns with Gibbs models and variable selection

We propose a method for detecting significant interactions in very large multivariate spatial point patterns. This methodology develops high dimensional data understanding in the point process setting. The method is based on modelling the patterns using a flexible Gibbs point process model to directly characterise point-to-point interactions at different spatial scales. By using the Gibbs framework significant interactions can also be captured at small scales. Subsequently, the Gibbs point process is fitted using a pseudo-likelihood approximation, and we select significant interactions automatically using the group lasso penalty with this likelihood approximation. Thus we estimate the multivariate interactions stably even in this setting. We demonstrate the feasibility of the method with a simulation study and show its power by applying it to a large and complex rainforest plant population data set of 83 species

A New Limit on the Antiproton Lifetime

Measurements of the cosmic ray pbar/p ratio are compared to predictions from an inhomogeneous disk-diffusion model of pbar production and propagation within the Galaxy, combined with a calculation of the modulation of the interstellar cosmic ray spectra as the particles propagate through the heliosphere to the Earth. The predictions agree with the observed pbar/p spectrum. Adding a finite pbar lifetime to the model, we obtain the limit tau_pbar > 0.8 Myr (90 % C.L.).Comment: 13 pages, 3 encapsulated Postscript figures, uses AASTeX; accepted by Astrophysical Journal; minor change

Large-scale groundwater modeling using global datasets: a test case for the Rhine-Meuse basin

The current generation of large-scale hydrological models does not include a groundwater flow component. Large-scale groundwater models, involving aquifers and basins of multiple countries, are still rare mainly due to a lack of hydro-geological data which are usually only available in developed countries. In this study, we propose a novel approach to construct large-scale groundwater models by using global datasets that are readily available. As the test-bed, we use the combined Rhine-Meuse basin that contains groundwater head data used to verify the model output. We start by building a distributed land surface model (30 arc-second resolution) to estimate groundwater recharge and river discharge. Subsequently, a MODFLOW transient groundwater model is built and forced by the recharge and surface water levels calculated by the land surface model. Results are promising despite the fact that we still use an offline procedure to couple the land surface and MODFLOW groundwater models (i.e. the simulations of both models are separately performed). The simulated river discharges compare well to the observations. Moreover, based on our sensitivity analysis, in which we run several groundwater model scenarios with various hydro-geological parameter settings, we observe that the model can reasonably well reproduce the observed groundwater head time series. However, we note that there are still some limitations in the current approach, specifically because the offline-coupling technique simplifies the dynamic feedbacks between surface water levels and groundwater heads, and between soil moisture states and groundwater heads. Also the current sensitivity analysis ignores the uncertainty of the land surface model output. Despite these limitations, we argue that the results of the current model show a promise for large-scale groundwater modeling practices, including for data-poor environments and at the global scale

Long-Baseline Study of the Leading Neutrino Oscillation at a Neutrino Factory

Within the framework of three-flavor neutrino oscillations, we consider the physics potential of \nu_e --> \nu_\mu appearance and \nu_\mu --> \nu_\mu survival measurements at a neutrino factory for a leading oscillation scale \delta m^2 ~ 3.5 \times 10^{-3} eV^2. Event rates are evaluated versus baseline and stored muon energy, and optimal values discussed. Over a sizeable region of oscillation parameter space, matter effects would enable the sign of \delta m^2 to be determined from a comparison of \nu_e --> \nu_\mu with \bar\nu_e --> \bar\nu_\mu event rates and energy distributions. It is important, therefore, that both positive and negative muons can be stored in the ring. Measurements of the \nu_\mu --> \nu_\mu survival spectrum could determine the magnitude of \delta m^2 and the leading oscillation amplitude with a precision of O(1%--2%).Comment: 33 pages, single-spaced Revtex, uses epsf.sty, 14 postscript figures. Added references, expanded conclusions, improved figs. 13 and 14. Version to be published in Phys. Rev.

Moment inversion problem for piecewise D-finite functions

We consider the problem of exact reconstruction of univariate functions with jump discontinuities at unknown positions from their moments. These functions are assumed to satisfy an a priori unknown linear homogeneous differential equation with polynomial coefficients on each continuity interval. Therefore, they may be specified by a finite amount of information. This reconstruction problem has practical importance in Signal Processing and other applications. It is somewhat of a ``folklore'' that the sequence of the moments of such ``piecewise D-finite''functions satisfies a linear recurrence relation of bounded order and degree. We derive this recurrence relation explicitly. It turns out that the coefficients of the differential operator which annihilates every piece of the function, as well as the locations of the discontinuities, appear in this recurrence in a precisely controlled manner. This leads to the formulation of a generic algorithm for reconstructing a piecewise D-finite function from its moments. We investigate the conditions for solvability of the resulting linear systems in the general case, as well as analyze a few particular examples. We provide results of numerical simulations for several types of signals, which test the sensitivity of the proposed algorithm to noise

COSMOGRAIL XVIII: time delays of the quadruply lensed quasar WFI2033-4723

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.Comment: Submitted to Astronomy and Astrophysic

Signals of R-parity violating supersymmetry in neutrino scattering at muon storage rings

Neutrino oscillation signals at muon storage rings can be faked by supersymmetric (SUSY) interactions in an R-parity violating scenario. We investigate the $\tau$-appearance signals for both long-baseline and near-site experiments, and conclude that the latter is of great use in distinguishing between oscillation and SUSY effects. On the other hand, SUSY can cause a manifold increase in the event rate for wrong-sign muons at a long-baseline setting, thereby providing us with signatures of new physics.Comment: 7 pages LaTeX, 4 ps figures, accepted for publication in Phys. Rev.

Neutrino tomography - Learning about the Earth's interior using the propagation of neutrinos

Because the propagation of neutrinos is affected by the presence of Earth matter, it opens new possibilities to probe the Earth's interior. Different approaches range from techniques based upon the interaction of high energy (above TeV) neutrinos with Earth matter, to methods using the MSW effect on the neutrino oscillations of low energy (MeV to GeV) neutrinos. In principle, neutrinos from many different sources (sun, atmosphere, supernovae, beams etc.) can be used. In this talk, we summarize and compare different approaches with an emphasis on more recent developments. In addition, we point out other geophysical aspects relevant for neutrino oscillations.Comment: 22 pages, 9 figures. Proceedings of ``Neutrino sciences 2005: Neutrino geophysics'', December 14-16, 2005, Honolulu, USA. Minor changes, some references added. Final version to appear in Earth, Moon, and Planet

An ultrashort pulse ultra-violet radiation undulator source driven by a laser plasma wakefield accelerator

Narrow band undulator radiation tuneable over the wavelength range of 150–260 nm has been produced by short electron bunches from a 2 mm long laser plasma wakefield accelerator based on a 20 TW femtosecond laser system. The number of photons measured is up to 9 × 106 per shot for a 100 period undulator, with a mean peak brilliance of 1 × 1018 photons/s/mrad2/mm2/0.1% bandwidth. Simulations estimate that the driving electron bunch r.m.s. duration is as short as 3 fs when the electron beam has energy of 120–130 MeV with the radiation pulse duration in the range of 50–100 fs

Neutrino Physics, Superbeams, and the Neutrino Factory

We summarize what has been learned about the neutrino mass spectrum and neutrino mixing, identify interesting open questions that can be answered by accelerator neutrino facilities of the future, and discuss the importance and physics of answering them.Comment: To appear in the Proceedings of the 4th International Workshop on Neutrino Factories (Nu Fact 02). LaTeX, 10 pages, 1 eps figur