Search

Quantum computers must achieve large-scale, fault-tolerant operation to deliver on their promise of transformational processing power [1-4]. This will require thousands or millions of high-fidelity quantum gates and similar numbers of qubits [5]. Demonstrations using neutral-atom qubits trapped and manipulated by lasers have shown that this modality can provide high two-qubit gate (CZ) fidelities and scalable operation [6-10]. However, the gates in these demonstrations are driven by lasers that do not resolve individual qubits, with universal computation enabled by physical mid-circuit shuttling of the qubits. This relatively slow operation will greatly extend runtimes for useful, large-scale computation. Here we demonstrate a universal neutral-atom quantum computer with gate rates limited by optical switching times, rather than shuttling, by individually addressing tightly focused laser beams at an array of single atoms. We achieve CZ fidelity of 99.35(4)% and local single qubit RZ gate fidelity of 99.902(8)%. Moreover, we demonstrate non-destructive readout of alkali-atom qubits with sub-percent loss, which boosts operational speed. This technique also enables us to measure 99.73(3)% CZ fidelity with atom-loss events excluded, which is a record among long lived neutral-atom qubits and highlights the path to higher fidelity and error correction. Our results represent a critical step towards large-scale, fault-tolerant neutral-atom quantum computers that can execute computations on practical timescales.

The $^{244}$Pu($^{50}$Ti,$xn$)$^{294-x}$Lv reaction was investigated at Lawrence Berkeley National Laboratory's 88-Inch Cyclotron facility. The experiment was aimed at the production of a superheavy element with $Z\ge 114$ by irradiating an actinide target with a beam heavier than $^{48}$Ca. Produced Lv ions were separated from the unwanted beam and nuclear reaction products using the Berkeley Gas-filled Separator and implanted into a newly commissioned focal plane detector system. Two decay chains were observed and assigned to the decay of $^{290}$Lv. The production cross section was measured to be $\sigma_{\rm prod}=0.44(^{+58}_{-28})$~pb at a center-of-target center-of-mass energy of 220(3)~MeV. This represents the first published measurement of the production of a superheavy element near the `Island-of-Stability', with a beam of $^{50}$Ti and is an essential precursor in the pursuit of searching for new elements beyond $Z=118$., Comment: Submitted to Physical Review Letters

Information on the $^{38}$S level scheme was expanded through experimental work utilizing a fusion-evaporation reaction and in-beam $\gamma$-ray spectroscopy. Prompt $\gamma$-ray transitions were detected by the Gamma-Ray Energy Tracking Array (GRETINA) and recoiling $^{38}$S residues were selected by the Fragment Mass Analayzer (FMA). Tools based on machine-learning techniques were developed and deployed for the first time in order to enhance the unique selection of $^{38}$S residues and identify any associated $\gamma$-ray transitions. The new level information, including the extension of the even-spin yrast sequence through $J^{\pi} = 8^{(+)}$, was interpreted in terms of a basic single-particle picture as well shell-model calculations which incorporated the empirically derived FSU interaction. A comparison between the properties of the yrast states in the even-$Z$ $N=22$ isotones from $Z=14$ to $20$, and for $^{36}$Si-$^{38}$S in particular, was also presented with an emphasis on the role and influence of the neutron $1p_{3/2}$ orbital on the structure in the region.

Excited-state spectroscopy from the first Facility for Rare Isotope Beams (FRIB) experiment is reported. A 24(2)-$\mu$s isomer was observed with the FRIB Decay Station initiator (FDSi) through a cascade of 224- and 401-keV $\gamma$ rays in coincidence with $^{32}\textrm{Na}$ nuclei. This is the only known microsecond isomer ($1{\text{ }\mu\text{s}}\leq T_{1/2} < 1\text{ ms}$) in the region. This nucleus is at the heart of the $N=20$ island of shape inversion and is at the crossroads of spherical shell-model, deformed shell-model, and ab initio theories. It can be represented as the coupling of a proton hole and neutron particle to $^{32}\textrm{Mg}$, $^{32}\textrm{Mg}+\pi^{-1} + \nu^{+1}$. This odd-odd coupling and isomer formation provides a sensitive measure of the underlying shape degrees of freedom of $^{32}\textrm{Mg}$, where the onset of spherical-to-deformed shape inversion begins with a low-lying deformed $2^+$ state at 885 keV and a low-lying shape-coexisting $0_2^+$ state at 1058 keV. We suggest two possible explanations for the 625-keV isomer in $^{32}$Na: a $6^-$ spherical shape isomer that decays by $E2$ or a $0^+$ deformed spin isomer that decays by $M2$. The present results and calculations are most consistent with the latter, indicating that the low-lying states are dominated by deformation., Comment: 7 pages, 5 figures, accepted by Physical Review Letters

