Proving Safety with Trace Automata and Bounded Model Checking

Loop under-approximation is a technique that enriches C programs with additional branches that represent the effect of a (limited) range of loop iterations. While this technique can speed up the detection of bugs significantly, it introduces redundant execution traces which may complicate the verification of the program. This holds particularly true for verification tools based on Bounded Model Checking, which incorporate simplistic heuristics to determine whether all feasible iterations of a loop have been considered. We present a technique that uses \emph{trace automata} to eliminate redundant executions after performing loop acceleration. The method reduces the diameter of the program under analysis, which is in certain cases sufficient to allow a safety proof using Bounded Model Checking. Our transformation is precise---it does not introduce false positives, nor does it mask any errors. We have implemented the analysis as a source-to-source transformation, and present experimental results showing the applicability of the technique

Phase Transition in a Random Fragmentation Problem with Applications to Computer Science

We study a fragmentation problem where an initial object of size x is broken into m random pieces provided x>x_0 where x_0 is an atomic cut-off. Subsequently the fragmentation process continues for each of those daughter pieces whose sizes are bigger than x_0. The process stops when all the fragments have sizes smaller than x_0. We show that the fluctuation of the total number of splitting events, characterized by the variance, generically undergoes a nontrivial phase transition as one tunes the branching number m through a critical value m=m_c. For m<m_c, the fluctuations are Gaussian where as for m>m_c they are anomalously large and non-Gaussian. We apply this general result to analyze two different search algorithms in computer science.Comment: 5 pages RevTeX, 3 figures (.eps

Whole lifespan microscopic observation of budding yeast aging through a microfluidic dissection platform

Important insights into aging have been generated with the genetically tractable and short-lived budding yeast. However, it is still impossible today to continuously track cells by high-resolution microscopic imaging (e.g., fluorescent imaging) throughout their entire lifespan. Instead, the field still needs to rely on a 50-y-old laborious and time-consuming method to assess the lifespan of yeast cells and to isolate differentially aged cells for microscopic snapshots via manual dissection of daughter cells from the larger mother cell. Here, we are unique in achieving continuous and high-resolution microscopic imaging of the entire replicative lifespan of single yeast cells. Our microfluidic dissection platform features an optically prealigned single focal plane and an integrated array of soft elastomer-based micropads, used together to allow for trapping of mother cells, removal of daughter cells, monitoring gradual changes in aging, and unprecedented microscopic imaging of the whole aging process. Using the platform, we found remarkable age-associated changes in phenotypes (e.g., that cells can show strikingly differential cell and vacuole morphologies at the moment of their deaths), indicating substantial heterogeneity in cell aging and death. We envision the microfluidic dissection platform to become a major tool in aging research.

On the expressiveness of forwarding in higher-order communication

Abstract. In higher-order process calculi the values exchanged in communications may contain processes. There are only two capabilities for received processes: execution and forwarding. Here we propose a limited form of forwarding: output actions can only communicate the parallel composition of statically known closed processes and processes received through previously executed input actions. We study the expressiveness of a higher-order process calculus featuring this style of communication. Our main result shows that in this calculus termination is decidable while convergence is undecidable.

Minimal size of a barchan dune

Barchans are dunes of high mobility which have a crescent shape and propagate under conditions of unidirectional wind. However, sand dunes only appear above a critical size, which scales with the saturation distance of the sand flux [P. Hersen, S. Douady, and B. Andreotti, Phys. Rev. Lett. {\bf{89,}} 264301 (2002); B. Andreotti, P. Claudin, and S. Douady, Eur. Phys. J. B {\bf{28,}} 321 (2002); G. Sauermann, K. Kroy, and H. J. Herrmann, Phys. Rev. E {\bf{64,}} 31305 (2001)]. It has been suggested by P. Hersen, S. Douady, and B. Andreotti, Phys. Rev. Lett. {\bf{89,}} 264301 (2002) that this flux fetch distance is itself constant. Indeed, this could not explain the proto size of barchan dunes, which often occur in coastal areas of high litoral drift, and the scale of dunes on Mars. In the present work, we show from three dimensional calculations of sand transport that the size and the shape of the minimal barchan dune depend on the wind friction speed and the sand flux on the area between dunes in a field. Our results explain the common appearance of barchans a few tens of centimeter high which are observed along coasts. Furthermore, we find that the rate at which grains enter saltation on Mars is one order of magnitude higher than on Earth, and is relevant to correctly obtain the minimal dune size on Mars.Comment: 11 pages, 10 figure

A class of Calogero type reductions of free motion on a simple Lie group

The reductions of the free geodesic motion on a non-compact simple Lie group G based on the $G_+ \times G_+$ symmetry given by left- and right multiplications for a maximal compact subgroup $G_+ \subset G$ are investigated. At generic values of the momentum map this leads to (new) spin Calogero type models. At some special values the `spin' degrees of freedom are absent and we obtain the standard $BC_n$ Sutherland model with three independent coupling constants from SU(n+1,n) and from SU(n,n). This generalization of the Olshanetsky-Perelomov derivation of the $BC_n$ model with two independent coupling constants from the geodesics on $G/G_+$ with G=SU(n+1,n) relies on fixing the right-handed momentum to a non-zero character of $G_+$. The reductions considered permit further generalizations and work at the quantized level, too, for non-compact as well as for compact G.Comment: shortened to 13 pages in v2 on request of Lett. Math. Phys. and corrected some spelling error

Population pharmacokinetics-based recommendations for a single delayed or missed dose of nusinersen

Nusinersen is an antisense oligonucleotide approved for the treatment of spinal muscular atrophy. The drug is given intrathecally at 12 mg, beginning with 3 loading doses at 2-week intervals, a fourth loading dose 30 days thereafter, and maintenance doses at 4-month intervals. This population pharmacokinetic model was developed to clarify how to maintain targeted nusinersen exposure after an unforeseen one-time delay or missed dose. Simulations demonstrated that the impact of a one-time delay in dosing or a missed dose on median cerebrospinal fluid exposures depended on duration of interruption and the regimen phase in which it occurred. Delays in loading doses delayed reaching the peak trough concentration by approximately the duration of the interruption. Resumption of the regimen as soon as possible resulted in achieving steady state trough concentration upon completion of the loading phase. A short delay (30-90 days) during the maintenance phase led to prolonged lower median cerebrospinal fluid concentration if all subsequent doses were shifted by the same 4-month interval. However, administration of the delayed dose, followed by the subsequent dose as originally scheduled, rapidly restored trough concentration. If a dose must be delayed, patients should return to the original dosing schedule as soon as possible

Calculation of Band Edge Eigenfunctions and Eigenvalues of Periodic Potentials through the Quantum Hamilton - Jacobi Formalism

We obtain the band edge eigenfunctions and the eigenvalues of solvable periodic potentials using the quantum Hamilton - Jacobi formalism. The potentials studied here are the Lam{\'e} and the associated Lam{\'e} which belong to the class of elliptic potentials. The formalism requires an assumption about the singularity structure of the quantum momentum function $p$, which satisfies the Riccati type quantum Hamilton - Jacobi equation, $p^{2} -i \hbar \frac{d}{dx}p = 2m(E- V(x))$ in the complex $x$ plane. Essential use is made of suitable conformal transformations, which leads to the eigenvalues and the eigenfunctions corresponding to the band edges in a simple and straightforward manner. Our study reveals interesting features about the singularity structure of $p$, responsible in yielding the band edge eigenfunctions and eigenvalues.Comment: 21 pages, 5 table

Classification of integrable Weingarten surfaces possessing an sl(2)-valued zero curvature representation

In this paper we classify Weingarten surfaces integrable in the sense of soliton theory. The criterion is that the associated Gauss equation possesses an sl(2)-valued zero curvature representation with a nonremovable parameter. Under certain restrictions on the jet order, the answer is given by a third order ordinary differential equation to govern the functional dependence of the principal curvatures. Employing the scaling and translation (offsetting) symmetry, we give a general solution of the governing equation in terms of elliptic integrals. We show that the instances when the elliptic integrals degenerate to elementary functions were known to nineteenth century geometers. Finally, we characterize the associated normal congruences