Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Physics

formerly Central European Journal of Physics

Editor-in-Chief: Seidel, Sally

Managing Editor: Lesna-Szreter, Paulina

1 Issue per year


IMPACT FACTOR 2017: 0.755
5-year IMPACT FACTOR: 0.820

CiteScore 2017: 0.83

SCImago Journal Rank (SJR) 2017: 0.241
Source Normalized Impact per Paper (SNIP) 2017: 0.537

ICV 2017: 162.45

Open Access
Online
ISSN
2391-5471
See all formats and pricing
More options …
Volume 16, Issue 1

Issues

Volume 13 (2015)

Super-structured photonic crystal fiber Bragg grating biosensor image model based on sparse matrix

Bo Li
  • Corresponding author
  • East China Jiaotong University, Information Engineering School, Nanchang, Jiangxi, China, 330013
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Jayanth Kanugo
  • Computer Science Department, University of Missouri-Columbia, Columbia, United States of America
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2018-11-30 | DOI: https://doi.org/10.1515/phys-2018-0091

Abstract

Fiber grating (FBG) is an important optical device of fiber, which is widely used in optical fiber communication and sensing. At the present stage, the fiber grating is almost prepared in the static state of the fiber, and then the grating is welded into the fiber grating array or network. With the continuous improvement of the application effect of fiber grating, has become one of the most promising, representative and fastest developing fiber passive devices. In this paper, fiber grating is classified according to the refractive index distribution of grating axis. The central wavelength of Bragg fiber grating is modulated by using external parameters (temperature or stress strain). The wave equation of insulating medium is obtained by the reflection characteristic and coupling mode theory of fiber grating. Using sparse matrix model of nerve action potential signal with wavelet decomposition layers, nerve action potential signal reconstruction of the relative error between the value and the original value contrast found that reconstructed and original signals are very close. Good results have been obtained for the sampling reconstruction of the filtered high signal-to-noise ratio neural action potential signal. Researchers have conducted extensive and in-depth research on fiber grating sensing technology, and achieved gratifying results. But with the of the engineering application of technical requirements, they need real-time monitoring. Due to the cross sensitivities of fiber grating, it became the bottleneck of multiple parameter measurement.

Keywords: Sparse matrix; Reflection characteristics; PCF Bragg grating sensor; Fiber grating

PACS: 02.30.Cj; 02.50.Fz; 42.70.Qs

1 Introduction

The fiber grating uses the photosensitivity of the fiber material to form a spatial phase grating in the core [1, 2, 3].

It has a set of vertical reflection peaks and it can be used as a comb filter. Besides, it has potential applications in multi-wavelength fiber lasers. The essence of the fiber grating is to form in the core (wavelength characteristics of Bragg scattering using a spatial phase grating) a narrow–band (projection or reflection) filter. Fiber gratings are made using photosensitivity in optical fibers [4, 5]. The so-called photosensitivity in an optical fiber means that the refractive index of the optical fiber will change with the spatial distribution of the intensity of the light when the laser passes through the doped fiber. The essence of the spatial phase grating formed in the core is the formation of a narrowband filter or mirror in the core [6]. With this feature, many unique fiber optic devices can be manufactured. These devices have a wide range of reflection bandwidth, small additional loss, small size, easy to couple with the optical fiber, compatible with other optical devices into one, without the impact of environmental dust and a series of excellent performance.

With the increasing application of fiber gratings, the types of fiber gratings are also increasing. According to the distribution of refractive index along the axial direction of the grating, the ultraviolet-written fiber grating can be divided into a uniform fiber grating and a non-uniform fiber grating [7]. The uniform fiber grating refers to a fiber grating in which the core refractive index variation amplitude and refractive index variation period (also referred to as the period of the fiber grating) are all constant along the fiber axis, such as a uniform fiber Bragg grating (period of refractive index change is generally of 0.1 um) and uniform long-period fiber grating (period of refractive index change is generally the order of 100 um); Non-uniform fiber grating refers to the core refractive index varies in amplitude or refractive index varies along the fiber’s axial direction. In 1978, K. O. Hill et al. first used the standing wave writing method to fabricate the first fiber grating in erbium-doped fiber. After more than 20 years of development, it has broad application prospects in fields such as optical fiber communication and optical fiber sensing [8, 9, 10]. With the continuous improvement of the manufacturing technology of fiber gratings, the photosensitivity of optical fibers has been gradually improved, and various kinds of special gratings have come out one after another. Some applications of fiber gratings have reached the commercialization level. With increasing application results, FBG has become one of the most promising, most representative and fastest-developing fiber-optic passive devices currently available. FBG-based sensing technologies have emerged in recent years, fully demonstrating advantages. However, FBG sensing technology is still in a stage of rapid development. It can be predicted that with the commercialization of FBG sensors, continuous improvement of performance, fiber gratings will show great vitality in the field of sensing, and play an important and irreplaceable role in national defense and national economic construction. Significance statement: This study the super-structured photonic crystal fiber bragg grating biosensor image model based on sparse matrix that can be beneficial for the field of optical fiber sensing technology with its unique advantages.

