Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FERMIONS”

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.

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Accelerated impurity solver for DMFT and its diagrammatic extensions

Here, we present ComCTQMC, a GPU accelerated quantum impurity solver. It uses the continuous-time quantum Monte Carlo (CTQMC) algorithm wherein the partition function is expanded in terms of the hybridisation function (CT-HYB). ComCTQMC supports both partition and worm-space measurements, and it uses improved estimators and the reduced density matrix to improve observable measurements whenever possible. ComCTQMC efficiently measures all one and two-particle Green's functions, all static observables which commute with the local Hamiltonian, and the occupation of each impurity orbital. ComCTQMC can solve complex-valued impurities with crystal fields that are hybridized to both fermionic and bosonic baths. Most importantly, ComCTQMC utilizes graphical processing units (GPUs), if available, to dramatically accelerate the CTQMC algorithm when the Hilbert space is sufficiently large. We demonstrate acceleration by a factor of over 600 (100) in a simulation of δ-Pu at 600 K with (without) crystal fields. In easier problems, the GPU offers less impressive acceleration or even decelerates the CTQMC. Here we describe the theory, algorithms, and structure used by ComCTQMC in order to achieve this set of features and level of acceleration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

SIMULATeQCD: A simple multi-GPU lattice code for QCD calculations

The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates to existing software. These challenges have to be balanced against constantly changing personnel. At the same time, there is a wide range of applications for HISQ fermions in QCD studies. This situation encourages the development of software featuring a HISQ action that is flexible, high-performing, open source, easy to use, and easy to adapt. In this technical paper, we explain the design strategy, provide implementation details, list available algorithms and modules, and show key performance indicators for SIMULATeQCD, a simple multi-GPU lattice code for large-scale QCD calculations, mainly developed and used by the HotQCD collaboration. The code is publicly available on GitHub.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Natural top quark condensation (a redux)

The Nambu--Jona-Lasinio (NJL) model involves a pointlike 4-fermion interaction. While it gives a useful description of chiral dynamics (mainly in QCD), it nonetheless omits the crucially important internal wave-function of a two-body bound state, $\phi(r)$. This becomes significant near critical coupling where $\phi(r)$ extends to large distance, leading to dilution and suppression of induced couplings $\propto \phi(0)$, such as the Yukawa and quartic couplings, as well as reduced fine-tuning of a hierarchy. In top quark condensation, where the BEH boson is a $\bar{t}t$ bound state and we have a UV completion such as topcolor, we must go beyond the NJL model and include effects of $\phi(r)$. We provide a formulation of this for the Brout-Englert-Higgs boson, and find that it leads to an extended $\phi(r)$, a significantly reduced and natural composite scale of $M_0 \sim 6$ TeV, a successful prediction for the quartic coupling, $\lambda$, and fine tuning that is reduced to a few percent, providing a compelling candidate solution to the naturalness problem of the BEH boson. The theory is testable and new physics should begin to emerge on the multi-TeV mass scales and possibly accessible to the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

From oversimplified to overlooked: The case for exploring rich dark sectors

The Standard Model (SM) of particle physics provides a very successful description of fundamental particles and their interactions but it is incomplete, as neutrino masses, dark matter and the baryon asymmetry of the Universe indicate. In addition, the origin of masses and of the approximate fundamental symmetries call out for deeper explanations. The quest for a New SM Theory, that extends the SM to a more general theory, is ongoing. For decades the main focus has been on the TeV scale, but despite an impressive theoretical and experimental effort, no hints of new physics at such scale has been found in experiments. Dark sectors provide an interesting alternative to TeV scale extensions of the SM to explain the open questions in particle and astroparticle physics. Going beyond minimal models, rich dark sectors extend the SM to a complex theory with multiple particles and interactions, in analogy to the SM itself. They have a wealth of theoretical and astrophysical/cosmological consequences and can lead to phenomenological signatures that can be markedly different to that of minimal ones. These include short-lived particles and semi-visible decay signatures, as opposed to minimal models where new states are typically long-lived and purely visible or invisible resonances. Given the experimental configurations and analysis strategies, current dark sector searches might miss such signatures. We advocate a dedicated programme of searches for rich dark sectors that overcomes the assumptions on minimality and on the long lifetime of particles and encompasses a broader range of possibilities. Here, we discuss a prototype model that includes a complex structure akin to the SM: multiple generations of fermions charged under a new spontaneously-broken gauge symmetry.

Dark matter↗
DOE OSTI2020

Dilute Fermi gas at fourth order in effective field theory

Using effective field theory methods, we calculate for the first time the complete fourth-order term in the Fermi-momentum or k F a s expansion for the ground-state energy of a dilute Fermi gas. The convergence behavior of the expansion is examined for the case of spin one-half fermions and compared against quantum Monte-Carlo results, showing that the Fermi-momentum expansion is well-converged at this order for |k F a s | ≲ 0.5.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Non-cancellation of the parity anomaly in the strong-field regime of QED 2+1

Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilson fermions in 2 + 1 space-time dimensions. In the static field limit, relevant at early times, we solve the problem analytically. Our results point out the fundamental role of quantum anomalies for strong-field phenomena, relevant for a wide range of condensed matter and high-energy applications, but also for the next generation of gauge theory quantum simulators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

