Search

We present a model of systems of cells with intracellular oscillators (`clocks'). This is motivated by examples from developmental biology and from the behavior of organisms on the threshold to multicellularity. Cells undergo random motion and adhere to each other. The adhesion strength between neighbors depends on their clock phases in addition to a constant baseline strength. The oscillators are linked via Kuramoto-type local interactions. The model is an advection-diffusion partial differential equation with nonlocal advection terms. We demonstrate that synchronized states correspond to Dirac-delta measure solutions of a weak version of the equation. To analyze the complex interplay of aggregation and synchronization, we then perform a linear stability analysis of the incoherent, spatially uniform state. This lets us classify possibly emerging patterns depending on model parameters. Combining these results with numerical simulations, we determine a range of possible far-from equilibrium patterns when baseline adhesion strength is zero: There is aggregation into separate synchronized clusters with or without global synchrony; global synchronization without aggregation; or unexpectedly a ``phase wave" pattern characterized by spatial gradients of clock phases. A 2D Lattice-Gas Cellular Automaton model confirms and illustrates these results.

We consider the geometric optics problem of finding a system of two reflectors that transform a spherical wavefront into a beam of parallel rays with prescribed intensity distribution. Using techniques from optimal transportation theory, it has been shown before that this problem is equivalent to an infinite dimensional linear programming (LP) problem. We investigate techniques for constructing the two reflectors numerically. A straightforward discretization of this problem has the disadvantage that the number of constraints increases rapidly with the mesh size. So with this technique only very coarse meshes are practical. To address this well-known issue we propose an iterative solution scheme. In each step an LP problem is solved. Information from the previous iteration step is used to reduce the number of constraints necessary. As a proof of concept we apply our proposed scheme to solve a problem with synthetic data. We give evidence that the scheme converges. We also show that it allows for much finer meshes than a simple discretization scheme. There exists a growing literature for the application of optimal transportation theory to other beam shaping problems. Our proposed scheme is easy to adapt for these problems as well., Comment: submitted

We consider the following geometric optics problem: Construct a system of two reflectors which transforms a spherical wavefront generated by a point source into a beam of parallel rays. This beam has a prescribed intensity distribution. We give a rigorous analysis of this problem. The reflectors we construct are (parts of) the boundaries of convex sets. We prove existence of solutions for a large class of input data and give a uniqueness result. To the author's knowledge, this is the first time that a rigorous mathematical analysis of this problem is given. The approach is based on optimal transportation theory. It yields a practical algorithm for finding the reflectors. Namely, the problem is equivalent to a constrained linear optimization problem., Comment: 5 Figures - pdf files attached to submission, but not shown in manuscript

The Cellular Potts Model (CPM) has been used for simulating various biological phenomena such as differential adhesion, fruiting body formation of the slime mold Dictyostelium discoideum, angiogenesis, cancer invasion, chondrogenesis in embryonic vertebrate limbs, and many others. In this paper, we derive continuous limit of discrete one dimensional CPM with the chemotactic interactions between cells in the form of a Fokker-Planck equation for the evolution of the cell probability density function. This equation is then reduced to the classical macroscopic Keller-Segel model. In particular, all coefficients of the Keller-Segel model are obtained from parameters of the CPM. Theoretical results are verified numerically by comparing Monte Carlo simulations for the CPM with numerics for the Keller-Segel model., Comment: 15 pages, 7 figures

It is shown that the problem of designing a two-reflector system transforming a plane wave front with given intensity into an output plane front with prescribed output intensity can be formulated and solved as the Monge-Kantorovich mass transfer problem., Comment: 25 pages, 2 figures