Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regularity”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bayesian optimization of the beam injection process into a storage ring

We have evaluated the data-efficient Bayesian optimization method for the specific task of injection tuning in a circular accelerator. In this paper, we describe the implementation of this method at the Karlsruhe Research Accelerator with up to nine tuning parameters, including the determination of the associated hyperparameters. We show that the Bayesian optimization method outperforms manual tuning and the commonly used Nelder-Mead optimization algorithm in both simulation and experiment. The algorithm was also successfully used to ease the commissioning phase after the installation of new injection magnets and is regularly used during accelerator operations. We demonstrate that the introduction of context variables that include intrabunch scattering effects, such as the Touschek effect, further improves the control and robustness of the injection process.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Approximate entropy analysis for nonlinear beam dynamics

In this paper, we applied an approximate entropy (ApEn) analysis to nonlinear beam dynamics. Due to the presence of strong nonlinear magnets, the degree of chaos in beam motion gradually increases proportionally with its amplitude. Such chaos can be quantitatively characterized with ApEn of beam turn-by-turn readings. ApEn analysis is a technique used to quantify the amount of regularity and the unpredictability of fluctuations with respect to time-series data. The ApEn, when applied as a chaos indicator, can then be used for nonlinear lattice optimization and analysis. The National Synchrotron Light Source II (NSLS-II) electron storage ring lattice was used as an example of a real-world application of our technique. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Inverse Method for Noninvasive Material Discrimination with Multienergy X-Ray Radiography

Active x-ray inspections have great potential for noninvasive detection of illicit materials. Dual-energy x-ray radiography, in particular, has been used to determine material composition by utilizing the energy-dependence of the x-ray attenuation coefficients. However, current implementations of this method are limited in their ability to determine the material composition of composite objects. Here, in this paper, we present an inverse algorithm that uses multienergy x-ray radiography data to noninvasively reconstruct the contents of a container. A critical feature of the current contribution is that material identification can be performed using conventional detectors and x-ray spectra with different endpoint energies produced by varying nothing more than the tube voltage. Adaptive regularization is used to increase the accuracy of material estimations from multienergy data sets. The utility of these methods is demonstrated with experimentally acquired radiographs obtained using a tunable x-ray source that produces spectra with endpoint energies of 100–450 keV. The object inspected is a scale model of a nuclear materials storage container composed of three-dimensional printed plastic and stainless-steel spheres inside a thin-walled steel container. Reconstructions of the steel sphere thicknesses are within a root-mean-square error of 0.37 cm.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Properties of La 0.7 Ca 0.3 MnO 3 under extreme tensile strain

The complex phase diagram of manganites with simultaneously active spin, charge, orbital, and lattice degrees of freedom continues providing surprises. In a recent groundbreaking experiment, membranes of the perovskite manganite La 0.7 Ca 0.3 MnO 3 (LCMO) deposited on a flexible polymer layer were strained up to 8% [S. S. Hong et al., Science 368, 71 (2020)], much more than achieved by regular strain induced by a rigid substrate. By increasing this strain, a metal-insulator transition was reported. Here we reproduce the results of the experiments using Monte Carlo simulations of the two-orbital double-exchange model including Jahn-Teller distortions at hole density x=1/3. The full phase diagram with varying temperature and Jahn-Teller coupling λ is presented. When the bandwidth W of mobile electrons is reduced, thus when the effective Jahn-Teller coupling λ/W is increased, a metal-insulator transition is found in our simulations, between a ferromagnetic metallic state with uniform charge distribution and an insulator with diagonal charge stripes that retains its ferromagnetic character. In between the hole-rich diagonals, staggered orbital order occurs. We also report resistivity and magnetization measurements alongside with spin correlations and charge structure factors. Our overall conclusions are in agreement with the recent experimental and density functional theory results by Hong et al., and we confirm much earlier ground state predictions of striped ferromagnetic order using energy optimization techniques by T. Hotta et al. [Phys. Rev. Lett. 86, 4922 (2001)]. The experimental observation of one of the states predicted by theory suggests that diagonal stripes could be achieved at other hole densities as well, such as x=1/4, if LCMO membranes with that hole doping were subject to similar strains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Magnetic properties of alternating Hubbard ladders

