SEARCH · Engineering Papers
Results for “DIRAC OPERATORS”
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.
Neutralizing topological obstructions to bubbles of nothing
Theories with compact extra dimensions can exhibit a vacuum instability known as a bubble of nothing. These decay modes can be obstructed if the internal manifold is stabilized by fluxes, or if it carries Wilson lines for background gauge fields, or if the instanton is incompatible with the spin structure. In each of these cases the decay can proceed by adding dynamical charged membranes or gauge fields. We give a general, bottom-up procedure for constructing approximate bubble of nothing solutions in models with internal spheres stabilized by flux and study the influence of the brane tension on the tunneling exponent, finding two branches of solutions that merge at a minimal superextremal value of the tension. In the case of Wilson operators and incompatible fermions, the relevant bubble is shown to be the Euclidean Reissner-Nordstrom black hole, and the ordinary decay exponent is modified by 1/g 2 effects. We examine the Dirac operator on this background and comment on the relevance for models of supergravity with gauged R-symmetry.
Microscopic origin of criticality at macroscale in QCD chiral phase transition
We reveal that the criticality of the chiral phase transition in QCD at the macroscale arises from the microscopic energy levels of its fundamental constituents, the quarks. We establish a novel relation between cumulants of the chiral order parameter (i.e., chiral condensate) and correlations among the energy levels of quarks (i.e., eigenspectra of the massless Dirac operator), which naturally leads to a generalization of the Banks-Casher relation. Based on this novel relation and through (2+1)-flavor lattice QCD calculations using the HISQ action with varying light quark masses in the vicinity of the chiral phase transition, we demonstrate that the correlations among the infrared part of the Dirac eigenspectra exhibit same universal scaling behaviors as expected of the cumulants of the chiral condensate. We find that these universal scaling behaviors extend up to the physical values of the up and down quark masses.
Optimizing a Multiple Right-Hand Side Dslash Kernel for Intel Knights Corner
Lattice Quantum Chromo-dynamics (LQCD) stands unique as the only computationally tractable, non-perturbative, and model-independent quantum field theory of the strong nuclear force. LQCD achieves tractability by discretizing our familiar 4-dimensional space-time onto a hypercubic lattice. The core component of the linear solver used in LQCD is the Wilson-Dirac operator whose finite difference component is called the Dslash operator. Dslash is a nearest neighbor stencil operator summing matrix-vector multiplications over lattice points. The arithmetic intensity of a regular Dslash operator acting on a single vector in single precision is between 0.88-2.29 (depending on the tricks employed) making it memory bandwidth bound on most current architectures, including Xeon Phi. Maximizing the performance of Dslash on Xeon Phi requires deliberate exploitation of its architectural features. One key feature of Xeon Phi is its powerful vector processing unit with lengthy vector registers. There is a class of LQCD computations which aims to solve the Dirac equation for the same lattice gauge fields multiple times. In this case one can gain additional reuse of the gauge field, and potentially simplified vectorization, through applying Dslash to several vectors at the same time; often referred to as Multiple Right Hand Side (MRHS) Dslash. Previous work in this area focused on applying this approach to the Highly Improved Staggered operator, and in cases to the multi-Grid algorithm for Wilson clover fermions. Here, we present two approaches to MRHS Dslash: a vector register blocking approach and one using the software package QPhiX with a custom Code Generator for low-level intrinsics. We explore the Knights Corner instruction set, memory hierarchy, and threading capabilities. For standard lattice sizes used in testing, we observed significant speedups, with sustained performance of over 700 GFLOPS (single precision) in one instance. We achieved up to 27% of theoretical peak performance compared to a maximum of 12% obtained by the previous single RHS method in QPhiX.
Spawning rings of exceptional points out of Dirac cones
The Dirac cone underlies many unique electronic properties of graphene and topological insulators, and its band structure—two conical bands touching at a single point—has also been realized for photons in waveguide arrays, atoms in optical lattices, and through accidental degeneracy. Deformation of the Dirac cone often reveals intriguing properties; an example is the quantum Hall effect, where a constant magnetic field breaks the Dirac cone into isolated Landau levels. A seemingly unrelated phenomenon is the exceptional point, also known as the parity–time symmetry breaking point, where two resonances coincide in both their positions and widths. Exceptional points lead to counter-intuitive phenomena such as loss-induced transparency, unidirectional transmission or reflection, and lasers with reversed pump dependence or single-mode operation. Dirac cones and exceptional points are connected: it was theoretically suggested that certain non-Hermitian perturbations can deform a Dirac cone and spawn a ring of exceptional points. Here we experimentally demonstrate such an ‘exceptional ring’ in a photonic crystal slab. Angle-resolved reflection measurements of the photonic crystal slab reveal that the peaks of reflectivity follow the conical band structure of a Dirac cone resulting from accidental degeneracy, whereas the complex eigenvalues of the system are deformed into a two-dimensional flat band enclosed by an exceptional ring. This deformation arises from the dissimilar radiation rates of dipole and quadrupole resonances, which play a role analogous to the loss and gain in parity–time symmetric systems. In this study, our results indicate that the radiation existing in any open system can fundamentally alter its physical properties in ways previously expected only in the presence of material loss and gain.
Multigrid deflation for Lattice QCD
Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.
Energy window stochastic density functional theory
We report that linear scaling density functional theory is important for understanding electronic structure properties of nanometer scale systems. Recently developed stochastic density functional theory can achieve linear or even sublinear scaling for various electronic properties without relying on the sparsity of the density matrix. The basic idea relies on projecting stochastic orbitals onto the occupied space by expanding the Fermi-Dirac operator and repeating this for N χ stochastic orbitals. Often, a large number of stochastic orbitals are required to reduce the statistical fluctuations (which scale as N$-1/2\atop{χ}$) below a tolerable threshold. In this work, we introduce a new stochastic density functional theory that can efficiently reduce the statistical fluctuations for certain observable and can also be integrated with an embedded fragmentation scheme. The approach is based on dividing the occupied space into energy windows and projecting the stochastic orbitals with a single expansion onto all windows simultaneously. This allows for a significant reduction of the noise as illustrated for bulk silicon with a large supercell. Finally, we also provide theoretical analysis to rationalize why the noise can be reduced only for a certain class of ground state properties, such as the forces and electron density.
Chiral anomalies and Wilson fermions
The Wilson formulation of fermions in lattice gauge theory provides a unified description of the chiral anomalies in the standard model. The discrete Dirac operator diagonalizes into a series of 2 × 2 blocks. In each block the possible eigenvalues either form a complex pair or separate into two real eigenvalues that have specific chirality. The collision of these pairs of eigenvalues occurs outside the perturbative region and provides a path between topological sectors.
Dirac spectral density in N f = 2 + 1 QCD at T = 230 MeV
We compute the renormalized Dirac spectral density in N f = 2 + 1 QCD at physical quark masses, temperature T = 230 MeV , and system size L s = 3.4 fm . To that end, we perform a pointwise continuum limit of the staggered density in lattice QCD with staggered quarks. We find, for the first time, that a clear infrared structure (IR peak) emerges in the density of the Dirac operator describing dynamical quarks. We also provide numerical evidence that a component of this peak, which becomes dominant in the thermodynamic limit, is due to a nontrivial accumulation of near-zero modes. Features of this structure are consistent with those previously attributed to the recently proposed IR phase of thermal QCD. Our results (i) provide the only complete first-principles evidence that these IR features exist and are physical, (ii) improve the upper bound for IR phase transition temperature T IR so that the new window is 200 < T IR < 230 MeV , and (iii) are consistent with the nonrestoration of anomalous U A ( 1 ) symmetry (chiral limit) below T = 230 MeV . Published by the American Physical Society 2024
Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors
In this study, we demonstrate the applicability of a recently proposed multiscale thermalization algorithm to two-color quantum chromodynamics (QCD) with two mass-degenerate fermion flavors. The algorithm involves refining an ensemble of gauge configurations that had been generated using a renormalization group (RG) matched coarse action, thereby producing a fine ensemble that is close to the thermalized distribution of a target fine action; the refined ensemble is subsequently rethermalized using conventional algorithms. Although the generalization of this algorithm from pure Yang-Mills theory to QCD with dynamical fermions is straightforward, we find that in the latter case, the method is susceptible to numerical instabilities during the initial stages of rethermalization when using the hybrid Monte Carlo algorithm. We find that these instabilities arise from large fermion forces in the evolution, which are attributed to an accumulation of spurious near-zero modes of the Dirac operator. We propose a simple strategy for curing this problem, and demonstrate that rapid thermalization--as probed by a variety of gluonic and fermionic operators--is possible with the use of this solution. Also, we study the sensitivity of rethermalization rates to the RG matching of the coarse and fine actions, and identify effective matching conditions based on a variety of measured scales.
Strange and charm quark spins from the anomalous Ward identity
In this article, we present a calculation of the strange and charm quark contributions to the nucleon spin from the anomalous Ward identity (AWI). This is performed with overlap valence quarks on 2 + 1 -flavor domain-wall fermion gauge configurations on a 2 4 3 × 64 lattice with lattice spacing a - 1 = 1.73 GeV and the light sea mass at m π = 330 MeV . To satisfy the AWI, the overlap fermion for the pseudoscalar density and the overlap Dirac operator for the topological density, which do not have multiplicative renormalization, are used to normalize the form factor of the local axial-vector current at finite q 2 . For the charm quark, we find that the negative pseudoscalar term almost cancels the positive topological term. For the strange quark, the pseudoscalar term is less negative than that of the charm. By imposing the AWI, the strange g A ( q 2 ) at q 2 = 0 is obtained by a global fit of the pseudoscalar and the topological form factors, together with g A ( q 2 ) and the induced pseudoscalar form factor h A ( q 2 ) at finite q 2 . The chiral extrapolation to the physical pion mass gives Δ s + Δ s ¯ = - 0.0403 ( 44 ) ( 78 ) .
The 't Hooft vertex for staggered fermions and flavor-singlet mesons
We derive the ’t Hooft vertex for staggered fermions and examine its symmetries for nonzero lattice spacing. We also derive a set of structural properties for the eigenvectors of the staggered Dirac operator, which should emerge in the continuum limit, if staggered fermions yield four species. This property also is needed for flavor-taste-singlet correlators to behave correctly. We then test numerically whether the needed structure arises: it does. This structure and symmetry of (unrooted) staggered fermions also imply that Creutz’s (latest) objections to the rooted determinant are without foundation.
Search on a hypercubic lattice using a quantum random walk. II. d=2
We investigate the spatial search problem on the two-dimensional square lattice, using the Dirac evolution operator discretized according to the staggered lattice fermion formalism. d=2 is the critical dimension for the spatial search problem, where infrared divergence of the evolution operator leads to logarithmic factors in the scaling behavior. As a result, the construction used in our accompanying article [A. Patel and M. A. Rahaman, Phys. Rev. A 82, 032330 (2010)] provides an O({radical}(N)lnN) algorithm, which is not optimal. The scaling behavior can be improved to O({radical}(NlnN)) by cleverly controlling the massless Dirac evolution operator by an ancilla qubit, as proposed by Tulsi [Phys. Rev. A 78, 012310 (2008)]. We reinterpret the ancilla control as introduction of an effective mass at the marked vertex, and optimize the proportionality constants of the scaling behavior of the algorithm by numerically tuning the parameters.
Hidden symmetries and killing tensors on curved spaces
Higher-order symmetries corresponding to Killing tensors are investigated. The intimate relation between Killing-Yano tensors and nonstandard supersymmetries is pointed out. In the Dirac theory on curved spaces, Killing-Yano tensors generate Dirac-type operators involved in interesting algebraic structures as dynamical algebras or even infinite dimensional algebras or superalgebras. The general results are applied to space-times which appear in modern studies. One presents the infinite dimensional superalgebra of Dirac type operators on the 4-dimensional Euclidean Taub-NUT space that can be seen as a twisted loop algebra. The existence of the conformal Killing-Yano tensors is investigated for some spaces with mixed 3-Sasakian structures.
Search on a hypercubic lattice using a quantum random walk. I. d>2
Random walks describe diffusion processes, where movement at every time step is restricted to only the neighboring locations. We construct a quantum random walk algorithm, based on discretization of the Dirac evolution operator inspired by staggered lattice fermions. We use it to investigate the spatial search problem, that is, to find a marked vertex on a d-dimensional hypercubic lattice. The restriction on movement hardly matters for d>2, and scaling behavior close to Grover's optimal algorithm (which has no restriction on movement) can be achieved. Using numerical simulations, we optimize the proportionality constants of the scaling behavior, and demonstrate the approach to that for Grover's algorithm (equivalent to the mean-field theory or the d{yields}{infinity} limit). In particular, the scaling behavior for d=3 is only about 25% higher than the optimal d{yields}{infinity} value.
On the relation between rigging inner product and master constraint direct integral decomposition
Canonical quantization of constrained systems with first-class constraints via Dirac's operator constraint method proceeds by the theory of Rigged Hilbert spaces, sometimes also called refined algebraic quantization. This method can work when the constraints form a Lie algebra. When the constraints only close with nontrivial structure functions, the Rigging map can no longer be defined. To overcome this obstacle, the master constraint method has been proposed which replaces the individual constraints by a weighted sum of absolute squares of the constraints. Now the direct integral decomposition (DID) methods, which are closely related to Rigged Hilbert spaces, become available and have been successfully tested in various situations. It is relatively straightforward to relate the rigging inner product to the path integral that one obtains via reduced phase space methods. However, for the master constraint, this is not at all obvious. In this paper we find sufficient conditions under which such a relation can be established. Key to our analysis is the possibility to pass to equivalent, Abelian constraints, at least locally in phase space. Then the master constraint DID for those Abelian constraints can be directly related to the rigging map and therefore has a path integral formulation.
Functional determinants of radial operators in AdS 2
We study the zeta-function regularization of functional determinants of Laplace and Dirac-type operators in two-dimensional Euclidean AdS 2 space. More specifically, we consider the ratio of determinants between an operator in the presence of background fields with circular symmetry and the free operator in which the background fields are absent. By Fourier-transforming the angular dependence, one obtains an infinite number of one-dimensional radial operators, the determinants of which are easy to compute. The summation over modes is then treated with care so as to guarantee that the result coincides with the two-dimensional zeta-function formalism. The method relies on some well-known techniques to compute functional determinants using contour integrals and the construction of the Jost function from scattering theory. Here, our work generalizes some known results in flat space. The extension to conformal AdS 2 geometries is also considered. We provide two examples, one bosonic and one fermionic, borrowed from the spectrum of fluctuations of the holographic 1/4-BPS latitude Wilson loop.
Upper bounds on the number of eigenvalues of stationary Schroedinger operators
In this article we provide several upper bounds on the number of eigenvalues, respectively, for Schroedinger operators of one-dimensional case, central potential case, and the case with point-interactions. Compared with Bargman's and Calogero's results, the new bounds can give finite estimates when the potential function has a heavy tail. Besides, a generalization of Calogero's bound is given, which can be applied to bell shaped potential functions rather than merely monotonous ones. With respect to Schroedinger operators containing Dirac functions, some bounds are also established.