"Universal" Distribution of Inter-Earthquake Times Explained

We propose a simple theory for the ``universal'' scaling law previously reported for the distributions of waiting times between earthquakes. It is based on a largely used benchmark model of seismicity, which just assumes no difference in the physics of foreshocks, mainshocks and aftershocks. Our theoretical calculations provide good fits to the data and show that universality is only approximate. We conclude that the distributions of inter-event times do not reveal more information than what is already known from the Gutenberg-Richter and the Omori power laws. Our results reinforces the view that triggering of earthquakes by other earthquakes is a key physical mechanism to understand seismicity.Comment: 4 pages with two figure

Spreading of Persistent Infections in Heterogeneous Populations

Up to now, the effects of having heterogeneous networks of contacts have been studied mostly for diseases which are not persistent in time, i.e., for diseases where the infectious period can be considered very small compared to the lifetime of an individual. Moreover, all these previous results have been obtained for closed populations, where the number of individuals does not change during the whole duration of the epidemics. Here, we go one step further and analyze, both analytically and numerically, a radically different kind of diseases: those that are persistent and can last for an individual's lifetime. To be more specific, we particularize to the case of Tuberculosis' (TB) infection dynamics, where the infection remains latent for a period of time before showing up and spreading to other individuals. We introduce an epidemiological model for TB-like persistent infections taking into account the heterogeneity inherent to the population structure. This sort of dynamics introduces new analytical and numerical challenges that we are able to sort out. Our results show that also for persistent diseases the epidemic threshold depends on the ratio of the first two moments of the degree distribution so that it goes to zero in a class of scale-free networks when the system approaches the thermodynamic limit.Comment: 12 pages and 2 figures. Revtex format. Submitted for publication

Mean-Field Interacting Boson Random Point Fields in Weak Harmonic Traps

A model of the mean-field interacting boson gas trapped by a weak harmonic potential is considered by the \textit{boson random point fields} methods. We prove that in the Weak Harmonic Trap (WHT) limit there are two phases distinguished by the boson condensation and by a different behaviour of the local particle density. For chemical potentials less than a certain critical value, the resulting Random Point Field (RPF) coincides with the usual boson RPF, which corresponds to a non-interacting (ideal) boson gas. For the chemical potentials greater than the critical value, the boson RPF describes a divergent (local) density, which is due to \textit{localization} of the macroscopic number of condensed particles. Notice that it is this kind of transition that observed in experiments producing the Bose-Einstein Condensation in traps

SLIMS—a user-friendly sample operations and inventory management system for genotyping labs

Summary: We present the Sample-based Laboratory Information Management System (SLIMS), a powerful and user-friendly open source web application that provides all members of a laboratory with an interface to view, edit and create sample information. SLIMS aims to simplify common laboratory tasks with tools such as a user-friendly shopping cart for subjects, samples and containers that easily generates reports, shareable lists and plate designs for genotyping. Further key features include customizable data views, database change-logging and dynamically filled pre-formatted reports. Along with being feature-rich, SLIMS' power comes from being able to handle longitudinal data from multiple time-points and biological sources. This type of data is increasingly common from studies searching for susceptibility genes for common complex diseases that collect thousands of samples generating millions of genotypes and overwhelming amounts of data. LIMSs provide an efficient way to deal with this data while increasing accessibility and reducing laboratory errors; however, professional LIMS are often too costly to be practical. SLIMS gives labs a feasible alternative that is easily accessible, user-centrically designed and feature-rich. To facilitate system customization, and utilization for other groups, manuals have been written for users and developers

A model for the distribution of aftershock waiting times

In this work the distribution of inter-occurrence times between earthquakes in aftershock sequences is analyzed and a model based on a non-homogeneous Poisson (NHP) process is proposed to quantify the observed scaling. In this model the generalized Omori's law for the decay of aftershocks is used as a time-dependent rate in the NHP process. The analytically derived distribution of inter-occurrence times is applied to several major aftershock sequences in California to confirm the validity of the proposed hypothesis.Comment: 4 pages, 3 figure

A point process framework for modeling electrical stimulation of the auditory nerve

Model-based studies of auditory nerve responses to electrical stimulation can provide insight into the functioning of cochlear implants. Ideally, these studies can identify limitations in sound processing strategies and lead to improved methods for providing sound information to cochlear implant users. To accomplish this, models must accurately describe auditory nerve spiking while avoiding excessive complexity that would preclude large-scale simulations of populations of auditory nerve fibers and obscure insight into the mechanisms that influence neural encoding of sound information. In this spirit, we develop a point process model of the auditory nerve that provides a compact and accurate description of neural responses to electric stimulation. Inspired by the framework of generalized linear models, the proposed model consists of a cascade of linear and nonlinear stages. We show how each of these stages can be associated with biophysical mechanisms and related to models of neuronal dynamics. Moreover, we derive a semi-analytical procedure that uniquely determines each parameter in the model on the basis of fundamental statistics from recordings of single fiber responses to electric stimulation, including threshold, relative spread, jitter, and chronaxie. The model also accounts for refractory and summation effects that influence the responses of auditory nerve fibers to high pulse rate stimulation. Throughout, we compare model predictions to published physiological data and explain differences in auditory nerve responses to high and low pulse rate stimulation. We close by performing an ideal observer analysis of simulated spike trains in response to sinusoidally amplitude modulated stimuli and find that carrier pulse rate does not affect modulation detection thresholds.Comment: 1 title page, 27 manuscript pages, 14 figures, 1 table, 1 appendi

Use of SMS texts for facilitating access to online alcohol interventions: a feasibility study

A41 Use of SMS texts for facilitating access to online alcohol interventions: a feasibility study In: Addiction Science & Clinical Practice 2017, 12(Suppl 1): A4

The iTEBD algorithm beyond unitary evolution

The infinite time-evolving block decimation (iTEBD) algorithm [Phys. Rev. Lett. 98, 070201 (2007)] allows to simulate unitary evolution and to compute the ground state of one-dimensional quantum lattice systems in the thermodynamic limit. Here we extend the algorithm to tackle a much broader class of problems, namely the simulation of arbitrary one-dimensional evolution operators that can be expressed as a (translationally invariant) tensor network. Relatedly, we also address the problem of finding the dominant eigenvalue and eigenvector of a one-dimensional transfer matrix that can be expressed in the same way. New applications include the simulation, in the thermodynamic limit, of open (i.e. master equation) dynamics and thermal states in 1D quantum systems, as well as calculations with partition functions in 2D classical systems, on which we elaborate. The present extension of the algorithm also plays a prominent role in the infinite projected entangled-pair states (iPEPS) approach to infinite 2D quantum lattice systems.Comment: 11 pages, 16 figures, 1 appendix with algorithms for specific types of evolution. A typo in the appendix figures has been corrected. Accepted in PR