TY - JOUR A1 - Ayanbayev, Birzhan A1 - Klebanov, Ilja A1 - Li, Han Cheng A1 - Sullivan, Tim J. T1 - Gamma-convergence of Onsager-Machlup functionals BT - I. With applications to maximum a posteriori estimation in Bayesian inverse problems JF - Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data N2 - The Bayesian solution to a statistical inverse problem can be summarised by a mode of the posterior distribution, i.e. a maximum a posteriori (MAP) estimator. The MAP estimator essentially coincides with the (regularised) variational solution to the inverse problem, seen as minimisation of the Onsager-Machlup (OM) functional of the posterior measure. An open problem in the stability analysis of inverse problems is to establish a relationship between the convergence properties of solutions obtained by the variational approach and by the Bayesian approach. To address this problem, we propose a general convergence theory for modes that is based on the Gamma-convergence of OM functionals, and apply this theory to Bayesian inverse problems with Gaussian and edge-preserving Besov priors. Part II of this paper considers more general prior distributions. KW - Bayesian inverse problems KW - Gamma-convergence KW - maximum a posteriori KW - estimation KW - Onsager-Machlup functional KW - small ball probabilities; KW - transition path theory Y1 - 2021 U6 - https://doi.org/10.1088/1361-6420/ac3f81 SN - 0266-5611 SN - 1361-6420 VL - 38 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Ayanbayev, Birzhan A1 - Klebanov, Ilja A1 - Lie, Han Cheng A1 - Sullivan, Tim J. T1 - Gamma-convergence of Onsager-Machlup functionals BT - II. Infinite product measures on Banach spaces JF - Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data N2 - We derive Onsager-Machlup functionals for countable product measures on weighted l(p) subspaces of the sequence space R-N. Each measure in the product is a shifted and scaled copy of a reference probability measure on R that admits a sufficiently regular Lebesgue density. We study the equicoercivity and Gamma-convergence of sequences of Onsager-Machlup functionals associated to convergent sequences of measures within this class. We use these results to establish analogous results for probability measures on separable Banach or Hilbert spaces, including Gaussian, Cauchy, and Besov measures with summability parameter 1 <= p <= 2. Together with part I of this paper, this provides a basis for analysis of the convergence of maximum a posteriori estimators in Bayesian inverse problems and most likely paths in transition path theory. KW - Bayesian inverse problems KW - Gamma-convergence KW - maximum a posteriori KW - estimation KW - Onsager-Machlup functional KW - small ball probabilities KW - transition path theory Y1 - 2021 U6 - https://doi.org/10.1088/1361-6420/ac3f82 SN - 0266-5611 SN - 1361-6420 VL - 38 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bandara, Lashi T1 - Functional calculus and harmonic analysis in geometry JF - São Paulo journal of mathematical sciences / Instituto de Matemática e Estatística da Universidade de São Paulo N2 - In this short survey article, we showcase a number of non-trivial geometric problems that have recently been resolved by marrying methods from functional calculus and real-variable harmonic analysis. We give a brief description of these methods as well as their interplay. This is a succinct survey that hopes to inspire geometers and analysts alike to study these methods so that they can be further developed to be potentially applied to a broader range of questions. KW - Functional calculus KW - Real-variable harmonic analysis KW - Elliptic boundary KW - value problems KW - Kato square root problem KW - Spectral flow KW - Riesz topology KW - Gigli-Mantegazza flow KW - Bisectorial operator Y1 - 2021 U6 - https://doi.org/10.1007/s40863-019-00149-0 SN - 1982-6907 SN - 2316-9028 VL - 15 IS - 1 SP - 20 EP - 53 PB - Springer CY - Cham ER - TY - JOUR A1 - Beckus, Siegfried A1 - Eliaz, Latif T1 - Eigenfunctions growth of R-limits on graphs JF - Journal of spectral theory / European Mathematical Society N2 - A characterization of the essential spectrum of Schrodinger operators on infinite graphs is derived involving the concept of R-limits. This concept, which was introduced previously for operators on N and Z(d) as "right-limits," captures the behaviour of the operator at infinity. For graphs with sub-exponential growth rate, we show that each point in sigma(ss)(H) corresponds to a bounded generalized eigenfunction of a corresponding R-limit of H. If, additionally, the graph is of uniform sub-exponential growth, also the converse inclusion holds. KW - Essential spectrum KW - Schrodinger operators KW - graphs KW - right limits KW - generalized eigenfunctions Y1 - 2021 U6 - https://doi.org/10.4171/JST/389 SN - 1664-039X SN - 1664-0403 VL - 11 IS - 4 SP - 1895 EP - 1933 PB - EMS Press, an imprint of the European Mathematical Society - EMS - Publishing House GmbH, Institut für Mathematik, Technische Universität CY - Berlin ER - TY - JOUR A1 - Bär, Christian T1 - The Faddeev-LeVerrier algorithm and the Pfaffian JF - Linear algebra and its applications N2 - We adapt the Faddeev-LeVerrier algorithm for the computation of characteristic polynomials to the computation of the Pfaffian of a skew-symmetric matrix. This yields a very simple, easy to implement and parallelize algorithm of computational cost O(n(beta+1)) where nis the size of the matrix and O(n(beta)) is the cost of multiplying n x n-matrices, beta is an element of [2, 2.37286). We compare its performance to that of other algorithms and show how it can be used to compute the Euler form of a Riemannian manifold using computer algebra. KW - Characteristic polynomial KW - Determinant KW - Pfaffian KW - Gauss-Bonnet-Chern KW - theorem Y1 - 2021 U6 - https://doi.org/10.1016/j.laa.2021.07.023 SN - 0024-3795 SN - 1873-1856 VL - 630 SP - 39 EP - 55 PB - Elsevier CY - New York ER - TY - JOUR A1 - Bär, Christian A1 - Mazzeo, Rafe T1 - Manifolds with many Rarita-Schwinger fields JF - Communications in mathematical physics N2 - The Rarita-Schwinger operator is the twisted Dirac operator restricted to 3/2-spinors. Rarita-Schwinger fields are solutions of this operator which are in addition divergence-free. This is an overdetermined problem and solutions are rare; it is even more unexpected for there to be large dimensional spaces of solutions. In this paper we prove the existence of a sequence of compact manifolds in any given dimension greater than or equal to 4 for which the dimension of the space of Rarita-Schwinger fields tends to infinity. These manifolds are either simply connected Kahler-Einstein spin with negative Einstein constant, or products of such spaces with flat tori. Moreover, we construct Calabi-Yau manifolds of even complex dimension with more linearly independent Rarita-Schwinger fields than flat tori of the same dimension. Y1 - 2021 U6 - https://doi.org/10.1007/s00220-021-04030-0 SN - 0010-3616 SN - 1432-0916 VL - 384 IS - 1 SP - 533 EP - 548 PB - Springer CY - Berlin ER - TY - JOUR A1 - Chang, Der-Chen A1 - Khalil, Sara A1 - Schulze, Bert-Wolfgang T1 - Analysis on regular corner spaces JF - The journal of geometric analysis N2 - We establish a new approach of treating elliptic boundary value problems (BVPs) on manifolds with boundary and regular corners, up to singularity order 2. Ellipticity and parametrices are obtained in terms of symbols taking values in algebras of BVPs on manifolds of corresponding lower singularity orders. Those refer to Boutet de Monvel's calculus of operators with the transmission property, see Boutet de Monvel (Acta Math 126:11-51, 1971) for the case of smooth boundary. On corner configuration operators act in spaces with multiple weights. We mainly study the case of upper left entries in the respective 2 x 2 operator block-matrices of such a calculus. Green operators in the sense of Boutet de Monvel (Acta Math 126:11-51, 1971) analogously appear in singular cases, and they are complemented by contributions of Mellin type. We formulate a result on ellipticity and the Fredholm property in weighted corner spaces, with parametrices of analogous kind. KW - Boutet de Monvel's calculus KW - Pseudo-differential operators KW - Singular cones KW - Mellin symbols with values in the edge calculus KW - Parametrices of elliptic operators KW - Kegel space Y1 - 2021 U6 - https://doi.org/10.1007/s12220-021-00614-3 SN - 1050-6926 SN - 1559-002X VL - 31 IS - 9 SP - 9199 EP - 9240 PB - Springer CY - New York ER - TY - JOUR A1 - Clavier, Pierre J. T1 - Borel-Écalle resummation of a two-point function JF - Annales Henri Poincaré : a journal of theoretical and mathematical physics / ed. jointly by the Institut Henri Poincaré and by the Swiss Physical Society N2 - We provide an overview of the tools and techniques of resurgence theory used in the Borel-ecalle resummation method, which we then apply to the massless Wess-Zumino model. Starting from already known results on the anomalous dimension of the Wess-Zumino model, we solve its renormalisation group equation for the two-point function in a space of formal series. We show that this solution is 1-Gevrey and that its Borel transform is resurgent. The Schwinger-Dyson equation of the model is then used to prove an asymptotic exponential bound for the Borel transformed two-point function on a star-shaped domain of a suitable ramified complex plane. This proves that the two-point function of the Wess-Zumino model is Borel-ecalle summable. Y1 - 2021 U6 - https://doi.org/10.1007/s00023-021-01057-w SN - 1424-0637 SN - 1424-0661 VL - 22 IS - 6 SP - 2103 EP - 2136 PB - Springer CY - Cham ER - TY - JOUR A1 - Cozzoni, Barbara A1 - Maibaum, Michael A1 - Hamm, Maximilian T1 - Thermal analysis and constraints for the MASCOT landing site selection on the asteroid Ryugu JF - Planetary and space science N2 - In June 2018, after 4 years of cruise, the Japanese space probe Hayabusa2 [1-Watanabe S. et al.: Hayabusa2 Mission Overview. (2017)] reached the Near-Earth Asteroid (162173) Ryugu. Hayabusa2 carried a small Lander named MASCOT (Mobile Asteroid Surface Scout) [2-Ho T. M. et al.: MASCOT-The Mobile Asteroid Surface Scout onboard the Hayabusa2 mission. (2017)], jointly developed by the German Aerospace Center (DLR) and the French Space Agency (CNES), to investigate Ryugu's surface structure, composition and physical properties including its thermal behaviour and magnetization in-situ. The Microgravity User Support Centre (DLR-MUSC) in Cologne was in charge of providing all thermal conditions and constraints necessary for the selection of the final landing site and for the final operations of the Lander MASCOT on the surface of the asteroid Ryugu. This article provides a comprehensive assessment of these thermal conditions and constraints, based on predictions performed with the Thermal Mathematical Model (TMM) of MASCOT using different asteroid surface thermal models, ephemeris data for approach as well as descent and hopping trajectories, the related operation sequences and scenarios and the possible environmental conditions driven by the Hayabusa2 spacecraft. A comparison with the real telemetry data confirms the analysis and provides further information about the asteroid characteristics. KW - MASCOT KW - Thermal mathematical model KW - Landing site selection KW - Small KW - spacecraft operations Y1 - 2021 U6 - https://doi.org/10.1016/j.pss.2021.105286 SN - 0032-0633 SN - 1873-5088 VL - 205 PB - Elsevier CY - Oxford ER - TY - JOUR A1 - Cvetković, Nada A1 - Conrad, Tim A1 - Lie, Han Cheng T1 - A convergent discretization method for transition path theory for diffusion processes JF - Multiscale modeling & simulation : a SIAM interdisciplinary journal N2 - Transition path theory (TPT) for diffusion processes is a framework for analyzing the transitions of multiscale ergodic diffusion processes between disjoint metastable subsets of state space. Most methods for applying TPT involve the construction of a Markov state model on a discretization of state space that approximates the underlying diffusion process. However, the assumption of Markovianity is difficult to verify in practice, and there are to date no known error bounds or convergence results for these methods. We propose a Monte Carlo method for approximating the forward committor, probability current, and streamlines from TPT for diffusion processes. Our method uses only sample trajectory data and partitions of state space based on Voronoi tessellations. It does not require the construction of a Markovian approximating process. We rigorously prove error bounds for the approximate TPT objects and use these bounds to show convergence to their exact counterparts in the limit of arbitrarily fine discretization. We illustrate some features of our method by application to a process that solves the Smoluchowski equation on a triple-well potential. KW - ergodic diffusion processes KW - transition paths KW - rare events KW - Monte Carlo KW - methods Y1 - 2021 U6 - https://doi.org/10.1137/20M1329354 SN - 1540-3459 SN - 1540-3467 VL - 19 IS - 1 SP - 242 EP - 266 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - THES A1 - Dahl, Dorothee Sophie T1 - Let's have FUN! Gamification im Mathematikunterricht N2 - Spiele und spieltypische Elemente wie das Sammeln von Treuepunkten sind aus dem Alltag kaum wegzudenken. Zudem werden sie zunehmend in Unternehmen oder in Lernumgebungen eingesetzt. Allerdings ist die Methode Gamification bisher für den pädagogischen Kontext wenig klassifiziert und für Lehrende kaum zugänglich gemacht worden. Daher zielt diese Bachelorarbeit darauf ab, eine systematische Strukturierung und Aufarbeitung von Gamification sowie innovative Ansätze für die Verwendung spieltypischer Elemente im Unterricht, konkret dem Mathematikunterricht, zu präsentieren. Dies kann eine Grundlage für andere Fachgebiete, aber auch andere Lehrformen bieten und so die Umsetzbarkeit von Gamification in eigenen Lehrveranstaltungen aufzeigen. In der Arbeit wird begründet, weshalb und mithilfe welcher Elemente Gamification die Motivation und Leistungsbereitschaft der Lernenden langfristig erhöhen, die Sozial- und Personalkompetenzen fördern sowie die Lernenden zu mehr Aktivität anregen kann. Zudem wird Gamification explizit mit grundlegenden mathematikdidaktischen Prinzipien in Verbindung gesetzt und somit die Relevanz für den Mathematikunterricht hervorgehoben. Anschließend werden die einzelnen Elemente von Gamification wie Punkte, Level, Abzeichen, Charaktere und Rahmengeschichte entlang einer eigens für den pädagogischen Kontext entwickelten Klassifikation „FUN“ (Feedback – User specific elements – Neutral elements) schematisch beschrieben, ihre Funktionen und Wirkung dargestellt sowie Einsatzmöglichkeiten im Unterricht aufgezeigt. Dies beinhaltet Ideen zu lernförderlichem Feedback, Differenzierungsmöglichkeiten und Unterrichtsrahmengestaltung, die in Lehrveranstaltungen aller Art umsetzbar sein können. Die Bachelorarbeit umfasst zudem ein spezifisches Beispiel, einen Unterrichtsentwurf einer gamifizierten Mathematikstunde inklusive des zugehörigen Arbeitsmaterials, anhand dessen die Verwendung von Gamification deutlich wird. Gamification offeriert oftmals Vorteile gegenüber dem traditionellen Unterricht, muss jedoch wie jede Methode an den Inhalt und die Zielgruppe angepasst werden. Weiterführende Forschung könnte sich mit konkreten motivationalen Strukturen, personenspezifischen Unterschieden sowie mit mathematischen Inhalten wie dem Problemlösen oder dem Wechsel zwischen verschiedenen Darstellungen hinsichtlich gamifizierter Lehrformen beschäftigen. N2 - Games and game-typical elements such as collecting points are an indispensable part of everyday life. In addition, they are used increasingly in companies or in learning environments. However, the method of gamification has been little classified for the pedagogical context and it has hardly been made accessible to teachers so far. Therefore, this bachelor’s thesis aims to present a systematic structure and reconditioning of gamification as well as innovative approaches for the implementation of game-typical elements in educational contexts, specifically in teaching mathematics. This thesis can provide a basis for other subject areas, but also for other forms of teaching and thus demonstrate the feasibility of gamification in own courses. The paper explains why and with which elements gamification can increase learners' motivation and willingness to perform in the long term, promote social and personal competences and encourage learners to become more active. Moreover, gamification is explicitly linked to basic mathematics didactic principles and thus emphasizes its relevance for mathematics teaching. Afterwards the individual elements of gamification such as points, levels, badges, characters and frame story are described schematically according to the classification “FUN” (Feedback – User specific elements – Neutral elements), developed especially for the educational context in the thesis. This includes ideas for learn-enhancing feedback, opportunities for differentiation and the design of teaching frameworks that can be implemented in courses of all kinds. The bachelor’s thesis also includes a specific example, a lesson plan for a gamified mathematics lesson including the associated working material, which illustrates the use of gamification. Gamification often offers advantages over traditional teaching, but like any method, it must be adapted to the content and the target group. Further research could focus on specific motivational structures, individual differences of students, and mathematical contents such as problem solving or changing representations regarding gamified teaching. KW - Gamification KW - Spiel KW - Motivation KW - Methode KW - Unterrichtsmethode KW - Feedback KW - Innovation KW - Lernen KW - Mathematikdidaktik KW - Mathematikunterricht KW - gamification KW - game KW - game-based KW - motivation KW - learning KW - feedback KW - method KW - teaching KW - teaching methods KW - didactics of mathematics Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-515937 ER - TY - JOUR A1 - De Oliveira Gomes, André A1 - Högele, Michael Anton T1 - The Kramers problem for SDEs driven by small, accelerated Lévy noise with exponentially light jumps JF - Stochastics and dynamics N2 - We establish Freidlin-Wentzell results for a nonlinear ordinary differential equation starting close to the stable state 0, say, subject to a perturbation by a stochastic integral which is driven by an epsilon-small and (1/epsilon)-accelerated Levy process with exponentially light jumps. For this purpose, we derive a large deviations principle for the stochastically perturbed system using the weak convergence approach developed by Budhiraja, Dupuis, Maroulas and collaborators in recent years. In the sequel, we solve the associated asymptotic first escape problem from the bounded neighborhood of 0 in the limit as epsilon -> 0 which is also known as the Kramers problem in the literature. KW - Freidlin-Wentzell theory KW - large deviations principle KW - accelerated small KW - noise Levy diffusions KW - first passage times KW - first exit location KW - strongly tempered stable Levy measure Y1 - 2021 U6 - https://doi.org/10.1142/S0219493721500192 SN - 0219-4937 SN - 1793-6799 VL - 21 IS - 04 PB - World Scientific CY - Singapore ER - TY - JOUR A1 - Denecke, Klaus-Dieter A1 - Hounnon, Hippolyte T1 - Partial Menger algebras of terms JF - Asian-European journal of mathematics N2 - The superposition operation S-n,S-A, n >= 1, n is an element of N, maps to each (n + 1)-tuple of n-ary operations on a set A an n-ary operation on A and satisfies the so-called superassociative law, a generalization of the associative law. The corresponding algebraic structures are Menger algebras of rank n. A partial algebra of type (n + 1) which satisfies the superassociative law as weak identity is said to be a partial Menger algebra of rank n. As a generalization of linear terms we define r-terms as terms where each variable occurs at most r-times. It will be proved that n-ary r-terms form partial Menger algebras of rank n. In this paper, some algebraic properties of partial Menger algebras such as generating systems, homomorphic images and freeness are investigated. As generalization of hypersubstitutions and linear hypersubstitutions we consider r-hypersubstitutions.U KW - n-ary operation KW - n-ary term KW - superposition of n-ary operations and n-ary KW - terms KW - linear term KW - r-term KW - Menger algebra of rank n KW - partial Menger KW - algebra of rank n KW - r-hypersubstitution Y1 - 2021 U6 - https://doi.org/10.1142/S1793557121500923 SN - 1793-5571 SN - 1793-7183 VL - 14 IS - 06 PB - World Scientific CY - Singapore ER - TY - JOUR A1 - Engbert, Ralf A1 - Rabe, Maximilian Michael A1 - Kliegl, Reinhold A1 - Reich, Sebastian T1 - Sequential data assimilation of the stochastic SEIR epidemic model for regional COVID-19 dynamics JF - Bulletin of mathematical biology : official journal of the Society for Mathematical Biology N2 - Newly emerging pandemics like COVID-19 call for predictive models to implement precisely tuned responses to limit their deep impact on society. Standard epidemic models provide a theoretically well-founded dynamical description of disease incidence. For COVID-19 with infectiousness peaking before and at symptom onset, the SEIR model explains the hidden build-up of exposed individuals which creates challenges for containment strategies. However, spatial heterogeneity raises questions about the adequacy of modeling epidemic outbreaks on the level of a whole country. Here, we show that by applying sequential data assimilation to the stochastic SEIR epidemic model, we can capture the dynamic behavior of outbreaks on a regional level. Regional modeling, with relatively low numbers of infected and demographic noise, accounts for both spatial heterogeneity and stochasticity. Based on adapted models, short-term predictions can be achieved. Thus, with the help of these sequential data assimilation methods, more realistic epidemic models are within reach. KW - Stochastic epidemic model KW - Sequential data assimilation KW - Ensemble Kalman KW - filter KW - COVID-19 Y1 - 2020 U6 - https://doi.org/10.1007/s11538-020-00834-8 SN - 0092-8240 SN - 1522-9602 VL - 83 IS - 1 PB - Springer CY - New York ER - TY - THES A1 - Engelhardt, Max Angel Ronan T1 - Zwischen Simulation und Beweis - eine mathematische Analyse des Bienaymé-Galton-Watson-Prozesses und sein Einsatz innerhalb des Mathematikunterrichts T1 - Between simulation and proof - a mathematical analysis of the Bienaymé-Galton-Watson-process and its application in mathematics lessons N2 - Die Bienaymé-Galton-Watson Prozesse können für die Untersuchung von speziellen und sich entwickelnden Populationen verwendet werden. Die Populationen umfassen Individuen, welche sich identisch, zufällig, selbstständig und unabhängig voneinander fortpflanzen und die jeweils nur eine Generation existieren. Die n-te Generation ergibt sich als zufällige Summe der Individuen der (n-1)-ten Generation. Die Relevanz dieser Prozesse begründet sich innerhalb der Historie und der inner- und außermathematischen Bedeutung. Die Geschichte der Bienaymé-Galton-Watson-Prozesse wird anhand der Entwicklung des Konzeptes bis heute dargestellt. Dabei werden die Wissenschaftler:innen verschiedener Disziplinen angeführt, die Erkenntnisse zu dem Themengebiet beigetragen und das Konzept in ihren Fachbereichen angeführt haben. Somit ergibt sich die außermathematische Signifikanz. Des Weiteren erhält man die innermathematische Bedeutsamkeit mittels des Konzeptes der Verzweigungsprozesse, welches auf die Bienaymé-Galton-Watson Prozesse zurückzuführen ist. Die Verzweigungsprozesse stellen eines der aussagekräftigsten Modelle für die Beschreibung des Populationswachstums dar. Darüber hinaus besteht die derzeitige Wichtigkeit durch die Anwendungsmöglichkeit der Verzweigungsprozesse und der Bienaymé-Galton-Watson Prozesse innerhalb der Epidemiologie. Es werden die Ebola- und die Corona-Pandemie als Anwendungsfelder angeführt. Die Prozesse dienen als Entscheidungsstütze für die Politik und ermöglichen Aussagen über die Auswirkungen von Maßnahmen bezüglich der Pandemien. Neben den Prozessen werden ebenfalls der bedingte Erwartungswert bezüglich diskreter Zufallsvariablen, die wahrscheinlichkeitserzeugende Funktion und die zufällige Summe eingeführt. Die Konzepte vereinfachen die Beschreibung der Prozesse und bilden somit die Grundlage der Betrachtungen. Außerdem werden die benötigten und weiterführenden Eigenschaften der grundlegenden Themengebiete und der Prozesse aufgeführt und bewiesen. Das Kapitel erreicht seinen Höhepunkt bei dem Beweis des Kritikalitätstheorems, wodurch eine Aussage über das Aussterben des Prozesses in verschiedenen Fällen und somit über die Aussterbewahrscheinlichkeit getätigt werden kann. Die Fälle werden anhand der zu erwartenden Anzahl an Nachkommen eines Individuums unterschieden. Es zeigt sich, dass ein Prozess bei einer zu erwartenden Anzahl kleiner gleich Eins mit Sicherheit ausstirbt und bei einer Anzahl größer als Eins, die Population nicht in jedem Fall aussterben muss. Danach werden einzelne Beispiele, wie der linear fractional case, die Population von Fibroblasten (Bindegewebszellen) von Mäusen und die Entstehungsfragestellung der Prozesse, angeführt. Diese werden mithilfe der erlangten Ergebnisse untersucht und einige ausgewählte zufällige Dynamiken werden im nachfolgenden Kapitel simuliert. Die Simulationen erfolgen durch ein in Python erstelltes Programm und werden mithilfe der Inversionsmethode realisiert. Die Simulationen stellen beispielhaft die Entwicklungen in den verschiedenen Kritikalitätsfällen der Prozesse dar. Zudem werden die Häufigkeiten der einzelnen Populationsgrößen in Form von Histogrammen angebracht. Dabei lässt sich der Unterschied zwischen den einzelnen Fällen bestätigen und es wird die Anwendungsmöglichkeit der Bienaymé-Galton-Watson Prozesse bei komplexeren Problemen deutlich. Histogramme bekräftigen, dass die einzelnen Populationsgrößen nur endlich oft vorkommen. Diese Aussage wurde von Galton aufgeworfen und in der Extinktions-Explosions-Dichotomie verwendet. Die dargestellten Erkenntnisse über das Themengebiet und die Betrachtung des Konzeptes werden mit einer didaktischen Analyse abgeschlossen. Die Untersuchung beinhaltet die Berücksichtigung der Fundamentalen Ideen, der Fundamentalen Ideen der Stochastik und der Leitidee „Daten und Zufall“. Dabei ergibt sich, dass in Abhängigkeit der gewählten Perspektive die Anwendung der Bienaymé-Galton-Watson Prozesse innerhalb der Schule plausibel ist und von Vorteil für die Schüler:innen sein kann. Für die Behandlung wird exemplarisch der Rahmenlehrplan für Berlin und Brandenburg analysiert und mit dem Kernlehrplan Nordrhein-Westfalens verglichen. Die Konzeption des Lehrplans aus Berlin und Brandenburg lässt nicht den Schluss zu, dass die Bienaymé-Galton-Watson Prozesse angewendet werden sollten. Es lässt sich feststellen, dass die zugrunde liegende Leitidee nicht vollumfänglich mit manchen Fundamentalen Ideen der Stochastik vereinbar ist. Somit würde eine Modifikation hinsichtlich einer stärkeren Orientierung des Lehrplans an den Fundamentalen Ideen die Anwendung der Prozesse ermöglichen. Die Aussage wird durch die Betrachtung und Übertragung eines nordrhein-westfälischen Unterrichtsentwurfes für stochastische Prozesse auf die Bienaymé-Galton-Watson Prozesse unterstützt. Darüber hinaus werden eine Concept Map und ein Vernetzungspentagraph nach von der Bank konzipiert um diesen Aspekt hervorzuheben. N2 - The Bienaymé-Galton-Watson processes can be used to study special and developing populations. These populations include individuals that reproduce identically, randomly, separately, independently of each other, and which exist only for one generation. The n-th generation is the random sum of the individuals of the (n-1)-th generation. The relevance of these processes is based on their history and their significance in mathematical and extra-mathematical contexts. The history of the Bienaymé-Galton-Watson processes is illustrated by the development of the concept to the present day. Various scientists from different disciplines who have contributed to the topic in their respective fields are listed. This illustrates moreover the significance in extra-mathematical contexts. Furthermore, the inner- mathematical magnitude is obtained by means of the superordinate concept of branching processes, which can be traced back to the Bienaymé-Galton-Watson processes. These branching processes are one of the most significant models for describing population growth. In addition, the current importance arises from the applicability of branching processes and the Bienaymé-Galton-Watson processes within epidemiology. The Ebola and Corona pandemics are mentioned as fields of application. The processes serve as a basis for political decision-making and enable statements made on the impact of pandemic measures. In addition to the processes, the conditional expectation value for discrete random variables, the probability generating function and the random sum are also introduced. These concepts simplify the description of the processes and thus form the basis of the considerations. Also, the required and further properties of the basic topics and processes are listed and demonstrated. The chapter reaches its climax with the proof of the criticality theorem, whereby a statement can be made about the extinction of the process in different cases and thus about the extinction probability. These cases are distinguished based on the expected number of offspring from the individuals. It turns out that a process with an expected number of less than one certainly becomes extinct. On the contrary, a process with a number greater than one does not necessarily has to die out. Individual examples are then given, such as the linear fractional case, the population of fibroblasts (connective tissue cells) of mice and the question of origin. These are investigated using the results obtained and some selected random dynamics are simulated in the following chapter. The simulations are carried out by a Python self-written program and are realized using the inversion method. These simulations exemplify the developments in the different criticality cases of the processes. Besides, the frequencies of the individual population sizes are displayed in the form of histograms. The difference between the individual cases can be confirmed and the analysis of the fibroblasts reveals the applicability of the Bienaymé-Galton-Watson processes to more complex problems. Histograms confirm that the individual population sizes occur only finitely often. This statement was raised by Galton and is used in the extinction-explosion dichotomy. The presented findings about the topic and the consideration of the concept are concluded with an analysis of didactic-background. This involves the fundamental ideas, the fundamental ideas of stochastics and the guiding idea of data and chance. Depending on the chosen perspective, the use of the Bienaymé-Galton-Watson processes within the school is plausible and may be beneficial for the students. For the treatment, the Rahmenlehrplan for Berlin and Brandenburg is analysed and compared with the core curriculum of Nord Rhine-Westphalia as an example. The design of the curriculum of Berlin and Brandenburg does not allow the conclusion of applying the Bienaymé-Galton-Watson processes. It can be seen that the underlying guiding idea is not fully compatible with some fundamental ideas of stochastics. Thus, a modification to the curriculum more oriented towards these fundamental ideas would allow the application of the processes. This statement is supported by the observation and transfer of a North Rhine-Westphalian teaching design for stochastic processes to the Bienaymé-Galton-Watson processes by means of chain letters. In addition, a concept map and a Vernetzungspentagraph by von der Bank are designed to highlight this aspect. KW - Bienaymé-Galton-Watson Prozess KW - Kritikalitätstheorem KW - Verzweigungsprozess KW - Populationen KW - linear fractional case KW - bedingter Erwartungswert KW - zufällige Summe KW - Simulation KW - wahrscheinlichkeitserzeugende Funktion KW - Historie der Verzweigungsprozesse KW - Instabilität des Prozesses KW - Aussterbewahrscheinlichkeit KW - Geometrische Reproduktionsverteilung KW - Fibroblasten KW - Entstehungsfragestellung KW - Fundamentale Ideen KW - Leitidee „Daten und Zufall“ KW - Rahmenlehrplan KW - Markov-Ketten KW - Corona KW - Bienaymé-Galton-Watson process KW - criticality theorem KW - branching process KW - populations KW - linear fractional case KW - conditional expectation value KW - random sum KW - simulation KW - probability generating function KW - history of branching processes KW - instability of the process KW - extinction probability KW - geometric reproduction distribution KW - fibroblasts KW - question of origin KW - fundamental ideas KW - guiding idea “Daten und Zufall” KW - Rahmenlehrplan KW - Markov chains KW - Corona Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-524474 ER - TY - JOUR A1 - Eshghi, Nasim A1 - Mach, Thomas A1 - Reichel, Lothar T1 - New matrix function approximations and quadrature rules based on the Arnoldi process JF - Journal of computational and applied mathematics N2 - The Arnoldi process can be applied to inexpensively approximate matrix functions of the form f (A)v and matrix functionals of the form v*(f (A))*g(A)v, where A is a large square non-Hermitian matrix, v is a vector, and the superscript * denotes transposition and complex conjugation. Here f and g are analytic functions that are defined in suitable regions in the complex plane. This paper reviews available approximation methods and describes new ones that provide higher accuracy for essentially the same computational effort by exploiting available, but generally not used, moment information. Numerical experiments show that in some cases the modifications of the Arnoldi decompositions proposed can improve the accuracy of v*(f (A))*g(A)v about as much as performing an additional step of the Arnoldi process. KW - Arnoldi process KW - Matrix function approximation KW - Quadrature rule Y1 - 2021 U6 - https://doi.org/10.1016/j.cam.2021.113442 SN - 0377-0427 SN - 1879-1778 VL - 391 PB - Elsevier CY - Amsterdam ER - TY - THES A1 - Etzold, Heiko T1 - Neue Zugänge zum Winkelbegriff T1 - New Ways to the Angle Concept BT - Fachdidaktische Entwicklungsforschung zur Ausbildung des Winkelfeldbegriffs bei Schülerinnen und Schülern der vierten Klassenstufe N2 - Die Vielfältigkeit des Winkelbegriffs ist gleichermaßen spannend wie herausfordernd in Hinblick auf seine Zugänge im Mathematikunterricht der Schule. Ausgehend von verschiedenen Vorstellungen zum Winkelbegriff wird in dieser Arbeit ein Lehrgang zur Vermittlung des Winkelbegriffs entwickelt und letztlich in konkrete Umsetzungen für den Schulunterricht überführt. Dabei erfolgt zunächst eine stoffdidaktische Auseinandersetzung mit dem Winkelbegriff, die von einer informationstheoretischen Winkeldefinition begleitet wird. In dieser wird eine Definition für den Winkelbegriff unter der Fragestellung entwickelt, welche Informationen man über einen Winkel benötigt, um ihn beschreiben zu können. So können die in der fachdidaktischen Literatur auftretenden Winkelvorstellungen aus fachmathematischer Perspektive erneut abgeleitet und validiert werden. Parallel dazu wird ein Verfahren beschrieben, wie Winkel – auch unter dynamischen Aspekten – informationstechnisch verarbeitet werden können, so dass Schlussfolgerungen aus der informationstheoretischen Winkeldefinition beispielsweise in dynamischen Geometriesystemen zur Verfügung stehen. Unter dem Gesichtspunkt, wie eine Abstraktion des Winkelbegriffs im Mathematikunterricht vonstatten gehen kann, werden die Grundvorstellungsidee sowie die Lehrstrategie des Aufsteigens vom Abstrakten zum Konkreten miteinander in Beziehung gesetzt. Aus der Verknüpfung der beiden Theorien wird ein grundsätzlicher Weg abgeleitet, wie im Rahmen der Lehrstrategie eine Ausgangsabstraktion zu einzelnen Winkelaspekten aufgebaut werden kann, was die Generierung von Grundvorstellungen zu den Bestandteilen des jeweiligen Winkelaspekts und zum Operieren mit diesen Begriffsbestandteilen ermöglichen soll. Hierfür wird die Lehrstrategie angepasst, um insbesondere den Übergang von Winkelsituationen zu Winkelkontexten zu realisieren. Explizit für den Aspekt des Winkelfeldes werden, anhand der Untersuchung der Sichtfelder von Tieren, Lernhandlungen und Forderungen an ein Lernmodell beschrieben, die Schülerinnen und Schüler bei der Begriffsaneignung unterstützen. Die Tätigkeitstheorie, der die genannte Lehrstrategie zuzuordnen ist, zieht sich als roter Faden durch die weitere Arbeit, wenn nun theoriebasiert Designprinzipien generiert werden, die in die Entwicklung einer interaktiven Lernumgebung münden. Hierzu wird u. a. das Modell der Artifact-Centric Activity Theory genutzt, das das Beziehungsgefüge aus Schülerinnen und Schülern, dem mathematischen Gegenstand und einer zu entwickelnden App als vermittelndes Medium beschreibt, wobei der Einsatz der App im Unterrichtskontext sowie deren regelgeleitete Entwicklung Bestandteil des Modells sind. Gemäß dem Ansatz der Fachdidaktischen Entwicklungsforschung wird die Lernumgebung anschließend in mehreren Zyklen erprobt, evaluiert und überarbeitet. Dabei wird ein qualitatives Setting angewandt, das sich der Semiotischen Vermittlung bedient und untersucht, inwiefern sich die Qualität der von den Schülerinnen und Schülern gezeigten Lernhandlungen durch die Designprinzipien und deren Umsetzung erklären lässt. Am Ende der Arbeit stehen eine finale Version der Designprinzipien und eine sich daraus ergebende Lernumgebung zur Einführung des Winkelfeldbegriffs in der vierten Klassenstufe. N2 - The diversity of the concept »angle« can be both exciting and challenging when looking at how to access it in mathematics education in schools. In this thesis, based on different ideas of the angle concept, a training course for conveying the concept will be developed and translated into concrete implementations for school teaching. First, there will be a didactical subject matter discussion of the angle concept, which will be accompanied by an angle definition from information theory. Through the didactical subject matter discussion, a definition for the angle concept will be developed which is guided by the question of what kind of information about an angle is needed in order to describe it. This way, the diverse ideas of the angle concept discussed in mathematics didactics literature can be once again derived and validated from a mathematical point of view. In parallel, a method will be described of how an angle - even one with dynamic aspects - can be handled in terms of information technology, so that conclusions can be drawn from a definition from information theory for dynamic geometry environments for instance. Considering how abstraction of the angle concept can take place in mathematics education, the Idea of Grundvorstellungen will then be connected to the structural principle of the Ascent From the Abstract to the Concrete. Based on the connection of these two theories, a training course will be developed that aims to construct an initial abstract of certain aspects of the angle concept which, in turn, aims at enabling the generating of Grundvorstellungen towards components of the angle concept and at operating with it. For this, the structural principle will be adapted – specifically to realize the transition from angle situations to angle contexts. For one aspect, the angular field, there will be a description of learning actions and demands on a learning model that supports students’ concept acquisition. The angular field, in this step, will be represented by vision fields of animals. Activity theory, on which the structural principle is based, depicts the recurring theme throughout this thesis when generating design principles that lead towards the development of an interactive learning environment. For this, the Artifact-Centric Activity Theory model will be used in order to describe connections between students, the mathematical topic and the to-be-created app. The use of the app in classroom situations, as well as its rule-governed development, are components of the model. Following a Design-Based Research approach, this learning environment will then go through several cycles of test, evaluation and revision. For this purpose, a qualitative setting will be applied using Semiotic Mediation. It will be used to investigate how far design principles, as well as their implementation, impacts on the quality of student’s learning actions. As an outcome of this thesis, a final version of the design principles and an ensuing learning environment that introduces the concept of »angular field« in grade four teaching will be created. KW - Winkel KW - Tätigkeitstheorie KW - Digitale Werkzeuge KW - Digital Tools KW - Activity Theory KW - Angle Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-504187 ER - TY - JOUR A1 - Fischer, Florian A1 - Keller, Matthias T1 - Riesz decompositions for Schrödinger operators on graphs JF - Journal of mathematical analysis and applications N2 - We study superharmonic functions for Schrodinger operators on general weighted graphs. Specifically, we prove two decompositions which both go under the name Riesz decomposition in the literature. The first one decomposes a superharmonic function into a harmonic and a potential part. The second one decomposes a superharmonic function into a sum of superharmonic functions with certain upper bounds given by prescribed superharmonic functions. As application we show a Brelot type theorem. KW - Potential theory KW - Green's function KW - Schrödinger operator KW - Weighted KW - graph KW - Subcritical KW - Greatest harmonic minorant Y1 - 2021 U6 - https://doi.org/10.1016/j.jmaa.2020.124674 SN - 0022-247X SN - 1096-0813 VL - 495 IS - 1 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Garmendia, Alfonso A1 - Zambon, Marco T1 - Quotients of singular foliations and Lie 2-group actions JF - Journal of noncommutative geometry N2 - Androulidakis-Skandalis (2009) showed that every singular foliation has an associated topological groupoid, called holonomy groupoid. In this note, we exhibit some functorial properties of this assignment: if a foliated manifold (M, FM ) is the quotient of a foliated manifold (P, FP ) along a surjective submersion with connected fibers, then the same is true for the corresponding holonomy groupoids. For quotients by a Lie group action, an analogue statement holds under suitable assumptions, yielding a Lie 2-group action on the holonomy groupoid. KW - Lie groupoid KW - singular foliation KW - fibration Y1 - 2021 U6 - https://doi.org/10.4171/JNCG/434 SN - 1661-6952 SN - 1661-6960 VL - 15 IS - 4 SP - 1251 EP - 1283 PB - EMS Press, an imprint of the European Mathematical Society - EMS - Publishing House GmbH, Institut für Mathematik, Technische Universität Berlin CY - Berlin ER - TY - JOUR A1 - Gottwald, Georg A. A1 - Reich, Sebastian T1 - Combining machine learning and data assimilation to forecast dynamical systems from noisy partial observations JF - Chaos : an interdisciplinary journal of nonlinear science N2 - We present a supervised learning method to learn the propagator map of a dynamical system from partial and noisy observations. In our computationally cheap and easy-to-implement framework, a neural network consisting of random feature maps is trained sequentially by incoming observations within a data assimilation procedure. By employing Takens's embedding theorem, the network is trained on delay coordinates. We show that the combination of random feature maps and data assimilation, called RAFDA, outperforms standard random feature maps for which the dynamics is learned using batch data. Y1 - 2021 U6 - https://doi.org/10.1063/5.0066080 SN - 1054-1500 SN - 1089-7682 VL - 31 IS - 10 PB - AIP CY - Melville ER - TY - JOUR A1 - Gottwald, Georg A. A1 - Reich, Sebastian T1 - Supervised learning from noisy observations BT - Combining machine-learning techniques with data assimilation JF - Physica : D, Nonlinear phenomena N2 - Data-driven prediction and physics-agnostic machine-learning methods have attracted increased interest in recent years achieving forecast horizons going well beyond those to be expected for chaotic dynamical systems. In a separate strand of research data-assimilation has been successfully used to optimally combine forecast models and their inherent uncertainty with incoming noisy observations. The key idea in our work here is to achieve increased forecast capabilities by judiciously combining machine-learning algorithms and data assimilation. We combine the physics-agnostic data -driven approach of random feature maps as a forecast model within an ensemble Kalman filter data assimilation procedure. The machine-learning model is learned sequentially by incorporating incoming noisy observations. We show that the obtained forecast model has remarkably good forecast skill while being computationally cheap once trained. Going beyond the task of forecasting, we show that our method can be used to generate reliable ensembles for probabilistic forecasting as well as to learn effective model closure in multi-scale systems. (C) 2021 Elsevier B.V. All rights reserved. KW - Data-driven modelling KW - Random feature maps KW - Data assimilation Y1 - 2021 U6 - https://doi.org/10.1016/j.physd.2021.132911 SN - 0167-2789 SN - 1872-8022 VL - 423 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Hartung, Niklas A1 - Wahl, Martin A1 - Rastogi, Abhishake A1 - Huisinga, Wilhelm T1 - Nonparametric goodness-of-fit testing for parametric covariate models in pharmacometric analyses JF - CPT: pharmacometrics & systems pharmacology N2 - The characterization of covariate effects on model parameters is a crucial step during pharmacokinetic/pharmacodynamic analyses. Although covariate selection criteria have been studied extensively, the choice of the functional relationship between covariates and parameters, however, has received much less attention. Often, a simple particular class of covariate-to-parameter relationships (linear, exponential, etc.) is chosen ad hoc or based on domain knowledge, and a statistical evaluation is limited to the comparison of a small number of such classes. Goodness-of-fit testing against a nonparametric alternative provides a more rigorous approach to covariate model evaluation, but no such test has been proposed so far. In this manuscript, we derive and evaluate nonparametric goodness-of-fit tests for parametric covariate models, the null hypothesis, against a kernelized Tikhonov regularized alternative, transferring concepts from statistical learning to the pharmacological setting. The approach is evaluated in a simulation study on the estimation of the age-dependent maturation effect on the clearance of a monoclonal antibody. Scenarios of varying data sparsity and residual error are considered. The goodness-of-fit test correctly identified misspecified parametric models with high power for relevant scenarios. The case study provides proof-of-concept of the feasibility of the proposed approach, which is envisioned to be beneficial for applications that lack well-founded covariate models. Y1 - 2021 U6 - https://doi.org/10.1002/psp4.12614 SN - 2163-8306 VL - 10 IS - 6 SP - 564 EP - 576 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Hastermann, Gottfried A1 - Reinhardt, Maria A1 - Klein, Rupert A1 - Reich, Sebastian T1 - Balanced data assimilation for highly oscillatory mechanical systems JF - Communications in applied mathematics and computational science : CAMCoS N2 - Data assimilation algorithms are used to estimate the states of a dynamical system using partial and noisy observations. The ensemble Kalman filter has become a popular data assimilation scheme due to its simplicity and robustness for a wide range of application areas. Nevertheless, this filter also has limitations due to its inherent assumptions of Gaussianity and linearity, which can manifest themselves in the form of dynamically inconsistent state estimates. This issue is investigated here for balanced, slowly evolving solutions to highly oscillatory Hamiltonian systems which are prototypical for applications in numerical weather prediction. It is demonstrated that the standard ensemble Kalman filter can lead to state estimates that do not satisfy the pertinent balance relations and ultimately lead to filter divergence. Two remedies are proposed, one in terms of blended asymptotically consistent time-stepping schemes, and one in terms of minimization-based postprocessing methods. The effects of these modifications to the standard ensemble Kalman filter are discussed and demonstrated numerically for balanced motions of two prototypical Hamiltonian reference systems. KW - data assimilation KW - ensemble Kalman filter KW - balanced dynamics KW - highly KW - oscillatory systems KW - Hamiltonian dynamics KW - geophysics Y1 - 2021 U6 - https://doi.org/10.2140/camcos.2021.16.119 SN - 1559-3940 SN - 2157-5452 VL - 16 IS - 1 SP - 119 EP - 154 PB - Mathematical Sciences Publishers CY - Berkeley ER - TY - JOUR A1 - Hethey, Christoph Philipp A1 - Hartung, Niklas A1 - Wangorsch, Gaby A1 - Weisser, Karin A1 - Huisinga, Wilhelm T1 - Physiology-based toxicokinetic modelling of aluminium in rat and man JF - Archives of toxicology : official journal of EUROTOX N2 - A sufficient quantitative understanding of aluminium (Al) toxicokinetics (TK) in man is still lacking, although highly desirable for risk assessment of Al exposure. Baseline exposure and the risk of contamination severely limit the feasibility of TK studies administering the naturally occurring isotope Al-27, both in animals and man. These limitations are absent in studies with Al-26 as a tracer, but tissue data are limited to animal studies. A TK model capable of inter-species translation to make valid predictions of Al levels in humans-especially in toxicological relevant tissues like bone and brain-is urgently needed. Here, we present: (i) a curated dataset which comprises all eligible studies with single doses of Al-26 tracer administered as citrate or chloride salts orally and/or intravenously to rats and humans, including ultra-long-term kinetic profiles for plasma, blood, liver, spleen, muscle, bone, brain, kidney, and urine up to 150 weeks; and (ii) the development of a physiology-based (PB) model for Al TK after intravenous and oral administration of aqueous Al citrate and Al chloride solutions in rats and humans. Based on the comprehensive curated Al-26 dataset, we estimated substance-dependent parameters within a non-linear mixed-effect modelling context. The model fitted the heterogeneous Al-26 data very well and was successfully validated against datasets in rats and humans. The presented PBTK model for Al, based on the most extensive and diverse dataset of Al exposure to date, constitutes a major advancement in the field, thereby paving the way towards a more quantitative risk assessment in humans. KW - PBTK KW - Toxicokinetics KW - Al-26 KW - Aluminium Y1 - 2021 U6 - https://doi.org/10.1007/s00204-021-03107-y SN - 0340-5761 SN - 1432-0738 VL - 95 IS - 9 SP - 2977 EP - 3000 PB - Springer CY - Berlin ; Heidelberg ER - TY - THES A1 - Hübner, Andrea T1 - Ein multityper Verzweigungsprozess als Modell zur Untersuchung der Ausbreitung von Covid-19 T1 - Modeling the spread of Covid-19 using a multitype branching process N2 - Im Zuge der Covid-19 Pandemie werden zwei Werte täglich diskutiert: Die zuletzt gemeldete Zahl der neu Infizierten und die sogenannte Reproduktionsrate. Sie gibt wieder, wie viele weitere Menschen ein an Corona erkranktes Individuum im Durchschnitt ansteckt. Für die Schätzung dieses Wertes gibt es viele Möglichkeiten - auch das Robert Koch-Institut gibt in seinem täglichen Situationsbericht stets zwei R-Werte an: Einen 4-Tage-R-Wert und einen weniger schwankenden 7-Tage-R-Wert. Diese Arbeit soll eine weitere Möglichkeit vorstellen, einige Aspekte der Pandemie zu modellieren und die Reproduktionsrate zu schätzen. In der ersten Hälfte der Arbeit werden die mathematischen Grundlagen vorgestellt, die man für die Modellierung benötigt. Hierbei wird davon ausgegangen, dass der Leser bereits ein Basisverständnis von stochastischen Prozessen hat. Im Abschnitt Grundlagen werden Verzweigungsprozesse mit einigen Beispielen eingeführt und die Ergebnisse aus diesem Themengebiet, die für diese Arbeit wichtig sind, präsentiert. Dabei gehen wir zuerst auf einfache Verzweigungsprozesse ein und erweitern diese dann auf Verzweigungsprozesse mit mehreren Typen. Um die Notation zu erleichtern, beschränken wir uns auf zwei Typen. Das Prinzip lässt sich aber auf eine beliebige Anzahl von Typen erweitern. Vor allem soll die Wichtigkeit des Parameters λ herausgestellt werden. Dieser Wert kann als durchschnittliche Zahl von Nachfahren eines Individuums interpretiert werden und bestimmt die Dynamik des Prozesses über einen längeren Zeitraum. In der Anwendung auf die Pandemie hat der Parameter λ die gleiche Rolle wie die Reproduktionsrate R. In der zweiten Hälfte dieser Arbeit stellen wir eine Anwendung der Theorie über Multitype Verzweigungsprozesse vor. Professor Yanev und seine Mitarbeiter modellieren in ihrer Veröffentlichung Branching stochastic processes as models of Covid-19 epidemic development die Ausbreitung des Corona Virus' über einen Verzweigungsprozess mit zwei Typen. Wir werden dieses Modell diskutieren und Schätzer daraus ableiten: Ziel ist es, die Reproduktionsrate zu ermitteln. Außerdem analysieren wir die Möglichkeiten, die Dunkelziffer (die Zahl nicht gemeldeter Krankheitsfälle) zu schätzen. Wir wenden die Schätzer auf die Zahlen von Deutschland an und werten diese schließlich aus. N2 - During the Covid-19 pandemic, the discussion about the situation has been dominated by two numbers: the number of daily new infected individuals and the reproduction rate. The latter is the average number of people, one infected individual will infect with the disease. Because the number of registered infected individuals is generally not equal to the actual number of people who carry the Corona virus, many facts about the pandemic have to be estimated and can not be known for certain. Since the reproduction rate is an important parameter to signify the course of the Pandemic, many ways to estimate it have been developed. The Institute of Robert Koch in Germany uses two reproduction rates R in their daily reports: The 4-days-R-value and the less fluctuating 7-days-Rvalue. This master thesis will develop another model to estimate the R-value and other interesting aspects of the pandemic. The first part of this thesis is dedicated to the mathematical foundations needed to understand the model. The reader is expected to already have basic understanding of stochastic processes. In the section Grundlagen we will discuss branching processes and present the results of their theory that are important for our work. We start by introducing simple branching processes and expand the results to multitype branching processes. In service of a simpler notation we will only consider twotype branching processes, but the results can be used for any number of types. The importance of the parameter λ shall be stressed. It can be seen as the average number of descendants of one individual and dictates the dynamic of the process over a long period of time. Applied to the modeling of the pandemic, λ plays the same role as the reproduction rate R. In the second part of this thesis will present an application of the previously developed theory about multitype branching processes. Prof. Yanev and his colleagues modeled in their publication Branching stochastic processes as models of Covid-19 epidemic development the spreading of the Corona virus by using a branching process with two types. We will discuss this model and deduce estimators from it. We want to estimate the reproduction rate and find a way to determine the number of not registered infected individuals. The estimators will be applied to the data from Germany and we will discuss the results. KW - Covid-19 KW - Corona KW - Reproduktionsrate KW - Verzweigungsprozess KW - Modellierung KW - Covid-19 KW - corona virus KW - reproduction rate KW - branching process KW - modeling Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-509225 ER - TY - JOUR A1 - Keller, Matthias A1 - Liu, Shiping A1 - Peyerimhoff, Norbert T1 - A note on eigenvalue bounds for non-compact manifolds JF - Mathematische Nachrichten N2 - In this article we prove upper bounds for the Laplace eigenvalues lambda(k) below the essential spectrum for strictly negatively curved Cartan-Hadamard manifolds. Our bound is given in terms of k(2) and specific geometric data of the manifold. This applies also to the particular case of non-compact manifolds whose sectional curvature tends to -infinity, where no essential spectrum is present due to a theorem of Donnelly/Li. The result stands in clear contrast to Laplacians on graphs where such a bound fails to be true in general. KW - Cheeger inequality KW - eigenvalues KW - Laplacian KW - negative curvature KW - Riemannian manifold Y1 - 2021 U6 - https://doi.org/10.1002/mana.201900209 SN - 0025-584X SN - 1522-2616 VL - 294 IS - 6 SP - 1134 EP - 1139 PB - Wiley-VCH CY - Weinheim ER - TY - JOUR A1 - Kempton, Mark A1 - Münch, Florentin A1 - Yau, Shing-Tung T1 - A homology vanishing theorem for graphs with positive curvature JF - Communications in analysis and geometry N2 - We prove a homology vanishing theorem for graphs with positive Bakry-' Emery curvature, analogous to a classic result of Bochner on manifolds [3]. Specifically, we prove that if a graph has positive curvature at every vertex, then its first homology group is trivial, where the notion of homology that we use for graphs is the path homology developed by Grigor'yan, Lin, Muranov, and Yau [11]. We moreover prove that the fundamental group is finite for graphs with positive Bakry-' Emery curvature, analogous to a classic result of Myers on manifolds [22]. The proofs draw on several separate areas of graph theory, including graph coverings, gain graphs, and cycle spaces, in addition to the Bakry-Emery curvature, path homology, and graph homotopy. The main results follow as a consequence of several different relationships developed among these different areas. Specifically, we show that a graph with positive curvature cannot have a non-trivial infinite cover preserving 3-cycles and 4-cycles, and give a combinatorial interpretation of the first path homology in terms of the cycle space of a graph. Furthermore, we relate gain graphs to graph homotopy and the fundamental group developed by Grigor'yan, Lin, Muranov, and Yau [12], and obtain an alternative proof of their result that the abelianization of the fundamental group of a graph is isomorphic to the first path homology over the integers. Y1 - 2021 UR - https://www.intlpress.com/site/pub/files/_fulltext/journals/cag/2021/0029/0006/CAG-2021-0029-0006-a005.pdf U6 - https://doi.org/10.4310/CAG.2021.v29.n6.a5 SN - 1019-8385 SN - 1944-9992 VL - 29 IS - 6 SP - 1449 EP - 1473 PB - International Press of Boston CY - Somerville ER - TY - JOUR A1 - Klein, Markus A1 - Rosenberger, Elke T1 - The tunneling effect for Schrödinger operators on a vector bundle JF - Analysis and mathematical physics N2 - In the semiclassical limit (h) over bar -> 0, we analyze a class of self-adjoint Schrodinger operators H-(h) over bar = (h) over bar L-2 + (h) over barW + V center dot id(E) acting on sections of a vector bundle E over an oriented Riemannian manifold M where L is a Laplace type operator, W is an endomorphism field and the potential energy V has non-degenerate minima at a finite number of points m(1),... m(r) is an element of M, called potential wells. Using quasimodes of WKB-type near m(j) for eigenfunctions associated with the low lying eigenvalues of H-(h) over bar, we analyze the tunneling effect, i.e. the splitting between low lying eigenvalues, which e.g. arises in certain symmetric configurations. Technically, we treat the coupling between different potential wells by an interaction matrix and we consider the case of a single minimal geodesic (with respect to the associated Agmon metric) connecting two potential wells and the case of a submanifold of minimal geodesics of dimension l + 1. This dimension l determines the polynomial prefactor for exponentially small eigenvalue splitting. KW - Laplace-type operator KW - Vector bundle KW - WKB-expansion KW - Quasimodes KW - Tunneling KW - Spectral gap KW - Complete asymptotics Y1 - 2021 U6 - https://doi.org/10.1007/s13324-021-00485-5 SN - 1664-2368 SN - 1664-235X VL - 11 IS - 2 PB - Springer International Publishing AG CY - Cham (ZG) ER - TY - JOUR A1 - Leung, Tsz Yan A1 - Leutbecher, Martin A1 - Reich, Sebastian A1 - Shepherd, Theodore G. T1 - Forecast verification BT - relating deterministic and probabilistic metrics JF - Quarterly journal of the Royal Meteorological Society N2 - The philosophy of forecast verification is rather different between deterministic and probabilistic verification metrics: generally speaking, deterministic metrics measure differences, whereas probabilistic metrics assess reliability and sharpness of predictive distributions. This article considers the root-mean-square error (RMSE), which can be seen as a deterministic metric, and the probabilistic metric Continuous Ranked Probability Score (CRPS), and demonstrates that under certain conditions, the CRPS can be mathematically expressed in terms of the RMSE when these metrics are aggregated. One of the required conditions is the normality of distributions. The other condition is that, while the forecast ensemble need not be calibrated, any bias or over/underdispersion cannot depend on the forecast distribution itself. Under these conditions, the CRPS is a fraction of the RMSE, and this fraction depends only on the heteroscedasticity of the ensemble spread and the measures of calibration. The derived CRPS-RMSE relationship for the case of perfect ensemble reliability is tested on simulations of idealised two-dimensional barotropic turbulence. Results suggest that the relationship holds approximately despite the normality condition not being met. KW - CRPS KW - ensembles KW - idealised turbulence KW - NWP KW - RMSE KW - verification Y1 - 2021 U6 - https://doi.org/10.1002/qj.4120 SN - 0035-9009 SN - 1477-870X VL - 147 IS - 739 SP - 3124 EP - 3134 PB - Wiley CY - Hoboken ER - TY - THES A1 - Maier, Corinna T1 - Bayesian data assimilation and reinforcement learning for model-informed precision dosing in oncology T1 - Bayes’sche Datenassimilation und Reinforcement Learning für die modellinformierte Präzisionsdosierung in der Onkologie N2 - While patients are known to respond differently to drug therapies, current clinical practice often still follows a standardized dosage regimen for all patients. For drugs with a narrow range of both effective and safe concentrations, this approach may lead to a high incidence of adverse events or subtherapeutic dosing in the presence of high patient variability. Model-informedprecision dosing (MIPD) is a quantitative approach towards dose individualization based on mathematical modeling of dose-response relationships integrating therapeutic drug/biomarker monitoring (TDM) data. MIPD may considerably improve the efficacy and safety of many drug therapies. Current MIPD approaches, however, rely either on pre-calculated dosing tables or on simple point predictions of the therapy outcome. These approaches lack a quantification of uncertainties and the ability to account for effects that are delayed. In addition, the underlying models are not improved while applied to patient data. Therefore, current approaches are not well suited for informed clinical decision-making based on a differentiated understanding of the individually predicted therapy outcome. The objective of this thesis is to develop mathematical approaches for MIPD, which (i) provide efficient fully Bayesian forecasting of the individual therapy outcome including associated uncertainties, (ii) integrate Markov decision processes via reinforcement learning (RL) for a comprehensive decision framework for dose individualization, (iii) allow for continuous learning across patients and hospitals. Cytotoxic anticancer chemotherapy with its major dose-limiting toxicity, neutropenia, serves as a therapeutically relevant application example. For more comprehensive therapy forecasting, we apply Bayesian data assimilation (DA) approaches, integrating patient-specific TDM data into mathematical models of chemotherapy-induced neutropenia that build on prior population analyses. The value of uncertainty quantification is demonstrated as it allows reliable computation of the patient-specific probabilities of relevant clinical quantities, e.g., the neutropenia grade. In view of novel home monitoring devices that increase the amount of TDM data available, the data processing of sequential DA methods proves to be more efficient and facilitates handling of the variability between dosing events. By transferring concepts from DA and RL we develop novel approaches for MIPD. While DA-guided dosing integrates individualized uncertainties into dose selection, RL-guided dosing provides a framework to consider delayed effects of dose selections. The combined DA-RL approach takes into account both aspects simultaneously and thus represents a holistic approach towards MIPD. Additionally, we show that RL can be used to gain insights into important patient characteristics for dose selection. The novel dosing strategies substantially reduce the occurrence of both subtherapeutic and life-threatening neutropenia grades in a simulation study based on a recent clinical study (CEPAC-TDM trial) compared to currently used MIPD approaches. If MIPD is to be implemented in routine clinical practice, a certain model bias with respect to the underlying model is inevitable, as the models are typically based on data from comparably small clinical trials that reflect only to a limited extent the diversity in real-world patient populations. We propose a sequential hierarchical Bayesian inference framework that enables continuous cross-patient learning to learn the underlying model parameters of the target patient population. It is important to note that the approach only requires summary information of the individual patient data to update the model. This separation of the individual inference from population inference enables implementation across different centers of care. The proposed approaches substantially improve current MIPD approaches, taking into account new trends in health care and aspects of practical applicability. They enable progress towards more informed clinical decision-making, ultimately increasing patient benefits beyond the current practice. N2 - Obwohl Patienten sehr unterschiedlich auf medikamentöse Therapien ansprechen, werden in der klinischen Praxis häufig noch standardisierte Dosierungsschemata angewendet. Bei Arzneimitteln mit engen therapeutischen Fenstern zwischen minimal wirksamen und toxischen Konzentrationen kann dieser Ansatz bei hoher interindividueller Variabilität zu häufigem Auftreten von Toxizitäten oder subtherapeutischen Konzentrationen führen. Die modellinformierte Präzisionsdosierung (MIPD) ist ein quantitativer Ansatz zur Dosisindividualisierung, der auf der mathematischen Modellierung von Dosis-Wirkungs-Beziehungen beruht und Daten aus dem therapeutischen Drug/Biomarker-Monitoring (TDM) einbezieht. Die derzeitigen MIPD-Ansätze verwenden entweder Dosierungstabellen oder einfache Punkt-Vorhersagen des Therapieverlaufs. Diesen Ansätzen fehlt eine Quantifizierung der Unsicherheiten, verzögerte Effekte werden nicht berücksichtigt und die zugrunde liegenden Modelle werden im Laufe der Anwendung nicht verbessert. Daher sind die derzeitigen Ansätze nicht ideal für eine fundierte klinische Entscheidungsfindung auf Grundlage eines differenzierten Verständnisses des individuell vorhergesagten Therapieverlaufs. Das Ziel dieser Arbeit ist es, mathematische Ansätze für das MIPD zu entwickeln, die (i) eine effiziente, vollständig Bayes’sche Vorhersage des individuellen Therapieverlaufs einschließlich der damit verbundenen Unsicherheiten ermöglichen, (ii) Markov-Entscheidungsprozesse mittels Reinforcement Learning (RL) in einen umfassenden Entscheidungsrahmen zur Dosisindividualisierung integrieren, und (iii) ein kontinuierliches Lernen zwischen Patienten erlauben. Die antineoplastische Chemotherapie mit ihrer wichtigen dosislimitierenden Toxizität, der Neutropenie, dient als therapeutisch relevantes Anwendungsbeispiel. Für eine umfassendere Therapievorhersage wenden wir Bayes’sche Datenassimilationsansätze (DA) an, um TDM-Daten in mathematische Modelle der Chemotherapie-induzierten Neutropenie zu integrieren. Wir zeigen, dass die Quantifizierung von Unsicherheiten einen großen Mehrwert bietet, da sie eine zuverlässige Berechnung der Wahrscheinlichkeiten relevanter klinischer Größen, z.B. des Neutropeniegrades, ermöglicht. Im Hinblick auf neue Home-Monitoring-Geräte, die die Anzahl der verfügbaren TDM-Daten erhöhen, erweisen sich sequenzielle DA-Methoden als effizienter und erleichtern den Umgang mit der Unsicherheit zwischen Dosierungsereignissen. Basierend auf Konzepten aus DA und RL, entwickeln wir neue Ansätze für MIPD. Während die DA-geleitete Dosierung individualisierte Unsicherheiten in die Dosisauswahl integriert, berücksichtigt die RL-geleitete Dosierung verzögerte Effekte der Dosisauswahl. Der kombinierte DA-RL-Ansatz vereint beide Aspekte und stellt somit einen ganzheitlichen Ansatz für MIPD dar. Zusätzlich zeigen wir, dass RL Informationen über die für die Dosisauswahl relevanten Patientencharakteristika liefert. Der Vergleich zu derzeit verwendeten MIPD Ansätzen in einer auf einer klinischen Studie (CEPAC-TDM-Studie) basierenden Simulationsstudie zeigt, dass die entwickelten Dosierungsstrategien das Auftreten subtherapeutischer Konzentrationen sowie lebensbedrohlicher Neutropenien drastisch reduzieren. Wird MIPD in der klinischen Routine eingesetzt, ist eine gewisse Modellverzerrung unvermeidlich. Die Modelle basieren in der Regel auf Daten aus vergleichsweise kleinen klinischen Studien, die die Heterogenität realer Patientenpopulationen nur begrenzt widerspiegeln. Wir schlagen einen sequenziellen hierarchischen Bayes’schen Inferenzrahmen vor, der ein kontinuierliches patientenübergreifendes Lernen ermöglicht, um die zugrunde liegenden Modellparameter der Ziel-Patientenpopulation zu erlernen. Zur Aktualisierung des Modells erfordert dieser Ansatz lediglich zusammenfassende Informationen der individuellen Patientendaten, was eine Umsetzung über verschiedene Versorgungszentren hinweg erlaubt. Die vorgeschlagenen Ansätze verbessern die derzeitigen MIPD-Ansätze erheblich, wobei neue Trends in der Gesundheitsversorgung und Aspekte der praktischen Anwendbarkeit berücksichtigt werden. Damit stellen sie einen Fortschritt in Richtung einer fundierteren klinischen Entscheidungsfindung dar. KW - data assimilation KW - Datenassimilation KW - reinforcement learning KW - model-informed precision dosing KW - pharmacometrics KW - oncology KW - modellinformierte Präzisionsdosierung KW - Onkologie KW - Pharmakometrie KW - Reinforcement Learning Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-515870 ER - TY - JOUR A1 - Matzka, Jürgen A1 - Stolle, Claudia A1 - Yamazaki, Yosuke A1 - Bronkalla, Oliver A1 - Morschhauser, Achim T1 - The geomagnetic Kp index and derived indices of geomagnetic activity JF - Space weather : the international journal of research and applications N2 - The geomagnetic Kp index is one of the most extensively used indices of geomagnetic activity, both for scientific and operational purposes. This article reviews the properties of the Kp index and provides a reference for users of the Kp index and associated data products as derived and distributed by the GFZ German Research Centre for Geosciences. The near real-time production of the nowcast Kp index is of particular interest for space weather services and here we describe and evaluate its current setup. Y1 - 2021 U6 - https://doi.org/10.1029/2020SW002641 SN - 1542-7390 VL - 19 IS - 5 PB - Wiley CY - New York ER - TY - THES A1 - Möhring, Jan T1 - Stochastic inversion for core field modeling using satellite data N2 - Magnetfeldmodellierung mit Kugelflächenfunktionen basiert auf der Inversion nach hunderten bis tausenden von Parametern. Dieses hochdimensionale Problem kann grundsätzlich als ein Optimierungsproblem formuliert werden, bei dem ein globales Minimum einer gewissen Zielfunktion berechnet werden soll. Um dieses Problem zu lösen, gibt es eine Reihe bekannter Ansätze, dazu zählen etwa gradientenbasierte Verfahren oder die Methode der kleinsten Quadrate und deren Varianten. Jede dieser Methoden hat verschiedene Vor- und Nachteile, beispielsweise bezüglich der Anwendbarkeit auf nicht-differenzierbare Funktionen oder der Laufzeit zugehöriger Algorithmen. In dieser Arbeit verfolgen wir das Ziel, einen Algorithmus zu finden, der schneller als die etablierten Verfahren ist und sich auch für nichtlineare Probleme anwenden lässt. Solche nichtlinearen Probleme treten beispielsweise bei der Abschätzung von Euler-Winkeln oder bei der Verwendung der robusteren L_1-Norm auf. Dazu untersuchen wir die Anwendbarkeit stochastischer Optimierungsverfahren aus der CMAES-Familie auf die Modellierung des geomagnetischen Feldes des Erdkerns. Es werden sowohl die Grundlagen der Kernfeldmodellierung und deren Parametrisierung anhand einiger Beispiele aus der Literatur besprochen, als auch die theoretischen Hintergründe der stochastischen Verfahren gegeben. Ein CMAES-Algorithmus wurde erfolgreich angewendet, um Daten der Swarm-Satellitenmission zu invertieren und daraus das Magnetfeldmodell EvoMag abzuleiten. EvoMag zeigt gute Übereinstimmung mit etablierten Modellen, sowie mit Observatoriumsdaten aus Niemegk. Wir thematisieren einige beobachtete Schwierigkeiten und präsentieren und diskutieren die Ergebnisse unserer Modellierung. N2 - Geomagnetic field modeling using spherical harmonics requires the inversion for hundreds to thousands of parameters. This large-scale problem can always be formulated as an optimization problem, where a global minimum of a certain cost function has to be calculated. A variety of approaches is known in order to solve this inverse problem, e.g. derivative-based methods or least-squares methods and their variants. Each of these methods has its own advantages and disadvantages, which affect for example the applicability to non-differentiable functions or the runtime of the corresponding algorithm. In this work, we pursue the goal to find an algorithm which is faster than the established methods and which is applicable to non-linear problems. Such non-linear problems occur for example when estimating Euler angles or when the more robust L_1 norm is applied. Therefore, we will investigate the usability of stochastic optimization methods from the CMAES family for modeling the geomagnetic field of Earth's core. On one hand, basics of core field modeling and their parameterization are discussed using some examples from the literature. On the other hand, the theoretical background of the stochastic methods are provided. A specific CMAES algorithm was successfully applied in order to invert data of the Swarm satellite mission and to derive the core field model EvoMag. The EvoMag model agrees well with established models and observatory data from Niemegk. Finally, we present some observed difficulties and discuss the results of our model. T2 - Stochastische Inversion für Kernfeldmodellierung mit Satellitendaten KW - Geomagnetismus KW - Kernfeldmodellierung KW - Optimierung KW - Evolutionsstrategien KW - Inverse Probleme KW - Geomagnetism KW - Core Field Modeling KW - Optimization KW - Evolution Strategies KW - Inverse Problems Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-498072 ER - TY - THES A1 - Oancea, Marius-Adrian T1 - Spin Hall effects in general relativity T1 - Spin Hall Effekte in der Allgemeinen Relativitätstheorie N2 - The propagation of test fields, such as electromagnetic, Dirac or linearized gravity, on a fixed spacetime manifold is often studied by using the geometrical optics approximation. In the limit of infinitely high frequencies, the geometrical optics approximation provides a conceptual transition between the test field and an effective point-particle description. The corresponding point-particles, or wave rays, coincide with the geodesics of the underlying spacetime. For most astrophysical applications of interest, such as the observation of celestial bodies, gravitational lensing, or the observation of cosmic rays, the geometrical optics approximation and the effective point-particle description represent a satisfactory theoretical model. However, the geometrical optics approximation gradually breaks down as test fields of finite frequency are considered. In this thesis, we consider the propagation of test fields on spacetime, beyond the leading-order geometrical optics approximation. By performing a covariant Wentzel-Kramers-Brillouin analysis for test fields, we show how higher-order corrections to the geometrical optics approximation can be considered. The higher-order corrections are related to the dynamics of the spin internal degree of freedom of the considered test field. We obtain an effective point-particle description, which contains spin-dependent corrections to the geodesic motion obtained using geometrical optics. This represents a covariant generalization of the well-known spin Hall effect, usually encountered in condensed matter physics and in optics. Our analysis is applied to electromagnetic and massive Dirac test fields, but it can easily be extended to other fields, such as linearized gravity. In the electromagnetic case, we present several examples where the gravitational spin Hall effect of light plays an important role. These include the propagation of polarized light rays on black hole spacetimes and cosmological spacetimes, as well as polarization-dependent effects on the shape of black hole shadows. Furthermore, we show that our effective point-particle equations for polarized light rays reproduce well-known results, such as the spin Hall effect of light in an inhomogeneous medium, and the relativistic Hall effect of polarized electromagnetic wave packets encountered in Minkowski spacetime. N2 - Unser grundlegendes Verständnis des Universums basiert auf Einsteins allgemeiner Relativitätstheorie, die eine Beschreibung in Form einer vierdimensional gekrümmten Raumzeit liefert, in der die Anziehungskraft der Gravitation in der Krümmung der Raumzeit kodiert ist. Die überwiegende Mehrheit der experimentellen Tests, die Einsteins allgemeine Relativitätstheorie bestätigt haben, basiert auf der Beobachtung elektromagnetischer Strahlung, die von entfernten astrophysikalischen Quellen wie Sternen oder Galaxien stammt. Daher ist ein tiefgreifendes Verständnis der Dynamik der sich in der Raumzeit ausbreitenden elektromagnetischen Strahlung von entscheidender Bedeutung. Elektromagnetische Phänomene werden durch Maxwell-Gleichungen beschrieben. Die Ausbreitung elektromagnetischer Strahlung in der Raumzeit ist jedoch sehr komplexe, und es ist im Allgemeinen nützlich, Näherungen zu betrachten, welche eine vereinfachte Beschreibung liefern. Auf diese Weise können die Haupteigenschaften des Systems in einem reduzierten Gleichungssystem codiert und die Gültigkeit der Näherung quantitativ kontrolliert werden. Beispielsweise kann die Ausbreitung elektromagnetischer Strahlung in der Raumzeit durch Anwendung der geometrischen Optik auf die Maxwell-Gleichungen beschrieben werden. Diese liefert ein Modell für die Ausbreitung elektromagnetischer Strahlung in Form von Lichtstrahlen, die sich auf dem kürzesten Weg zwischen zwei Punkten ausbreiten. Im Kontext von Einsteins allgemeiner Relativitätstheorie entsprechen dise Lichtstrahlen den Nullgeodäten der zugrunde liegenden gekrümmten Raumzeit. Für die meisten astrophysikalischen Anwendungen von Interesse, wie die Beobachtung von Himmelskörpern oder Gravitationslinsen, stellen die Näherungen der geometrischen Optik und damit die Beschreibung der Ausbreitung elektromagnetischer Strahlung durch Lichtstrahlen ein zufriedenstellendes theoretisches Modell dar. In dieser Arbeit untersuchen wir mögliche Korrekturen der Ausbreitung elektromagnetischer Strahlung in der Raumzeit, welche durch die Näherung der geometrischen Optik nicht erfasst werden. Solche Korrekturen sind aus der Optik bekannt, wo beobachtet wurde, dass die Ausbreitung von Lichtstrahlen in bestimmten Materialien durch die Polarisation des Lichts beeinflusst werden kann. Diese Korrekturen sind als Spin-Hall-Effekt von Licht bekannt. In dieser Arbeit wird gezeigt, dass ein ähnlicher Effekt für elektromagnetische Strahlung auftreten kann, welche sich in gekrümmter Raumzeit in der Nähe massiver astrophysikalischer Objekte wie Schwarzer Löcher oder Sterne ausbreitet. Darüber hinaus präsentieren wir, basierend auf der Dirac-Gleichung, eine ähnliche Analyse für die Bewegung von Elektronen in gekrümmten Raumzeiten. KW - spin Hall effect KW - gravitation KW - black hole KW - Schwarzes Loch KW - Gravitation KW - Spin Hall effekte Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-502293 ER - TY - JOUR A1 - Oster, Mathias A1 - Dias, Marcelo A. A1 - Wolff, Timo de A1 - Evans, Myfanwy T1 - Reentrant tensegrity BT - a three-periodic, chiral, tensegrity structure that is auxetic JF - Science advances / American Association for the Advancement of Science N2 - We present a three-periodic, chiral, tensegrity structure and demonstrate that it is auxetic. Our tensegrity structure is constructed using the chiral symmetry Pi(+) cylinder packing, transforming cylinders to elastic elements and cylinder contacts to incompressible rods. The resulting structure displays local reentrant geometry at its vertices and is shown to be auxetic when modeled as an equilibrium configuration of spatial constraints subject to a quasi-static deformation. When the structure is subsequently modeled as a lattice material with elastic elements, the auxetic behavior is again confirmed through finite element modeling. The cubic symmetry of the original structure means that the auxetic behavior is observed in both perpendicular directions and is close to isotropic in magnitude. This structure could be the simplest three-dimensional analog to the two-dimensional reentrant honeycomb. This, alongside the chirality of the structure, makes it an interesting design target for multifunctional materials. Y1 - 2021 U6 - https://doi.org/10.1126/sciadv.abj6737 SN - 2375-2548 VL - 7 IS - 50 PB - American Association for the Advancement of Science CY - Washington ER - TY - JOUR A1 - Pathiraja, Sahani Darschika A1 - Reich, Sebastian A1 - Stannat, Wilhelm T1 - McKean-Vlasov SDEs in nonlinear filtering JF - SIAM journal on control and optimization : a publication of the Society for Industrial and Applied Mathematics N2 - Various particle filters have been proposed over the last couple of decades with the common feature that the update step is governed by a type of control law. This feature makes them an attractive alternative to traditional sequential Monte Carlo which scales poorly with the state dimension due to weight degeneracy. This article proposes a unifying framework that allows us to systematically derive the McKean-Vlasov representations of these filters for the discrete time and continuous time observation case, taking inspiration from the smooth approximation of the data considered in [D. Crisan and J. Xiong, Stochastics, 82 (2010), pp. 53-68; J. M. Clark and D. Crisan, Probab. Theory Related Fields, 133 (2005), pp. 43-56]. We consider three filters that have been proposed in the literature and use this framework to derive Ito representations of their limiting forms as the approximation parameter delta -> 0. All filters require the solution of a Poisson equation defined on R-d, for which existence and uniqueness of solutions can be a nontrivial issue. We additionally establish conditions on the signal-observation system that ensures well-posedness of the weighted Poisson equation arising in one of the filters. KW - data assimilation KW - feedback particle filter KW - Poincare inequality KW - well-posedness KW - nonlinear filtering KW - McKean-Vlasov KW - mean-field equations Y1 - 2022 U6 - https://doi.org/10.1137/20M1355197 SN - 0363-0129 SN - 1095-7138 VL - 59 IS - 6 SP - 4188 EP - 4215 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - THES A1 - Perera, Upeksha T1 - Solutions of direct and inverse Sturm–Liouville problems T1 - Lösungen von direkten und inversen Sturm-Liouville-Problemen N2 - Lie group method in combination with Magnus expansion is utilized to develop a universal method applicable to solving a Sturm–Liouville Problem (SLP) of any order with arbitrary boundary conditions. It is shown that the method has ability to solve direct regular and some singular SLPs of even orders (tested up to order eight), with a mix of boundary conditions (including non-separable and finite singular endpoints), accurately and efficiently. The present technique is successfully applied to overcome the difficulties in finding suitable sets of eigenvalues so that the inverse SLP problem can be effectively solved. Next, a concrete implementation to the inverse Sturm–Liouville problem algorithm proposed by Barcilon (1974) is provided. Furthermore, computational feasibility and applicability of this algorithm to solve inverse Sturm–Liouville problems of order n=2,4 is verified successfully. It is observed that the method is successful even in the presence of significant noise, provided that the assumptions of the algorithm are satisfied. In conclusion, this work provides methods that can be adapted successfully for solving a direct (regular/singular) or inverse SLP of an arbitrary order with arbitrary boundary conditions. N2 - Die Lie-Gruppen-Methode in Kombination mit der Magnus-Expansion wird verwendet, um eine universelle Methode zu entwickeln, die zur Lösung eines Sturm-Liouville-Problems (SLP) beliebiger Ordnung mit beliebigen Randbedingungen anwendbar ist. Es wird gezeigt, dass die Methode in der Lage ist, direkte reguläre und einige singuläre SLPs gerader Ordnung (getestet bis zur 8. Ordnung) mit einer Mischung von Randbedingungen (einschließlich nicht trennbarer und endlicher singulärer Endpunkte) genau und effizient zu lösen. Die vorliegende Technik wird erfolgreich angewendet, um die Schwierigkeiten beim Finden geeigneter Sätze von Eigenwerten zu überwinden, so dass das inverse SLP-Problem effektiv gelöst werden kann. Als nächstes wird eine konkrete Implementierung des von Barcilon (1974) vorgeschlagenen inversen Sturm-Liouville-Problemalgorithmus bereitgestellt. Weiterhin wird die rechnerische Durchführbarkeit und Anwendbarkeit dieses Algorithmus zur Lösung inverser Sturm-Liouville-Probleme der Ordnung n=2,4 erfolgreich verifiziert. Es wird beobachtet, dass das Verfahren selbst bei Vorhandensein von signifikantem Rauschen erfolgreich ist, vorausgesetzt, dass die Annahmen des Algorithmus erfüllt sind. Zusammenfassend stellt diese Arbeit Methoden zur Verfügung, die erfolgreich zur Lösung eines direkten (regulär/singulären) oder inversen SLP beliebiger Ordnung mit beliebigen Randbedingungen angepasst werden können. KW - Sturm-Liouville problem KW - Inverse Sturm-Liouville problem KW - Higher-order Sturm-Liouville problem KW - Sturm-Liouville-Problem höherer Ordnung KW - Inverses Sturm-Liouville-Problem KW - Sturm-Liouville-Problem Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-530064 ER - TY - JOUR A1 - Redmann, Martin A1 - Freitag, Melina A. T1 - Optimization based model order reduction for stochastic systems JF - Applied mathematics and computation N2 - In this paper, we bring together the worlds of model order reduction for stochastic linear systems and H-2-optimal model order reduction for deterministic systems. In particular, we supplement and complete the theory of error bounds for model order reduction of stochastic differential equations. With these error bounds, we establish a link between the output error for stochastic systems (with additive and multiplicative noise) and modified versions of the H-2-norm for both linear and bilinear deterministic systems. When deriving the respective optimality conditions for minimizing the error bounds, we see that model order reduction techniques related to iterative rational Krylov algorithms (IRKA) are very natural and effective methods for reducing the dimension of large-scale stochastic systems with additive and/or multiplicative noise. We apply modified versions of (linear and bilinear) IRKA to stochastic linear systems and show their efficiency in numerical experiments. KW - Model order reduction KW - Stochastic systems KW - Optimality conditions KW - Sylvester equations KW - Levy process Y1 - 2021 U6 - https://doi.org/10.1016/j.amc.2020.125783 SN - 0096-3003 SN - 1873-5649 VL - 398 PB - Elsevier CY - New York ER - TY - JOUR A1 - Reich, Sebastian A1 - Weissmann, Simon T1 - Fokker-Planck particle systems for Bayesian inference: computational approaches JF - SIAM ASA journal on uncertainty quantification N2 - Bayesian inference can be embedded into an appropriately defined dynamics in the space of probability measures. In this paper, we take Brownian motion and its associated Fokker-Planck equation as a starting point for such embeddings and explore several interacting particle approximations. More specifically, we consider both deterministic and stochastic interacting particle systems and combine them with the idea of preconditioning by the empirical covariance matrix. In addition to leading to affine invariant formulations which asymptotically speed up convergence, preconditioning allows for gradient-free implementations in the spirit of the ensemble Kalman filter. While such gradient-free implementations have been demonstrated to work well for posterior measures that are nearly Gaussian, we extend their scope of applicability to multimodal measures by introducing localized gradient-free approximations. Numerical results demonstrate the effectiveness of the considered methodologies. KW - Bayesian inverse problems KW - Fokker-Planck equation KW - gradient flow KW - affine KW - invariance KW - gradient-free sampling methods KW - localization Y1 - 2021 U6 - https://doi.org/10.1137/19M1303162 SN - 2166-2525 VL - 9 IS - 2 SP - 446 EP - 482 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - JOUR A1 - Rodríguez Zuluaga, Juan A1 - Stolle, Claudia A1 - Yamazaki, Yosuke A1 - Xiong, Chao A1 - England, Scott L. T1 - A synoptic-scale wavelike structure in the nighttime equatorial ionization anomaly JF - Earth and Space Science : ESS N2 - Both ground- and satellite-based airglow imaging have significantly contributed to understanding the low-latitude ionosphere, especially the morphology and dynamics of the equatorial ionization anomaly (EIA). The NASA Global-scale Observations of the Limb and Disk (GOLD) mission focuses on far-ultraviolet airglow images from a geostationary orbit at 47.5 degrees W. This region is of particular interest at low magnetic latitudes because of the high magnetic declination (i.e., about -20 degrees) and proximity of the South Atlantic magnetic anomaly. In this study, we characterize an exciting feature of the nighttime EIA using GOLD observations from October 5, 2018 to June 30, 2020. It consists of a wavelike structure of a few thousand kilometers seen as poleward and equatorward displacements of the EIA-crests. Initial analyses show that the synoptic-scale structure is symmetric about the dip equator and appears nearly stationary with time over the night. In quasi-dipole coordinates, maxima poleward displacements of the EIA-crests are seen at about +/- 12 degrees latitude and around 20 and 60 degrees longitude (i.e., in geographic longitude at the dip equator, about 53 degrees W and 14 degrees W). The wavelike structure presents typical zonal wavelengths of about 6.7 x 10(3) km and 3.3 x 10(3) km. The structure's occurrence and wavelength are highly variable on a day-to-day basis with no apparent dependence on geomagnetic activity. In addition, a cluster or quasi-periodic wave train of equatorial plasma depletions (EPDs) is often detected within the synoptic-scale structure. We further outline the difference in observing these EPDs from FUV images and in situ measurements during a GOLD and Swarm mission conjunction. KW - equatorial ionization anomaly KW - equatorial ionosphere KW - equatorial plasma bubbles KW - wave structure KW - forcing from below Y1 - 2021 U6 - https://doi.org/10.1029/2020EA001529 SN - 2333-5084 VL - 8 IS - 2 PB - American Geophysical Union CY - Malden, Mass. ER - TY - JOUR A1 - Roos, Saskia A1 - Otoba, Nobuhiko T1 - Scalar curvature and the multiconformal class of a direct product Riemannian manifold JF - Geometriae dedicata N2 - For a closed, connected direct product Riemannian manifold (M, g) = (M-1, g(1)) x ... x (M-l, g(l)), we define its multiconformal class [[g]] as the totality {integral(2)(1)g(1) circle plus center dot center dot center dot integral(2)(l)g(l)} of all Riemannian metrics obtained from multiplying the metric gi of each factor Mi by a positive function fi on the total space M. A multiconformal class [[ g]] contains not only all warped product type deformations of g but also the whole conformal class [(g) over tilde] of every (g) over tilde is an element of[[ g]]. In this article, we prove that [[g]] contains a metric of positive scalar curvature if and only if the conformal class of some factor (Mi, gi) does, under the technical assumption dim M-i = 2. We also show that, even in the case where every factor (M-i, g(i)) has positive scalar curvature, [[g]] contains a metric of scalar curvature constantly equal to -1 and with arbitrarily large volume, provided l = 2 and dim M = 3. KW - Positive scalar curvature KW - Constant scalar curvature KW - The Yamabe KW - problem KW - Warped product KW - Umbilic product KW - Twisted product Y1 - 2021 U6 - https://doi.org/10.1007/s10711-021-00636-9 SN - 0046-5755 SN - 1572-9168 VL - 214 IS - 1 SP - 801 EP - 829 PB - Springer CY - Dordrecht ER - TY - JOUR A1 - Ruchi, Sangeetika A1 - Dubinkina, Svetlana A1 - Wiljes, Jana de T1 - Fast hybrid tempered ensemble transform filter formulation for Bayesian elliptical problems via Sinkhorn approximation JF - Nonlinear processes in geophysics / European Geosciences Union ; American Geophysical Union N2 - Identification of unknown parameters on the basis of partial and noisy data is a challenging task, in particular in high dimensional and non-linear settings. Gaussian approximations to the problem, such as ensemble Kalman inversion, tend to be robust and computationally cheap and often produce astonishingly accurate estimations despite the simplifying underlying assumptions. Yet there is a lot of room for improvement, specifically regarding a correct approximation of a non-Gaussian posterior distribution. The tempered ensemble transform particle filter is an adaptive Sequential Monte Carlo (SMC) method, whereby resampling is based on optimal transport mapping. Unlike ensemble Kalman inversion, it does not require any assumptions regarding the posterior distribution and hence has shown to provide promising results for non-linear non-Gaussian inverse problems. However, the improved accuracy comes with the price of much higher computational complexity, and the method is not as robust as ensemble Kalman inversion in high dimensional problems. In this work, we add an entropy-inspired regularisation factor to the underlying optimal transport problem that allows the high computational cost to be considerably reduced via Sinkhorn iterations. Further, the robustness of the method is increased via an ensemble Kalman inversion proposal step before each update of the samples, which is also referred to as a hybrid approach. The promising performance of the introduced method is numerically verified by testing it on a steady-state single-phase Darcy flow model with two different permeability configurations. The results are compared to the output of ensemble Kalman inversion, and Markov chain Monte Carlo methods results are computed as a benchmark. Y1 - 2021 U6 - https://doi.org/10.5194/npg-28-23-2021 SN - 1023-5809 SN - 1607-7946 VL - 28 IS - 1 SP - 23 EP - 41 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Saynisch-Wagner, Jan A1 - Bärenzung, Julien A1 - Hornschild, Aaron A1 - Irrgang, Christopher A1 - Thomas, Maik T1 - Tide-induced magnetic signals and their errors derived from CHAMP and Swarm satellite magnetometer observations JF - Earth, planets and space : EPS N2 - Satellite-measured tidal magnetic signals are of growing importance. These fields are mainly used to infer Earth's mantle conductivity, but also to derive changes in the oceanic heat content. We present a new Kalman filter-based method to derive tidal magnetic fields from satellite magnetometers: KALMAG. The method's advantage is that it allows to study a precisely estimated posterior error covariance matrix. We present the results of a simultaneous estimation of the magnetic signals of 8 major tides from 17 years of Swarm and CHAMP data. For the first time, robustly derived posterior error distributions are reported along with the reported tidal magnetic fields. The results are compared to other estimates that are either based on numerical forward models or on satellite inversions of the same data. For all comparisons, maximal differences and the corresponding globally averaged RMSE are reported. We found that the inter-product differences are comparable with the KALMAG-based errors only in a global mean sense. Here, all approaches give values of the same order, e.g., 0.09 nT-0.14 nT for M2. Locally, the KALMAG posterior errors are up to one order smaller than the inter-product differences, e.g., 0.12 nT vs. 0.96 nT for M2. KW - Tides KW - Electromagnetic induction KW - Error covariance KW - Satellite magnetometer observations Y1 - 2021 U6 - https://doi.org/10.1186/s40623-021-01557-3 SN - 1880-5981 VL - 73 IS - 1 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Schanner, Maximilian Arthus A1 - Mauerberger, Stefan A1 - Korte, Monika A1 - Holschneider, Matthias T1 - Correlation based time evolution of the archeomagnetic field JF - Journal of geophysical research : JGR ; an international quarterly. B, Solid earth N2 - In a previous study, a new snapshot modeling concept for the archeomagnetic field was introduced (Mauerberger et al., 2020, ). By assuming a Gaussian process for the geomagnetic potential, a correlation-based algorithm was presented, which incorporates a closed-form spatial correlation function. This work extends the suggested modeling strategy to the temporal domain. A space-time correlation kernel is constructed from the tensor product of the closed-form spatial correlation kernel with a squared exponential kernel in time. Dating uncertainties are incorporated into the modeling concept using a noisy input Gaussian process. All but one modeling hyperparameters are marginalized, to reduce their influence on the outcome and to translate their variability to the posterior variance. The resulting distribution incorporates uncertainties related to dating, measurement and modeling process. Results from application to archeomagnetic data show less variation in the dipole than comparable models, but are in general agreement with previous findings. Y1 - 2021 U6 - https://doi.org/10.1029/2020JB021548 SN - 2169-9313 SN - 2169-9356 VL - 126 IS - 7 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Schick, Thomas A1 - Seyedhosseini, Mehran T1 - On an index theorem of Chang, Weinberger and Yu JF - Münster journal of mathematics N2 - In this paper we prove a strengthening of a theorem of Chang, Weinberger and Yu on obstructions to the existence of positive scalar curvature metrics on compact manifolds with boundary. They construct a relative index for the Dirac operator, which lives in a relative K-theory group, measuring the difference between the fundamental group of the boundary and of the full manifold. Whenever the Riemannian metric has product structure and positive scalar curvature near the boundary, one can define an absolute index of the Dirac operator taking value in the K-theory of the C*-algebra of fundamental group of the full manifold. This index depends on the metric near the boundary. We prove that (a slight variation of) the relative index of Chang, Weinberger and Yu is the image of this absolute index under the canonical map of K-theory groups. This has the immediate corollary that positive scalar curvature on the whole manifold implies vanishing of the relative index, giving a conceptual and direct proof of the vanishing theorem of Chang, Weinberger and Yu (rather: a slight variation). To take the fundamental groups of the manifold and its boundary into account requires working with maximal C*-completions of the involved *-algebras. A significant part of this paper is devoted to foundational results regarding these completions. On the other hand, we introduce and propose a more conceptual and more geometric completion, which still has all the required functoriality. Y1 - 2021 U6 - https://doi.org/10.17879/59019522628 SN - 1867-5778 SN - 1867-5786 VL - 14 IS - 1 SP - 123 EP - 154 PB - WWU, Fachbereich Mathematik und Informatik CY - Münster ER - TY - JOUR A1 - Schindler, Daniel A1 - Moldenhawer, Ted A1 - Stange, Maike A1 - Lepro, Valentino A1 - Beta, Carsten A1 - Holschneider, Matthias A1 - Huisinga, Wilhelm T1 - Analysis of protrusion dynamics in amoeboid cell motility by means of regularized contour flows JF - PLoS Computational Biology : a new community journal N2 - Amoeboid cell motility is essential for a wide range of biological processes including wound healing, embryonic morphogenesis, and cancer metastasis. It relies on complex dynamical patterns of cell shape changes that pose long-standing challenges to mathematical modeling and raise a need for automated and reproducible approaches to extract quantitative morphological features from image sequences. Here, we introduce a theoretical framework and a computational method for obtaining smooth representations of the spatiotemporal contour dynamics from stacks of segmented microscopy images. Based on a Gaussian process regression we propose a one-parameter family of regularized contour flows that allows us to continuously track reference points (virtual markers) between successive cell contours. We use this approach to define a coordinate system on the moving cell boundary and to represent different local geometric quantities in this frame of reference. In particular, we introduce the local marker dispersion as a measure to identify localized membrane expansions and provide a fully automated way to extract the properties of such expansions, including their area and growth time. The methods are available as an open-source software package called AmoePy, a Python-based toolbox for analyzing amoeboid cell motility (based on time-lapse microscopy data), including a graphical user interface and detailed documentation. Due to the mathematical rigor of our framework, we envision it to be of use for the development of novel cell motility models. We mainly use experimental data of the social amoeba Dictyostelium discoideum to illustrate and validate our approach.
Author summary Amoeboid motion is a crawling-like cell migration that plays an important key role in multiple biological processes such as wound healing and cancer metastasis. This type of cell motility results from expanding and simultaneously contracting parts of the cell membrane. From fluorescence images, we obtain a sequence of points, representing the cell membrane, for each time step. By using regression analysis on these sequences, we derive smooth representations, so-called contours, of the membrane. Since the number of measurements is discrete and often limited, the question is raised of how to link consecutive contours with each other. In this work, we present a novel mathematical framework in which these links are described by regularized flows allowing a certain degree of concentration or stretching of neighboring reference points on the same contour. This stretching rate, the so-called local dispersion, is used to identify expansions and contractions of the cell membrane providing a fully automated way of extracting properties of these cell shape changes. We applied our methods to time-lapse microscopy data of the social amoeba Dictyostelium discoideum. Y1 - 2021 U6 - https://doi.org/10.1371/journal.pcbi.1009268 SN - 1553-734X SN - 1553-7358 VL - 17 IS - 8 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Shlapunov, Alexander A1 - Tarchanov, Nikolaj Nikolaevič T1 - An open mapping theorem for the Navier-Stokes type equations associated with the de Rham complex over R-n JF - Siberian electronic mathematical reports = Sibirskie ėlektronnye matematičeskie izvestija N2 - We consider an initial problem for the Navier-Stokes type equations associated with the de Rham complex over R-n x[0, T], n >= 3, with a positive time T. We prove that the problem induces an open injective mappings on the scales of specially constructed function spaces of Bochner-Sobolev type. In particular, the corresponding statement on the intersection of these classes gives an open mapping theorem for smooth solutions to the Navier-Stokes equations. KW - Navier-Stokes equations KW - de Rham complex KW - open mapping theorem Y1 - 2021 U6 - https://doi.org/10.33048/semi.2021.18.108 SN - 1813-3304 VL - 18 IS - 2 SP - 1433 EP - 1466 PB - Institut Matematiki Imeni S. L. Soboleva CY - Novosibirsk ER - TY - JOUR A1 - Silverii, Francesca A1 - Maccaferri, Francesco A1 - Richter, Gudrun A1 - Gonzalez Cansado, Borja A1 - Wang, Rongjiang A1 - Hainzl, Sebastian A1 - Dahm, Torsten T1 - Poroelastic model in a vertically sealed gas storage BT - a case study from cyclic injection/production in a carbonate aquifer JF - Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society N2 - Natural gas can be temporarily stored in a variety of underground facilities, such as depleted gas and oil fields, natural aquifers and caverns in salt rocks. Being extensively monitored during operations, these systems provide a favourable opportunity to investigate how pressure varies in time and space and possibly induces/triggers earthquakes on nearby faults. Elaborate and detailed numerical modelling techniques are often applied to study gas reservoirs. Here we show the possibilities and discuss the limitations of a flexible and easily formulated tool that can be straightforwardly applied to simulate temporal pore-pressure variations and study the relation with recorded microseismic events. We use the software POEL (POroELastic diffusion and deformation) which computes the poroelastic response to fluid injection/extraction in a horizontally layered poroelastic structure. We further develop its application to address the presence of vertical impermeable faults bounding the reservoir and of multiple injection/extraction sources. Exploiting available information on the reservoir geometry and physical parameters, and records of injection/extraction rates for a gas reservoir in southern Europe, we perform an extensive parametric study considering different model configurations. Comparing modelled spatiotemporal pore-pressure variations with in situ measurements, we show that the inclusion of vertical impermeable faults provides an improvement in reproducing the observations and results in pore-pressure accumulation near the faults and in a variation of the temporal pore-pressure diffusion pattern. To study the relation between gas storage activity and recorded local microseismicity, we applied different seismicity models based on the estimated porepressure distribution. This analysis helps to understand the spatial distribution of seismicity and its temporal modulation. The results show that the observed microseismicity could be partly linked to the storage activity, but the contribution of tectonic background seismicity cannot be excluded. KW - Permeability and porosity KW - Gas and hydrate systems KW - Europe KW - Induced KW - seismicity Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab268 SN - 0956-540X SN - 1365-246X VL - 227 IS - 2 SP - 1322 EP - 1338 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Wormell, Caroline L. A1 - Reich, Sebastian T1 - Spectral convergence of diffusion maps BT - Improved error bounds and an alternative normalization JF - SIAM journal on numerical analysis / Society for Industrial and Applied Mathematics N2 - Diffusion maps is a manifold learning algorithm widely used for dimensionality reduction. Using a sample from a distribution, it approximates the eigenvalues and eigenfunctions of associated Laplace-Beltrami operators. Theoretical bounds on the approximation error are, however, generally much weaker than the rates that are seen in practice. This paper uses new approaches to improve the error bounds in the model case where the distribution is supported on a hypertorus. For the data sampling (variance) component of the error we make spatially localized compact embedding estimates on certain Hardy spaces; we study the deterministic (bias) component as a perturbation of the Laplace-Beltrami operator's associated PDE and apply relevant spectral stability results. Using these approaches, we match long-standing pointwise error bounds for both the spectral data and the norm convergence of the operator discretization. We also introduce an alternative normalization for diffusion maps based on Sinkhorn weights. This normalization approximates a Langevin diffusion on the sample and yields a symmetric operator approximation. We prove that it has better convergence compared with the standard normalization on flat domains, and we present a highly efficient rigorous algorithm to compute the Sinkhorn weights. KW - diffusion maps KW - graph Laplacian KW - Sinkhorn problem KW - kernel methods Y1 - 2021 U6 - https://doi.org/10.1137/20M1344093 SN - 0036-1429 SN - 1095-7170 VL - 59 IS - 3 SP - 1687 EP - 1734 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - THES A1 - Zadorozhnyi, Oleksandr T1 - Contributions to the theoretical analysis of the algorithms with adversarial and dependent data N2 - In this work I present the concentration inequalities of Bernstein's type for the norms of Banach-valued random sums under a general functional weak-dependency assumption (the so-called $\cC-$mixing). The latter is then used to prove, in the asymptotic framework, excess risk upper bounds of the regularised Hilbert valued statistical learning rules under the τ-mixing assumption on the underlying training sample. These results (of the batch statistical setting) are then supplemented with the regret analysis over the classes of Sobolev balls of the type of kernel ridge regression algorithm in the setting of online nonparametric regression with arbitrary data sequences. Here, in particular, a question of robustness of the kernel-based forecaster is investigated. Afterwards, in the framework of sequential learning, the multi-armed bandit problem under $\cC-$mixing assumption on the arm's outputs is considered and the complete regret analysis of a version of Improved UCB algorithm is given. Lastly, probabilistic inequalities of the first part are extended to the case of deviations (both of Azuma-Hoeffding's and of Burkholder's type) to the partial sums of real-valued weakly dependent random fields (under the type of projective dependence condition). KW - Machine learning KW - nonparametric regression KW - kernel methods KW - regularisation KW - concentration inequalities KW - learning rates KW - sequential learning KW - multi-armed bandits KW - Sobolev spaces Y1 - 2021 ER - TY - THES A1 - Zass, Alexander T1 - A multifaceted study of marked Gibbs point processes T1 - Facetten von markierten Gibbsschen Punktprozessen N2 - This thesis focuses on the study of marked Gibbs point processes, in particular presenting some results on their existence and uniqueness, with ideas and techniques drawn from different areas of statistical mechanics: the entropy method from large deviations theory, cluster expansion and the Kirkwood--Salsburg equations, the Dobrushin contraction principle and disagreement percolation. We first present an existence result for infinite-volume marked Gibbs point processes. More precisely, we use the so-called entropy method (and large-deviation tools) to construct marked Gibbs point processes in R^d under quite general assumptions. In particular, the random marks belong to a general normed space S and are not bounded. Moreover, we allow for interaction functionals that may be unbounded and whose range is finite but random. The entropy method relies on showing that a family of finite-volume Gibbs point processes belongs to sequentially compact entropy level sets, and is therefore tight. We then present infinite-dimensional Langevin diffusions, that we put in interaction via a Gibbsian description. In this setting, we are able to adapt the general result above to show the existence of the associated infinite-volume measure. We also study its correlation functions via cluster expansion techniques, and obtain the uniqueness of the Gibbs process for all inverse temperatures β and activities z below a certain threshold. This method relies in first showing that the correlation functions of the process satisfy a so-called Ruelle bound, and then using it to solve a fixed point problem in an appropriate Banach space. The uniqueness domain we obtain consists then of the model parameters z and β for which such a problem has exactly one solution. Finally, we explore further the question of uniqueness of infinite-volume Gibbs point processes on R^d, in the unmarked setting. We present, in the context of repulsive interactions with a hard-core component, a novel approach to uniqueness by applying the discrete Dobrushin criterion to the continuum framework. We first fix a discretisation parameter a>0 and then study the behaviour of the uniqueness domain as a goes to 0. With this technique we are able to obtain explicit thresholds for the parameters z and β, which we then compare to existing results coming from the different methods of cluster expansion and disagreement percolation. Throughout this thesis, we illustrate our theoretical results with various examples both from classical statistical mechanics and stochastic geometry. N2 - Diese Arbeit konzentriert sich auf die Untersuchung von markierten Gibbs-Punkt-Prozessen und stellt insbesondere einige Ergebnisse zu deren Existenz und Eindeutigkeit vor. Dabei werden Ideen und Techniken aus verschiedenen Bereichen der statistischen Mechanik verwendet: die Entropie-Methode aus der Theorie der großen Abweichungen, die Cluster-Expansion und die Kirkwood-Salsburg-Gleichungen, das Dobrushin-Kontraktionsprinzip und die Disagreement-Perkolation. Wir präsentieren zunächst ein Existenzergebnis für unendlich-volumige markierte Gibbs-Punkt-Prozesse. Genauer gesagt verwenden wir die sogenannte Entropie-Methode (und Werkzeuge der großen Abweichung), um markierte Gibbs-Punkt-Prozesse in R^d unter möglichst allgemeinen Annahmen zu konstruieren. Insbesondere gehören die zufälligen Markierungen zu einem allgemeinen normierten Raum und sind nicht beschränkt. Außerdem lassen wir Interaktionsfunktionale zu, die unbeschränkt sein können und deren Reichweite endlich, aber zufällig ist. Die Entropie-Methode beruht darauf, zu zeigen, dass eine Familie von endlich-volumigen Gibbs-Punkt-Prozessen zu sequentiell kompakten Entropie-Niveau-Mengen gehört, und daher dicht ist. Wir stellen dann unendlich-dimensionale Langevin-Diffusionen vor, die wir über eine Gibbssche Beschreibung in Wechselwirkung setzen. In dieser Umgebung sind wir in der Lage, das vorangehend vorgestellte allgemeine Ergebnis anzupassen, um die Existenz des zugehörigen unendlich-dimensionalen Maßes zu zeigen. Wir untersuchen auch seine Korrelationsfunktionen über Cluster-Expansions Techniken und erhalten die Eindeutigkeit des Gibbs-Prozesses für alle inversen Temperaturen β und Aktivitäten z unterhalb einer bestimmten Schwelle. Diese Methode beruht darauf, zunächst zu zeigen, dass die Korrelationsfunktionen des Prozesses eine so genannte Ruelle-Schranke erfüllen, um diese dann zur Lösung eines Fixpunktproblems in einem geeigneten Banach-Raum zu verwenden. Der Eindeutigkeitsbereich, den wir erhalten, wird dann aus den Modellparametern z und β definiert, für die ein solches Problem genau eine Lösung hat. Schließlich untersuchen wir die Frage nach der Eindeutigkeit von unendlich-volumigen Gibbs-Punkt-Prozessen auf R^d im unmarkierten Fall weiter. Im Zusammenhang mit repulsiven Wechselwirkungen basierend auf einer Hartkernkomponente stellen wir einen neuen Ansatz zur Eindeutigkeit vor, indem wir das diskrete Dobrushin-Kriterium im kontinuierlichen Rahmen anwenden. Wir legen zunächst einen Diskretisierungsparameter a>0 fest und untersuchen dann das Verhalten des Bereichs der Eindeutigkeit, wenn a gegen 0 geht. Mit dieser Technik sind wir in der Lage, explizite Schwellenwerte für die Parameter z und β zu erhalten, die wir dann mit bestehenden Ergebnissen aus den verschiedenen Methoden der Cluster-Expansion und der Disagreement-Perkolation vergleichen. In dieser Arbeit illustrieren wir unsere theoretischen Ergebnisse mit verschiedenen Beispielen sowohl aus der klassischen statistischen Mechanik als auch aus der stochastischen Geometrie. KW - marked Gibbs point processes KW - Langevin diffusions KW - Dobrushin criterion KW - Entropy method KW - Cluster expansion KW - Kirkwood--Salsburg equations KW - DLR equations KW - Markierte Gibbs-Punkt-Prozesse KW - Entropiemethode KW - Cluster-Expansion KW - DLR-Gleichungen KW - Dobrushin-Kriterium KW - Kirkwood-Salsburg-Gleichungen KW - Langevin-Diffusions Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-512775 ER -