Engineering topics

Wirth, R.

Publications and source records attributed to Wirth, R..

7 records

DOE OSTI2022

Ab initio studies of the double–Gamow-Teller transition and its correlation with neutrinoless double- β decay

Here, we use chiral interactions and several ab initio methods to compute the nuclear matrix elements (NMEs) for ground-state-to-ground-state double Gamow-Teller transitions in a range of isotopes and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the ab initio calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double beta decay.

20 ≤ A ≤ 38↗

DOE OSTI2021

Ab Initio Calculation of the Contact Operator Contribution in the Standard Mechanism for Neutrinoless Double Beta Decay

Starting from chiral nuclear interactions, we evaluate the contribution of the leading-order contact transition operator to the nuclear matrix element (NME) of neutrinoless double-beta decay, assuming a light Majorana neutrino-exchange mechanism. The corresponding low-energy constant (LEC) is determined by fitting the transition amplitude of the nn → ppe – e – process to a recently proposed synthetic datum. We examine the dependence of the amplitude on similarity renormalization group (SRG) scale and chiral expansion order of the nuclear interaction, finding that both dependences can be compensated to a large extent by readjusting the LEC. We evaluate the contribution of both the leading-order contact operator and standard long-range operator to the neutrinoless double-beta decays in the light nuclei 6,8 He and the candidate nucleus 48 Ca. Furthermore, our results provide the first clear demonstration that the contact term enhances the NME by 43(7)% in 48 Ca, where the uncertainty is propagated from the synthetic datum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Singular value decomposition and similarity renormalization group evolution of nuclear interactions

One of the main challenges for ab initio nuclear many-body theory is the growth of computational and storage costs as calculations are extended to heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform Singular Value Decompositions of nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization scheme and resolution scales. We develop and implement the Similarity Renormalization Group (SRG) evolution of the factorized interaction, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves, and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Lastly, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the In-Medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Symmetry reduction of tensor networks in many-body theory: I. Automated symbolic evaluation of SU(2) algebra

Abstract The ongoing progress in (nuclear) many-body theory is accompanied by an ever-rising increase in complexity of the underlying formalisms used to solve the stationary Schrödinger equation. The associated working equations at play in state-of-the-art ab initio nuclear many-body methods can be analytically reduced with respect to angular-momentum, i.e. SU (2), quantum numbers whenever they are effectively employed in a symmetry-restricted context. The corresponding procedure constitutes a tedious and error-prone but yet an integral part of the implementation of those many-body frameworks. Indeed, this symmetry reduction is a key step to advance modern simulations to higher accuracy since the use of symmetry-adapted tensors can decrease the computational complexity by orders of magnitude. While attempts have been made in the past to automate the (anti-) commutation rules linked to Fermionic and Bosonic algebras at play in the derivation of the working equations, there is no systematic account to achieve the same goal for their symmetry reduction. In this work, the first version of an automated tool performing graph-theory-based angular-momentum reduction is presented. Taking the symmetry-unrestricted expressions of a generic tensor network as an input, the code provides their angular-momentum-reduced form in an error-safe way in a matter of seconds. Several state-of-the-art many-body methods serve as examples to demonstrate the generality of the approach and to highlight the potential impact on the many-body community.

Tichai, A.↗

DOE OSTI2020

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48 Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. In this work, we use our scheme to compute the matrix element that governs the neutrinoless double beta decay of 48 Ca to 48 Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76 Ge, 130 Te, and 136 Xe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

The transformation of diamond to graphite: Experiments reveal the presence of an intermediate linear carbon phase

Natural diamonds that have been partially replaced by graphite have been observed to occur in natural rocks. While the graphite-to-diamond phase transition has been extensively studied the opposite of this (diamond to graphite) remains poorly understood. We performed high-pressure and temperature hydrous and anhydrous experiments up to 1.0 GPa and 1300 °C using Amplex premium virgin synthetic diamonds (20–40 μm size) as the starting material mixed with Mg(OH) 2 as a source of H 2 O for the hydrous experiments. The experiments revealed that the diamond-to-graphite transformation at P = 1.0 GPa and T = 1300 °C was triggered by the presence of H 2 O and was accomplished through a three-stage process. Stage 1: diamond reacts with a supercritical H 2 O producing an intermediate 200–500 nm size “globular carbon” phase. This phase is a linear carbon chain; i.e. a polyyne or carbyne. Stage 2: the linear carbon chains are unstable and highly reactive, and they decompose by zigzagging and cross-linking to form sp 2 -hybridized structures. Stage 3: normal, disordered, and onion-like graphite is produced by the decomposition of the sp-hybridized carbon chains which are re-organized into sp 2 bonds. Our experiments show that there is no direct transformation from sp 3 C-bonds into sp 2 C-bonds. Our hydrous high-pressure and high-temperature experiments show that the diamond-to-graphite transformation requires an intermediate metastable phase of a linear hydrocarbon. This process also provides a simple mechanism for the substitution of other elements into the graphite structure (e.g. H, S, O).

36 MATERIALS SCIENCE↗