Lattice QCD calculations for nuclear physics
I present some recent results regarding the use of numerical Lattice QCD techniques to describe light nuclear systems from the underlying theory of the strong interaction, QCD.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
I present some recent results regarding the use of numerical Lattice QCD techniques to describe light nuclear systems from the underlying theory of the strong interaction, QCD.
The relation between the hadronic short-distance constituent quark and gluon particle limit and the long-range confining domain is yet one of the most challenging aspects of particle physics due to the strong coupling nature of Quantum Chromodynamics, the fundamental theory of the strong interactions. The central question is how one can compute hadronic properties from first principles; i.e., directly from the QCD Lagrangian. The most successful theoretical approach thus far has been to quantize QCD on discrete lattices in Euclidean space-time. Lattice numerical results follow from computation of frame-dependent moments of distributions in Euclidean space and dynamical observables in Minkowski spacetime, such as the time-like hadronic form factors, are not amenable to Euclidean lattice computations. The Dyson-Schwinger methods have led to many important insights, such as the infrared fixed point behavior of the strong coupling constant, but in practice, the analyses are limited to ladder approximation in Landau gauge. Baryon spectroscopy and the excitation dynamics of nucleon resonances encoded in the nucleon transition form factors can provide fundamental insight into the strong-coupling dynamics of QCD. New theoretical tools are thus of primary interest for the interpretation of the results expected at the new mass scale and kinematic regions accessible to the JLab 12 GeV Upgrade Project. The AdS/CFT correspondence between gravity or string theory on a higher-dimensional anti-de Sitter (AdS) space and conformal field theories in physical space-time has led to a semiclassical approximation for strongly-coupled QCD, which provides physical insights into its nonperturbative dynamics. The correspondence is holographic in the sense that it determines a duality between theories in different number of space-time dimensions. This geometric approach leads in fact to a simple analytical and phenomenologically compelling nonperturbative approximation to the full light-front QCD Hamiltonian 'Light-Front Holography'. Light-Front Holography is in fact one of the most remarkable features of the AdS/CFT correspondence. The Hamiltonian equation of motion in the light-front (LF) is frame independent and has a structure similar to eigenmode equations in AdS space. This makes a direct connection of QCD with AdS/CFT methods possible. Remarkably, the AdS equations correspond to the kinetic energy terms of the partons inside a hadron, whereas the interaction terms build confinement and correspond to the truncation of AdS space in an effective dual gravity approximation. One can also study the gauge/gravity duality starting from the bound-state structure of hadrons in QCD quantized in the light-front. The LF Lorentz-invariant Hamiltonian equation for the relativistic bound-state system is P{sub {mu}}P{sup {mu}}|{psi}(P)> = (P{sup +}P{sup -} - P{sub {perpendicular}}{sup 2})|{psi}(P)> = M{sup 2}|{psi}(P)>, P{sup {+-}} = P{sup 0} {+-} P{sup 3}, where the LF time evolution operator P{sup -} is determined canonically from the QCD Lagrangian. To a first semiclassical approximation, where quantum loops and quark masses are not included, this leads to a LF Hamiltonian equation which describes the bound-state dynamics of light hadrons in terms of an invariant impact variable {zeta} which measures the separation of the partons within the hadron at equal light-front time {tau} = x{sup 0} + x{sup 3}. This allows us to identify the holographic variable z in AdS space with an impact variable {zeta}. The resulting Lorentz-invariant Schroedinger equation for general spin incorporates color confinement and is systematically improvable. Light-front holographic methods were originally introduced by matching the electromagnetic current matrix elements in AdS space with the corresponding expression using LF theory in physical space time. It was also shown that one obtains identical holographic mapping using the matrix elements of the energy-momentum tensor by perturbing the AdS metric around its static solution. A gravity dual to QCD is not known, but the mechanisms of confinement can be incorporated in the gauge/gravity correspondence by modifying the AdS geometry in the large infrared (IR) domain z {approx} 1 = {Lambda}{sub QCD}, which also sets the scale of the strong interactions. In this simplified approach we consider the propagation of hadronic modes in a fixed effective gravitational background asymptotic to AdS space, which encodes salient properties of the QCD dual theory, such as the ultraviolet (UV) conformal limit at the AdS boundary, as well as modifications of the background geometry in the large z IR region to describe confinement. The modified theory generates the point-like hard behavior expected from QCD, instead of the soft behavior characteristic of extended objects.
Type IIn supernovae (SNe) show spectral evidence for strong interaction between their blast wave and dense circumstellar material (CSM) around the progenitor star. SN 2010jl was the brightest core-collapse supernova in 2010, and it was a Type IIn explosion with strong CSM interaction. Andrews et al. recently reported evidence for an infrared (IR) excess in SN 2010jl, indicating either new dust formation or the heating of CSM dust in an IR echo. Here we report multi-epoch spectra of SN 2010jl that reveal the tell-tale signature of new dust formation: emission-line profiles becoming systematically more blueshifted as the red side of the line is blocked by increasing extinction. The effect is seen clearly in the intermediate-width (400-4000 km s{sup -1}) component of H{alpha} beginning roughly 30 days after explosion. Moreover, we present near-IR spectra demonstrating that the asymmetry in the hydrogen-line profiles is wavelength dependent, appearing more pronounced at shorter wavelengths. This evidence suggests that new dust grains had formed quickly in the post-shock shell of SN 2010jl arising from CSM interaction. Since the observed dust temperature has been attributed to an IR echo and not to new dust, either (1) IR excess emission at {lambda} < 5 {mu}m is not a particularly sensitive tracer of new dust formation in SNe, or (2) some assumptions about expected dust temperatures might require further study. Lastly, we discuss one possible mechanism other than dust that might lead to increasingly blueshifted line profiles in SNe IIn, although the wavelength dependence of the asymmetry argues against this hypothesis in the case of SN 2010jl.
Abstract—Superconducting heterostructures interlayered with a material with strong spin–orbit interaction are studied. As materials with strong spin–orbit interaction, we choose strontium iridates featuring the Ruddlesden–Popper structure (Sr{sub n} {sub + 1}Ir{sub n}O{sub 3n} {sub + 1}; n = 1, 2, ∞): paramagnetic metal SrIrO{sub 3} and magnetic insulator Sr{sub 2}IrO{sub 4}. Using laser ablation, epitaxial films of these materials are obtained, as well as heterostructures based on epitaxial films of cuprate superconducting YBa{sub 2}Cu{sub 3}O{sub x}. For electrical measurements, mesa-structures of micrometer sizes with top electrodes made of Au/Nb bilayer are fabricated.
We analyze the observed properties of dwarf galaxies, which are dark matter dominated astrophysical objects, by assuming that dark matter is in the form of a strongly-coupled, dilute Bose-Einstein condensate. The basic astrophysical properties of the condensate (density profile, rotational velocity, and mass profile, respectively), are derived from a variational principle. To test the validity of the model we compare first the tangential velocity equation of the model with a sample of eight rotation curves of dwarf galaxies. We find a good agreement between the theoretically predicted rotation curves (without any baryonic component) and the observational data. The mean value of the logarithmic inner slope of the mass density profile of dwarf galaxies is also obtained, and it is shown that the observed value of this parameter is in agreement with the theoretical results. The predictions of the Bose-Einstein condensate model are also systematically compared with the predictions of the standard Cold Dark Matter model. The non-singular density profiles of the Bose-Einstein condensed dark matter generally show the presence of an extended core, whose presence is due to the strong interaction between dark matter particles.
A two-dimensional Fourier transform of hadron form factors allows to determine their charge density in transverse space. We show that this method can be applied to any virtual photon induced transition, such as {gamma}{sup *}(q)+N{yields}{pi}N. Only Fock states that are common to the initial and final states contribute to the amplitudes, which are determined by the overlap of the corresponding light-front wave functions. Their transverse extent may be studied as a function of the final state configuration, allowing qualitatively new insight into strong interaction dynamics. Fourier transforming the cross section (rather than the amplitude) gives the distribution of the transverse distance between the virtual photon interaction vertices in the scattering amplitude and its complex conjugate. While the measurement of parton distributions in longitudinal momentum depends on the leading twist approximation (-q{sup 2}{yields}{infinity} limit), all q{sup 2}<0 values contribute to the Fourier transform, with the transverse resolution increasing with the available range in q{sup 2}. We illustrate the method using QED amplitudes.
Unsteady laminar flow with convective heat transfer through a curved square duct rotating at a constant angular velocity about the center of curvature is investigated numerically by using a spectral method, and covering a wide range of the Taylor number −300≤Tr≤1000 for the Dean number Dn = 1000. A temperature difference is applied across the vertical sidewalls for the Grashof number Gr = 100, where the outer wall is heated and the inner wall cooled, the top and bottom walls being adiabatic. Flow characteristics are investigated with the effects of rotational parameter, Tr, and the pressure-driven parameter, Dn, for the constant curvature 0.001. Time evolution calculations as well as their phase spaces show that the unsteady flow undergoes through various flow instabilities in the scenario ‘multi-periodic → chaotic → steady-state → periodic → multi-periodic → chaotic’, if Tr is increased in the positive direction. For negative rotation, however, time evolution calculations show that the flow undergoes in the scenario ‘multi-periodic → periodic → steady-state’, if Tr is increased in the negative direction. Typical contours of secondary flow patterns and temperature profiles are obtained at several values of Tr, and it is found that the unsteady flow consists of two- to six-vortex solutions if the duct rotation is involved. External heating is shown to generate a significant temperature gradient at the outer wall of the duct. This study also shows that there is a strong interaction between the heating-induced buoyancy force and the centrifugal-Coriolis instability in the curved channel that stimulates fluid mixing and consequently enhances heat transfer in the fluid.
We investigate all pure quantum-electrodynamics corrections to the np{yields}1s, n=2-4 transition energies of pionic hydrogen larger than 1 meV, which requires an accurate evaluation of all relevant contributions up to order {alpha}{sup 5}. These values are needed to extract an accurate strong interaction shift from experiment. Many small effects, such as second-order and double vacuum polarization contribution, proton and pion self-energies, finite size and recoil effects are included with exact mass dependence. Our final value differs from previous calculations by up to {approx_equal}11 ppm for the 1s state, while a recent experiment aims at a 4 ppm accuracy.
We investigate a single impurity interacting with a free two-dimensional atomic Fermi gas. The interaction between the impurity and the gas is characterized by an arbitrary attractive short-range potential, which, in two dimensions, always admits a two-particle bound state. We provide analytical expressions for the energy and the effective mass of the dressed impurity by including the two-body bound state, which is crucial for strong interactions, in the integral equation for the effective interaction. Using the same method, we also give the results for the polaron parameters in one and three dimensions and find good agreement with previous results. Thus, our relations can be used as a simple way to estimate the polaron parameters once the two-body bound state of the interaction potential is known.
The precision spectroscopy of kaonic hydrogen and kaonic helium x-rays offers a unique possibility to precisely determine the K{sup -}N strong interaction in the low energy limit. Results and status of modern experiments using silicon drift x-ray detectors (SDDs), namely, KEK E570, SIDDHARTA at DA{Phi}NE and E17 at J-PARC, are discussed.
The study of baryon resonances provides a deeper understanding of the strong interaction because the dynamics and relevant degrees of freedom hidden within them are reflected by the properties of the excited states of baryons. Higher-lying excited states at and above 1.9 GeV/c{sup 2} are generally predicted to have strong couplings to the {pi}{pi}N final states via {pi}{Delta} or {rho}N intermediate states. Double-pion photoproduction is therefore important to find and investigate properties of highmass resonances. The CLAS g9a (FROST) experiment, as part of the N* spectroscopy program at Jefferson Laboratory (JLab), has accumulated photoproduction data using linearly- and circularly-polarized photons incident on a longitudinally-polarized butanol target in the photon energy range 0.3 to 2.4 GeV. In this contribution, the extraction of the target asymmetry for the reaction {gamma}p {yields} p{pi}{sup +}{pi}{sup -} will be described and preliminary results will be presented.
Mach-Zehnder atom interferometry requires hold-time phase squeezing to attain readout accuracy below the standard quantum limit. This increases its sensitivity to phase diffusion, restoring shot-noise scaling of the optimal signal-to-noise ratio in the presence of interactions. The contradiction between the preparations required for readout accuracy and robustness to interactions is removed by monitoring Rabi-Josephson oscillations instead of relative-phase oscillations during signal acquisition. Optimizing the signal-to-noise ratio with a Gaussian squeezed input, we find that hold-time number squeezing satisfies both demands and that sub-shot-noise scaling is retained even for strong interactions.
The condensation of half-light half-matter exciton polaritons in semiconductor optical cavities is a striking example of macroscopic quantum coherence in a solid-state platform. Quantum coherence is possible only when there are strong interactions between the exciton polaritons provided by their excitonic constituents. Rydberg excitons with high principal value exhibit strong dipole–dipole interactions in cold atoms. However, polaritons with the excitonic constituent that is an excited state, namely Rydberg exciton polaritons (REPs), have not yet been experimentally observed. Here, we observe the formation of REPs in a single crystal CsPbBr 3 perovskite cavity without any external fields. These polaritons exhibit strong nonlinear behavior that leads to a coherent polariton condensate with a prominent blue shift. Furthermore, the REPs in CsPbBr 3 are highly anisotropic and have a large extinction ratio, arising from the perovskite’s orthorhombic crystal structure. Our observation not only sheds light on the importance of many-body physics in coherent polariton systems involving higher-order excited states, but also paves the way for exploring these coherent interactions for solid-state quantum optical information processing.
Single top quark production at the International Linear Collider (ILC) can be used to obtain high precision measurements of the V{sub tb} Cabibbo-Kobayashi-Maskawa quark-mixing matrix (CKM) element as well as the effective tbW coupling. We have calculated the QCD correction for the cross section in the context of an effective vector boson approximation. Our results show a {approx}10% increase due to the strong interaction.
We employ the QCD sum rules method for description of nucleons in nuclear matter. We show that this approach provides a consistent formalism for solving various problems of nuclear physics. Such nucleon characteristics as the Dirac effective mass m* and the vector self-energy {Sigma}{sub V} are expressed in terms of the in-medium values of QCD condensates. The values of these parameters at saturation density and the dependence on the baryon density and on the neutron-to-proton density ratio is in agreement with the results, obtained by conventional nuclear physics methods. The contributions to m* and {Sigma}{sub V} are related to observables and do not require phenomenological parameters. The scalar interaction is shown to be determined by the pion-nucleon {sigma} term. The nonlinear behavior of the scalar condensate may appear to provide a possible mechanism of the saturation. The approach provided reasonable results for renormalization of the axial coupling constant, for the contribution of the strong interactions to the neutron-proton mass difference and for the behavior of the structure functions of the in-medium nucleon. The approach enables to solve the problems which are difficult or unaccessible for conventional nuclear physics methods. The method provides guidelines for building the nuclear forces. The three-body interactions emerge within the method in a natural way. Their rigorous calculation will be possible in the framework of self-consistent calculation in nuclear matter of the scalar condensate and of the nucleon effective mass m*.
A long-standing goal of nuclear theory is to determine the properties of atomic nuclei based on the fundamental interactions among the protons and neutrons (i.e., nucleons). By adopting nucleon-nucleon (NN), three-nucleon (NNN) and higher-nucleon interactions determined from either meson-exchange theory or QCD, with couplings fixed by few-body systems, we preserve the predictive power of nuclear theory. This foundation enables tests of nature's fundamental symmetries and offers new vistas for the full range of complex nuclear phenomena. Basic questions that drive our quest for a microscopic predictive theory of nuclear phenomena include: (1) What controls nuclear saturation; (2) How the nuclear shell model emerges from the underlying theory; (3) What are the properties of nuclei with extreme neutron/proton ratios; (4) Can we predict useful cross sections that cannot be measured; (5) Can nuclei provide precision tests of the fundamental laws of nature; and (6) Under what conditions do we need QCD to describe nuclear structure, among others. Along with other ab initio nuclear theory groups, we have pursued these questions with meson-theoretical NN interactions, such as CD-Bonn and Argonne V18, that were tuned to provide high-quality descriptions of the NN scattering phase shifts and deuteron properties. We then add meson-theoretic NNN interactions such as the Tucson-Melbourne or Urbana IX interactions. More recently, we have adopted realistic NN and NNN interactions with ties to QCD. Chiral perturbation theory within effective field theory ({chi}EFT) provides us with a promising bridge between QCD and hadronic systems. In this approach one works consistently with systems of increasing nucleon number and makes use of the explicit and spontaneous breaking of chiral symmetry to expand the strong interaction in terms of a dimensionless constant, the ratio of a generic small momentum divided by the chiral symmetry breaking scale taken to be about 1 GeV/c. The resulting NN and NNN interactions, characterized by the order of the expansion retained (e.g. 'next-to-next-to leading order' is NNLO), provide a high-quality fit to the NN data and the A = 3 ground-state (g.s.) properties. The derivations of NN, NNN, etc. interactions within meson-exchange and {chi}EFT are well-established but are not subjects of this review. Our focus is solution of the non-relativistic quantum many-body Hamiltonian that includes these interactions using our no core shell model (NCSM) formalism. In the next section we will briefly outline the NCSM formalism and then present applications, results and extensions in later sections.
A useful definition of orbital degeneracy—form-degeneracy—is introduced, which is distinct from the usual energy-degeneracy: Two canonical spatial orbitals are form-degenerate when the energy expectation value in the restricted Hartree–Fock (RHF) wave function is unaltered upon a two-electron excitation from one of these orbitals to the other. Form-degenerate orbitals tend to have isomorphic electron densities and occur in the highest-occupied and lowest-unoccupied molecular orbitals (HOMOs and LUMOs) of strongly correlated systems. Here, we present a mathematical proof of the existence of a triplet instability in a real or complex RHF wave function of a finite system in the space of real or complex unrestricted Hartree–Fock wave functions when HOMO and LUMO are energy- or form-degenerate. We also show that a singlet instability always exists in a real RHF wave function of a finite system in the space of complex RHF wave functions, when HOMO and LUMO are form-degenerate, but have nonidentical electron densities, or are energy-degenerate. These theorems provide Hartree–Fock-theory-based explanations of Hund’s rule, a singlet instability in Jahn–Teller systems, biradicaloid electronic structures, and a triplet instability during some covalent bond breaking. They also suggest (but not guarantee) the spontaneous formation of a spin density wave (SDW) in a metallic solid. The stability theory underlying these theorems extended to a continuous orbital-energy spectrum proves the existence of an oscillating (nonspiral) SDW instability in one- and three-dimensional homogeneous electron gases, but only at low densities or for strong interactions.
Quantum many-body systems with coherent and controllable interactions enable the realization of novel quantum materials, the quantum simulation of problems that cannot be simulated on classical computers, computational systems that are exponentially faster than existing classical algorithms, and the direct investigation of the conjectured duality between gravitational theories and quantum field theories. This project focusses on quantum phase transitions in strongly interacting many-body spin systems, and in particular on the dynamical evolution near the critical point and across the phase transition. The dynamics of a large quantum system containing more than 50 qubits is investigated in various critical regimes. In particular, the Kibble-Zurek transition has been characterized and critical exponents have been measured. Furthermore, quantum many-body scars have been observed and characterized for the first time. We are investigating the Kibble-Zurek mechanism in the quantum phase transition in one dimension from a disordered system to a system with ordering induced by Rydberg interactions between neighboring atoms. The Kibble-Zurek mechanism describes the formation of topological defects as a system evolves across a phase transitions, and is expected to apply to a large number of problems, including cosmology, where it should describe the structure of domain formation in the early universe.