We introduce SupercheQ, a family of quantum protocols that achieves asymptotic advantage over classical protocols for checking the equivalence of files, a task also known as fingerprinting. The first variant, SupercheQ-EE (Efficient Encoding), uses n qubits to verify files with 2^O(n) bits -- an exponential advantage in communication complexity (i.e. bandwidth, often the limiting factor in networked applications) over the best possible classical protocol in the simultaneous message passing setting. Moreover, SupercheQ-EE can be gracefully scaled down for implementation on circuits with poly(n^l) depth to enable verification for files with O(n^l) bits for arbitrary constant l. The quantum advantage is achieved by random circuit sampling, thereby endowing circuits from recent quantum supremacy and quantum volume experiments with a practical application. We validate SupercheQ-EE's performance at scale through GPU simulation. The second variant, SupercheQ-IE (Incremental Encoding), uses n qubits to verify files with O(n^2) bits while supporting constant-time incremental updates to the fingerprint. Moreover, SupercheQ-IE only requires Clifford gates, ensuring relatively modest overheads for error-corrected implementation. We experimentally demonstrate proof-of-concepts through Qiskit Runtime on IBM quantum hardware. We envision SupercheQ could be deployed in distributed data settings, accompanying replicas of important databases.

The collective structure of $^{106}$Cd is elucidated by multi-step Coulomb excitation of a 3.849 MeV/$A$ beam of $^{106}$Cd on a 1.1 mg/cm$^2$ $^{208}$Pb target using GRETINA-CHICO2 at ATLAS. Fourteen $E2$ matrix elements were obtained. The nucleus $^{106}$Cd is a prime example of emergent collectivity that possesses a simple structure: it is free of complexity caused by shape coexistence and has a small, but collectively active number of valence nucleons. This work follows in a long and currently active quest to answer the fundamental question of the origin of nuclear collectivity and deformation, notably in the cadmium isotopes. The results are discussed in terms of phenomenological models, the shell model, and Kumar-Cline sums of $E2$ matrix elements. The ${\langle 0_2^+ ||E2||2_1^+ \rangle}$ matrix element is determined for the first time, providing a total, converged measure of the electric quadrupole strength, $\langle Q^2 \rangle$, of the first-excited $2_1^+$ level relative to the $0_1^+$ ground state, which does not show an increase as expected of harmonic and anharmonic vibrations. Strong evidence for triaxial shapes in weakly collective nuclei is indicated; collective vibrations are excluded. This is contrary to the only other cadmium result of this kind in $^{114}$Cd by C. Fahlander et al., Nucl. Phys. A485, 327 (1988), which is complicated by low-lying shape coexistence near midshell., Comment: 7 pages, 4 figures

The $^{25}$F($5/2^+) (-1p) ^{24}$O reaction was studied at the NSCL using the S800 spectrometer. The experimental spectroscopic factor for the ground-state to ground-state transition indicates a substantial depletion of the proton $d_{5/2}$ strength compared to shell-model expectations. Our result supports the findings reported by Tang \textit{et al.}, from their study of the $(p,2p)$ reaction at RIBF. The overlap between the $^{25}$F and $^{24}$O ground-states is considerably less than anticipated if $^{24}$O acted as a robust and rigid doubly-magic core in $^{25}$F. We interpret the results within the framework of the Particle-Vibration Coupling (PVC) of a $d_{5/2}$ proton coupled to a quadrupole phonon of an effective core. This approach provides a good description of the experimental data by requiring an effective $^{24}$O* core with a phonon energy of $\hbar\omega_2$= 3.2 MeV, and a $B(E2) ~ 2.7$ W.u., softer and more collective than a bare $^{24}$O. Both the Nilsson deformed mean field and the PVC models appear to capture the properties of the effective core of $^{25}$F, suggesting that the additional proton tends to polarize the free, doubly magic $^{24}$O in such a way that it becomes either slightly deformed or a quadrupole vibrator., Comment: 6 pages, 2 figures

