Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forcing fields”

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 OSTI2025

Competing conduction mechanisms in high performance carbon nanotube fibers

The performance of carbon nanotube (CNT) cables, a contender for copper-wire replacement, is tied to its metallic and semi-conducting-like conductivity responses with temperature; the origin of the semi-conducting-like response however is an underappreciated incongruity in literature. With controlled aspect-ratio and doping-degree, over 61 unique cryogenic experiments including anisotropy and Hall measurements, CNT cable performance is explored at extreme temperatures (65 mK) and magnetic-fields (60 T). A semi-conducting-like conductivity response with temperature becomes temperature-independent approaching absolute-zero, uniquely demonstrating for the first time the necessity of heterogeneous fluctuation induced tunneling; complete de-doping leads to localized hopping, contrasting graphite's pure metallic-like response. High-field magneto-resistance (including novel +22 % longitudinal magneto-resistance near room-temperature) is analyzed with hopping and classical two-band models, both yielding a similar parameter useful for conductor development. Varying field-orientation angle uncovers significant two- and four-fold symmetries that are shown to be from Aharonov-Bohm-like corrections to the curvature-induced bandgap, a first for macroscale CNT fibers. Tight-binding calculations using Green's Function formalism model the largest, coherent transport to-date in commensurate CNT bundles in magnetic-field, revealing non-uniform transmission across bundle cross-sections with doping restoring uniformity; independent of doping, transport in bundle-junction-bundle systems are predominantly from CNTs adjacent to the other bundle—demonstrating that smaller bundles are more efficient for electronic transport. The final impact is predicting the ultimate conductivity of heterogeneous CNT cables using temperature and field-dependent transport, surpassing conductivity of traditional metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Composite pressure cell for pulsed magnets

We report that extreme pressures and high magnetic fields can affect materials in profound and fascinating ways. However, large pressures and fields are often mutually incompatible; the rapidly changing fields provided by pulsed magnets induce eddy currents in the metallic components used in conventional pressure cells, causing serious heating, forces, and vibration. Here, we report a diamond-anvil-cell made mainly out of insulating composites that minimizes inductive heating while retaining sufficient strength to apply pressures of up to 8 GPa. Any residual metallic component is made of low-conductivity metals and patterned to reduce eddy currents. The simple design enables rapid sample or pressure changes, desired by pulsed-magnetic-field-facility users. The pressure cell has been used in pulsed magnetic fields of up to 65 T with no noticeable heating at cryogenic temperatures. Several measurement techniques are possible inside the cell at temperatures as low as 500 mK.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Combining Eddy Covariance Towers, Field Measurements, and the MEMS 2 Ecosystem Model Improves Confidence in the Climate Impacts of Bioenergy With Carbon Capture and Storage

ABSTRACT Carbon dioxide removal technologies such as bioenergy with carbon capture and storage (BECCS) are required if the effects of climate change are to be reversed over the next century. However, BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty. We place these net ecosystem carbon balance forcings in the context of associated radiative forcings from changes to soil N 2 O and CH 4 fluxes, land surface albedo, embedded fossil fuel use, and geologically stored carbon. Results from long term eddy covariance measurements, a soil and plant carbon inventory, and the MEMS 2 process‐based ecosystem model all agree that establishing perennials such as switchgrass or mixed prairie on former cropland resulted in net negative radiative forcing (i.e., global cooling) of −26.5 to −39.6 fW m −2 over 100 years. Establishing these perennials on former grassland sites had similar climate mitigation impacts of −19.3 to −42.5 fW m −2 . However, the largest climate mitigation came from establishing corn for BECCS on former cropland or grassland, with radiative forcings from −38.4 to −50.5 fW m −2 , due to its higher plant productivity and therefore more geologically stored carbon. Our results highlight the strengths and limitations of each method for quantifying the field scale climate impacts of BECCS and show that utilizing multiple methods can increase confidence in the final radiative forcing estimates.

