Optical fiber waveguide having minimum modal dispersion
BACKGROUND OF THE INVENTION 
     This invention relates to multimode, circular-symmetric, isotropic, optical fiber waveguides and, more particularly, to optical fiber waveguides in which it is desirable to achieve minimum modal dispersion. 
     Optical fiber waveguides are now well recognized in the art as desirable mediums for transmitting optical information. Initially these optical fibers were constructed with a core having a uniform index of refraction surrounded by a cladding with a lower valued index of refraction. In this type of optical fiber waveguide with a stepped index profile and optical energy is coupled into the core, and the energy is transmitted to the far end of the optical fiber through a process of multiple reflections from the core-cladding interface. 
     One difficulty found with the multimode optical fiber having a stepped index profile is related to the fact that the various modes take widely different transit times. The modes which encounter very few reflections appear at the receiving end of the fiber much sooner than the modes which encounter many reflections. Since the latter modes are caused to travel through a longer length of the medium before reaching the receiving end of the fiber, this multimode effect causes any optical pulse transmitted through the fiber to encounter pulse dispersion. 
     A technique for reducing the effect of this multimode dispersion was disclosed in the article entitled &#34;Multimode Theory of Gradient Core Fibers&#34; by D. Gloge and E. A. J. Marcatili published in the November, 1973 issue of the Bell System Technical Journal pg. 1653-1678. In accordance with this technique, the index of refraction is caused to change along the radius of the fiber. The index of refraction at the core center has the highest value, and the index is changed in a roughly parabolic shape so as to decrease to the value of the index in the cladding at the core-cladding interface. The index profile in this type of fiber waveguide is given by the equation ##EQU4## where n 1  is the on-axis refractive index, n 2  is the refractive index of the cladding and of the fiber core at radius a, ##EQU5## and a is the core radius. For α = 2(1-Δ) the fiber has an almost parabolic index profile, and the modes do not differ in transit time much from each other. 
     In the Gloge-Marcatili analysis, it was assumed that a parameter identified as profile dispersion is negligible. This parameter, profile dispersion, will be identified more completely hereinafter. At this point, it need not only be said that the profile dispersion is a function of the rate of change of the index with respect to wavelength. 
     The above-identified analysis by Gloge and Marcatili was extended in a very important way by D. B. Keck and R. Olshansky to optical fibers, wherein the profile dispersion is constant throughout the radius of the fiber core. See, for example, U.S. Pat. No. 3,904,268 entitled &#34;Optical Waveguide Having Optimal Index Gradient&#34; issued Sept. 9, 1975 to D. B. Keck and R. Olshansky. Where the profile dispersion is constant, it was determined by Keck and Olshansky that the index profile still follows a power law for minimum modal dispersion, but the exponent α in the index profile equation requires a value other than 2(1-Δ). Specifically, the exponent α should satisfy the following equation for optimal modal dispersion: ##EQU6## where ##EQU7## 
     It has recently been determined though that the profile dispersion in graded index profiles is not constant with respect to radius for some of the dopants that are presently being used to shape the index profile. See, for example, the article entitled &#34;Pulse Broadening in Multimode Optical Fiber With Large Δn/n&#34; by J. A. Arnaud and J. W. Fleming to be published in Electronics Letters. In accordance with the mathematical technique presented by Arnaud and Fleming, the RMS impulse response width of an optical fiber can be determined using measured values of dn/dλ. As pointed out in their article, the RMS impulse response width for fibers with germanium dopant do not achieve the optimum impulse response dictated by the theory established by Olshansky and Keck. As is further pointed out as in the Arnaud and Fleming article, this lack of agreement is due primarily to the fact that (1/n )dn/dλ is not a constant as assumed in Olshansky and Keck&#39;s theory. In terms of the analysis to be presented hereinafter (1/n)dn/dλ is assumed to be an arbitrary function of λ and r, and consequently the results of Gloge et al and Olshansky et al are extended to cover a vast class of fibers. 
     SUMMARY OF THE INVENTION 
     It is an object of the present invention to increase the information carrying capacity of an optical fiber communication system by employing a graded index optical waveguide having an optimal index gradient for reducing dispersion even in those fibers which use dopants that have profile dispersions other than constant. Briefly, the invention relates to optical fiber waveguides in which the index of refraction is graded in the core in accordance with the following equation 
     
         n.sup.2 = n.sub.1.sup.2 (1-F) 
    
     where n 1  is the index on the axis and the function F is the profile function of the fiber. This function F is 0 on the axis and is a function of a radial coordinate r and wavelength λ within the core. At the cladding-core interface the function F is equal to 2Δ where ##EQU8## and n 2  is the index of refraction in the cladding. The profile dispersion p in this type of optical waveguide can be expressed in terms of the following equation ##EQU9## where N 1  is the group index on axis. As pointed out hereinafter, the parameter p, the profile dispersion, can be obtained by taking measurements of n and dn/dλ on the various glasses with different concentration of dopants to be used in the construction of the fiber. 
     Using these parameters the profile function (and consequently the index profile) that achieves minimum modal dispersion in a fiber having any arbitrary profile dispersion can be obtained in accordance with the present invention by solving the following general equation ##EQU10## where D is an arbitrary function of λ having a value approximately equal to 2 for minimum modal dispersion. This parameter D, designated hereinafter as the dispersion parameter, is more precisely set equal to 1 + √1 - 2Δ to achieve minimum impulse response width. 
     As mentioned before, by taking measurements of the index of refraction and its derivatives with respect to λ for the doped glasses that are to be used to provide the index profile, the profile dispersion p can be determined for all values of F. Accordingly, the above general equation can then be solved in order to determine the function of F with respect to radius. This determination of F dictates the index profile that is required to provide the minimum impulse response width. 
     In cases where only a single wavelength is to be used and, therefore, D is a constant, the profile dispersion is only a function of the radius r; therefore the partial differential equation is reduced to an ordinary differential equation whose solution yields the optimum index profile in terms of the profile dispersion p and the dispersion parameter D: ##EQU11## where the profile dispersion p is a function of F determined by measurement of n and its derivatives with respect to wavelength. This profile dispersion and the constant D can be inserted into this integral equation and the value of the radius r for any particular value of the index profile F can be determined from the equation. The resulting specification of the radial position for each value of F completely determines the dielectric profile of n versus r. 
     Up to now we have assumed p known (from measurements) and have calculated F via either one of above two equations I or II. 
     In accordance with a second important feature of the present invention, the index profile can be chosen to satisfy another criteria besides that of minimum modal dispersion. This criteria might, for example, be the establishment of an optical fiber with larger numerical aperture. The index profile under these circumstances is known as a function of radius and can, therefore, be utilized in the above general equation (I) to find the profile dispersion which is necessary to provide minimum modal dispersion. As a result of this additional feature, an optical fiber with minimum modal dispersion and a larger numerical aperture than was heretofore available in the prior art can be manufactured. To construct this fiber a maximum amount of index-increasing dopant can be utilized on the axis of the core and this index-increasing dopant can be graded to 0 at a point in the radius between the axis and the core-cladding interface. From this point outward toward the core-cladding interface an index-decreasing dopant can be utilized grading out to a maximum amount of this dopant at the core-cladding interface. This fiber will be shown later to have a discontinuous index profile and the minimum modal dispersion is achieved by providing a profile dispersion which is dictated by the above equation (I). 
    
    
     BRIEF DESCRIPTION OF THE DRAWING 
     FIG. 1 shows a segment of an optical fiber waveguide used to practice the present invention; 
     FIG. 2 is a graph of transit time versus mode parameter for a fiber constructed in accordance with the present invention; 
     FIG. 3 is a graph of the profile dispersion p versus F obtained from measurements on bulk samples having various concentrations of the dopants GeO 2  and B 2  O 3  in a SiO 2  optical fiber; 
     FIG. 4 is a plot of the correct deviation of the profile function F from the F for a square law versus a normalized radius obtained by using the exact theory of the present invention as compared with the simplified theory of Olshansky and Keck which assumes a constant profile dispersion; and 
     FIG. 5A and B are graphs of the dielectric profile and the profile dispersion versus radius respectively for a fiber designed to provide minimum modal dispersion and increased numerical aperture. 
    
    
     DETAILED DESCRIPTION 
     The present invention is useful in optical fiber waveguides of the type illustrated in FIG. 1. In this type of waveguide the core 100 is circular-symmetric, having a cross-sectional radius equal to a and the index of refraction n is graded from the on axis value of n 1  to the value of n 2  at the core-cladding interface 101. With n 1  &gt; n 2  optical energy coupled into the core is confined to the core and this energy is propagated through the optical fiber to a receiving location. 
     In my above-identified article with D. C. Gloge entitled &#34;Multimode Theory of Graded-Core Fibers,&#34; a WKB approximation is utilized to obtain an equation which relates the propagation constant β of a mode characterized by the radial and azimuthal wave numbers μ and ν to the free space propagation constant k, where k = 2π/λ. The refractive index is designated by n, and the radial coordinate is designated by r. This equation is as follows ##EQU12## where 
     
         ρ = √(k.sup.2 n.sup.2 -β.sup.2)r.sup.2 - ν.sup.2. (2) 
    
     In these equations, r 1  and r 2  are the two turning points in the radial coordinate that cause the radical to be equal to zero (ρ=0), and between which points most of the field of the mode is concentrated. In further manipulation of this equation, it is useful to define the index of refraction with the following equation: 
     
         n.sup.2 = n.sub.1.sup.2 (1-F)                              (3) 
    
     and to define the propagation constant β by the following equation 
     
         β.sup.2 = k.sup.2 n.sub.1.sup.2 (1-B)                 (4) 
    
     where F is a profile function which is equal to zero on the axis and is an arbitrary function of r and wavelength λ within the core, while it is equal to 2Δ in the cladding. Similarly, B is a mode parameter which varies between zero for the lowest order mode and is equal to 2Δ for the modes whose phase velocities coincide with that of a plain wave in the cladding. Using these definitions of β and n in the radical representing ρ in equation 2 becomes: 
     
         ρ = √(kn.sub.1 r).sup.2 (B-F) - ν.sup.2      (5) 
    
     the group velocity of a mode can be introduced by taking the derivative with respect to the free space wavelength, λ, of both sides of equation (1). Inasmuch as the radial mode number μ should not be dependent on wavelength, the derivative of μ with respect to wavelength can be set equal to zero thereby resulting in the following equation: ##EQU13## where ##EQU14## N 1  is the group index on the axis, and p is a generalized version of the profile dispersion. A similar profile dispersion is defined and measured in the article entitled &#34;Profile Dispersion in Multi-Mode Fibers: Measurement and Analysis&#34;, by D. Gloge, I. P. Kaminow, H. M. Presby, Electronics Letters, Sept. 18, 1976, Vol. 11, No. 19 pp. 469-471. Since F is both a function of r and λ, the profile dispersion as given in (8) as in general also a function of λ and r. 
     The derivative of the mode parameter B with respect to λ can be expressed in terms of the group delay t of a particular mode by taking the square roots and derivatives of both sides of equation (1). The following equation results: ##EQU15## in which T, the flight time for the ray on the core axis, is related to the delay of a plain wave in a medium having a group index N 1  and a length L by the following equation: ##EQU16## where c is the velocity of light in free space. Substituting the value of dB/dλ obtained from equation (9) into equation (6) yields the following integral: ##EQU17## This integral equation was solved in prior art references by assuming that the profile dispersion, p, is a constant with respect to radius. As pointed out in the above-identified article by J. A. Arnaud and J. W. Fleming, this assumption cannot be made where large concentrations of chemical dopants are utilized in order to achieve a desired index profile. 
     Taking note of the fact that the radical ρ = o at the radial turning points of r 1  and r 2 , the following self-evident expression may be written: ##EQU18## This integral expression becomes very useful if the (δρ/δr)is obtained by operating upon equation (5) and this result is then inserted into equation (12). The following equation then results: ##EQU19## In this latter equation, it should be noted that B, the mode parameter, is not a function of the radial coordinate r, and therefore equation (13) can be rewritten in the following form: ##EQU20## In addition to the mode parameter B, the group delay t of the mode is also independent of the radial coordinate r, as T, the flight time on axis, and therefore equation 11 can be expressed in the following form: ##EQU21## 
     Equation 15 can then be divided by equation 14 to yield the following general expression: ##EQU22## This expression is valid for any circular-symmetric optical fiber with isotropic dielectric and in general, is still difficult to solve exactly. If, however, a particular special family of fibers is considered wherein the index profile and the profile dispersion are related in accordance with the following equation: ##EQU23## where D is an arbitrary function of the wavelength λ, the above general expression (16) can then be solved and reduced to the following form:  In this latter expression, t, the group delay of a particular mode is characterized by the mode parameter B, the flight time of a ray on axis T, and this newly introduced dispersion parameter D. This expression indicates that the group delay of a mode is only a function of the mode parameter B and the dispersion parameter D. Most importantly, this group delay is independent of the mode numbers (which means that modes with the same propagation constant have the same delay) and is also independent of the profile function and of the profile dispersion. 
     The impulse response width can be determined from equation (18) by finding the differences between the flight time of the slowest and fastest modes or rays for any given value of the dispersion parameter D. It can easily be determined that the minimum differences between the fastest and slowest modes occurs if D is chosen in accordance with the following equation: 
     
         D.sub.o = 1 + √1 - 2Δor D.sub.o ≅ 2 - Δ for Δ &lt;&lt; 1                                              (19) 
    
     With this value of D equal to D o  the lowest and highest order modes characterized by B = O and B = 2Δ respectively are the slowest modes and they arrive at the end of the fiber after a flight time equal to T. The modes characterized by B = 1 - √1 - 2Δ are the fastest modes and they arrive at the end of the fiber at a minimum time dictated by the following equation: ##EQU24## A plot of equation (18) for the optimum value of the dispersion parameter D O  is given in FIG. 2. 
     As pointed out hereinabove, the dispersion parameter D is ordinarily a function of wavelength λ. If, however, the operation of the fiber is restricted to a single wavelength or to a narrow band, the dispersion parameter D is a constant value. Under these circumstances the generalized equation (17) can be simplified. The profile dispersion is a function of the radial coordinate r, and the partial derivative of the profile function F is reduced to an ordinary derivative. Under these conditions equation (17) takes on the simplified form of ##EQU25## 
     This latter equation can be used two ways. If the profile dispersion p is known as a function of the profile function F, the equation can be solved to determine the index profile that will give the minimum impulse width. Secondly, the index profile can be specified in an arbitrary manner, and the equation can be solved to determine what profile dispersion versus radius is necessary to obtain minimum impulse width. 
     In accordance with the first usage of the equation, that is where the profile dispersion is known as a function of F or r, a most convenient form of equation (21) is as follows: ##EQU26## To use this latter equation, the profile dispersion p can be measured for a multimode fiber in the manner set forth in the article entitled &#34;Profile Dispersion in Multimode Fibers: Measurement and Analysis&#34; by D. Gloge, I. P. Kaminow and H. M. Presby, Electronics Letters, Sept. 18, 1975, Vol. 11, No. 19, pp. 469-471. The Gloge et al article presents a desirable technique for measuring the index of refraction of optical fibers since it measures that index directly on a fiber after the fiber has undergone its total history of temperature changes. As is well known to those skilled in the art, the index of refraction is a function not only of the wavelength and dopants used but also of the numerous temperature changes that the glassy material has encountered. 
     Still others in the art measure index of refraction and the derivatives of this index by taking measurements on bulk samples. Mr. J. W. Fleming has taken measurements on bulk samples using GeO 2  -B 2  O 2  -SiO 2  glasses and presented these results at the fall meeting of the American Ceramic Society in October, 1975. The results of these measurements are also presented on pp. 418-428 of the text entitled &#34;Beam and Fiber Optics&#34; by J. A. Arnaud, Academic Press, Feb. 1, 1976. In all of these measurements both the index of refraction n (or n 2 ) and the derivative of n 2   with respect to the wavelength λ can be determined for the various dopants that are to be used in an optical fiber. 
     Thus, using the definition of the profile function F given in equation (3), we can easily establish an equation which gives us the values of F for each of the various indices n generated with different quantities of dopant to be utilized. This equation is as follows ##EQU27## A simple derivative of the latter equation with respect to wavelength also provides us with the derivative of the profile function with respect to λ ##EQU28## Using equations (23), (24) and equation (8), which defines the profile dispersion, we can obtain the value of profile dispersion for each and every value of F that can be achieved by the various quantities of dopant to be utilized in the construction of an optical fiber. Hence, measurements of the type already performed in the art can be utilized to determine the dispersion parameter as a function of F. These values can be utilized in equation (22) along with the desired value of the dispersion parameter D, and this equation can then be solved to determine the radial location r for each and every value of the profile function F. In summary, the dispersion parameter p can be obtained by measurement techniques well known in the art and equation (22) can be utilized in order to determine the index profile F that is required in order to obtain minimum modal dispersion. 
     FIGS. 3 and 4 illustrate an application of the above technique in the design of an optical fiber having minimum modal dispersion using germanium oxide (GeO 2 ) and boron oxide dopants (B 2  O 3 ) dopants in a silicon dioxide (SiO 2 ) fiber. As pointed out hereinabove, measurements on bulk samples of silicon dioxide using germanium oxide and boron oxide as dopants have been performed by Fleming and previously reported in the literature. Using the results of these measurements as reported on page 419 of the above-identified text by Arnaud, values of n 2  and the derivative of n 2  with respect to wavelength can be computed for each of the various concentrations of GeO 2  and B 2  O 3 . As is well known to those skilled in the art, germanium oxide increases refractive index and boron oxide decreases the refractive index. Arbitrarily choosing a maximum value of 17 percent molar concentration of germanium oxide as the material to be used on the fiber axis, values of the profile function F can be determined for all of the other molar concentrations to be utilized in the fiber by using equation (23). Using both the determined values of n 2  and dn 2  /dλ the derivative of the profile function F with respect to wavelength can be computed using equation (24). At this point, the profile dispersion p can be computed using equation (8) for each of the molar concentrations to be utilized in the fiber. A plot of the profile dispersion versus the profile function for the selected values of germanium oxide and boron oxide to be used in the fiber is presented in FIG. 3. In FIG. 3 each of the computed points is associated with a two-numeral description. The first numeral indicates the molar concentration of germanium oxide (GeO 2 ) and the second numeral indicates the molar concentration of boron oxide (B 2  O 3 ). 
     Using the values of profile dispersion obtained by computation and plotted in FIG. 3, equation (22) can then be solved to determine the radial location r for each value F of the profile function. For the dopants under consideration, solution of equation (22) yields a profile function F that is illustrated in FIG. 4. In FIG. 4 the deviation of the profile function from a square law profile function is plotted versus (r/a) 2 , normalized radial coordinate. Square law profile function refers to the case where the dielectric profile has a parabolic shape. Horizontal line 400 in FIG. 4 is the plot which would result from a parabolic dielectric profile. The difference between the profile function that is obtained using the exact theory of the present invention and the square law profile is shown in FIG. 4 as the solid line 401. As indicated in FIG. 4, the profile function slightly exceeds the equivalent square law profile function for values of the radial coordinate between 0 and √0.4a. For radial coordinates in excess of this value, the profile function F is less than the square law profile function. Recalling the fundamental definition of F given in equation (3), line 401 in FIG. 4 indicates the dielectric profile should be smaller than a parabolic function up to a radial coordinate of about √0.4a and larger than the same parabolic function beyond this point. 
     Also plotted in FIG. 4 as line 402 is the difference between the profile function obtained using the above-identified Keck-Olshansky theory and the square law profile function. As clearly indicated in FIG. 4, there is a significant difference between the profile function obtained in accordance with the present exact theory and the one which is obtained in accordance with the Keck-Olshansky theory. As further indicated by line 402 in FIG. 4, the Keck-Olshansky theory requires a dielectric profile which for all radial coordinates has a dielectric constant greater than the dielectric constant in a parabolic shape. As pointed out hereinabove, the Keck-Olshansky theory assumes a profile dispersion which is constant with respect to radius and this is clearly not in keeping with the concentration of dopants under consideration. 
     As indicated hereinabove, it is only in those cases where the dispersion parameter p is a constant that the profile function will follow a power law. Setting p equal to a constant value of P o , the integral equation (22) can be solved to provide the following function for the index profile: ##EQU29## where 
     
         α = D(2-P.sub.o) -2.                                 (26) 
    
     This last equation establishes the relation between the dispersion parameter D of the fiber and the alpha value so widely used in the prior art literature for fibers with constant profile dispersion. This value of α can most easily be compared to the value presented in the above-identified patent of Keck and Olshansky by recalling from equation (19) hereinabove that the dispersion parameter D has a value approximately equal to 2 for minimum modal dispersion. 
     As pointed out hereinabove, equation (21) can also be utilized to solve for the profile dispersion which is necessary in order to provide minimum modal dispersion with an arbitrary profile function. This usage can be illustrated in the case of an optical fiber that is designed to provide an increased numerical aperture. As is well known to those skilled in the art, the numerical aperture of a fiber is a function of Δ which in turn can be increased by increasing the difference between the index of refraction at the axis of the core and the index of refraction at the core-cladding interface. 
     To obtain an increased numerical aperture one could use a maximum amount of index-increasing dopant such as germanium oxide on the fiber axis and this dopant could be decreased to a concentration of O at some point in the core between the core axis and the core-cladding interface. At this point in the core, increasing amounts of an index-decreasing dopant could be added up to the core-cladding interface at which point a constant index could be maintained throughout the cladding. 
     As an illustration let us choose the profile function of the fiber shown in FIG. 5A and represented by the following equations: ##EQU30## To achieve this dielectric profile a maximum amount of index increasing dopant is utilized at r = 0 and this dopant is decreased in molar concentration to the point in the core where r = a o . At this point the index-increasing dopant has a molar concentration of O and an index-decreasing dopant is added up to the point where the radial coordinate r = a. At this latter point the index-decreasing dopant has a maximum molar concentration. 
     Substituting equations (27) and (28) into equation (21) and assuming that the dispersion parameter D has its optimum value D o , the required profile dispersion needed to achieve minimum modal dispersion is as follows: ##EQU31## These equations for the profile dispersion are plotted in FIG. 5B. A first constant value of profile dispersion is needed between r = o and r = a o  and a second constant value is required between r = a o  and r = a. If the materials chosen to achieve this profile dispersion have a measured profile dispersion that is other than constant, this deviation from a constant value must be taken into acount by solving equation (21) for the correct measured profile dispersion. 
     Up to this point in determining the index profile which provides minimum modal dispersion we have assumed that the dispersion parameter D is a constant. This assumption is very accurate when the optical fiber is operated over a small range of wavelengths. As pointed out hereinabove, however, the dispersion parameter D is a function of λ and variations in this parameter can result in variations in the impulse width unless additional criteria to be described hereinbelow are satisfied. 
     As pointed out hereinabove in connection with FIG. 2, the optimum dispersion parameter results in a group delay time for the slowest modes equal to T and a group delay time for the fastest modes equal to the value provided by equation (20). If τ represents the time spread between the fastest and slowest modes, the minimum value for τ corresponding to the value of D utilized in FIG. 2 can be expressed by the following equation: ##EQU32## From this equation it can be seen that fibers with the same Δ have the same minimum impulse response width and this minimum impulse response width is independent of the index profile and profile dispersion provided the index profile and profile dispersion are related in accordance with the design equation (19). If Δ &lt;&lt; 1, the minimum impulse response width can be expressed by the following equation: ##EQU33## If the dispersion parameter D is allowed to deviate from its optimum value of D o  by a parameter designated by δ in the following equation: 
     
         D = (1+δ)D.sub.o.                                    (34) 
    
     The ratio of τ/τ min  can be expressed where δ &lt;&lt; 1 and where Δ &lt;&lt; 1 by the following equation: ##EQU34## This latter equation (35) indicates that a small fractional departure of the dispersion parameter from its optimum value results in a substantial change in the impulse response width. In fact, if δ is equal to Δ, then the impulse response width encountered is nine times larger than the minimum impulse response width. 
     In view of this sensitivity of the impulse response width to variations in the dispersion parameter, fibers that are to be operated over large bands of wavelengths require additional design criteria to increase the band over which a fiber operates with low modal dispersion. In fact, it is desirable to have the dispersion parameter equal to its optimum value provided in equation (19) and in addition to have the derivative of the dispersion parameter with respect to wavelengths equal to 0. This latter criteria can be met by satisfying the following equation: ##EQU35## which is derived from equation (17) by taking the derivative of both sides with respect to λ and setting the derivative of D with respect to λ equal to 0. Inasmuch as this last equation involves a derivative of the profile dispersion with respect to λ, this last equation requires control over the second derivative F with respect to λ. If this second derivative of the profile function is controlled by the designer and both of the equations (17) and (36) are satisfied, minimum modal dispersion can be controlled over substantially broader frequency bands. 
     It can easily be extrapolated that control over higher derivatives permits even further demands on D. In fact, if all the higher derivatives are controllable, D(λ) can be chosen arbitrarily and the profile F is the solution of equation (17) subject to the conditions of being zero at r = 0 and 2Δ(λ) at r = a.