Gate model quantum computers promise to solve currently intractable computational problems if they can be operated at scale with long coherence times and high fidelity logic. Neutral atom hyperfine qubits provide inherent scalability due to their identical characteristics, long coherence times, and ability to be trapped in dense multi-dimensional arrays\cite{Saffman2010}. Combined with the strong entangling interactions provided by Rydberg states\cite{Jaksch2000,Gaetan2009,Urban2009}, all the necessary characteristics for quantum computation are available. Here we demonstrate several quantum algorithms on a programmable gate model neutral atom quantum computer in an architecture based on individual addressing of single atoms with tightly focused optical beams scanned across a two-dimensional array of qubits. Preparation of entangled Greenberger-Horne-Zeilinger (GHZ) states\cite{Greenberger1989} with up to 6 qubits, quantum phase estimation for a chemistry problem\cite{Aspuru-Guzik2005}, and the Quantum Approximate Optimization Algorithm (QAOA)\cite{Farhi2014} for the MaxCut graph problem are demonstrated. These results highlight the emergent capability of neutral atom qubit arrays for universal, programmable quantum computation, as well as preparation of non-classical states of use for quantum enhanced sensing., Comment: 4 figures and supplemental material. V2 details on phase estimation and RB data included

Exclusive cross sections for the $^{43}$P$(-1p)^{42}$Si reaction to the lowest $0^+$ and $2^+$ states, measured at NSCL with GRETINA and the S800, are interpreted in terms of a two-level mixing (collective) model of oblate and prolate co-existing shapes. Using the formalism developed for deformed nuclei we calculate the spectroscopic amplitudes and exclusive cross-sections in the strong coupling limit, where for $^{43}$P the schematic wavefunction includes the coupling of the Nilsson [211]$\frac{1}{2}$ proton orbit. Good agreement with the experimental data is obtained when the amplitude of the oblate configuration is $\gtrsim$ 80\%, suggesting that both nuclei are predominantly oblate, in line with theoretical expectations., Comment: 5 pages, 3 figures

UCGretina, a GEANT4 simulation of the GRETINA gamma-ray tracking array of highly-segmented high-purity germanium detectors is described. We have developed a model of the array, in particular of the Quad Module and the capsules, that gives good agreement between simulated and measured photopeak efficiencies over a broad range of gamma-ray energies and reproduces the shape of the measured Compton continuum. Both of these features are needed in order to accurately extract gamma-ray yields from spectra collected in in-beam gamma-ray spectroscopy measurements with beams traveling at $v/c \gtrsim 0.3$ at the National Superconducting Cyclotron Laboratory and the Facility for Rare Isotope Beams. In the process of developing the model, we determined that millimeter-scale layers of passive germanium surrounding the active volumes of the simulated crystals must be included in order to reproduce measured photopeak efficiencies. We adopted a simple model of effective passive layers and developed heuristic methods of determining passive-layer thicknesses by comparison of simulations and measurements for a single crystal and for the full array. Prospects for future development of the model are discussed., Comment: 26 pages, 12 figures

Spectroscopic factors to low-lying negative-parity states in $^{11}$Be extracted from the $^{12}$B($d$,$^3$He)$^{11}$Be proton-removal reaction are interpreted within the rotational model. Earlier predictions of the $p$-wave proton removal strengths in the strong coupling limit of the Nilsson model underestimated the spectroscopic factors to the $3/2^-_1$ and $5/2^-_1$ states and suggested that deviations in the $1^+$ ground state of the odd-odd $^{12}$B due to Coriolis coupling should be further explored. In this work we use the Particle Rotor Model to take into account these effects and obtain a good description of the level scheme in $^{11}$B, with a moderate $K$-mixing of the proton Nilsson levels [110]1/2 and [101]3/2. This mixing, present in the $1^+$ bandhead of $^{12}$B, is key to explaining the proton pickup data., Comment: 5 Pages, 2 Figures, to be submitted to Phys. Rev. C

The low-spin structure of the semimagic 64Ni nucleus has been considerably expanded: combining four experiments, several 0+ and 2+ excited states were identified below 4.5 MeV, and their properties established. The Monte Carlo shell model accounts for the results and unveils an unexpectedly complex landscape of coexisting shapes: a prolate 0+ excitation is located at a surprisingly high energy (3463 keV), with a collective 2+ state 286 keV above it, the first such observation in Ni isotopes. The evolution in excitation energy of the prolate minimum across the neutron N = 40 subshell gap highlights the impact of the monopole interaction and its variation in strength with N.