We explore the Hubbard Hamiltonian on ladders where the number of sites per rung alternates between two and three. These geometries are bipartite with nonequal or equal number of sites on the two sublattices. Thus they share a key feature of the Hubbard model in a class of lattices which Lieb has shown analytically to exhibit long-range ferrimagnetic order while being amenable to powerful numeric approaches developed for quasi-one-dimensional geometries. The density matrix renormalization group (DMRG) method is used to obtain the groundstate properties, e.g., excitation gaps, charge and spin densities as well as their correlation functions at half filling. We show the existence of long-range ferrimagnetic order in the one-dimensional ladder geometries. Our work provides detailed quantitative results which complement the general theorem of Lieb for generalized bipartite lattices. It also addresses the issue of how the alternation between quasi-long-range order and spin liquid behavior for uniform ladders with odd and even numbers of legs might be affected by a regular alternation pattern.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dislocation content of grain boundary phase junctions and its relation to grain boundary excess properties

In this work, we analyze the dislocation content of grain boundary (GB) phase junctions, i.e., line defects separating two different GB phases coexisting on the same GB plane. While regular GB disconnections have been characterized for a variety of interfaces, GB phase junctions formed by GBs with different structures and different numbers of excess atoms have not been previously studied. We apply a general Burgers circuit analysis to calculate the Burgers vectors b of junctions in two Σ5 Cu boundaries previously simulated with molecular dynamics. The Burgers vectors of these junctions cannot be described by the displacement shift complete (DSC) lattice alone. We show that, in general, the normal component of b is not equal to the difference in the GB excess volumes but contains another contribution from the numbers of GB atoms per unit area ΔN* required to transform one GB phase into another. In the boundaries studied, the latter component dominates and even changes the sign of b. We derive expressions for the normal and tangential components of b in terms of the DSC lattice vectors and the non-DSC part due to ΔN* and additional GB excess properties, including excess volume and shears. These expressions provide a connection between GB phase transformations driven by the GB free energy difference and the motion of GB junctions under applied normal and shear stresses. The proposed analysis quantifies b and therefore makes it possible to calculate the elastic part of the energy of these defects, evaluate their contribution to the nucleation barrier during GB phase transformations, and treat elastic interactions with other defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Reaching the equilibrium state of the frustrated triangular Ising magnet Ca 3 Co 2 O 6

Ca 3 Co 2 O 6 is a frustrated magnet consisting of a triangular arrangement of chains of Ising spins. It shows regular magnetization steps versus magnetic field every 1.2 T that are metastable with very slow dynamics. This has puzzled the community for many years and given rise to numerous potential theories. Here we approach the problem by seeking the elusive magnetic equilibrium state at T=2 K. To this end, we explore two approaches: (1) bypassing the slow dynamics produced by changing fields by instead field-cooling directly to the target temperature and (2) quantum annealing in transverse magnetic fields. While we observe no measurable effect of the quantum annealing in fields up to 7 T, which is likely due to the large Ising anisotropy of Co spins in this material, we find that for the field cooling in longitudinal fields we achieve the predicted equilibrium 1/3 magnetization. We perform Monte Carlo simulations of the ground state phase diagram and we also simulate the quantum annealing process and find good agreement between experiment and theory. Thus we present an investigation of the elusive ground state properties of the canonical frustrated triangular system Ca 3 Co 2 O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Drude weights in one-dimensional systems with a single defect

Ballistic transport of a quantum system can be characterized by Drude weight, which quantifies the response of the system to a uniform electric field in the infinitely long timescale. The Drude weight is often discussed in terms of the Kohn formula, which gives the Drude weight by the derivative of the energy eigenvalue of a finite-size system with the periodic boundary condition in terms of the Aharonov-Bohm flux. Recently, the Kohn formula is generalized to nonlinear responses. However, the nonlinear Drude weight determined by the Kohn formula often diverges in the thermodynamic limit. In order to elucidate the issue, in this work we examine a simple example of a one-dimensional tight-binding model in the presence of a single defect at zero temperature. We find that its linear and nonlinear Drude weights given by the Kohn formula (i) depend on the Aharonov-Bohm flux and (ii) diverge proportionally to a power of the system size. Here, we argue that the problem can be attributed to different order of limits. The Drude weight according to the Kohn formula (“Kohn-Drude weight”) indicates the response of a finite-size system to an adiabatic insertion of the Aharonov-Bohm flux. While it is a well-defined physical quantity for a finite-size system, its thermodynamic limit does not always describe the ballistic transport of the bulk. The latter should be rather characterized by a “bulk Drude weight” defined by taking the thermodynamic limit first before the zero-frequency limit. While the potential issue of the order of limits has been sometimes discussed within the linear response, the discrepancy between the two limits is amplified in nonlinear Drude weights. We demonstrate the importance of the low-energy excitations of O(1/L), which are excluded from the Kohn-Drude weight, in regularizing the bulk Drude weight.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Anharmonic phonon behavior via irreducible derivatives: Self-consistent perturbation theory and molecular dynamics

Cubic phonon interactions are now regularly computed from first principles, and the quartic interactions have begun to receive more attention. Given this realistic anharmonic vibrational Hamiltonian, the classical phonon Green's function can be precisely measured using molecular dynamics, which can then be used to rigorously assess the range of validity for self-consistent diagrammatic approaches in the classical limit. Here we use the bundled irreducible derivative approach to efficiently and precisely compute the cubic and quartic phonon interactions of CaF 2 , systematically obtaining the vibrational Hamiltonian purely in terms of irreducible derivatives. We demonstrate that the 4-phonon sunset diagram has an important contribution to the optical phonon linewidths beyond T = 500 K. Reasonable results are obtained even at T = 900 K when performing self-consistency using the 4-phonon loop diagram and evaluating the 3-phonon bubble and 4-phonon sunset diagrams post-self-consistency. Further improvements are obtained by performing quasiparticle perturbation theory, where both the 4-phonon loop and the real part of the 3-phonon bubble are employed during self-consistency. Finally, our irreducible derivative approach to self-consistent perturbation theory is a robust tool for studying anharmonic phonons in both the quantum and classical regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

(2+1)-dimensional compact Lifshitz theory, tensor gauge theory, and fractons

The (2+1)-dimensional continuum Lifshitz theory of a free compact scalar field plays a prominent role in a variety of quantum systems in condensed matter physics and high energy physics. It is known that in compact space, it has an infinite ground-state degeneracy. In order to understand this theory better, we consider two candidate lattice regularizations of it using the modified Villain formalism. We show that these two lattice theories have significantly different global symmetries (including a dipole global symmetry), anomalies, ground-state degeneracies, and dualities. In particular, one of them is self-dual. Given these theories and their global symmetries, we can couple them to corresponding gauge theories. Here, these are two different U(1) tensor gauge theories. The resulting models have excitations with restricted mobility, i.e., fractons. Finally, we give an exact lattice realization of the fracton and lineon-elasticity dualities for the Lifshitz theory, and scalar and vector charge gauge theories.

2-dimensional systems↗
DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

New, low-energy excitations in 107 Mo and 109 Mo

New ground-state level with spin-parity 1/2 + is established in 109 Mo, 69.8 keV below the previously reported 5/2(+) ground state in this nucleus, based on precise spectroscopy measurements of gamma radiation following spontaneous fission of 248 Cm performed using the Eurogam2 array of anti-Compton spectrometers. Analogous measurement of gamma radiation following spontaneous fission of 252 Cf, performed using the Gammasphere array, confirms the new 1/2 + ground state of 107 Mo, proposed recently and establishes the isomeric character of the 5/2 + first excited state in 107 Mo. Two beta-decaying isomers are suggested in 111 Mo nucleus based on regular energy systematics, supporting previous predictions. Low-energy excitations in Mo isotopes are interpreted and compared to calculations reported in the literature. The results suggest shape transition from prolate to oblate deformation at N >= 67.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Operator evolution from the similarity renormalization group and the Magnus expansion

The Magnus expansion is an efficient alternative to solving similarity renormalization group (SRG) flow equations with high-order, memory-intensive ordinary differential equation solvers. The numerical simplifications it offers for operator evolution are particularly valuable for in-medium SRG calculations, though challenges remain for difficult problems involving intruder states. Here we test the Magnus approach in an analogous but more accessible situation, which is the free-space SRG treatment of the spurious bound states arising from a leading-order chiral effective field theory (EFT) potential with very high cutoffs. In this work, we show that the Magnus expansion passes these tests and then use the investigations as a springboard to address various aspects of operator evolution that have renewed relevance in the context of the scale and scheme dependence of nuclear processes. These aspects include SRG operator flow with band- versus block-diagonal generators, universality for chiral EFT Hamiltonians and associated operators with different regularization schemes, and the impact of factorization arising from scale separation. Implications for short-range correlations physics and the possibilities for reconciling high- and low-resolution treatments of nuclear structure and reactions are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Hyperheavy spherical and toroidal nuclei: The role of shell structure

