Robustness and applicability of Markov chain Monte Carlo algorithms for eigenvalue problems

In this paper we analyse applicability and robustness of Markov chain Monte Carlo algorithms for eigenvalue problems. We restrict our consideration to real symmetric matrices. Almost Optimal Monte Carlo (MAO) algorithms for solving eigenvalue problems are formulated. Results for the structure of both - systematic and probability error are presented. It is shown that the values of both errors can be controlled independently by different algorithmic parameters. The results present how the systematic error depends on the matrix spectrum. The analysis of the probability error is presented. It shows that the close (in some sense) the matrix under consideration is to the stochastic matrix the smaller is this error. Sufficient conditions for constructing robust and interpolation Monte Carlo algorithms are obtained. For stochastic matrices an interpolation Monte Carlo algorithm is constructed. A number of numerical tests for large symmetric dense matrices are performed in order to study experimentally the dependence of the systematic error from the structure of matrix spectrum. We also study how the probability error depends on the balancing of the matrix. (c) 2007 Elsevier Inc. All rights reserved.

A neurofuzzy network knowledge extraction and extended gram-Schmidt algorithm for model subspace decomposition

This paper introduces a new neurofuzzy model construction and parameter estimation algorithm from observed finite data sets, based on a Takagi and Sugeno (T-S) inference mechanism and a new extended Gram-Schmidt orthogonal decomposition algorithm, for the modeling of a priori unknown dynamical systems in the form of a set of fuzzy rules. The first contribution of the paper is the introduction of a one to one mapping between a fuzzy rule-base and a model matrix feature subspace using the T-S inference mechanism. This link enables the numerical properties associated with a rule-based matrix subspace, the relationships amongst these matrix subspaces, and the correlation between the output vector and a rule-base matrix subspace, to be investigated and extracted as rule-based knowledge to enhance model transparency. The matrix subspace spanned by a fuzzy rule is initially derived as the input regression matrix multiplied by a weighting matrix that consists of the corresponding fuzzy membership functions over the training data set. Model transparency is explored by the derivation of an equivalence between an A-optimality experimental design criterion of the weighting matrix and the average model output sensitivity to the fuzzy rule, so that rule-bases can be effectively measured by their identifiability via the A-optimality experimental design criterion. The A-optimality experimental design criterion of the weighting matrices of fuzzy rules is used to construct an initial model rule-base. An extended Gram-Schmidt algorithm is then developed to estimate the parameter vector for each rule. This new algorithm decomposes the model rule-bases via an orthogonal subspace decomposition approach, so as to enhance model transparency with the capability of interpreting the derived rule-base energy level. This new approach is computationally simpler than the conventional Gram-Schmidt algorithm for resolving high dimensional regression problems, whereby it is computationally desirable to decompose complex models into a few submodels rather than a single model with large number of input variables and the associated curse of dimensionality problem. Numerical examples are included to demonstrate the effectiveness of the proposed new algorithm.

Investigation of liquid MALDI and optimization for instrument tuning and quantitative measurements

The success of Matrix-assisted laser desorption / ionisation (MALDI) in fields such as proteomics has partially but not exclusively been due to the development of improved data acquisition and sample preparation techniques. This has been required to overcome some of the short comings of the commonly used solid-state MALDI matrices such as - cyano-4-hydroxycinnamic acid (CHCA) and 2,5-dihydroxybenzoic acid (DHB). Solid state matrices form crystalline samples with highly inhomogeneous topography and morphology which results in large fluctuations in analyte signal intensity from spot to spot and positions within the spot. This means that efficient tuning of the mass spectrometer can be impeded and the use of MALDI MS for quantitative measurements is severely impeded. Recently new MALDI liquid matrices have been introduced which promise to be an effective alternative to crystalline matrices. Generally the liquid matrices comprise either ionic liquid matrices (ILMs) or a usually viscous liquid matrix which is doped with a UV lightabsorbing chromophore [1-3]. The advantages are that the droplet surface is smooth and relatively uniform with the analyte homogeneously distributed within. They have the ability to replenish a sampling position between shots negating the need to search for sample hot-spots. Also the liquid nature of the matrix allows for the use of additional additives to change the environment to which the analyte is added.

