Using the entire BABAR gamma(4S) dataset, the first two-dimensional unbinned angular analysis of the semileptonic decay (B) over bar -> Dl(-)(v) over bar (l) is performed, employing hadronic reconstruction of the tag-side B meson from gamma(4S) -> B (B) over bar. Here, l denotes the light charged leptons e and mu. Anovel data-driven signal-background separation procedure with minimal dependence on simulation is developed. This procedure preserves all multidimensional correlations present in the data. The expected sin(2) theta(l) dependence of the differential decay rate in the Standard Model is demonstrated, where theta(l) is the lepton helicity angle. Including input from the latest lattice QCD calculations and previously available experimental data, the underlying form factors are extracted using both model-independent (BGL) and dependent (CLN) methods. Comparisons with lattice calculations show flavor SU(3) symmetry to be a good approximation in the B-(s) -> D-(s) sector. Using the BGL results, the CKM matrix element vertical bar V-cb vertical bar = (41.09 +/- 1.16) x 10(-3) and the Standard Model prediction of the lepton-flavor universality violation variable R(D) = 0.300 +/- 0.004, are extracted. The value of vertical bar V-cb vertical bar from (B) over bar -> Dl(-)(v) over bar (l) tends to be higher than that extracted using (B) over bar -> Dl(-)(v) over bar (l). The Standard Model R(D) calculation is at a 1.97 sigma tension with the latest HFLAV experimental average.

We study the processes e+e- -> K+K-a degrees a degrees a degrees y, K degrees SK +/- a -/+ a degrees a degrees y, and K degrees SK +/- a -/+ a+a-y in which an energetic photon is radiated from the initial state. The data were collected with the BABAR detector at the SLAC National Accelerator Laboratory. About 1200, 2600, and 6000 events, respectively, are selected from a data sample corresponding to an integrated luminosity of 469 fb-1. The invariant mass of the hadronic final state defines the effective e+e- center-of-mass energy. The center-of-mass energies range from threshold to 4.5 GeV. From the mass spectra, the first ever measurements of the e+e- -> K+K-a degrees a degrees a degrees, e+e- -> K degrees SK +/- a -/+ a degrees a degrees, and e+e- -> K degrees SK +/- a -/+ a+a- cross sections are performed. The contributions from the intermediate states that include eta, phi, rho, K*(892), and other resonances are presented. We observe the J=psi and psi(2S) in most of these final states and measure the corresponding branching fractions, many of them for the first time.

A new mechanism has been proposed to simultaneously explain the presence of dark matter and the matter-antimatter asymmetry in the universe. This scenario predicts exotic $B$ meson decays into a baryon and a dark sector anti-baryon ($\psi_D$) with branching fractions accessible at $B$ factories. We present a search for $B \rightarrow \Lambda \psi_D$ decays using data collected by the $BABAR$ experiment at SLAC. This reaction is identified by fully reconstructing the accompanying $B$ meson and requiring the presence of a single $\Lambda$ baryon in the remaining particles. No significant signal is observed, and bounds on the $B \rightarrow \Lambda \psi_D$ branching fraction are derived in the range $0.13 - 5.2\times 10^{-5}$ for $1.0 < m_{\psi_D} < 4.2$ GeV/$c^{2}$. These results set strong constraints on the parameter space allowed by the theory.