The properties of toroidal hyperheavy even-even nuclei and the role of toroidal shell structure are extensively studied within covariant density functional theory. The general trends in the evolution of toroidal shapes in the Z ≈ 130–180 region of the nuclear chart are established for the first time. These nuclei are stable with respect to breathing deformations. The most compact fat toroidal nuclei are located in the Z ≈ 136, N ≈ 206 region of the nuclear chart, but thin toroidal nuclei become dominant with increasing proton number and toward proton and neutron drip lines. The roles of toroidal shell structure, its regularity, supershell structure, and shell gaps as well as the role of different groups of the pairs of the orbitals in its formation are investigated in detail. The lowest in energy solutions at axial symmetry are characterized either by large shell gaps or low density of the single-particle states in the vicinity of the Fermi level in at least one of the subsystems (proton or neutron). Related quantum shell effects are expected to act against the instabilities in breathing and sausage deformations for these subsystems. The investigation with a large set of covariant energy density functionals reveals that substantial proton Z = 154 and 186 and neutron N = 228, 308, and 406 spherical shell gaps exist in all functionals. The nuclei in the vicinity of the combination of these particle numbers form the islands of stability of spherical hyperheavy nuclei. Furthermore, the study suggests that the N = 210 toroidal shell gap plays a substantial role in the stabilization of fat toroidal nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Two- and three-nucleon contact interactions and ground-state energies of light- and medium-mass nuclei

Classes of two-nucleon (2N) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet np scattering length and deuteron binding energy at LO, and np and pp scattering data in the laboratory-energy ranges 0–15 MeV at NLO and 0–25 MeV at N3LO. These interactions are regularized by including two Gaussian cutoffs, one for T = 0 and the other for T = 1 channels. The cutoffs are taken to vary in the ranges R 0 = (1.5 –2.3) fm and R 1 = (1.5 –3.0) fm. The 780 (1100) data points up to 15 (25) MeV energy, primarily differential cross sections, are fitted by the NLO (N3LO) models with a χ 2 / datum about 1.7 or less (well below 1.5), when harder cutoff values are adopted. As a first application, we report results for the binding energies of nuclei with mass numbers A = 3 –6 and 16 obtained with selected LO and NLO 2N models both by themselves as well as in combination with a LO three-nucleon (3N) contact interaction. The latter is characterized by a single low-energy constant that is fixed to reproduce the experimental 3 H binding energy. The inclusion of the 3N interaction largely removes the sensitivity to cutoff variations in the few-nucleon systems and leads to predictions for the 3 He and 4 He binding energies that cluster around 7.8 and 30 MeV, respectively. However, in 16 O this cutoff sensitivity remains rather strong. Finally, predictions at LO only are also reported for medium-mass nuclei with A = 40 , 48, and 90.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Comparative study of He 6 β -decay based on different similarity-renormalization-group evolved chiral interactions

We report on a study of the Gamow-Teller matrix element contributing to ${}^6$He $\beta$-decay with similarity renormalization group (SRG) versions of momentum- and configuration-space two-nucleon interactions. These interactions are derived from two different formulations of chiral effective field theory ($\chi$EFT) -- without and with the explicit inclusion of $\Delta$-isobars. We consider evolution parameters $\Lambda_{\rm SRG}$ in the range between 1.2 and 2.0 fm$^{-1}$ and, for the $\Delta$-less case, also the unevolved (bare) interaction. The axial current contains one- and two-body terms, consistently derived at tree level (no loops) in the two distinct $\chi$EFT formulations we have adopted here. The ${}^6$He and ${}^6$Li ground-state wave functions are obtained from hyperspherical-harmonics (HH) solutions of the nuclear many-body problem. In $A\,$=$\,6$ systems, the HH method is limited at present to treat only two-body interactions and non-SRG evolved currents. Our results exhibit a significant dependence on $\Lambda_{\text{SRG}}$ of the contributions associated with two-body currents, suggesting that a consistent SRG-evolution of these is needed in order to obtain reliable estimates. We also show that the contributions from one-pion-exchange currents depend strongly on the model (chiral) interactions and on the momentum- or configuration-space cutoffs used to regularize them. These results might prove helpful in clarifying the origin of the sign difference recently found in No-Core-Shell-Model and Quantum Monte Carlo calculations of the ${}^6$He Gamow-Teller matrix element.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