A pseudo-time EnKF incorporating shape based reconstruction for diffuse optical tomography

Purpose: The authors aim at developing a pseudo-time, sub-optimal stochastic filtering approach based on a derivative free variant of the ensemble Kalman filter (EnKF) for solving the inverse problem of diffuse optical tomography (DOT) while making use of a shape based reconstruction strategy that enables representing a cross section of an inhomogeneous tumor boundary by a general closed curve. Methods: The optical parameter fields to be recovered are approximated via an expansion based on the circular harmonics (CH) (Fourier basis functions) and the EnKF is used to recover the coefficients in the expansion with both simulated and experimentally obtained photon fluence data on phantoms with inhomogeneous inclusions. The process and measurement equations in the pseudo-dynamic EnKF (PD-EnKF) presently yield a parsimonious representation of the filter variables, which consist of only the Fourier coefficients and the constant scalar parameter value within the inclusion. Using fictitious, low-intensity Wiener noise processes in suitably constructed ``measurement'' equations, the filter variables are treated as pseudo-stochastic processes so that their recovery within a stochastic filtering framework is made possible. Results: In our numerical simulations, we have considered both elliptical inclusions (two inhomogeneities) and those with more complex shapes (such as an annular ring and a dumbbell) in 2-D objects which are cross-sections of a cylinder with background absorption and (reduced) scattering coefficient chosen as mu(b)(a)=0.01mm(-1) and mu('b)(s)=1.0mm(-1), respectively. We also assume mu(a) = 0.02 mm(-1) within the inhomogeneity (for the single inhomogeneity case) and mu(a) = 0.02 and 0.03 mm(-1) (for the two inhomogeneities case). The reconstruction results by the PD-EnKF are shown to be consistently superior to those through a deterministic and explicitly regularized Gauss-Newton algorithm. We have also estimated the unknown mu(a) from experimentally gathered fluence data and verified the reconstruction by matching the experimental data with the computed one. Conclusions: The PD-EnKF, which exhibits little sensitivity against variations in the fictitiously introduced noise processes, is also proven to be accurate and robust in recovering a spatial map of the absorption coefficient from DOT data. With the help of shape based representation of the inhomogeneities and an appropriate scaling of the CH expansion coefficients representing the boundary, we have been able to recover inhomogeneities representative of the shape of malignancies in medical diagnostic imaging. (C) 2012 American Association of Physicists in Medicine. [DOI: 10.1118/1.3679855]

Statistical properties of entropy-consuming fluctuations in jammed states of laponite suspensions: Fluctuation relations and generalized Gumbel distribution

We report a universal large deviation behavior of spatially averaged global injected power just before the rejuvenation of the jammed state formed by an aging suspension of laponite clay under an applied stress. The probability distribution function (PDF) of these entropy consuming strongly non-Gaussian fluctuations follow an universal large deviation functional form described by the generalized Gumbel (GG) distribution like many other equilibrium and nonequilibrium systems with high degree of correlations but do not obey the Gallavotti-Cohen steady-state fluctuation relation (SSFR). However, far from the unjamming transition (for smaller applied stresses) SSFR is satisfied for both Gaussian as well as non-Gaussian PDF. The observed slow variation of the mean shear rate with system size supports a recent theoretical prediction for observing GG distribution.

Protein-Protein Interactions in Clathrin Vesicular Assembly: Radial Distribution of Evolutionary Constraints in Interfaces

In eukaryotic organisms clathrin-coated vesicles are instrumental in the processes of endocytosis as well as intracellular protein trafficking. Hence, it is important to understand how these vesicles have evolved across eukaryotes, to carry cargo molecules of varied shapes and sizes. The intricate nature and functional diversity of the vesicles are maintained by numerous interacting protein partners of the vesicle system. However, to delineate functionally important residues participating in protein-protein interactions of the assembly is a daunting task as there are no high-resolution structures of the intact assembly available. The two cryoEM structures closely representing intact assembly were determined at very low resolution and provide positions of C alpha atoms alone. In the present study, using the method developed by us earlier, we predict the protein-protein interface residues in clathrin assembly, taking guidance from the available low-resolution structures. The conservation status of these interfaces when investigated across eukaryotes, revealed a radial distribution of evolutionary constraints, i.e., if the members of the clathrin vesicular assembly can be imagined to be arranged in spherical manner, the cargo being at the center and clathrins being at the periphery, the detailed phylogenetic analysis of these members of the assembly indicated high-residue variation in the members of the assembly closer to the cargo while high conservation was noted in clathrins and in other proteins at the periphery of the vesicle. This points to the strategy adopted by the nature to package diverse proteins but transport them through a highly conserved mechanism.

On the solution of bivariate population balance equations for aggregation: Pivotwise expansion of solution domain

The solution of a bivariate population balance equation (PBE) for aggregation of particles necessitates a large 2-d domain to be covered. A correspondingly large number of discretized equations for particle populations on pivots (representative sizes for bins) are solved, although at the end only a relatively small number of pivots are found to participate in the evolution process. In the present work, we initiate solution of the governing PBE on a small set of pivots that can represent the initial size distribution. New pivots are added to expand the computational domain in directions in which the evolving size distribution advances. A self-sufficient set of rules is developed to automate the addition of pivots, taken from an underlying X-grid formed by intersection of the lines of constant composition and constant particle mass. In order to test the robustness of the rule-set, simulations carried out with pivotwise expansion of X-grid are compared with those obtained using sufficiently large fixed X-grids for a number of composition independent and composition dependent aggregation kernels and initial conditions. The two techniques lead to identical predictions, with the former requiring only a fraction of the computational effort. The rule-set automatically reduces aggregation of particles of same composition to a 1-d problem. A midway change in the direction of expansion of domain, effected by the addition of particles of different mean composition, is captured correctly by the rule-set. The evolving shape of a computational domain carries with it the signature of the aggregation process, which can be insightful in complex and time dependent aggregation conditions. (c) 2012 Elsevier Ltd. All rights reserved.

Distribution of transverse chain fluctuations in harmonically confined semiflexible polymers

Two different experimental studies of polymer dynamics based on single-molecule fluorescence imaging have recently found evidence of heterogeneities in the widths of the putative tubes that surround filaments of F-actin during their motion in concentrated solution. In one J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchesner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)], the observations were explained in terms of the statistics of a worm-like chain confined to a potential determined self-consistently by a binary collision approximation, and in the other B. Wang, J. Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010)], they were explained in terms of the scaling properties of a random fluid of thin rods. In this paper, we show, using an exact path integral calculation, that the distribution of the length-averaged transverse fluctuations of a harmonically confined weakly bendable rod (one possible realization of a semiflexible chain in a tube), is in good qualitative agreement with the experimental data, although it is qualitatively different in analytic structure from the earlier theoretical predictions. We also show that similar path integral techniques can be used to obtain an exact expression for the time correlation function of fluctuations in the tube cross section. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4712306]

