Anomalous radiative transitions

Anomalous transitions involving photons derived by many-body interaction of the form, $\partial_{\mu} G^{\mu}$, in the standard model are studied. This does not affect the equation of motion in the bulk, but makes wave functions modified, and causes the unusual transition characterized by the time-independent probability. In the transition probability at a time-interval $T$ expressed generally in the form $P=T \Gamma_0 +P^{(d)}$, now with $ P^{(d)} \neq 0 $. The diffractive term $P^{(d)}$ has the origin in the overlap of waves of the initial and final states, and reveals the characteristics of waves. In particular, the processes of the neutrino-photon interaction ordinarily forbidden by Landau-Yang's theorem ($\Gamma_0=0$) manifests itself through the boundary interaction. The new term leads to physical processes over a wide energy range to have finite probabilities. New methods of detecting neutrinos using laser are proposed that are based on this difractive term, which enhance the detectability of neutrinos by many orders of magnitude.


Matter wave and S[T]
In modern science and technology, quantum mechanics plays fundamental roles. Despite the fact that stationary phenomena and method have been well developed, those of non-stationary phenomena have not. In the former, the de Broglie wave length p determines a typical length and is of a smicroscopic size, and scatterings or reactions in macroscopic scale are considered independent, and successive reactions have been treated under the independent scattering hypothesis. The probability of the event that they occur are computed by the incoherent sum of each value. In the latter, time and space variables vary simultaneously, and a new scale, which can be much larger than the de Broglie wave length, emerges. They appear in overlapping regions of the initial and final waves, and show unique properties of intriguing quantum mechanical waves.
A transition rate computed with a method for stationary waves with initial and final states defined at the infinite time-interval T = ∞ is independent of the details of the wave functions. They hold characteristics of particles and preserve symmetry of the system. Transitions occurring at a finite T , however, reveal characteristics of waves, the dependence on the boundary conditions [1,2] 1 , and the probability, where P (d) is the diffractive term which has often escaped attention by researchers. The rate Γ 0 is computed with Fermi's golden rule [6,7,8,9,10], and preserves the internal and space-time symmetries, including the kineticenergy conservation. Γ 0 holds the characteristic properties of particles and the hypothesis of independent scatterings is valid. For a particle of small mass, m s , Γ 0 (p i , m s ) behaves because the characteristic length, the de Broglie wave length, is determined with p i . The region where P (d) is ignorable is called the particle-zone.
Overlapping waves in the initial and final states have the finite interaction energy, and reveal unique properties of waves [1,2]. Because the interaction energy is part of the total energy, sharing with the kinetic energy, the conservation law of the kinetic-energy is violated. Consequently, the state becomes non-uniform in time and the transition probability has a new component P (d) , showing characteristics of waves. The term P (d) T was shown to behave with a new scale of length, ( msc ) · ( E i msc 2 ). Accordingly the correction is proportional to the ratio of two small quantities and does not follow Eq. (2). P (d) reflects the non-stationary waves and is not computed with the stationary waves. In the region where P (d) is important, the hypothesis of independent scattering is invalid, and interference unique to waves manifests. This region is called the wave-zone, and extends to a large area for light particles. P (d) has been ignored, but gives important contributions to the probability in various processes. Especially P (d) is inevitable for the process of Γ 0 = 0 and P (d) = 0, which often appears. Furthermore, P (d) can be enhanced drastically, if the overlap of the waves is constructive in wide area. This happens for small m s , large E i , even for large T , and reveals macroscopic quantum phenomena. Processes of large P (d) involving photon and neutrino in the standard model are studied in the present paper.
An example of showing Γ 0 = 0, P (d) = 0 is a system of fields described by a free part L 0 and interaction part L int of total derivative, where G is a polynomial of fields φ l (x). φ l (x) follows free equation, L int decouples from the equation and does not modify the equation of motion in classical and quantum mechanics. Nevertheless, a wave function |Ψ(t) follows a Schrödinger equation in the interaction picture, where the free part, H 0 , and the interaction part, H int , are derived from the previous Lagrangian, and G int stands for G of the interaction picture. A solution at t, is expressed with G(t), and the state at t > 0 is modified by the interaction. The initial state |Ψ(0) int prepared at t = 0 is transformed to the other state of t-independent weight. Hence, Eq. (7), is like stationary, and Γ 0 = 0 , and P (d) = 0. Physical observables are expressed by the probability of the events, which are specified by the initial and final states. For those at finite T , normal S-matrix, S[∞], which satisfies the boundary condition at T = ∞ instead of those at finite T , is useless. S[T ] that satisfies the boundary condition at T [1,2] is necessary and was constructed. S[T ] is applied to the system described by Eq. (4).
hence a matrix element of S[T ] between two eigenstates of H 0 , |α and |β of eigenvalues E α and E β , is decomposed into two components where β|S (n) |α and β|S (d) |α get contributions from cases E β = E α and E β = E α , and give Γ 0 T and P (d) respectively. The deviation of the kineticenergies, E β − E α , in the latter is due to the interaction energy of the overlapping waves, which depends on the coordinate system. Therefore, it is understood that H int is not Lorentz invariant. Thus the kinetic-energy nonconserving term, which was mentioned by Pierls and Landau [4] so as to give negligibly small correction, yields P (d) [5] 2 . Because H 0 is an generator of the Poincare group, Eq. (8) shows that S (d) [T] and P (d) violate the Poincare invariance. In the system described by Eq. (4), the first term disappears but the second term does not, Γ 0 = 0, and P (d) = 0.
S[T ] is expressed with the boundary conditions for the scalar field φ(x) [12,13], where φ in (x) and φ out (x) satisfy the free wave equation, and φ f , φ f in , and φ f out are the expansion coefficient of φ(x), φ in (x), and φ out (x), with the normalized wave functions f (x) of the form The function f (x) indicates the wave function that the out-going wave interacts in a successive reaction of the process. The out-going photon studied in the following section interacts with atom or nucleus and their wave functions are used for f (x  [12,13,8,9,10] can be replaced with plane waves for a practical computation for S[∞] [18,19,20,21,22,23,24,25], but can not be done so for P (d) [14,15,16]. P (d) is derived from S (d) [T ; f ], and depends on f (x).
Photon is massless in vacuum and has a small effective mass determined by the plasma frequency in matter, and neutrino is nearly massless. Thus they have the large wave-zone of revealing wave phenomena caused by P (d) . These small masses make P (d) appear in a macroscopic scale and significantly affect physical reactions. In this small (or zero) mass region the effects of diffractive term P (d) are pronounced. This is our interest in the present paper. The produced photon interacts with matter with the electromagnetic interaction, which leads to macroscopic observables. The term P (d) of the processes of Γ 0 = 0 such as 2γ decays of 1 + meson and γ and ν reactions are shown to be relevant to many physical processes, including possible experimental observation of relic neutrino. The enhancement of the probability for light particles with intense photons based on the normal component Γ 0 was proposed in Refs. [27,28], and the collective interaction between electrons and neutrino derived from the normal component Γ 0 was considered in Ref. [29]. Our theory is based on the probability P (d) , hence differs from the previous ones in many respects.
This paper is organized in the following manner. In Section 2, the couplings of two photon with 1 + state through triangle diagram is obtained. In section 3 and 4, positronium and heavy quarkonium are studied and their P (d) are computed. Based on these studies we go on to investigate the interaction of photons and neutrinos. In section 5, neutrino-photon interaction of the order αG F , and various implications to high energy neutrino phenomena are presented. In Section 6 we explore the implication of the photon-neutrino coupling on experimental settings. Summary is given in section 7.

Coupling of 1 + meson with two photons
The coupling of γγ with axial vector states, 1 + meson composed of e − e + , QQ, and νν are studied using an effective Lagrangian expressed by local fields. From symmetry considerations, an effective interaction of the 1 + state φ µ 1 with two photons has the form, where F αβ is the electromagnetic field, and the coupling strength g is computed later. In a transition of plane waves in infinite-time interval, the spacetime boundary is at the infinity, and the transition amplitude is computed with the plane waves, in the form, and vanishes, where p i and p f are four-dimensional momenta of the initial and final states andM is the invariant amplitude. This shows that the amplitude proportional to δ 4 (p i − p f ) and the transition rate Γ 0 vanish. The rate of 1 + → γγ decay vanishes in general systems, because the state of two photons of momenta ( p, − p) does not couple with a massive 1 + particle. Hence which is known as Landau-Yang's theorem [30,31]. The term S int is written as a surface term in four-dimensional space-time, which is determined by the wave functions of the initial and final states. Accordingly, the transition amplitude derived from this surface action is not proportional to T , but has a weaker T dependence. Thus P (d) comes from the surface term, and does not have the delta function of kinetic-energy conservation. Kinetic energy of the final states deviates from that of the initial state due to the finite interaction energy between them. The deviation becomes larger and P (d) is expected to increase with larger overlap. We find P (d) in the following.

Triangle diagram
The interaction of the form Eq. (13) is generated by one loop effect in the standard model. The scalar and axial vector current in QED, we have where A µ (x) is the photon field and l(x) is the electron field, coupled through two photons in the bulk through the triangle diagram Fig. 1. The matrix elements are where ǫ µ (k) is the polarization vector for the photon. The triangle diagram for the axial vector current Eq. (21) has been studied in connection with axial anomaly and π 0 → γγ decay [32,33,34,35,36] and now is applied to P (d) for two photon transitions of the axial vector meson and neutrino. The triangle diagram Fig. 1 shows that the interaction occurs localy in space and time, but the transition amplitude is the integral over the coordinates and receives the large diffractive contribution if the neutrino and photon are spatially spread waves. In Fig. 1, the in-coming and out-going waves are expressed by lines, but they are in fact the spread waves, which is obvious in the figures in Ref. [37] and in Fig. 2. Γ 5,α is expressed also with f 1 in the form The coefficients f 0 and f 1 are given by the integral over the Feynman parameters, γ (detected) ν γ ν Figure 2: Diagram of the neutrino photon scattering, ν + γ → ν + γ of spatially spread waves. where f 1 in various kinematical regions is

Positronium
The bound states of a positronium with the orbital angular momentum L = 1, S = 1 have total angular momentum J = 2, 1, 0. These states at rest of P = 0 are expressed with non-relativistic wave functions and creation operators of l + and l − of momentum p, spin ±1/2, and P-wave wave functions p i F (p) as are given in Appendix A, where Others are defined in the same manner. These bound states couple with the lepton pairs with effective local interactions which are expressed as where the coupling strengths are computed from where the ǫ µ , ǫ µν are polarization vector or tensor for the massive vector and tensor mesons. We have The decay rates for 0 + and 2 + were studied in [38,39], so that here we concentrate on 1 + and 0 + as a reference for 1 + . Fields φ 0 (x), φ 1 (x), and φ 2 (x) couple with two photons through the triangle diagrams. Their interactions with two photons are summarized in the following effective Lagrangian,

Axial vector positronium
Here we study the two photon decay of axial vector positronium, which is governed by the second term of the right-hand side of Eq. (34). The matrix element of the axial current between the vacuum and two photon state was computed by Refs. [32,33,34,35,36]. Because Γ 0 of two photon decay of axial vector meson vanishes due to Landau-Yang's theorem, but P (d) does not, we give the detailed derivation of P (d) .
From the effective interaction, Eq. (34), the probability amplitude of the event that one of the photons of k γ from the decay of c µ of p c is detected at The initial state is normalized, and the coupling in Eq. (33) has (2π) 3 2 for the initial state, and The state k γ , X γ , T γ | is normalized and, Integration over k 2 is made prior to the integration over x, in order for M to satisfy the boundary condition of S[T ], and we have, where The wavepacket expands in the transverse direction, and σ T in large x 0 − T γ is given by We later use The stationary phase for large x 0 − T γ exists in the time-like region λ ≥ 0 [14], and the function and its derivative is proportional to θ(λ). Thus the integration in Eq. (35) is made over the region λ ≥ 0, which has the boundary, λ = 0. Consequently the transition amplitude does not vanish if the integrand at the boundary is finite. It is shown that this is the case in fact. σ γ is the size of the nucleus or atomic wave function that the photon interacts with and is estimated later. For the sake of simplicity, we use the Gaussian form for the main part in this paper.
Substituting Eq. (21), we have the amplitude which depends on the momenta and coordinates (T γ , X γ ) of the final state and T m . Although M is written as the integral over the surface, λ = 0, the expression Eq. (44) is useful and applied for computing the probability per particle where V is a normalization volume for the initial state, and the momentum of the non-observed final state is integrated over the whole positive energy region and the position of the observed particle are integrated inside the detector.
Following the method of the previous works [1,2] and Appendix B, we write the probability with a correlation function. In the integral With variables the integral is written as where λ + = x 2 + and x + is integrated over the region λ + ≥ 0, and the integral is computed with the value at the boundary, λ + = 0. δx are integrated over the whole region, and the second term in the second line vanishes.
Using the formulas [1,2], and the integrals given in Appendix E, we compute the probability. It is worthwhile to note that the right-cone singularity exists in the kinematical region θ(phase-space) and the probability becomes finite in this region [1,2]. The natural unit, c = = 1, is taken in the majority of places, but c and are written explicitly when it is necessary, and MKSA unit is taken in later parts. After tedious calculations, we have the probability where I i is given in Appendix E, and and Integrating over the gamma's coordinate X γ , we obtain the total volume, which is canceled by the factor V −1 from the normalization of the initial state. The total probability, for the high-energy gamma rays, is then expressed as where L = cT is the length of the decay region. The kinematical region of the final states is expressed by the step function [2], which is different from the on-mass shell condition. Γ 0 vanishes and P is composed of log T and constant terms. The constant is inversely proportional to m 2 γ , and becomes large for the small m γ .

Transition probability
A photon γ is massless in vacuum and has an effective mass in matter, which is given by plasma frequency in matter [40] as where n e is the density of electron, and for the value, Substituting we have, In the air, n e = 3 × 10 25 /m 3 , the mass agrees with which is comparable to the neutrino mass. At a macroscopic T , ωT is much larger than 1.
We have the probability at the system of p c = (E c , 0, 0, p c ), and the total probability

Heavy quarkonium
The decay of axial vector mesons composed of heavy quarks exhibits the same phenomena. The heavy quark mesons composed of charm and bottom quarks are observed, and show rich decay properties in the two photon decay, radiative transitions, and two gluon decays. Because quarks interact via electromagnetic and strong interactions, non-perturbative effects are not negligible, but the symmetry consideration is valid. Moreover, the nonrelativistic representations are good for these bound states, because quark masses are much greater than the confinement scale. Furthermore they have small spatial sizes. Accordingly we represent them with local fields and find their interactions using the coupling strengths of Eq. (33) and the values of the triangle diagrams, Eqs. (20), and (21).

qq → γ + γ
Up and down quarks have charge 2e/3 or e/3, and color triplet. Hence the probabilities of two photon decays are obtained by the expression Eq. (67) with charges of quarks 2e/3, e/3, and color factor. It is highly desired to obtain the experimental value for 1 + .

qq → gluon + gluon
A meson composed of heavy quarks decays to light hadrons through gluons. Color singlet two gluon states are equivalent to two photon states. Accordingly the two gluon decays are calculated in the equivalent way to that of two photon decays as far as the perturbative calculations are concerned. The transition rates for 0 + and 2 + may be calculated in this manner. The total rates for L = 0 charmonium, J/Ψ and η c agree with the values obtained by the perturbative calculations. J/Ψ is C = −1 and decays to three gluons and the latter is C = +1 and decays to two gluons. The former rate is of α 3 s and the latter rate is of α 2 s , where α s is the coupling strength of gluon. Their widths are Γ = 93 keV or Γ = 26.7 MeV, and are consistent with small coupling strength α s ≈ 0.2. Now L = 1 states have C = +1. Hence a meson of J = 2 and that of J = 0 decay to two gluons, whereas the decay rate of a J = 1 meson vanishes by Landau-Yang's theorem.
A cc meson of the quantum number 1 + is slightly different from that of positronium, because a gluon hadronizes by a non-perturbative effect, which is peculiar to the gluon. The hadronization length is not rigorously known, but it would be reasonable to assume that the length is on the order of the size of pion. The time interval T is then a microscopic value. P (d) for this time T is estimated in the following.
A gluon also hadronizes in the interval of the lightest hadron size, The gluon plasma frequency is estimated from quark density, and Then we have At small ωT ,g(ωT ) varies with T as is shown in Fig. 2 of Ref. [2] and Higher order corrections also modify the rates for qq. The values to light hadrons are estimated by [41,42,43,44,45]. The total values for light hadrons are expressed with a singlet and octet components H 1 and H 8 as Using values for χ 0 and χ 2 from [49] Γ(χ 0 → light hadrons) = 1.8 MeV, we have the rate for χ 1 from H 8 , The experimental value for χ 1 is The large discrepancy among experiments for Γ χ 1 →light hadrons may suggest that the value depends on the experimental situation, which may be a feature of P (d) . We have

E1 transition:
which is expressed by the effective interaction in the local limit. The action Eq. (84) does not take the form of the total derivative, but is written as The second term shows the interaction of the local electromagnetic coupling of the current j µ = ∂ ν O µν and the first term shows the surface term. This leads to the constant probability at a finite T , P (d) . The induced P (d) for muon decay was computed in [1], and large probability from P (d) compared with T Γ 0 was found in the region cT = 1 m. The radiative transitions of heavy quarkonium are deeply connected with other radiative transitions and the detailed analysis will be presented elsewhere.
Spin 0 and 2 mesons , φ 0 and φ 2 , show the same E 1 transitions and photons show the same behavior from P (d) [46,47,48,49]. A pair of photons of the continuous energy spectrum are produced in the wave zone, and are correlated. On the other hand, a pair of photons of the discrete energy spectrum are produced in the particle zones and are not correlated. Thus the photons in the continuous spectrum are different from the simple background, and it would be possible to confirm the correlation by measuring the time coincident of the two photons.

Neutrino-photon interaction
The neutrino photon interaction of the strength αG F induced from higher order effects vanishes due to the Landau-Yang-Gell-Mann theorem. But P (d) is free from the theorem and gives observable effects. Moreover, although the strength seems much weaker than the normal weak process, that is enhanced drastically if the photon's effective mass is extremely small. Because P (d) does not preserve Lorentz invariance, a careful treatment is required. The triangle diagram of Fig. 1 is expressed in terms of the action , where the axial vector meson in Eq. (34) was replaced by the neutrino current. The mass of the neutrino is extremely small, and was neglected. The S νγ leads to the neutrino gamma reactions with Γ 0 = 0, P (d) = 0 on the order αG F .

ν + γ → ν + γ
The rates Γ 0 of the events vanishes on the order αG F [50] due to Landau-Yang's theorem. The higher order effects were also shown to be extremely small [34] and these processes have been ignored. The theorem is derived from the rigorous conservation law of the kinetic-energy and angular momentum. However these do not hold in P (d) , due to the interaction energy caused by the overlap of the initial and final wave functions. Consequently, neither P (d) nor the transition probability vanish. These processes are reconsidered with P (d) .
From Eq. (86), we have the probability amplitude of the event in which one of the photon of k γ interacts with another object or is detected at X γ as The amplitude is expressed in the same manner as Eq. (44), where the corrections due to higher order in δ k are ignored in the following calculations as in Section (3.1). The probability averaged over the initial spin per unit of particles of initial state, is where V is the normalization volume of the initial state, and where I 0 and I i are given in Appendix E. The probability P is not Lorentz invariant and the values in the CM frame of the initial neutrino and photon and those of the general frame do not agree generally.

Center of mass frame
The phase space integral over the momentum in the CM frame Thus we have the probability ; high-energy, (95) ; low-energy, whereω in the log term is the average ω, and ǫ is the deviation of the index of refraction from unity given in Appendix C. The log term was ignored in the right-hand side of the low energy region. We found that in the majority of region, 1 ωE √ σγ becomes much larger than √ σ γ log(ωT ). Here we assume that the medium is not ionized.

Moving frame
At the frame p ν 1 = (0, 0, p ν 1 ), p γ 1 = (0, 0, −p γ 1 ), p ν 1 > p γ 1 the probability in high energy region is In the low energy region, the second term of P is given in the form, and the first term proportional to p 6 ν 1 was ignored. Numerical constants C 0 and C ′ 0 are proportional to ( G F √ 2 2α π ) 2 . The photon effective mass at high energy region, m γ , and the deviation of the refraction constant from unity at low energy region, ǫ, are extremely small in the dilute gas, and P (d) becomes large in these situations.

Neutrino interaction with uniform magnetic field
ν + B → ν + γ The action Eq. (86) leads to the coherent interactions of neutrinos with the macroscopic electric or magnetic fields. These fields are expressed in MKSA unit. Accordingly we express the Lagrangian in MKSA unit, which is summarized inAppendix D, and compute the probabilities. The magnetic field B in the z-direction is expressed by the field strength and we have the action where g B is given in Appendix D. Because S νγ (B) is reduced to the surface term, the rate vanishes, Γ 0 = 0, but P (d) = 0. Furthermore, S νγ (B) is not Lorentz invariant, and P (d) for ν i → ν j + γ becomes proportional to m 2 ν of much larger magnitude than the naive expectation.
The amplitude is We have the probability from Eq. (46), The convergence condition on the light-cone singularity is satisfied in the kinematical region, θ(phase-space). Thus the momentum satisfies Solving k γ , we have the condition for fraction x = kγ where α ± = O(1). The process ν i → ν j + γ occurs with the probability Eq.
which is satisfied in dilute gas. If the inequality Eq. (107) is not satisfied, this probability vanishes. The probablity P reflects the large overlap of initial and final states and is not Lorentz invariant. Consequently, although the integration region in Eq. (103) is narrow in phase space determined by θ(phase-space), which is proportional to the mass-squared difference, δω 2 12 , the integrand is as large as p 3 ν . P becomes much larger than the value obtained from Fermi's golden rule.

High energy neutrino
At high energy, Eq. (A.49) are substituted. For the case that p ν 1 is nonparallel to B, we have where ζ is the angle between p ν 1 and B, The integral I 0 is almost independent of k γ . Ignoring the dependence, we integrate the photon's momentum, and have The probability P of Eq. (110) is proportional to (ω 0 ν ) 2 k 3 ν 1 , which is very different from the rate of the normal neutrino radiative decay, Γ = G 2 F m 5 ν (m ν 1 /E ν 1 )× (numerial factor), especially for high energy neutrino. Moreover, I 0 = 1 m 2 γ , can be extremely large in dilute gas, thus P is enormously enhanced.
If the momentum of initial neutrino is parallel to the magnetic field, ζ = 0, we have P in Eq. (111) is proportional to (ω 0 ν ) 4 k ν 1 , and is negligibly small compared to that of Eq. (110).

Low energy neutrino
At low energy, Eq. (A.52) are substituted. I 0 is inversely proportional to k 2 γ and we have dP dx and Using α ± , C(m) and D(m) are computed easily. 25

Neutrino interaction with uniform electric field
For a uniform electric field in the z-direction, we have ν i → ν j + γ in the electric field is almost the same as that in the magnetic field.

High energy neutrino
We have the probability in the high energy region,

Low energy neutrino
The probability in the low energy region is

Neutrino interaction with nucleus electric field
In space-time near a nucleus, there is the Coulombic electric field E nucl due to the nucleus, and one of F µν in the action is replaced with E nucl . The rate estimated in Ref. [34] was much smaller by a factor 10 −4 or more than the value of the normal process due to charged current interaction. Here we estimate P (d) for the same process. The action becomes which causes the unusual radiative interaction of neutrino in matter. The probability P (d) of the high energy neutrino, where 2k 1 · k 2 ≫ m 2 e , and we substitute the value f 1 = 1 8m 2 l . Since F µν nucl due to a bound nucleus is shortrange, the probability is not enhanced.

Neutrino interaction with laser wave
In the scattering of neutrino with a classical electromagnetic wave due to laser, one of F µν in the action is replaced with the electromagnetic field E laser of laser of the form The probability P (d) of this process is computed with the action where 2k 1 · k 2 ≈ 0, and we substitute the value f 1 = 1 8m 2 l . The amplitude and probability are almost equivalent to those of the uniform electric field.

Implications to neutrino reactions in matter and fields
An initial neutrino is transformed to another neutrino and a photon following the probability P (d) . The photon in the final state interacts with a microscopic object in matter with the electromagnetic interaction, and loses the energy. Thus the size σ γ in P (d) is determined by its wave function, and the probability P (d) × σ γA , where σ γA is the cross section of the photon and is much larger than that of weak reactions, determines the effective cross section of the whole process. Hence the effective cross section can be as large as that of the normal weak process caused by the charged current interaction.

Effective cross section
The probability of the event that the photon reacts on another object is expressed by P (d) in Eqs. (96), (97), and (98), and that of the final photon.
If a system initially has photons of density n γ (E γ ), the number of photons are multiplied, and the probability of the event that the initial neutrino is transformed is given by P (d) × n γ . In the system of electric or magnetic fields, the initial neutrino is transformed to the final neutrino and photon. Hence P (d) ×n γ for the former case, and P (d) for the latter case are important parameters to be compared with the experiments. The effective cross section, for the process where the photon in the final state interacts with atoms A of the cross section σ γA , is for the former case, and σ (d) for the latter cases. The cross sections Eqs. (122) and (123) are compared with that of the charged current weak process Since σ γA is much larger than σ weak νA , by a factor 10 14 or more, σ Accordingly, 10 −14 or 10 −15 is the critical value for the photon neutrino process to be relevant and important. If the value is larger, then the reaction that is dectated to vanish due to Landau-Yang's theorem manifests with a sizable probability.

6.2
The probability P (d) is of the order of αG F and almost independent of time. The probability in this order has been considered vanishing, and this process has not been studied. If the magnitude is sizable, these neutrino processes should be included in astronomy and others. The process ν +ν → γ + γ is almost equivalent to ν + γ → ν + γ, and we do not study in this paper.
A system of high temperature has many photons, and a neutrino makes a transition through its collision with the photons. The probability is determined by the product between the number of photons n γ and each probability In a thermal equilibrium of higher temperature, the density is about and we have the product for a head-on collision

The sun
In the core of the sun, R = 10 9 meters, (127) kT ≈ 2 keV, and the solar neutrino has the energy around 1 − 10 MeV. The photon's energy distribution is given by the Planck distribution, and the mean free path for the head-on collision is where σ γA = 10 −24 cm 2 , m γ = 1 eV and n A = 10 29 /cm 3 are used. The value is much longer than the sun's radius. For the neutrino of higher energy, we use Eq. (97), and have p 0 ν = 10 GeV, thus the length exceeds the sun's radius for p ν > 25 GeV. The high energy neutrino does not escape from the core if the energy is higher than around 25 GeV.
The mean free path becomes smaller in the lower matter density region. Thus in the region of small photon's effective mass, the neutrino does not escape but loses its majority of energy. This is totally different from the standard behavior of the supernova neutrino.

Neutron star
If magnetic field is as high as 10 9 [Tesla], then the probability becomes large. The energy of the neutrino is transfered to the photon's energy.

Low energy reaction
In low energy region, P (d) behaves as Eq.(96) and The effective transition probability P (d) n γ and the cross section depend on the photon density.

ν + (E, B) → ν ′ + γ
The radiative transition of one neutrino to another lighter neutrino and photon in the electromagnetic field occurs with the probability P (d) . Because P (d) is not proportional to T but almost constant, the number of parents decreases fast at small T , and remains the same afterward without decreasing. Now the photon in the final state reacts with matter with sizable magnitude and the probability of whole process is expressed with the effective cross section.

High energy neutrinos
The transition probability of the high energy neutrino in the magnetic field B [Tesla] and electric field E [V/meter] are For the parameters, m γ c 2 = 10 −9 eV, √ σ γ = 10 −13 meters, they become in the magnetic field

Low energy neutrinos
The transition probability of the low energy neutrino are They become in the typical situations ǫ = 10 −20 , √ σ γ = 10 −13 m, p ν = 1 eV, E 0 = 10 3 GV/m, (140) The probabilities Eqs. (136), (137), and (140) are caused by the overlap of waves on the parent and daughters. The phases of waves become cancelled at small m γ or ǫ, and more waves are added constructively then. Thus the effect become larger as they become smaller. The probability shows this and is inversely proportional to m 2 γ or ǫ. They become extremely large as m γ → 0 or ǫ → 0.

Table of processes
The P (d) could be tested in various neutrino processes of wide energy regions. The diffractive term P (d) becomes substantial in magnitude at high energy or fields. So this process may be relevant to the neutrino of high energy or at high fields, and may give new insights into or measurability of the following processes.
1. The accelerator neutrino has high energy and total intensity of the order 10 20 . At B = 10 Tesla and P ν = 1(10) GeV, we have P ≈ 10 −18 (10 −15 ), and using the laser of E 0 = 10 3 GV/m at P ν = 1(10) GeV, we have P ≈ 10 −16 (10 −13 ). The neutrino from the accelerator may be probed by detecting the final gamma rays. For example, at the beam dump of LHC we may set up laser to detect the neutrino interaction.
2. For the reactor neutrino, of the flux ≈ 10 20 /s per reactor, a detector of high magnetic field of the order 10 Tesla, or of high intensity laser, may be able to detect the neutrino.
3. Direct observations of solar neutrinos, which have the energy in 0.5 MeV to 10 MeV and flux around 10 15 /m 2 sec, using the detector with strong magnetic field similar to that for axion search [51] would be possible. 4. In Supernova or neutron star, the neutrino photon reaction would give a new important process, because the final photon interacts with matter strongly. As to the detectability of neutrinos from the sun, 10 GeV is the threshold from the sun, while they from SN interact with 10 MeV photons. 5. We find that P (d) becomes maximum at ζ = π 2 from Eq.(118), and may apply this effect to enhancing the neutrino flux. The term also has the momentum dependence, which may be exploites for "optical" effect of neutrino through the photon interaction.
6. The photon neutrino reaction may be useful for the relic neutrino detection. The reaction rate may be enhanced with such method as the neutrinos mirrors that collect them. [52] We may take advantage of the above effect.
7. The probability becomes huge at extremely high energy. So, this process may be relevant to ultra-high energy neutrino process. [53,54,55]), 8. Neutrino may interact with electromagnetic fields in Cumulonimbus cloud.
8-1. Lightening has total energy of the order 900 MJ ≈ 10 9 CV and current 10 6 A in a short period. cB at a radius r = 1 cm is 1.5 × 10 9 N/C. Assuming E = cB = 1.5 × 10 9 N/C, and m γ c 2 = 10 −11 eV, P ν = 10 MeV we have P E+B ≈ 10 −15 . The neutrino inevitably loses its energy, and the photons of the continuous spectrum are emitted. This may be related with the upper-atmospheric lightening [56] .
8-2. The gamma rays observed in Cumulonimbus cloud, [57], may be connected with the diffractive component.
9. Primordial magnetic fluctuations with zero frequency [58] may have interacted with neutrinos before neutrino detachment from the hot neutrino plasma(>GeV temperature) during the Big Bang. Such signatures may be carried by neutrinos (which is now relic neutrinos).