2 Methods and principles

2.1 Sparse matrix model of neural action potential signals

Nerve action potential is a record of when a neuron is active [11]. It contains the most critical information about neuronal activity. Similar to other bioelectrical signals, neural action potentials can be sparse under orthogonal wavelet basis. Using this feature, a sparse matrix of neural signals can be constructed in the wavelet domain.

The construction of the sparse matrix is based on the Mallat algorithm of the discrete wavelet transform (DWT). The wavelet filter is used to decompose the signal [12]. The specific process can be expressed by the following formula (1):

A0f(t)=f(t)Ajf(t)=kH(2tk)Aj1f(t)Djf(t)=kG(2tk)Aj1f(t)(1)

In equation (1), t is the discrete time sequence number, t = 1, 2, L, N; f(t) is the original signal; j is the number of layers, j = 1, 2, L, J, where J = log2N. H, G is the wavelet decomposition filter in the time domain, which is actually the filter coefficient; Aj Ai is the wavelet coefficient of the signal f(t) at the approximate part (i.e., the low frequency part) of the jth layer, and D is the detail of the signal f(t) at the jth layer.

The specific decomposition process can be represented by Figure 1 and Figure 2.

One-layer wavelet decomposition
Figure 1

One-layer wavelet decomposition

Multilayer wavelet decomposition
Figure 2

Multilayer wavelet decomposition

2.2 Optical mask distribution expression

Light through the mask, the light field distribution can be expressed as the following formula (2):

E(x)=1t1(x)t2(x)(2)

Here, t1(x) and t2(x) are the transmittances of light transmitted through the amplitude mask and the phase mask, respectively. They can be expressed as the following formula (3):

t1(x)=0d/2xa/21a/2xa/20a/2xd/2t2(x)=expjϕ1x+(2J+1)Λ/2<Λ/2expjϕ2xJΛ/2<Λ/2(3)

Where d is the period of the amplitude mask; a is the width of the tooth of the amplitude mask; L is the total length of the grating; J = 0, ±1, ±2, · · · , 𝛬 is the period of the phase mask. When the phase is produced, expanding the formula (1) with Fourier series, the formula (1) can be expressed as the following formula (4):

E(x)=n=+Anexpj2πnxdn=+Cmexpj2πmxΛ(4)

Where, An=d/2d/2t1(x)expj2πnxddx;Cm=Λ/2Λ/2t2(x)expj2πmxddx

Equation (4) is the expression of the light field distribution of the ultraviolet laser after the periodic square wave modulation phase mask.

2.3 Principle and calculation of fiber Bragg grating reflection characteristics

Fiber Bragg grating is a waveguide with periodically varying refractive index. From coupled-mode theory, it can be seen that when the waveguide is disturbed periodically, it will lead to the coupling of propagation modes. Fiber Bragg grating sensing principle diagram was shown in Figure 3.

Fiber Bragg grating sensing principle diagram
Figure 3

Fiber Bragg grating sensing principle diagram

Maxwell’s equations can be expressed as the following Equations (5) and (6):

×H¯=J¯+tε0E¯+P¯(5)×E¯=t(μH¯)(6)

For optical waveguides, there is no free charge, so J¯= 0. Equations (5) and (6), the wave equation of insulating medium can be, the expression can be written as the following Equation (7):

2E¯(r¯,t)=με02E¯(r¯,t)t2+μ2t2P¯(r¯,t)(7)

Due to the non-uniformity of the medium fluctuation, it can be considered as perturbation, and the polarization intensity can be expressed as the following Equation (8):