The force distribution function of an oscillator model of polymer stretching at constant velocity

Recent simulations of the stretching of tethered biopolymers at a constant speed v (Ponmurugan and Vemparala, 2011 Phys. Rev. E 84 060101(R)) have suggested that for any time t, the distribution of the fluctuating forces f responsible for chain deformation is governed by a relation of the form P(+ f)/ P(- f) = expgamma f], gamma being a coefficient that is solely a function of v and the temperature T. This result, which is reminiscent of the fluctuation theorems applicable to stochastic trajectories involving thermodynamic variables, is derived in this paper from an analytical calculation based on a generalization of Mazonka and Jarzynski's classic model of dragged particle dynamics Mazonka and Jarzynski, 1999 arXiv:cond-\textbackslashmat/9912121v1]. However, the analytical calculations suggest that the result holds only if t >> 1 and the force fluctuations are driven by white rather than colored noise; they further suggest that the coefficient gamma in the purported theorem varies not as v(0.15)T-(0.7), as indicated by the simulations, but as vT(-1).

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.

Bimodal distribution of motility and cell fate in Dictyostelium discoideum

Pre-starvation amoebae of Dictyostelium discoideum exhibit random movements. Starved cells aggregate by directed movements (chemotaxis) towards cyclic AMP and differentiate into live spores or dead stalk cells. Many differences between presumptive spore and stalk cells precede differentiation. We have examined whether cell motility-related factors are also among them. Cell speeds and localisation of motility-related signalling molecules were monitored by live cell imaging and immunostaining (a) in nutrient medium during growth, (b) immediately following transfer to starvation medium and (c) in nutrient medium that was re-introduced after a brief period of starvation. Cells moved randomly under all three conditions but mean speeds increased following transfer from nutrient medium to starvation medium; the transition occurred within 15 min. The distribution of speeds in starvation medium was bimodal: about 20% of the cells moved significantly faster than the remaining 80%. The motility-related molecules F-actin, PTEN and PI3 kinase were distributed differently in slow and fast cells. Among starved cells, the calcium content of slower cells was lower than that of the faster cells. All differences reverted within 15 min after restoration of the nutrient medium. The slow/fast distinction was missing in Polysphondylium pallidum, a cellular slime mould that lacks the presumptive stalk and spore cell classes, and in the trishanku (triA(center dot)) mutant of D. discoideum, in which the classes exist but are unstable. The transition from growth to starvation triggers a spontaneous and reversible switch in the distribution of D. discoideum cell speeds. Cells whose calcium content is relatively low (known to be presumptive spore cells) move slower than those whose calcium levels are higher (known to be presumptive stalk cells). Slow and fast cells show different distributions of motility-related proteins. The switch is indicative of a bistable mechanism underlying cell motility.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

