What can be inferred from surrogate data testing?

Surrogate data testing for linearity is frequently applied to confirm the results of nonlinear time series analysis. We argue that this, in general, is not possible.Comment: 1 pag

Global helioseismic evidence for a deeply penetrating Solar meridional flow consisting of multiple flow cells

We use a novel global helioseismic analysis method to infer the meridional flow in the deep Solar interior. The method is based on the perturbation of eigenfunctions of Solar p modes due to meridional flow. We apply this method to time series obtained from Dopplergrams measured by the Michelson Doppler Imager aboard the Solar and Heliospheric Observatory (SOHO) covering the observation period 2004-2010. Our results show evidence that the meridional flow reaches down to the base of the convection zone. The flow profile has a complex spatial structure consisting of multiple flow cells distributed in depth and latitude. Toward the Solar surface, our results are in good agreement with flow measurements from local helioseismology

Increasing insightful thinking in analytic geometry

Elsewhere in this issue Ferdinand Verhulst described the discussion of the interaction of anal- ysis and geometry in the 19th century. In modern times such discussions come up again and again. As of 2014, synthetic geometry will not be part of the Dutch `vwo ¿ mathematics B¿ programme anymore. Instead, the focus will be more on analytic geometry. Mark Timmer and Nellie Verhoef explored possibilities to connect the two disciplines in order to have students look at analytical exercises from a more synthetic point of vie

Electron-hole spectra created by adsorption on metals from density-functional theory

Non-adiabaticity in adsorption on metal surfaces gives rise to a number of measurable effects, such as chemicurrents and exo-electron emission. Here we present a quantitative theory of chemicurrents on the basis of ground-state density-functional theory (DFT) calculations of the effective electronic potential and the Kohn-Sham band structure. Excitation probabilities are calculated both for electron-hole pairs and for electrons and holes separately from first-order time-dependent perturbation theory. This is accomplished by evaluating the matrix elements (between Kohn-Sham states) of the rate of change of the effective electronic potential between subsequent (static) DFT calculations. Our approach is related to the theory of electronic friction, but allows for direct access to the excitation spectra. The method is applied to adsorption of atomic hydrogen isotopes on the Al(111) surface. The results are compatible with the available experimental data (for noble metal surfaces); in particular, the observed isotope effect in H versus D adsorption is described by the present theory. Moreover, the results are in qualitative agreement with computationally elaborate calculations of the full dynamics within time-dependent density-functional theory, with the notable exception of effects due to the spin dynamics. Being a perturbational approach, the method proposed here is simple enough to be applied to a wide class of adsorbates and surfaces, while at the same time allowing us to extract system-specific information.Comment: 23 pages, 9 figures, accepted for publication in Phys. Rev. B, http://prb.aps.org/, v2: some major improvements, plus correction of minor error

Cross-spectral analysis of physiological tremor and muscle activity. I. Theory and application to unsynchronized EMG

We investigate the relationship between the extensor electromyogram (EMG) and tremor time series in physiological hand tremor by cross-spectral analysis. Special attention is directed to the phase spectrum and the effects of observational noise. We calculate the theoretical phase spectrum for a second order linear stochastic process and compare the results to measured tremor data recorded from subjects who did not show a synchronized EMG activity in the corresponding extensor muscle. The results show that physiological tremor is well described by the proposed model and that the measured EMG represents a Newtonian force by which the muscle acts on the hand.Comment: 9 pages, 6 figures, to appear in Biological Cybernetic

Confluence reduction for Markov automata

Markov automata are a novel formalism for specifying systems exhibiting nondeterminism, probabilistic choices and Markovian rates. Recently, the process algebra MAPA was introduced to efficiently model such systems. As always, the state space explosion threatens the analysability of the models generated by such specifications. We therefore introduce confluence reduction for Markov automata, a powerful reduction technique to keep these models small. We define the notion of confluence directly on Markov automata, and discuss how to syntactically detect confluence on the MAPA language as well. That way, Markov automata generated by MAPA specifications can be reduced on-the-fly while preserving divergence-sensitive branching bisimulation. Three case studies demonstrate the significance of our approach, with reductions in analysis time up to an order of magnitude

Selective low concentration ammonia sensing in a microfluidic lab-on-a-chip

In the medical community, there is a considerable interest in a diagnostic breath analyzer for ammonia that is selectively enough to measure in exhaled air and small enough for the small volumes available in such an application. An indirect measurement system for low gaseous ammonia concentrations has been miniaturized and integrated on a chip in order to reach this goal. The detection limit of the system was calculated to be 1.1 parts per billion (ppb). The response time was determined to be 1.6 min with a gas How of 50 ml/min. The required gas volume for one measurement is therefore sufficiently small, although sampling assistance is required for breath analysis. The selectivity of the system is sufficient to measure ammonia concentrations in the low-ppb range. The system is even sufficiently selective to be used in environments that contain elevated carbon dioxide levels, like exhaled air. The lower ammonia concentration expected in diagnostic breath analysis applications, 50 ppb, was demonstrated to be detectable

Cessation of X-ray Pulsation of GX 1+4

We report results from our weekly monitoring campaign on the X-ray pulsar GX 1+4 with the {\em Rossi X-ray Timing Explorer} satellite. The spin-down trend of GX 1+4 was continuing, with the pulsar being at its longest period ever measured (about 138.7 s). At the late stage of the campaign, the source entered an extended faint state, when its X-ray (2-60 keV) flux decreased significantly to an average level of $\sim 3 \times 10^{-10} ergs cm^{-2} s^{-1}$. It was highly variable in the faint state; the flux dropped to as low as $\sim 3 \times 10^{-11} ergs cm^{-2} s^{-1}$. In several observations during this period, the X-ray pulsation became undetectable. We can, therefore, conclude conservatively that the pulsed fraction, which is normally $\gtrsim$ 70% (peak-to-peak), must have decreased drastically in those cases. This is very similar to what was observed of GX 1+4 in 1996 when it became similarly faint in X-ray. In fact, the flux at which the cessation of X-ray pulsation first occurred is nearly the same as it was in 1996. We suggest that we have, once again, observed the propeller effect in GX 1+4, a phenomenon that is predicted by theoretical models of accreting X-ray pulsars.Comment: 13 pages, 9 figures (available at http://www.physics.purdue.edu/~cui/ftp/cuifigs.tar.gz). To appear in Ap

Linear and nonlinear time series analysis of the black hole candidate Cygnus X-1

We analyze the variability in the X-ray lightcurves of the black hole candidate Cygnus X-1 by linear and nonlinear time series analysis methods. While a linear model describes the over-all second order properties of the observed data well, surrogate data analysis reveals a significant deviation from linearity. We discuss the relation between shot noise models usually applied to analyze these data and linear stochastic autoregressive models. We debate statistical and interpretational issues of surrogate data testing for the present context. Finally, we suggest a combination of tools from linear andnonlinear time series analysis methods as a procedure to test the predictions of astrophysical models on observed data.Comment: 15 pages, to appear in Phys. Rev.

Distinguishing Direct from Indirect Interactions in Oscillatory Networks with Multiple Time Scales

Peer reviewedPublisher PD