TY - JOUR A1 - Menz, Stephan A1 - Latorre, Juan C. A1 - Schütte, Christof A1 - Huisinga, Wilhelm T1 - Hybrid stochastic-deterministic solution of the chemical master equation JF - Multiscale modeling & simulation : a SIAM interdisciplinary journal N2 - The chemical master equation (CME) is the fundamental evolution equation of the stochastic description of biochemical reaction kinetics. In most applications it is impossible to solve the CME directly due to its high dimensionality. Instead, indirect approaches based on realizations of the underlying Markov jump process are used, such as the stochastic simulation algorithm (SSA). In the SSA, however, every reaction event has to be resolved explicitly such that it becomes numerically inefficient when the system's dynamics include fast reaction processes or species with high population levels. In many hybrid approaches, such fast reactions are approximated as continuous processes or replaced by quasi-stationary distributions in either a stochastic or a deterministic context. Current hybrid approaches, however, almost exclusively rely on the computation of ensembles of stochastic realizations. We present a novel hybrid stochastic-deterministic approach to solve the CME directly. Our starting point is a partitioning of the molecular species into discrete and continuous species that induces a partitioning of the reactions into discrete-stochastic and continuous-deterministic processes. The approach is based on a WKB (Wentzel-Kramers-Brillouin) ansatz for the conditional probability distribution function (PDF) of the continuous species (given a discrete state) in combination with Laplace's method of integral approximation. The resulting hybrid stochastic-deterministic evolution equations comprise a CME with averaged propensities for the PDF of the discrete species that is coupled to an evolution equation of the related expected levels of the continuous species for each discrete state. In contrast to indirect hybrid methods, the impact of the evolution of discrete species on the dynamics of the continuous species has to be taken into account explicitly. The proposed approach is efficient whenever the number of discrete molecular species is small. We illustrate the performance of the new hybrid stochastic-deterministic approach in an application to model systems of biological interest. KW - chemical master equation KW - hybrid model KW - multiscale analysis KW - partial averaging KW - asymptotic approximation KW - WKB ansatz Y1 - 2012 U6 - https://doi.org/10.1137/110825716 SN - 1540-3459 VL - 10 IS - 4 SP - 1232 EP - 1262 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - JOUR A1 - Zaourar, N. A1 - Hamoudi, M. A1 - Mandea, M. A1 - Balasis, G. A1 - Holschneider, Matthias T1 - Wavelet-based multiscale analysis of geomagnetic disturbance JF - EARTH PLANETS AND SPACE N2 - The dynamics of external contributions to the geomagnetic field is investigated by applying time-frequency methods to magnetic observatory data. Fractal models and multiscale analysis enable obtaining maximum quantitative information related to the short-term dynamics of the geomagnetic field activity. The stochastic properties of the horizontal component of the transient external field are determined by searching for scaling laws in the power spectra. The spectrum fits a power law with a scaling exponent beta, a typical characteristic of self-affine time-series. Local variations in the power-law exponent are investigated by applying wavelet analysis to the same time-series. These analyses highlight the self-affine properties of geomagnetic perturbations and their persistence. Moreover, they show that the main phases of sudden storm disturbances are uniquely characterized by a scaling exponent varying between 1 and 3, possibly related to the energy contained in the external field. These new findings suggest the existence of a long-range dependence, the scaling exponent being an efficient indicator of geomagnetic activity and singularity detection. These results show that by using magnetogram regularity to reflect the magnetosphere activity, a theoretical analysis of the external geomagnetic field based on local power-law exponents is possible. KW - Geomagnetic field KW - magnetosphere KW - geomagnetic storm KW - multiscale analysis KW - spectral exponent Y1 - 2013 U6 - https://doi.org/10.5047/eps.2013.05.001 SN - 1343-8832 SN - 1880-5981 VL - 65 IS - 12 SP - 1525 EP - 1540 PB - TERRA SCIENTIFIC PUBL CO CY - TOKYO ER - TY - GEN A1 - Zaourar, Naima A1 - Hamoudi, Mohamed A1 - Mandea, Mioara A1 - Balasis, Georgios A1 - Holschneider, Matthias T1 - Wavelet-based multiscale analysis of geomagnetic disturbance T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - The dynamics of external contributions to the geomagnetic field is investigated by applying time-frequency methods to magnetic observatory data. Fractal models and multiscale analysis enable obtaining maximum quantitative information related to the short-term dynamics of the geomagnetic field activity. The stochastic properties of the horizontal component of the transient external field are determined by searching for scaling laws in the power spectra. The spectrum fits a power law with a scaling exponent β, a typical characteristic of self-affine time-series. Local variations in the power-law exponent are investigated by applying wavelet analysis to the same time-series. These analyses highlight the self-affine properties of geomagnetic perturbations and their persistence. Moreover, they show that the main phases of sudden storm disturbances are uniquely characterized by a scaling exponent varying between 1 and 3, possibly related to the energy contained in the external field. These new findings suggest the existence of a long-range dependence, the scaling exponent being an efficient indicator of geomagnetic activity and singularity detection. These results show that by using magnetogram regularity to reflect the magnetosphere activity, a theoretical analysis of the external geomagnetic field based on local power-law exponents is possible. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 947 KW - geomagnetic field KW - magnetosphere KW - geomagnetic storm KW - multiscale analysis KW - spectral exponent Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-436912 SN - 1866-8372 VL - 65 IS - 12 SP - 1525 EP - 1540 ER -