P¯(r¯,t)=[ε(r¯)ε0]E¯(r¯,t)+P¯Perturbation(r¯,t)(8)

Substituting (8) into (7), we can obtain the following equation (9):

2E¯yμε2E¯yt2=μ2t2[P¯Perturbation(r¯,t)](9)

E x and E z are similar. For the TE mode, if you omit the guided mode and the radiation mode coupling, there is the following equation (10):

Ey(r,t)=12mAm(z)Ey(m)(x)ei(ωtβmz)+c.c.(10)

where m is the number of the inherent mode, c.c. is the complex conjugate, and the inherent mode field satisfies the non-disrupting sex wave equation. The non-disrupting sex wave equation can be expressed as following:

(2x2βm2)Ey(m)(r)+ω2με(r)Ey(m)(r)=0(11)

where ε(r) = ε0n2(r), n(r) is the medium refractive index. Substituting (10) into (11), there is the following Equation (12):

eiωtmAm2βm2Ey(m)+2Ey(m)x2+ω2με(r)Ey(m)eiβmz+122iβmdAmdz+d2Amdz2Ey(m)eiβmz+c.c.=μ2t2[PPerturbation(r,t)](12)

The sum of the first three terms of the equation is equal to zero, and the slow amplitude changes are similar. Omitting the second derivative, because the second derivative satisfies the following Equation (13), we get:

d2Amdz2βmdAmdz(13)

Substituting (13) into (12), there is the following Equation (14):

dAsdzei(ωt+βmz)dAs+dzei(ωt+βmz)c.c.=i2ω2t2+[P¯Perturbation(r¯,t)]Ey(x)dx(14)

Where superscript − and + respectively represent propagation in the -z direction and the +z direction. The perturbation of the dielectric constant can be expressed as the periodic fluctuation of refractive index, which can be expressed as the following Equation (15):

P¯Perturbation(r¯,t)=Δn2(r¯)ε0E¯(r¯,t)(15)

Since Δn2(r¯)is a scalar, as can be seen from Equation (8), the periodic structure can only couple the TE mode to the TE mode or the TM mode to the TM mode. However, it cannot couple the TE mode.

For the TE mode propagation, substitute (10) into (15) and (14), there is the following Equation (16):

dAsdzei(ωt+βmz)dAs+dzei(ωt+βmz)c.c.(16)=iε04ω2t2m[Am+Δn2(x,z)]Ey(m)(x)Ey(s)(x)dxei(ωt+βmz)+c.c.]

Suppose the period of perturbation Δn2(x, z) is Δ, and lπΛβs,where l is the integer. So the Equation (17) can be obtained, which can be expressed as the following:

Δn2(x,z)Δn2(x)j=+ajei2jπΛz(17)

Substituting (17) into (16), there is:

dAsdz=iε04ωAs++Δn2(x)]Ey(s)(x)2dxei(2lπΛβs)(18)

Then the coupling equation by backward wave and forward wave at l-th harmonic can be expressed as the following Equation (19):

dAsdz=KcAs+ei2(Δβ)z(19)

Similarly, there is the following Equation (20)

dAs+dz=KcAsei2(Δβ)z(20)

In the Equation (20), Kc and Δβ can be expressed as the following forms:

Kc=iωε04at+Δn2(x)]Ey(s)(x)2dx(21)ΔββslπΛβsβ0(22)

Where Kc is called as the coupling coefficient. Equations (19) and (20) are the coupled mode equations for the forward and backward modes propagating in the periodic waveguide. Due to the agreement of two modes, so the power of these two modes is conserved.

For fiber Bragg grating, the waveguide structure is shown in Figure 4. Grating length is L, the grating period is Δ. The amplitude of light from the incident fiber, the amplitude of the backward wave at the boundary is 0. Assume that the fiber Bragg grating distribution is a strict sine function, the refractive index can be expressed as

Fiber Bragg grating waveguide structure diagram
Figure 4

Fiber Bragg grating waveguide structure diagram

n(z)=n¯+Δn(z)=n¯+δnsin(2πzΛ)(23)

The coupling coefficient of the grating is K=πδnηλBwhere η is the coefficient related to the mode energy in the remaining core. Here η = −V2 is approximated as the value of the fiber, V characterizes the modulus of the fiber transmission. Solving for Equations (19) and (20), the Equations (24) and (25) can be obtained:

As(z)eiβz=As+(0)iKeiβzΔβsinh(SL)+iScosh(SL)sinhS(zl)(24)As(z)eiβz=As+(0)eiβzΔβsinh(SL)+iScosh(SL)ΔβsinhS(zl)+iScoshS(zl)(25)

Where S=K2Δβ2so the fiber Bragg grating reflectivity is obtained from (24) and (25), which can be expressed as:

R=K2sinh2(SL)Δβ2sinh2(SL)+S2cosh2(SL)(K2>Δβ2)K2sin2(SL)Δβ2K2cos2(SL)(K2<Δβ2)(26)

When the incident light wavelength is equal to the center wavelength of the fiber Bragg grating [Bragg wavelength]. Δβ = 0, there is:

λB=2neffΛ(27)R(λB)=tanh2(KL)(28)

Where neff is called as the core effective refractive index.

Results and discussion

2.4 Nodes / edge mixed triangular units

The basic idea of finite element method is to separate the continuous solution area into a set of finite combinations, which are connected together in a certain way. Because the elements can be combined in different ways of joining and the elements themselves can have different shapes, the complex solution domain of geometric shapes can be modeled.

Figure 5(a) and (b) a typical triangular unit, expressed with e. The coordinates of the i-th vertex is (xi , yi)(i =

1, 2, 3). The longitudinal field component is represented by Ezie. The edge length of the edge is expressed as l. The magnitude of the tangential field in the direction of the edge. The area of the triangle is expressed as Δe, which can be expressed as:

Δe=121x1ey1e1x2ey2e1x3ey3e(29)

Define the area coordinates of the triangle as Lie=Δie/Δe where Δieis the value of the determinant obtained after the i-th row of the determinant. It has been replaced. Among them, the longitudinal field component is expressed as:

Eze(x,y)=i=13LieEzie=LeTEze=EzeTLe(30)

is the column vectors; superscript of T represents vector transpose. The horizontal field component can be expressed as:

Ete=i=13NtieEtie=NteTEte=EteTNte(31)

Where, N1te=(L1eL2eL2eL1e)l1e,N2te=(L2eL3eL3eL2e)l2e.

On the basis of Figure 5, it keeps the air filling rate constant. When the hole spacing, the change tendency of waveguide dispersion curve in the wave band of interest, just as shown in Figure 6. As can be seen from Figure 4, not only the slope of the curve has changed, but also the curve has shifted. When the ratio of a/Δ is constant and the inter-hole spacing is increased, the curve (chromatic dispersion zero point) moves toward the long wavelength while the slope becomes smaller. When the ratio of a/Δ is constant and the inter-hole spacing decreases, the curve shifts toward the short wavelength move. The accuracy of this method is much higher than that of conventional fiber optic biosensors and photonic crystal fiber biosensors [13, 14]. This is also the main reason for the application of superstructure photonic crystal fibers based on biosensor image models [15, 16].

(a) Hybrid edge/ node elements with triangular shape (b) The vector expansion function corresponding to any edge (e.g. edge 1) has tangential element only this edge
Figure 5

(a) Hybrid edge/ node elements with triangular shape (b) The vector expansion function corresponding to any edge (e.g. edge 1) has tangential element only this edge

The effect of the pitch on waveguide dispersion while keeping air filling coefflcient unchanged
Figure 6

The effect of the pitch on waveguide dispersion while keeping air filling coefflcient unchanged

Due to the periodic change the refractive index of fiber core area, waveguide conditions changed, coupling phenomenon also occurred in a certain wavelength range, fiber Bragg grating is a narrow–band filter transmission (reflection).

According to the above rules, you can first adjust the aperture so that the 1st and 2nd zero waveguide dispersion curve shifted to 1.55 μm wavelength window. Since the material dispersion curve is approximately linear on this window, if the ratio of is kept constant and then the hole spacing is adjusted, the slope of the waveguide dispersion will increase or decrease. As a result, it is necessary to find suitable structural parameters. The material dispersion and the waveguide dispersion slope are approximately equal and the signs are opposite, so that the total dispersion is approximately zero. The numerical simulation structure is shown in Figure 7. In this case, the dispersion is ± 0.2 ps/km/nm in the 1550 ± 100 nm band with a dispersion slope of −0.003 ps/km/nm2. The mode field 𝛬 = 3.24 μm, the nonlinear coefficient γ = 15.8 W−1 · km−1, which is much larger than that of conventional SMF. The zero-dispersion point is also located near the pump laser wavelength (1.55 μm). These data ensure that the parametric amplifier has a high gain and a large bandwidth.