We study the processes ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{K}^{+}{K}^{\ensuremath{-}}{\ensuremath{\pi}}^{0}{\ensuremath{\pi}}^{0}{\ensuremath{\pi}}^{0}\ensuremath{\gamma}$, ${K}_{S}^{0}{K}^{\ifmmode\pm\else\textpm\fi{}}{\ensuremath{\pi}}^{\ensuremath{\mp}}{\ensuremath{\pi}}^{0}{\ensuremath{\pi}}^{0}\ensuremath{\gamma}$, and ${K}_{S}^{0}{K}^{\ifmmode\pm\else\textpm\fi{}}{\ensuremath{\pi}}^{\ensuremath{\mp}}{\ensuremath{\pi}}^{+}{\ensuremath{\pi}}^{\ensuremath{-}}\ensuremath{\gamma}$ in which an energetic photon is radiated from the initial state. The data were collected with the BABAR detector at the SLAC National Accelerator Laboratory. About 1200, 2600, and 6000 events, respectively, are selected from a data sample corresponding to an integrated luminosity of $469\text{ }\text{ }{\mathrm{fb}}^{\ensuremath{-}1}$. The invariant mass of the hadronic final state defines the effective ${e}^{+}{e}^{\ensuremath{-}}$ center-of-mass energy. The center-of-mass energies range from threshold to 4.5 GeV. From the mass spectra, the first ever measurements of the ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{K}^{+}{K}^{\ensuremath{-}}{\ensuremath{\pi}}^{0}{\ensuremath{\pi}}^{0}{\ensuremath{\pi}}^{0}$, ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{K}_{S}^{0}{K}^{\ifmmode\pm\else\textpm\fi{}}{\ensuremath{\pi}}^{\ensuremath{\mp}}{\ensuremath{\pi}}^{0}{\ensuremath{\pi}}^{0}$, and ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{K}_{S}^{0}{K}^{\ifmmode\pm\else\textpm\fi{}}{\ensuremath{\pi}}^{\ensuremath{\mp}}{\ensuremath{\pi}}^{+}{\ensuremath{\pi}}^{\ensuremath{-}}$ cross sections are performed. The contributions from the intermediate states that include $\ensuremath{\eta}$, $\ensuremath{\phi}$, $\ensuremath{\rho}$, ${K}^{*}(892)$, and other resonances are presented. We observe the $J/\ensuremath{\psi}$ and $\ensuremath{\psi}(2S)$ in most of these final states and measure the corresponding branching fractions, many of them for the first time.

Using the entire BaBar Υ(4S) data set, the first two-dimensional unbinned angular analysis of the semileptonic decay B→ D ℓ^- ν_ℓ is performed, employing hadronic reconstruction of the tag-side B meson from Υ(4S)→ BB. Here, ℓ denotes the light charged leptons e and μ. A novel data-driven signal-background separation procedure with minimal dependence on simulation is developed. This procedure preserves all multi-dimensional correlations present in the data. The expected sin^2θ_ℓ dependence of the differential decay rate in the Standard Model is demonstrated, where θ_ℓ is the lepton helicity angle. Including input from the latest lattice QCD calculations and previously available experimental data, the underlying form factors are extracted using both model-independent (BGL) and dependent (CLN) methods. Comparisons with lattice calculations show flavor SU(3) symmetry to be a good approximation in the B_(s)→ D_(s) sector. Using the BGL results, the CKM matrix element |V_cb|=(41.09± 1.16)× 10^-3 and the Standard Model prediction of the lepton-flavor universality violation variable ℛ(D)=0.300± 0.004, are extracted. The value of |V_cb| from B→ D ℓ^- ν_ℓ tends to be higher than that extracted using B→ D ℓ^- ν_ℓ. The Standard Model ℛ(D) calculation is at a 1.97σ tension with the latest HFLAV experimental average.

A dedicated measurement of additional radiation in e+e- -> mu+mu- gamma and e+e- -> pi+ pi- gamma initial-state radiation events is presented using the full BABAR data sample. For the first time results are presented at next -to-and next-to-next-to-leading order, with one and two additional photons, respectively, for radiation from the initial and final states. Comparison with predictions from PHOKHARA and AFKQED Monte Carlo generators is performed, revealing discrepancies in the one-photon rates and angular distributions for the former. This disagreement has a negligible effect on the BABAR measurement of the e+e- -> pi+ pi-(gamma) cross section, but could affect other measurements significantly. This study sheds a new light on the longstanding discrepancy in this channel that affects the theoretical prediction of hadronic vacuum polarization contributions to the muon magnetic moment anomaly.

A new dark sector antibaryon, denoted ψ_{D}, could be produced in decays of B mesons. This Letter presents a search for B^{+}→ψ_{D}+p (and the charge conjugate) decays in e^{+}e^{-} annihilations at 10.58 GeV, using data collected in the BABAR experiment. Data corresponding to an integrated luminosity of 398  fb^{-1} are analyzed. No evidence for a signal is observed. Branching fraction upper limits in the range from 10^{-7}-10^{-5} are obtained at 90% confidence level for masses of 1.0<m_{ψ_{D}}<4.3  GeV/c^{2}. The result is also reinterpreted to provide the first limits on a supersymmetric model with R-parity violation and a light neutralino.

