Abstract
The possibility of second-harmonic generation based on surface dipole and bulk multipole nonlinearities in silica nanowires is investigated numerically. Both circular and microstructured nanowires are considered. Phase matching is provided by propagating the pump field in the fundamental mode, while generating the second harmonic in one of the modes of the multiplet. This is shown to work in both circular and microstructured nanowires, although only one of the modes can be phase-matched in the microstructure. The prospect of obtaining large conversion efficiencies in silica-based nanowires is critically discussed, based on simulations of second-harmonic generation in nanowires with a fluctuating phase-matching wavelength. It is concluded that efficient wavelength conversion will require strong improvements in the nanowire uniformity, peak powers well in excess of , increase of the second-order nonlinearity by an order of magnitude by use of a different base material, or highly polarizable surface coatings.
© 2010 Optical Society of America
1. INTRODUCTION
Second-harmonic generation (SHG) in doped-silica fibers was discovered in the early 1980s and intensively studied in the following years. In 1986 conversion efficiencies as high as 3–5% were reported [1]. A subsequent theoretical study [2] showed that such high efficiencies were orders of magnitude above what could reasonably be explained by interface or bulk multipole contributions to the second-order susceptibility, . Since the isotropic nature of amorphous silica forbids the presence of a (bulk) dipole contribution to , the origin of the experimentally observed phenomena was initially a puzzle, but was eventually explained by the formation of a bulk dipole through multiphoton processes involving both pump and SHG light [3]. In this model, the multiphoton inscription will naturally lead to the formation of a grating that compensates the phase mismatch arising from waveguide and material dispersion in the fiber. It also explains why the efficient SHG reported in [1] was only found after exposing the fiber to pump radiation for a certain incubation period, which was subsequently found to be shortened by the presence of seed SHG radiation [4]. However, in spite of much effort it has proved impossible to scale up the SHG efficiency beyond the level of a few percent. This can be explained as a self-saturation effect caused by the interference of the SHG light itself with the grating formed in the fiber [3].
Recent experiments on second- and third-harmonic generation in pure-silica nanowires found that a noticeable SHG signal could be observed in a nanowire with an estimated interaction length of only [5]. The authors attributed this finding to a surface effect. Indeed, the advent of nanowire technology [6] calls for a reexamination of the surface and bulk multipole effects in fibers. This is because the nanowires offer a much larger index contrast than the values assumed in [2], allowing SHG phase matching between the fundamental and second-order modes, and greatly increasing the achievable intensities at the surface. Furthermore, the microstructured optical fibers developed during the last decade offer a possibility of manufacturing nanoscale fiber cores embedded in a fiber with a sufficient outer diameter to allow for practical handling outside of the laboratory. However, the complex geometries of microstructured fibers imply that the analytical expressions valid for circular step-index fibers cannot be used, and must be replaced by numerical simulations.
The purpose of this work is to investigate the prospects for surface and bulk multipole SHG by phase matching between fundamental and second-order modes in circular as well as structured nanowires. While circular nanowires are studied by analytical methods, the structured fibers will be studied numerically using the finite-element method. It will be shown that significant power conversion can in principle be achieved in nanowire lengths of a few cm, but that structural fluctuations along the nanowire are likely to be a severe limitation, at least for silica-based nanowires.
2. THEORY AND NUMERICAL METHODS
In the small-signal limit where pump depletion is negligible the SHG process can be described by the equation
where , are the field amplitudes of the fundamental and second-harmonic signals, respectively; is the phase mismatch between the guided waves at the fundamental and second-harmonic frequencies, and is the overlap integralwhere is the second-harmonic frequency. The fields are normalized to unit power, i.e., the pump and signal powers, and , are given by and . The electric and magnetic fields of the guided modes are expressed as for . In the following, the shorthand notation will be used. is the second-order nonlinear polarization, which will have both dipole contributions from the fiber surface and multipole contributions from the bulk. The bulk contributions can be written asA third term, proportional to , vanishes (or can be included in the surface term) for fibers with piecewise constant material profiles, such as the ones considered here. For the surface contributions, it is convenient to express the results in terms of the field and polarization components parallel and perpendicular to the fiber surface. Three distinct terms then contribute to : where is the unit vector normal to the surface. The surface components appearing in Eqs. (7, 8, 9), as well as the bulk multipole parameters γ and δ have recently been measured by Rodriguez et al. [7]. Two sets of values, differing by about a factor of 1.5 were reported, resulting from the use of two different calibration methods. In the experiments, the contribution from the bulk multipole term proportional to γ in Eq. (5) cannot be separated from the surface contributions in Eqs. (7, 8). However, various indirect arguments support the notion that the bulk susceptibility is predominantly of magnetic dipole character [8], which implies . Under this assumption, and taking the lowest estimates put forward in [7], one arrives at the components given in Table 1 .In the case of a circular silica nanowire, analytic solutions for the guided-mode electromagnetic fields are available. Throughout this work, the input (pump) field will be assumed to propagate in the fundamental mode of the fiber, while SHG (signal) field generation in either the or the modes will be considered. Working in cylindrical coordinates, the fields of the modes inside the silica nanowire can be expressed as [9]
where are Bessel functions of the first kind, are modified Bessel functions of the second kind, is the vacuum wave vector of the light, is the silica refractive index, is the modal effective index, a is the nanowire radius, and φ is a phase that determines the polarization. Two orthogonal polarization states can be obtained by, e.g., setting . The harmonic dependence on z and t, , has been omitted in these formulas. The propagation constant, β can be determined from the implicit equation [9]From Eqs. (10, 11, 12, 13, 14, 15, 5, 6, 7, 8, 9) it may be shown that pump light in the mode will only generate SHG signal in the mode with the same value of φ. Alternatively, it may generate signal in the mode, whose field distribution can be expressed as [9] and whose propagation constant is determined by the equation [9] This mode is often referred to as the “radially polarized” mode because of the strictly radial nature of the transverse electric-field component. It should be noted that because of the large index contrast and small dimensions of silica nanowires, the z components of the fields may be significant, and full-vectorial expressions for the overlap integrals are essential. In the limit of small index contrasts, the and modes are both part of the multiplet, which also comprises the mode, whose only nonvanishing electric-field component is . The overlap between this mode and arising from pump light can be shown to vanish. The nonvanishing surface contributions to then become: Here and are the amplitudes of the fields when normalized to unit power. The expressions for the bulk polarization terms are somewhat more involved and will be omitted here. A discussion of these terms can be found in [2].Whereas the case of a nanowire with a circular transverse profile can be treated analytically, numerical solutions are required for the microstructured nanowire. In the present work, guided-mode fields and propagation constants were determined using the finite-element modeling tool COMSOL3.5. The generic fiber structure is depicted in Fig. 1 . The three parameters characterizing the structure are the outer radius of the microstructure R, the bridge width , and the curvature radius of the curved core surfaces . In the calculations, R was set to , and absorbing perfectly-matched-layer boundary conditions were applied beyond that radius, to emulate a fiber with a large ring of massive silica surrounding the microstructure. The parameters and were taken as variables. The material dispersion of pure silica was described by a three-term Sellmeier polynomial, with the coefficients given by Okamoto [9].
In the limit of a small SHG signal, and assuming perfect phase matching, , the evolution of the SHG intensity is given by:
Even in the case of perfect phase matching, Eq. (26) will become inaccurate when the SHG power becomes large enough that pump depletion must be taken into account. A complete description of SHG dynamics will not be attempted here, since the undepleted-pump approximation is sufficient to reveal whether significant (e.g. ) SHG conversion efficiencies are realistic. As is evident from Eq. (26), this requires where is the pump power and L the length of the nanowire.Away from perfect phase matching, but still assuming an undepleted pump, Eq. (1) can be integrated to
For a spectrally narrow pump, can be found by a first-order expansion around the perfect phase matching wavelength aswhere , are the group indices of the fundamental and second-harmonic modes, respectively. From this expression, the spectral FWHM (at the SHG wavelength) can be evaluated toThe bandwidth thus scales inversely with the interaction length L.3. NUMERICAL RESULTS
3A. Circular Nanowire
Surface and bulk nonlinear coefficients were obtained for pump radiation in the fundamental mode and SHG radiation in the and modes. The phase-matched SHG wavelengths, , for and modes are plotted against the nanowire diameter in Fig. 2 . is seen to be roughly similar to the fiber diameter, which implies that phase matching happens at a V-parameter value of for the fundamental wavelength, where is the usual V-parameter as defined in the theory of step-index fibers. The difference in between and modes is about 10%.
In Fig. 3 , the sum of all SHG nonlinear coefficients, is plotted together with the individual terms as calculated from Eqs. (20, 21, 22, 23, 24, 25). For the mode, the term is clearly dominant, whereas for the mode both the and give significant contributions to the total sum. The total SHG coefficient roughly scales with . This is readily understood from Eq. (2) and Fig. 2. The fiber diameter scales with λ, which implies that the electric field strength for a fixed total power must scale with . The scaling of the prefactor in Eq. (2) is then cancelled by the λ scaling of the surface area over which the fields are integrated.
3B. Microstructured Nanowire
Replacing the perfectly circular nanowire with the microstructure depicted in Fig. 1 has significant implications for the phase matching conditions. The general modal structure known from the circular nanowire is preserved, and one can find second-order modes with field distributions resembling both the and modes. The symmetry arguments that imply that only two second-order modes are coupled to a given fundamental mode by the surface and bulk SHG terms also hold up because of the mirror symmetry of the fiber in a plane running parallel to any of the three core-supporting struts through its center. However, if the struts are regarded as slab waveguides coupled to the core region of the microstructure, the transverse field of the -like core mode matches well onto the TE modes in any of these slab waveguides. The mode, on the other hand, can only match well onto the TE mode in one of the struts, and the mode always matches onto the TM modes of the struts. As a result, the -like core mode effectively sees a higher cladding index in the microstructure than the and modes, and therefore ends up having the highest effective index. Furthermore, the effective cladding index increases with increasing frequency, which makes SHG phase matching more problematic, because both the increasing material index and the increasing cladding index must now be compensated by the effective-index difference between fundamental and second-order modes.
In the present work, only values of were studied. Since the values where phase matching occurs turn out to be of the order of , this corresponds to a strut width of a few tens of nm. It is unclear whether still narrower struts can be fabricated with a reasonable accuracy. With these values of , it was found that the frequency-doubled -like core mode could not be phase matched to the fundamental mode, at least not in a regime where both modes were reasonably well confined in the microstructure. On the other hand, the -like mode could be phase-matched in geometries where both the fundamental and second-order modes were reasonably well confined.
In Fig. 4 the phase matching curves and the sum of surface contributions to the SHG coefficients are reported. The bulk contributions were ignored, since the results for the silica nanowires indicated them to be of minor importance. Results are shown as a function of the structural parameter (see Fig. 1) for different values of the strut width relative to . For , phase matching to the mode could no longer be found. However, for the lower bridge widths, the variations with are within . Comparing to the ideal silica nanowire, the nonlinear coefficients for a given SHG wavelength are reduced by a factor of roughly 0.6–0.8 in the microstructures. This can be understood as a consequence of the phase matching appearing in a regime where the guided modes are less well confined because of the above mentioned influence of the struts on the SHG phase matching. Indeed, the effective areas of the phase matched modes in the microstructure are found to be larger than in the nanowire for the same phase matching wavelength. Also the effective index of the modes at phase matching is slightly lower in the microstructure, in spite of the presence of the struts, which should raise the effective cladding index. It must therefore be expected that structures with a larger number of struts, such as the well-known sixfold symmetric microstructures, will be more difficult to phase match, and will have lower values of the surface SHG coefficients.
3C. Impact of Structural Fluctuations
The central question of the present work is whether the nonlinear coefficients calculated in the previous subsections will allow significant power conversion fractions in realistic nanowires. Using Eq. (26) for a quick estimate, the SHG coefficients calculated for the microstructured nanowires lead to power conversion at the 10% level in fiber lengths of a few centimeters at a pump power level of . For nanosecond pulses, such peak power levels have been shown to be feasible in silica nanowires [10]. Also, nanowires with lengths of many centimeters have been demonstrated, with propagation losses of the order of or smaller [10]. However, a propagation length of the order of wavelengths sets rather severe restrictions on the tolerable phase mismatch . This implies that efficient power conversion to the second-harmonic frequency will only occur over a very narrow bandwidth. In Fig. 5 , the product of defined in Eq. (29) with the fiber length L is plotted as a function of the fiber diameter. It is evident that fiber lengths of the order of will lead to bandwidths of the order of . Therefore, the efficiency of the process will be highly sensitive to the uniformity of the nanowire. In this subsection, some simple propagation calculations will serve to illustrate the level of uniformity needed if serious power conversion fractions are to be achieved.
To describe propagation of pulses in a fluctuating fiber structure, Eq. (1) is generalized to
The pulse power is distributed on discrete frequency components , whose spacing is proportional to the inverse of the repetition rate (this is a different notation than the one used in Eqs. (3, 4)). Note thatas long as , are expanded to linear order in the frequency arguments. In the numerical simulations, the spectral distribution of the pump pulse is assumed to be an unchirped Gaussian. As a mathematical model for the structural fluctuations, the position-dependent phase matching wavelength is expressed as a Gaussian Fourier series with random-phase coefficientswhere is a random phase factor. One can then show thatthat is, can be regarded as a correlation length for the structural fluctuations. A typical realization with is shown in Fig. 6 . The function can be found from Eq. (28).In the numerical simulations, Eq. (30) is discretized in steps much smaller than . In each step is assumed constant, and Eq. (30) is propagated according to
where is the value of in the mth discretization segment. The calculations were done for an SHG wavelength of , with parameters for and appropriate for the microstructured nanowire with . Transform-limited Gaussian pump pulses with a temporal width of and a peak power of were used.For short propagation distances, fulfilling that
one has coherent SHG, and the second-harmonic power increases with . For longer distances, a transition to an incoherent regime occurs, where the SHG becomes an incoherent sum of contributions from different fiber segments, and the SHG power in this case scales with z [11]. In this limit, there are strong fluctuations in the SHG efficiency between different realizations of the random structure. An example is shown in Fig. 7 . The SHG power as a function of propagation length, averaged over 200 different realizations, is compared to the results for the particular realizations which lead to maximal and minimal SHG power at the output end. In Fig. 8 , curves averaged over 200 realizations are shown for different levels of structural fluctuations, and for fiber lengths of either 10 or . It is evident that the SHG power initially follows a dependence, which at some point breaks off , and at longer propagation distance becomes a linear dependence. The transition point depends on the magnitude of the structural fluctuations. This magnitude is quantified by , where is the maximal deviation of from the nominal value, and L is the fiber length. It can be seen that a fiber length of allows a conversion efficiency of (in the undepleted-pump approximation) for perfect phase matching, but that the conversion efficiency begins to deteriorate for , corresponding to . Since the SHG phase matching wavelength was found to scale almost linearly with the nanowire diameter, this corresponds to an atomic-scale roughness, which is unrealistic to maintain over longer fiber segments. For a fiber length of , a similar conversion efficiency is obtained for , but in this case must be lower. This is difficult, because one must expect that in a longer fiber the peak will also be larger.Experimentally, substantial efforts have gone into improving the uniformity of tapered nanowires. Shi et al. reported diameter fluctuations of in a nanowire with a diameter of [12], which would roughly translate into a figure of . Vukovic et al. studied tapering of microstructured fibers, and were able to maintain an outer-diameter fluctuation of over a length of [13]. Assuming that this level of accuracy carries over to the core dimensions and the phase matching wavelength, this would correspond to a fluctuation of for the fiber considered here, and thus . Calculations done at this level of fluctuations yielded averaged conversion efficiencies of for a fiber length of , and for a fiber length of , i.e., linear scaling in the incoherent limit. These figures indicate that power conversion factors in excess of 10% would require very long fibers while maintaining a very precise diameter control.
Since the lack of phase matching is due to random fluctuations rather than an inherent phase mismatch in the waveguide structure, the use of quasi-phase matching techniques is of little use to increase the conversion efficiency in this case. Instead, the conversion efficiency may be enhanced either by increasing the peak power of the pump pulses, or by increasing the surface nonlinear coefficient. The conversion efficiency will scale linearly with pulse power until the breakdown threshold is reached. Conversion efficiencies in excess of 10% in fibers shorter than 1 m will therefore require power levels beyond . It is unclear whether such high powers in sub-micrometer cores will be tolerable, even in very short pulses. On the other hand, if the nonlinear coefficient can be increased, it follows from Eq. (1) that the power conversion efficiency will scale with . Thus, increasing the surface nonlinearity by an order of magnitude compared to the values assumed for silica in this work would be enough to provide useful conversion efficiencies in realistic fibers. Such an increase in surface nonlinearity might be achieved by use of soft glasses with a higher polarizability than silica as a base material. Direct drawing of lead silicate microstructured fibers with core dimensions of the order of was recently demonstrated experimentally [14]. Another approach would be to keep silica as a base material, but coat the surface of the nanowire with a film of highly polarizable material, as was recently suggested by Xu et al. [15]. The potential of this idea can be appreciated by the following argument: The dominant term in the silica surface tensor has a magnitude of . One can think of this as if a surface layer with a thickness of had a dipole polarizability of . Using more polarizable materials, this polarizability could be enhanced by at least a factor of 2 [15], and the surface layer could be extended to a thickness of 10 nm or more without drastically changing the guiding properties of the fiber, especially for the longer pump wavelengths. Of course, the uniformity of the polarizable layer would be a crucial issue. Furthermore, an increase of the surface coefficient will in itself not be very useful if it happens at the expense of a similar reduction in the power threshold, a caveat which applies for both the soft-glass and the coated-silica approach.
An alternative application of nanowire SHG may be in sensing applications, since the surface second-order polarizability will be highly sensitive to adsorbed molecules on the silica surface. Indeed, surface SHG on bulk samples has long been used to characterize adsorbates, and a recent paper has suggested a similar application for silica microspheres [16, 17]. Finally, it should be mentioned that a very recent paper demonstrates the use of both second- and third-harmonic generation to characterize the diameter of a nanowire with an accuracy of about 2% [18].
4. CONCLUSION
In conclusion, the surface dipole and bulk multipole contributions to the second-order nonlinear coefficient in silica nanowires have been studied numerically, along with phase matching conditions and bandwidths. It is concluded that the length and loss figures of present-day nanowires allow for significant conversion efficiencies, but that nanowire uniformity is likely to be a severe limitation. Changing the fiber material or applying a polarizable surface coating may have the potential to overcome this limitation by providing an increase in the nonlinear coefficient of an order of magnitude. The use of microstructured fibers instead of circular nanowires is shown to be possible, although the nonlinear coefficients are slightly reduced, and -like modes are difficult to phase match.
1. U. Osterberg and W. Margulis, “Dye laser pumped by Nd:YAG laser pulses frequency doubled in a glass optical fiber,” Opt. Lett. 11, 516–518 (1986). [CrossRef] [PubMed]
2. R. Terhune and D. Weinberger, “Second-harmonic generation in fibers,” J. Opt. Soc. Am. B 4, 661–674 (1987). [CrossRef]
3. D. Anderson, V. Mizrahi, and J. Sipe, “Model for second-harmonic generation in glass optical fibers based on asymmetric photoelectron emission from defect sites,” Opt. Lett. 16, 796–798 (1991). [CrossRef] [PubMed]
4. R. Stolen and H. Tom, “Self-organized phase-matched harmonic generation in optical fibers,” Opt. Lett. 12, 585–587 (1987). [CrossRef] [PubMed]
5. V. Grubsky and J. Feinberg, “Phase-matched third-harmonic uv generation using low-order modes in a glass micro-fiber,” Opt. Commun. 274, 447–450 (2007). [CrossRef]
6. G. Brambilla, F. Xu, P. Horak, Y. Jung, F. Koizumi, N. Sessions, E. Koukharenko, X. Feng, G. Murugan, J. Wilkinson, and D. Richardson, “Optical fiber nanowires and microwires: fabrication and applications,” Adv. Opt. Photon. 1, 107–161 (2009). [CrossRef]
7. F. Rodriguez, F. X. Wang, and M. Kauranen, “Calibration of the second-order nonlinear optical susceptibility of surface and bulk of glass,” Opt. Express 16, 8704–8710 (2008). [CrossRef] [PubMed]
8. F. Rodriguez, F. X. Wang, B. Canfield, S. Cattaneo, and M. Kauranen, “Multipolar tensor analysis of second-order nonlinear optical response of surface and bulk of glass,” Opt. Express 15, 8695–8701 (2007). [CrossRef] [PubMed]
9. K. Okamoto, Fundamentals of Optical Waveguides (Academic, 2000).
10. S. Leon-Saval, T. Birks, W. Wadsworth, P. S. J. Russell, and M. Mason, “Supercontinuum generation in submicron fibre waveguides,” Opt. Express 12, 2864–2869 (2004). [CrossRef] [PubMed]
11. X. Vidal and J. Martorell, “Generation of light in media with a random distribution of nonlinear domains,” Phys. Rev. Lett. 97, 013902 (2006). [CrossRef] [PubMed]
12. L. Shi, X. Chen, H. Liu, Y. Chen, Z. Ye, W. Liao, and Y. Xia, “Fabrication of submicron-diameter silica fibers using electric strip heater,” Opt. Express 14, 5055–5060 (2006). [CrossRef] [PubMed]
13. N. Vukovic, N. G. R. Broderick, M. Petrovich, and G. Brambilla, “Novel method for the fabrication of long optical fiber tapers,” IEEE Photon. Technol. Lett. 20, 1264–1266 (2008). [CrossRef]
14. H. Ebendorff-Heidepriem, S. C. Warren-Smith, and T. M. Monro, “Suspended nanowires: fabrication, design and characterization of fibers with nanoscale cores,” Opt. Express 17, 2646–2657 (2009). [CrossRef] [PubMed]
15. Y. Xu, A. Wang, J. Heflin, and Z. Liu, “Proposal and analysis of a silica fiber with large and thermodynamically stable second order nonlinearity,” Appl. Phys. Lett. 90, 211110-1–3 (2007). [CrossRef]
16. J. Dominguez-Juarez, G. Kozyreff, and J. Martorell, “Surface nonlinear light generation in microresonators,” CLEO/Europe—EQEC 2009—European Conference on Lasers and Electro-Optics and the European Quantum Electronics Conference (IEEE, 2009), p. 1.
17. G. Kozyreff, J. Dominguez Juarez, J. Martorell, and J. Martorell, “Whispering-gallery-mode phase matching for surface second-order nonlinear optical processes in spherical microresonators,” Phys. Rev. A 77043817 (2008). [CrossRef]
18. U. Wiedemann, K. Karapetyan, C. Dan, D. Pritzkau, W. Alt, D. Meschede, and S. Irsen, “Measurement of submicrometre diameters of tapered optical fibres using harmonic generation,” Opt. Express 18, 7693–7704 (2010). [CrossRef] [PubMed]