Study of the excited $1^-$ charm and charm-strange mesons

We give a systematical study on the recently reported excited charm and charm-strange mesons with potential $1^-$ spin-parity, including the $D^*_{s1}(2700)^+$, $D^*_{s1}(2860)^+$, $D^*(2600)^0$, $D^*(2650)^0$, $D^*_1(2680)^0$ and $D^*_1(2760)^0$. The main strong decay properties are obtained by the framework of Bethe-Salpeter (BS) methods. Our results reveal that the two $1^-$ charm-strange mesons can be well described by the further $2^3\!S_1$-$1^3\!D_1$ mixing scheme with a mixing angle of $8.7^{+3.9}_{-3.2}$ degrees. The predicted decay ratio $\frac{\mathcal{B}(D^*K)}{\mathcal{B}(D~K)}$ for $D^*_{s1}(2860)$ is $0.62^{+0.22}_{-0.12}$.~$D^*(2600)^0$ can also be explained as the $2^3\!S_1$ predominant state with a mixing angle of $-(7.5^{+4.0}_{-3.3})$ degrees. Considering the mass range, $D^*(2650)^0$ and $D^*_1(2680)^0$ are more likely to be the $2^3\!S_1$ predominant states, although the total widths under both the $2^3\!S_1$ and $1^3\!D_1$ assignments have no great conflict with the current experimental data. The calculated width for LHCb $D^*_1(2760)^0$ seems about 100 \si{MeV} larger than experimental measurement if taking it as $1^3\!D_1$ or $1^3\!D_1$ dominant state $c\bar u$. The comparisons with other calculations and several important decay ratios are also present. For the identification of these $1^-$ charm mesons, further experimental information, such as $\frac{\mathcal{B}(D\pi)}{\mathcal{B}(D^*\pi)}$ are necessary.


