60-XX PROBABILITY THEORY AND STOCHASTIC PROCESSES (For additional applications, see 11Kxx, 62-XX, 90-XX, 91-XX, 92-XX, 93-XX, 94-XX)
Refine
Year of publication
Document Type
- ZIB-Report (13)
- Article (1)
- Doctoral Thesis (1)
- Master's Thesis (1)
Language
- English (16)
Is part of the Bibliography
- no (16)
Keywords
- Markov State Models (3)
- Mixed-Integer Programming (3)
- NESS (3)
- Non-reversible Markov Processes (3)
- 2-photon microscopy (1)
- Adaptive Importance Sampling (1)
- Bayesian inverse problems (1)
- Girsanov (1)
- Karhunen–Loève expansion (1)
- Metadynamics (1)
Institute
- Numerical Mathematics (9)
- Mathematical Optimization (7)
- Computational Molecular Design (3)
- Mathematical Optimization Methods (2)
- Modeling and Simulation of Complex Processes (2)
- Visual Data Analysis (2)
- Visual and Data-centric Computing (2)
- Image Analysis in Biology and Materials Science (1)
- Uncertainty Quantification (1)
The paper deals with the workload and busy period for the M/GI/1
system under FCFS discipline, where the customers may become
impatient during their waiting for service with generally
distributed maximal waiting times and also during their service
with generally distributed maximal service times depending
on the time waited for service. This general impatience mechanism,
originally introduced by Kovalenko (1961) and considered
by Daley (1965), too, covers the special cases of impatience
on waiting times as well as impatience on sojourn times,
for which Boxma et al. (2010), (2011) gave new results
and outlined special cases recently. Our unified approach
bases on the vector process of workload and busy time.
Explicit representations for the LSTs of workload and
busy period are given in case of phase-type distributed
impatience.
The topic of this thesis is the examination of an optimization model
which stems from the clustering process of non-reversible markov processes.
We introduce the cycle clustering problem und formulate it as a mixed
integer program (MIP).
We prove that this problem is N P-hard and discuss polytopal aspects
such as facets and dimension. The focus of this thesis is the development of
solving methods for this clustering problem. We develop problem specific
primal heuristics, as well as separation methods and an approximation
algorithm. These techniques are implemented in practice as an application
for the MIP solver SCIP.
Our computational experiments show that these solving methods result
in an average speedup of ×4 compared to generic solvers and that our
application is able to solve more instances to optimality within the given
time limit of one hour.
Neurotransmission at chemical synapses relies on the calcium-induced fusion of synaptic vesicles with the presynaptic membrane. The distance to the calcium channels determines the release probability and thereby the postsynaptic signal. Suitable models of the process need to capture both the mean and the variance observed in electrophysiological measurements of the postsynaptic current. In this work, we propose a method to directly compute the exact first- and second-order moments for signals generated by a linear reaction network under convolution with an impulse response function, rendering computationally expensive numerical simulations of the underlying stochastic counting process obsolete. We show that the autocorrelation of the process is central for the calculation of the filtered signal’s second-order moments, and derive a system of PDEs for the cross-correlation functions (including the autocorrelations) of linear reaction networks with time-dependent rates. Finally, we employ our method to efficiently compare different spatial coarse graining approaches for a specific model of synaptic vesicle fusion. Beyond the application to neurotransmission processes, the developed theory can be applied to any linear reaction system that produces a filtered stochastic signal.
Despite the success of constraint programming (CP) for scheduling, the much wider penetration of mixed integer programming (MIP) technology into business applications means that many practical scheduling problems are being addressed with MIP, at least as an initial approach. Furthermore, there has been impressive and well-documented improvements in the power of generic MIP solvers over the past decade.
We empirically demonstrate that on an existing set of resource allocation and scheduling problems standard MIP and CP models are now competitive with the state-of-the-art manual decomposition approach. Motivated by this result, we formulate two tightly coupled hybrid models based on constraint integer programming (CIP) and demonstrate that these models,
which embody advances in CP and MIP, are able to out-perform the CP, MIP, and decomposition models. We conclude that both MIP and CIP
are technologies that should be considered along with CP for solving scheduling problems.
This paper proposes a new method for probabilistic analysis of online algorithms. It is based on the notion of stochastic dominance. We develop the method for
the online bin coloring problem introduced by Krumke et al (2008). Using methods for the stochastic
comparison of Markov chains we establish the result that the performance of the online algorithm GreedyFit is stochastically better than the performance of the algorithm OneBin for any number of items processed. This result gives a more realistic
picture than competitive analysis and explains the behavior observed in simulations.
Given a sequence of Cauchy-distributed random variables defined by a sequence of location parameters and a sequence of scale parameters, we consider another sequence of random variables that is obtained by perturbing the location or scale parameter sequences. Using a result of Kakutani on equivalence of infinite product measures, we provide sufficient conditions for the equivalence of laws of the two sequences.
We deal with an infinite-server system where the
service speed is governed by a stationary and ergodic
process with countably many states. Applying a random
time transformation such that the service speed
becomes one, the sojourn time of a class of virtual
requests with given required service time is equal
in distribution to an additive functional defined
via a stationary version of the time-changed process.
Thus bounds for the expectation of functions of additive
functionals yield bounds for the expectation
of functions of virtual sojourn times, in particular
bounds for fractional moments and the distribution
function. Interpreting the $GI(n)/GI(n)/\infty$ system or
equivalently the $GI(n)/GI$ system under state-dependent
processor sharing as an infinite-server system with
random states given by the number $n$ of requests
in the system provides results for sojourn times
of virtual requests. In case of $M(n)/GI(n)/\infty$,
the sojourn times of arriving and added requests are
equal in distribution to sojourn times of virtual
requests in modified systems, which yields many results
for the sojourn times of arriving and added requests.
In case of integer moments, the bounds generalize
earlier results for $M/GI(n)/\infty$. In particular,
the mean sojourn times of arriving and added requests
in $M(n)/GI(n)/\infty$ are proportional to the required
service time, generalizing Cohen's famous result
for $M/GI(n)/\infty$.
An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates
(2017)
In this article we propose an adaptive importance sampling scheme for dynamical quantities of high dimensional complex systems which are metastable. The main idea of this article is to combine a method coming from Molecular Dynamics Simulation, Metadynamics, with a theorem from stochastic analysis, Girsanov's theorem. The proposed algorithm has two advantages compared to a standard estimator of dynamic quantities: firstly, it is possible to produce estimators with a lower variance and, secondly, we can speed up the sampling. One of the main problems for building importance sampling schemes for metastable systems is to find the metastable region in order to manipulate the potential accordingly. Our method circumvents this problem by using an assimilated version of the Metadynamics algorithm and thus creates a non-equilibrium dynamics which is used to sample the equilibrium quantities.
Some connections between importance sampling and enhanced sampling methods in molecular dynamics
(2017)
Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Schütte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling.
This article extends the framework of Bayesian inverse problems in infinite-dimensional parameter spaces, as advocated by Stuart (Acta Numer. 19:451–559, 2010) and others, to the case of a heavy-tailed prior measure in the family of stable distributions, such as an infinite-dimensional Cauchy distribution, for which polynomial moments are infinite or undefined. It is shown that analogues of the Karhunen–Loève expansion for square-integrable random variables can be used to sample such measures. Furthermore, under weaker regularity assumptions than those used to date, the Bayesian posterior measure is shown to depend Lipschitz continuously in the Hellinger metric upon perturbations of the misfit function and observed data.