CellExcite: an efficient simulation environment for excitable cells

Background Brain, heart and skeletal muscle share similar properties of excitable tissue, featuring both discrete behavior (all-or-nothing response to electrical activation) and continuous behavior (recovery to rest follows a temporal path, determined by multiple competing ion flows). Classical mathematical models of excitable cells involve complex systems of nonlinear differential equations. Such models not only impair formal analysis but also impose high computational demands on simulations, especially in large-scale 2-D and 3-D cell networks. In this paper, we show that by choosing Hybrid Automata as the modeling formalism, it is possible to construct a more abstract model of excitable cells that preserves the properties of interest while reducing the computational effort, thereby admitting the possibility of formal analysis and efficient simulation. Results We have developed CellExcite, a sophisticated simulation environment for excitable-cell networks. CellExcite allows the user to sketch a tissue of excitable cells, plan the stimuli to be applied during simulation, and customize the diffusion model. CellExcite adopts Hybrid Automata (HA) as the computational model in order to efficiently capture both discrete and continuous excitable-cell behavior. Conclusions The CellExcite simulation framework for multicellular HA arrays exhibits significantly improved computational efficiency in large-scale simulations, thus opening the possibility for formal analysis based on HA theory. A demo of CellExcite is available at http://www.cs.sunysb.edu/~eha/ webcite