On the effective theory of neutrino-electron and neutrino-quark interactions

We determine the four Fermi effective theory of neutrino interactions within the Standard Model including one-loop electroweak radiative corrections, in combination with the measured muon lifetime and precision electroweak data. Including two-loop matching and three-loop running corrections, we determine lepton coefficients accounting for all large logarithms through relative order $\cal{O}(\alpha \alpha_s)$ and quark coefficients accounting for all large logarithms through ${\cal{O}}(\alpha)$. We present four-fermion coefficients valid in $n_f=3$ and $n_f=4$ flavor quark theories, as well as in the extreme low-energy limit. We relate the coefficients in this limit to neutrino charge radii governing matter effects via forward neutrino scattering on charged particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Gluon parton distribution of the pion from lattice QCD

We present the first determination of the x-dependent pion gluon distribution from lattice QCD using the pseudo-PDF approach. We use lattice ensembles with 2+1+1 flavors of highly improved staggered quarks (HISQ), generated by MILC Collaboration, at two lattice spacings a ≈ 0.12 and 0.15 fm and three pion masses M π ≈ 220, 310 and 690 MeV. We use clover fermions for the valence action and momentum smearing to achieve pion boost momentum up to 2.29 GeV. We find that the dependence of the pion gluon parton distribution on lattice spacing and pion mass is mild. We compare our results from the lightest pion mass ensemble with the determination by JAM and xFitter global fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Emergence of 4 H $J^π$ = 1 - resonance in contact theories

We obtain the s- and p-wave low-energy scattering parameters for n 3 H elastic scattering and the position of the 4 H J π = 1 - resonance using the pionless effective field theory at leading order. Results are extracted with three numerical techniques: confining the system in a harmonic oscillator trap, solving the Faddeev-Yakubovsky equations in configuration space, and using an effective two-body cluster approach. The renormalization of the theory for the relevant amplitudes is assessed in a cutoff-regulator range between 1 fm -1 and 10 fm -1 . Most remarkably, we find a cutoff-stable/RG-invariant resonance in the 4 H J π = 1 - system. This p-wave resonance is a universal consequence of a shallow two-body state and the introduction of a three-body s-wave scale set by the triton binding energy. The stabilization of a resonant state in a few-fermion system through pure contact interactions has a significant consequence for the powercounting of the pionless theory. Specifically, it suggests the appearance of similar resonant states also in larger nuclei, like 16-oxygen, in which the theory’s leading order does not predict stable states. Those resonances would provide a starting state to be moved to the correct physical position by the perturbative insertion of sub-leading orders, possibly resolving the discrepancy between data and contact EFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum state preparation with resolution refinement

We introduce a method called resolution refinement that allows one to bootstrap eigenstate preparation on a quantum computer. We first prepare an eigenstate of a low-resolution Hamiltonian using any method of choice. The eigenstate is then lifted to higher resolution and adiabatically evolved to produce the corresponding eigenstate of a higher-fidelity Hamiltonian. We give examples of resolution refinement applied to both single-particle basis states as well as a spatial lattice grid. For basis refinement, we compute few-body ground states of the Busch model for interacting particles in a harmonic trap in one dimension. For lattice refinement, we compute Hartree-Fock nuclear states for a central Woods-Saxon potential in three dimensions, and we compute bound states and continuum states in a multi-species Hubbard model of fermions in one dimension. In all cases, the method is efficient and requires an adiabatic evolution time that scales with the inverse of the energy gap times the square root of the system size. We show that this very favorable scaling arises from the fact that resolution refinement does not make large changes to the structure or energies of the low-energy eigenstates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Lattice QCD Inputs for nuclear double beta decay

Second order β -decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double- β decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double- β decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double- β decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double- β decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn → ppe^-e^-\bar{v}_e\bar{v}_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double- β decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double- β decay processes themselves or using inputs from LQCD. Finally, future prospects for direct calculations in larger nuclei are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Unexpected Hydride: Ce 4 B 2 C 2 H 2.42 , a Stuffed Variant of the Nd 2 BC Structure Type

Ce 4 B 2 C 2 H 2.42 was grown as large crystals from a cerium/copper eutectic flux. The structure was characterized by single-crystal X-ray and neutron diffraction and was found to be a stuffed variant of Nd 2 BC with the addition of two interstitial hydrogen positions. The tetrahedral hydrogen position is fully occupied, while the octahedral position has an occupancy of 42(3)%. Initial synthesis was due to hydrogen contamination of the cerium metal but has been successfully repeated using anthracene as a carbon and hydrogen source. Density of states calculations suggest that the incorporation of hydrogen stabilizes the compound with respect to the nonhydrided model. Magnetic susceptibility data show a complex magnetic ordering at 7.7 K that originates from the localized electron on the Ce 3+ in the structure. The trivalent state is also supported by X-ray photoelectron spectroscopy measurements. Heat capacity and electrical resistivity data show that the phase transition is broad in temperature, which may be due to structural disorder. Furthermore, the large low temperature value of C/T also indicates possible heavy fermion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Study of Polycrystalline Bulk Sr 3 OsO 6 Double-Perovskite Insulator: Comparison with 1000 K Ferromagnetic Epitaxial Films

