Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the system's symmetries.

Full S matrix calculation via a single real-symmetric Lanczos recursion: The Lanczos artificial boundary inhomogeneity method

We present an efficient and robust method for the calculation of all S matrix elements (elastic, inelastic, and reactive) over an arbitrary energy range from a single real-symmetric Lanczos recursion. Our new method transforms the fundamental equations associated with Light's artificial boundary inhomogeneity approach [J. Chem. Phys. 102, 3262 (1995)] from the primary representation (original grid or basis representation of the Hamiltonian or its function) into a single tridiagonal Lanczos representation, thereby affording an iterative version of the original algorithm with greatly superior scaling properties. The method has important advantages over existing iterative quantum dynamical scattering methods: (a) the numerically intensive matrix propagation proceeds with real symmetric algebra, which is inherently more stable than its complex symmetric counterpart; (b) no complex absorbing potential or real damping operator is required, saving much of the exterior grid space which is commonly needed to support these operators and also removing the associated parameter dependence. Test calculations are presented for the collinear H+H-2 reaction, revealing excellent performance characteristics. (C) 2004 American Institute of Physics.

On pairwise comparison matrices that can be made consistent by the modification of a few elements

Pairwise comparison matrices are often used in Multi-attribute Decision Making forweighting the attributes or for the evaluation of the alternatives with respect to a criteria. Matrices provided by the decision makers are rarely consistent and it is important to index the degree of inconsistency. In the paper, the minimal number of matrix elements by the modification of which the pairwise comparison matrix can be made consistent is examined. From practical point of view, the modification of 1, 2, or, for larger matrices, 3 elements seems to be relevant. These cases are characterized by using the graph representation of the matrices. Empirical examples illustrate that pairwise comparison matrices that can be made consistent by the modification of a few elements are present in the applications.

Spin alignment measurements of the K*(0)(892) and phi(1020) vector mesons in heavy ion collisions at root S(NN)=200 GeV

We present the first spin alignment measurements for the K*(0)(892) and phi(1020) vector mesons produced at midrapidity with transverse momenta up to 5 GeV/c at root s(NN) = 200 GeV at RHIC. The diagonal spin-density matrix elements with respect to the reaction plane in Au+Au collisions are rho(00) = 0.32 +/- 0.04 (stat) +/- 0.09 (syst) for the K*(0) (0.8 < p(T) < 5.0 GeV/c) and rho(00) = 0.34 +/- 0.02 (stat) +/- 0.03 (syst) for the phi (0.4 < p(T) < 5.0 GeV/c) and are constant with transverse momentum and collision centrality. The data are consistent with the unpolarized expectation of 1/3 and thus no evidence is found for the transfer of the orbital angular momentum of the colliding system to the vector-meson spins. Spin alignments for K(*0) and phi in Au+Au collisions were also measured with respect to the particle's production plane. The phi result, rho(00) = 0.41 +/- 0.02 (stat) +/- 0.04 (syst), is consistent with that in p+p collisions, rho(00) = 0.39 +/- 0.03 (stat) +/- 0.06 (syst), also measured in this work. The measurements thus constrain the possible size of polarization phenomena in the production dynamics of vector mesons.

rho(0) photoproduction in ultraperipheral relativistic heavy ion collisions at root s(NN)=200 GeV

Photoproduction reactions occur when the electromagnetic field of a relativistic heavy ion interacts with another heavy ion. The STAR Collaboration presents a measurement of rho(0) and direct pi(+)pi(-) photoproduction in ultraperipheral relativistic heavy ion collisions at root s(NN) = 200 GeV. We observe both exclusive photoproduction and photoproduction accompanied by mutual Coulomb excitation. We find a coherent cross section of sigma(AuAu -> Au*Au*rho(0)) = 530 +/- 19(stat.) +/- 57(syst.) mb, in accord with theoretical calculations based on a Glauber approach, but considerably below the predictions of a color dipole model. The rho 0 transverse momentum spectrum (p(T)(2)) is fit by a double exponential curve including both coherent and incoherent coupling to the target nucleus; we find sigma(inc)/sigma(coh) = 0.29 +/- 0.03 (stat.) +/- 0.08 (syst.). The ratio of direct pi(+)pi(-) to rho(0) production is comparable to that observed in gamma(p) collisions at HERA and appears to be independent of photon energy. Finally, the measured rho(0) spin helicity matrix elements agree within errors with the expected s-channel helicity conservation.