The near zero flatten dispersion characteristics of PCF near 1.55 μm
Figure 7

The near zero flatten dispersion characteristics of PCF near 1.55 μm

It can be seen that the relative error calculation values of the highest point (peak), lowest point (peak) and time course of the original signal and the reconstructed signal are relatively small, and basically less than 10%, all less than 15%. This shows that the reconstructed signal is very close to the original signal, that is, the method used in this paper achieves better results by sampling and reconstructing the neural action potential signal with high signal-to-noise ratio after filtering.

The fabrication process of chemical component fiber Bragg grating generally includes the steps of hydrogen loading, hydrogen diffusion and high temperature annealing. The degradation and regeneration of grating occur at high temperature annealing stage. The refractive index modulation reduced to 20%. The high temperature test process of chemical component fiber grating is shown in Figure 8. This kind of grating has excellent high temperature stability, able to work under the high temperature of environment. Similar to previous experimental studies [17, 18, 19], the amount of change in the measurement results of this study, the sensitivity of the refractive index. Both the curve variation law and the feature distribution have reached a high degree of agreement, thus ensuring the accuracy and effectiveness of the research results. This lays a theoretical foundation for the application of superstructure photonic crystal fiber Bragg grating in biosensor image [20].

Chemical composition of optical fiber grating under refractive index modulation normalized decline curve
Figure 8

Chemical composition of optical fiber grating under refractive index modulation normalized decline curve

3 Conclusions

This study finds that the application of fiber grating sensor is a new subject and has a broad development prospect. That can be beneficial for large structure engineering, power engineering, geotechnical engineering, traffic engineering, mining engineering, blasting engineering, petrochemical industry, aerospace, biomedical, ships, oil exploration, and military weapons and equipment, etc. In this paper, the characteristics and basic principles of fiber grating are briefly introduced, and the application status of fiber grating sensing technology is analyzed. Then, according to the feature that the neural action potential can be sparse, the sparse matrix of neural signals can be constructed in the wavelet domain, and the signal can be decomposed by the wavelet filter. Furthermore, according to the principle of reflection characteristics of fiber grating, the wave equation of insulating medium is calculated from the coupled mode theory. Finally, under the condition of invariable in keeping the air filling factor, a comprehensive analysis of the influence of pitch of waveguide dispersion, it is concluded that the hole spacing increases, the curve is zero to the long wavelength direction, and the slope has shrunk, and hole spacing and small curve move in the direction of short wavelength. Contrast nerve action potential signal reconstruction based on sparse matrix values and original values, that the smaller values of the relative error between them, that was very close to reconstruct signal and original signal, the method adopted by the after filtering by a high signal-to-noise ratio of nerve action potential signal sampling refactoring got good results. The study will help researchers identify key areas that many have been unable to explore. Although the research on fiber grating sensing technology in China is relatively late, it can be predicted that with the deepening of the research on fiber grating sensing technology, there will be a new breakthrough in its application. Therefore, a number of major projects will be launched in succession, and fiber grating sensing technology will provide an effective way for long-term stable and real-time online security monitoring of these projects.

Acknowledgement

This work was supported by Natural Science Foundation of China (No. 61863013), Key R & D projects of Jiangxi science and Technology Department of China (No. 20161BBE50091) and Science and Technology Foundation of Jiangxi Educational Committee of China (No. 150529).

