Tania function

Explicit plot of $y=\mathrm {Tania} (x)$

Tania function is the particular solution f = Tania of the equation

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (1)~~~f'(z)={\frac {f(z)}{1+f(z)}}$

satisfying the initial condition

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (2)~~~f(0)=1\ .$

Relation with other special functions

f= ArcTania(x + i y) in the x,y-plane with levels u = ℜ(f) = constant, and v = ℑ(f) = constant.

f= Tania(x + i y) in the x,y-plane with levels u = ℜ(f) = constant, and v = ℑ(f) = constant.

The inverse function of Tania can be called "ArcTania"; ArcTania = Tania⁻¹. The ArcTania is an elementary function:

$\mathrm {(3)} \ \ \mathrm {ArcTania} (z)=z+\ln(z)-1$

The complex maps of f = ArcTania(z) and f = Tania(z) are shown in figures at right in the plane x = ℜ(z), y = ℑ(z) with lines u = ℜ(f) = constant, and lines v = ℑ(f) = constant. Symbol ℜ denotes "real part" and symbol ℑ denotes "imaginary part".

At least for |ℑ(z)| < π, the relation below holds

$\mathrm {ArcTania} (\mathrm {Tania} (z))=z$

In certain range of values of the argument, the Tania function can be expressed through the Lambert W-function (known also as the ProductLog function) as follows:

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (4)~~~\mathrm {Tania} (z)=\mathrm {LambertW} (\exp(1\!+\!z))~$, $~~|\Im (z)|<\pi$

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (5)~~~\mathrm {LambertW} (z)=\mathrm {Tania} {\big (}\ln(z)-1{\big )}~$

At least in vicinity of the real axis, the Tania function can be expressed through the Lambert W-function; while the principal branch of the Lambert W-function can be expressed through the Tania function for all complex values of the argument.

After to use, implement and describe the Tania function, it was revealed, that very similar function, WrightOmega is already described :

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (6)~~~\mathrm {WrightOmega} (z)=\mathrm {Tania} (z\!-\!1)$

History

The Tania function is used in the Laser science for the description of intensity of light in the idealized gain medium.

The Tania function appears as the intensity of a wave in the uniformly pumped gain medium, measured in unit of the saturation intensity. Then, its argument has sense of length of propagation, measured in units of the inverse of unsaturated gain. Such a model describes the idealized laser amplifier with simple kinetics, and can be applied, for example, to the Yb-doped ceramics pumped in vicinity of wavelength of order of 940nm; the maximum gain corresponds to wavelength in vicinity of 1030 nm. For this wavelength, the saturation intensity may have value of order of $I_{\rm {sat}}=100{\rm {W/mm}}^{2}$, and the unsaturated gain (that depends on the concentration of dopant and the intensity of pump) may be of order of $g_{0}\approx 1/\mathrm {mm}$. This determines the orders of magnitude of values of the physical parameters, while the intensity along the propagation at distance $x$ can be described with intensity $I(x)$ satisfying the equation

$\displaystyle (7)~~~I'(x)={\frac {g_{0}I(x)}{1+I(x)/I_{\rm {sat}}}}$

The solution $I$ can be expressed as

$\displaystyle (8)~~~I(x)=I_{\rm {sat}}\mathrm {Tania} (g_{0}x-g_{0}x_{1})$

where $x_{1}$ is distance at which the intensity saturates the gain, id est, $I(x_{1})=I_{\rm {sat}}$.

The phenomenology of waves with such a saturated gain has been considered in century 20; some results are published in the Soviet Journal of Quantum Electronics.

The function "Tania" is named after the short name of T.I. Kuznetsova.

For the applications, the simple and robust implementation of the Tania function is an important tools of the scientific research. For this reason, the description of the Tania function is loaded at TORI (Tools for Outstanding Research and Investigations). Only then the relation with WrightOmega had been revealed. The WrightOmega function is implemented in many programming languages including Mathematica.

Properties of Tania and ArcTania

The $\mathrm {Tania} (z)$ has two cut lines, $\Im (z)=\pm \pi ~$, $~\Re (z)\leq -2$. In the figure, these lines are marked with dashed half-lines. The branch points are determined by the singularity of the expression in the right hand side of equation (1), they correspond to values of ArcTania equal to $-1=\exp(\pm \mathrm {i} \pi )$. Using (3), these points can be expressed as follows:

$\!\!\!\!\!\!\!\!\!(9)~~~z_{\pm }=\exp(\pm \mathrm {i} \pi )+\ln(\exp(\pm \mathrm {i} \pi ))-1=-2\pm \mathrm {i} \pi$.

In the left hand side of the figure, between the cut lines lines, the function exponentially approaches zero. In other directions, the modulus of the function is smoothly growing; the growth is similar to that of the linear function.

Tania is real–holomorhic;

$\!\!\!\!\!\!\!\!\!(10)~~~\mathrm {Tania} (z^{*})=\mathrm {Tania} (z)^{*}$

