High-power laser fiber system

Abstract

The invention relates a power fiber laser system including at least one single-mode fiber laser, emitting at a signal wavelength, the fiber including at least one outer cladding and a core, in which the core of the fiber has a radially graded index. The fiber includes, at least over a part of its length, a geometrical section having a graded fiber-core radius that decreases between an input end of the section and an output end of the section, the core radius and the index variation between the cladding and the fiber at the input end being such that the normalized frequency at the signal wavelength is less than the normalized cutoff frequency at which the fiber becomes unimodal.

Claims

1. A power fiber laser system comprising: at least one single-mode fiber laser, emitting at a signal wavelength, said fiber laser comprising at least one outer cladding and at least two cores, one of the at least two cores being configured as a central core and further configured such that a laser emission beam emerges, the central core being arranged in a tapered section of the fiber, a single reflector directly coupled to only one end of the central core and configured to provide feedback at said only one end of the central core, and remaining ones of the at least two cores are configured without a reflector, wherein each of the at least two cores of said fiber laser has a radially graded index and comprises, at least over a part of its length, a geometrical section, having a graded fiber-core radius that decreases between an input end of the geometrical section and an output end of the geometrical section, the graded fiber-core radius and an index variation between the cladding and the core at the input end being configured such that a normalized frequency of the signal wavelength is less than a normalized cutoff frequency at which the fiber laser becomes unimodal.

2. The power fiber laser system as claimed in claim 1, wherein the graded index has a parabolic profile.

3. The power fiber laser system as claimed in claim 1, wherein the graded index has a triangular profile.

4. The power fiber laser system as claimed in claim 1, wherein the fiber laser comprises a section of conical geometry, a first, cladding radius (r.sub.1) and a second, core radius (r.sub.2), the first and second radii being defined as a function of length z, respectively, by the following equations:
r.sub.1(z)=r.sub.1c(1m.sub.1z)
r.sub.2(z)=r.sub.2c(1m.sub.2z), wherein r.sub.1c and r.sub.2c are the initial radius; and wherein m.sub.1 and m.sub.2 are the conicity.

5. The power fiber laser system as claimed in claim 1, wherein the at least two cores are distributed and configured so as to couple modes in said at least two cores.

6. The power fiber laser system as claimed in claim 5, wherein the at least two cores are distributed and configured in a geometrical array in a plane perpendicular to the length of said at least one single-mode fiber laser.

7. The power fiber laser system as claimed in claim 5, wherein the at least two cores are distributed in a circle.

8. The power fiber laser system as claimed in claim 1, further comprising a preform comprising an assembly of individual fiber claddings and cores.

9. The power fiber laser system as claimed in claim 1, further comprising: a fiber pigtail configured to receive pump power from a laser diode; and a bundle of radially coupled fibers comprising a stretched construction so as to optimize a coupling therebetween, wherein the single mode fiber laser is configured such that the signal wavelength of the laser emission beam emerges from the central core.

10. The power fiber laser system as claimed in claim 5, wherein the single reflector is configured to implement a mirror function.

11. The power fiber laser system as claimed in claim 10, wherein the single reflector includes a Bragg reflector.

12. The power fiber laser system as claimed in claim 1, further comprising axicons coupled to opposite ends of the at least one single-mode fiber laser comprising conical sections.

13. The power fiber laser system as claimed in claim 1, further comprising at least one fiber having a silica cladding.

14. The power fiber laser system as claimed in claim 1, further comprising at least one fiber comprising a GeO.sub.2-doped silica fiber core.

15. The power fiber laser system as claimed in claim 14, wherein a fiber core of the at least one fiber comprises peripheral layers containing varying amounts of GeO.sub.2 so as to ensure a radially graded index within the fiber core.

16. The power fiber laser system as claimed in claim 1, further comprising at least one fiber having a fiber core made of silica doped with rare-earth ions.

17. The power fiber laser system as claimed in claim 6, wherein the single reflector is configured to implement a mirror function.

18. The power fiber laser system as claimed in claim 17, wherein the single reflector is configured to implement a Bragg reflector.

19. A power fiber laser system comprising: a single-mode fiber laser configured to emit a signal wavelength comprising: a fiber including at least one outer cladding and at least two cores, one of the at least two cores being configured and arranged as a central core from which a laser emission beam emerges, the central core being arranged in a tapered section of the fiber, a single reflector directly coupled to only one end of the central core and configured to provide feedback at said only one end of the central core and remaining ones of the at least two cores are configured to operate without a reflector, each said core of the fiber comprising a radially graded index that comprises, at least over a part of its length, a geometrical section having a graded fiber-core radius that decreases homothetically between an input end and an output end, and for each said core, a core radius and the index variation between the cladding and the fiber at the input end being configured such that a normalized frequency of the signal wavelength is less than the normalized cutoff frequency at which the fiber becomes unimodal.

