Search NASA⌕ Search

SEARCH · Search NASA

Results for “FINITE-RANGE INTERACTIONS”

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 37 records · Page 2

Quantum spin liquids bootstrapped from Ising criticality in Rydberg arrays

Arrays of Rydberg atoms constitute a highly tunable, strongly interacting venue for the pursuit of exotic states of matter. We develop a strategy for accessing a family of fractionalized phases known as quantum spin liquids in two-dimensional Rydberg arrays. We specifically use effective field theory methods to study arrays assembled from Rydberg chains tuned to an Ising phase transition that famously hosts emergent fermions propagating within each chain. This highly entangled starting point allows us to naturally access spin liquids familiar from Kitaev's honeycomb model—albeit from an entirely different framework. In particular, we argue that finite-range repulsive Rydberg interactions, which frustrate nearby symmetry-breaking orders, can enable coherent propagation of emergent fermions between the chains in which they were born. Delocalization of emergent fermions across the full two-dimensional Rydberg array yields a gapless Z 2 spin liquid with a single massless Dirac cone. Here, the Rydberg occupation numbers exhibit universal power-law correlations that provide a straightforward experimental diagnostic of this phase. We further show that explicitly breaking symmetries perturbs the gapless spin liquid into gapped, topologically ordered descendants: Breaking lattice symmetries generates toric-code topological order, whereas introducing Floquet-mediated chirality generates non-Abelian Ising topological order. In the toric-code phase, we analytically construct microscopic incarnations of non-Abelian defects, which can be created and transported by dynamically controlling the atom positions in the array. Our work suggests that appropriately tuned Rydberg arrays provide a cold-atoms counterpart of solid-state “Kitaev materials” and, more generally, it spotlights a different angle for pursuing experimental platforms for Abelian and non-Abelian fractionalization.

36 MATERIALS SCIENCE↗

Solving scattering problems in the half-line using methods developed for scattering in the full line

We reduce the solution of the scattering problem defined on the half-line [0,∞) by a real or complex potential v(x) and a general homogeneous boundary condition at x=0 to that of the extension of v(x) to the full line that vanishes for x<0. We find an explicit expression for the reflection amplitude of the former problem in terms of the reflection and transmission amplitudes of the latter. We obtain a set of conditions on these amplitudes under which the potential in the half-line develops bound states, spectral singularities, and time-reversed spectral singularities where the potential acts as a perfect absorber. We examine the application of these results in the study of the scattering properties of a δ-function potential and a finite barrier potential defined in [0,∞), discuss optical systems modeled by these potentials, and explore the configurations in which these systems act as a laser or perfect absorber. In particular, we derive an explicit formula for the laser threshold condition for a slab laser with a single mirror and establish the surprising fact that a nearly perfect mirror gives rise to a lower threshold gain than a perfect mirror. We also offer a nonlinear extension of our approach which allows for utilizing a recently developed nonlinear transfer matrix method in the full line to deal with finite-range nonlinear scattering interactions defined in the half-line.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Circuit Quantum Electrodynamics in Hyperbolic Space: From Photon Bound States to Frustrated Spin Models

Circuit quantum electrodynamics is one of the most promising platforms for efficient quantum simulation and computation. In recent groundbreaking experiments, the immense flexibility of superconducting microwave resonators was utilized to realize hyperbolic lattices that emulate quantum physics in negatively curved space. Here we investigate experimentally feasible settings in which a few superconducting qubits are coupled to a bath of photons evolving on the hyperbolic lattice. We compare our numerical results for finite lattices with analytical results for continuous hyperbolic space on the Poincaré disk. We find good agreement between the two descriptions in the long-wavelength regime. We show that photon-qubit bound states have a curvature-limited size. We propose to use a qubit as a local probe of the hyperbolic bath, for example, by measuring the relaxation dynamics of the qubit. We find that, although the boundary effects strongly impact the photonic density of states, the spectral density is well described by the continuum theory. We show that interactions between qubits are mediated by photons propagating along geodesics. We demonstrate that the photonic bath can give rise to geometrically frustrated hyperbolic quantum spin models with finite-range or exponentially decaying interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantified Gamow shell model interaction for p s d -shell nuclei

