The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section of $\phi$(1020)-meson production at forward rapidity in $p$$+$$p$ collisions at $\sqrt{s}=$510 GeV via the dimuon decay channel. The partial cross section in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $2<p_T<7$ GeV/$c$ is $\sigma_\phi=[2.28 \pm 0.09\,{\rm (stat)} \pm 0.14\,{\rm (syst)} \pm 0.27\, {\rm (norm)}] \times 10^{-2}$~mb. The energy dependence of $\sigma_\phi$ ($1.2<|y|<2.2, \; 2<p_T<5$ GeV/$c$) is studied using the PHENIX measurements at $\sqrt{s}=$200 and 510 GeV and the Large-Hadron-Collider measurements at $\sqrt{s}=$2.76 and 7 TeV. The experimental results are compared to various event generator predictions ({\sc pythia6, pythia8, phojet, ampt, epos3,} and {\sc epos-lhc}).

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.

The experimental study of the collisions of heavy nuclei at relativistic energies has established the properties of the quark-gluon plasma (QGP), a state of hot, dense nuclear matter in which quarks and gluons are not bound into hadrons. In this state, matter behaves as a nearly inviscid fluid that efficiently translates initial spatial anisotropies into correlated momentum anisotropies among the produced particles, producing a common velocity field pattern known as collective flow. In recent years, comparable momentum anisotropies have been measured in small-system proton-proton ($p$$+$$p$) and proton-nucleus ($p$$+$$A$) collisions, despite expectations that the volume and lifetime of the medium produced would be too small to form a QGP. Here, we report on the observation of elliptic and triangular flow patterns of charged particles produced in proton-gold ($p$$+$Au), deuteron-gold ($d$$+$Au), and helium-gold ($^3$He$+$Au) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_{NN}}}$~=~200 GeV. The unique combination of three distinct initial geometries and two flow patterns provides unprecedented model discrimination. Hydrodynamical models, which include the formation of a short-lived QGP droplet, provide a simultaneous description of these measurements.

We conducted an improved search for the simultaneous capture of two $K$-shell electrons on the $^{124}$Xe and $^{126}$Xe nuclei with emission of two neutrinos using 800.0 days of data from the XMASS-I detector. A novel method to discriminate $\gamma$-ray/$X$-ray or double electron capture signals from $\beta$-ray background using scintillation time profiles was developed for this search. No significant signal was found when fitting the observed energy spectra with the expected signal and background. Therefore, we set the most stringent lower limits on the half-lives at $2.1 \times 10^{22}$ and $1.9 \times 10^{22}$ years for $^{124}$Xe and $^{126}$Xe, respectively, with 90% confidence level. These limits improve upon previously reported values by a factor of 4.5.

An annual modulation signal due to the Earth's rotation around the Sun would be one of the strongest indications of the direct detection of dark matter. In 2016, we reported a search for dark matter by looking for an annual modulation with the single-phase liquid xenon XMASS-I detector. That analysis indicated a slightly negative amplitude at low energy. In this work, we included more than one year of additional data, which more than doubles the exposure to 800 live days with our 832 kg target mass. When we assume Weakly Interacting Massive Particle (WIMP) dark matter elastically scattering on the xenon target, the exclusion upper limit for the WIMP-nucleon cross section was improved by a factor of two to 1.9$\times$10$^{-41}$cm$^2$ at 8 GeV/c$^2$ with our newly implemented data selection through a likelihood method. For the model-independent case, without assuming any specific dark matter model, we obtained more consistency with the null hypothesis than before with a $p$-value of 0.11 in the 1--20 keV energy region. This search is probing this region with an exposure that is larger than that of DAMA/LIBRA. We also did not find any significant amplitude in the data for periodicity with periods between 50 and 600 days.