20. The power fiber laser system of claim 19, wherein the single reflector is a mirror or a Bragg reflector.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The invention will be better understood and other advantages will become clear on reading the following nonlimiting description and by virtue of the appended figures in which:

(2) FIG. 1 shows the various optical and geometrical parameters that characterize an optical fiber;

(3) FIG. 2 shows schematically a power fiber laser, used in a laser system according to the invention;

(4) FIG. 3a shows the variation in the radius of the HE.sub.11 mode as a function of the distance in a tapered fiber structure;

(5) FIG. 3b shows the irradiance distribution as a function of the radial distance and of the propagation distance through the tapered fiber.

(6) FIG. 3c shows the irradiance at z=1 m for a fiber having a 50% conicity from a 100 m initial diameter;

(7) FIG. 4 shows the irradiance distribution as a function of distance in three 100 m initial diameter fibers comprising tapered sections with a 50% conicity over 1 m;

(8) FIG. 5 shows a graph of the variation in the radius of the mode normalized to the core radius of a fiber as a function of the normalized frequency parameter and as a function of the index gradient;

(9) FIG. 6 shows a graph of the variation in the radius of the HE.sub.11 as a function of the structure in a parabolic-profile tapered fiber;

(10) FIG. 7 shows a graph of the variation in the irradiance as a function of the distance in the conical part of a graded-index fiber;

(11) FIG. 8 shows a graph of the variation in the irradiance
|a.sub.1(z)|.sup.2 at a given wavelength;

(12) FIG. 9 shows a two-fiber linear coupler;

(13) FIG. 10 shows the variation in the coupling coefficient between the two fibers as a function of the normalized distance separating the two waveguides in the spectral range of the (laser) signal;

(14) FIGS. 11a, 11b and 11c show the behavior of the fiber coupler as a function of the wavelength, the propagation distance and the normalized spacing between the two optical-fiber cores;

(15) FIG. 12a shows the variation in the local coupling coefficient as a function of the distance z;

(16) FIG. 12b shows the effective coupling factor
C.sub.effz as a function of the distance z into the tapered part;

(17) FIG. 12c shows the normalized irradiance distribution of the excited waveguide in a linear coupler configuration with a conical structure at 1.07 m (curve 12c.sub.0) and at 1.09 m (curve 12c.sub.1) as a function of the propagation distance;

(18) FIG. 13a shows a embodiment of the invention in which the fiber lasers are assembled with a square symmetry;

(19) FIGS. 13b and 13c show coupling between waveguides as a function of the propagation distance and the wavelength;

(20) FIG. 13d shows the irradiance distribution in the waveguide 1 as a function of z for various wavelengths;

(21) FIG. 14a shows a symmetrical geometry consisting of a stack of cylindrically configured fibers;

(22) FIG. 14b shows the variation in the irradiance in the central fiber and the irradiance obtained in one of the peripheral coupling fibers as a function of the propagation distance with
s=5 and at 1.07 m and FIG. 14c gives a general idea of the spectral distribution at the output of the fibers (i.e. central and peripheral fibers);

(23) FIG. 15a shows a perspective view of an array of fibers, distributed cylindrically, used in a embodiment of the invention;

(24) FIG. 15b shows a cross section of another embodiment of the invention, in which a preform comprising an assembly of individual fiber claddings and cores;

(25) FIG. 16a shows a cross section of an embodiment of a laser system according to the invention;

(26) FIG. 16b shows another embodiment of the invention comprising an axicon;

(27) FIG. 17 shows the variation in the refractive index of GeO.sub.2-doped silica; and

(28) FIG. 18 shows a fiber production step implemented in the invention.

DETAILED DESCRIPTION

(29) The essence of the invention lies in the use of fiber, the core of which has a radial index profile, to produce a power fiber laser system. This profile may notably be parabolic or triangulardetails of the behavior obtained with these configurations will be given in the description below.

(30) Moreover, it may be of particular benefit to combine such a radial index profile with a longitudinal fiber-core gradient.

(31) FIG. 2 thus shows a schematic of a power fiber laser used in a laser system according to the invention. This fiber is intended to be supplied with laser power PL so as, at the out of the fiber, to generate a laser beam emitting at a signal wavelength .

(32) This fiber comprises a section T having an input end Exe and an output end Exs, a longitudinal core-radius gradient and a radial index profile

(33) The analysis below pertains to the properties obtained with such a decreasing fiber core radius gradient:

(34) As is known, although the solution of Maxwell's equations applied to a cylindrically-symmetric two-dielectric structure leads to an expression for the radial field distribution in the form of Bessel functions and modified Bessel functions of the first kind, this field distribution may be approximated by a Gaussian function with a 1/e mode radius (for the field but 1/e.sup.2 for the irradiance) having the form:

(35) $w_0\left(\right)=a\left(0.65+\frac{1.619}{V^{3/2}\left(\right)}+\frac{2.879}{V^6\left(\right)}\right).$

