Improved linear response for stochastically driven systems

The recently developed short-time linear response algorithm, which predicts the average response of a nonlinear chaotic system with forcing and dissipation to small external perturbation, generally yields high precision of the response prediction, although suffers from numerical instability for long response times due to positive Lyapunov exponents. However, in the case of stochastically driven dynamics, one typically resorts to the classical fluctuation-dissipation formula, which has the drawback of explicitly requiring the probability density of the statistical state together with its derivative for computation, which might not be available with sufficient precision in the case of complex dynamics (usually a Gaussian approximation is used). Here we adapt the short-time linear response formula for stochastically driven dynamics, and observe that, for short and moderate response times before numerical instability develops, it is generally superior to the classical formula with Gaussian approximation for both the additive and multiplicative stochastic forcing. Additionally, a suitable blending with classical formula for longer response times eliminates numerical instability and provides an improved response prediction even for long response times

Cross sections for geodesic flows and \alpha-continued fractions

We adjust Arnoux's coding, in terms of regular continued fractions, of the geodesic flow on the modular surface to give a cross section on which the return map is a double cover of the natural extension for the \alpha-continued fractions, for each $\alpha$ in (0,1]. The argument is sufficiently robust to apply to the Rosen continued fractions and their recently introduced \alpha-variants.Comment: 20 pages, 2 figure

Large closed queueing networks in semi-Markov environment and its application

The paper studies closed queueing networks containing a server station and $k$ client stations. The server station is an infinite server queueing system, and client stations are single-server queueing systems with autonomous service, i.e. every client station serves customers (units) only at random instants generated by a strictly stationary and ergodic sequence of random variables. The total number of units in the network is $N$. The expected times between departures in client stations are $(N\mu_j)^{-1}$. After a service completion in the server station, a unit is transmitted to the $j$th client station with probability $p_{j}$ $(j=1,2,...,k)$, and being processed in the $j$th client station, the unit returns to the server station. The network is assumed to be in a semi-Markov environment. A semi-Markov environment is defined by a finite or countable infinite Markov chain and by sequences of independent and identically distributed random variables. Then the routing probabilities $p_{j}$ $(j=1,2,...,k)$ and transmission rates (which are expressed via parameters of the network) depend on a Markov state of the environment. The paper studies the queue-length processes in client stations of this network and is aimed to the analysis of performance measures associated with this network. The questions risen in this paper have immediate relation to quality control of complex telecommunication networks, and the obtained results are expected to lead to the solutions to many practical problems of this area of research.Comment: 35 pages, 1 figure, 12pt, accepted: Acta Appl. Mat

Multilevel Analysis of Oscillation Motions in Active Regions of the Sun

We present a new method that combines the results of an oscillation study made in optical and radio observations. The optical spectral measurements in photospheric and chromospheric lines of the line-of-sight velocity were carried out at the Sayan Solar Observatory. The radio maps of the Sun were obtained with the Nobeyama Radioheliograph at 1.76 cm. Radio sources associated with the sunspots were analyzed to study the oscillation processes in the chromosphere-corona transition region in the layer with magnetic field B=2000 G. A high level of instability of the oscillations in the optical and radio data was found. We used a wavelet analysis for the spectra. The best similarities of the spectra of oscillations obtained by the two methods were detected in the three-minute oscillations inside the sunspot umbra for the dates when the active regions were situated near the center of the solar disk. A comparison of the wavelet spectra for optical and radio observations showed a time delay of about 50 seconds of the radio results with respect to optical ones. This implies a MHD wave traveling upward inside the umbral magnetic tube of the sunspot. Besides three-minute and five-minute ones, oscillations with longer periods (8 and 15 minutes) were detected in optical and radio records.Comment: 17 pages, 9 figures, accepted to Solar Physics (18 Jan 2011). The final publication is available at http://www.springerlink.co

Multi-level Dynamical Systems: Connecting the Ruelle Response Theory and the Mori-Zwanzig Approach

In this paper we consider the problem of deriving approximate autonomous dynamics for a number of variables of a dynamical system, which are weakly coupled to the remaining variables. In a previous paper we have used the Ruelle response theory on such a weakly coupled system to construct a surrogate dynamics, such that the expectation value of any observable agrees, up to second order in the coupling strength, to its expectation evaluated on the full dynamics. We show here that such surrogate dynamics agree up to second order to an expansion of the Mori-Zwanzig projected dynamics. This implies that the parametrizations of unresolved processes suited for prediction and for the representation of long term statistical properties are closely related, if one takes into account, in addition to the widely adopted stochastic forcing, the often neglected memory effects.Comment: 14 pages, 1 figur

Melatonin prevents cytosolic calcium overload, mitochondrial damage and cell death due to toxically high doses of dexamethasone-induced oxidative stress in human neuroblastoma SH-SY5Y cells

Stressor exposure activates the hypothalamic-pituitary-adrenal (HPA) axis and causes elevations in the levels of glucocorticoids (GC) from the adrenal glands. Increasing evidence has demonstrated that prolonged exposure to high GC levels can lead to oxidative stress, calcium deregulation, mitochondrial dysfunction and apoptosis in a number of cell types. However, melatonin, via its antioxidant activity, exhibits a neuroprotective effect against oxidative stress-induced cell death. Therefore, in the present study, we explored the protective effect of melatonin in GC-induced toxicity in human neuroblastoma SH-SY5Y cells. Cellular treatment with the toxically high doses of the synthetic GC receptor agonist, dexamethasone (DEX) elicited marked decreases in the levels of glutathione and increases in ROS production, lipid peroxidation and cell death. DEX toxicity also induced increases in the levels of cytosolic calcium and mitochondrial fusion proteins (Mfn1 and Opa1) but decreases in the levels of mitochondrial fission proteins (Fis1 and Drp1). Mitochondrial damage was observed in large proportions of the DEX-treated cells. Pretreatment of the cells with melatonin substantially prevented the DEX-induced toxicity. These results suggest that melatonin might exert protective effects against oxidative stress, cytosolic calcium overload and mitochondrial damage in DEX-induced neurotoxicity