We established a method to assay $^{210}$Pb and $^{210}$Po contaminations in the bulk of copper samples using a low-background alpha particle counter. The achieved sensitivity for the $^{210}$Pb and $^{210}$Po contaminations reaches a few mBq/kg. Due to this high sensitivity, the $^{210}$Pb and $^{210}$Po contaminations in oxygen free copper bulk were identified and measured for the first time. The $^{210}$Pb contaminations of our oxygen free copper samples were 17-40 mBq/kg. Based on our investigation of copper samples in each production step, the $^{210}$Pb in oxygen free copper was understood to be a small residual of an electrolysis process. This method to measure bulk contaminations of $^{210}$Pb and $^{210}$Po could be applied to other materials.

Recently, multiparticle-correlation measurements of relativistic $p/d/^3$He$+$Au, $p$$+$Pb, and even $p$$+$$p$ collisions have shown surprising collective signatures. Here we present beam-energy-scan measurements of 2-, 4-, and 6-particle angular correlations in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200, 62.4, 39, and 19.6 GeV. We also present measurements of 2- and 4-particle angular correlations in $p$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. We find the 4-particle cumulant to be real-valued for $d$$+$Au collisions at all four energies. We also find that the 4-particle cumulant in $p$$+$Au has the opposite sign as that in $d$$+$Au, indicating that the correlations are geometrical and therefore collective in origin. Further we find that the 6-particle cumulant agrees with the 4-particle cumulant in $d$$+$Au collisions at 200 GeV, which indicates that nonflow effects are subdominant and provides strong evidence of collectivity.

A search for the exotic meson $X(5568)$ decaying into the $B^0_s \pi^{\pm}$ final state is performed using data corresponding to $9.6 \textrm{fb}^{-1}$ from $p{\bar p}$ collisions at $\sqrt{s} = 1960$ GeV recorded by the Collider Detector at Fermilab. No evidence for this state is found and an upper limit of 6.7\% at the 95\% confidence level is set on the fraction of $B^0_s$ produced through the $X(5568) \rightarrow B^0_s \, \pi^{\pm}$ process.

We present measurements of the elliptic flow, $v_2$, as a function of transverse momentum ($p_T$), pseudorapidity ($\eta$), and centrality in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200, 62.4, 39, and 19.6 GeV. The beam-energy scan of $d$$+$Au collisions provides a testing ground for the onset of flow signatures in small collision systems. We measure a nonzero $v_2$ signal at all four collision energies, which, at midrapidity and low $p_T$, is consistent with predictions from viscous hydrodynamic models. Comparisons with calculations from parton transport models (based on the {\sc ampt} Monte Carlo generator) show good agreement with the data at midrapidity to forward ($d$-going) rapidities and low $p_T$. At backward (Au-going) rapidities and $p_T>1.5$ GeV/$c$, the data diverges from {\sc ampt} calculations of $v_2$ relative to the initial geometry, indicating the possible dominance of nongeometry related correlations, referred to as nonflow. We also present measurements of the charged-particle multiplicity ($dN_{\rm ch}/d\eta$) as a function of $\eta$ in central $d$$+$Au collisions at the same energies. We find that in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200 GeV the $v_2$ scales with $dN_{\rm ch}/d\eta$ over all $\eta$ in the PHENIX acceptance. At $\sqrt{s_{_{NN}}}=$62.4, and 39 GeV, $v_2$ scales with $dN_{\rm ch}/d\eta$ at midrapidity and forward rapidity, but falls off at backward rapidity. This departure from the $dN_{\rm ch}/d\eta$ scaling may be a further indication of nonflow effects dominating at backward rapidity.

During 2015 the Relativistic Heavy Ion Collider (RHIC) provided collisions of transversely polarized protons with Au and Al nuclei for the first time, enabling the exploration of transverse-single-spin asymmetries with heavy nuclei. Large single-spin asymmetries in very forward neutron production have been previously observed in transversely polarized $p$$+$$p$ collisions at RHIC, and the existing theoretical framework that was successful in describing the single-spin asymmetry in $p$$+$$p$ collisions predicts only a moderate atomic-mass-number ($A$) dependence. In contrast, the asymmetries observed at RHIC in $p$$+$$A$ collisions showed a surprisingly strong $A$ dependence in inclusive forward neutron production. The observed asymmetry in $p$$+$Al collisions is much smaller, while the asymmetry in $p$$+$Au collisions is a factor of three larger in absolute value and of opposite sign. The interplay of different neutron production mechanisms is discussed as a possible explanation of the observed $A$ dependence.