(36) This function, formulated by D. Marcuse and notably described in the article Loss analysis of single-mode fiber splices, Bell Syst. Tech. J., Vol. 56, No. 5, 1977, p. 703-719, is applicable in the range
0.8<V<2.8.
A variant allows the propagation of the mode to be directly connected to the signal wavelength (the laser in the case considered) for a given cutoff wavelength .sub.c.

(37) This wavelength .sub.c defines the spectral range in which only the groups of HE11 modes with their two polarization eigenstates (LP01) can propagate without loss (in the electromagnetical sense). For smaller wavelengths, the groups of LP.sub.11 modes may exist and propagate with the group of LP.sub.01 modes (multimodal regime). This relationship is the following:

(38) $w_0\left(\right)a\left(0.65+0.434{\left(\frac{}{{}_c}\right)}^{3/2}+0.0149{\left(\frac{}{{}_c}\right)}^6\right).$

(39) By way of example, it is possible to calculate the variation in the mode radius from the optical and geometrical parameters that characterize an optical fiber.

(40) In the presence of a linear taper, the size of the mode increases as a function of the ratio /c whereas the core radius decreases.

(41) Thus, by way of example, it is possible to consider a fiber laser emitting a signal wave at =1.07 m and having a core radius and index variation =103 such that the initial cutoff wavelength is equal to C=0.75 m. This cutoff wavelength C=0.75 m corresponds to a normalized frequency parameter V=2.4048 and to a mode radius w0(0.75 m)=4.85 m.

(42) At the signal wavelength, a parameter V equal to V=1.68 and a mode size (radius) w0(1.07 m)=6.715 m are obtained.

(43) These parameters lead to an HE11 mode weakly confined to the region covered by the core radius of the fiber.

(44) Typically, the fiber has a cladding diameter of 100 m. It should be noted that, generally, fiber lasers called double-cladding fiber lasers are characterized by a cladding diameter of about 400 m, so as to couple pump diodes via an assembly of multimode couplers.

(45) According to a first embodiment proposed by the invention, the fiber-core radius gradient may be a conical profile, that is to say that, typically, when tapered over an arbitrary distance of 1 m, so as to obtain adiabatic coupling with a conicity of 50%, corresponding to variation from an initial cladding diameter of 100 m to 50 m, it is possible to calculate the variation in the various radii (of the core and cladding) homothetically:

