Winding Numbers, Unwinding Numbers, and the Lambert W Function

The unwinding number of a complex number was introduced to process automatic computations involving complex numbers and multi-valued complex functions, and has been successfully applied to computations involving branches of the Lambert W function. In this partly expository note we discuss the unwinding number from a purely topological perspective, and link it to the classical winding number of a curve in the complex plane. We also use the unwinding number to give a representation of the branches Wk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_k$$\end{document} of the Lambert W function as a line integral.


Introduction
The function E : z → z exp z is holomorphic throughout the complex plane, E (z) = 0 if and only if z = −1, and E(z) = 0 if and only if z = 0. However, for any non-zero complex number a, the equation E(z) = a has infinitely many solutions, and the (multi-valued) inverse W of E is known as the Lambert W function. We refer the reader to [3][4][5] for the basic properties of W , and many examples of its use in a variety of different problems. There is a standard construction, and labelling, of a particular set . . . , W −1 , W 0 , W 1 , . . . of branches of W which are defined in [4,5], and illustrated below, and which we now describe. First, C is the complex plane, and Fig. 1 The codomains of the branches W k so that C is the complex plane cut along the interval (−∞, 0], and E is the complex plane cut along the interval (−∞, −1/e]. The branch W 0 is a conformal map of E onto the region labelled W 0 (E ) in Fig. 1 and, for each non-zero integer k, the branch W k is a conformal map of C onto the region labelled W k (C ) in Fig. 1.
Except for a representation of W 0 as an infinite series that is valid throughout E (see [3]), no explict formula for W k is known, and the purpose of [5] was to present the following result which can be used in computer algebra systems to make computations involving the various branches W k . As usual, ln(z) is the principal branch of the complex logarithm defined for all non-zero z by ln(z) = ln |z| + iθ(z), where here, ln |z| is the real logarithm, and θ (z) is the unique choice of the argument of z in the interval (−π, π].
Theorem 1 [5] Let W k be the branches of the Lambert function as defined in [4,5]. Then if k = −1 and z ∈ [−1/e, 0); ln(z) + 2π ki otherwise. (1.1) The idea of the unwinding number of a complex variable was introduced in [5,6] in order to accommodate computations involving complex numbers in various computer algebra systems. The discussions there focussed on the problem of computing multi-valued functions such as the complex logarithm and the inverse trigonometric functions and, in particular, on providing automatic computations near the discontinuities that occur on the branch cuts of these functions. In particular, the authors of [5,6] used the idea of the unwinding number in their proof of Theorem 1. In Sect. 2 we consider Theorem 1 from the perspective of conformal mapping, and we provide a self-contained, topological, proof which does not use the unwinding number.
Although the reason for the choice of the name unwinding number is hinted at in [6], its connection to the fundamental and well established notion of the winding number of a curve in the complex plane seems not to have been described explicitly in the literature, and it is even claimed in one publication that the concept of the unwinding number has no connection with the winding number of a curve. This is plainly false (the unwinding number actually is the winding number of a certain curve); indeed, any complete treatment of multi-valued functions, holomorphic functions, Cauchy's theorem, and so on, necessarily, and inevitably, involves the topological notion of the winding number of a curve in one (perhaps hidden) form or another, and the entire subject of complex analysis can even be developed topologically from the single concept of the winding number of a curve without ever mentioning an integral [1,2]. In Sect. 3 we give a brief description of the unwinding number from a topological perspective.