In theories with the large extra dimensions beyond the standard 4-dimensional spacetime, axions could propagate in such extra dimensions, and acquire Kaluza-Klein (KK) excitations. These KK axions are produced in the Sun and could solve unexplained heating of the solar corona. While most of the solar KK axions escape from the solar system, a small fraction is gravitationally trapped in orbits around the Sun. They would decay into two photons inside a terrestrial detector. The event rate is expected to modulate annually depending on the distance from the Sun. We have searched for the annual modulation signature using $832\times 359$ kg$\cdot$days of XMASS-I data. No significant event rate modulation is found, and hence we set the first experimental constraint on the KK axion-photon coupling of $4.8 \times 10^{-12}\, \mathrm{GeV}^{-1}$ at 90% confidence level for a KK axion number density of $\bar{n}_\mathrm{a} = 4.07 \times 10^{13}\, \mathrm{m}^{-3}$, the total number of extra dimensions $n = 2$, and the number of extra dimensions $\delta = 2$ that axions can propagate in.

We report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $p$+$p$ collisions at $\sqrt{s_{NN}}=200$~GeV. The $e^+e^-$ pair yield from $b\bar{b}$ and $c\bar{c}$ is separated by exploiting a double differential fit done simultaneously in dielectron invariant mass and $p_T$. We used three different event generators, {\sc pythia}, {\sc mc@nlo}, and {\sc powheg}, to simulate the $e^+e^-$ spectra from $c\bar{c}$ and $b\bar{b}$ production. The data can be well described by all three generators within the detector acceptance. However, when using the generators to extrapolate to $4\pi$, significant differences are observed for the total cross section. These difference are less pronounced for $b\bar{b}$ than for $c\bar{c}$. The same model dependence was observed in already published $d$+$A$ data. The $p$+$p$ data are also directly compared with $d$+$A$ data in mass and $p_T$, and within the statistical accuracy no nuclear modification is seen.

We observe variation in the resistively-detected nuclear magnetic resonance (RDNMR) lineshapes in quantum Hall breakdown. The breakdown is locally occurred in a gate-defined quantum point contact (QPC) region. Of particular interest is the observation of a dispersive lineshape occured when the bulk 2D electron gas (2DEG) is set to $\nu_{\rm{b}} = 2$ and the QPC filling factor to the vicinity of $\nu_{\rm{QPC}} = 1$, strikingly resemble the dispersive lineshape observed on a 2D quantum Hall state. This previously unobserved lineshape in a QPC points to simultaneous occurrence of two hyperfine-mediated spin flip-flop processes within the QPC. Those events give rise to two different sets of nuclei polarized in the opposite direction and positioned at a separate region with different degree of electronic polarizations.

We report on the results of a $Suzaku$ observation of the plasma in the filament located between the two massive clusters of galaxies Abell 399 and Abell 401. Abell 399 ($z$=0.0724) and Abell 401 ($z$=0.0737) are expected to be in the initial phase of a cluster merger. In the region between the two clusters, we find a clear enhancement in the temperature of the filament plasma from 4 keV (expected value from a typical cluster temperature profile) to $kT\sim$6.5 keV. Our analysis also shows that filament plasma is present out to a radial distance of 15' (1.3 Mpc) from a line connecting the two clusters. The temperature profile is characterized by an almost flat radial shape with $kT\sim$6-7 keV within 10' or $\sim$0.8 Mpc. Across $r$=8'~from the axis, the temperature of the filament plasma shows a drop from 6.3 keV to 5.1 keV, indicating the presence of a shock front. The Mach number based on the temperature drop is estimated to be ${\cal M}\sim$1.3. We also successfully determined the abundance profile up to 15' (1.3 Mpc), showing an almost constant value ($Z$=0.3 solar) at the cluster outskirt. We estimated the Compton $y$-parameter to be $\sim$14.5$\pm1.3\times10^{-6}$, which is in agreement with $Planck$'s results (14-17$\times10^{-6}$ on the filament). The line of sight depth of the filament is $l\sim$1.1 Mpc, indicating that the geometry of filament is likely a pancake shape rather than cylindrical. The total mass of the filamentary structure is $\sim$7.7$\times10^{13}~\rm M_{\odot}$. We discuss a possible interpretation of the drop of X-ray emission at the rim of the filament, which was pushed out by the merging activity and formed by the accretion flow induced by the gravitational force of the filament.