References

  • [1]

    Liu Y., Zhang J., Model study of the influence of ambient temperature and installation types on surface temperature measurement by using a fiber bragg grating sensor, Sensors, 2016, 16.7, 975 Google Scholar

  • [2]

    Sengupta S., Ghorai I., Swapan K., Biswas P., Design of superstructure fiber Bragg grating with efflcient mode coupling for simultaneous strain and temperature measurement with low cross-sensitivity, IEEE Sensors Journal, 2016, 16.22, 7941–7949. Google Scholar

  • [3]

    Chen T., Experimental study on cross-sensitivity of temperature and vibration of embedded fiber Bragg grating sensors, Opto–electr. Lett., 2018, 14.2, 92–97. Google Scholar

  • [4]

    Zhou L., Liang D., Zeng J., The influence of temperature on the deuteration mechanism and experimental research of fiber grating sensor performance, China Laser, 2012, 39.4, 1–6. Google Scholar

  • [5]

    Chong S., Dye concentration determination with cross-sensitivity compensation, Sensors and Actuators B: Chemical, 2016, 226, 450–456. Google Scholar

  • [6]

    Rajan G., Experimental study and analysis of a polymer fiber Bragg grating embedded in a composite material, J. Lightwave Techn., 2014, 32.9, 1726-1733. Google Scholar

  • [7]

    Zhang D., Zheng Y., Research on cross-sensitivity problem of FBG temperature sensor, Journal of Atmospheric and Environmental Optics, 2016, 11.3, 226-233. Google Scholar

  • [8]

    Carlo M., Michele A., Rosaria D., Fiber Bragg Grating Measuring System for Simultaneous Monitoring of Temperature and Humidity in Mechanical Ventilation, Sensors, 2017, 17.4, 749. Google Scholar

  • [9]

    Hu J., Yang Y., Liu X., Based on gap fiber grating micro-gap and temperature measurement technology, China Laser, 2014, 41.11, 201-206. Google Scholar

  • [10]

    Liu S., Chen T., Li R., Study on the Effect of Paste Effect of Substrate FBG Temperature Sensor on Its Performance, Optoelectr. Laser, 2016, 27.7, 692-698. Google Scholar

  • [11]

    Liu Y., Guo Z., Zhang Y., Simultaneous measurement of single fiber grating pressure and temperature, China Laser, 2010, 27. 11, 1002-1006. Google Scholar

  • [12]

    Huang J., Zhou Z., Tan Y., Design and Experimental Study of a Fiber Bragg Grating Pressure Sensor, International Conference on Innovative Design and Manufacturing, August 13-15, 2014, Montreal Canada, Piscataway: IEEE, 2014, 217-221. Google Scholar

  • [13]

    Zhu D., Qin H., Analysis of Sampled Fiber Bragg Grating Based on the Photonic Crystal Theory, Semicond. Photonics & Techn., 2009, 39.3, 184-188. Google Scholar

  • [14]

    Geng D., Yang D., Zhang X., All-optical generation of microwave using a photonic crystal fiber Brillouin laser based on Bragg grating Fabry-Perot cavity, Microwave & Optical Techn. Lett., 2010, 50.3, 809-814. Google Scholar

  • [15]

    Liu H., Tan C., Zhu C., Multi-parameters measurement based on cascaded Bragg gratings in magnetic fluid-infiltrated photonic crystal fibre, Optica Acta Int. J. Optics, 2016, 64.8, 887-894. Google Scholar

  • [16]

    Rindorf L., Bang O., Sensitivity of photonic crystal fiber grating sensors: biosensing, refractive index, strain, and temperature sensing, J. Optic. Soc. Amer. B, 2007, 25.3, 310-324. Google Scholar

  • [17]

    Tan X., Geng Y., Li X., High-birefringence photonic crystal fiber Michelson interferometer with cascaded fiber Bragg grating for pressure and temperature discrimination, Optic. Engin., 2016, 55.9, 090508. Google Scholar

  • [18]

    Wang B., Huang Y., Research on characteristics of the twelve polygon photonic crystal fiber grating, Semicond. Optoelectr., 2014, 35.6, 1000-1003. Google Scholar

  • [19]

    Yiping W., Wang M., Huang X., In fiber Bragg grating twist sensor based on analysis of polarization dependent loss, Optics Express, 2013, 21.10, 11913-11920. Google Scholar

  • [20]

    Carpintero G., Balakier K., Yang Z., Microwave Photonic Integrated Circuits for Millimeter-Wave Wireless Communications, Lightwave Technology Journal, 2014, 32.20, 3495-3501. Google Scholar

About the article

Received: 2018-08-26

Accepted: 2018-09-23

Published Online: 2018-11-30


Citation Information: Open Physics, Volume 16, Issue 1, Pages 727–733, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0091.

Export Citation

© 2018 B. Li and J. Kanugo, published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0

Comments (0)

Please log in or register to comment.
Log in