(46) $\begin{array}{c}b\left(z\right)\\ z\left(z\right)\end{array}\}=\{\begin{array}{c}{b}_C\left(1-0.5z\right)\\ a_C\left(1-0.5z\right)\end{array}$

(47) The variation in the mode radius as a function of distance is then shown in FIG. 3a, the lines 3a.sub.1 and 3a.sub.2 showing, respectively, the variation in the core radius and the cladding radius as a function of the distance z.

(48) At one end of the tapered part, the normalized frequency parameter is equal to V=0.84, a value close to the limit of validity of the formula defining the Gaussian model of the modal distribution of the HE.sub.11 mode or group of LP.sub.01 modes. The curves of FIG. 3a clearly show that over a z1 m length of section of conical profile, the 1/e field radius of the mode is approximately equal to the radius of the cladding. Beyond this, the mode is subjected to diffraction and refraction effects leading to beam losses.

(49) It is known that the optical damage threshold of a dielectric medium (in terms of fluence (J/m2) or in terms of irradiance (W/m2)) is always less at the interface compared to that characterizing the bulk behavior. The conical geometry thus allows this limitation to be partly overcome, in particular, if it desired to maximize the optical power that a laser or an optical-fiber amplifier system can deliver.

(50) Thus, two criteria need to be fulfilled: single-mode optical fibers having a low cutoff wavelength for the group of LP11 modes compared to the signal wavelength must be developed that also have a conicity at the end of each fiber laser.

(51) It is then possible to evaluate the effects of a taper on the field, assuming an irradiance proportional to the field modulus squared:
I(r,z)|E|.sup.2.

(52) To calculate the radial extension of the mode, the irradiance distribution is calculated. To do this, the power transmitted by the mode is defined by the formula:

(53) $P_T=I_0\left(z\right){}_0^2{}_0^{}{e}^{-2{\left(\frac{r}{w_0\left(z\right)}\right)}^2}rdrd=I_0\frac{w_0^2\left(z\right)}{2}$

(54) and the distribution of the irradiance (W/m.sup.2) is thus:

(55) $I\left(r,z\right)=\frac{2P_T}{w_0^2\left(z\right)}{e}^{-2{\left(\frac{r}{w_0\left(z\right)}\right)}^2}.$

(56) The irradiance distribution as a function of the radial distance and the propagation distance through the tapered fiber is shown in FIG. 3b, taking the transmitted optical power to be 1 W.

(57) At the distance z=1 m it is observed that the power distribution is not contained within the cladding diameter (50 m in the case shown).

(58) FIG. 3c shows the irradiance distribution at this point and for the proposed conicity.

(59) A straight-forward calculation gives the percentage of the power contained in a centered circular opening of radius b:

(60) ${}_g=\frac{{}_0^b{e}^{-2{\left(\frac{r}{w_0}\right)}^2}\mathrm{rdr}}{{}_0^{}{e}^{-2{\left(\frac{r}{w_0}\right)}^2}\mathrm{rdr}}=1-e^{-2{\left(\frac{b}{w_0}\right)}^2}.$

(61) On the basis of these preliminary data, it is possible to analyze a combination of several fibers in tape form, for example three fibers that undergo collectively a tapering of their ends. Considering the case of uncoupled fibers, i.e. there is no mode coupling, a description of the irradiance distribution can be obtained, as shown in FIG. 4.

(62) Thus, an overall low power density and a reduced effective emission area are obtained at the fiber output. Such an architecture therefore allows the optical flux restrictions to be met whilst increasing the radiance of the emission source.

(63) The essence of the invention lies in the use of fiber, the core of which has an index profile. This profile may notably be parabolic or triangular.

(64) It is possible to consider a Gloge-Marcatili index gradient distribution, described notably in Multimode theory of graded core fibers by D. Gloge and E. A. J. Marcatili, Bell Sys. Tech. J., Vol. 52, 1973, p. 1563-1578:

(65) $n^2\left(r\right)=\{\begin{array}{cc}{n}_1^2\left(1-2{\left(r/a\right)}^a\right)& ra\\ n_1^2\left(1-2\right)& .Math.r.Math.>a;\end{array}$

(66) For =2 (parabolic profile), the normalized cutoff frequency parameter of the mode group adjacent to the LP01 mode (unimodal) group is equal to Vc=3.518, as described in the article Cut-off frequency in radially inhomogeneous single-mode fiber by W. A. Gambling, D. N. Payne and H. Matsumura, Electr. Letters 13 (5) Mar. 3, 1977, p. 139, as compared to a parameter value of Vc=2.4048 for a step-index fiber
(=).

(67) The core radius is thus increased by the ratio Vc(=2)/Vc(=)=1.463. This ratio is greater still if a triangular profile is considered for which Vc=4.381 and Vc (=1)/Vc(=)=1.8218.

(68) A more general formula giving the variation in the mode radius as a function of the parameter V for a fiber with a parabolic profile is given by D. Marcuse in Gaussian approximation of the fundamental modes of graded-index fibers, JOSA 68 (1) January 1978, p. 103.

(69) 0$\frac{w_0}{a}=\frac{A}{V^{2/+2}}+\frac{B}{V^{3/2}}+\frac{C}{V^6}$$A={\left(\frac{2}{5}\left(1+4{\left(\frac{2}{}\right)}^{5/6}\right)\right)}^{1/2}$$B=e^{0.298/}-1+1.478\left(1-e^{-0.077}\right)$$C=3.76+e^{4.19/{}^{0.418}}$

(70) The above set of parameters and equations validates the concept of the present invention, proposing to produce laser sources with this type of index gradient.

(71) Not only does this type of structure enable the threshold at which nonlinear effects appear and optical damage occurs to be increased, but it also enables losses related to bending or microbending to be reduced by increasing the index difference (i.e. the optical guiding), since Vc is increased.

(72) FIG. 5 shows the variation in the mode radius normalized to the core radius of a fiber as a function of the normalized frequency parameter and as a function of the index gradient: curve 5a is for a coefficient =1, the curve 5b for =2, and curve 5c is for an infinite a coefficient.

(73) A formula similar to that developed above for a step-index fiber but for a tapered parabolically-graded fiber gives the variation in the mode radius as a function of z.

(74) In this case, the relative index difference is =4103 (four times larger than for a step-index fiber) and the cutoff wavelength is taken to be Vc=3.518 for a wavelength of 0.8 m.

(75) FIG. 6 shows the variation in the HE11 mode radius, in microns, as a function of the distance into a tapered structure with a parabolic profile; also shown are the core radius (curve a) and the cladding radius (curve b).

(76) FIG. 7 shows the variation in the irradiance as a function of the distance into the conical part of a fiber with a parabolic profile, the parameters of which are the following: a radius a=3.446 m; a variation =4103; and a cutoff wavelength .sub.c=0.8 m. At the signal wavelength of the laser (=1.07 m) the initial normalized frequency parameter is equal to V=2.623 for an identical conicity to that used in the analysis of a tapered step-index fiber.

(77) Thus it may be seen that, using a fiber laser with a parabolic index-profile leads not only to strong confinement at the laser wavelength, reducing the threshold for coherent emission, but also to coupling with a conical part that enables the optical damage thresholds to be increased at the ends.

(78) System according to the invention comprising several coupled fiber lasers

(79) In the case of two parallel waveguides with identical optical and geometrical parameters (modes characterized by one and the same propagation constant ) spaced such that the evanescent parts are able to partially intercept the adjacent waveguide, in the absence of differential losses, the coupling equations may be written as:

(80) $\frac{}{z}.Math.\begin{array}{c}{a}_1\left(z\right)\\ a_2\left(z\right)\end{array}.Math.+i.Math.\begin{array}{cc}& C\\ C& \end{array}.Math..Math.\begin{array}{c}{a}_1\left(z\right)\\ a_2\left(z\right)\end{array}.Math.=.Math.\begin{array}{c}0\\ 0\end{array}.Math.$

(81) where a.sub.j(z) corresponds to the amplitude of the complex fields transmitted by the waveguides,
j=1,2

(82) and C is the coupling coefficient between waveguides.

(83) The solution of these equations for incident in-phase amplitudes .sub.1(0) and .sub.2(0) is:
a.sub.1(z)=a.sub.1(0)cos(Cz)ia.sub.2(0)sin(Cz)
a.sub.2(z)=a.sub.2(0)cos(Cz)ia.sub.1(0)sin(Cz)

(84) giving a modulus |a.sub.j|.sup.2 corresponding to the irradiance of:
|a.sub.1(z)|.sup.2=a.sub.1.sup.2(0)cos.sup.2(Cz)+a.sub.2.sup.2(0)sin.sup.2(Cz)
|a.sub.2(z)|.sup.2=a.sub.1.sup.2(0)cos.sup.2(Cz)+a.sub.1.sup.2(0)sin 2(Cz)

(85) When a single field is incident, the equation for the linear coupler are:
|a.sub.1(z)|.sup.2=a.sub.1.sup.2(0)cos.sup.2(Cz)
|a.sub.2(z)|.sup.2=a.sub.1.sup.2(0)sin.sup.2(Cz)

(86) When there are two incident fields, but with a phase difference of between them:
a.sub.1(0)=a.sub.1
a.sub.2(0)=a.sub.2e.sup.i
the following solutions are obtained:
|a.sub.1(z)|.sup.2=a.sub.1.sup.2(0)cos.sup.2(Cz)a.sub.1(0)a.sub.2(0)sin(2Cz)sin()+a.sub.2.sup.2(0)sin.sup.2(Cz)
|a.sub.2(z)|.sup.2=a.sub.2.sup.2(0)cos.sup.2(Cz)+a.sub.1(0)a.sub.2(0)sin(2Cz)sin()+a.sub.1.sup.2(0)sin.sup.2(Cz)

(87) A graph of the variation in the irradiance
|a.sub.1(z)|.sup.2
at a given wavelength is plotted in FIG. 8 with the following parameter settings: .sub.1(0)=1, .sub.2(0)=1 and C=1.

(88) Generally, the linear coupler makes use of the ability to address the output power as a function of the excitation of one of the two waveguides. For this application, only the phase relationships between two waves simultaneously coupled at both ends are of interest. On the other hand, if it were desired to obtain coherent transfer, the analysis of the behavior of the coupler would need to take into account the relative initial phase between the two fields, explaining the above analysis.

(89) Still more generally, in the presence of two coupled waveguides, characterized by a difference
=.sub.1.sub.2
between the propagation constants of the structures, the following are obtained:

(90) ${.Math.a_1\left(z\right).Math.}^2=\frac{1}{{\left(\frac{}{2}\right)}^2+C^2}\left[\begin{array}{c}{C}^2{\cos }^2\left(a_1^2-a_2^2\right)-a_1{a}_{2}C\cos {\sin }^2+{\left(\frac{}{2}\right)}^2{a}_1^2{a}_2^2{C}^2\\ -2a_1{a}_{2}C\sin \sqrt{C^2+{\left(\frac{}{2}\right)}^2}\cos \sin \end{array}\right]$${.Math.a_2\left(z\right).Math.}^2=\frac{1}{C^2+{\left(\frac{}{2}\right)}^2}[\begin{array}{c}-C^2{\cos }^2\left(a_1^2-a_2^2\right)+a_1{a}_{2}C\cos {\sin }^2+{\left(\frac{}{2}\right)}^2{a}_2^2+a_1^2{C}^2\\ +2a_1{a}_{2}C\sin \sqrt{C^2+{\left(\frac{}{2}\right)}^2}\cos \sin \end{array}]$$\mathrm{where}$$=z\sqrt{C^2+{\left(\frac{}{2}\right)}^2}.$

(91) For =0 and in the presence of a single field incident on the coupler (a.sub.2(0)=0), the conventional equation for a directional coupler with different propagation constants:

(92) $\frac{{.Math.a_2\left(z\right).Math.}^2}{{.Math.a_1\left(0\right).Math.}^2}=\frac{C^2}{C^2+{\left(\frac{}{2}\right)}^2}{\sin }^2\left(\sqrt{{C^2\left(\frac{}{2}\right)}^2}z\right).$

(93) For two identical cylindrical axisymmetric dielectric waveguides, shown in FIG. 9, the coupling coefficient C is given by:

(94) $C=\frac{\sqrt{2}}{a}\frac{U^2}{V^3}\frac{K_0\left(\mathrm{Wd}/a\right)}{K_1^2\left(W\right)}$
where d is the distance between the two cores of radius a and U and W are the transverse propagation constants in the core and in the cladding respectively of the optical fiber.
K.sub.0(x)
and
K.sub.1(x)
are 0-order and 1.sup.st-order modified Bessel functions of the second kind that have x as their argument.

(95) The parameters U and W result from the solution of Maxwell's equations applied to a cylindrical waveguide structure with two dielectrics.

(96) $\frac{U^2}{a^2}=k^2{n}_1^2-{}^2$$\frac{W^2}{a^2}={}^2-k^2{n}_2^2$

(97) where is the propagation constant associated with the electromagnetic mode able to propagate in the waveguided regime. It is from these equations that the normalized frequency parameter is defined:
V.sup.2=U.sup.2+W.sup.2=k.sup.2a.sup.2n.sub.1.sup.22.

(98) From the optical and geometrical properties of the step-index fibers it is possible to calculate the amplitude of the coupling coefficient as a function of the normalized distance s=d/a between the two longitudinal guiding axes. In the case of a fiber laser for which .sub.c=0.75 m, =10.sup.3 and a=4.41 m, the variation in the coupling coefficient between the two fibers shown in FIG. 10 is obtained, shown as a function of the normalized distance between the two waveguides in the spectral range of the (laser) signal.

(99) FIGS. 11a, 11b and 11c show the behavior of the fiber coupler as a function of the wavelength, the propagation distance and the normalized distance between the two cores of the optical fibers.

(100) More precisely, FIG. 11a shows the variation in the power transmitted by the excited waveguide of the linear multifiber coupler as a function of s=d/a and of the propagation distance at 1.07 m, FIG. 11b shows the power at the point z=0.2 m as a function of the wavelength and as a function of the parameter s=d/a and FIG. 11c shows the spectral distribution at the end of the excited waveguide as a function of the distance: z=0.2 m (curve 11c4) at 1 m (curve 11c0) for s=4, with the following intermediate curves 11c3, 11c2, 11c1.

(101) In the case shown, initially an optical wave is coupled to only one of the two waveguides and the coupling is given by:

(102) $\frac{{.Math.a_1\left(z\right).Math.}^2}{{.Math.a_1\left(0\right).Math.}^2}={\cos }^2\left(\mathrm{Cz}\right).$

(103) In the context of the present invention, it is convenient to consider a linear coupler with a tapered structure. Such an architecture has a coupling coefficient that varies with the propagation distance.

(104) For two coupled fibers with a conicity of 50%, the coupling coefficient depends on z through the z-dependence of the various parameters:

(105) $C\left(z\right)=\frac{\sqrt{2}}{a\left(z\right)}\frac{{U\left(z\right)}^2}{{V\left(z\right)}^3}\frac{K_0\left(W\left(z\right)d\left(z\right)/a\left(z\right)\right)}{K_1^2\left(W\left(z\right)\right)}.$

(106) Since
V(z)=ka(z)n.sub.1{square root over (2)}
it follows that the transverse propagation constants are also functions of z:
V(z).sup.2=U(z).sup.2+W(z).sup.2=k.sup.2a(z).sup.2n.sub.1.sup.22.

(107) It turns out that the ratio
d(z)/a(z)
is constant whatever z for a linear, and therefore homothetic (constant conicity), stretch.

(108) FIG. 12a shows the variation in the local coupling coefficient (in m1) as a function of distance z and the signal wavelength in m.

(109) To evaluate the effect of coupling one waveguide to another, it is convenient to calculate the effective coupling factor as a function of z, i.e. the integral:

(110) $C_{\mathrm{eff}}z={}_0^{z}C\left(z^{}\right){\mathrm{dz}}^{}.$

(111) FIG. 12b shows the effective coupling factor C.sub.eff.sup.z as a function of the distance z into the tapered part.

(112) With these longitudinal variations in the coupling coefficient, the irradiance of the excited waveguide varies as shown in FIG. 12c and according to the following formula:

(113) $\frac{{.Math.a_1\left(z\right).Math.}^2}{{.Math.a_1\left(0\right).Math.}^2}={\cos }^2\left({}_0^{z}C\left(z^{}\right)\mathrm{dz}\right).$

(114) FIG. 12c shows the normalized irradiance distribution in the excited waveguide in a linear coupler configuration with a conical structure at 1.07 m (curve 12c0) and at 1.09 m (curve 12c1) as a function of the propagation distance.

(115) The concept of multichannel optical fibers may be extended and the coupling equations in the absence of losses may be written as:

(116) 0$\frac{{\mathrm{da}}_m}{\mathrm{dz}}+{}_m{a}_m=-\underset{nm}{\underset{n=1}{\overset{N}{.Math.}}}{C}_{\mathrm{mn}}{.Math.}_n$

(117) Exemplary embodiments of a multichannel multifiber laser system according to the invention:

(118) According to one embodiment of the invention, the fiber lasers are bundled with a square symmetry, as shown in FIG. 13a. In this configuration it is convenient to distinguish two inter-fiber coupling coefficients: a direct coupling coefficient C.sub.L and a transverse coupling coefficient C.sub.T.

(119) The coupling equations are written as:

(120) $\frac{{\mathrm{da}}_1}{\mathrm{dz}}+i{a}_1=-i\left(a_2+a_4\right)C_L-{\mathrm{ia}}_3{C}_T$$\frac{{\mathrm{da}}_2}{\mathrm{dz}}+i{a}_2=-i\left(a_1+a_3\right)C_L-{\mathrm{ia}}_4{C}_T$$\frac{{\mathrm{da}}_3}{\mathrm{dz}}+i{a}_3=-i\left(a_2+a_4\right)C_L-{\mathrm{ia}}_1{C}_T$$\frac{{\mathrm{da}}_4}{\mathrm{dz}}+i{a}_4=-i\left(a_1+a_3\right)C_L-{\mathrm{ia}}_2{C}_T$

(121) The waveguides may be considered to be identical and the coupling coefficient independent of the polarization state. The characteristic equation is obtained by substituting the function:
a.sub.j=B.sub.je.sup.ize.sup.iz

(122) where B.sub.j and .sub.j are the elements of the eigenvectors and the eigenvalues respectively.

(123) The following matrix equation is obtained:

(124) $i{e}^{-i\left(+\right)z}.Math.\begin{array}{cccc}-& C_L& C_T& C_L\\ C_L& -& C_L& C_T\\ C_T& C_L& -& C_L\\ C_L& C_T& C_L& -\end{array}.Math..Math.\begin{array}{c}{B}_1\\ B_2\\ B_3\\ B_4\end{array}.Math.=0$
the characteristic equation of which is:

(125) $.Math.\begin{array}{cccc}-& C_L& C_T& C_L\\ C_L& -& C_L& C_T\\ C_T& C_L& -& C_L\\ C_L& C_T& C_L& -\end{array}.Math.=0$

(126) The four eigenvalues are:
.sub.1=C.sub.T2C.sub.L
.sub.2=C.sub.T+2C.sub.L
.sub.3=.sub.4=C.sub.T
and the associated eigenvectors are:

(127) ${\stackrel{.fwdarw.}{B}}_1=.Math.\begin{array}{c}-1\\ 1\\ -1\\ 1\end{array}.Math.{\stackrel{.fwdarw.}{B}}_2=.Math.\begin{array}{c}1\\ 1\\ 1\\ 1\end{array}.Math.{\stackrel{.fwdarw.}{B}}_3=.Math.\begin{array}{c}-1\\ 0\\ 1\\ 0\end{array}.Math.{\stackrel{.fwdarw.}{B}}_3=.Math.\begin{array}{c}0\\ -1\\ 0\\ 1\end{array}.Math.$

(128) Thus, four possible modes are defined.

(129) In the presence of a single incident field the irradiance distributions are given by the formulae:

(130) $I_1\left(z\right)=\frac{1}{4}\left(1+2\cos \left(2C_{L}z\right)\cos \left(2C_{T}z\right)+{\cos }^2\left(2C_{L}z\right)\right)$$I_2\left(z\right)=\frac{1}{4}{\sin }^2\left(2C_{L}z\right)$$I_3\left(z\right)=\frac{1}{4}\left(1-2\cos \left(2C_{L}z\right)\cos \left(2C_{T}z\right)+{\cos }^2\left(2C_{L}z\right)\right)$$I_4\left(z\right)=\frac{1}{4}{\sin }^2\left(2C_{L}z\right)$$\mathrm{with}$$C_T\left(s\right)=C_L\left(s/\cos \left(/4\right)\right).$
The irradiance distributions were calculated as a function of the propagation distance in this example for a linear multifiber coupler with a ratio s=5. FIGS. 13b and 13c show the coupling between waveguides as a function of the propagation distance and the wavelength. More precisely, FIG. 13b shows the variation in the irradiances at 1.06 m as a function of z for a unit of power injected into the waveguide 1 (curve 13b.sub.1). The irradiances in 2 and 4 are given by the curve 13b.sub.2 and the irradiance in the waveguide 3 is shown by the curve 13b.sub.3. FIG. 13c shows the variation in the irradiances at 1.09 m as a function of z for a unit of power injected into the waveguide 1 (curve 13c.sub.1). The irradiances in 2 and 4 are given by the curve 13c.sub.24 and the irradiance in the waveguide 3 is shown by the curve 13c.sub.3.

(131) The transfer behavior between waveguides 1 and 3 requires two coupling coefficients (C.sub.T and C.sub.L)

(132) FIG. 13d shows more clearly the behavior of the irradiance distributions in the excited waveguide as a function of the wavelength, the irradiance distribution in waveguide 1 as a function of z being given both for =1.06 m (curve 13d.sub.1) and =1.09 m (curve 13d.sub.2).

(133) According to another embodiment of the invention, the laser system may have a symmetrical geometry consisting of a stack of cylindrically configured fibers.

(134) This configuration is shown in FIG. 14a. Typically, the initial irradiance may be carried by the central fiber. In this case, the irradiances transmitted by fibers 1 to 6 result from coupling between fibers separated by a common distance s (arranged in an isosceles triangle). The irradiances are given by:

(135) $I_0\left(z\right)=1-\frac{6}{7}{\sin }^2\left(\sqrt{7}{C}_{L}z\right)$$I_n\left(z\right)=\frac{1}{7}{\sin }^2\left(\sqrt{7}{C}_{L}z\right)$$n=1,\mathrm{.Math.},6.$

(136) FIG. 14b shows the variation in the irradiance in the central fiber and the variation in the irradiance in one of the coupled peripheral fibers as a function of the propagation distance with s=5 and at 1.07 m and FIG. 14b also shows the irradiance distribution in the central fiber and in one of the six fibers on the periphery with s=5 at 1.07 m as a function of the propagation distance.

(137) FIG. 14c shows the spectral distribution of the irradiance in the central fiber and in one of the six fibers on the periphery with s=5 as a function of the propagation distance: z0=0.2 m, z1=0.4 m, z2=0.6 m, z3=0.8 m and z4=1 m. The pairs of curves, respectively (14cz0, 14cz0), (14cz1, 14cz1), (14cz2, 14cz2), (14cz3, 14cz3) and (14cz4, 14cz4), give the spectral distribution at the output of the central fiber, and at the output of an adjacent fiber.

(138) To produce this type of configuration, one possible geometry consists in taking a bundle of fiber lasers that are stretched at the end. To do this, it is possible to assemble N fibers fi around a central fiber fic, one end of the fibers having been stretched so as to obtain a structure with multiple couplings. FIG. 15a shows a perspective view of such an assembly.

(139) It is also possible to produce a preform made of several rods; typically they may be made of silica locally doped with a rare earth as shown in FIG. 15b which shows a cross section of such a preform P and its peripheral fibers fi.sub.j around a central fiber fi.sub.c. The central dual-index rod may be doped or undoped so as to provide the output coupler channel for the laser emission. The bundle is then drawn into fibers, the distance between cores being homothetically maintained, and then one of the ends is stretched.

(140) Such a multiple-core fiber is well suited to efficient pumping of the doped regions as it makes use of the modes of the pump waveguiding-structure which have a nonzero azimuthal dependence.

(141) The geometry of the multifiber laser may take the form shown in FIG. 16a, having: a fiber pigtail intended to receive pump power from a laser diode; a bundle of radially coupled fibers stretched so as to optimize the coupling; and a central fiber from which the laser beam emerges.

(142) At the end of this fiber may be placed a Bragg reflector providing a feed-back mirror function or an external mirror. A self-organization occurs depending on the forward and return paths in the bundle of laser cavities

(143) According to another embodiment, the laser system comprises: a parabolic or triangular graded-index core structure; a bundle of multicore couplers (monolithic in the form of a multifiber or specific preform), the bundle being stretched at one end; and an optical axicon at the other end.

(144) Such an architecture enables the polarization state to be controlled, this being the equivalent of providing specific couplers and multiplexing the polarization states, as shown in FIG. 16b.

(145) Generally, to produce fibers having, at least over one section, a tapered region, it is possible to use a chemical etch (for example employing hydrofluoric acid HF).

(146) It is also possible to envisage welding a fiber with a dedicated undoped end to obtain the desired profiles (provided that the two guiding structures are optimized so as to achieve mode matching).

(147) According to one embodiment of the invention, the graded-index core fibers may be produced using the method described below.

(148) It is known that the index of a fiber core made of silica (SiO2) can be increased by incorporating into it molar percentages of, for example, germanium oxide. FIG. 17 thus shows the change in the refractive index of the silica composite caused by a gas stream carrying Ge, P or B cations.

(149) Thus, by implementing chemical vapor deposition, it is possible to deposit, within a rod, internal layers, the index of which varies depending on the concentration of the stream carrying a varying percentage of Ge cations, enabling thereby the desired radially graded-index of the present invention to be obtained, as shown in FIG. 18.