Lattice Boltzmann method for simulating the viscous flow in large distensible blood vessels

A lattice Boltzmann method for simulating the viscous flow in large distensible blood vessels is presented by introducing a boundary condition for elastic and moving boundaries. The mass conservation for the boundary condition is tested in detail. The viscous flow in elastic vessels is simulated with a pressure-radius relationship similar to that of the Pulmonary blood vessels. The numerical results for steady flow agree with the analytical prediction to very high accuracy, and the simulation results for pulsatile flow are comparable with those of the aortic flows observed experimentally. The model is expected to find many applications for studying blood flows in large distensible arteries, especially in those suffering from atherosclerosis. stenosis. aneurysm, etc.

Studying the contact point and interface moving in a sinusoidal tube with lattice Boltzmann method

The multicomponent nonideal gas lattice Boltzmann model by Shan and Chen (S-C) is used to study the immiscible displacement in a sinusoidal tube. The movement of interface and the contact point (contact line in three-dimension) is studied. Due to the roughness of the boundary, the contact point shows "stick-slip" mechanics. The "stick-slip" effect decreases as the speed of the interface increases. For fluids that are nonwetting, the interface is almost perpendicular to the boundaries at most time, although its shapes at different position of the tube are rather different. When the tube becomes narrow, the interface turns a complex curves rather than remains simple menisci. The velocity is found to vary considerably between the neighbor nodes close to the contact point, consistent with the experimental observation that the velocity is multi-values on the contact line. Finally, the effect of three boundary conditions is discussed. The average speed is found different for different boundary conditions. The simple bounce-back rule makes the contact point move fastest. Both the simple bounce-back and the no-slip bounce-back rules are more sensitive to the roughness of the boundary in comparison with the half-way bounce-back rule. The simulation results suggest that the S-C model may be a promising tool in simulating the displacement behaviour of two immiscible fluids in complex geometry.

Lattice dynamics study of microstructures of electrorheological fluids

The lattice dynamics method is used to study the stability of the chain structures formed in electrorheological (ER) fluids. The appearance of the soft modes in the phonon dispersion of the structures indicates that the chains tend to distort and aggregate into thicker columns due to the electrostatic attractive forces and thermal generated forces between them. The results show that the stability of the chains relies on their width and the separation between them. The complete chain structures are more stable than the chains with defects. The results can be used to elucidate the densification phenomenon of the chains in the structuring process of ER fluids in the quiescent state.

Lattice Boltzmann simulation of viscous fluid systems with elastic boundaries

A lattice Boltzmann model able to simulate viscous fluid systems with elastic and movable boundaries is proposed. By introducing the virtual distribution function at the boundary, the Galilean invariance is recovered for the full system. As examples of application, the how in elastic vessels is simulated with the pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady how are in good agreement with the analytical prediction, while the simulation results for pulsative how agree with the experimental observation of the aortic flows qualitatively. The approach has potential application in the study of the complex fluid systems such as the suspension system as well as the arterial blood flow.

Reply to “Comment on ‘Temperature-dependent orientational ordering on a spherical surface modeled with a lattice spin model’ ”

A reply to the comment of S. Romano, Phys. Rev. E 2015 on our previous paper is provided.

Thiol and Sulfenic Acid Oxidation of AhpE, the One-Cysteine Peroxiredoxin from Mycobacterium tuberculosis: Kinetics, Acidity Constants, and Conformational Dynamics

