Calculating the Fickian diffusivity for a lattice-based random walk with agents and obstacles of different shapes and sizes

Random walk models are often used to interpret experimental observations of the motion of biological cells and molecules. A key aim in applying a random walk model to mimic an in vitro experiment is to estimate the Fickian diffusivity (or Fickian diffusion coefficient),D. However, many in vivo experiments are complicated by the fact that the motion of cells and molecules is hindered by the presence of obstacles. Crowded transport processes have been modeled using repeated stochastic simulations in which a motile agent undergoes a random walk on a lattice that is populated by immobile obstacles. Early studies considered the most straightforward case in which the motile agent and the obstacles are the same size. More recent studies considered stochastic random walk simulations describing the motion of an agent through an environment populated by obstacles of different shapes and sizes. Here, we build on previous simulation studies by analyzing a general class of lattice-based random walk models with agents and obstacles of various shapes and sizes. Our analysis provides exact calculations of the Fickian diffusivity, allowing us to draw conclusions about the role of the size, shape and density of the obstacles, as well as examining the role of the size and shape of the motile agent. Since our analysis is exact, we calculateDdirectly without the need for random walk simulations. In summary, we find that the shape, size and density of obstacles has a major influence on the exact Fickian diffusivity. Furthermore, our results indicate that the difference in diffusivity for symmetric and asymmetric obstacles is significant.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Size effect on the lattice parameter of KCl during mechanical milling

The size effect on the lattice parameter of ionic KCl nanocrystals was studied systematically during mechanical milling of Pure KCl powder under vacuum. The results suggest anomalous lattice expansion, with the lattice parameter increasing from 6.278 angstrom at d = 6 mu m to 6.30307 angstrom at d = 85 mn. The defects generated during ball milling of KCl and surface stress are deemed to be responsible for this lattice parameter expansion. (C) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Proton Spin-Lattice Relaxation in Sodium Ammonium Selenate Dihydrate:NaNHaSeOa.2He0

Proton spin-lattice relaxation studies in sodium ammonium selenate dihydrate carried out in the temperature range 130 to 300 K at 10 MHz show a continuous change in T, at T, indicating a second order phase transition. This compound is a typical case of a highly hindered solid wherein the thermally activated reorientations of ammonium ions freeze well above 77 K, as seen by NMR.Untersuchimgen der Protonen-Spin-Gitter-Relaxation in Natriuni-Ammoniumselenat-Dihydrat bei 10 MHz im Temperaturbereich 130 bis 300 K zeigen eine kontinuierliche Andernng in TI bei T, und ergeben einen Phasenubergang zweiter Art. Diese Verbindung ist ein typischer Fall eines stark ,,behinderten" Festkarpers, in dein die thermisch aktivierten Reorientierungen der Ammoniumionen weit oberhalb 77 H einfrieren, wie die NMR-Ergebnisse zeigen.

Distortion of Dimension under Quasiconformal Mappings