Correlation between MMPs and their inhibitors in breast cancer tumor tissue specimens and in cell lines with different metastatic potential

Background: The metastatic disease rather than the primary tumor itself is responsible for death in most solid tumors, including breast cancer. The role of matrix metalloproteinases ( MMPs), tissue inhibitors of MMPs (TIMPs) and Reversion-inducing cysteine-rich protein with Kazal motifs ( RECK) in the metastatic process has previously been established. However, in all published studies only a limited number of MMPs/MMP inhibitors was analyzed in a limited number of cell lines. Here, we propose a more comprehensive approach by analyzing the expression levels of several MMPs (MMP-2, MMP-9 and MMP-14) and MMP inhibitors (TIMP-1, TIMP-2 and RECK) in different models ( five human breast cancer cell lines, 72 primary breast tumors and 30 adjacent normal tissues). Methods: We analyzed the expression levels of MMP-2, MMP-9 and MMP-14 and their inhibitors (TIMP-1, TIMP-2 and RECK) by quantitative RT-PCR (qRT-PCR) in five human breast cancer cell lines presenting increased invasiveness and metastatic potential, 72 primary breast tumors and 30 adjacent normal tissues. Moreover, the role of cell-extracellular matrix elements interactions in the regulation of expression and activity of MMPs and their inhibitors was analyzed by culturing these cell lines on plastic or on artificial ECM (Matrigel). Results: The results demonstrated that MMPs mRNA expression levels displayed a positive and statistically significant correlation with the transcriptional expression levels of their inhibitors both in the cell line models and in the tumor tissue samples. Furthermore, the expression of all MMP inhibitors was modulated by cell-Matrigel contact only in highly invasive and metastatic cell lines. The enzyme/inhibitor balance at the transcriptional level significantly favors the enzyme which is more evident in tumor than in adjacent non-tumor tissue samples. Conclusion: Our results suggest that the expression of MMPs and their inhibitors, at least at the transcriptional level, might be regulated by common factors and signaling pathways. Therefore, the multi-factorial analysis of these molecules could provide new and independent prognostic information contributing to the determination of more adequate therapy strategies for each patient.`

TL and OSL properties of KAlSi(3)O(8):Mn, obtained by sol-gel process

Thermoluminescence (TL) and Optically Stimulated Luminescence (OSL) properties of KAlSi(3)O(8):Mn glasses obtained through the sol gel technique were investigated. Samples were obtained with five different molar concentrations of 0.25, 0.5, 1, 2 and 5 mol% of manganese. Transmission Electronic Microscopy (TEM) indicated the occurrence of nanoparticles composed by glass matrix elements with Mn. Best results for TL response were obtained with 0.5 mol% Mn doped sample, which exhibits a TL peak at 180 degrees C. The TL spectrum of this sample presents a broad emission band from 450 to 700 nm with a peak at 575 nm approximately. The emission band fits very well with the characteristic lines of the Mn(2+) emission features. According to this fact, the band at 410 nm can be ascribed to (6)A(1)(S) -> (4)A(1)(G), (4)E(G) transition, while the 545 nm band can be attributed to the superposition of the transitions (6)A(1)(S) -> (4)T(2)(G) and (6)A(1)(S) -> (4)T(1)(G). The dependence of the TL response with the energy of X-rays (27-41 keV) showed a small decrease of the TL intensity in the high energy region. Excitation with blue LEDs showed OSL in the UV region with a fast decay component. (C) 2011 Elsevier Ltd. All rights reserved.

Algebraic Bethe ansatz method for the exact calculation of energy spectra and form factors: applications to models of Bose-Einstein condensates and metallic nanograins

In this review we demonstrate how the algebraic Bethe ansatz is used for the calculation of the-energy spectra and form factors (operator matrix elements in the basis of Hamiltonian eigenstates) in exactly solvable quantum systems. As examples we apply the theory to several models of current interest in the study of Bose-Einstein condensates, which have been successfully created using ultracold dilute atomic gases. The first model we introduce describes Josephson tunnelling between two coupled Bose-Einstein condensates. It can be used not only for the study of tunnelling between condensates of atomic gases, but for solid state Josephson junctions and coupled Cooper pair boxes. The theory is also applicable to models of atomic-molecular Bose-Einstein condensates, with two examples given and analysed. Additionally, these same two models are relevant to studies in quantum optics; Finally, we discuss the model of Bardeen, Cooper and Schrieffer in this framework, which is appropriate for systems of ultracold fermionic atomic gases, as well as being applicable for the description of superconducting correlations in metallic grains with nanoscale dimensions.; In applying all the above models to. physical situations, the need for an exact analysis of small-scale systems is established due to large quantum fluctuations which render mean-field approaches inaccurate.

Casimir invariants and characteristic identities for gl(infinity)

A full set of (higher-order) Casimir invariants for the Lie algebra gl(infinity) is constructed and shown to be well defined in the category O-FS generated by the highest weight (unitarizable) irreducible representations with only a finite number of nonzero weight components. Moreover, the eigenvalues of these Casimir invariants are determined explicitly in terms of the highest weight. Characteristic identities satisfied by certain (infinite) matrices with entries from gl(infinity) are also determined and generalize those previously obtained for gl(n) by Bracken and Green [A. J. Bracken and H. S. Green, J. Math. Phys. 12, 2099 (1971); H. S. Green, ibid. 12, 2106 (1971)]. (C) 1997 American Institute of Physics.

Combined chemical separation of Lu, Hf, Sm, Nd, and REEs from a single rock digest: Precise and accurate isotope Determinations of Lu-Hf and Sm-Nd using multicollector-ICPMS

A combined procedure for separating Lu, Hf, Sm, Nd, and rare earth elements (REEs) from a single sample digest is presented. The procedure consists of the following five steps: (1) sample dissolution via sodium peroxide sintering; (2) separation of the high field strength elements from the REEs and other matrix elements by a HF-free anion-exchange column procedure; (3) purification of Hf on a cation-exchange resin; (4) separation of REEs from other matrix elements by cation exchange; (5) Lu, Sm, and Nd separation from the other REEs by reversed-phase ion chromatography. Analytical reproducibilities of Sm-Nd and Lu-Hf isotope systematics are demonstrated for standard solutions and international rock reference materials. Results show overall good reproducibilities for Sm-Nd systematics independent of the rock type analyzed. For the Lu-Hf systematics, the reproducibility of the parent/daughter ratio is much better for JB-1 (basalt) than for two analyzed felsic crustal rocks (DR-N and an Archaean granitoid). It is demonstrated that this poorer reproducibility of the Lu/Hf ratio is truly caused by sample heterogeneity; thus, results are geologically reasonable.

New textures for the lepton mass matrices

We study predictive textures for the lepton mass matrices in which the charged-lepton mass matrix has either four or five zero matrix elements while the neutrino Majorana mass matrix has, respectively, either four or three zero matrix elements. We find that all the viable textures of these two kinds share many predictions: the neutrino mass spectrum is inverted, the sum of the light-neutrino masses is close to 0.1 eV, the Dirac phase delta in the lepton mixing matrix is close to either 0 or pi, and the mass term responsible for neutrinoless double-beta decay lies in between 12 and 22 meV. (C) 2014 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license.

Measurement of the inclusive jet cross-section in proton--proton collisions at s√=7 TeV using 4.5 fb−1 of data with the ATLAS detector

The inclusive jet cross-section is measured in proton--proton collisions at a centre-of-mass energy of 7 TeV using a data set corresponding to an integrated luminosity of 4.5 fb−1 collected with the ATLAS detector at the Large Hadron Collider in 2011. Jets are identified using the anti-kt algorithm with radius parameter values of 0.4 and 0.6. The double-differential cross-sections are presented as a function of the jet transverse momentum and the jet rapidity, covering jet transverse momenta from 100 GeV to 2 TeV. Next-to-leading-order QCD calculations corrected for non-perturbative effects and electroweak effects, as well as Monte Carlo simulations with next-to-leading-order matrix elements interfaced to parton showering, are compared to the measured cross-sections. A quantitative comparison of the measured cross-sections to the QCD calculations using several sets of parton distribution functions is performed.

Straightening out the concept of direct and indirect input requirements

Gim & Kim (1998) proposed a generalization of Jeong (1982, 1984) reinterpretation of the Hawkins-Simon condition for macroeconomic stability to off-diagonal matrix elements. This generalization is conceptually relevant for it offers a complementary view of interindustry linkages beyond final or net output influence. The extension is completely similar to the 'total flow' idea introduced by Szyrmer (1992) or the 'output-to-output' multiplier of Miller & Blair (2009). However the practical implementation of Gim & Kim is actually faulty since it confuses the appropriate order of output normalization. We provide a new and elementary solution for the correct formalization using standard interindustry accounting concepts.

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA