Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear structure & decay”

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

Model-Independent Study of Structure in B + → D + D - K + Decays

The only anticipated resonant contributions to B + → D + D - K + decays are charmonium states in the D + D - channel. A model-independent analysis, using LHCb proton-proton collision data taken at center of-mass energies of $\sqrt{s}$ = 7, 8, and 13 TeV, corresponding to a total integrated luminosity of 9 fb -1 , is carried out to test this hypothesis. The description of the data assuming that resonances only manifest in decays to the D + D - pair is shown to be incomplete. This constitutes evidence for a new contribution to the decay, potentially one or more new charm-strange resonances in the D + K - channel with masses around 2.9 GeV/c 2

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Scalable Circuits for Preparing Ground States on Digital Quantum Computers: The Schwinger Model Vacuum on 100 Qubits

The vacuum of the lattice Schwinger model is prepared on up to 100 qubits of IBM’s Eagle-processor quantum computers. A new algorithm to prepare the ground state of a gapped translationally invariant system on a quantum computer is presented, which we call “scalable circuits ADAPT-VQE” (SC-ADAPT-VQE). This algorithm uses the exponential decay of correlations between distant regions of the ground state, together with ADAPT-VQE, to construct quantum circuits for state preparation that can be scaled to arbitrarily large systems. These scalable circuits can be determined with use of classical computers, avoiding the challenging task of optimizing parameterized circuits on a quantum computer. SC-ADAPT-VQE is applied to the Schwinger model, and is shown to be systematically improvable, with an accuracy that converges exponentially with circuit depth. Both the structure of the circuits and the deviations of prepared wave functions are found to become independent of the number of spatial sites, L . This allows a controlled extrapolation of the circuits, determined with use of small or modest-sized systems, to arbitrarily large L . The circuits for the Schwinger model are determined on lattices up to L = 14 (28 qubits) with the Qiskit classical simulator, and are subsequently scaled up to prepare the L = 50 (100 qubits) vacuum on IBM’s 127-superconducting-qubit quantum computers ibm_brisbane and ibm_cusco. After introduction of an improved error-mitigation technique, which we call “operator decoherence renormalization”, the chiral condensate and charge-charge correlators obtained from the quantum computers are found to be in good agreement with classical matrix product state simulations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Search for heavy pseudoscalar and scalar bosons decaying to a top quark pair in proton–proton collisions at $\sqrt{s}$ = 13 TeV

A search for pseudoscalar or scalar bosons decaying to a top quark pair ($\textrm{t}$$\bar{\textrm{t}}$) in final states with one or two charged leptons is presented. The analyzed proton–proton collision data was recorded at $\sqrt{s}$ = 13 TeV by the CMS experiment at the CERN LHC and corresponds to an integrated luminosity of 138 fb -1 . The invariant mass $m_{{\textrm{t}}{\bar{\textrm{t}}}}$ of the reconstructed $\textrm{t}$$\bar{\textrm{t}}$ system and variables sensitive to its spin and parity are used to discriminate against the standard model background. Interference between pseudoscalar or scalar boson production and the standard model $\textrm{t}$$\bar{\textrm{t}}$ continuum is included, leading to peak-dip structures in the $m_{{\textrm{t}}{\bar{\textrm{t}}}}$ distribution. An excess of the data above the background prediction, based on perturbative quantum chromodynamics (QCD) calculations, is observed near the kinematic $\textrm{t}$$\bar{\textrm{t}}$ production threshold, while good agreement is found for high $m_{{\textrm{t}}{\bar{\textrm{t}}}}$. The data are consistent with the background prediction if the contribution from a simplified model of a color-singlet $^1\textrm{S}^{|1|}_0$ ${{\textrm{t}}{\bar{\textrm{t}}}}$ quasi-bound state $η_{\textrm{t}}$, inspired by nonrelativistic QCD, is added. Upper limits at 95% confidence level are set on the coupling between the pseudoscalar or scalar bosons and the top quark for boson masses in the range 365–1000 GeV, relative widths between 0.5% and 25%, and two background scenarios with or without $η_{\textrm{t}}$ contribution.