In a recent experiment, carried out at RIBF/RIKEN, the $^{25}$F$(p,2p)$$^{24}$O reaction was studied at 270 MeV/A in inverse kinematics. Derived spectroscopic factors suggest that the effective core of $^{25}$F significantly differs from a free $^{24}$O nucleus. We interpret these results within the Particle-Rotor Model and show that the experimental level scheme of $^{25}$F can be understood in the rotation-aligned coupling scheme, with its $5/2^+_1$ ground state as the band-head of a decoupled band. The excitation energies of the observed $1/2_1^+$ and $9/2_1^+$ states correlate strongly with the rotational energy of the effective core, seen by the odd proton, and allow us to estimate its $2^+$ energy at $\approx$ 3.2 MeV and a moderate quadrupole deformation, $\epsilon_2 \approx 0.15$. The measured fragmentation of the $\pi d_{5/2}$ single-particle strength is discussed and some further experiments suggested., Comment: 5 pages, 3 figures. Submitted to Physical Review C

The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number $Z$ along $N=20$ is of interest both in terms of better understanding the evolution of nuclear structure between spherical $^{40}$Ca and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the $N=20$ Island of Inversion, the npnh (where n=2,4,6) neutron excitations across a diminished $N=20$ gap result in deformed and collective ground states, as observed in $^{32}$Mg. In heavier isotones, npnh excitations do not dominate in the ground states, but are present in the relatively low-lying level schemes. With the aim of identifying the expected 2p2h$\otimes\mathrm{s}_{1/2^+}$ state in $^{35}$P, the only $N=20$ isotone for which the neutron 2p2h excitation bandhead has not yet been identified, the $^{36}$S(d,$^3$He)$^{35}$P reaction has been revisited in inverse kinematics with the HELical Orbit Spectrometer (HELIOS) at the Argonne Tandem Linac Accelerator System (ATLAS). While a candidate state has not been located, an upper limit for the transfer reaction cross-section to populate such a configuration within a 2.5 to 3.6\,MeV energy range, provides a stringent constraint on the wavefunction compositions in both $^{36}$S and $^{35}$P.

New experimental data obtained from $\gamma$-ray tagged one-neutron and one-proton knockout from $^{55}$Co is presented. A candidate for the sought-after $T=1, T_z = 0, J^{\pi} = 6^+$ state in $^{54}$Co is proposed based on a comparison to the new data on $^{54}$Fe, the corresponding observables predicted by large-scale-shell-model (LSSM) calculations in the full $fp$-model space employing charge-dependent contributions, and isospin-symmetry arguments. Furthermore, possible isospin-symmetry breaking in the $A=54$, $T=1$ triplet is studied by calculating the experimental $c$ coefficients of the isobaric mass multiplet equation (IMME) up to the maximum possible spin $J=6$ expected for the $(1f_{7/2})^{-2}$ two-hole configuration relative to the doubly-magic nucleus $^{56}$Ni. The experimental quantities are compared to the theoretically predicted $c$ coefficients from LSSM calculations using two-body matrix elements obtained from a realistic chiral effective field theory potential at next-to-next-to-next-to-leading order (N$^3$LO)., Comment: 6 pages, 5 figures. Work has been published

Imperfections in analog-to-digital conversion (ADC) cannot be ignored when signal digitization requirements demand both wide dynamic range and high resolution, as is the case for the Majorana Demonstrator 76Ge neutrinoless double beta decay search. Enabling the experiment's high-resolution spectral analysis and efficient pulse shape discrimination required careful measurement and correction of ADC nonlinearites. A simple measurement protocol was developed that did not require sophisticated equipment or lengthy data taking campaigns. A slope-dependent hysteresis was observed and characterized. A correction applied to digitized waveforms prior to signal processing reduced the differential and integral nonlinearites by an order of magnitude, eliminating these as dominant contributions to the systematic energy uncertainty at the double-beta decay Q value.