Background: The structure of weakly bound and unbound nuclei close to particle drip lines is one of the major science drivers of nuclear physics. A comprehensive understanding of these systems goes beyond the traditional configuration interaction approach formulated in the Hilbert space of localized states (nuclear shell model) and requires an open quantum system description. The complex-energy Gamow shell model (GSM) provides such a framework as it is capable of describing resonant and nonresonant many-body states on equal footing. Purpose: To make reliable predictions, quality input is needed that allows for the full uncertainty quantification of theoretical results. In this study, we carry out the optimization of an effective GSM (one-body and two-body) interaction in the psdf-shell-model space. The resulting interaction is expected to describe nuclei with 5 ≤ A ≲ 12 at the p-sd-shell interface. Method: The one-body potential of the 4 He core is modeled by a Woods-Saxon + spin-orbit + Coulomb potential, and the finite-range nucleon-nucleon interaction between the valence nucleons consists of central, spin-orbit, tensor, and Coulomb terms. The GSM is used to compute key fit observables. The χ 2 optimization is performed using the Gauss-Newton algorithm augmented by the singular value decomposition technique. The resulting covariance matrix enables quantification of statistical errors within the linear regression approach. Results: The optimized one-body potential reproduces nucleon- 4 He scattering phase shifts up to an excitation energy of 20 MeV. The two-body interaction built on top of the optimized one-body field is adjusted to the bound and unbound ground-state binding energies and selected excited states of the helium, lithium, and beryllium isotopes up to A = 9 . A very good agreement with experimental results was obtained for binding energies. First applications of the optimized interaction include predictions for two-nucleon correlation densities and excitation spectra of light nuclei with quantified uncertainties. In conclusion: The new interaction will enable comprehensive and fully quantified studies of structure and reactions aspects of nuclei from the psd region of the nuclear chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Search for excited states in O 25

Theoretical calculations suggest the presence of low-lying excited states in 25 O. Previous experimental searches by means of proton knockout on 26 F produced no evidence for such excitations. We search for excited states in 25 O using the 24 O(d,p) 25 O reaction. The theoretical analysis of excited states in unbound 25,27 O is based on the configuration interaction approach that accounts for couplings to the scattering continuum. We use invariant-mass spectroscopy to measure neutron-unbound states in 25 O. For the theoretical approach, we use the complex-energy Gamow Shell Model and Density Matrix Renormalization Group method with a finite-range two-body interaction optimized to the bound states and resonances of 23–26 O, assuming a core of 22 O. Here we predict energies, decay widths, and asymptotic normalization coefficients. Our calculations in a large spdf space predict several low-lying excited states in 25 O of positive and negative parity, and we obtain an experimental limit on the relative cross section of a possible J π = 1/2 + state with respect to the ground state of 25 O at σ 1/2+ /σ g.s. = 0.25$^{+1.0}_{-0.25}$. We also discuss how the observation of negative parity states in 25 O could guide the search for the low-lying negative parity states in 27 O. Previous experiments based on the proton knockout of 26 F suffered from the low cross sections for the population of excited states in 25 O because of low spectroscopic factors. In this respect, neutron transfer reactions carry more promise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Saturating the one-axis twisting quantum Cramér-Rao bound with a total spin readout

We show that the lowest quantum Cramér-Rao bound achievable in interferometry with a one-axis twisted spin coherent state is saturated by the asymptotic method of moments error of a protocol that uses one call to the one-axis twisting, one call to time-reversed one-axis twisting, and a final total spin measurement (i.e., a twist-untwist protocol). The result is derived by first showing that the metrological phase diagram for one-axis twisting is asymptotically characterized by a single quantum Fisher information value N(N + 1)/2 for all times, then constructing a twist-untwist protocol having a method of moments error that saturates this value. The case of finite-range one-axis twisting is similarly analyzed, and a simple functional form for the metrological phase diagram is found in both the short-range and long-range interaction regimes. Numerical evidence suggests that the finite-range analogues of twist-untwist protocols can exhibit a method of moments error that asymptotically saturates the lowest quantum Cramér-Rao bound achievable in interferometry with finite-range one-axis twisted spin coherent states for all interaction times.

36 MATERIALS SCIENCE↗

Operator product expansion beyond leading order for two-component fermions

Here, we consider a homogeneous, balanced gas of strongly interacting fermions in two spin states interacting through a large scattering length. Finite-range corrections are needed for a quantitative description of data which experiments and numerical simulations have provided. We use a perturbative field-theoretical framework and a tool called the operator product expansion (OPE), which together allow for the expression of finite-range corrections to the universal relations and momentum distribution. By using the OPE, we derive the 1/k 6 part of the momentum tail, which is related to the sum of the derivative of the energy with respect to the finite range and the averaged kinetic energy of opposite spin pairs. By comparing the 1/k 4 term and the 1/k 6 correction in the momentum distribution to provided quantum Monte Carlo (QMC) data, we show that including the 1/k 6 part offers marked improvements. Our field-theoretical approach enables for a clear understanding of the role of the scattering length and finite effective range in the universal relations and the momentum distribution.