The P_33(1232) resonance contribution into the amplitudes M_{1+}^{3/2},E_{1+}^{3/2},S_{1+}^{3/2} from an analysis of the p(e,e'p)\pi^0 data at Q^2 = 2.8, 3.2, and 4 (GeV/c)^2 within dispersion relation approach

Within the fixed-t dispersion relation approach we have analysed the TJNAF and DESY data on the exclusive p(e,e'p)\pi^0 reaction in order to find the P_{33}(1232) resonance contribution into the multipole amplitudes M_{1+}^{3/2},E_{1+}^{3/2},S_{1+}^{3/2}. As an input for the resonance and nonresonance contributions into these amplitudes the earlier obtained solutions of the integral equations which follow from dispersion relations are used. The obtained values of the ratio E2/M1 for the \gamma^* N \to P_{33}(1232) transition are: 0.039\pm 0.029, 0.121\pm 0.032, 0.04\pm 0.031 for Q^2= 2.8, 3.2, and 4 (GeV/c)^2, respectively. The comparison with the data at low Q^2 shows that there is no evidence for the presence of the visible pQCD contribution into the transition \gamma N \to P_{33}(1232) at Q^2=3-4 GeV^2. The ratio S_{1+}^{3/2}/M_{1+}^{3/2} for the resonance parts of multipoles is: -0.049\pm 0.029, -0.099\pm 0.041, -0.085\pm 0.021 for Q^2= 2.8, 3.2, and 4 (GeV/c)^2, respectively. Our results for the transverse form factor G_T(Q^2) of the \gamma^* N \to P_{33}(1232) transition are lower than the values obtained from the inclusive data. With increasing Q^2, Q^4G_T(Q^2) decreases, so there is no evidence for the presence of the pQCD contribution here too

Transmission Characteristics of Primate Vocalizations: Implications for Acoustic Analyses

Acoustic analyses have become a staple method in field studies of animal vocal communication, with nearly all investigations using computer-based approaches to extract specific features from sounds. Various algorithms can be used to extract acoustic variables that may then be related to variables such as individual identity, context or reproductive state. Habitat structure and recording conditions, however, have strong effects on the acoustic structure of sound signals. The purpose of this study was to identify which acoustic parameters reliably describe features of propagated sounds. We conducted broadcast experiments and examined the influence of habitat type, transmission height, and re-recording distance on the validity (deviation from the original sound) and reliability (variation within identical recording conditions) of acoustic features of different primate call types. Validity and reliability varied independently of each other in relation to habitat, transmission height, and re-recording distance, and depended strongly on the call type. The smallest deviations from the original sounds were obtained by a visually-controlled calculation of the fundamental frequency. Start- and end parameters of a sound were most susceptible to degradation in the environment. Because the recording conditions can have appreciable effects on acoustic parameters, it is advisable to validate the extraction method of acoustic variables from recordings over longer distances before using them in acoustic analyses

The ADAMTS (A Disintegrin and Metalloproteinase with Thrombospondin motifs) family

The ADAMTS (A Disintegrin and Metalloproteinase with Thrombospondin motifs) enzymes are secreted, multi-domain matrix-associated zinc metalloendopeptidases that have diverse roles in tissue morphogenesis and patho-physiological remodeling, in inflammation and in vascular biology. The human family includes 19 members that can be sub-grouped on the basis of their known substrates, namely the aggrecanases or proteoglycanases (ADAMTS1, 4, 5, 8, 9, 15 and 20), the procollagen N-propeptidases (ADAMTS2, 3 and 14), the cartilage oligomeric matrix protein-cleaving enzymes (ADAMTS7 and 12), the von-Willebrand Factor proteinase (ADAMTS13) and a group of orphan enzymes (ADAMTS6, 10, 16, 17, 18 and 19). Control of the structure and function of the extracellular matrix (ECM) is a central theme of the biology of the ADAMTS, as exemplified by the actions of the procollagen-N-propeptidases in collagen fibril assembly and of the aggrecanases in the cleavage or modification of ECM proteoglycans. Defects in certain family members give rise to inherited genetic disorders, while the aberrant expression or function of others is associated with arthritis, cancer and cardiovascular disease. In particular, ADAMTS4 and 5 have emerged as therapeutic targets in arthritis. Multiple ADAMTSs from different sub-groupings exert either positive or negative effects on tumorigenesis and metastasis, with both metalloproteinase-dependent and -independent actions known to occur. The basic ADAMTS structure comprises a metalloproteinase catalytic domain and a carboxy-terminal ancillary domain, the latter determining substrate specificity and the localization of the protease and its interaction partners; ancillary domains probably also have independent biological functions. Focusing primarily on the aggrecanases and proteoglycanases, this review provides a perspective on the evolution of the ADAMTS family, their links with developmental and disease mechanisms, and key questions for the future

An improved method for measuring muon energy using the truncated mean of dE/dx

The measurement of muon energy is critical for many analyses in large Cherenkov detectors, particularly those that involve separating extraterrestrial neutrinos from the atmospheric neutrino background. Muon energy has traditionally been determined by measuring the specific energy loss (dE/dx) along the muon's path and relating the dE/dx to the muon energy. Because high-energy muons (E_mu > 1 TeV) lose energy randomly, the spread in dE/dx values is quite large, leading to a typical energy resolution of 0.29 in log10(E_mu) for a muon observed over a 1 km path length in the IceCube detector. In this paper, we present an improved method that uses a truncated mean and other techniques to determine the muon energy. The muon track is divided into separate segments with individual dE/dx values. The elimination of segments with the highest dE/dx results in an overall dE/dx that is more closely correlated to the muon energy. This method results in an energy resolution of 0.22 in log10(E_mu), which gives a 26% improvement. This technique is applicable to any large water or ice detector and potentially to large scintillator or liquid argon detectors.Comment: 12 pages, 16 figure

All-particle cosmic ray energy spectrum measured with 26 IceTop stations

We report on a measurement of the cosmic ray energy spectrum with the IceTop air shower array, the surface component of the IceCube Neutrino Observatory at the South Pole. The data used in this analysis were taken between June and October, 2007, with 26 surface stations operational at that time, corresponding to about one third of the final array. The fiducial area used in this analysis was 0.122 km^2. The analysis investigated the energy spectrum from 1 to 100 PeV measured for three different zenith angle ranges between 0{\deg} and 46{\deg}. Because of the isotropy of cosmic rays in this energy range the spectra from all zenith angle intervals have to agree. The cosmic-ray energy spectrum was determined under different assumptions on the primary mass composition. Good agreement of spectra in the three zenith angle ranges was found for the assumption of pure proton and a simple two-component model. For zenith angles {\theta} < 30{\deg}, where the mass dependence is smallest, the knee in the cosmic ray energy spectrum was observed between 3.5 and 4.32 PeV, depending on composition assumption. Spectral indices above the knee range from -3.08 to -3.11 depending on primary mass composition assumption. Moreover, an indication of a flattening of the spectrum above 22 PeV were observed.Comment: 38 pages, 17 figure

The search for transient astrophysical neutrino emission with IceCube-DeepCore

We present the results of a search for astrophysical sources of brief transient neutrino emission using IceCube and DeepCore data acquired between 2012 May 15 and 2013 April 30. While the search methods employed in this analysis are similar to those used in previous IceCube point source searches, the data set being examined consists of a sample of predominantly sub-TeV muon-neutrinos from the Northern Sky (-5 degrees < delta < 90 degrees) obtained through a novel event selection method. This search represents a first attempt by IceCube to identify astrophysical neutrino sources in this relatively unexplored energy range. The reconstructed direction and time of arrival of neutrino events are used to search for any significant self-correlation in the data set. The data revealed no significant source of transient neutrino emission. This result has been used to construct limits at timescales ranging from roughly 1 s to 10 days for generic soft-spectra transients. We also present limits on a specific model of neutrino emission from soft jets in core-collapse supernovae

Search for Relativistic Magnetic Monopoles with IceCube

We present the first results in the search for relativistic magnetic monopoles with the IceCube detector, a subsurface neutrino telescope located in the South Polar ice cap containing a volume of 1 km$^{3}$. This analysis searches data taken on the partially completed detector during 2007 when roughly 0.2 km$^{3}$ of ice was instrumented. The lack of candidate events leads to an upper limit on the flux of relativistic magnetic monopoles of \Phi_{\mathrm{90%C.L.}}\sim 3\e{-18}\fluxunits for $\beta\geq0.8$. This is a factor of 4 improvement over the previous best experimental flux limits up to a Lorentz boost $\gamma$ below $10^{7}$. This result is then interpreted for a wide range of mass and kinetic energy values.Comment: 11 pages, 11 figures. v2 is minor text edits, no changes to resul

Search for non-relativistic Magnetic Monopoles with IceCube

The IceCube Neutrino Observatory is a large Cherenkov detector instrumenting $1\,\mathrm{km}^3$ of Antarctic ice. The detector can be used to search for signatures of particle physics beyond the Standard Model. Here, we describe the search for non-relativistic, magnetic monopoles as remnants of the GUT (Grand Unified Theory) era shortly after the Big Bang. These monopoles may catalyze the decay of nucleons via the Rubakov-Callan effect with a cross section suggested to be in the range of $10^{-27}\,\mathrm{cm^2}$ to $10^{-21}\,\mathrm{cm^2}$. In IceCube, the Cherenkov light from nucleon decays along the monopole trajectory would produce a characteristic hit pattern. This paper presents the results of an analysis of first data taken from May 2011 until May 2012 with a dedicated slow-particle trigger for DeepCore, a subdetector of IceCube. A second analysis provides better sensitivity for the brightest non-relativistic monopoles using data taken from May 2009 until May 2010. In both analyses no monopole signal was observed. For catalysis cross sections of $10^{-22}\,(10^{-24})\,\mathrm{cm^2}$ the flux of non-relativistic GUT monopoles is constrained up to a level of $\Phi_{90} \le 10^{-18}\,(10^{-17})\,\mathrm{cm^{-2}s^{-1}sr^{-1}}$ at a 90% confidence level, which is three orders of magnitude below the Parker bound. The limits assume a dominant decay of the proton into a positron and a neutral pion. These results improve the current best experimental limits by one to two orders of magnitude, for a wide range of assumed speeds and catalysis cross sections.Comment: 20 pages, 20 figure