I. Introduction
Recently lots of natural parity charm and charm-strange mesons are observed in experiments [1][2][3][4][5][6][7][8][9][10], which are summarized in Tab. I, where we have combined the statistical, systematic and model errors in quadrature for simplicity. These new resonances have great importance in improving our knowledge of the radial and orbital charmed excitations. Especially for the spin-parity 1 − charm and charm-strange states, there may exist the 2 3 S 1 -1 3 D 1 mixing , which makes the assignments more complicated.
For the corresponding charm mesons, the GI model [11,12] predicts the 2 3 S 1 and 1 3 D 1 states cū locate in the mass range of about 2.64 and 2.82 GeV, respectively, while the mass of 1 3 D 3 state is predicted to be quite close to the 1 3 D 1 state. BaBar in 2010 [3] reported two natural parity resonances D * (2600) 0 and D * (2760) 0 . Further more, they measured the following ratio of branching fraction as, Again we have introduced an abbreviation R D + for the sake of simplicity. Later in 2013 LHCb [5] discovered two natural parity charmed particle D * (2650) 0 and D * J (2760) 0 . Then in 2015 LHCb [9] reported the 1 − state D * 1 (2760) 0 , which has a large width of 177 ± 38 MeV. Very recently, in 2016 by using the amplitude analysis, LHCb collaboration measured a 1 − state D * 1 (2680) 0 and a 3 − state D * 3 (2760) 0 [10]. The later one's mass and total width seems consistent with the BaBar D * (2760) 0 and LHCb D * J (2760) 0 . These experimental data are also summarized in Tab. I, where the isospin partners of these neutral charm mesons are also listed in the bottom of Tab. I for comparison. The D * (2760) 0 , D * J (2760) 0 and D * 3 (2760) 0 can be interpreted as the same particle, namely, the 3 − state cū, while this interpretation is favored by Refs. [15,[26][27][28]. Then there are still four natural parity resonances, D * (2600) 0 , D * (2650) 0 , D * 1 (2680) 0 and D * 1 (2760) 0 in the mass range of 2.6 2.8 GeV. In the traditional conventions of charm meson spectroscopy, these four resonances should correspond the 2 3 S 1 and 1 3 D 1 states cū or the mixtures of them.
These newly observed charm resonances have also been studied with the 2 3 S 1 , 1 3 D 1 assignments or the 2S-1D mixing scheme in theory by several models, including the non-relativistic quark model [21,27], the heavy quark effective theory [28,29], the effective Lagrangian approach based on heavy quark chiral symmetry [30], the EHQ decay formula [15,24] and the QPC model [23,[31][32][33][34]. However, the current theoretical calculations for these higher mass charmed mesons can not be well consistent with the experimental data. We find the calculated ratio R D + [D * (2600) 0 ] for taking it as the 2 3 S 1 state is usually greater than the experimental value Eq. (3) [27][28][29][30]32], while Ref. [30] argues that no quantum number assignments for pure state at mass 2600 MeV is able to reproduce the experimental ratio.
Generally, all the physical mesons have definite J P spin-parity or J P C for quarkonia. In the relativistic situations, the spin S and orbital angular momentum L are no longer the good quantum numbers, and the physical states are not always located in the definite 2S+1 L J states. This situations become obvious in the 1 + and 1 − mesons, for the 1 + states we always have to make the 1 P 1 -3 P 1 mixing to fit the physical states [35,36], while for the 1 − states the 2 3 S 1 -1 3 D 1 mixing is needed to fit the experimental measurements [37]. So in a more effective and appropriate method to describe the bound state, we should focus on the J P (C) , which are the good quantum numbers in any case. In principal, if we use a full relativistic method to solve the eigenstate problem of the bound mesons with definite J P (C) , we do not need mixing to fit the data. We have tried this by BS method in a previous work to study the state D * s1 (2700) [38], based on the BS wave function constructed directly from the quantum number J P = 1 − . The Salpeter wave functions of 1 − states were given, and by solving the full Salpeter equations, we obtain the eigenstates for cs and find that all the states include both S and D-wave components. The first state is 1S dominant while D-wave components can be ignored. The second state is 2S dominant, which is the first radial excited state. The third state, which is the second radial excited state, is predominant by 1D components. But our previous results, including the mass spectra and decays can not fit the data very well. The reason is that, we also make some approximations. The first is the instantaneous approximation, which assumes the potential is static, since the four-dimensional BS equation with non-static potential is quite difficult to solve. The second is the interaction kernel, where we choose the Coulomb-like plus linear potential. The Coulomb potential comes from the single-gluon-exchange, where we only keep the first order of QCD interaction. Also the linear confinement potential is introduced by phenomenological analysis. Since BS equation is a integral equation, then the kernel include all the ladder diagrams contributions but not the cross diagrams and the annihilation diagrams. These approximations have some effects in diagonalizing the mass matrix. So our method is not a full theory and not a full relativistic method, and can not exactly fit the experimental measurements. To overcome this discrepancy, we will make a further mixing to fit the physical states. In this study, we will give a continuous study of these 1 − states open charm mesons. We make a further mixing by the second and third radial excited states. Our mixing angle may be smaller than other non-relativistic methods since some relativistic corrections have already been kept in.
The manuscript is organized as below: In Section II we give the theoretical formalisms of the strong decays by BS methods; then in Section III the numerical results and detailed discusses are present; finally, we give a brief summary and conclusion about this work.

II. Theoretic calculations
In this section first we give a brief review on the calculations of transition matrix element and BS methods; then the 1 − states Salpeter wave functions are present. The Feynman diagram for strong decays of charmed meson is showed in Fig. 1, where we use subscript 1 and 2 to denote the final charmed meson and light meson, respectively. By using the reduction formula, the transition matrix element for decay D * s → D ( * ) K can be written as [41],

II.1. Transition matrix element
where P , P 1 and P 2 denote the momenta of initial state D * s , final charmed meson D ( * ) and final light meson K, respectively (see Fig. 1); M 2 is the mass of final light meson. Φ 2 (x) is used to describe the light scalar meson field. The PCAC relation reads where F 2 is the decay constant of the light scalar meson; q = u or d corresponds to the K + and K 0 respectively; the abbreviation Γ µ = γ µ γ 5 is used. Inserting the PCAC relation into Eq. (4), with the low energy theorem, Eq. (4) can be expressed as Then the decay amplitude M can be described as, where the transition matrix element D ( * ) (P 1 )|sΓ µ q|D * s (P ) can be calculated by Salpeter method and will be derived in next subsection. The decay width Γ is then expressed as, where is used; |M| 2 stands for the average over initial spins and sum over final spins; When the light meson is η, the η − η ′ mixing should be considered. In this work we use the following mixing conventions, η 8 and η 1 are the SU(3) octet and singlet states, respectively. We use the mixing angle θ η = 19 • . To include this mixing effect, the PCAC relation reads where in the last step, we have only remained thesΓ µ s part since others have no contribution here; f η 8 and f η 1 are the corresponding decay constants of η 8 and η 1 , respectively. When the π 0 is involved in the final states, the PCAC relation reads Again we have only kept the contributory part.