CMS↗
DOE OSTI2022

Discovery of the new isotope 251 Lr and the α-decay fine structure in 253 L

The products of the 203, 205 Tl ( 50 Ti, 2n) fusion-evaporation reactions were studied using the recently commissioned Argonne Gas-Filled Analyzer at Argonne National Laboratory. Two α-decay activities with energies of 9210(19) and 9246(19) keV and half-lives of $42^{+42}_{-14}$ and $24.4^{+7.0}_{-4.5}$ ms were observed which were followed by the known α decays of 247 Md and 243 Es . They are interpreted as originating from the 1/2 ⁻ [521] and 7/2 ⁻ [514] single-proton Nilsson states in the hitherto unknown isotope 251 Lr . From the measured Q α values the 1/2 ⁻ level was placed 117(27) keV above the 7/2 - level in 251 Lr in contrast to 255 Lr where the 1/2 ⁻ level is the lowest. Also, the α decay of 253 Lr was studied in more detail and a new α line at 8660(20) keV was found and a new half-life value of 2.46(32) s for an isomeric state in 253 Lr was measured. The 251, 253, 255 Lr Q α values were compared with predictions of various mass models. The relative energies of the 1/2 ⁻ [521] and 7/2 ⁻ [514] single-proton Nilsson states in 251, 253, 255 Lr isotopes were compared with results of the cranking shell model with pairing treated using the particle-number-conserving method. The level separation and, in particular, the level order change between 251 Lr and 255 Lr was reproduced only when the hexacontetrapole deformation ϵ 6 was included in the calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Almost medium-free measurement of the Hoyle state direct-decay component with a TPC

The structure of the Hoyle state, a highly α -clustered state at 7.65 MeV in 12 C, has long been the subject of debate. Understanding if the system comprises of three weakly interacting α particles in the 0s orbital, known as an α-condensate state, is possible by studying the decay branches of the Hoyle state. The direct decay of the Hoyle state into three α particles, rather than through the 8 Be ground state, can be identified by studying the energy partition of the three α particles arising from the decay. This paper provides details on the breakup mechanism of the Hoyle stating using a new experimental technique. Method: By using β-delayed charged-particle spectroscopy of 12 N using the Texas active target time-projection chamber, a high-sensitivity measurement of the direct 3α decay ratio can be performed without contributions from pileup events. A Bayesian approach to understanding the contribution of the direct components via a likelihood function shows that the direct component is <0.043% at the 95% confidence level. This value is in agreement with several other studies, and, here, we can demonstrate that a small nonsequential component with a decay fraction of about 10 –4 is most likely. Here, the measurement of the non-sequential component of the Hoyle state decay is performed in an almost medium-free reaction for the first time. The derived upper limit is in agreement with previous studies and demonstrates sensitivity to the absolute branching ratio. Further experimental studies would need to be combined with robust microscopic theoretical understanding of the decay dynamics to provide additional insight into the idea of the Hoyle state as an α condensate.

6 ≤ A ≤ 19↗
DOE OSTI2020

The CUORE Detector and Results

The cryogenic underground observatory for rare events (CUORE) is a cryogenic experiment searching for neutrinoless double beta decay (0νββ) of 130 Te. The detector consists of an array of 988TeO 2 crystals arranged in a compact cylindrical structure of 19 towers. We report the CUORE initial operations and optimization campaigns. We then present the CUORE results on 0νββ and 2νββ decay of 130 Te obtained from the analysis of the physics data acquired in 2017.

0vββ↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

