Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum mechanics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Shock compression response of an insensitive high explosive single crystal: 1,1-diamino-2,2-dinitroethene (FOX-7)

Insensitive high explosives (IHEs) – that do not compromise performance – are of immense interest as a safer alternative to conventional high explosives, such as pentaerythritol tetranitrate (PETN) and cyclotrimethylenetrinitramine (RDX). Despite the strong interest in using IHEs, shock compression experiments on IHE single crystals have not been reported. To address this need, plate impact experiments were conducted to measure wave profiles in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shocked to 21 GPa longitudinal stress. Particle velocity histories, measured using laser interferometry, show a clear two-wave structure (elastic-inelastic response) at modest stresses (<3.8 GPa). Wave profiles at higher stresses show a single (overdriven) wave. Measured shock velocities and wave profiles provide accurate Hugoniot data to 21 GPa. The measured wave profiles to 21 GPa show no sign of energy release due to chemical decomposition and constitute the first demonstration of an IHE single crystal insensitivity under plane shock compression. Numerical simulations using a phenomenological material model developed for FOX-7 showed good agreement with the measured wave profiles. The experimental results and continuum simulations presented here constitute a significant first step in gaining insight into the shock compression response of IHE single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A state-of-the-art review of experimental and computational studies of granular materials: Properties, advances, challenges, and future directions

Modeling of heterogeneous materials and media is a problem of fundamental importance to a wide class of phenomena and systems, ranging from condensed matter physics, soft materials, and composite media to porous media, biological systems, geosystems, ceramic engineering, pharmaceutical science and even in space discoveries. Among the most important materials are granular systems, which have received intense interest from the engineering, physics, and mathematics communities. In this review paper, the recent developments and new advances in experimental, and computational methods on a variety of scales and physics that extend understanding to a wide range of materials and phenomena are reviewed. Experimental advances include computed neutron and nanometer-scale tomography, magnetic resonance imaging, refractive index matching, digital image correlation, acoustic emission analysis, and the most recent 4D techniques. Furthermore, a tremendous shift has occurred from the continuum scale to micro-scale and developing multiscale approaches. As such, various computational methods, including, constitutive modeling, discrete modeling, and multiscale approaches, have been developed. In conclusion, aside from all these evolutions, more complicated modeling called coupled, or multiphysics, systems representing a simultaneous presence of heat, fluid, chemical variation, and mechanical effect are also explored.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Real-time correlation of crystallization and segmental order in conjugated polymers

Modulating the segmental order in the morphology of conjugated polymers is widely recognized as a crucial factor for achieving optimal electronic properties and mechanical deformability. However, it is worth noting that the segmental order is typically associated with the crystallization process, which can result in rigid and brittle long-range ordered crystalline domains. To precisely control the morphology, a comprehensive understanding of how highly anisotropic conjugated polymers form segmentally ordered structures with ongoing crystallization is essential, yet currently elusive. To fill this knowledge gap, we developed a novel approach with a combination of stage-type fast scanning calorimetry and micro-Raman spectroscopy to capture the series of specimens with a continuum in the polymer percent crystallinity and detect the segmental order in real-time. Through the investigation of conjugated polymers with different backbones and side-chain structures, we observed a generally existing phenomenon that the degree of segmental order saturates before the maximum crystallinity is achieved. This disparity allows the conjugated polymers to achieve good charge carrier mobility while retaining good segmental dynamic mobility through the tailored treatment. Moreover, the crystallization temperature to obtain optimal segmental order can be predicted based on T g and T m of conjugated polymers. In conclusion, this in-depth characterization study provides fundamental insights into the evolution of segmental order during crystallization, which can aid in designing and controlling the optoelectronic and mechanical properties of conjugated polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The statistical mechanics of near-BPS black holes

Due to the failure of thermodynamics for low temperature near-extremal black holes, it has long been conjectured that a 'thermodynamic mass gap' exists between an extremal black hole and the lightest near-extremal state. For non-supersymmetric near-extremal black holes in Einstein gravity with an AdS2 throat, no such gap was found. Rather, at that energy scale, the spectrum exhibits a continuum of states, up to non-perturbative corrections. In this paper, we compute the partition function of near-BPS black holes in supergravity where the emergent, broken, symmetry is PSU(1, 1|2). To reliably compute this partition function, we show that the gravitational path integral can be reduced to that of a $\mathcal{N}=4$ supersymmetric extension of the Schwarzian theory, which we define and exactly quantize. In contrast to the non-supersymmetric case, we find that black holes in supergravity have a mass gap and a large extremal black hole degeneracy consistent with the Bekenstein–Hawking area. Furthermore, our results verify a plethora of string theory conjectures, concerning the scale of the mass gap and the counting of extremal micro-states.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Preseismic Fault Creep and Elastic Wave Amplitude Precursors Scale With Lab Earthquake Magnitude for the Continuum of Tectonic Failure Modes