A dedicated measurement of additional radiation in $e^+e^-\to\mu^+\mu^-\gamma$ and $e^+e^-\to\pi^+\pi^-\gamma$ initial-state-radiation events is presented using the full BABAR data sample. For the first time results are presented at next-to- and next-to-next-to-leading order, with one and two additional photons, respectively, for radiation from the initial and final states. Comparison with predictions from Phokhara and AfkQed Monte Carlo generators is performed, revealing discrepancies in the one-photon rates and angular distributions for the former. This disagreement has a negligible effect on the BABAR measurement of the $e^+e^-\to\pi^+\pi^-(\gamma)$ cross section, but could affect other measurements significantly. This study sheds a new light on the longstanding discrepancy in this channel that affects the theoretical prediction of hadronic vacuum polarization contributions to the muon magnetic moment anomaly.

We present a new technique, "tau polarimetry," for measuring the longitudinal beam polarization present in an e(+)e(-) collider through the analysis of e(+)e(-) -> tau(+)tau(-) events. By exploiting the sensitivity of tau decay kinematics to the longitudinal polarization of the beams, we demonstrate that the longitudinal polarization can be measured with a 3 per mil systematic uncertainty at the interaction point using a technique that is independent of spin and beam transport modeling. Using 424.2 +/- 1.8 fb(-1) of BABAR data at root s = 10.58 GeV, the average longitudinal polarization of the PEP-II e(+)e(-) collider has been measured to be < P > = 0.0035 +/- 0.0024(stat) +/- 0.0029(sys). The systematic uncertainty studies are described in detail, which can serve as a guide for future applications of tau polarimetry. A proposed e(-) beam longitudinal polarization upgrade to the SuperKEKB e(+)e(-) collider would benefit from this technique.

We present a new technique, ``tau polarimetry,'' for measuring the longitudinal beam polarization present in an ${e}^{+}{e}^{\ensuremath{-}}$ collider through the analysis of ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{\ensuremath{\tau}}^{+}{\ensuremath{\tau}}^{\ensuremath{-}}$ events. By exploiting the sensitivity of $\ensuremath{\tau}$ decay kinematics to the longitudinal polarization of the beams, we demonstrate that the longitudinal polarization can be measured with a 3 per mil systematic uncertainty at the interaction point using a technique that is independent of spin and beam transport modeling. Using $424.2\ifmmode\pm\else\textpm\fi{}1.8\text{ }\text{ }{\mathrm{fb}}^{\ensuremath{-}1}$ of BABAR data at $\sqrt{s}=10.58\text{ }\text{ }\mathrm{GeV}$, the average longitudinal polarization of the PEP-II ${e}^{+}{e}^{\ensuremath{-}}$ collider has been measured to be $⟨P⟩=0.0035\ifmmode\pm\else\textpm\fi{}0.002{4}_{\mathrm{stat}}\ifmmode\pm\else\textpm\fi{}0.002{9}_{\mathrm{sys}}$. The systematic uncertainty studies are described in detail, which can serve as a guide for future applications of tau polarimetry. A proposed ${e}^{\ensuremath{-}}$ beam longitudinal polarization upgrade to the SuperKEKB ${e}^{+}{e}^{\ensuremath{-}}$ collider would benefit from this technique.

This article presents a model-independent search for an additional, mostly sterile, heavy neutral lepton (HNL), that is capable of mixing with the Standard Model tau neutrino with a mixing strength of jU tau 4j2, corresponding to the absolute square of the extended Pontecorvo-Maki-Nakagawa-Sakata matrix element. Data from the BABAR experiment, with a total integrated luminosity of 424 fb-1, are analyzed using a kinematic approach that makes no assumptions on the model behind the origins of the HNL, its lifetime or decay modes. No significant signal is found. Upper limits on jU tau 4j2 at the 95% confidence level, depend on the HNL mass hypothesis and vary from 2.31 x 10-2 to 5.04 x 10-6 (with all uncertainties considered), across the mass range 100 < m4 < 1300 MeV/c2; the more stringent limits being placed at higher masses.

