Scientific Research

An Academic Publisher

**On Radiative Decay of Heavy Vector Mesons V → Pγ in the Scalar Strong Interaction Hadron Theory** ()

^{}

*V*→

*P*

*γ*is treated semi-classically in the scalar strong interaction hadron theory. The treatment successfully employs the new wave function of the vector meson updated here. The ratio of the available

*J/*

*ψ*and decay rates agrees with prediction. The values of the predicted rates are also in order of magnitude agreement with measurements. These agreements are the only ones directly computed from a first principles’ theory.

Share and Cite:

*V*→

*Pγ*in the Scalar Strong Interaction Hadron Theory.

*Journal of Modern Physics*,

**9**, 584-595. doi: 10.4236/jmp.2018.94040.

1. Introduction

Radiative decay of heavy ground state vector meson V → Pγ has been treated earlier [1] and Section 6.4 of [2] , hereafter denoted by I, in the scalar strong interaction hadron theory [3] . The dimensional approximation used was unsatisfactory and the results disagreed with later data [4] . Further, the meson wave functions used in these earlier works have recently been updated [5] . The purpose of this paper is to provide a revised treatment of the decay V → Pγ taking into account these two developments as well as including a new gauge field component ignored earlier.

This paper is divided into two following sections. In Sec. 2, the wave function for pseudoscalar mesons at rest I (4.3.2) is replaced by [5 (8a)] and that of the vector meson I (4.3.3) by a corrected wave function in (1.6) below. The wave functions of heavy, slowly moving pseudoscalar mesons treated in I Section 3.5 are replaced by new ones in §2.3 below. These results are then used in the new approach in Sec. 3 to evaluate the V → Pγ rates semi-classically.

2. Meson Wave Functions

The scalar strong interaction hadron theory allows inherently for two types of quark-antiquark interactions, a Coulomb type and a harmonic oscillator type [3, §7], I (3.2.8). When the theory was initially being developed in the early 1990’s, data suggested that the confining potential is of Coulomb plus linear type [6] [7] . Since a linear confinement at large quark separations is inherent in the theory [3, §7], I (3.2.19) the Coulomb type of potential was, somewhat indiscriminatingly, chosen.

Such a potential could rather successfully account for ground state meson spectra [8] , I Section 5.3-4, but failed to account for the excited meson spectra I Section 5.5-7. This was largely remedied by replacing the Coulomb potential by the harmonic oscillator type [5] . Here, a correction to the ground state vector meson wave function in [5] will be made.

In §2.1, the updated ground state meson wave function is given. The wave equations of a nonrelativistic pseudoscalar meson are given in §2.2 and the wave functions for heavy mesons are derived in §2.3.

2.1. Wave Functions of Ground State Mesons at Rest

For ground state singlet mesons, I (3.4.1) or [5 (4)] with J = 0 holds. The updated and normalized pseudoscalar meson wave function using the harmonic type of potential reads [5 (8b) with J = 0].

${\psi}_{00}\left(r\right)=\frac{1}{\sqrt{\Omega}}{\alpha}_{00}\mathrm{exp}\left(-\frac{{d}_{h}}{2}{r}^{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\alpha}_{00}={\left(\frac{{d}_{h}}{\text{\pi}}\right)}^{3/4}=0.0577\text{\hspace{0.17em}}{\text{GeV}}^{3/2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}J=0$ (1.1)

in which d_{h} = 0.07 GeV^{2} of [5 (11)] has been used.

For ground state triplet mesons, I (3.4.3) or [5 (4)] with J = 1 derived from I (3.2.11a) is incorrect in two signs and is replaced by

$\left(\frac{{\partial}^{2}}{\partial {r}^{2}}-\frac{2}{r}\frac{\partial}{\partial r}+\frac{2}{{r}^{2}}+{d}_{m0}-{d}_{h}^{2}{r}^{2}+\frac{{E}_{10}^{2}}{4}-{M}_{m}^{2}\right){\psi}_{1}\left(r\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\psi}_{1}\left(r\right)={\chi}_{1}\left(r\right)$ (1.2)

The nonlinear potential Φ_{c1} there has been dropped in the linearized (1.2) which, like (1.1), has a harmonic oscillator type of solution [5 (8a)]

${\psi}_{10}\left(r\right)=\frac{1}{\sqrt{\Omega}}\mathrm{exp}\left(-\frac{{d}_{h}}{2}{r}^{2}\right){\displaystyle \underset{\nu =0}{\overset{n}{\sum}}{a}_{\nu}{r}^{s+\nu}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{d}_{h}>0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{a}_{0}\ne 0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{a}_{odd}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=2{n}_{r}$ (1.3)

where n_{r} is the radial quantum number. An extra subscript 0 has been added to indicate that ψ_{10} is of zeroth order in later perturbational calculations. Following [5 (9)], the series terminates when

${a}_{n+2}=\frac{2{d}_{h}\left(s+n-\frac{1}{2}\right)-\frac{1}{4}\left({E}_{10}^{2}-{\left({m}_{p}+{m}_{r}\right)}^{2}\right)-{d}_{m0}}{\left(s+n+2\right)\left(s+n-1\right)+2}{a}_{n}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}s=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}\text{\hspace{0.17em}}2$ (1.4)

These gives the vector meson mass E_{10},

$\frac{1}{4}{E}_{10}^{2}=\frac{1}{4}{\left({m}_{p}+{m}_{r}\right)}^{2}-{d}_{m0}+5{d}_{h},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}s=1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=2$ (1.5)

Here, the lowest allowed s and n values have been chosen. The choice n = 0 would lead to a mass E_{10} < E_{00}, the mass of the corresponding pseudoscalar meson, contrary to data. Therefore, the next lowest value n = 2 is chosen. This relation is the same as [5 (10)] with J = 1 and n = 0 for vector meson so that the spectra found in [5] remain unchanged. The only difference is that the wave function [5 (8b) with J = 1] is replaced by the updated, normalized wave function for a vector meson obtained from (1.3)

${\psi}_{10}\left(r\right)=\frac{1}{\sqrt{\Omega}}{\alpha}_{10}r\left(1-2{d}_{h}{r}^{2}\right)\mathrm{exp}\left(-\frac{{d}_{h}}{2}{r}^{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\alpha}_{10}=2.44\times {10}^{-3}\text{\hspace{0.17em}}{\text{GeV}}^{5/2}$ (1.6)

2.2. Wave Equations for Nonrelativistic Pseudoscalar Mesons

The starting point is I §3.5.1 which introduces the small parameter

${\epsilon}_{0}=\left|\underset{\_}{K}\right|/{E}_{00}\ll 1$ (1.7)

where K is the momentum of the nonrelativistic pseudoscalar meson. Expand the meson energy E_{K} in I (3.1.6) and the wave functions in I (3.1.7a), suppressing the relative time factor in I (3.1.9), in the form

${\psi}_{0}\left(\underset{\_}{x}\right)={\displaystyle \underset{i}{\sum}{\psi}_{0i}\left(\underset{\_}{x}\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{\_}{\psi}\left(\underset{\_}{x}\right)={\displaystyle \underset{i}{\sum}{\underset{\_}{\psi}}_{i}\left(\underset{\_}{x}\right)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\psi \to \chi ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{E}_{K}={\displaystyle \underset{i}{\sum}{E}_{J\text{\hspace{0.17em}}i}}$ (1.8)

where i denotes the ith order in ε_{0}. Only freely moving mesons are considered so that the nonlinear potential Φ_{c} in I (3.2.8) drops out according to (1.1, 3) where Ω → ¥. Therefore, Φ_{m} of I (3.2.8a) is independent of the meson wave functions (1.1, 3).

To zeroth order in ε_{0}, I (3.5.8) goes over to I (3.2.10b) and

${\chi}_{00}\left(r\right)=-{\psi}_{00}\left(r\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\underset{\_}{\chi}}_{0}\left(r\right)={\underset{\_}{\psi}}_{0}\left(r\right)=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{E}_{K}={E}_{0}={E}_{00}$ (1.9)

To first order in ε_{0}, I (3.5.8), (3.2.8) with the potentials replaced by those in (1.2) leads to

${\chi}_{01}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\psi}_{01}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{E}_{01}=0$ (1.10a)

$\begin{array}{l}\left({E}_{00}^{2}/4-\Delta \right){\underset{\_}{\chi}}_{1}+2\underset{\_}{\partial}\left(\underset{\_}{\partial}{\underset{\_}{\chi}}_{1}\right)+{E}_{00}\underset{\_}{\partial}\times {\underset{\_}{\chi}}_{1}+\left({d}_{m0}-{d}_{h}^{2}{r}^{2}-{M}_{m}^{2}\right){\underset{\_}{\psi}}_{1}\\ =\frac{1}{2}{E}_{00}\underset{\_}{K}{\chi}_{00}-\underset{\_}{K}\times \underset{\_}{\partial}{\chi}_{00}\end{array}$ (1.10b)

$\begin{array}{l}\left({E}_{00}^{2}/4-\Delta \right){\underset{\_}{\psi}}_{1}+2\underset{\_}{\partial}\left(\underset{\_}{\partial}{\underset{\_}{\psi}}_{1}\right)-{E}_{00}\underset{\_}{\partial}\times {\underset{\_}{\psi}}_{1}+\left({d}_{m0}-{d}_{h}^{2}{r}^{2}-{M}_{m}^{2}\right){\underset{\_}{\chi}}_{1}\\ =\frac{1}{2}{E}_{00}\underset{\_}{K}{\psi}_{00}+\underset{\_}{K}\times \underset{\_}{\partial}{\psi}_{00}\end{array}$ (1.10c)

To second order in ε_{0}, the singlet part of I (3.5.8) becomes

$\begin{array}{l}\left({E}_{00}^{2}/4+\Delta \right){\chi}_{02}-\left({d}_{m0}-{d}_{h}^{2}{r}^{2}-{M}_{m}^{2}\right){\psi}_{02}\\ =-\underset{\_}{\partial}\left(\underset{\_}{K}\times {\underset{\_}{\chi}}_{1}\right)-\left(\frac{1}{2}{E}_{00}{E}_{02}+\frac{1}{4}{\underset{\_}{K}}^{2}\right){\chi}_{00}+\frac{1}{2}{E}_{00}\underset{\_}{K}{\underset{\_}{\chi}}_{1}\end{array}$ (1.11a)

$\begin{array}{l}\left({E}_{00}^{2}/4+\Delta \right){\psi}_{02}-\left({d}_{m0}-{d}_{h}^{2}{r}^{2}-{M}_{m}^{2}\right){\chi}_{02}\\ =\underset{\_}{\partial}\left(\underset{\_}{K}\times {\underset{\_}{\psi}}_{1}\right)-\left(\frac{1}{2}{E}_{00}{E}_{02}+\frac{1}{4}{\underset{\_}{K}}^{2}\right){\psi}_{00}+\frac{1}{2}{E}_{00}\underset{\_}{K}{\underset{\_}{\psi}}_{1}\end{array}$ (1.11b)

The spherical symmetry present in the ε_{0} = 0 limit is broken by the momentum K so that separation of variables in the relative space x cannot be carried out. This renders that each of (1.10, 11) consists of two coupled second order partial differential equations containing eight dependent variables and cannot be readily solved analytically.

2.3. Approximative Heavy Meson Wave Functions

For heavy mesons, the last term is small next to the first term on the right sides of (1.10b, 10c) when

${E}_{00}\gg 2\sqrt{{d}_{h}}\approx 0.529\text{\hspace{0.17em}}\text{GeV}$ (1.12)

Here, the magnitude of
$\underset{\_}{\partial}$ in (1.10) is d_{h}r according to (1.1) where r has been replaced by some mean value
${r}_{0}\approx 1/\sqrt{{d}_{h}}$ . This inequality holds roughly for mesons containing b or c quarks but not for kaon and pion. In these cases, an approximate solution can be found. Let

$\underset{\_}{K}=\left(0,0,K\right)$ (1.13)

The approximation consists of making the simplifying ansatz

${\underset{\_}{\chi}}_{1}\left(\underset{\_}{x}\right)=-{\underset{\_}{\psi}}_{1}\left(\underset{\_}{x}\right)$ (1.14)

in (1.10b, 10c). Addition and subtraction of thee two equations yields

$\left({E}_{00}^{2}/4-\Delta +{d}_{h}^{2}{r}^{2}-{d}_{m0}+{M}_{m}^{2}\right){\underset{\_}{\psi}}_{1}=\frac{1}{2}{E}_{00}\underset{\_}{K}{\psi}_{00}-2\underset{\_}{\partial}\left(\underset{\_}{\partial}{\underset{\_}{\psi}}_{1}\right)$ (1.15a)

$-{E}_{00}\underset{\_}{\partial}\times {\underset{\_}{\psi}}_{1}=\underset{\_}{K}\times \underset{\_}{\partial}{\psi}_{00}$ (1.15b)

The ansatz (1.14) thus leads to that these two equations determine only one unknown ${\underset{\_}{\psi}}_{1}$ ; there has to be an inconsistency. Solving (1.15b) with (1.13) gives

${\underset{\_}{\psi}}_{1}\left(\underset{\_}{x}\right)=\frac{\underset{\_}{K}}{{E}_{00}}{\psi}_{00}\left(r\right)=\left(0,0,{\psi}_{1z0}\left(\underset{\_}{x}\right)\right)$ (1.16)

which shows that
${\underset{\_}{\psi}}_{1}$ is of order ε_{0}. Inserting this expression into (1.15a) using (1.1) gives

$\frac{1}{2}{E}_{00}K{\psi}_{00}\left(r\right)=\frac{1}{2}{E}_{00}K{\psi}_{00}\left(r\right)+2\frac{K}{{E}_{00}}{d}_{h}\left(1-\stackrel{\u2322}{z}\stackrel{\u2322}{z}{d}_{h}{r}^{2}\right){\psi}_{00}\left(r\right)$ (1.17)

The last term comes from the last term in (1.15a) and causes that (1.17) cannot be satisfied and is inconsistent with (1.16).

However, the last term in (1.17) will be small for heavy mesons having large E_{00} so that (1.17) and the solution (1.16) both hold approximately. Replacing r by the mean value r_{0} above and putting
$\stackrel{\u2322}{z}\stackrel{\u2322}{z}=1/3$ in anticipation of later angular integration, the criterion is that the ratio between the two terms on the right side of (1.17) be small. This gives another requirement for heavy meson approximation

$\frac{8{d}_{h}}{3{E}_{00}^{2}}\approx \frac{0.187}{{E}_{00}^{2}}\ll 1$ (1.18)

This ration is 9.6 for pion, 0.762 for kaon, 0.054 for D, and 0.0067 for B meson. Thus, (1.18) is well satisfied by mesons with b or d quark but again not by kaon or pion. This inequality agrees roughly with that given by (1.12). There is however no great loss here; pions as decay products often move relativistically so that the first order (1.10) no longer holds in the first place.

The last terms in (1.15a) will also introduce a correction to (1.16). As an estimate, let the corrected (1.16) be an average of (1.16) and (1.17) multiplied by $2/{E}_{00}^{2}$ and generalized to include the x and y components of the last term in (1.15a);

${\underset{\_}{\psi}}_{1}\left(\underset{\_}{x}\right)=\left({\psi}_{1x1}\left(\underset{\_}{x}\right),{\psi}_{1y1}\left(\underset{\_}{x}\right),{\psi}_{1z0}\left(\underset{\_}{x}\right)+{\psi}_{1z1}\left(\underset{\_}{x}\right)\right)$ (1.19)

$\begin{array}{l}{\psi}_{1x1}\left(\underset{\_}{x}\right)=-\frac{2K}{{E}_{00}^{3}}\stackrel{\u2322}{x}\stackrel{\u2322}{z}{d}_{h}^{2}{r}^{2}{\psi}_{00}\left(r\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\psi}_{1y1}\left(\underset{\_}{x}\right)=-\frac{2K}{{E}_{00}^{3}}\stackrel{\u2322}{y}\stackrel{\u2322}{z}{d}_{h}^{2}{r}^{2}{\psi}_{00}\left(r\right)\\ {\psi}_{1z1}\left(\underset{\_}{x}\right)=\frac{2K}{{E}_{00}^{3}}{d}_{h}\left(1-\stackrel{\u2322}{z}\stackrel{\u2322}{z}{d}_{h}{r}^{2}\right){\psi}_{00}\left(r\right)\end{array}$ (1.20)

The second order Equations (1.11) are treated analogously. Subtraction and addition of (1.11a) and (1.11b) leads to

${\chi}_{02}\left(\underset{\_}{x}\right)=-{\psi}_{02}\left(\underset{\_}{x}\right)$ (1.21a)

$\left(\frac{{E}_{00}^{2}}{4}+\Delta -{d}_{h}^{2}{r}^{2}+{d}_{m0}-{M}_{m}^{2}\right){\psi}_{02}=-\left(\frac{{E}_{00}{E}_{02}}{2}+\frac{{K}^{2}}{4}\right){\psi}_{00}+\frac{{E}_{00}\underset{\_}{K}{\underset{\_}{\psi}}_{1}}{2}$ (1.21b)

$\underset{\_}{\partial}\left(\underset{\_}{K}\times {\underset{\_}{\psi}}_{1}\right)=0$ (1.21c)

The last equation is satisfied for
${\underset{\_}{\psi}}_{1}$ given by (1.16). It is slightly violated by (1.19, 20) but the violations contain odd powers of the angles
$\stackrel{\u2322}{r}$ I (3.1.7b) and will vanish upon integration over the angles later. The left operator of (1.21b) is the same as the linearized ones in I (3.4.1) and vanishes for ψ_{02} µ ψ_{00}; ψ_{02} can be absorbed into the zeroth order ψ_{00} and put to 0. Thus,

$\frac{1}{2}{E}_{00}\underset{\_}{K}{\underset{\_}{\psi}}_{1}=-\left(\frac{1}{2}{E}_{00}{E}_{02}+\frac{1}{4}{K}^{2}\right){\psi}_{00},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\psi}_{02}={\chi}_{02}=0$ (1.22)

Inserting (1.16) into the first of (1.22) leads to

${E}_{02}={K}^{2}/2{E}_{00}$ (1.23)

3. Radiative Decay of Heavy Vector Meson V à Pγ

In this section, the wave functions found in Sec. 2 above are applied to the decay of a heavy ground state vector meson V into the corresponding pseudoscalar meson P and a photon γ. The treatment is semi-classical; the electromagnetic field is not quantized. This is in accord with the quantum mechanical nature of the present theory, which cannot be quantized.

Such decays have been treated earlier [1] and I Section 6.4. The estimated D*^{±} decay rate turned out to be too small compared to subsequent measurement. This section follows basically [1] with two main differences. Firstly, the assumption that the U(1) gauge field is limited to the two component photon field is removed. Secondly, the Coulomb form of wave functions of the ground state mesons adopted in [1] and I (4.3.2-3) is replaced by harmonic oscillator form (1.1.6). Consistently, the approximative wave functions of I §3.5.3 are replaced by those in §2.3 above.

Eugene Wigner once said: “Once the equation of motion is known, the rest is engineering”. In this sense, the present treatment may be such an “engineering” and assumes no “model”, as in the literature [9] .

The formalism of I Section 6.4 is largely taken over here.

3.1. Wave Functions of Decaying Meson

For a free meson, the meson Equations I (2.3.22) with Φ_{c} = 0 in I (3.2.8a) hold. If an electromagnetic field A is introduced on the quark level, I (2.3.22) is replaced by I (6.1.10). The magnitude of the difference between these two sets of equations is small, of the order of quark charges. Therefore, the associated decay V → Pγ can be formulated as a first order perturbational problem. The wave function of the decaying meson is taken to be a modified form of I (3.1.5, 6, 9) with I (3, 5.6),

$\begin{array}{l}{\psi}^{a\stackrel{\dot{}}{b}}\left(X,x\right)={\displaystyle \underset{K}{\sum}{b}_{JK}\left({\delta}^{a\stackrel{\dot{}}{b}}{\psi}_{JK}\left(\underset{\_}{x}\right)-{\underset{\_}{\sigma}}^{a\stackrel{\dot{}}{b}}{\underset{\_}{\psi}}_{JK}\left(\underset{\_}{x}\right)\right)\mathrm{exp}\left(-i{E}_{JK}{X}^{0}+i\text{\hspace{0.05em}}{\underset{\_}{K}}_{J}\underset{\_}{X}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\psi \to \chi}\\ {b}_{JK}={a}_{JK}+{a}_{JK}^{\left(1\right)}\left({X}^{0}\right)\end{array}$ (2.1)

The subscript J = 0, 1 refers to pseudoscalar and vector mesons, respectively. a_{JK} is unity here but is in a quantized case in §3.4 below to be elevated to an annihilation operator annihilating an initial meson of spin J and momentum K.
${a}_{JK}^{\left(1\right)}\left({X}^{0}\right)$ is a small first order amplitude that varies slowly with time and, in the quantized case in §3.4, becomes an operator that “slowly” transforms the initial vector meson to some intermediate state. It is zero at X^{0} = -¥. Similarly,
${a}_{JK}^{\ast}$ enters ψ*; the complex conjugate of ψ, is also unit and is to be elevated to a creation operator creating a final state with the same J and K.
${a}_{JK}^{\left(1\right)\ast}\left({X}^{0}\right)$ is the complex conjugate of
${a}_{JK}^{\left(1\right)}\left({X}^{0}\right)$ and, in the quantized case, becomes an operator that “slowly” creates the same final state as that created by
${a}_{JK}^{\ast}$ . These quantization assignments are phenomenological and rudimentary.

In V → Pγ, ${a}_{JK}^{\left(1\right)}\left({X}^{0}\right)$ is caused by the qA terms in (2.3) below and is hence of first order in quark charge. Picking out only one J value and one K value in the summations of (2.1), it leads to the simplified form

${\chi}^{a\stackrel{\dot{}}{b}}\left(X,x\right)={\chi}_{0}^{a\stackrel{\dot{}}{b}}+{\chi}_{1}^{a\stackrel{\dot{}}{b}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\chi}_{1}^{a\stackrel{\dot{}}{b}}=\left({a}_{JK}^{\left(1\right)}\left({X}^{0}\right)/{a}_{JK}\right){\chi}_{0}^{a\stackrel{\dot{}}{b}}$ (2.2)

where the subscripts 0 and 1 denote orders in quark charge. ${\chi}_{0}^{a\stackrel{\dot{}}{b}}$ is simply (2.1) with ${a}_{JK}^{\left(1\right)}\left({X}^{0}\right)=0$ .

3.2. U(1) Gauge Field

The photon is observable and the electromagnetic field is thus introduced on the meson level similar to that in [1 (2.1)] and to the U(1) gauge field in I (6.2.3-4). The equivalent of I (6.1.10a) becomes

$\begin{array}{l}\left({\partial}_{I}^{a\stackrel{\dot{}}{b}}+i\frac{1}{2}{q}_{p}{A}^{a\stackrel{\dot{}}{b}}\left(X\right)\right)\left({\partial}_{II}^{f\stackrel{\dot{}}{e}}-i\frac{1}{2}{q}_{r}{A}^{f\stackrel{\dot{}}{e}}\left(X\right)\right){\chi}_{\left(pr\right)\stackrel{\dot{}}{b}f}\left(X,x\right)\\ -\left({M}_{m}^{2}-{\Phi}_{m}\left(\underset{\_}{x}\right)\right){\psi}_{\left(pr\right)}^{a\stackrel{\dot{}}{e}}\left(X,x\right)=0\end{array}$ (2.3)

where I (3.1.3a) and the transition of I (2.3.23) to (3.1.11) have been used. Here,
${\partial}_{I}$ and
${\partial}_{II}$ operate on χ only. Putting the time component A_{0} to 0 and expand A(X) in plane waves, as in [1 (3.3)],

$\underset{\_}{A}\left(X\right)={\displaystyle {\sum}_{{K}_{\gamma}}\frac{1}{\sqrt{2{E}_{r}\Omega}}{\displaystyle {\sum}_{T}\underset{\_}{e}{a}_{T}\left({\underset{\_}{K}}_{\gamma}\right)\mathrm{exp}\left(-i{E}_{r}{X}^{0}+i{\underset{\_}{K}}_{r}\underset{\_}{X}\right)}}+c.c.\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{\_}{e}=\left({e}_{X},{e}_{Y},{e}_{Z}\right)$

(2.4)

where E_{γ} and K_{γ} are the energy and momentum of the photon and Ω a large normalization box. a_{T}(K) is the analog of a_{JK} in (2.1) and is set to a unity here but is elevated to an annihilation operator in the quantized case below.

In [1] and I, the unit vector
$\underset{\_}{e}=\left({e}_{T},{e}_{Z}=0\right)=\left({e}_{X},{e}_{Y},{e}_{Z}=0\right)$ for
${\underset{\_}{K}}_{\gamma}=\left(0,0,{K}_{\gamma}\right)$ representing the photon field associated with the both transverse modes T = 1, 2. The assumption is removed here and e_{Z} ≠ 0 is allowed here even if it is not part of the photon field.

3.3. First Order Relations

Inserting (2.2) into (2.3), multiplying it by ${\chi}_{\stackrel{\dot{}}{e}a}^{\ast}$ and integrating over X and x, the first order part reads

$\begin{array}{l}{\displaystyle \int {\text{d}}^{4}X{\text{d}}^{4}x{\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}[{\partial}_{I}^{a\stackrel{\dot{}}{b}}{\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{1\stackrel{\dot{}}{b}f}-\left({M}_{m}^{2}-{\Phi}_{m}\right){\psi}_{1}^{a\stackrel{\dot{}}{e}}+i\frac{1}{2}{q}_{p}{A}^{a\stackrel{\dot{}}{b}}\left(X\right){\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{0\stackrel{\dot{}}{b}f}}\\ -i\frac{1}{2}{q}_{r}{\partial}_{I}^{a\stackrel{\dot{}}{b}}{A}^{f\stackrel{\dot{}}{e}}\left(X\right){\chi}_{0\stackrel{\dot{}}{b}f}]=0\end{array}$ (2.5)

Applying I (6.1.3) and the second of (2.2) and noting that
${a}_{JK}^{\left(1\right)}\left({X}^{0}\right)$ can with good approximation be moved to the left of
${\partial}_{I}$ or
${\partial}_{II}$ because it varies slowly over X^{0}, (2.5) becomes

$\begin{array}{l}{\displaystyle \int {\text{d}}^{4}X{\text{d}}^{4}x}\{{\partial}_{I}^{a\stackrel{\dot{}}{b}}\left({a}_{JK}^{\left(1\right)}\left({X}^{0}\right)/{a}_{JK}\right){\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}{\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{0\stackrel{\dot{}}{b}f}+i\frac{1}{2}{q}_{p}{A}^{a\stackrel{\dot{}}{b}}\left(X\right){\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}\left({\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{0\stackrel{\dot{}}{b}f}\right)\\ +\text{\hspace{0.17em}}i\frac{1}{2}{q}_{r}\left({\partial}_{I}^{a\stackrel{\dot{}}{b}}{\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}\right){A}^{f\stackrel{\dot{}}{e}}\left(X\right){\chi}_{0\stackrel{\dot{}}{b}f}-i\frac{1}{2}{q}_{r}{\partial}_{I}^{a\stackrel{\dot{}}{b}}{\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}{A}^{f\stackrel{\dot{}}{e}}\left(X\right){\chi}_{0\stackrel{\dot{}}{b}f}\\ -\left({a}_{JK}^{\left(1\right)}\left({X}^{0}\right)/{a}_{JK}\right)\left[\left({\partial}_{I}^{a\stackrel{\dot{}}{b}}{\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}\right)\left({\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{0\stackrel{\dot{}}{b}f}\right)+\left({M}_{m}^{2}-{\Phi}_{m}\right){\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}{\psi}_{0}^{a\stackrel{\dot{}}{e}}\right]\}=0\end{array}$ (2.6)

The
${a}_{JK}^{\left(1\right)}\left({X}^{0}\right)$ factor multiplying the brackets in (2.6) can with good approximation be moved outside the integral sign because it varies slowly with X^{0}. Substituting I (6.1.3) once more into (2.6), it is seen that the surface term associated with the first term on the right of I (6.1.3) vanishes upon integration. The remaining terms in the bracket also drop out by virtue of the equation of motion I (2.3.22) for steady state mesons. The next to last terms in (2.6), the
$-i{q}_{IIr}{\partial}_{I}$ term, is also a surface term and vanishes after integration, noting that A(X) is real and satisfies periodic boundary conditions at large X^{μ} as in (2.4). Equation (2.6) can now be written as

${{S}^{\prime}}_{mAd}={{S}^{\prime}}_{mAs}$ (2.7)

${{S}^{\prime}}_{mAd}={\displaystyle \int {\text{d}}^{4}X{\text{d}}^{4}x\left[{\partial}_{I}^{a\stackrel{\dot{}}{b}}\left({a}_{JK}^{\left(1\right)}\left({X}^{0}\right)/{a}_{JK}\right){\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}{\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{0\stackrel{\dot{}}{b}f}\right]}$ (2.8)

${{S}^{\prime}}_{mAs}=i{\displaystyle \int {\text{d}}^{4}X{\text{d}}^{4}x\left\{\frac{1}{2}{q}_{p}\left[{A}^{a\stackrel{\dot{}}{b}}\left(X\right){\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}\left({\partial}_{II}^{f\stackrel{\dot{}}{e}}{\chi}_{0\stackrel{\dot{}}{b}f}\right)\right]-\frac{1}{2}{q}_{r}\left[{A}_{\stackrel{\dot{}}{f}e}\left(X\right){\chi}_{0}^{\ast b\stackrel{\dot{}}{f}}\left({\partial}_{I\stackrel{\dot{}}{a}b}{\chi}_{0}^{e\stackrel{\dot{}}{a}}\right)\right]\right\}}$ (2.9)

The last relation has been extracted from (2.5) directly.

3.4. Rudimentary Quantization and Decay Amplitude

The semi-classical treatment mentioned in the beginning of this section is analogous to the treatment of time-dependent problems in quantum mechanics. A justification is that the energies involved here are low so that typical field-theoretical effects such as vacuum polarization and self energy are small.

The following rudimentary quantization procedures are accordingly adopted. Let $|0\rangle $ and $\langle 0|$ denote vacuum states, one has conventionally,

$\langle f|i\rangle =\langle 0|i\rangle =\langle f|0\rangle =0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\langle 0|0\rangle =\langle f|f\rangle =1$ (2.10)

$|i\rangle =|V\left({\underset{\_}{K}}_{1}=0\right)\rangle ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\langle f|=\langle P\left(\underset{\_}{K}\right),{\gamma}_{T}\left({\underset{\_}{K}}_{\gamma}\right)|$ (2.11)

where V denotes a vector meson at rest which decays into a pseudoscalar meson P with momentum K and photon γ_{T} with momentum K_{g}.

Insert (2.2) into (2.8) and sandwich it between
$\langle f|$ and
$|i\rangle $ . a_{10} in χ operated on the initial state
$|i\rangle $ picks out the zeroth order initial vector meson at rest,

${\chi}_{0\stackrel{\dot{}}{b}f}|i\rangle \to {\underset{\_}{\sigma}}_{\stackrel{\dot{}}{b}f}\stackrel{^}{r}{\chi}_{10}\left(r\right)\mathrm{exp}\left(-i{E}_{10}{X}^{0}\right)$ (2.12)

where (2.1) and I (3.2.4a, 5b) have been used and χ_{10} = ψ_{10} is given in (1.6). Since (2.8) does not contain the final state photon A(X), χ* operating on the final state
$\langle f|$ becomes equivalent to the complex conjugate of the initial state (2.12). Using I (3.1.4, 10a) and (3.5.6), letting the a’s considered below (2.1) be elevated to operators according to the interpretations there and integrating over X^{0} leads to

$\langle f|{{S}^{\prime}}_{mAd}|i\rangle =-i\frac{1}{4}{E}_{10}{S}_{fi}{\displaystyle \int {\text{d}}^{3}\underset{\_}{X}{\displaystyle \int {\text{d}}^{4}x{\chi}_{10}^{2}\left(r\right)}}$ (2.13)

${S}_{fi}=\langle f|{a}_{JK}^{\ast}{a}_{JK}^{\left(1\right)}\left({X}^{0}\to \infty \right)|i\rangle $ (2.14)

S_{fi} corresponds to the conventional S-matrix element and is interpreted as the decay amplitude via the assignments of the a’s below (2.1).

Next, place (2.9) between
$\langle f|$ and
$|i\rangle $ of (2.11) and elevate a_{JK} and
${a}_{JK}^{\ast}$ to annihilation and creation operators, respectively, so that they are on the same level as a_{T}(K_{γ}) in (2.4). Here, (2.12) is applicable.
${a}_{{T}^{\prime}}^{\ast}\left({\underset{\_}{K}}_{\gamma}\right)$ in c.c. part of (2.4) picks out a final state photon of polarization T with momentum K_{γ}. The final state meson operator
${\chi}_{0}^{\ast}$ operating on
$\langle f|$ picks out a pseudoscalar meson having a momentum K with the wave function

$\langle f|{\chi}_{0\stackrel{\dot{}}{e}a}^{\ast}\to \left({\delta}_{\stackrel{\dot{}}{e}a}{\chi}_{0K}^{*}\left(\underset{\_}{x}\right)+{\underset{\_}{\sigma}}_{\stackrel{\dot{}}{e}a}{\underset{\_}{\chi}}_{0K}^{*}\left(\underset{\_}{x}\right)\right)\mathrm{exp}\left(i{E}_{0K}{X}^{0}-i\underset{\_}{K}\underset{\_}{X}\right)$ (2.15)

where K stands for K_{J} in (2.1) for J = 0. Insert (2.4, 12, 15) into (2.9), apply I (3.1.4, 10a) and (3.5.6) and integrate over X. The result reads

$\langle f|{{S}^{\prime}}_{mAs}|i\rangle =\frac{i}{\sqrt{2{E}_{r}\Omega}}{\left(2\text{\pi}\right)}^{4}\delta \left({E}_{0K}+{E}_{r}-{E}_{10}\right)\delta \left(\underset{\_}{K}+{\underset{\_}{K}}_{\gamma}\right){\displaystyle \int \text{d}{x}^{0}}{I}_{q}$ (2.16)

$\begin{array}{l}{I}_{q}={\displaystyle \int {\text{d}}^{3}\underset{\_}{x}}[\left({q}_{p}+{q}_{r}\right)\left(\underset{\_}{e}\text{\hspace{0.05em}}{\underset{\_}{\chi}}_{0K}^{*}\left(\underset{\_}{x}\right)\left(\underset{\_}{\partial}\stackrel{^}{r}{\chi}_{10}\left(r\right)\right)+\frac{1}{2}{E}_{10}\left({\underset{\_}{\chi}}_{0K}^{*}\left(\underset{\_}{x}\right)\times \underset{\_}{e}\right)\stackrel{^}{r}{\chi}_{10}\left(r\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\left({q}_{p}-{q}_{r}\right)i\frac{1}{2}{E}_{10}\underset{\_}{e}\stackrel{^}{r}{\chi}_{0K}^{*}\left(\underset{\_}{x}\right){\chi}_{10}\left(r\right)]\end{array}$ (2.17)

in which only one of T = 1 or 2 in (2.4) has been included. The
${\chi}_{0}^{\ast}$ ^{ }wave functions obey I (3.5.8) which has not been solved, as was indicated in I §3.5.3.

For the heavy mesons containing b or c quark, I Table 5.1 shows that the masses of the vector meson and the associated pseudoscalar meson are close to each other so that the latter moves slowly in a radiative decay of the former. For these mesons, the criteria (1.7, 12, 18) are satisfied and approximative solutions have been given by (1.9, 10a, 14, 16, 19, 20) where χ_{00} stands for
${\chi}_{0K}^{\ast}$ and
${\underset{\_}{\chi}}_{1}$ _{ }for
${\underset{\_}{\chi}}_{0K}^{\ast}$ in (2.17). Inserting these results together with (1.1, 6) into (2.17) and carrying out the angular integrations, terms with odd powers of
$\stackrel{^}{r}$ vanish, and

$\begin{array}{l}{I}_{q}=\left({q}_{p}+{q}_{r}\right)\frac{K}{{E}_{0K}}{e}_{Z}{I}_{dh}\\ {I}_{dh}=4\text{\pi}{\alpha}_{00}{\alpha}_{10}{\displaystyle \int \text{d}}r\text{\hspace{0.05em}}{r}^{2}\left(1+\frac{2{d}_{h}}{{E}_{0K}^{2}}\left(1-\frac{1}{3}{d}_{h}{r}^{2}\right)\right)\left(1-7{d}_{h}{r}^{2}+2{d}_{h}^{2}{r}^{4}\right)\mathrm{exp}\left(-{d}_{h}{r}^{2}\right)\end{array}$ (2.18)

Let the photon be directed along the Z axis consistent with (1.13) and put

${\underset{\_}{K}}_{\gamma}=\left(0,0,{K}_{\gamma}\right)=-\underset{\_}{K}=\left(0,0,-K\right)$ (2.19)

the iK_{γ}X term in (2.4) and the −iKX term in (2.15) cancel out in (2.9) and the last d function in (2.16) drops out. Equating (2.13) to (2.16) with (2.18, 19) according to (2.7), the decay amplitude becomes

${S}_{fi}=-\frac{4}{{E}_{10}}\frac{1}{\sqrt{2{E}_{r}\Omega}}2\text{\pi}\delta \left({E}_{0K}+{E}_{r}-{E}_{10}\right)\left({q}_{p}+{q}_{r}\right){e}_{Z}{I}_{dh}$ (2.20)

in which the integrals over the relative time x^{0} and the laboratory space X have been cancelled out and I (4.2.8) has been consulted.

3.5. Decay Rate

The decay rate is

$\Gamma \left(V\to P\gamma \right)={\displaystyle {\sum}_{\text{finalstates}}{\left|{S}_{fi}\right|}^{2}/{T}_{d}}={\displaystyle {\sum}_{T}{\displaystyle {\sum}_{K}{\left|{S}_{fi}\right|}^{2}/{T}_{d}}}$ (2.21)

where T_{d} is a long time during which decay takes place and Σ_{T} = 2 for the both photon polarizations mentioned below (2.4). Further,

${{\displaystyle \sum}}_{K}^{}=\frac{\Omega}{{\left(2\text{\pi}\right)}^{3}}{\displaystyle \int {\text{d}}^{3}\underset{\_}{K}}\to \frac{\Omega}{{\left(2\text{\pi}\right)}^{3}}4\text{\pi}{\displaystyle \int \text{d}}K{K}^{2}$ (2.22)

With (2.19) and E_{r} = K_{γ}, one finds

$\delta \left({E}_{0K}+{E}_{r}-{E}_{10}\right)=\delta \left(\sqrt{{E}_{00}^{2}+{K}^{2}}+K-{E}_{10}\right)=\delta \left(K-{K}_{0}\right)\left(1-{K}_{0}/{E}_{10}\right)$ (2.23a)

${K}_{0}=\left({E}_{10}^{2}-{E}_{00}^{2}\right)/2{E}_{10}$ (2.23b)

This δ function in (2.20) will be squared in (2.21). Using (2.23a), this square becomes

${\left[\delta \left({E}_{0K}+{E}_{r}-{E}_{10}\right)\right]}^{2}=\delta \left(K-{K}_{0}\right){\left(1-\frac{{K}_{0}}{{E}_{10}}\right)}^{2}\frac{{T}_{d}}{\text{2\pi}}$ (2.24)

Combining (2.20-22, 24) leads to the decay rate

$\Gamma \left(V\to P\gamma \right)=\frac{16}{\text{\pi}}\frac{{K}_{0}^{3}}{{E}_{10}^{2}\left({E}_{00}^{2}+{K}_{0}^{2}\right)}{\left(1-\frac{{K}_{0}}{{E}_{10}}\right)}^{2}{\left({q}_{p}+{q}_{r}\right)}^{2}{e}_{Z}^{2}{I}_{dh}^{2}$ (2.25)

Here, e_{Z} is not fixed here but is a parameter < 1. This is as far as the present semi-classical treatment can carry. Perhaps a quantized version of the present treatment, if possible to be devised, can fix this e_{Z}.

e_{Z} is the third, longitudinal, component of the unit verctor e in the photon field (2.4) and balances off the momentum K_{0} of the pseudoscalar meson. Its absence in the earlier treatment [1] and I led to the assumption that higher order effects were responsible for such radiative decays. This is now seen not to the case. e_{Z}, being separate from the photon associated with the transverse components e_{X}, e_{Y} in (2.4) is responsible for the decays.

The first order (2.9) can be complemented to include second order terms containing the square of the photon field A(X). The coefficient for this square representing the mass of the photon contains only odd powered $\stackrel{^}{r}$ terms and vanishes after integration over the relative space x. This verifies that the photon remains massless.

Finally, it is remarked that (2.25) remains unchanged if the employed equation for χ (2.2) were replaced by an equivalent equation for ψ analogous to I (6.1.10b). This is due to that I (6.1.10a, 10b) originate from I (2.2.4a, 4b) which are invariant under the interchanges χ « ψ and I « II. The last one leads to a sign change of the relative coordinate x in I (3.1.3a) but does not affect the decay amplitude (2.20).

3.6. Comparison with Data

Table 1 gives the V → Pγ decay rate for some heavy mesons most of them listed in I Table 5.1.

Table 1. Decay rate Γ(V → Pγ) for some heavy mesons. The meson momentum K_{0} is given by (2.23b). The expansion parameter ε_{0} is defined in (1.7). q_{p} and q_{r} are the quark charges. The predicted decay rate Γ(V → Pγ) is given by (2.25) which holds only for nonrelativistic and heavy pseudoscalar mesons satisfying
${e}_{0}\ll \text{1}$ of (1.7) and (1.12, 18). Only the mesons with b or c quark qualify. Equation (2.25) does not apply to K*. This is indicated by the parentheses around the so-obtained decay rates. e_{Z} is the longitudinal polarization vector in the unit vector e in (2.4) and should be < 1 but is not fixed in the present nonquantized treatment. The underlined entries are the only ones available for direct comparison with predictions.

For mesons with b quark, there is no data available. The K* decay rates are put inside parentheses to indicate that Γof (2.25) is not applicable. This is due to thate_{0} is too large and does not satisfy the criteria (1.7) and the kaons are too light so that (1.12, 18) are violated. Further, the effect of the singularity associated with lighter mesons mentioned above §5.7.2 has not been investigated and may impact upon the validity of the approximations leading to (2.25). Still the predictions of (2.25) for kaons are included to indicate the trend that, as the mesons get lighter and move faster; the decay rates increase. Here, they are 10 - 15 times too high.

For mesons containing c quark, J/ψ and D*^{±} rates from (2.25) can be compared to data. The ratio between these two measured rates is 1.58/1.33 = 1.19 and can be 0.73 up to 1.95 within error limits. The predicted value 1.02/0.69 = 1.48 lies well within these limits. This agreement can hardly be any coincidence inasmuch as the ratio between their (q_{p} + q_{r})^{2} values is 16.

Putting e_{Z} = 1, the measured rates are greater than the predicted ones by a factor of 1.5 - 1.9. Since e_{Z} < 1 in the present semi-classical treatment, this factor is actually still greater. Nevertheless, the predicted values and data are of the same magnitude and may be regarded to be in basic agreement with each other. Consider that the mass ratios of K*, D* and B* are 0.44:1:2.86 while the ratios of predicted decay rates are 632:1:4.3 × 10^{−5}, the above discrepancy of a factor of 1.5 - 1.9 or greater is negligible. The approximations introduced in (1.19, 20) can however only lead to a correction of ±10% for D* in (2.25).

These results indicate that the harmonic oscillator type of meson wave functions (1.1, 6) is useful not only for meson spectra but also for decay problems. They provide further support to the scalar strong interaction hadron theory. No other first principles’ theory can make a prediction of this kind.

If the present semi-classical treatment can be quantized and e_{Z} can be fixed, a wholly new problem, a more precise prediction may be found.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] |
Hoh, F.C. (1999) International Journal of Theoretical Physics, 38, 2647-2664. https://doi.org/10.1023/A:1026617512188 |

[2] |
Hoh, F.C. (2011) Scalar Strong Interaction Hadron Theory. Nova Science Publishers, Hauppauge, 2017, Revised Second Edition, Denoted by I. https://www.novapublishers.com/catalog/product_info.php?products_id=27069 |

[3] |
Hoh, F.C. (1993) International Journal of Theoretical Physics, 32, 1111-1133. https://doi.org/10.1007/BF00671793 |

[4] | Patrignani, C., et al. (Particle Data Group) (2016) Chinese Physics C, 40, Article ID: 100001. |

[5] |
Hoh, F.C. (2017) Journal of Modern Physics, 8, 1127-1133. https://doi.org/10.4236/jmp.2017.88073 |

[6] |
Lichtenberg, D.B. (1987) International Journal of Theoretical Physics, A2, 1669. https://doi.org/10.1142/S0217751X87000879 |

[7] | Ludwig, W. and Falter, C. (1988) Symmetry in Physics. Springer, Berlin. |

[8] |
Hoh, F.C. (1996) Journal of Physics G: Nuclear and Particle Physics, G22, 85. https://doi.org/10.1088/0954-3899/22/1/007 |

[9] |
Eichten, E., et al. (2008) Reviews of Modern Physics, 80, 1161. https://doi.org/10.1103/RevModPhys.80.1161 |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.