In Silico Prediction of the Mechanobiological Response of Arterial Tissue: Application to Angioplasty and Stenting

One way to restore physiological blood flow to occluded arteries involves the deformation of plaque using an intravascular balloon and preventing elastic recoil using a stent. Angioplasty and stent implantation cause unphysiological loading of the arterial tissue, which may lead to tissue in-growth and reblockage; termed “restenosis.” In this paper, a computational methodology for predicting the time-course of restenosis is presented. Stress-induced damage, computed using a remaining life approach, stimulates inflammation (production of matrix degrading factors and growth stimuli). This, in turn, induces a change in smooth muscle cell phenotype from contractile (as exists in the quiescent tissue) to synthetic (as exists in the growing tissue). In this paper, smooth muscle cell activity (migration, proliferation, and differentiation) is simulated in a lattice using a stochastic approach to model individual cell activity. The inflammation equations are examined under simplified loading cases. The mechanobiological parameters of the model were estimated by calibrating the model response to the results of a balloon angioplasty study in humans. The simulation method was then used to simulate restenosis in a two dimensional model of a stented artery. Cell activity predictions were similar to those observed during neointimal hyperplasia, culminating in the growth of restenosis. Similar to experiment, the amount of neointima produced increased with the degree of expansion of the stent, and this relationship was found to be highly dependant on the prescribed inflammatory response. It was found that the duration of inflammation affected the amount of restenosis produced, and that this effect was most pronounced with large stent expansions. In conclusion, the paper shows that the arterial tissue response to mechanical stimulation can be predicted using a stochastic cell modeling approach, and that the simulation captures features of restenosis development observed with real stents. The modeling approach is proposed for application in three dimensional models of cardiovascular stenting procedures.

Received signal strength in large-scale wireless relay sensor network: a stochastic ray approach

The authors consider a point percolation lattice representation of a large-scale wireless relay sensor network (WRSN) deployed in a cluttered environment. Each relay sensor corresponds to a grid point in the random lattice and the signal sent by the source is modelled as an ensemble of photons that spread in the space, which may 'hit' other sensors and are 'scattered' around. At each hit, the relay node forwards the received signal to its nearest neighbour through direction-selective relaying. The authors first derive the distribution that a relay path reaches a prescribed location after undergoing certain number of hops. Subsequently, a closed-form expression of the average received signal strength (RSS) at the destination can be computed as the summation of all signal echoes' energy. Finally, the effect of the anomalous diffusion exponent ß on the mean RSS in a WRSN is studied, for which it is found that the RSS scaling exponent e is given by (3ß-1)/ß. The results would provide useful insight into the design and deployment of large-scale WRSNs in future. © 2011 The Institution of Engineering and Technology.

Critical assessment of two-qubit post-Markovian master equations

A post-Markovian master equation has been recently proposed as a tool to describe the evolution of a system coupled to a memory-keeping environment [A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101 ( R) ( 2005)]. For a single qubit affected by appropriately chosen environmental conditions, the corresponding dynamics is always legitimate and physical. Here we extend such a situation to the case of two qubits, only one of which experiences the environmental effects. We show how, despite the innocence of such an extension, the introduction of the second qubit should be done cum grano salis to avoid consequences such as the breaking of the positivity of the associated dynamical map. This hints at the necessity of using care when adopting phenomenologically derived models for evolutions occurring outside the Markovian framework.

Cold-formed steel sections with web openings subjected to web crippling under two-flange loading conditions - Part II: Parametric study and proposed design equations

A parametric study of cold-formed steel sections with web openings subjected to web crippling was undertaken using finite element analysis, to investigate the effects of web holes and cross-section sizes on the web crippling strengths of channel sections subjected to web crippling under both interior-two-flange (ITF) and end-two-flange (ETF) loading conditions. In both loading conditions, the hole was centred beneath the bearing plate. It was demonstrated that the main factors influencing the web crippling strength are the ratio of the hole depth to the flat depth of the web, and the ratio of the length of bearing plates to the flat depth of the web. In this paper, design recommendations in the form of web crippling strength reduction factors are proposed, that are conservative to both the experimental and finite element results.

Reliability Assessment of 3-D Space Frame Structures Applying Stochastic Finite Element Analysis

Throughout design development of satellite structure, stress engineer is usually challenged with randomness in applied loads and material properties. To overcome such problem, a risk-based design is applied which estimates satellite structure probability of failure under static and thermal loads. Determining probability of failure can help to update initially applied factors of safety that were used during structure preliminary design phase. These factors of safety are related to the satellite mission objective. Sensitivity-based analysis is to be implemented in the context of finite element analysis (probabilistic finite element method or stochastic finite element method (SFEM)) to determine the probability of failure for satellite structure or one of its components.

R-matrix-Floquet theory of multiphoton processes:VIII. A linear equations method

A new linear equations method for calculating the R-matrix, which arises in the R-matrix-Floquet theory of multiphoton processes, is introduced. This method replaces the diagonalization of the Floquet Hamiltonian matrix by the solution of a set of linear simultaneous equations which are solved, in the present work, by the conjugate gradient method. This approach uses considerably less computer memory and can be readily ported onto parallel computers. It will thus enable much larger problems of current interest to be treated. This new method is tested by applying it to three-photon ionization of helium at frequencies where double resonances with a bound state and autoionizing states are important. Finally, an alternative linear equations method, which avoids the explicit calculation of the R-matrix by incorporating the boundary conditions directly, is described in an appendix.

Regularized Online Mixture of Gaussians for Background Subtraction