A Proof of Theorem 1
We present a topological proof of Theorem 1 (which is different from the proof in [6]); then we consider the result from the perspective of conformal maps.
Proof Consider a branch W k , where k is a non-zero integer. Since the functions z and W k (z) are continuous and non-zero in the simply connected region C , it is a topological fact that the functions ln(z) and ln W k (z) are defined, single-valued, and continuous, on in C , and it follows from this that each of the two functions ln(z) − ln W k (z) and W k (z) are a single-valued, continuous, choice of log exp W k (z) in C . It follows that there is some integer p such that As E(2πik) = 2πik, we see that 2πik is a fixed point of W k , and this shows that p = −k. This completes the proof for k = 0. The proof for W 0 follows in a similar way because z/W 0 (z) (which has a removable singularity at z = 0) is non-zero throughout E. Finally, it is sufficient to establish (1.1) with k = 0 when z ∈ C , and with k = 0 when z ∈ E, for the values of W k on the boundary of its domain are uniquely determined by continuity from within its domain. Fig. 2 The conformal map F : z → z + ln(z) Fig. 3 The conformal map of W 2 (C ) onto the strip = {x + iy : 3π < y < 5π } Theorem 1 can be viewed entirely in terms of conformal maps, and this provides a more visual way of thinking about the different branches W k . To begin, the left-hand side of (1.1) is the map z → F W k (z) , where F(z) = z + ln z, and the codomain of W k lies in the domain of F. Thus, to understand (1.1) we need to know about F, and it is known that F maps C onto the complex plane cut along the two half-lines given by {x + iy : x ≤ −1, y = π } and {x + iy : x ≤ −1, y = −π } (see [8] and Fig. 2) For brevity, we shall only consider W 2 , and the following result (illustrated in Fig. 3) is simply a reformulation of Theorem 1.
The point of Theorem 2 is this: from the perspective of complex analysis, there is essentially no difference between the cut plane C and a (conformally equivalent) infinite strip ; thus it is more natural to study a conformal map between a horizontal strip and W 2 (C ) (than a conformal map between C and W 2 (C )) because and W 2 (C ) have a similar geometric shape. Theorem 2 shows that the map from W 2 (C ) to the given strip is a perturbation of the identity map by an amount ln(z) which, for z in W 2 (C ), has imaginary part between −π and 0. Thus this gives a visual explanation of the fact that W 2 (C ) approximates the strip {2π < y < 4π } at its left-hand end, and the strip {3π < y < 5π } at its right-hand end.

The Unwinding Number
In order to create a single-valued choice of the logarithm it is necessary to introduce a branch cut, but this has undesirable consequences for computations since, in general, familiar algebraic identities will fail near these branch cuts. For example, in many cases ln(z 1 z 2 ) = ln z 1 + ln z 2 , but in some cases (for example, when z 1 = z 2 = −1) we have ln(z 1 z 2 ) = ln(z 1 ) + ln(z 2 ). The unwinding number K (z) of a complex number z was introduced in [5,6] in order to resolve these difficulties in a way that makes it possible to engage automatic computing of multi-valued complex functions in various computer algebra systems. Indeed, with the unwinding number K (z) available, we have, for example, for all non-zero z 1 and z 2 . The discussions in [5,6] focussed on the problem of providing automatic computations near the discontinuities that occur on the branch cuts of these functions. In particular, the authors used their idea of the unwinding to give a proof of Theorem 1. Here, we give a brief description of the unwinding number from a topological perspective. We now consider, with p an integer, the three regions = C\(−∞, 0]; S = {x + iy : − π < y < π}; in C which are illustrated in Fig. 4, and three conformal maps between these regions, namely where E(z) = exp z, L(z) = ln(z) and T (z) = z + 2πi p. Trivially, E T L(z) = z on , where we denote the composition of functions by juxtaposition so that f g(z) = f g(z) . It follows that T L E(z) = z on S p , so that According to [6], this means that − p = K (z), so that which is the definition of K (z) in [6]. Note that if we put z = ln z 1 + ln z 2 then (3.1) follows immediately. Now suppose z ∈ S p , and let γ be the (vertical) straight line segment from z −2πi p in S to z in S p , parametrized by γ (t) = (z − 2πi p) + 2πi pt, where t ∈ [0, 1]. As the image of γ under the exponential map is a circle with centre 0 that is traversed It follows that various formulae involving the unwinding number (such as those found in [4][5][6]) can be considered to be algebraic descriptions of topological properties of curves in the plane.

An Integral Formula for the Function W 2 (z)
We end with an expression for W 2 (z) as a line integral, and for this we shall need the following (known) general result.

Theorem 3 Let F be a bijective conformal map of a simply connected region D onto a simple connected region D . Let γ be a simple closed curve in D, and F(γ ) its image in D .
Suppose that z 0 is inside γ , and w 0 = F(z 0 ); then w 0 is inside , and The proof is by the Residue theorem, for z 0 is the only singularity of the integrand that is inside γ , and this is a simple pole with residue R, where
We now consider Fig. 5 below: this represents the top part of Fig. 3 in which we have drawn z 0 , w 0 , γ and F(γ ) (and where we are taking D and D in Theorem 3 to be W 2 (C ) and , respectively). According to Theorem 3, we now find that, for z in C , we have which is the promised formula for W 2 (z) as a line integral. If we take γ to be a small circle of radius r centred at the fixed point 4πi of W 2 , so that ζ = 4πi + re iθ , then we obtain the formula ire iθ dθ. and it would be interesting to see whether this can be derived from (4.1).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.