The expansions of Tania in diffrent cases, and, in particular, at the branchpoints $z_{\pm }$ is described below.

Expansion of Tania at infinity

The complex maps of Tania and ArcTania above show that asymptotically each of these functions becomes similar to its argument. Being far from the branch points (and outside the cut lines), the grid of the lines of constant real and constant imaginary part looks similar to that of the identity function. This indicates that the leading term of the expansion at infinity of $\mathrm {Tania} (z)$ should be just $z$. Such an expansion is suggested in this section.

At large values of the argument, Tania shows almost linear growth and can be expanded as follows:

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (11)~~~\mathrm {Tania} (z)=$ $z+1\!-\!\ln(z)+{\frac {\ln(z)\!-\!1}{z}}$ $+{\frac {\ln(z)^{2}-4\ln(z)+3}{2z^{2}}}$ $+{\frac {2\ln(z)^{2}-15\ln(z)+30\ln(z)-17}{6z^{3}}}$ $+O{\Big (}{\frac {\ln(z)}{z}}{\Big )}^{4}$

At $|z|\gg 1$, the truncated series (7) can be used for the evaluation of $\mathrm {Tania} (z)$, while the argument is not between the cut lines in figure. For the evaluation, the asymptotic can be written also as follows:

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!(12)~~~\mathrm {Tania} (z)=$ $(z\!+\!1)\!-\!\ln(z\!+\!1)+{\frac {\ln(z\!+\!1)}{z+1}}$ $+{\Big (}{\frac {\ln(z\!+\!1)}{z+1}}{\Big )}^{\!2}{\Big (}{\frac {1}{2}}-\ln(z\!+\!1)^{-1}{\Big )}+{\Big (}{\frac {\ln(z\!+\!1)}{z+1}}{\Big )}^{\!3}{\Big (}{\frac {1}{3}}-{\frac {3}{4}}\ln(z\!+\!1)^{-1}+\ln(z\!+\!1)^{-2}{\Big )}$ $+{\Big (}{\frac {\ln(z\!+\!1)}{z+1}}{\Big )}^{\!4}{\Big (}{\frac {1}{4}}-{\frac {11}{\!5}}\ln(z\!+\!1)^{-1}+3\ln(z\!+\!1)^{-2}-\ln(z\!+\!1)^{-3}{\Big )}$ $+O{\Big (}{\frac {\ln(z\!+\!1)}{z+1}}{\Big )}^{\!5}$

Asymptotic expansion of Tania between the cut lines

The expansions (11) and (12) are not valid between the cut lines. For this range, for large values of the argument, the special expansion is suggested in this section.

For $-\Re (z)\!\gg \!1$ at $|\Im (z)|\!<\!\pi$

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (13)~~~\mathrm {Tania} (z)~\sim ~$ $\varepsilon -\varepsilon ^{2}+{\frac {3}{2}}\varepsilon ^{3}-{\frac {7}{2}}\varepsilon ^{4}+O(\varepsilon ^{5})$

where $\varepsilon =\exp(1\!+\!z)$.

Such expansion can be obtained iterating assignment

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (14)~~~f~=~\varepsilon ~\exp(-f)$

with the initial approach $f=\varepsilon +O(\varepsilon ^{2})$.

Together with the expansion (11), the expansion (13) covers the most of the complex plane, leaving only a finite-size domain that includes the branch points and the origin of the coordinates. The expansion at these points are considered below.

Expansion of Tania at singularities

$\mathrm {Tania} (z)$ has two cut lines at $\Re (z)\!\leq \!-2$, $\Im (z)\!=\!\pm \pi$. The expansions above do not reproduce the behavior of Tania in vicinity of the branch points. In this section, the expansion at these branch points is suggested.

For the expansion of $\mathrm {Tania} (z)$ at the upper branch point $z=-2+\pi \mathrm {i}$, the convenient small parameter is

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!(11)~~~t=\mathrm {i} {\sqrt {{\frac {2}{9}}(z+2-\pi \mathrm {i} )~}}$

then, the expansion of Tania can be written as follows:

$\displaystyle \!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!(15)~~~\mathrm {Tania} (z)=$ $-1+3t-3t^{2}$ $+{\frac {3}{4}}t^{3}$ $+{\frac {3}{10}}t^{4}$ $+{\frac {9}{160}}t^{5}$ $-{\frac {3}{70}}t^{6}$ $-{\frac {1251}{22400}}t^{7}$ $-{\frac {9}{280}}t^{8}+...$

The Mathematica Series routine easy calculates some tens of the first coefficients of the expansion. The series seems to converge while $|z\!+2\!-\!\pi \mathrm {i} |\!<\!2\pi$; in vicinity of another branch point, id est, at $z\approx -2\!-\!\mathrm {i}$, such an approximation is not valid.

Similar expansion of Tania at point $-2-\pi \mathrm {i}$ can be obtained by complex conjugation of expressions (9) and (10).

