Monte Carlo simulations of fluid vesicles with in plane orientational ordering

We present a method for simulating fluid vesicles with in-plane orientational ordering. The method involves computation of local curvature tensor and parallel transport of the orientational field on a randomly triangulated surface. It is shown that the model reproduces the known equilibrium conformation of fluid membranes and work well for a large range of bending rigidities. Introduction of nematic ordering leads to stiffening of the membrane. Nematic ordering can also result in anisotropic rigidity on the surface leading to formation of membrane tubes.Comment: 11 Pages, 12 Figures, To appear in Phys. Rev.

Self-organized stable pacemakers near the onset of birhythmicity

General amplitude equations for reaction-diffusion systems near to the soft onset of birhythmicity described by a supercritical pitchfork-Hopf bifurcation are derived. Using these equations and applying singular perturbation theory, we show that stable autonomous pacemakers represent a generic kind of spatiotemporal patterns in such systems. This is verified by numerical simulations, which also show the existence of breathing and swinging pacemaker solutions. The drift of self-organized pacemakers in media with spatial parameter gradients is analytically and numerically investigated.Comment: 4 pages, 4 figure

Role of disclinations in determining the morphology of deformable fluid interfaces

We study the equilibrium shapes of vesicles, with an in-plane nematic order, using a Monte-Carlo scheme and show that highly curved shapes, like tubes and discs, with a striking similarity to the structures engendered by certain curvature sensing peripheral membrane proteins, can be spontaneously generated by anisotropic directional curvature with nematic disclinations playing and important role. We show that the coupling between nematic order and local curvature could lead to like defects moving towards each other and unlike defects moving away, in turn leading to tube formation. Thermally induced defect pair production lead to branched tubular structures. It is also shown that helical arrangement of the membrane tubes, with nematic field spiraling around it, is a dominant soft mode of the system.Comment: 6 Figures; Soft Matter, Advance Article 201

Exact Solutions for Domain Walls in Coupled Complex Ginzburg - Landau Equations

The complex Ginzburg Landau equation (CGLE) is a ubiquitous model for the evolution of slowly varying wave packets in nonlinear dissipative media. A front (shock) is a transient layer between a plane-wave state and a zero background. We report exact solutions for domain walls, i.e., pairs of fronts with opposite polarities, in a system of two coupled CGLEs, which describe transient layers between semi-infinite domains occupied by each component in the absence of the other one. For this purpose, a modified Hirota bilinear operator, first proposed by Bekki and Nozaki, is employed. A novel factorization procedure is applied to reduce the intermediate calculations considerably. The ensuing system of equations for the amplitudes and frequencies is solved by means of computer-assisted algebra. Exact solutions for mutually-locked front pairs of opposite polarities, with one or several free parameters, are thus generated. The signs of the cubic gain/loss, linear amplification/attenuation, and velocity of the coupled-front complex can be adjusted in a variety of configurations. Numerical simulations are performed to study the stability properties of such fronts.Comment: Journal of the Physical Society of Japan, in pres

Spectral plots and the representation and interpretation of biological data

It is basic question in biology and other fields to identify the char- acteristic properties that on one hand are shared by structures from a particular realm, like gene regulation, protein-protein interaction or neu- ral networks or foodwebs, and that on the other hand distinguish them from other structures. We introduce and apply a general method, based on the spectrum of the normalized graph Laplacian, that yields repre- sentations, the spectral plots, that allow us to find and visualize such properties systematically. We present such visualizations for a wide range of biological networks and compare them with those for networks derived from theoretical schemes. The differences that we find are quite striking and suggest that the search for universal properties of biological networks should be complemented by an understanding of more specific features of biological organization principles at different scales.Comment: 15 pages, 7 figure

Three-dimensional pattern formation, multiple homogeneous soft modes, and nonlinear dielectric electroconvection

Patterns forming spontaneously in extended, three-dimensional, dissipative systems are likely to excite several homogeneous soft modes ($\approx$ hydrodynamic modes) of the underlying physical system, much more than quasi one- and two-dimensional patterns are. The reason is the lack of damping boundaries. This paper compares two analytic techniques to derive the patten dynamics from hydrodynamics, which are usually equivalent but lead to different results when applied to multiple homogeneous soft modes. Dielectric electroconvection in nematic liquid crystals is introduced as a model for three-dimensional pattern formation. The 3D pattern dynamics including soft modes are derived. For slabs of large but finite thickness the description is reduced further to a two-dimensional one. It is argued that the range of validity of 2D descriptions is limited to a very small region above threshold. The transition from 2D to 3D pattern dynamics is discussed. Experimentally testable predictions for the stable range of ideal patterns and the electric Nusselt numbers are made. For most results analytic approximations in terms of material parameters are given.Comment: 29 pages, 2 figure

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

For a given subspace, the Rayleigh-Ritz method projects the large quadratic eigenvalue problem (QEP) onto it and produces a small sized dense QEP. Similar to the Rayleigh-Ritz method for the linear eigenvalue problem, the Rayleigh-Ritz method defines the Ritz values and the Ritz vectors of the QEP with respect to the projection subspace. We analyze the convergence of the method when the angle between the subspace and the desired eigenvector converges to zero. We prove that there is a Ritz value that converges to the desired eigenvalue unconditionally but the Ritz vector converges conditionally and may fail to converge. To remedy the drawback of possible non-convergence of the Ritz vector, we propose a refined Ritz vector that is mathematically different from the Ritz vector and is proved to converge unconditionally. We construct examples to illustrate our theory.Comment: 20 page

Instances and connectors : issues for a second generation process language

This work is supported by UK EPSRC grants GR/L34433 and GR/L32699Over the past decade a variety of process languages have been defined, used and evaluated. It is now possible to consider second generation languages based on this experience. Rather than develop a second generation wish list this position paper explores two issues: instances and connectors. Instances relate to the relationship between a process model as a description and the, possibly multiple, enacting instances which are created from it. Connectors refers to the issue of concurrency control and achieving a higher level of abstraction in how parts of a model interact. We believe that these issues are key to developing systems which can effectively support business processes, and that they have not received sufficient attention within the process modelling community. Through exploring these issues we also illustrate our approach to designing a second generation process language.Postprin