Czochralski Growth and Characterization of a Compositionally Complex Rare Earth Aluminum Garnet Scintillator: (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce

Compositionally complex oxides have garnered increasing interest for their enhanced phase stability and tunable functional properties, yet their development as bulk single crystal scintillators remains limited. Herein, we report the Czochralski growth and characterization of (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce (GYTLAG), a compositionally complex garnet incorporating four dodecahedrally coordinated principal rare earth elements. The garnet phase was confirmed by powder and single crystal X-ray diffraction, and macroscopic defects are described. X-ray absorption near-edge structure measurements confirm the 3+ oxidation state of all rare earths and support their occupation of the same crystallographic site; white line intensity variations correlate with the anticipated segregation behavior. Elemental segregation is quantified by SEM/EDS and ICP-OES, and a linear trend was established between the segregation coefficient and the difference between each rare earth’s ionic radius (r) and the average ionic radius (AIR) of the dodecahedral site. This trend offers a predictive framework for compositional control in future REAG crystals grown by the Czochralski method. Photoluminescence and radioluminescence measurements reveal both Ce 3+ and Tb 3+ emission. Scintillation pulses exhibit four-component decay with dominant ~230 µs and ~1.2 ms components, and the light yield is estimated to be up to 43,000 ph/MeV under 137Cs γ-ray excitation. GYTLAG also demonstrates a strong radioluminescence efficiency and 50% lower afterglow at 20 ms compared to a LuAG:Ce reference, underscoring its promise for scintillator applications.

Compositionally complex oxide, high entropy oxide,↗
DOE OSTI2025

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Effective field theory for radiative corrections to charged-current processes. II. Axial-vector coupling

Here we discuss the hadronic structure-dependent radiative corrections to the axial-vector coupling that controls single-nucleon weak charged-current processes—commonly denoted by 𝑔 𝐴 . We match the Standard Model at the GeV scale onto chiral perturbation theory at next-to-leading order in the one-nucleon sector, in the presence of electromagnetic and weak interactions. As a result, we provide a representation for the corrections to 𝑔 𝐴 in terms of infrared finite convolutions of simple kernels with the single-nucleon matrix elements of time-ordered products of two and three quark bilinears (vector, axial-vector, and pseudoscalar). We discuss strategies to determine the required nonperturbative input from data, lattice-QCD (+QED), and possibly hadronic models. This work paves the way for a precise comparison of the values of the ratio 𝑔 𝐴 ⁡/𝑔 𝑉 extracted from experiment and from lattice QCD, which constrain physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Shear Behavior of Artificial Clay Seams within Bedded Salt Structures

Bedded salt contains thin layers of clay, also known as clay seams, in-between far thicker layers of salt. These inhomogeneities are thought to have first-order effects on the closure of nearby drifts and potential roof collapses. Despite their importance, characterizations of the peak shear strength and residual shear strength of clay seams in salt are extremely rare in the published literature. A previous paper reported results from laboratory direct shear experiments on clay seam samples from the Permian Basin in New Mexico. These clay seams behaved similar to intact salt, which was attributed to the abundance of salt crystals intersecting the clay seams. None of those specimens contained a distinct ¼" - ½" (6 -12 mm) thick clay seam, as has been observed in drifts at the Waste Isolation Pilot Plant (WIPP). Due to the difficulty in obtaining WIPP samples with these types of clay seams, artificial seams of bentonite and brine sandwiched between sections of salt were created and shear tested. Eight 4" diameter samples were created with either a ¼" or ½" a thick seam and then consolidated at 3000 psi prior to shear testing. The direct shear tests on these samples were performed at nominal normal stresses representative of expected WIPP in-situ conditions (500 to 1500 psi). The resulting shear stress vs. shear displacement curves exhibited a peak followed by a gradual decay of shear strength. The shear stress never transitioned to a true residual shear stress plateau, so the final shear strength at the end of each test (0.75" of shear displacement) was analyzed instead. Both the peak shear strength and the final shear strength conformed to Mohr- Coulomb behavior with friction angles and cohesion strengths consistent with a saturated, highly consolidated, clay. These new artificial clay seam results and the previous clay-interspersed-with-salt results likely bound the expected shear behavior of WIPP clay seams.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Second Target Station High-Fidelity Target Activation Comparison

The development of the Second Target Station (STS) target system at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) is well underway. The target system at STS consists of a rotating target disk that contains 21 segments of tungsten clad in tantalum clad in steel. A key aspect of the design of the target system is to account for the delayed heating and material damage caused by the delayed dose from decaying radionuclides. These radionuclides are a product of either spallation reactions or transmutation of the nuclei in the target system. These radionuclides build up in the target system components over the lifetime of the facility, and the radiation that is emitted can deposit energy in the components causing significant component heating and material damage. Monte Carlo N-Particle (MCNP) Version 6.2 transports the various particle species and calculates the spallation products and neutron fluxes throughout the target system. These spallation products and neutron fluxes along with the material definition of each component are relayed to the CINDER2008 transmutation code to calculate the radionuclide inventories and the corresponding decay gamma emission spectra. MCNP6.2 coupled with CINDER2008 is the computational method-of-choice for the analysis discussed in the following sections of this report. The analysis focuses on validating major assumptions in calculating the radionuclide inventory in the STS target system: all of the target segments are fresh, unirradiated material when the protons are incident on the segment, the average of the 21 segments of the target is sufficient to represent a single segment, and that averaging the proton pulse structure over time does not significantly affect the radionuclide inventory. The position-averaged, high-fidelity, and single-tally computational methods are used to validate the assumptions and provide a point of comparison to evaluate how the assumptions impact the radionuclide inventories. A more detailed explanation of the three computational methods is provided in Section 2. The position-averaged and single-tally methods are less computationally expensive when compared with the high-fidelity method where 54,000 individual calculations are needed to calculate 1 hr of STS operation. Section 3 details the comparison of the three methods to show that the assumptions made in the position-averaged method do not significantly impact the radionuclide inventory after 1 hr of operation. The discussions and results in this report are for 1 hr of operation. Due to the computational cost associated with calculating the transmutation and activation using the high-fidelity method, only 1 hr of operation has been calculated. The discrepancies observed after 1 hr of operation are not extrapolated out to longer operational times, and this report does not address how the discrepancies between the computational methods may manifest for longer operational periods.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