74 ATOMIC AND MOLECULAR PHYSICS↗

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Axial deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator Basis

HFBTHO is a physics computer code that is used to model the structure of the nucleus. It is an implementation of the nuclear energy Density Functional Theory (DFT), where the energy of the nucleus is obtained by integration over space of some phenomenological energy density, which is itself a functional of the neutron and proton densities. In HFBTHO, the energy density derives either from the zero-range Dkyrme or the finite-range Gogny effective two-body interaction between nucleons. Nuclear superfluidity is treated at the Hartree-Fock-Bogoliubov (HFB) approximation, and axial-symmetry of the nuclear shape is assumed. This version is the 3rd release of the program; the two previous versions were published in Computer Physics Communications [1,2]. The previous version was released at LLNL under GPL 3 Open Source License and was given release code LLNL-CODE-573953.

Perez, R. Navarro↗

Numerical analysis for finite-range multitype stochastic contact financial market dynamic systems

In an attempt to reproduce and study the dynamics of financial markets, a random agent-based financial price model is developed and investigated by the finite-range multitype contact dynamic system, in which the interaction and dispersal of different types of investment attitudes in a stock market are imitated by viruses spreading. With different parameters of birth rates and finite-range, the normalized return series are simulated by Monte Carlo simulation method and numerical studied by power-law distribution analysis and autocorrelation analysis. To better understand the nonlinear dynamics of the return series, a q-order autocorrelation function and a multi-autocorrelation function are also defined in this work. The comparisons of statistical behaviors of return series from the agent-based model and the daily historical market returns of Shanghai Composite Index and Shenzhen Component Index indicate that the proposed model is a reasonable qualitative explanation for the price formation process of stock market systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron-Phonon Systems on a Universal Quantum Computer

We present an algorithm that extends existing quantum algorithms for simulating fermion systems in quantum chemistry and condensed matter physics to include bosons in general and phonons in particular. We introduce a qubit representation for the low-energy subspace of phonons which allows an efficient simulation of the evolution operator of the electron-phonon systems. As a consequence of the Nyquist-Shannon sampling theorem, the phonons are represented with exponential accuracy on a discretized Hilbert space with a size that increases linearly with the cutoff of the maximum phonon number. The additional number of qubits required by the presence of phonons scales linearly with the size of the system. The additional circuit depth is constant for systems with finite-range electron-phonon and phonon-phonon interactions and linear for long-range electron-phonon interactions. Our algorithm for a Holstein polaron problem was implemented on an Atos quantum learning machine quantum simulator employing the quantum phase estimation method. The energy and the phonon number distribution of the polaron state agree with exact diagonalization results for weak, intermediate, and strong electron-phonon coupling regimes.

97 MATHEMATICS AND COMPUTING↗

One-dimensional Hubbard-Holstein model with finite-range electron-phonon coupling

The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling g to the phonons drives the system towards a Peierls charge density wave state, whereas increasing the electron-electron interaction U drives the fermions into a Mott antiferromagnet. At low g and U, or when doped, the system is metallic. In one dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long-range coupling to phonons, with characteristic range ξ, interpolating between the Holstein and Fröhlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range ξ and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction U will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. Furthermore, the phase separation region is once again present in the U ≠ 0 case, even for small values of the coupling range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of choosing the Q{sub α}-values and daughter density distributions on the magic numbers predicted by α decays

We investigated the effects of choosing the α-decay energy (Q{sub α}) and particular density distributions of the daughter nuclei on predicting unknown magic numbers, or confirming some of them. The Q{sub α} values were calculated using three models, namely the finite-range droplet model (FRDM), the Weizsäcker–Skyrme (WS4) mass model, and a recent fitting formula. The α-decay half-lives (T{sub α}) of superheavy nuclei of Z = 118–124 have been calculated within the generalized density-dependent cluster model. The double-folding α+core potential was determined based on the M3Y-Reid nucleon–nucleon interaction. We have used the two-parameter Fermi distribution to represent the nuclear density, with four sets of parameters that were determined based on electron scattering data and different microscopic calculations. Upon investigating the relative stability of the mentioned superheavy isotopic chains, we found that the Q{sub α}-values deduced from FRDM and WS4 models predict different neutron magic numbers. The examined fitting formula failed to predict any magic number larger than N = 162. Under normalization to the proton and neutron (N) numbers, the different choices of the density distributions showed almost no effect on the general behavior of T{sub α} with N, but changed its values up to one order of magnitude. Even though the maxima and minima of T{sub α} clearly follow the characteristic minima and maxima of Q{sub α}, respectively, the results showed that precise density distribution and deformations, in addition to a correct preformation probability, are necessarily required for estimating an accurate value of the half-life.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Beta-Delayed Neutron Emission of 𝑁=84 132Cd