Power properties if invariant tests for spatial autocorrelation in linear regression

This paper derives some exact power properties of tests for spatial autocorrelation in the context of a linear regression model. In particular, we characterize the circumstances in which the power vanishes as the autocorrelation increases, thus extending the work of Krämer (2005). More generally, the analysis in the paper sheds new light on how the power of tests for spatial autocorrelation is affected by the matrix of regressors and by the spatial structure. We mainly focus on the problem of residual spatial autocorrelation, in which case it is appropriate to restrict attention to the class of invariant tests, but we also consider the case when the autocorrelation is due to the presence of a spatially lagged dependent variable among the regressors. A numerical study aimed at assessing the practical relevance of the theoretical results is included

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

Monte Carlo phase diagram for a polydisperse diblock copolymer melt

The phase diagram for an AB diblock copolymer melt with polydisperse A blocks and monodisperse B blocks is evaluated using lattice-based Monte Carlo simulations. Experiments on this system have shown that the A-block polydispersity shifts the order-order transitions (OOTs) towards higher A-monomer content, while the order-disorder transition (ODT) moves towards higher temperatures when the A blocks form the minority domains and lower temperatures when the A blocks form the matrix. Although self-consistent field theory (SCFT) correctly accounts for the change in the OOTs, it incorrectly predicts the ODT to shift towards higher temperatures at all diblock copolymer compositions. In contrast, our simulations predict the correct shifts for both the OOTs and the ODT. This implies that polydispersity amplifies the fluctuation-induced correction to the mean-field ODT, which we attribute to a reduction in packing frustration. Consistent with this explanation, polydispersity is found to enhance the stability of the perforated-lamellar phase.

High-strength, healable, supramolecular polymer nanocomposites

A supramolecular polymer blend, formed via π-π interactions between a π-electron rich pyrenyl endcapped oligomer and a chain-folding oligomer containing pairs of π-electron poor naphthalene-diimide (NDI) units, has been reinforced with cellulose nanocrystals (CNCs) to afford a healable nanocomposite material. Nanocomposites with varying weight percentage of CNCs (from 1.25 to 20.0 wt.%) within the healable supramolecular polymeric matrix have been prepared via solvent casting followed by compression molding, and their mechanical properties and healing behavior have been evaluated. It is found that homogeneously dispersed films can be formed with CNCs at less than 10 wt.%. Above 10 wt.% CNC heterogeneous nanocomposites were obtained. All the nanocomposites formed could be re-healed upon exposure to elevated temperatures although, for the homogeneous films, it was found that the healing rate was reduced with increasing CNC content. The best combination of healing efficiency and mechanical properties was obtained with the 7.5 wt.% CNC nanocomposite which exhibited a tensile modulus enhanced by as much as a factor of 20 over the matrix material alone and could be fully re-healed at 85 °C within 30 minutes. Thus it is demonstrated that supramolecular nanocomposites can afford greatly enhanced mechanical properties relative to the unreinforced polymer, while still allowing efficient thermal healing.

A Lab-in-a-briefcase for rapid prostate specific antigen (PSA) screening from whole blood

We present a new concept for rapid and fully portable Prostate Specific Antigen (PSA) measurement, termed “Lab-in-a-Briefcase”, which integrates an affordable microfluidic ELISA platform utilising a melt-extruded fluoropolymer Micro Capillary Film (MCF) containing 10 bore, 200 μm internal diameter capillaries, a disposable multi-syringe aspirator (MSA) plus a sample tray pre-loaded with all required immunoassay reagents, and a portable film scanner for colorimetric signal digital quantitation. Each MSA can perform 10 replicate microfluidic immunoassays on 8 samples, allowing 80measurements to be made in less than 15 minutes based on semi-automated operation and norequirement of additional fluid handling equipment. An assay was optimised for measurement of a clinically relevant range of PSA from 0.9 to 60.0 ng/ml in 15 minutes with CVs in the order of 5% based on intra-assay variability when read using a consumer flatbed film scanner. The PSA assay performance in the MSA remained robust in the presence of undiluted or 1:2 diluted human serum or whole blood, and the matrix effect could simply be overcome by extending sample incubation times. The PSA "Lab-in-a-briefcase" is particularly suited to a low-resource health setting where diagnostic labs and automated immunoassay systems are not accessible, by allowing PSA measurement outside the laboratory using affordable equipment.

