Research Article Volume 2 Issue 5
The High Energy Accelerator Organization (KEK), Japan
Correspondence: Yoshimasa Kurihara, The High Energy Accelerator Organization (KEK), Oho 1?1, Tsukuba, Ibaraki 305?0801, Japan, Tel 8129 8796 088
Received: August 28, 2018 | Published: October 12, 2018
Citation: Kurihara Y. Electroweak symmetry braking and quantization of the higgs field in early universe. Phys Astron Int J. 2018;2(5):468-474. DOI: 10.15406/paij.2018.02.00126
In this paper, we propose a novel method to treat the electroweak symmetry braking. In this method, a conformal metric is employed and the Higgs field is scaled owing to the conformal function and the mass parameters of the quadratic term of the Higgs potential has a time dependence through the conformal function, and it induces the phase transition. Quantization of the Higgs field is induced associated with the canonical quantization of general relativity. The cosmic inflation and the electroweak phase–transition are discussed in a framework of the scaled field. The Friedmann equations are numerically solved and an example of a possible solution to match with the cosmic inflation scenario is given in this research.
Keywords: electroweak symmetry, conformal metric, higgs potential, cosmic inflation, phase transition, friedmann equations
After the discovery of the Higgs boson1,2 in 2012, the standard theory of particle physics (SM) is established as the canon of a fundamental physics. According to the SM, the electroweak symmetry is spontaneously broken owing to the Higgs mechanism,3–5 and the current universe is considered to be filled with the Higgs filed which has a finite vacuum expectation–value. On the other hand, the electroweak symmetry is expected to be in the unbroken phase in the early universe before the cosmic inflation. A standard scenario of the big–bang cosmology is that the electroweak phase–transition from the unbroken to the broken phases might occur during (or at the beginning of) the cosmic inflation, and the universe was re–heated after the cosmic inflation, and then the big–bang started.
The cosmic inflation was proposed to solve the flatness and horizon problems by several authors independently6–9 in 1980. In these “old” models, the cosmic inflation is induced by the Higgs filed (or some scalar filed which is referred to as the inflaton field), and it is terminated when the electroweak phase–transition of the first–kind occurred. These models are suffered by the vacuum–bubble problem which destroyed isotropy of the universe. In 1982, the “new” inflation models10–12 are proposed which utilize the phase–transition of second–kind to avoid the vacuum–bubble problem. The inflation terminated moderately in the “new” models. Although these “new’ models can solve the vacuum–bubble problem, they require a fine tunning of initial parameters to realize the cosmic inflation for a enough time duration. Yet other models of the cosmic inflation was proposed such as the chaotic inflation,13 the Higgs inflation with non–minimal coupling to gravity14 and so on.
In any inflation scenario, the electroweak phase–transition is critical to induce and terminate the cosmic inflation. However, a widely accepted mechanism to induce the symmetry braking is not established yet. The early scenario15 of the symmetry braking using higher–order radiative corrections is now excluded in precise measurements of the SM parameters. The radiative braking scenario is re–examined and concluded that it is still viable if an additional scalar field is introduced.16 This scenario is intensively investigated in literatures.17,18 Recently, this scenario is extended19 using the classically conformal ℬ−ℒ extension of the SM.20
We propose a new and novel mechanism of the electroweak symmetry braking in this study. In the proposed method, which is referred to as the “scaled scalar–field method”, a conformal metric is employed, and the Higgs field is scaled owing to the conformal function. Quantization of the Higgs field is induced due to the canonical quantization of general relativity. In consequence, a mass parameter of the Higgs field (a quadratic term of the Higgs potential) has the time dependence through the conformal function, and it causes the phase transition. The scaled scalar–field method is introduced in section 2 after a brief explanation of our geometrical setups. The Friedmann equation which governs the cosmic inflation is formulated using the scaled scalar–field method in section 3.1. After some appropriate approximation, the Friedmann equations are numerically solved and an existence of a possible solution to match with the cosmic inflation scenario is shown in section 3.2. A summary of the method and consequences on the inflation scenario is provided in section 4.
Geometrical setups
A scalar field defined on a four dimensional space time manifold ℳ with a GL(1,3) symmetry is considered in this study. First, classical general relativity and the scalar field defined on ℳ are summarized using a vierbein formalism. The formalism and terminology in this study follow our previous works.21–23 At each point on ℳ , a local Lorentz manifold M with a Poincaé symmetry ISO(1,3)=SO(1,3)⋉T4 is associated, where T4 is a four–dimensional translation group. On a open neighborhood around any point x∈U⊂ℳ , a trivial frame vector is expressed as xμ , and a trivial vector bundle (frame bundle) can be introduced. An orthonormal basis of ∂μ in Tℳ and dxμ in T*ℳ are also introduced. A short–hand notations of ∂μ=∂/∂xμ are used through out this study. The Einstein’s equivalent theorem insists an existence of an isomorphism ℳ→M at any point in ℳ . A metric tensor g•• in ℳ is mapped to η•• in M using a vierbein ℰaμ as gμν=ηabℰaμℰbν . An orthogonal basis in T*M and TM are respectively expressed as dxa=ℰaμdxμ and ∂a=ℰμa∂μ using a vierbein and its inverse. The Einstein convention for repeated indices is used though out this study. In addition, Greek and Roman indices are used for a coordinate on ℳ and M , respectively. A local Lorentz metric and the Levi Civita tensor are respectively defined as η••=diag(1,−1,−1,−1) and ε•••• with ε0123=1 . Dummy Roman–indices are often abbreviate to dots (or asterisks), when the index pairing of the Einstein convention is obvious, such as ηabℰaμℰbν=η⋅⋅ℰ⋅μℰ⋅ν . When multiple dots appear in an expression, pairing must be a left–to–right order at both upper and lower indices, e.g. a⋅⋅∧b⋅⋅=aab∧bab . A principal connection of the fiber bundle so(1,3)→M is represented as ωaμb , which is referred to as the spin connection. The spin connection satisfies a metric compatibility condition as ωaμ⋅η⋅b=ωabμ=−ωbaμ . A vierbein for ea=ℰaμdxμ and a GL(1,3) invariant volume form v=ε⋅⋅⋅⋅e⋅∧e⋅∧e⋅∧e⋅/4! are introduced. Similarly, the three–dimensional volume form and two–dimensional surface form are also introduced as Va=εa⋯e⋅∧e⋅∧e⋅/3! and Sab=εab⋅⋅e⋅∧e⋅/2 , respectively. The volume form Va is a three–dimensional volume perpendicular to ea , and the surface form Sab is a two–dimensional plane perpendicular to both ea and eb . Fraktur letters are used for differential forms. A unit of c=1 is used while the reduced Planck constant ℏ and Newtonian gravitational constant G (or the Einstein constant κ=4πG in our convention) written explicitly. In this units, there are two physical dimensions, the length and mass dimensions, which are denoted as L and M , respectively.
The Lagrangian for pure gravity without the cosmological term and matter fields is expressed as
LG=12S⋅⋅∧ℜ⋅⋅, (1)
ℜab=dwab+wa⋅∧w⋅b, (2)
where w•• is the spin one–form, which is defined as wab=ωaμ⋅η⋅bdxμ . A two form ℜ•• is referred to as the curvature form, that is a rank– 2 Lorentz tensor on M .
A Lagrangian of a scalar field φ(x) can be expressed in the vierbein formalism as
LS=1213!S⋅⋅∧(η∗∗ι∗s⋅∧ι∗s⋅−V(φ)e⋅∧e⋅), (3)
where V(φ) is a potential energy. A scalar–field two–form s• is defined as
sa=dφ∧ea=(∂⋅φ)e⋅∧ea. (4)
Here ιa=ιξa is a contraction operator with respect to the trivial coordinate–vector field ξa=ηa⋅ℰμ⋅∂μ . A physical dimension of the scalar field is set as L−1 in this study. In consequence, a Lagrangian form has null physical dimension as
[φ]=L−1→{[ι•s•]=L−1[V]=L−4→[LS]=1, (5)
where [•] shows the physical dimension of a quantity • . On the other hand, the gravitational Lagrangian has a length square dimension [LG]=L2 . An action integral is defined as
ℐ=∫(1κℏLG+2LS). (6)
Physical constants in front of the gravitational and scalar–field Lagrangian are chosen to adjust the physical dimension of the action integral to be null. Owing to require a stationary condition on a variation of the action integral with respect to the vierbein form δeℐ=0 , one can obtain the Euler–Lagrange equation as
12εa⋅⋅⋅ℜ⋅⋅∧e⋅=−κℏta, (7)
where ta is the energy–momentum three–form of the scalar field, which can be represented as
ta=−13!εa⋯(ι∗s⋅)∧(ι∗s⋅)∧e⋅+V(φ)Va. (8)
Here, rising and lowering Roman indices are performed using a Lorentz metric tensor, e.g. ιasb=ηa⋅ι⋅eb . A torsionless condition and Klein–Gordon equation can be obtained from δ(w,dw)ℐ=0 and δ(φ,dφ)ℐ=0 , respectively.
The scalar–filed Lagrangian LS given in (3) can be expressed using a trivial frame vector in M or ℳ as
(3)=v[12η⋅⋅∂⋅φ∂⋅φ−V(φ)],
=√−gdx0∧dx1∧dx2∧dx3[12gμ1μ2∂μ1φ∂μ2φ−V(φ)], (9)
where g is a determinant of a metric tensor g•• . Here a relation ιaeb=δba is used. The Einstein equation (7) can be expressed using components of a trivial basis on M as
Rab−12ηabR=−2κℏTab, (10)
Tab=∂aφ∂bφ−12ηab∂⋅φ∂⋅φ+ηabV, (11)
where R•• and R are the Ricci and scalar curvature, respectively. An energy–momentum tensor T•• is defined through the relation ta=V⋅T⋅a .
Conformal metric and scaled field method
The Friedmann–Lemaître–Robertson–Walker (FLRW) metric24–27 is considered in the homogeneous and isotropic universe using a coordinate (t,r,θ,ϕ) on ℳ such as;
ds2=dt2−Ω2(t)(f−2(r)dr2+r2dθ2+r2sin2θdϕ2), (12)
where f2(r)=1−Kr2 , and K is a constant related the space–time curvature. This metric can be expressed using a conformal time τ(t)=∫d˜t/Ω(˜t) as
(12)=Ω2(τ)(dτ2−f−2(r)dr2−r2dθ2−r2sin2θdϕ2). (13)
From the metric (13) and the torsion–less condition, the spin form can be obtained as
w••=(0Hf−1drHrdθHrsinθdϕ0−fdθ−fsinθdϕ0−cosθdϕ0), (14)
where H=((dΩ/dτ)/Ω)(τ) . The lower half is omitted because it is obvious due to antisymmetry of the spin form. From (14) and the surface form, we can obtain the classical Hamiltonian as22
ℌFLRW=−12w⋅∗∧w∗⋅∧S⋅⋅=Ω2[3H2−(fr)2]dτ∧drf∧(rdθ)∧(rsinθdϕ), (15)
and the Liouville form as
12w⋅⋅∧S⋅⋅=Ω2[−cotθrdτ∧drf∧(rsinθdϕ)
+2frdτ∧(rdθ)∧(rsinθdϕ)+3Hdrf∧(rdθ)∧(rsinθdϕ)]. (16)
A canonical quantization of general relativity requires the commutation relation between the spin form and surface from as22
[ˆwa1a2(x),ˆSb1b2(y)]=−iℏGδ(4)(x−y)δ[a1b1δa2]b2, (17)
and otherwise zero, where δ[a1b1δa2]b2=δa1b1δa2b2−δa2b1δa1b2 . Here, ˆw•• and ˆS•• are operators respectively corresponding to the spin and surface forms, which are formally described as
{ˆw••=w••,ˆS••=iℏGδδw••. (18)
The commutation relation (17) can be represented using terms of the FLRM metric. The metric tensor is a functional of two functions Ω(τ) and f(r) other than the integration measure. On the other hand, the Liouville–form (16) includes a derivative ˙Ω , and it does not have a term df(r)/dr . Therefore, quantization of the system can be performed by replacing a scale faction Ω by the corresponding operator as ˆΩ . The conformal–time derivative1 ˙Ω=dΩ/dτ is included only in the last term of (16). Thus, a non–zero component of the commutation relation can be obtained from (16) and (17) as
[ˆw⋅⋅(r),ˆS⋅⋅(r′)]|τ=τ′=2⋅3[ˆΩ(τ),ˆ˙Ω(τ)]d3r' (19)
where r=(τ,r) , r′=(τ′,r') and r , r' are three–dimensional spacial vectors. The three–dimensional integration can be expressed as
d3r'=dr′f(r′)∧(r′dθ)∧(r′sinθdϕ). (20)
On the other hand, the tight hand side of (17) is represented as
(17)=−i(2ℏG)δ3(r−r'). (21)
Therefore, the commutation relation is obtained as
[ˆΩ(τ),ˆ˙Ω(τ)]d3r'=−iℏG3δ3(r−r'). (22)
This can be understood as the equal–time commutation relation of the conformal metric. While this commutation relation is obtained from the commutation relation (17), one can obtain the same result based on quantization by Nakanishi28 as shown in Appendix A.
The scalar field can be defined in the conformal metric with K=0 in the Cartesian coordinate for three–dimensional space
ds2=Ω2(τ)(dτ2−(dx1)2−(dx2)2−(dx3)2) (23)
Instead of a polar–coordinate because virbeins are now independent of r . The action integral of a scalar Lagrangian (9) in the conformal FLRW metric can be obtained as29
ℐs=∫LS=∫dτ∧dx1∧dx2∧dx3[12η⋅⋅(∂⋅χ)(∂⋅χ)+12¨ΩΩχ2−V(χ)] (24)
using a local conformal coordinate, where χ(x)=Ω(τ)φ(x) . We note that √−det{g••}=Ω4 for the conformal metric (23). For further discussions, we specify the potential energy as the Higgs–type field as;
V(χ)=Ω4V(φ)=−12(Ωμℏ)2χ2+λ4χ4, (25)
where μ and λ>0 is real constants. Here the physical dimension of each term is given as;
[μ]=M,[ℏ]=LM→[μℏ]=L−1,[Ωμℏχ]=L−2. (26)
This field χ is provided owing to scale the scalar field φ by the conformal function Ω(τ) , which is referred to as the scaled scalar field (SSF). While a quartic term has no corrections, a quadratic term receives corrections from the conformal function and its second derivative. According to a standard procedure for field quantization (canonical quantization), the field and its conjugate momentum are replaced by corresponding operators. Here, the non–zero component of equal–time commutation relations of the SSF are naturally introduced from (21) as;
[ˆχ(τ,x),ˆ˙χ(τ,x')]=ˆφ(τ,x)ˆφ(τ,x')[ˆΩ(τ),ˆ˙Ω(τ)],
=iℏG3|ˆφ(x)|2δ(3)(x−x')→iδ(3)(x−x'), (27)
where ˆφ(x) is a scaler field operator. The overall factor ℏG/3 can absorbed by re–definition of the scalar field. Although the ˆφ is an operator, it is assumed to commutes each other as [ˆφ,ˆ˙φ]=0 . Non–commutativity of χ and ˙χ is induced by that of Ω and ˙Ω .
Since the commutation relation (27) for the scaled–field operator ˆχ is the same as the standard Klein–Gordon field–operator, the standard procedure of the field quantization in a momentum space using the Fock Hilbert–space can be performed as usual. We note that, in the SSF formalism, the commutation relation is required only on the space–time metric (vierbein). Any observers in the local space time cannot observe the scalar field independently form the space time metric. This is one of a realization of the concept given in Kurihara30 such that “Classical mechanics in the stochastic space is equivalent to quantum mechanics on the standard space time manifold”. According to this concept, only the vierbein is quantized using the quantum commutation relations with keeping the scalar field classical. Even if the local observer is in the flat space–time, one may observe the SSF which is quantized due to the quantum space time. From equations (24) and (25), the effective potential can be obtained as;
˜V(χ)=−12(˜μℏ)2χ2+λ4χ4, (28)
where
(˜μℏ)2=¨ΩΩ+(Ωμℏ)2. (29)
An effective mass of the SSF has time dependence through the relation (29), and the electroweak phase–transition from unbroken to broken phases can occur through it. Although the SSF mass was very small, which corresponds to the unbroken phase, owing to a small value of Ω and ¨Ω/Ω at the early universe, it can be large as the same as that in the current universe after the inflation due to the second term with a large value of Ω . In the current universe, the first term of (29) is negligibly small compared with the second term, and the time dependence of the mass term through Ω(t) is much smaller than the current accuracy of Higgs mass measurements.
A vacuum expectation value v0=˜μ/√λ and higgs mass mh=√2˜μ can be extracted from the effective potential. The Klein–Gordon equation obtained from the action (24) with respect to the scaled field χ can be expressed as
¨χ−Δχ+∂˜V∂χ=0, (30)
where Δ=∑i=1,2,3∂i∂i . For the uniform and isotoropic universe, one can set Δχ=0 . In this case, the energy–momentum tensor (11) is represented as
T00=12˙χ2+˜V, (31)
Tii=12˙χ2−˜V,i=1,2,3(not summed), (32)
and other components are zero. This energy–momentum tensor can be interpreted as a density ( ρ ) and pressure ( p ) of a perfect fluid as T00=ρ and Tii=p , respectively.
Friedmann equation
The Einstein equation (10) for perfect fluid with the conformal metric (13) can be expressed as follows:
3(H2+K)=+κℏρ, (33)
2˙H+H2+K=−κℏp, (34)
where the curvature K is put buck in this subsection. These equations, refereed to as the Friedmann equations, can be rearranged as
H2=2κℏ3ρ−K, (35)
˙H=−κℏ3(ρ+3p). (36)
When the density and pressure are provided from the SSF, the Friedmann equations have a form
H2=+2κℏ3(12˙χ2+˜V)−K, (37)
˙H=−2κℏ3(˙χ2−˜V). (38)
If appropriate initial conditions are given, the history of the universe can be obtained by solving equations (30), (37) and (38), simultaneously. Further simplification is possible in this case as follows: A conformal–time derivative of (37) provides a equation
dH2dτ=2κℏ3(¨χ+∂˜V∂χ)˙χ−κ3ℏd˜μ2dτχ2,
=2κℏ3Δχ˙χ−κ3ℏd˜μ2dτχ2, (39)
where the Klein–Gordon equation (30) is used. Under the assumption that the SSF is the uniform and isotoropic, a spacial derivative are set to zero again as Δχ=0 . Therefore, the equation (39) can be expressed as
dH2dτ=−κ3ℏd˜μ2dτχ2. (40)
Due to definitions of ˜μ2 in (29), its conformal time derivative is expressed as
d˜μ2dτ=ℏ2ddτ(¨ΩΩ)+2Ω˙Ωμ2. (41)
Under the assume that higher derivatives of the scale function Ω are much smaller than both the scale function itself and first derivative of that, the first term on the righthand side can be ignored compared with the second term. The validity of this assumption will be confirmed later in this section. In this case, the equation (40) can be approximated as
˙ΩΩ(τ)=√κ3ℏμΩ(τ)χ(τ), (42)
where H(τ)>0 is assumed through any τ . This equation can be express using the proper time t using a relation Ω(τ)dτ=dt as
dΩ(t)dt=√κ3ℏμΩ(t)χ(t). (43)
The second Friedmann equation (38) can be approximated under the same assumptions as
dΩ(t)dt=√κ3[Ω(t)2(2(dχ(t)dt)2+(μℏχ(t))2)−λ2χ(t)4]12. (44)
Above two Friemann equations (43) and (44) must be consistent each other within the approximation. It leads a following differential equation
˙χ(τ)=±√λ2χ(τ)2, (45)
with respect to the conformal time τ . This equation can be solved as
χ(τ)=(∓√λ2τ+1χ0)−1, (46)
where χ0=χ(τ=0) . This solution does not give any oscillating fields because of the requirement Δχ=0 .
Numerical calculations
Next the scaling function for the cosmic inflation era is treated in this section. During the cosmic inflation, a condition dΩ/dt≫dχ/dt is expected naturally. Due to the definition of the SSF, this condition can be satisfied if φ≪1 and dφ/dt≪(dΩ/dt)/Ω during the cosmic inflation. A validity of this assumption will be discussed later in this section. The Friedmann equation (43) can be solved under the assumption as;
Ω(t)=Ω0exp(H0t), (47)
where H0=√κ/3ℏμχ0 with an initial condition of Ω(t=0)=Ω0 . Due to the assumption above, the SSF stays constant at χ(0<t<te)=χ0 during the inflation. An inflation starting time is set to t=0 . An initial value of the SSF may be given by a quantum fluctuation of the field, which is naturally expected to be an order unity. Here the initial value χ0=1 in the Planck units is assumed. On the other hand, the Higgs potential parameters, μ and λ , at the beginning of the universe are set to be the same as the current universe at the tree level for simplicity. From the recent measurement of the Higgs mass,31 the quadratic term can be obtained as Ωμ/ℏ=90 and . Therefore, the initial value is obtained in the Planck units as , where is the Planck mass in the natural units. In order to make the inflation scenario working well, the e–folding number defined as must be greater than or around 70. When the normalization and the e–folding number are required, the initial value of the scaling function can be obtained as . The inflation termination–time can be obtained by solving the such as second. In our scenario, a speed of changing a vacuum expectation value was very high, and thus, the SSF vacuum is staying at the origin for a short period. Then, delayed explosion of the SSF field caught the vacuum expectation value up and the inflation was terminated. Therefore, the electroweak phase–transition must be the second kind.
The evolution of the scaling function can be fixed completely using above parameters under the approximations. Before a starting time of the inflation, the first term of (29) was much smaller than the second term. At that period, the potential energy of the SSF had a minimum at as shown in Figure 1 . During the inflation, the SSF had a finite vacuum expectation value (Figure 1 ). At the end of the inflation, the vacuum expectation value arrived at the same value as that in the current universe (Figure 1 ). Although the vacuum expectation value stayed almost constant at the beginning inflation, it grew very rapidly after the second term of (29) became dominant in the potential energy as shown in Figure 2. The inflation was terminated when the SSF arrived at the vacuum expectation value, and it fixed the e–folding number. A precise value of the inflation ending time must be evaluated by solving equations (43) and (45) (or equivalently (46)), simultaneously. We note that the relation between the conformal and proper times can be fixed only after the solution of the scaling function is obtained. One cannot solve equations analytically without the approximation that the SFF is constant during the inflation. When the inflation duration changed around its nominal value, the e–folding number varied about as shown in Figure 3. On the other hand, the same variation for the duration of the inflation affects the vacuum expectation value about to , as shown in Figure 2. A tolerance of the inflation duration is rather narrow to realize a current observed value of the vacuum expectation value.
Figure 1 The potential energy of the SSF before the inflation (solid line), during the inflation (dotted line), and at the end of the inflation (dashed line).
Figure 2 The potential–energy term of the SSF normalized using the SM value in the current universe. If we require that the cosmic inflation terminated when the potential–energy term of the SSF arrived at the current value , a duration of the in ation is second.
Figure 3 The e–folding number different from nominal value is shown with respect to the duration time of the ination. When the ination duration changed around its nominal value, the e–folding number varied about .
The solution (47) is obtained under the assumption of . The validity of this assumption during the cosmic inflation must be confirmed. If the time dependence of the SSF is put back in the solution as , the time derivative of the scaling function becomes
(48)
In calculations for the cosmic inflation, the time derivative term in the right hand side is neglected. The validity of this approximation can be examined using the equation of motion. The time evolution of the SSF is governed by the equation (46) with respect to the conformal time. For a conversion from the conformal time to the proper time, the approximated solution of (47) is used for a numerical calculation. A numerical result of a time evolution of the tern is shown in Figure 4. It is shown that that term is less than unity during the cosmic inflation, and thus the effect on the result is a factor of two on at most.
We propose a novel method to treat the electroweak symmetry braking, which is named the scaled scalar–field method. In this method, a conformal metric is employed and the Higgs field is scaled owing to the conformal function. In consequence, a mass parameter of the Higgs field (a quadratic term of the Higgs potential) has the time dependence through the conformal function, and it causes the phase transition. Quantization of the Higgs field is induced associated with the canonical quantization of general relativity. In the context of the scaled field, only the vierbein is quantized owing to the quantum commutation relations with keeping the scalar field classical.
The cosmic inflation and the electroweak phase–transition are investigated in a framework of the scaled field. The Friedmann equations and their appropriate approximations are provided using the scaled field method. The Friedmann equations are numerically solved and an example of a possible solution to match with the cosmic inflation scenario is shown. The electroweak phase–transition induced by the scaled field is the second kind, and thus, the fine–tuning problem is still exists.
I appreciate the kind hospitality of all members of the theory group of Nikhef, particularly Prof. J. Vermaseren and Prof. E. Laenen. A major part of this study has been conducted during my stay at Nikhef in 2017. In addition, I would like to thank Dr. Y. Sugiyama for his continuous encouragement and fruitful discussion.
Author declares there is no conflict of interest.
©2018 Kurihara. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.