We report on the first search for electron-muon lepton flavor violation (LFV) in the decay of a b quark and b antiquark bound state. We look for the LFV decay ϒ(3S)→e^{±}μ^{∓} in a sample of 118 million ϒ(3S) mesons from 27 fb^{-1} of data collected with the BABAR detector at the SLAC PEP-II e^{+}e^{-} collider operating with a 10.36 GeV center-of-mass energy. No evidence for a signal is found, and we set a limit on the branching fraction B[ϒ(3S)→e^{±}μ^{∓}]<3.6×10^{-7} at 90% C. L. This result can be interpreted as a limit Λ_{NP}/g_{NP}^{2}>80 TeV on the energy scale Λ_{NP} divided by the coupling-squared g_{NP}^{2} of relevant new physics (NP).

Axionlike particles (ALPs) are predicted in many extensions of the standard model, and their masses can naturally be well below the electroweak scale. In the presence of couplings to electroweak bosons, these particles could be emitted in flavor-changing B meson decays. We report herein a search for an ALP, a, in the reaction B^{±}→K^{±}a, a→γγ using data collected by the BABAR experiment at SLAC. No significant signal is observed, and 90% confidence level upper limits on the ALP coupling to electroweak bosons are derived as a function of ALP mass, improving current constraints by several orders of magnitude in the range 0.175  GeV<m_{a}<4.78  GeV.

Collider searches for dark sectors, new particles interacting only feebly with ordinary matter, have largely focused on identifying signatures of new mediators, leaving much of dark sector structures unexplored. In particular, the existence of dark matter bound states (darkonia) remains to be investigated. This possibility could arise in a simple model in which a dark photon (A^{'}) is light enough to generate an attractive force between dark fermions. We report herein a search for a J^{PC}=1^{--} darkonium state, the ϒ_{D}, produced in the reaction e^{+}e^{-}→γϒ_{D}, ϒ_{D}→A^{'}A^{'}A^{'}, where the dark photons subsequently decay into pairs of leptons or pions, using 514  fb^{-1} of data collected with the BABAR detector. No significant signal is observed, and we set bounds on the γ-A^{'} kinetic mixing as a function of the dark sector coupling constant for 0.001<m_{A^{'}}<3.16  GeV and 0.05<m_{ϒ_{D}}<9.5  GeV.