Efficient animal-serum free 3D cultivation method for adult human neural crest-derived stem cell therapeutics

Due to their broad differentiation potential and their persistence into adulthood, human neural crest-derived stem cells (NCSCs) harbour great potential for autologous cellular therapies, which include the treatment of neurodegenerative diseases and replacement of complex tissues containing various cell types, as in the case of musculoskeletal injuries. The use of serum-free approaches often results in insufficient proliferation of stem cells and foetal calf serum implicates the use of xenogenic medium components. Thus, there is much need for alternative cultivation strategies. In this study we describe for the first time a novel, human blood plasma based semi-solid medium for cultivation of human NCSCs. We cultivated human neural crest-derived inferior turbinate stem cells (ITSCs) within a blood plasma matrix, where they revealed higher proliferation rates compared to a standard serum-free approach. Three-dimensionality of the matrix was investigated using helium ion microscopy. ITSCs grew within the matrix as revealed by laser scanning microscopy. Genetic stability and maintenance of stemness characteristics were assured in 3D cultivated ITSCs, as demonstrated by unchanged expression profile and the capability for self-renewal. ITSCs pre-cultivated in the 3D matrix differentiated efficiently into ectodermal and mesodermal cell types, particularly including osteogenic cell types. Furthermore, ITSCs cultivated as described here could be easily infected with lentiviruses directly in substrate for potential tracing or gene therapeutic approaches. Taken together, the use of human blood plasma as an additive for a completely defined medium points towards a personalisable and autologous cultivation of human neural crest-derived stem cells under clinical grade conditions.

Effect of fat replacement by inulin on textural and structural properties of short dough biscuits

The aim of this study was to evaluate the effects of inulin as fat replacer on short dough biscuits and their corresponding doughs. A control formulation, with no replacement, and four formulations in which 10, 20, 30, and 40 % of shortening was replaced by inulin were studied. In the dough, shortening was observed surrounding flour components. At higher fat replacement levels, flour was more available for hydration leading to significant (P<0.05) harder doughs: from 2.76 (0.12)N in 10 % fat-replaced biscuits to 5.81 (1.56)N in 30 % fat-replaced ones. Biscuit structure was more continuous than dough structure. A continuous fat layer coated the matrix surface, where starch granules were embedded. In general, weight loss during baking and water activity decreased significantly (P<0.05) as fat replacement increased. Biscuit dimensions and aeration decreased when fat replacement increased, e.g., width gain was +1.20 mm in 10 fat-replaced biscuits and only +0.32 mm in 40 % fat-replaced ones. Panelist found biscuits with 20 % of fat replacement slightly harder than control biscuits. It can be concluded that shortening may be partially replaced, up to 20 %, with inulin. These low fat biscuits are similar than the control biscuits, and they can have additional health benefits derived from inulin presence.

Functionality of several cake ingredients: a comprehensive approach

The roles of some cake ingredients – oil, a leavening agent, and inulin – in the structure and physicochemical properties of batter and cakes were studied in four different formulations. Oil played an important role in the batter stability, due to its contribution to increasing batter viscosity and occluding air during mixing. The addition of the leavening agent was crucial to the final height and sponginess of the cakes. When inulin was used as a fat replacer, the absence of oil caused a decrease in the stability of the batter, where larger air bubbles were occluded. Inulin dispersed uniformly in the batter could create a competition for water with the flour components: gluten was not properly hydrated and some starch granules were not fully incorporated into the matrix. Thus, the development of a continuous network was disrupted and the cake was shorter and softer; it contained interconnected air cells in the crumb, and was easily crumbled. The structure studies were decisive to understand the physicochemical properties.

Commercial thickeners used by patients with dysphagia: rheological and structural behaviour in different food matrices

