Stochastic Galerkin methods for non-affine coefficient representations are known to cause major difficulties from theoretical and numerical points of view. In this work, an adaptive Galerkin FE method for linear parametric PDEs with lognormal coefficients discretized in Hermite chaos polynomials is derived. It employs problem-adapted function spaces to ensure solvability of the variational formulation. The inherently high computational complexity of the parametric operator is made tractable by using hierarchical tensor representations. For this, a new tensor train format of the lognormal coefficient is derived and verified numerically. The central novelty is the derivation of a reliable residual-based a posteriori error estimator. This can be regarded as a unique feature of stochastic Galerkin methods. It allows for an adaptive algorithm to steer the refinements of the physical mesh and the anisotropic Wiener chaos polynomial degrees. For the evaluation of the error estimator to become feasible, a numerically efficient tensor format discretization is developed. Benchmark examples with unbounded lognormal coefficient fields illustrate the performance of the proposed Galerkin discretization and the fully adaptive algorithm.

The Guide to the Expression of Uncertainty in Measurement (GUM) has led to a harmonization of uncertainty evaluation throughout metrology. The simplicity of its employed methodology has fostered the broad acceptance of the GUM among metrologists. However, this simplicity also compromises best practice and does not provide state-of-the-art data analysis. Specifically, metrologists often possess useful prior knowledge about the measurand which cannot be accounted for by the GUM. Bayesian uncertainty evaluation, on the other hand, takes prior knowledge about the measurand as its starting point. While this technique has been successfully applied in several instances in metrology, its broad access for metrologists is still lacking. One reason for this is that application of Bayesian methods and their computation, e.g. by means of Markov chain Monte Carlo (MCMC) methods, usually requires statistical skills.In this work, we propose a simple method for Bayesian uncertainty evaluation applicable to measurement models that depend linearly on a single input quantity for which Type A information is available, and several input quantities for which Type B information is given. This category of measurement models covers many uncertainty evaluations in metrology. The approach utilizes a specific class of prior distributions to encode the prior information about the measurand. Explicit guidance is provided on how to choose the parameters of the prior in the light of one's prior knowledge. MCMC methods are not required for the calculation of results; instead, a simple Monte Carlo procedure is developed that allows to draw uncorrelated samples from the posterior distribution, which greatly simplifies convergence assessments. The proposed Bayesian approach allows the treatment of small numbers of observations for a Type A uncertainty evaluation, including the case of only a single observation. Examples are provided that illustrate the proposed approach and corresponding open source Python software is made available.

Our systematic study shows that vimentin forms individual 10 nm diameter filaments in the presence of monovalent ions and thicker filaments that aggregate into dense networks in the presence of multivalent ions.

A sampling-free approach to Bayesian inversion with an explicit polynomial representation of the parameter densities is developed, based on an affineparametric representation of a linear forward model. This becomes feasible due to the complete treatment in function spaces, which requires an efficient model reduction technique for numerical computations.The advocated perspective yields the crucial benefit that error bounds can be derived for all occuring approximations, leading to provable convergence subject to the discretization parameters. Moreover, it enables a fully adaptive a posteriori control with automatic problem-dependent adjustments of the employed discretizations. The method is discussed in the context of modern hierarchical tensor representations, which are used for the evaluation of a random PDE (the forward model) and the subsequent high-dimensional quadrature of the log-likelihood, alleviating the 'curse of dimensionality'. Numerical experiments demonstrate the performance and confirm the theoretical results.

Fourier transform infrared (FTIR) spectroscopy is a powerful technique in analytical chemistry. Typically, spatially distributed spectra of the substance of interest are conducted simultaneously using FTIR spectrometers equipped with array detectors. Scanning-based methods such as near-field FTIR spectroscopy, on the other hand, are a promising alternative providing higher spatial resolution. However, serial recording severely limits their application due to the long acquisition times involved and the resulting stability issues. We demonstrate that it is possible to significantly reduce the measurement time of scanning methods by applying the mathematical technique of low-rank matrix reconstruction. Data from a previous pilot study of Leishmania strains are analyzed by randomly selecting 5% of the interferometer samples. The results obtained for bioanalytical fingerprinting using the proposed approach are shown to be essentially the same as those obtained from the full set of data. This finding can significantly foster the practical applicability of high-resolution serial scanning techniques in analytical chemistry and is also expected to improve other applications of FTIR spectroscopy and spectromicroscopy.