Mixture of Gaussians (MoG) modelling [13] is a popular approach to background subtraction in video sequences. Although the algorithm shows good empirical performance, it lacks theoretical justification. In this paper, we give a justification for it from an online stochastic expectation maximization (EM) viewpoint and extend it to a general framework of regularized online classification EM for MoG with guaranteed convergence. By choosing a special regularization function, l1 norm, we derived a new set of updating equations for l1 regularized online MoG. It is shown empirically that l1 regularized online MoG converge faster than the original online MoG .

Relative role of deterministic and stochastic determinants of soil animal community: a spatially explicit analysis of oribatid mites

1. Ecologists are debating the relative role of deterministic and stochastic determinants of community structure. Although the high diversity and strong spatial structure of soil animal assemblages could provide ecologists with an ideal ecological scenario, surprisingly little information is available on these assemblages.
2. We studied species-rich soil oribatid mite assemblages from a Mediterranean beech forest and a grassland. We applied multivariate regression approaches and analysed spatial autocorrelation at multiple spatial scales using Moran's eigenvectors. Results were used to partition community variance in terms of the amount of variation uniquely accounted for by environmental correlates (e.g. organic matter) and geographical position. Estimated neutral diversity and immigration parameters were also applied to a soil animal group for the first time to simulate patterns of community dissimilarity expected under neutrality, thereby testing neutral predictions.
3. After accounting for spatial autocorrelation, the correlation between community structure and key environmental parameters disappeared: about 40% of community variation consisted of spatial patterns independent of measured environmental variables such as organic matter. Environmentally independent spatial patterns encompassed the entire range of scales accounted for by the sampling design (from tens of cm to 100 m). This spatial variation could be due to either unmeasured but spatially structured variables or stochastic drift mediated by dispersal. Observed levels of community dissimilarity were significantly different from those predicted by neutral models.
4. Oribatid mite assemblages are dominated by processes involving both deterministic and stochastic components and operating at multiple scales. Spatial patterns independent of the measured environmental variables are a prominent feature of the targeted assemblages, but patterns of community dissimilarity do not match neutral predictions. This suggests that either niche-mediated competition or environmental filtering or both are contributing to the core structure of the community. This study indicates new lines of investigation for understanding the mechanisms that determine the signature of the deterministic component of animal community assembly.

Stochastic and deterministic processes interact in the assembly of desert microbial communities on a global scale

Extreme arid regions in the worlds' major deserts are typified by quartz pavement terrain. Cryptic hypolithic communities colonize the ventral surface of quartz rocks and this habitat is characterized by a relative lack of environmental and trophic complexity. Combined with readily identifiable major environmental stressors this provides a tractable model system for determining the relative role of stochastic and deterministic drivers in community assembly. Through analyzing an original, worldwide data set of 16S rRNA-gene defined bacterial communities from the most extreme deserts on the Earth, we show that functional assemblages within the communities were subject to different assembly influences. Null models applied to the photosynthetic assemblage revealed that stochastic processes exerted most effect on the assemblage, although the level of community dissimilarity varied between continents in a manner not always consistent with neutral models. The heterotrophic assemblages displayed signatures of niche processes across four continents, whereas in other cases they conformed to neutral predictions. Importantly, for continents where neutrality was either rejected or accepted, assembly drivers differed between the two functional groups. This study demonstrates that multi-trophic microbial systems may not be fully described by a single set of niche or neutral assembly rules and that stochasticity is likely a major determinant of such systems, with significant variation in the influence of these determinants on a global scale.

Stochastic particle acceleration in the lobes of giant radio galaxies

We investigate the acceleration of particles by Alfven waves via the second-order Fermi process in the lobes of giant radio galaxies. Such sites are candidates for the accelerators of ultra-high-energy cosmic rays (UHECR). We focus on the nearby Fanaroff-Riley type I radio galaxy Centaurus A. This is motivated by the coincidence of its position with the arrival direction of several of the highest energy Auger events. The conditions necessary for consistency with the acceleration time-scales predicted by quasi-linear theory are reviewed. Test particle calculations are performed in fields which guarantee electric fields with no component parallel to the local magnetic field. The results of quasi-linear theory are, to an order of magnitude, found to be accurate at low turbulence levels for non-relativistic Alfven waves and at both low and high turbulence levels in the mildly relativistic case. We conclude that for pure stochastic acceleration via Alfven waves to be plausible as the generator of UHECR in Cen A, the baryon number density would need to be several orders of magnitude below currently held upper limits.

Pressure effect on vibrational frequency and dephasing of 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids

Raman spectra in the range of the totally symmetric stretching mode of the [PF6]− anion, νs(PF6), have been measured for 1-alkyl-3-methylimidazolium ionic liquids [CnC1im][PF6], for n = 4, 6, and 8, as a function of pressure at room temperature. The ionic liquids [C6C1im][PF6] and [C8C1im][PF6] remain in an amorphous phase up to 3.5 GPa, in contrast to [C4C1im][PF6], whichcrystallizes above ∼0.5 GPa. Equations of state based either on a group contribution model or Carnahan-Starling-van der Waals model have been used to estimate the densities of the ionic liquids at high pressures. The shifts of the vibrational frequency of νs(PF6) with density observed in [C6C1im][PF6] and in [C8C1im][PF6] have been calculated by a hard-sphere model of a pseudo-diatomic solute under short-range repulsive interactions with the neighboring particles. The stochastic model of Kubo for vibrational dephasing has been used to obtain the amplitude of vibrational frequency fluctuation, ⟨Δω 2⟩, and the relaxation time of frequency fluctuation, τ c , as a function of density by Raman band shape analysis of the νs(PF6) mode of [C6C1im][PF6] and [C8C1im][PF6].