Refine
Year of publication
- 2019 (59) (remove)
Document Type
- Article (44)
- Doctoral Thesis (8)
- Other (3)
- Postprint (2)
- Monograph/Edited Volume (1)
- Review (1)
Language
- English (59) (remove)
Is part of the Bibliography
- yes (59) (remove)
Keywords
- Dirac operator (2)
- Magnus expansion (2)
- asymptotic method (2)
- clone (2)
- discrepancy principle (2)
- higher-order Sturm–Liouville problems (2)
- inverse Sturm–Liouville problems (2)
- linear hypersubstitution (2)
- linear term (2)
- nonlinear operator (2)
Institute
- Institut für Mathematik (59) (remove)
Continuous insight into biological processes has led to the development of large-scale, mechanistic systems biology models of pharmacologically relevant networks. While these models are typically designed to study the impact of diverse stimuli or perturbations on multiple system variables, the focus in pharmacological research is often on a specific input, e.g., the dose of a drug, and a specific output related to the drug effect or response in terms of some surrogate marker.
To study a chosen input-output pair, the complexity of the interactions as well as the size of the models hinders easy access and understanding of the details of the input-output relationship.
The objective of this thesis is the development of a mathematical approach, in specific a model reduction technique, that allows (i) to quantify the importance of the different state variables for a given input-output relationship, and (ii) to reduce the dynamics to its essential features -- allowing for a physiological interpretation of state variables as well as parameter estimation in the statistical analysis of clinical data. We develop a model reduction technique using a control theoretic setting by first defining a novel type of time-limited controllability and observability gramians for nonlinear systems. We then show the superiority of the time-limited generalised gramians for nonlinear systems in the context of balanced truncation for a benchmark system from control theory.
The concept of time-limited controllability and observability gramians is subsequently used to introduce a state and time-dependent quantity called the input-response (ir) index that quantifies the importance of state variables for a given input-response relationship at a particular time.
We subsequently link our approach to sensitivity analysis, thus, enabling for the first time the use of sensitivity coefficients for state space reduction. The sensitivity based ir-indices are given as a product of two sensitivity coefficients. This allows not only for a computational more efficient calculation but also for a clear distinction of the extent to which the input impacts a state variable and the extent to which a state variable impacts the output.
The ir-indices give insight into the coordinated action of specific state variables for a chosen input-response relationship.
Our developed model reduction technique results in reduced models that still allow for a mechanistic interpretation in terms of the quantities/state variables of the original system, which is a key requirement in the field of systems pharmacology and systems biology and distinguished the reduced models from so-called empirical drug effect models. The ir-indices are explicitly defined with respect to a reference trajectory and thereby dependent on the initial state (this is an important feature of the measure). This is demonstrated for an example from the field of systems pharmacology, showing that the reduced models are very informative in their ability to detect (genetic) deficiencies in certain physiological entities. Comparing our novel model reduction technique to the already existing techniques shows its superiority.
The novel input-response index as a measure of the importance of state variables provides a powerful tool for understanding the complex dynamics of large-scale systems in the context of a specific drug-response relationship. Furthermore, the indices provide a means for a very efficient model order reduction and, thus, an important step towards translating insight from biological processes incorporated in detailed systems pharmacology models into the population analysis of clinical data.
Low thermal conductivity boulder with high porosity identified on C-type asteroid (162173) Ryugu
(2019)
C-type asteroids are among the most pristine objects in the Solar System, but little is known about their interior structure and surface properties. Telescopic thermal infrared observations have so far been interpreted in terms of a regolith-covered surface with low thermal conductivity and particle sizes in the centimetre range. This includes observations of C-type asteroid (162173) Ryugu1,2,3. However, on arrival of the Hayabusa2 spacecraft at Ryugu, a regolith cover of sand- to pebble-sized particles was found to be absent4,5 (R.J. et al., manuscript in preparation). Rather, the surface is largely covered by cobbles and boulders, seemingly incompatible with the remote-sensing infrared observations. Here we report on in situ thermal infrared observations of a boulder on the C-type asteroid Ryugu. We found that the boulder’s thermal inertia was much lower than anticipated based on laboratory measurements of meteorites, and that a surface covered by such low-conductivity boulders would be consistent with remote-sensing observations. Our results furthermore indicate high boulder porosities as well as a low tensile strength in the few hundred kilopascal range. The predicted low tensile strength confirms the suspected observational bias6 in our meteorite collections, as such asteroidal material would be too frail to survive atmospheric entry7
Low thermal conductivity boulder with high porosity identified on C-type asteroid (162173) Ryugu
(2019)
C-type asteroids are among the most pristine objects in the Solar System, but little is known about their interior structure and surface properties. Telescopic thermal infrared observations have so far been interpreted in terms of a regolith-covered surface with low thermal conductivity and particle sizes in the centimetre range. This includes observations of C-type asteroid (162173) Ryugu1,2,3. However, on arrival of the Hayabusa2 spacecraft at Ryugu, a regolith cover of sand- to pebble-sized particles was found to be absent4,5 (R.J. et al., manuscript in preparation). Rather, the surface is largely covered by cobbles and boulders, seemingly incompatible with the remote-sensing infrared observations. Here we report on in situ thermal infrared observations of a boulder on the C-type asteroid Ryugu. We found that the boulder’s thermal inertia was much lower than anticipated based on laboratory measurements of meteorites, and that a surface covered by such low-conductivity boulders would be consistent with remote-sensing observations. Our results furthermore indicate high boulder porosities as well as a low tensile strength in the few hundred kilopascal range. The predicted low tensile strength confirms the suspected observational bias6 in our meteorite collections, as such asteroidal material would be too frail to survive atmospheric entry7.
An efficient immunosurveillance of CD8(+) T cells in the periphery depends on positive/negative selection of thymocytes and thus on the dynamics of antigen degradation and epitope production by thymoproteasome and immunoproteasome in the thymus. Although studies in mouse systems have shown how thymoproteasome activity differs from that of immunoproteasome and strongly impacts the T cell repertoire, the proteolytic dynamics and the regulation of human thymoproteasome are unknown. By combining biochemical and computational modeling approaches, we show here that human 20S thymoproteasome and immunoproteasome differ not only in the proteolytic activity of the catalytic sites but also in the peptide transport. These differences impinge upon the quantity of peptide products rather than where the substrates are cleaved. The comparison of the two human 20S proteasome isoforms depicts different processing of antigens that are associated to tumors and autoimmune diseases.
We develop a technique for the multivariate data analysis of perturbed self-sustained oscillators. The approach is based on the reconstruction of the phase dynamics model from observations and on a subsequent exploration of this model. For the system, driven by several inputs, we suggest a dynamical disentanglement procedure, allowing us to reconstruct the variability of the system's output that is due to a particular observed input, or, alternatively, to reconstruct the variability which is caused by all the inputs except for the observed one. We focus on the application of the method to the vagal component of the heart rate variability caused by a respiratory influence. We develop an algorithm that extracts purely respiratory-related variability, using a respiratory trace and times of R-peaks in the electrocardiogram. The algorithm can be applied to other systems where the observed bivariate data can be represented as a point process and a slow continuous signal, e.g. for the analysis of neuronal spiking. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.
For a singularly perturbed parabolic - ODE system we construct the asymptotic expansion in the small parameter in the case, when the degenerate equation has a double root. Such systems, which are called partly dissipative reaction-diffusion systems, are used to model various natural processes, including the signal transmission along axons, solid combustion and the kinetics of some chemical reactions. It turns out that the algorithm of the construction of the boundary layer functions and the behavior of the solution in the boundary layers essentially differ from that ones in case of a simple root. The multizonal initial and boundary layers behaviour was stated.
Probabilistic integration of a continuous dynamical system is a way of systematically introducing discretisation error, at scales no larger than errors introduced by standard numerical discretisation, in order to enable thorough exploration of possible responses of the system to inputs. It is thus a potentially useful approach in a number of applications such as forward uncertainty quantification, inverse problems, and data assimilation. We extend the convergence analysis of probabilistic integrators for deterministic ordinary differential equations, as proposed by Conrad et al. (Stat Comput 27(4):1065-1082, 2017. ), to establish mean-square convergence in the uniform norm on discrete- or continuous-time solutions under relaxed regularity assumptions on the driving vector fields and their induced flows. Specifically, we show that randomised high-order integrators for globally Lipschitz flows and randomised Euler integrators for dissipative vector fields with polynomially bounded local Lipschitz constants all have the same mean-square convergence rate as their deterministic counterparts, provided that the variance of the integration noise is not of higher order than the corresponding deterministic integrator. These and similar results are proven for probabilistic integrators where the random perturbations may be state-dependent, non-Gaussian, or non-centred random variables.
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.
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.
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.
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.
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.
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.
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.
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.
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).