The cross section and transverse single-spin asymmetries of $\mu^{-}$ and $\mu^{+}$ from open heavy-flavor decays in polarized $p$+$p$ collisions at $\sqrt{s}=200$ GeV were measured by the PHENIX experiment during 2012 at the Relativistic Heavy Ion Collider. Because heavy-flavor production is dominated by gluon-gluon interactions at $\sqrt{s}=200$ GeV, these measurements offer a unique opportunity to obtain information on the trigluon correlation functions. The measurements are performed at forward and backward rapidity ($1.4<|y|<2.0$) over the transverse momentum range of $1.25<p_T<7$ GeV/$c$ for the cross section and $1.25<p_T<5$ GeV/$c$ for the asymmetry measurements. The obtained cross section is compared to a fixed-order-plus-next-to-leading-log perturbative-quantum-chromodynamics calculation. The asymmetry results are consistent with zero within uncertainties, and a model calculation based on twist-3 three-gluon correlations agrees with the data.

We report the first measurement of the full angular distribution for inclusive $J/\psi\rightarrow\mu^{+}\mu^{-}$ decays in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. The measurements are made for $J/\psi$ transverse momentum $2<p_{T}<10$ GeV/$c$ and rapidity $1.2<y<2.2$ in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient $\lambda_{\theta}$ is strongly negative at low $p_{T}$ and becomes close to zero at high $p_{T}$, while the azimuthal coefficient $\lambda_{\phi}$ is close to zero at low $p_{T}$, and becomes slightly negative at higher $p_{T}$. The frame-independent coefficient $\tilde{\lambda}$ is strongly negative at all $p_{T}$ in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

We report the first measurement of the fraction of $J/\psi$ mesons coming from $B$-meson decay ($F_{B{\rightarrow}J/\psi}$) in $p$+$p$ collisions at $\sqrt{s}=$ 510 GeV. The measurement is performed using the forward silicon vertex detector and central vertex detector at PHENIX, which provide precise tracking and distance-of-closest-approach determinations, enabling the statistical separation of $J/\psi$ due to $B$-meson decays from prompt $J/\psi$. The measured value of $F_{B{\rightarrow}J/\psi}$ is 8.1\%$\pm$2.3\% (stat)$\pm$1.9\% (syst) for $J/\psi$ with transverse momenta $0<p_T<5$ GeV/$c$ and rapidity $1.2<|y|<2.2$. The measured fraction $F_{B{\rightarrow}J/\psi}$ at PHENIX is compared to values measured by other experiments at higher center of mass energies and to fixed-order-next-to-leading-logarithm and color-evaporation-model predictions. The $b\bar{b}$ cross section per unit rapidity ($d\sigma/dy(pp{\rightarrow}b\bar{b})$) extracted from the obtained $F_{B{\rightarrow}J/\psi}$ and the PHENIX inclusive $J/\psi$ cross section measured at 200 GeV scaled with color-evaporation-model calculations, at the mean $B$ hadron rapidity $y={\pm}1.7$ in 510 GeV $p$$+$$p$ collisions, is $3.63^{+1.92}_{-1.70}\mu$b, and it is consistent with the fixed-order-next-to-leading-logarithm calculations.