The neutron-rich Ni isotopes have attracted attention in recent years due to the occurrence of shape or configuration coexistence. We report on the difference in population of excited final states in 70Ni following gamma-ray tagged one-proton, one-neutron, and two-proton knockout from 71Cu, 71Ni, and 72Zn rare-isotope beams, respectively. Using variations observed in the relative transition intensities, signaling the changed population of specific final states in the different reactions, the role of neutron and proton configurations in excited states of 70Ni is probed schematically, with the goal of identifying those that carry, as leading configuration, proton excitations across the Z = 28 shell closure. Such states are suggested in the literature to form a collective structure associated with prolate deformation. Adding to the body of knowledge for 70Ni, 29 new transitions are reported, of which 15 are placed in its level scheme., Comment: accepted for publication in Phys. Rev. C

We present a high-resolution in-beam $\gamma$-ray spectroscopy study of excited states in the mirror nuclei $^{55}$Co and $^{55}$Ni following one-nucleon knockout from a projectile beam of $^{56}$Ni. The newly determined partial cross sections and the $\gamma$-decay properties of excited states provide a test of state-of-the-art nuclear structure models and probe mirror symmetry in unique ways. A mirror asymmetry for the partial cross sections leading to the two lowest $3/2^-$ states in the $A = 55$ mirror pair was identified as well as a significant difference in the $E1$ decays from the $1/2^+_1$ state to the same two $3/2^-$ states. The mirror asymmetry in the partial cross sections cannot be reconciled with the present shell-model picture or small mixing introduced in a two-state model. The observed mirror asymmetry in the $E1$ decay pattern, however, points at stronger mixing between the two lowest $3/2^-$ states in $^{55}$Co than in its mirror $^{55}$Ni.

In the late stages of stellar core-collapse, prior to core bounce, electron captures on medium-heavy nuclei drive deleptonization and simulations require the use of accurate reaction rates. Nuclei with neutron number near $N=50$, just above atomic number $Z=28$, play an important role, but rates used in astrophysical simulations rely primarily on a relatively simple single-state approximation. In order to improve the accuracy of astrophysical simulations, experimental data are needed to test the electron-capture rates and to guide the development of better theoretical models. This work presents the results of the $^{86}$Kr($t$,$^{3}$He+$\gamma$) experiment at the NSCL, from which an upper limit for the Gamow-Teller strength up to an excitation energy in $^{86}$Br of 5 MeV is extracted. The derived upper limit for the electron-capture rate on $^{86}$Kr indicates that the rate estimated through the single-state approximation is too high and that rates based on Gamow-Teller strengths estimated in shell-model and QRPA calculations are more accurate. The QRPA calculations tested in this manner were used for estimating the electron capture rates for 78 isotopes near $N=50$ and above $Z=28$. The impact of using these new electron-capture rates in simulations of supernovae instead of the rates based on the single-state approximation is investigated, indicating a significant reduction in the deleptonization that affects multi-messenger signals, such as the emission of neutrinos and gravitational waves., Comment: 15 pages, 10 figures

A more detailed test of the implementation of nuclear forces that drive shell evolution in the pivotal nucleus \nuc{42}{Si} -- going beyond earlier comparisons of excited-state energies -- is important. The two leading shell-model effective interactions, SDPF-MU and SDPF-U-Si, both of which reproduce the low-lying \nuc{42}{Si}($2^+_1$) energy, but whose predictions for other observables differ significantly, are interrogated by the population of states in neutron-rich \nuc{42}{Si} with a one-proton removal reaction from \nuc{43}{P} projectiles at 81~MeV/nucleon. The measured cross sections to the individual \nuc{42}{Si} final states are compared to calculations that combine eikonal reaction dynamics with these shell-model nuclear structure overlaps. The differences in the two shell-model descriptions are examined and linked to predicted low-lying excited $0^+$ states and shape coexistence. Based on the present data, which are in better agreement with the SDPF-MU calculations, the state observed at 2150(13)~keV in \nuc{42}{Si} is proposed to be the ($0^+_2$) level., Comment: accepted in Physical Review Letters

