Reversibility of Red blood Cell deformation

The ability of cells to undergo reversible shape changes is often crucial to their survival. For Red Blood Cells (RBCs), irreversible alteration of the cell shape and flexibility often causes anemia. Here we show theoretically that RBCs may react irreversibly to mechanical perturbations because of tensile stress in their cytoskeleton. The transient polymerization of protein fibers inside the cell seen in sickle cell anemia or a transient external force can trigger the formation of a cytoskeleton-free membrane protrusion of micrometer dimensions. The complex relaxation kinetics of the cell shape is shown to be responsible for selecting the final state once the perturbation is removed, thereby controlling the reversibility of the deformation. In some case, tubular protrusion are expected to relax via a peculiar "pearling instability".Comment: 4 pages, 3 figure

Limiting shapes of confined lipid vesicles

We theoretically study the shapes of lipid vesicles confined to a spherical cavity, elaborating a framework based on the so-called limiting shapes constructed from geometrically simple structural elements such as double-membrane walls and edges. Partly inspired by numerical results, the proposed non-compartmentalized and compartmentalized limiting shapes are arranged in the bilayer-couple phase diagram which is then compared to its free-vesicle counterpart. We also compute the area-difference-elasticity phase diagram of the limiting shapes and we use it to interpret shape transitions experimentally observed in vesicles confined within another vesicle. The limiting-shape framework may be generalized to theoretically investigate the structure of certain cell organelles such as the mitochondrion

Dynamics of Fluid Vesicles in Oscillatory Shear Flow

The dynamics of fluid vesicles in oscillatory shear flow was studied using differential equations of two variables: the Taylor deformation parameter and inclination angle $\theta$. In a steady shear flow with a low viscosity $\eta_{\rm {in}}$ of internal fluid, the vesicles exhibit steady tank-treading motion with a constant inclination angle $\theta_0$. In the oscillatory flow with a low shear frequency, $\theta$ oscillates between $\pm \theta_0$ or around $\theta_0$ for zero or finite mean shear rate $\dot\gamma_{\rm m}$, respectively. As shear frequency $f_{\gamma}$ increases, the vesicle oscillation becomes delayed with respect to the shear oscillation, and the oscillation amplitude decreases. At high $f_{\gamma}$ with $\dot\gamma_{\rm m}=0$, another limit-cycle oscillation between $\theta_0-\pi$ and $-\theta_0$ is found to appear. In the steady flow, $\theta$ periodically rotates (tumbling) at high $\eta_{\rm {in}}$, and $\theta$ and the vesicle shape oscillate (swinging) at middle $\eta_{\rm {in}}$ and high shear rate. In the oscillatory flow, the coexistence of two or more limit-cycle oscillations can occur for low $f_{\gamma}$ in these phases. For the vesicle with a fixed shape, the angle $\theta$ rotates back to the original position after an oscillation period. However, it is found that a preferred angle can be induced by small thermal fluctuations.Comment: 11 pages, 13 figure

Transport coefficients of off-lattice mesoscale-hydrodynamics simulation techniques

The viscosity and self-diffusion constant of particle-based mesoscale hydrodynamic methods, multi-particle collision dynamics (MPC) and dissipative particle dynamics (DPD), are investigated, both with and without angular-momentum conservation. Analytical results are derived for fluids with an ideal-gas equation of state and a finite-time-step dynamics, and compared with simulation data. In particular, the viscosity is derived in a general form for all variants of the MPC method. In general, very good agreement between theory and simulations is obtained.Comment: 12 pages, 9 figure

Magnetorotational Instability in Liquid Metal Couette Flow

Despite the importance of the magnetorotational instability (MRI) as a fundamental mechanism for angular momentum transport in magnetized accretion disks, it has yet to be demonstrated in the laboratory. A liquid sodium alpha-omega dynamo experiment at the New Mexico Institute of Mining and Technology provides an ideal environment to study the MRI in a rotating metal annulus (Couette flow). A local stability analysis is performed as a function of shear, magnetic field strength, magnetic Reynolds number, and turbulent Prandtl number. The later takes into account the minimum turbulence induced by the formation of an Ekman layer against the rigidly rotating end walls of a cylindrical vessel. Stability conditions are presented and unstable conditions for the sodium experiment are compared with another proposed MRI experiment with liquid gallium. Due to the relatively large magnetic Reynolds number achievable in the sodium experiment, it should be possible to observe the excitation of the MRI for a wide range of wavenumbers and further to observe the transition to the turbulent state.Comment: 12 pages, 22 figures, 1 table. To appear in the Astrophysical Journa