In order to achieve a safe swallowing in patients with dysphagia, liquids must be thickened. In this work, two commercial starch based thickeners dissolved in water, whole milk, apple juice and tomato juice were studied. The thickeners were Resource®, composed of modified maize starch and Nutilis®, composed of modified maize starch and gums. They were formulated at two different concentrations corresponding to nectar- and pudding-like consistencies. Influence of composition, concentration and food matrix on rheological properties and structure of the resulting pastes were analysed. Viscoelastic measurements and microscopic observations of the thickeners dissolved in water revealed structural differences due to the presence of gums. When the thickeners were dissolved in the other food matrices significant statistical interactions were found between the matrix and the thickener-type in both the viscoelastic and flow parameters. The most relevant differences were observed for the nectar-like consistency with Nutilis® thickener in milk and apple juice. These samples had lower zero viscosity values and higher loss tangent values, that corresponded to weaker structured systems. Light microscopy images showed that the matrix formed by swollen starch granules was interrupted by the presence of gums. The structure of the matrices in pudding-like formulations became more continuous irrespectively of the matrix employed, and also differences in viscoelasticity among samples diminished. Although differences were observed in zero shear viscosity values among samples, the viscosity of the beverages at 50 s−1 – commonly used as a reference for swallowing – was similar for all samples regardless of the matrix used.

Strangely dual orbifold equivalence I

In this brief note we prove orbifold equivalence between two potentials described by strangely dual exceptional unimodular singularities of type K14 and Q10 in two different ways. The matrix factorizations proving the orbifold equivalence give rise to equations whose solutions are permuted by Galois groups which differ for different expressions of the same singularity.

An Alternative Section Method for Casting and Posterior Laser Welding of Metallic Frameworks for an Implant-Supported Prosthesis

Purpose: The aim of this study was to compare the accuracy of fit of three types of implant-supported frameworks cast in Ni-Cr alloy: specifically, a framework cast as one piece compared to frameworks cast separately in sections to the transverse or the diagonal axis, and later laser welded. Materials and Methods: Three sets of similar implant-supported frameworks were constructed. The first group of six 3-unit implant-supported frameworks were cast as one piece, the second group of six were sectioned in the transverse axis of the pontic region prior to casting, and the last group of six were sectioned in the diagonal axis of the pontic region prior to casting. The sectioned frameworks were positioned in the matrix (10 N(.)cm torque) and laser welded. To evaluate passive fit, readings were made with an optical microscope with both screws tightened and with only one-screw tightened. Data were submitted to ANOVA and Tukey-Kramer`s test (p < 0.05). Results: When both screws were tightened, no differences were found between the three groups (p > 0.05). In the single-screw-tightened test, with readings made opposite to the tightened side, the group cast as one piece (57.02 +/- 33.48 mu m) was significantly different (p < 0.05) from the group sectioned diagonally (18.92 +/- 4.75 mu m) but no different (p > 0.05) from the group transversally sectioned (31.42 +/- 20.68 mu m). On the tightened side, no significant differences were found between the groups (p > 0.05). Conclusions: Results of this study showed that casting diagonally sectioned frameworks lowers misfit levels of prosthetic implant-supported frameworks and also improves the levels of passivity to the same frameworks when compared to structures cast as one piece.

Wear of Ball Attachments After 1 to 8 Years of Clinical Use: A Qualitative Analysis

The purpose of this study was to analyze and compare scanning electron microscopy (SEM) observations of ball attachments that had been worn by patients during three periods of clinical use. One hundred forty-four specimens of ball anchor attachments (gold alloy matrix and titanium patrix) were studied by SEM after periods of approximately 1, 3.5, and 8 years of clinical use. Twenty new attachment components were examined as controls. SEM images revealed signs of mechanical wear for the ball attachments studied. The surfaces of the titanium patrix were associated primarily with roughening after short-term use, whereas surfaces of the gold alloy matrix showed wear, roughening, and loss of microscopic material in the form of flakes. Severe mechanical wear on both surfaces was noted after longer periods of use. The mechanical changes were not correlated with patient-mediated observations regarding the time-dependent retentive efficacy of the attachments. One year of clinical wear appeared to have limited effect on the ball attachment tested. Conversely, longer periods of use led to marked modifications in shape of the matrix and patrix components. Int J Prosthodont 2011;24:270-272.