We investigate electronic excitations in La2-x(Br,Sr)xCuO4 using resonant inelastic x-ray scattering (RIXS) at the oxygen K edge. RIXS spectra of the hole-doped cuprates show clear momentum dependence below 1 eV. The spectral weight exhibits positive dispersion and shifts to higher energy with increasing hole concentration. Theoretical calculation of the dynamical charge structure factor on oxygen orbitals in a three-band Hubbard model is consistent with the experimental observation of the momentum and doping dependence, and therefore the dispersive mode is ascribed to intraband charge excitations which have been observed in electron-doped cuprates.

A measurement of the inclusive production cross section of isolated prompt photons in proton-antiproton collisions at center-of-mass energy $\sqrt{s}$=1.96TeV is presented. The results are obtained using the full Run II data sample collected with the Collider Detector at the Fermilab Tevatron, which corresponds to an integrated luminosity of 9.5fb$^{-1}$. The cross section is measured as a function of photon transverse energy, $E_T^{\gamma}$, in the range 30$ < E_T^{\gamma} <$500GeV and in the pseudorapidity region $|\eta^{\gamma}|<$1.0. The results are compared with predictions from parton-shower Monte Carlo models at leading order in quantum chromodynamics (QCD) and from next-to-leading order perturbative QCD calculations. The latter show good agreement with the measured cross section.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. Correlations of charged hadrons of $0.7<p_T<10$ GeV/$c$ with $\pi^0$ mesons of $4<p_T<15$ GeV/$c$ or isolated direct photons of $7<p_T<15$ GeV/$c$ are used to study nonperturbative effects generated by initial-state partonic transverse momentum and final-state transverse momentum from fragmentation. The nonperturbative behavior is characterized by measuring the out-of-plane transverse momentum component $p_{\rm out}$ perpendicular to the axis of the trigger particle, which is the high-$p_T$ direct photon or $\pi^0$. Nonperturbative evolution effects are extracted from Gaussian fits to the away-side inclusive-charged-hadron yields for different trigger-particle transverse momenta ($p_T^{\rm trig}$). The Gaussian widths and root mean square of $p_{\rm out}$ are reported as a function of the interaction hard scale $p_T^{\rm trig}$ to investigate possible transverse-momentum-dependent evolution differences between the $\pi^0$-h$^\pm$ and direct photon-h$^\pm$ correlations and factorization breaking effects. The widths are found to decrease with $p_T^{\rm trig}$, which indicates that the Collins-Soper-Sterman soft factor is not driving the evolution with the hard scale in nearly back-to-back dihadron and direct photon-hadron production in $p$$+$$p$ collisions. This behavior is in contrast to Drell-Yan and semi-inclusive deep-inelastic scattering measurements.

The fraction of $J/\psi$ mesons which come from B-meson decay, $\textrm{F}_{B{\rightarrow}J/\psi}$, is measured for J/$\psi$ rapidity \mbox{$1.2<|y|<2.2$} and $p_T>0$ in $p$+$p$ and Cu+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV with the PHENIX detector. The extracted fraction is $\textrm{F}_{B{\rightarrow}J/\psi}$ = 0.025 $\pm$ 0.006(stat) $\pm$ 0.010(syst) for $p$+$p$ collisions. For Cu+Au collisions, $\textrm{F}_{B{\rightarrow}J/\psi}$ is 0.094 $\pm$ 0.028(stat) $\pm$ 0.037(syst) in the Au-going direction ($-2.2<y<-1.2$) and 0.089 $\pm$ 0.026(stat) $\pm$ 0.040(syst) in the Cu-going direction ($1.2<y<2.2$). The nuclear modification factor, $R_{\rm CuAu}$, of B mesons in Cu+Au collisions is consistent with binary scaling of measured yields in $p$+$p$ at both forward and backward rapidity.