A model-independent technique was used to determine the $\gamma$-ray Strength Function ($\gamma$SF) of $^{56}$Fe down to $\gamma$-ray energies less than 1 MeV for the first time with GRETINA using the $(p,p')$ reaction at 16 MeV. No difference was observed in the energy dependence of the $\gamma$SF built on $2^{+}$ and $4^{+}$ final states, supporting the Brink hypothesis. In addition, angular distribution and polarization measurements were performed. The angular distributions are consistent with dipole radiation. The polarization results show a small bias towards magnetic character in the region of the enhancement.

The underrepresentation of African-American faculty in the US professoriate has persisted for some time. Relatedly, adjunct faculty remain a fast-growing sector of the professoriate. Adjunct faculty include "experts" and "specialists" who teach postsecondary courses with a narrow focus and with content tailored to their full-time employment. Using a qualitative narrative approach, I delineate ways I construct meaning for myself as a part-timer. I develop a cultural interpretation of adjunct teaching that provides alternative view of professional socialization to counter the ongoing challenge of increasing the number of Black faculty in higher education. In doing so, three themes from the data (ideal worker as adjunct, historical role models, and clinical value of course content) suggest the identity of part-time faculty is situated in personal, professional, and cultural experiences. [For the complete volume, "The Beauty and the Burden of Being a Black Professor. Diversity in Higher Education. Volume 24," see ED614149.]

Recent results from RIKEN/RIBF on the low-lying level structure of 29F are interpreted within the Particle-Rotor Model. We show that the experimental data can be understood in the Rotation-aligned Coupling Scheme, with the 5/2+ ground state as the bandhead of a decoupled band. In this picture, the energy of the observed 1/2+1 state correlates strongly with the rotational energy of the core and provides an estimate of the 2+ energy in 28O. Our analysis suggest a moderate deformation, epsilon_2 ~ 0.17, and places the 2+ in 28O at ~ 2.4 MeV., Comment: 3 pages, 3 figures, To be submitted to Phys. Lett. B

Spectroscopic factors in 10Be, 11Be and 12Be, extracted from (d,p), one neutron knockout, and (p,d) reactions are interpreted within the rotational model. Assuming that the ground state and first excited state of 11Be can be associated with the 1/2[220] and 1/2[101] Nilsson levels, the strong coupling limit gives simple expressions that relate the amplitudes of these wavefunctions (in the spherical basis) with the measured cross-sections and derived spectroscopic factors. We obtain good agreement with both the measured magnetic moment of the ground state in 11Be and the reaction data., Comment: 4 pages, 1 figure

The calcium isotopes have emerged as an important testing ground for new microscopically derived shell-model interactions, and a great deal of focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Ca following one-neutron knockout reactions from $^{48,50}$Ca. The observed reduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca, as compared to $^{47}$Ca, is consistent with the description given by shell-model calculations based on two- and three-nucleon forces in the neutron $pf$ model space, implying a fragmentation of the $l$=3 strength to higher-lying states. The experimental result is inconsistent with both the GXPF1 interaction routinely used in this region of the nuclear chart and with microscopic calculations in an extended model space including the $\nu1g_{9/2}$ orbital., Comment: 6 pages, 2 figures

The dynamics of superfluid 4He at and above the Landau quasiparticle regime is investigated by high precision inelastic neutron scattering measurements of the dynamic structure factor. A highly structured response is observed above the familiar phonon-maxon-roton spectrum, characterized by sharp thresholds for phonon-phonon, maxon-roton and roton-roton coupling processes. The experimental dynamic structure factor is compared to the calculation of the same physical quantity by a Dynamic Many-body theory including three-phonon processes self-consistently. The theory is found to provide a quantitative description of the dynamics of the correlated bosons for energies up to about three times that of the Landau quasiparticles., Comment: 5 pages, 3 figures

The neutron-rich nucleus $^{144}$Ba ($t_{1/2}$=11.5 s) is expected to exhibit some of the strongest octupole correlations among nuclei with mass numbers $A$ less than 200. Until now, indirect evidence for such strong correlations has been inferred from observations such as enhanced $E1$ transitions and interleaving positive- and negative-parity levels in the ground-state band. In this experiment, the octupole strength was measured directly by sub-barrier, multi-step Coulomb excitation of a post-accelerated 650-MeV $^{144}$Ba beam on a 1.0-mg/cm$^2$ $^{208}$Pb target. The measured value of the matrix element, $\langle 3_1^- \| \mathcal{M}(E3) \| 0_1^+ \rangle=0.65(^{+17}_{-23})$ $e$b$^{3/2}$, corresponds to a reduced $B(E3)$ transition probability of 48($^{+25}_{-34}$) W.u. This result represents an unambiguous determination of the octupole collectivity, is larger than any available theoretical prediction, and is consistent with octupole deformation., Comment: 5 pages, 2 figures, 1 table, Physical Review Letters 2016