Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Impact of plastic deformation on magnetoacoustic properties of Fe-2%Si alloy

The paper presents the results of a complementary study including magnetic hysteresis loops B(H), magnetic Barkhausen noise (MBN) and magnetoacoustic emission (MAE) signals measurements for plastically deformed Fe-2%Si samples. The investigated samples had been plastically deformed with plastic strain level (epsilon(p)) up to 8%. The properties of B(H) loops are quantified using the coercivity H(C) and maximum differential permeability mu(rmax) as parameters. The MBN and MAE voltage signals were analysed by means of rms-like voltage (Ub and Ua, respectively) envelopes, plotted as a function of applied field strength. Integrals of the Ub and Ua voltages over half of a period of magnetization were then calculated. It has been found that He and integrals of Ub increase, while mu(rmax) decreases monotonically with increasing epsilon(p). The MAE (Ua) peak voltage at first decreases, then peaks at epsilon(p) approximate to 1.5% and finally decreases again. The integral of the Ua voltage at first increases for low epsilon(p) and then decreases for epsilon(p) > 1.5%. All those various dependence types suggest the possibility of detection of various stages of microstructure change. The above-mentioned results are discussed qualitatively in the paper. Some modelling of the discussed dependency is also presented. (C) 2008 Elsevier Ltd. All rights reserved.

Light scattering in polycrystalline alumina with bi-dimensionally large surface grains

Alumina ceramics with high in-line transmittance at 0.5-1.0 mm-thickness were prepared with different doping additives by sintering at 1850 degrees C in vacuum for 1-8 h. Depending on the additive contents and sintering variables bi-dimensionally large surface grains, caused by surface evaporation of MgO, had grown parallel to the surface with similar to 100 mu m thickness and lateral sizes up to the millimeter range. The abnormal grain-growth process also resulted in the formation of pores entrapped inside the large surface grains within a narrow zone at 10-20 mu m distance from the surface. The fraction of these pores is thickness-invariant. Scattering factors associated to the pores entrapped inside the bi-dimensionally large surface grains, second-phase particles, grain-boundaries, and microstructural surface defects are derived from the results of in-line transmission (at 600 nm) and are used together with microstructural characteristics to explain the light transmittance in these materials. (C) 2008 Elsevier Ltd. All rights reserved.

Rotation-discriminating template matching based on Fourier coefficients of radial projections with robustness to scaling and partial occlusion

We consider brightness/contrast-invariant and rotation-discriminating template matching that searches an image to analyze A for a query image Q We propose to use the complex coefficients of the discrete Fourier transform of the radial projections to compute new rotation-invariant local features. These coefficients can be efficiently obtained via FFT. We classify templates in ""stable"" and ""unstable"" ones and argue that any local feature-based template matching may fail to find unstable templates. We extract several stable sub-templates of Q and find them in A by comparing the features. The matchings of the sub-templates are combined using the Hough transform. As the features of A are computed only once, the algorithm can find quickly many different sub-templates in A, and it is Suitable for finding many query images in A, multi-scale searching and partial occlusion-robust template matching. (C) 2009 Elsevier Ltd. All rights reserved.

A Separation Principle for the Continuous-Time LQ-Problem With Markovian Jump Parameters

In this paper, we devise a separation principle for the finite horizon quadratic optimal control problem of continuous-time Markovian jump linear systems driven by a Wiener process and with partial observations. We assume that the output variable and the jump parameters are available to the controller. It is desired to design a dynamic Markovian jump controller such that the closed loop system minimizes the quadratic functional cost of the system over a finite horizon period of time. As in the case with no jumps, we show that an optimal controller can be obtained from two coupled Riccati differential equations, one associated to the optimal control problem when the state variable is available, and the other one associated to the optimal filtering problem. This is a separation principle for the finite horizon quadratic optimal control problem for continuous-time Markovian jump linear systems. For the case in which the matrices are all time-invariant we analyze the asymptotic behavior of the solution of the derived interconnected Riccati differential equations to the solution of the associated set of coupled algebraic Riccati equations as well as the mean square stabilizing property of this limiting solution. When there is only one mode of operation our results coincide with the traditional ones for the LQG control of continuous-time linear systems.

Flatness-based control of a single qubit gate

This work considers the open-loop control problem of steering a two-level quantum system from any initial to any final condition. The model of this system evolves on the state space X = SU(2), having two inputs that correspond to the complex amplitude of a resonant laser field. A symmetry preserving flat output is constructed using a fully geometric construction and quaternion computations. Simulation results of this flatness-based open-loop control are provided.

