@misc{Deuflhard, author = {Deuflhard, Peter}, title = {Adaptive Pseudo-transient Continuation for Nonlinear Steady State Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6814}, number = {02-14}, abstract = {Pseudo--transient continuation methods are quite popular for the numerical solution of steady state problems, typically in PDEs. They are based on an embedding into a time dependent initial value problem. In the presence of dynamical invariants the Jacobian matrix of the nonlinear equation system is bound to be singular. The paper presents a convergence analysis which takes this property into account -- in contrast to known approaches. On the basis of the new analysis adaptive algorithms are suggested in detail. These include a variant with Jacobian approximations as well as inexact pseudo--transient continuation, both of which play an important role in discretized PDEs. Numerical experiments are left to future work.}, language = {en} } @misc{Deuflhard, author = {Deuflhard, Peter}, title = {From Molecular Dynamics to Conformational Dynamics in Drug Design}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6878}, number = {02-20}, abstract = {Computational drug design studies molecular recognition in the {\em virtual lab}. The arising Hamiltonian dynamics is known to be chaotic and ill-conditioned already after picoseconds, whereas times are \$msec\$ up to \$min\$. Classical molecular dynamics with long term trajectory computation gives, at best, information about time and statistical ensemble averages. The present paper surveys a recent new modeling approach called {\em conformational dynamics}, which is due to the author and Ch. Sch{\"u}tte. This approach achieves information about the dy time scales by telescoping a short term deterministic model with a statistical model. Examples of small biomolecules are included.}, language = {en} } @misc{Deuflhard, author = {Deuflhard, Peter}, title = {A Comparison of Related Concepts in Computational Chemistry and Mathematics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7138}, number = {02-46}, abstract = {This article studies the relation of the two scientific languages Chemistry and Mathematics via three selected comparisons: (a) QSSA versus dynamic ILDM in reaction kinetics, (b) lumping versus discrete Galerkin methods in polymer chemistry, and (c) geometrical conformations versus metastable conformations in drug design. The common clear message from these comparisons is that chemical intuition may pave the way for mathematical concepts just as chemical concepts may gain from mathematical precisioning. Along this line, significant improvements in chemical research and engineering have already been possible -- and can be further expected in the future from the dialogue between the two scientific languages.}, language = {en} } @misc{Steinbach, author = {Steinbach, Marc}, title = {Topological Index Criteria in DAE for Water Networks}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8821}, number = {05-49}, abstract = {The dynamics of pressurized water distribution networks are naturally modeled by differential algebraic equations (DAE). This paper investigates fundamental structural properties of such a DAE model under weak regularity assumptions. The usual partial derivative-based index-1 condition is shown to be necessary and sufficient for several index concepts, as well as sufficient for solvability in a strong sense. Using the physical properties of nonlinear network elements and the inherent saddle point structure of network hydraulics, we then derive purely topological index criteria based on the network graph and the choice of control variables. Several examples illustrate the theoretical results and explore different non-index-1 situations. A brief discussion of the implications for operative planning by discrete time DAE boundary value problems concludes the paper.}, language = {en} } @misc{Deuflhard, author = {Deuflhard, Peter}, title = {Differential Equations in Technology and Medicine. Computational Concepts, Adaptive Algorithms, and Virtual Labs}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4223}, number = {SC-99-34}, abstract = {This series of lectures has been given to a class of mathematics postdocs at a European summer school on Computational Mathematics Driven by Industrial Applications in Martina Franca, Italy (organized by CIME). It deals with a variety of challenging real life problems selected from clinical cancer therapy, communication technology, polymer production, and pharmaceutical drug design. All of these problems from rather diverse application areas share two common features: (a) they have been modelled by various differential equations -- elliptic, parabolic, or Schr{\"o}dinger--type partial differential equations, countable ordinary diffential equations, or Hamiltonian systems, (b) their numerical solution has turned out to be real challenge to computational mathematics.}, language = {en} } @misc{HuisingaPesceKosloffetal., author = {Huisinga, Wilhelm and Pesce, Lorenzo and Kosloff, Ronnie and Saalfrank, Peter}, title = {Faber and Newton Polynomial Integrators for Open-System Density Matrix Propagation}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3725}, number = {SC-98-29}, abstract = {Two polynomial expansions of the time-evolution superoperator to directly integrate Markovian Liouville-von Neumann (LvN) equations for quantum open systems, namely the Newton interpolation and the Faber approximation, are presented and critically compared. Details on the numerical implementation including error control, and on the performance of either method are given. In a first physical application, a damped harmonic oscillator is considered. Then, the Faber approximation is applied to compute a condensed phase absorption spectrum, for which a semi--analytical expression is derived. Finally, even more general applications are discussed. In all applications considered here it is found that both the Newton and Faber integrators are fast, general, stable, and accurate.}, language = {en} } @misc{LeimkuhlerReich, author = {Leimkuhler, Benedict and Reich, Sebastian}, title = {The Numerical Solution of Constrained Hamiltonian Systems.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-868}, number = {SC-92-16}, abstract = {A Hamiltonian system subject to smooth constraints can typically be viewed as a Hamiltonian system on a manifold. Numerical computations, however, must be performed in \$ R^n\$. In this paper, canonical transformations from ``Hamiltonian differential--algebraic equations'' to ODEs in Euclidean space are considered. In \S2, canonical parameterizations or local charts are developed and it is shown how these can be computed in a practical framework. In \S3 we consider the construction of unconstrained Hamiltonian ODE systems in the space in which the constraint manifold is embedded which preserve the constraint manifold as an integral invariant and whose flow reduces to the flow of the constrained system along the manifold. It is shown that certain of these unconstrained Hamiltonian systems force Lyapunov stability of the constraint--invariants, while others lead to an unstable invariant. In \S4, we compare various projection techniques which might be incorporated to better insure preservation of the constraint--invariants in the context of numerical discretization. Numerical experiments illustrate the degree to which the constraint and symplectic invariants are maintained under discretization of various formulations. {\bf Keywords:} differential--algebraic equations, Hamiltonian systems, canonical discretization schemes. {\bf AMS(MOS):} subject classification 65L05.}, language = {en} } @misc{FischerGoetschelWeiser, author = {Fischer, Lisa and G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Lossy data compression reduces communication time in hybrid time-parallel integrators}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63961}, abstract = {Parallel in time methods for solving initial value problems are a means to increase the parallelism of numerical simulations. Hybrid parareal schemes interleaving the parallel in time iteration with an iterative solution of the individual time steps are among the most efficient methods for general nonlinear problems. Despite the hiding of communication time behind computation, communication has in certain situations a significant impact on the total runtime. Here we present strict, yet no sharp, error bounds for hybrid parareal methods with inexact communication due to lossy data compression, and derive theoretical estimates of the impact of compression on parallel efficiency of the algorithms. These and some computational experiments suggest that compression is a viable method to make hybrid parareal schemes robust with respect to low bandwidth setups.}, language = {en} } @misc{WeiserGhosh, author = {Weiser, Martin and Ghosh, Sunayana}, title = {Theoretically optimal inexact SDC methods}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53140}, abstract = {In several inital value problems with particularly expensive right hand side computation, there is a trade-off between accuracy and computational effort in evaluating the right hand sides. We consider inexact spectral deferred correction (SDC) methods for solving such non-stiff initial value problems. SDC methods are interpreted as fixed point iterations and, due to their corrective iterative nature, allow to exploit the accuracy-work-tradeoff for a reduction of the total computational effort. On one hand we derive an error model bounding the total error in terms of the right hand side evaluation errors. On the other hand, we define work models describing the computational effort in terms of the evaluation accuracy. Combining both, a theoretically optimal tolerance selection is worked out by minimizing the total work subject to achieving the requested tolerance.}, language = {en} } @phdthesis{Reich, author = {Reich, Sebastian}, title = {Dynamical Systems, Numerical Integration, and Exponentially Small Estimates}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3622}, number = {SC-98-19}, abstract = {Dynamical systems with two well-separated time-scales are investigated using normal form theory. Exponential estimates for the normal form truncation error are derived and applied to the numerical integration of differential equations (backward error analysis) and the reduction of highly oscillatory Hamiltonian systems (constrained dynamics and correcting potentials). The theoretical results are used to formulate new algorithms for the time integration of conservative Hamiltonian systems (projected multiple time stepping, soft constraints, rigid bodies, symplectic variable step-size methods).}, language = {en} } @misc{HolderLeimkuhlerReich, author = {Holder, Thomas and Leimkuhler, Benedict and Reich, Sebastian}, title = {Explicit Variable Step-Size and Time-Reversible Integration}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3607}, number = {SC-98-17}, abstract = {A variable step-size, semi-explicit variant of the explicit St{\"o}rmer-Verlet method has been proposed for the time-reversible integration of Newton's equations of motion by Huang \& Leimkuhler. Here we propose a fully explicit version of this approach applicable to explicit and symmetric integration methods for general time-reversible differential equations. As applications, we discuss the variable step-size, time-reversible, and fully explicit integration of rigid body motion and reversible Nos\'e-Hoover dynamics.}, language = {en} } @misc{GelssMateraSchuette, author = {Gelß, Patrick and Matera, Sebastian and Sch{\"u}tte, Christof}, title = {Solving the master equation without kinetic Monte Carlo: tensor train approximations for a CO oxidation model}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55743}, abstract = {In multiscale models of heterogeneous catalysis, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. This usually is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO_2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased stiffness.}, language = {en} }