We study the processes gamma gamma -> eta(c )-> eta'K+K-, eta'pi(+)pi(-), and eta pi(+)pi(-) using a data sample of 519 fb(-1) recorded with the BABAR detector operating at the SLAC PEP-II asymmetric-energy e(+)e(-) collider at center-of-mass energies at and near the Upsilon(nS) (n = 2, 3, 4) resonances. This is the first observation of the decay eta(c)-> eta'K+K- and we measure the branching fraction Gamma (eta(c) -> eta' K+K- )/(Gamma(eta(c) -> eta'pi(+)pi(-)) = 0.644 +/- 0.039(stat) +/- 0.032(sys). Significant interference is observed between gamma gamma ->eta(c)->eta pi(+)pi(-) and the nonresonant two-photon process gamma gamma -> eta pi(+)pi(-). A Dalitz plot analysis is performed of eta(c) decays to eta'K+K-, eta'pi(+)pi(-) and eta pi(+)pi(-). Combined with our previous analysis of eta(c) -> K (K) over bar pi, we measure the K-0*(1430) parameters and the ratio between its eta'K and pi K couplings. The decay eta(c) -> eta'pi(+)pi(-) is dominated by the f(0)(2100) resonance, also observed in J/psi radiative decays. A new a(0) (1700) -> eta pi resonance is observed in the eta(c) -> eta pi(+)pi(-) channel. We also compare eta(c) decays to eta and eta' final states in association with scalar mesons as they relate to the identification of the scalar glueball.

We study the processes e(+)e(-) -> pi+ pi-pi(0)pi(0)pi(0)pi(0)gamma and pi(+) pi(-) pi(0)pi(0)pi eta gamma in which an energetic photon is radiated from the initial state. The data were collected with the BABAR detector at the SLAC National Accelerator Laboratory. About 7300 and 870 events, respectively, are selected from a data sample corresponding to an integrated luminosity of 469 fb(-1). The invariant mass of the hadronic final state defines the effective e(+)e(-) center-of-mass energy. The center-of-mass energies range from threshold to 4.5 GeV. From the mass spectra, the first ever measurements of the e(+)e(-) -> pi(+)pi(-) pi(0)pi(0)pi(0)pi(0) and the e(+)e(-) -> pi(+)pi(-) pi(0)pi(0)pi(0)eta cross sections are performed. The contributions from omega pi(0)pi(0)pi(0), eta pi(+)pi(-)pi(0), omega eta, and other intermediate states are presented. We observe the J=psi and psi(2S) in most of these final states and measure the corresponding branching fractions, many of them for the first time.

We study the processes $e^+e^-\to 2(\pi^+\pi^-)\pi^0\pi^0\pi^0\gamma$ and $2(\pi^+\pi^-)\pi^0\pi^0\eta\gamma$ in which an energetic photon is radiated from the initial state. The data were collected with the \babar~ detector at SLAC. About 14\,000 and 4700 events, respectively, are selected from a data sample corresponding to an integrated luminosity of 469 fb^{-1}. The invariant mass of the hadronic final state defines the effective \epem center-of-mass energy. The center-of-mass energies range from threshold to 4.5 GeV. From the mass spectra, the first ever measurements of the $e^+e^-\to 2(\pi^+\pi^-)\pi^0\pi^0\pi^0$ and the $e^+e^-\to2(\pi^+\pi^-)\pi^0\pi^0\eta$ cross sections are performed. The contributions from $\omega\pi^+\pi^-\pi^0\pi^0$, $\eta2(\pi^+\pi^-)$, and other intermediate states are presented. We observe the $J/\psi$ and $\psi(2S)$ in most of these final states and measure the corresponding branching fractions, many of them for the first time.

We study the processes e(+)e(-) -> 2(pi(+)pi(-))pi(0)pi(0)pi(0)gamma and 2(pi(+)pi(-))pi(0)pi(0)eta gamma in which an energetic photon is radiated from the initial state. The data were collected with the BABAR detector at SLAC. About 14 000 and 4700 events, respectively, are selected from a data sample corresponding to an integrated luminosity of 469 fb(-1). The invariant mass of the hadronic final state defines the effective e(+)e(-) center-of-mass energy. The center-of-mass energies range from threshold to 4.5 GeV. From the mass spectra, the first ever measurements of the e(+)e(-) -> 2(pi(+)pi(-))pi(0)pi(0)pi(0) and the e(+)e(-) -> 2(pi(+)pi(-))pi(0)pi(0)eta cross sections are performed. The contributions from omega pi(+)pi(-)pi(0)pi(0), eta 2(pi(+)pi(-)), and other intermediate states are presented. We observe the J/psi and psi(2S) in most of these final states and measure the corresponding branching fractions, many of them for the first time.

J. P. Lees, V. Poireau, V. Tisserand, E. Grauges, A. Palano, G. Eigen, D. N. Brown, Yu. G. Kolomensky, M. Fritsch, H. Koch, T. Schroeder, R. Cheaib, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, E. A. Kozyrev, E. A. Kravchenko, A. P. Onuchin, S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, K. Yu. Todyshev, A. J. Lankford, B. Dey, J. W. Gary, O. Long, A. M. Eisner, W. S. Lockman, W. Panduro Vazquez, D. S. Chao, C. H. Cheng, B. Echenard, K. T. Flood, D. G. Hitlin, J. Kim, Y. Li, D. X. Lin, T. S. Miyashita, P. Ongmongkolkul, J. Oyang, F. C. Porter, M. Röhrken, Z. Huard, B. T. Meadows, B. G. Pushpawela, M. D. Sokoloff, L. Sun, ∗ J. G. Smith, S. R. Wagner, D. Bernard, M. Verderi, D. Bettoni, C. Bozzi, R. Calabrese, G. Cibinetto, E. Fioravanti, I. Garzia, E. Luppi, V. Santoro, A. Calcaterra, R. de Sangro, G. Finocchiaro, S. Martellotti, P. Patteri, I. M. Peruzzi, M. Piccolo, M. Rotondo, A. Zallo, S. Passaggio, C. Patrignani, † B. J. Shuve, H. M. Lacker, B. Bhuyan, U. Mallik, C. Chen, J. Cochran, S. Prell, A. V. Gritsan, N. Arnaud, M. Davier, F. Le Diberder, A. M. Lutz, G. Wormser, D. J. Lange, D. M. Wright, J. P. Coleman, E. Gabathuler, ‡ D. E. Hutchcroft, D. J. Payne, C. Touramanis, A. J. Bevan, F. Di Lodovico, § R. Sacco, G. Cowan, Sw. Banerjee, D. N. Brown, C. L. Davis, A. G. Denig, W. Gradl, K. Griessinger, A. Hafner, K. R. Schubert, R. J. Barlow, ¶ G. D. Lafferty, R. Cenci, A. Jawahery, D. A. Roberts, R. Cowan, S. H. Robertson, R. M. Seddon, N. Neri, F. Palombo, L. Cremaldi, R. Godang, ∗∗ D. J. Summers, P. Taras, G. De Nardo, C. Sciacca, G. Raven, C. P. Jessop, J. M. LoSecco, K. Honscheid, R. Kass, A. Gaz, M. Margoni, M. Posocco, G. Simi, F. Simonetto, R. Stroili, S. Akar, E. Ben-Haim, M. Bomben, G. R. Bonneaud, G. Calderini, J. Chauveau, G. Marchiori, J. Ocariz, M. Biasini, E. Manoni, A. Rossi, G. Batignani, S. Bettarini, M. Carpinelli, †† G. Casarosa, M. Chrzaszcz, F. Forti, M. A. Giorgi, A. Lusiani, B. Oberhof, E. Paoloni, M. Rama, G. Rizzo, J. J. Walsh, L. Zani, A. J. S. Smith, F. Anulli, R. Faccini, F. Ferrarotto, F. Ferroni, ‡‡ A. Pilloni, G. Piredda, ‡ C. Bünger, S. Dittrich, O. Grünberg, M. Heß, T. Leddig, C. Voß, R. Waldi, T. Adye, F. F. Wilson, S. Emery, G. Vasseur, D. Aston, C. Cartaro, M. R. Convery, J. Dorfan, W. Dunwoodie, M. Ebert, R. C. Field, B. G. Fulsom, M. T. Graham, C. Hast, W. R. Innes, ‡ P. Kim, D. W. G. S. Leith, ‡ S. Luitz, D. B. MacFarlane, D. R. Muller, H. Neal, B. N. Ratcliff, A. Roodman, M. K. Sullivan, J. Va’vra, W. J. Wisniewski, M. V. Purohit, J. R. Wilson, A. Randle-Conde, S. J. Sekula, H. Ahmed, M. Bellis, P. R. Burchat, E. M. T. Puccio, M. S. Alam, J. A. Ernst, R. Gorodeisky, N. Guttman, D. R. Peimer, A. Soffer, S. M. Spanier, J. L. Ritchie, R. F. Schwitters, J. M. Izen, X. C. Lou, F. Bianchi, F. De Mori, A. Filippi, D. Gamba, L. Lanceri, L. Vitale, F. Martinez-Vidal, A. Oyanguren, J. Albert, A. Beaulieu, F. U. Bernlochner, G. J. King, R. Kowalewski, T. Lueck, I. M. Nugent, J. M. Roney, A. Sibidanov, R. J. Sobie, N. Tasneem, T. J. Gershon, P. F. Harrison, T. E. Latham, R. Prepost, and S. L. Wu

Many scenarios of physics beyond the standard model predict the existence of new gauge singlets, which might be substantially lighter than the weak scale. The experimental constraints on additional scalars with masses in the MeV to GeV range could be significantly weakened if they interact predominantly with leptons rather than quarks. At an e^{+}e^{-} collider, such a leptophilic scalar (ϕ_{L}) would be produced predominantly through radiation from a τ lepton. We report herein a search for e^{+}e^{-}→τ^{+}τ^{-}ϕ_{L}, ϕ_{L}→ℓ^{+}ℓ^{-} (ℓ=e, μ) using data collected by the BABAR experiment at SLAC. No significant signal is observed, and we set limits on the ϕ_{L} coupling to leptons in the range 0.04<m_{ϕ_{L}}<7.0  GeV. These bounds significantly improve upon the current constraints, excluding almost entirely the parameter space favored by the observed discrepancy in the muon anomalous magnetic moment below 4 GeV at 90% confidence level.