The PHENIX Collaboration has measured the ratio of the yields of $\psi(2S)$ to $\psi(1S)$ mesons produced in $p$$+$$p$, $p$$+$Al, $p$$+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV over the forward and backward rapidity intervals $1.2<|y|<2.2$. We find that the ratio in $p$$+$$p$ collisions is consistent with measurements at other collision energies. In collisions with nuclei, we find that in the forward ($p$-going or $^{3}$He-going) direction, the relative yield of $\psi(2S)$ mesons to $\psi(1S)$ mesons is consistent with the value measured in \pp collisions. However, in the backward (nucleus-going) direction, the $\psi(2S)$ is preferentially suppressed by a factor of $\sim$2. This suppression is attributed in some models to breakup of the weakly-bound $\psi(2S)$ through final state interactions with comoving particles, which have a higher density in the nucleus-going direction. These breakup effects may compete with color screening in a deconfined quark-gluon plasma to produce sequential suppression of excited quarkonia states.

We searched for the $CP$-violating rare decay of neutral kaon, $K_{L} \to \pi^0 \nu \overline{\nu}$, in data from the first 100 hours of physics running in 2013 of the J-PARC KOTO experiment. One candidate event was observed while $0.34\pm0.16$ background events were expected. We set an upper limit of $5.1\times10^{-8}$ for the branching fraction at the 90\% confidence level (C.L.). An upper limit of $3.7\times10^{-8}$ at the 90\% C.L. for the $K_{L} \to \pi^{0} X^{0}$decay was also set for the first time, where $X^{0}$ is an invisible particle with a mass of 135 MeV/$c^{2}$.

We describe the new version (v2.73y) of the code HFODD which solves the nuclear Skyrme Hartree-Fock or Skyrme Hartree-Fock-Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we have implemented the following new features: (i) full proton-neutron mixing in the particle-hole channel for Skyrme functionals, (ii) the Gogny force in both particle-hole and particle-particle channels, (iii) linear multi-constraint method at finite temperature, (iv) fission toolkit including the constraint on the number of particles in the neck between two fragments, calculation of the interaction energy between fragments, and calculation of the nuclear and Coulomb energy of each fragment, (v) the new version 200d of the code HFBTHO, together with an enhanced interface between HFBTHO and HFODD, (vi) parallel capabilities, significantly extended by adding several restart options for large-scale jobs, (vii) the Lipkin translational energy correction method with pairing, (viii) higher-order Lipkin particle-number corrections, (ix) interface to a program plotting single-particle energies or Routhians, (x) strong-force isospin-symmetry-breaking terms, and (xi) the Augmented Lagrangian Method for calculations with 3D constraints on angular momentum and isospin. Finally, an important bug related to the calculation of the entropy at finite temperature and several other little significant errors of the previous published version were corrected.

We report on a measurement of the $D^{+}$-meson production cross section as a function of transverse momentum ($p_T$) in proton-antiproton ($p\bar{p}$) collisions at 1.96 TeV center-of-mass energy, using the full data set collected by the Collider Detector at Fermilab in Tevatron Run II and corresponding to 10 fb$^{-1}$ of integrated luminosity. We use $D^{+} \to K^-\pi^+\pi^+$ decays fully reconstructed in the central rapidity region $|y|<1$ with transverse momentum down to 1.5 GeV/$c$, a range previously unexplored in $p\bar{p}$ collisions. Inelastic $p\bar{p}$-scattering events are selected online using minimally-biasing requirements followed by an optimized offline selection. The $K^-\pi^+\pi^+$ mass distribution is used to identify the $D^+$ signal, and the $D^+$ transverse impact-parameter distribution is used to separate prompt production, occurring directly in the hard scattering process, from secondary production from $b$-hadron decays. We obtain a prompt $D^+$ signal of 2950 candidates corresponding to a total cross section $\sigma(D^+, 1.5 < p_T < 14.5~\mbox{GeV/}c, |y|<1) = 71.9 \pm 6.8 (\mbox{stat}) \pm 9.3 (\mbox{syst})~\mu$b. While the measured cross sections are consistent with theoretical estimates in each $p_T$ bin, the shape of the observed $p_T$ spectrum is softer than the expectation from quantum chromodynamics. The results are unique in $p\bar{p}$ collisions and can improve the shape and uncertainties of future predictions.