An integral conservative gridding--algorithm using Hermitian curve interpolation

The problem of re-sampling spatially distributed data organized into regular or irregular grids to finer or coarser resolution is a common task in data processing. This procedure is known as 'gridding' or 're-binning'. Depending on the quantity the data represents, the gridding-algorithm has to meet different requirements. For example, histogrammed physical quantities such as mass or energy have to be re-binned in order to conserve the overall integral. Moreover, if the quantity is positive definite, negative sampling values should be avoided. The gridding process requires a re-distribution of the original data set to a user-requested grid according to a distribution function. The distribution function can be determined on the basis of the given data by interpolation methods. In general, accurate interpolation with respect to multiple boundary conditions of heavily fluctuating data requires polynomial interpolation functions of second or even higher order. However, this may result in unrealistic deviations (overshoots or undershoots) of the interpolation function from the data. Accordingly, the re-sampled data may overestimate or underestimate the given data by a significant amount. The gridding-algorithm presented in this work was developed in order to overcome these problems. Instead of a straightforward interpolation of the given data using high-order polynomials, a parametrized Hermitian interpolation curve was used to approximate the integrated data set. A single parameter is determined by which the user can control the behavior of the interpolation function, i.e. the amount of overshoot and undershoot. Furthermore, it is shown how the algorithm can be extended to multidimensional grids. The algorithm was compared to commonly used gridding-algorithms using linear and cubic interpolation functions. It is shown that such interpolation functions may overestimate or underestimate the source data by about 10-20%, while the new algorithm can be tuned to significantly reduce these interpolation errors. The accuracy of the new algorithm was tested on a series of x-ray CT-images (head and neck, lung, pelvis). The new algorithm significantly improves the accuracy of the sampled images in terms of the mean square error and a quality index introduced by Wang and Bovik (2002 IEEE Signal Process. Lett. 9 81-4).

Monitoring protein-protein interactions between the mammalian integral membrane transporters and PDZ-interacting partners using a modified split-ubiquitin membrane yeast two-hybrid system

PDZ-binding motifs are found in the C-terminal tails of numerous integral membrane proteins where they mediate specific protein-protein interactions by binding to PDZ-containing proteins. Conventional yeast two-hybrid screens have been used to probe protein-protein interactions of these soluble C termini. However, to date no in vivo technology has been available to study interactions between the full-length integral membrane proteins and their cognate PDZ-interacting partners. We previously developed a split-ubiquitin membrane yeast two-hybrid (MYTH) system to test interactions between such integral membrane proteins by using a transcriptional output based on cleavage of a transcription factor from the C terminus of membrane-inserted baits. Here we modified MYTH to permit detection of C-terminal PDZ domain interactions by redirecting the transcription factor moiety from the C to the N terminus of a given integral membrane protein thus liberating their native C termini. We successfully applied this "MYTH 2.0" system to five different mammalian full-length renal transporters and identified novel PDZ domain-containing partners of the phosphate (NaPi-IIa) and sulfate (NaS1) transporters that would have otherwise not been detectable. Furthermore this assay was applied to locate the PDZ-binding domain on the NaS1 protein. We showed that the PDZ-binding domain for PDZK1 on NaS1 is upstream of its C terminus, whereas the two interacting proteins, NHERF-1 and NHERF-2, bind at a location closer to the N terminus of NaS1. Moreover NHERF-1 and NHERF-2 increased functional sulfate uptake in Xenopus oocytes when co-expressed with NaS1. Finally we used MYTH 2.0 to demonstrate that the NaPi-IIa transporter homodimerizes via protein-protein interactions within the lipid bilayer. In summary, our study establishes the MYTH 2.0 system as a novel tool for interactive proteomics studies of membrane protein complexes.

Evolutionary demography of long-lived monocarpic perennials: a time-lagged integral projection model

1. The evolution of flowering strategies (when and at what size to flower) in monocarpic perennials is determined by balancing current reproduction with expected future reproduction, and these are largely determined by size-specific patterns of growth and survival. However, because of the difficulty in following long-lived individuals throughout their lives, this theory has largely been tested using short-lived species (< 5 years). 2. Here, we tested this theory using the long-lived monocarpic perennial Campanula thyrsoides which can live up to 16 years. We used a novel approach that combined permanent plot and herb chronology data from a 3-year field study to parameterize and validate integral projection models (IPMs). 3. Similar to other monocarpic species, the rosette leaves of C. thyrsoides wither over winter and so size cannot be measured in the year of flowering. We therefore extended the existing IPM framework to incorporate an additional time delay that arises because flowering demography must be predicted from rosette size in the year before flowering. 4. We found that all main demographic functions (growth, survival probability, flowering probability and fecundity) were strongly size-dependent and there was a pronounced threshold size of flowering. There was good agreement between the predicted distribution of flowering ages obtained from the IPMs and that estimated in the field. Mostly, there was good agreement between the IPM predictions and the direct quantitative field measurements regarding the demographic parameters lambda, R-0 and T. We therefore conclude that the model captures the main demographic features of the field populations. 5. Elasticity analysis indicated that changes in the survival and growth function had the largest effect (c. 80%) on lambda and this was considerably larger than in short-lived monocarps. We found only weak selection pressure operating on the observed flowering strategy which was close to the predicted evolutionary stable strategy. 6. Synthesis. The extended IPM accurately described the demography of a long-lived monocarpic perennial using data collected over a relatively short period. We could show that the evolution of flowering strategies in short- and long-lived monocarps seem to follow the same general rules but with a longevity-related emphasis on survival over fecundity.