Non-linear modal analysis for beams subjected to axial loads: Analytical and finite-element solutions

A rigorous derivation of non-linear equations governing the dynamics of an axially loaded beam is given with a clear focus to develop robust low-dimensional models. Two important loading scenarios were considered, where a structure is subjected to a uniformly distributed axial and a thrust force. These loads are to mimic the main forces acting on an offshore riser, for which an analytical methodology has been developed and applied. In particular, non-linear normal modes (NNMs) and non-linear multi-modes (NMMs) have been constructed by using the method of multiple scales. This is to effectively analyse the transversal vibration responses by monitoring the modal responses and mode interactions. The developed analytical models have been crosschecked against the results from FEM simulation. The FEM model having 26 elements and 77 degrees-of-freedom gave similar results as the low-dimensional (one degree-of-freedom) non-linear oscillator, which was developed by constructing a so-called invariant manifold. The comparisons of the dynamical responses were made in terms of time histories, phase portraits and mode shapes. (C) 2008 Elsevier Ltd. All rights reserved.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

One-dimensional point interaction with three complex parameters

For a pair of non-Hermitian Hamiltonian H and its Hermitian adjoint H(dagger), there are situations in which their eigenfunctions form a biorthogonal system. We illustrate such a situation by means of a one-particle system with a one-dimensional point interaction in the form of the Fermi pseudo-potential. The interaction consists of three terms with three strength parameters g(i) (i = 1, 2 and 3), which are all complex. This complex point interaction is neither Hermitian nor PT-invariant in general. The S-matrix for the transmission reflection problem constructed with H (or with H(dagger)) in the usual manner is not unitary, but it conforms to the pseudo-unitarity that we define. The pseudounitarity is closely related to the biorthogonality of the eigenfunctions. The eigenvalue spectrum of H with the complex interaction is generally complex but there are cases where the spectrum is real. In such a case H and H(dagger) form a pseudo-Hermitian pair.

Hypomania: a transcultural perspective

This study examined the transcultural robustness of a screening instrument for hypomania, the Hypomania Checklist-32. first revised version(HCL-32 R1) It was carried out in 2606 patients from twelve countries in five geographic regions (Northern, Southern and Eastern Europe, South America and East Asia) In addition, GAMIAN Europe contributed data from its members Exploratory and confirmatory factor analyses were used to examine the transregional stability of the measurement properties of the HCL-32 R 1, including the influence of sex and age as covariates A cross cultures, a two-factor structure was confirmed the first factor (F 1) reflected the more positive aspects of hypomania (being more active, elated, self-confident, and cognitively, enhanced), the second factor (F2) reflected the more negative aspects (being irritable, impulsive, careless. more substance use). The measurement properties of the HCL-32 R1 were largely invariant across cultures Only few items showed transcultural differences in their relation to hypomania as measured by the test F2 was higher among men and in more severe manic syndromes, F1 was highest in North and East Europe and lowest in South America Vie scores decreased slightly with age The frequency of the 32 items showed remarkable similarities across geographic areas, with two exceptions. South Europeans had lower symptom frequencies in general and East Europeans higher rates of substance use These findings support the international applicability of the HCL-32 R1 as a screening instrument for hypomania

Lower cytokine secretion ex vivo by natural killer T cells in HIV-infected individuals is associated with higher CD161 expression

Objective: Natural killer T (NKT) cells are efficiently targeted by HIV and severely reduced in numbers in the circulation of infected individuals. The functional capacity of the remaining NKT cells in HIV-infected individuals is poorly characterized. This study measured NKT cell cytokine production directly ex vivo and compared these responses with both the disease status and NKT subset distribution of individual patients. Methods: NKT cell frequencies, subsets, and ex-vivo effector functions were measured in the peripheral blood mononuclear cells of HIV-infected patients and healthy controls by flow cytometry. We measured cytokines from NKT cells after stimulation with either a-galactosyl ceramide-loaded CD1d dimers (DimerX-alpha GalCer) or phorbol myristate acetate and ionomycin. Results: The frequencies of NKT cells secreting interferon-gamma and tumor necrosis factor-alpha were significantly lower in HIV-infected patients than healthy controls after DimerX-alpha GalCer treatment, but responses were similar after treatment with phorbol myristate acetate and ionomycin. The magnitude of the interferon-gamma response to DimerX-alpha GalCer correlated inversely with the number of years of infection. Both interferon-gamma and tumor necrosis factor-alpha production in response to DimerX-alpha GalCer correlated inversely with CD161 expression. Conclusion: The ex-vivo Th1 responses of circulating NKT cells to CD1d-glycolipid complexes are impaired in HIV-infected patients. NKT cell functions may be progressively lost over time in HIV infection, and CD161 is implicated in the regulation of NKT cell responsiveness. (C) 2009 Wolters Kluwer Health vertical bar Lippincott Williams & Wilkins