Falvo, Grant [Department of Plant, Soil and Microb↗
DOE OSTI2024

Gravitational scattering and beyond from extreme mass ratio effective field theory

We explore a recently proposed effective field theory describing electromagnetically or gravitationally interacting massive particles in an expansion about their mass ratio, also known as the self-force (SF) expansion. By integrating out the deviation of the heavy particle about its inertial trajectory, we obtain an effective action whose only degrees of freedom are the lighter particle together with the photon or graviton, all propagating in a Coulomb or Schwarzschild background. The 0SF dynamics are described by the usual background field method, which at 1SF is supplemented by a “recoil operator” that encodes the wobble of the heavy particle, and similarly computable corrections appearing at 2SF and higher. Our formalism exploits the fact that the analytic expressions for classical backgrounds and particle trajectories encode dynamical information to all orders in the couplings, and from them we extract multiloop integrands for perturbative scattering. As a check, we study the two-loop classical scattering of scalar particles in electromagnetism and gravity, verifying known results. We then present new calculations for the two-loop classical scattering of dyons, and of particles interacting with an additional scalar or vector field coupling directly to the lighter particle but only gravitationally to the heavier particle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN

This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.

42 ENGINEERING↗
DOE OSTI2026

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗
DOE OSTI2023

Engineering of Self‐Aggregation‐Resistant MnO 2 Heterostructure with A Built‐in Field for Enhanced High‐Mass‐Loading Energy Storage

Abstract Although MnO 2 has been intensively investigated for energy storage, further applications are limited by van der Waals force‐triggered self‐aggregation that always leads to poorly exposed active sites and compromised reaction dynamics, especially under high‐mass‐loading conditions. Herein, by synergistically coupling interfacial modulation with the Kirkendall effect, this work achieves in situ topological structure reorganization of MnOOH toward the high‐aspect‐ratio MnO 2 heterostructure (Heter‐MnO 2 ) with fully exposed active sites, which is ready to assemble into self‐supporting high‐mass‐loading film (30 mg cm −2 ) with restrained self‐aggregation. Theoretical calculation and dynamics analysis results demonstrate the generation of the built‐in field within the heterostructure, thus enhancing the electronic‐transfer and ionic‐adsorption/transport rates. As such, the 30 mg cm −2 Heter‐MnO 2 electrode achieves a superior areal capacitance of 4762 mF cm −2 at 1 mA cm −2 and a sound rate performance (79% at 100 mA cm −2 ) comparable to those of low‐mass‐loading/thin‐film electrodes. As a proof of concept, the fabricated planar interdigital quasi‐solid‐state symmetric micro‐supercapacitor (MSC) based on the Heter‐MnO 2 electrode can deliver a remarkable areal capacitance of 181 mF cm −2 and a considerable volumetric energy density of 10.3 mWh cm −3 . This methodology highlights the promise of surface/interface chemistry modulation for the configuration of easy‐to‐integrate hierarchical nanostructures to better meet practical energy applications.

Chemistry↗
DOE OSTI2021

Integrated atomic force microscopy and x-ray irradiation for in situ characterization of radiation-induced processes

Understanding radiation-induced chemical and physical transformations at material interfaces is important across diverse fields, but experimental approaches are often limited to either ex situ observations or in situ electron microscopy or synchrotron-based methods in which cases the radiation type and dose are inextricably tied to the imaging basis itself. In this work we overcome this limitation by demonstrating integration of an X-ray source with an atomic force microscope to directly monitor radiolytically-driven interfacial chemistry at the nanoscale. We illustrate the value of in situ observations by examining effects of radiolysis on material adhesion forces in aqueous solution, as well as examining the production of alkali nitrates at the interface between an alkali halide crystal surface and air. For the examined salt-air interface, direct visualization under flexible experimental conditions greatly extends prior observations by enabling the transformation process to be followed comprehensively from source-to-sink with mass balance quantitation. Overall, our novel rad-AFM opens doors into understanding the dynamics of radiolytically-driven mass transfer and surface alteration at the nanoscale in real-time.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Near-cancellation of up- and down-gradient momentum transport in forced magnetized shear-flow turbulence

Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Furthermore, even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to 98%. These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the falloff rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Entanglement-enhanced optomechanical sensor array with application to dark matter searches

Abstract Squeezed light has long been used to enhance the precision of a single optomechanical sensor. An emerging set of proposals seeks to use arrays of optomechanical sensors to detect weak distributed forces, for applications ranging from gravity-based subterranean imaging to dark matter searches; however, a detailed investigation into the quantum-enhancement of this approach remains outstanding. Here, we propose an array of entanglement-enhanced optomechanical sensors to improve the broadband sensitivity of distributed force sensing. By coherently operating the optomechanical sensor array and distributing squeezing to entangle the optical fields, the array of sensors has a scaling advantage over independent sensors (i.e., $$\sqrt{M}\to M$$ M → M , where M is the number of sensors) due to coherence as well as joint noise suppression due to multi-partite entanglement. As an illustration, we consider entanglement-enhancement of an optomechanical accelerometer array to search for dark matter, and elucidate the challenge of realizing a quantum advantage in this context.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Gravitational vortex disruption in mesofluidic separators

Typically, bump array devices separate submicrometer particles where gravitational effects remain small. However, as particle size increases from a few micrometers to hundreds of micrometers, gravitational effects become important and settling may induce vortex formation. Here, this article shows under laminar flow that these vortices suffice to disrupt flow fields, frustrating particle separation. The limiting diameter, above which vortices do form, is determined, and the dependence of the limiting diameter on viscosity, particle properties, and separator orientation is quantified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗
DOE OSTI2025

Magnetic and Mechanical Analysis of Bi-2212 Rutherford Cable in a Cos-Theta Sub-Scale Dipole Coil

The U.S. Magnet Development Program (US-MDP) explores high-field accelerator magnets compatible with operational conditions beyond the limits of Nb$_3$Sn technology. The ongoing R&D High-Temperature Superconductors (HTS) suggests using Bi$_2$Sr$_2$CaCu$_2$O$_{8-x}$ (Bi-2212) as superconducting element. Bi-2212 Rutherford cables maintain a high critical current (I$_C$) when exposed to a large external magnetic field. However, Bi-2212 exhibits an oversensitive stress-strain response when subject to large Lorentz forces. This paper reports on the magnetic and mechanical analysis of the Bi-2212 cosine-theta insert being developed at Fermilab for a hybrid magnet composed of two external layers of Nb$_3$Sn and two internal layers of Bi-2212. We performed a FEM analysis of the insert to estimate the HTS stress state in the coil's strands under magnetic and mechanical loads.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Light mesons within the basis light-front quantization framework

In this work, we study the light-unflavored mesons as relativistic bound states in the nonperturbative Hamiltonian formalism of the basis light-front quantization approach. The dynamics for the valence quarks of these mesons is specified by an effective Hamiltonian containing the one-gluon exchange interaction and the confining potentials both introduced in our previous work on heavy quarkonia, supplemented additionally by a pseudoscalar contact interaction. We diagonalize this Hamiltonian in our basis function representation to obtain the mass spectrum and the light-front wave functions (LFWFs). Based on these LFWFs, we then study the structure of these mesons by computing the electromagnetic form factors, the decay constants, the parton distribution amplitudes, and the parton distribution functions. Our results are comparable to those from experiments and other theoretical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Bridging hadronic and vacuum structure by heavy quarkonia

We discuss the central and, mostly, spin-dependent potentials in heavy quarkonia $\overline{b}$$b$,$\overline{c}$$c$, with two goals in mind. The first is phenomenological: using the splitting between the 1S and 2S pairs, as well as the 1P and 2P quartet masses, we obtain pertinent matrix elements of all these potentials. The second is theoretical: using standard wave functions, we compute these matrix elements, quantifying the perturbative, nonperturbative and “string” contributions. The model for nonperturbative effects is a “dense instanton liquid model,” in which the QCD vacuum is made mostly of (strongly color correlated) instanton-anti-instanton pairs or “molecules.” Their contributions are calculated via standard Wilson lines, with or without extra powers of gauge fields. We find a reasonable description of all central, spin-spin, and spin-orbit forces at distances 𝑟 = 0–0.7 fm, relevant for $\overline{b}$$b$ and $\overline{c}$$c$ quarkonia.

Bound states↗