A simple closure approximation for slow dynamics of a multiscale system: nonlinear and multiplicative coupling

Multiscale dynamics are ubiquitous in applications of modern science. Because of time scale separation between relatively small set of slowly evolving variables and (typically) much larger set of rapidly changing variables, direct numerical simulations of such systems often require relatively small time discretization step to resolve fast dynamics, which, in turn, increases computational expense. As a result, it became a popular approach in applications to develop a closed approximate model for slow variables alone, which both effectively reduces the dimension of the phase space of dynamics, as well as allows for a longer time discretization step. In this work we develop a new method for approximate reduced model, based on the linear fluctuation-dissipation theorem applied to statistical states of the fast variables. The method is suitable for situations with quadratically nonlinear and multiplicative coupling. We show that, with complex quadratically nonlinear and multiplicative coupling in both slow and fast variables, this method produces comparable statistics to what is exhibited by an original multiscale model. In contrast, it is observed that the results from the simplified closed model with a constant coupling term parameterization are consistently less precise

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

The Hopf algebra structure of the Z3_3-graded quantum supergroup GLq,j(1∣1)_{q,j}(1|1)

In this work, we give some features of the Z$_3$-graded quantum supergroup

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

Photo-Induced Oxidative Stress Impairs Mitochondrial Metabolism in Neurons and Astrocytes

Photodynamic therapy is selective destruction of cells stained with a photosensitizer upon irradiation with light at a specific wavelength in the presence of oxygen. Cell death upon photodynamic treatment is known to occur mainly due to free radical production and subsequent development of oxidative stress. During photodynamic therapy of brain tumors, healthy cells are also damaged; considering this, it is important to investigate the effect of the treatment on normal neurons and glia. We employed live-cell imaging technique to investigate the cellular mechanism of photodynamic action of radachlorin (200 nM) on neurons and astrocytes in primary rat cell culture. We found that the photodynamic effect of radachlorin increases production of reactive oxygen species measured by dihydroethidium and significantly decrease mitochondrial membrane potential. Mitochondrial depolarization was independent of opening of mitochondrial permeability transition pore and was insensitive to blocker of this pore cyclosporine A. However, irradiation of cells with radachlorin dramatically decreased NADH autofluorescence and also reduced mitochondrial NADH pool suggesting inhibition of mitochondrial respiration by limitation of substrate. This effect could be prevented by inhibition of poly (ADP-ribose) polymerase (PARP) with DPQ. Thus, irradiation of neurons and astrocytes in the presence of radachlorin leads to activation of PARP and decrease in NADH that leads to mitochondrial dysfunction