Thanks to visit codestin.com
Credit goes to dlmf.nist.gov

About the Project
4 Elementary FunctionsLogarithm, Exponential, Powers

§4.13 Lambert W-Function

The Lambert W-function W(z) is the solution of the equation

4.13.1 WeW=z.

On the z-interval [0,) there is one real solution, and it is nonnegative and increasing. On the z-interval (e1,0) there are two real solutions, one increasing and the other decreasing. We call the increasing solution for which W(z)W(e1)=1 the principal branch and denote it by W0(z). See Figure 4.13.1.

See accompanying text
Figure 4.13.1: Branches W0(x), W±1(x0i) of the Lambert W-function. Magnify

The decreasing solution can be identified as W±1(x0i). Other solutions of (4.13.1) are other branches of W(z). They are denoted by Wk(z), k, and have the property

4.13.1_1 Wk(z)=lnk(z)ln(lnk(z))+o(1),
|z|,

where lnk(z)=ln(z)+2πik. W0(z) is a single-valued analytic function on (,e1], real-valued when z>e1, and has a square root branch point at z=e1. See (4.13.6) and (4.13.9_1). The other branches Wk(z) are single-valued analytic functions on (,0], have a logarithmic branch point at z=0, and, in the case k=±1, have a square root branch point at z=e10i respectively. See Figure 4.13.2.

See accompanying text
Figure 4.13.2: The W(z) function on the first 5 Riemann sheets. W(z) maps the first Riemann sheet |ph(z+e1)|<π in the middle of the left-hand side to the region enclosed by the green curve on the right-hand side; it maps the Riemann sheet π<phz<3π on the left-hand side to the region enclosed by the pink, green and orange curves on the right-hand side, etc. Magnify

Alternative notations are Wp(x) for W0(x), Wm(x) for W1(x+0i), both previously used in this section, the Wright ω-function ω(z)=W(ez), which is single-valued, satisfies

4.13.1_2 ω(z)+ln(ω(z))=z,

and has several advantages over the Lambert W-function (see Lawrence et al. (2012)), and the tree T-function T(z)=W(z), which is a solution of

4.13.1_3 TeT=z.

Properties include:

4.13.2 W0(e1) =W±1(e10i)=1,
W0(0) =0,
W0(e) =1.
4.13.3 Moved to (4.13.1_2).
4.13.3_1 W0(xex)={x,1x,(no simpler form),x<1.
4.13.3_2 W±1(xex0i)={(no simpler form),1x,x,x<1.
4.13.4 dWdz=eW1+W=Wz(1+W).
4.13.4_1 dnWdzn=enWpn1(W)(1+W)2n1,
n=1,2,3,,

in which the pn(x) are polynomials of degree n with

4.13.4_2 p0(x) =1,
pn(x) =(1+x)pn1(x)+(1n(x+3))pn1(x),
n=1,2,3,.

Explicit representations for the pn(x) are given in Kalugin and Jeffrey (2011).

4.13.5 W0(z)=n=1(n)n1n!zn,
|z|<e1.
4.13.5_1 (W0(z)z)a=eaW0(z)=n=0a(n+a)n1n!(z)n,
|z|<e1, a.
4.13.5_2 11+W0(z)=n=0nnn!zn,
|z|<e1.
4.13.5_3 (1+W0(z))2=12n=1nn2n!(z)n,
|z|<e1.
4.13.6 W(e1(t2/2))=n=0(1)n1cntn,
|t|<2π,

where t0 for W0, t0 for W±1 on the relevant branch cuts,

4.13.7 c0=1,c1=1,c2=13,c3=136,c4=1270,
4.13.8 cn=cn1n+112k=2n1ckcn+1k,
n=2,3,4,,

and

4.13.9 135(2n+1)c2n+1=gn,

where gn is defined in §5.11(i). See Jeffrey and Murdoch (2017) for an explicit representation for the cn in terms of associated Stirling numbers.

4.13.9_1 W0(z)=n=0dn(ez+1)n/2,
|ez+1|<1, |ph(z+e1)|<π,

where

4.13.9_2 d0 =1,d1=2,d2=23,d3=11362,d4=43135,
(n+2)d1dn+1 =2dn+n2k=1n1dkdnkn+22k=1n1dk+1dnk+1,
n=1,2,3,.

For the definition of Stirling cycle numbers of the first kind [nk] see (26.13.3). As |z|

4.13.10 Wk(z)ξklnξk+n=1(1)nξknm=1n[nnm+1](lnξk)mm!,

where ξk=ln(z)+2πik. For large enough |z| the series on the right-hand side of (4.13.10) is absolutely convergent to its left-hand side. In the case of k=0 and real z the series converges for ze. As x0

4.13.11 W±1(x0i)ηlnη+n=11ηnm=1n[nnm+1](lnη)mm!,

where η=ln(1/x). For these results and other asymptotic expansions see Corless et al. (1997).

For integrals of W(z) use the substitution w=W(z), z=wew and dz=(w+1)ewdw. Examples are

4.13.12 W(z)dz=zW(z)+zW(z)z,
4.13.13 W(z)zdz=12W(z)2+W(z),
4.13.14 2sin(W(z))dz=z(1+1W(z))sin(W(z))zcos(W(z)).
4.13.15 W0(z)=zπ0π(1tcott)2+t2z+tetcottcsctdt.
4.13.16 W0(z)=1π0πln(1+zsinttetcott)dt.

For these and other integral representations of the Lambert W-function see Kheyfits (2004), Kalugin et al. (2012) and Mező (2020).

For the foregoing results and further information see Borwein and Corless (1999), Corless et al. (1996), de Bruijn (1961, pp. 25–28), Olver (1997b, pp. 12–13), and Siewert and Burniston (1973).

For a generalization of the Lambert W-function connected to the three-body problem see Scott et al. (2006), Scott et al. (2013) and Scott et al. (2014).