New method for level-lifetime measurements with thick scintillators

Level lifetimes provide key insight into the structure of atomic nuclei and serve as stringent tests of theoretical descriptions. Though several methods for determining level lifetimes exist for both reaction measurements and decay studies, here the focus is on techniques involving the direct measurement of time differences between population and subsequent depopulation of excited states. The techniques presented herein are broadly applicable across multiple timing ranges, but the approach is specifically described for the βγ timing method. A multi-step, amplitude-dependent time walk correction was employed to address the data analysis complications that arise from using thick scintillators for electron detection. Further, additional corrections for the depth of the interaction must also be performed when the parent isotope has been implanted into the detector. Techniques are presented for performing these time walk and depth of interaction corrections. Subsequently, a new Monte Carlo method utilizing measured detector responses obtained from the data, coupled with chi-square minimization, is presented for extracting excited state lifetimes ≳ 100 picoseconds. The framework of this Monte Carlo method is developed for the decay of a state in 68 Zn with a known 1.6 ps half life, which is considered prompt given the detection sensitivity of the technique, and then benchmarked using two other excited states in neutron-rich Ni isotopes with 120(34) ps and 1.05(3) ns half lives. Using this new method which takes into account the thick scintillator used, these same half lives were measured to be 135(10) ps and 1.04(24) ns, respectively. The overall good agreement demonstrates the validity of the technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Time-resolved vacuum-ultraviolet photoelectron spectroscopy of the à 1 A u state of acetylene

