TY - JOUR
A1 - Ayanbayev, Birzhan
A1 - Klebanov, Ilja
A1 - Li, Han Cheng
A1 - Sullivan, Tim J.
T1 - Gamma-convergence of Onsager-Machlup functionals
BT - I. With applications to maximum a posteriori estimation in Bayesian inverse problems
JF - Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data
N2 - The Bayesian solution to a statistical inverse problem can be summarised by a mode of the posterior distribution, i.e. a maximum a posteriori (MAP) estimator. The MAP estimator essentially coincides with the (regularised) variational solution to the inverse problem, seen as minimisation of the Onsager-Machlup (OM) functional of the posterior measure. An open problem in the stability analysis of inverse problems is to establish a relationship between the convergence properties of solutions obtained by the variational approach and by the Bayesian approach. To address this problem, we propose a general convergence theory for modes that is based on the Gamma-convergence of OM functionals, and apply this theory to Bayesian inverse problems with Gaussian and edge-preserving Besov priors. Part II of this paper considers more general prior distributions.
KW - Bayesian inverse problems
KW - Gamma-convergence
KW - maximum a posteriori
KW - estimation
KW - Onsager-Machlup functional
KW - small ball probabilities;
KW - transition path theory
Y1 - 2021
U6 - https://doi.org/10.1088/1361-6420/ac3f81
SN - 0266-5611
SN - 1361-6420
VL - 38
IS - 2
PB - IOP Publ. Ltd.
CY - Bristol
ER -
TY - JOUR
A1 - Lie, Han Cheng
A1 - Stahn, Martin
A1 - Sullivan, Tim J.
T1 - Randomised one-step time integration methods for deterministic operator differential equations
JF - Calcolo
N2 - 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.
KW - Time integration
KW - Operator differential equations
KW - Randomisation
KW - Uncertainty quantification
Y1 - 2022
U6 - https://doi.org/10.1007/s10092-022-00457-6
SN - 0008-0624
SN - 1126-5434
VL - 59
IS - 1
PB - Springer
CY - Milano
ER -
TY - JOUR
A1 - Malem-Shinitski, Noa
A1 - Ojeda, Cesar
A1 - Opper, Manfred
T1 - Variational bayesian inference for nonlinear hawkes process with gaussian process self-effects
JF - Entropy
N2 - 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.
KW - Bayesian inference
KW - point process
KW - Gaussian process
Y1 - 2022
U6 - https://doi.org/10.3390/e24030356
SN - 1099-4300
VL - 24
IS - 3
PB - MDPI
CY - Basel
ER -
TY - JOUR
A1 - Ayanbayev, Birzhan
A1 - Klebanov, Ilja
A1 - Lie, Han Cheng
A1 - Sullivan, Tim J.
T1 - Gamma-convergence of Onsager-Machlup functionals
BT - II. Infinite product measures on Banach spaces
JF - Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data
N2 - We derive Onsager-Machlup functionals for countable product measures on weighted l(p) subspaces of the sequence space R-N. Each measure in the product is a shifted and scaled copy of a reference probability measure on R that admits a sufficiently regular Lebesgue density. We study the equicoercivity and Gamma-convergence of sequences of Onsager-Machlup functionals associated to convergent sequences of measures within this class. We use these results to establish analogous results for probability measures on separable Banach or Hilbert spaces, including Gaussian, Cauchy, and Besov measures with summability parameter 1 <= p <= 2. Together with part I of this paper, this provides a basis for analysis of the convergence of maximum a posteriori estimators in Bayesian inverse problems and most likely paths in transition path theory.
KW - Bayesian inverse problems
KW - Gamma-convergence
KW - maximum a posteriori
KW - estimation
KW - Onsager-Machlup functional
KW - small ball probabilities
KW - transition path theory
Y1 - 2021
U6 - https://doi.org/10.1088/1361-6420/ac3f82
SN - 0266-5611
SN - 1361-6420
VL - 38
IS - 2
PB - IOP Publ. Ltd.
CY - Bristol
ER -
TY - JOUR
A1 - Lazarides, Rebecca
A1 - Rubach, Charlott
T1 - Instructional characteristics in mathematics classrooms
BT - relationships to achievement goal orientation and student engagement
JF - Mathematics Education Research Journal
N2 - This longitudinal study examined relationships between student-perceived teaching for meaning, support for autonomy, and competence in mathematic classrooms (Time 1), and students’ achievement goal orientations and engagement in mathematics 6 months later (Time 2). We tested whether student-perceived instructional characteristics at Time 1 indirectly related to student engagement at Time 2, via their achievement goal orientations (Time 2), and, whether student gender moderated these relationships. Participants were ninth and tenth graders (55.2% girls) from 46 classrooms in ten secondary schools in Berlin, Germany. Only data from students who participated at both timepoints were included (N = 746 out of total at Time 1 1118; dropout 33.27%). Longitudinal structural equation modeling showed that student-perceived teaching for meaning and support for competence indirectly predicted intrinsic motivation and effort, via students’ mastery goal orientation. These paths were equivalent for girls and boys. The findings are significant for mathematics education, in identifying motivational processes that partly explain the relationships between student-perceived teaching for meaning and competence support and intrinsic motivation and effort in mathematics.
KW - Mathematics classrooms
KW - Instruction
KW - Intrinsicmotivation
KW - Achievement goal orientation
KW - Effort
KW - Gender
Y1 - 2017
U6 - https://doi.org/10.1007/s13394-017-0196-4
SN - 1033-2170
SN - 2211-050X
VL - 29
SP - 201
EP - 217
PB - Springer
CY - Dordrecht
ER -
TY - JOUR
A1 - Eshghi, Nasim
A1 - Mach, Thomas
A1 - Reichel, Lothar
T1 - New matrix function approximations and quadrature rules based on the Arnoldi process
JF - Journal of computational and applied mathematics
N2 - The Arnoldi process can be applied to inexpensively approximate matrix functions of the form f (A)v and matrix functionals of the form v*(f (A))*g(A)v, where A is a large square non-Hermitian matrix, v is a vector, and the superscript * denotes transposition and complex conjugation. Here f and g are analytic functions that are defined in suitable regions in the complex plane. This paper reviews available approximation methods and describes new ones that provide higher accuracy for essentially the same computational effort by exploiting available, but generally not used, moment information. Numerical experiments show that in some cases the modifications of the Arnoldi decompositions proposed can improve the accuracy of v*(f (A))*g(A)v about as much as performing an additional step of the Arnoldi process.
KW - Arnoldi process
KW - Matrix function approximation
KW - Quadrature rule
Y1 - 2021
U6 - https://doi.org/10.1016/j.cam.2021.113442
SN - 0377-0427
SN - 1879-1778
VL - 391
PB - Elsevier
CY - Amsterdam
ER -
TY - JOUR
A1 - Omel'chenko, Oleh
A1 - Laing, Carlo R.
T1 - Collective states in a ring network of theta neurons
JF - Proceedings of the Royal Society of London. Series A, Mathematical, physical and engineering sciences
N2 - We consider a ring network of theta neurons with non-local homogeneous coupling. We analyse the corresponding continuum evolution equation, analytically describing all possible steady states and their stability. By considering a number of different parameter sets, we determine the typical bifurcation scenarios of the network, and put on a rigorous footing some previously observed numerical results.
KW - theta neurons
KW - neural networks
KW - bumps
Y1 - 2022
U6 - https://doi.org/10.1098/rspa.2021.0817
SN - 1364-5021
SN - 1471-2946
VL - 478
IS - 2259
PB - Royal Society
CY - London
ER -
TY - JOUR
A1 - Kolbe, Benedikt Maximilian
A1 - Evans, Myfanwy E.
T1 - Enumerating isotopy classes of tilings guided by the symmetry of triply
JF - Siam journal on applied algebra and geometry
N2 - 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.
KW - isotopic tiling theory
KW - mapping class group
KW - orbifolds
KW - group
KW - presentations
KW - representations of groups as automorphism groups of
KW - algebraic systems
KW - triply periodic minimal surface
KW - Delaney--Dress
KW - tiling theory
KW - hyperbolic tilings
KW - two-dimensional topology
Y1 - 2022
U6 - https://doi.org/10.1137/20M1358943
SN - 2470-6566
VL - 6
IS - 1
SP - 1
EP - 40
PB - Society for Industrial and Applied Mathematics
CY - Philadelphia
ER -
TY - JOUR
A1 - Pasemann, Gregor
A1 - Flemming, Sven
A1 - Alonso, Sergio
A1 - Beta, Carsten
A1 - Stannat, Wilhelm
T1 - Diffusivity estimation for activator-inhibitor models
BT - theory and application to intracellular dynamics of the actin cytoskeleton
JF - Journal of nonlinear science
N2 - A theory for diffusivity estimation for spatially extended activator-inhibitor dynamics modeling the evolution of intracellular signaling networks is developed in the mathematical framework of stochastic reaction-diffusion systems. In order to account for model uncertainties, we extend the results for parameter estimation for semilinear stochastic partial differential equations, as developed in Pasemann and Stannat (Electron J Stat 14(1):547-579, 2020), to the problem of joint estimation of diffusivity and parametrized reaction terms. Our theoretical findings are applied to the estimation of effective diffusivity of signaling components contributing to intracellular dynamics of the actin cytoskeleton in the model organism Dictyostelium discoideum.
KW - Parametric drift estimation
KW - Stochastic reaction– diffusion
KW - systems
KW - Maximum likelihood estimation
KW - Actin cytoskeleton dynamics
Y1 - 2021
U6 - https://doi.org/10.1007/s00332-021-09714-4
SN - 0938-8974
SN - 1432-1467
VL - 31
IS - 3
PB - Springer
CY - New York
ER -
TY - JOUR
A1 - Kürschner, Patrick
A1 - Freitag, Melina A.
T1 - Inexact methods for the low rank solution to large scale Lyapunov equations
JF - BIT : numerical mathematics ; the leading applied mathematics journal for all computational mathematicians
N2 - 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.
KW - Lyapunov equation
KW - alternating direction implicit
KW - rational Krylov
KW - subspaces
KW - low-rank approximations
Y1 - 2019
U6 - https://doi.org/10.1007/s10543-020-00813-4
SN - 0006-3835
SN - 1572-9125
VL - 60
IS - 4
SP - 1221
EP - 1259
PB - Springer
CY - Dordrecht
ER -
TY - JOUR
A1 - Hehn, Jennifer
A1 - Mendez, Daniel
A1 - Uebernickel, Falk
A1 - Brenner, Walter
A1 - Broy, Manfred
T1 - On integrating design thinking for human-centered requirements engineering
JF - IEEE software
N2 - We elaborate on the possibilities and needs to integrate design thinking into requirements engineering, drawing from our research and project experiences. We suggest three approaches for tailoring and integrating design thinking and requirements engineering with complementary synergies and point at open challenges for research and practice.
KW - requirements engineering
KW - prototypes
KW - software
KW - electronic mail
KW - tools
KW - organizations
KW - design thinking
Y1 - 2019
U6 - https://doi.org/10.1109/MS.2019.2957715
SN - 0740-7459
SN - 1937-4194
VL - 37
IS - 2
SP - 25
EP - 31
PB - Inst. of Electr. and Electronics Engineers
CY - Los Alamitos
ER -
TY - JOUR
A1 - Rosenau, Philip
A1 - Pikovskij, Arkadij
T1 - Waves in strongly nonlinear Gardner-like equations on a lattice
JF - Nonlinearity / the Institute of Physics and the London Mathematical Society
N2 - We introduce and study a family of lattice equations which may be viewed either as a strongly nonlinear discrete extension of the Gardner equation, or a non-convex variant of the Lotka-Volterra chain. Their deceptively simple form supports a very rich family of complex solitary patterns. Some of these patterns are also found in the quasi-continuum rendition, but the more intriguing ones, like interlaced pairs of solitary waves, or waves which may reverse their direction either spontaneously or due a collision, are an intrinsic feature of the discrete realm.
KW - nonlinear lattice
KW - solitary wave
KW - Gardner equation
KW - compacton
Y1 - 2021
U6 - https://doi.org/10.1088/1361-6544/ac0f51
SN - 0951-7715
SN - 1361-6544
VL - 34
IS - 8
SP - 5872
EP - 5896
PB - IOP Publ. Ltd.
CY - Bristol
ER -
TY - JOUR
A1 - Cvetković, Nada
A1 - Conrad, Tim
A1 - Lie, Han Cheng
T1 - A convergent discretization method for transition path theory for diffusion processes
JF - Multiscale modeling & simulation : a SIAM interdisciplinary journal
N2 - Transition path theory (TPT) for diffusion processes is a framework for analyzing the transitions of multiscale ergodic diffusion processes between disjoint metastable subsets of state space. Most methods for applying TPT involve the construction of a Markov state model on a discretization of state space that approximates the underlying diffusion process. However, the assumption of Markovianity is difficult to verify in practice, and there are to date no known error bounds or convergence results for these methods. We propose a Monte Carlo method for approximating the forward committor, probability current, and streamlines from TPT for diffusion processes. Our method uses only sample trajectory data and partitions of state space based on Voronoi tessellations. It does not require the construction of a Markovian approximating process. We rigorously prove error bounds for the approximate TPT objects and use these bounds to show convergence to their exact counterparts in the limit of arbitrarily fine discretization. We illustrate some features of our method by application to a process that solves the Smoluchowski equation on a triple-well potential.
KW - ergodic diffusion processes
KW - transition paths
KW - rare events
KW - Monte Carlo
KW - methods
Y1 - 2021
U6 - https://doi.org/10.1137/20M1329354
SN - 1540-3459
SN - 1540-3467
VL - 19
IS - 1
SP - 242
EP - 266
PB - Society for Industrial and Applied Mathematics
CY - Philadelphia
ER -
TY - JOUR
A1 - Perera, Upeksha
A1 - Böckmann, Christine
T1 - Solutions of Direct and Inverse Even-Order Sturm-Liouville Problems Using Magnus Expansion
JF - Mathematics
N2 - 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.
KW - higher-order Sturm–Liouville problems
KW - inverse Sturm–Liouville problems
KW - Magnus expansion
Y1 - 2019
U6 - https://doi.org/10.3390/math7060544
SN - 2227-7390
VL - 7
IS - 6
PB - MDPI
CY - Basel, Schweiz
ER -
TY - JOUR
A1 - Cseh, Ágnes
A1 - Kavitha, Telikepalli
T1 - Popular matchings in complete graphs
JF - Algorithmica : an international journal in computer science
N2 - Our input is a complete graph G on n vertices where each vertex has a strict ranking of all other vertices in G. The goal is to construct a matching in G that is popular. A matching M is popular if M does not lose a head-to-head election against any matching M ': here each vertex casts a vote for the matching in {M,M '} in which it gets a better assignment. Popular matchings need not exist in the given instance G and the popular matching problem is to decide whether one exists or not. The popular matching problem in G is easy to solve for odd n. Surprisingly, the problem becomes NP-complete for even n, as we show here. This is one of the few graph theoretic problems efficiently solvable when n has one parity and NP-complete when n has the other parity.
KW - Popular matching
KW - Complexity
KW - Stable matching
Y1 - 2021
U6 - https://doi.org/10.1007/s00453-020-00791-7
SN - 0178-4617
SN - 1432-0541
VL - 83
IS - 5
SP - 1493
EP - 1523
PB - Springer
CY - New York
ER -
TY - JOUR
A1 - Chang, Der-Chen
A1 - Mahmoudi, Mahdi Hedayat
A1 - Schulze, Bert-Wolfgang
T1 - Volterra operators in the edge-calculus
JF - Analysis and Mathematical Physics
N2 - We study the Volterra property of a class of anisotropic pseudo-differential operators on R x B for a manifold B with edge Y and time-variable t. This exposition belongs to a program for studying parabolicity in such a situation. In the present consideration we establish non-smoothing elements in a subalgebra with anisotropic operator-valued symbols of Mellin type with holomorphic symbols in the complex Mellin covariable from the cone theory, where the covariable t of t extends to symbolswith respect to t to the lower complex v half-plane. The resulting space ofVolterra operators enlarges an approach of Buchholz (Parabolische Pseudodifferentialoperatoren mit operatorwertigen Symbolen. Ph. D. thesis, Universitat Potsdam, 1996) by necessary elements to a new operator algebra containing Volterra parametrices under an appropriate condition of anisotropic ellipticity. Our approach avoids some difficulty in choosing Volterra quantizations in the edge case by generalizing specific achievements from the isotropic edge-calculus, obtained by Seiler (Pseudodifferential calculus on manifolds with non-compact edges, Ph. D. thesis, University of Potsdam, 1997), see also Gil et al. (in: Demuth et al (eds) Mathematical research, vol 100. Akademic Verlag, Berlin, pp 113-137, 1997; Osaka J Math 37: 221-260, 2000).
KW - Volterra operator
KW - Anisotropic pseudo-differential operators
KW - Edge calculus
KW - Operator-valued symbols of Mellin type
Y1 - 2018
U6 - https://doi.org/10.1007/s13324-018-0238-4
SN - 1664-2368
SN - 1664-235X
VL - 8
IS - 4
SP - 551
EP - 570
PB - Springer
CY - Basel
ER -
TY - JOUR
A1 - Flad, Heinz-Jürgen
A1 - Flad-Harutyunyan, Gohar
A1 - Schulze, Bert-Wolfgang
T1 - Explicit Green operators for quantum mechanical Hamiltonians
BT - II. edge-type singularities of the helium atom
JF - Asian-European journal of mathematics : AEJM
N2 - 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.
KW - Singular analysis
KW - Schrodinger equation
KW - many-electron systems
KW - asymptotic properties of eigenfunctions
Y1 - 2020
U6 - https://doi.org/10.1142/S1793557120501223
SN - 1793-5571
SN - 1793-7183
VL - 13
IS - 7
PB - World Scientific
CY - Singapore
ER -
TY - JOUR
A1 - Rungrottheera, Wannarut
A1 - Lyu, Xiaojing
A1 - Schulze, Bert-Wolfgang
T1 - Parameter-dependent edge calculus and corner parametrices
JF - Journal of nonlinear and convex analysis : an international journal
N2 - Let B be a compact manifold with smooth edge of dimension > 0. We study the interplay between parameter-dependent edge algebra algebra on B and operator families belonging to the corner calculus, and we characterize parametrices in the corner case.
KW - Edge calculus
KW - corner parametrices
Y1 - 2018
SN - 1345-4773
SN - 1880-5221
VL - 19
IS - 12
SP - 2021
EP - 2051
PB - Yokohama Publishers
CY - Yokohama
ER -
TY - JOUR
A1 - Stauffer, Maxime
A1 - Mengesha, Isaak
A1 - Seifert, Konrad
A1 - Krawczuk, Igor
A1 - Fischer, Jens
A1 - Serugendo, Giovanna Di Marzo
T1 - A computational turn in policy process studies
BT - coevolving network dynamics of policy change
JF - Complexity
N2 - The past three decades of policy process studies have seen the emergence of a clear intellectual lineage with regard to complexity. Implicitly or explicitly, scholars have employed complexity theory to examine the intricate dynamics of collective action in political contexts. However, the methodological counterparts to complexity theory, such as computational methods, are rarely used and, even if they are, they are often detached from established policy process theory. Building on a critical review of the application of complexity theory to policy process studies, we present and implement a baseline model of policy processes using the logic of coevolving networks. Our model suggests that an actor's influence depends on their environment and on exogenous events facilitating dialogue and consensus-building. Our results validate previous opinion dynamics models and generate novel patterns. Our discussion provides ground for further research and outlines the path for the field to achieve a computational turn.
Y1 - 2022
U6 - https://doi.org/10.1155/2022/8210732
SN - 1076-2787
SN - 1099-0526
VL - 2022
PB - Wiley-Hindawi
CY - London
ER -
TY - JOUR
A1 - Hedayat Mahmoudi, Mahdi
A1 - Schulze, Bert-Wolfgang
T1 - Corner boundary value problems
JF - Asian-European journal of mathematics
N2 - The paper develops some crucial steps in extending the first-order cone or edge calculus to higher singularity orders. We focus here on order 2, but the ideas are motivated by an iterative approach for higher singularities.
KW - Mellin operators
KW - Mellin oscillatory integrals
KW - exit calculus
KW - weighted Sobolev spaces
Y1 - 2016
U6 - https://doi.org/10.1142/S1793557117500541
SN - 1793-5571
SN - 1793-7183
VL - 10
IS - 1
PB - World Scientific
CY - Singapore
ER -
TY - JOUR
A1 - Cseh, Agnes
A1 - Faenza, Yuri
A1 - Kavitha, Telikepalli
A1 - Powers, Vladlena
T1 - Understanding popular matchings via stable matchings
JF - SIAM journal on discrete mathematics
N2 - An instance of the marriage problem is given by a graph G = (A boolean OR B, E), together with, for each vertex of G, a strict preference order over its neighbors. A matching M of G is popular in the marriage instance if M does not lose a head-to-head election against any matching where vertices are voters. Every stable matching is a min-size popular matching; another subclass of popular matchings that always exists and can be easily computed is the set of dominant matchings. A popular matching M is dominant if M wins the head-to-head election against any larger matching. Thus, every dominant matching is a max-size popular matching, and it is known that the set of dominant matchings is the linear image of the set of stable matchings in an auxiliary graph. Results from the literature seem to suggest that stable and dominant matchings behave, from a complexity theory point of view, in a very similar manner within the class of popular matchings. The goal of this paper is to show that there are instead differences in the tractability of stable and dominant matchings and to investigate further their importance for popular matchings. First, we show that it is easy to check if all popular matchings are also stable; however, it is co-NP hard to check if all popular matchings are also dominant. Second, we show how some new and recent hardness results on popular matching problems can be deduced from the NP-hardness of certain problems on stable matchings, also studied in this paper, thus showing that stable matchings can be employed to show not only positive results on popular matchings (as is known) but also most negative ones. Problems for which we show new hardness results include finding a min-size (resp., max-size) popular matching that is not stable (resp., dominant). A known result for which we give a new and simple proof is the NP-hardness of finding a popular matching when G is nonbipartite.
KW - popular matching
KW - stable matching
KW - complexity
KW - dominant matching
Y1 - 2022
U6 - https://doi.org/10.1137/19M124770X
SN - 0895-4801
SN - 1095-7146
VL - 36
IS - 1
SP - 188
EP - 213
PB - Society for Industrial and Applied Mathematics
CY - Philadelphia
ER -
TY - JOUR
A1 - Chang, Der-Chen
A1 - Khalil, Sara
A1 - Schulze, Bert-Wolfgang
T1 - Analysis on regular corner spaces
JF - The journal of geometric analysis
N2 - We establish a new approach of treating elliptic boundary value problems (BVPs) on manifolds with boundary and regular corners, up to singularity order 2. Ellipticity and parametrices are obtained in terms of symbols taking values in algebras of BVPs on manifolds of corresponding lower singularity orders. Those refer to Boutet de Monvel's calculus of operators with the transmission property, see Boutet de Monvel (Acta Math 126:11-51, 1971) for the case of smooth boundary. On corner configuration operators act in spaces with multiple weights. We mainly study the case of upper left entries in the respective 2 x 2 operator block-matrices of such a calculus. Green operators in the sense of Boutet de Monvel (Acta Math 126:11-51, 1971) analogously appear in singular cases, and they are complemented by contributions of Mellin type. We formulate a result on ellipticity and the Fredholm property in weighted corner spaces, with parametrices of analogous kind.
KW - Boutet de Monvel's calculus
KW - Pseudo-differential operators
KW - Singular cones
KW - Mellin symbols with values in the edge calculus
KW - Parametrices of elliptic operators
KW - Kegel space
Y1 - 2021
U6 - https://doi.org/10.1007/s12220-021-00614-3
SN - 1050-6926
SN - 1559-002X
VL - 31
IS - 9
SP - 9199
EP - 9240
PB - Springer
CY - New York
ER -
TY - JOUR
A1 - Eckert, Silvia
A1 - Herden, Jasmin
A1 - Stift, Marc
A1 - Durka, Walter
A1 - Kleunen, Mark van
A1 - Joshi, Jasmin Radha
T1 - Traces of genetic but not epigenetic adaptation in the invasive goldenrod Solidago canadensis despite the absence of population structure
JF - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe
N2 - Biological invasions may result from multiple introductions, which might compensate for reduced gene pools caused by bottleneck events, but could also dilute adaptive processes. A previous common-garden experiment showed heritable latitudinal clines in fitness-related traits in the invasive goldenrod Solidago canadensis in Central Europe. These latitudinal clines remained stable even in plants chemically treated with zebularine to reduce epigenetic variation. However, despite the heritability of traits investigated, genetic isolation-by-distance was non-significant. Utilizing the same specimens, we applied a molecular analysis of (epi)genetic differentiation with standard and methylation-sensitive (MSAP) AFLPs. We tested whether this variation was spatially structured among populations and whether zebularine had altered epigenetic variation. Additionally, we used genome scans to mine for putative outlier loci susceptible to selection processes in the invaded range. Despite the absence of isolation-by-distance, we found spatial genetic neighborhoods among populations and two AFLP clusters differentiating northern and southern Solidago populations. Genetic and epigenetic diversity were significantly correlated, but not linked to phenotypic variation. Hence, no spatial epigenetic patterns were detected along the latitudinal gradient sampled. Applying genome-scan approaches (BAYESCAN, BAYESCENV, RDA, and LFMM), we found 51 genetic and epigenetic loci putatively responding to selection. One of these genetic loci was significantly more frequent in populations at the northern range. Also, one epigenetic locus was more frequent in populations in the southern range, but this pattern was lost under zebularine treatment. Our results point to some genetic, but not epigenetic adaptation processes along a large-scale latitudinal gradient of S. canadensis in its invasive range.
T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1271
KW - AFLP
KW - MSAP
KW - cytosine methylation
KW - spatial autocorrelation
KW - genome scan
Y1 - 2022
U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-566758
SN - 1866-8372
SP - 1
EP - 17
PB - Universitätsverlag Potsdam
CY - Potsdam
ER -
TY - JOUR
A1 - Dube, Jonas
A1 - Böckmann, Christine
A1 - Ritter, Christoph
T1 - Lidar-Derived Aerosol Properties from Ny-Ålesund, Svalbard during the MOSAiC Spring 2020
JF - Remote sensing / Molecular Diversity Preservation International (MDPI)
N2 - 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.
KW - aerosol
KW - Arctic haze
KW - lidar
KW - microphysical properties
KW - backtrajectories;
KW - Ny-Alesund
KW - Svalbard
KW - MOSAiC
KW - aerosol-boundary layer interactions
Y1 - 2022
U6 - https://doi.org/10.3390/rs14112578
SN - 2072-4292
VL - 14
IS - 11
PB - MDPI
CY - Basel
ER -
TY - JOUR
A1 - Pornsawad, Pornsarp
A1 - Sapsakul, Nantawan
A1 - Böckmann, Christine
T1 - A modified asymptotical regularization of nonlinear ill-posed problems
JF - Mathematics
N2 - 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.
KW - nonlinear operator
KW - regularization
KW - discrepancy principle
KW - asymptotic method
KW - optimal rate
Y1 - 2019
U6 - https://doi.org/10.3390/math7050419
SN - 2227-7390
VL - 7
PB - MDPI
CY - Basel, Schweiz
ET - 5
ER -
TY - JOUR
A1 - Pornsawad, Pornsarp
A1 - Böckmann, Christine
A1 - Panitsupakamon, Wannapa
T1 - The Levenberg–Marquardt regularization for the backward heat equation with fractional derivative
JF - Electronic transactions on numerical analysis - ETNA
N2 - The backward heat problem with time-fractional derivative in Caputo's sense is studied. The inverse problem is severely ill-posed in the case when the fractional order is close to unity. A Levenberg-Marquardt method with a new a posteriori stopping rule is investigated. We show that optimal order can be obtained for the proposed method under a Hölder-type source condition. Numerical examples for one and two dimensions are provided.
KW - ill-posed problems
KW - time-fractional derivative
KW - backward heat problem
KW - Levenberg-Marquardt method
KW - a posteriori stopping rule
KW - optimal order
Y1 - 2022
SN - 978-3-7001-8258-0
U6 - https://doi.org/10.1553/etna_vol57s67
SN - 1068-9613
VL - 57
SP - 67
EP - 79
PB - Kent State University
CY - Kent
ER -
TY - JOUR
A1 - Wormell, Caroline L.
A1 - Reich, Sebastian
T1 - Spectral convergence of diffusion maps
BT - Improved error bounds and an alternative normalization
JF - SIAM journal on numerical analysis / Society for Industrial and Applied Mathematics
N2 - Diffusion maps is a manifold learning algorithm widely used for dimensionality reduction. Using a sample from a distribution, it approximates the eigenvalues and eigenfunctions of associated Laplace-Beltrami operators. Theoretical bounds on the approximation error are, however, generally much weaker than the rates that are seen in practice. This paper uses new approaches to improve the error bounds in the model case where the distribution is supported on a hypertorus. For the data sampling (variance) component of the error we make spatially localized compact embedding estimates on certain Hardy spaces; we study the deterministic (bias) component as a perturbation of the Laplace-Beltrami operator's associated PDE and apply relevant spectral stability results. Using these approaches, we match long-standing pointwise error bounds for both the spectral data and the norm convergence of the operator discretization. We also introduce an alternative normalization for diffusion maps based on Sinkhorn weights. This normalization approximates a Langevin diffusion on the sample and yields a symmetric operator approximation. We prove that it has better convergence compared with the standard normalization on flat domains, and we present a highly efficient rigorous algorithm to compute the Sinkhorn weights.
KW - diffusion maps
KW - graph Laplacian
KW - Sinkhorn problem
KW - kernel methods
Y1 - 2021
U6 - https://doi.org/10.1137/20M1344093
SN - 0036-1429
SN - 1095-7170
VL - 59
IS - 3
SP - 1687
EP - 1734
PB - Society for Industrial and Applied Mathematics
CY - Philadelphia
ER -
TY - JOUR
A1 - Keller, Matthias
A1 - Liu, Shiping
A1 - Peyerimhoff, Norbert
T1 - A note on eigenvalue bounds for non-compact manifolds
JF - Mathematische Nachrichten
N2 - In this article we prove upper bounds for the Laplace eigenvalues lambda(k) below the essential spectrum for strictly negatively curved Cartan-Hadamard manifolds. Our bound is given in terms of k(2) and specific geometric data of the manifold. This applies also to the particular case of non-compact manifolds whose sectional curvature tends to -infinity, where no essential spectrum is present due to a theorem of Donnelly/Li. The result stands in clear contrast to Laplacians on graphs where such a bound fails to be true in general.
KW - Cheeger inequality
KW - eigenvalues
KW - Laplacian
KW - negative curvature
KW - Riemannian manifold
Y1 - 2021
U6 - https://doi.org/10.1002/mana.201900209
SN - 0025-584X
SN - 1522-2616
VL - 294
IS - 6
SP - 1134
EP - 1139
PB - Wiley-VCH
CY - Weinheim
ER -
TY - JOUR
A1 - Peng, Junhao
A1 - Sandev, Trifce
A1 - Kocarev, Ljupco
T1 - First encounters on Bethe lattices and Cayley trees
JF - Communications in nonlinear science & numerical simulation
N2 - In this work we consider the first encounter problems between a fixed and/or mobile target A and a moving trap B on Bethe lattices and Cayley trees. The survival probabilities (SPs) of the target A on the both kinds of structures are considered analytically and compared. On Bethe lattices, the results show that the fixed target will still prolong its survival time, whereas, on Cayley trees, there are some initial positions where the target should move to prolong its survival time. The mean first encounter time (MFET) for mobile target A is evaluated numerically and compared with the mean first passage time (MFPT) for the fixed target A. Different initial settings are addressed and clear boundaries are obtained. These findings are helpful for optimizing the strategy to prolong the survival time of the target or to speed up the search process on Cayley trees, in relation to the target's movement and the initial position configuration of the two walkers. We also present a new method, which uses a small amount of memory, for simulating random walks on Cayley trees. (C) 2020 Elsevier B.V. All rights reserved.
KW - Random walks
KW - Survival probability
KW - Mean first encounter time
KW - Bethe
KW - lattices
KW - Cayley trees
Y1 - 2021
U6 - https://doi.org/10.1016/j.cnsns.2020.105594
SN - 1007-5704
SN - 1878-7274
VL - 95
PB - Elsevier
CY - Amsterdam
ER -
TY - JOUR
A1 - Andjelkovic, Marko
A1 - Simevski, Aleksandar
A1 - Chen, Junchao
A1 - Schrape, Oliver
A1 - Stamenkovic, Zoran
A1 - Krstić, Miloš
A1 - Ilic, Stefan
A1 - Ristic, Goran
A1 - Jaksic, Aleksandar
A1 - Vasovic, Nikola
A1 - Duane, Russell
A1 - Palma, Alberto J.
A1 - Lallena, Antonio M.
A1 - Carvajal, Miguel A.
T1 - A design concept for radiation hardened RADFET readout system for space applications
JF - Microprocessors and microsystems
N2 - Instruments for measuring the absorbed dose and dose rate under radiation exposure, known as radiation dosimeters, are indispensable in space missions. They are composed of radiation sensors that generate current or voltage response when exposed to ionizing radiation, and processing electronics for computing the absorbed dose and dose rate. Among a wide range of existing radiation sensors, the Radiation Sensitive Field Effect Transistors (RADFETs) have unique advantages for absorbed dose measurement, and a proven record of successful exploitation in space missions. It has been shown that the RADFETs may be also used for the dose rate monitoring. In that regard, we propose a unique design concept that supports the simultaneous operation of a single RADFET as absorbed dose and dose rate monitor. This enables to reduce the cost of implementation, since the need for other types of radiation sensors can be minimized or eliminated. For processing the RADFET's response we propose a readout system composed of analog signal conditioner (ASC) and a self-adaptive multiprocessing system-on-chip (MPSoC). The soft error rate of MPSoC is monitored in real time with embedded sensors, allowing the autonomous switching between three operating modes (high-performance, de-stress and fault-tolerant), according to the application requirements and radiation conditions.
KW - RADFET
KW - Radiation hardness
KW - Absorbed dose
KW - Dose rate
KW - Self-adaptive MPSoC
Y1 - 2022
U6 - https://doi.org/10.1016/j.micpro.2022.104486
SN - 0141-9331
SN - 1872-9436
VL - 90
PB - Elsevier
CY - Amsterdam
ER -
TY - JOUR
A1 - Chang, Der-Chen
A1 - Schulze, Bert-Wolfgang
T1 - Corner spaces and Mellin quantization
JF - Journal of nonlinear and convex analysis : an international journal
N2 - Manifolds with corners in the present investigation are non-smooth configurations - specific stratified spaces - with an incomplete metric such as cones, manifolds with edges, or corners of piecewise smooth domains in Euclidean space. We focus here on operators on such "corner manifolds" of singularity order <= 2, acting in weighted corner Sobolev spaces. The corresponding corner degenerate pseudo-differential operators are formulated via Mellin quantizations, and they also make sense on infinite singular cones.
KW - Mellin quantizations
KW - operator-valued symbols
KW - weighted edge and corner spaces
Y1 - 2018
SN - 1345-4773
SN - 1880-5221
VL - 19
IS - 2
SP - 179
EP - 195
PB - Yokohama Publishers
CY - Yokohama
ER -
TY - JOUR
A1 - Chang, Der-Chen
A1 - Schulze, Bert-Wolfgang
T1 - Ellipticity on spaces with higher singularities
JF - Science China Mathematics
N2 - We study corner-degenerate pseudo-differential operators of any singularity order and develop ellipticity based on the principal symbolic hierarchy, associated with the stratification of the underlying space. We construct parametrices within the calculus and discuss the aspect of additional trace and potential conditions along lower-dimensional strata.
KW - pseudo-differential operators
KW - operator-valued symbols
KW - Fourier and Mellin transforms
Y1 - 2017
U6 - https://doi.org/10.1007/s11425-016-0519-9
SN - 1674-7283
SN - 1869-1862
VL - 60
IS - 11
SP - 2053
EP - 2076
PB - Science China Press
CY - Beijing
ER -
TY - JOUR
A1 - Khalil, Sara
A1 - Schulze, Bert-Wolfgang
T1 - Boundary problems on a manifold with edge
JF - Asian-European Journal of Mathematics
N2 - We establish a calculus of boundary value problems (BVPs) on a manifold N with boundary and edge, based on Boutet de Monvel’s theory of BVPs in the case of a smooth boundary and on the edge calculus, where in the present case the model cone has a base which is a compact manifold with boundary. The corresponding calculus with boundary and edge is a unification of both structures and controls different operator-valued symbolic structures, in order to obtain ellipticity and parametrices.
KW - manifolds with edge and boundary
KW - distribution with asymptotics
KW - ellipticity
KW - Fredholm property
Y1 - 2017
U6 - https://doi.org/10.1142/S1793557117500875
SN - 1793-5571
SN - 1793-7183
VL - 10
IS - 2
PB - World Scientific
CY - Singapore
ER -
TY - JOUR
A1 - Chang, Der-Chen
A1 - Hedayat Mahmoudi, Mahdi
A1 - Schulze, Bert-Wolfgang
T1 - Singular degenerate operators
JF - Applicable analysis : an international journal
N2 - We outline some simplified and more general method for constructing parametrices on higher singular spaces. We also outline basic ideas on operators on manifolds with conical or edge singularities.
KW - Operators on singular cones
KW - Mellin symbols with values in the edge calculus
KW - parametrices of elliptic operators
Y1 - 2017
U6 - https://doi.org/10.1080/00036811.2017.1336546
SN - 0003-6811
SN - 1563-504X
VL - 96
IS - 14
SP - 2434
EP - 2456
PB - Routledge, Taylor & Francis Group
CY - Abingdon
ER -
TY - JOUR
A1 - Harutyunyan, Gohar
A1 - Schulze, Bert-Wolfgang
T1 - The Zaremba problem with singular interfaces as a corner boundary value problem
JF - Potential analysis : an international journal devoted to the interactions between potential theory, probability theory, geometry and functional analysis
N2 - We study mixed boundary value problems for an elliptic operator A on a manifold X with boundary Y, i.e., Au = f in int X, T (+/-) u = g(+/-) on int Y+/-, where Y is subdivided into subsets Y+/- with an interface Z and boundary conditions T+/- on Y+/- that are Shapiro-Lopatinskij elliptic up to Z from the respective sides. We assume that Z subset of Y is a manifold with conical singularity v. As an example we consider the Zaremba problem, where A is the Laplacian and T- Dirichlet, T+ Neumann conditions. The problem is treated as a corner boundary value problem near v which is the new point and the main difficulty in this paper. Outside v the problem belongs to the edge calculus as is shown in Bull. Sci. Math. ( to appear). With a mixed problem we associate Fredholm operators in weighted corner Sobolev spaces with double weights, under suitable edge conditions along Z {v} of trace and potential type. We construct parametrices within the calculus and establish the regularity of solutions.
KW - Zaremba problem
KW - corner Sobolev spaces with double weights
KW - pseudo-differential boundary value problems
Y1 - 2006
U6 - https://doi.org/10.1007/s11118-006-9020-6
SN - 0926-2601
VL - 25
SP - 327
EP - 369
PB - Springer
CY - Dordrecht
ER -
TY - JOUR
A1 - Hedayat Mahmoudi, Mahdi
A1 - Schulze, Bert-Wolfgang
T1 - A new approach to the second order edge calculus
JF - Journal of pseudo-differential operators and applications
N2 - We establish essential steps of an iterative approach to operator algebras, ellipticity and Fredholm property on stratified spaces with singularities of second order. We cover, in particular, corner-degenerate differential operators. Our constructions are focused on the case where no additional conditions of trace and potential type are posed, but this case works well and will be considered in a forthcoming paper as a conclusion of the present calculus.
KW - Operators on singular manifolds
KW - Mellin transform
KW - Stratified spaces
KW - Ellipticity and parametrices
Y1 - 2018
U6 - https://doi.org/10.1007/s11868-017-0191-2
SN - 1662-9981
SN - 1662-999X
VL - 9
IS - 2
SP - 265
EP - 300
PB - Springer
CY - Basel
ER -
TY - JOUR
A1 - Khalil, Sara
A1 - Schulze, Bert-Wolfgang
T1 - Calculus on a Manifold with Edge and Boundary
JF - Complex analysis and operator theory
N2 - We study elements of the calculus of boundary value problems in a variant of Boutet de Monvel’s algebra (Acta Math 126:11–51, 1971) on a manifold N with edge and boundary. If the boundary is empty then the approach corresponds to Schulze (Symposium on partial differential equations (Holzhau, 1988), BSB Teubner, Leipzig, 1989) and other papers from the subsequent development. For non-trivial boundary we study Mellin-edge quantizations and compositions within the structure in terms a new Mellin-edge quantization, compared with a more traditional technique. Similar structures in the closed case have been studied in Gil et al.
KW - algebra
KW - Mellin quantization
Y1 - 2019
U6 - https://doi.org/10.1007/s11785-018-0800-y
SN - 1661-8254
SN - 1661-8262
VL - 13
IS - 6
SP - 2627
EP - 2670
PB - Springer
CY - Basel
ER -
TY - JOUR
A1 - Schulze, Bert-Wolfgang
A1 - Seiler, Jörg
T1 - Elliptic complexes on manifolds with boundary
JF - The journal of geometric analysis
N2 - We show that elliptic complexes of (pseudo) differential operators on smooth compact manifolds with boundary can always be complemented to a Fredholm problem by boundary conditions involving global pseudodifferential projections on the boundary (similarly as the spectral boundary conditions of Atiyah, Patodi, and Singer for a single operator). We prove that boundary conditions without projections can be chosen if, and only if, the topological Atiyah-Bott obstruction vanishes. These results make use of a Fredholm theory for complexes of operators in algebras of generalized pseudodifferential operators of Toeplitz type which we also develop in the present paper.
KW - Elliptic complexes
KW - Manifolds with boundary
KW - Atiyah-Bott obstruction
KW - Toeplitz-type pseudodifferential operators
Y1 - 2018
U6 - https://doi.org/10.1007/s12220-018-0014-6
SN - 1050-6926
SN - 1559-002X
VL - 29
IS - 1
SP - 656
EP - 706
PB - Springer
CY - New York
ER -
TY - JOUR
A1 - Zöller, Gert
A1 - Hainzl, Sebastian
A1 - Tilmann, Frederik
A1 - Woith, Heiko
A1 - Dahm, Torsten
T1 - Comment on: Wikelski, Martin; Mü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
JF - Ethology
N2 - 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.
KW - animal behavior
KW - earthquake precursor
KW - error diagram
KW - prediction
KW - randomness
KW - statistics
Y1 - 2020
U6 - https://doi.org/10.1111/eth.13105
SN - 0179-1613
SN - 1439-0310
VL - 127
IS - 3
SP - 302
EP - 306
PB - Wiley
CY - Hoboken
ER -
TY - JOUR
A1 - Pohle, Jennifer
A1 - Adam, Timo
A1 - Beumer, Larissa
T1 - Flexible estimation of the state dwell-time distribution in hidden semi-Markov models
JF - Computational statistics & data analysis
N2 - 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.
KW - Penalized likelihood
KW - Smoothing
KW - Time series
KW - Animal movement modeling
Y1 - 2022
U6 - https://doi.org/10.1016/j.csda.2022.107479
SN - 0167-9473
SN - 1872-7352
VL - 172
PB - Elsevier
CY - Amsterdam
ER -
TY - JOUR
A1 - Biskaborn, Boris
A1 - Smith, Sharon L.
A1 - Noetzli, Jeannette
A1 - Matthes, Heidrun
A1 - Vieira, Goncalo
A1 - Streletskiy, Dmitry A.
A1 - Schoeneich, Philippe
A1 - Romanovsky, Vladimir E.
A1 - Lewkowicz, Antoni G.
A1 - Abramov, Andrey
A1 - Allard, Michel
A1 - Boike, Julia
A1 - Cable, William L.
A1 - Christiansen, Hanne H.
A1 - Delaloye, Reynald
A1 - Diekmann, Bernhard
A1 - Drozdov, Dmitry
A1 - Etzelmueller, Bernd
A1 - Grosse, Guido
A1 - Guglielmin, Mauro
A1 - Ingeman-Nielsen, Thomas
A1 - Isaksen, Ketil
A1 - Ishikawa, Mamoru
A1 - Johansson, Margareta
A1 - Johannsson, Halldor
A1 - Joo, Anseok
A1 - Kaverin, Dmitry
A1 - Kholodov, Alexander
A1 - Konstantinov, Pavel
A1 - Kroeger, Tim
A1 - Lambiel, Christophe
A1 - Lanckman, Jean-Pierre
A1 - Luo, Dongliang
A1 - Malkova, Galina
A1 - Meiklejohn, Ian
A1 - Moskalenko, Natalia
A1 - Oliva, Marc
A1 - Phillips, Marcia
A1 - Ramos, Miguel
A1 - Sannel, A. Britta K.
A1 - Sergeev, Dmitrii
A1 - Seybold, Cathy
A1 - Skryabin, Pavel
A1 - Vasiliev, Alexander
A1 - Wu, Qingbai
A1 - Yoshikawa, Kenji
A1 - Zheleznyak, Mikhail
A1 - Lantuit, Hugues
T1 - Permafrost is warming at a global scale
JF - Nature Communications
N2 - Permafrost warming has the potential to amplify global climate change, because when frozen sediments thaw it unlocks soil organic carbon. Yet to date, no globally consistent assessment of permafrost temperature change has been compiled. Here we use a global data set of permafrost temperature time series from the Global Terrestrial Network for Permafrost to evaluate temperature change across permafrost regions for the period since the International Polar Year (2007-2009). During the reference decade between 2007 and 2016, ground temperature near the depth of zero annual amplitude in the continuous permafrost zone increased by 0.39 +/- 0.15 degrees C. Over the same period, discontinuous permafrost warmed by 0.20 +/- 0.10 degrees C. Permafrost in mountains warmed by 0.19 +/- 0.05 degrees C and in Antarctica by 0.37 +/- 0.10 degrees C. Globally, permafrost temperature increased by 0.29 +/- 0.12 degrees C. The observed trend follows the Arctic amplification of air temperature increase in the Northern Hemisphere. In the discontinuous zone, however, ground warming occurred due to increased snow thickness while air temperature remained statistically unchanged.
Y1 - 2019
U6 - https://doi.org/10.1038/s41467-018-08240-4
SN - 2041-1723
VL - 10
PB - Nature Publ. Group
CY - London
ER -
TY - JOUR
A1 - Mueller-Schoell, Anna
A1 - Groenland, Stefanie L.
A1 - Scherf-Clavel, Oliver
A1 - van Dyk, Madele
A1 - Huisinga, Wilhelm
A1 - Michelet, Robin
A1 - Jaehde, Ulrich
A1 - Steeghs, Neeltje
A1 - Huitema, Alwin D. R.
A1 - Kloft, Charlotte
T1 - Therapeutic drug monitoring of oral targeted antineoplastic drugs
JF - European journal of clinical pharmacology
N2 - Purpose This review provides an overview of the current challenges in oral targeted antineoplastic drug (OAD) dosing and outlines the unexploited value of therapeutic drug monitoring (TDM). Factors influencing the pharmacokinetic exposure in OAD therapy are depicted together with an overview of different TDM approaches. Finally, current evidence for TDM for all approved OADs is reviewed. Methods A comprehensive literature search (covering literature published until April 2020), including primary and secondary scientific literature on pharmacokinetics and dose individualisation strategies for OADs, together with US FDA Clinical Pharmacology and Biopharmaceutics Reviews and the Committee for Medicinal Products for Human Use European Public Assessment Reports was conducted. Results OADs are highly potent drugs, which have substantially changed treatment options for cancer patients. Nevertheless, high pharmacokinetic variability and low treatment adherence are risk factors for treatment failure. TDM is a powerful tool to individualise drug dosing, ensure drug concentrations within the therapeutic window and increase treatment success rates. After reviewing the literature for 71 approved OADs, we show that exposure-response and/or exposure-toxicity relationships have been established for the majority. Moreover, TDM has been proven to be feasible for individualised dosing of abiraterone, everolimus, imatinib, pazopanib, sunitinib and tamoxifen in prospective studies. There is a lack of experience in how to best implement TDM as part of clinical routine in OAD cancer therapy. Conclusion Sub-therapeutic concentrations and severe adverse events are current challenges in OAD treatment, which can both be addressed by the application of TDM-guided dosing, ensuring concentrations within the therapeutic window.
KW - targeted antineoplastic drugs
KW - tyrosine kinase inhibitors
KW - therapeutic
KW - drug monitoring
KW - oral anticancer drugs
KW - personalised medicine
Y1 - 2020
U6 - https://doi.org/10.1007/s00228-020-03014-8
SN - 0031-6970
SN - 1432-1041
VL - 77
IS - 4
SP - 441
EP - 464
PB - Springer
CY - Heidelberg
ER -
TY - JOUR
A1 - Grisic, Ana-Marija
A1 - Eser, Alexander
A1 - Huisinga, Wilhelm
A1 - Reinisch, Walter
A1 - Kloft, Charlotte
T1 - Quantitative relationship between infliximab exposure and inhibition of C-reactive protein synthesis to support inflammatory bowel disease management
JF - British journal of clinical pharmacology
N2 - Aim Quantitative and kinetic insights into the drug exposure-disease response relationship might enhance our knowledge on loss of response and support more effective monitoring of inflammatory activity by biomarkers in patients with inflammatory bowel disease (IBD) treated with infliximab (IFX). This study aimed to derive recommendations for dose adjustment and treatment optimisation based on mechanistic characterisation of the relationship between IFX serum concentration and C-reactive protein (CRP) concentration.
Methods Data from an investigator-initiated trial included 121 patients with IBD during IFX maintenance treatment. Serum concentrations of IFX, antidrug antibodies (ADA), CRP, and disease-related covariates were determined at the mid-term and end of a dosing interval. Data were analysed using a pharmacometric nonlinear mixed-effects modelling approach. An IFX exposure-CRP model was generated and applied to evaluate dosing regimens to achieve CRP remission.
Results The generated quantitative model showed that IFX has the potential to inhibit up to 72% (9% relative standard error [RSE]) of CRP synthesis in a patient. IFX concentration leading to 90% of the maximum CRP synthesis inhibition was 18.4 mu g/mL (43% RSE). Presence of ADA was the most influential factor on IFX exposure. With standard dosing strategy, >= 55% of ADA+ patients experienced CRP nonremission. Shortening the dosing interval and co-therapy with immunomodulators were found to be the most beneficial strategies to maintain CRP remission.
Conclusions With the generated model we could for the first time establish a robust relationship between IFX exposure and CRP synthesis inhibition, which could be utilised for treatment optimisation in IBD patients.
KW - C‐ reactive protein remission
KW - inflammatory bowel disease
KW - infliximab dosing
Y1 - 2020
U6 - https://doi.org/10.1111/bcp.14648
SN - 0306-5251
SN - 1365-2125
VL - 87
IS - 5
SP - 2374
EP - 2384
PB - Wiley
CY - Hoboken
ER -
TY - JOUR
A1 - Kluwe, Franziska
A1 - Michelet, Robin
A1 - Müller-Schöll, Anna
A1 - Maier, Corinna
A1 - Klopp-Schulze, Lena
A1 - van Dyk, Madele
A1 - Mikus, Gerd
A1 - Huisinga, Wilhelm
A1 - Kloft, Charlotte
T1 - Perspectives on model-informed precision dosing in the digital health era
BT - challenges, opportunities, and recommendations
JF - Clinical pharmacology & therapeutics
Y1 - 2020
U6 - https://doi.org/10.1002/cpt.2049
SN - 0009-9236
SN - 1532-6535
VL - 109
IS - 1
SP - 29
EP - 36
PB - Wiley
CY - Hoboken
ER -
TY - JOUR
A1 - Nassar, Yomna M.
A1 - Hohmann, Nicolas
A1 - Michelet, Robin
A1 - Gottwalt, Katharina
A1 - Meid, Andreas D.
A1 - Burhenne, Jürgen
A1 - Huisinga, Wilhelm
A1 - Haefeli, Walter E.
A1 - Mikus, Gerd
A1 - Kloft, Charlotte
T1 - Quantification of the Time Course of CYP3A Inhibition, Activation, and Induction Using a Population Pharmacokinetic Model of Microdosed Midazolam Continuous Infusion
JF - Clinical Pharmacokinetics
N2 - 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.
Y1 - 2022
U6 - https://doi.org/10.1007/s40262-022-01175-6
SN - 0312-5963
SN - 1179-1926
VL - 61
IS - 11
SP - 1595
EP - 1607
PB - Springer
CY - Northcote
ER -
TY - JOUR
A1 - Wicha, Sebastian G.
A1 - Huisinga, Wilhelm
A1 - Kloft, Charlotte
T1 - Translational pharmacometric evaluation of typical antibiotic broad-spectrum combination therapies against staphylococcus aureus exploiting in vitro information
JF - CPT: pharmacometrics & systems pharmacology
N2 - Broad-spectrum antibiotic combination therapy is frequently applied due to increasing resistance development of infective pathogens. The objective of the present study was to evaluate two common empiric broad-spectrum combination therapies consisting of either linezolid (LZD) or vancomycin (VAN) combined with meropenem (MER) against Staphylococcus aureus (S. aureus) as the most frequent causative pathogen of severe infections. A semimechanistic pharmacokinetic-pharmacodynamic (PK-PD) model mimicking a simplified bacterial life-cycle of S. aureus was developed upon time-kill curve data to describe the effects of LZD, VAN, and MER alone and in dual combinations. The PK-PD model was successfully (i) evaluated with external data from two clinical S. aureus isolates and further drug combinations and (ii) challenged to predict common clinical PK-PD indices and breakpoints. Finally, clinical trial simulations were performed that revealed that the combination of VAN-MER might be favorable over LZD-MER due to an unfavorable antagonistic interaction between LZD and MER.
Y1 - 2017
U6 - https://doi.org/10.1002/psp4.12197
SN - 2163-8306
VL - 6
SP - 512
EP - 522
PB - Wiley
CY - Hoboken
ER -
TY - JOUR
A1 - Edlund, Helena
A1 - Grisic, Ana-Marija
A1 - Steenholdt, Casper
A1 - Ainsworth, Mark Andrew
A1 - Brynskov, Torn
A1 - Huisinga, Wilhelm
A1 - Kloft, Charlotte
T1 - Absence of Relationship Between Crohn's Disease Activity Index or C-Reactive Protein and Infliximab Exposure Calls for Objective Crohn's Disease Activity Measures for the Evaluation of Treatment Effects at Treatment Failure
JF - Therapeutic drug monitoring : official journal of the International Association of Therapeutic Drug Monitoring and Clinical Toxicology
N2 - Background: Circulating infliximab (IFX) concentrations correlate with clinical outcomes, forming the basis of the IFX concentration monitoring in patients with Crohn's disease. This study aims to investigate and refine the exposure-response relationship by linking the disease activity markers "Crohn's disease activity index" (CDAI) and C-reactive protein (CRP) to IFX exposure. In addition, we aim to explore the correlations between different disease markers and exposure metrics.
Methods: Data from 47 Crohn's disease patients of a randomized controlled trial were analyzed post hoc. All patients had secondary treatment failure at inclusion and had received intensified IFX of 5 mg/kg every 4 weeks for up to 20 weeks. Graphical analyses were performed to explore exposure-response relationships. Metrics of exposure included area under the concentration-time curve (AUC) and trough concentrations (Cmin). Disease activity was measured by CDAI and CRP values, their change from baseline/last visit, and response/remission outcomes at week 12.
Results: Although trends toward lower Cmin and lower AUC in nonresponders were observed, neither CDAI nor CRP showed consistent trends of lower disease activity with higher IFX exposure across the 30 evaluated relationships. As can be expected, Cmin and AUC were strongly correlated with each other. Contrarily, the disease activity markers were only weakly correlated with each other.
Conclusions: No significant relationship between disease activity, as evaluated by CDAI or CRP, and IFX exposure was identified. AUC did not add benefit compared with Cmin. These findings support the continued use of Cmin and call for stringent objective disease activity (bio-)markers (eg, endoscopy) to form the basis of personalized IFX therapy for Crohn's disease patients with IFX treatment failure.
Y1 - 2019
U6 - https://doi.org/10.1097/FTD.0000000000000590
SN - 0163-4356
SN - 1536-3694
VL - 41
IS - 2
SP - 235
EP - 242
PB - Lippincott Williams & Wilkins
CY - Philadelphia
ER -
TY - JOUR
A1 - Knöchel, Jane
A1 - Kloft, Charlotte
A1 - Huisinga, Wilhelm
T1 - Understanding and reducing complex systems pharmacology models based on a novel input-response index
JF - Journal of pharmacokinetics and pharmacodynamics
N2 - A growing understanding of complex processes in biology has led to large-scale mechanistic models of pharmacologically relevant processes. These models are increasingly used to study the response of the system to a given input or stimulus, e.g., after drug administration. Understanding the input–response relationship, however, is often a challenging task due to the complexity of the interactions between its constituents as well as the size of the models. An approach that quantifies the importance of the different constituents for a given input–output relationship and allows to reduce the dynamics to its essential features is therefore highly desirable. In this article, we present a novel state- and time-dependent quantity called the input–response index that quantifies the importance of state variables for a given input–response relationship at a particular time. It is based on the concept of time-bounded controllability and observability, and defined with respect to a reference dynamics. In application to the brown snake venom–fibrinogen (Fg) network, the input–response indices give insight into the coordinated action of specific coagulation factors and about those factors that contribute only little to the response. We demonstrate how the indices can be used to reduce large-scale models in a two-step procedure: (i) elimination of states whose dynamics have only minor impact on the input–response relationship, and (ii) proper lumping of the remaining (lower order) model. In application to the brown snake venom–fibrinogen network, this resulted in a reduction from 62 to 8 state variables in the first step, and a further reduction to 5 state variables in the second step. We further illustrate that the sequence, in which a recursive algorithm eliminates and/or lumps state variables, has an impact on the final reduced model. The input–response indices are particularly suited to determine an informed sequence, since they are based on the dynamics of the original system. In summary, the novel measure of importance provides a powerful tool for analysing the complex dynamics of large-scale systems and a means for very efficient model order reduction of nonlinear systems.
KW - Control theory
KW - Model order reduction
KW - Blood coagulation network
KW - Nonlinear systems
Y1 - 2017
U6 - https://doi.org/10.1007/s10928-017-9561-x
SN - 1567-567X
SN - 1573-8744
VL - 45
IS - 1
SP - 139
EP - 157
PB - Springer Science + Business Media B.V.
CY - New York
ER -
TY - JOUR
A1 - Stachanow, Viktoria
A1 - Neumann, Uta
A1 - Blankenstein, Oliver
A1 - Bindellini, Davide
A1 - Melin, Johanna
A1 - Ross, Richard
A1 - Whitaker, Martin J. J.
A1 - Huisinga, Wilhelm
A1 - Michelet, Robin
A1 - Kloft, Charlotte
T1 - Exploring dried blood spot cortisol concentrations as an alternative for monitoring pediatric adrenal insufficiency patients
BT - a model-based analysis
JF - Frontiers in pharmacology
N2 - Congenital adrenal hyperplasia (CAH) is the most common form of adrenal insufficiency in childhood; it requires cortisol replacement therapy with hydrocortisone (HC, synthetic cortisol) from birth and therapy monitoring for successful treatment. In children, the less invasive dried blood spot (DBS) sampling with whole blood including red blood cells (RBCs) provides an advantageous alternative to plasma sampling.
Potential differences in binding/association processes between plasma and DBS however need to be considered to correctly interpret DBS measurements for therapy monitoring. While capillary DBS samples would be used in clinical practice, venous cortisol DBS samples from children with adrenal insufficiency were analyzed due to data availability and to directly compare and thus understand potential differences between venous DBS and plasma. A previously published HC plasma pharmacokinetic (PK) model was extended by leveraging these DBS concentrations.
In addition to previously characterized binding of cortisol to albumin (linear process) and corticosteroid-binding globulin (CBG; saturable process), DBS data enabled the characterization of a linear cortisol association with RBCs, and thereby providing a quantitative link between DBS and plasma cortisol concentrations. The ratio between the observed cortisol plasma and DBS concentrations varies highly from 2 to 8. Deterministic simulations of the different cortisol binding/association fractions demonstrated that with higher blood cortisol concentrations, saturation of cortisol binding to CBG was observed, leading to an increase in all other cortisol binding fractions.
In conclusion, a mathematical PK model was developed which links DBS measurements to plasma exposure and thus allows for quantitative interpretation of measurements of DBS samples.
KW - adrenal insufficiency
KW - cortisol
KW - dried blood spots
KW - pediatrics
KW - pharmacokinetics
KW - binding
KW - association
KW - red blood cells
Y1 - 2022
U6 - https://doi.org/10.3389/fphar.2022.819590
SN - 1663-9812
VL - 13
PB - Frontiers Media
CY - Lausanne
ER -
TY - JOUR
A1 - Hijazi, Saddam
A1 - Freitag, Melina A.
A1 - Landwehr, Niels
T1 - POD-Galerkin reduced order models and physics-informed neural networks for solving inverse problems for the Navier-Stokes equations
JF - Advanced modeling and simulation in engineering sciences : AMSES
N2 - We present a Reduced Order Model (ROM) which exploits recent developments in Physics Informed Neural Networks (PINNs) for solving inverse problems for the Navier-Stokes equations (NSE). In the proposed approach, the presence of simulated data for the fluid dynamics fields is assumed. A POD-Galerkin ROM is then constructed by applying POD on the snapshots matrices of the fluid fields and performing a Galerkin projection of the NSE (or the modified equations in case of turbulence modeling) onto the POD reduced basis. A POD-Galerkin PINN ROM is then derived by introducing deep neural networks which approximate the reduced outputs with the input being time and/or parameters of the model. The neural networks incorporate the physical equations (the POD-Galerkin reduced equations) into their structure as part of the loss function. Using this approach, the reduced model is able to approximate unknown parameters such as physical constants or the boundary conditions. A demonstration of the applicability of the proposed ROM is illustrated by three cases which are the steady flow around a backward step, the flow around a circular cylinder and the unsteady turbulent flow around a surface mounted cubic obstacle.
KW - Proper orthogonal decomposition
KW - Inverse problems
KW - Physics-based machine learning
KW - Navier-Stokes equations
Y1 - 2023
U6 - https://doi.org/10.1186/s40323-023-00242-2
SN - 2213-7467
VL - 10
IS - 1
PB - SpringerOpen
CY - Berlin
ER -
TY - JOUR
A1 - Molkenthin, Christian
A1 - Donner, Christian
A1 - Reich, Sebastian
A1 - Zöller, Gert
A1 - Hainzl, Sebastian
A1 - Holschneider, Matthias
A1 - Opper, Manfred
T1 - GP-ETAS: semiparametric Bayesian inference for the spatio-temporal epidemic type aftershock sequence model
JF - Statistics and Computing
N2 - The spatio-temporal epidemic type aftershock sequence (ETAS) model is widely used to describe the self-exciting nature of earthquake occurrences. While traditional inference methods provide only point estimates of the model parameters, we aim at a fully Bayesian treatment of model inference, allowing naturally to incorporate prior knowledge and uncertainty quantification of the resulting estimates. Therefore, we introduce a highly flexible, non-parametric representation for the spatially varying ETAS background intensity through a Gaussian process (GP) prior. Combined with classical triggering functions this results in a new model formulation, namely the GP-ETAS model. We enable tractable and efficient Gibbs sampling by deriving an augmented form of the GP-ETAS inference problem. This novel sampling approach allows us to assess the posterior model variables conditioned on observed earthquake catalogues, i.e., the spatial background intensity and the parameters of the triggering function. Empirical results on two synthetic data sets indicate that GP-ETAS outperforms standard models and thus demonstrate the predictive power for observed earthquake catalogues including uncertainty quantification for the estimated parameters. Finally, a case study for the l'Aquila region, Italy, with the devastating event on 6 April 2009, is presented.
KW - Self-exciting point process
KW - Hawkes process
KW - Spatio-temporal ETAS model
KW - Bayesian inference
KW - Sampling
KW - Earthquake modeling
KW - Gaussian process
KW - Data augmentation
Y1 - 2022
U6 - https://doi.org/10.1007/s11222-022-10085-3
SN - 0960-3174
SN - 1573-1375
VL - 32
IS - 2
PB - Springer
CY - Dordrecht
ER -
TY - JOUR
A1 - Kucharski, Maciej
A1 - Ergintav, Arzu
A1 - Ahmad, Wael Abdullah
A1 - Krstić, Miloš
A1 - Ng, Herman Jalli
A1 - Kissinger, Dietmar
T1 - A Scalable 79-GHz Radar Platform Based on Single-Channel Transceivers
JF - IEEE Transactions on Microwave Theory and Techniques
N2 - This paper presents a scalable E-band radar platform based on single-channel fully integrated transceivers (TRX) manufactured using 130-nm silicon-germanium (SiGe) BiCMOS technology. The TRX is suitable for flexible radar systems exploiting massive multiple-input-multipleoutput (MIMO) techniques for multidimensional sensing. A fully integrated fractional-N phase-locked loop (PLL) comprising a 39.5-GHz voltage-controlled oscillator is used to generate wideband frequency-modulated continuous-wave (FMCW) chirp for E-band radar front ends. The TRX is equipped with a vector modulator (VM) for high-speed carrier modulation and beam-forming techniques. A single TRX achieves 19.2-dBm maximum output power and 27.5-dB total conversion gain with input-referred 1-dB compression point of -10 dBm. It consumes 220 mA from 3.3-V supply and occupies 3.96 mm(2) silicon area. A two-channel radar platform based on full-custom TRXs and PLL was fabricated to demonstrate high-precision and high-resolution FMCW sensing. The radar enables up to 10-GHz frequency ramp generation in 74-84-GHz range, which results in 1.5-cm spatial resolution. Due to high output power, thus high signal-to-noise ratio (SNR), a ranging precision of 7.5 mu m for a target at 2 m was achieved. The proposed architecture supports scalable multichannel applications for automotive FMCW using a single local oscillator (LO).
KW - Automotive
KW - E-band
KW - frequency-modulated continuous-wave (FMCW)
KW - patch antenna
KW - phase-locked loop (PLL)
KW - power amplifier (PA)
KW - radar
KW - scalable
KW - transceiver (TRX)
Y1 - 2019
U6 - https://doi.org/10.1109/TMTT.2019.2914104
SN - 0018-9480
SN - 1557-9670
VL - 67
IS - 9
SP - 3882
EP - 3896
PB - Inst. of Electr. and Electronics Engineers
CY - Piscataway
ER -
TY - JOUR
A1 - Sharma, Shubham
A1 - Hainzl, Sebastian
A1 - Zöller, Gert
A1 - Holschneider, Matthias
T1 - Is Coulomb stress the best choice for aftershock forecasting?
JF - Journal of geophysical research : Solid earth
N2 - The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model.
Y1 - 2020
U6 - https://doi.org/10.1029/2020JB019553
SN - 2169-9313
SN - 2169-9356
VL - 125
IS - 9
PB - American Geophysical Union
CY - Washington
ER -
TY - JOUR
A1 - Engbert, Ralf
A1 - Rabe, Maximilian Michael
A1 - Kliegl, Reinhold
A1 - Reich, Sebastian
T1 - Sequential data assimilation of the stochastic SEIR epidemic model for regional COVID-19 dynamics
JF - Bulletin of mathematical biology : official journal of the Society for Mathematical Biology
N2 - Newly emerging pandemics like COVID-19 call for predictive models to implement precisely tuned responses to limit their deep impact on society. Standard epidemic models provide a theoretically well-founded dynamical description of disease incidence. For COVID-19 with infectiousness peaking before and at symptom onset, the SEIR model explains the hidden build-up of exposed individuals which creates challenges for containment strategies. However, spatial heterogeneity raises questions about the adequacy of modeling epidemic outbreaks on the level of a whole country. Here, we show that by applying sequential data assimilation to the stochastic SEIR epidemic model, we can capture the dynamic behavior of outbreaks on a regional level. Regional modeling, with relatively low numbers of infected and demographic noise, accounts for both spatial heterogeneity and stochasticity. Based on adapted models, short-term predictions can be achieved. Thus, with the help of these sequential data assimilation methods, more realistic epidemic models are within reach.
KW - Stochastic epidemic model
KW - Sequential data assimilation
KW - Ensemble Kalman
KW - filter
KW - COVID-19
Y1 - 2020
U6 - https://doi.org/10.1007/s11538-020-00834-8
SN - 0092-8240
SN - 1522-9602
VL - 83
IS - 1
PB - Springer
CY - New York
ER -
TY - JOUR
A1 - Gerlach, Moritz Reinhardt
A1 - Glück, Jochen
T1 - On a convergence theorem for semigroups of positive integral operators
JF - Comptes Rendus Mathematique
N2 - We give a new and very short proof of a theorem of Greiner asserting that a positive and contractive -semigroup on an -space is strongly convergent in case it has a strictly positive fixed point and contains an integral operator. Our proof is a streamlined version of a much more general approach to the asymptotic theory of positive semigroups developed recently by the authors. Under the assumptions of Greiner's theorem, this approach becomes particularly elegant and simple. We also give an outlook on several generalisations of this result.
Y1 - 2017
U6 - https://doi.org/10.1016/j.crma.2017.07.017
SN - 1631-073X
SN - 1778-3569
VL - 355
SP - 973
EP - 976
PB - Elsevier
CY - Paris
ER -
TY - JOUR
A1 - Gerlach, Moritz Reinhardt
T1 - Convergence of dynamics and the Perron-Frobenius operator
JF - Israel Journal of Mathematics
N2 - We complete the picture how the asymptotic behavior of a dynamical system is reflected by properties of the associated Perron-Frobenius operator. Our main result states that strong convergence of the powers of the Perron-Frobenius operator is equivalent to setwise convergence of the underlying dynamic in the measure algebra. This situation is furthermore characterized by uniform mixing-like properties of the system.
Y1 - 2018
U6 - https://doi.org/10.1007/s11856-018-1671-7
SN - 0021-2172
SN - 1565-8511
VL - 225
IS - 1
SP - 451
EP - 463
PB - Hebrew univ magnes press
CY - Jerusalem
ER -
TY - JOUR
A1 - Gerlach, Moritz Reinhardt
A1 - Glück, Jochen
T1 - Convergence of positive operator semigroups
JF - Transactions of the American Mathematical Society
N2 - We present new conditions for semigroups of positive operators to converge strongly as time tends to infinity. Our proofs are based on a novel approach combining the well-known splitting theorem by Jacobs, de Leeuw, and Glicksberg with a purely algebraic result about positive group representations. Thus, we obtain convergence theorems not only for one-parameter semigroups but also for a much larger class of semigroup representations. Our results allow for a unified treatment of various theorems from the literature that, under technical assumptions, a bounded positive C-0-semigroup containing or dominating a kernel operator converges strongly as t ->infinity. We gain new insights into the structure theoretical background of those theorems and generalize them in several respects; especially we drop any kind of continuity or regularity assumption with respect to the time parameter.
KW - Positive semigroups
KW - semigroup representations
KW - asymptotic behavior
KW - kernel operator
Y1 - 2019
U6 - https://doi.org/10.1090/tran/7836
SN - 0002-9947
SN - 1088-6850
VL - 372
IS - 9
SP - 6603
EP - 6627
PB - American Mathematical Soc.
CY - Providence
ER -
TY - JOUR
A1 - Edeko, Nikolai
A1 - Gerlach, Moritz Reinhardt
A1 - Kühner, Viktoria
T1 - Measure-preserving semiflows and one-parameter Koopman semigroups
JF - Semigroup forum
N2 - For a finite measure space X, we characterize strongly continuous Markov lattice semigroups on Lp(X) by showing that their generator A acts as a derivation on the dense subspace D(A)L(X). We then use this to characterize Koopman semigroups on Lp(X) if X is a standard probability space. In addition, we show that every measurable and measure-preserving flow on a standard probability space is isomorphic to a continuous flow on a compact Borel probability space.
KW - Measure-preserving semiflow
KW - Koopman semigroup
KW - Derivation
KW - Topological model
Y1 - 2019
U6 - https://doi.org/10.1007/s00233-018-9960-3
SN - 0037-1912
SN - 1432-2137
VL - 98
IS - 1
SP - 48
EP - 63
PB - Springer
CY - New York
ER -
TY - JOUR
A1 - Gerlach, Moritz Reinhardt
A1 - Glück, Jochen
T1 - Lower bounds and the asymptotic behaviour of positive operator semigroups
JF - Ergodic theory and dynamical systems
N2 - If (T-t) is a semigroup of Markov operators on an L-1-space that admits a nontrivial lower bound, then a well-known theorem of Lasota and Yorke asserts that the semigroup is strongly convergent as t -> infinity. In this article we generalize and improve this result in several respects. First, we give a new and very simple proof for the fact that the same conclusion also holds if the semigroup is merely assumed to be bounded instead of Markov. As a main result, we then prove a version of this theorem for semigroups which only admit certain individual lower bounds. Moreover, we generalize a theorem of Ding on semigroups of Frobenius-Perron operators. We also demonstrate how our results can be adapted to the setting of general Banach lattices and we give some counterexamples to show optimality of our results. Our methods combine some rather concrete estimates and approximation arguments with abstract functional analytical tools. One of these tools is a theorem which relates the convergence of a time-continuous operator semigroup to the convergence of embedded discrete semigroups.
Y1 - 2017
U6 - https://doi.org/10.1017/etds.2017.9
SN - 0143-3857
SN - 1469-4417
VL - 38
SP - 3012
EP - 3041
PB - Cambridge Univ. Press
CY - New York
ER -
TY - JOUR
A1 - Gerlach, Moritz Reinhardt
A1 - Glück, Jochen
T1 - Mean ergodicity vs weak almost periodicity
JF - Studia mathematica
N2 - We provide explicit examples of positive and power-bounded operators on c(0) and l(infinity) which are mean ergodic but not weakly almost periodic. As a consequence we prove that a countably order complete Banach lattice on which every positive and power-bounded mean ergodic operator is weakly almost periodic is necessarily a KB-space. This answers several open questions from the literature. Finally, we prove that if T is a positive mean ergodic operator with zero fixed space on an arbitrary Banach lattice, then so is every power of T .
KW - positive operators
KW - weakly almost periodic
KW - order continuous norm
KW - KB-space
KW - mean ergodic
Y1 - 2019
U6 - https://doi.org/10.4064/sm170918-20-3
SN - 0039-3223
SN - 1730-6337
VL - 248
IS - 1
SP - 45
EP - 56
PB - Polska Akademia Nauk, Instytut Matematyczny
CY - Warszawa
ER -
TY - JOUR
A1 - Gerlach, Moritz
A1 - Glück, Jochen
A1 - Kunze, Markus
T1 - Stability of transition semigroups and applications to parabolic equations
JF - Transactions of the American Mathematical Society
N2 - This paper deals with the long-term behavior of positive operator semigroups on spaces of bounded functions and of signed measures, which have applications to parabolic equations with unbounded coefficients and to stochas-tic analysis. The main results are a Tauberian type theorem characterizing the convergence to equilibrium of strongly Feller semigroups and a generalization of a classical convergence theorem of Doob. None of these results requires any kind of time regularity of the semigroup.
KW - Transition probabilities
KW - strong Feller property
KW - asymptotic
KW - behavior
KW - invariant measure
KW - parabolic equations
Y1 - 2023
U6 - https://doi.org/10.1090/tran/8620
SN - 0002-9947
SN - 1088-6850
VL - 376
IS - 1
SP - 153
EP - 180
PB - American Mathematical Soc.
CY - Providence
ER -
TY - JOUR
A1 - Dimitrova, Ilinka
A1 - Koppitz, Jörg
T1 - On relative ranks of the semigroup of orientation-preserving transformations on infinite chain with restricted range
JF - Communications in algebra
N2 - 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.
KW - Order-preserving transformations
KW - orientation-preserving
KW - transformations
KW - relative rank
KW - restricted range
KW - transformation
KW - semigroups on infinite chain
Y1 - 2022
U6 - https://doi.org/10.1080/00927872.2021.2000998
SN - 0092-7872
SN - 1532-4125
VL - 50
IS - 5
SP - 2157
EP - 2168
PB - Taylor & Francis Group
CY - Philadelphia
ER -
TY - JOUR
A1 - Dimitrova, Ilinka
A1 - Koppitz, Jörg
T1 - On relative ranks of the semigroup of orientation-preserving transformations on infinite chains
JF - Asian-European journal of mathematics
N2 - 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.
KW - Transformation semigroups on infinite chains
KW - order-preserving
KW - transformations
KW - orientation-preserving transformations
KW - relative rank
Y1 - 2020
U6 - https://doi.org/10.1142/S1793557121501461
SN - 1793-5571
SN - 1793-7183
VL - 14
IS - 08
PB - World Scientific
CY - Singapore
ER -
TY - JOUR
A1 - Kretzschmar, Mirjam E.
A1 - Ashby, Ben
A1 - Fearon, Elizabeth
A1 - Overton, Christopher E.
A1 - Panovska-Griffiths, Jasmina
A1 - Pellis, Lorenzo
A1 - Quaife, Matthew
A1 - Rozhnova, Ganna
A1 - Scarabel, Francesca
A1 - Stage, Helena B.
A1 - Swallow, Ben
A1 - Thompson, Robin N.
A1 - Tildesley, Michael J.
A1 - Villela, Daniel Campos
T1 - Challenges for modelling interventions for future pandemics
JF - Epidemics
N2 - Mathematical modelling and statistical inference provide a framework to evaluate different non-pharmaceutical and pharmaceutical interventions for the control of epidemics that has been widely used during the COVID-19 pandemic. In this paper, lessons learned from this and previous epidemics are used to highlight the challenges for future pandemic control. We consider the availability and use of data, as well as the need for correct parameterisation and calibration for different model frameworks. We discuss challenges that arise in describing and distinguishing between different interventions, within different modelling structures, and allowing both within and between host dynamics. We also highlight challenges in modelling the health economic and political aspects of interventions. Given the diversity of these challenges, a broad variety of interdisciplinary expertise is needed to address them, combining mathematical knowledge with biological and social insights, and including health economics and communication skills. Addressing these challenges for the future requires strong cross disciplinary collaboration together with close communication between scientists and policy makers.
KW - Mathematical models
KW - Pandemics
KW - Pharmaceutical interventions
KW - Non-pharmaceutical interventions
KW - Policy support
Y1 - 2022
U6 - https://doi.org/10.1016/j.epidem.2022.100546
SN - 1755-4365
SN - 1878-0067
VL - 38
PB - Elsevier
CY - Amsterdam
ER -