Polycrystalline Sr 3 OsO 6 , which is an ordered double-perovskite insulator, is synthesized via solid-state reaction under high-temperature and high-pressure conditions of 1200 ° C and 6 GPa. The synthesis enables us to conduct a comparative study of the bulk form of Sr3OsO6 toward revealing the driving mechanism of 1000 K ferromagnetism, which has recently been discovered for epitaxially grown Sr 3 OsO 6 films. Unlike the film, the bulk is dominated by antiferromagnetism rather than ferromagnetism. Thus, robust ferromagnetic order appears only when Sr 3 OsO 6 is under the influence of interfaces. A specific heat capacity of 39.6(9) X 10 -3 J mol -1 K -2 is found at low temperatures (<17 K). This value is remarkably high, suggesting the presence of possible Fermionic-like excitations at the magnetic ground state. While the bulk and film forms of Sr 3 OsO 6 share the same lattice basis and electrically insulating state, the magnetism is entirely different between them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Flux Growth of Cerium Nickel Gallides Studied by In Situ Neutron Diffraction

In this work, reactions of cerium and nickel in excess molten gallium were monitored by neutron diffraction during heating and cooling. The formation of binary intermediates CeGa 2 and Ni 2 Ga 3 was observed during heating. During cooling of the molten mixture from 900 °C, precipitation of BaAl 4 -type CeNi 0.74 Ga 3.26 occurred at 850 °C. Upon cooling to 650 °C, this compound reacted in the flux to form Ce 2 NiGa 10 and then Ce 2 NiGa 12 , the latter of which persisted to room temperature. Making use of this information, subsequent reactions were quenched at 750 °C to isolate crystals of CeNi 0.74 Ga 3.26 for further study. Similar reactions replacing Ce with La and quenching above 750 °C yielded LaNi 0.35 Ga 3.65 crystals. Magnetic susceptibility studies on CeNi 0.74 Ga 3.26 indicate that the cerium is trivalent; the Ce 3+ moments undergo a strongly anisotropic ferromagnetic ordering with moment perpendicular to the c axis below 7 K. Heat capacity data show little evidence of heavy fermion behavior. Resistivity measurements show that both LaNi0.35Ga3.65 and CeNi 0.74 Ga 3.26 exhibit metallic behavior. Density of states calculations support this and indicate that Ni/Ga mixing in the compound stabilizes the structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Synthesis and Characterization of Metastable Cobalt Honeycomb KCoAsO 4

The Kitaev model has served as a long-sought-after target in the realization of a quantum spin liquid that could host Majorana Fermions. Such non-Abelian anyons could revolutionize quantum computing if properly implemented to overcome decoherence. A 3d 7 electronic configuration, like Co 2+ , has been explored by theory and experimental work to design Kitaev materials. Here, in this study, we report the synthesis of a new cobaltate honeycomb material KCoAsO 4 . The compound is synthesized through a low-temperature solution route and crystallized in space group R¯3 with lattice parameters a = 5.0394(1) and c = 28.6790(1) as determined by neutron powder diffraction. The crystal structure follows motifs similar to those of the honeycomb compound BaCo2(AsO 4 ) 2 but presents differing magnetic behavior. Magnetization/heat capacity measurements on the powder show antiferromagnetic transition T N = 14 K. Two lower-temperature transitions are present in susceptibility at low field that resemble spin reorientations. Magnetization data as a function of field have curvature indicative of metamagnetic behavior below the magnetic ordering temperature, with the magnetic ordering suppressed upon application of a higher magnetic field. Computational studies suggest the presence of a weak nearest-neighbor Kitaev term, K 1 , consistent with related honeycomb cobaltates. Together, the data suggest that this material should present a new platform for developing Kitaev quantum spin liquids.

cobaltate↗
DOE OSTI2022

QForte: An Efficient State-Vector Emulator and Quantum Algorithms Library for Molecular Electronic Structure

Here, we introduce a novel open-source software package QForte, a comprehensive development tool for new quantum simulation algorithms. QForte incorporates functionality for handling molecular Hamiltonians, fermionic encoding, ansatz construction, time evolution, and state-vector emulation, requiring only a classical electronic structure package as a dependency. QForte also contains black-box implementations of a wide variety of quantum algorithms, including variational and projective quantum eigensolvers, adaptive eigensolvers, quantum imaginary time evolution, and quantum Krylov methods. We highlight two features of QForte: (i) how the Python class structure of QForte enables the facile implementation of new algorithms, and (ii) how existing algorithms can be executed in just a few lines of code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