Ultrafast time-resolved photoelectron spectra are reported for the vacuum-ultraviolet (VUV) photoionization of acetylene following excitation to the à 1 A u state via UV absorption at 200 nm. The excitation energy lies above the lowest dissociation threshold to C 2 H X̃ 2 Σ + + H, as well as above the threshold for adiabatic dissociation of the à 1 A u state to form C 2 H (à 2 Π) + H. The time-dependent mass spectra and photoelectron spectra provide insight into the intramolecular decay processes of the à 1 A u state. In addition, photoelectron spectra of the à 1 A u state with VUV light access both the X̃ 2 Π u and à 2 Σ g + states of the ion, as well as the predicted, but previously unobserved, 1 2 Π g state, which corresponds to a two-hole, one-particle configuration that lies in close proximity to the à 2 Σ g + state. The 1 2 Π g state is split into 2 A 2 + 2 B 2 and 2 A g + 2 B g states in the cis and trans configurations, respectively. In conclusion, electronic structure calculations, along with trajectory calculations, reproduce the principal features of the experimental data and confirm the assignment of the 1 2 Π g state.

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

CUORE: The first bolometric experiment at the ton scale for the search for neutrino-less double beta decay

The Cryogenic Underground Observatory for Rare Events (CUORE) is the most massive bolometric experiment searching for neutrino-less double beta (0νββ) decay. The detector consists of an array of 988 TeO 2 crystals (742 kg) arranged in a compact cylindrical structure of 19 towers. Here, this paper will describe the CUORE experiment, including the cryostat, and present the detector performance during the first year of running. Additional detail will describe the effort made in improving the energy resolution in the 130 Te 0νββ decay region of interest (ROI) and the suppression of backgrounds. A description of work to lower the energy threshold in order to give CUORE the sensitivity to search for other rare events, such as dark matter, will also be provided.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Photoproduction of the f 2 ( 1270 ) Meson Using the CLAS Detector

The quark structure of the f 2 ( 1270 ) meson has, for many years, been assumed to be a pure quark-antiquark (q$\bar{q}$) resonance with quantum numbers J PC = 2 ++ . Recently, it was proposed that the f 2 ( 1270 ) is a molecular state made from the attractive interaction of two ρ mesons. Such a state would be expected to decay strongly to final states with charged pions due to the dominant decay ρ → π + π - , whereas decay to two neutral pions would likely be suppressed. Here, we measure for the first time the reaction γp → π 0 π 0 p , using the CEBAF Large Acceptance Spectrometer detector at Jefferson Lab for incident beam energies between 3.6 and 5.4 GeV. Differential cross sections, dσ/dt , for f 2 ( 1270 ) photoproduction are extracted with good precision due to low backgrounds and are compared to theoretical calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Observation of new excited ${B} ^0_{s}$ states

A structure is observed in the ${B} ^{\pm }{K} ^{\mp }$ mass spectrum in a sample of proton–proton collisions at centre-of-mass energies of 7, 8, and 13 TeV, collected with the LHCb detector and corresponding to a total integrated luminosity of 9$\,\text {fb} ^{-1}$. The structure is interpreted as the result of overlapping excited ${B} ^0_{s}$ states. With high significance, a two-peak hypothesis provides a better description of the data than a single resonance. Under this hypothesis the masses and widths of the two states, assuming they decay directly to ${B} ^{\pm }{K} ^{\mp }$, are determined to be $\begin{aligned} m_1&= 6063.5 \pm 1.2 \text { (stat)} \pm 0.8\text { (syst)}\,\text {Me}\text {V}, \\ \Gamma _1&= 26 \pm 4 \text { (stat)} \pm 4\text { (syst)}\,\text {Me}\text {V}, \\ m_2&= 6114 \pm 3 \text { (stat)} \pm 5\text { (syst)}\,\text {Me}\text {V}, \\ \Gamma _2&= 66 \pm 18 \text { (stat)} \pm 21\text { (syst)}\,\text {Me}\text {V}. \end{aligned}$ Alternative values assuming a decay through ${B} ^{*\pm }{K} ^{\mp }$, with a missing photon from the $B^{*\pm } \rightarrow B^{\pm }\gamma $ decay, which are shifted by approximately 45 $\,\text {Me}$V, are also determined. The possibility of a single state decaying in both channels is also considered. The ratio of the total production cross-section times branching fraction of the new states relative to the previously observed $B_{s2}^{*0}$ state is determined to be $0.87 \pm 0.15 \text { (stat)} \pm 0.19 \text { (syst)}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