Shape Effect on Electronic and Photovoltaic Properties of CdS Nanocrystals

Changes in electronic and photovoltaic properties of semiconductor nanocrystals predominantly due to changes in shape are discussed here. Cadmium sulfide (CdS) semiconductor nanocrystals of various shapes (tetrapod, tetrahedron, sphere and rod) obtained using an optimized solvothermal process exhibited a mixed cubic (zinc blende) and hexagonal (wurtzite) crystal structure. The simultaneous presence of the two crystal phases in varying amounts is observed to play a pivotal role in determining both the electronic and photovoltaic properties of the CdS nanocrystals. Light to electrical energy conversion efficiencies (measured in two-electrode configuration laboratory solar cells) remarkably decreased by one order in magnitude from tetrapod -> tetrahedron -> sphere -> rod. The tetrapod-CdS nanocrystals, which displayed the highest light to electrical energy conversion efficiency, showed a favorable shift in position of the conduction band edge leading to highest rate of electron injection (from CdS nanocrystal to the wide band gap semiconductor viz, titanium dioxide, TiO2) and lowest rate of electron-hole recombination (higher free electron lifetimes).

A mixture model approach for formant tracking and the robustness of student's-t distribution

We address the problem of robust formant tracking in continuous speech in the presence of additive noise. We propose a new approach based on mixture modeling of the formant contours. Our approach consists of two main steps: (i) Computation of a pyknogram based on multiband amplitude-modulation/frequency-modulation (AM/FM) decomposition of the input speech; and (ii) Statistical modeling of the pyknogram using mixture models. We experiment with both Gaussian mixture model (GMM) and Student's-t mixture model (tMM) and show that the latter is robust with respect to handling outliers in the pyknogram data, parameter selection, accuracy, and smoothness of the estimated formant contours. Experimental results on simulated data as well as noisy speech data show that the proposed tMM-based approach is also robust to additive noise. We present performance comparisons with a recently developed adaptive filterbank technique proposed in the literature and the classical Burg's spectral estimator technique, which show that the proposed technique is more robust to noise.

A broad pore size distribution mesoporous SnO2 as anode for lithium-ion batteries