Quasiconformal mappings are natural generalizations of conformal mappings. They are homeomorphisms with 'bounded distortion' of which there exist several approaches. In this work we study dimension distortion properties of quasiconformal mappings both in the plane and in higher dimensional Euclidean setting. The thesis consists of a summary and three research articles. A basic property of quasiconformal mappings is the local Hölder continuity. It has long been conjectured that this regularity holds at the Sobolev level (Gehring's higher integrabilty conjecture). Optimal regularity would also provide sharp bounds for the distortion of Hausdorff dimension. The higher integrability conjecture was solved in the plane by Astala in 1994 and it is still open in higher dimensions. Thus in the plane we have a precise description how Hausdorff dimension changes under quasiconformal deformations for general sets. The first two articles contribute to two remaining issues in the planar theory. The first one concerns distortion of more special sets, for rectifiable sets we expect improved bounds to hold. The second issue consists of understanding distortion of dimension on a finer level, namely on the level of Hausdorff measures. In the third article we study flatness properties of quasiconformal images of spheres in a quantitative way. These also lead to nontrivial bounds for their Hausdorff dimension even in the n-dimensional case.

Characterization of Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease, reveals a unique mode of DNA-binding, helical distortion and cleavage, compared to a canonical LAGLIDADG homing endonuclease

Mycobacterium leprae, which has undergone reductive evolution leaving behind a minimal set of essential genes, has retained intervening sequences in four of its genes implicating a vital role for them in the survival of the leprosy bacillus. A single in-frame intervening sequence has been found embedded within its recA gene. Comparison of M. leprae recA intervening sequence with the known intervening sequences indicated that it has the consensus amino acid sequence necessary for being a LAGLIDADG-type homing endonuclease. In light of massive gene decay and function loss in the leprosy bacillus, we sought to investigate whether its recA intervening sequence encodes a catalytically active homing endonuclease. Here we show that the purified M. leprae RecA intein (PI-MleI) binds to cognate DNA and displays endonuclease activity in the presence of alternative divalent cations, Mg2+ or Mn2+. A combination of approaches including four complementary footprinting assays such as DNase I, Cu/phenanthroline, methylation protection and KMnO4, enhancement of 2-aminopurine fluorescence and mapping of the cleavage site revealed that PI-MleI binds to cognate DNA flanking its insertion site, induces helical distortion at the cleavage site and generates two staggered double-strand breaks. Taken together, these results implicate that PI-MleI possess a modular structure with separate domains for DNA target recognition and cleavage, each with distinct sequence preferences. From a biological standpoint, it is tempting to speculate that our findings have implications for understanding the evolution of LAGLIDADG family of homing endonucleases

Electron spin resonance study of copper(II) doped ferroelectric ammonium sulphate

Single crystal electron spin resonance studies of Cu2+ doped ferroelectric ammonium sulphate ((NH4)2SO4, Tc = 223 K) are reported at 300 and 77 K. The Cu2+ ion is found to enter the lattice interstitially with a trigonal bipyramidal coordination. Proton superhyperfine interaction is found for magnetic field directions close to the a-axis. Changes are observed in the 77 K recordings indicating a distortion of the trigonal bipyramid consistent with crystal structure data. An increase of the proton superhyperfine constant in the ferroelectric phase is indicative of stronger hydrogen bonding. The Cu2+ ion doped as an impurity in a trigonal bipyramid environment in a diamagnetic host lattice is reported for the first time.

Pyrrolidine ring conformations in prolyl peptides from 13C spin-lattice relaxation times

A method was developed in the framework of a bistable jump model to obtain the pyrrolidine ring conformations in proline peptides from 13C spin-lattice relaxation times. Equations are presented expressing the ring torsions in terms of the 13C spin-lattice relaxation times of the ring carbons. This method was applied to 26 pyrrolidine ring systems and acceptable conformations were obtained.

Rapid distortion of axisymmetric turbulence

A generalization of the isotropic theory of Batchelor & Proudman (1954) is developed to estimate the effect of sudden but arbitrary three-dimensional distortion on homogeneous, initially axisymmetric turbulence. The energy changes due to distortion are expressed in terms of the Fourier coefficients of an expansion in zonal harmonics of the two independent scalar functions that describe the axisymmetric spectral tensor. However, for two special but non-trivial forms of this tensor, which represent possibly the simplest kinds of non-isotropic turbulence and specify the angular distribution but not the wavenumber dependence, the energy ratios have been determined in closed form. The deviation of the ratio from its isotropic value is the product of a factor containing R, the initial value of the ratio of the longitudinal to the transverse energy component, and another factor depending only on the geometry of the distortion. It is found that, in axisymmetric and large two-dimensional contractions, the isotropic theory gives nearly the correct longitudinal energy, but (when R > 1) over-estimates the increase in the transverse energy; the product of the two intensities varies little unless the distortion is very large, thus accounting for the stress-freezing observed in rapidly accelerated shear flows.Comparisons with available experimental data for the spectra and for the energy ratios show reasonable agreement. The different ansatzes predict results in broad qualitative agreement with a simple strategem suggested by Reynolds & Tucker (1975), but the quantitative differences are not always negligible.

Large Low-Frequency Fluctuations in the Velocity of a Driven Vortex Lattice in a Single Crystal of 2H-NbSe2 Superconductor

The driven state of a well-ordered flux line lattice in a single crystal of 2H-NbSe2 in the time domain has revealed the presence of substantial fluctuations in velocity, with large and distinct time periods (similar to seconds). A superposition of a periodic drive in the driven vortex lattice causes distinct changes in these fluctuations. We propose that prior to the onset of the peak effect there exists a heretofore unexplored regime of coherent dynamics, with unexpected behavior in velocity fluctuations.

Origin of activation of Lattice Oxygen and Synergistic Interaction in Bimetal-Ionic Ce0.89Fe0.1Pd0.01O2-delta Catalyst

Flourite-type nanocrystalline Ce0.9Fe0.1O2-delta and Ce0.89Fe0.1Pd0.01O2-delta solid solutions have been synthesized by solution combustion method,'.which show higher oxygen storage/release property (OSC) compared to CeO2 and Ce0.8Zr0.2O2. Temperature programmed reduction an XPS study reveal that the presence of Pd ion in Ce0.9Fe0.1O2-delta facilitates complete reduction of Fe3+ to Fe2+ state and partial reduction of Ce4+ to Ce3+ state at.temperatures as low as 105 degrees C compared to 400 degrees C for monometal-ionic Ce0.9Fe0.1O2-delta. Fe3+ ion is reduced to Fe2+ and not to Feo due to favorable redox potential for Ce4+ + Fe2+ -> Ce3+ + Fe3+ reaction. Using first-principles density functional theory calculation we determine M-O (M = Pd, Fe, Ce) bond lengths, and find that bond lengths vary from shorter (2.16 angstrom) to longer (2.9 angstrom) bond distances compared to mean Ce-O bond distance of 2.34 angstrom. for CeO2. Using these results in bond valence analysis, we show that oxygen with bond valences as low as -1.55 are created, leading to activation of lattice oxygen in the bimetal ionic catalyst. Temperatures of CO oxidation and NO reduction by CO/H-2 are lower with the bimetalionic Ce0.89Fe0.1Pd0.01O2-delta catalyst compared to monometal-ionic Ce0.9Fe0.1O2-delta and Ce0.99Pd0.01O2-delta catalysts. From XPS studies of Pd impregnated on CeO2 and Fe2O3 oxides, we show that the synergism leading to low temperature activation of lattice oxygen in bimetal-ionic catalyst Ce0.89Fe0.1Pd0.01O2-delta is due to low-temperature reduction of Pd2+ to Pd-0, followed by Pd-0 + 2Fe(3+) -> Pd2+ + 2Fe(2+), Pd-0 + 2Ce(4+) -> Pd2+ + 2Ce(3+) redox reaction.

Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

Structural Investigation of Activated Lattice Oxygen in Ce1-xSnxO2 and Ce1-x-ySnxPdyO2-delta by EXAFS and DFT calculation

Substitution of Sn4+ ion in CeO2 creates activated oxygen in Ce0.8Sn0.2O2 leading to higher oxygen storage capacity compared to Ce0.8Zr0.2O2. With Pd ion substitution in Ce0.8Sn0.2O2,activation of oxygen is further enhanced as observed from the H-2/TPR study. Both EXAFS analysis and DFT calculation reveal that in the solid solution Ceexhibits 4 + 4 coordination, Sri exhibits 4 + 2 + 2 coordination and Pd has 4 + 3 coordination. While the oxygen in the First four coordination with short M-O bonds are strongly held in the lattice, the oxygens in the second and higher coordinations with long M-O bonds are weakly bound, and they are the activated oxygen ill the lattice. Bond valence analysis shows that oxygen with valencies as low its 1.65 are created by the Sn and Pd ion Substitution. Another interesting observation is that H-2/TPR experiment of Ce1-xSnxO2 shows a broad peak starting from 200 to 500 degrees C, while the same reduction is achieved in a single step at similar to 110 degrees C in presence Pd2+ on. Substitution of Pd2+ ion thus facilitates synergistic reduction of the catalyst at lower temperature. We have shown that simultaneous reduction of the Ce4+ and Sr4+ ions by Pd-0 is the synergistic interaction leading to high oxygen storage capacity at low temperature.

Quantum phase analysis of 1D superconducting quantum dot lattice using extended Bose-Hubbard model

We have obtained the quantum phase diagram of a one-dimensional superconducting quantum dot lattice using the extended Bose-Hubbard model for different commensurabilities. We describe the nature of different quantum phases at the charge degeneracy point. We find a direct phase transition from the Mott insulating phase to the superconducting phase for integer band fillings of Cooper pairs. We predict explicitly the presence of two kinds of repulsive Luttinger liquid phases, besides the charge density wave and superconducting phases for half-integer band fillings. We also predict that extended range interactions are necessary to obtain the correct phase boundary of a one-dimensional interacting Cooper system. We have used the density matrix renormalization group method and Abelian bosonization to study our system.