Summary and the future prospect
We have found that the photon interaction expressed by the total derivative Eqs. (13) or (86), which are derived from the triangle diagram in the standard model, causes the unusual transitions characterized by the time-independent probability. The interaction Lagrangian of this form does not give rise to any physical effect in classical physics, because the equation of motion is not affected. In quantum mechanics, this assertion is correct for the transition rate Γ 0 . However, this does not apply to the diffractive term P (d) , which manifests the wave characteristics of the initial and final waves. Our results show that P (d) is relevant to experiments and important in understanding many phenomena in nature.
The neutral particles do not interact with the photon in classical mechanics. In quantum field theory, the vacuum fluctuation expressed by triangle diagram gives the effective interaction to the neutral particles such as 1 + meson → γγ and ν +γ(B, E) → ν +γ. However, they have vanishing rates due to Landau-Yang-Gell-Mann's theorem. P (d) does not vanish, nevertheless, and holds unusual properties such as the violation of the kinetic-energy conservation and that of Lorentz invariance. Furthermore the magnitudes become comparable to or even larger than the normal weak processes. Accordingly, the two photon or two gluon decays of the neutral axial vector mesons composed of a pair of electrons or quarks come to have finite decay probabilities. They will be tested in experiments. The neutrino photon processes, which have been ignored, also have finite probabilities from P (d) . It will be interesting to observe the neutrino photon processes directly using electric or magnetic fields, or laser and neutrino beams in various energy regions. The diffractive probability P (d) would be also important for understanding the wide neutrino processes in earth, star, astronomy, and cosmology.
The diffractive probability P (d) is caused by the overlap of wave functions of the parent and decay products, which makes the interaction energy finite and the kinetic energy vary. Consequently the final state has continuous spectrum of the kinetic energy and possesses a wave nature unique to the waves. The unique feature of P (d) , i.e., independence of on the time-interval T , shows that the number of parents, which decrease like e −Γt in the normal decay, is constant now. The state of parent and daughter is expressed by the quasi-stationary states, which is expressed by the superposition of different energies and different from the normal stationary state of the form e −iEt ψ(x i ). The probability of the events that the neutrinos or photons are detected is computed with S[T ] that satisfies the boundary condition of the physical processes. Applying S[T ] we obtained the results that can be compared with experiments. The pattern of the probability is determined by the difference of angular velocities, ω = ω E − ω dB , where ω E = E ν / and ω dB = c| p ν |/ . The quantity ω takes the extremely small value m 2 c 4 /(2E ) for light particles such as neutrinos or photons [49] in matter. Consequently, the diffraction term becomes finite in the macroscopic spatial region of r ≤ 2πE c m 2 c 4 . This allows us to introduce a new class of experimental measurement possibilities of the deployment of photons to detect weakly interacting particles such as neutrinos. Because the modern technology on the electric and magnetic fields and laser, a large number of coherent photons are possible and the effects we have derived may have important implications in detecting and enhancing the measurement of neutrinos with photons. We see a variety of detectability opportunities that have eluded attention till now. These happen either with high energy neutrinos such as from cosmic rays and accelerators; or in high fields (such as intense laser and strong magnetic fields ). In the latter examples, neutrinos from reactors, accelerators, the sun, supernovae, thunder clouds and even polar ice may be detected with enhanced probabilities using intense lasers. We also mention the probability estimate for the primordial relic neutrinos and embedded information in them.
where the integrands satisfy g(t + , x + = ±∞) = 0, f (t − , x − = ±∞) = 0 (A.5) is made with the change of variables. Due to Eq. (A.5), I would vanish if the integration region were from −∞ to +∞. We write I with variables x + = x 1 +x 2 2 , x − = x 1 − x 2 , and have I = I 1 − I 2 , (A.6) g(x + )), where λ + = x 2 + . It is noted that x + is integrated in the restricted region but x − is integrated in whole region. Thus g(x + ))| x + =x +,min , (A.7) where the functional form g(x + ) = g(x + − vt + ) was used. I is computed with the slope of g(x + ) at the boundary x + = t + . We apply this method for computing four dimensional integrals.

Appendix C
In the case of the low energy region the photon has no effective mass, but is expressed by the index of refraction very close to unity in dilute gas, n = 1 + ǫ. Consequently the integrand in P (d) is proportional to 2 Pγ ǫ . and the Lagrangian density of electric and magnetic fields The Lagrangian density of electronic fields is, .18) and that of QED is L QED = L e + L EM + +ecA 0ψ γ 0 ψ − ecA lψ γ l ψ. The spinor is normalized as s u( k, s)ū( k, s) = γ · k +ω 0 2ω 0 , (A. 27) and the light-cone singularity is where the less-singular and regular terms are in · · · . The action S = dx(L QED +ψ(x)γ µ (1 − γ 5 )ψ(x)J µ (x)), (A.29) governs the dynamics of electron, photon and neutrino. Integrating ψ(x) andψ(x), we have Det(D + J µ γ µ (1 − γ 5 )) and the effective action between neutrino and photon