Using the time-of-flight technique, we measured the beta-delayed neutron emission of 132 Cd . From our large-scale shell model (LSSM) calculation using the N3⁢LO interaction [, Phys. Rev. Lett. 131, 022501 (2023)], we suggest the decay is dominated by the transformation of a neutron in the 𝑔7/2 orbital, deep below the Fermi surface, into a proton in the 𝑔9/2 orbital. We compare the beta-decay half-lives and neutron branching ratios of nuclei with 𝑍<50 and 𝑁 ≥82 obtained with our LSSM with those of leading “global” models such as finite-range droplet model (FRDM). Our calculations match known half-lives and neutron branching ratios well and suggest that current leading models overestimate the yet-to-be-measured half-lives. Our model, backed by the 132 Cd decay data presented here, offers robust predictive power for nuclei of astrophysical interest such as 𝑟-process waiting points.

Madurga, M [Department of Physics and Astronomy, U↗

Finite-range pairing in nuclear density functional theory

Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. Here, we show that a folding radius of about 1 fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.

Nuclear density functional theory↗

Nuclear properties for astrophysical and radioactive-ion-beam applications (II)

We tabulate the ground-state odd-proton and odd-neutron spins and parities, proton and neutron pairing gaps, one- and two-neutron separation energies, quantities related to β-delayed one- and two-neutron emission probabilities, average energy and average number of emitted neutrons, β-decay energy release and half-life with respect to Gamow–Teller decay with a phenomenological treatment of first-forbidden decays, one- and two-proton separation energies, and α-decay energy release and half-life for 9318 nuclei ranging from {sup 16}O to {sup 339}136 and extending from the proton drip line to the neutron drip line. This paper is a new and improved version of Atomic Data And Nuclear Data Tables [66 131 (1997)]. The starting point of our present work is the new study (FRDM(2012)) of nuclear ground-state masses and deformations based on the finite-range droplet model and folded-Yukawa single-particle potential published in a previous issue of Atomic Data And Nuclear Data Tables [109–110, 1 (2016)]. The β-delayed neutron-emission probabilities and Gamow–Teller β-decay rates are obtained from a quasi-particle random-phase approximation with single-particle levels and wave functions at the calculated nuclear ground-state shapes as input quantities. A development since 1997 is we now use a Hauser–Feshbach approach to account for (n, γ) competition and treat first-forbidden decay in a phenomenological approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Renormalization of a finite-range inverse-cube potential

We study the regularization and renormalization of a finite-range inverse-cube potential in the two- and three-body sectors. Specifically, we compare and contrast three different regulation schemes frequently used to study few-body systems as well as the associated renormalization group (RG) flows. We also calculate bound-state and scattering observables over a wide range of cutoffs, demonstrating the sufficiency of a two-body contact interaction to renormalize two- and three-body observables. Here, we supplement these plots with quantified analyses of the observables' residual cutoff dependences.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gogny-force-derived effective shell-model Hamiltonian

The density-dependent finite-range Gogny force has been used to derive the effective Hamiltonian for the shell-model calculations of nuclei. The density dependence simulates an equivalent three-body force, while the finite range gives a Gaussian distribution of the interaction in the momentum space and hence leads to an automatic smooth decoupling between low-momentum and high-momentum components of the interaction, which is important for finite-space shell-model calculations. Two-body interaction matrix elements, single-particle energies, and the core energy of the shell model can be determined by the unified Gogny force. The analytical form of the Gogny force is advantageous to treat cross-shell cases, while it is difficult to determine the cross-shell matrix elements and single-particle energies using an empirical Hamiltonian by fitting experimental data with a large number of matrix elements. In this paper, we have applied the Gogny-force effective shell-model Hamiltonian to the p- and sd-shell nuclei. The results show good agreements with experimental data and other calculations using empirical Hamiltonians. The experimentally known neutron drip line of oxygen isotopes and the ground states of typical nuclei 10B and 18N can be reproduced, in which the role of three-body force is non-negligible. Here, the Gogny-force-derived effective Hamiltonian has also been applied to the cross-shell calculations of the sd–pf shell.

20 ≤ A ≤ 38↗