Rotational integral geometry of tensor valuations

We derive a new rotational Crofton formula for Minkowski tensors. In special cases, this formula gives (1) the rotational average of Minkowski tensors defined on linear subspaces and (2) the functional defined on linear subspaces with rotational average equal to a Minkowski tensor. Earlier results obtained for intrinsic volumes appear now as special cases.

Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

An Integral View of Fast Shocks Around Supernova 1006

Supernova remnants are among the most spectacular examples of astrophysical pistons in our cosmic neighborhood. The gas expelled by the supernova explosion is launched with velocities ~1000 kilometers per second into the ambient, tenuous interstellar medium, producing shocks that excite hydrogen lines. We have used an optical integral-field spectrograph to obtain high-resolution spatial-spectral maps that allow us to study in detail the shocks in the northwestern rim of supernova 1006. The two-component Hα line is detected at 133 sky locations. Variations in the broad line widths and the broad-to-narrow line intensity ratios across tens of atomic mean free paths suggest the presence of suprathermal protons, the potential seed particles for generating high-energy cosmic rays.

Quasi-Self-Dual Exponential Lévy Processes

The important application of semistatic hedging in financial markets naturally leads to the notion of quasi--self-dual processes. The focus of our study is to give new characterizations of quasi--self-duality. We analyze quasi--self-dual Lévy driven markets which do not admit arbitrage opportunities and derive a set of equivalent conditions for the stochastic logarithm of quasi--self-dual martingale models. Since for nonvanishing order parameter two martingale properties have to be satisfied simultaneously, there is a nontrivial relation between the order and shift parameter representing carrying costs in financial applications. This leads to an equation containing an integral term which has to be inverted in applications. We first discuss several important properties of this equation and, for some well-known Lévy-driven models, we derive a family of closed-form inversion formulae.

Justifying induction on modal μ-formulae

We define a rank function for formulae of the propositional modal μ-calculus such that the rank of a fixed point is strictly bigger than the rank of any of its finite approximations. A rank function of this kind is needed, for instance, to establish the collapse of the modal μ-hierarchy over transitive transition systems. We show that the range of the rank function is ωω. Further we establish that the rank is computable by primitive recursion, which gives us a uniform method to generate formulae of arbitrary rank below ωω.

On the similarity of Sturm-Liouville operators with non-Hermitian boundary conditions to self-adjoint and normal operators

We consider one-dimensional Schrödinger-type operators in a bounded interval with non-self-adjoint Robin-type boundary conditions. It is well known that such operators are generically conjugate to normal operators via a similarity transformation. Motivated by recent interests in quasi-Hermitian Hamiltonians in quantum mechanics, we study properties of the transformations and similar operators in detail. In the case of parity and time reversal boundary conditions, we establish closed integral-type formulae for the similarity transformations, derive a non-local self-adjoint operator similar to the Schrödinger operator and also find the associated “charge conjugation” operator, which plays the role of fundamental symmetry in a Krein-space reformulation of the problem.

Finite volume effects in chiral perturbation theory with twisted boundary conditions

We study the effects of a finite cubic volume with twisted boundary conditions on pseudoscalar mesons. We apply Chiral Perturbation Theory in the p-regime and introduce the twist by means of a constant vector field. The corrections of masses, decay constants, pseudoscalar coupling constants and form factors are calculated at next-to-leading order. We detail the derivations and compare with results available in the literature. In some case there is disagreement due to a different treatment of new extra terms generated from the breaking of the cubic invariance. We advocate to treat such terms as renormalization terms of the twisting angles and reabsorb them in the on-shell conditions. We confirm that the corrections of masses, decay constants, pseudoscalar coupling constants are related by means of chiral Ward identities. Furthermore, we show that the matrix elements of the scalar (resp. vector) form factor satisfies the Feynman–Hellman Theorem (resp. the Ward–Takahashi identity). To show the Ward–Takahashi identity we construct an effective field theory for charged pions which is invariant under electromagnetic gauge transformations and which reproduces the results obtained with Chiral Perturbation Theory at a vanishing momentum transfer. This generalizes considerations previously published for periodic boundary conditions to twisted boundary conditions. Another method to estimate the corrections in finite volume are asymptotic formulae. Asymptotic formulae were introduced by Lüscher and relate the corrections of a given physical quantity to an integral of a specific amplitude, evaluated in infinite volume. Here, we revise the original derivation of Lüscher and generalize it to finite volume with twisted boundary conditions. In some cases, the derivation involves complications due to extra terms generated from the breaking of the cubic invariance. We isolate such terms and treat them as renormalization terms just as done before. In that way, we derive asymptotic formulae for masses, decay constants, pseudoscalar coupling constants and scalar form factors. At the same time, we derive also asymptotic formulae for renormalization terms. We apply all these formulae in combination with Chiral Perturbation Theory and estimate the corrections beyond next-to-leading order. We show that asymptotic formulae for masses, decay constants, pseudoscalar coupling constants are related by means of chiral Ward identities. A similar relation connects in an independent way asymptotic formulae for renormalization terms. We check these relations for charged pions through a direct calculation. To conclude, a numerical analysis quantifies the importance of finite volume corrections at next-to-leading order and beyond. We perform a generic Analysis and illustrate two possible applications to real simulations.