Taylor expansions of Tania

The expansions above, in principle, allow to cover all the complex plane with the approximations of the Tania function. But the expansion at the branch points may be not convenient while dealing with Tania of a real argument. For this case, the Taylor expansions of Tania at the real axis are suggested in this section.

The Taylor expansion of Tania at zero

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (16)~~~\mathrm {Tania} (z)=$ $1$ $+{\frac {z}{2}}$ $+{\frac {z^{2}}{16}}$ $-{\frac {z^{3}}{192}}$ $-{\frac {z^{4}}{3072}}$ $+{\frac {13z^{5}}{61440}}$ $-{\frac {47z^{6}}{1474560}}$ $+{\frac {73z^{7}}{41287680}}$ $+...$

converges at $|z|<{\sqrt {4+\pi ^{2}}}\approx 3.724$; this series can be obtained with the Mathematica's operator InverseSeries, converting the expansion

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (17)~~~\mathrm {ArcTania} (z)=$ $2(z\!-\!1)$ $-{\frac {(z\!-\!1)^{2}}{2}}$ $+{\frac {(z\!-\!1)^{3}}{3}}$ $-{\frac {(z\!-\!1)^{4}}{4}}$ $+{\frac {(z\!-\!1)^{5}}{5}}$ $+{\frac {(z\!-\!1)^{6}}{6}}$ $+~{O}(z\!-\!1)^{7}$

However, more terms can be added to the series in (16) and (17). At moderate values of $|z|\lesssim 1$, the exansion (11) can be used for the efficient (quick and precise) evaluation of $\mathrm {Tania} (z)$; in order to get 14 significant decimal digits, it is sufficient to take 20 terms.

In the similar way, the expansion of $\mathrm {Tania} (z)$ at $z=-2+1/\mathrm {e} \approx -1.632$ can be written as follows:

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\displaystyle (18)~~~\mathrm {Tania} (z)=$ ${\frac {1}{\mathrm {e} }}$ $+{\frac {1/\mathrm {e} }{1+1/\mathrm {e} }}$ ${\big (}z+2-{\frac {1}{\mathrm {e} }}{\big )}$ $+{\frac {1/\mathrm {e} }{2(1+1/\mathrm {e} )^{3}}}{\big (}z+2-{\frac {1}{\mathrm {e} }}{\big )}^{2}$ $+{\frac {(1-2/\mathrm {e} )/\mathrm {e} }{6(1+1/\mathrm {e} )^{5}}}{\big (}z+2-{\frac {1}{\mathrm {e} }}{\big )}^{3}+...$

The series converges at ${\big |}z+2-{\frac {1}{\mathrm {e} }}{\big |}$ $<$ ${\sqrt {\pi ^{2}+1/\mathrm {e} ^{2}}}\approx 3.1445~$.

The truncated expansions above cover all the complex plane with at least 3–digit approximations.

Evaluation of Tania function

The series above can be used for the evaluation of $\mathrm {Tania} (z)$, dependently on value of $z$. The increase of the precision of the evaluation is easier to achieve not with the increase of number of terms in the expansion, but applying the Newton method. Let $s_{0}$ be initial approximation, obtained with the appropriate formulas from the set above; let

$\displaystyle \!\!\!\!\!\!\!\!\!\!(19)~~~s_{n+1}=s_{n}+{\frac {z-\mathrm {ArcTania} (s_{n})}{\mathrm {ArcTania} '(s_{n})}}$

where

$\displaystyle \!\!\!\!\!\!\!\!\!\!(20)~~~\mathrm {ArcTania} '(z)=1+1/z$

is derivative of the ArcTania function by equation (3). The primary approximations considered in the previous sections allow to get good initial approximation for any point of the complex plane; then the sequence $s$ converges to $\mathrm {Tania} (z)$.

The sequence $s$ converges faster than exponentially. Practically, of order of 5 terms in some of the primary approximations are sufficient, and then the third iteration $s_{3}$ approximates $\mathrm {Tania} (z)$ with 15 decimal digits. Few tens operations are sufficient to get the maximal precision available with the complex(double) variables; the graphics can be plotted in the real time with the C++ generators supplied. (Click the figure to see the implementation.)

Conclusions

The Tania function came from the Laser science describing evolution of the light intensity in the gain medium with simple kinetics. The inverse of Tania function is elementary function. The Tania function can be expressed through the function WrightOmega , but the C++ complex(double) implementation Tania.cin is already prepared; it is used for generators of the figures. The Lambert W function (known also as ProductLog function) can be expressed through the Tania function as well as through the WrightOmega function.

Note

Some of the content in this article previously has appeared on the TORI (Tools for Outstanding Research and Investigation) wiki as Tania function. In 2013 February-March, TORI had been attacked and moved from http://tori.ils.uec.ac.jp/TORI to http://mizugadro.mydns.jp/t