Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: since its discovery in 1974 @xcite , the charmonium system has become the prototypical ` hydrogen atom ' of meson spectroscopy @xcite . the experimentally clear spectrum of relatively narrow states below the open - charm dd threshold of 3.73  gev can be identified with the 1s , 1p and 2s levels predicted by potential models , which incorporate a color coulomb term at short distances and a linear scalar confining term at large distances . spin - dependent interquark forces are evident in the splittings of states within these multiplets , and the observed splittings are consistent with the predictions of a one gluon exchange ( oge ) breit - fermi hamiltonian , combined with a linear scalar confining interaction . discussions of the theoretical importance and experimental status of heavy quarkonium , including recent experimental results for charmonium , have been given by quigg @xcite , galik @xcite , the cern quarkonium working group @xcite , seth @xcite and swarnicki @xcite . recently there has been a resurgence of interest in charmonium , due to the realization that b factories can contribute to the study of the missing states @xcite , and to high - statistics experiments at bes @xcite and cleo @xcite and the planned gsi @xmath3 facility @xcite . the possibility of contributions from b factories was dramatically illustrated by the recent discovery of the long missing @xmath4 @xmath5 state by the belle collaboration @xcite , which has since been confirmed by babar @xcite , and has also been observed by cleo in @xmath6 collisions @xcite . additional interest in spectroscopy has followed the discovery of the remarkable x@xmath7 by belle @xcite and cdf @xcite in b decays to @xmath8 ; assuming that this is a real resonance rather than a threshold effect , the x@xmath7 is presumably either a dd@xmath9 charmed meson molecule @xcite or a narrow @xmath10 d - wave state @xcite . very recent observations of the x@xmath7 in @xmath11 and @xmath12 by belle support a @xmath13 dd@xmath9 molecule assignment @xcite . there has also been experimental activity in the spin - singlet p - wave sector , with recent reports of the observation of the elusive @xmath14 @xmath15 state by cleo @xcite . finally , the surprisingly large cross sections for double charmonium production in @xmath16 reported by belle @xcite suggest that it may be possible to study c @xmath17 @xmath0 states in @xmath16 without using the higher - order @xmath18 two - photon annihilation process . one open topic of great current interest in spectroscopy is the search for the @xmath19 and @xmath20 states , which are expected to be quite narrow due to the absence of open - charm decay modes . a second topic is the lorentz nature of confinement ; in pure models this is tested by the multiplet splittings of orbitally - excited states . for example , with pure scalar confinement as is normally assumed there is no spin - spin hyperfine interaction at o@xmath21 , so the masses of spin - singlets ( such as the @xmath22 @xmath15 ) are degenerate with the corresponding triplet c.o.g . ( center of gravity ) ; here this is the @xmath1p@xmath23 c.o.g . , at 3525  mev . in the original cornell model @xcite it was assumed that confinement acts as the time component of a lorentz vector , which lifts the degeneracy of the @xmath15 and the @xmath1p@xmath23 c.o.g . another possibility is that confinement may be a more complicated mix of scalar and timelike vector @xcite . of course these simple potential model considerations may be complicated by mass shifts due to other effects , such as couplings to open - flavor channels @xcite . a third topic is the search for exotica such as hybrids ; the level of mixing between conventional quarkonium and hybrid basis states falls rapidly with increasing quark mass , which suggests that nonexotic hybrids may be more easily distinguished from conventional quarkonia in charmonium than in the light quark sectors . since lgt ( lattice gauge theory ) predicts that the lightest hybrids lie near 4.4  gev @xcite , there is a strong incentive to establish the  background " spectrum of conventional states up to and somewhat beyond this mass . a final topic of current interest is the importance of mixing between quark model @xmath24 basis states and two - meson continua , which has been cited as a possible reason for the low masses of the recently discovered d@xmath25 states @xcite . the effects of  unquenching the quark model " by including meson loops can presumably be studied effectively in the system , in which the experimental spectrum of states is relatively unambiguous . the success of the @xmath24 quark model is surprising , in view of the probable importance of corrections to the valence approximation ; the range of validity of the naive  quenched " @xmath24 quark model is an interesting and open question @xcite . motivated by this revived interest in spectroscopy , we have carried out a theoretical study of the expected properties of charmonium states , notably the poorly understood higher - mass levels above dd threshold . two variants of potential models are used in this study , a conventional nonrelativistic model based on the schrdinger equation with a coulomb plus linear potential , and the godfrey - isgur relativized potential model . we give results for all states in the multiplets 1 - 4s , 1 - 3p , 1 - 2d , 1 - 2f and 1 g , comprising 40 resonances in total . predictions are given for quantities which are likely to be of greatest experimental interest , which are the spectrum of states , e1 ( and some m1 ) electromagnetic transition rates , and strong partial and total widths of states above open - charm threshold . similar results for many of the electromagnetic transition rates have recently been reported by ebert _ the @xmath26 leptonic and two - photon widths are not discussed in detail here , as they have been considered extensively elsewhere ; see for example @xcite and @xcite and references cited therein . as a minimal model of the charmonium system we use a nonrelativistic potential model , with wavefunctions determined by the schrdinger equation with a conventional quarkonium potential . we use the standard color coulomb plus linear scalar form , and also include a gaussian - smeared contact hyperfine interaction in the zeroth - order potential . the central potential is v_0^(c|c)(r ) = - + br + _ ( r ) _ c _ |c where @xmath27 . the four parameters @xmath28 are determined by fitting the spectrum . the spin - spin contact hyperfine interaction is one of the spin - dependent terms predicted by one gluon exchange ( oge ) forces . the contact form @xmath29 is actually an artifact of an o@xmath30 expansion of the t - matrix @xcite , so replacing it by an interaction with a range @xmath31 comparable to @xmath32 is not an unwarranted modification . we treat the remaining spin - dependent terms as mass shifts using leading - order perturbation theory . these are the oge spin - orbit and tensor and a longer - ranged inverted spin - orbit term , which arises from the assumed lorentz scalar confinement . these are explicitly v_spin - dep = . the spin - orbit operator is diagonal in a @xmath33 basis , with the matrix elements @xmath34/2 $ ] . the tensor operator t has nonvanishing diagonal matrix elements only between @xmath35 spin - triplet states , which are    ^3l _ | | ^3l _ =  . for experimental input we use the masses of the 11 reasonably well established @xmath0 states , which are given in table  [ table_spectrum ] ( rounded to 1  mev ) . the parameters that follow from fitting these masses are @xmath36 . given these values , we can predict the masses and matrix elements of the currently unknown states ; table  [ table_spectrum ] and fig . [ spectrum ] show the predicted spectrum . the godfrey - isgur model is a `` relativized '' extension of the nonrelativistic model of the previous section . this model assumes a relativistic dispersion relation for the quark kinetic energy , a qcd - motivated running coupling @xmath37 , a flavor - dependent potential smearing parameter @xmath38 , and replaces factors of quark mass with quark kinetic energy . details of the model and the method of solution may be found in ref.@xcite . the hamiltonian consists of a relativistic kinetic term and a generalized quark - antiquark potential , = h_0 + _ q|q ( , ) where _ 0 = +  . just as in the nonrelativistic model , the quark - antiquark potential v@xmath39 assumed here incorporates the lorentz vector one gluon exchange interaction at short distances and a lorentz scalar linear confining interaction . to first order in @xmath40 , v@xmath41 reduces to the standard non - relativistic result given by eqns.(1 ) and ( 2 ) ( with @xmath42 replaced by a running coupling constant , @xmath37 ) . the full set of model parameters is given in ref.@xcite . note that the string tension and quark mass ( @xmath43 gev@xmath44 and @xmath45 gev ) are significantly larger than the values used in our nonrelativistic model . one important aspect of this model is that it gives reasonably accurate results for the spectrum and matrix elements of quarkonia of all @xmath46 quark flavors , whereas the nonrelativistic model of the previous section is only fitted to the system . the spectra predicted by the nr and gi models ( table  [ table_spectrum ] and fig .  [ spectrum ] ) are quite similar for s- and p - wave states , largely because of the constraints provided by the experimental candidates for these multiplets . we note in passing that these potential model results are very similar to the most recent predictions of the charmonium spectrum from lgt @xcite . at higher l we have only the l=2 @xmath47d@xmath48 and @xmath49d@xmath48 states @xmath50 and @xmath51 to constrain the models , and the predicted mean d - wave multiplet masses differ by _ ca . _ 50  mev . for @xmath52 the absence of experimental states allows a relatively large scatter of predicted mean masses , which differ by as much as @xmath53  mev in the 1 g multiplet . ( the splittings within higher l multiplets in contrast are rather similar . ) the mean multiplet masses predicted by the two models differ largely because of the values assumed for the string tension @xmath54 , which is @xmath55  gev@xmath44 in the gi model but is a rather smaller @xmath56  gev@xmath44 in the nr model . identification of any l=3 or l=4 state would be very useful as a constraint on the spectrum of higher - l states generally . this is unfortunately a difficult problem , since these states are not easily accessible experimentally . as we shall see , one possibility is to produce the @xmath1f@xmath57 @xmath58 state , which is formed in e1 radiative transitions from the @xmath51 and decays dominantly to dd and dd@xmath9 . ( for simplicity in this paper we abbreviate the final state d@xmath59 as  dd " , the state d@xmath60d@xmath9 as  dd@xmath9 " , and so forth . )    1.0 cm    ) . the solid lines are experiment , and the broken lines are theory ( nr model left , gi right ) . spin triplet levels are dashed , and spin singlets are dotted . the dd open - charm threshold at 3.73 gev is also shown.,width=321 ] radiative transitions of higher - mass charmonium states are of interest largely because they provide one of the few pathways between states with different quantum numbers . since typical e1 radiative partial widths of charmonia are 10s to 100s of kev , corresponding to significant branching fractions of @xmath61 to @xmath62 , large event samples of radially excited s - wave states produced in @xmath16 annihilation could be used to identify radially excited p - wave states , which are not otherwise easily produced . similarly , the e1 radiative transition of the nominally 2@xmath1d@xmath48 @xmath51 can be used to produce an f - wave state ; this multiplet would likely be very difficult to reach using other mechanisms . we evaluate these e1 radiative partial widths using    _ e1 ( n ^2+1l _  ^2+1l_ + ) = c_fi _    e_c^2 |_f | r | _ i |^2 e_^3  ,    where @xmath63 is the @xmath64-quark charge in units of @xmath65 , @xmath66 is the fine - structure constant , e@xmath67 is the final photon energy , @xmath68 is the total energy of the final state , @xmath69 is the mass of the initial state , the spatial matrix element @xmath70 involves the initial and final radial wavefunctions , and the angular matrix element @xmath71 is c_fi=(, ) ( 2 + 1 ) \ {    1 } ^2 . this is the result of ref.@xcite , except for our inclusion of the relativistic phase space factor of @xmath72 , which is usually not far from unity . ( the gi results do not include this phase space factor . ) we evaluate these radiative partial widths in both the nr and gi models . for the nr model the matrix elements latexmath:[$\langle { n'}^{2{\s}'+1}{\l}'_{{\j } ' } |\ ; r \ ;     linear plus smeared hyperfine wavefunctions of the potential model described in sec.iia , and for the gi model they were evaluated using the wavefunctions of ref.@xcite . since the masses predicted for unknown states differ in the two models , the assumed photon energy e@xmath67 differs as well ; these photon energies are given in the e1 and m1 transition tables ( tables  [ e1rad_stop]-[m1rad ] ) together with the radiative partial widths . some e1 transitions that are of special importance for the study of higher charmonium states are discussed in the text . transitions from initial @xmath74 states are of greatest interest in this regard , since these can be studied with high statistics at @xmath16 machines . these can provide access to the spin - triplet members of the 2p and 3p multiplets in particular , starting from the @xmath75 and @xmath76 . e1 radiative transitions may also be useful in identifying the narrow 1d @xmath77 and @xmath78 states , since they are all predicted to have large partial widths ( _ ca . _ 300  kev ) to the 1p @xmath79 and @xmath15 states . although m1 rates are typically rather weaker than e1 rates , they are nonetheless interesting because they may allow access to spin - singlet states that are very difficult to produce otherwise . it is also interesting that the known m1 rates show serious disagreement between theory and experiment . this is in part due to the fact that m1 transitions between different spatial multiplets , such as @xmath80 ( 2s @xmath81 1s ) , are nonzero only due to small relativistic corrections to a vanishing lowest - order m1 matrix element . the m1 radiative partial widths are evaluated using    _ m1 ( n ^2+1l _  ^2+1l_ + ) = 2j+12l+1 _  _  ,  1 e_c^2 |_f | _ i     ( see the previous e1 formula for definitions . ) the gi m1 radiative rates do not incorporate the @xmath72 phase space factor , but do include a @xmath82 recoil factor . we quote nr results both with and without this recoil factor . the photon energies e@xmath67 depend on the model in most cases , since we have assumed theoretical masses for unknown states .    as with the e1 rates , we quote m1 radiative widths in both the nr potential model of sec.iia and the gi model @xcite of sec.iib . although the m1 decay rates involve a relatively simple off - diagonal matrix element of the magnetic moment operator , we note that our predicted rates show considerable variation with the model assumptions , and agreement with the two known m1 rates is not good . in the off - diagonal ( 2s @xmath81 1s ) transition @xmath80 this is because the zeroth - order m1 matrix element vanishes due to orthogonality of the spatial wavefunctions , and the nonzero predicted rate results from rather model - dependent corrections to the wavefunctions and recoil effects . evidently the corrections we have included do not accurately predict the observed partial width ; for this reason the other predicted m1 rates between different spatial multiplets are suspect , and in any case are evidently strongly dependent on recoil factors . better experimental data will be very important for improving our description of these apparently simple but evidently poorly understood m1 radiative transitions . -0.2 cm      -0.2 cm the dominant strong decays ( when allowed by phase space ) are transitions to open - flavor final states . in these open - flavor decays the initial meson decays through production of a light @xmath24 quark - antiquark pair @xmath83 , followed by separation into two open - charm mesons . remarkably , the qcd mechanism underlying this dominant decay process is still poorly understood . in quark model calculations this decay process is modelled by a simple phenomenological pair production amplitude . the pair is usually assumed to be produced with vacuum @xmath84 quantum numbers , and variants of the decay model make different assumptions regarding the spatial dependence of the pair production amplitude relative to the initial pair . the simplest of these models is the @xmath1p@xmath2 model , originally introduced by micu  @xcite , which assumes that the new @xmath24 pair is produced with vacuum quantum numbers ( @xmath1p@xmath2 ) by a spatially constant pair - production amplitude @xmath85 . leyaouanc  _ et al . _ subsequently applied the @xmath1p@xmath2 model to meson @xcite and baryon @xcite open flavor strong decays in a series of publications in the 1970s . they also evaluated strong decay partial widths of the three states @xmath50 , @xmath75 and @xmath76 in the @xmath1p@xmath2 model @xcite ; the relation between this early work and our contribution will be discussed . this decay model , which has since been applied extensively to the decays of light mesons and baryons , was originally adopted largely because of the success in predicting the d / s amplitude ratio in the decay @xmath86 @xcite . another stringent test of strong decay models is provided by the very tight limit on the spin singlet to spin singlets transition @xmath87 from the ves collaboration @xcite , @xmath88 this branching fraction is predicted to be zero in the @xmath1p@xmath2 model , but would not necessarily be negligible in a different decay model or if final state interactions are important . ( final state interactions combined with the @xmath1p@xmath2 model allow two - stage transitions such as @xmath89 , although the direct decay @xmath90 is forbidden in the @xmath1p@xmath2 model . ) recent variants of the @xmath1p@xmath2 model modify the pair production vertex  @xcite or modulate the spatial dependence of the pair production amplitude to simulate a gluonic flux tube @xcite . ( the latter is the  flux - tube decay model " , which gives very similar predictions to the @xmath1p@xmath2 model in practice . ) another class of decay models assumes that the pair production amplitude transforms as the time component of a lorentz vector . the cornell group used a decay model of this type for charmonium @xcite . this model appears to describe the partial and total widths of some charmonium states reasonably well , although it is known to disagree with experimental amplitude ratios in the light quark sector @xcite .    in this work we employ a formalism that is equivalent to the original constant - amplitude version of the @xmath1p@xmath2 decay model , although we have simplified the calculations by deriving momentum - space feynman rules @xcite instead of using the real - space convolution integrals of leyaouanc  _ et al . _ @xcite .    in our formalism the @xmath1p@xmath2 model describes decay matrix elements using a pair production hamiltonian which is the nonrelativistic limit of h_decay = _ q 2 m_q d^ 3 x  |_q _ q  . here @xmath91 is a dirac quark field , @xmath92 is the constituent quark mass , and @xmath85 is the dimensionless pair production amplitude , which is fitted to data . the decay arises from the matrix element of this model hamiltonian between an initial meson state @xmath93 and a final meson pair @xmath94 , which is nonzero due to the @xmath95 term . here we will use sho ( simple harmonic oscillator ) wavefunctions for the mesons , with a universal width parameter @xmath96 . the pair - production strength parameter @xmath85 is found to be roughly flavor - independent in light meson decays involving pair production of @xmath97 , @xmath98 and @xmath99 pairs . in our recent extensive studies of light ( @xmath100 ) meson decays @xcite we found that the value = 0.4 gives a reasonably accurate description of the overall scale of decay widths . the open - flavor decay amplitudes of the three states @xmath50 , @xmath75 and @xmath76 were evaluated in the @xmath1p@xmath2 model by leyaouanc  _ et al . _ in the late 1970s @xcite . their calculations are generally quite similar to the approach used here , although we note three important differences :    0.2 cm 1 ) leyaouanc  _ et al . _ fitted the pair production strength @xmath85 and wavefunction length scale @xmath101 to the charmonium decay data , whereas we used the same pair production amplitude @xmath85 as in light meson decays , and our sho wavefunction length scale @xmath102 is taken from potential models . 0.2 cm 2 ) we use the correct reduced mass coordinates for the unsymmetric open - charm meson wavefunctions , whereas leyaouanc  _ et al . _ assumed symmetric sho wavefunctions for all mesons . this approximation proves to be quite important numerically , since @xmath103 . 0.2 cm 3 ) we consider all energetically allowed open - flavor decay modes . leyaouanc  _ et al . _ did not consider some s+p decays of the @xmath76 , nor did they consider the important j@xmath104 d meson singlet - triplet mixing angle . ( the masses of the p - wave charmed mesons were unknown at that time . ) 0.2 cm although we have improved on the earlier calculations of leyaouanc  _ et al . _ , and evaluate a much more extensive set of decay amplitudes , we do concur with their important conclusion that the dd mode of the @xmath75 is naturally suppressed given a @xmath105s@xmath48 assignment , and also that s+p modes are important for the @xmath76 . other more recent applications of the @xmath1p@xmath2 decay model to charmonium include a study of the strong decays of various @xmath0 x(3872 ) candidates to dd @xcite , and an examination of the hypothesis that the @xmath75 and @xmath106 are linear superpositions of a hybrid and a conventional @xmath107 state @xcite . the only other detailed studies of the open - flavor strong decays of charmonia of which we are aware are the early work by the cornell group @xcite , and recent more detailed applications of this model to other charmonium states @xcite . the decay model of the cornell group assumes a pair creation operator which is the time component of a lorentz vector ( @xmath108 ) , rather than the lorentz scalar assumed by the @xmath1p@xmath2 model . as we noted previously , this @xmath108 model does not agree well with some light meson decay amplitude ratios . in the original cornell model decay calculations only the combined channel contributions to @xmath101 were derived , and only channels containing pairs of s - wave open - charm mesons were included . the very interesting relative amplitudes in the d@xmath9d@xmath9 channel have not been evaluated in this model . tables  [ strong_3s4s]-[strong_1 g ] present the partial widths and strong decay amplitudes we find for all kinematically allowed open - flavor decay modes of all the charmonium states listed in table  [ table_spectrum ] . we have evaluated these in the @xmath1p@xmath2 model , assuming sho wavefunctions with a width parameter of @xmath109  gev and a pair production amplitude of @xmath110 . sho wavefunctions have the advantage that the decay amplitudes can be determined analytically , and we have found that the numerical results are usually not strongly dependent on the details of the spatial wavefunctions @xcite , unless they are near a node . this width parameter was chosen by comparing the overlap of nr and gi potential model wavefunctions with sho wavefunctions , which was found to be largest for both charmonium and open - charm mesons for a @xmath96 near 0.5  gev .    as a test of the accuracy of these parameters , especially the value assumed for the pair production amplitude @xmath85 ( which is taken from light meson decays ) , in fig . [ gamma ] we compare the predicted total widths for @xmath109  gev and variable @xmath85 to the values observed for the four known states above open - charm threshold , @xmath50 , @xmath75 , @xmath51 and @xmath76 . although there is some scatter in the value of @xmath85 specified by these widths , a choice of @xmath111 is evidently near optimum , and yields an average error of only 29% , which is very reasonable for this phenomenological model . in the decay tables we will quote numerical results for our standard light - meson value of @xmath110 , which is rather close to this value . a subsequent change in @xmath85 will simply scale the widths as @xmath112 .    for this paper the strong decay amplitudes were first evaluated analytically , following which numerical values were determined given our parameters @xmath109  gev , @xmath110 , and our constituent quark masses @xmath113  gev , @xmath114  gev and @xmath115  gev . ( quark mass _ ratios _ are required to specify the final open - charm meson wavefunctions . ) we determine the decay kinematics using external meson masses taken from the particle data group @xcite ( charge averaged where appropriate ) and from recent belle results @xcite . these masses are m@xmath116 gev , m@xmath117 gev , m@xmath118 gev @xcite , m@xmath119(narrow)@xmath120 gev , m@xmath121(broad)@xmath122 gev @xcite , m@xmath123 gev , m@xmath124 gev , m@xmath125 gev , m@xmath126 gev , m@xmath127 gev . the focus collaboration @xcite estimate a somewhat higher mass for the scalar @xmath128 , however their result is complicated by the presence of both scalar and broad axial vector ( @xmath129 ) contributions @xcite . the @xmath130 axial vector @xmath131 and @xmath132 mesons @xmath133 and @xmath129 are assumed to be coherent superpositions of quark model spin - singlet and spin - triplet states ,    @xmath134    we define this singlet - triplet mixing angle in the ls coupling scheme , in accord with the conventions of ref.@xcite . this reference also discusses other conventions for this mixing angle that have appeared in the literature . in the heavy - quark limit we expect to find a  magic " mixing angle , due to the quark mass dependence of the spin - orbit and tensor terms , which is @xmath135 if the expectation of the heavy - quark spin - orbit interaction is negative ( positive ) @xcite . in the following we assume the first case , which is supported by the reported widths of the @xmath136 p - wave @xmath131 mesons . we note however that finite quark mass effects and mixing induced by higher - order fock states can substantially modify this mixing angle , so it should more generally be treated as a free parameter . we will discuss the dependence of our strong decay amplitudes on this mixing angle for a few sensitive cases . the four known states above dd threshold , @xmath50 , @xmath75 , @xmath51 and @xmath76 are of special interest because they are easily produced at @xmath16 machines . accordingly we will first discuss our predicted strong decay amplitudes and widths for these states . this will be followed by a more general discussion of strong decays by multiplet . the @xmath50 is generally assumed to be the @xmath1d@xmath48 state , perhaps with a significant @xmath49s@xmath48 component @xcite . ( this additional component can explain the leptonic width , which is much larger than expected for a pure @xmath1d@xmath48 state @xcite . )    for a pure @xmath1d@xmath48 state at the @xmath50 mass we predict a dd width of 43  mev with our parameters , which is rather larger than the experimental value of @xmath137  mev . the partial width of a mixed 2s - d state    function of the mixing angle @xmath138 . in this simple model , fitting the experimental @xmath50 width requires a mixing angle of @xmath139  @xmath140  @xmath141 . assuming that the leptonic widths scale as the s - wave component squared , this mixing angle predicts an @xmath16 partial width ratio of |_thy . = 0.10 0.03  . it is interesting that this is consistent with experiment , |_expt . = 0.12 0.02  , so the small @xmath50 total width may indeed be an effect of @xmath1d@xmath48 - 2@xmath1s@xmath48 mixing .     as a function of the @xmath1d@xmath48 - 2@xmath1s@xmath48 mixing angle @xmath138 . the experimental @xmath142 limits are shown as horizontal dashed lines.,width=302 ]      next in mass is the @xmath75 , which is a very interesting case for the study of strong decays . four open - charm modes are energetically allowed , dd , dd@xmath9 , d@xmath9d@xmath9 and d@xmath143d@xmath143 , although d@xmath9d@xmath9 has little phase space . there is some experimental evidence for the three nonstrange modes from mark  i at slac @xcite . remarkably , the reported relative branching fractions ( scaled by @xmath144 ) show a very strong preference for d@xmath9 final states , @xmath145 . this motivated suggestions that the @xmath75 might be a @xmath146 molecule @xcite . the mark  i results were ( ^*0^*0   :  ^0^*0  :   ^0 ^ 0 )  = -0.8 cm 12840 : 4.0 0.8 : 0.20.1  . 1.0 cm this is much less striking once the @xmath144 factors are restored , which gives ( ^*0^*0   :  ^0^*0   :   ^0 ^ 0 )  = -0.8 cm 10.31 : 0.95 0.19 : 0.120.06  , 0.8 cm where we have used the @xmath147 mode to set the scale . these relative branching fractions should be compared to our theoretical @xmath75 branching fractions , not the more often quoted @xmath148 ratios . comparison with table  [ strong_3s4s ] shows that the predicted @xmath1p@xmath2 model partial and total widths for a 3@xmath1s@xmath48 @xmath75 state are actually in accord with the mark  i results . the total width is predicted to be 74  mev , which is reasonably close to the pdg experimental average of @xmath149  mev . seth @xcite estimates a rather larger width of @xmath150  mev for this state from crystal ball and new bes data , which is closer to our predicted width . the @xmath151 and @xmath146 branching fractions are predicted to be approximately equal , consistent with mark  i. the  dd mode is predicted to be much smaller , despite its larger phase space , as it is accidentally near a node in the @xmath1p@xmath2 model decay amplitude . ( this was previously noted by leyaouanc _ _ @xcite . ) this set of predictions constitutes a very nontrivial success of the @xmath1p@xmath2 decay model . the actual width to the suppressed dd mode is probably a few mev ; the 0.1  mev quoted in the table depends strongly on the location of the node , and varying the sho width parameter @xmath96 by @xmath152 suggests that a dd partial width of a few mev is plausible . this is again consistent with the mark  i result . the @xmath146 mode is especially interesting for strong decay studies , since there are three independent decay amplitudes for @xmath153 , @xmath154p@xmath48 , @xmath155p@xmath48 and @xmath155f@xmath48 . if the @xmath75 is a pure s - wave state , we predict a zero @xmath155f@xmath48 @xmath146 decay amplitude , and the ratio of the nonzero amplitudes is @xmath155p@xmath156p@xmath157 ( independent of the radial wavefunction ) . as we shall see when we discuss the @xmath51 , very different @xmath146 amplitude ratios are predicted for a @xmath1d@xmath48 initial state .    the unobserved mode d@xmath143d@xmath143 is predicted to have a branching fraction of about @xmath158 ; this prediction is fortunately not especially sensitive to the location of the decay amplitude node for @xmath75 to two pseudoscalars . this branching fraction is of special interest because it determines the event rates available for studies of @xmath159 weak decays at @xmath16 colliders such as cleo . finally , we note that a large sample of @xmath75 events might be used to access the 2p multiplet , since the e1 radiative branching fractions for @xmath160 are expected to be _ ca . _ @xmath161 ( table  [ e1rad_stop ] ) . comparison of the mass of the @xmath74 @xmath51 with potential model predictions ( fig . [ spectrum ] ) immediately suggests a 2@xmath1d@xmath48 assignment . there may also be a significant s - wave component , since the @xmath51 has a much larger @xmath16 width than one would expect for a pure d - wave state @xcite . ( this was also the case for the @xmath50 . ) there are five open - flavor decay modes available to a @xmath74 vector at this mass , dd , dd@xmath9 , d@xmath9d@xmath9 , d@xmath143d@xmath143 and d@xmath162d@xmath163 . the predicted decay amplitudes and partial widths for a 2@xmath1d@xmath48 @xmath51 are given in table  [ strong_1d2d ] . the theoretical total width of 74  mev is in good agreement with the experimental value of @xmath164  mev . seth @xcite quotes a similar width of @xmath165  mev for this state from crystal ball and recent bes data . the leading mode is predicted to be d@xmath9d@xmath9 with a branching fraction of @xmath166 , followed by comparable dd and d@xmath162d@xmath163 modes ( both @xmath167 ) , and a somewhat weaker d@xmath143d@xmath143 ( @xmath168 ) . the dd@xmath9 branching fraction is predicted to be very small , since it is suppressed by a decay amplitude node near the physical point . the mode d@xmath9d@xmath9 is again especially interesting , due to the three decay amplitudes allowed for this final state . ( the other modes have only single amplitudes . ) for a pure d - wave assignment the ratio of the two @xmath146 p - wave amplitudes is independent of the radial wavefunction , and is @xmath155p@xmath156p@xmath169 . the @xmath155f@xmath48 amplitude is predicted to be the largest given this 2@xmath1d@xmath48 assignment , whereas it is zero for an s - wave state . clearly , a determination of these d@xmath9d@xmath9 decay amplitude ratios would be an extremely interesting test of decay models . experimentally , to date nothing has been reported regarding the exclusive hadronic decay modes of the @xmath51 . the final resonance known above dd threshold is the @xmath74 @xmath76 . again , the mass of this resonance relative to potential model predictions suggests a assignment , in this case 4@xmath1s@xmath48 . of course this identification requires independent confirmation , since hybrids are predicted to first appear near this mass by lgt simulations @xcite , and the lightest hybrid multiplet in the flux tube model includes a @xmath74 state @xcite . ten open - charm strong decay modes are allowed for the @xmath76 , seven with @xmath131 meson final states ( @xmath170 ) , and three with @xmath132 . the predicted branching fractions of a 4@xmath1s@xmath48 state at this mass are quite characteristic , and ( if the @xmath1p@xmath2 decay model is accurate ) may be useful in confirming this assignment . as a note of caution , a 4s state has three radial nodes , and some of the smaller predicted branching fractions are sensitive to the locations of the nodes . the predicted total width of a 4@xmath1s@xmath48 meson at this mass with our parameters is 77  mev , somewhat larger than the experimental pdg average width of @xmath171  mev . seth @xcite notes however that the width of this state in crystal ball and recent bes data is rather larger , and quotes an average of @xmath172  mev , on the opposite side of the @xmath1p@xmath2 decay model prediction . the largest exclusive mode is predicted to be the s+p combination dd@xmath48 , where d@xmath48 is the narrower of the two @xmath136 @xmath131 axial mesons near 2.42 - 2.43  gev . since the d@xmath48 is rather narrow ( @xmath173-@xmath174  mev ) and decays dominantly to d@xmath175 , there should be a strong @xmath76 signal in dd@xmath175 final states . although both s - wave and d - wave dd@xmath48 final states are allowed , in the @xmath1p@xmath2 model the hqet ( heavy quark effective theory ) d@xmath48 mixing angle @xmath138 is just the value needed to give a zero s - wave @xmath176 amplitude . thus we have the striking prediction that the dominant @xmath76 decay mode is dd@xmath48 , in d - wave rather than s - wave . the second - largest @xmath76 branching fraction is predicted to be another s+p mode , dd@xmath177 . the d@xmath177 is also moderately narrow , so this isobar can also be isolated from the observed final state . ( the pdg quotes a neutral d@xmath177 total width of @xmath178  mev , but belle @xcite and focus @xcite find somewhat broader values of @xmath179  mev and @xmath180  mev respectively ; see @xcite for a recent discussion . the d@xmath177 has significant branching fractions to both d@xmath175 and d@xmath181 , so the dd@xmath177 mode of the @xmath76 should be observable in both dd@xmath181 and dd@xmath175 .    a final important @xmath76 mode is predicted to be d@xmath9d@xmath9 , which should be comparable in strength to dd@xmath177 . as noted in the previous discussions of @xmath75 and @xmath51 decays , d@xmath9d@xmath9 is an especially interesting decay mode because it has three amplitudes , and the @xmath155p@xmath48/@xmath154p@xmath48 amplitude ratio is independent of the radial wavefunction for pure s - wave or d - wave states . if the @xmath76 is indeed an s - wave ( 4@xmath1s@xmath48 ) state to a good approximation , we expect this ratio to be @xmath155p@xmath48/@xmath154p@xmath157 , and the @xmath155f@xmath48 amplitude should be zero . it is interesting to note that @xmath76 decays may provide access to the recently discovered d@xmath182 . although the channel d@xmath183d@xmath184 has a threshold of 4429  mev , 14  mev above the nominal mass of the @xmath76 , the width of the @xmath76 and the fact that the decay @xmath185d@xmath183d@xmath184 is purely s - wave implies that one may observe significant d@xmath184 production just above threshold , near e@xmath186  mev . this is illustrated in fig . [ 4415_scan ] , in which we show theoretical @xmath187 decay model partial widths to d@xmath188d@xmath189 and d@xmath143d@xmath190 as functions of the @xmath191s@xmath48 mass . ( this calculation assumes pure @xmath132 d@xmath189 and d@xmath190 states , and should be modified accordingly if they have significant non-@xmath132 components . ) unfortunately , @xmath76 decays are not expected to be similarly effective in producing the d@xmath190 , because the assumed  magic - mixed " hqet state @xmath192 predicts a vanishing @xmath193 s - wave decay amplitude .     and @xmath194 as a function of the assumed @xmath195 mass . , width=321 ]      the two unknown states in the 3s and 4s multiplets are the 3@xmath154s@xmath2 and 4@xmath154s@xmath2 pseudoscalars . to evaluate their strong decays we have assigned potential model masses to these states , specifically 4043  mev and 4384  mev . the resulting total widths are predicted to be moderate , 80  mev and 61  mev respectively , and both states should be observable in dd@xmath9 and d@xmath9d@xmath9 ( see table  [ strong_3s4s ] ) . no other open - charm strong decay modes are allowed for a 3@xmath154s@xmath2 @xmath196 . although these modes are also important for the 4@xmath154s@xmath2 @xmath197 , its largest branching fraction is predicted to be to the s+p combination dd@xmath177 ( @xmath198 ) . these states may be observable in @xmath6 and hadronic production ( through @xmath199 fusion ) , analogous to the @xmath200 and @xmath5 . m1 transitions from the higher vectors to these states are unfortunately predicted to be rather weak , since they are either hindered or have little phase space . ( see table  [ m1rad ] . )      the mean 2p multiplet mass is predicted to be near 3.95  gev . there is an interesting disagreement between the nr and gi models regarding the scale of mass splittings within this multiplet , as well as in the 3p states . ( see table  [ table_spectrum ] and fig . [ spectrum ] . ) there are few open - flavor strong modes available to the 2p states ( table  [ strong_2p3p ] ) . one result of this restricted phase space is the prediction of a fairly narrow 2@xmath1p@xmath2 @xmath201 ( assuming the nr model mass ) that decays only to dd , with a total width of just 30  mev . this small width is due in part to a decay amplitude node , and a variation of the total width by a factor of two could easily be accommodated . the two axial states 2@xmath1p@xmath48 @xmath202 and 2@xmath154p@xmath48 @xmath203 can only decay to dd@xmath9 , and both have nonzero s- and d - wave decay amplitudes ; their d / s ratios are @xmath204 and @xmath205 respectively . observation of the d / s ratio in a dd@xmath9 enhancement could be a useful check of the assumption of a resonant @xmath136 contribution . the relatively narrow 2@xmath1p@xmath2 state @xmath201 can be produced in @xmath6 , as can its somewhat wider 2@xmath1p@xmath57 @xmath206 partner . additional possible production mechanisms include @xmath199 fusion ( for the 2@xmath1p@xmath2 and 2@xmath1p@xmath57 ) , e1 radiative transitions from higher @xmath74 states ( to all 2p @xmath79 states ) , and b decays . the e1 branching fractions of the @xmath75 and @xmath51 to 2p @xmath79 states are predicted to be sufficiently large ( _ ca . _  @xmath161 , see tables  [ e1rad_stop],[e1rad_1d2d ] ) to allow identification of these states at bes and cleo . although no 2p @xmath0 state has been clearly established experimentally , there are two recent reports of enhancements which may be due to states in this multiplet . the belle collaboration has reported evidence of an enhancement in @xmath12 with a mass and width of @xmath207  mev and @xmath208  mev @xcite ; this is compatible with expectations for 2p @xmath79 @xmath0 states . ( although 2p @xmath0 states should strongly favor open - charm strong decay modes , inelastic fsis ( final state interactions ) will allow weaker transitions to closed - charm final states such as @xmath12 . these fsis should be most important in s - wave , as in the virtual second - order processes @xmath209 . ) a search for evidence of this state in open - charm final states , with much larger branching fractions than to @xmath12 , will be a crucial test of a charmonium assignment . assuming that the enhancement is indeed due to a 2p @xmath79 state , if it is the @xmath210 state it will populate only dd , if @xmath13 only dd@xmath9 , and if @xmath211 comparable branching fractions to dd and dd@xmath9 are expected .    there is a second report of a possible 2p state , also from the belle collaboration @xcite , in double charmonium production . an enhancement is seen in @xmath212 at a mass of @xmath213  mev , with a limit on the total width of @xmath214  mev ( @xmath215 ) . as this enhancement is not seen in @xmath12 , it has been suggested @xcite that it is distinct from the 3940  mev signal discussed above . consideration of possible 2p candidates ( table  [ strong_2p3p ] ) suggests that this enhancement may be due to the @xmath216p@xmath217 , which is predicted to have a similar mass of 3925  mev . ( the predicted total width however is a somewhat larger 165  mev . ) the 3p states have an expected mean multiplet mass of about 4.3  gev . as with the 2p states there is a significant difference between the nr and gi models in the scale of mass splittings predicted within the multiplet . this leads for example to a 90  mev difference between the models for the predicted mass of the 3@xmath1p@xmath2 scalar . many open - charm strong decay channels are open at the 3p mass scale . the individual decay amplitudes however are typically somewhat smaller than for the 2p states , due in part to nodes in the decay amplitudes . this has the interesting result that the predicted mean total width of 3p states is smaller than for 2p states ( 58  mev for 3p versus 91  mev for 2p ) . the predicted branching fraction to the mode d@xmath9d@xmath9 is large for all four 3p states , and is the largest for all except the spin - singlet 3@xmath154p@xmath48 @xmath218 ; in this case the leading mode is dd@xmath219 . for all 3p states except the 3@xmath1p@xmath48 the d@xmath9d@xmath9 decay mode has multiple amplitudes , which are predicted to be comparable in strength . an experimental determination of any of these amplitude ratios would provide an interesting test of decay models . the decays of the tensor 3@xmath1p@xmath57 @xmath220 to dd@xmath48 and dd@xmath221 are also interesting , in that the prediction that the broad d@xmath221 mode dominates is a test of the axial - d mixing angle .    as with the 2p states , the most important experimental problem is to find an effective production mechanism . again , @xmath6 and @xmath199 fusion can produce the @xmath211 and @xmath210 states ( here 3@xmath1p@xmath57 and 3@xmath1p@xmath2 ) , although there should be moderate suppression of the production amplitudes relative to the 1p ground states due to the higher masses and smaller short - distance wavefunctions . e1 transitions from the @xmath76 could lead to the identification of the 3p states , since radiative partial widths of 50 - 100  kev are expected to the 3@xmath1p@xmath57 and 3@xmath1p@xmath48 ; these correspond to branching fractions of _ ca . _  @xmath161 . production of the 3@xmath1p@xmath2 state however is expected to be suppressed , due to a radiative decay amplitude node . since the @xmath1d@xmath57 and @xmath154d@xmath57 states do not have allowed open - flavor decay modes , and we have already considered @xmath1d@xmath48 decays in the section on the @xmath50 , we restrict our discussion to the @xmath1d@xmath222 . this currently unknown state is predicted to have a mass of 3806  mev in the nr potential model . it is expected to be quite narrow simply because it has little phase space to the only open - flavor mode dd , which has an f - wave centrifugal barrier ; with our parameters the predicted total strong width is just 0.5  mev . the radiative width for the e1 transition @xmath223 is predicted to be about 0.3  mev ( table  [ e1rad_1d2d ] ) , so the total width of this @xmath1d@xmath222 state should be near 1  mev , with comparable strong and radiative branching fractions .    finding an effective production mechanism for this interesting state is again a crucial problem . as it has c  @xmath224 , it can not be made in @xmath225 or @xmath199 fusion . one possibility is that the @xmath1d@xmath222 may be observable in e1 transitions from the radial @xmath49p@xmath57 tensor , since the radiative partial width is predicted to be a relatively large @xmath53  kev . similarly , the narrow @xmath1d@xmath57 and @xmath154d@xmath57 states may be observable in e1 transitions from 2p states , assuming that a production mechanism can be found for 2p states . all three of these states can also be produced in b - decays @xcite . the states in the 2d multiplet are predicted to be essentially degenerate in both the nr and gi potential models . the @xmath49d@xmath48 candidate @xmath51 suggests a 2d multiplet mass near 4.16  gev , which is consistent with potential model expectations . since the @xmath51 was discussed previously , here we will consider only the unknown states @xmath49d@xmath222 , @xmath49d@xmath57 and @xmath226d@xmath57 . these states are predicted to be rather broad , with theoretical total widths ranging from 78  mev for the @xmath49d@xmath48 ( consistent with the @xmath51 ) to 148  mev for the @xmath49d@xmath222 . the energetically allowed open - flavor modes are all of s+s type , since the first s+p mode is dd@xmath227 , at a nominal mass of 4175  mev . the decays of 2d states are dominated by dd , dd@xmath9 and d@xmath9d@xmath9 , with smaller contributions from the charm - strange final states d@xmath143d@xmath143 and d@xmath162d@xmath163 . if these states ( besides the @xmath51 ) can be produced , there are several interesting tests of their strong decay amplitudes . the @xmath49d@xmath222 @xmath228 , for example , is predicted to have d@xmath9d@xmath9 as the dominant mode , with both @xmath154f@xmath222 and @xmath155f@xmath222 amplitudes present , in the ratio @xmath155f@xmath222/@xmath154f@xmath229 . the remaining 2d states have many multiamplitude decays , in which the higher - l partial waves have large or dominant amplitudes . ( compare the p- and f - waves in the decays in table  [ strong_1d2d ] . )      the four states in the 1f multiplet are predicted to be almost degenerate in both the nr and gi models , although the models disagree somewhat regarding this mass . in the nr potential model the states are expected near 4025  mev , whereas in the gi model the expected mass is near 4095  mev . this discrepancy is mainly due to the different values assumed for the string tension . it is unfortunate that production amplitudes for these higher - l states are expected to be quite weak , since the 1f multiplet is predicted to contain a rather narrow state , the @xmath1f@xmath230 @xmath231 ( table  [ strong_1f ] ) . assuming our nr model @xmath1f@xmath230 mass of 4021  mev , we predict a strong width for this state of just 8.3  mev . this is partly a result of the centrifugal barrier ( both open modes , dd and dd@xmath9 , are g - waves ) , and is also because there is essentially no phase space for d@xmath9d@xmath9 . if the gi masses are more accurate we would no longer expect a very narrow @xmath1f@xmath230 state , since the d@xmath9d@xmath9 partial width increases rapidly above threshold . at the gi mass of 4095  mev the @xmath1f@xmath230 is predicted to have a partial width to d@xmath9d@xmath9 of 38  mev , and a total width of 60  mev . the single well - established mechanism that might be exploited to produce a 1f state is an e1 radiative transition from the @xmath51 to the @xmath1f@xmath57 @xmath232 . this transition however has a rather weak , model - dependent partial width of 20 - 50  kev ( table  [ e1rad_1d2d ] ) , and the @xmath1f@xmath57 is predicted to be the broadest 1f state , with a total width of about 160  mev . if an effective production mechanism for 1f states is identified , the e1 decays of these states can be used to populate the narrow 1d states ; the e1 transitions @xmath1f@xmath233d@xmath222 , @xmath1f@xmath234d@xmath57 and @xmath154f@xmath235d@xmath57 are predicted to have rather large radiative widths of _ ca . _  300 - 400  kev ( see table  [ e1rad_1f ] ) .      to complete our study we have evaluated the spectrum and decays of all states in the 2f and 1 g multiplets ( tables  [ strong_2fa ] and [ strong_1 g ] ) , although the lack of a clear experimental route to these states suggests that it may be quite difficult to test these predictions . the multiplets are again predicted to be nearly degenerate in both the nr and gi potential models , with somewhat lower masses predicted by the nr model . the mean 2f and 1 g masses are approximately 4350  mev and 4225  mev respectively in the nr model , and 4425  mev and 4320  mev in the gi model . assuming nr model masses , these states are predicted to have total widths that range from a minimum of 58  mev for the @xmath1g@xmath236 @xmath237 to a maximum of 180  mev for the @xmath49f@xmath57 @xmath238 . the individual decay modes do show some interesting features ; although the dominant modes are usually of s+s type , some s+p modes are important in 2f decays . for example , the largest decay modes of the @xmath49f@xmath222 @xmath239 and @xmath49f@xmath57 @xmath238 are both s+p , with branching fractions of b(@xmath240 dd@xmath241 and b(@xmath242 dd@xmath243 respectively . we have computed the spectrum , all allowed e1 and some m1 electromagnetic partial widths , and all allowed open - charm strong decay amplitudes of the 40 charmonium states expected to _ ca . _ 4.4  gev . these were evaluated using two potential models , the quark model formulation of electromagnetic transition amplitudes , and the @xmath187 strong decay model . the predictions of the spectrum should be useful in the identification of new states , and in tests of the lorentz structure of confinement and the nature of spin - dependent forces . the transition amplitudes will be useful in searches for new higher - mass states , as they suggest which final states should be populated preferentially by the decays of a given initial state , as well as predicting the often characteristic strong decay amplitude ratios within a given final state . these predictions may also be useful in distinguishing between conventional charmonia and exotica such as charmed meson molecules ( perhaps including the x(3872 ) ) and hybrids . these results all implicitly test the accuracy of the pure - assumption made in these models , and can serve as benchmarks for calculations that relax this assumption . future experimental measurements of strong partial widths and decay amplitude ratios can provide very important tests of strong decay models , the @xmath1p@xmath2 model in particular . we have discussed many specific examples of these tests in the previous sections ; here we will summarize some of our most important observations . the @xmath50 is the lightest known state above open - charm threshold . we noted that the relatively large @xmath16 width of the @xmath50 is difficult to understand if it is assumed to be a pure @xmath1d@xmath48 state , but this problem is solved if the @xmath50 has a significant admixture of the @xmath244 @xmath0 basis state . this mixing also solves the @xmath50 total width problem , since the mixing angle required to fit the @xmath16 width is consistent with the value required to fit the @xmath1p@xmath2 model dd strong width to experiment . this mixing angle has been discussed in earlier references ( see for example refs.@xcite ) , although the observation that the same mixing angle resolves both the @xmath16 and total width discrepancies has not been noted previously . the relative branching fractions of the @xmath75 to dd , dd@xmath9 and d@xmath9d@xmath9 , although often cited as anomalous , are naturally explained in the @xmath187 model by a node in the dd decay amplitude . a very important test of this strong decay model and the nature of the @xmath75 follows from a simple measurement of the ratios of the three independent d@xmath9d@xmath9 decay amplitudes . the @xmath187 model predicts that the @xmath245 amplitude is zero , and the ratio of the two p - wave amplitudes is @xmath246 , independent of the radial wavefunction , provided that the @xmath75 is a pure s - wave state . similarly , the d@xmath9d@xmath9 amplitude ratios in @xmath51 decays test the decay model and the assumed @xmath247 assignment for this state . the @xmath245 d@xmath9d@xmath9 amplitude is predicted to be the largest for a pure @xmath247 @xmath51 , and the p - wave amplitude ratio is predicted to be @xmath248 , again independent of the radial wavefunction . finally , we noted that it may be possible to identify a higher - l state ( the @xmath249 @xmath232 member of the f - wave multiplet ) in e1 radiative @xmath51 decays . the @xmath76 has ten open - flavor modes , and their branching fractions have never been determined . we predict that the largest is the unusual s+p combination dd@xmath48 . it is especially notable that strong decays of the @xmath76 may provide a novel production mechanism for the enigmatic @xmath250 , through the above - threshold decay @xmath251 , as shown in fig . [ 4415_scan ] . one may also produce higher - mass charmonium states from the @xmath76 through electromagnetic transitions . for example , @xmath76 e1 decays are predicted to produce 3p states with branching fractions at the per mil level . these states could then be identified through their subsequent open - charm strong decays .    as discussed in section v.c , two new states near 3940  mev have recently been reported by belle . 2p @xmath79 states are natural assignments , as they have theoretical quark model masses ranging from 3852 to 3972 mev and widths of 50 - 100 mev , similar to the reported values . a  measurement of the relative branching fractions to @xmath252 , @xmath151 and @xmath146 and comparisons to our @xmath1p@xmath2 decay model predictions should allow the determination of the quantum numbers of these states , and will show whether they are indeed consistent with 2p @xmath0 assignments . a better determination of the properties of conventional charmonium states at higher masses is important not only because of the improved understanding of qcd that will follow ( especially aspects of confinement and strong decays ) , but also because qualitatively different types of resonances are expected at these masses . these new states include charmed meson molecules and charmonium hybrids , and the identification of these novel excitations will obviously be easier if the conventional charmonium spectrum is well established . the recent results from b factories and new programs at bes , cleo and gsi have led to a resurgence of interest in the physics of charmonium . we argue that a detailed experimental investigation of the spectrum of excited charmonium states and their decay properties will considerably improve our understanding of the nonperturbative aspects of qcd . this research was supported in part by the natural sciences and engineering research council of canada , the u.s . national science foundation through grant nsf - phy-0244786 at the university of tennessee , and the u.s . department of energy under contracts de - ac05 - 00or22725 at oak ridge national laboratory and de - fg02 - 00er41135 at the university of pittsburgh , and by pparc grant pp / b500607 . j.  j.  aubert _ et al . _ [ e598 collaboration ] , phys .  rev . * 33 * , 1404 ( 1974 ) . j.  e.  augustin _ et al . _ [ slac - sp-017 collaboration ] , phys . * 33 * , 1406 ( 1974 ) . t.  appelquist and h.  d.  politzer , phys . lett .   * 34 * , 43 ( 1975 ) . a.  de rujula and s.  l.  glashow , phys . * 34 * , 46 ( 1975 ) . t.  appelquist , a.  de rujula , h.  d.  politzer and s.  l.  glashow , phys . lett .   * 34 * , 365 ( 1975 ) . e.  eichten , k.  gottfried , t.  kinoshita , j.  b.  kogut , k.  d.  lane and t.  m.  yan , phys . lett .   * 34 * , 369 ( 1975 ) [ erratum - ibid .   * 36 * , 1276 ( 1976 ) ] . c.  quigg , arxiv : hep - ph/0403187 . r.  s.  galik , arxiv : hep - ph/0408190 . n.  brambilla _ et al . _ , arxiv : hep - ph/0412158 . k.  k.  seth , arxiv : hep - ex/0501022 . this reference in effect quotes a combined average mass for the @xmath15 of @xmath253  mev from new cleo measurements .    k.  k.  seth , arxiv : hep - ex/0504050 . k.  k.  seth , arxiv : hep - ex/0504052 . t.  skwarnicki , arxiv : hep - ex/0505050 . e.j.eichten , k.lane and c.quigg , phys . 89 , 162002 ( 2002 ) [ arxiv : hep - ph/0206018 ] . the besiii detector . preliminary design report , ihep - bepch - sb-13 ( jan . 2004 ) . s .- k.choi _ et al . _ ( belle collaboration ) , phys . 89 , 102001 ( 2002 ) , err . 89 , 129901 ( 2002 ) . [ arxiv : hep - ex/0206002 ] . b.  aubert _ et al . _ [ babar collaboration ] , phys . 92 , 142002 ( 2004 ) [ arxiv : hep - ex/0311038 ] . ( cleo collaboration ) , phys . lett.92 , 142001 ( 2004 ) [ arxiv : hep - ex/0312058 ] . s.  k.  choi _ et al . _ [ belle collaboration ] , phys . 91 , 262001 ( 2003 ) [ arxiv : hep - ex/0309032 ] . d.  acosta _ et al . _ [ cdf ii collaboration ] , phys . rev . 93 , 072001 ( 2004 ) [ arxiv : hep - ex/0312021 ] . f.  e.  close and p.  r.  page , phys . lett . b578 , 119 ( 2004 ) [ arxiv : hep - ph/0309253 ] . e.  s.  swanson , phys . b588 , 189 ( 2004 ) [ arxiv : hep - ph/0311229 ] . n.  a.  tornqvist , phys . b * 590 * , 209 ( 2004 ) [ arxiv : hep - ph/0402237 ] . t.  barnes and s.  godfrey , phys . d69 , 054008 ( 2004 ) [ arxiv : hep - ph/0311162 ] . e.  j.  eichten , k.  lane and c.  quigg , phys . d69 , 094019 ( 2004 ) [ arxiv : hep - ph/0401210 ] . k.  abe , arxiv : hep - ex/0505037 . k.  abe , arxiv : hep - ex/0505038 . a.  tomaradze , j.  phys . * 9 * , 119 ( 2005 ) [ arxiv : hep - ex/0410090 ] . the reported @xmath15 signal is in the decay chain @xmath254 , @xmath255 . the masses found in two different inclusive analyses were @xmath256  mev and @xmath256  mev ( with an estimated systematic error of @xmath257  mev , and @xmath258  mev in exclusive decays ( with six different identified @xmath200 final states ) . k.  abe _ et al . _ [ belle collaboration ] , phys . 89 , 142001 ( 2002 ) [ arxiv : hep - ex/0205104 ] . k.  abe _ et al . _ [ belle collaboration ] , phys . d70 , 071102 ( 2004 ) [ arxiv : hep - ex/0407009 ] . p.  pakhlov , [ arxiv : hep - ex/0412041 ] . p.  pakhlov ( private communication ) . e.  eichten , k.  gottfried , t.  kinoshita , k.  d.  lane and t.  m.  yan , phys . d * 17 * , 3090 ( 1978 ) [ erratum - ibid . d * 21 * , 313 ( 1980 ) ] . d.  ebert , r.  n.  faustov and v.  o.  galkin , arxiv : hep - ph/0503012 . c.  w.  bernard _ et al . _ [ milc collaboration ] , phys . d * 56 * , 7039 ( 1997 ) [ arxiv : hep - lat/9707008 ] . x.liao and t.manke , [ arxiv : hep - lat/0210030 ] . z.  h.  mei and x.  q.  luo , int . j.  mod .  phys . a * 18 * , 5713 ( 2003 ) [ arxiv : hep - lat/0206012 ] . g.  s.  bali , eur . j.  a * 19 * , 1 ( 2004 ) [ arxiv : hep - lat/0308015 ] . b.  aubert _ et al . _ [ babar collaboration ] , phys . 90 , 242001 ( 2003 ) [ arxiv : hep - ex/0304021 ] . d.  besson _ et al . _ [ cleo collaboration ] , aip conf . 698 , 497 ( 2004 ) [ arxiv : hep - ex/0305017 ] . e.  s.  swanson , phys . b * 582 * , 167 ( 2004 ) [ arxiv : hep - ph/0309296 ] . d.  ebert , r.  n.  faustov and v.  o.  galkin , phys . d * 67 * , 014027 ( 2003 ) [ arxiv : hep - ph/0210381 ] . t.  barnes , [ arxiv : hep - ph/0406327 ] . li , f.  e.  close and t.  barnes , phys . d * 43 * , 2161 ( 1991 ) . e.  s.  ackleh and t.  barnes , phys . d * 45 * , 232 ( 1992 ) . e.  s.  ackleh , t.  barnes and f.  e.  close , phys .  rev . d * 46 * , 2257 ( 1992 ) . t.  barnes and g.  i.  ghandour , phys . lett . b118 , 411 ( 1982 ) . s.godfrey and n.isgur , phys . d32 , 189 ( 1985 ) . m.  okamoto _ et al . _ [ cp - pacs collaboration ] , phys . d * 65 * , 094508 ( 2002 ) [ arxiv : hep - lat/0112020 ] . s.  hashimoto and t.  onogi , ann . rev .  nucl .  part . sci .   * 54 * , 451 ( 2004 ) [ arxiv : hep - ph/0407221 ] . w.kwong and j.l.rosner , phys . d38 , 279 ( 1988 ) . l. micu , nucl . b10 , 521 ( 1969 ) . a.  le yaouanc , l.  oliver , o.  pne and j .- c . raynal , phys . d8 , ( 1973 ) 2223 . a.  le yaouanc , l.  oliver , o.  pne and j .- c . raynal , phys . d9 , 1415 ( 1974 ) . a.  le yaouanc , o.  pne , j .- c . raynal and l.  oliver , nucl . b149 , 321 ( 1979 ) . a.  le yaouanc , l.  oliver , o.  pne and j .- c . raynal , phys . b71 , 397 ( 1977 ) . le yaouanc , l.  oliver , o.  pne and j .- c . raynal , phys . b72 , 57 ( 1977 ) . r. kokoski and n. isgur , phys . d35 , 907 ( 1987 ) . e.s.ackleh , t.barnes and e.s . swanson , phys . d54 , 6811 ( 1996 ) . [ arxiv : hep - ph/9604355 ] . m.  nozar _ et al . _ [ e852 collaboration ] , phys . lett . b541 , 35 ( 2002 ) [ arxiv : hep - ex/0206026 ] . d.  v.  amelin _ et al . _ [ ves collaboration ] , [ arxiv : hep - ex/9810013 ] . ( particle data group ) , phys . b592 , 1 ( 2004 ) . w.roberts and b.silvestre-brac , phys . d57 , 1694 ( 1998 ) . [ arxiv : hep - ph/9708235 ] . t.barnes , f.e.close , p.r.page and e.s.swanson , phys . d55 , 4157 ( 1997 ) . [ arxiv : hep - ph/9609339 ] . t.barnes , n.black and p.r.page , phys . d68 , 054014 ( 2003 ) [ arxiv : nucl - th/0208072 ] . f.  e.  close and p.  r.  page , phys . b * 366 * , 323 ( 1996 ) [ arxiv : hep - ph/9507407 ] . p.  geiger and e.  s.  swanson , phys . d50 , 6855 ( 1994 ) [ arxiv : hep - ph/9405238 ] . h.  g.  blundell and s.  godfrey , phys . d53 , 3700 ( 1996 ) [ arxiv : hep - ph/9508264 ] . et al . _ [ belle collaboration ] , phys . d69 , 112002 ( 2004 ) [ arxiv : hep - ex/0307021 ] . e.  w.  vaandering [ focus collaboration ] , [ arxiv : hep - ex/0406044 ] . r.kutschke , private communication .    s.godfrey and r.kokoski , phys . d43 , 1679 ( 1991 ) . j.  l.  rosner , [ arxiv : hep - ph/0411003 ] . j.  l.  rosner , [ arxiv : hep - ph/0405196 ] . et al . _ , b69 , 503 ( 1977 ) . m.b.voloshin and l.b.okun , jetp lett . 23 , 333 ( 1976 ) [ pisma zh . fiz . 3 , 369 ( 1976 ) ] . a.derjula , h.georgi and s.l.glashow , phys . 38 , 317 ( 1977 ) . s.iwao , lett . nuovo cim . 28 , 305 ( 1980 ) . k.  k.  seth , [ arxiv : hep - ex/0405007 ] . n.isgur , r.kokoski and j.paton , phys . 54 , 869 ( 1985 ) . j.  m.  link _ et al . _ [ focus collaboration ] , phys . b586 , 11 ( 2004 ) [ arxiv : hep - ex/0312060 ] . p.s.cooper ,  new results in charm meson spectroscopy from focus and selex " , presented at ghp2004 ( fermilab , 24 - 26 oct 2004 ) . k.  abe _ et al . _ [ belle collaboration ] , [ arxiv : hep - ex/0408126 ] . t. skwarnicki , int j. mod . a * 19 * , 1030 ( 2004 ) . y.zhu , ph.d . thesis , california institute of technology , 1988 , caltech report calt-68 - 1513 ( unpublished ) .
Abstract:

Answer: this paper gives results for the spectrum , all allowed e1 radiative partial widths ( and some important m1 widths ) and all open - charm strong decay amplitudes of all 40 @xmath0 states expected up to the mass of the 4s multiplet , just above 4.4  gev . the spectrum and radiative widths are evaluated using two models , the relativized godfrey - isgur model and a nonrelativistic potential model . the electromagnetic transitions are evaluated using coulomb plus linear plus smeared hyperfine wavefunctions , both in a nonrelativistic potential model and in the godfrey - isgur model . the open - flavor strong decay amplitudes are determined assuming harmonic oscillator wavefunctions and the @xmath1p@xmath2 decay model . this work is intended to motivate future experimental studies of higher - mass charmonia , and may be useful for the analysis of high - statistics data sets to be accumulated by the bes , cleo and gsi facilities . h l p  s