Drug resistance and virulence of Mycobacterium tuberculosis are partially related to the pathogen`s antioxidant systems. Peroxide detoxification in this bacterium is achieved by the heme-containing catalase peroxidase and different two-cysteine peroxiredoxins. M. tuberculosis genome also codifies for a putative one-cysteine peroxiredoxin, alkyl hydroperoxide reductase E (MtAhpE). Its expression was previously demonstrated at a transcriptional level, and the crystallographic structure of the recombinant protein was resolved under reduced and oxidized states. Herein, we report that the conformation of MtAhpE changed depending on its single cysteine redox state, as reflected by different tryptophan fluorescence properties and changes in quaternary structure. Dynamics of fluorescence changes, complemented by competition kinetic assays, were used to perform protein functional studies. MtAhE reduced peroxynitrite 2 orders of magnitude faster than hydrogen peroxide (1.9 x 10(7) M(-1) s(-1) vs 8.2 x 10(4) M(-1) s(-1) at pH 7.4 and 25 degrees C, respectively). The latter also caused cysteine overoxidation to sulfinic acid, but at much slower rate constant (40 M(-1) s(-1)). The pK(a) of the thiol in the reduced enzyme was 5.2, more than one unit lower than that of the sulfenic acid in the oxidized enzyme. The pH profile of hydrogen peroxide-mediated thiol and sulfenic acid oxidations indicated thiolate and sulfenate as the reacting species. The formation of sulfenic acid as well as the catalytic peroxidase activity of MtAhpE was demonstrated using the artificial reducing substrate thionitrobenzoate. Taken together, our results indicate that MtAhpE is a relevant component in the antioxidant repertoire of M. tuberculosis probably involved in peroxide and specially peroxynitrite detoxification.

Chaotic synchronization in 2D lattice for scene segmentation

Synchronization and chaos play important roles in neural activities and have been applied in oscillatory correlation modeling for scene and data analysis. Although it is an extensively studied topic, there are still few results regarding synchrony in locally coupled systems. In this paper we give a rigorous proof to show that large numbers of coupled chaotic oscillators with parameter mismatch in a 2D lattice can be synchronized by providing a sufficiently large coupling strength. We demonstrate how the obtained result can be applied to construct an oscillatory network for scene segmentation. (C) 2007 Elsevier B.V. All rights reserved.

Clustering and diffusion in a symplectic map lattice with non-local coupling

We study a symplectic chain with a non-local form of coupling by means of a standard map lattice where the interaction strength decreases with the lattice distance as a power-law, in Such a way that one can pass continuously from a local (nearest-neighbor) to a global (mean-field) type of coupling. We investigate the formation of map clusters, or spatially coherent structures generated by the system dynamics. Such clusters are found to be related to stickiness of chaotic phase-space trajectories near periodic island remnants, and also to the behavior of the diffusion coefficient. An approximate two-dimensional map is derived to explain some of the features of this connection. (C) 2008 Elsevier Ltd. All rights reserved.

A strongly attractive Fermi gas in an optical lattice

We study strongly attractive fermions in an optical lattice superimposed by a trapping potential. We calculate the densities of fermions and condensed bound molecules at zero temperature. There is a competition between dissociated fermions and molecules leading to a reduction of the density of fermions at the trap center. (C) 2010 Elsevier B.V. All rights reserved.

Stochastic lattice gas model describing the dynamics of the SIRS epidemic process

We study a stochastic process describing the onset of spreading dynamics of an epidemic in a population composed of individuals of three classes: susceptible (S), infected (I), and recovered (R). The stochastic process is defined by local rules and involves the following cyclic process: S -> I -> R -> S (SIRS). The open process S -> I -> R (SIR) is studied as a particular case of the SIRS process. The epidemic process is analyzed at different levels of description: by a stochastic lattice gas model and by a birth and death process. By means of Monte Carlo simulations and dynamical mean-field approximations we show that the SIRS stochastic lattice gas model exhibit a line of critical points separating the two phases: an absorbing phase where the lattice is completely full of S individuals and an active phase where S, I and R individuals coexist, which may or may not present population cycles. The critical line, that corresponds to the onset of epidemic spreading, is shown to belong in the directed percolation universality class. By considering the birth and death process we analyze the role of noise in stabilizing the oscillations. (C) 2009 Elsevier B.V. All rights reserved.

CRITICAL BEHAVIOR AND THRESHOLD OF COEXISTENCE OF A PREDATOR-PREY STOCHASTIC MODEL IN A 2D LATTICE

We investigate the critical behavior of a stochastic lattice model describing a predator-prey system. By means of Monte Carlo procedure we simulate the model defined on a regular square lattice and determine the threshold of species coexistence, that is, the critical phase boundaries related to the transition between an active state, where both species coexist and an absorbing state where one of the species is extinct. A finite size scaling analysis is employed to determine the order parameter, order parameter fluctuations, correlation length and the critical exponents. Our numerical results for the critical exponents agree with those of the directed percolation universality class. We also check the validity of the hyperscaling relation and present the data collapse curves.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Conservative ensembles for nonequilibrium lattice-gas systems

We show how to set up a constant particle ensemble for the steady state of nonequilibrium lattice-gas systems which originally are defined on a constant rate ensemble. We focus on nonequilibrium systems in which particles are created and annihilated on the sites of a lattice and described by a master equation. We consider also the case in which a quantity other than the number of particle is conserved. The conservative ensembles can be useful in the study of phase transitions and critical phenomena particularly discontinuous phase transitions.

Ab initio study of the electronic and optical properties of sillimanite (Al(2)SiO(5)) crystal

Ab initio calculations based on the density functional theory (DFT) are used to investigate the electronic and optical properties of sillimanite. The geometrical parameters of the unit cell, which contain 32 atoms, have been fully optimized and are in good agreement with the experimental data. The electronic structure shows that sillimanite has an indirect band gap of 5.18 eV. The complex dielectric function and optical constants, such as extinction coefficient, refractive index, reflectivity and energy-loss spectrum, are calculated. The optical properties of sillimanite are discussed based on the band structure calculations. It is shown that the O-2p states and Al-3s, Si-3s states play the major role in optical transitions as initial and final states, respectively. (C) 2011 Elsevier B.V. All rights reserved.

Linear covariant gauges on the lattice

Linear covariant gauges, such as Feynman gauge, are very useful in perturbative calculations. Their non-perturbative formulation is, however, highly non-trivial. In particular, it is a challenge to define linear covariant gauges on a lattice. We consider a class of gauges in lattice gauge theory that coincides with the perturbative definition of linear covariant gauges in the formal continuum limit. The corresponding gauge-fixing procedure is described and analyzed in detail, with an application to the pure SU(2) case. In addition, results for the gluon propagator in the two-dimensional case are given. (C) 2008 Elsevier B.V. All rights reserved.