II.2. Transition matrix element with Salpeter wave function
In this subsection we briefly review the BS methods. The BS equation is an four-dimensional integral equation, which reads in momentum space as [39] ( where Ψ(q) is the four dimensional BS wave function; V (q−k) stands for the BS interaction kernel; p 1 and p 2 are the quark and anti-quark momentum respectively, while m 1 and m 2 are the corresponding masses (see Fig. 1). It is more convenient to express the p 1 and p 2 with the total momentum P and inner relative momentum q as where q P = P ·q M and q ⊥ = q − P M q P , in rest frame of initial meson they correspond to the q 0 and q respectively; the 3-dimensional integration η(q ⊥ ) can be understood as the BS vertex for bound states; V (|q ⊥ − k ⊥ |) denotes the instantaneous interaction kernel, namely, the inner interaction are assumed to be a static potential. As usual, in this work, the specific interaction kernel V (r) we use are the Coulomb-like potential plus the unquenched scalar confinement one [47], where λ is the string constant, α s (r) is the running strong coupling constant, and V 0 is a free constant fixed by fitting the data. By Fourier transformation, the potential V ( q ) in momentum space reads, where the running coupling constant α s ( q ) = 12π 27 1 log(a+ q 2 /Λ 2 QCD ) . Then under the instantaneous approximation, the BS equation (12) can be written as S(p 1 ) and S(−p 2 ) are the propagators for the quark and anti-quark respectively, and can be decomposed as where are the projection operators, which have the following forms, It can be easily check that, the projection operators satisfy the following relations: Since the BS kernel is assumed to be instantaneous, we can perform a contour integration over q P on both sides of Eq. (17), then we achieve the Salpeter equation as To make further simplification, we introduce four new wave functions ϕ ±± (q ⊥ ) with the definitions as where ϕ ++ is then called the positive Salpeter wave function, while ϕ −− is called the negative Salpeter wave function.
Then with the help of Eqs. (20), the Salpeter equation (21) can be further expressed as the following 4 coupled equations [40] From above equations, we can see that in the weak binding condition, namely, M ∼ (ω 1 + ω 2 ), ϕ −− is much smaller compared with ϕ ++ and can be ignored in the calculations. However, these four equations play equivalent roles in solving the eigenstate problem. The normalization condition for Salpeter wave function reads According to Mandelstam formalism [50], the transition matrix element D ( * ) (P 1 )|sΓ µ q|D * s (P ) can be expressed as whereφ ′++ P 1 is defined as γ 0 (ϕ ′++ P 1 ) † γ 0 , and ϕ ′++ P 1 is the positive Salpeter wave function of the final state; and m ′ 2 are the constituent quark and anti-quark masses in the final charmed meson (see Fig. 1). To achieve a final result, we still need to know the specific form of the corresponding Salpeter wave functions.

II.3. Salpeter wave function
The Salpeter wave functions involved in this calculations include the 0 − , 1 − and 1 + states, which corresponds to the 1 S 0 , 3 S 1 ( 3 D 1 ) and 1 P 1 ( 3 P 1 ) states within the non-relativistic models. The Salpeter wave function of 0 − state can be seen in Ref. [51]. Here we only give the 1 − state Salpeter wave function. The 3 S 1 and 3 D 1 states share the same spin-parity J P = 1 − . We rewrite the Salpeter wave function for 1 − states [47] as below, where f i (i = 1, 2, · · · , 8) are the radial wave functions; ǫ µP q ⊥ ξ = ǫ µναβ P ν q α ⊥ ξ β and ǫ µναβ is the totally antisymmetric Levi-Civita tensor; ξ denotes the polarization vector for initial state and fulfills P · ξ = 0, ξ By using Salpeter equations (23), we obtain the following 4 constraint conditions, which left us 4 independent wave functions f 3 , f 4 , f 5 and f 6 , only depending on | q | directly. In the paper, ω i is defined as m 2 i + q 2 (i = 1, 2). It can be easily check that, with above Salpeter wave function form, every item in Eq. (28) has the same quantum number J P = 1 − . Noticed that this wave function form and constraint conditions for 1 − state are not exactly the same with that in Ref. [47], however, it can be proved that the two forms are totally equivalent.
According to the definitions Eq. (22), the positive Salpter wave function ϕ ++ for 1 − state is then expressed as And the corresponding coefficients A i are The negative Salpeter wave function ϕ −− can be obtained similarly or by ϕ −− = ϕ − ϕ ++ . Then inserting the expressions of ϕ ++ and ϕ −− into the coupled Salpeter equations (24) and (25), we achieve the radial eigenvalue equations, which can be solved numerically. The normalization condition for 1 − states Salpeter wave functions now becomes, Interested readers can see a more detailed procedures on solving the full Salpeter equations in our previous works [45,[51][52][53]. By solving the Salpeter equations, finally we achieve these eight radial wave functions numerically, which are showed in Fig. 2. Fig. 2(a) shows the 8 radial wave functions of the first radial excited state, and Fig. 2(b) shows the radial wave functions of the second radial excited state. From the two diagrams, also considering that in Eq. (28), the direction of momentum / q ⊥ has contribution to the S or D wave [54], we can conclude that both the first and second radial excited states have S and D wave components, while the first radial excited state is 2 3 S 1 predominant and the second radial excited state is a 1 3 D 1 dominant state as has been stated in Ref. [55].
Since the decay final states include 1 + meson, we show our treatment of their wave functions. Generally, the bound state mesons consisting of unequal masses of quark and anti-quark do not have definite charge conjugation parity. So the physical 1 + states D 1 and D ′ 1 can be considered as the admixtures of 1 ++ ( 3 P 1 ) and 1 +− ( 1 P 1 ) states. Here we will follow the mixing conventions in  Refs. [56,57], where the mixing form for 1 + states is defined by the mixing angle α 1P as The heavy quark effective predicts that, in the limit m Q → ∞ the mixing angle for 1 + states are expressed as α 1P = arctan 1/2 = 35.3 • . This result will be used in the strong decay calculations when D (′) 1 mesons are involved in the final states. The Salpeter wave functions for 1 +− and 1 ++ states can be found in Ref. [58].
Having these numerical Salpeter wave functions, we can calculate the 3-dimensional integral of the transition matrix element D ( * ) (P 1 )|sΓ µ q|D * s (P ) in Eq. (27). The detailed information on performing this integral can be found in our previous work Refs. [43,46].

III. Numerical Results and Discussions
First we specify the corresponding parameters used in this work. The constituent quark masses and other parameters to characterize the model are before [43] a = e = 2. The free parameter V 0 is fixed by fitting the mass eigenvalue to experimental value. The decay constants we used are f π = 130.4 MeV, f K = 156 MeV [7], f η 8 = 1.26f π and f η 1 = 1.07f π . The mixing angle θ between η − η ′ we choose is θ η = 19 • with the mixing convention in Eq. (9). Other involved parameters are from PDG data [7] unless otherwise specified.

III.1. 1 − charm-strange mesons
Both the D * s1 (2700) and D * s1 (2860) share spin-parity J P = 1 − determined by experiments. In the first place, we take them as the pure first and second radial excited states, which are dominant by 2 3 S 1 and 1 3 D 1 respectively. So in this work we still label the first radial excited state as 2 3 S 1 and the second excited state as 1 3 D 1 . The calculated main strong decays properties are listed in Tab. II. From Tab. II we can see that, for D * s1 (2700), neither the 2 3 S 1 nor 1 3 D 1 assignment can produce the Tab. II: Decay widths of D * s1 (2710) and D * s1 (2860) as the 2 3 S 1 and 1 3 D 1 dominant cs states in MeV.  (2860) respectively are roughly consistent with the experimental measurements. Also notice that the total width is only the half of experimental value when taking D * s1 (2700) as the 1 3 D 1 state. So these assignments can not produce experimental data. One also notes that the predicted total decay widths in this paper are much larger than our previous calculation [38], the reason is that we have chosen different D * s1 (2700) mass as input, besides the difference of phase space, the node structure of 2S state also has sensitive effect due to the variance of phase space.
Then we introduce the further 2 3 S 1 -1 3 D 1 mixing scheme in the 1 − charm-strange system. The mixing form we used is defined as, and the mixing states D * sa and D * sb correspond to D * s1 (2700) and D * s1 (2860), respectively. To get the experimental branching ratio R K [D * s1 (2700)] = 0.91 −0.18 +0. 18 , we obtain the mixing angle θ s = (8.7 +3.9 −3.2 ) • . The errors in experimental ratio R K [D * s1 (2700)] are combined in quadrature for simplicity. The uncertainties in our mixing angle θ s and other obtained results are induced by varying the experimental ratio R K [D * s1 (2700)] in 1σ range. The decay properties with 2S-1D mixing are listed in Tab. III. The total width for D * s1 (2700) is about 100.8 MeV, which agrees well with the experimental measurement Γ = 117 ± 13 MeV [7]. Our results are also consistent with that in Ref. [23], where the mixing angle is about 6.8 • 11.2 • and the calculated Γ[D * s1 (2700)] is about 100 MeV. With the obtained mixing angle, the total width for D * s1 (2860) we obtain is 108.8 MeV, which is also comparable with the LHCb result Γ[D * s (2860)] = 159 ± 80.3 MeV, but less than the result 300 MeV in Ref. [23]. Furthermore, the predicted ratio R K [D * s1 (2860)] = 0.62, which is also consistent with the result 0.6 0.8 in Refs. [17,23,24] and could be used to test this 2 3 S 1 -1 3 D 1 mixing scheme in the future measurements. The comparisons of our results with other predictions can be seen in Tab. IV.
If we do not restrict the mass of 2 3 S 1 state is less than that of 1 3 D 1 state, we can obtain another large mixing angle θ s ≃ −77 • , which could also reproduce the experiment ratio R K [D * s1 (2710)] = 0.91. The strong decay properties are also listed in Tab. III. In such case, the Γ[D * s1 (2700)] is 63 MeV, which is about the half of experimental value; the total width for its orthogonal partner D * s1 (2860) is about 150 MeV; the ratio R K [D * s1 (2860)] is 1.09. Refs. [19][20][21] also achieved a large mixing angle −(57 77) • . We noticed that in Ref. [18] the large mixing angle is also obtained, although in the Tab. III: Decay widths in MeV for D * s1 (2700) and D * s1 (2860) under the further 2 3 S 1 -1 3 D 1 mixing. The mixing angle θ s is in unit of degree.  Tab. IV: Comparisons with other Refs. when taking D * s1 (2700) and D * s1 (2860) as the mixtures of 2 3 S 1 -1 3 D 1 cs. Decay width Γ is in unit of MeV and the mixing angle θ s is in unit of degree.
Above all, we calculate the strong decays properties by taking all these four resonances as the 2 3 S 1 or 1 3 D 1 state cū. The obtained results are listed in Tab. V. We can see that, when taking D * (2600) as the 2 3 S 1 state, both the total widths and ratio R D + [D * (2600)] are comparable with the BaBar Tab. V: Decay properties of D * (2600), D * (2650), D * 1 (2680) and D * 1 (2760) as the 2 3 S 1 or 1 3 D 1 dominant cū states in unit of MeV.  [3]. Since for the other three resonances, only the total widths can be used to compare with experiments. From the calculated total widths, we can only make rough judgments, both the 2 3 S 1 and 1 3 D 1 assignments seem reasonable for D * (2650) and D * 1 (2680), while also considering the mass predictions [11,12], they are more likely to be the 2 3 S 1 states. Taking D * 1 (2760) as 1 3 D 1 state cū gives the width 290 MeV, which is about 100 MeV larger than the LHCb measurement 180 MeV [9]. We also find that the decay channel D + 1 π − becomes quite important in the decay of 1 3 D 1 state, hence we define ratio . This ratio is quite sensitive to the assignments of 2 3 S 1 or 1 3 D 1 . All in all, for identification of these excited 1 − resonances, the consistent measurements from experiments are necessary and pivotal. D * (2600) seems consistent with the 2 3 S 1 assignment, while the predicted ratio Γ(D + π − ) Γ(D * + π − ) = 0.54 is a little larger than BaBar measurement 0.32 [3]. This small discrepancy between theoretical and experimental results hints, there exists a small mixing between the 2 3 S 1 and 1 3 D 1 states. The physical quantity R D + [D * (2600)] can behave as a good restriction to the mixing angle, just as what we have done in the 1 − charm-strange systems. Then again we introduce the 2 3 S 1 -1 3 D 1 mixing scheme as At first, we take D * (2600) as the 1 − state cū dominant by 2 3 S 1 components, while D * (2650) as the orthogonal partner of D * (2600). To fix the ratio R D + [D * (2600)] at BaBar's measurement 0.32 −0.09 +0.09 [3], we obtain the mixing angle θ u = −(7.5 +4.0 −3.3 ) • . The theoretical uncertainties are induced by varying the experimental ratio R D + [D * (2600)] in 1σ range of its central value. Our results reveal that the mixing angle θ u is not sensitive to the mass of D * b . When m D * b ranges from 2.65 to 2.78 GeV, the variation of mixing angle is about 0.1 degree. So we will ignore this tiny difference in the following statements. The partial decay widths are listed in Tab. VI, where D * (2650), D * (2680) or D * 1 (2760) is taken as the orthogonal partner of D * (2600). The dependence of Γ D * b and ratio R D + (D * b ) over the mass of D * b can be seen in Fig. 3, where we let m D * b range from 2.65 to 2.78 GeV. It can be seen clearly in Fig. 3(a) that the corresponding Γ D * b increases from about 140 MeV to 290 MeV. The predicted ratio R D + [D * b ] goes down from 9.4 to 6.8, which is displayed in Fig. 3(b). The calculated total width Γ D * (2600) = 66 MeV is comparable with BaBar's measurement 93 MeV [3], while Γ D *     Of course, we can still take D * 1 (2680) as the |D * a state while D * 1 (2760) as the |D * b state, the results should behave similar with above tests, the mixing angle will be even smaller, and the corresponding properties will behave almost the same with that under the assignments without this further mixing.
The comparisons of our results with others can be found in Tab 4. D * (2650) is more likely to be the 2 3 S 1 predominant state. There is no great conflict in the total widths of 2 3 S 1 and 1 3 D 1 assignments, while its mass is more consistent with 2 3 S 1 state. The ratio Γ(D + π − ) Γ(D * + π − ) behaves quite differently for the two assignments, which is 0.5 for 2 3 S 1 assignment and 4.4 for 1 3 D 1 assignment. Hence this ratio can be used to test the essence of D * (2650). 5. For D * 1 (2680), the situation is similar with D * (2650). There exists no great conflict in the total widths between 2 3 S 1 and 1 3 D 1 assignments compared with experimental measurements. The ratios Γ(D + π − ) Γ(D * + π − ) are 0.46 and 4.15 for 2 3 S 1 and 1 3 D 1 assignment, respectively, therefore can be used to discriminate the essence of D * 1 (2680).
Tab. VIII: Comparison with other Refs. when taking D * (2600) as the mixture of 2 3 S 1 -1 3 D 1 cū. Γ Tot is in unit of MeV and the mixing angle θ u is in unit of degree.

IV. Summary
In this work, we carried out a systematical research on the potential 1 − open charm mesons, including the charm-strange mesons D * s1 (2700) and D * s1 (2860), charm mesons D * (2600), D * (2650), D * 1 (2680) and D * 1 (2760). The main strong decay properties by taking these natural spin-parity resonances as the 2 3 S 1 or 1 3 D 1 states are achieved by using the Bethe-Salpeter methods. In particularly, the further 2S-1D mixing scheme is used to explain both the 1 − charm and charm-strange mesons. The obtained results and predicted properties can be tested in the near future experiments.
Our results reveal that, D * s1 (2700) and D * s1 (2860) can be well described by the further 2 3 S 1 -1 3 D 1 mixing scheme with a small mixing angle (8.7 +3.9 −3.2 ) • . Both the total widths and ratio of corresponding partial decay widths are consistent with the experimental measurements. Our predicted ratio Γ(D * K) Γ(D * K) for D * s1 (2860) is 0.62 +0.22 −0.12 , which could be used to test this 2S-1D mixing scheme in future. For the corresponding charm mesons, since the experimental measurements are not consistent with each other, the identification and assignments are much more difficult. Based on our results, the BaBar D * (2600) [3] can be explained by the same mixing scheme with a mixing angle of −(7.5 +4.0 −3.3 ) • . D * (2650) [5] and D * 1 (2680) [10] are more likely to be 2 3 S 1 predominant states, since their masses are consistent with the 2 3 S 1 predictions, while our calculated total widths are both comparable with the experimental measurements under the 2 3 S 1 or 1 3 D 1 assignments. Our results also show that, the measured width of D * 1 (2760) is much smaller than the theoretical calculations under the 1 3 D 1 assignment. This would be an obstacle to identify D * 1 (2760) as the 1 3 D 1 predominant cū. There still exit puzzles and difficulties in identifications of these new excited charmed mesons. Further precise measurements of their properties are needed and important.