We demonstrate here that mesoporous tin dioxide (abbreviated M-SnO2) with a broad pore size distribution can be a prospective anode in lithium-ion batteries. M-SnO2 with pore size ranging between 2 and 7.5 nm was synthesized using a hydrothermal procedure involving two different surfactants of slightly different sizes, and characterized. The irreversible capacity loss that occurs during the first discharge and charge cycle is 890 mAh g(-1), which is smaller than the 1,010-mAh g(-1) loss recorded for mesoporous SnO2 (abbreviated S-SnO2) synthesized using a single surfactant. After 50 cycles, the discharge capacity of M-SnO2 (504 mAh g(-1)) is higher than that of S-SnO2 (401 mAh g(-1)) and solid nanoparticles of SnO2 (abbreviated nano-SnO2 < 4 mAh g(-1)) and nano-SnO2. Transmission electron microscopy revealed higher disorder in the pore arrangement in M-SnO2. This, in turn imparts lower stiffness to M-SnO2 (elastic modulus, E (R) a parts per thousand aEuro parts per thousand 14.5 GPa) vis-a-vis S-SnO2 (E (R) a parts per thousand aEuro parts per thousand 20.5 GPa), as obtained using the nanoindentation technique. Thus, the superior battery performance of M-SnO2 is attributed to its intrinsic material mechanical property. The fluidity of the internal microstructure of M-SnO2 resulted in a lower degree of aggregation of Sn particles compared to S-SnO2 and nano-SnO2 structural stabilization and long-term cyclability.

Potential and electric field distribution in a ceramic disc insulator string with faulty insulators

Ceramic/Porcelain suspension disc insulators are widely used in power systems to provide electrical insulation and mechanically support for high-voltage transmission lines. These insulators are subjected to a variety of stresses, including mechanical, electrical and environmental. These stresses act in unison. The exact nature and magnitude of these stresses vary significantly and depends on insulator design, application and its location. Due to various reasons the insulator disc can lose its electrical insulation properties without any noticeable mechanical failure. Such a condition while difficult to recognize, can enhance the stress on remaining healthy insulator discs in the string further may lead to a flashover. To understand the stress enhancement due to faulty discs in a string, attempt has been made to simulate the potential and electric field profiles for various disc insulators presently used in the country. The results of potential and electric filed stress obtained for normal and strings with faulty insulator discs are presented.

A note on the inverse mode shape problem for bars, beams, and plates

Motivated by the idea of designing a structure for a desired mode shape, intended towards applications such as resonant sensors, actuators and vibration confinement, we present the inverse mode shape problem for bars, beams and plates in this work. The objective is to determine the cross-sectional profile of these structures, given a mode shape, boundary condition and the mass. The contribution of this article is twofold: (i) A numerical method to solve this problem when a valid mode shape is provided in the finite element framework for both linear and nonlinear versions of the problem. (ii) An analytical result to prove the uniqueness and existence of the solution in the case of bars. This article also highlights a very important question of the validity of a mode shape for any structure of given boundary conditions.

A Model for the Temperature Distribution in Resin Impregnated Paper Bushings

Resin impregnated paper (RIP) is a relatively new insulation system recommended for the use in transformer bushings. In the recent past, RIP has acquired prominence as insulation in bushings, over conventional oil impregnated paper (OIP), in view of its overwhelming advantages the more important among them being low dielectric loss and possibility for positioning the bushing at any desired angle over the transformer. In addition, the fact that such systems do not pose problems of fire hazard is counted as a very important consideration. The disadvantage of RIP compared to OIP, however, is its much higher cost and involved manufacturing process. The temperature rise in RIP bushings under normal operating conditions is seen to be a difficult parameter to control in view of the limited options for effective cooling. It is therefore essential to take serious note of this aspect, to arrest rapid deterioration of bushing. The degradation of dry-type insulation such as RIP is often due to thermal stress. The long time performance thereof, depends strongly, on the maximum operating temperature. With this in view, the Authors have developed a theoretical model and computational method to study the temperature distribution in the body of insulation. The Authors consider that the basis for the model as being the temperature and electric stress aided AC conductivity. The ensuing heat balance (continuity) equations in 2-D cylindrical geometry are treated as a Dirichelet-Neumann boundary value problem.