One-Dimensional Confinement and Enhanced Jahn-Teller Instability in LaVO3_3

Ordering and quantum fluctuations of orbital degrees of freedom are studied theoretically for LaVO$_3$ in spin-C-type antiferromagnetic state. The effective Hamiltonian for the orbital pseudospin shows strong one-dimensional anisotropy due to the negative interference among various exchange processes. This significantly enhances the instability toward lattice distortions for the realistic estimate of the Jahn-Teller coupling by first-principle LDA+$U$ calculations, instead of favoring the orbital singlet formation. This explains well the experimental results on the anisotropic optical spectra as well as the proximity of the two transition temperatures for spin and orbital orderings.Comment: 4 pages including 4 figure

Particle-Based Mesoscale Hydrodynamic Techniques

Dissipative particle dynamics (DPD) and multi-particle collision (MPC) dynamics are powerful tools to study mesoscale hydrodynamic phenomena accompanied by thermal fluctuations. To understand the advantages of these types of mesoscale simulation techniques in more detail, we propose new two methods, which are intermediate between DPD and MPC -- DPD with a multibody thermostat (DPD-MT), and MPC-Langevin dynamics (MPC-LD). The key features are applying a Langevin thermostat to the relative velocities of pairs of particles or multi-particle collisions, and whether or not to employ collision cells. The viscosity of MPC-LD is derived analytically, in very good agreement with the results of numerical simulations.Comment: 7 pages, 2 figures, 1 tabl

Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids

In this review, we describe and analyze a mesoscale simulation method for fluid flow, which was introduced by Malevanets and Kapral in 1999, and is now called multi-particle collision dynamics (MPC) or stochastic rotation dynamics (SRD). The method consists of alternating streaming and collision steps in an ensemble of point particles. The multi-particle collisions are performed by grouping particles in collision cells, and mass, momentum, and energy are locally conserved. This simulation technique captures both full hydrodynamic interactions and thermal fluctuations. The first part of the review begins with a description of several widely used MPC algorithms and then discusses important features of the original SRD algorithm and frequently used variations. Two complementary approaches for deriving the hydrodynamic equations and evaluating the transport coefficients are reviewed. It is then shown how MPC algorithms can be generalized to model non-ideal fluids, and binary mixtures with a consolute point. The importance of angular-momentum conservation for systems like phase-separated liquids with different viscosities is discussed. The second part of the review describes a number of recent applications of MPC algorithms to study colloid and polymer dynamics, the behavior of vesicles and cells in hydrodynamic flows, and the dynamics of viscoelastic fluids

On particle acceleration and trapping by Poynting flux dominated flows

Using particle-in-cell (PIC) simulations, we study the evolution of a strongly magnetized plasma slab propagating into a finite density ambient medium. Like previous work, we find that the slab breaks into discrete magnetic pulses. The subsequent evolution is consistent with diamagnetic relativistic pulse acceleration of \cite{liangetal2003}. Unlike previous work, we use the actual electron to proton mass ratio and focus on understanding trapping vs. transmission of the ambient plasma by the pulses and on the particle acceleration spectra. We find that the accelerated electron distribution internal to the slab develops a double-power law. We predict that emission from reflected/trapped external electrons will peak after that of the internal electrons. We also find that the thin discrete pulses trap ambient electrons but allow protons to pass through, resulting in less drag on the pulse than in the case of trapping of both species. Poynting flux dominated scenarios have been proposed as the driver of relativistic outflows and particle acceleration in the most powerful astrophysical jets.Comment: 25 pages, Accepted by Plasma Physics and Controlled Fusio

On the nature of the FBS blue stellar objects and the completeness of the Bright Quasar Survey. II

In Paper I (Mickaelian et al. 1999), we compared the surface density of QSOs in the Bright Quasar Survey (BQS) and in the First Byurakan Survey (FBS) and concluded that the completeness of the BQS is of the order of 70% rather than 30-50% as suggested by several authors. A number of new observations recently became available, allowing a re-evaluation of this completeness. We now obtain a surface density of QSOs brighter than B = 16.16 in a subarea of the FBS covering ~2250 deg^2, equal to 0.012 deg^-2 (26 QSOs), implying a completeness of 53+/-10%.Comment: LaTeX 2e, 11 pages, 3 tables and 3 figures (included in text). To appear in Astrophysics. Uses a modified aaspp4.sty (my_aaspp4.sty), included in packag