Tectonic faults fail in a continuum of modes from slow earthquakes to elastodynamic rupture. Precursory variations in elastic wavespeed and amplitude, interpreted as indicators of imminent failure, have been observed in limited natural settings and lab experiments where they are thought to arise from contact rejuvenation and microcracking within and around the fault zone. However, the physical mechanisms and connections to fault creep are poorly understood. Here we vary loading stiffness during frictional shear to generate a range of slip modes and measure fault zone properties using transmitted elastic waves. We find that elastic wave amplitudes show clear changes before fault failure. The temporal onset of amplitude reduction scales with lab earthquake magnitude and the magnitude of this reduction varies with fault slip. Our data provide clear evidence of precursors to lab earthquakes and suggest that continuous seismic monitoring could be useful for assessing fault state and seismic hazard potential.

58 GEOSCIENCES↗
DOE OSTI2020

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Signatures of the attractive interaction in spin spectra of one-dimensional cuprate chains

Identifying the minimal model for cuprates is crucial for explaining the high- T c pairing mechanism. Recent photoemission experiments have suggested a significant near-neighbor attractive interaction V in cuprate chains, favoring pairing instability. To determine its strength, we systematically investigate the dynamical spin structure factors S ( q , ω ) using the density matrix renormalization group. Our analysis quantitatively reveals a notable softening in the two-spinon continuum, particularly evident in the intense spectrum at large momentum. This softening is primarily driven by the renormalization of the superexchange interaction, as determined by a comparison with the slave-boson theory. We also demonstrate the feasibility of detecting this spectral shift in thin-film samples using resonant inelastic x-ray scattering. Therefore, this provides a distinctive fingerprint for the attractive interaction, motivating future experiments to unveil essential ingredients in cuprates. Published by the American Physical Society 2024

Shen, Zecheng↗
DOE OSTI2021

Electromechanics of lipid-modulated gating of potassium channels

Experimental studies reveal that the anionic lipid phosphatidic acid (POPA), non-phospholipid cholesterol, and cationic lipid DOTAP inhibit the gating of voltage-sensitive potassium (Kv) channels. In this work, we develop a continuum electromechanical model to investigate the interaction of these lipids with the ion channel. Our model suggests that: (i) POPA lipids may restrict the vertical motion of the voltage-sensor domain through direct electrostatic interactions; (ii) cholesterol may oppose the radial motion of the pore domain of the channel by increasing the mechanical rigidity of the membrane; and (iii) DOTAP can reduce the effect of electrostatic forces by regulating the dielectric constant at the channel–lipid interface. The electromechanical model predictions for the three lipid types match well with the experimental observations and provide mechanistic insights into lipid-dependent gating of Kv channels.

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

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Infrared spectroscopy of nearby radio active early-type galaxies - II: spectral atlas

ABSTRACT We present a near-infrared spectroscopic atlas of nearby, bright early-type galaxies with radio emission, containing 163 galaxies observed by the Palomar 200″ TripleSpec instrument, measuring the emission line fluxes for H, He, [Fe ii], and H2 and determined the nuclear excitation mechanisms. By stacking spectra, we deduced the H2 excitation temperature (1957 ± 182 K) and dominant excitation mechanism (thermal and shock heating combined) from the K-band emission line sequence. Stacking also produces an ‘average’ spectrum of absorption features and spectral indices from the literature; the CO12 absorption line index versus J − K colour shows a trend of stronger nuclear activity producing a weaker CO12 index and a redder (flatter) continuum. The correlations between the radio and the emission-line luminosities finds a trend with radio power; however, the large scatter in the upper limits shows that the two are not directly coupled and the duty cycles of SF and AGN activity are not synchronized.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Pushing the limits of the periodic table — A review on atomic relativistic electronic structure theory and calculations for the superheavy elements

