@article{EshghiMachReichel2021, author = {Eshghi, Nasim and Mach, Thomas and Reichel, Lothar}, title = {New matrix function approximations and quadrature rules based on the Arnoldi process}, series = {Journal of computational and applied mathematics}, volume = {391}, journal = {Journal of computational and applied mathematics}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0377-0427}, doi = {10.1016/j.cam.2021.113442}, pages = {12}, year = {2021}, abstract = {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.}, language = {en} } @article{PohleAdamBeumer2022, author = {Pohle, Jennifer and Adam, Timo and Beumer, Larissa}, title = {Flexible estimation of the state dwell-time distribution in hidden semi-Markov models}, series = {Computational statistics \& data analysis}, volume = {172}, journal = {Computational statistics \& data analysis}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0167-9473}, doi = {10.1016/j.csda.2022.107479}, pages = {15}, year = {2022}, abstract = {Hidden semi-Markov models generalise hidden Markov models by explicitly modelling the time spent in a given state, the so-called dwell time, using some distribution defined on the natural numbers. While the (shifted) Poisson and negative binomial distribution provide natural choices for such distributions, in practice, parametric distributions can lack the flexibility to adequately model the dwell times. To overcome this problem, a penalised maximum likelihood approach is proposed that allows for a flexible and data-driven estimation of the dwell-time distributions without the need to make any distributional assumption. This approach is suitable for direct modelling purposes or as an exploratory tool to investigate the latent state dynamics. The feasibility and potential of the suggested approach is illustrated in a simulation study and by modelling muskox movements in northeast Greenland using GPS tracking data. The proposed method is implemented in the R-package PHSMM which is available on CRAN.}, language = {en} } @misc{KolbeEvans2020, author = {Kolbe, Benedikt Maximilian and Evans, Myfanwy E.}, title = {Isotopic tiling theory for hyperbolic surfaces}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1}, issn = {1866-8372}, doi = {10.25932/publishup-54428}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-544285}, pages = {30}, year = {2020}, abstract = {In this paper, we develop the mathematical tools needed to explore isotopy classes of tilings on hyperbolic surfaces of finite genus, possibly nonorientable, with boundary, and punctured. More specifically, we generalize results on Delaney-Dress combinatorial tiling theory using an extension of mapping class groups to orbifolds, in turn using this to study tilings of covering spaces of orbifolds. Moreover, we study finite subgroups of these mapping class groups. Our results can be used to extend the Delaney-Dress combinatorial encoding of a tiling to yield a finite symbol encoding the complexity of an isotopy class of tilings. The results of this paper provide the basis for a complete and unambiguous enumeration of isotopically distinct tilings of hyperbolic surfaces.}, language = {en} } @article{PathirajaLeeuwen2022, author = {Pathiraja, Sahani Darschika and Leeuwen, Peter Jan van}, title = {Multiplicative Non-Gaussian model error estimation in data assimilation}, series = {Journal of advances in modeling earth systems : JAMES}, volume = {14}, journal = {Journal of advances in modeling earth systems : JAMES}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1942-2466}, doi = {10.1029/2021MS002564}, pages = {23}, year = {2022}, abstract = {Model uncertainty quantification is an essential component of effective data assimilation. Model errors associated with sub-grid scale processes are often represented through stochastic parameterizations of the unresolved process. Many existing Stochastic Parameterization schemes are only applicable when knowledge of the true sub-grid scale process or full observations of the coarse scale process are available, which is typically not the case in real applications. We present a methodology for estimating the statistics of sub-grid scale processes for the more realistic case that only partial observations of the coarse scale process are available. Model error realizations are estimated over a training period by minimizing their conditional sum of squared deviations given some informative covariates (e.g., state of the system), constrained by available observations and assuming that the observation errors are smaller than the model errors. From these realizations a conditional probability distribution of additive model errors given these covariates is obtained, allowing for complex non-Gaussian error structures. Random draws from this density are then used in actual ensemble data assimilation experiments. We demonstrate the efficacy of the approach through numerical experiments with the multi-scale Lorenz 96 system using both small and large time scale separations between slow (coarse scale) and fast (fine scale) variables. The resulting error estimates and forecasts obtained with this new method are superior to those from two existing methods.}, language = {en} } @misc{KellerPinchoverPogorzelski2020, author = {Keller, Matthias and Pinchover, Yehuda and Pogorzelski, Felix}, title = {From hardy to rellich inequalities on graphs}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {3}, issn = {1866-8372}, doi = {10.25932/publishup-54214}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-542140}, pages = {22}, year = {2020}, abstract = {We show how to deduce Rellich inequalities from Hardy inequalities on infinite graphs. Specifically, the obtained Rellich inequality gives an upper bound on a function by the Laplacian of the function in terms of weighted norms. These weights involve the Hardy weight and a function which satisfies an eikonal inequality. The results are proven first for Laplacians and are extended to Schrodinger operators afterwards.}, language = {en} } @misc{MazzonettoSalimova2020, author = {Mazzonetto, Sara and Salimova, Diyora}, title = {Existence, uniqueness, and numerical approximations for stochastic burgers equations}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {4}, issn = {1866-8372}, doi = {10.25932/publishup-51579}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-515796}, pages = {26}, year = {2020}, abstract = {In this article, we propose an all-in-one statement which includes existence, uniqueness, regularity, and numerical approximations of mild solutions for a class of stochastic partial differential equations (SPDEs) with non-globally monotone nonlinearities. The proof of this result exploits the properties of an existing fully explicit space-time discrete approximation scheme, in particular the fact that it satisfies suitable a priori estimates. We also obtain almost sure and strong convergence of the approximation scheme to the mild solutions of the considered SPDEs. We conclude by applying the main result of the article to the stochastic Burgers equations with additive space-time white noise.}, language = {en} } @article{HetheyHartungWangorschetal.2021, author = {Hethey, Christoph Philipp and Hartung, Niklas and Wangorsch, Gaby and Weisser, Karin and Huisinga, Wilhelm}, title = {Physiology-based toxicokinetic modelling of aluminium in rat and man}, series = {Archives of toxicology : official journal of EUROTOX}, volume = {95}, journal = {Archives of toxicology : official journal of EUROTOX}, number = {9}, publisher = {Springer}, address = {Berlin ; Heidelberg}, issn = {0340-5761}, doi = {10.1007/s00204-021-03107-y}, pages = {2977 -- 3000}, year = {2021}, abstract = {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.}, language = {en} } @article{LilienkampvonSpechtWeatherilletal.2022, author = {Lilienkamp, Henning and von Specht, Sebastian and Weatherill, Graeme and Caire, Giuseppe and Cotton, Fabrice}, title = {Ground-Motion modeling as an image processing task}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120220008}, pages = {1565 -- 1582}, year = {2022}, abstract = {We construct and examine the prototype of a deep learning-based ground-motion model (GMM) that is both fully data driven and nonergodic. We formulate ground-motion modeling as an image processing task, in which a specific type of neural network, the U-Net, relates continuous, horizontal maps of earthquake predictive parameters to sparse observations of a ground-motion intensity measure (IM). The processing of map-shaped data allows the natural incorporation of absolute earthquake source and observation site coordinates, and is, therefore, well suited to include site-, source-, and path-specific amplification effects in a nonergodic GMM. Data-driven interpolation of the IM between observation points is an inherent feature of the U-Net and requires no a priori assumptions. We evaluate our model using both a synthetic dataset and a subset of observations from the KiK-net strong motion network in the Kanto basin in Japan. We find that the U-Net model is capable of learning the magnitude???distance scaling, as well as site-, source-, and path-specific amplification effects from a strong motion dataset. The interpolation scheme is evaluated using a fivefold cross validation and is found to provide on average unbiased predictions. The magnitude???distance scaling as well as the site amplification of response spectral acceleration at a period of 1 s obtained for the Kanto basin are comparable to previous regional studies.}, language = {en} } @article{CvetkovićConradLie2021, author = {Cvetković, Nada and Conrad, Tim and Lie, Han Cheng}, title = {A convergent discretization method for transition path theory for diffusion processes}, series = {Multiscale modeling \& simulation : a SIAM interdisciplinary journal}, volume = {19}, journal = {Multiscale modeling \& simulation : a SIAM interdisciplinary journal}, number = {1}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {1540-3459}, doi = {10.1137/20M1329354}, pages = {242 -- 266}, year = {2021}, abstract = {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.}, language = {en} } @article{DubeBoeckmannRitter2022, author = {Dube, Jonas and B{\"o}ckmann, Christine and Ritter, Christoph}, title = {Lidar-Derived Aerosol Properties from Ny-{\AA}lesund, Svalbard during the MOSAiC Spring 2020}, series = {Remote sensing / Molecular Diversity Preservation International (MDPI)}, volume = {14}, journal = {Remote sensing / Molecular Diversity Preservation International (MDPI)}, number = {11}, publisher = {MDPI}, address = {Basel}, issn = {2072-4292}, doi = {10.3390/rs14112578}, pages = {17}, year = {2022}, abstract = {In this work, we present Raman lidar data (from a Nd:YAG operating at 355 nm, 532 nm and 1064 nm) from the international research village Ny-Alesund for the time period of January to April 2020 during the Arctic haze season of the MOSAiC winter. We present values of the aerosol backscatter, the lidar ratio and the backscatter Angstrom exponent, though the latter depends on wavelength. The aerosol polarization was generally below 2\%, indicating mostly spherical particles. We observed that events with high backscatter and high lidar ratio did not coincide. In fact, the highest lidar ratios (LR > 75 sr at 532 nm) were already found by January and may have been caused by hygroscopic growth, rather than by advection of more continental aerosol. Further, we performed an inversion of the lidar data to retrieve a refractive index and a size distribution of the aerosol. Our results suggest that in the free troposphere (above approximate to 2500 m) the aerosol size distribution is quite constant in time, with dominance of small particles with a modal radius well below 100 nm. On the contrary, below approximate to 2000 m in altitude, we frequently found gradients in aerosol backscatter and even size distribution, sometimes in accordance with gradients of wind speed, humidity or elevated temperature inversions, as if the aerosol was strongly modified by vertical displacement in what we call the "mechanical boundary layer". Finally, we present an indication that additional meteorological soundings during MOSAiC campaign did not necessarily improve the fidelity of air backtrajectories.}, language = {en} } @article{KemptonMuenchYau2021, author = {Kempton, Mark and M{\"u}nch, Florentin and Yau, Shing-Tung}, title = {A homology vanishing theorem for graphs with positive curvature}, series = {Communications in analysis and geometry}, volume = {29}, journal = {Communications in analysis and geometry}, number = {6}, publisher = {International Press of Boston}, address = {Somerville}, issn = {1019-8385}, doi = {10.4310/CAG.2021.v29.n6.a5}, pages = {1449 -- 1473}, year = {2021}, abstract = {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.}, language = {en} } @article{DimitrovaKoppitz2020, author = {Dimitrova, Ilinka and Koppitz, J{\"o}rg}, title = {On relative ranks of the semigroup of orientation-preserving transformations on infinite chains}, series = {Asian-European journal of mathematics}, volume = {14}, journal = {Asian-European journal of mathematics}, number = {08}, publisher = {World Scientific}, address = {Singapore}, issn = {1793-5571}, doi = {10.1142/S1793557121501461}, pages = {15}, year = {2020}, abstract = {In this paper, we determine the relative rank of the semigroup OP(X) of all orientation-preserving transformations on infinite chains modulo the semigroup O(X) of all order-preserving transformations.}, language = {en} } @article{MalemShinitskiOjedaOpper2022, author = {Malem-Shinitski, Noa and Ojeda, Cesar and Opper, Manfred}, title = {Variational bayesian inference for nonlinear hawkes process with gaussian process self-effects}, series = {Entropy}, volume = {24}, journal = {Entropy}, number = {3}, publisher = {MDPI}, address = {Basel}, issn = {1099-4300}, doi = {10.3390/e24030356}, pages = {22}, year = {2022}, abstract = {Traditionally, Hawkes processes are used to model time-continuous point processes with history dependence. Here, we propose an extended model where the self-effects are of both excitatory and inhibitory types and follow a Gaussian Process. Whereas previous work either relies on a less flexible parameterization of the model, or requires a large amount of data, our formulation allows for both a flexible model and learning when data are scarce. We continue the line of work of Bayesian inference for Hawkes processes, and derive an inference algorithm by performing inference on an aggregated sum of Gaussian Processes. Approximate Bayesian inference is achieved via data augmentation, and we describe a mean-field variational inference approach to learn the model parameters. To demonstrate the flexibility of the model we apply our methodology on data from different domains and compare it to previously reported results.}, language = {en} } @article{AyanbayevKlebanovLieetal.2021, author = {Ayanbayev, Birzhan and Klebanov, Ilja and Lie, Han Cheng and Sullivan, Tim J.}, title = {Gamma-convergence of Onsager-Machlup functionals}, series = {Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data}, volume = {38}, journal = {Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data}, number = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0266-5611}, doi = {10.1088/1361-6420/ac3f82}, pages = {35}, year = {2021}, abstract = {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.}, language = {en} } @article{BellingeriFrizPaychaetal.2022, author = {Bellingeri, Carlo and Friz, Peter and Paycha, Sylvie and Preiß, Rosa Lili Dora}, title = {Smooth rough paths, their geometry and algebraic renormalization}, series = {Vietnam journal of mathematics}, volume = {50}, journal = {Vietnam journal of mathematics}, number = {3}, publisher = {Springer}, address = {Singapore}, issn = {2305-221X}, doi = {10.1007/s10013-022-00570-7}, pages = {719 -- 761}, year = {2022}, abstract = {We introduce the class of "smooth rough paths" and study their main properties. Working in a smooth setting allows us to discard sewing arguments and focus on algebraic and geometric aspects. Specifically, a Maurer-Cartan perspective is the key to a purely algebraic form of Lyons' extension theorem, the renormalization of rough paths following up on [Bruned et al.: A rough path perspective on renormalization, J. Funct. Anal. 277(11), 2019], as well as a related notion of "sum of rough paths". We first develop our ideas in a geometric rough path setting, as this best resonates with recent works on signature varieties, as well as with the renormalization of geometric rough paths. We then explore extensions to the quasi-geometric and the more general Hopf algebraic setting.}, language = {en} } @article{AyanbayevKlebanovLietal.2021, author = {Ayanbayev, Birzhan and Klebanov, Ilja and Li, Han Cheng and Sullivan, Tim J.}, title = {Gamma-convergence of Onsager-Machlup functionals}, series = {Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data}, volume = {38}, journal = {Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data}, number = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0266-5611}, doi = {10.1088/1361-6420/ac3f81}, pages = {32}, year = {2021}, abstract = {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.}, language = {en} } @article{DimitrovaKoppitz2022, author = {Dimitrova, Ilinka and Koppitz, J{\"o}rg}, title = {On relative ranks of the semigroup of orientation-preserving transformations on infinite chain with restricted range}, series = {Communications in algebra}, volume = {50}, journal = {Communications in algebra}, number = {5}, publisher = {Taylor \& Francis Group}, address = {Philadelphia}, issn = {0092-7872}, doi = {10.1080/00927872.2021.2000998}, pages = {2157 -- 2168}, year = {2022}, abstract = {Let X be an infinite linearly ordered set and let Y be a nonempty subset of X. We calculate the relative rank of the semigroup OP(X,Y) of all orientation-preserving transformations on X with restricted range Y modulo the semigroup O(X,Y) of all order-preserving transformations on X with restricted range Y. For Y = X, we characterize the relative generating sets of minimal size.}, language = {en} } @article{LeungLeutbecherReichetal.2021, author = {Leung, Tsz Yan and Leutbecher, Martin and Reich, Sebastian and Shepherd, Theodore G.}, title = {Forecast verification}, series = {Quarterly journal of the Royal Meteorological Society}, volume = {147}, journal = {Quarterly journal of the Royal Meteorological Society}, number = {739}, publisher = {Wiley}, address = {Hoboken}, issn = {0035-9009}, doi = {10.1002/qj.4120}, pages = {3124 -- 3134}, year = {2021}, abstract = {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.}, language = {en} } @inproceedings{RungrottheeraChangSchulze2020, author = {Rungrottheera, Wannarut and Chang, Der-Chen and Schulze, Bert-Wolfgang}, title = {The edge calculus of singularity order >3}, series = {Journal of nonlinear and convex analysis : an international journal}, volume = {21}, booktitle = {Journal of nonlinear and convex analysis : an international journal}, number = {2}, publisher = {Yokohama Publishers}, address = {Yokohama}, issn = {1345-4773}, pages = {387 -- 401}, year = {2020}, abstract = {We study Mellin pseudo-differential algebras on singular straight cones and manifolds with singularity of order >= 3. Those are necessary to express parametrices of elliptic differential operators with a corresponding cornerdegenerate behavior, and we obtain regularity in weighted spaces.}, language = {en} } @article{PornsawadSungcharoenBoeckmann2020, author = {Pornsawad, Pornsarp and Sungcharoen, Parada and B{\"o}ckmann, Christine}, title = {Convergence rate of the modified Landweber method for solving inverse potential problems}, series = {Mathematics : open access journal}, volume = {8}, journal = {Mathematics : open access journal}, number = {4}, publisher = {MDPI}, address = {Basel}, issn = {2227-7390}, doi = {10.3390/math8040608}, pages = {22}, year = {2020}, abstract = {In this paper, we present the convergence rate analysis of the modified Landweber method under logarithmic source condition for nonlinear ill-posed problems. The regularization parameter is chosen according to the discrepancy principle. The reconstructions of the shape of an unknown domain for an inverse potential problem by using the modified Landweber method are exhibited.}, language = {en} } @article{KuerschnerFreitag2020, author = {K{\"u}rschner, Patrick and Freitag, Melina A.}, title = {Inexact methods for the low rank solution to large scale Lyapunov equations}, series = {BIT : numerical mathematics ; the leading applied mathematics journal for all computational mathematicians}, volume = {60}, journal = {BIT : numerical mathematics ; the leading applied mathematics journal for all computational mathematicians}, number = {4}, publisher = {Springer}, address = {Dordrecht}, issn = {0006-3835}, doi = {10.1007/s10543-020-00813-4}, pages = {1221 -- 1259}, year = {2020}, abstract = {The rational Krylov subspace method (RKSM) and the low-rank alternating directions implicit (LR-ADI) iteration are established numerical tools for computing low-rank solution factors of large-scale Lyapunov equations. In order to generate the basis vectors for the RKSM, or extend the low-rank factors within the LR-ADI method, the repeated solution to a shifted linear system of equations is necessary. For very large systems this solve is usually implemented using iterative methods, leading to inexact solves within this inner iteration (and therefore to "inexact methods"). We will show that one can terminate this inner iteration before full precision has been reached and still obtain very good accuracy in the final solution to the Lyapunov equation. In particular, for both the RKSM and the LR-ADI method we derive theory for a relaxation strategy (e.g. increasing the solve tolerance of the inner iteration, as the outer iteration proceeds) within the iterative methods for solving the large linear systems. These theoretical choices involve unknown quantities, therefore practical criteria for relaxing the solution tolerance within the inner linear system are then provided. The theory is supported by several numerical examples, which show that the total amount of work for solving Lyapunov equations can be reduced significantly.}, language = {en} } @article{HartungWahlRastogietal.2021, author = {Hartung, Niklas and Wahl, Martin and Rastogi, Abhishake and Huisinga, Wilhelm}, title = {Nonparametric goodness-of-fit testing for parametric covariate models in pharmacometric analyses}, series = {CPT: pharmacometrics \& systems pharmacology}, volume = {10}, journal = {CPT: pharmacometrics \& systems pharmacology}, number = {6}, publisher = {Nature Publ. Group}, address = {London}, issn = {2163-8306}, doi = {10.1002/psp4.12614}, pages = {564 -- 576}, year = {2021}, abstract = {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.}, language = {en} } @article{JiaAnslanChenetal.2022, author = {Jia, Weihan and Anslan, Sten and Chen, Fahu and Cao, Xianyong and Dong, Hailiang and Dulias, Katharina and Gu, Zhengquan and Heinecke, Liv and Jiang, Hongchen and Kruse, Stefan and Kang, Wengang and Li, Kai and Liu, Sisi and Liu, Xingqi and Liu, Ying and Ni, Jian and Schwalb, Antje and Stoof-Leichsenring, Kathleen R. and Shen, Wei and Tian, Fang and Wang, Jing and Wang, Yongbo and Wang, Yucheng and Xu, Hai and Yang, Xiaoyan and Zhang, Dongju and Herzschuh, Ulrike}, title = {Sedimentary ancient DNA reveals past ecosystem and biodiversity changes on the Tibetan Plateau: overview and prospects}, series = {Quaternary science reviews : the international multidisciplinary research and review journal}, volume = {293}, journal = {Quaternary science reviews : the international multidisciplinary research and review journal}, publisher = {Elsevier}, address = {Oxford}, issn = {0277-3791}, doi = {10.1016/j.quascirev.2022.107703}, pages = {14}, year = {2022}, abstract = {Alpine ecosystems on the Tibetan Plateau are being threatened by ongoing climate warming and intensified human activities. Ecological time-series obtained from sedimentary ancient DNA (sedaDNA) are essential for understanding past ecosystem and biodiversity dynamics on the Tibetan Plateau and their responses to climate change at a high taxonomic resolution. Hitherto only few but promising studies have been published on this topic. The potential and limitations of using sedaDNA on the Tibetan Plateau are not fully understood. Here, we (i) provide updated knowledge of and a brief introduction to the suitable archives, region-specific taphonomy, state-of-the-art methodologies, and research questions of sedaDNA on the Tibetan Plateau; (ii) review published and ongoing sedaDNA studies from the Tibetan Plateau; and (iii) give some recommendations for future sedaDNA study designs. Based on the current knowledge of taphonomy, we infer that deep glacial lakes with freshwater and high clay sediment input, such as those from the southern and southeastern Tibetan Plateau, may have a high potential for sedaDNA studies. Metabarcoding (for microorganisms and plants), metagenomics (for ecosystems), and hybridization capture (for prehistoric humans) are three primary sedaDNA approaches which have been successfully applied on the Tibetan Plateau, but their power is still limited by several technical issues, such as PCR bias and incompleteness of taxonomic reference databases. Setting up high-quality and open-access regional taxonomic reference databases for the Tibetan Plateau should be given priority in the future. To conclude, the archival, taphonomic, and methodological conditions of the Tibetan Plateau are favorable for performing sedaDNA studies. More research should be encouraged to address questions about long-term ecological dynamics at ecosystem scale and to bring the paleoecology of the Tibetan Plateau into a new era.}, language = {en} } @article{SchannerKorteHolschneider2022, author = {Schanner, Maximilian and Korte, Monika and Holschneider, Matthias}, title = {ArchKalmag14k: A kalman-filter based global geomagnetic model for the holocene}, series = {Journal of geophysical research : Solid earth}, volume = {127}, journal = {Journal of geophysical research : Solid earth}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2021JB023166}, pages = {17}, year = {2022}, abstract = {We propose a global geomagnetic field model for the last 14 thousand years, based on thermoremanent records. We call the model ArchKalmag14k. ArchKalmag14k is constructed by modifying recently proposed algorithms, based on space-time correlations. Due to the amount of data and complexity of the model, the full Bayesian posterior is numerically intractable. To tackle this, we sequentialize the inversion by implementing a Kalman-filter with a fixed time step. Every step consists of a prediction, based on a degree dependent temporal covariance, and a correction via Gaussian process regression. Dating errors are treated via a noisy input formulation. Cross correlations are reintroduced by a smoothing algorithm and model parameters are inferred from the data. Due to the specific statistical nature of the proposed algorithms, the model comes with space and time-dependent uncertainty estimates. The new model ArchKalmag14k shows less variation in the large-scale degrees than comparable models. Local predictions represent the underlying data and agree with comparable models, if the location is sampled well. Uncertainties are bigger for earlier times and in regions of sparse data coverage. We also use ArchKalmag14k to analyze the appearance and evolution of the South Atlantic anomaly together with reverse flux patches at the core-mantle boundary, considering the model uncertainties. While we find good agreement with earlier models for recent times, our model suggests a different evolution of intensity minima prior to 1650 CE. In general, our results suggest that prior to 6000 BCE the data is not sufficient to support global models.}, language = {en} } @article{HanischLudewig2022, author = {Hanisch, Florian and Ludewig, Matthias}, title = {A rigorous construction of the supersymmetric path integral associated to a compact spin manifold}, series = {Communications in mathematical physics}, volume = {391}, journal = {Communications in mathematical physics}, number = {3}, publisher = {Springer}, address = {Berlin ; Heidelberg}, issn = {0010-3616}, doi = {10.1007/s00220-022-04336-7}, pages = {1209 -- 1239}, year = {2022}, abstract = {We give a rigorous construction of the path integral in N = 1/2 supersymmetry as an integral map for differential forms on the loop space of a compact spin manifold. It is defined on the space of differential forms which can be represented by extended iterated integrals in the sense of Chen and Getzler-Jones-Petrack. Via the iterated integral map, we compare our path integral to the non-commutative loop space Chern character of Guneysu and the second author. Our theory provides a rigorous background to various formal proofs of the Atiyah-Singer index theorem for twisted Dirac operators using supersymmetric path integrals, as investigated by Alvarez-Gaume, Atiyah, Bismut and Witten.}, language = {en} } @article{DeOliveiraGomesHoegele2021, author = {De Oliveira Gomes, Andr{\´e} and H{\"o}gele, Michael Anton}, title = {The Kramers problem for SDEs driven by small, accelerated L{\´e}vy noise with exponentially light jumps}, series = {Stochastics and dynamics}, volume = {21}, journal = {Stochastics and dynamics}, number = {04}, publisher = {World Scientific}, address = {Singapore}, issn = {0219-4937}, doi = {10.1142/S0219493721500192}, pages = {44}, year = {2021}, abstract = {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.}, language = {en} } @article{LieStahnSullivan2022, author = {Lie, Han Cheng and Stahn, Martin and Sullivan, Tim J.}, title = {Randomised one-step time integration methods for deterministic operator differential equations}, series = {Calcolo}, volume = {59}, journal = {Calcolo}, number = {1}, publisher = {Springer}, address = {Milano}, issn = {0008-0624}, doi = {10.1007/s10092-022-00457-6}, pages = {33}, year = {2022}, abstract = {Uncertainty quantification plays an important role in problems that involve inferring a parameter of an initial value problem from observations of the solution. Conrad et al. (Stat Comput 27(4):1065-1082, 2017) proposed randomisation of deterministic time integration methods as a strategy for quantifying uncertainty due to the unknown time discretisation error. We consider this strategy for systems that are described by deterministic, possibly time-dependent operator differential equations defined on a Banach space or a Gelfand triple. Our main results are strong error bounds on the random trajectories measured in Orlicz norms, proven under a weaker assumption on the local truncation error of the underlying deterministic time integration method. Our analysis establishes the theoretical validity of randomised time integration for differential equations in infinite-dimensional settings.}, language = {en} } @article{WiljesTong2020, author = {Wiljes, Jana de and Tong, Xin T.}, title = {Analysis of a localised nonlinear ensemble Kalman Bucy filter with complete and accurate observations}, series = {Nonlinearity}, volume = {33}, journal = {Nonlinearity}, number = {9}, publisher = {IOP Publ.}, address = {Bristol}, issn = {0951-7715}, doi = {10.1088/1361-6544/ab8d14}, pages = {4752 -- 4782}, year = {2020}, abstract = {Concurrent observation technologies have made high-precision real-time data available in large quantities. Data assimilation (DA) is concerned with how to combine this data with physical models to produce accurate predictions. For spatial-temporal models, the ensemble Kalman filter with proper localisation techniques is considered to be a state-of-the-art DA methodology. This article proposes and investigates a localised ensemble Kalman Bucy filter for nonlinear models with short-range interactions. We derive dimension-independent and component-wise error bounds and show the long time path-wise error only has logarithmic dependence on the time range. The theoretical results are verified through some simple numerical tests.}, language = {en} } @article{SchickSeyedhosseini2021, author = {Schick, Thomas and Seyedhosseini, Mehran}, title = {On an index theorem of Chang, Weinberger and Yu}, series = {M{\"u}nster journal of mathematics}, volume = {14}, journal = {M{\"u}nster journal of mathematics}, number = {1}, publisher = {WWU, Fachbereich Mathematik und Informatik}, address = {M{\"u}nster}, issn = {1867-5778}, doi = {10.17879/59019522628}, pages = {123 -- 154}, year = {2021}, abstract = {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.}, language = {en} } @article{DereudreHoudebert2018, author = {Dereudre, David and Houdebert, Pierre}, title = {Sharp phase transition for the continuum Widom-Rowlinson model}, series = {Annales de l'Institut Henri Poincar{\´e}. B, Probability and statistics}, volume = {57}, journal = {Annales de l'Institut Henri Poincar{\´e}. B, Probability and statistics}, number = {1}, publisher = {Association des Publications de l'Institut Henri Poincar{\´e}}, address = {Bethesda, Md.}, issn = {0246-0203}, doi = {10.1214/20-AIHP1082}, pages = {387 -- 407}, year = {2018}, abstract = {The Widom-Rowlinson model (or the Area-interaction model) is a Gibbs point process in R-d with the formal Hamiltonian defined as the volume of Ux epsilon omega B1(x), where. is a locally finite configuration of points and B-1(x) denotes the unit closed ball centred at x. The model is also tuned by two other parameters: the activity z > 0 related to the intensity of the process and the inverse temperature beta >= 0 related to the strength of the interaction. In the present paper we investigate the phase transition of the model in the point of view of percolation theory and the liquid-gas transition. First, considering the graph connecting points with distance smaller than 2r > 0, we show that for any beta >= 0, there exists 0 <(similar to a)(zc) (beta, r) < +infinity such that an exponential decay of connectivity at distance n occurs in the subcritical phase (i.e. z <(similar to a)(zc) (beta, r)) and a linear lower bound of the connection at infinity holds in the supercritical case (i.e. z >(similar to a)(zc) (beta, r)). These results are in the spirit of recent works using the theory of randomised tree algorithms (Probab. Theory Related Fields 173 (2019) 479-490, Ann. of Math. 189 (2019) 75-99, Duminil-Copin, Raoufi and Tassion (2018)). Secondly we study a standard liquid-gas phase transition related to the uniqueness/non-uniqueness of Gibbs states depending on the parameters z, beta. Old results (Phys. Rev. Lett. 27 (1971) 1040-1041, J. Chem. Phys. 52 (1970) 1670-1684) claim that a non-uniqueness regime occurs for z = beta large enough and it is conjectured that the uniqueness should hold outside such an half line ( z = beta >= beta(c) > 0). We solve partially this conjecture in any dimension by showing that for beta large enough the non-uniqueness holds if and only if z = beta. We show also that this critical value z = beta corresponds to the percolation threshold (similar to a)(zc) (beta, r) = beta for beta large enough, providing a straight connection between these two notions of phase transition.}, language = {en} } @article{CozzoniMaibaumHamm2021, author = {Cozzoni, Barbara and Maibaum, Michael and Hamm, Maximilian}, title = {Thermal analysis and constraints for the MASCOT landing site selection on the asteroid Ryugu}, series = {Planetary and space science}, volume = {205}, journal = {Planetary and space science}, publisher = {Elsevier}, address = {Oxford}, issn = {0032-0633}, doi = {10.1016/j.pss.2021.105286}, pages = {11}, year = {2021}, abstract = {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.}, language = {en} } @article{HinzSchwarz2022, author = {Hinz, Michael and Schwarz, Michael}, title = {A note on Neumann problems on graphs}, series = {Positivity}, volume = {26}, journal = {Positivity}, number = {4}, publisher = {Springer}, address = {Dordrecht}, issn = {1385-1292}, doi = {10.1007/s11117-022-00930-0}, pages = {23}, year = {2022}, abstract = {We discuss Neumann problems for self-adjoint Laplacians on (possibly infinite) graphs. Under the assumption that the heat semigroup is ultracontractive we discuss the unique solvability for non-empty subgraphs with respect to the vertex boundary and provide analytic and probabilistic representations for Neumann solutions. A second result deals with Neumann problems on canonically compactifiable graphs with respect to the Royden boundary and provides conditions for unique solvability and analytic and probabilistic representations.}, language = {en} } @article{MeraTarkhanov2022, author = {Mera, Azal Jaafar Musa and Tarkhanov, Nikolai}, title = {An elliptic equation of finite index in a domain}, series = {Boletin de la Sociedad Matem{\´a}tica Mexicana}, volume = {28}, journal = {Boletin de la Sociedad Matem{\´a}tica Mexicana}, number = {2}, publisher = {Springer International}, address = {New York [u.a.]}, issn = {1405-213X}, doi = {10.1007/s40590-022-00442-7}, pages = {10}, year = {2022}, abstract = {We give an example of first order elliptic equation for a complex-valued function in a plane domain which has a finite number of linearly independent solutions for any right-hand side. No boundary value conditions are thus required.}, language = {en} } @article{KortenkampKuzleReitzKoncebovski2023, author = {Kortenkamp, Ulrich and Kuzle, Ana and Reitz-Koncebovski, Karen}, title = {Fachdidaktisches Wissen aus dem Fachwissen generieren}, series = {PSI-Potsdam: Ergebnisbericht zu den Aktivit{\"a}ten im Rahmen der Qualit{\"a}tsoffensive Lehrerbildung (2019-2023) (Potsdamer Beitr{\"a}ge zur Lehrerbildung und Bildungsforschung ; 3)}, journal = {PSI-Potsdam: Ergebnisbericht zu den Aktivit{\"a}ten im Rahmen der Qualit{\"a}tsoffensive Lehrerbildung (2019-2023) (Potsdamer Beitr{\"a}ge zur Lehrerbildung und Bildungsforschung ; 3)}, number = {3}, publisher = {Universit{\"a}tsverlag Potsdam}, address = {Potsdam}, isbn = {978-3-86956-568-2}, issn = {2626-3556}, doi = {10.25932/publishup-61760}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-617602}, pages = {171 -- 191}, year = {2023}, abstract = {Das Mathematik-Teilprojekt SPIES-M zielt auf eine st{\"a}rkere Professionsorientierung und die Verkn{\"u}pfung von Fachwissenschaft und Fachdidaktik in der universit{\"a}ren Lehrkr{\"a}ftebildung. Zu allen großen Inhaltsgebieten der Mathematik wurden neue Lehrveranstaltungen konzipiert und in den Studienordnungen s{\"a}mtlicher Lehr{\"a}mter Mathematik an der Universit{\"a}t Potsdam implementiert. F{\"u}r die Konzeption wurden theoriebasiert Gestaltungsprinzipien herausgearbeitet, die sowohl f{\"u}r das Design als auch f{\"u}r die Evaluation und Weiterentwicklung der Lehrveranstaltungen nach dem Design-Research-Ansatz genutzt werden k{\"o}nnen. Die Umsetzung der Gestaltungsprinzipien wird am Beispiel der Fundamentalen Idee der Proportionalit{\"a}t verdeutlicht und dabei aufgezeigt, wie Studierende dazu bef{\"a}higt werden k{\"o}nnen, fachdidaktisches Wissen aus fachmathematischen Inhalten zu generieren. Die Entwicklung des Professionswissens der Studierenden wird mithilfe unterschiedlicher Instrumente untersucht, um R{\"u}ckschl{\"u}sse auf die Wirksamkeit der neu konzipierten Lehrveranstaltungen zu ziehen. F{\"u}r die Untersuchungen im Mixed-Methods-Design werden neben Beobachtungen in Lehrveranstaltungen eigens konzipierte Wissenstests, Gruppeninterviews, Unterrichtsentw{\"u}rfe aus Praxisphasen und Lerntageb{\"u}cher genutzt. Die Studierendenperspektive wird durch Befragungen zur wahrgenommenen (Berufs-)Relevanz der Lehrveranstaltungen erhoben. Weiteres wesentliches Element der Begleitforschung ist die kollegiale Supervision durch sogenannte „Spies" (Spione), die die Veranstaltungen kriteriengeleitet beobachten und anschließend gemeinsam mit den Dozierenden reflektieren. Die bisherigen Ergebnisse werden hier pr{\"a}sentiert und hinsichtlich ihrer Implikationen diskutiert. Die im Projekt entwickelten Gestaltungsprinzipien als Werkzeug f{\"u}r Design und Evaluation sowie das Spies-Konzept der kollegialen Supervision werden f{\"u}r die Qualit{\"a}tsentwicklung von Lehrveranstaltungen zum Transfer vorgeschlagen.}, language = {de} } @article{BaerBandara2022, author = {B{\"a}r, Christian and Bandara, Lashi}, title = {Boundary value problems for general first-order elliptic differential operators}, series = {Journal of functional analysis}, volume = {282}, journal = {Journal of functional analysis}, number = {12}, publisher = {Elsevier}, address = {Amsterdam [u.a.]}, issn = {0022-1236}, doi = {10.1016/j.jfa.2022.109445}, pages = {69}, year = {2022}, abstract = {We study boundary value problems for first-order elliptic differential operators on manifolds with compact boundary. The adapted boundary operator need not be selfadjoint and the boundary condition need not be pseudo-local.We show the equivalence of various characterisations of elliptic boundary conditions and demonstrate how the boundary conditions traditionally considered in the literature fit in our framework. The regularity of the solutions up to the boundary is proven. We show that imposing elliptic boundary conditions yields a Fredholm operator if the manifold is compact. We provide examples which are conveniently treated by our methods.}, language = {en} } @phdthesis{MalemShinitski2023, author = {Malem-Shinitski, Noa}, title = {Bayesian inference and modeling for point processes with applications from neuronal activity to scene viewing}, doi = {10.25932/publishup-61495}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-614952}, school = {Universit{\"a}t Potsdam}, pages = {vii, 129}, year = {2023}, abstract = {Point processes are a common methodology to model sets of events. From earthquakes to social media posts, from the arrival times of neuronal spikes to the timing of crimes, from stock prices to disease spreading -- these phenomena can be reduced to the occurrences of events concentrated in points. Often, these events happen one after the other defining a time--series. Models of point processes can be used to deepen our understanding of such events and for classification and prediction. Such models include an underlying random process that generates the events. This work uses Bayesian methodology to infer the underlying generative process from observed data. Our contribution is twofold -- we develop new models and new inference methods for these processes. We propose a model that extends the family of point processes where the occurrence of an event depends on the previous events. This family is known as Hawkes processes. Whereas in most existing models of such processes, past events are assumed to have only an excitatory effect on future events, we focus on the newly developed nonlinear Hawkes process, where past events could have excitatory and inhibitory effects. After defining the model, we present its inference method and apply it to data from different fields, among others, to neuronal activity. The second model described in the thesis concerns a specific instance of point processes --- the decision process underlying human gaze control. This process results in a series of fixated locations in an image. We developed a new model to describe this process, motivated by the known Exploration--Exploitation dilemma. Alongside the model, we present a Bayesian inference algorithm to infer the model parameters. Remaining in the realm of human scene viewing, we identify the lack of best practices for Bayesian inference in this field. We survey four popular algorithms and compare their performances for parameter inference in two scan path models. The novel models and inference algorithms presented in this dissertation enrich the understanding of point process data and allow us to uncover meaningful insights.}, language = {en} } @article{ReichWeissmann2021, author = {Reich, Sebastian and Weissmann, Simon}, title = {Fokker-Planck particle systems for Bayesian inference: computational approaches}, series = {SIAM ASA journal on uncertainty quantification}, volume = {9}, journal = {SIAM ASA journal on uncertainty quantification}, number = {2}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {2166-2525}, doi = {10.1137/19M1303162}, pages = {446 -- 482}, year = {2021}, abstract = {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.}, language = {en} } @phdthesis{Schindler2023, author = {Schindler, Daniel}, title = {Mathematical modeling and simulation of protrusion-driven cell dynamics}, doi = {10.25932/publishup-61327}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-613275}, school = {Universit{\"a}t Potsdam}, pages = {VI, 161}, year = {2023}, abstract = {Amoeboid cell motility takes place in a variety of biomedical processes such as cancer metastasis, embryonic morphogenesis, and wound healing. In contrast to other forms of cell motility, it is mainly driven by substantial cell shape changes. Based on the interplay of explorative membrane protrusions at the front and a slower-acting membrane retraction at the rear, the cell moves in a crawling kind of way. Underlying these protrusions and retractions are multiple physiological processes resulting in changes of the cytoskeleton, a meshwork of different multi-functional proteins. The complexity and versatility of amoeboid cell motility raise the need for novel computational models based on a profound theoretical framework to analyze and simulate the dynamics of the cell shape. The objective of this thesis is the development of (i) a mathematical framework to describe contour dynamics in time and space, (ii) a computational model to infer expansion and retraction characteristics of individual cell tracks and to produce realistic contour dynamics, (iii) and a complementing Open Science approach to make the above methods fully accessible and easy to use. In this work, we mainly used single-cell recordings of the model organism Dictyostelium discoideum. Based on stacks of segmented microscopy images, we apply a Bayesian approach to obtain smooth representations of the cell membrane, so-called cell contours. We introduce a one-parameter family of regularized contour flows to track reference points on the contour (virtual markers) in time and space. This way, we define a coordinate system to visualize local geometric and dynamic quantities of individual contour dynamics in so-called kymograph plots. In particular, we introduce the local marker dispersion as a measure to identify membrane protrusions and retractions in a fully automated way. This mathematical framework is the basis of a novel contour dynamics model, which consists of three biophysiologically motivated components: one stochastic term, accounting for membrane protrusions, and two deterministic terms to control the shape and area of the contour, which account for membrane retractions. Our model provides a fully automated approach to infer protrusion and retraction characteristics from experimental cell tracks while being also capable of simulating realistic and qualitatively different contour dynamics. Furthermore, the model is used to classify two different locomotion types: the amoeboid and a so-called fan-shaped type. With the complementing Open Science approach, we ensure a high standard regarding the usability of our methods and the reproducibility of our research. In this context, we introduce our software publication named AmoePy, an open-source Python package to segment, analyze, and simulate amoeboid cell motility. Furthermore, we describe measures to improve its usability and extensibility, e.g., by detailed run instructions and an automatically generated source code documentation, and to ensure its functionality and stability, e.g., by automatic software tests, data validation, and a hierarchical package structure. The mathematical approaches of this work provide substantial improvements regarding the modeling and analysis of amoeboid cell motility. We deem the above methods, due to their generalized nature, to be of greater value for other scientific applications, e.g., varying organisms and experimental setups or the transition from unicellular to multicellular movement. Furthermore, we enable other researchers from different fields, i.e., mathematics, biophysics, and medicine, to apply our mathematical methods. By following Open Science standards, this work is of greater value for the cell migration community and a potential role model for other Open Science contributions.}, language = {en} } @article{ShlapunovTarchanov2021, author = {Shlapunov, Alexander and Tarchanov, Nikolaj Nikolaevič}, title = {An open mapping theorem for the Navier-Stokes type equations associated with the de Rham complex over R-n}, series = {Siberian electronic mathematical reports = Sibirskie ėlektronnye matematičeskie izvestija}, volume = {18}, journal = {Siberian electronic mathematical reports = Sibirskie ėlektronnye matematičeskie izvestija}, number = {2}, publisher = {Institut Matematiki Imeni S. L. Soboleva}, address = {Novosibirsk}, issn = {1813-3304}, doi = {10.33048/semi.2021.18.108}, pages = {1433 -- 1466}, year = {2021}, abstract = {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.}, language = {en} } @article{BeckusEliaz2021, author = {Beckus, Siegfried and Eliaz, Latif}, title = {Eigenfunctions growth of R-limits on graphs}, series = {Journal of spectral theory / European Mathematical Society}, volume = {11}, journal = {Journal of spectral theory / European Mathematical Society}, number = {4}, publisher = {EMS Press, an imprint of the European Mathematical Society - EMS - Publishing House GmbH, Institut f{\"u}r Mathematik, Technische Universit{\"a}t}, address = {Berlin}, issn = {1664-039X}, doi = {10.4171/JST/389}, pages = {1895 -- 1933}, year = {2021}, abstract = {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.}, language = {en} } @article{RoppLesurBaerenzungetal.2020, author = {Ropp, Guillaume and Lesur, Vincent and B{\"a}renzung, Julien and Holschneider, Matthias}, title = {Sequential modelling of the Earth's core magnetic field}, series = {Earth, Planets and Space}, volume = {72}, journal = {Earth, Planets and Space}, number = {1}, publisher = {Springer}, address = {New York}, issn = {1880-5981}, doi = {10.1186/s40623-020-01230-1}, pages = {15}, year = {2020}, abstract = {We describe a new, original approach to the modelling of the Earth's magnetic field. The overall objective of this study is to reliably render fast variations of the core field and its secular variation. This method combines a sequential modelling approach, a Kalman filter, and a correlation-based modelling step. Sources that most significantly contribute to the field measured at the surface of the Earth are modelled. Their separation is based on strong prior information on their spatial and temporal behaviours. We obtain a time series of model distributions which display behaviours similar to those of recent models based on more classic approaches, particularly at large temporal and spatial scales. Interesting new features and periodicities are visible in our models at smaller time and spatial scales. An important aspect of our method is to yield reliable error bars for all model parameters. These errors, however, are only as reliable as the description of the different sources and the prior information used are realistic. Finally, we used a slightly different version of our method to produce candidate models for the thirteenth edition of the International Geomagnetic Reference Field.}, language = {en} } @article{MatzkaStolleYamazakietal.2021, author = {Matzka, J{\"u}rgen and Stolle, Claudia and Yamazaki, Yosuke and Bronkalla, Oliver and Morschhauser, Achim}, title = {The geomagnetic Kp index and derived indices of geomagnetic activity}, series = {Space weather : the international journal of research and applications}, volume = {19}, journal = {Space weather : the international journal of research and applications}, number = {5}, publisher = {Wiley}, address = {New York}, issn = {1542-7390}, doi = {10.1029/2020SW002641}, pages = {21}, year = {2021}, abstract = {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.}, language = {en} } @article{PathirajaReichStannat2021, author = {Pathiraja, Sahani Darschika and Reich, Sebastian and Stannat, Wilhelm}, title = {McKean-Vlasov SDEs in nonlinear filtering}, series = {SIAM journal on control and optimization : a publication of the Society for Industrial and Applied Mathematics}, volume = {59}, journal = {SIAM journal on control and optimization : a publication of the Society for Industrial and Applied Mathematics}, number = {6}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {0363-0129}, doi = {10.1137/20M1355197}, pages = {4188 -- 4215}, year = {2021}, abstract = {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.}, language = {en} } @article{GarmendiaZambon2021, author = {Garmendia, Alfonso and Zambon, Marco}, title = {Quotients of singular foliations and Lie 2-group actions}, series = {Journal of noncommutative geometry}, volume = {15}, journal = {Journal of noncommutative geometry}, number = {4}, publisher = {EMS Press, an imprint of the European Mathematical Society - EMS - Publishing House GmbH, Institut f{\"u}r Mathematik, Technische Universit{\"a}t Berlin}, address = {Berlin}, issn = {1661-6952}, doi = {10.4171/JNCG/434}, pages = {1251 -- 1283}, year = {2021}, abstract = {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.}, language = {en} } @article{RuchiDubinkinaWiljes2021, author = {Ruchi, Sangeetika and Dubinkina, Svetlana and Wiljes, Jana de}, title = {Fast hybrid tempered ensemble transform filter formulation for Bayesian elliptical problems via Sinkhorn approximation}, series = {Nonlinear processes in geophysics / European Geosciences Union ; American Geophysical Union}, volume = {28}, journal = {Nonlinear processes in geophysics / European Geosciences Union ; American Geophysical Union}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1023-5809}, doi = {10.5194/npg-28-23-2021}, pages = {23 -- 41}, year = {2021}, abstract = {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.}, language = {en} } @article{SchannerMauerbergerKorteetal.2021, author = {Schanner, Maximilian Arthus and Mauerberger, Stefan and Korte, Monika and Holschneider, Matthias}, title = {Correlation based time evolution of the archeomagnetic field}, series = {Journal of geophysical research : JGR ; an international quarterly. B, Solid earth}, volume = {126}, journal = {Journal of geophysical research : JGR ; an international quarterly. B, Solid earth}, number = {7}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2020JB021548}, pages = {22}, year = {2021}, abstract = {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.}, language = {en} } @article{DemarisWidigsonIlvemarketal.2022, author = {D{\´e}maris, Alix and Widigson, Ella S. K. and Ilvemark, Johan F. K. F. and Steenholdt, Casper and Seidelin, Jakob B. and Huisinga, Wilhelm and Michelet, Robin and Aulin, Linda B. S. and Kloft, Charlotte}, title = {Ulcerative colitis and acute severe ulcerative colitis patients are overlooked in infliximab population pharmacokinetic models}, series = {Pharmaceutics / Molecular Diversity Preservation International}, volume = {14}, journal = {Pharmaceutics / Molecular Diversity Preservation International}, number = {10}, publisher = {MDPI}, address = {Basel}, issn = {1999-4923}, doi = {10.3390/pharmaceutics14102095}, pages = {32}, year = {2022}, abstract = {Ulcerative colitis (UC) is part of the inflammatory bowels diseases, and moderate to severe UC patients can be treated with anti-tumour necrosis alpha monoclonal antibodies, including infliximab (IFX). Even though treatment of UC patients by IFX has been in place for over a decade, many gaps in modelling of IFX PK in this population remain. This is even more true for acute severe UC (ASUC) patients for which early prediction of IFX pharmacokinetic (PK) could highly improve treatment outcome. Thus, this review aims to compile and analyse published population PK models of IFX in UC and ASUC patients, and to assess the current knowledge on disease activity impact on IFX PK. For this, a semi-systematic literature search was conducted, from which 26 publications including a population PK model analysis of UC patients receiving IFX therapy were selected. Amongst those, only four developed a model specifically for UC patients, and only three populations included severe UC patients. Investigations of disease activity impact on PK were reported in only 4 of the 14 models selected. In addition, the lack of reported model codes and assessment of predictive performance make the use of published models in a clinical setting challenging. Thus, more comprehensive investigation of PK in UC and ASUC is needed as well as more adequate reports on developed models and their evaluation in order to apply them in a clinical setting.}, language = {en} } @article{BaerenzungHolschneiderWichtetal.2020, author = {Baerenzung, Julien and Holschneider, Matthias and Wicht, Johannes and Lesur, Vincent and Sanchez, Sabrina}, title = {The Kalmag model as a candidate for IGRF-13}, series = {Earth, planets and space}, volume = {72}, journal = {Earth, planets and space}, number = {1}, publisher = {Springer}, address = {New York}, issn = {1880-5981}, doi = {10.1186/s40623-020-01295-y}, pages = {13}, year = {2020}, abstract = {We present a new model of the geomagnetic field spanning the last 20 years and called Kalmag. Deriving from the assimilation of CHAMP and Swarm vector field measurements, it separates the different contributions to the observable field through parameterized prior covariance matrices. To make the inverse problem numerically feasible, it has been sequentialized in time through the combination of a Kalman filter and a smoothing algorithm. The model provides reliable estimates of past, present and future mean fields and associated uncertainties. The version presented here is an update of our IGRF candidates; the amount of assimilated data has been doubled and the considered time window has been extended from [2000.5, 2019.74] to [2000.5, 2020.33].}, language = {en} } @article{HastermannReinhardtKleinetal.2021, author = {Hastermann, Gottfried and Reinhardt, Maria and Klein, Rupert and Reich, Sebastian}, title = {Balanced data assimilation for highly oscillatory mechanical systems}, series = {Communications in applied mathematics and computational science : CAMCoS}, volume = {16}, journal = {Communications in applied mathematics and computational science : CAMCoS}, number = {1}, publisher = {Mathematical Sciences Publishers}, address = {Berkeley}, issn = {1559-3940}, doi = {10.2140/camcos.2021.16.119}, pages = {119 -- 154}, year = {2021}, abstract = {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.}, language = {en} } @article{NassarHohmannMicheletetal.2022, author = {Nassar, Yomna M. and Hohmann, Nicolas and Michelet, Robin and Gottwalt, Katharina and Meid, Andreas D. and Burhenne, J{\"u}rgen and Huisinga, Wilhelm and Haefeli, Walter E. and Mikus, Gerd and Kloft, Charlotte}, title = {Quantification of the Time Course of CYP3A Inhibition, Activation, and Induction Using a Population Pharmacokinetic Model of Microdosed Midazolam Continuous Infusion}, series = {Clinical Pharmacokinetics}, volume = {61}, journal = {Clinical Pharmacokinetics}, number = {11}, publisher = {Springer}, address = {Northcote}, issn = {0312-5963}, doi = {10.1007/s40262-022-01175-6}, pages = {1595 -- 1607}, year = {2022}, abstract = {Background Cytochrome P450 (CYP) 3A contributes to the metabolism of many approved drugs. CYP3A perpetrator drugs can profoundly alter the exposure of CYP3A substrates. However, effects of such drug-drug interactions are usually reported as maximum effects rather than studied as time-dependent processes. Identification of the time course of CYP3A modulation can provide insight into when significant changes to CYP3A activity occurs, help better design drug-drug interaction studies, and manage drug-drug interactions in clinical practice. Objective We aimed to quantify the time course and extent of the in vivo modulation of different CYP3A perpetrator drugs on hepatic CYP3A activity and distinguish different modulatory mechanisms by their time of onset, using pharmacologically inactive intravenous microgram doses of the CYP3A-specific substrate midazolam, as a marker of CYP3A activity. Methods Twenty-four healthy individuals received an intravenous midazolam bolus followed by a continuous infusion for 10 or 36 h. Individuals were randomized into four arms: within each arm, two individuals served as a placebo control and, 2 h after start of the midazolam infusion, four individuals received the CYP3A perpetrator drug: voriconazole (inhibitor, orally or intravenously), rifampicin (inducer, orally), or efavirenz (activator, orally). After midazolam bolus administration, blood samples were taken every hour (rifampicin arm) or every 15 min (remaining study arms) until the end of midazolam infusion. A total of 1858 concentrations were equally divided between midazolam and its metabolite, 1'-hydroxymidazolam. A nonlinear mixed-effects population pharmacokinetic model of both compounds was developed using NONMEM (R). CYP3A activity modulation was quantified over time, as the relative change of midazolam clearance encountered by the perpetrator drug, compared to the corresponding clearance value in the placebo arm. Results Time course of CYP3A modulation and magnitude of maximum effect were identified for each perpetrator drug. While efavirenz CYP3A activation was relatively fast and short, reaching a maximum after approximately 2-3 h, the induction effect of rifampicin could only be observed after 22 h, with a maximum after approximately 28-30 h followed by a steep drop to almost baseline within 1-2 h. In contrast, the inhibitory impact of both oral and intravenous voriconazole was prolonged with a steady inhibition of CYP3A activity followed by a gradual increase in the inhibitory effect until the end of sampling at 8 h. Relative maximum clearance changes were +59.1\%, +46.7\%, -70.6\%, and -61.1\% for efavirenz, rifampicin, oral voriconazole, and intravenous voriconazole, respectively. Conclusions We could distinguish between different mechanisms of CYP3A modulation by the time of onset. Identification of the time at which clearance significantly changes, per perpetrator drug, can guide the design of an optimal sampling schedule for future drug-drug interaction studies. The impact of a short-term combination of different perpetrator drugs on the paradigm CYP3A substrate midazolam was characterized and can define combination intervals in which no relevant interaction is to be expected.}, language = {en} } @article{ClavierGuoPaychaetal.2020, author = {Clavier, Pierre and Guo, Li and Paycha, Sylvie and Zhang, Bin}, title = {Locality and renormalization: universal properties and integrals on trees}, series = {Journal of mathematical physics}, volume = {61}, journal = {Journal of mathematical physics}, number = {2}, publisher = {American Institute of Physics}, address = {College Park, Md.}, issn = {0022-2488}, doi = {10.1063/1.5116381}, pages = {19}, year = {2020}, abstract = {The purpose of this paper is to build an algebraic framework suited to regularize branched structures emanating from rooted forests and which encodes the locality principle. This is achieved by means of the universal properties in the locality framework of properly decorated rooted forests. These universal properties are then applied to derive the multivariate regularization of integrals indexed by rooted forests. We study their renormalization, along the lines of Kreimer's toy model for Feynman integrals.}, language = {en} } @misc{EhlenFloegeGoebeletal.2023, author = {Ehlen, Tobias and Fl{\"o}ge, Annie and G{\"o}bel, Franziska and Keller, Peter and Rœlly, Sylvie}, title = {{\"U}bungsbuch zur Stochastik}, editor = {Keller, Peter and Rœlly, Sylvie}, publisher = {Universit{\"a}tsverlag Potsdam}, address = {Potsdam}, isbn = {978-3-86956-563-7}, doi = {10.25932/publishup-59593}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-595939}, pages = {306}, year = {2023}, abstract = {Dieses Buch stellt {\"U}bungen zu den Grundbegriffen und Grunds{\"a}tzen der Stochastik und ihre L{\"o}sungen zur Verf{\"u}gung. So wie man Tonleitern in der Musik trainiert, so berechnet man {\"U}bungsaufgaben in der Mathematik. In diesem Sinne soll dieses {\"U}bungsbuch vor allem als Vorlage dienen f{\"u}r das eigenst{\"a}ndige, eigenverantwortliche Lernen und {\"U}ben. Die Sch{\"o}nheit und Einzigartigkeit der Wahrscheinlichkeitstheorie besteht darin, dass sie eine Vielzahl von realen Ph{\"a}nomenen modellieren kann. Daher findet man hier Aufgaben mit Verbindungen zur Geometrie, zu Gl{\"u}cksspielen, zur Versicherungsmathematik, zur Demographie und vielen anderen Themen.}, language = {de} } @article{BaerHanke2022, author = {B{\"a}r, Christian and Hanke, Bernhard}, title = {Local flexibility for open partial differential relations}, series = {Communications on pure and applied mathematics / issued by the Courant Institute of Mathematical Sciences, New York Univ.}, volume = {75}, journal = {Communications on pure and applied mathematics / issued by the Courant Institute of Mathematical Sciences, New York Univ.}, number = {6}, publisher = {Wiley}, address = {Hoboken}, issn = {0010-3640}, doi = {10.1002/cpa.21982}, pages = {1377 -- 1415}, year = {2022}, abstract = {We show that local deformations, near closed subsets, of solutions to open partial differential relations can be extended to global deformations, provided all but the highest derivatives stay constant along the subset. The applicability of this general result is illustrated by a number of examples, dealing with convex embeddings of hypersurfaces, differential forms, and lapse functions in Lorentzian geometry. The main application is a general approximation result by sections that have very restrictive local properties on open dense subsets. This shows, for instance, that given any K is an element of Double-struck capital R every manifold of dimension at least 2 carries a complete C-1,C- 1-metric which, on a dense open subset, is smooth with constant sectional curvature K. Of course, this is impossible for C-2-metrics in general.}, language = {en} } @article{RoosOtoba2021, author = {Roos, Saskia and Otoba, Nobuhiko}, title = {Scalar curvature and the multiconformal class of a direct product Riemannian manifold}, series = {Geometriae dedicata}, volume = {214}, journal = {Geometriae dedicata}, number = {1}, publisher = {Springer}, address = {Dordrecht}, issn = {0046-5755}, doi = {10.1007/s10711-021-00636-9}, pages = {801 -- 829}, year = {2021}, abstract = {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.}, language = {en} } @article{MauerbergerSchannerKorteetal.2020, author = {Mauerberger, Stefan and Schanner, Maximilian Arthus and Korte, Monika and Holschneider, Matthias}, title = {Correlation based snapshot models of the archeomagnetic field}, series = {Geophysical journal international}, volume = {223}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggaa336}, pages = {648 -- 665}, year = {2020}, abstract = {For the time stationary global geomagnetic field, a new modelling concept is presented. A Bayesian non-parametric approach provides realistic location dependent uncertainty estimates. Modelling related variabilities are dealt with systematically by making little subjective apriori assumptions. Rather than parametrizing the model by Gauss coefficients, a functional analytic approach is applied. The geomagnetic potential is assumed a Gaussian process to describe a distribution over functions. Apriori correlations are given by an explicit kernel function with non-informative dipole contribution. A refined modelling strategy is proposed that accommodates non-linearities of archeomagnetic observables: First, a rough field estimate is obtained considering only sites that provide full field vector records. Subsequently, this estimate supports the linearization that incorporates the remaining incomplete records. The comparison of results for the archeomagnetic field over the past 1000 yr is in general agreement with previous models while improved model uncertainty estimates are provided.}, language = {en} } @article{Clavier2021, author = {Clavier, Pierre J.}, title = {Borel-{\´E}calle resummation of a two-point function}, series = {Annales Henri Poincar{\´e} : a journal of theoretical and mathematical physics / ed. jointly by the Institut Henri Poincar{\´e} and by the Swiss Physical Society}, volume = {22}, journal = {Annales Henri Poincar{\´e} : a journal of theoretical and mathematical physics / ed. jointly by the Institut Henri Poincar{\´e} and by the Swiss Physical Society}, number = {6}, publisher = {Springer}, address = {Cham}, issn = {1424-0637}, doi = {10.1007/s00023-021-01057-w}, pages = {2103 -- 2136}, year = {2021}, abstract = {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.}, language = {en} } @misc{Dahl2023, type = {Master Thesis}, author = {Dahl, Dorothee Sophie}, title = {Zahlen in den Fingern}, doi = {10.25932/publishup-60762}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-607629}, school = {Universit{\"a}t Potsdam}, pages = {118}, year = {2023}, abstract = {Die Debatte {\"u}ber den Einsatz von digitalen Werkzeugen in der mathematischen Fr{\"u}hf{\"o}rderung ist hoch aktuell. Lernspiele werden konstruiert, mit dem Ziel, mathematisches, informelles Wissen aufzubauen und so einen besseren Schulstart zu erm{\"o}glichen. Doch allein die digitale und spielerische Aufarbeitung f{\"u}hrt nicht zwingend zu einem Lernerfolg. Daher ist es umso wichtiger, die konkrete Implementation der theoretischen Konstrukte und Interaktionsm{\"o}glichkeiten mit den Werkzeugen zu analysieren und passend aufzubereiten. In dieser Masterarbeit wird dazu exemplarisch ein mathematisches Lernspiel namens „Fingu" f{\"u}r den Einsatz im vorschulischen Bereich theoretisch und empirisch im Rahmen der Artifact-Centric Activity Theory (ACAT) untersucht. Dazu werden zun{\"a}chst die theoretischen Hintergr{\"u}nde zum Zahlensinn, Zahlbegriffserwerb, Teil-Ganze-Verst{\"a}ndnis, der Anzahlwahrnehmung und -bestimmung, den Anzahlvergleichen und der Anzahldarstellung mithilfe von Fingern gem{\"a}ß der Embodied Cognition sowie der Verwendung von digitalen Werkzeugen und Multi-Touch-Ger{\"a}ten umfassend beschrieben. Anschließend wird die App Fingu erkl{\"a}rt und dann theoretisch entlang des ACAT-Review-Guides analysiert. Zuletzt wird die selbstst{\"a}ndig durchgef{\"u}hrte Studie mit zehn Vorschulkindern erl{\"a}utert und darauf aufbauend Verbesserungs- und Entwicklungsm{\"o}glichkeiten der App auf wissenschaftlicher Grundlage beigetragen. F{\"u}r Fingu l{\"a}sst sich abschließend festhalten, dass viele Prozesse wie die (Quasi-)Simultanerfassung oder das Z{\"a}hlen gef{\"o}rdert werden k{\"o}nnen, f{\"u}r andere wie das Teil-Ganze-Verst{\"a}ndnis aber noch Anpassungen und/oder die Begleitung durch Erwachsene n{\"o}tig ist.}, language = {de} } @article{GottwaldReich2021, author = {Gottwald, Georg A. and Reich, Sebastian}, title = {Supervised learning from noisy observations}, series = {Physica : D, Nonlinear phenomena}, volume = {423}, journal = {Physica : D, Nonlinear phenomena}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0167-2789}, doi = {10.1016/j.physd.2021.132911}, pages = {15}, year = {2021}, abstract = {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.}, language = {en} } @article{RedmannFreitag2021, author = {Redmann, Martin and Freitag, Melina A.}, title = {Optimization based model order reduction for stochastic systems}, series = {Applied mathematics and computation}, volume = {398}, journal = {Applied mathematics and computation}, publisher = {Elsevier}, address = {New York}, issn = {0096-3003}, doi = {10.1016/j.amc.2020.125783}, pages = {18}, year = {2021}, abstract = {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.}, language = {en} } @misc{PereraBoeckmann2019, author = {Perera, Upeksha and B{\"o}ckmann, Christine}, title = {Solutions of direct and inverse even-order Sturm-Liouville problems using Magnus expansion}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1336}, issn = {1866-8372}, doi = {10.25932/publishup-47341}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-473414}, pages = {24}, year = {2019}, abstract = {In this paper 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 for up to eight), with a mix of (including non-separable and finite singular endpoints) boundary conditions, 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. The inverse SLP algorithm proposed by Barcilon (1974) is utilized in combination with the Magnus method so that a direct SLP of any (even) order and an inverse SLP of order two can be solved effectively.}, language = {en} } @article{PereraBoeckmann2019, author = {Perera, Upeksha and B{\"o}ckmann, Christine}, title = {Solutions of Direct and Inverse Even-Order Sturm-Liouville Problems Using Magnus Expansion}, series = {Mathematics}, volume = {7}, journal = {Mathematics}, number = {6}, publisher = {MDPI}, address = {Basel, Schweiz}, issn = {2227-7390}, doi = {10.3390/math7060544}, pages = {24}, year = {2019}, abstract = {In this paper 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 for up to eight), with a mix of (including non-separable and finite singular endpoints) boundary conditions, 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. The inverse SLP algorithm proposed by Barcilon (1974) is utilized in combination with the Magnus method so that a direct SLP of any (even) order and an inverse SLP of order two can be solved effectively.}, language = {en} } @misc{PornsawadSapsakulBoeckmann2019, author = {Pornsawad, Pornsarp and Sapsakul, Nantawan and B{\"o}ckmann, Christine}, title = {A modified asymptotical regularization of nonlinear ill-posed problems}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1335}, issn = {1866-8372}, doi = {10.25932/publishup-47343}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-473433}, pages = {19}, year = {2019}, abstract = {In this paper, we investigate the continuous version of modified iterative Runge-Kutta-type methods for nonlinear inverse ill-posed problems proposed in a previous work. The convergence analysis is proved under the tangential cone condition, a modified discrepancy principle, i.e., the stopping time T is a solution of ∥𝐹(𝑥𝛿(𝑇))-𝑦𝛿∥=𝜏𝛿+ for some 𝛿+>𝛿, and an appropriate source condition. We yield the optimal rate of convergence.}, language = {en} } @article{PornsawadSapsakulBoeckmann2019, author = {Pornsawad, Pornsarp and Sapsakul, Nantawan and B{\"o}ckmann, Christine}, title = {A modified asymptotical regularization of nonlinear ill-posed problems}, series = {Mathematics}, volume = {7}, journal = {Mathematics}, edition = {5}, publisher = {MDPI}, address = {Basel, Schweiz}, issn = {2227-7390}, doi = {10.3390/math7050419}, pages = {19}, year = {2019}, abstract = {In this paper, we investigate the continuous version of modified iterative Runge-Kutta-type methods for nonlinear inverse ill-posed problems proposed in a previous work. The convergence analysis is proved under the tangential cone condition, a modified discrepancy principle, i.e., the stopping time T is a solution of ∥𝐹(𝑥𝛿(𝑇))-𝑦𝛿∥=𝜏𝛿+ for some 𝛿+>𝛿, and an appropriate source condition. We yield the optimal rate of convergence.}, language = {en} } @article{LeungLeutbecherReichetal.2020, author = {Leung, Tsz Yan and Leutbecher, Martin and Reich, Sebastian and Shepherd, Theodore G.}, title = {Impact of the mesoscale range on error growth and the limits to atmospheric predictability}, series = {Journal of the atmospheric sciences}, volume = {77}, journal = {Journal of the atmospheric sciences}, number = {11}, publisher = {American Meteorological Soc.}, address = {Boston}, issn = {0022-4928}, doi = {10.1175/JAS-D-19-0346.1}, pages = {3769 -- 3779}, year = {2020}, abstract = {Global numerical weather prediction (NWP) models have begun to resolve the mesoscale k(-5/3) range of the energy spectrum, which is known to impose an inherently finite range of deterministic predictability per se as errors develop more rapidly on these scales than on the larger scales. However, the dynamics of these errors under the influence of the synoptic-scale k(-3) range is little studied. Within a perfect-model context, the present work examines the error growth behavior under such a hybrid spectrum in Lorenz's original model of 1969, and in a series of identical-twin perturbation experiments using an idealized two-dimensional barotropic turbulence model at a range of resolutions. With the typical resolution of today's global NWP ensembles, error growth remains largely uniform across scales. The theoretically expected fast error growth characteristic of a k(-5/3) spectrum is seen to be largely suppressed in the first decade of the mesoscale range by the synoptic-scale k(-3) range. However, it emerges once models become fully able to resolve features on something like a 20-km scale, which corresponds to a grid resolution on the order of a few kilometers.}, language = {en} } @article{BeckusBellissardDeNittis2020, author = {Beckus, Siegfried and Bellissard, Jean and De Nittis, Giuseppe}, title = {Spectral continuity for aperiodic quantum systems}, series = {Journal of mathematical physics}, volume = {61}, journal = {Journal of mathematical physics}, number = {12}, publisher = {American Institute of Physics}, address = {Melville, NY}, issn = {0022-2488}, doi = {10.1063/5.0011488}, pages = {19}, year = {2020}, abstract = {This work provides a necessary and sufficient condition for a symbolic dynamical system to admit a sequence of periodic approximations in the Hausdorff topology. The key result proved and applied here uses graphs that are called De Bruijn graphs, Rauzy graphs, or Anderson-Putnam complex, depending on the community. Combining this with a previous result, the present work justifies rigorously the accuracy and reliability of algorithmic methods used to compute numerically the spectra of a large class of self-adjoint operators. The so-called Hamiltonians describe the effective dynamic of a quantum particle in aperiodic media. No restrictions on the structure of these operators other than general regularity assumptions are imposed. In particular, nearest-neighbor correlation is not necessary. Examples for the Fibonacci and the Golay-Rudin-Shapiro sequences are explicitly provided illustrating this discussion. While the first sequence has been thoroughly studied by physicists and mathematicians alike, a shroud of mystery still surrounds the latter when it comes to spectral properties. In light of this, the present paper gives a new result here that might help uncovering a solution.}, language = {en} } @phdthesis{LopezValencia2023, author = {Lopez Valencia, Diego Andres}, title = {The Milnor-Moore and Poincar{\´e}-Birkhoff-Witt theorems in the locality set up and the polar structure of Shintani zeta functions}, doi = {10.25932/publishup-59421}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-594213}, school = {Universit{\"a}t Potsdam}, pages = {147}, year = {2023}, abstract = {This thesis bridges two areas of mathematics, algebra on the one hand with the Milnor-Moore theorem (also called Cartier-Quillen-Milnor-Moore theorem) as well as the Poincar{\´e}-Birkhoff-Witt theorem, and analysis on the other hand with Shintani zeta functions which generalise multiple zeta functions. The first part is devoted to an algebraic formulation of the locality principle in physics and generalisations of classification theorems such as Milnor-Moore and Poincar{\´e}-Birkhoff-Witt theorems to the locality framework. The locality principle roughly says that events that take place far apart in spacetime do not infuence each other. The algebraic formulation of this principle discussed here is useful when analysing singularities which arise from events located far apart in space, in order to renormalise them while keeping a memory of the fact that they do not influence each other. We start by endowing a vector space with a symmetric relation, named the locality relation, which keeps track of elements that are "locally independent". The pair of a vector space together with such relation is called a pre-locality vector space. This concept is extended to tensor products allowing only tensors made of locally independent elements. We extend this concept to the locality tensor algebra, and locality symmetric algebra of a pre-locality vector space and prove the universal properties of each of such structures. We also introduce the pre-locality Lie algebras, together with their associated locality universal enveloping algebras and prove their universal property. We later upgrade all such structures and results from the pre-locality to the locality context, requiring the locality relation to be compatible with the linear structure of the vector space. This allows us to define locality coalgebras, locality bialgebras, and locality Hopf algebras. Finally, all the previous results are used to prove the locality version of the Milnor-Moore and the Poincar{\´e}-Birkhoff-Witt theorems. It is worth noticing that the proofs presented, not only generalise the results in the usual (non-locality) setup, but also often use less tools than their counterparts in their non-locality counterparts. The second part is devoted to study the polar structure of the Shintani zeta functions. Such functions, which generalise the Riemman zeta function, multiple zeta functions, Mordell-Tornheim zeta functions, among others, are parametrised by matrices with real non-negative arguments. It is known that Shintani zeta functions extend to meromorphic functions with poles on afine hyperplanes. We refine this result in showing that the poles lie on hyperplanes parallel to the facets of certain convex polyhedra associated to the defining matrix for the Shintani zeta function. Explicitly, the latter are the Newton polytopes of the polynomials induced by the columns of the underlying matrix. We then prove that the coeficients of the equation which describes the hyperplanes in the canonical basis are either zero or one, similar to the poles arising when renormalising generic Feynman amplitudes. For that purpose, we introduce an algorithm to distribute weight over a graph such that the weight at each vertex satisfies a given lower bound.}, language = {en} } @article{SchindlerMoldenhawerStangeetal.2021, author = {Schindler, Daniel and Moldenhawer, Ted and Stange, Maike and Lepro, Valentino and Beta, Carsten and Holschneider, Matthias and Huisinga, Wilhelm}, title = {Analysis of protrusion dynamics in amoeboid cell motility by means of regularized contour flows}, series = {PLoS Computational Biology : a new community journal}, volume = {17}, journal = {PLoS Computational Biology : a new community journal}, number = {8}, publisher = {PLoS}, address = {San Fransisco}, issn = {1553-734X}, doi = {10.1371/journal.pcbi.1009268}, pages = {33}, year = {2021}, abstract = {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.}, language = {en} } @article{KellerPinchoverPogorzelski2019, author = {Keller, Matthias and Pinchover, Yehuda and Pogorzelski, Felix}, title = {Criticality theory for Schr{\"o}dinger operators on graphs}, series = {Journal of spectral theory}, volume = {10}, journal = {Journal of spectral theory}, number = {1}, publisher = {European Mathematical Society}, address = {Z{\"u}rich}, issn = {1664-039X}, doi = {10.4171/JST/286}, pages = {73 -- 114}, year = {2019}, abstract = {We study Schrodinger operators given by positive quadratic forms on infinite graphs. From there, we develop a criticality theory for Schrodinger operators on general weighted graphs.}, language = {en} } @article{KellerPinchoverPogorzelski2020, author = {Keller, Matthias and Pinchover, Yehuda and Pogorzelski, Felix}, title = {From hardy to rellich inequalities on graphs}, series = {Proceedings of the London Mathematical Society}, volume = {122}, journal = {Proceedings of the London Mathematical Society}, number = {3}, publisher = {Wiley}, address = {Hoboken}, issn = {0024-6115}, doi = {10.1112/plms.12376}, pages = {458 -- 477}, year = {2020}, abstract = {We show how to deduce Rellich inequalities from Hardy inequalities on infinite graphs. Specifically, the obtained Rellich inequality gives an upper bound on a function by the Laplacian of the function in terms of weighted norms. These weights involve the Hardy weight and a function which satisfies an eikonal inequality. The results are proven first for Laplacians and are extended to Schrodinger operators afterwards.}, language = {en} } @article{MaierWiljesHartungetal.2022, author = {Maier, Corinna Sabrina and Wiljes, Jana de and Hartung, Niklas and Kloft, Charlotte and Huisinga, Wilhelm}, title = {A continued learning approach for model-informed precision dosing}, series = {CPT: pharmacometrics \& systems pharmacology}, volume = {11}, journal = {CPT: pharmacometrics \& systems pharmacology}, number = {2}, publisher = {London}, address = {Nature Publ. Group}, issn = {2163-8306}, doi = {10.1002/psp4.12745}, pages = {185 -- 198}, year = {2022}, abstract = {Model-informed precision dosing (MIPD) is a quantitative dosing framework that combines prior knowledge on the drug-disease-patient system with patient data from therapeutic drug/ biomarker monitoring (TDM) to support individualized dosing in ongoing treatment. Structural models and prior parameter distributions used in MIPD approaches typically build on prior clinical trials that involve only a limited number of patients selected according to some exclusion/inclusion criteria. Compared to the prior clinical trial population, the patient population in clinical practice can be expected to also include altered behavior and/or increased interindividual variability, the extent of which, however, is typically unknown. Here, we address the question of how to adapt and refine models on the level of the model parameters to better reflect this real-world diversity. We propose an approach for continued learning across patients during MIPD using a sequential hierarchical Bayesian framework. The approach builds on two stages to separate the update of the individual patient parameters from updating the population parameters. Consequently, it enables continued learning across hospitals or study centers, because only summary patient data (on the level of model parameters) need to be shared, but no individual TDM data. We illustrate this continued learning approach with neutrophil-guided dosing of paclitaxel. The present study constitutes an important step toward building confidence in MIPD and eventually establishing MIPD increasingly in everyday therapeutic use.}, language = {en} } @article{GottwaldReich2021, author = {Gottwald, Georg A. and Reich, Sebastian}, title = {Combining machine learning and data assimilation to forecast dynamical systems from noisy partial observations}, series = {Chaos : an interdisciplinary journal of nonlinear science}, volume = {31}, journal = {Chaos : an interdisciplinary journal of nonlinear science}, number = {10}, publisher = {AIP}, address = {Melville}, issn = {1054-1500}, doi = {10.1063/5.0066080}, pages = {8}, year = {2021}, abstract = {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.}, language = {en} } @article{DeneckeHounnon2021, author = {Denecke, Klaus-Dieter and Hounnon, Hippolyte}, title = {Partial Menger algebras of terms}, series = {Asian-European journal of mathematics}, volume = {14}, journal = {Asian-European journal of mathematics}, number = {06}, publisher = {World Scientific}, address = {Singapore}, issn = {1793-5571}, doi = {10.1142/S1793557121500923}, pages = {14}, year = {2021}, abstract = {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}, language = {en} } @article{GaidzikPathirajaSaalfeldetal.2020, author = {Gaidzik, Franziska and Pathiraja, Sahani Darschika and Saalfeld, Sylvia and Stucht, Daniel and Speck, Oliver and Thevenin, Dominique and Janiga, Gabor}, title = {Hemodynamic data assimilation in a subject-specific circle of Willis geometry}, series = {Clinical Neuroradiology}, volume = {31}, journal = {Clinical Neuroradiology}, number = {3}, publisher = {Springer}, address = {Heidelberg}, issn = {1869-1439}, doi = {10.1007/s00062-020-00959-2}, pages = {643 -- 651}, year = {2020}, abstract = {Purpose The anatomy of the circle of Willis (CoW), the brain's main arterial blood supply system, strongly differs between individuals, resulting in highly variable flow fields and intracranial vascularization patterns. To predict subject-specific hemodynamics with high certainty, we propose a data assimilation (DA) approach that merges fully 4D phase-contrast magnetic resonance imaging (PC-MRI) data with a numerical model in the form of computational fluid dynamics (CFD) simulations. Methods To the best of our knowledge, this study is the first to provide a transient state estimate for the three-dimensional velocity field in a subject-specific CoW geometry using DA. High-resolution velocity state estimates are obtained using the local ensemble transform Kalman filter (LETKF). Results Quantitative evaluation shows a considerable reduction (up to 90\%) in the uncertainty of the velocity field state estimate after the data assimilation step. Velocity values in vessel areas that are below the resolution of the PC-MRI data (e.g., in posterior communicating arteries) are provided. Furthermore, the uncertainty of the analysis-based wall shear stress distribution is reduced by a factor of 2 for the data assimilation approach when compared to the CFD model alone. Conclusion This study demonstrates the potential of data assimilation to provide detailed information on vascular flow, and to reduce the uncertainty in such estimates by combining various sources of data in a statistically appropriate fashion.}, language = {en} } @article{BaerMazzeo2021, author = {B{\"a}r, Christian and Mazzeo, Rafe}, title = {Manifolds with many Rarita-Schwinger fields}, series = {Communications in mathematical physics}, volume = {384}, journal = {Communications in mathematical physics}, number = {1}, publisher = {Springer}, address = {Berlin}, issn = {0010-3616}, doi = {10.1007/s00220-021-04030-0}, pages = {533 -- 548}, year = {2021}, abstract = {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.}, language = {en} } @article{MariucciRaySzabo2020, author = {Mariucci, Ester and Ray, Kolyan and Szabo, Botond}, title = {A Bayesian nonparametric approach to log-concave density estimation}, series = {Bernoulli : official journal of the Bernoulli Society for Mathematical Statistics and Probability}, volume = {26}, journal = {Bernoulli : official journal of the Bernoulli Society for Mathematical Statistics and Probability}, number = {2}, publisher = {International Statistical Institute}, address = {The Hague}, issn = {1350-7265}, doi = {10.3150/19-BEJ1139}, pages = {1070 -- 1097}, year = {2020}, abstract = {The estimation of a log-concave density on R is a canonical problem in the area of shape-constrained nonparametric inference. We present a Bayesian nonparametric approach to this problem based on an exponentiated Dirichlet process mixture prior and show that the posterior distribution converges to the log-concave truth at the (near-) minimax rate in Hellinger distance. Our proof proceeds by establishing a general contraction result based on the log-concave maximum likelihood estimator that prevents the need for further metric entropy calculations. We further present computationally more feasible approximations and both an empirical and hierarchical Bayes approach. All priors are illustrated numerically via simulations.}, language = {en} } @article{HydeEvans2022, author = {Hyde, Stephen T. and Evans, Myfanwy E.}, title = {Symmetric tangled Platonic polyhedra}, series = {Proceedings of the National Academy of Sciences of the United States of America}, volume = {119}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {1}, publisher = {National Acad. of Sciences}, address = {Washington}, issn = {0027-8424}, doi = {10.1073/pnas.2110345118}, pages = {10}, year = {2022}, abstract = {Conventional embeddings of the edge-graphs of Platonic polyhedra, {f,z}, where f,z denote the number of edges in each face and the edge-valence at each vertex, respectively, are untangled in that they can be placed on a sphere (S-2) such that distinct edges do not intersect, analogous to unknotted loops, which allow crossing-free drawings of S-1 on the sphere. The most symmetric (flag-transitive) realizations of those polyhedral graphs are those of the classical Platonic polyhedra, whose symmetries are *2fz, according to Conway's two-dimensional (2D) orbifold notation (equivalent to Schonflies symbols I-h, O-h, and T-d). Tangled Platonic {f,z} polyhedra-which cannot lie on the sphere without edge-crossings-are constructed as windings of helices with three, five, seven,... strands on multigenus surfaces formed by tubifying the edges of conventional Platonic polyhedra, have (chiral) symmetries 2fz (I, O, and T), whose vertices, edges, and faces are symmetrically identical, realized with two flags. The analysis extends to the "theta(z)" polyhedra, {2,z}. The vertices of these symmetric tangled polyhedra overlap with those of the Platonic polyhedra; however, their helicity requires curvilinear (or kinked) edges in all but one case. We show that these 2fz polyhedral tangles are maximally symmetric; more symmetric embeddings are necessarily untangled. On one hand, their topologies are very constrained: They are either self-entangled graphs (analogous to knots) or mutually catenated entangled compound polyhedra (analogous to links). On the other hand, an endless variety of entanglements can be realized for each topology. Simpler examples resemble patterns observed in synthetic organometallic materials and clathrin coats in vivo.}, language = {en} } @article{HermannHumbert2020, author = {Hermann, Andreas and Humbert, Emmanuel}, title = {Mass functions of a compact manifold}, series = {Journal of geometry and physics : JGP}, volume = {154}, journal = {Journal of geometry and physics : JGP}, publisher = {Elsevier}, address = {Amsterdam [u.a.]}, issn = {0393-0440}, doi = {10.1016/j.geomphys.2020.103650}, pages = {14}, year = {2020}, abstract = {Let M be a compact manifold of dimension n. In this paper, we introduce the Mass Function a >= 0 bar right arrow X-+(M)(a) (resp. a >= 0 bar right arrow X--(M)(a)) which is defined as the supremum (resp. infimum) of the masses of all metrics on M whose Yamabe constant is larger than a and which are flat on a ball of radius 1 and centered at a point p is an element of M. Here, the mass of a metric flat around p is the constant term in the expansion of the Green function of the conformal Laplacian at p. We show that these functions are well defined and have many properties which allow to obtain applications to the Yamabe invariant (i.e. the supremum of Yamabe constants over the set of all metrics on M).}, language = {en} } @article{KolbeEvans2022, author = {Kolbe, Benedikt Maximilian and Evans, Myfanwy E.}, title = {Enumerating isotopy classes of tilings guided by the symmetry of triply}, series = {Siam journal on applied algebra and geometry}, volume = {6}, journal = {Siam journal on applied algebra and geometry}, number = {1}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {2470-6566}, doi = {10.1137/20M1358943}, pages = {1 -- 40}, year = {2022}, abstract = {We present a technique for the enumeration of all isotopically distinct ways of tiling a hyperbolic surface of finite genus, possibly nonorientable and with punctures and boundary. This generalizes the enumeration using Delaney--Dress combinatorial tiling theory of combinatorial classes of tilings to isotopy classes of tilings. To accomplish this, we derive an action of the mapping class group of the orbifold associated to the symmetry group of a tiling on the set of tilings. We explicitly give descriptions and presentations of semipure mapping class groups and of tilings as decorations on orbifolds. We apply this enumerative result to generate an array of isotopically distinct tilings of the hyperbolic plane with symmetries generated by rotations that are commensurate with the threedimensional symmetries of the primitive, diamond, and gyroid triply periodic minimal surfaces, which have relevance to a variety of physical systems.}, language = {en} } @article{GueneysuKeller2020, author = {G{\"u}neysu, Batu and Keller, Matthias}, title = {Feynman path integrals for magnetic Schr{\"o}dinger operators on infinite weighted graphs}, series = {Journal d'analyse math{\´e}matique}, volume = {141}, journal = {Journal d'analyse math{\´e}matique}, number = {2}, publisher = {The Magnes Press, the Hebrew Univ.}, address = {Jerusalem}, issn = {0021-7670}, doi = {10.1007/s11854-020-0110-y}, pages = {751 -- 770}, year = {2020}, abstract = {We prove a Feynman path integral formula for the unitary group exp(-itL(nu,theta)), t >= 0, associated with a discrete magnetic Schrodinger operator L-nu,L-theta on a large class of weighted infinite graphs. As a consequence, we get a new Kato-Simon estimate vertical bar exp(- itL(nu,theta))(x,y)vertical bar <= exp( -tL(-deg,0))(x,y), which controls the unitary group uniformly in the potentials in terms of a Schrodinger semigroup, where the potential deg is the weighted degree function of the graph.}, language = {en} } @article{AzzaliPaycha2020, author = {Azzali, Sara and Paycha, Sylvie}, title = {Spectral zeta-invariants lifted to coverings}, series = {Transactions of the American Mathematical Society}, volume = {373}, journal = {Transactions of the American Mathematical Society}, number = {9}, publisher = {American Mathematical Society}, address = {Providence, RI}, issn = {0002-9947}, doi = {10.1090/tran/8067}, pages = {6185 -- 6226}, year = {2020}, abstract = {The canonical trace and the Wodzicki residue on classical pseudo-differential operators on a closed manifold are characterised by their locality and shown to be preserved under lifting to the universal covering as a result of their local feature. As a consequence, we lift a class of spectral zeta-invariants using lifted defect formulae which express discrepancies of zeta-regularised traces in terms of Wodzicki residues. We derive Atiyah's L-2-index theorem as an instance of the Z(2)-graded generalisation of the canonical lift of spectral zeta-invariants and we show that certain lifted spectral zeta-invariants for geometric operators are integrals of Pontryagin and Chern forms.}, language = {en} } @article{HoudebertZass2022, author = {Houdebert, Pierre and Zass, Alexander}, title = {An explicit Dobrushin uniqueness region for Gibbs point processes with repulsive interactions}, series = {Journal of applied probability / Applied Probability Trust}, volume = {59}, journal = {Journal of applied probability / Applied Probability Trust}, number = {2}, publisher = {Cambridge Univ. Press}, address = {Cambridge}, issn = {0021-9002}, doi = {10.1017/jpr.2021.70}, pages = {541 -- 555}, year = {2022}, abstract = {We present a uniqueness result for Gibbs point processes with interactions that come from a non-negative pair potential; in particular, we provide an explicit uniqueness region in terms of activity z and inverse temperature beta. The technique used relies on applying to the continuous setting the classical Dobrushin criterion. We also present a comparison to the two other uniqueness methods of cluster expansion and disagreement percolation, which can also be applied for this type of interaction.}, language = {en} } @article{KluweMicheletMuellerSchoelletal.2020, author = {Kluwe, Franziska and Michelet, Robin and M{\"u}ller-Sch{\"o}ll, Anna and Maier, Corinna and Klopp-Schulze, Lena and van Dyk, Madele and Mikus, Gerd and Huisinga, Wilhelm and Kloft, Charlotte}, title = {Perspectives on model-informed precision dosing in the digital health era}, series = {Clinical pharmacology \& therapeutics}, volume = {109}, journal = {Clinical pharmacology \& therapeutics}, number = {1}, publisher = {Wiley}, address = {Hoboken}, issn = {0009-9236}, doi = {10.1002/cpt.2049}, pages = {29 -- 36}, year = {2020}, language = {en} } @article{MaoutsaReichOpper2020, author = {Maoutsa, Dimitra and Reich, Sebastian and Opper, Manfred}, title = {Interacting particle solutions of Fokker-Planck equations through gradient-log-density estimation}, series = {Entropy}, volume = {22}, journal = {Entropy}, number = {8}, publisher = {MDPI}, address = {Basel}, issn = {1099-4300}, doi = {10.3390/e22080802}, pages = {35}, year = {2020}, abstract = {Fokker-Planck equations are extensively employed in various scientific fields as they characterise the behaviour of stochastic systems at the level of probability density functions. Although broadly used, they allow for analytical treatment only in limited settings, and often it is inevitable to resort to numerical solutions. Here, we develop a computational approach for simulating the time evolution of Fokker-Planck solutions in terms of a mean field limit of an interacting particle system. The interactions between particles are determined by the gradient of the logarithm of the particle density, approximated here by a novel statistical estimator. The performance of our method shows promising results, with more accurate and less fluctuating statistics compared to direct stochastic simulations of comparable particle number. Taken together, our framework allows for effortless and reliable particle-based simulations of Fokker-Planck equations in low and moderate dimensions. The proposed gradient-log-density estimator is also of independent interest, for example, in the context of optimal control.}, language = {en} } @article{LudewigRosenberger2020, author = {Ludewig, Matthias and Rosenberger, Elke}, title = {Asymptotic eigenfunctions for Schr{\"o}dinger operators on a vector bundle}, series = {Reviews in mathematical physics}, volume = {32}, journal = {Reviews in mathematical physics}, number = {7}, publisher = {World Scientific}, address = {Singapore}, issn = {0129-055X}, doi = {10.1142/S0129055X20500208}, pages = {28}, year = {2020}, abstract = {In the limit (h) over bar -> 0, we analyze a class of Schr{\"o}dinger operators H-(h) over bar = (h) over bar L-2 + (h) over barW + V .id(epsilon) acting on sections of a vector bundle epsilon over a Riemannian manifold M where L is a Laplace type operator, W is an endomorphism field and the potential energy V has a non-degenerate minimum at some point p is an element of M. We construct quasimodes of WKB-type near p for eigenfunctions associated with the low-lying eigenvalues of H-(h) over bar. These are obtained from eigenfunctions of the associated harmonic oscillator H-p,H-(h) over bar at p, acting on smooth functions on the tangent space.}, language = {en} } @article{HartungBorghardt2020, author = {Hartung, Niklas and Borghardt, Jens Markus}, title = {A mechanistic framework for a priori pharmacokinetic predictions of orally inhaled drugs}, series = {PLoS Computational Biology : a new community journal}, volume = {16}, journal = {PLoS Computational Biology : a new community journal}, number = {12}, publisher = {PLoS}, address = {San Fransisco}, issn = {1553-734X}, doi = {10.1371/journal.pcbi.1008466}, pages = {24}, year = {2020}, abstract = {Author summary
The use of orally inhaled drugs for treating lung diseases is appealing since they have the potential for lung selectivity, i.e. high exposure at the site of action -the lung- without excessive side effects. However, the degree of lung selectivity depends on a large number of factors, including physiochemical properties of drug molecules, patient disease state, and inhalation devices. To predict the impact of these factors on drug exposure and thereby to understand the characteristics of an optimal drug for inhalation, we develop a predictive mathematical framework (a "pharmacokinetic model"). In contrast to previous approaches, our model allows combining knowledge from different sources appropriately and its predictions were able to adequately predict different sets of clinical data. Finally, we compare the impact of different factors and find that the most important factors are the size of the inhaled particles, the affinity of the drug to the lung tissue, as well as the rate of drug dissolution in the lung. In contrast to the common belief, the solubility of a drug in the lining fluids is not found to be relevant. These findings are important to understand how inhaled drugs should be designed to achieve best treatment results in patients.
The fate of orally inhaled drugs is determined by pulmonary pharmacokinetic processes such as particle deposition, pulmonary drug dissolution, and mucociliary clearance. Even though each single process has been systematically investigated, a quantitative understanding on the interaction of processes remains limited and therefore identifying optimal drug and formulation characteristics for orally inhaled drugs is still challenging. To investigate this complex interplay, the pulmonary processes can be integrated into mathematical models. However, existing modeling attempts considerably simplify these processes or are not systematically evaluated against (clinical) data. In this work, we developed a mathematical framework based on physiologically-structured population equations to integrate all relevant pulmonary processes mechanistically. A tailored numerical resolution strategy was chosen and the mechanistic model was evaluated systematically against data from different clinical studies. Without adapting the mechanistic model or estimating kinetic parameters based on individual study data, the developed model was able to predict simultaneously (i) lung retention profiles of inhaled insoluble particles, (ii) particle size-dependent pharmacokinetics of inhaled monodisperse particles, (iii) pharmacokinetic differences between inhaled fluticasone propionate and budesonide, as well as (iv) pharmacokinetic differences between healthy volunteers and asthmatic patients. Finally, to identify the most impactful optimization criteria for orally inhaled drugs, the developed mechanistic model was applied to investigate the impact of input parameters on both the pulmonary and systemic exposure. Interestingly, the solubility of the inhaled drug did not have any relevant impact on the local and systemic pharmacokinetics. Instead, the pulmonary dissolution rate, the particle size, the tissue affinity, and the systemic clearance were the most impactful potential optimization parameters. In the future, the developed prediction framework should be considered a powerful tool for identifying optimal drug and formulation characteristics.}, language = {en} } @article{FladFladHarutyunyanSchulze2020, author = {Flad, Heinz-J{\"u}rgen and Flad-Harutyunyan, Gohar and Schulze, Bert-Wolfgang}, title = {Explicit Green operators for quantum mechanical Hamiltonians}, series = {Asian-European journal of mathematics : AEJM}, volume = {13}, journal = {Asian-European journal of mathematics : AEJM}, number = {7}, publisher = {World Scientific}, address = {Singapore}, issn = {1793-5571}, doi = {10.1142/S1793557120501223}, pages = {64}, year = {2020}, abstract = {We extend our approach of asymptotic parametrix construction for Hamiltonian operators from conical to edge-type singularities which is applicable to coalescence points of two particles of the helium atom and related two electron systems including the hydrogen molecule. Up to second-order, we have calculated the symbols of an asymptotic parametrix of the nonrelativistic Hamiltonian of the helium atom within the Born-Oppenheimer approximation and provide explicit formulas for the corresponding Green operators which encode the asymptotic behavior of the eigenfunctions near an edge.}, language = {en} } @article{LudewigRoos2020, author = {Ludewig, Matthias and Roos, Saskia}, title = {The chiral anomaly of the free fermion in functorial field theory}, series = {Annales Henri Poincar{\´e} : a journal of theoretical and mathematical physics}, volume = {21}, journal = {Annales Henri Poincar{\´e} : a journal of theoretical and mathematical physics}, number = {4}, publisher = {Springer International Publishing AG}, address = {Cham (ZG)}, issn = {1424-0637}, doi = {10.1007/s00023-020-00893-6}, pages = {1191 -- 1233}, year = {2020}, abstract = {When trying to cast the free fermion in the framework of functorial field theory, its chiral anomaly manifests in the fact that it assigns the determinant of the Dirac operator to a top-dimensional closed spin manifold, which is not a number as expected, but an element of a complex line. In functorial field theory language, this means that the theory is twisted, which gives rise to an anomaly theory. In this paper, we give a detailed construction of this anomaly theory, as a functor that sends manifolds to infinite-dimensional Clifford algebras and bordisms to bimodules.}, language = {en} } @article{BandaraBryan2020, author = {Bandara, Lashi and Bryan, Paul}, title = {Heat kernels and regularity for rough metrics on smooth manifolds}, series = {Mathematische Nachrichten}, volume = {293}, journal = {Mathematische Nachrichten}, number = {12}, publisher = {Wiley-VCH}, address = {Weinheim}, issn = {0025-584X}, doi = {10.1002/mana.201800459}, pages = {2255 -- 2270}, year = {2020}, abstract = {We consider rough metrics on smooth manifolds and corresponding Laplacians induced by such metrics. We demonstrate that globally continuous heat kernels exist and are Holder continuous locally in space and time. This is done via local parabolic Harnack estimates for weak solutions of operators in divergence form with bounded measurable coefficients in weighted Sobolev spaces.}, language = {en} } @article{KolbeEvans2020, author = {Kolbe, Benedikt Maximilian and Evans, Myfanwy E.}, title = {Isotopic tiling theory for hyperbolic surfaces}, series = {Geometriae dedicata}, volume = {212}, journal = {Geometriae dedicata}, number = {1}, publisher = {Springer}, address = {Dordrecht}, issn = {0046-5755}, doi = {10.1007/s10711-020-00554-2}, pages = {177 -- 204}, year = {2020}, abstract = {In this paper, we develop the mathematical tools needed to explore isotopy classes of tilings on hyperbolic surfaces of finite genus, possibly nonorientable, with boundary, and punctured. More specifically, we generalize results on Delaney-Dress combinatorial tiling theory using an extension of mapping class groups to orbifolds, in turn using this to study tilings of covering spaces of orbifolds. Moreover, we study finite subgroups of these mapping class groups. Our results can be used to extend the Delaney-Dress combinatorial encoding of a tiling to yield a finite symbol encoding the complexity of an isotopy class of tilings. The results of this paper provide the basis for a complete and unambiguous enumeration of isotopically distinct tilings of hyperbolic surfaces.}, language = {en} } @article{PereraBoeckmann2020, author = {Perera, Upeksha and B{\"o}ckmann, Christine}, title = {Solutions of Sturm-Liouville problems}, series = {Mathematics}, volume = {8}, journal = {Mathematics}, number = {11}, publisher = {MDPI}, address = {Basel}, issn = {2227-7390}, doi = {10.3390/math8112074}, pages = {14}, year = {2020}, abstract = {This paper further improves the Lie group method with Magnus expansion proposed in a previous paper by the authors, to solve some types of direct singular Sturm-Liouville problems. 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 higher order (for n=2,4) are 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 a method that can be adapted successfully for solving a direct (regular/singular) or inverse Sturm-Liouville problem (SLP) of an arbitrary order with arbitrary boundary conditions.}, language = {en} } @article{BeckusPinchover2020, author = {Beckus, Siegfried and Pinchover, Yehuda}, title = {Shnol-type theorem for the Agmon ground state}, series = {Journal of spectral theory}, volume = {10}, journal = {Journal of spectral theory}, number = {2}, publisher = {EMS Publishing House}, address = {Z{\"u}rich}, issn = {1664-039X}, doi = {10.4171/JST/296}, pages = {355 -- 377}, year = {2020}, abstract = {LetH be a Schrodinger operator defined on a noncompact Riemannianmanifold Omega, and let W is an element of L-infinity (Omega; R). Suppose that the operator H + W is critical in Omega, and let phi be the corresponding Agmon ground state. We prove that if u is a generalized eigenfunction ofH satisfying vertical bar u vertical bar <= C-phi in Omega for some constant C > 0, then the corresponding eigenvalue is in the spectrum of H. The conclusion also holds true if for some K is an element of Omega the operator H admits a positive solution in (Omega) over bar = Omega \ K, and vertical bar u vertical bar <= C psi in (Omega) over bar for some constant C > 0, where psi is a positive solution of minimal growth in a neighborhood of infinity in Omega. Under natural assumptions, this result holds also in the context of infinite graphs, and Dirichlet forms.}, language = {en} } @article{BlanchardMuecke2020, author = {Blanchard, Gilles and M{\"u}cke, Nicole}, title = {Kernel regression, minimax rates and effective dimensionality}, series = {Analysis and applications}, volume = {18}, journal = {Analysis and applications}, number = {4}, publisher = {World Scientific}, address = {New Jersey}, issn = {0219-5305}, doi = {10.1142/S0219530519500258}, pages = {683 -- 696}, year = {2020}, abstract = {We investigate if kernel regularization methods can achieve minimax convergence rates over a source condition regularity assumption for the target function. These questions have been considered in past literature, but only under specific assumptions about the decay, typically polynomial, of the spectrum of the the kernel mapping covariance operator. In the perspective of distribution-free results, we investigate this issue under much weaker assumption on the eigenvalue decay, allowing for more complex behavior that can reflect different structure of the data at different scales.}, language = {en} } @article{MazzonettoSalimova2020, author = {Mazzonetto, Sara and Salimova, Diyora}, title = {Existence, uniqueness, and numerical approximations for stochastic burgers equations}, series = {Stochastic analysis and applications}, volume = {38}, journal = {Stochastic analysis and applications}, number = {4}, publisher = {Taylor \& Francis Group}, address = {Philadelphia}, issn = {0736-2994}, doi = {10.1080/07362994.2019.1709503}, pages = {623 -- 646}, year = {2020}, abstract = {In this article, we propose an all-in-one statement which includes existence, uniqueness, regularity, and numerical approximations of mild solutions for a class of stochastic partial differential equations (SPDEs) with non-globally monotone nonlinearities. The proof of this result exploits the properties of an existing fully explicit space-time discrete approximation scheme, in particular the fact that it satisfies suitable a priori estimates. We also obtain almost sure and strong convergence of the approximation scheme to the mild solutions of the considered SPDEs. We conclude by applying the main result of the article to the stochastic Burgers equations with additive space-time white noise.}, language = {en} } @phdthesis{Gehring2023, author = {Gehring, Penelope}, title = {Non-local boundary conditions for the spin Dirac operator on spacetimes with timelike boundary}, doi = {10.25932/publishup-57775}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-577755}, school = {Universit{\"a}t Potsdam}, pages = {100}, year = {2023}, abstract = {Non-local boundary conditions - for example the Atiyah-Patodi-Singer (APS) conditions - for Dirac operators on Riemannian manifolds are rather well-understood, while not much is known for such operators on Lorentzian manifolds. Recently, B{\"a}r and Strohmaier [15] and Drago, Große, and Murro [27] introduced APS-like conditions for the spin Dirac operator on Lorentzian manifolds with spacelike and timelike boundary, respectively. While B{\"a}r and Strohmaier [15] showed the Fredholmness of the Dirac operator with these boundary conditions, Drago, Große, and Murro [27] proved the well-posedness of the corresponding initial boundary value problem under certain geometric assumptions. In this thesis, we will follow the footsteps of the latter authors and discuss whether the APS-like conditions for Dirac operators on Lorentzian manifolds with timelike boundary can be replaced by more general conditions such that the associated initial boundary value problems are still wellposed. We consider boundary conditions that are local in time and non-local in the spatial directions. More precisely, we use the spacetime foliation arising from the Cauchy temporal function and split the Dirac operator along this foliation. This gives rise to a family of elliptic operators each acting on spinors of the spin bundle over the corresponding timeslice. The theory of elliptic operators then ensures that we can find families of non-local boundary conditions with respect to this family of operators. Proceeding, we use such a family of boundary conditions to define a Lorentzian boundary condition on the whole timelike boundary. By analyzing the properties of the Lorentzian boundary conditions, we then find sufficient conditions on the family of non-local boundary conditions that lead to the well-posedness of the corresponding Cauchy problems. The well-posedness itself will then be proven by using classical tools including energy estimates and approximation by solutions of the regularized problems. Moreover, we use this theory to construct explicit boundary conditions for the Lorentzian Dirac operator. More precisely, we will discuss two examples of boundary conditions - the analogue of the Atiyah-Patodi-Singer and the chirality conditions, respectively, in our setting. For doing this, we will have a closer look at the theory of non-local boundary conditions for elliptic operators and analyze the requirements on the family of non-local boundary conditions for these specific examples.}, language = {en} } @article{deWiljesPathirajaReich2020, author = {de Wiljes, Jana and Pathiraja, Sahani Darschika and Reich, Sebastian}, title = {Ensemble transform algorithms for nonlinear smoothing problems}, series = {SIAM journal on scientific computing}, volume = {42}, journal = {SIAM journal on scientific computing}, number = {1}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {1064-8275}, doi = {10.1137/19M1239544}, pages = {A87 -- A114}, year = {2020}, abstract = {Several numerical tools designed to overcome the challenges of smoothing in a non-linear and non-Gaussian setting are investigated for a class of particle smoothers. The considered family of smoothers is induced by the class of linear ensemble transform filters which contains classical filters such as the stochastic ensemble Kalman filter, the ensemble square root filter, and the recently introduced nonlinear ensemble transform filter. Further the ensemble transform particle smoother is introduced and particularly highlighted as it is consistent in the particle limit and does not require assumptions with respect to the family of the posterior distribution. The linear update pattern of the considered class of linear ensemble transform smoothers allows one to implement important supplementary techniques such as adaptive spread corrections, hybrid formulations, and localization in order to facilitate their application to complex estimation problems. These additional features are derived and numerically investigated for a sequence of increasingly challenging test problems.}, language = {en} } @phdthesis{Günther2023, author = {G{\"u}nther, Claudia-Susanne}, title = {Das Eigene und das Fremde}, school = {Universit{\"a}t Potsdam}, pages = {245}, year = {2023}, abstract = {Die vorliegende Arbeit stellt eine Untersuchung des Fremdverstehens von Lehrkr{\"a}ften im Mathematikunterricht dar. Mit ‚Fremdverstehen' soll dabei - in Anlehnung an den Soziologen Alfred Sch{\"u}tz - der Prozess bezeichnet werden, in welchem eine Lehrkraft versucht, das Verhalten einer Sch{\"u}lerin oder eines Sch{\"u}lers zu verstehen, indem sie dieses Verhalten auf ein Erleben zur{\"u}ckf{\"u}hrt, das ihm zugrunde gelegen haben k{\"o}nnte. Als ein wesentliches Merkmal des Prozesses stellt Sch{\"u}tz in seiner Theorie des Fremdverstehens heraus, dass das Fremdverstehen eines Menschen immer auch auf seinen eigenen Erlebnissen basiert. Aus diesem Grund wird in der Arbeit ein methodischer Zweischritt vorgenommen: Es werden zun{\"a}chst die mathematikbezogenen Erlebnisse zweier Lehrkr{\"a}fte nachgezeichnet, bevor dann ihr Fremdverstehen in konkreten Situationen im Mathematikunterricht rekonstruiert wird. In der ersten Teiluntersuchung (= der Rekonstruktion eigener Erlebnisse der untersuchten Lehrkr{\"a}fte) erfolgt die Datenerhebung mit Hilfe biographisch-narrativer Interviews, in denen die untersuchten Lehrkr{\"a}fte angeregt werden, ihre mathematikbezogene Lebensgeschichte zu erz{\"a}hlen. Die Analyse dieser Interviews wird im Sinne der rekonstruktiven Fallanalyse vorgenommen. Insgesamt f{\"u}hrt die erste Teiluntersuchung zu textlichen Darstellungen der rekonstruierten mathematikbezogenen Lebensgeschichte der untersuchten Mathematiklehrkr{\"a}fte. In der zweiten Teiluntersuchung (= der Rekonstruktion des Fremdverstehens der untersuchten Lehrkr{\"a}fte) werden dann narrative Interviews gef{\"u}hrt, in denen die untersuchten Lehrkr{\"a}fte von ihrem Fremdverstehen in konkreten Situationen im Mathematikunterricht erz{\"a}hlen. Die Analyse dieser Interviews erfolgt mit Hilfe eines dreischrittigen Analyseverfahrens, welches die Autorin eigens zum Zweck der Rekonstruktion von Fremdverstehen entwickelte. Am Ende dieser zweiten Teiluntersuchung werden sowohl das rekonstruierte Fremdverstehen der Lehrkr{\"a}fte in verschiedenen Unterrichtssituationen dargestellt als auch Strukturen, die sich in ihrem Fremdverstehen abzeichnen. Mit Hilfe einer theoretischen Verallgemeinerung werden schließlich - auf Basis der Ergebnisse der zweiten Teiluntersuchung - Aussagen {\"u}ber f{\"u}nf Merkmale des Fremdverstehens von Lehrkr{\"a}ften im Mathematikunterricht im Allgemeinen gewonnen. Mit diesen Aussagen vermag die Arbeit eine erste Beschreibung davon hervorzubringen, wie sich das Ph{\"a}nomen des Fremdverstehens von Lehrkr{\"a}ften im Mathematikunterricht ausgestalten kann.}, language = {de} } @article{FischerKeller2021, author = {Fischer, Florian and Keller, Matthias}, title = {Riesz decompositions for Schr{\"o}dinger operators on graphs}, series = {Journal of mathematical analysis and applications}, volume = {495}, journal = {Journal of mathematical analysis and applications}, number = {1}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0022-247X}, doi = {10.1016/j.jmaa.2020.124674}, pages = {22}, year = {2021}, abstract = {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.}, language = {en} } @article{ZoellerHainzlTilmannetal.2020, author = {Z{\"o}ller, Gert and Hainzl, Sebastian and Tilmann, Frederik and Woith, Heiko and Dahm, Torsten}, title = {Comment on: Wikelski, Martin; M{\"u}ller, Uschi; Scocco, Paola; Catorci, Andrea; Desinov, Lev V.; Belyaev, Mikhail Y.; Keim, Daniel A.; Pohlmeier, Winfried; Fechteler, Gerhard; Mai, Martin P. : Potential short-term earthquake forecasting by farm animal monitoring. - Ethology. - 126 (2020), 9. - S. 931 - 941. -ISSN 0179-1613. - eISSN 1439-0310. - doi 10.1111/eth.13078}, series = {Ethology}, volume = {127}, journal = {Ethology}, number = {3}, publisher = {Wiley}, address = {Hoboken}, issn = {0179-1613}, doi = {10.1111/eth.13105}, pages = {302 -- 306}, year = {2020}, abstract = {Based on an analysis of continuous monitoring of farm animal behavior in the region of the 2016 M6.6 Norcia earthquake in Italy, Wikelski et al., 2020; (Seismol Res Lett, 89, 2020, 1238) conclude that animal activity can be anticipated with subsequent seismic activity and that this finding might help to design a "short-term earthquake forecasting method." We show that this result is based on an incomplete analysis and misleading interpretations. Applying state-of-the-art methods of statistics, we demonstrate that the proposed anticipatory patterns cannot be distinguished from random patterns, and consequently, the observed anomalies in animal activity do not have any forecasting power.}, language = {en} } @article{Clavier2020, author = {Clavier, Pierre J.}, title = {Double shuffle relations for arborified zeta values}, series = {Journal of algebra}, volume = {543}, journal = {Journal of algebra}, publisher = {Elsevier}, address = {San Diego}, issn = {0021-8693}, doi = {10.1016/j.jalgebra.2019.10.015}, pages = {111 -- 155}, year = {2020}, abstract = {Arborified zeta values are defined as iterated series and integrals using the universal properties of rooted trees. This approach allows to study their convergence domain and to relate them to multiple zeta values. Generalisations to rooted trees of the stuffle and shuffle products are defined and studied. It is further shown that arborified zeta values are algebra morphisms for these new products on trees.}, language = {en} } @article{KayaFreitag2022, author = {Kaya, Adem and Freitag, Melina A.}, title = {Conditioning analysis for discrete Helmholtz problems}, series = {Computers and mathematics with applications : an international journal}, volume = {118}, journal = {Computers and mathematics with applications : an international journal}, publisher = {Elsevier Science}, address = {Amsterdam}, issn = {0898-1221}, doi = {10.1016/j.camwa.2022.05.016}, pages = {171 -- 182}, year = {2022}, abstract = {In this paper, we examine conditioning of the discretization of the Helmholtz problem. Although the discrete Helmholtz problem has been studied from different perspectives, to the best of our knowledge, there is no conditioning analysis for it. We aim to fill this gap in the literature. We propose a novel method in 1D to observe the near-zero eigenvalues of a symmetric indefinite matrix. Standard classification of ill-conditioning based on the matrix condition number is not true for the discrete Helmholtz problem. We relate the ill-conditioning of the discretization of the Helmholtz problem with the condition number of the matrix. We carry out analytical conditioning analysis in 1D and extend our observations to 2D with numerical observations. We examine several discretizations. We find different regions in which the condition number of the problem shows different characteristics. We also explain the general behavior of the solutions in these regions.}, language = {en} }