Fluctuations of Matrix Elements of Regular Functions of Gaussian Random Matrices
NASA Astrophysics Data System (ADS)
Lytova, A.; Pastur, L.
2009-01-01
We find the limit of the variance and prove the Central Limit Theorem (CLT) for the matrix elements φ jk ( M), j, k=1,…, n of a regular function φ of the Gaussian matrix M (GOE and GUE) as its size n tends to infinity. We show that unlike the linear eigenvalue statistics Tr φ( M), a traditional object of random matrix theory, whose variance is bounded as n→∞ and the CLT is valid for Tr φ( M)- E{Tr φ( M)}, the variance of φ jk ( M) is O(1/ n), and the CLT is valid for sqrt{n}(\\varphi _{jk}(M)-E\\{\\varphi _{jk}(M)\\}) . This shows the role of eigenvectors in the forming of the asymptotic regime of various functions (statistics) of random matrices. Our proof is based on the use of the Fourier transform as a basic characteristic function, unlike the Stieltjes transform and moments, used in majority of works of the field. We also comment on the validity of analogous results for other random matrices.
Extracellular matrix fluctuations during early embryogenesis
Szabó, A; Rupp, P A; Rongish, B J; Little, C D; Czirók, A
2011-01-01
Extracellular matrix (ECM) movements and rearrangements were studied in avian embryos during early stages of development. We show that the ECM moves as a composite material, whereby distinct molecular components as well as spatially separated layers exhibit similar displacements. Using scanning wide field and confocal microscopy we show that the velocity field of ECM displacement is smooth in space and that ECM movements are correlated even at locations separated by several hundred micrometers. Velocity vectors, however, strongly fluctuate in time. The autocorrelation time of the velocity fluctuations is less than a minute. Suppression of the fluctuations yields a persistent movement pattern that is shared among embryos at equivalent stages of development. The high resolution of the velocity fields allows a detailed spatio-temporal characterization of important morphogenetic processes, especially tissue dynamics surrounding the embryonic organizer (Hensen’s node). PMID:21750366
Matrix Elements for Hylleraas CI
NASA Astrophysics Data System (ADS)
Harris, Frank E.
The limitation to at most a single interelectron distance in individual configurations of a Hylleraas-type multiconfiguration wave function restricts significantly the types of integrals occurring in matrix elements for energy calculations, but even then if the formulation is not handled efficiently the angular parts of these integrals escalate to create expressions of great complexity. This presentation reviews ways in which the angular-momentum calculus can be employed to systematize and simplify the matrix element formulas, particularly those for the kinetic-energy matrix elements.
Relativistic Dipole Matrix Element Zeros
NASA Astrophysics Data System (ADS)
Lajohn, L. A.; Pratt, R. H.
2002-05-01
There is a special class of relativistic high energy dipole matrix element zeros (RZ), whose positions with respect to photon energy ω , only depend on the bound state l quantum number according to ω^0=mc^2/(l_b+1) (independent of primary quantum number n, nuclear charge Z, central potential V and dipole retardation). These RZ only occur in (n,l_b,j_b)arrow (ɛ , l_b+1,j_b) transitions such as ns_1/2arrow ɛ p_1/2; np_3/2arrow ɛ d_3/2: nd_5/2arrow ɛ f_5/2 etc. The nonrelativistic limit of these matrix elements can be established explicitly in the Coulomb case. Within the general matrix element formalism (such as that in [1]); when |κ | is substituted for γ in analytic expressions for matrix elements, the zeros remain, but ω^0 now becomes dependent on n and Z. When the reduction to nonrelativistic form is completed by application of the low energy approximation ω mc^2 mc^2, the zeros disappear. This nonzero behavior was noted in nonrelativistic dipole Coulomb matrix elements by Fano and Cooper [2] and later proven by Oh and Pratt[3]. (J. H. Scofield, Phys. Rev. A 40), 3054 (1989 (U. Fano and J. W. Cooper, Rev. Mod. Phys. 40), 441 (1968). (D. Oh and R. H. Pratt, Phys. Rev. A 34), 2486 (1986); 37, 1524 (1988); 45, 1583 (1992).
Baykara, N. A.; Guervit, Ercan; Demiralp, Metin
2012-12-10
In this work a study on finite dimensional matrix approximations to products of quantum mechanical operators is conducted. It is emphasized that the matrix representation of the product of two operators is equal to the product of the matrix representation of each of the operators when all the fluctuation terms are ignored. The calculation of the elements of the matrices corresponding to the matrix representation of various operators, based on three terms recursive relation is defined. Finally it is shown that the approximation quality depends on the choice of higher values of n, namely the dimension of Hilbert space.
Weak matrix elements for CP violation.
Lee, W.; Gupta, R.; Christ, N.; Fleming, G. T.; Kilcup, G.; Liu, G.; Mawhinney, R.; Sharpe, S.; Wu, L.; Bhattacharya, T.
2001-01-01
We present preliminary results of matrix elements of four fermion operators relevant to the determination of e and E ' / E using staggered fermions. To calculate the matrix elements relevant to CP violation in Kaon decays it is important to use a lattice formulation which preserves (some) chiral symmetry.
Graphical evaluation of relativistic matrix elements
NASA Technical Reports Server (NTRS)
Huang, K. N.
1978-01-01
A graphical representation of angular momentum was used to evaluate relativistic matrix elements between antisymmetrized states of many particle configurations having any number of open shells. The antisymmetrized matrix element was expanded as a sum of semisymmetrized matrix elements. The diagram representing a semisymmetrized matrix element was composed of four diagram blocks; the bra block, the ket block, the spectator block, and the interaction block. The first three blocks indicate the couplings of the two interacting configurations while the last depends on the interaction and is the replaceable component. Interaction blocks for relativistic operators and commonly used potentials were summarized in ready to use forms. A simple step by step procedure was prescribed generally for calculating antisymmetrized matrix elements of one and two particle operators.
Rolling Element Bearing Stiffness Matrix Determination (Presentation)
Guo, Y.; Parker, R.
2014-01-01
Current theoretical bearing models differ in their stiffness estimates because of different model assumptions. In this study, a finite element/contact mechanics model is developed for rolling element bearings with the focus of obtaining accurate bearing stiffness for a wide range of bearing types and parameters. A combined surface integral and finite element method is used to solve for the contact mechanics between the rolling elements and races. This model captures the time-dependent characteristics of the bearing contact due to the orbital motion of the rolling elements. A numerical method is developed to determine the full bearing stiffness matrix corresponding to two radial, one axial, and two angular coordinates; the rotation about the shaft axis is free by design. This proposed stiffness determination method is validated against experiments in the literature and compared to existing analytical models and widely used advanced computational methods. The fully-populated stiffness matrix demonstrates the coupling between bearing radial, axial, and tilting bearing deflections.
Distributions of off-diagonal scattering matrix elements: Exact results
Nock, A. Kumar, S. Sommers, H.-J. Guhr, T.
2014-03-15
Scattering is a ubiquitous phenomenon which is observed in a variety of physical systems which span a wide range of length scales. The scattering matrix is the key quantity which provides a complete description of the scattering process. The universal features of scattering in chaotic systems is most generally modeled by the Heidelberg approach which introduces stochasticity to the scattering matrix at the level of the Hamiltonian describing the scattering center. The statistics of the scattering matrix is obtained by averaging over the ensemble of random Hamiltonians of appropriate symmetry. We derive exact results for the distributions of the real and imaginary parts of the off-diagonal scattering matrix elements applicable to orthogonally-invariant and unitarily-invariant Hamiltonians, thereby solving a long standing problem. -- Highlights: •Scattering problem in complex or chaotic systems. •Heidelberg approach to model the chaotic nature of the scattering center. •A novel route to the nonlinear sigma model based on the characteristic function. •Exact results for the distributions of off-diagonal scattering-matrix elements. •Universal aspects of the scattering-matrix fluctuations.
Fermi matrix element with isospin breaking
NASA Astrophysics Data System (ADS)
Guichon, P. A. M.; Thomas, A. W.; Saito, K.
2011-02-01
Prompted by the level of accuracy now being achieved in tests of the unitarity of the CKM matrix, we consider the possible modification of the Fermi matrix element for the β-decay of a neutron, including possible in-medium and isospin violating corrections. While the nuclear modifications lead to very small corrections once the Behrends-Sirlin-Ademollo-Gatto theorem is respected, the effect of the u-d mass difference on the conclusion concerning Vud is no longer insignificant. Indeed, we suggest that the correction to the value of |+|+| is at the level of 10.
NLO matrix elements and truncated showers
NASA Astrophysics Data System (ADS)
Höche, Stefan; Krauss, Frank; Schönherr, Marek; Siegert, Frank
2011-08-01
In this publication, an algorithm is presented that combines the ME+PS approach to merge sequences of tree-level matrix elements into inclusive event samples [1] with the P owheg method, which combines exact next-to-leading order matrix element results with the parton shower [2, 3]. It was developed in parallel to the ME nloPS technique discussed in [4] and has been implemented in the event generator S herpa [5, 6]. The benefits of this approach are exemplified by some first predictions for a number of processes, namely the production of jets in e + e --annihilation, in deep-inelastic ep scattering, in association with single W, Z or Higgs bosons, and with vector boson pairs at hadron colliders.
Proton decay matrix elements from lattice QCD
Aoki, Yasumichi; Shintani, Eigo; Collaboration: RBC Collaboration; UKQCD Collaboration
2012-07-27
We report on the calculation of the matrix elements of nucleon to pseudoscalar decay through a three quark operator, a part of the low-energy, four-fermion, baryon-number-violating operator originating from grand unified theories. The direct calculation of the form factors using domain-wall fermions on the lattice, incorporating the u, d and s sea-quarks effects yields the results with all the relevant systematic uncertainties controlled for the first time.
NASA Astrophysics Data System (ADS)
De Corato, M.; Slot, J. J. M.; Hütter, M.; D'Avino, G.; Maffettone, P. L.; Hulsen, M. A.
2016-07-01
In this paper, we present a finite element implementation of fluctuating hydrodynamics with a moving boundary fitted mesh for treating the suspended particles. The thermal fluctuations are incorporated into the continuum equations using the Landau and Lifshitz approach [1]. The proposed implementation fulfills the fluctuation-dissipation theorem exactly at the discrete level. Since we restrict the equations to the creeping flow case, this takes the form of a relation between the diffusion coefficient matrix and friction matrix both at the particle and nodal level of the finite elements. Brownian motion of arbitrarily shaped particles in complex confinements can be considered within the present formulation. A multi-step time integration scheme is developed to correctly capture the drift term required in the stochastic differential equation (SDE) describing the evolution of the positions of the particles. The proposed approach is validated by simulating the Brownian motion of a sphere between two parallel plates and the motion of a spherical particle in a cylindrical cavity. The time integration algorithm and the fluctuating hydrodynamics implementation are then applied to study the diffusion and the equilibrium probability distribution of a confined circle under an external harmonic potential.
Measuring Sparticles with the Matrix Element
Alwall, Johan; Freitas, Ayres; Mattelaer, Olivier; /INFN, Rome3 /Rome III U. /Louvain U.
2012-04-10
We apply the Matrix Element Method (MEM) to mass determination of squark pair production with direct decay to quarks and LSP at the LHC, showing that simultaneous mass determination of squarks and LSP is possible. We furthermore propose methods for inclusion of QCD radiation effects in the MEM. The goal of the LHC at CERN, scheduled to start this year, is to discover new physics through deviations from the Standard Model (SM) predictions. After discovery of deviations from the SM, the next step will be classification of the new physics. An important first goal in this process will be establishing a mass spectrum of the new particles. One of the most challenging scenarios is pair-production of new particles which decay to invisible massive particles, giving missing energy signals. Many methods have been proposed for mass determination in such scenarios (for a recent list of references, see e.g. [1]). In this proceeding, we report the first steps in applying the Matrix Element Method (MEM) in the context of supersymmetric scenarios giving missing energy signals. After a quick review of the MEM, we will focus on squark pair production, a process where other mass determination techniques have difficulties to simultaneously determine the LSP and squark masses. Finally, we will introduce methods to extend the range of validity of the MEM, by taking into account initial state radiation (ISR) in the method.
Diagonal multisoliton matrix elements in finite volume
NASA Astrophysics Data System (ADS)
Pálmai, T.; Takács, G.
2013-02-01
We consider diagonal matrix elements of local operators between multisoliton states in finite volume in the sine-Gordon model and formulate a conjecture regarding their finite size dependence which is valid up to corrections exponential in the volume. This conjecture extends the results of Pozsgay and Takács which were only valid for diagonal scattering. In order to test the conjecture, we implement a numerical renormalization group improved truncated conformal space approach. The numerical comparisons confirm the conjecture, which is expected to be valid for general integrable field theories. The conjectured formula can be used to evaluate finite temperature one-point and two-point functions using recently developed methods.
Precision measurement of transition matrix elements via light shift cancellation.
Herold, C D; Vaidya, V D; Li, X; Rolston, S L; Porto, J V; Safronova, M S
2012-12-14
We present a method for accurate determination of atomic transition matrix elements at the 10(-3) level. Measurements of the ac Stark (light) shift around "magic-zero" wavelengths, where the light shift vanishes, provide precise constraints on the matrix elements. We make the first measurement of the 5s - 6p matrix elements in rubidium by measuring the light shift around the 421 and 423 nm zeros through diffraction of a condensate off a sequence of standing wave pulses. In conjunction with existing theoretical and experimental data, we find 0.3235(9)ea(0) and 0.5230(8)ea(0) for the 5s - 6p(1/2) and 5s - 6p(3/2) elements, respectively, an order of magnitude more accurate than the best theoretical values. This technique can provide needed, accurate matrix elements for many atoms, including those used in atomic clocks, tests of fundamental symmetries, and quantum information. PMID:23368314
Excited State Effects in Nucleon Matrix Element Calculations
Constantia Alexandrou, Martha Constantinou, Simon Dinter, Vincent Drach, Karl Jansen, Theodoros Leontiou, Dru B Renner
2011-12-01
We perform a high-statistics precision calculation of nucleon matrix elements using an open sink method allowing us to explore a wide range of sink-source time separations. In this way the influence of excited states of nucleon matrix elements can be studied. As particular examples we present results for the nucleon axial charge g{sub A} and for the first moment of the isovector unpolarized parton distribution x{sub u-d}. In addition, we report on preliminary results using the generalized eigenvalue method for nucleon matrix elements. All calculations are performed using N{sub f} = 2+1+1 maximally twisted mass Wilson fermions.
Analytical expressions for vibrational matrix elements of Morse oscillators
Zuniga, J.; Hidalgo, A.; Frances, J.M.; Requena, A.; Lopez Pineiro, A.; Olivares del Valle, F.J.
1988-10-15
Several exact recursion relations connecting different Morse oscillator matrix elements associated with the operators q/sup ..cap alpha../e/sup -//sup ..beta..//sup aq/ and q/sup ..cap alpha../e/sup -//sup ..beta..//sup aq/(d/dr) are derived. Matrix elements of the other useful operators may then be obtained easily. In particular, analytical expressions for (y/sup k/d/dr) and (y/sup k/d/dr+(d/dr)y/sup k/), matrix elements of interest in the study of the internuclear motion in polyatomic molecules, are obtained.
Analytic formula for quadrupole-quadrupole matrix elements
NASA Astrophysics Data System (ADS)
Rosensteel, G.
1990-12-01
An analytic formula is reported for general matrix elements of the microscopic quadrupole-quadrupole operator in the U(3)-boson approximation. The complete infinite-dimensional basis of A-fermion wave functions is compatible with the harmonic-oscillator shell model and consists of np-nh configurations, with spurious center-of-mass excitations removed, which are symmetry adapted to the Elliott U(3) and symplectic Sp(3,R) models. The formula expresses the general Q2.Q2 matrix element with respect to this complete orthonormal basis as a Racah SU(3) U coefficient times a closed-shell matrix element. An oscillator closed-shell matrix element of Q2.Q2 is a square root of a rational function of the integer quantum numbers of the U(3) basis.
[Hematopoietic microenvironment: cellular and extracellular matrix elements].
Minguell, J J; Fernández, M; Tetas, M; Martínez, J; Bruzzone, M; Rodríguez, J P
1988-06-01
In bone marrow, cellular stroma together with extracellular matrix (EM) provide an adequate microenvironment for the proliferation and differentiation of hemopoietic progenitor cells. In this article we describe studies on the cell characteristics of a main stromal phenotype, a fibroblast-like cell and its ability to produce in vitro EM components. Comparative studies were performed in fibroblast cultures derived from normal and acute lymphoblastic leukemic (ALL) bone marrow. The grow characteristics of fibroblasts from ALL marrow as well as its capacity to synthetize collagen, fibronectin and GAGs are impaired when compared to fibroblast from normal marrow. Thus, in ALL the impaired production of EM biomolecules by a transient damaged population of stromal cells, may contribute to the development of a defective microenvironment for hemopoiesis. PMID:3154858
Neutrinoless double-β decay and nuclear transition matrix elements
Rath, P. K.
2015-10-28
Within mechanisms involving the light Majorana neutrinos, squark-neutrino, Majorons, sterile neutrinos and heavy Majorana neutrino, nuclear transition matrix elements for the neutrinoless (β{sup −}β{sup −}){sub 0ν} decay of {sup 96}Zr, {sup 100}Mo, {sup 128,130}Te and {sup 150}Nd nuclei are calculated by employing the PHFB approach. Effects due to finite size of nucleons, higher order currents, short range correlations, and deformations of parent as well as daughter nuclei on the calculated matrix elements are estimated. Uncertainties in nuclear transition matrix elements within long-ranged mechanisms but for double Majoron accompanied (β{sup −}β{sup −}ϕϕ){sub 0ν} decay modes are 9%–15%. In the case of short ranged heavy Majorona neutrino exchange mechanism, the maximum uncertainty is about 35%. The maximum systematic error within the mechanism involving the exchange of light Majorana neutrino is about 46%.
Acceleration of matrix element computations for precision measurements
Brandt, Oleg; Gutierrez, Gaston; Wang, M. H.L.S.; Ye, Zhenyu
2014-11-25
The matrix element technique provides a superior statistical sensitivity for precision measurements of important parameters at hadron colliders, such as the mass of the top quark or the cross-section for the production of Higgs bosons. The main practical limitation of the technique is its high computational demand. Using the example of the top quark mass, we present two approaches to reduce the computation time of the technique by a factor of 90. First, we utilize low-discrepancy sequences for numerical Monte Carlo integration in conjunction with a dedicated estimator of numerical uncertainty, a novelty in the context of the matrix element technique. We then utilize a new approach that factorizes the overall jet energy scale from the matrix element computation, a novelty in the context of top quark mass measurements. The utilization of low-discrepancy sequences is of particular general interest, as it is universally applicable to Monte Carlo integration, and independent of the computing environment.
Neutrinoless double-β decay and nuclear transition matrix elements
NASA Astrophysics Data System (ADS)
Rath, P. K.
2015-10-01
Within mechanisms involving the light Majorana neutrinos, squark-neutrino, Majorons, sterile neutrinos and heavy Majorana neutrino, nuclear transition matrix elements for the neutrinoless (β-β-)0ν decay of 96Zr, 100Mo, 128,130Te and 150Nd nuclei are calculated by employing the PHFB approach. Effects due to finite size of nucleons, higher order currents, short range correlations, and deformations of parent as well as daughter nuclei on the calculated matrix elements are estimated. Uncertainties in nuclear transition matrix elements within long-ranged mechanisms but for double Majoron accompanied (β-β-ϕϕ)0ν decay modes are 9%-15%. In the case of short ranged heavy Majorona neutrino exchange mechanism, the maximum uncertainty is about 35%. The maximum systematic error within the mechanism involving the exchange of light Majorana neutrino is about 46%.
Some measurements for determining strangeness matrix elements in the nucleon
Henley, E.M.; Pollock, S.J.; Ying, S. ); Frederico, T. , Sao Jose dos Campos, SP . Inst. de Estudos Avancados); Krein, . Inst. de Fisica Teorica); Williams, A.G. )
1991-01-01
Some experiments to measure strangeness matrix elements of the proton are proposed. Two of these suggestions are described in some detail, namely electro-production of phi mesons and the difference between neutrino and antineutrino scattering for isospin zero targets such as deuterium.
Some measurements for determining strangeness matrix elements in the nucleon
Henley, E.M.; Pollock, S.J.; Ying, S.; Frederico, T.; Krein,; Williams, A.G.
1991-12-31
Some experiments to measure strangeness matrix elements of the proton are proposed. Two of these suggestions are described in some detail, namely electro-production of phi mesons and the difference between neutrino and antineutrino scattering for isospin zero targets such as deuterium.
Expectation Matrix Based Quantum Dynamics of a Univariate System at the Zero Fluctuation Limit
Demiralp, Metin
2007-12-26
The variation of the expectation matrix of position and momentum operator in time can serve us to investigate the evolution of a quantum system in time. This brings the utilization of the ODEs instead of Schroedinger's equation at the expense of incapability for the calculation of the wave function. As long as we deal with the observables which can be expressed in terms of position and momentum operators this may be quite practical to know about the quantum dynamics of the system under consideration. Expectation matrix of an operator becomes a function of the expectation matrices of the position and momentum operator when the fluctuations diminish to zero. At this limit, the time-variant ODEs for the expectation matrices of the position and momentum operator can be handled by using the matrix algebraic tools even in the case of nonlinearities in the potential function. This work presents certain details about these points.
Aggelen, Helen van; Department of Chemistry, Duke University, Durham, North Carolina 27708 ; Yang, Yang; Yang, Weitao
2014-05-14
Despite their unmatched success for many applications, commonly used local, semi-local, and hybrid density functionals still face challenges when it comes to describing long-range interactions, static correlation, and electron delocalization. Density functionals of both the occupied and virtual orbitals are able to address these problems. The particle-hole (ph-) Random Phase Approximation (RPA), a functional of occupied and virtual orbitals, has recently known a revival within the density functional theory community. Following up on an idea introduced in our recent communication [H. van Aggelen, Y. Yang, and W. Yang, Phys. Rev. A 88, 030501 (2013)], we formulate more general adiabatic connections for the correlation energy in terms of pairing matrix fluctuations described by the particle-particle (pp-) propagator. With numerical examples of the pp-RPA, the lowest-order approximation to the pp-propagator, we illustrate the potential of density functional approximations based on pairing matrix fluctuations. The pp-RPA is size-extensive, self-interaction free, fully anti-symmetric, describes the strong static correlation limit in H{sub 2}, and eliminates delocalization errors in H{sub 2}{sup +} and other single-bond systems. It gives surprisingly good non-bonded interaction energies – competitive with the ph-RPA – with the correct R{sup −6} asymptotic decay as a function of the separation R, which we argue is mainly attributable to its correct second-order energy term. While the pp-RPA tends to underestimate absolute correlation energies, it gives good relative energies: much better atomization energies than the ph-RPA, as it has no tendency to underbind, and reaction energies of similar quality. The adiabatic connection in terms of pairing matrix fluctuation paves the way for promising new density functional approximations.
Magic wavelengths, matrix elements, polarizabilities, and lifetimes of Cs
NASA Astrophysics Data System (ADS)
Safronova, M. S.; Safronova, U. I.; Clark, Charles W.
2016-07-01
Motivated by recent interest in their applications, we report a systematic study of Cs atomic properties calculated by a high-precision relativistic all-order method. Excitation energies, reduced matrix elements, transition rates, and lifetimes are determined for levels with principal quantum numbers n ≤12 and orbital angular momentum quantum numbers l ≤3 . Recommended values and estimates of uncertainties are provided for a number of electric-dipole transitions and the electric dipole polarizabilities of the n s , n p , and n d states. We also report a calculation of the electric quadrupole polarizability of the ground state. We display the dynamic polarizabilities of the 6 s and 7 p states for optical wavelengths between 1160 and 1800 nm and identify corresponding magic wavelengths for the 6 s -7 p1 /2 and 6 s -7 p3 /2 transitions. The values of relevant matrix elements needed for polarizability calculations at other wavelengths are provided.
Importance of Matrix Elements in the ARPES Spectra of BISCO
Bansil, A.; Lindroos, M.
1999-12-13
We have carried out extensive first-principles angle-resolved photointensity (ARPES) simulations in Bi2212 wherein the photoemission process is modeled realistically by taking into account the full crystal wave functions of the initial and final states in the presence of the surface. The spectral weight of the ARPES feature associated with the CuO{sub 2} plane bands is found to undergo large and systematic variations with k{sub (parallel} {sub sign)} as well as the energy and polarization of the incident photons. These theoretical predictions are in good accord with the corresponding measurements, indicating that the remarkable observed changes in the spectral weights in Bi2212 are essentially a matrix element effect and that the importance of matrix elements should be kept in mind in analyzing the ARPES spectra in the high T{sub c} 's. (c) 1999 The American Physical Society.
Importance of Matrix Elements in the ARPES Spectra of BISCO
NASA Astrophysics Data System (ADS)
Bansil, A.; Lindroos, M.
1999-12-01
We have carried out extensive first-principles angle-resolved photointensity (ARPES) simulations in Bi2212 wherein the photoemission process is modeled realistically by taking into account the full crystal wave functions of the initial and final states in the presence of the surface. The spectral weight of the ARPES feature associated with the CuO2 plane bands is found to undergo large and systematic variations with k∥ as well as the energy and polarization of the incident photons. These theoretical predictions are in good accord with the corresponding measurements, indicating that the remarkable observed changes in the spectral weights in Bi2212 are essentially a matrix element effect and that the importance of matrix elements should be kept in mind in analyzing the ARPES spectra in the high Tc's.
Glueball Spectrum and Matrix Elements on Anisotropic Lattices
Y. Chen; A. Alexandru; S.J. Dong; T. Draper; I. Horvath; F.X. Lee; K.F. Liu; N. Mathur; C. Morningstar; M. Peardon; S. Tamhankar; B.L. Young; J.B. Zhang
2006-01-01
The glueball-to-vacuum matrix elements of local gluonic operators in scalar, tensor, and pseudoscalar channels are investigated numerically on several anisotropic lattices with the spatial lattice spacing ranging from 0.1fm - 0.2fm. These matrix elements are needed to predict the glueball branching ratios in J/{psi} radiative decays which will help identify the glueball states in experiments. Two types of improved local gluonic operators are constructed for a self-consistent check and the finite volume effects are studied. We find that lattice spacing dependence of our results is very weak and the continuum limits are reliably extrapolated, as a result of improvement of the lattice gauge action and local operators. We also give updated glueball masses with various quantum numbers.
A stochastic method for computing hadronic matrix elements
Alexandrou, Constantia; Constantinou, Martha; Dinter, Simon; Drach, Vincent; Jansen, Karl; Hadjiyiannakou, Kyriakos; Renner, Dru B.
2014-01-24
In this study, we present a stochastic method for the calculation of baryon 3-point functions which is an alternative to the typically used sequential method offering more versatility. We analyze the scaling of the error of the stochastically evaluated 3-point function with the lattice volume and find a favorable signal to noise ratio suggesting that the stochastic method can be extended to large volumes providing an efficient approach to compute hadronic matrix elements and form factors.
Algebraic evaluation of matrix elements in the Laguerre function basis
NASA Astrophysics Data System (ADS)
McCoy, A. E.; Caprio, M. A.
2016-02-01
The Laguerre functions constitute one of the fundamental basis sets for calculations in atomic and molecular electron-structure theory, with applications in hadronic and nuclear theory as well. While similar in form to the Coulomb bound-state eigenfunctions (from the Schrödinger eigenproblem) or the Coulomb-Sturmian functions (from a related Sturm-Liouville problem), the Laguerre functions, unlike these former functions, constitute a complete, discrete, orthonormal set for square-integrable functions in three dimensions. We construct the SU(1, 1) × SO(3) dynamical algebra for the Laguerre functions and apply the ideas of factorization (or supersymmetric quantum mechanics) to derive shift operators for these functions. We use the resulting algebraic framework to derive analytic expressions for matrix elements of several basic radial operators (involving powers of the radial coordinate and radial derivative) in the Laguerre function basis. We illustrate how matrix elements for more general spherical tensor operators in three dimensional space, such as the gradient, may then be constructed from these radial matrix elements.
Weak matrix elements on the lattice - Circa 1995
Soni, A.
1995-10-03
Status of weak matrix elements is reviewed. In particular, e{prime}/e, B {yields} K*{gamma}, B{sub B} and B{sub B}, are discussed and the overall situation with respect to the lattice effort and some of its phenomenological implications are summarised. For e{prime}/e the need for the relevant matrix elements is stressed in view of the forthcoming improved experiments. For some of the operators, (e.g. O{sub 6}), even bound on their matrix elements would be very helpful. On B {yields} K{degrees}{gamma}, a constant behavior of T{sub 2} appears disfavored although dependence of T{sub 2} could, of course, be milder than a simple pole. Improved data is badly needed to settle this important issue firmly, especially in view of its ramification for extractions of V{sub td} from B {yields} {rho}{gamma}. On B{sub {kappa}}, the preliminary result from JLQCD appears to contradict Sharpe et al. JLQCD data seems to fit very well to linear {alpha} dependence and leads to an appreciably lower value of B{sub {kappa}}. Four studies of B{sub {kappa}} in the {open_quotes}full{close_quotes} (n{sub f} = 2) theory indicate very little quenching effects on B{sub {kappa}}; the full theory value seems to be just a little less than the quenched result. Based on expectations from HQET, analysis of B-parameter (B{sub h}{ell}) for the heavy-light mesons via B{sub h}{ell}) = constant + constants{prime}/m{sub h}{ell} is suggested. A summary of an illustrative sample of hadron matrix elements is given and constraints on CKM parameters (e.g. V{sub td}/V{sub ts}, on the unitarity triangle and on x{sub s}/x{sub d}, emerging from the lattice calculations along with experimental results are briefly discussed. In quite a few cases, for the first time, some indication of quenching errors on weak matrix elements are now becoming available.
Multi-jet Merging with NLO Matrix Elements
Siegert, Frank; Hoche, Stefan; Krauss, Frank; Schonherr, Marek; /Dresden, Tech. U.
2011-08-18
In the algorithm presented here, the ME+PS approach to merge samples of tree-level matrix elements into inclusive event samples is combined with the POWHEG method, which includes exact next-to-leading order matrix elements in the parton shower. The advantages of the method are discussed and the quality of its implementation in SHERPA is exemplified by results for e{sup +}e{sup -} annihilation into hadrons at LEP, for deep-inelastic lepton-nucleon scattering at HERA, for Drell-Yan lepton-pair production at the Tevatron and for W{sup +}W{sup -}-production at LHC energies. The simulation of hard QCD radiation in parton-shower Monte Carlos has seen tremendous progress over the last years. It was largely stimulated by the need for more precise predictions at LHC energies where the large available phase space allows additional hard QCD radiation alongside known Standard Model processes or even signals from new physics. Two types of algorithms have been developed, which allow to improve upon the soft-collinear approximations made in the parton shower, such that hard radiation is simulated according to exact matrix elements. In the ME+PS approach [1] higher-order tree-level matrix elements for different final-state jet multiplicity are merged with each other and with subsequent parton shower emissions to generate an inclusive sample. Such a prescription is invaluable for analyses which are sensitive to final states with a large jet multiplicity. The only remaining deficiency of such tree-level calculations is the large uncertainty stemming from scale variations. The POWHEG method [2] solves this problem for the lowest multiplicity subprocess by combining full NLO matrix elements with the parton shower. While this leads to NLO accuracy in the inclusive cross section and the exact radiation pattern for the first emission, it fails to describe higher-order emissions with improved accuracy. Thus it is not sufficient if final states with high jet multiplicities are considered
Acceleration of matrix element computations for precision measurements
Brandt, Oleg; Gutierrez, Gaston; Wang, M. H.L.S.; Ye, Zhenyu
2014-11-25
The matrix element technique provides a superior statistical sensitivity for precision measurements of important parameters at hadron colliders, such as the mass of the top quark or the cross-section for the production of Higgs bosons. The main practical limitation of the technique is its high computational demand. Using the example of the top quark mass, we present two approaches to reduce the computation time of the technique by a factor of 90. First, we utilize low-discrepancy sequences for numerical Monte Carlo integration in conjunction with a dedicated estimator of numerical uncertainty, a novelty in the context of the matrix elementmore » technique. We then utilize a new approach that factorizes the overall jet energy scale from the matrix element computation, a novelty in the context of top quark mass measurements. The utilization of low-discrepancy sequences is of particular general interest, as it is universally applicable to Monte Carlo integration, and independent of the computing environment.« less
SYMBMAT: Symbolic computation of quantum transition matrix elements
NASA Astrophysics Data System (ADS)
Ciappina, M. F.; Kirchner, T.
2012-08-01
We have developed a set of Mathematica notebooks to compute symbolically quantum transition matrices relevant for atomic ionization processes. The utilization of a symbolic language allows us to obtain analytical expressions for the transition matrix elements required in charged-particle and laser induced ionization of atoms. Additionally, by using a few simple commands, it is possible to export these symbolic expressions to standard programming languages, such as Fortran or C, for the subsequent computation of differential cross sections or other observables. One of the main drawbacks in the calculation of transition matrices is the tedious algebraic work required when initial states other than the simple hydrogenic 1s state need to be considered. Using these notebooks the work is dramatically reduced and it is possible to generate exact expressions for a large set of bound states. We present explicit examples of atomic collisions (in First Born Approximation and Distorted Wave Theory) and laser-matter interactions (within the Dipole and Strong Field Approximations and different gauges) using both hydrogenic wavefunctions and Slater-Type Orbitals with arbitrary nlm quantum numbers as initial states. Catalogue identifier: AEMI_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEMI_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC license, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 71 628 No. of bytes in distributed program, including test data, etc.: 444 195 Distribution format: tar.gz Programming language: Mathematica Computer: Single machines using Linux or Windows (with cores with any clock speed, cache memory and bits in a word) Operating system: Any OS that supports Mathematica. The notebooks have been tested under Windows and Linux and with versions 6.x, 7.x and 8.x Classification: 2.6 Nature of problem
Calculation of hadronic matrix elements using lattice QCD
Gupta, R.
1993-08-01
The author gives a brief introduction to the scope of lattice QCD calculations in his effort to extract the fundamental parameters of the standard model. This goal is illustrated by two examples. First the author discusses the extraction of CKM matrix elements from measurements of form factors for semileptonic decays of heavy-light pseudoscalar mesons such as D {yields} Ke{nu}. Second, he presents the status of results for the kaon B parameter relevant to CP violation. He concludes the talk with a short outline of his experiences with optimizing QCD codes on the CM5.
Matrix element analyses of dark matter scattering and annihilation
NASA Astrophysics Data System (ADS)
Kumar, Jason; Marfatia, Danny
2013-07-01
We provide a compendium of results at the level of matrix elements for a systematic study of dark matter scattering and annihilation. We identify interactions that yield spin-dependent and spin-independent scattering and specify whether the interactions are velocity and/or momentum suppressed. We identify the interactions that lead to s-wave or p-wave annihilation, and those that are chirality suppressed. We also list the interaction structures that can interfere in scattering and annihilation processes. Using these results, we point out situations in which deviations from the standard lore are obtained.
NASA Astrophysics Data System (ADS)
Spohn, Herbert
2011-03-01
In 1986 Kardar, Parisi, and Zhang (KPZ) proposed a stochastic evolution equation for growing interfaces, thereby triggering an intense study of growth processes with local growth rules. Specifically we have in mind the recent spectacular experiment of Takeuchi and Sano on droplet growth in a thin film of turbulent liquid crystal. Over the last ten years one has studied universal probability density functions on the basis of simplified lattice growth models. Surprisingly enough the one-point shape fluctuations are governed by the same statistical laws as the largest eigenvalue of a random matrix, Gaussian Unitary Ensemble (GUE) in case of a curved front and Gaussian Orthogonal Ensemble (GOE) for a flat front. Recently we obtained the first exact solution of the KPZ equation for initial conditions corresponding to droplet growth, thereby providing the probability density function for the height at any time. For long times we recover the universal statistical properties as computed from lattice growth models.
Gorissen, Mieke; Hooyberghs, Jef; Vanderzande, Carlo
2009-02-01
Cumulants of a fluctuating current can be obtained from a free-energy-like generating function, which for Markov processes equals the largest eigenvalue of a generalized generator. We determine this eigenvalue with the density-matrix renormalization group for stochastic systems. We calculate the variance of the current in the different phases, and at the phase transitions, of the totally asymmetric exclusion process. Our results can be described in the terms of a scaling ansatz that involves the dynamical exponent z . We also calculate the generating function of the dynamical activity (total number of configuration changes) near the absorbing-state transition of the contact process. Its scaling properties can be expressed in terms of known critical exponents. PMID:19391693
NASA Astrophysics Data System (ADS)
Gligor, M.; Ausloos, M.
2007-05-01
The statistical distances between countries, calculated for various moving average time windows, are mapped into the ultrametric subdominant space as in classical Minimal Spanning Tree methods. The Moving Average Minimal Length Path (MAMLP) algorithm allows a decoupling of fluctuations with respect to the mass center of the system from the movement of the mass center itself. A Hamiltonian representation given by a factor graph is used and plays the role of cost function. The present analysis pertains to 11 macroeconomic (ME) indicators, namely the GDP (x1), Final Consumption Expenditure (x2), Gross Capital Formation (x3), Net Exports (x4), Consumer Price Index (y1), Rates of Interest of the Central Banks (y2), Labour Force (z1), Unemployment (z2), GDP/hour worked (z3), GDP/capita (w1) and Gini coefficient (w2). The target group of countries is composed of 15 EU countries, data taken between 1995 and 2004. By two different methods (the Bipartite Factor Graph Analysis and the Correlation Matrix Eigensystem Analysis) it is found that the strongly correlated countries with respect to the macroeconomic indicators fluctuations can be partitioned into stable clusters.
Importance of Matrix Elements in the ARPES Spectra of BISCO
NASA Astrophysics Data System (ADS)
Bansil, A.; Lindroos, M.
2000-03-01
We have carried out extensive first-principles angle-resolved photointensity (ARPES) simulations in Bi2212 wherein the photoemission process is modelled realistically by taking into account the full crystal wavefunctions of the initial and final states in the presence of the surface.(A. Bansil and M. Lindroos, Phys. Rev. Letters (Dec 13, 1999)) The spectral weight of the ARPES feature associated with the CuO2 plane bands is found to undergo large and systematic variations with k_allel as well as the energy and polarization of the incident photons. These theoretical predictions are in good accord with the corresponding measurements, indicating that the remarkable observed changes in the spectral weights in Bi2212 are essentially a matrix element effect and that the importance of matrix elements should be kept in mind in analyzing the ARPES spectra in the high-Tc's. Another notable implication of this work is that the integral (over energy) of the ARPES intensity does not yield the momentum density of the electron gas. We will also discuss some of our simulations aimed at gaining insight into the connectivity of the Fermi surface in Bi2212 around the M-point, the effects of modulations, and related issues. Work supported in part by the U.S.D.O.E.
SYMBMAT: Symbolic computation of quantum transition matrix elements
NASA Astrophysics Data System (ADS)
Ciappina, M. F.; Kirchner, T.
2012-08-01
We have developed a set of Mathematica notebooks to compute symbolically quantum transition matrices relevant for atomic ionization processes. The utilization of a symbolic language allows us to obtain analytical expressions for the transition matrix elements required in charged-particle and laser induced ionization of atoms. Additionally, by using a few simple commands, it is possible to export these symbolic expressions to standard programming languages, such as Fortran or C, for the subsequent computation of differential cross sections or other observables. One of the main drawbacks in the calculation of transition matrices is the tedious algebraic work required when initial states other than the simple hydrogenic 1s state need to be considered. Using these notebooks the work is dramatically reduced and it is possible to generate exact expressions for a large set of bound states. We present explicit examples of atomic collisions (in First Born Approximation and Distorted Wave Theory) and laser-matter interactions (within the Dipole and Strong Field Approximations and different gauges) using both hydrogenic wavefunctions and Slater-Type Orbitals with arbitrary nlm quantum numbers as initial states. Catalogue identifier: AEMI_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEMI_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC license, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 71 628 No. of bytes in distributed program, including test data, etc.: 444 195 Distribution format: tar.gz Programming language: Mathematica Computer: Single machines using Linux or Windows (with cores with any clock speed, cache memory and bits in a word) Operating system: Any OS that supports Mathematica. The notebooks have been tested under Windows and Linux and with versions 6.x, 7.x and 8.x Classification: 2.6 Nature of problem
NASA Astrophysics Data System (ADS)
Li, L.; Wang, K.; Li, H.; Eibert, T. F.
2014-11-01
A hybrid higher-order finite element boundary integral (FE-BI) technique is discussed where the higher-order FE matrix elements are computed by a fully analytical procedure and where the gobal matrix assembly is organized by a self-identifying procedure of the local to global transformation. This assembly procedure applys to both, the FE part as well as the BI part of the algorithm. The geometry is meshed into three-dimensional tetrahedra as finite elements and nearly orthogonal hierarchical basis functions are employed. The boundary conditions are implemented in a strong sense such that the boundary values of the volume basis functions are directly utilized within the BI, either for the tangential electric and magnetic fields or for the asssociated equivalent surface current densities by applying a cross product with the unit surface normals. The self-identified method for the global matrix assembly automatically discerns the global order of the basis functions for generating the matrix elements. Higher order basis functions do need more unknowns for each single FE, however, fewer FEs are needed to achieve the same satisfiable accuracy. This improvement provides a lot more flexibility for meshing and allows the mesh size to raise up to λ/3. The performance of the implemented system is evaluated in terms of computation time, accuracy and memory occupation, where excellent results with respect to precision and computation times of large scale simulations are found.
Precision Study of Excited State Effects in Nucleon Matrix Elements
Simon Dinter, Constantia Alexandrou, Martha Constantinou, Vincent Drach, Karl Jansen, Dru B. Renner
2011-10-01
We present a dedicated precision analysis of the influence of excited states on the calculation of several nucleon matrix elements. This calculation is performed at fixed values of the lattice spacing, volume and pion mass that are typical of contemporary lattice computations. We focus on the nucleon axial charge, g{sub A}, for which we use 7,500 measurements, and on the average momentum of the unpolarized isovector parton distribution, x{sub u-d}, for which we use 23,000 measurements. All computations are done employing N{sub f}=2+1+1 maximally-twisted-mass Wilson fermions and non-perturbatively calculated renormalization factors. We find that excited state effects are negligible for g{sub A} and lead to a O(10%) downward shift for x{sub u-d}.
Fabrication of synthetic diffractive elements using advanced matrix laser lithography
NASA Astrophysics Data System (ADS)
Škereň, M.; Svoboda, J.; Květoň, M.; Fiala, P.
2013-02-01
In this paper we present a matrix laser writing device based on a demagnified projection of a micro-structure from a computer driven spatial light modulator. The device is capable of writing completely aperiodic micro-structures with resolution higher than 200 000 DPI. An optical system is combined with ultra high precision piezoelectric stages with an elementary step ~ 4 nm. The device operates in a normal environment, which significantly decreases the costs compared to competitive technologies. Simultaneously, large areas can be exposed up to 100 cm2. The capabilities of the constructed device will be demonstrated on particular elements fabricated for real applications. The optical document security is the first interesting field, where the synthetic image holograms are often combined with sophisticated aperiodic micro-structures. The proposed technology can easily write simple micro-gratings creating the color and kinetic visual effects, but also the diffractive cryptograms, waveguide couplers, and other structures recently used in the field of optical security. A general beam shaping elements and special photonic micro-structures are another important applications which will be discussed in this paper.
Depolarizing differential Mueller matrix of homogeneous media under Gaussian fluctuation hypothesis.
Devlaminck, Vincent
2015-10-01
In this paper, we address the issue of the existence of a solution of depolarizing differential Mueller matrix for a homogeneous medium. Such a medium is characterized by linear changes of its differential optical properties with z the thickness of the medium. We show that, under a short correlation distance assumption, it is possible to derive such linear solution, and we clarify this solution in the particular case where the random fluctuation processes associated to the optical properties are Gaussian white noise-like. A solution to the problem of noncommutativity of a previously proposed model [J. Opt. Soc. Am.30, 2196 (2013)JOSAAH0030-394110.1364/JOSAA.30.002196] is given by assuming a random permutation of the order of the layers and by averaging all the differential matrices resulting from these permutations. It is shown that the underlying assumption in this case is exactly the Gaussian white noise assumption. Finally, a recently proposed approach [Opt. Lett.39, 4470 (2014)OPLEDP0146-959210.1364/OL.39.004470] for analysis of the statistical properties related to changes in optical properties is revisited, and the experimental conditions of application of these results are specified. PMID:26479926
NASA Technical Reports Server (NTRS)
Newman, M. B.; Filstrup, A. W.
1973-01-01
Linear (8 node), parabolic (20 node), cubic (32 node) and mixed (some edges linear, some parabolic and some cubic) have been inserted into NASTRAN, level 15.1. First the dummy element feature was used to check out the stiffness matrix generation routines for the linear element in NASTRAN. Then, the necessary modules of NASTRAN were modified to include the new family of elements. The matrix assembly was changed so that the stiffness matrix of each isoparametric element is only generated once as the time to generate these higher order elements tends to be much longer than the other elements in NASTRAN. This paper presents some of the experiences and difficulties of inserting a new element or family of elements into NASTRAN.
Shell model nuclear matrix elements for competing mechanisms contributing to double beta decay
Horoi, Mihai
2013-12-30
Recent progress in the shell model approach to the nuclear matrix elements for the double beta decay process are presented. This includes nuclear matrix elements for competing mechanisms to neutrionless double beta decay, a comparison between closure and non-closure approximation for {sup 48}Ca, and an updated shell model analysis of nuclear matrix elements for the double beta decay of {sup 136}Xe.
Matrix element method for high performance computing platforms
NASA Astrophysics Data System (ADS)
Grasseau, G.; Chamont, D.; Beaudette, F.; Bianchini, L.; Davignon, O.; Mastrolorenzo, L.; Ochando, C.; Paganini, P.; Strebler, T.
2015-12-01
Lot of efforts have been devoted by ATLAS and CMS teams to improve the quality of LHC events analysis with the Matrix Element Method (MEM). Up to now, very few implementations try to face up the huge computing resources required by this method. We propose here a highly parallel version, combining MPI and OpenCL, which makes the MEM exploitation reachable for the whole CMS datasets with a moderate cost. In the article, we describe the status of two software projects under development, one focused on physics and one focused on computing. We also showcase their preliminary performance obtained with classical multi-core processors, CUDA accelerators and MIC co-processors. This let us extrapolate that with the help of 6 high-end accelerators, we should be able to reprocess the whole LHC run 1 within 10 days, and that we have a satisfying metric for the upcoming run 2. The future work will consist in finalizing a single merged system including all the physics and all the parallelism infrastructure, thus optimizing implementation for best hardware platforms.
A top quark mass measurement using a matrix element method
Linacre, Jacob Thomas
2009-01-01
A measurement of the mass of the top quark is presented, using top-antitop pair (t$\\bar{t}$) candidate events for the lepton+jets decay channel. The measurement makes use of Tevatron p$\\bar{p}$ collision data at centre-of-mass energy √s = 1.96 TeV, collected at the CDF detector. The top quark mass is measured by employing an unbinned maximum likelihood method where the event probability density functions are calculated using signal (t$\\bar{t}$) and background (W+jets) matrix elements, as well as a set of parameterised jet-to-parton mapping functions. The likelihood function is maximised with respect to the top quark mass, the fraction of signal events, and a correction to the jet energy scale (JES) of the calorimeter jets. The simultaneous measurement of the JES correction (Δ_{JES}) provides an in situ jet energy calibration based on the known mass of the hadronically decaying W boson. Using 578 lepton+jets candidate events corresponding to 3.2 fb ^{-1} of integrated luminosity, the top quark mass is measured to be m_{t} = 172.4± 1.4 (stat+Δ_{JES}) ±1.3 (syst) GeV=c^{2}, one of the most precise single measurements to date.
Controlling excited-state contamination in nucleon matrix elements
NASA Astrophysics Data System (ADS)
Yoon, Boram; Gupta, Rajan; Bhattacharya, Tanmoy; Engelhardt, Michael; Green, Jeremy; Joó, Bálint; Lin, Huey-Wen; Negele, John; Orginos, Kostas; Pochinsky, Andrew; Richards, David; Syritsyn, Sergey; Winter, Frank; Nucleon Matrix Elements NME Collaboration
2016-06-01
We present a detailed analysis of methods to reduce statistical errors and excited-state contamination in the calculation of matrix elements of quark bilinear operators in nucleon states. All the calculations were done on a 2 +1 -flavor ensemble with lattices of size 323×64 generated using the rational hybrid Monte Carlo algorithm at a =0.081 fm and with Mπ=312 MeV . The statistical precision of the data is improved using the all-mode-averaging method. We compare two methods for reducing excited-state contamination: a variational analysis and a 2-state fit to data at multiple values of the source-sink separation tsep. We show that both methods can be tuned to significantly reduce excited-state contamination and discuss their relative advantages and cost effectiveness. A detailed analysis of the size of source smearing used in the calculation of quark propagators and the range of values of tsep needed to demonstrate convergence of the isovector charges of the nucleon to the tsep→∞ estimates is presented.
Configuration interaction matrix elements for the quantum Hall effect
NASA Astrophysics Data System (ADS)
Wooten, Rachel; Macek, Joseph
2015-03-01
In the spherical model of the quantum Hall system, the two-body matrix elements and pseudopotentials can be found analytically in terms of a general scalar pair interaction potential by expressing the pair interaction as a weighted sum over Legendre polynomials. For non-infinite systems, only a finite set of terms in the potential expansion contribute to the interactions; the contributing terms define an effective spatial potential for the system. The connection between the effective spatial potential and the pseudopotential is one-to-one for finite systems, and any completely defined model pseudopotential can be analytically inverted to give a unique corresponding spatial potential. This technique of inverting the pseudopotential to derive effective spatial potentials may be of use for developing accurate model spatial potentials for quantum Monte Carlo simulations. We demonstrate the technique and the corresponding spatial potentials for a few example model pseudopotentials. Supported by Office of Basic Energy Sciences, U.S. DOE, Grant DE-FG02-02ER15283 to the University of Tennessee.
Spin correlations and velocity scaling in color-octet nonrelativistic QCD matrix elements
NASA Astrophysics Data System (ADS)
Bodwin, Geoffrey T.; Lee, Jungil; Sinclair, D. K.
2005-07-01
We compute spin-dependent decay matrix elements for S-wave charmonium and bottomonium in lattice nonrelativistic quantum chromodynamics (NRQCD). Particular emphasis is placed upon the color-octet matrix elements, since the corresponding production matrix elements are expected to appear in the dominant contributions to the production cross sections at large transverse momenta. We use three slightly different versions of the heavy-quark lattice Green’s functions in order to minimize the contributions that scale as powers of the ultraviolet cutoff. The lattice matrix elements that we calculate obey the hierarchy that is suggested by the velocity-scaling rules of NRQCD.
Analytical Derivations of Single-Particle Matrix Elements in Nuclear Shell Model
NASA Astrophysics Data System (ADS)
Fatah, Aziz H.; Radhi, R. A.; Abdullah, Nzar R.
2016-07-01
We present analytical method to calculate single particle matrix elements used in atomic and nuclear physics. We show seven different formulas of matrix elements of the operator f(r)dr m where f(r) = rμ, rμ jJ(qr), V(r) corresponding to the Gaussian and the Yukawa potentials used in nuclear shell models and nuclear structure. In addition, we take into account a general integral formula of the matrix element
NASA Astrophysics Data System (ADS)
Yang, Yang; van Aggelen, Helen; Yang, Weitao
2014-03-01
Double, Rydberg and charge transfer (CT) excitations have been great challenges for time-dependent density functional theory (TDDFT). Starting from an (N +/- 2) -electron single-determinant reference, we investigate excitations for the N-electron system through the pairing matrix fluctuation, which contains information on two-electron addition/removal processes. We adopt the particle-particle random phase approximation (pp-RPA) and the particle-particle Tamm-Dancoff approximation (pp-TDA) to approximate the pairing matrix fluctuation and then determine excitation energies by the differences of two-electron addition/removal energies. This approach captures all types of interesting excitations: single and double excitations are described accurately, Rydberg excitations are in good agreement with experimental data and CT excitations display correct 1/R dependence. Furthermore, the pp-RPA and the pp-TDA have a computational cost similar to TDDFT and consequently are promising for practical calculations.
Yang, Yang; Aggelen, Helen van; Yang, Weitao
2013-12-14
Double, Rydberg, and charge transfer (CT) excitations have been great challenges for time-dependent density functional theory (TDDFT). Starting from an (N ± 2)-electron single-determinant reference, we investigate excitations for the N-electron system through the pairing matrix fluctuation, which contains information on two-electron addition/removal processes. We adopt the particle-particle random phase approximation (pp-RPA) and the particle-particle Tamm-Dancoff approximation (pp-TDA) to approximate the pairing matrix fluctuation and then determine excitation energies by the differences of two-electron addition/removal energies. This approach captures all types of interesting excitations: single and double excitations are described accurately, Rydberg excitations are in good agreement with experimental data and CT excitations display correct 1/R dependence. Furthermore, the pp-RPA and the pp-TDA have a computational cost similar to TDDFT and consequently are promising for practical calculations.
NASA Astrophysics Data System (ADS)
Yang, Yang; van Aggelen, Helen; Yang, Weitao
2013-12-01
Double, Rydberg, and charge transfer (CT) excitations have been great challenges for time-dependent density functional theory (TDDFT). Starting from an (N ± 2)-electron single-determinant reference, we investigate excitations for the N-electron system through the pairing matrix fluctuation, which contains information on two-electron addition/removal processes. We adopt the particle-particle random phase approximation (pp-RPA) and the particle-particle Tamm-Dancoff approximation (pp-TDA) to approximate the pairing matrix fluctuation and then determine excitation energies by the differences of two-electron addition/removal energies. This approach captures all types of interesting excitations: single and double excitations are described accurately, Rydberg excitations are in good agreement with experimental data and CT excitations display correct 1/R dependence. Furthermore, the pp-RPA and the pp-TDA have a computational cost similar to TDDFT and consequently are promising for practical calculations.
NASA Astrophysics Data System (ADS)
Kinghorn, Donald Bruce
The matrix differential calculus is introduced to the quantum chemistry community via new matrix derivations of integral formulas and gradients for Hamiltonian matrix elements in a basis of correlated Gaussian functions. Requisite mathematical background material on Kronecker products, Hadamard products, the vec and vech operators, linear structures, and matrix differential calculus is presented. New matrix forms for the kinetic and potential energy operators are presented. Integrals for overlap, kinetic energy and potential energy matrix elements are derived in matrix form using matrix calculus. The gradient of the energy functional with respect to the correlated Gaussian exponent matrices is derived. Burdensome summation notation is entirely replaced with a compact matrix notation that is both theoretically and computationally insightful. These new formulas in the basis of explicitly correlated Gaussian basis functions, are implemented and applied to find variational upper bounds for non-relativistic ground states of ^4He, ^{infty}He, Ps_2, ^9Be, and ^ {infty}Be. Analytic gradients of the energy are included to speed optimization of the exponential variational parameters. Five different nonlinear optimization subroutines (algorithms) are compared: TN, truncated Newton; DUMING, quasi-Newton; DUMIDH, modified Newton; DUMCGG, conjugate gradient; POWELL, direction set (non-gradient). The new analytic gradient formulas are found to significantly accelerate optimizations that require gradients. The truncated Newton algorithm is found to outperform the other optimizers for the selected test cases. Computer timings and energy bounds are reported. The new TN bounds surpass previously reported bounds with the same basis size.
Application of the Finite-Element Z-Matrix Method to e-H2 Collisions
NASA Technical Reports Server (NTRS)
Huo, Winifred M.; Brown, David; Langhoff, Stephen R. (Technical Monitor)
1997-01-01
The present study adapts the Z-matrix formulation using a mixed basis of finite elements and Gaussians. This is a energy-independent basis which allows flexible boundary conditions and is amenable to efficient algorithms for evaluating the necessary matrix elements with molecular targets.
Destructive interference of E2 matrix elements in a triaxial rotor model
Allmond, James M; Wood, J. L.; Kulp, W. D.
2010-01-01
A triaxial rotor model with independent inertia and electric quadrupole tensors is applied to nuclei that have certain E2 matrix elements equal to zero. It is shown that such vanishing E2 matrix elements are explained by the model as a destructive interference effect. The example of 196Pt is considered.
Wilkins, Michael J; Wrighton, Kelly C; Nicora, Carrie D; Williams, Kenneth H; McCue, Lee Ann; Handley, Kim M; Miller, Chris S; Giloteaux, Ludovic; Montgomery, Alison P; Lovley, Derek R; Banfield, Jillian F; Long, Philip E; Lipton, Mary S
2013-01-01
While microbial activities in environmental systems play a key role in the utilization and cycling of essential elements and compounds, microbial activity and growth frequently fluctuates in response to environmental stimuli and perturbations. To investigate these fluctuations within a saturated aquifer system, we monitored a carbon-stimulated in situ Geobacter population while iron reduction was occurring, using 16S rRNA abundances and high-resolution tandem mass spectrometry proteome measurements. Following carbon amendment, 16S rRNA analysis of temporally separated samples revealed the rapid enrichment of Geobacter-like environmental strains with strong similarity to G. bemidjiensis. Tandem mass spectrometry proteomics measurements suggest high carbon flux through Geobacter respiratory pathways, and the synthesis of anapleurotic four carbon compounds from acetyl-CoA via pyruvate ferredoxin oxidoreductase activity. Across a 40-day period where Fe(III) reduction was occurring, fluctuations in protein expression reflected changes in anabolic versus catabolic reactions, with increased levels of biosynthesis occurring soon after acetate arrival in the aquifer. In addition, localized shifts in nutrient limitation were inferred based on expression of nitrogenase enzymes and phosphate uptake proteins. These temporal data offer the first example of differing microbial protein expression associated with changing geochemical conditions in a subsurface environment. PMID:23472107
Wilkins, Michael J.; Wrighton, Kelly C.; Nicora, Carrie D.; Williams, Kenneth H.; McCue, Lee Ann; Handley, Kim M.; Miller, Chris S.; Giloteaux, Ludovic; Montgomery, Alison P.; Lovley, Derek R.; Banfield, Jillian F.; Long, Philip E.; Lipton, Mary S.
2013-01-01
While microbial activities in environmental systems play a key role in the utilization and cycling of essential elements and compounds, microbial activity and growth frequently fluctuates in response to environmental stimuli and perturbations. To investigate these fluctuations within a saturated aquifer system, we monitored a carbon-stimulated in situ Geobacter population while iron reduction was occurring, using 16S rRNA abundances and high-resolution tandem mass spectrometry proteome measurements. Following carbon amendment, 16S rRNA analysis of temporally separated samples revealed the rapid enrichment of Geobacter-like environmental strains with strong similarity to G. bemidjiensis. Tandem mass spectrometry proteomics measurements suggest high carbon flux through Geobacter respiratory pathways, and the synthesis of anapleurotic four carbon compounds from acetyl-CoA via pyruvate ferredoxin oxidoreductase activity. Across a 40-day period where Fe(III) reduction was occurring, fluctuations in protein expression reflected changes in anabolic versus catabolic reactions, with increased levels of biosynthesis occurring soon after acetate arrival in the aquifer. In addition, localized shifts in nutrient limitation were inferred based on expression of nitrogenase enzymes and phosphate uptake proteins. These temporal data offer the first example of differing microbial protein expression associated with changing geochemical conditions in a subsurface environment. PMID:23472107
Wilkins, Michael J.; Wrighton, Kelly C.; Nicora, Carrie D.; Williams, Kenneth H.; McCue, Lee Ann; Handley, Kim M.; Miller, C. S.; Giloteaux, L.; Montgomery, A. P.; Lovley, Derek R.; Banfield, Jillian F.; Long, Philip E.; Lipton, Mary S.
2013-03-05
While microbial activities in environmental systems play a key role in the utilization and cycling of essential elements and compounds, microbial activity and growth frequently fluctuates in response to environmental stimuli and perturbations. To investigate these fluctuations within a saturated aquifer system, we monitored a carbon-stimulated in situ Geobacter population while iron reduction was occurring, using 16S rRNA abundances and high-resolution tandem mass spectrometry proteome measurements. Following carbon amendment, 16S rRNA analysis of temporally separated samples revealed the rapid enrichment of Geobacter-like environmental strains with strong similarity to G. bemidjiensis. Tandem mass spectrometry proteomics measurements suggest high carbon flux through Geobacter respiratory pathways, and the synthesis of anapleurotic four carbon compounds from acetyl-CoA via pyruvate ferredoxin oxidoreductase activity. Across a 40-day period where Fe(III) reduction was occurring, fluctuations in protein expression reflected changes in anabolic versus catabolic reactions, with increased levels of biosynthesis occurring soon after acetate arrival in the aquifer. In addition, localized shifts in nutrient limitation were inferred based on expression of nitrogenase enzymes and phosphate uptake proteins. These temporal data offer the first example of differing microbial protein expression associated with changing geochemical conditions in a subsurface environment.
EH3 matrix mineralogy with major and trace element composition compared to chondrules
NASA Astrophysics Data System (ADS)
Lehner, S. W.; McDonough, W. F.; NéMeth, P.
2014-12-01
We investigated the matrix mineralogy in primitive EH3 chondrites Sahara 97072, ALH 84170, and LAR 06252 with transmission electron microscopy; measured the trace and major element compositions of Sahara 97072 matrix and ferromagnesian chondrules with laser-ablation, inductively coupled, plasma mass spectrometry (LA-ICPMS); and analyzed the bulk composition of Sahara 97072 with LA-ICPMS, solution ICPMS, and inductively coupled plasma atomic emission spectroscopy. The fine-grained matrix of EH3 chondrites is unlike that in other chondrite groups, consisting primarily of enstatite, cristobalite, troilite, and kamacite with a notable absence of olivine. Matrix and pyroxene-rich chondrule compositions differ from one another and are distinct from the bulk meteorite. Refractory lithophile elements are enriched by a factor of 1.5-3 in chondrules relative to matrix, whereas the matrix is enriched in moderately volatile elements. The compositional relation between the chondrules and matrix is reminiscent of the difference between EH3 pyroxene-rich chondrules and EH3 Si-rich, highly sulfidized chondrules. Similar refractory element ratios between the matrix and the pyroxene-rich chondrules suggest the fine-grained material primarily consists of the shattered, sulfidized remains of the formerly pyroxene-rich chondrules with the minor addition of metal clasts. The matrix, chondrule, and metal-sulfide nodule compositions are probably complementary, suggesting all the components of the EH3 chondrites came from the same nebular reservoir.
NASA Astrophysics Data System (ADS)
Daehnick, W. W.
1983-07-01
Matrix elements for the effective two-nucleon interaction have been deduced from the population of multiplets near closed shells as observed in direct transfer reactions. In the evaluation, the limited purity of such multiples was taken into consideration, typically by weighting the observed fractions of the two-nucleon configurations by their spectroscopic strenghts and by using the resulting energy centroids. In a few cases, off-diagonal matrix elements are available from empirical wave funcitons. The systematic errors for particle-particle matrix elements extracted directly and those obtained from Pandya transformations were found to go in opposite directions. In some cases, this feautre of the empirical mehtod could be used to suggest upper and lower “bounds” for the extracted matrix elements. Diagonal matrix elements for the empirical residual interaction show a number of features suggestive of an underlying simplicity in the interaction of bound nucleons. Within experimental uncertainties (of about 10% for T=0 matrix elements) the monopole parts of the matrix elements are fit well with a simple A-0.75 dependence, and the data available to date do not reveal any significant monopole dependence on the quantum numbers of the interacting nucleons. The usefulness of scaling is suggested. Generally, diagonal matrix elements EJ( j1, j2) normalized by the extracted A-dependent monopole strength agree within expected experimental uncertainties whether derived from particle-particle or particle-hole multiples and whether extracted from the beginning or the end of a major shell. For values J≠0, the diagonal EJ( j2) matrix elements seem to follow two universal functions which depend on the semi-classical coupling angles θ 12, but are otherwise independent on j. For j1≠ j2 several “typical” functions ƒ(θ 12) can be constructed which fit subsets of the data and differ in a predictable way. The general features of the bound-nucleon interaction appear
Measurement of the CKM Matrix Elements |Vcb| and |Vub| at the B-factories
Menges, Wolfgang
2006-08-01
Recent results on inclusive and exclusive semileptonic B decays from B-factories are presented. The impact of these measurements on the determination of the CKM matrix elements |V{sub ub}| and |V{sub cb}| is discussed.
Can one measure nuclear matrix elements of neutrinoless double {beta} decay?
Rodin, Vadim; Faessler, Amand
2009-10-15
By making use of the isospin conservation by strong interaction, the Fermi 0{nu}{beta}{beta} nuclear matrix element M{sub F}{sup 0{nu}} is transformed to acquire the form of an energy-weighted double Fermi transition matrix element. This useful representation allows reconstruction of the total M{sub F}{sup 0{nu}} provided a small isospin-breaking Fermi matrix element between the isobaric analog state in the intermediate nucleus and the ground state of the daughter nucleus could be measured, e.g., by charge-exchange reactions. Such a measurement could set a scale for the 0{nu}{beta}{beta} nuclear matrix elements and help to discriminate between the different nuclear structure models in which calculated M{sub F}{sup 0{nu}} may differ by as much as a factor of 5 (that translates to about 20% difference in the total M{sup 0{nu}})
NASA Astrophysics Data System (ADS)
Moynier, F.; Pringle, E.; Hezel, D.
2015-07-01
The variations of Zn isotope ratios among carbonaceous chondrites show that the volatile element depletion in solar system material occurred in the solar nebula. We will also present the Zn isotopic composition of chondrules and matrix from carbonaceous chondrites.
Calculation of Radiative Corrections to E1 matrix elements in the Neutral Alkalis
Sapirstein, J; Cheng, K T
2004-09-28
Radiative corrections to E1 matrix elements for ns-np transitions in the alkali metal atoms lithium through francium are evaluated. They are found to be small for the lighter alkalis but significantly larger for the heavier alkalis, and in the case of cesium much larger than the experimental accuracy. The relation of the matrix element calculation to a recent decay rate calculation for hydrogenic ions is discussed, and application of the method to parity nonconservation in cesium is described.
Search for Off-Diagonal Density Matrix Elements for Atoms in a Supersonic Beam
NASA Astrophysics Data System (ADS)
Rubenstein, Richard A.; Dhirani, Al-Amin; Kokorowski, David A.; Roberts, Tony D.; Smith, Edward T.; Smith, Winthrop W.; Bernstein, Herbert J.; Lehner, Jana; Gupta, Subhadeep; Pritchard, David E.
1999-03-01
We demonstrate the absence of off-diagonal elements for the density matrix of a supersonic Na atomic beam, thus showing that there are no coherent wave packets emerging from this source. We used a differentially detuned separated oscillatory field longitudinal interferometer to search for off-diagonal density matrix elements in the longitudinal energy/momentum basis. Our study places a stringent lower bound on their possible size over an off-diagonal energy range from 0 to 100 kHz.
Elemental: a new framework for distributed memory dense matrix computations.
Romero, N.; Poulson, J.; Marker, B.; Hammond, J.; Van de Geijn, R.
2012-02-14
Parallelizing dense matrix computations to distributed memory architectures is a well-studied subject and generally considered to be among the best understood domains of parallel computing. Two packages, developed in the mid 1990s, still enjoy regular use: ScaLAPACK and PLAPACK. With the advent of many-core architectures, which may very well take the shape of distributed memory architectures within a single processor, these packages must be revisited since the traditional MPI-based approaches will likely need to be extended. Thus, this is a good time to review lessons learned since the introduction of these two packages and to propose a simple yet effective alternative. Preliminary performance results show the new solution achieves competitive, if not superior, performance on large clusters.
Matrix elements and duality for type 2 unitary representations of the Lie superalgebra gl(m|n)
Werry, Jason L.; Gould, Mark D.; Isaac, Phillip S.
2015-12-15
The characteristic identity formalism discussed in our recent articles is further utilized to derive matrix elements of type 2 unitary irreducible gl(m|n) modules. In particular, we give matrix element formulae for all gl(m|n) generators, including the non-elementary generators, together with their phases on finite dimensional type 2 unitary irreducible representations which include the contravariant tensor representations and an additional class of essentially typical representations. Remarkably, we find that the type 2 unitary matrix element equations coincide with the type 1 unitary matrix element equations for non-vanishing matrix elements up to a phase.
Matrix elements and duality for type 2 unitary representations of the Lie superalgebra gl(m|n)
NASA Astrophysics Data System (ADS)
Werry, Jason L.; Gould, Mark D.; Isaac, Phillip S.
2015-12-01
The characteristic identity formalism discussed in our recent articles is further utilized to derive matrix elements of type 2 unitary irreducible gl(m|n) modules. In particular, we give matrix element formulae for all gl(m|n) generators, including the non-elementary generators, together with their phases on finite dimensional type 2 unitary irreducible representations which include the contravariant tensor representations and an additional class of essentially typical representations. Remarkably, we find that the type 2 unitary matrix element equations coincide with the type 1 unitary matrix element equations for non-vanishing matrix elements up to a phase.
Use of R-matrix theory in light element evaluations
Hale, G.M.
1992-12-31
R-matrix theory is a general framework for describing nuclear reactions (neutron- and charged-particle-induced) that is particularly well-suited for including resonances. We will review some unconventional properties of resonances that arise from this theory, including non-Breit-Wigner (BW) resonances and shadow poles, and discuss their physical consequences. Examples will be given from the analyses of several light systems that have been used in ENDF evaluations, including {sup 4}He, {sup 5}He, {sup 15}N, and {sup 17}O. The resonances in the helium systems tend to be significantly non-BW in character, while almost all the resonances in {sup 15}N and {sup 17}O are Breit-Wigner. An interesting exception in the case of {sup 15}N indicates that some of the sub-threshold levels that have been assumed to be bound might be virtual. We find that fitting data from all possible reactions simultaneously results in level schemes for the compound systems that differ in some cases significantly from the ``accepted`` published level information.
Matrix exponentials, SU(N) group elements, and real polynomial roots
NASA Astrophysics Data System (ADS)
Van Kortryk, T. S.
2016-02-01
The exponential of an N × N matrix can always be expressed as a matrix polynomial of order N - 1. In particular, a general group element for the fundamental representation of SU(N) can be expressed as a matrix polynomial of order N - 1 in a traceless N × N hermitian generating matrix, with polynomial coefficients consisting of elementary trigonometric functions dependent on N - 2 invariants in addition to the group parameter. These invariants are just angles determined by the direction of a real N-vector whose components are the eigenvalues of the hermitian matrix. Equivalently, the eigenvalues are given by projecting the vertices of an (" separators=" N - 1 ) -simplex onto a particular axis passing through the center of the simplex. The orientation of the simplex relative to this axis determines the angular invariants and hence the real eigenvalues of the matrix.
A Matrix Model for Reliability of a Cold-Standby system with Identical Repairable Elements
NASA Astrophysics Data System (ADS)
Farahpour, Peyman; Mahshid, Kamrouz; Sharifi, Mani; Palizban, Aidin
2011-09-01
In this paper we studied a cold standby system with n identical constant failure rate repairable elements. The system has m repairmen and each repairman only works on the one failed element. After failings one element, another element replace immediately. The failure and repair rate of each element is constant as λ, μ. At first a matrix model presented to determine the state of the system. Then we establish the differential equations between the states of the system and finally with a numerical example, we illustrate the method of solving the equations. This paper divided to five main parts, we present some studies about the redundancy allocation and the marcovian models in the introduction. In the second part introduce the system description. In the third part differential equations of the system have been presented in a matrix. A numerical example presented in the 4th part to illustrated how to work with these equations. Last parts we deal with conclusion and future studies.
Efficient computation of Hamiltonian matrix elements between non-orthogonal Slater determinants
NASA Astrophysics Data System (ADS)
Utsuno, Yutaka; Shimizu, Noritaka; Otsuka, Takaharu; Abe, Takashi
2013-01-01
We present an efficient numerical method for computing Hamiltonian matrix elements between non-orthogonal Slater determinants, focusing on the most time-consuming component of the calculation that involves a sparse array. In the usual case where many matrix elements should be calculated, this computation can be transformed into a multiplication of dense matrices. It is demonstrated that the present method based on the matrix-matrix multiplication attains ˜80% of the theoretical peak performance measured on systems equipped with modern microprocessors, a factor of 5-10 better than the normal method using indirectly indexed arrays to treat a sparse array. The reason for such different performances is discussed from the viewpoint of memory access.
Quenched domain wall QCD with DBW2 gauge action toward nucleon decay matrix element calculation
NASA Astrophysics Data System (ADS)
Aoki, Yasumichi
2001-10-01
The domain wall fermion action is a promising way to control chiral symmetry in lattice gauge theory. By the good chiral symmetry of this approach even at finite lattice spacing, one is able to extract hadronic matrix elements, like kaon weak matrix elements, for which the symmetry is extremely important. Ordinary fermions with poor chiral symmetry make calculation difficult because of the large mixing of operators with different chiral structure. Even though the domain wall fermion action with the simple Wilson gauge action has a good chiral symmetry, one can further improve the symmetry by using a different gauge action. We take a non-perturbatively improved action, the DBW2 action of the QCD Taro group. Hadron masses are systematically examined for a range of parameters. Application to nucleon decay matrix element is also discussed.
Uncertainties in nuclear transition matrix elements of neutrinoless ββ decay
Rath, P. K.
2013-12-30
To estimate the uncertainties associated with the nuclear transition matrix elements M{sup (K)} (K=0ν/0N) for the 0{sup +} → 0{sup +} transitions of electron and positron emitting modes of the neutrinoless ββ decay, a statistical analysis has been performed by calculating sets of eight (twelve) different nuclear transition matrix elements M{sup (K)} in the PHFB model by employing four different parameterizations of a Hamiltonian with pairing plus multipolar effective two-body interaction and two (three) different parameterizations of Jastrow short range correlations. The averages in conjunction with their standard deviations provide an estimate of the uncertainties associated the nuclear transition matrix elements M{sup (K)} calculated within the PHFB model, the maximum of which turn out to be 13% and 19% owing to the exchange of light and heavy Majorana neutrinos, respectively.
B(s) 0-mixing matrix elements from lattice QCD for the Standard Model and beyond
NASA Astrophysics Data System (ADS)
Bazavov, A.; Bernard, C.; Bouchard, C. M.; Chang, C. C.; DeTar, C.; Du, Daping; El-Khadra, A. X.; Freeland, E. D.; Gámiz, E.; Gottlieb, Steven; Heller, U. M.; Kronfeld, A. S.; Laiho, J.; Mackenzie, P. B.; Neil, E. T.; Simone, J.; Sugar, R.; Toussaint, D.; Van de Water, R. S.; Zhou, Ran; Fermilab Lattice; MILC Collaborations
2016-06-01
We calculate—for the first time in three-flavor lattice QCD—the hadronic matrix elements of all five local operators that contribute to neutral B0- and Bs-meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral B -meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio ξ =1.206 (18 )(6 ), where the second error stems from the omission of charm-sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from B mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light-valence quarks and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral B -meson oscillation frequencies to determine the CKM matrix elements |Vt d|=8.00 (34 )(8 )×10-3, |Vt s|=39.0 (1.2 )(0.4 )×10-3, and |Vt d/Vt s|=0.2052 (31 )(10 ), which differ from CKM-unitarity expectations by about 2 σ . These results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels.
Calculation of radiative corrections to E1 matrix elements in the neutral alkali metals
Sapirstein, J.; Cheng, K.T.
2005-02-01
Radiative corrections to E1 matrix elements for ns-np transitions in the alkali-metal atoms lithium through francium are evaluated. They are found to be small for the lighter alkali metals but significantly larger for the heavier alkali metals, and in the case of cesium much larger than the experimental accuracy. The relation of the matrix element calculation to a recent decay rate calculation for hydrogenic ions is discussed, and application of the method to parity nonconservation in cesium is described.
Matrix elements for type 1 unitary irreducible representations of the Lie superalgebra gl(m|n)
Gould, Mark D.; Isaac, Phillip S.; Werry, Jason L.
2014-01-15
Using our recent results on eigenvalues of invariants associated to the Lie superalgebra gl(m|n), we use characteristic identities to derive explicit matrix element formulae for all gl(m|n) generators, particularly non-elementary generators, on finite dimensional type 1 unitary irreducible representations. We compare our results with existing works that deal with only subsets of the class of type 1 unitary representations, all of which only present explicit matrix elements for elementary generators. Our work therefore provides an important extension to existing methods, and thus highlights the strength of our techniques which exploit the characteristic identities.
Quantum Monte Carlo calculations of electroweak transition matrix elements in A=6,7 nuclei
Pervin, Muslema; Pieper, Steven C.; Wiringa, R. B.
2007-12-15
Green's function Monte Carlo (GFMC) calculations of magnetic dipole, electric quadrupole, Fermi, and Gamow-Teller transition matrix elements are reported for A=6,7 nuclei. The matrix elements are extrapolated from mixed estimates that bracket the relevant electroweak operator between variational Monte Carlo (VMC) and GFMC propagated wave functions. Because they are off-diagonal terms, two mixed estimates are required for each transition, with a VMC initial (final) state paired with a GFMC final (initial) state. The realistic Argonne v{sub 18} two-nucleon and Illinois-2 three-nucleon interactions are used to generate the nuclear states. In most cases we find good agreement with experimental data.
Analytic matrix elements for the two-electron atomic basis with logarithmic terms
Liverts, Evgeny Z.; Barnea, Nir
2014-08-01
The two-electron problem for the helium-like atoms in S-state is considered. The basis containing the integer powers of ln r, where r is a radial variable of the Fock expansion, is studied. In this basis, the analytic expressions for the matrix elements of the corresponding Hamiltonian are presented. These expressions include only elementary and special functions, what enables very fast and accurate computation of the matrix elements. The decisive contribution of the correct logarithmic terms to the behavior of the two-electron wave function in the vicinity of the triple-coalescence point is reaffirmed.
LIBS Detection of Heavy Metal Elements in Liquid Solutions by Using Wood Pellet as Sample Matrix
NASA Astrophysics Data System (ADS)
Wen, Guanhong; Sun, Duixiong; Su, Maogen; Dong, Chenzhong
2014-06-01
Laser-induced breakdown spectroscopy (LIBS) has been applied to the analysis of heavy metals in liquid samples. A new approach was presented to lower the limit of detection (LOD) and minimize the sample matrix effects, in which dried wood pellets absorbed the given amounts of Cr standard solutions and then were baked because they have stronger and rapid absorption properties for liquid samples as well as simple elemental compositions. In this work, we have taken a typical heavy metal Cr element as an example, and investigated the spectral feasibility of Cr solutions and dried wood pellets before and after absorbing Cr solutions at the same experimental conditions. The results were demonstrated to successfully produce a superior analytical response for heavy metal elements by using wood pellet as sample matrix according to the obtained LOD of 0.07 ppm for Cr element in solutions.
NASA Astrophysics Data System (ADS)
Kumar, Deepak; Roy, Rene; Kweon, Jin-Hwe; Choi, Jin-ho
2015-10-01
Sub-laminate damage in the form of matrix cracking and delamination was simulated by using interface cohesive elements in the finite element (FE) software ABAQUS. Interface cohesive elements were inserted parallel to the fiber orientation in the transverse ply with equal spacing (matrix cracking) and between the interfaces (delamination). Matrix cracking initiation in the cohesive elements was based on stress traction separation laws and propagated under mixed-mode loading. We expanded the work of Shi et al. (Appl. Compos. Mater. 21, 57-70 2014) to include delamination and simulated additional [45/-45/0/90]s and [02/90n]s {n = 1,2,3} CFRP laminates and a [0/903]s GFRP laminate. Delamination damage was quantified numerically in terms of damage dissipative energy. We observed that transverse matrix cracks can propagate to the ply interface and initiate delamination. We also observed for [0/90n/0] laminates that as the number of 90° ply increases past n = 2, the crack density decreases. The predicted crack density evolution compared well with experimental results and the equivalent constraint model (ECM) theory. Empirical relationships were established between crack density and applied stress by linear curve fitting. The reduction of laminate elastic modulus due to cracking was also computed numerically and it is in accordance with reported experimental measurements.
NASA Astrophysics Data System (ADS)
Kumar, Deepak; Roy, Rene; Kweon, Jin-Hwe; Choi, Jin-ho
2016-06-01
Sub-laminate damage in the form of matrix cracking and delamination was simulated by using interface cohesive elements in the finite element (FE) software ABAQUS. Interface cohesive elements were inserted parallel to the fiber orientation in the transverse ply with equal spacing (matrix cracking) and between the interfaces (delamination). Matrix cracking initiation in the cohesive elements was based on stress traction separation laws and propagated under mixed-mode loading. We expanded the work of Shi et al. (Appl. Compos. Mater. 21, 57-70 2014) to include delamination and simulated additional [45/-45/0/90]s and [02/90n]s { n = 1,2,3} CFRP laminates and a [0/903]s GFRP laminate. Delamination damage was quantified numerically in terms of damage dissipative energy. We observed that transverse matrix cracks can propagate to the ply interface and initiate delamination. We also observed for [0/90n/0] laminates that as the number of 90° ply increases past n = 2, the crack density decreases. The predicted crack density evolution compared well with experimental results and the equivalent constraint model (ECM) theory. Empirical relationships were established between crack density and applied stress by linear curve fitting. The reduction of laminate elastic modulus due to cracking was also computed numerically and it is in accordance with reported experimental measurements.
Finite-Element Z-Matrix Calculation of Electron-N2 Collisions
NASA Technical Reports Server (NTRS)
Huo, Winifred M.; Dateo, Christopher E.
1999-01-01
The finite element Z-matrix method has been applied in a multichannel study of e-N2 Collisions for electron energies from threshold to 30 eV. General agreement is obtained comparing with existing experimental and theoretical data. Some discrepancies are also found.
Matrix elements of explicitly correlated Gaussian basis functions with arbitrary angular momentum
NASA Astrophysics Data System (ADS)
Joyce, Tennesse; Varga, Kálmán
2016-05-01
A new algorithm for calculating the Hamiltonian matrix elements with all-electron explicitly correlated Gaussian functions for quantum-mechanical calculations of atoms with arbitrary angular momentum is presented. The calculations are checked on several excited states of three and four electron systems. The presented formalism can be used as unified framework for high accuracy calculations of properties of small atoms and molecules.
Matrix-Assisted Plasma Atomization Emission Spectrometry for Surface Sampling Elemental Analysis
Yuan, Xin; Zhan, Xuefang; Li, Xuemei; Zhao, Zhongjun; Duan, Yixiang
2016-01-01
An innovative technology has been developed involving a simple and sensitive optical spectrometric method termed matrix-assisted plasma atomization emission spectrometry (MAPAES) for surface sampling elemental analysis using a piece of filter paper (FP) for sample introduction. MAPAES was carried out by direct interaction of the plasma tail plume with the matrix surface. The FP absorbs energy from the plasma source and releases combustion heating to the analytes originally present on its surface, thus to promote the atomization and excitation process. The matrix-assisted plasma atomization excitation phenomenon was observed for multiple elements. The FP matrix served as the partial energy producer and also the sample substrate to adsorb sample solution. Qualitative and quantitative determinations of metal ions were achieved by atomic emission measurements for elements Ba, Cu, Eu, In, Mn, Ni, Rh and Y. The detection limits were down to pg level with linear correlation coefficients better than 0.99. The proposed MAPAES provides a new way for atomic spectrometry which offers advantages of fast analysis speed, little sample consumption, less sample pretreatment, small size, and cost-effective. PMID:26762972
Matrix-Assisted Plasma Atomization Emission Spectrometry for Surface Sampling Elemental Analysis
NASA Astrophysics Data System (ADS)
Yuan, Xin; Zhan, Xuefang; Li, Xuemei; Zhao, Zhongjun; Duan, Yixiang
2016-01-01
An innovative technology has been developed involving a simple and sensitive optical spectrometric method termed matrix-assisted plasma atomization emission spectrometry (MAPAES) for surface sampling elemental analysis using a piece of filter paper (FP) for sample introduction. MAPAES was carried out by direct interaction of the plasma tail plume with the matrix surface. The FP absorbs energy from the plasma source and releases combustion heating to the analytes originally present on its surface, thus to promote the atomization and excitation process. The matrix-assisted plasma atomization excitation phenomenon was observed for multiple elements. The FP matrix served as the partial energy producer and also the sample substrate to adsorb sample solution. Qualitative and quantitative determinations of metal ions were achieved by atomic emission measurements for elements Ba, Cu, Eu, In, Mn, Ni, Rh and Y. The detection limits were down to pg level with linear correlation coefficients better than 0.99. The proposed MAPAES provides a new way for atomic spectrometry which offers advantages of fast analysis speed, little sample consumption, less sample pretreatment, small size, and cost-effective.
Matrix-Assisted Plasma Atomization Emission Spectrometry for Surface Sampling Elemental Analysis.
Yuan, Xin; Zhan, Xuefang; Li, Xuemei; Zhao, Zhongjun; Duan, Yixiang
2016-01-01
An innovative technology has been developed involving a simple and sensitive optical spectrometric method termed matrix-assisted plasma atomization emission spectrometry (MAPAES) for surface sampling elemental analysis using a piece of filter paper (FP) for sample introduction. MAPAES was carried out by direct interaction of the plasma tail plume with the matrix surface. The FP absorbs energy from the plasma source and releases combustion heating to the analytes originally present on its surface, thus to promote the atomization and excitation process. The matrix-assisted plasma atomization excitation phenomenon was observed for multiple elements. The FP matrix served as the partial energy producer and also the sample substrate to adsorb sample solution. Qualitative and quantitative determinations of metal ions were achieved by atomic emission measurements for elements Ba, Cu, Eu, In, Mn, Ni, Rh and Y. The detection limits were down to pg level with linear correlation coefficients better than 0.99. The proposed MAPAES provides a new way for atomic spectrometry which offers advantages of fast analysis speed, little sample consumption, less sample pretreatment, small size, and cost-effective. PMID:26762972
$$B^0_{(s)}$$-mixing matrix elements from lattice QCD for the Standard Model and beyond
Bazavov, A.; Bernard, C.; Bouchard, C. M.; Chang, C. C.; DeTar, C.; Du, Daping; El-Khadra, A. X.; Freeland, E. D.; Gamiz, E.; Gottlieb, Steven; et al
2016-06-28
We calculate—for the first time in three-flavor lattice QCD—the hadronic matrix elements of all five local operators that contribute to neutral B0- and Bs-meson mixing in and beyond the Standard Model. We present a complete error budget for each matrix element and also provide the full set of correlations among the matrix elements. We also present the corresponding bag parameters and their correlations, as well as specific combinations of the mixing matrix elements that enter the expression for the neutral B-meson width difference. We obtain the most precise determination to date of the SU(3)-breaking ratio ξ=1.206(18)(6), where the second errormore » stems from the omission of charm-sea quarks, while the first encompasses all other uncertainties. The threefold reduction in total uncertainty, relative to the 2013 Flavor Lattice Averaging Group results, tightens the constraint from B mixing on the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle. Our calculation employs gauge-field ensembles generated by the MILC Collaboration with four lattice spacings and pion masses close to the physical value. We use the asqtad-improved staggered action for the light-valence quarks and the Fermilab method for the bottom quark. We use heavy-light meson chiral perturbation theory modified to include lattice-spacing effects to extrapolate the five matrix elements to the physical point. We combine our results with experimental measurements of the neutral B-meson oscillation frequencies to determine the CKM matrix elements |Vtd| = 8.00(34)(8)×10-3, |Vts| = 39.0(1.2)(0.4)×10-3, and |Vtd/Vts| = 0.2052(31)(10), which differ from CKM-unitarity expectations by about 2σ. In addition, these results and others from flavor-changing-neutral currents point towards an emerging tension between weak processes that are mediated at the loop and tree levels.« less
Bubin, Sergiy; Adamowicz, Ludwik
2006-06-14
In this work we present analytical expressions for Hamiltonian matrix elements with spherically symmetric, explicitly correlated Gaussian basis functions with complex exponential parameters for an arbitrary number of particles. The expressions are derived using the formalism of matrix differential calculus. In addition, we present expressions for the energy gradient that includes derivatives of the Hamiltonian integrals with respect to the exponential parameters. The gradient is used in the variational optimization of the parameters. All the expressions are presented in the matrix form suitable for both numerical implementation and theoretical analysis. The energy and gradient formulas have been programmed and used to calculate ground and excited states of the He atom using an approach that does not involve the Born-Oppenheimer approximation. PMID:16784284
NASA Technical Reports Server (NTRS)
Caruso, J. J.; Trowbridge, D.; Chamis, C. C.
1989-01-01
The mechanics of materials approach (definition of E, G, Nu, and Alpha) and the finite element method are used to explore the effects of partial bonding and fiber fracture on the behavior of high temperature metal matrix composites. Composite ply properties are calculated for various degrees of disbonding to evaluate the sensitivity of these properties to the presence of fiber/matrix disbonding and fiber fracture. The mechanics of materials approach allows for the determination of the basic ply material properties needed for design/analysis of composites. The finite element method provides the necessary structural response (forces and displacements) for the mechanics of materials equations. Results show that disbonding of fractured fibers affect only E sub (111) and alpha sub (111) significantly.
NASA Technical Reports Server (NTRS)
Caruso, J. J.; Chamis, C. C.; Trowbridge, D.
1989-01-01
The mechanics of materials approach (definition of E, G, nu, and alpha) and the finite element method are used to explore the effects of partial bonding and fiber fracture on the behavior of high temperature metal matrix composites. Composite ply properties are calculated for various degrees of disbonding to evaluate the sensitivity of these properties to the presence of fiber/matrix disbonding and fiber fracture. The mechanics of materials approach allows for the determination of the basic ply material properties needed for design/analysis of composites. The finite element method provides the necessary structural response (forces and displacements) for the mechanics of materials equations. Results show that disbonding of fractured fibers affect only E-l(11) and alpha-l(11) significantly.
Kinetic-energy matrix elements for atomic Hylleraas-CI wave functions.
Harris, Frank E
2016-05-28
Hylleraas-CI is a superposition-of-configurations method in which each configuration is constructed from a Slater-type orbital (STO) product to which is appended (linearly) at most one interelectron distance rij. Computations of the kinetic energy for atoms by this method have been difficult due to the lack of formulas expressing these matrix elements for general angular momentum in terms of overlap and potential-energy integrals. It is shown here that a strategic application of angular-momentum theory, including the use of vector spherical harmonics, enables the reduction of all atomic kinetic-energy integrals to overlap and potential-energy matrix elements. The new formulas are validated by showing that they yield correct results for a large number of integrals published by other investigators. PMID:27250282
Kinetic-energy matrix elements for atomic Hylleraas-CI wave functions
NASA Astrophysics Data System (ADS)
Harris, Frank E.
2016-05-01
Hylleraas-CI is a superposition-of-configurations method in which each configuration is constructed from a Slater-type orbital (STO) product to which is appended (linearly) at most one interelectron distance rij. Computations of the kinetic energy for atoms by this method have been difficult due to the lack of formulas expressing these matrix elements for general angular momentum in terms of overlap and potential-energy integrals. It is shown here that a strategic application of angular-momentum theory, including the use of vector spherical harmonics, enables the reduction of all atomic kinetic-energy integrals to overlap and potential-energy matrix elements. The new formulas are validated by showing that they yield correct results for a large number of integrals published by other investigators.
T -dual Ramond-Ramond couplings on D-branes from S -matrix elements
NASA Astrophysics Data System (ADS)
Babaei Velni, Komeil
2016-03-01
Using the linear T -dual Ward identity associated with the NS-NS B-field or two Neveu-Schwarz B-fields (NSNS) gauge transformations, some Ramond-Ramond (RR) couplings on Dp -branes have been found at order O (α'2) . We examine the C(p -1 ) couplings with the S -matrix elements of one RR, one graviton and one antisymmetric B-field vertex operator. We find the consistency of T -dual S -matrix elements and explicit results of scattering string amplitude and show that the string amplitude reproduces these couplings as well as some other couplings. This illustration is found for C(p -3 ) couplings in the literature, which is extended to the C(p -1 ) couplings in this paper.
Short-distance matrix elements for $D$-meson mixing for 2+1 lattice QCD
Chang, Chia Cheng
2015-01-01
We study the short-distance hadronic matrix elements for D-meson mixing with partially quenched N_{f} = 2+1 lattice QCD. We use a large set of the MIMD Lattice Computation Collaboration's gauge configurations with a^{2} tadpole-improved staggered sea quarks and tadpole-improved Lüscher-Weisz gluons. We use the a^{2} tadpole-improved action for valence light quarks and the Sheikoleslami-Wohlert action with the Fermilab interpretation for the valence charm quark. Our calculation covers the complete set of five operators needed to constrain new physics models for D-meson mixing. We match our matrix elements to the MS-NDR scheme evaluated at 3 GeV. We report values for the Beneke-Buchalla-Greub-Lenz-Nierste choice of evanescent operators.
Fluctuation phenomena in layered superconductors
Klemm, R.A.
1996-10-01
Gaussian fluctuations in layered superconductors have been the subject of study for many years. Although the FD was studied in detail long ago, the FC (fluctuation conductivity) was studied only recently, since the MT and DOS diagrams were previously neglected. Recent comparisons with experiment on YBCO have shown that the DOS diagrams are important and can lead to qualitatively different behaviors for the FC parallel and perpendicular to the layers. In both cases, Gaussian fluctuations fit the data above {Tc} very well, even for YBCO. To date, nearly all calculations of fluctuation quantities were for B{parallel}{cflx c}. Nevertheless, it should be possible to treat an arbitrary B, but the evaluation of the required matrix elements for the fluctuation quantities will be more complicated.
NASA Technical Reports Server (NTRS)
Buehler, Martin G. (Inventor)
1988-01-01
A set of addressable test structures, each of which uses addressing schemes to access individual elements of the structure in a matrix, is used to test the quality of a wafer before integrated circuits produced thereon are diced, packaged and subjected to final testing. The electrical characteristic of each element is checked and compared to the electrical characteristic of all other like elements in the matrix. The effectiveness of the addressable test matrix is in readily analyzing the electrical characteristics of the test elements and in providing diagnostic information.
Finite-element grid improvement by minimization of stiffness matrix trace
NASA Technical Reports Server (NTRS)
Kittur, Madan G.; Huston, Ronald L.; Oswald, Fred B.
1987-01-01
A new and simple method of finite-element grid improvement is presented. The objective is to improve the accuracy of the analysis. The procedure is based on a minimization of the trace of the stiffness matrix. For a broad class of problems this minimization is seen to be equivalent to minimizing the potential energy. The method is illustrated with the classical tapered bar problem examined earlier by Prager and Masur. Identical results are obtained.
Finite-element grid improvement by minimization of stiffness matrix trace
NASA Technical Reports Server (NTRS)
Kittur, Madan G.; Huston, Ronald L.; Oswald, Fred B.
1989-01-01
A new and simple method of finite-element grid improvement is presented. The objective is to improve the accuracy of the analysis. The procedure is based on a minimization of the trace of the stiffness matrix. For a broad class of problems this minimization is seen to be equivalent to minimizing the potential energy. The method is illustrated with the classical tapered bar problem examined earlier by Prager and Masur. Identical results are obtained.
An improved method for extracting matrix elements from lattice three-point functions
C. Aubin, K. Orginos
2011-12-01
The extraction of matrix elements from baryon three-point functions is complicated by the fact that the signal-to-noise drops rapidly as a function of time. Using a previously discussed method to improve the signal-to-noise for lattice two-point functions, we use this technique to do so for lattice three-point functions, using electromagnetic form factors for the nucleon and Delta as an example.
Useful extremum principle for the variational calculation of matrix elements. II
NASA Technical Reports Server (NTRS)
Gerjuoy, E.; Rosenberg, L.; Spruch, L.
1975-01-01
Recent work (Gerjuoy et al., 1974) on variational principles for diagonal bound state matrix elements of arbitrary Hermitian operators is extended. In particular, it is shown that the previously derived minimum principle for the trial auxiliary function appearing in such variational principles can be constructed using a modified Hamiltonian possessing not heretofore recognized positive definite properties. Thus there is at least one alternative to the particular modified Hamiltonian on which the results of Gerjuoy et al. (1974) originally were based.
Closed expressions for matrix elements of the trigonometric Pöschl-Teller potential
NASA Astrophysics Data System (ADS)
Rey, Michaël; Michelot, Françoise
2010-10-01
Analytical matrix elements of the x ( n>0) and [[d/dx operators are derived using the eigenfunctions of the symmetric trigonometric Pöschl-Teller potential. The closed formulas are written in terms of Gauss hypergeometric functions and could be used in variational calculations to describe vibrational energy levels associated with bending modes. Multiprecision computational packages are considered in order to obtain an arbitrary level of precision.
D'Ariano, G M; Lo Presti, P
2001-05-01
Quantum operations describe any state change allowed in quantum mechanics, including the evolution of an open system or the state change due to a measurement. We present a general method based on quantum tomography for measuring experimentally the matrix elements of an arbitrary quantum operation. As input the method needs only a single entangled state. The feasibility of the technique for the electromagnetic field is shown, and the experimental setup is illustrated based on homodyne tomography of a twin beam. PMID:11328133
Nuclear Matrix Elements for two-neutrino DBD in Te isotopes
Bes, D. R.; Civitarese, O.
2009-11-09
Theoretical matrix elements, for the ground-state to ground-state two-neutrino double-beta-decay mode of {sup 128,130}Te isotopes, are calculated within a formalism which describes interactions between neutrons in a superfluid phase and protons in a normal phase. The model is basically a parameter-free one, since all relevant parameters are fixed from phenomenology. A comparison with the available experimental data is presented.
Dipole transition-matrix elements of the one-electron heterodiatomic quasimolecules
Devdariani, A.; Kereselidze, T.M.; Noselidze, I. L.; Dalimier, E.; Angelo, P.; Schott, R.; Sauvan, P.
2005-02-01
The problem of dipole transition-matrix element calculation for optical transitions in multiply charged one-electron diatomic quasimolecules with unequal nuclear charges Z{sub 1} and Z{sub 2} has been stated and solved. The quasimolecule Z{sub 1}eZ{sub 2} is a unique example of a two-center system for which the energy terms and dipole transition moments have been calculated precisely in the frame of a nonrelativistic approach. Particular examples for the optical transitions with Z{sub 1}=1.5,2,2.5,3 and Z{sub 2}=1 and with the principal quantum number of the united ion n{sub u}=1,2,3,4 have been tabulated. The scaling rules make it possible to determine the matrix elements for quasimolecules having nuclear charge ratios such as 2:1, 3:1, 3:2, and 5:2. Zeros at intermediate R and zero limiting values at large R are the highlighted features of the matrix elements. The heteronucleus case generates a large number of asymptotically forbidden transitions corresponding to transitions of an electron from one ion to another.
MOON for neutrino-less ββ decays and ββ nuclear matrix elements
NASA Astrophysics Data System (ADS)
Ejiri, H.
2009-11-01
The MOON project aims at spectroscopic 0vββ studies with the v-mass sensitivity of 100-30 meV by measuring two beta rays from 100Mo and/or 82Se. The detector is a compact super-module of multi-layer PL scintillator plates. R&D works made by the pro to-type MOON-1 and the small PL plate show the possible energy resolution of around σ~2.2%, as required for the mass sensitivity. Nuclear matrix elements M2v for 2vββ are shown to be given by the sum ΣLMk of the 2vββ matrix elements Mk through intermediate quasi-particle states in the Fermi-surface, where Mi is obtained experimentally by using the GT(Jπ = 1+) matrix elements of Mi(k) and Mf(k) for the successive single-β transitions through the k-th intermediate state.
Matrix elements and few-body calculations within the unitary correlation operator method
Roth, R.; Hergert, H.; Papakonstantinou, P.; Neff, T.; Feldmeier, H.
2005-09-01
We employ the unitary correlation operator method (UCOM) to construct correlated, low-momentum matrix elements of realistic nucleon-nucleon interactions. The dominant short-range central and tensor correlations induced by the interaction are included explicitly by an unitary transformation. Using correlated momentum-space matrix elements of the Argonne V18 potential, we show that the unitary transformation eliminates the strong off-diagonal contributions caused by the short-range repulsion and the tensor interaction and leaves a correlated interaction dominated by low-momentum contributions. We use correlated harmonic oscillator matrix elements as input for no-core shell model calculations for few-nucleon systems. Compared to the bare interaction, the convergence properties are dramatically improved. The bulk of the binding energy can already be obtained in very small model spaces or even with a single Slater determinant. Residual long-range correlations, not treated explicitly by the unitary transformation, can easily be described in model spaces of moderate size allowing for fast convergence. By varying the range of the tensor correlator we are able to map out the Tjon line and can in turn constrain the optimal correlator ranges.
Urinary stones as a novel matrix for human biomonitoring of toxic and essential elements.
Kuta, J; Smetanová, S; Benová, D; Kořistková, T; Machát, J
2016-02-01
Monitoring of body burden of toxic elements is usually based on analysis of concentration of particular elements in blood, urine and/or hair. Analysis of these matrices, however, predominantly reflects short- or medium-term exposure to trace elements or pollutants. In this work, urinary stones were investigated as a matrix for monitoring long-term exposure to toxic and essential elements. A total of 431 samples of urinary calculi were subjected to mineralogical and elemental analysis by infrared spectroscopy and inductively coupled plasma mass spectrometry. The effect of mineralogical composition of the stones and other parameters such as sex, age and geographical location on contents of trace and minor elements is presented. Our results demonstrate the applicability of such approach and confirm that the analysis of urinary calculi can be helpful in providing complementary information on human exposure to trace metals and their excretion. Analysis of whewellite stones (calcium oxalate monohydrate) with content of phosphorus <0.6 % has been proved to be a promising tool for biomonitoring of trace and minor elements. PMID:25736734
A Data Matrix Method for Improving the Quantification of Element Percentages of SEM/EDX Analysis
NASA Technical Reports Server (NTRS)
Lane, John
2009-01-01
A simple 2D M N matrix involving sample preparation enables the microanalyst to peer below the noise floor of element percentages reported by the SEM/EDX (scanning electron microscopy/ energy dispersive x-ray) analysis, thus yielding more meaningful data. Using the example of a 2 3 sample set, there are M = 2 concentration levels of the original mix under test: 10 percent ilmenite (90 percent silica) and 20 percent ilmenite (80 percent silica). For each of these M samples, N = 3 separate SEM/EDX samples were drawn. In this test, ilmenite is the element of interest. By plotting the linear trend of the M sample s known concentration versus the average of the N samples, a much higher resolution of elemental analysis can be performed. The resulting trend also shows how the noise is affecting the data, and at what point (of smaller concentrations) is it impractical to try to extract any further useful data.
Evaluation of Solid Modeling Software for Finite Element Analysis of Woven Ceramic Matrix Composites
NASA Technical Reports Server (NTRS)
Nemeth, Noel N.; Mital, Subodh; Lang, Jerry
2010-01-01
Three computer programs, used for the purpose of generating 3-D finite element models of the Repeating Unit Cell (RUC) of a textile, were examined for suitability to model woven Ceramic Matrix Composites (CMCs). The programs evaluated were the open-source available TexGen, the commercially available WiseTex, and the proprietary Composite Material Evaluator (COMATE). A five-harness-satin (5HS) weave for a melt-infiltrated (MI) silicon carbide matrix and silicon carbide fiber was selected as an example problem and the programs were tested for their ability to generate a finite element model of the RUC. The programs were also evaluated for ease-of-use and capability, particularly for the capability to introduce various defect types such as porosity, ply shifting, and nesting of a laminate. Overall, it was found that TexGen and WiseTex were useful for generating solid models of the tow geometry; however, there was a lack of consistency in generating well-conditioned finite element meshes of the tows and matrix. TexGen and WiseTex were both capable of allowing collective and individual shifting of tows within a ply and WiseTex also had a ply nesting capability. TexGen and WiseTex were sufficiently userfriendly and both included a Graphical User Interface (GUI). COMATE was satisfactory in generating a 5HS finite element mesh of an idealized weave geometry but COMATE lacked a GUI and was limited to only 5HS and 8HS weaves compared to the larger amount of weave selections available with TexGen and WiseTex.
Experimental studies of nuclear matrix elements for neutrino-less ββ decays
Ejiri, H.
2013-12-30
Nuclear matrix elements M{sup 0ν} for neutrino less double beta decays (0νββ) are crucial for neutrino studies in 0νββ experiments. The neutrino mass to be studied is sensitive to M{sup 0ν}, while theoretical calculations for M{sup 0ν} are hard. Thus experimental studies of nuclear structures and single β matrix elements M{sub β} associated with 0νββ are useful to confirm and help the theoretical calculations. This reports briefly experimental methods and recent charge exchange reaction studies for M{sub β}. The single β elements for M{sup ±}(2{sup −}) associated with M{sup 0ν}(2{sup −}), which is the major component of M{sup 0ν}, are found to be reduced (quenched) much by the spin isospin correlation and the nuclear medium (non-nucleonic isobar) effect. The present result suggests that the spin isospin components of M{sup 0ν} is fairly reduced (quenched)
Differential cross sections and spin density matrix elements for the reaction gamma p -> p omega
M. Williams, D. Applegate, M. Bellis, C.A. Meyer
2009-12-01
High-statistics differential cross sections and spin density matrix elements for the reaction gamma p -> p omega have been measured using the CLAS at Jefferson Lab for center-of-mass (CM) energies from threshold up to 2.84 GeV. Results are reported in 112 10-MeV wide CM energy bins, each subdivided into cos(theta_CM) bins of width 0.1. These are the most precise and extensive omega photoproduction measurements to date. A number of prominent structures are clearly present in the data. Many of these have not previously been observed due to limited statistics in earlier measurements.
NASA Astrophysics Data System (ADS)
Genest, Vincent X.; Miki, Hiroshi; Vinet, Luc; Zhedanov, Alexei
2014-01-01
The multivariate Meixner polynomials are shown to arise as matrix elements of unitary representations of the SO(d, 1) group on oscillator states. These polynomials depend on d discrete variables and are orthogonal with respect to the negative multinomial distribution. The emphasis is put on the bivariate case for which the SO(2, 1) connection is used to derive the main properties of the polynomials: orthogonality relation, raising/lowering relations, generating function, recurrence relations and difference equations as well as explicit expressions in terms of standard (univariate) Krawtchouk and Meixner polynomials. It is explained how these results generalize directly to d variables.
Spin Density Matrix Elements from {rho}{sup 0} and {phi} Meson Electroproduction at HERMES
Borissov, A.
2009-03-23
Exclusive production of {rho}{sup 0} and {phi} mesons on hydrogen and deuterium targets is studied in the HERMES kinematic region 1matrix elements are presented. Violation of s-Channel Helicity Conservation is observed through several non-zero values of SDMEs for {rho}{sup 0}, but not for {phi}. In exclusive {rho}{sup 0} production on the proton an indication is observed of a contribution of unnatural-parity exchange amplitudes, for which the dependence on Q{sup 2} and t' is shown.
Electron-H2 Collisions Studied Using the Finite Element Z-Matrix Method
NASA Technical Reports Server (NTRS)
Huo, Winifred M.; Brown, David; Langhoff, Stephen R. (Technical Monitor)
1997-01-01
We have applied the Z-matrix method, using a mixed basis of finite elements and Gaussians, to study e-H2 elastic and inelastic collisions. Special attention is paid to the quality of the basis set and the treatment of electron correlation. The calculated cross sections are invariant, to machine accuracy, with respect to the choice of parameters a, b, d, e as long as they satisfy Equation (3). However, the log derivative approach, i.e., the choice a = -e = 1, b = d = 0 appears to converge slightly faster than other choices. The cross sections agree well with previous theoretical results. Comparison will be made with available experimental data.
Reduced matrix elements of spin–spin interactions for the atomic f-electron configurations
Yeung, Y.Y.
2014-03-15
A re-examination of some major references on the intra-atomic magnetic interactions over the last six decades reveals that there exist some gaps or puzzles concerning the previous studies of the spin–spin interactions for the atomic f-shell electrons. Hence, tables are provided for the relevant reduced matrix elements of the four double-tensor operators z{sub r} (r=1,2,3, and 4) of rank 2 in both the orbital and spin spaces. The range of the tables covers all states of the configurations from f{sup 4} to f{sup 7}.
Spin Density Matrix Elements in exclusive production of ω mesons at Hermes
NASA Astrophysics Data System (ADS)
Marianski, B.; Terkulov, A.
2014-03-01
Spin density matrix elements have been determined for exclusive ω meson production on hydrogen and deuterium targets, in the kinematic region of 1.0 < Q2 < 10.0 GeV2, 3.0 < W < 6.3 GeV and -t' < 0.2 GeV2. The data, from which SDMEs are determined, were accumulated with the HERMES forward spectrometer during the running period of 1996 to 2007 using the 27.6 GeV electron or positron beam of HERA. A sizable contribution of unnatural parity exchange amplitudes is found for exclusive ω meson production.
Stochastic method with low mode substitution for nucleon isovector matrix elements
NASA Astrophysics Data System (ADS)
Yang, Yi-Bo; Alexandru, Andrei; Draper, Terrence; Gong, Ming; Liu, Keh-Fei; χ QCD Collaboration
2016-02-01
We introduce a stochastic method with low-mode substitution to evaluate the connected three-point functions. The isovector matrix elements of the nucleon for the axial-vector coupling gA3, scalar couplings gS3 and the quark momentum fraction ⟨x ⟩u -d are calculated with overlap fermion on 2 +1 flavor domain-wall configurations on a 243×64 lattice at mπ=330 MeV with lattice spacing a =0.114 fm .
Many-body correlations of QRPA in nuclear matrix elements of double-beta decay
Terasaki, J.
2015-10-28
We present two new ideas on the quasiparticle random-phase approximation (QRPA) approach for calculating nuclear matrix elements of double-beta decay. First, it is necessary to calculate overlaps of the QRPA states obtained on the basis of the ground states of different nuclei. We calculate this overlap using quasiboson vacua as the QRPA ground states. Second, we show that two-particle transfer paths are possible to use for the calculation under the closure approximation. A calculation is shown for {sup 150}Nd→{sup 150}Sm using these two new ideas, and their implication is discussed.
Nucleon matrix elements with Nf=2+1+1 maximally twisted fermions
Simon Dinter, Constantia Alexandrou, Martha Constantinou, Vincent Drach, Karl Jansen, Dru Renner
2010-06-01
We present the first lattice calculation of nucleon matrix elements using four dynamical flavors. We use the Nf=2+1+1 maximally twisted mass formulation. The renormalization is performed non-perturbatively in the RI'-MOM scheme and results are given for the vector and axial vector operators with up to one-derivative. Our calculation of the average momentum of the unpolarized non-singlet parton distribution is presented and compared to our previous results obtained from the Nf=2 case.
Matrix elements of scalar three-electron operators for the atomic f shell
Hansen, J.E.; Judd, B.R.; Crosswhite H.
1996-01-01
Tables are provided for the matrix elements of an orthogonal set of Hermitian three-electron operators t{sub i} for the states of the f shell. The t{sub i} are scalar with respect to the total spin S and total orbital angular momentum L, and they are among the effective operators needed to be included in an f-electron Hamiltonian in order to represent the coupling of the ground configuration f{sup N} to excited configurations via the interelectronic Coulomb interaction. 15 refs., 2 tabs.
Benes, Petr; Simkovic, Fedor
2009-11-09
The nuclear matrix elements M{sup 0v} of the neutrinoless double beta decay (0v{beta}{beta}-decay) are systematically evaluated using the self-consistent renormalized quasiparticle random phase approximation (SRQRPA). The residual interaction and the two-nucleon short-range correlations are derived from the charge-dependent Bonn (CD-Bonn) potential. The importance of further progress in the calculation of the 0v{beta}{beta}-decay nuclear matrix elements is stressed.
Measurement of the top quark mass in the lepton+jets final state with the matrix element method
Abazov, V.M.; Abbott, B.; Abolins, M.; Acharya, B.S.; Adams, M.; Adams, T.; Agelou, M.; Aguilo, E.; Ahn, S.H.; Ahsan, M.; Alexeev, G.D.; /Buenos Aires U. /Rio de Janeiro, CBPF /Sao Paulo, IFT /Alberta U. /Simon Fraser U. /York U., Canada /McGill U. /Hefei, CUST /Andes U., Bogota /Charles U. /Prague, Tech. U.
2006-09-01
We present a measurement of the top quark mass with the Matrix Element method in the lepton+jets final state. As the energy scale for calorimeter jets represents the dominant source of systematic uncertainty, the Matrix Element likelihood is extended by an additional parameter, which is defined as a global multiplicative factor applied to the standard energy scale. The top quark mass is obtained from a fit that yields the combined statistical and systematic jet energy scale uncertainty.
NASA Astrophysics Data System (ADS)
CYKLIS, P.
2001-07-01
In positive-displacement compressor manifolds there are pressure pulsations due to their cyclic operation. The analysis of pressure pulsations in the compressor manifolds is important for various reasons: they directly affect the quantity of energy required for medium compression due to dynamic pressure charging, or inversely, dynamic suppression of suction and discharge processes; they cause mechanical vibrations of compressed gas piping network, they cause aerodynamic and mechanical noise; they affect the dynamics of working valves in valve compressors, they intensify the process of heat convection in heat exchangers in the gas network. The Helmholtz model used so far, which is the basis for users, who deal with pressure pulsation damping, contains many simplifying assumptions. This is because; a straight pipe segment substitutes each element of the piping system. In many cases this model is insufficient. An attempt of the analysis of other shapes was presented in references [1-3] but only simple geometry elements were considered. In other papers [4-8] the influence of the mean flow velocity caused problems. In the presented method, on the basis of pressure pulsation measurement results, firstly a division into the forward and backward going wave is determined, then the elements of the scattering (transmittance) matrix are calculated defining the installation element. This allows introducing the correction for gas mean velocity. The results of the method using correction for the gas mean velocity have been compared with the results without correction and Helmholtz model showing better accuracy.
Precision Measurement of the Neutron Twist-3 Matrix Element d2n: Probing Color Forces
NASA Astrophysics Data System (ADS)
Posik, M.; Flay, D.; Parno, D. S.; Allada, K.; Armstrong, W.; Averett, T.; Benmokhtar, F.; Bertozzi, W.; Camsonne, A.; Canan, M.; Cates, G. D.; Chen, C.; Chen, J.-P.; Choi, S.; Chudakov, E.; Cusanno, F.; Dalton, M. M.; Deconinck, W.; de Jager, C. W.; Deng, X.; Deur, A.; Dutta, C.; El Fassi, L.; Franklin, G. B.; Friend, M.; Gao, H.; Garibaldi, F.; Gilad, S.; Gilman, R.; Glamazdin, O.; Golge, S.; Gomez, J.; Guo, L.; Hansen, O.; Higinbotham, D. W.; Holmstrom, T.; Huang, J.; Hyde, C.; Ibrahim, H. F.; Jiang, X.; Jin, G.; Katich, J.; Kelleher, A.; Kolarkar, A.; Korsch, W.; Kumbartzki, G.; LeRose, J. J.; Lindgren, R.; Liyanage, N.; Long, E.; Lukhanin, A.; Mamyan, V.; McNulty, D.; Meziani, Z.-E.; Michaels, R.; Mihovilovič, M.; Moffit, B.; Muangma, N.; Nanda, S.; Narayan, A.; Nelyubin, V.; Norum, B.; Nuruzzaman; Oh, Y.; Peng, J. C.; Qian, X.; Qiang, Y.; Rakhman, A.; Riordan, S.; Saha, A.; Sawatzky, B.; Shabestari, M. H.; Shahinyan, A.; Širca, S.; Solvignon, P.; Subedi, R.; Sulkosky, V.; Tobias, W. A.; Troth, W.; Wang, D.; Wang, Y.; Wojtsekhowski, B.; Yan, X.; Yao, H.; Ye, Y.; Ye, Z.; Yuan, L.; Zhan, X.; Zhang, Y.; Zhang, Y.-W.; Zhao, B.; Zheng, X.; Jefferson Lab Hall A Collaboration
2014-07-01
Double-spin asymmetries and absolute cross sections were measured at large Bjorken x (0.25≤x ≤0.90), in both the deep-inelastic and resonance regions, by scattering longitudinally polarized electrons at beam energies of 4.7 and 5.9 GeV from a transversely and longitudinally polarized He3 target. In this dedicated experiment, the spin structure function g2He3 was determined with precision at large x, and the neutron twist-3 matrix element d2n was measured at ⟨Q2⟩ of 3.21 and 4.32 GeV2/c2, with an absolute precision of about 10-5. Our results are found to be in agreement with lattice QCD calculations and resolve the disagreement found with previous data at ⟨Q2⟩=5 GeV2/c2. Combining d2n and a newly extracted twist-4 matrix element f2n, the average neutron color electric and magnetic forces were extracted and found to be of opposite sign and about 30 MeV /fm in magnitude.
Precision Measurement of the Neutron Twist-3 Matrix Element dn2: Probing Color Forces
Posik, Matthew; Flay, David; Parno, Diana; Allada, Kalyan; Armstrong, Whitney; Averett, Todd; Benmokhtar, Fatiha; Bertozzi, William; Camsonne, Alexandre; Canan, Mustafa; Cates, Gordon; Chen, Chunhua; Chen, Jian-Ping; Choi, Seonho; Chudakov, Eugene; Cusanno, Francesco; Dalton, Mark; Deconinck, Wouter; De Jager, Cornelis; Deng, Xiaoyan; Deur, Alexandre; Dutta, Chiranjib; El Fassi, Lamiaa; Franklin, Gregg; Friend, Megan; Gao, Haiyan; Garibaldi, Franco; Gilad, Shalev; Gilman, Ronald; Glamazdin, Oleksandr; Golge, Serkan; Gomez, Javier; Guo, Lei; Hansen, Jens-Ole; Higinbotham, Douglas; Holmstrom, Timothy; Huang, J; Hyde, Charles; Ibrahim Abdalla, Hassan; Jiang, Xiaodong; Jin, Ge; Katich, Joseph; Kelleher, Aidan; Kolarkar, Ameya; Korsch, Wolfgang; Kumbartzki, Gerfried; LeRose, John; Lindgren, Richard; Liyanage, Nilanga; Long, Elena; Lukhanin, Oleksandr; Mamyan, Vahe; McNulty, Dustin; Meziani, Zein-Eddine; Michaels, Robert; Mihovilovic, Miha; Moffit, Bryan; Muangma, Navaphon; Nanda, Sirish; Narayan, Amrendra; Nelyubin, Vladimir; Norum, Blaine; Nuruzzaman, nfn; Oh, Yongseok; Peng, Jen-chieh; Qian, Xin; Qiang, Yi; Rakhman, Abdurahim; Riordan, Seamus; Saha, Arunava; Sawatzky, Bradley; Hashemi Shabestari, Mitra; Shahinyan, Albert; Sirca, Simon; Solvignon-Slifer, Patricia; Subedi, Ramesh; Sulkosky, Vincent; Tobias, William; Troth, Wolfgang; Wang, Diancheng; Wang, Y; Wojtsekhowski, Bogdan; Yan, X; Yao, Huan; Ye, Yunxiu; Ye, Zhihong; Yuan, Lulin; Zhan, X; Zhang, Y; Zhang, Y -W; Zhao, Bo; Zheng, Xiaochao
2014-07-01
Double-spin asymmetries and absolute cross sections were measured at large Bjorken x (0.25 lte x lte 0.90), in both the deep-inelastic and resonance regions, by scattering longitudinally polarized electrons at beam energies of 4.7 and 5.9 GeV from a transversely and longitudinally polarized 3He target. In this dedicated experiment, the spin structure function g2 on 3He was determined with precision at large x, and the neutron twist-three matrix element dn2 was measured at ?Q2? of 3.21 and 4.32 GeV2/c2, with an absolute precision of about 10?5. Our results are found to be in agreement with lattice QCD calculations and resolve the disagreement found with previous data at ?Q2?= 5 GeV2/c2. Combining dn2 and a newly extracted twist-four matrix element, fn2, the average neutron color electric and magnetic forces were extracted and found to be of opposite sign and about 60 MeV/fm in magnitude.
NASA Astrophysics Data System (ADS)
Monthus, Cécile
2016-07-01
For short-ranged disordered quantum models in one dimension, the many-body-localization is analyzed via the adaptation to the many-body context (Serbyn et al 2015 Phys. Rev. X 5 041047) of the Thouless point of view on the Anderson transition: the question is whether a local interaction between two long chains is able to reshuffle completely the eigenstates (delocalized phase with a volume-law entanglement) or whether the hybridization between tensor states remains limited (many-body-localized phase with an area-law entanglement). The central object is thus the level of hybridization induced by the matrix elements of local operators, as compared with the difference of diagonal energies. The multifractal analysis of these matrix elements of local operators is used to analyze the corresponding statistics of resonances. Our main conclusion is that the critical point is characterized by the strong-multifractality spectrum f(0≤slant α ≤slant 2)=\\fracα{2} , well known in the context of Anderson localization in spaces of effective infinite dimensionality, where the size of the Hilbert space grows exponentially with the volume. Finally, the possibility of a delocalized non-ergodic phase near criticality is discussed.
Symbolic algorithms for the computation of Moshinsky brackets and nuclear matrix elements
NASA Astrophysics Data System (ADS)
Ursescu, D.; Tomaselli, M.; Kuehl, T.; Fritzsche, S.
2005-12-01
To facilitate the use of the extended nuclear shell model (NSM), a FERMI module for calculating some of its basic quantities in the framework of MAPLE is provided. The Moshinsky brackets, the matrix elements for several central and non-central interactions between nuclear two-particle states as well as their expansion in terms of Talmi integrals are easily given within a symbolic formulation. All of these quantities are available for interactive work. Program summaryTitle of program:Fermi Catalogue identifier:ADVO Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVO Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Licensing provisions:None Computer for which the program is designed and others on which is has been tested:All computers with a licence for the computer algebra package MAPLE [Maple is a registered trademark of Waterloo Maple Inc., produced by MapleSoft division of Waterloo Maple Inc.] Instalations:GSI-Darmstadt; University of Kassel (Germany) Operating systems or monitors under which the program has beentested: WindowsXP, Linux 2.4 Programming language used:MAPLE 8 and 9.5 from MapleSoft division of Waterloo Maple Inc. Memory required to execute with typical data:30 MB No. of lines in distributed program including test data etc.:5742 No. of bytes in distributed program including test data etc.:288 939 Distribution program:tar.gz Nature of the physical problem:In order to perform calculations within the nuclear shell model (NSM), a quick and reliable access to the nuclear matrix elements is required. These matrix elements, which arise from various types of forces among the nucleons, can be calculated using Moshinsky's transformation brackets between relative and center-of-mass coordinates [T.A. Brody, M. Moshinsky, Tables of Transformation Brackets, Monografias del Instituto de Fisica, Universidad Nacional Autonoma de Mexico, 1960] and by the proper use of the nuclear states in different coupling notations
Top quark mass measurement from dilepton events at CDF II with the matrix-element method
Abulencia, A.; Acosta, D.; Adelman, Jahred A.; Affolder, T.; Akimoto, T.; Albrow, M.G.; Ambrose, D.; Amerio, S.; Amidei, D.; Anastassov, A.; Anikeev, K.; /Taiwan, Inst. Phys. /Argonne /Barcelona, IFAE /Baylor U. /INFN, Bologna /Bologna U. /Brandeis U. /UC, Davis /UCLA /UC, San Diego /UC, Santa Barbara
2006-05-01
We describe a measurement of the top quark mass using events with two charged leptons collected by the CDF II detector from p{bar p} collisions with {radical}s = 1.96 TeV at the Fermilab Tevatron. The likelihood in top mass is calculated for each event by convoluting the leading order matrix element describing q{bar q} {yields} t{bar t} {yields} b{ell}{nu}{sub {ell}}{bar b}{ell}{prime} {nu}{sub {ell}}, with detector resolution functions. The presence of background events in the data sample is modeled using similar calculations involving the matrix elements for major background processes. In a data sample with integrated luminosity of 340 pb{sup -1}, we observe 33 candidate events and measure M{sub top} = 165.2 {+-} 6.1(stat.) {+-} 3.4(syst.) GeV/c{sup 2}. This measurement represents the first application of this method to events with two charged leptons and is the most precise single measurement of the top quark mass in this channel.
Top Quark Mass Measurement in the Lepton plus Jets Channel Using a Modified Matrix Element Method
Aaltonen, T.; Adelman, J.; Akimoto, T.; Alvarez Gonzalez, B.; Amerio, S.; Amidei, D.; Anastassov, A.; Annovi, A.; Antos, J.; Apollinari, G.; Apresyan, A.; /Purdue U. /Waseda U.
2008-12-01
The authors report a measurement of the top quark mass, m{sub t}, obtained from p{bar p} collisions at {radical}s = 1.96 TeV at the Fermilab Tevatron using the CDF II detector. They analyze a sample corresponding to an integrated luminosity of 1.9 rfb{sup -1}. They select events with an electron or muon, large missing transverse energy, and exactly four high-energy jets in the central region of the detector, at least one of which is tagged as coming from a b quark. They calculate a signal likelihood using a matrix element integration method, where the matrix element is modified by using effective propagators to take into account assumptions on event kinematics. The event likelihood is a function of m{sub t} and a parameter JES that determines in situ the calibration of the jet energies. They use a neural network discriminant to distinguish signal from background events. They also apply a cut on the peak value of each event likelihood curve to reduce the contribution of background and badly reconstructed events. Using the 318 events that pass all selection criteria, they find m{sub t} = 172.7 {+-} 1.8 (stat. + JES) {+-} 1.2(syst.) GeV/c{sup 2}.
Precision measurement of the neutron twist-3 matrix element d(2)(n): probing color forces.
Posik, M; Flay, D; Parno, D S; Allada, K; Armstrong, W; Averett, T; Benmokhtar, F; Bertozzi, W; Camsonne, A; Canan, M; Cates, G D; Chen, C; Chen, J-P; Choi, S; Chudakov, E; Cusanno, F; Dalton, M M; Deconinck, W; de Jager, C W; Deng, X; Deur, A; Dutta, C; El Fassi, L; Franklin, G B; Friend, M; Gao, H; Garibaldi, F; Gilad, S; Gilman, R; Glamazdin, O; Golge, S; Gomez, J; Guo, L; Hansen, O; Higinbotham, D W; Holmstrom, T; Huang, J; Hyde, C; Ibrahim, H F; Jiang, X; Jin, G; Katich, J; Kelleher, A; Kolarkar, A; Korsch, W; Kumbartzki, G; LeRose, J J; Lindgren, R; Liyanage, N; Long, E; Lukhanin, A; Mamyan, V; McNulty, D; Meziani, Z-E; Michaels, R; Mihovilovič, M; Moffit, B; Muangma, N; Nanda, S; Narayan, A; Nelyubin, V; Norum, B; Nuruzzaman; Oh, Y; Peng, J C; Qian, X; Qiang, Y; Rakhman, A; Riordan, S; Saha, A; Sawatzky, B; Shabestari, M H; Shahinyan, A; Širca, S; Solvignon, P; Subedi, R; Sulkosky, V; Tobias, W A; Troth, W; Wang, D; Wang, Y; Wojtsekhowski, B; Yan, X; Yao, H; Ye, Y; Ye, Z; Yuan, L; Zhan, X; Zhang, Y; Zhang, Y-W; Zhao, B; Zheng, X
2014-07-11
Double-spin asymmetries and absolute cross sections were measured at large Bjorken x (0.25≤x≤0.90), in both the deep-inelastic and resonance regions, by scattering longitudinally polarized electrons at beam energies of 4.7 and 5.9 GeV from a transversely and longitudinally polarized (3)He target. In this dedicated experiment, the spin structure function g(2)((3)He) was determined with precision at large x, and the neutron twist-3 matrix element d(2)(n) was measured at ⟨Q(2)⟩ of 3.21 and 4.32 GeV(2)/c(2), with an absolute precision of about 10(-5). Our results are found to be in agreement with lattice QCD calculations and resolve the disagreement found with previous data at ⟨Q(2)⟩=5 GeV(2)/c(2). Combining d(2)(n) and a newly extracted twist-4 matrix element f(2)(n), the average neutron color electric and magnetic forces were extracted and found to be of opposite sign and about 30 MeV/fm in magnitude. PMID:25062166
Measurement of single top quark production at D0 using a matrix element method
Mitrevski, Jovan Pavle; /Columbia U.
2007-07-01
Until now, the top quark has only been observed produced in pairs, by the strong force. According to the standard model, it can also be produced singly, via an electroweak interaction. Top quarks produced this way provide powerful ways to test the charged-current electroweak interactions of the top quark, to measure |V{sub tb}|, and to search for physics beyond the standard model. This thesis describes the application of the matrix element analysis technique to the search for single top quark production with the D0 detector using 0.9 fb{sup -1} of Run II data. From a comparison of the matrix element discriminants between data and the background model, assuming a Standard Model s-channel to t-channel cross section ratio of {sigma}{sub s}/{sigma}{sub t} = 0.44, we measure the single top quark production cross section: {sigma}(p{bar p} {yields} tb + X, tqb + X) = 4.8{sub -1.4}{sup +1.6} pb. This result has a p-value of 0.08%, corresponding to a 3.2 standard deviation Gaussian equivalent significance.
Analytical techniques for the evaluation of asymptotic matrix elements in electromagnetic problems
NASA Astrophysics Data System (ADS)
Park, Seong-Ook
1997-07-01
Analytical solutions of the asymptotic part of impedance matrix elements, whose integrand contains either a singularity or exhibits a strongly oscillatory behavior, are developed for the analysis of various geometries. These developed formulas are applied using the Method of Moments (MoM), to the analysis of wire antennas, planar transmission lines, and printed circuits including microstrip dipoles, asymmetric gap discontinuities, and arbitrary shapes of planar circuits, where computational speed and accuracy are important. The first example considered is the singularity treatment in integrands of wire-type antenna problems. To overcome the problem, first the singularity is subtracted out, and then this part is integrated analytically in the transform domain. Asymptotic matrix elements of the multilayer planar transmission lines are solved analytically by using Chebyshev polynomial basis functions in the spectral domain. These results are applied to open and coupled microstrip lines. The analysis of printed circuits involves Sommerfeld-type integrals which are extremely difficult to evaluate. To enhance the speed and accuracy of Sommerfeld-type integrals, this dissertation, for the first time, presents an analytical transformation technique. This formulation allows the infinite double integral of the asymptotic part of the impedance matrix to be transformed into a finite one-dimensional integral. Using this approach, the asymptotic part of self and mutual interactions between the triangular edge mode basis functions, along an electrically narrow strip, is solved analytically. These results are applied to microstrip dipoles and asymmetric gap discontinuities. This formula provides highly accurate results with minimal computational effort. This work is also extended to roof- top subdomain basis functions, and the obtained results can be used to solve arbitrarily shaped planar geometries. As an example, computed results of the Radar Cross Section (RCS) of a microstrip
Measurement of the top quark mass in the dilepton final state using the matrix element method
Grohsjean, Alexander; /Munich U.
2008-12-01
The top quark, discovered in 1995 by the CDF and D0 experiments at the Fermilab Tevatron Collider, is the heaviest known fundamental particle. The precise knowledge of its mass yields important constraints on the mass of the yet-unobserved Higgs boson and allows to probe for physics beyond the Standard Model. The first measurement of the top quark mass in the dilepton channel with the Matrix Element method at the D0 experiment is presented. After a short description of the experimental environment and the reconstruction chain from hits in the detector to physical objects, a detailed review of the Matrix Element method is given. The Matrix Element method is based on the likelihood to observe a given event under the assumption of the quantity to be measured, e.g. the mass of the top quark. The method has undergone significant modifications and improvements compared to previous measurements in the lepton+jets channel: the two undetected neutrinos require a new reconstruction scheme for the four-momenta of the final state particles, the small event sample demands the modeling of additional jets in the signal likelihood, and a new likelihood is designed to account for the main source of background containing tauonic Z decay. The Matrix Element method is validated on Monte Carlo simulated events at the generator level. For the measurement, calibration curves are derived from events that are run through the full D0 detector simulation. The analysis makes use of the Run II data set recorded between April 2002 and May 2008 corresponding to an integrated luminosity of 2.8 fb{sup -1}. A total of 107 t{bar t} candidate events with one electron and one muon in the final state are selected. Applying the Matrix Element method to this data set, the top quark mass is measured to be m{sub top}{sup Run IIa} = 170.6 {+-} 6.1(stat.){sub -1.5}{sup +2.1}(syst.)GeV; m{sub top}{sup Run IIb} = 174.1 {+-} 4.4(stat.){sub -1.8}{sup +2.5}(syst.)GeV; m{sub top}{sup comb} = 172.9 {+-} 3.6(stat
ERIC Educational Resources Information Center
Arnold, Randy J.; Arndt, Brett; Blaser, Emilia; Blosser, Chris; Caulton, Dana; Chung, Won Sog; Fiorenza, Garrett; Heath, Wyatt; Jacobs, Alex; Kahng, Eunice; Koh, Eun; Le, Thao; Mandla, Kyle; McCory, Chelsey; Newman, Laura; Pithadia, Amit; Reckelhoff, Anna; Rheinhardt, Joseph; Skljarevski, Sonja; Stuart, Jordyn; Taylor, Cassie; Thomas, Scott; Tse, Kyle; Wall, Rachel; Warkentien, Chad
2011-01-01
A multivitamin tablet and liquid are analyzed for the elements calcium, magnesium, iron, zinc, copper, and manganese using atomic absorption spectrometry. Linear calibration and standard addition are used for all elements except calcium, allowing for an estimate of the matrix effects encountered for this complex sample. Sample preparation using…
Measurement of the RMS Parity Violating Matrix Element in URANIUM-239
NASA Astrophysics Data System (ADS)
Zhu, Xianzhou (Joe).
We report the first determination of the Root -Mean-Square (RMS) parity violating matrix element in a compound nucleus (CN) system, ^{239 }U. The experiment was performed using the intense pulsed epithermal neutron beam available at the Los Alamos Neutron Scattering Center (LANSCE). The helicity dependence of neutron transmission through a spin zero target (^{238}U) is measured for neutron energies from 6 eV to 300 eV. Parity violation is analyzed on 17 p-wave resonances among which five show 2sigma or larger effects. The largest is a 7sigma effect at the 63.5 eV resonance which shows a parity violating asymmetry of p = 2.6%. A likelihood analysis is performed on these 17 parity violating asymmetries, and the RMS parity violating matrix element is determined for the first time to be M = 0.59_sp{-0.25}{+0.50} meV which corresponds to a parity violating spreading width ofGamma^{PV} = (1.0 {+1.7atop -0.8} ) times 10^{-7} {rm eV}.Using statistical nuclear spectroscopy, we are able to relate M to the effective nucleon-nucleon (NN) interaction. The result is | alpha_{p}| ~ (4 {+4atop -2} ) times 10^{-7} where alpha_{p} is the ratio of the parity violating strength to the parity conserving strength in the effective NN interaction. This agrees qualitatively with the estimate of free NN interaction. The consistency of the experimental measurement with expectation suggests that the manifestation of parity violating NN interaction in CN is understood. It is a challenging problem for the theorists to relate the RMS matrix element in the CN to the underlying NN interaction, therefore providing alternative ways to determine the Desplanques -Donoghue-Holstein (DDH) parameters of the NN interaction. The success of the parity violation study also validates the proposed experiment of studying the time reversal symmetry violation utilizing the large enhancement in the CN.
NASA Astrophysics Data System (ADS)
Yao, J. M.; Song, L. S.; Hagino, K.; Ring, P.; Meng, J.
2015-02-01
We report a systematic study of nuclear matrix elements (NMEs) in neutrinoless double-β decays with a state-of-the-art beyond-mean-field covariant density functional theory. The dynamic effects of particle-number and angular-momentum conservations as well as quadrupole shape fluctuations are taken into account with projections and generator coordinate method for both initial and final nuclei. The full relativistic transition operator is adopted to calculate the NMEs. The present systematic studies show that in most of the cases there is a much better agreement with the previous nonrelativistic calculation based on the Gogny force than in the case of the nucleus 150Nd found by Song et al. [Phys. Rev. C 90, 054309 (2014), 10.1103/PhysRevC.90.054309]. In particular, we find that the total NMEs can be well approximated by the pure axial-vector coupling term with a considerable reduction of the computational effort.
OMC studies for the matrix elements in ββ decay
Zinatulina, D.; Brudanin, V.; Egorov, V.; Shirchenko, M.; Vasiliev, R.; Yyutlandov, I.; Briançon, Ch.; Petitjean, C.
2013-12-30
Energy and time spectra of gamma-rays following μ-capture in natural Kr, Se, Cd and Sm, as well as isotopic enriched {sup 82}Kr, {sup 76}Se, {sup 106}Cd and {sup 150}Sm, have been measured. Total life-times of muons in different isotopes, as well as partial μ-capture rates to the excited states of {sup 48}Sc, {sup 76}As and {sup 106}Ag, were extracted. These results are discussed in the context of the double-beta decay matrix elements. The data are also compared with data from theoretical calculations and with data from charge-exchange reactions on {sup 48}Ti. It is the first time that μ-capture and charge-exchange reaction data are being compared in the context of ββ decay.
Nucleon distribution amplitudes and proton decay matrix elements on the lattice
Braun, Vladimir M.; Goeckeler, Meinulf; Kaltenbrunner, Thomas; Warkentin, Nikolaus; Horsley, Roger; Zanotti, James M.; Nakamura, Yoshifumi; Pleiter, Dirk; Rakow, Paul E. L.; Schaefer, Andreas; Schierholz, Gerrit; Stueben, Hinnerk
2009-02-01
Baryon distribution amplitudes (DAs) are crucial for the theory of hard exclusive reactions. We present a calculation of the first few moments of the leading-twist nucleon DA within lattice QCD. In addition we deal with the normalization of the next-to-leading (twist-four) DAs. The matrix elements determining the latter quantities are also responsible for proton decay in grand unified theories. Our lattice evaluation makes use of gauge field configurations generated with two flavors of clover fermions. The relevant operators are renormalized nonperturbatively with the final results given in the MS scheme. We find that the deviation of the leading-twist nucleon DA from its asymptotic form is less pronounced than sometimes claimed in the literature.
Thermoplastic matrix composites - Finite-element analysis of mode I and mode II failure specimens
NASA Technical Reports Server (NTRS)
Bankert, Ray J.; Lambropoulos, Nicholas D.; Shephard, Mark S.; Sternstein, Sanford S.
1989-01-01
A finite-element analysis was conducted to evaluate the stress distributions within mode I and mode II failure specimens, assuming both isotropic and orthotropic elastic material properties. The effects of anisotropy on both the magnitude and the location of the highest stress concentration at the vicinity of the crack tip are significant. The results from modeling realistic blunt crack tip geometry and resin-rich zones imply that local variations in the microstructure strongly influence the stress state near the crack tip and therefore the measured fracture properties. In addition, the features of a viscoelastic model for thermoplastic matrices are described. This model will be used in future investigations of matrix-dominated failure phenomena.
Influence of Pairing on the Nuclear Matrix Elements of the Neutrinoless ββ Decays
NASA Astrophysics Data System (ADS)
Caurier, E.; Menéndez, J.; Nowacki, F.; Poves, A.
2008-02-01
We study in this Letter the neutrinoless double beta decay nuclear matrix elements (NME’s) in the framework of the interacting shell model. We analyze them in terms of the total angular momentum of the decaying neutron pair and as a function of the seniority truncations in the nuclear wave functions. This point of view turns out to be very adequate to gauge the accuracy of the NME’s predicted by different nuclear models. In addition, it gives back the protagonist role in this process to the pairing interaction, the one which is responsible for the very existence of double beta decay emitters. We show that low seniority approximations, comparable to those implicit in the quasiparticle RPA in a spherical basis, tend to overestimate the NME’s in several decays.
Measurement of spin correlation in tt production using a matrix element approach.
Abazov, V M; Abbott, B; Acharya, B S; Adams, M; Adams, T; Alexeev, G D; Alkhazov, G; Alton, A; Alverson, G; Alves, G A; Ancu, L S; Aoki, M; Arov, M; Askew, A; Åsman, B; Atramentov, O; Avila, C; BackusMayes, J; Badaud, F; Bagby, L; Baldin, B; Bandurin, D V; Banerjee, S; Barberis, E; Baringer, P; Barreto, J; Bartlett, J F; Bassler, U; Bazterra, V; Beale, S; Bean, A; Begalli, M; Begel, M; Belanger-Champagne, C; Bellantoni, L; Beri, S B; Bernardi, G; Bernhard, R; Bertram, I; Besançon, M; Beuselinck, R; Bezzubov, V A; Bhat, P C; Bhatnagar, V; Blazey, G; Blessing, S; Bloom, K; Boehnlein, A; Boline, D; Boos, E E; Borissov, G; Bose, T; Brandt, A; Brandt, O; Brock, R; Brooijmans, G; Bross, A; Brown, D; Brown, J; Bu, X B; Buehler, M; Buescher, V; Bunichev, V; Burdin, S; Burnett, T H; Buszello, C P; Calpas, B; Camacho-Pérez, E; Carrasco-Lizarraga, M A; Casey, B C K; Castilla-Valdez, H; Chakrabarti, S; Chakraborty, D; Chan, K M; Chandra, A; Chen, G; Chevalier-Théry, S; Cho, D K; Cho, S W; Choi, S; Choudhary, B; Cihangir, S; Claes, D; Clutter, J; Cooke, M; Cooper, W E; Corcoran, M; Couderc, F; Cousinou, M-C; Croc, A; Cutts, D; Das, A; Davies, G; De, K; de Jong, S J; De la Cruz-Burelo, E; Déliot, F; Demarteau, M; Demina, R; Denisov, D; Denisov, S P; Desai, S; Deterre, C; DeVaughan, K; Diehl, H T; Diesburg, M; Dominguez, A; Dorland, T; Dubey, A; Dudko, L V; Duggan, D; Duperrin, A; Dutt, S; Dyshkant, A; Eads, M; Edmunds, D; Ellison, J; Elvira, V D; Enari, Y; Evans, H; Evdokimov, A; Evdokimov, V N; Facini, G; Ferbel, T; Fiedler, F; Filthaut, F; Fisher, W; Fisk, H E; Fortner, M; Fox, H; Fuess, S; Garcia-Bellido, A; Gavrilov, V; Gay, P; Geng, W; Gerbaudo, D; Gerber, C E; Gershtein, Y; Ginther, G; Golovanov, G; Goussiou, A; Grannis, P D; Greder, S; Greenlee, H; Greenwood, Z D; Gregores, E M; Grenier, G; Gris, Ph; Grivaz, J-F; Grohsjean, A; Grünendahl, S; Grünewald, M W; Guillemin, T; Guo, F; Gutierrez, G; Gutierrez, P; Haas, A; Hagopian, S; Haley, J; Han, L; Harder, K; Harel, A; Hauptman, J M; Hays, J; Head, T; Hebbeker, T; Hedin, D; Hegab, H; Heinson, A P; Heintz, U; Hensel, C; Heredia-De la Cruz, I; Herner, K; Hesketh, G; Hildreth, M D; Hirosky, R; Hoang, T; Hobbs, J D; Hoeneisen, B; Hohlfeld, M; Hubacek, Z; Huske, N; Hynek, V; Iashvili, I; Illingworth, R; Ito, A S; Jabeen, S; Jaffré, M; Jamin, D; Jayasinghe, A; Jesik, R; Johns, K; Johnson, M; Johnston, D; Jonckheere, A; Jonsson, P; Joshi, J; Jung, A W; Juste, A; Kaadze, K; Kajfasz, E; Karmanov, D; Kasper, P A; Katsanos, I; Kehoe, R; Kermiche, S; Khalatyan, N; Khanov, A; Kharchilava, A; Kharzheev, Y N; Khatidze, D; Kirby, M H; Kohli, J M; Kozelov, A V; Kraus, J; Kulikov, S; Kumar, A; Kupco, A; Kurča, T; Kuzmin, V A; Kvita, J; Lammers, S; Landsberg, G; Lebrun, P; Lee, H S; Lee, S W; Lee, W M; Lellouch, J; Li, L; Li, Q Z; Lietti, S M; Lim, J K; Lincoln, D; Linnemann, J; Lipaev, V V; Lipton, R; Liu, Y; Liu, Z; Lobodenko, A; Lokajicek, M; Lopes de Sa, R; Lubatti, H J; Luna-Garcia, R; Lyon, A L; Maciel, A K A; Mackin, D; Madar, R; Magaña-Villalba, R; Malik, S; Malyshev, V L; Maravin, Y; Martínez-Ortega, J; McCarthy, R; McGivern, C L; Meijer, M M; Melnitchouk, A; Menezes, D; Mercadante, P G; Merkin, M; Meyer, A; Meyer, J; Miconi, F; Mondal, N K; Muanza, G S; Mulhearn, M; Nagy, E; Naimuddin, M; Narain, M; Nayyar, R; Neal, H A; Negret, J P; Neustroev, P; Novaes, S F; Nunnemann, T; Obrant, G; Orduna, J; Osman, N; Osta, J; Otero y Garzón, G J; Padilla, M; Pal, A; Parashar, N; Parihar, V; Park, S K; Parsons, J; Partridge, R; Parua, N; Patwa, A; Penning, B; Perfilov, M; Peters, K; Peters, Y; Petridis, K; Petrillo, G; Pétroff, P; Piegaia, R; Piper, J; Pleier, M-A; Podesta-Lerma, P L M; Podstavkov, V M; Polozov, P; Popov, A V; Prewitt, M; Price, D; Prokopenko, N; Protopopescu, S; Qian, J; Quadt, A; Quinn, B; Rangel, M S; Ranjan, K; Ratoff, P N; Razumov, I; Renkel, P; Rijssenbeek, M; Ripp-Baudot, I; Rizatdinova, F; Rominsky, M; Ross, A; Royon, C; Rubinov, P; Ruchti, R; Safronov, G; Sajot, G; Salcido, P; Sánchez-Hernández, A; Sanders, M P; Sanghi, B; Santos, A S; Savage, G; Sawyer, L; Scanlon, T; Schamberger, R D; Scheglov, Y; Schellman, H; Schliephake, T; Schlobohm, S; Schwanenberger, C; Schwienhorst, R; Sekaric, J; Severini, H; Shabalina, E; Shary, V; Shchukin, A A; Shivpuri, R K; Simak, V; Sirotenko, V; Skubic, P; Slattery, P; Smirnov, D; Smith, K J; Snow, G R; Snow, J; Snyder, S; Söldner-Rembold, S; Sonnenschein, L; Soustruznik, K; Stark, J; Stolin, V; Stoyanova, D A; Strauss, M; Strom, D; Stutte, L; Suter, L; Svoisky, P; Takahashi, M; Tanasijczuk, A; Taylor, W; Titov, M; Tokmenin, V V; Tsai, Y-T; Tsybychev, D; Tuchming, B; Tully, C; Uvarov, L; Uvarov, S; Uzunyan, S; Van Kooten, R; van Leeuwen, W M; Varelas, N; Varnes, E W; Vasilyev, I A; Verdier, P; Vertogradov, L S; Verzocchi, M; Vesterinen, M; Vilanova, D; Vokac, P; Wahl, H D; Wang, M H L S; Warchol, J; Watts, G; Wayne, M; Weber, M; Welty-Rieger, L; White, A; Wicke, D; Williams, M R J; Wilson, G W; Wobisch, M; Wood, D R; Wyatt, T R; Xie, Y; Xu, C; Yacoob, S; Yamada, R; Yang, W-C; Yasuda, T; Yatsunenko, Y A; Ye, Z; Yin, H; Yip, K; Youn, S W; Yu, J; Zelitch, S; Zhao, T; Zhou, B; Zhu, J; Zielinski, M; Zieminska, D; Zivkovic, L
2011-07-15
We determine the fraction of tt events with spin correlation, assuming that the spin of the top quark is either correlated with the spin of the top antiquark as predicted by the standard model or is uncorrelated. For the first time we use a matrix-element-based approach to study tt spin correlation. We use tt → W+ b W- b → ℓ+ νbℓ- ν b final states produced in pp collisions at a center-of-mass energy sqrt(s)=1.96 TeV, where ℓ denotes an electron or a muon. The data correspond to an integrated luminosity of 5.4 fb(-1) and were collected with the D0 detector at the Fermilab Tevatron collider. The result agrees with the standard model prediction. We exclude the hypothesis that the spins of the tt are uncorrelated at the 97.7% C.L. PMID:21838349
Influence of Pairing on the Nuclear Matrix Elements of the Neutrinoless {beta}{beta} Decays
Caurier, E.; Nowacki, F.
2008-02-08
We study in this Letter the neutrinoless double beta decay nuclear matrix elements (NME's) in the framework of the interacting shell model. We analyze them in terms of the total angular momentum of the decaying neutron pair and as a function of the seniority truncations in the nuclear wave functions. This point of view turns out to be very adequate to gauge the accuracy of the NME's predicted by different nuclear models. In addition, it gives back the protagonist role in this process to the pairing interaction, the one which is responsible for the very existence of double beta decay emitters. We show that low seniority approximations, comparable to those implicit in the quasiparticle RPA in a spherical basis, tend to overestimate the NME's in several decays.
Measurement of the top quark mass using the matrix element technique in dilepton final states
Abazov, V. M.; Abbott, B.; Acharya, B. S.; Adams, M.; Adams, T.; Agnew, J. P.; Alexeev, G. D.; Alkhazov, G.; Alton, A.; Askew, A.; et al
2016-08-18
Here, we present a measurement of the top quark mass in pp collisions at a center-of-mass energy of 1.96 TeV at the Fermilab Tevatron collider. The data were collected by the D0 experiment corresponding to an integrated luminosity of 9.7 fb-1. The matrix element technique is applied to tt events in the final state containing leptons (electrons or muons) with high transverse momenta and at least two jets. The calibration of the jet energy scale determined in the lepton+jets final state of tt decays is applied to jet energies. This correction provides a substantial reduction in systematic uncertainties. We obtain amore » top quark mass of mt = 173.93±1.84 GeV.« less
Differential cross sections and spin density matrix elements for the reaction γp→pω
NASA Astrophysics Data System (ADS)
Williams, M.; Applegate, D.; Bellis, M.; Meyer, C. A.; Adhikari, K. P.; Anghinolfi, M.; Baghdasaryan, H.; Ball, J.; Battaglieri, M.; Bedlinskiy, I.; Berman, B. L.; Biselli, A. S.; Bookwalter, C.; Briscoe, W. J.; Brooks, W. K.; Burkert, V. D.; Careccia, S. L.; Carman, D. S.; Cole, P. L.; Collins, P.; Crede, V.; D'Angelo, A.; Daniel, A.; Vita, R. De; Sanctis, E. De; Deur, A.; Dey, B.; Dhamija, S.; Dickson, R.; Djalali, C.; Dodge, G. E.; Doughty, D.; Dugger, M.; Dupre, R.; Alaoui, A. El; Elouadrhiri, L.; Eugenio, P.; Fedotov, G.; Fegan, S.; Fradi, A.; Gabrielyan, M. Y.; Garçon, M.; Gevorgyan, N.; Gilfoyle, G. P.; Giovanetti, K. L.; Girod, F. X.; Gohn, W.; Golovatch, E.; Gothe, R. W.; Griffioen, K. A.; Guidal, M.; Guo, L.; Hafidi, K.; Hakobyan, H.; Hanretty, C.; Hassall, N.; Hicks, K.; Holtrop, M.; Ilieva, Y.; Ireland, D. G.; Ishkhanov, B. S.; Isupov, E. L.; Jawalkar, S. S.; Jo, H. S.; Johnstone, J. R.; Joo, K.; Keller, D.; Khandaker, M.; Khetarpal, P.; Kim, W.; Klein, A.; Klein, F. J.; Krahn, Z.; Kubarovsky, V.; Kuleshov, S. V.; Kuznetsov, V.; Livingston, K.; Lu, H. Y.; Mayer, M.; McAndrew, J.; McCracken, M. E.; McKinnon, B.; Mikhailov, K.; Mirazita, M.; Mokeev, V.; Moreno, B.; Moriya, K.; Morrison, B.; Moutarde, H.; Munevar, E.; Nadel-Turonski, P.; Nepali, C. S.; Niccolai, S.; Niculescu, G.; Niculescu, I.; Niroula, M. R.; Niyazov, R. A.; Osipenko, M.; Ostrovidov, A. I.; Paris, M.; Park, K.; Park, S.; Pasyuk, E.; Pereira, S. Anefalos; Perrin, Y.; Pisano, S.; Pogorelko, O.; Pozdniakov, S.; Price, J. W.; Procureur, S.; Protopopescu, D.; Raue, B. A.; Ricco, G.; Ripani, M.; Ritchie, B. G.; Rosner, G.; Rossi, P.; Sabatié, F.; Saini, M. S.; Salamanca, J.; Salgado, C.; Schott, D.; Schumacher, R. A.; Seraydaryan, H.; Sharabian, Y. G.; Smith, E. S.; Sober, D. I.; Sokhan, D.; Stepanyan, S. S.; Stoler, P.; Strakovsky, I. I.; Strauch, S.; Taiuti, M.; Tedeschi, D. J.; Tkachenko, S.; Ungaro, M.; Vineyard, M. F.; Voutier, E.; Watts, D. P.; Weinstein, L. B.; Weygand, D. P.; Wood, M. H.; Zhang, J.; Zhao, B.
2009-12-01
High-statistics differential cross sections and spin-density matrix elements for the reaction γp→pω have been measured using the CEBAF large acceptance spectrometer (CLAS) at Jefferson Lab for center-of-mass (c.m.) energies from threshold up to 2.84 GeV. Results are reported in 11210-MeV wide c.m. energy bins, each subdivided into cosθc.m.ω bins of width 0.1. These are the most precise and extensive ω photoproduction measurements to date. A number of prominent structures are clearly present in the data. Many of these have not previously been observed due to limited statistics in earlier measurements.
NASA Astrophysics Data System (ADS)
Nordin, Gregory P.; Jones, Michael W.; Kulick, Jeffrey H.; Lindquist, Robert G.; Kowel, Stephen T.
1996-12-01
We describe the design, construction, and performance of the first real-time autostereoscopic 3D display based on the partial pixel 3D display architecture. The primary optical components of the 3D display are an active-matrix liquid crystal display and a diffractive optical element (DOE). The display operates at video frame rates and is driven with a conventional VGA signal. 3D animations with horizontal motion parallax are readily viewable as sets of stereo images. Formation of the virtual viewing slits by diffraction from the partial pixel apertures is experimentally verified. The measured contrast and perceived brightness of the display are excellent, but there are minor flaws in image quality due to secondary images. The source of these images and how they may be eliminated is discussed. The effects of manufacturing-related systematic errors in the DOE are also analyzed.
HELAC-Onia: An automatic matrix element generator for heavy quarkonium physics
NASA Astrophysics Data System (ADS)
Shao, Hua-Sheng
2013-11-01
By the virtues of the Dyson-Schwinger equations, we upgrade the published code HELAC to be capable to calculate the heavy quarkonium helicity amplitudes in the framework of NRQCD factorization, which we dub HELAC-Onia. We rewrote the original HELAC to make the new program be able to calculate helicity amplitudes of multi P-wave quarkonium states production at hadron colliders and electron-positron colliders by including new P-wave off-shell currents. Therefore, besides the high efficiencies in computation of multi-leg processes within the Standard Model, HELAC-Onia is also sufficiently numerical stable in dealing with P-wave quarkonia (e.g. h,χ) and P-wave color-octet intermediate states. To the best of our knowledge, it is a first general-purpose automatic quarkonium matrix elements generator based on recursion relations on the market.
A modified Finite Element-Transfer Matrix for control design of space structures
NASA Technical Reports Server (NTRS)
Tan, T.-M.; Yousuff, A.; Bahar, L. Y.; Konstandinidis, M.
1990-01-01
The Finite Element-Transfer Matrix (FETM) method was developed for reducing the computational efforts involved in structural analysis. While being widely used by structural analysts, this method does, however, have certain limitations, particularly when used for the control design of large flexible structures. In this paper, a new formulation based on the FETM method is presented. The new method effectively overcomes the limitations in the original FETM method, and also allows an easy construction of reduced models that are tailored for the control design. Other advantages of this new method include the ability to extract open loop frequencies and mode shapes with less computation, and simplification of the design procedures for output feedback, constrained compensation, and decentralized control. The development of this new method and the procedures for generating reduced models using this method are described in detail and the role of the reduced models in control design is discussed through an illustrative example.
Mesh refinement in finite element analysis by minimization of the stiffness matrix trace
NASA Technical Reports Server (NTRS)
Kittur, Madan G.; Huston, Ronald L.
1989-01-01
Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.
Cwik, T.; Jamnejad, V.; Zuffada, C.
1994-12-31
The usefulness of finite element modeling follows from the ability to accurately simulate the geometry and three-dimensional fields on the scale of a fraction of a wavelength. To make this modeling practical for engineering design, it is necessary to integrate the stages of geometry modeling and mesh generation, numerical solution of the fields-a stage heavily dependent on the efficient use of a sparse matrix equation solver, and display of field information. The stages of geometry modeling, mesh generation, and field display are commonly completed using commercially available software packages. Algorithms for the numerical solution of the fields need to be written for the specific class of problems considered. Interior problems, i.e. simulating fields in waveguides and cavities, have been successfully solved using finite element methods. Exterior problems, i.e. simulating fields scattered or radiated from structures, are more difficult to model because of the need to numerically truncate the finite element mesh. To practically compute a solution to exterior problems, the domain must be truncated at some finite surface where the Sommerfeld radiation condition is enforced, either approximately or exactly. Approximate methods attempt to truncate the mesh using only local field information at each grid point, whereas exact methods are global, needing information from the entire mesh boundary. In this work, a method that couples three-dimensional finite element (FE) solutions interior to the bounding surface, with an efficient integral equation (IE) solution that exactly enforces the Sommerfeld radiation condition is developed. The bounding surface is taken to be a surface of revolution (SOR) to greatly reduce computational expense in the IE portion of the modeling.
Nuclear-Structure Data Relevant to Neutinoless-Double-Beta-Decay Matrix Elements
NASA Astrophysics Data System (ADS)
Kay, Benjamin
2015-10-01
An observation of neutrinoless double beta decay is one of the most exciting prospects in contemporary physics. It follows that calculations of the nuclear matrix elements for this process are of high priority. The change in the wave functions between the initial and final states of the neutrinoless-double-beta-decay candidates 76Ge-->76Se, 100Mo-->100Ru, 130Te-->130Xe, and 136Xe-->136Ba have been studied with transfer reactions. The data are focused on the change in the occupancies of the valence orbitals in the ground states as two neutrons decay into two protons. The results set a strict constraint on any theoretical calculations describing this rearrangement and thus on the magnitude of the nuclear matrix elements for this process, which currently exhibit uncertainties at the factor of 2-4 level. Prior to these measurements there were limited experimental data were available A = 76 and 100 systems, and very limited data for the A = 130 and 136 systems, in a large part due to the gaseous Xe isotopes involved. The uncertainties on most of these data are estimated to range from 0.1-0.3 nucleons. The program started with the A = 76 system, with subsequent calculations, modified to reproduce the experimental occupancies, exhibiting a significant reduction in the discrepancy between various models. New data are available for the A = 100 , 130, and 136 systems. I review the program, making detailed comparisons between the latest theoretical calculations and the experimental data where available. This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics, under Contract Number DE-AC02-06CH11357.
On-shell Delta I= 3/2 kaon weak matrix elements with nonzero total momentum
Yamazaki, T.
2009-05-20
We present our results for the on-shell {Delta}I = 3/2 kaon decay matrix elements using domain wall fermions and the DBW2 gauge action at one coarse lattice spacing corresponding to a{sup -1} = 1.31 GeV in the quenched approximation. The on-shell matrix elements are evaluated in two different frames: the center-of-mass frame and nonzero total-momentum frame. We employ the formula proposed by Lellouch and Luescher in the center-of-mass frame, and its extension for a nonzero total-momentum frame to extract the infinite volume, on-shell, center-of-mass frame decay amplitudes. We determine the decay amplitude at the physical pion mass and momentum from the chiral extrapolation and an interpolation of the relative momentum using the results calculated in the two frames. We have obtained ReA{sub 2} = 1.66(23)(+48/-03)(+53/-0) x 10{sup -8} GeV and ImA{sub 2} = -1.181(26)(+141/-014)(+44/-0) x 10{sup -12} GeV at the physical point, using the data at the relatively large pion mass, m{sub {pi}} > 0.35 GeV. The first error is statistic, and the second and third are systematic. The second error is estimated with several fits of the chiral extrapolation including the (quenched) chiral perturbation formula at next to leading order using only lighter pion masses. The third one is estimated with an analysis using the lattice dispersion relation. The result of ReA{sub 2} is reasonably consistent with experiment.
FRODO: a MuPAD program to calculate matrix elements between contracted wavefunctions
NASA Astrophysics Data System (ADS)
Angeli, C.; Cimiraglia, R.
2005-09-01
A symbolic program performing the Formal Reduction of Density Operators (FRODO) has been developed in the MuPAD computer algebra system with the purpose of evaluating the matrix elements of the electronic Hamiltonian between internally contracted functions in a complete active space (CAS) scheme. The program is illustrated making use of two meaningful examples. Program summaryTitle of program:FRODO Catalogue identifier:ADVY Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVY Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computer:Any computer on which the MuPAD computer algebra system can be installed Operating systems under which the program has been tested:Linux Programming language used:MuPAD vs. 2.5.3 for Linux No. of lines in distributed program, including test data, etc.:3939 No. of bytes in distributed program, including test data, etc.:19 661 Distribution format:tar.gz Nature of physical problem: In order to improve on the CAS-SCF wavefunction one can resort to multireference perturbation theory or configuration interaction based on internally contracted functions (ICF) which are obtained by application of the excitation operators to the reference CAS-SCF wavefunction. The formulation of such matrix elements is quite cumbersome and a computer algebra system like MuPAD appears ideally suited to perform such a task. Method of solution: The method adopted consists in successively eliminating all occurrences of inactive orbital indices (core and virtual) from the products of excitation operators which appear in the definition of the ICF's and in the electronic Hamiltonian expressed in the second quantization formalism. Restrictions due to the complexity of the problem: The program is limited to no more than doubly excited ICF's.
Spin Density Matrix Elements in Exclusive Production of Omega Mesons at HERMES
NASA Astrophysics Data System (ADS)
Marukyan, Hrachya
2016-02-01
Exclusive electroproduction of ω mesons on unpolarized hydrogen and deuterium targets is studied at HERMES in the kinematic region of Q2 > 1.0GeV2, 3.0GeV < W < 6.3GeV, and ‑ t‧ < 0.2GeV2. The data were accumulated during the 1996-2007 running period using the 27.6GeV longitudinally polarized electron or positron beams at HERA. The determination of the virtual-photon longitudinal-to-transverse cross-section ratio shows that a considerable part of the cross section arises from transversely polarized photons. Spin density matrix elements are derived and presented in projections of Q2 or ‑ t‧. Violation of s-channel helicity conservation is observed for some of these elements. A sizable contribution from unnatural-parity-exchange amplitudes is found and the phase shift between those amplitudes that describe transverse ω production by longitudinal and transverse virtual photons is determined for the first time. Good agreement is found between the HERMES proton data and results of a pQCD-inspired phenomenological model that includes pion-pole contributions.
Heavy-ion double charge exchange reactions: A tool toward 0 νββ nuclear matrix elements
NASA Astrophysics Data System (ADS)
Cappuzzello, F.; Cavallaro, M.; Agodi, C.; Bondì, M.; Carbone, D.; Cunsolo, A.; Foti, A.
2015-11-01
The knowledge of the nuclear matrix elements for the neutrinoless double beta decay is fundamental for neutrino physics. In this paper, an innovative technique to extract information on the nuclear matrix elements by measuring the cross section of a double charge exchange nuclear reaction is proposed. The basic point is that the initial- and final-state wave functions in the two processes are the same and the transition operators are similar. The double charge exchange cross sections can be factorized in a nuclear structure term containing the matrix elements and a nuclear reaction factor. First pioneering experimental results for the 40Ca(18O,18Ne)40Ar reaction at 270 MeV incident energy show that such cross section factorization reasonably holds for the crucial 0+ → 0+ transition to 40Args, at least at very forward angles.
Smith, F.G.; Wiederin, D.R.; Mortlock, R.
1994-12-31
Determination of the rare earth elements is important in the study of sedimentary processes. Geological and environmental samples often contain very low levels of these elements, and detection by plasma spectroscopy (ICP-AES, ICP-MS) is difficult unless a preconcentration and/or matrix elimination procedure is performed prior to analysis.; An automated batch preconcentration/matrix elimination system offers rapid, off-line sample preparation for a variety of sample types. A chelating form of a solid suspended reagent is added to a pH-adjusted sample. The suspended reagent with any bound elements are trapped in a hollow fiber membrane filter while unbound matrix components are washed to waste. The reagent with bound analytes are then released in a small volume. The system works in concert with an autosampler for unattended operation. Application to a variety of geological and environmental samples will be described.
0 ν β β and 2 ν β β nuclear matrix elements in the interacting boson model with isospin restoration
NASA Astrophysics Data System (ADS)
Barea, J.; Kotila, J.; Iachello, F.
2015-03-01
We introduce a method for isospin restoration in the calculation of nuclear matrix elements (NMEs) for 0 ν β β and 2 ν β β decay within the framework of the microscopic interacting boson model (IBM-2). With this method, we calculate the NMEs for all processes of interest in 0 ν β-β- and 2 ν β-β- and in 0 ν β+β+ , 0 ν EC β+ , R 0 ν ECEC , 2 ν β+β+ , 2 ν EC β+ , and 2 ν ECEC . With this method, the Fermi matrix elements for 2 ν β β vanish, and those for 0 ν β β are considerably reduced.
Matrix elements in the coupled-cluster approach - With application to low-lying states in Li
NASA Technical Reports Server (NTRS)
Martensson-Pendrill, Ann-Marie; Ynnerman, Anders
1990-01-01
A procedure is suggested for evaluating matrix elements of an operator between wavefunctions in the coupled-cluster form. The use of the exponential ansatz leads to compact exponential expressions also for matrix elements. Algorithms are developed for summing all effects of one-particle clusters and certain chains of two-particle clusters (containing the well-known random-phase approximation as a subset). The treatment of one-particle perturbations in single valence states is investigated in detail. As examples the oscillator strength for the 2s-2p transition in Li as well as the hyperfine structure for the two states are studied and compared to earlier work.
NASA Technical Reports Server (NTRS)
Sanfeliz, Jose G.
1993-01-01
Micromechanical modeling via elastic-plastic finite element analyses were performed to investigate the effects that the residual stresses and the degree of matrix work hardening (i.e., cold-worked, annealed) have upon the behavior of a 9 vol percent, unidirectional W/Cu composite, undergoing tensile loading. The inclusion of the residual stress-containing state as well as the simulated matrix material conditions proved to be significant since the Cu matrix material exhibited plastic deformation, which affected the subsequent tensile response of the composite system. The stresses generated during cooldown to room temperature from the manufacturing temperature were more of a factor on the annealed-matrix composite, since they induced the softened matrix to plastically flow. This event limited the total load-carrying capacity of this matrix-dominated, ductile-ductile type material system. Plastic deformation of the hardened-matrix composite during the thermal cooldown stage was not considerable, therefore, the composite was able to sustain a higher stress before showing any appreciable matrix plasticity. The predicted room temperature, stress-strain response, and deformation stages under both material conditions represented upper and lower bounds characteristic of the composite's tensile behavior. The initial deformation stage for the hardened material condition showed negligible matrix plastic deformation while for the annealed state, its initial deformation stage showed extensive matrix plasticity. Both material conditions exhibited a final deformation stage where the fiber and matrix were straining plastically. The predicted stress-strain results were compared to the experimental, room temperature, tensile stress-strain curve generated from this particular composite system. The analyses indicated that the actual thermal-mechanical state of the composite's Cu matrix, represented by the experimental data, followed the annealed material condition.
LATTICE MATRIX ELEMENTS AND CP VIOLATION IN B AND KA PHYSICS: STATUS AND OUTLOOK.
SONI,A.
2003-01-03
Status of lattice calculations of hadron matrix elements along with CP violation in B and in K systems is reviewed. Lattice has provided useful input which, in conjunction with experimental data, leads to the conclusion that CP-odd phase in the CKM matrix plays the dominant role in the observed asymmetry in B {yields} {psi}K{sub s}. It is now quite likely that any beyond the SM, CP-odd, phase will cause only small deviations in B-physics. Search for the effects of the new phase(s) will consequently require very large data samples as well as very precise theoretical predictions. Clean determination of all the angles of the unitarity triangle therefore becomes essential. In this regard B {yields} KD{sup 0} processes play a unique role. Regarding K-decays, remarkable progress made by theory with regard to maintenance of chiral symmetry on the lattice is briefly discussed. First application already provide quantitative information on B{sub K} and the {Delta}I = 1/2 rule. The enhancement in ReA{sub 0} appears to arise solely from tree operators, esp. Q{sub 2}; penguin contribution to ReA{sub 0} appears to be very small. However, improved calculations are necessary for {epsilon}{prime}/{epsilon} as there the contributions of QCD penguins and electroweak penguins largely seem to cancel. There are good reasons, though, to believe that these cancellations will not survive improvements that are now underway. Importance of determining the unitarity triangle purely from K-decays is also emphasized.
NASA Astrophysics Data System (ADS)
Ebel, Denton S.; Brunner, Chelsea; Konrad, Kevin; Leftwich, Kristin; Erb, Isabelle; Lu, Muzhou; Rodriguez, Hugo; Crapster-Pregont, Ellen J.; Friedrich, Jon M.; Weisberg, Michael K.
2016-01-01
The relative abundances and chemical compositions of the macroscopic components or "inclusions" (chondrules and refractory inclusions) and fine-grained mineral matrix in chondritic meteorites provide constraints on astrophysical theories of inclusion formation and chondrite accretion. We present new techniques for analysis of low count/pixel Si, Mg, Ca, Al, Ti and Fe X-ray intensity maps of rock sections, and apply them to large areas of CO and CV chondrites, and the ungrouped Acfer 094 chondrite. For many thousands of manually segmented and type-identified inclusions, we are able to assess, pixel-by-pixel, the major element content of each inclusion. We quantify the total fraction of refractory elements accounted for by various types of inclusion and matrix. Among CO chondrites, both matrix and inclusion Mg/Si ratios approach the solar (and bulk CO) ratio with increasing petrologic grade, but Si remains enriched in inclusions relative to matrix. The oxidized CV chondrites with higher matrix/inclusion ratios exhibit more severe aqueous alteration (oxidation), and their excess matrix accounts for their higher porosity relative to reduced CV chondrites. Porosity could accommodate an original ice component of matrix as the direct cause of local alteration of oxidized CV chondrites. We confirm that major element abundances among inclusions differ greatly, across a wide range of CO and CV chondrites. These abundances in all cases add up to near-chondritic (solar) bulk abundance ratios in these chondrites, despite wide variations in matrix/inclusion ratios and inclusion sizes: chondrite components are complementary. This complementarity provides a robust meteoritic constraint for astrophysical disk models.
Bubin, Sergiy; Adamowicz, Ludwik
2008-03-21
In this work we consider explicitly correlated complex Gaussian basis functions for expanding the wave function of an N-particle system with the L=1 total orbital angular momentum. We derive analytical expressions for various matrix elements with these basis functions including the overlap, kinetic energy, and potential energy (Coulomb interaction) matrix elements, as well as matrix elements of other quantities. The derivatives of the overlap, kinetic, and potential energy integrals with respect to the Gaussian exponential parameters are also derived and used to calculate the energy gradient. All the derivations are performed using the formalism of the matrix differential calculus that facilitates a way of expressing the integrals in an elegant matrix form, which is convenient for the theoretical analysis and the computer implementation. The new method is tested in calculations of two systems: the lowest P state of the beryllium atom and the bound P state of the positronium molecule (with the negative parity). Both calculations yielded new, lowest-to-date, variational upper bounds, while the number of basis functions used was significantly smaller than in previous studies. It was possible to accomplish this due to the use of the analytic energy gradient in the minimization of the variational energy. PMID:18361554
NASA Astrophysics Data System (ADS)
Angeli, C.; Cimiraglia, R.
2013-02-01
A symbolic program performing the Formal Reduction of Density Operators (FRODO), formerly developed in the MuPAD computer algebra system with the purpose of evaluating the matrix elements of the electronic Hamiltonian between internally contracted functions in a complete active space (CAS) scheme, has been rewritten in Mathematica. New version : A program summaryProgram title: FRODO Catalogue identifier: ADV Y _v2_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVY_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 3878 No. of bytes in distributed program, including test data, etc.: 170729 Distribution format: tar.gz Programming language: Mathematica Computer: Any computer on which the Mathematica computer algebra system can be installed Operating system: Linux Classification: 5 Catalogue identifier of previous version: ADV Y _v1_0 Journal reference of previous version: Comput. Phys. Comm. 171(2005)63 Does the new version supersede the previous version?: No Nature of problem. In order to improve on the CAS-SCF wavefunction one can resort to multireference perturbation theory or configuration interaction based on internally contracted functions (ICFs) which are obtained by application of the excitation operators to the reference CAS-SCF wavefunction. The previous formulation of such matrix elements in the MuPAD computer algebra system, has been rewritten using Mathematica. Solution method: The method adopted consists in successively eliminating all occurrences of inactive orbital indices (core and virtual) from the products of excitation operators which appear in the definition of the ICFs and in the electronic Hamiltonian expressed in the second quantization formalism. Reasons for new version: Some years ago we published in this journal a couple of papers [1, 2
Symmetry of Isoscalar Matrix Elements and Systematics in the sd and beginning of fp shells
NASA Astrophysics Data System (ADS)
Orce, J. N.; Petkov, P.; Velázquez, V.; McKay, C. J.; Lesher, S. R.; Choudry, S.; Mynk, M.; Linnemann, A.; Jolie, J.; von Brentano, P.; Werner, V.; Yates, S. W.; McEllistrem, M. T.
2006-03-01
A careful determination of the lifetime and measurement of the branching ratio for decay of the first 2T=1+ state in 42Sc has allowed an accurate experimental test of charge independence in the A = 42 isobaric triplet. A lifetime of 69(17) fs was measured at the University of Kentucky, while relative intensities for the 975 keV and 1586 keV transitions depopulating the first 2T=1+ state have been determined at the University of Cologne as 100(1) and 8(1), respectively. Both measurements give an isoscalar matrix element, M0, of 6.4(9) (W.u.)1/2. This result confirms charge independence for the A=42 isobaric triplet. Shell model calculations have been carried out for understanding the global trend of M0 values for A = 4n + 2 isobaric triplets ranging from A = 18 to A = 42. The 21 (T=1)+ → 01 (T=1)+ transition energies, reduced transition probabilities and M0 values are reproduced to a high degree of accuracy. The trend of M0 strength along the sd shell is interpreted in terms of the shell structure. Certain discrepancies arise at the extremes of the sd shell, for the A = 18 and A = 38 isobaric triplets, which might be explained in terms of the low valence space at the extremes of the sd shell.
Electron scattering from large molecules: a 3d finite element R-matrix approach
NASA Astrophysics Data System (ADS)
Tonzani, Stefano; Greene, Chris H.
2005-05-01
To solve the Schr"odinger equation for scattering of a low energy electron from a molecule, we present a three-dimensional finite element R-matrix method [S. Tonzani and C. H. Greene, J. Chem. Phys. 122 01411, (2005)]. Using the static exchange and local density approximations, we can use directly the molecular potentials extracted from ab initio codes (GAUSSIAN 98 in the work described here). A local polarization potential based on density functional theory [F. A. Gianturco and A. Rodriguez-Ruiz, Phys. Rev. A 47, 1075 (1993)] approximately describes the long range attraction to the molecular target induced by the scattering electron without adjustable parameters. We have used this approach successfully in calculations of cross sections for small and medium sized molecules (like SF6, XeF6, C60 and Uracil). This method will be useful to treat the electron-induced dynamics of extended molecular systems, possibly of biological interest, where oth er more complex ab initio methods are difficult to apply.
NASA Astrophysics Data System (ADS)
Fatchurrohman, N.; Marini, C. D.; Suraya, S.; Iqbal, AKM Asif
2016-02-01
The increasing demand of fuel efficiency and light weight components in automobile sectors have led to the development of advanced material parts with improved performance. A specific class of MMCs which has gained a lot of attention due to its potential is aluminium metal matrix composites (Al-MMCs). Product performance investigation of Al- MMCs is presented in this article, where an Al-MMCs brake disc is analyzed using finite element analysis. The objective is to identify the potentiality of replacing the conventional iron brake disc with Al-MMCs brake disc. The simulation results suggested that the MMCs brake disc provided better thermal and mechanical performance as compared to the conventional cast iron brake disc. Although, the Al-MMCs brake disc dissipated higher maximum temperature compared to cast iron brake disc's maximum temperature. The Al-MMCs brake disc showed a well distributed temperature than the cast iron brake disc. The high temperature developed at the ring of the disc and heat was dissipated in circumferential direction. Moreover, better thermal dissipation and conduction at brake disc rotor surface played a major influence on the stress. As a comparison, the maximum stress and strain of Al-MMCs brake disc was lower than that induced on the cast iron brake disc.
NASA Astrophysics Data System (ADS)
Ablinger, J.; Behring, A.; Blümlein, J.; De Freitas, A.; von Manteuffel, A.; Schneider, C.
2016-05-01
Three loop ladder and V-topology diagrams contributing to the massive operator matrix element AQg are calculated. The corresponding objects can all be expressed in terms of nested sums and recurrences depending on the Mellin variable N and the dimensional parameter ε. Given these representations, the desired Laurent series expansions in ε can be obtained with the help of our computer algebra toolbox. Here we rely on generalized hypergeometric functions and Mellin-Barnes representations, on difference ring algorithms for symbolic summation, on an optimized version of the multivariate Almkvist-Zeilberger algorithm for symbolic integration, and on new methods to calculate Laurent series solutions of coupled systems of differential equations. The solutions can be computed for general coefficient matrices directly for any basis also performing the expansion in the dimensional parameter in case it is expressible in terms of indefinite nested product-sum expressions. This structural result is based on new results of our difference ring theory. In the cases discussed we deal with iterative sum- and integral-solutions over general alphabets. The final results are expressed in terms of special sums, forming quasi-shuffle algebras, such as nested harmonic sums, generalized harmonic sums, and nested binomially weighted (cyclotomic) sums. Analytic continuations to complex values of N are possible through the recursion relations obeyed by these quantities and their analytic asymptotic expansions. The latter lead to a host of new constants beyond the multiple zeta values, the infinite generalized harmonic and cyclotomic sums in the case of V-topologies.
Characterization of metal matrix composites by linear ultrasonics and finite element modeling.
Chen, Xuesheng; Sharples, Steve D; Clark, Matt; Wright, David
2013-02-01
Titanium metal matrix composites (TiMMCs) offer advantages over traditional materials for aerospace applications due to the increased mechanical strength of the materials. But the non-destructive inspection of these materials, especially with ultrasound, is in an infancy stage. If the manufacturing process of TiMMC is not correctly controlled, then disbonds and voids between the fibers can result. The effective microstructure of the composite makes difficulty to interpret results from traditional ultrasound techniques because of the scattering caused by fibers; the scattering prevents the ultrasound from penetrating far into the composite region and produces a background signal masking any reflections from voids. In this paper, relatively low frequency ultrasound is used to probe the composite region, and the state of the composite (porosity) is inferred from the velocity of the ultrasound traversing the composite. The relationship between the velocity and porosity is complex in this regime, so finite element (FE) analysis is used to model the composite regions and relate the velocity to the porosity. The FE simulated results are validated by ultrasound velocity measurements. PMID:23363095
Menéndez, Javier
2013-12-30
We explore the theoretical uncertainties related to the transition operator of neutrinoless double-beta (0νββ) decay. The transition operator used in standard calculations is a product of one-body currents, that can be obtained phenomenologically as in Tomoda [1] or Šimkovic et al. [2]. However, corrections to the operator are hard to obtain in the phenomenological approach. Instead, we calculate the 0νββ decay operator in the framework of chiral effective theory (EFT), which gives a systematic order-by-order expansion of the transition currents. At leading orders in chiral EFT we reproduce the standard one-body currents of Refs. [1] and [2]. Corrections appear as two-body (2b) currents predicted by chiral EFT. We compute the effects of the leading 2b currents to the nuclear matrix elements of 0νββ decay for several transition candidates. The 2b current contributions are related to the quenching of Gamow-Teller transitions found in nuclear structure calculations.
A measurement of the top quark mass with a matrix element method
Gibson, Adam Paul; /UC, Berkeley
2006-12-01
The authors present a measurement of the mass of the top quark. The event sample is selected from proton-antiproton collisions, at 1.96 TeV center-of-mass energy, observed with the CDF detector at Fermilab's Tevatron. They consider a 318 pb{sup -1} dataset collected between March 2002 and August 2004. They select events that contain one energetic lepton, large missing transverse energy, exactly four energetic jets, and at least one displaced vertex b tag. The analysis uses leading-order t{bar t} and background matrix elements along with parameterized parton showering to construct event-by-event likelihoods as a function of top quark mass. From the 63 events observed with the 318 pb{sup -1} dataset they extract a top quark mass of 172.0 {+-} 2.6(stat) {+-} 3.3(syst) GeV/c{sup 2} from the joint likelihood. The mean expected statistical uncertainty is 3.2 GeV/c{sup 2} for m{sub t} = 178 GTeV/c{sup 2} and 3.1 GeV/c{sup 2} for m{sub t} = 172.5 GeV/c{sup 2}. The systematic error is dominated by the uncertainty of the jet energy scale.
NASA Astrophysics Data System (ADS)
Angeli, C.; Cimiraglia, R.
2013-02-01
A symbolic program performing the Formal Reduction of Density Operators (FRODO), formerly developed in the MuPAD computer algebra system with the purpose of evaluating the matrix elements of the electronic Hamiltonian between internally contracted functions in a complete active space (CAS) scheme, has been rewritten in Mathematica. New version : A program summaryProgram title: FRODO Catalogue identifier: ADV Y _v2_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVY_v2_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 3878 No. of bytes in distributed program, including test data, etc.: 170729 Distribution format: tar.gz Programming language: Mathematica Computer: Any computer on which the Mathematica computer algebra system can be installed Operating system: Linux Classification: 5 Catalogue identifier of previous version: ADV Y _v1_0 Journal reference of previous version: Comput. Phys. Comm. 171(2005)63 Does the new version supersede the previous version?: No Nature of problem. In order to improve on the CAS-SCF wavefunction one can resort to multireference perturbation theory or configuration interaction based on internally contracted functions (ICFs) which are obtained by application of the excitation operators to the reference CAS-SCF wavefunction. The previous formulation of such matrix elements in the MuPAD computer algebra system, has been rewritten using Mathematica. Solution method: The method adopted consists in successively eliminating all occurrences of inactive orbital indices (core and virtual) from the products of excitation operators which appear in the definition of the ICFs and in the electronic Hamiltonian expressed in the second quantization formalism. Reasons for new version: Some years ago we published in this journal a couple of papers [1, 2
Symmetry of Isoscalar Matrix Elements and Systematics in the sd and beginning of fp shells
Orce, J. N.; McKay, C. J.; Lesher, S. R.; Choudry, S.; Mynk, M.; McEllistrem, M. T.; Petkov, P.; Velazquez, V.; Linnemann, A.; Jolie, J.; Brentano, P. von; Werner, V.; Yates, S. W.
2006-03-13
A careful determination of the lifetime and measurement of the branching ratio for decay of the first 2{sub T=1}{sup +} state in 42Sc has allowed an accurate experimental test of charge independence in the A = 42 isobaric triplet. A lifetime of 69(17) fs was measured at the University of Kentucky, while relative intensities for the 975 keV and 1586 keV transitions depopulating the first 2{sub T=1}{sup +} state have been determined at the University of Cologne as 100(1) and 8(1), respectively. Both measurements give an isoscalar matrix element, M0, of 6.4(9) (W.u.)1/2. This result confirms charge independence for the A=42 isobaric triplet. Shell model calculations have been carried out for understanding the global trend of M0 values for A = 4n + 2 isobaric triplets ranging from A = 18 to A = 42. The 2{sub 1(T=1)}{sup +} {yields} 0{sub 1(T=1)}{sup +} transition energies, reduced transition probabilities and M0 values are reproduced to a high degree of accuracy. The trend of M0 strength along the sd shell is interpreted in terms of the shell structure. Certain discrepancies arise at the extremes of the sd shell, for the A = 18 and A 38 isobaric triplets, which might be explained in terms of the low valence space at the extremes of the sd shell.
Bijker, R.; Dieperink, A.E.L.
1982-12-01
It is shown that the dominance of ..beta --> gamma.. and ..gamma -->..g over ..beta -->..gE2 transitions in the SU(3) limit of the interacting-boson-approximation model, reported by Warner and Casten can be explained simply in terms of properties of the intrinsic E2 matrix elements.
NASA Astrophysics Data System (ADS)
Gallup, G. A.
1986-07-01
In a recent article [Phys. Rev. A 31, 2107 (1985)] Leasure and Balint-Kurti claim to give a more efficient algorithm than any previously available for determining matrix elements of the Hamiltonian in valence-bond calculations. Actually, an algorithm of no significant difference and the same efficiency has been available since 1972 and has been applied to valence-bond calculations.
Hollaus, K; Magele, C; Merwa, R; Scharfetter, H
2004-02-01
Magnetic induction tomography of biological tissue is used to reconstruct the changes in the complex conductivity distribution by measuring the perturbation of an alternating primary magnetic field. To facilitate the sensitivity analysis and the solution of the inverse problem a fast calculation of the sensitivity matrix, i.e. the Jacobian matrix, which maps the changes of the conductivity distribution onto the changes of the voltage induced in a receiver coil, is needed. The use of finite differences to determine the entries of the sensitivity matrix does not represent a feasible solution because of the high computational costs of the basic eddy current problem. Therefore, the reciprocity theorem was exploited. The basic eddy current problem was simulated by the finite element method using symmetric tetrahedral edge elements of second order. To test the method various simulations were carried out and discussed. PMID:15005313
NASA Technical Reports Server (NTRS)
Jenkins, J. M.; Taylor, A. H.; Sakata, I. F.
1985-01-01
A hybrid spar of titanium with an integrally brazed composite, consisting of an aluminum matrix reinforced with boron-carbide-coated fibers, was heated in an oven and the resulting thermal stresses were measured. Uniform heating of the spar in an oven resulted in thermal stresses arising from the effects of dissimilar materials and anisotropy of the metal matrix composite. Thermal stresses were calculated from a finite element structural model using anisotropic material properties deduced from constituent properties and rules of mixtures. Comparisons of calculated thermal stresses with measured thermal stresses on the spar are presented. It was shown that failure to account for anisotropy in the metal matrix composite elements would result in large errors in correlating measured and calculated thermal stresses. It was concluded that very strong material characterization efforts are required to predict accurate thermal stresses in anisotropic composite structures.
NASA Astrophysics Data System (ADS)
Hinlein, E. S.; Ostendorf, D. W.
2010-12-01
A matrix of mini-wells was installed and instrumented to monitor the hyporheic zone linking a river and its floodplain deposit in Eastern Massachusetts. The Neponset River is of small to moderate size with seasonal approximate depths ranging from 1-3 m, widths of 5-15 m, and flows varying from 0.3-28 m3/s in the study area adjacent to a major interstate highway. Transport in the hyporheic zone is governed by the floodplain deposit and the river which combine to induce diurnal timescale dispersive mixing. A low steady groundwater gradient drives transport toward the river at a rate of approximately 10 meters per year. The floodplain deposit is made up of silty sand to a depth of 21 m underlain by 9 m of low permeability medium silt. The mini-well installation is in the medium silt river bank. Because of the presence of conductive ions in the floodplain deposit, it is possible to use the mini-well matrix to document a front where groundwater and river water meet both through changes in hydraulic head and groundwater conductivity. Specific conductivity values are in the range of 500 uS/cm in the river and 1500 uS/cm in the adjacent 4.5 m deep well 15 m away. Four clusters of existing monitoring wells currently measure head and conductivity outside the hyporheic zone in 4.5 m deep wells along a line from 15-300 m away from the river. A monitoring station at the river measures river level and conductivity. The mini-well matrix captures the final 15 m between the last well cluster and the river. Data from the existing well cluster adjacent to the hyporheic zone shows response to high river levels on the order of 12-24 hrs, with a predicted excursion amplitude of 1 m in response to a hydrograph amplitude of 2 m from a storm of Feb, 2010. The mini-well matrix will measure the horizontal excursions as well as vertical gradients of head and conductivity. Due to complications associated with river bank topography and the sensitivity of the area, the mini-wells were installed
Birkholzer, J.; Karasaki, K.
1996-07-01
Fracture network simulators have extensively been used in the past for obtaining a better understanding of flow and transport processes in fractured rock. However, most of these models do not account for fluid or solute exchange between the fractures and the porous matrix, although diffusion into the matrix pores can have a major impact on the spreading of contaminants. In the present paper a new finite element code TRIPOLY is introduced which combines a powerful fracture network simulator with an efficient method to account for the diffusive interaction between the fractures and the adjacent matrix blocks. The fracture network simulator used in TRIPOLY features a mixed Lagrangian-Eulerian solution scheme for the transport in fractures, combined with an adaptive gridding technique to account for sharp concentration fronts. The fracture-matrix interaction is calculated with an efficient method which has been successfully used in the past for dual-porosity models. Discrete fractures and matrix blocks are treated as two different systems, and the interaction is modeled by introducing sink/source terms in both systems. It is assumed that diffusive transport in the matrix can be approximated as a one-dimensional process, perpendicular to the adjacent fracture surfaces. A direct solution scheme is employed to solve the coupled fracture and matrix equations. The newly developed combination of the fracture network simulator and the fracture-matrix interaction module allows for detailed studies of spreading processes in fractured porous rock. The authors present a sample application which demonstrate the codes ability of handling large-scale fracture-matrix systems comprising individual fractures and matrix blocks of arbitrary size and shape.
Yeh, Rong-Guan; Lin, Chung-Wu; Abbod, Maysam F.; Shieh, Jiann-Shing
2012-01-01
A detrended fluctuation analysis (DFA) method is applied to image analysis. The 2-dimensional (2D) DFA algorithms is proposed for recharacterizing images of lymph sections. Due to Burkitt lymphoma (BL) and diffuse large B-cell lymphoma (DLBCL), there is a significant different 5-year survival rates after multiagent chemotherapy. Therefore, distinguishing the difference between BL and DLBCL is very important. In this study, eighteen BL images were classified as group A, which have one to five cytogenetic changes. Ten BL images were classified as group B, which have more than five cytogenetic changes. Both groups A and B BLs are aggressive lymphomas, which grow very fast and require more intensive chemotherapy. Finally, ten DLBCL images were classified as group C. The short-term correlation exponent α1 values of DFA of groups A, B, and C were 0.370 ± 0.033, 0.382 ± 0.022, and 0.435 ± 0.053, respectively. It was found that α1 value of BL image was significantly lower (P < 0.05) than DLBCL. However, there is no difference between the groups A and B BLs. Hence, it can be concluded that α1 value based on DFA statistics concept can clearly distinguish BL and DLBCL image. PMID:23365623
NASA Astrophysics Data System (ADS)
Sturtz, Timothy M.
Source apportionment models attempt to untangle the relationship between pollution sources and the impacts at downwind receptors. Two frameworks of source apportionment models exist: source-oriented and receptor-oriented. Source based apportionment models use presumed emissions and atmospheric processes to estimate the downwind source contributions. Conversely, receptor based models leverage speciated concentration data from downwind receptors and apply statistical methods to predict source contributions. Integration of both source-oriented and receptor-oriented models could lead to a better understanding of the implications sources have on the environment and society. The research presented here investigated three different types of constraints applied to the Positive Matrix Factorization (PMF) receptor model within the framework of the Multilinear Engine (ME-2): element ratio constraints, spatial separation constraints, and chemical transport model (CTM) source attribution constraints. PM10-2.5 mass and trace element concentrations were measured in Winston-Salem, Chicago, and St. Paul at up to 60 sites per city during two different seasons in 2010. PMF was used to explore the underlying sources of variability. Information on previously reported PM10-2.5 tire and brake wear profiles were used to constrain these features in PMF by prior specification of selected species ratios. We also modified PMF to allow for combining the measurements from all three cities into a single model while preserving city-specific soil features. Relatively minor differences were observed between model predictions with and without the prior ratio constraints, increasing confidence in our ability to identify separate brake wear and tire wear features. Using separate data, source contributions to total fine particle carbon predicted by a CTM were incorporated into the PMF receptor model to form a receptor-oriented hybrid model. The level of influence of the CTM versus traditional PMF was
Matrix elements of the electromagnetic operator between kaon and pion states
Baum, I.; Lubicz, V.; Martinelli, G.; Orifici, L.; Simula, S.
2011-10-01
We compute the matrix elements of the electromagnetic operator sF{sub {mu}{nu}}{sigma}{sup {mu}{nu}}d between kaon and pion states, using lattice QCD with maximally twisted-mass fermions and two flavors of dynamical quarks (N{sub f}=2). The operator is renormalized nonperturbatively in the RI'/MOM scheme and our simulations cover pion masses as light as 270 MeV and three values of the lattice spacing from {approx_equal}0.07 up to {approx_equal}0.1 fm. At the physical point our result for the corresponding tensor form factor at zero-momentum transfer is f{sub T}{sup K{pi}}(0)=0.417(14{sub stat})(5{sub syst}), where the systematic error does not include the effect of quenching the strange and charm quarks. Our result differs significantly from the old quenched result f{sub T}{sup K{pi}}(0)=0.78(6) obtained by the SPQ{sub cd}R Collaboration with pion masses above 500 MeV. We investigate the source of this difference and conclude that it is mainly related to the chiral extrapolation. We also study the tensor charge of the pion and obtain the value f{sub T}{sup {pi}{pi}}(0)=0.195(8{sub stat})(6{sub syst}) in good agreement with, but more accurate than the result f{sub T}{sup {pi}{pi}}(0)=0.216(34) obtained by the QCDSF Collaboration using higher pion masses.
A finite element model of the effects of primary creep in an Al-SiC metal matrix composite
NASA Astrophysics Data System (ADS)
Atkins, Steven L.; Gibeling, Jeffery C.
1995-12-01
A two dimensional axisymmetric finite element model has been developed to study the creep behavior of a high-temperature aluminum alloy matrix (alloy 8009) reinforced with 11 vol pct silicon carbide paniculate. Because primary creep represents a significant portion of the total creep strain for this matrix alloy, the emphasis of the present investigation is on the influence of primary creep on the high-temperature behavior of the composite. The base alloy and composite are prepared by rapid solidification processing, resulting in a very fine grain size and the absence of precipitates that may complicate modeling of the composite. Because the matrix microstructure is unaffected by the presence of the SiC paniculate, this material is particularly well suited to continuum finite element modeling. Stress contours, strain contours, and creep curves are presented for the model. While the final distribution of stresses and strains is unaffected by the inclusion of primary creep, the overall creep response of the model reveals a significant primary strain transient. The effects of true primary creep are more significant than the primary-like transient introduced by the redistribution of stresses after loading. Examination of the stress contours indicates that the matrix axial and shear components become less uniform while the effective stress becomes more homogeneous as creep progresses and that the distribution of stresses do not change significantly with time after the strain rate reaches a steady state. These results also confirm that load transfer from the matrix to reinforcement occurs primarily through the shear stress. It is concluded that inclusion of matrix primary creep is essential to obtaining accurate representations of the creep response of metal matrix composites.
Mahfuz, H.; Ahsan Mian, A.K.M.; Vaidya, U.K.; Brown, T.; Jeelani, S.
1995-10-01
A 3D-unit cell for 0.90 laminated composites has been developed to predict the composite behavior under longitudinal tensile loading condition. 3D contact element has been used to model the fiber matrix interface. Two interface conditions, namely, infinitely strong and weakly bonded, are considered in the analysis. Both large displacement and plastic strain behavior for the matrix are considered to account for the geometric and material non-linearities. Investigations were carried out at three temperatures to compare the composite response obtained from mechanical tests at those temperatures. Stress-strain behavior and the local stress distributions at the fiber as well as at the matrix are presented, and their effects on the failure of the interface are discussed in the paper. The material under investigation was SiC{sub f}/Si{sub 3}N{sub 4}.
NASA Astrophysics Data System (ADS)
Whiting, Daniel J.; Keaveney, James; Adams, Charles S.; Hughes, Ifan G.
2016-04-01
Applying large magnetic fields to gain access to the hyperfine Paschen-Back regime can isolate three-level systems in a hot alkali metal vapors, thereby simplifying usually complex atom-light interactions. We use this method to make the first direct measurement of the |<5 P ||e r ||5 D >| matrix element in 87Rb. An analytic model with only three levels accurately models the experimental electromagnetically induced transparency spectra and extracted Rabi frequencies are used to determine the dipole matrix element. We measure |<5 P3 /2||e r ||5 D5 /2>| =(2.290 ±0 .002stat±0 .04syst) e a0 , which is in excellent agreement with the theoretical calculations of Safronova, Williams, and Clark [Phys. Rev. A 69, 022509 (2004), 10.1103/PhysRevA.69.022509].
NASA Astrophysics Data System (ADS)
Demić, Aleksandar; Radovanović, Jelena; Milanović, Vitomir
2016-08-01
We present a method for modeling nonparabolicity effects (NPE) in quantum nanostructures in presence of external electric and magnetic field by using second order perturbation theory. The method is applied to analysis of quantum well structure and active region of a quantum cascade laser (QCL). This model will allow us to examine the influence of magnetic field on dipole matrix element in QCL structures, which will provide a better insight to how NPE can affect the gain of QCL structures.
0{nu}{beta}{beta}-decay nuclear matrix elements with self-consistent short-range correlations
Simkovic, Fedor; Faessler, Amand; Muether, Herbert; Rodin, Vadim; Stauf, Markus
2009-05-15
A self-consistent calculation of nuclear matrix elements of the neutrinoless double-beta decays (0{nu}{beta}{beta}) of {sup 76}Ge, {sup 82}Se, {sup 96}Zr, {sup 100}Mo, {sup 116}Cd, {sup 128}Te, {sup 130}Te, and {sup 136}Xe is presented in the framework of the renormalized quasiparticle random phase approximation (RQRPA) and the standard QRPA. The pairing and residual interactions as well as the two-nucleon short-range correlations are for the first time derived from the same modern realistic nucleon-nucleon potentials, namely, from the charge-dependent Bonn potential (CD-Bonn) and the Argonne V18 potential. In a comparison with the traditional approach of using the Miller-Spencer Jastrow correlations, matrix elements for the 0{nu}{beta}{beta} decay are obtained that are larger in magnitude. We analyze the differences among various two-nucleon correlations including those of the unitary correlation operator method (UCOM) and quantify the uncertainties in the calculated 0{nu}{beta}{beta}-decay matrix elements.
Energy levels and transition probability matrix elements of ruby for maser applications
NASA Technical Reports Server (NTRS)
Berwin, R. W.
1971-01-01
Program computes fine structure energy levels of ruby as a function of magnetic field. Included in program is matrix formulation, each row of which contains a magnetic field and four corresponding energy levels.
NASA Astrophysics Data System (ADS)
Franco, E.; Maiani, L.; Martinelli, G.; Morelli, A.
1988-10-01
The K-pi and K-pi-pi elements of left-right four fermion operators in quenched lattice QCD at beta=6 are computed. The soft-pion relations derived from the chiral structure of the operators are checked. A large enhancement of matrix elements is observed and interpreted as the effect of a scalar octet pole in the pi-K channel. This observation has implications for the related calculation of weak matrix elements.
Shekhar, R.; Karunasagar, D.; Ranjit, M.; Arunachalam, J.
2009-10-15
An open-to-air type electrolyte cathode discharge (ELCAD) has been developed with a new design. The present configuration leads to a stable plasma even at low flow rates (0.96 mL/min). Plasma fluctuations arising from the variations in the gap between solid anode and liquid cathode were eliminated by providing a V-groove to the liquid glass-capillary. Cathode (ground) connection is given to the solution at the V-groove itself. Interfaced to atomic emission spectrometry (AES), its analytical performance is evaluated. The optimized molarity of the solution is 0.2 M. The analytical response curves for Ca, Cu, Cd, Pb, Hg, Fe, and Zn demonstrated good linearity. The limit of detections of Ca, Cu, Cd, Pb, Hg, Fe, and Zn are determined to be 17, 11, 5, 45, 15, 28, and 3 ng mL{sup -1}. At an integration time of 0.3 s, the relative standard deviation (RSD) values of the acid blank solutions are found to be less than 10% for the elements Ca, Cu, Cd, Hg, Fe, and Zn and 18% for Pb. The method is applied for the determination of the elemental constituents in different matrix materials such as tuna fish (IAEA-350), oyster tissue (NIST SRM 1566a), and coal fly ash (CFA SRM 1633b). The obtained results are in good agreement with the certified values. The accuracy is found to be between 7% and 0.6% for major to trace levels of constituent elements and the precision between 11% and 0.6%. For the injection of 100 {mu} L of 200 ng mL{sup -1} mercury solution at the flow rate of 0.8 mL/min, the flow injection studies resulted in the relative standard deviation (RSD) of 8%, concentration detection limit of 10 ng/mL, and mass detection limit of 1 ng for mercury.
Shekhar, R; Karunasagar, D; Ranjit, Manjusha; Arunachalam, J
2009-10-01
An open-to-air type electrolyte cathode discharge (ELCAD) has been developed with a new design. The present configuration leads to a stable plasma even at low flow rates (0.96 mL/min). Plasma fluctuations arising from the variations in the gap between solid anode and liquid cathode were eliminated by providing a V-groove to the liquid glass-capillary. Cathode (ground) connection is given to the solution at the V-groove itself. Interfaced to atomic emission spectrometry (AES), its analytical performance is evaluated. The optimized molarity of the solution is 0.2 M. The analytical response curves for Ca, Cu, Cd, Pb, Hg, Fe, and Zn demonstrated good linearity. The limit of detections of Ca, Cu, Cd, Pb, Hg, Fe, and Zn are determined to be 17, 11, 5, 45, 15, 28, and 3 ng mL(-1). At an integration time of 0.3 s, the relative standard deviation (RSD) values of the acid blank solutions are found to be less than 10% for the elements Ca, Cu, Cd, Hg, Fe, and Zn and 18% for Pb. The method is applied for the determination of the elemental constituents in different matrix materials such as tuna fish (IAEA-350), oyster tissue (NIST SRM 1566a), and coal fly ash (CFA SRM 1633b). The obtained results are in good agreement with the certified values. The accuracy is found to be between 7% and 0.6% for major to trace levels of constituent elements and the precision between 11% and 0.6%. For the injection of 100 microL of 200 ng mL(-1) mercury solution at the flow rate of 0.8 mL/min, the flow injection studies resulted in the relative standard deviation (RSD) of 8%, concentration detection limit of 10 ng/mL, and mass detection limit of 1 ng for mercury. PMID:19715301
Lloyd, J W; Rook, J S; Braselton, E; Shea, M E
2000-02-01
A study was designed to model the fluctuations of nine specific element concentrations in mammary secretions from periparturient mares over time. During the 1992 foaling season, serial samples of mammary secretions were collected from all 18 pregnant Arabian mares at the Michigan State University equine teaching and research center. Non-linear regression techniques were used to model the relationship between element concentration in mammary secretions and days from foaling (which connected two separate sigmoid curves with a spline function); indicator variables were included for mare and mare parity. Element concentrations in mammary secretions varied significantly during the periparturient period in mares. Both time trends and individual variability explained a significant portion of the variation in these element concentrations. Multiparous mares had lower concentrations of K and Zn, but higher concentrations of Na. Substantial serial and spatial correlation were detected in spite of modeling efforts to avoid the problem. As a result, p-values obtained for parameter estimates were likely biased toward zero. Nonetheless, results of this analysis indicate that monitoring changes in mammary-secretion element concentrations might reasonably be used as a predictor of impending parturition in the mare. In addition, these results suggest that element concentrations warrant attention in the development of neonatal milk-replacement therapies. This study demonstrates that non-linear regression can be used successfully to model time-series data in animal-health management. This approach should be considered by investigators facing similar analytical challenges. PMID:10782599
NASA Astrophysics Data System (ADS)
Haxton, Wick; Lunardini, Cecilia
2008-09-01
Semi-leptonic electroweak interactions in nuclei—such as β decay, μ capture, charged- and neutral-current neutrino reactions, and electron scattering—are described by a set of multipole operators carrying definite parity and angular momentum, obtained by projection from the underlying nuclear charge and three-current operators. If these nuclear operators are approximated by their one-body forms and expanded in the nucleon velocity through order |p→|/M, where p→ and M are the nucleon momentum and mass, a set of seven multipole operators is obtained. Nuclear structure calculations are often performed in a basis of Slater determinants formed from harmonic oscillator orbitals, a choice that allows translational invariance to be preserved. Harmonic-oscillator single-particle matrix elements of the multipole operators can be evaluated analytically and expressed in terms of finite polynomials in q, where q is the magnitude of the three-momentum transfer. While results for such matrix elements are available in tabular form, with certain restriction on quantum numbers, the task of determining the analytic form of a response function can still be quite tedious, requiring the folding of the tabulated matrix elements with the nuclear density matrix, and subsequent algebra to evaluate products of operators. Here we provide a Mathematica script for generating these matrix elements, which will allow users to carry out all such calculations by symbolic manipulation. This will eliminate the errors that may accompany hand calculations and speed the calculation of electroweak nuclear cross sections and rates. We illustrate the use of the new script by calculating the cross sections for charged- and neutral-current neutrino scattering in 12C. Program summaryProgram title: SevenOperators Catalogue identifier: AEAY_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAY_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Gazmeh, Meisam; Bahreini, Maryam; Tavassoli, Seyed Hassan; Asnaashari, Mohammad
2015-01-01
Introduction: In this study, laser induced breakdown spectroscopy (LIBS) is used for qualitative analysis of healthy and carious teeth. The technique of laser ablation is receiving increasing attention for applications in dentistry, specifically for the treatment of teeth such as drilling of micro-holes and plaque removal. Methods: A quality-switched (Q-switched) Neodymium-Doped Yttrium Aluminium Garnet (Nd:YAG) laser operating at wavelength of 1064 nm, pulse energy of 90 mJ/pulse, repetition rate of 2Hz and pulse duration of 6 ns was used in this analysis. In the process of ablation a luminous micro-plasma is normally generated which may be exploited for on-line elemental analysis via laser induced breakdown spectroscopy technique. We propose laser induced breakdown spectroscopy as a rapid, in situ and easy method for monitoring drilling process. Results: The results of elemental analysis show the presence of some trace elements in teeth including P, Ca, Mg, Zn, K, Sr, C, Na, H, O and the permeability of some amalgam (teeth filling materials) elements including Hg, Ag, Cu and Sn into dental matrix. Conclusion: This study addresses the ability of LIBS in elemental analysis of teeth and its feasibility in acute identification of healthy and carious teeth during drilling process for future clinical applications. PMID:25987971
NASA Astrophysics Data System (ADS)
Wu, Yueqian; Yang, Minglin; Sheng, Xinqing; Ren, Kuan Fang
2015-05-01
Light scattering properties of absorbing particles, such as the mineral dusts, attract a wide attention due to its importance in geophysical and environment researches. Due to the absorbing effect, light scattering properties of particles with absorption differ from those without absorption. Simple shaped absorbing particles such as spheres and spheroids have been well studied with different methods but little work on large complex shaped particles has been reported. In this paper, the surface Integral Equation (SIE) with Multilevel Fast Multipole Algorithm (MLFMA) is applied to study scattering properties of large non-spherical absorbing particles. SIEs are carefully discretized with piecewise linear basis functions on triangle patches to model whole surface of the particle, hence computation resource needs increase much more slowly with the particle size parameter than the volume discretized methods. To improve further its capability, MLFMA is well parallelized with Message Passing Interface (MPI) on distributed memory computer platform. Without loss of generality, we choose the computation of scattering matrix elements of absorbing dust particles as an example. The comparison of the scattering matrix elements computed by our method and the discrete dipole approximation method (DDA) for an ellipsoid dust particle shows that the precision of our method is very good. The scattering matrix elements of large ellipsoid dusts with different aspect ratios and size parameters are computed. To show the capability of the presented algorithm for complex shaped particles, scattering by asymmetry Chebyshev particle with size parameter larger than 600 of complex refractive index m = 1.555 + 0.004 i and different orientations are studied.
The nuclear matrix elements of 0νββ decay and the NUMEN project at INFN-LNS
NASA Astrophysics Data System (ADS)
Cappuzzello, F.; Agodi, C.; Balestra, F.; Bijker, R.; Bonanno, D.; Bongiovanni, D.; Branchina, V.; Calabrese, S.; Calabretta, L.; Calanna, A.; Calvo, D.; Carbone, D.; Cavallaro, M.; Colonna, M.; Ferrero, S.; Foti, A.; Finocchiaro, P.; Giraudo, G.; Greco, V.; Iazzi, F.; Introzzi, R.; Lanzalone, G.; Lavagno, A.; Lo Presti, D.; Longhitano, F.; Muoio, A.; Pandola, L.; Rifuggiato, D.; Ruslan, M. V.; Santopinto, E.; Scaltrito, L.; Tudisco, S.; Zagatto, V.
2016-05-01
An innovative technique to access the nuclear matrix elements entering the expression of the life time of the double beta decay by relevant cross sections measurements of double charge exchange reactions is proposed. A key aspect of the project is the use of the MAGNEX large acceptance magnetic spectrometer, for the detection of the ejectiles, and of the LNS K800 Superconducting Cyclotron (CS), for the acceleration of the required high resolution and low emittance heavyion beams, already in operation at INFN Laboratory Nazionali del Sud in Catania (Italy).
Agodi, C. Calabretta, L.; Calanna, A.; Carbone, D.; Cavallaro, M.; Colonna, M.; Cuttone, G.; Finocchiaro, P.; Pandola, L.; Rifuggiato, D.; Tudisco, S.; Cappuzzello, F.; Greco, V.; Bonanno, D. L.; Bongiovanni, D. G.; Longhitano, F.; Branchina, V.; Foti, A.; Lo Presti, D.; Lanzalone, G.; and others
2015-10-28
In the NUMEN Project it is proposed an innovative technique to access the nuclear matrix elements entering in the expression of the life-time of the neutrinoless double beta decay, using relevant cross sections of double charge exchange reactions. A key aspect is the use of MAGNEX large acceptance magnetic spectrometer, for the detection of the ejectiles, and of the INFN Laboratori Nazionali del Sud (LNS) K800 Superconducting Cyclotron (CS), for the acceleration of the required high resolution and low emittance heavy-ion beams.
NASA Astrophysics Data System (ADS)
Agodi, C.; Cappuzzello, F.; Bonanno, D. L.; Bongiovanni, D. G.; Branchina, V.; Calabretta, L.; Calanna, A.; Carbone, D.; Cavallaro, M.; Colonna, M.; Cuttone, G.; Foti, A.; Finocchiaro, P.; Greco, V.; Lanzalone, G.; Lo Presti, D.; Longhitano, F.; Muoio, A.; Pandola, L.; Rifuggiato, D.; Tudisco, S.
2015-10-01
In the NUMEN Project it is proposed an innovative technique to access the nuclear matrix elements entering in the expression of the life-time of the neutrinoless double beta decay, using relevant cross sections of double charge exchange reactions. A key aspect is the use of MAGNEX large acceptance magnetic spectrometer, for the detection of the ejectiles, and of the INFN Laboratori Nazionali del Sud (LNS) K800 Superconducting Cyclotron (CS), for the acceleration of the required high resolution and low emittance heavy-ion beams.
Investigation of the E2 and E3 matrix elements in 200Hg using direct nuclear reactions
NASA Astrophysics Data System (ADS)
Rand, Evan; Bildstein, Vinzenz; Diaz Varela, Alejandra; Garrett, Paul; Hadinia, Baharak; Jamieson, Drew; Jigmeddorj, Badamsambuu; Laffoley, Alex; Leach, Kyle; Maclean, Andrew; Svensson, Carl; Ball, Gordon; Faestermann, Thomas; Hertenberger, Ralf; Wirth, Hans-Friedrich
2014-09-01
To date, 199Hg provides the most stringent limit on an atomic electric dipole moment (EDM). The existence of a permanent EDM would be a clear signal of CP violation from new physics beyond the Standard Model. Theoretical nuclear-structure calculations for 199Hg are challenging, and give varied predictions for the excited-state spectrum. Understanding the E2 and E3 strengths in 198 , 199 , 200Hg will make it possible to develop a nuclear structure model for the Schiff strength based on these matrix elements, and thereby constrain present models that predict the contribution of octupole collectivity to the Schiff moment of the nucleus. This work comprises two experiments using the Q3D magnetic spectrograph at the Maier-Leibnitz Laboratory. These experiments utilized a 22 MeV deuteron beam incident on a target of 200Hg32S. The first experiment accesses the E 2 and E 3 matrix elements in 200Hg via inelastic deuteron scattering. The second experiment, 200Hg(d , t) 199Hg, yields important information on the single-particle nature of 199Hg. Preliminary results will be presented. To date, 199Hg provides the most stringent limit on an atomic electric dipole moment (EDM). The existence of a permanent EDM would be a clear signal of CP violation from new physics beyond the Standard Model. Theoretical nuclear-structure calculations for 199Hg are challenging, and give varied predictions for the excited-state spectrum. Understanding the E2 and E3 strengths in 198 , 199 , 200Hg will make it possible to develop a nuclear structure model for the Schiff strength based on these matrix elements, and thereby constrain present models that predict the contribution of octupole collectivity to the Schiff moment of the nucleus. This work comprises two experiments using the Q3D magnetic spectrograph at the Maier-Leibnitz Laboratory. These experiments utilized a 22 MeV deuteron beam incident on a target of 200Hg32S. The first experiment accesses the E 2 and E 3 matrix elements in 200Hg via
Wang, X Y; Zhang, J H; Sun, Q L; Yao, Z Y; Deng, B G; Guo, W Y; Wang, L; Dong, W H; Wang, F; Zhao, C P; Wang, T Y
2015-01-01
Preliminary studies have suggested that a characteristic element of the matrix attachment region (MAR) in human interferon-β mediates the adhesion of vectors to Chinese hamster ovary (CHO) cells. In this study, we investigated if vector adhesion increased nerve growth factor (NGF) expression in CHO cells. The MAR characteristic element sequence of human interferon-β was inserted into the multiple-cloning site of the pEGFP-C1 vector. The target NGF gene was inserted upstream of the MAR characteristic element sequence to construct the MAR/NGF expression vector. The recombinant plasmid was transfected into CHO cells and stable monoclonal cells were selected using G418. NGF mRNA and protein expression was detected by reverse transcriptase-polymerase chain reaction and enzyme-linked immunosorbent assay, respectively. Plasmid reduction experiments were used to determine the state of transfected plasmid in mammalian cells. The insertion of MAR into the vector increased NGF expression levels in CHO cells (1.93- fold) compared to the control. The recombinant plasmid expressing the MAR sequence was digested into a linear space vector. The inserted MAR and NGF sequences were consistent with those inserted into the plasmid before recombination. Therefore, we concluded that the MAR characteristic element mediates vector adhesion to CHO cells and enhances the stability and efficiency of the target gene expression. PMID:26345852
Kolb, Florian; Schmoltner, Kerstin; Huth, Michael; Hohenau, Andreas; Krenn, Joachim; Klug, Andreas; List, Emil J W; Plank, Harald
2013-08-01
The development of simple gas sensing concepts is still of great interest for science and technology. The demands on an ideal device would be a single-step fabrication method providing a device which is sensitive, analyte-selective, quantitative, and reversible without special operating/reformation conditions such as high temperatures or special environments. In this study we demonstrate a new gas sensing concept based on a nanosized PtC metal-matrix system fabricated in a single step via focused electron beam induced deposition (FEBID). The sensors react selectively on polar H2O molecules quantitatively and reversibly without any special reformation conditions after detection events, whereas non-polar species (O2, CO2, N2) produce no response. The key elements are isolated Pt nanograins (2-3 nm) which are embedded in a dielectric carbon matrix. The electrical transport in such materials is based on tunneling effects in the correlated variable range hopping regime, where the dielectric carbon matrix screens the electric field between the particles, which governs the final conductivity. The specific change of these dielectric properties by the physisorption of polar gas molecules (H2O) can change the tunneling probability and thus the overall conductivity, allowing their application as a simple and straightforward sensing concept. PMID:23818049
Study on Hankel matrix-based SVD and its application in rolling element bearing fault diagnosis
NASA Astrophysics Data System (ADS)
Jiang, Huiming; Chen, Jin; Dong, Guangming; Liu, Tao; Chen, Gang
2015-02-01
Based on the traditional theory of singular value decomposition (SVD), singular values (SVs) and ratios of neighboring singular values (NSVRs) are introduced to the feature extraction of vibration signals. The proposed feature extraction method is called SV-NSVR. Combined with selected SV-NSVR features, continuous hidden Markov model (CHMM) is used to realize the automatic classification. Then the SV-NSVR and CHMM based method is applied in fault diagnosis and performance assessment of rolling element bearings. The simulation and experimental results show that this method has a higher accuracy for the bearing fault diagnosis compared with those using other SVD features, and it is effective for the performance assessment of rolling element bearings.
Cwik, T.; Katz, D.S.
1996-12-31
Finite element modeling has proven useful for accurately simulating scattered or radiated electromagnetic fields from complex three-dimensional objects whose geometry varies on the scale of a fraction of an electrical wavelength. An unstructured finite element model of realistic objects leads to a large, sparse, system of equations that needs to be solved efficiently with regard to machine memory and execution time. Both factorization and iterative solvers can be used to produce solutions to these systems of equations. Factorization leads to high memory requirements that limit the electrical problem size of three-dimensional objects that can be modeled. An iterative solver can be used to efficiently solve the system without excessive memory use and in a minimal amount of time if the convergence rate is controlled.
Direct finite element solution on an optical laboratory matrix-vector processor
NASA Technical Reports Server (NTRS)
Casasent, David; Riedl, Steven
1988-01-01
The first optical laboratory system results employing a direct LU decomposition solution of a system of linear algebraic equations are presented for a finite element problem solution. This also represents the first laboratory demonstration of the use of sign-magnitude negative number representation as well as new bit partitioning techniques to increase the accuracy of an optical encoded processor beyond the number of bit channels available.
NASA Technical Reports Server (NTRS)
Duffy, Stephen F.
1998-01-01
The development of modeling approaches for the failure analysis of ceramic-based material systems used in high temperature environments was the primary objective of this research effort. These materials have the potential to support many key engineering technologies related to the design of aeropropulsion systems. Monolithic ceramics exhibit a number of useful properties such as retention of strength at high temperatures, chemical inertness, and low density. However, the use of monolithic ceramics has been limited by their inherent brittleness and a large variation in strength. This behavior has motivated material scientists to reinforce the monolithic material with a ceramic fiber. The addition of a second ceramic phase with an optimized interface increases toughness and marginally increases strength. The primary purpose of the fiber is to arrest crack growth, not to increase strength. The material systems of interest in this research effort were laminated ceramic matrix composites, as well as two- and three- dimensional fabric reinforced ceramic composites. These emerging composite systems can compete with metals in many demanding applications. However, the ongoing metamorphosis of ceramic composite material systems, and the lack of standardized design data has in the past tended to minimize research efforts related to structural analysis. Many structural components fabricated from ceramic matrix composites (CMC) have been designed by "trial and error." The justification for this approach lies in the fact that during the initial developmental phases for a material system fabrication issues are paramount. Emphasis is placed on demonstrating feasibility rather than fully understanding the processes controlling mechanical behavior. This is understandable during periods of rapid improvements in material properties for any composite system. But to avoid the ad hoc approach, the analytical methods developed under this effort can be used to develop rational structural
NASA Technical Reports Server (NTRS)
Haque, A.; Ahmed, L.; Ware, T.; Jeelani, S.; Verrilli, Michael J. (Technical Monitor)
2001-01-01
The stress concentrations associated with circular notches and subjected to uniform tensile loading in woven ceramic matrix composites (CMCs) have been investigated for high-efficient turbine engine applications. The CMC's were composed of Nicalon silicon carbide woven fabric in SiNC matrix manufactured through polymer impregnation process (PIP). Several combinations of hole diameter/plate width ratios and ply orientations were considered in this study. In the first part, the stress concentrations were calculated measuring strain distributions surrounding the hole using strain gages at different locations of the specimens during the initial portion of the stress-strain curve before any microdamage developed. The stress concentration was also calculated analytically using Lekhnitskii's solution for orthotropic plates. A finite-width correction factor for anisotropic and orthotropic composite plate was considered. The stress distributions surrounding the circular hole of a CMC's plate were further studied using finite element analysis. Both solid and shell elements were considered. The experimental results were compared with both the analytical and finite element solutions. Extensive optical and scanning electron microscopic examinations were carried out for identifying the fracture behavior and failure mechanisms of both the notched and notched specimens. The stress concentration factors (SCF) determined by analytical method overpredicted the experimental results. But the numerical solution underpredicted the experimental SCF. Stress concentration factors are shown to increase with enlarged hole size and the effects of ply orientations on stress concentration factors are observed to be negligible. In all the cases, the crack initiated at the notch edge and propagated along the width towards the edge of the specimens.
Gardiner, Bruce S.; Wong, Kelvin K. L.; Joldes, Grand R.; Rich, Addison J.; Tan, Chin Wee; Burgess, Antony W.; Smith, David W.
2015-01-01
This paper presents a framework for modelling biological tissues based on discrete particles. Cell components (e.g. cell membranes, cell cytoskeleton, cell nucleus) and extracellular matrix (e.g. collagen) are represented using collections of particles. Simple particle to particle interaction laws are used to simulate and control complex physical interaction types (e.g. cell-cell adhesion via cadherins, integrin basement membrane attachment, cytoskeletal mechanical properties). Particles may be given the capacity to change their properties and behaviours in response to changes in the cellular microenvironment (e.g., in response to cell-cell signalling or mechanical loadings). Each particle is in effect an ‘agent’, meaning that the agent can sense local environmental information and respond according to pre-determined or stochastic events. The behaviour of the proposed framework is exemplified through several biological problems of ongoing interest. These examples illustrate how the modelling framework allows enormous flexibility for representing the mechanical behaviour of different tissues, and we argue this is a more intuitive approach than perhaps offered by traditional continuum methods. Because of this flexibility, we believe the discrete modelling framework provides an avenue for biologists and bioengineers to explore the behaviour of tissue systems in a computational laboratory. PMID:26452000
NASA Astrophysics Data System (ADS)
Maury, Cécile; Sirven, Jean-Baptiste; Tabarant, Michel; L'Hermite, Daniel; Courouau, Jean-Louis; Gallou, Catherine; Caron, Nadège; Moutiers, Gilles; Cabuil, Valérie
2013-04-01
Liquid sodium is used as coolant in sodium-cooled fast nuclear reactors. Among many parameters to monitor to ensure the safe operation of the reactor, the coolant chemical purity is a relevant indicator of several undesirable situations, like corrosion of structural materials or sodium contamination, which may release different elements in the coolant. Several techniques have already been implemented to measure the sodium purity, but their response time is long and not suited for continuous monitoring. Therefore, as a complement to them, laser-induced breakdown spectroscopy (LIBS) is considered as a promising technique for real-time analysis of the coolant purity. In this paper we report on the first LIBS quantitative measurements performed in liquid sodium at 150 °C. Calibration curves were traced for lead and indium using the standard addition method. Important intensity drifts and fluctuations were observed, mostly due to pressure variations in the sodium oven. Background subtraction and/or normalization was used to compensate for those intensity fluctuations. To describe the effect of these corrections on the analytical signal noise, a simple model was proposed and its results were found to satisfactorily fit the experimental data. Using this approach, the best detection limits were obtained for the background-subtracted and normalized data, and were found to be 6 ppm for lead and 5 ppm for indium.
Williams, M.; Applegate, D.; Bellis, M.; Meyer, C. A.; Dey, B.; Dickson, R.; Krahn, Z.; McCracken, M. E.; Moriya, K.; Schumacher, R. A.; Adhikari, K. P.; Careccia, S. L.; Dodge, G. E.; Klein, A.; Mayer, M.; Nepali, C. S.; Niroula, M. R.; Seraydaryan, H.; Tkachenko, S.; Weinstein, L. B.
2009-12-15
High-statistics differential cross sections and spin-density matrix elements for the reaction {gamma}p{yields}p{omega} have been measured using the CEBAF large acceptance spectrometer (CLAS) at Jefferson Lab for center-of-mass (c.m.) energies from threshold up to 2.84 GeV. Results are reported in 112 10-MeV wide c.m. energy bins, each subdivided into cos{theta}{sub c.m.}{sup {omega}} bins of width 0.1. These are the most precise and extensive {omega} photoproduction measurements to date. A number of prominent structures are clearly present in the data. Many of these have not previously been observed due to limited statistics in earlier measurements.
The (d,2He) reaction on Mo96 and the double-β decay matrix elements for Zr96
NASA Astrophysics Data System (ADS)
Dohmann, H.; Bäumer, C.; Frekers, D.; Grewe, E.-W.; Harakeh, M. N.; Hollstein, S.; Johansson, H.; Popescu, L.; Rakers, S.; Savran, D.; Simon, H.; Thies, J. H.; van den Berg, A. M.; Wörtche, H. J.; Zilges, A.
2008-10-01
The 96Mo(d,2He)96Nb charge-exchange reaction was investigated at an incident energy of Ed=183.5 MeV. An excitation-energy resolution of 110 keV was achieved. The experiment was performed at KVI, Groningen, using the magnetic spectrometer BBS at three angular positions: 0°,2.5°, and 6°. We found that below 6 MeV almost the entire Gamow-Teller (GT+) strength is concentrated in a single state at 0.69 MeV excitation energy. As Mo96 is the daughter of the ββ decay nucleus Zr96, the present result provides information about the nuclear matrix elements active in the 2νββ decay of Zr96.
Menendez, J.; Poves, A.
2009-10-15
We discuss the variation of the nuclear matrix element (NME) for the neutrinoless double beta (0{nu}{beta}{beta}) decay of {sup 76}Ge when the wave functions are constrained to reproduce the experimental occupancies of the two nuclei involved in the transition. In the interacting shell model description the value of the NME is enhanced about 15% compared to previous calculations, whereas in the QRPA the NME's are reduced by 20%-30%. This diminishes the discrepancies between both approaches. In addition, we discuss the effect of the short-range correlations on the NME in light of the recently proposed parametrizations based on a consistent renormalization of the 0{nu}{beta}{beta} transition operator.
NASA Astrophysics Data System (ADS)
Turovtsev, V. V.; Orlov, Yu. D.; Tsirulev, A. N.
2015-08-01
The advantages of the orthonormal basis set of 2π-periodic Mathieu functions compared to the trigonometric basis set in calculations of torsional states of molecules are substantiated. Explicit expressions are derived for calculating the Hamiltonian matrix elements of a one-dimensional torsional Schrödinger equation with a periodic potential of the general form in the basis set of Mathieu functions. It is shown that variation of a parameter of Mathieu functions allows the rotation potential and the structural function to be approximated with a good accuracy by a small number of series terms. The conditions for the best choice of this parameter are specified, and approximations are obtained for torsional potentials of n-butane upon rotation about the central C-C bond and of its univalent radical n-butyl C2H5C·H2 upon rotation of the C·H2 group. All algorithms are implemented in the Maple package.
The (d,2He) reaction on Se76 and the double-β-decay matrix elements for A=76
NASA Astrophysics Data System (ADS)
Grewe, E.-W.; Bäumer, C.; Dohmann, H.; Frekers, D.; Harakeh, M. N.; Hollstein, S.; Johansson, H.; Popescu, L.; Rakers, S.; Savran, D.; Simon, H.; Thies, J. H.; van den Berg, A. M.; Wörtche, H. J.; Zilges, A.
2008-10-01
The (d,2He) charge-exchange reaction on Se76 was studied at an incident energy of 183 MeV. The outgoing two protons in the 1S0 state, referred to as He2, were both momentum analyzed and detected by the same spectrometer and detector. The experiment was performed at KVI, Groningen, using the magnetic spectrometer BBS at three angular positions: 0°,2.5°, and 5°. Excitation-energy spectra of the residual nucleus As76 were obtained with an energy resolution of about 120 keV (FWHM). Gamow-Teller (GT+) transition strengths were extracted up to 5 MeV and compared with those from an (n,p) experiment at low resolution. Together with the GT- transition strengths from the 76Ge(p,n) experiment leading to the same intermediate nucleus, the nuclear matrix element of the two-neutrino double-β decay of Ge76 was evaluated.
NASA Astrophysics Data System (ADS)
Bierenbaum, Isabella; Blümlein, Johannes; Klein, Sebastian
2007-09-01
We calculate the O(αs2) massive operator matrix elements for the twist-2 operators, which contribute to the heavy flavor Wilson coefficients in unpolarized deeply inelastic scattering in the region Q≫m. The calculation has been performed using light-cone expansion techniques. We confirm an earlier result obtained in [M. Buza, Y. Matiounine, J. Smith, R. Migneron, W.L. van Neerven, Nucl. Phys. B 472 (1996) 611, arxiv:/hep-ph/9601302]. The calculation is carried out without using the integration-by-parts method and in Mellin space using harmonic sums, which lead to a significant compactification of the analytic results derived previously. The results allow to determine the heavy flavor Wilson coefficients for F(x,Q) to O(αs2) and for F(x,Q) to O(αs3) for all but the power suppressed terms ∝(/Q)k,k⩾1.
The nuclear matrix elements of 0vββ decay and the NUMEN project at INFN-LNS
NASA Astrophysics Data System (ADS)
Cappuzzello, F.; Agodi, C.; Aciksoz, E.; Acosta, L.; Aslanouglou, X.; Auerbach, N.; Bijker, R.; Bonanno, D.; Bongiovanni, D.; Borello, T.; Boudhaim, S.; Bouhssa, M. L.; Boztosun, I.; Calabretta, L.; Calanna, A.; Carbone, D.; Cavallaro, M.; Calvo, D.; Chávez Lomelí, E. R.; Colonna, M.; D'Agostino, G.; Deshmukh, N.; de Faria, P. N.; Ferrero, A.; Foti, A.; Finocchiaro, P.; Gomes, P. R. S.; Greco, V.; Hacisalihoglu, A.; Housni, Z.; Khouaja, A.; Inchaou, J.; Lanzalone, G.; La Via, F.; Lay, J. A.; Lenske, H.; Linares, R.; Lubian, J.; Iazzi, F.; Introzzi, R.; Lavagno, A.; Lo Presti, D.; Medina, N.; Mendes, D. R.; Muoio, A.; Oliveira, J. R. B.; Pakou, A.; Pandola, L.; Rifuggiato, D.; Rodrigues, M. R. D.; Santagati, G.; Santopinto, E.; Scaltrito, L.; Sgouros, O.; Solakcı, S. O.; Soukeras, V.; Tudisco, S.; Vsevolodovna, R. I. M.; Zagatto, V.
2016-07-01
An innovative technique to access the nuclear matrix elements entering the expression of the life time of the double beta decay by relevant cross section measurements of double charge exchange reactions is proposed. A key aspect of the project is the use of the MAGNEX large acceptance magnetic spectrometer, for the detection of the ejectiles, and of the LNS K800 Superconducting Cyclotron (CS), for the acceleration of the required high resolution and low emittance heavy-ion beams, already in operation at INFN Laboratory Nazionali del Sud in Catania (Italy). However, a major upgrade is foreseen for the INFN-LNS research infrastructure to cope with beam currents as high as several ppA required by the project.
Am phases in the matrix of a U–Pu–Zr alloy with Np, Am, and rare-earth elements
Janney, Dawn E.; Kennedy, J. Rory; Madden, James W.; O’Holleran, Thomas P.
2015-01-01
Phases and microstructures in the matrix of an as-cast U-Pu-Zr alloy with 3 wt% Am, 2% Np, and 8% rare-earth elements were characterized by scanning and transmission electron microscopy. The matrix consists primarily of two phases, both of which contain Am: ζ-(U, Np, Pu, Am) (~70 at% U, 5% Np, 14% Pu, 1% Am, and 10% Zr) and δ-(U, Np, Pu, Am)Zr_{2} (~25% U, 2% Np, 10-15% Pu, 1-2% Am, and 55-60 at% Zr). These phases are similar to those in U-Pu-Zr alloys, although the Zr content in ζ-(U, Np, Pu, Am) is higher than that in ζ-(U, Pu) and the Zr content in δ-(U, Np, Pu, Am)Zr_{2} is lower than that in δ-UZr_{2}. Nanocrystalline actinide oxides with structures similar to UO2 occurred in some areas, but may have formed by reactions with the atmosphere during sample handling. Planar features consisting of a central zone of ζ-(U, Np, Pu, Am) bracketed by zones of δ-(U, Np, Pu, Am)Zr_{2} bound irregular polygons ranging in size from a few micrometers to a few tens of micrometers across. The rest of the matrix consists of elongated domains of ζ-(U, Np, Pu, Am) and δ-(U, Np, Pu, Am)Zr_{2}. Each of these domains is a few tens of nanometers across and a few hundred nanometers long. The domains display strong preferred orientations involving areas a few hundred nanometers to a few micrometers across.
Spin density matrix elements for radiative decays of the omega meson in photoproduction at 5 GeV
NASA Astrophysics Data System (ADS)
Mokaya, Fridah
2016-03-01
The photoproduction of ω(782) meson on the nucleon at high energies is well described by a sum of t-channel exchanges. In the high energy limit of diffractive scattering, where Pomeron exchange dominates the total cross section, the helicity of the incident photon is transferred directly to the vector meson. At intermediate energies, other Regge exchanges compete with the Pomeron, leading to a complex energy dependence in the spin density matrix for vector mesons like the omega. High statistics measurements of the spin density matrix elements for the reaction γp ωp, ω π0 γ are presented based on data taken with the Radphi experiment at Jefferson Lab in the energy range 4.4 - 5.5 GeV. The results binned in Eγ and |t | are analysed in both the Gottfried Jackson and s-channel helicity frames and compared to a model with the Pomeron and other Regge exchanges contributing to the omega meson photoproduction amplitude.
NASA Astrophysics Data System (ADS)
Stafilov, Trajče; Zendelovska, Dragica; Pavlovska, Gorica; Čundeva, Katarina
2002-05-01
The interferences of Ca and Mg as matrix elements in dolomite and gypsum on Ag, Cd, Cr, Mn, Tl and Zn absorbances during their electrothermal atomic absorption spectrometric (ETAAS) determination are investigated. The results reveal that Ca and Mg do not interfere on Zn and Mn, tend to decrease absorbances of Ag, Cd and Cr, while Tl suffers the most significant influence. A flotation separation method is proposed to eliminate matrix interferences. Hydrated iron(III) oxide, Fe 2O 3· xH 2O, and iron(III) hexamethylenedithiocarbamate, Fe(HMDTC) 3, are applied as flotation collectors. The influence of hydrophobic dithiocarbamate anion, HMDTC, on flotation recoveries of each analyte is studied. The most suitable concentrations of dolomite and gypsum solutions for flotation are determined. To avoid flotation suppression due to the reaction of Ca 2+ and Mg 2+ with surfactant ions, a fit foaming agent was selected. The elements present in dolomite and gypsum as traces have been analyzed by ETAAS. Their ETAAS limits of detection following flotation are found to be 0.021 μg·g -1 for Ag, 0.019 μg·g -1 for Cd, 0.014 μg·g -1 for Cr and 0.11 μg·g -1 for Tl. The determination of Mn and Zn can be performed by flame AAS (FAAS). The limit of detection for Mn is 1.5 μg·g -1, while for Zn 0.8 μg·g -1.
NASA Astrophysics Data System (ADS)
Upcraft, L. M.
2010-09-01
There is an ongoing need for numerically efficient algorithms that are capable of calculating the radiative and collisional rates of arbitrarily complex ions that are present in hot plasmas to a level of accuracy that surpasses that available in many existing approximations. Hydrogen-like solutions for determining these rates in more general ions by use of an effective (and generally non-integer) atomic number frequently give poor results and are of limited validity. This paper illustrates that results accurate to of order 20% can be obtained for matrix elements of both rates for arbitrarily complex ions by use of hydrogenic wavefunctions that use different effective atomic numbers for the initial and final sub-shells. Not only does this allow for the realistic modelling of inner shell transitions, it naturally allows for the physical effect of orbital relaxation. It is shown that the integral of the generalised oscillator strength used by the Plane-wave Born approximation has an analytic solution that can be reduced to a form suitable for efficient numerical integration over an arbitrary electron distribution. Extensive use of the computer algebra package Mathematica ® has generated a unique formula for each transition and the results have been transformed to efficient fortran 90 code for all transitions between non-relativistic sub-shells with principal quantum numbers n ≤ 10. In the case of the collisional matrix elements these are typically two to three orders of magnitude faster to calculate than by direct numerical evaluation. The fortran code is available upon request from the author.
On-shell {delta}I=3/2 kaon weak matrix elements with nonzero total momentum
Yamazaki, Takeshi
2009-05-01
We present our results for the on-shell {delta}I=3/2 kaon decay matrix elements using domain wall fermions and the DBW2 gauge action at one coarse lattice spacing corresponding to a{sup -1}=1.31 GeV in the quenched approximation. The on-shell matrix elements are evaluated in two different frames: the center-of-mass frame and nonzero total-momentum frame. We employ the formula proposed by Lellouch and Luescher in the center-of-mass frame, and its extension for a nonzero total-momentum frame to extract the infinite volume, on-shell, center-of-mass frame decay amplitudes. We determine the decay amplitude at the physical pion mass and momentum from the chiral extrapolation and an interpolation of the relative momentum using the results calculated in the two frames. We have obtained ReA{sub 2}=1.66(23)((+48/-03))((+53/-0))x10{sup -8} GeV and ImA{sub 2}=-1.181(26)((+141/-014))((+44/-0))x10{sup -12} GeV at the physical point, using the data at the relatively large pion mass, m{sub {pi}}>0.35 GeV. The first error is statistic, and the second and third are systematic. The second error is estimated with several fits of the chiral extrapolation including the (quenched) chiral perturbation formula at next to leading order using only lighter pion masses. The third one is estimated with an analysis using the lattice dispersion relation. The result of ReA{sub 2} is reasonably consistent with experiment.
FERM3D: A finite element R-matrix electron molecule scattering code
NASA Astrophysics Data System (ADS)
Tonzani, Stefano
2007-01-01
FERM3D is a three-dimensional finite element program, for the elastic scattering of a low energy electron from a general polyatomic molecule, which is converted to a potential scattering problem. The code is based on tricubic polynomials in spherical coordinates. The electron-molecule interaction is treated as a sum of three terms: electrostatic, exchange, and polarization. The electrostatic term can be extracted directly from ab initio codes ( GAUSSIAN 98 in the work described here), while the exchange term is approximated using a local density functional. A local polarization potential based on density functional theory [C. Lee, W. Yang, R.G. Parr, Phys. Rev. B 37 (1988) 785] describes the long range attraction to the molecular target induced by the scattering electron. Photoionization calculations are also possible and illustrated in the present work. The generality and simplicity of the approach is important in extending electron-scattering calculations to more complex targets than it is possible with other methods. Program summaryTitle of program:FERM3D Catalogue identifier:ADYL_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYL_v1_0 Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computer for which the program is designed and others on which it has been tested:Intel Xeon, AMD Opteron 64 bit, Compaq Alpha Operating systems or monitors under which the program has been tested:HP Tru64 Unix v5.1, Red Hat Linux Enterprise 3 Programming language used:Fortran 90 Memory required to execute with typical data:900 MB (neutral CO 2), 2.3 GB (ionic CO 2), 1.4 GB (benzene) No. of bits in a word:32 No. of processors used:1 Has the code been vectorized?:No No. of lines in distributed program, including test data, etc.:58 383 No. of bytes in distributed program, including test data, etc.:561 653 Distribution format:tar.gzip file CPC Program library subprograms used:ADDA, ACDP Nature of physical problem:Scattering of an
NASA Astrophysics Data System (ADS)
Regueiro, R. A.; Yu, S.
2010-12-01
The paper models grain-scale micro-cracking in shale at grain-matrix interfaces, assuming constituents are composed of quart silt grains and compacted clay matrix for a typical shale. The influence of grain-matrix-grain interaction on micro-crack patterns is investigated. Elasto-plastic pressure-sensitive cohesive-surface models are inserted at grain-matrix interfaces and intra-clay-matrix finite element facets, while a bulk elasto-plasticity model with bifurcation is employed for the clay matrix to compare to the intra-clay-matrix cohesive-surface model. Numerical examples are presented under two-dimensional plane strain condition at small strains. A procedure is proposed to upscale grain-scale micro-cracking to predict macro-fracture nucleation and propagation in shale and other bound particulate materials. It is shown that using cohesive surface elements (CSEs) at all finite element facets in the clay matrix mesh to simulate micro-cracking in the clay matrix leads to mesh-dependent results. Using CSEs at grain-clay-matrix interfaces is physical and not mesh dependent. We also considered using bulk pressure-sensitive elasto-plasticity with bifurcation condition within the clay matrix to attempt to predict onset of localization around grains in the simulations. It was encouraging to see that for both the single grain and multiple grain simulations, the finite element region in the clay matrix meshes where bifurcation was first detected around the grains was nearly the same. This gives us confidence that once a proper post-bifurcation constitutive model is implemented within an embedded discontinuity formulation, micro-cracking nucleation and propagation at the grain-scale in shale can be properly simulated, which will provide the basis for up-scaling to macro-cracks within a multiscale method for fracture in shale. Other items to address in future research are: (i) include transverse isotropy (elastic and plastic) for the bulk clay matrix elasto-plasticity model
Saucedo-Mora, L; Marrow, T J
2016-07-13
The problem of multi-scale modelling of damage development in a SiC ceramic fibre-reinforced SiC matrix ceramic composite tube is addressed, with the objective of demonstrating the ability of the finite-element microstructure meshfree (FEMME) model to introduce important aspects of the microstructure into a larger scale model of the component. These are particularly the location, orientation and geometry of significant porosity and the load-carrying capability and quasi-brittle failure behaviour of the fibre tows. The FEMME model uses finite-element and cellular automata layers, connected by a meshfree layer, to efficiently couple the damage in the microstructure with the strain field at the component level. Comparison is made with experimental observations of damage development in an axially loaded composite tube, studied by X-ray computed tomography and digital volume correlation. Recommendations are made for further development of the model to achieve greater fidelity to the microstructure. This article is part of the themed issue 'Multiscale modelling of the structural integrity of composite materials'. PMID:27242308
NASA Technical Reports Server (NTRS)
Abdul-Aziz, Ali; Baaklini, George Y.; Bhatt, Ramakrishna T.
2003-01-01
Most reverse engineering approaches involve imaging or digitizing an object and then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. The rapid prototyping technique builds high-quality physical prototypes directly from computer-aided design files. This fundamental technique for interpreting and interacting with large data sets is being used here via Velocity2 (an integrated image-processing software, ref. 1) using computed tomography (CT) data to produce a prototype three-dimensional test specimen model for analyses. A study at the NASA Glenn Research Center proposes to use these capabilities to conduct a combined nondestructive evaluation (NDE) and finite element analysis (FEA) to screen pretest and posttest structural anomalies in structural components. A tensile specimen made of silicon nitrite (Si3N4) ceramic matrix composite was considered to evaluate structural durability and deformity. Ceramic matrix composites are being sought as candidate materials to replace nickel-base superalloys for turbine engine applications. They have the unique characteristics of being able to withstand higher operating temperatures and harsh combustion environments. In addition, their low densities relative to metals help reduce component mass (ref. 2). Detailed three-dimensional volume rendering of the tensile test specimen was successfully carried out with Velocity2 (ref. 1) using two-dimensional images that were generated via computed tomography. Subsequent, three-dimensional finite element analyses were performed, and the results obtained were compared with those predicted by NDE-based calculations and experimental tests. It was shown that Velocity2 software can be used to render a three-dimensional object from a series of CT scan images with a minimum level of complexity. The analytical results (ref. 3) show that the high-stress regions correlated well with the damage sites identified by the CT scans
Judd, B.R.; Lo, E.
1996-01-01
Single-electron excitations of the atomic f shell are usually taken into account in analyses of lanthanide and actinide spectra by including six three-electron scalar operators t{sub i} in the Hamiltonian. Their matrix elements have been factorized for all f{sup N} by using the Wigner-Eckart theorem applied to Racah`s groups SO(7) and G{sub 2}. The two component parts, namely a reduced matrix element in G{sub 2} and an isoscalar factor, are tabulated. This provides a compact representation of the values of the matrix elements and also enables unusual selection rules and proportionalities to be exposed. 16 refs., 2 tabs.
Parrish, Robert M; Hohenstein, Edward G; Schunck, Nicolas F; Sherrill, C David; Martínez, Todd J
2013-09-27
Configuration-space matrix elements of N-body potentials arise naturally and ubiquitously in the Ritz-Galerkin solution of many-body quantum problems. For the common specialization of local, finite-range potentials, we develop the exact tensor hypercontraction method, which provides a quantized renormalization of the coordinate-space form of the N-body potential, allowing for a highly separable tensor factorization of the configuration-space matrix elements. This representation allows for substantial computational savings in chemical, atomic, and nuclear physics simulations, particularly with respect to difficult "exchangelike" contractions. PMID:24116775
NASA Astrophysics Data System (ADS)
Rath, P. K.; Chandra, R.; Chaturvedi, K.; Lohani, P.; Raina, P. K.
2016-02-01
The model-dependent uncertainties in the nuclear transition matrix elements for the Majoron-accompanied neutrinoless double-β decay (0+→0+transition) of Zr,9694, 100Mo, Te,130128, and 150Nd isotopes are calculated by employing the projected-Hartree-Fock-Bogoliubov formalism with four different parametrizations of the pairing plus multipolar two-body interactions and three different parametrizations of the Jastrow short-range correlations. Uncertainties in the nuclear transition matrix elements turn out to be less than 15% and 21% for decays involving the emission of single and double Majorons, respectively.
NASA Astrophysics Data System (ADS)
Guseinov, I. I.; Mamedov, B. A.
2011-04-01
In this study, a new method is proposed for evaluating electric multipole transition (radial) matrix elements of the generalized type Hnl,n'l'k in hydrogenic atom and ions using the Slater type orbitals (STOs). The formula obtained allows the determination of all multipole transition matrix elements between two different nonrelativistic radial wave functions Rnl and R. A comparative study carried out between the results of analytical computations and other numerical simulations shows that the methods agree well and emphasizing thus the effectiveness and accuracy of the proposed analytical expressions. The simple equation thus obtained has been found to be remarkable accurate and has shown a wide range of applicability.
Kozono, Shuji; Haraguchi, Hiroki
2007-07-31
The determination of 52 impurity elements in niobium materials (niobium metal, niobium oxide (V), and niobium pentaethoxide) was performed by inductively coupled plasma mass spectrometry (ICP-MS) with on-line anion exchange matrix separation as well as direct nebulization. Niobium material samples were decomposed with a mixture of hydrofluoric acid and nitric acid to prepare 10% niobium solutions. In the on-line anion exchange matrix separation/ICP-MS, the niobium and hydrofluoric acid concentrations in sample solution were adjusted to 5% and ca. 8M, respectively. The solution was then injected into the carrier stream from the sample loop of injection valve to pass through an anion exchange resin column. In the anion exchange separation, niobium in the fluoro-complex form was adsorbed on the resin, while impurity elements were eluted. The eluted elements were introduced into ICP-MS for the determination of 25 impurity elements. On the other hand, 27 impurity elements could not be separated well from niobium matrix under the above anion exchange conditions, and then the sample solution with the niobium concentration of max. 0.2% containing internal standard elements was injected from the sample loop of injection valve directly to introduce into ICP-MS. As a result, 52 impurity elements in three kinds of niobium materials could be determined at the ng g(-1) level. PMID:19071834
Kaon matrix elements and CP violation from lattice QCD with 2+1 flavors of domain wall fermions
NASA Astrophysics Data System (ADS)
Li, Shu
Low energy constants describing the weak, two-pion decays of K mesons in chiral perturbation theory are computed using 2+1 flavors of domain wall fermions in a finite volume with spatial extent 2.74 fm and for a single inverse lattice spacing 1/a = 1.73 GeV. Partially quenched perturbation theory is used in both leading and next-to-leading order. The non-perturbative regularization independent RI/MOM renormalization scheme is employed to determine these low energy constants in the continuum, RI normalization scheme with 20% statistical errors but systematic errors which are estimated to lie between 50 and 100% depending on the operator. These low energy constants are then used to estimate the DeltaI = 1/2 and DeltaI = 3/2 K → pipi decay matrix elements and epsilon'/epsilon. The poor convergence of chiral perturbation theory for quark masses as large as that of the strange quark severely limits the accuracy of these results.
Sahoo, Bijaya K.; Gopakumar, Geetha; Chaudhuri, Rajat K.; Das, B.P.; Merlitz, Holger; Mahapatra, Uttam Sinha; Mukherjee, Debashis
2003-10-01
The relativistic coupled-cluster method is applied to calculate the magnetic dipole hyperfine constant 'A' of the 6s{sub 1/2}, 6p{sub 1/2}, 6p{sub 3/2}, and 5d{sub 3/2} states of singly ionized barium. After the inclusion of two-body correlation effects into the computation of the hyperfine matrix elements, the accuracy of the obtained values was significantly increased compared to earlier computations. Based on these numbers and earlier calculations of the electric dipole transitions and excitation energies, an estimate for the accuracy of the vertical bar [5p{sup 6}]6s{sub 1/2}>{yields} vertical bar [5p{sup 6}]5d{sub 3/2}> parity-nonconserving electric dipole transition amplitude is carried out. The results suggest that for the first time, to our knowledge, a precision of better than 1% is feasible for this transition amplitude.
YAMAZAKI, T.
2006-07-23
We calculate {Delta}I = 3/2 kaon decay matrix elements using domain wall fermions and the DBW2 gauge action at one coarse lattice spacing corresponding to a{sup -1} = 1.3 GeV. We employ the Lellouch and Luescher formula and its extension for non-zero total momentum to extract the infinite volume, center-of-mass frame decay amplitudes. The decay amplitudes obtained from the methods correspond to those from the indirect method with full order chiral perturbation theory. We confirm that the result is consistent with the previous result calculated with H-parity (anti-periodic) boundary condition by investigating the relative momentum dependence. We evaluate the decay amplitude ReA{sub 2} at the physical point by a chiral extrapolation with a polynomial function of m{sub {pi}}{sup 2} and the relative momentum as well as the {Delta}I = 3/2 electroweak penguin contributions to {var_epsilon}{prime}/{var_epsilon}. We found that the result of ReA{sub 2} reasonably agrees with the experiment.
Freeman, John
2007-01-01
A measurement of the top quark mass in t$\\bar{t}$ → l + jets candidate events, obtained from p$\\bar{p}$ collisions at √s = 1.96 TeV at the Fermilab Tevatron using the CDF II detector, is presented. The measurement approach is that of a matrix element method. For each candidate event, a two dimensional likelihood is calculated in the top pole mass and a constant scale factor, 'JES', where JES multiplies the input particle jet momenta and is designed to account for the systematic uncertainty of the jet momentum reconstruction. As with all matrix element techniques, the method involves an integration using the Standard Model matrix element for t$\\bar{t}$ production and decay. However, the technique presented is unique in that the matrix element is modified to compensate for kinematic assumptions which are made to reduce computation time. Background events are dealt with through use of an event observable which distinguishes signal from background, as well as through a cut on the value of an event's maximum likelihood. Results are based on a 955 pb^{-1} data sample, using events with a high-p_{T} lepton and exactly four high-energy jets, at least one of which is tagged as coming from a b quark; 149 events pass all the selection requirements. They find M_{meas} = 169.8 ± 2.3(stat.) ± 1.4(syst.) GeV/c^{2}.
Freeman, John C
2007-01-01
A measurement of the top quark mass in t$\\bar{t}$ → l + jets candidate events, obtained from p$\\bar{p}$ collisions at √s = 1.96 TeV at the Fermilab Tevatron using the CDF II detector, is presented. The measurement approach is that of a matrix element method. For each candidate event, a two dimensional likelihood is calculated in the top pole mass and a constant scale factor, 'JES', where JES multiplies the input particle jet momenta and is designed to account for the systematic uncertainty of the jet momentum reconstruction. As with all matrix elements techniques, the method involves an integration using the Standard Model matrix element for tt production and decay. however, the technique presented is unique in that the matrix element is modified to compensate for kinematic assumptions which are made to reduce computation time. Background events are dealt with through use of an event observable which distinguishes signal from background, as well as through a cut on the value of an event's maximum likelihood. Results are based on a 955 pb^{-1 }data sample, using events with a high-p_{T} lepton and exactly four high-energy jets, at least one of which is tagged as coming from a b quark; 149 events pass all the selection requirements. They find M_{meas} = 169.8 ± 2.3(stat.) ± 1.4(syst.) GeV/c^{2}.
Wang, Xiao-Ping; Zhang, Ji-Long
2007-07-01
Twelve camphor (cinnamomum camphora) tree bark samples were collected from Hiroshima and Kyoto, and the matrix element composition and morphology of the outer surface of these camphor tree bark samples were studied by EDXS and SEM respectively. After a dry decomposition, DOWEX 1-X8 anion exchange resin was used to separate uranium from matrix elements in these camphor tree bark samples. Finally, 235U/238 U isotope ratios in purified uranium solutions were determined by MC-ICP-MS. It was demonstrated that the outer surface of these camphor tree bark samples is porous and rough, with Al, Ca, Fe, K, Mg, Si, C, O and S as its matrix element composition. Uranium in these camphor tree bark samples can be efficiently separated and quantitatively recovered from the matrix element composition. Compared with those collected from Kyoto, the camphor tree bark samples collected from Hiroshima have significantly higher uranium contents, which may be due to the increased aerosol mass concentration during the city reconstruction. Moreover, the 235 U/23.U isotope ratios in a few camphor tree bark samples collected from Hiroshima are slightly higher than 0.007 25. PMID:17944430
Knupp, P.M.
1999-01-18
Structured mesh quality optimization methods are extended to optimization of unstructured triangular, quadrilateral, and mixed finite element meshes. N"ew interpretations of well-known nodally-bssed objective functions are made possible using matrices and matrix norms. The matrix perspective also suggests several new objective functions. Particularly significant is the interpretation of the Oddy metric and the Smoothness objective functions in terms of the condition number of the metric tensor and Jacobian matrix, respectively. Objective functions are grouped according to dimensionality to form weighted combinations. A simple unconstrained local optimum is computed using a modiiied N-ewton iteration. The optimization approach was implemented in the CUBIT mesh generation code and tested on several problems. Results were compared against several standard element-based quaIity measures to demonstrate that good mesh quality can be achieved with nodally-based objective functions.
NASA Technical Reports Server (NTRS)
Ray, Sujit K.; Utku, Senol; Wada, Ben K.
1986-01-01
The stiffness-matrix formulation for the rectangular finite element described by Melosh (1963) and Weaver and Johnston (1984) is generalized to orthotropic materials with material axes not necessarily coincident with the x and y axes; i.e., the condition d(13) = d(23) = 0 is removed. Also included are explicit expressions for the element load vector associated with nonuniform temperature increase in the element. Applications to the analysis of thermal stresses in thin Si-crystal ribbons subjected to temperature changes with highly nonuniform lengthwise and transverse gradients (Utku et al., 1986) and to the simulation of the thermoviscoelastic behavior of growing Si ribbons (Utku and Ray, 1986) are indicated.
Knupp, P.M.
1999-03-26
Three-dimensional unstructured tetrahedral and hexahedral finite element mesh optimization is studied from a theoretical perspective and by computer experiments to determine what objective functions are most effective in attaining valid, high quality meshes. The approach uses matrices and matrix norms to extend the work in Part I to build suitable 3D objective functions. Because certain matrix norm identities which hold for 2 x 2 matrices do not hold for 3 x 3 matrices. significant differences arise between surface and volume mesh optimization objective functions. It is shown, for example, that the equivalence in two-dimensions of the Smoothness and Condition Number of the Jacobian matrix objective functions does not extend to three dimensions and further. that the equivalence of the Oddy and Condition Number of the Metric Tensor objective functions in two-dimensions also fails to extend to three-dimensions. Matrix norm identities are used to systematically construct dimensionally homogeneous groups of objective functions. The concept of an ideal minimizing matrix is introduced for both hexahedral and tetrahedral elements. Non-dimensional objective functions having barriers are emphasized as the most logical choice for mesh optimization. The performance of a number of objective functions in improving mesh quality was assessed on a suite of realistic test problems, focusing particularly on all-hexahedral ''whisker-weaved'' meshes. Performance is investigated on both structured and unstructured meshes and on both hexahedral and tetrahedral meshes. Although several objective functions are competitive, the condition number objective function is particularly attractive. The objective functions are closely related to mesh quality measures. To illustrate, it is shown that the condition number metric can be viewed as a new tetrahedral element quality measure.
NASA Astrophysics Data System (ADS)
Amato, Fulvio; Nava, Silvia; Lucarelli, Franco; Querol, Xavier; Alastuey, Andrés.
2010-05-01
Despite the high environmental and health burden of road dust emissions in urban environments, there is still a dearth of knowledge on the effectiveness of some possible remediation measures such as street cleaning activities. As a consequence of the recent notification from the EU commission for the exceedances of PM limit values (1999/30/EC), several EU countries (Austria, Germany and UK among others) have introduced street cleaning as one of the main measures to be taken in order to meet these limits in the next future. Nevertheless, the effectiveness of street cleaning is still far from being definitively determined since only few tests have been carried out so far and with very different local conditions. An intensive campaign was carried out during spring 2009 in the city centre of Barcelona (NE of Spain) by means of the application of innovative techniques: i) the hourly elemental composition of size segregated PM was used to investigate short term variability of specific tracers of road dust resuspension; ii) a Positive Matrix Factorization was applied in order to identify the daily pattern of each PM source. Hourly elemental concentrations were obtained by a two-stage streaker sampler, where particles are separated on different stages: an impactor deposits the aerosol coarse fraction (aerodynamic diameter between 2.5-10 µm) on a Kapton foil while the fine fraction (<2.5 µm) is collected on a Nuclepore filter having 0.4 µm pores. The two collecting plates (Kapton and Nuclepore) are paired on a cartridge which rotates at constant speed for a week: this produces a circular continuous deposition of particulate matter (streak) on both stages. Totally 349 samples were collected onto three pairs of stages and analyzed by Proton Induced X-Ray Emission (PIXE) external beam facility in Florence, based on a Van de Graaff accelerator. This facility has been used several times in the past for aerosol studies. A Positive Matrix Factorization (PMF) model was applied
Chen, Zhenhua; Chen, Xun; Wu, Wei
2013-04-28
In this series, the n-body reduced density matrix (n-RDM) approach for nonorthogonal orbitals and their applications to ab initio valence bond (VB) methods are presented. As the first paper of this series, Hamiltonian matrix elements between internally contracted VB wave functions are explicitly provided by means of nonorthogonal orbital based RDM approach. To this end, a more generalized Wick's theorem, called enhanced Wick's theorem, is presented both in arithmetical and in graphical forms, by which the deduction of expressions for the matrix elements between internally contracted VB wave functions is dramatically simplified, and the matrix elements are finally expressed in terms of tensor contractions of electronic integrals and n-RDMs of the reference VB self-consistent field wave function. A string-based algorithm is developed for the purpose of evaluating n-RDMs in an efficient way. Using the techniques presented in this paper, one is able to develop new methods and efficient algorithms for nonorthogonal orbital based many-electron theory much easier than by use of the first quantized formulism. PMID:23635123
NASA Astrophysics Data System (ADS)
Chen, Zhenhua; Chen, Xun; Wu, Wei
2013-04-01
In this series, the n-body reduced density matrix (n-RDM) approach for nonorthogonal orbitals and their applications to ab initio valence bond (VB) methods are presented. As the first paper of this series, Hamiltonian matrix elements between internally contracted VB wave functions are explicitly provided by means of nonorthogonal orbital based RDM approach. To this end, a more generalized Wick's theorem, called enhanced Wick's theorem, is presented both in arithmetical and in graphical forms, by which the deduction of expressions for the matrix elements between internally contracted VB wave functions is dramatically simplified, and the matrix elements are finally expressed in terms of tensor contractions of electronic integrals and n-RDMs of the reference VB self-consistent field wave function. A string-based algorithm is developed for the purpose of evaluating n-RDMs in an efficient way. Using the techniques presented in this paper, one is able to develop new methods and efficient algorithms for nonorthogonal orbital based many-electron theory much easier than by use of the first quantized formulism.
Vernarsky, Brian J.
2014-01-01
In an effort towards a ''complete'' experiment for the ω meson, we present studies from an experiment with an unpolarized target and a circularly polarized photon beam (g1c), carried out using the CEBAF Large Acceptance Spectrometer (CLAS) at Jefferson Lab. The experiment was analyzed using an extended maximum likelihood fit with partial wave amplitudes. New likelihood functions were calculated to account for the polarization of the photon beam. Both circular and linear polarizations are explored. The results of these fits are then used to project out the spin density matrix for the {omega}. First measurements of the {rho}{sup 3} spin density matrix elements will be presented using this method.
NASA Astrophysics Data System (ADS)
Newell, Kenneth James
A three year research effort is completed with the development of the Discrete Element Consolidation Analyzer (DECA) for process modeling the formation of titanium composites from powder-fiber monotapes. The primary goal of the DECA process model is to provide a statistically realistic analysis of the various physical processes necessary to achieve higher quality composites from the powder-fiber technique. Over the course of this effort, research and code development was conducted in three distinct stages. The first stage focused on the simulation of initial geometry of the powder and fibers as well as the evolution of tape configuration during the pre-consolidation processing steps. The second stage developed the mechanics of the discrete element powder consolidation and the material characterization methods necessary to model the viscoplastic response of the powder to transient thermal and mechanical boundary conditions. The final stage incorporated the presence of fibers to evaluate the interaction mechanics and possible fibers damage resulting from discrete powder-fiber contacts. As a conclusion to the research, DECA model predictions of density versus time for various consolidation profiles are directly compared to actual consolidation test results and a DECA prescribed process profile is used to fabricate a 6sp{''} × 6sp{''} composite panel of Ti-6242/SCS-6. In completing this research, the discrete element modeling technique has proven to be a powerful tool for the analysis and simulation of metal powder consolidation as well as the consolidation of metal matrix composites. The DECA code orchestrates the use of particle kinetics, some simple aspects of gas dynamics, elasticity, plasticity, creep and various innovative material characterization methods to produce a seamless analysis for powder metallurgy processing of composites. Through the application of the DECA capability, many aspects of the processing stages have been elucidated for further
Jaus, W.; Rasche, G. )
1990-01-01
We calculate nuclear-structure corrections to the {ital ft} values of the eight accurately measured superallowed {beta}{sup +} decays. The statistical fit for the average {ital ft} value is very good. The resulting new value for the matrix element of the Kobayashi-Maskawa (KM) matrix is {vert bar}{ital V}{sub {ital ud}}{vert bar}=0.9735(5). The error in {vert bar}{ital V}{sub {ital ud}}{vert bar} has thus been reduced by 50%. Combining this value for {vert bar}{ital V}{sub {ital ud}}{vert bar} with the presently accepted results from kaon-, hyperon-, and {ital B}-decay constraints, the unitarity of the KM matrix for three generations of quarks seems to be violated.