We review the progress in atomic structure theory with a focus on superheavy elements and their predicted ground state configurations important for an element’s placement in the periodic table. To understand the electronic structure and correlations in the regime of large atomic numbers, it is essential to correctly solve the Dirac equation in strong Coulomb fields, and to take into account quantum electrodynamic effects. We specifically focus on the fundamental difficulties encountered when dealing with the many-particle Dirac equation. We further discuss the possibility for future many-electron atomic structure calculations going beyond the critical nuclear charge $Z_{crit}≈170$, where levels such as the $1s$ shell dive into the negative energy continuum ($E_{nκ} < –m_{e}c^2$). The nature of the resulting Gamow states within a rigged Hilbert space formalism is highlighted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Sensitivity of Thermal Predictions to Uncertain Surface Tension Data in Laser Additive Manufacturing

To understand the process-microstructure relationships in additive manufacturing (AM), it is necessary to predict the solidification characteristics in the melt pool. This study investigates the influence of Marangoni driven fluid flow on the predicted melt pool geometry and solidification conditions using a continuum finite volume model. A calibrated laser absorptivity was determined by comparing the model predictions (neglecting fluid flow) against melt pool dimensions obtained from single laser melt experiments on a nickel super alloy 625 (IN625) plate. Using this calibrated efficiency, predicted melt pool geometries agree well with experiments across a range of process conditions. When fluid mechanics is considered, a surface tension gradient recommended for IN625 tends to overpredict the influence of convective heat transfer, but the use of an intermediate value reported from experimental measurements of a similar nickel super alloy produces excellent experimental agreement. Despite its significant effect on the melt pool geometry predictions, fluid flow was found to have a small effect on the predicted solidification conditions compared to processing conditions. We find that this result suggests that under certain circumstances, a model only considering conductive heat transfer is sufficient for approximating process-microstructure relationships in laser AM. Extending the model to multiple laser passes further showed that fluid flow also has a small effect on the solidification conditions compared to the transient variations in the process. Furthermore, limitations of the current model and areas of improvement, including uncertainties associated with the phenomenological model inputs are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Linear gyrokinetic simulations of toroidal Alfvén eigenmodes in the Mega-Amp Spherical Tokamak

Linear gyrokinetic (GK) simulations using the Gyrokinetic Toroidal Code (GTC) have been performed to investigate Toroidicity-driven Alfvén Eigenmodes (TAEs) driven by the neutral beam injection (NBI) induced fast ions in the Mega-Amp Spherical Tokamak (MAST) to identify the non-perturbative and kinetic effects of thermal plasma. A specific TAE in MAST discharge 26887, with an on-axis NBI power of approximately 1.5 MW and plasma current around 800 kA, exhibited frequency chirping, and the tangential soft x-ray camera array resolved the radial mode structure peaked near |q|=1.5. Various excitation methods were used in the GTC linear simulations, illustrating this code's capability to realistically represent the mechanisms and behaviors of fast ion-driven TAEs in spherical tokamaks. The radial structures from these GK simulations closely match measurements and calculations performed using the NOVA ideal MHD code, though with the frequencies approximately 10 kHz lower, likely due to various kinetic and non-perturbative effects. The simulations measured the damping rates due to continuum damping, radiative damping, and ion Landau damping, revealing that ion Landau damping has the most significant contribution to the total damping rate of the TAE. A comparison of growth rates of TAEs excited by fast ion Maxwellian and slowing-down distributions shows that the TAEs excited by a fast ion anisotropic pitch distribution (as part of the slowing-down distributions) are more unstable compared to those excited by a Maxwellian distribution with an equivalent fast ion beta. This shows that the use of fast ion anisotropy alters the number of fast ions to be in shear Alfvén resonance, and hence, it can greatly affect the stability of TAEs. These tests can be performed with the GTC but impossible with ideal MHD simulations, highlighting the necessity of kinetic simulations such as the GTC for a precise prediction of the TAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2025

Bottomonium meson spectrum with quenched and unquenched quark models

An open question in hadronic phenomenology concerns the “unquenching” effects of higher Fock space components on the leading Fock space description of hadrons. Here, we address this by making a comparison of the bottomonium spectrum as computed with the relativized Godfrey–Isgur quark model and an unquenched coupled channel model driven by the “$^3P_0$” mechanism of hadronic decay. Our results show that both models can describe the spectrum well, indicating that the influence of coupled channel effects can be largely absorbed into the parameters of the quenched quark model. This conclusion is reinforced by a perturbative calculation that shows that the spin-dependence of mass splittings due to mixing with the continuum recapitulates quenched quark model spin-dependent interactions. We also show that softening of the quark-antiquark wavefunction due to continuum mixing improves the description of vector bottomonium decay constants. Together, these results illustrate and substantiate the surprising robustness of simple constituent quark model descriptions of hadrons.

Sultan, M. Atif [Nankai Univ., Tianjin (China); Un↗