Abstract
The space-time properties of the expressions of sub-cycle and single-cycle focused optical pulses with radially polarized electromagnetic field based on the Sink-Source model are studied. The self-induced blue shift of the center frequency of spectrum in the center of the pulse field is found to have an important impact on the electrons acceleration. When the electrons approach to the center of pulse, the electrons will obtain a large kinetic energy gain in a short time. The effect of radiation-reaction force can’t be ignored if the net kinetic energy gain of electrons is more than GeVs. The electrons will deviate from the original acceleration channel and the gain of kinetic energy that electrons may gain will be greatly reduced if the radiation-reaction effect is considered. In contrast to the few-cycle laser pulse accelerating electrons, the gain of kinetic energy obtained by electrons is a few times higher and the corresponding peak optical power is one order of magnitude lower in the case of the sub-cycle laser pulses accelerating electrons. The maximal kinetic energy gain of electrons is robust against the variation of the incident angles.
© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Electron acceleration by ultrashort optical pulse is an interesting and important topic of research. Many works have been done on this topic [1–8]. The electron acceleration by few-cycle and sub-cycle laser pulses are also investigated [9–15]. In the most cases, the radial polarized optical beams were focused to obtain the optical field with strong axial electric field. Two methods are used to describe the radial polarized optical beams. One method is based on Lax series where the electric fields are expanded as series of diffraction angle , where is the beam waist and is the Rayleigh range [11]. The other method is based on plane wave angular spectrum representation [16–18]. But for the sub-cycle pulses, the temporal and spatial parameters are coupled, which causes the special properties of sub-cycle pulse [19]. The expressions of sub-cycle focused optical pulse with radially polarized electromagnetic (EM) field should be deduced directly from Maxwell’s equation. The expressions of sub-cycle pulsed focused vector beams can be derived from Sink-Source Model [14,19–21]. Many new properties of sub-cycle pulse, such as the self-induced blue-shift of the central frequency have been found and proved by the experiment [22]. In our previous works, the expressions of sub-cycle pulsed focused beams with the transverse polarization of electric field on the beam waist plane are obtained. The electron can be accelerated to sub-GeV energies by the sub-cycle pulses when the peak intensity is [10]. The expression of Poisson-like ultrashort pulse based on Sink-Source Model is given in [14]. The spectrum of the pulse can be reduced to a Gaussian function if its duration is sufficiently long. In this paper, both of the expressions of Lorentz type and Poisson type sub-cycle and single-cycle focused optical pulses with radially polarized EM field are derived from the Sink-Source Model. The acceleration of electrons by the two types of pulses is studied after accounting for the radiation-reaction force.
The process of radiation is fundamentally quantum. Such both of the electrons and the pulses should be processed by the quantum methods [14]. In this paper, the cases of sub-cycle laser pulses are considered. The space and time scales of pulse are extremely limited. Such classical methods can be good approximations. The electronic trajectories can also be conveniently studied.
The effect of the spin force on the acceleration of electrons is also researched by some papers [23,24]. For the sub-cycle pulse cases, the effects of the spin force are a wide range of issues and need to be discussed carefully in the future. In this paper, the topics are focused on the effect of the time-space structure of sub-cycle pulse on the classical radiation reaction force and the radiation reaction force on the motion of atom.
2. Model and equations
The expressions of sub-cycle and single-cycle focused pulses propagating along a certain direction can be deduced from dipole model using the complex sink-source source theory [14]. The complex sink-source source theory is traditionally considered to overcome the singularity problem on the beam waist plane of the expressions of the complex point source theory. But recent research indicates that the discontinuities still arise in phasors generated from the complex source-sink model for all odd orbital angular momentum modes [25]. In order to study the effect of different pulse configurations on the electron acceleration, the Lorentz type and Poisson type pulses are studied. The temporal shape of Poisson pulse vanishes more quickly at large time than that of the Lorentz pulse. The space-time coupling property of ultrashort pulse beams are discussed in the Ref [26]. and the Lorentz sub-cycle pulse is discussed in the Ref [27].
The expressions of sub-cycle focused pulse with radially polarized EM field can be deduced by using the Sink-Source method and the dipole moment model [14]. The dipole moment can be described as an oscillating dipole source located at the origin of the coordinate. The upper arrow indicates the vectors and the upper dot indicates derivative.
can be written as for the Lorentz type pulse and for the Poisson like pulse. is the pulse width parameter, is the carrier frequency and is a real positive parameter.The focused pulse propagating in the z direction can be obtained by introducing the spatiotemporal translation
where. is the wave vector. is the waist spot size, and is a constant value in the simulation. The carrier frequency , together with the waist spot size, determines the parameter. The complex distance is introduced as
According to the Sink-Source method, the optical field can be obtained by folding the source with the function,The complex retarded times are introduced as and . Such time domain pulse functions can be written as and.The detailed derivation process of the expressions of electromagnetic field can be found in the Ref [14]. and they can be written asTo remove from the expression of electromagnetic field, the peak value of the electric field is assumed to be . The peak value of the electric field appears at the position and in the time can be written asSuch expressions of can be deduced asThen can be replaced by .The electron acceleration by the Lorentzian focused sub-cycle pulse can be studied by solving the relativistic Newton-Lorentz equation. After the radiation-reaction force of electrons is considered, the modified relativistic Newton-Lorentz equation is given as:
where is the velocity of particle scaled by . The classical radiation reaction model is discussed in Ref [28]. In this paper, the Landau-Lifshitz equation is selected because it is concise and clear, and is convenient for numerical calculations [29]. The component equations of Eq. (17) and the displacement equations can be written asIn the SI system, the general form of the radiation reaction term can be written as:Equation (8) can be solved by some high precision numerical algorithms. The initial speed of the electron is and the initial kinetic energy of the electron is .The radiation losses in the process of electron acceleration can be calculated from the Lienard-Wiechert potential. The radiation power of an electron in the acceleration process can be written as
the radiation energy loss can be written as . The velocity and the acceleration of an electron in the spherical coordinate system can be written aswhere . The expressions of are the same forms.In the numerical computations, we use the fifth-order prediction-correction Runge-Kutta algorithm. In order to avoid too long calculation time, the algorithm must be adaptive variable step. We use the matlab function ode45 to do the numerical computations. In order to ensure the accuracy of the result, the relative error parameter ‘RelTol’ is set toand the absolute error parameter ‘AbsTol’ is set to . The numerical results can be confirmed by our self-written fifth-order Runge-Kutta Cash-Karp algorithm. The computing speed of ode45 is faster.
3. Space-time characteristics of sub-cycle and single-cycle pulses
The space-time characteristics in the center of the non-paraxial pulse field can be analyzed by setting . The temporal waveform, power and the frequency spectrum of a 0.45-cycle Poisson-like pulse are shown in Fig. 1. The frequency is and the period is . The parameter is 4.12.
From Fig. 1(a), we can define the number of cycles of the pulse by the ratio of the full width at half maximum(FWHM) of the power to the period . Figure 1(b) shows the self-induced blueshift of the center frequency of spectrum in the center of the pulse field. The center frequency of spectrum is in Fig. 1(b), which is much higher than. The reason is that the sub-cycle pulse has an ultra-wide spectrum. For the tightly focused sub-cycle radially polarized EM field, there are strong axial electric field distributions on the waist plane. According to the principle of the Rayleigh diffraction limit, the focus spot of the high frequency component is much smaller than that of the low frequency component on the waist plane. Figure 2 is given to clarify this question.
From Fig. 2, one can see that the focus spots of the short wavelength components are far smaller than that of the long wavelength components. Therefore, the high-frequency components are more concentrated in the center of the pulse than that of the low frequency component. That is called as the central frequency blue shift of pulse. This has a significant impact on the acceleration of electrons due to the large electric field gradient in the center of the pulse.
4. Results and discussions
The initial electron velocity components in thedirection are and . The initial velocities can be easily prepared by a simple electron accelerator. The choice of the initial velocities has no effect on the main conclusions in the paper. These cases that the electrons initially enter the light beam at a wide range angle are considered. If the time of the pulse peak is zero, the initial interaction time is negative for the 0.45-cycle pulse and negative for the 1.37-cycle pulse. Zero time is the time when the pulse center arrives. The initial positions of the electron are, where is the initial velocity component in thedirection and is the initial incidence angle. The initial electron velocity is .The frequencyis . The peak intensity of pulse is . For reading, the different input parameters of most figures are summarized in the Table 1.
In Table 1, c indicates the velocity of light and fs indicates femtosecond.
4.1 The 0.45-cycle pulse case
For the Poisson type 0.45-cycle pulse, the exit kinetic energy of the electron versus the carrier-envelope phase (CEP) and the beam waist is shown in Fig. 3. The initial positions of an electron are and the velocity is . The incidence angle is 10 degrees in Fig. 3.
Figure 3 shows that there is a big difference for the kinetic energy of electrons when the radiation-reaction effect is considered in Fig. 3(a) and is neglected in Fig. 3(b). The ideal working point in the waist-CEP plane is modified when the radiation reaction force is considered. In Fig. 3(a), the maximal exit kinetic energy gain of an electron is when the CEP is and the beam waist is . In Fig. 3(b), the maximal exit kinetic energy gain of an electron is when the CEP is and the beam waist is . This can be explained from Fig. 3(c)-3(d). In order to clearly demonstrate the effects of radiation reaction, only the case of the initial acceleration phase is shown in Fig. 3(d). The significant differences of kinetic energy of electrons occur from the moment of 600 femtoseconds. When the electrons enter the center of the laser pulse, the strong electric field gradients cause significant acceleration of electrons, the radiation reaction force also increases significantly at the same time. The cumulative effect of the radiation reaction force causes that the electron can’t always stay in the acceleration phase [30]. From the embedded graph of Fig. 3(c) and Fig. 3(d), one can see that the electron experiences more acceleration and deceleration processes when the radiation reaction force is considered than that the radiation reaction force is neglected. The self-induced blue shift of the frequency in the center of the sub-cycle pulse field significantly increases the effect of the radiation reaction force on the electron’s motion. From Fig. 3(d), one can see that the electron leaves the beam with a larger ejection angle when the radiation reaction force is considered than that is neglected. This indicates that the electron leaves the beam in a relatively short period of time. The electron can obtain more than exit kinetic energy gain in the ranges of CEP varying fromtoand beam waist varying from to when the radiation reaction force is considered. This indicates that the electron can obtain the high kinetic energy gain in the wide range of the beam waist and CEP parameters. The parameter ranges in which the electron obtains high exit kinetic energy are far smaller in the CEP interval of tothan that of 0 to. This indicates that the initial direction of the electromagnetic field plays an important role in the electron acceleration. The radiation loss of the acceleration electron in the case of Fig. 3(c) is .
The exit kinetic energy of the electron versus the incidence angle for the Poisson type 0.45-cycle pulse is shown in Fig. 4. The exit kinetic energy of the electron for the case of the Lorentz type 0.45-cycle pulse is also shown in Fig. 4. The peak intensities of the two types of pulses are.
Figure 4(a) shows the exit kinetic energy of the electron versus the beam waist and CEP for the initial incident angle of 90 degrees. Comparing Fig. 3(a) and Fig. 4(a), one can see that there are only some minor changes for the ideal working point when the incident angle changes. From Fig. 4(b), we can see that there are only some little changes for the maximal kinetic energy of an electron when the incident angle varies from 10 degrees to 90 degrees. This shows that the maximal kinetic energy gain of electrons is robust against the variation of the incident angles. The reason is that the Poisson type sub-cycle pulses have large time domain electric field gradients and the fringe electric fields disappear quickly. Such an electron can enter the center of the pulse field and interacts with the peak field. In Fig. 4(c), the maximal exit kinetic energy gain of an electron is when the CEP is and the beam waist is . The parameter ranges to win a high kinetic gain are also much smaller in Fig. 4(c) than that in Fig. 4(a). Figure 4(d) also shows that the exit kinetic energy of electrons will decrease as the incidence angle increases. The reason is that the electric field amplitude of the Lorentz sub-cycle pulse decreases slowly with time. Such electrons are easily driven away from the beam by the fringe electric field of the pulse. The exit kinetic energy of electrons is also sensitive to the Lorentz sub-cycle pulse parameters.
Figure 4 also shows that the electron can obtain more than 9 GeV energy gain in a very wide range of incident angles by the Poisson sub-cycle pulse when the peak intensity is , which is far higher than that of the Lorentz sub-cycle pulse in the same peak intensity of light. Compared with the case of few-cycle laser pulse accelerating electrons, the gain of kinetic energy of electrons is a few times higher and the corresponding peak optical power is one order of magnitude lower in the case of the sub-cycle laser pulses accelerating electrons [13]. Such a Poisson type sub-cycle pulse is more suitable for accelerating electrons than a Lorentz type sub-cycle pulse.
The case that is shown in Fig. 5. The temporal configuration of the pulse is Poisson type.
The maximal exit kinetic energy gain of electrons is in Fig. 5(a). Figure 5 shows that the exit kinetic energy of the electron is robust against the variation of the incident angles. From Fig. 5(a), one can see that the electron can obtain more than exit kinetic energy gain in the wide ranges of CEP varying fromtoand beam waist varying from to. Figure 5(c) shows that the maximal kinetic energy of an electron remains approximately constant when the incident angle varies from 10 degrees to 90 degrees.
4.2 The 1.37-cycle pulse case
The exit kinetic energy of the electron versus the CEP and the beam waist is shown in Fig. 6 for the Poisson type 1.37-cycle pulse. is and the parameter is 28.4.
The maximal exit kinetic energy gain of electrons is in Fig. 6(a), which is lower than the maximal kinetic energy in Fig. 5(a). This shows that the sub-cycle pulses are more suitable for accelerating electrons than the single cycle pulses. The reason is that the electrons are easier to be driven away from the beam by the fringe electric field of the pulse if the pulse duration is shorter. The electron can obtain more than exit kinetic energy gain in the wide ranges of CEP varying from to and beam waist varying from to .The maximal kinetic energy gain of electrons is also robust against the variation of the incident angles. The radiation loss of the acceleration electron in the case of Fig. 6(b) is .
In Fig. 7, is .The temporal configuration of the pulse is Poisson type.
Compared with Fig. 6, the initial velocity of electrons is higher in Fig. 7. We can see that the parameter ranges for obtaining the high exit kinetic energy gain is wider in Fig. 7(b) than that in Fig. 7(a). If there is a large angle between the direction of the velocity of electrons and the propagation velocity of the beam, electrons easily enter the center of the pulse and are not easily driven away by the fringe electric field of the pulse. The maximal exit kinetic energy gain of electrons is in Fig. 7(b). Figure 7(c) shows that the maximal kinetic energy gain of electrons is also robust against the variation of the incident angles. Figure 5 and Fig. 7 show that the electron can obtain more net exit kinetic energy gain with shorter FWHM of pulses.
4.3 The different initial positions and misaligned cases
In this section, these cases that the electrons are initially located at different locations are studied. Two cases are considered. The first case is that the initial interaction time is changed to negative for the 0.45-cycle pulse and negative for the1.37-cycle pulse. The second case is that the initial positions of electrons are changed to for both the 0.45-cycle and 1.37-cycle pulses. The incidence angle is 10 degrees. The first case is shown in Fig. 8.
The maximal exit kinetic energy gain of electrons is in Fig. 9(a) and in Fig. 9(b). Figure 9 shows that the electron can obtain the high exit kinetic energy in the wide CEP and beam waist parameter ranges. This indicates that with some specific CEP and beam waist parameters, the electron can enter the capture acceleration channel when interacting with the fringe electric fields. Such an electron can meet the center of the pulse field and obtain a high exit kinetic energy gain.
The low velocity electrons can be easily prepared by a simple, commonly used electron accelerator. In the experiment, the most concern is the alignment of electrons and laser pulses. The alignment means that the electrons and the center of pulse will arrive at the same point at the same time if the interaction of electrons with pulses isn’t considered. In this paper, the analysis points out that this is not necessarily required. The arrival time can be staggered by dozens of femtoseconds. If the electron is captured by the edge field of pulse, the electron has the opportunity to interact with the center field of pulse and obtain a high exit kinetic energy.
Sub-cycle laser pulses can already be produced in the laboratory by the synthetic methods [31]. With the development of technology, sub-cycle laser pulses will be more convenient to produced and used for the electron acceleration experiments.
5. Conclusions
Using the Poisson-like and Lorentz type sub-cycle and single cycle pulse models, the electrons acceleration problem is studied. The space-time characteristics of sub-cycle and single-cycle pulses are found to have an important impact on the electrons acceleration. Compared with the Lorentz type pulse, the Poisson-like pulses can accelerate electrons to the high exit kinetic energy states in the wide CEP and beam waist parameter ranges. The maximal kinetic energy gain of electrons is robust against the variation of the incident angles for the Poisson-like pulse cases. The electron can obtain more net exit kinetic energy gain with the shorter FWHM of pulses. Compared with the few-cycle laser pulse accelerating electrons case, the gain of kinetic energy obtained by electrons is a few times higher and the corresponding peak optical power is one order of magnitude lower in the sub-cycle laser pulses accelerating electrons case. The radiation-reaction effect is also important in the sub-cycle and single cycle laser pulses accelerating electrons case.
Funding
National Natural Science Foundation of China (61727821, 61475139, 11564005, 11764008); Key Project for Innovation Research Groups of Guizhou Provincial Department of Education, China (KY[2016]030 and 028, KY[2017]035). The Thousand Levels of Innovative Talents of Guizhou Province ([2016]016).
References
1. C. N. Harvey, A. Gonoskov, A. Ilderton, and M. Marklund, “Quantum Quenching of Radiation Losses in Short Laser Pulses,” Phys. Rev. Lett. 118(10), 105004 (2017). [CrossRef] [PubMed]
2. Z. Zhang, L. Yan, Y. Du, W. Huang, C. Tang, and Z. Huang, “Generation of high-power, tunable terahertz radiation from laser interaction with a relativistic electron beam,” Phys. Rev. Accel. Beams 20(5), 050701 (2017). [CrossRef]
3. S. Carbajo, E. A. Nanni, L. J. Wong, G. Moriena, P. D. Keathley, G. Laurent, and R. J. Dwayne Miller, “Direct longitudinal laser acceleration of electrons in free space,” Phys. Rev. Accel. Beams 19(2), 021303 (2016). [CrossRef]
4. M. Dunning, E. Hemsing, C. Hast, T. O. Raubenheimer, S. Weathersby, D. Xiang, and F. Fu, “Demonstration of Cascaded Optical Inverse Free-Electron Laser Accelerator,” Phys. Rev. Lett. 110(24), 244801 (2013). [CrossRef] [PubMed]
5. Y. I. Salamin and F. H. M. Faisal, “Ultrahigh electron acceleration and Compton emission spectra in a superintense laser pulse and a uniform axial magnetic field,” Phys. Rev. A 61(4), 043801 (2000). [CrossRef]
6. E. N. Nerush and I. Yu. Kostyukov, “Carrier-Envelope Phase Effects in Plasma-Based Electron Acceleration with Few-Cycle Laser Pulses,” Phys. Rev. Lett. 103(3), 035001 (2009). [CrossRef] [PubMed]
7. G. I. Dudnikova, V. Yu. Bychenkov, A. Maksimchuk, G. Mourou, J. Nees, S. G. Bochkarev, and V. A. Vshivkov, “Electron acceleration by few-cycle laser pulses with single-wavelength spot size,” Phys. Rev. E. 67(2), 026416 (2003). [CrossRef] [PubMed]
8. K. Schmid, L. Veisz, F. Tavella, S. Benavides, R. Tautz, D. Herrmann, A. Buck, B. Hidding, A. Marcinkevicius, U. Schramm, M. Geissler, J. Meyer-Ter-Vehn, D. Habs, and F. Krausz, “Few-Cycle Laser-Driven Electron Acceleration,” Phys. Rev. Lett. 102(12), 124801 (2009). [CrossRef] [PubMed]
9. F. Terranova, “Particle acceleration in subcycle optical cells,” Phys. Rev. Spec. Top. Accel. Beams 17(7), 071301 (2014). [CrossRef]
10. X. M. Cai, J. Y. Zhao, Q. Lin, and J. L. Luo, “Electron acceleration by subcycle pulsed focused vector beams,” J. Opt. Soc. Am. B 33(2), 158–164 (2016). [CrossRef]
11. B. Rau, T. Tajima, and H. Hojo, “Coherent Electron Acceleration by Subcycle Laser Pulses,” Phys. Rev. Lett. 78(17), 3310–3313 (1997). [CrossRef]
12. V. Marceau, A. April, and M. Piché, “Electron acceleration driven by ultrashort and nonparaxial radially polarized laser pulses,” Opt. Lett. 37(13), 2442–2444 (2012). [CrossRef] [PubMed]
13. C. Varin, M. Piché, and M. A. Porras, “Acceleration of electrons from rest to GeV energies by ultrashort transverse magnetic laser pulses in free space,” Phys. Rev. E. 71(2), 026603 (2005). [CrossRef] [PubMed]
14. A. April, Coherence and Ultrashort Pulse Laser Emission (InTech, pp.355–382, 2010).
15. S. Payeur, S. Fourmaux, B. E. Schmidt, J. P. MacLean, C. Tchervenkov, F. Légaré, M. Piche, and J. C. Kieffer, “Generation of a beam of fast electrons by tightly focusing a radially polarized ultrashort laser pulse,” Appl. Phys. Lett. 101(4), 041105 (2012). [CrossRef]
16. J. L. Liu, Z. M. Sheng, and J. Zheng, “Exact solution of the vectorial field structure of a light beam,” Opt. Commun. 284(19), 4646–4649 (2011). [CrossRef]
17. R. Martínez-Herrero, P. M. Mejías, S. Bosch, and A. Carnicer, “Vectorial structure of nonparaxial electromagnetic beams,” J. Opt. Soc. Am. A 18(7), 1678–1680 (2001). [CrossRef] [PubMed]
18. R. Martinez-Herrero, P. M. Mejias, and A. Manjavacas, “On the longitudinal polarization of non-paraxial electromagnetic fields,” Appl. Phys. B 99(3), 579–584 (2010). [CrossRef]
19. Q. Lin, J. Zheng, and W. Becker, “Subcycle Pulsed Focused Vector Beams,” Phys. Rev. Lett. 97(25), 253902 (2006). [CrossRef] [PubMed]
20. C. J. R. Sheppard and S. Saghafi, “Beam modes beyond the paraxial approximation: A scalar treatment,” Phys. Rev. A 57(4), 2971–2979 (1998). [CrossRef]
21. E. Heyman, B. Z. Steinberg, and L. B. Felsen, “Spectral analysis of focus wave modes,” J. Opt. Soc. Am. A 4(11), 2081–2091 (1987). [CrossRef]
22. Q. Lin, J. Zheng, J. Dai, I. Chen Ho, and X. C. Zhang, “Intrinsic chirp of single-cycle pulses,” Phys. Rev. A 81(4), 043821 (2010). [CrossRef]
23. S. M. Mahajan, F. A. Asenjo, and R. D. Hazeltine, “Comparison of the electron-spin force and radiation reaction force,” Mon. Not. R. Astron. Soc. 446(4), 4112–4115 (2015). [CrossRef]
24. J. P. Costella and B. H. J. McKellar, “Electromagnetic Deflection of Spinning Particles,” Int. J. Mod. Phys. A 9(3), 461–473 (1994). [CrossRef]
25. A. Vikartofsky, L. W. Pi, and A. F. Starace, “Discontinuities in the electromagnetic fields of vortex beams in the complex source-sink model,” Phys. Rev. A. 95(5), 053826 (2017). [CrossRef]
26. C. F. R. Caron and R. M. Potvliege, “Free-space propagation of ultrashort pulses: Space-time couplings in Gaussian pulse beams,” J. Mod. Opt. 46(13), 1881–1891 (1999). [CrossRef]
27. J. Zheng, E. M. Qiu, and Q. Lin, “High harmonic generation with sub-cycle pulses,” J. Opt. 13(7), 075206 (2011). [CrossRef]
28. H. Spohn, “The critical manifold of the Lorentz-Dirac equation,” Europhys. Lett. 50(3), 287–292 (2000). [CrossRef]
29. L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, 4th ed. (Pergamon Oxford, 1989).
30. P. X. Wang, Y. K. Ho, X. Q. Yuan, Q. Kong, N. Cao, A. M. Sessler, E. Esarey, and Y. Nishida, “Vacuum electron acceleration by an intense laser,” Appl. Phys. Lett. 78(15), 2253–2255 (2001). [CrossRef]
31. A. A. Voronin, J. M. Mikhailova, M. Gorjan, Z. Major, and A. M. Zheltikov, “Pulse compression to subcycle field waveforms with split-dispersion cascaded hollow fibers,” Opt. Lett. 38(21), 4354–4357 (2013). [CrossRef] [PubMed]