Study of the discovery potential of light custodians in a clean decay channel at the LHC

Models of warped extra dimensions with custodial symmetry usually predict the existence of a light Kaluza-Klein fermion arising as a partner of the right-handed top quark, sometimes called light custodians which we will denote (b) over tilde (R). The production of these particles at the LHC can give rise to multi-W events which could be observed in same-sign dilepton channels, but its mass reconstruction is challenging. In this paper we study the possibility of finding a signal for the pair production of this new particle at the LHC focusing on a rarer, but cleaner decay mode of a light custodian into a Z boson and a b-quark. In this mode it would be possible to reconstruct the light custodian mass. In addition to the dominant standard model QCD production processes, we include the contribution of a Kaluza-Klein gluon first mode. We find that (b) over tilde (R) stands out from the background as a peak in the bZ invariant mass. However, when taking into account only the electronic and muonic decay modes of the Z boson and b-tagging efficiencies, the LHC will have access only to the very light range of masses, m((b) over tilde) = O(500) GeV.

Quantized fields and gravitational particle creation in f (R) expanding universes

The problem of cosmological particle creation for a spatially flat, homogeneous and isotropic universes is discussed in the context of f (R) theories of gravity. Different from cosmological models based on general relativity theory, it is found that a conformal invariant metric does not forbid the creation of massless particles during the early stages (radiation era) of the universe. (C) 2010 Elsevier B.V. All rights reserved.

Dynamic portrait of the planetary 2/1 mean-motion resonance - II. Systems with a more massive inner planet

This paper presents the second part in our study of the global structure of the planar phase space of the planetary three-body problem, when both planets lie in the vicinity of a 2/1 mean-motion resonance. While Paper I was devoted to cases where the outer planet is the more massive body, the present work is devoted to the cases where the more massive body is the inner planet. As before, outside the well-known Apsidal Corotation Resonances (ACR), the phase space shows a complex picture marked by the presence of several distinct regimes of resonant and non-resonant motion, crossed by families of periodic orbits and separated by chaotic zones. When the chosen values of the integrals of motion lead to symmetric ACR, the global dynamics are generally similar to the structure presented in Paper I. However, for asymmetric ACR the resonant phase space is strikingly different and shows a galore of distinct dynamical states. This structure is shown with the help of dynamical maps constructed on two different representative planes, one centred on the unstable symmetric ACR and the other on the stable asymmetric equilibrium solution. Although the study described in the work may be applied to any mass ratio, we present a detailed analysis for mass values similar to the Jupiter-Saturn case. Results give a global view of the different dynamical states available to resonant planets with these characteristics. Some of these dynamical paths could have marked the evolution of the giant planets of our Solar system, assuming they suffered a temporary capture in the 2/1 resonance during the latest stages of the formation of our Solar system.

The Immunomodulatory Action of Sialostatin L on Dendritic Cells Reveals Its Potential to Interfere with Autoimmunity

Sialostatin L (SialoL) is a secreted cysteine protease inhibitor identified in the salivary glands of the Lyme disease vector Ixodes scapularis. In this study, we reveal the mechanisms of SialoL immunomodulatory actions on the vertebrate host. LPS-induced maturation of dendritic cells from C57BL/6 mice was significantly reduced in the presence of SialoL. Although OVA degradation was not affected by the presence of SialoL in dendritic cell cultures, cathepsin S activity was partially inhibited, leading to an accumulation of a 10-kDa invariant chain intermediate in these cells. As a consequence, in vitro Ag-specific CD4(+) T cell proliferation was inhibited in a time-dependent manner by SialoL, and further studies engaging cathepsin S(-/-) or cathepsin L(-/-) dendritic cells confirmed that the immunomodulatory actions of SialoL are mediated by inhibition of cathepsin S. Moreover, mice treated with SialoL displayed decreased early T cell expansion and recall response upon antigenic stimulation. Finally, SialoL administration during the immunization phase of experimental autoimmune encephalomyelitis in mice significantly prevented disease symptoms, which was associated with impaired IFN-gamma and IL-17 production and specific T cell proliferation. These results illuminate the dual mechanism by which a human disease vector protein modulates vertebrate host immunity and reveals its potential in prevention of an autoimmune disease. The Journal of Immunology, 2009, 182: 7422-7429.