Grunsky matrix

In mathematics, the Grunsky matrices, or Grunsky operators, are matrices introduced by Grunsky (1939) in complex analysis and geometric function theory. They correspond to either a single holomorphic function on the unit disk or a pair of holomorphic functions on the unit disk and its complement. The Grunsky inequalities express boundedness properties of these matrices, which in general are contraction operators or in important special cases unitary operators. As Grunsky showed, these inequalities hold if and only if the holomorphic function is univalent. The inequalities are equivalent to the inequalities of Goluzin, discovered in 1947. Roughly speaking, the Grunsky inequalities give information on the coefficients of the logarithm of a univalent function; later generalizations by Milin, starting from the Lebedev–Milin inequality, succeeded in exponentiating the inequalities to obtain inequalities for the coefficients of the univalent function itself. Historically the inequalities were used in proving special cases of the Bieberbach conjecture up to the sixth coefficient; the exponentiated inequalities of Milin were used by de Branges in the final solution. The Grunsky operators and their Fredholm determinants are related to spectral properties of bounded domains in the complex plane. The operators have further applications in conformal mapping, Teichmüller theory and conformal field theory.

Contents

Grunsky matrix

If f(z) is a holomorphic univalent function on the unit disk, normalized so that f(0) = 0 and f'(0) = 1, the function

 g(z) = f(z^{-1})^{-1}

is a non-vanishing univalent function on |z| > 1 having a simple pole at ∞ with residue 1:

 g(z) = z %2B b_0 %2B b_1z^{-1} %2B b_2 z^{-2} %2B \cdots

The same inversion formula applied to g gives back f and establishes a one-one correspondence between these two classes of function.

The Grunsky matrix of (cnm) of g is defined by the equation

\log\frac{g(z)-g(\zeta)}{z-\zeta} = -\sum_{m,n>0}c_{nm}z^{-m}\zeta^{-n}

It is a symmetric matrix. Its entries are called the Grunsky coefficients of g.

Note that

\log{g(z^{-1})-g(\zeta^{-1})\over z^{-1}-\zeta^{-1}} = \log{f(z)-f(\zeta)\over z -\zeta} -\log {f(z)\over z} -\log {f(\zeta)\over \zeta},

so that that the coefficients can be expressed directly in terms of f. Indeed if

\log{f(z)-f(\zeta)\over z -\zeta} = -\sum_{m,n\ge 0} d_{mn} z^n \zeta^n,

then for m, n > 0

 \displaystyle{d_{mn} = c_{mn}}

and d0n = dn0 is given by

 \log {f(z)\over z} = \sum_{n>0} d_{0n} z^n

with

 \displaystyle{d_{00}=0.}

Grunsky inequalities

If f is a holomorphic function on the unit disk with Grunsky matrix (cnm), the Grunsky inequalities state that

 \left|\sum_{1\le m,n \le N} c_{mn}\lambda_m \lambda_n \right|\le \sum_{1\le n\le N} |\lambda_n|^2/n

for any finite sequence of complex numbers λ1, ..., λN.

Faber polynomials

The Grunsky coefficients of a normalized univalent function in |z| > 1

g(z)=z%2Bb_0 %2B b_1z^{-1} %2B b_2 z^{-2} %2B \cdots

are polynomials in the coefficients bi which can be computed recursively in terms of the Faber polynomials Φn, a monic polynomial of degree n depending on g.

Taking the derivative in z of the defining relation of the Grunsky coefficients and multiplying by z gives

 {z g^\prime(z)\over g(z) -g(\zeta)} - {z\over z-\zeta} = \sum_{m,n>0} m c_{mn} z^{-m} \zeta^{-n}.

The Faber polynomials are defined by the relation

 {z g^\prime(z)\over g(z) -w}=\sum_{n\ge 0} \Phi_n(w) z^{-n}.

Dividing this relation by z and integrating between z and ∞ gives

 \log {g(z) - w\over z} =-\sum_{n\ge 1}{1\over n}\Phi_n(w) z^{-n}.

This gives the recurrence relations for n > 0

\Phi_n(w)= (w-b_0) \Phi_{n-1}(w) -nb_n -\sum_{0\le i \le n-1} b_{n-i} \Phi_{i}(w)

with

\Phi_0(w)\equiv 1.

Thus

\sum_{n\ge 0} \Phi_n(g(z))\zeta^{-n} =1%2B\sum_{n\ge 1} \left(z^n %2B\sum_{m\ge 1} c_{nm}z^{-m}\right)\zeta^{-n},

so that for n ≥ 1

\Phi_n(g(z))=z^n %2B \sum_{m\ge 1} c_{nm}z^{-m}.

The latter property uniquely determines the Faber polynomial of g.

Milin's area theorem

Let g(z) be a univalent function on |z| > 1 normalized so that

 g(z) = z %2B b_1 z^{-1} %2B b_2 z^{-2} %2B \cdots

and let f(z) be a non-constant holomorphic function on C.

If

 f(g(z))=\sum_{-\infty}^\infty c_n z^n

is the Laurent expansion on z > 1, then

\sum_{n>0}  n|c_n|^2 \le \sum_{n>0} n|c_{-n}|^2.

Proof

If Ω is a bounded open region with smooth boundary ∂Ω and h is a differentiable function on Ω extending to a continuous function on the closure, then, by Stokes theorem applied to the differential 1-form ω = h(z)dz,

 \int_{\partial \Omega} h(z) \, dz = \int_{\partial \Omega} \omega = \iint_\Omega d\omega  = \iint_\Omega (i\partial_x-\partial_y)h \, dx dy= 2i \iint_\Omega \partial_{\overline{z}} h \,dx dy.

For r > 1, let Ωr be the complement of the image of |z|> r under g(z), a bounded domain. Then, by the above identity with h = f', the area of fr) is given by

 A(r)=\iint_{\Omega_r} |f'(z)|^2\, dx dy = {1\over 2i}\int_{\partial\Omega_r} \overline{f(z)}f^\prime(z) \,dz={1\over 2i}\int_{|w|=r} \overline{f(g(w)))} f^\prime(g(w)) g^\prime(w)\, dw.

Hence

\displaystyle{A(r)=\pi \sum_n n|c_{-n}|^2 r^{2n}.}

Since the area is non-negative

 \sum_{n>0} n|c_n|^2 r^{-2n} \le \sum_{n>0} n |c_{-n}|^2 r^{2n}.

The result follows by letting r decrease to 1.

Milin's proof of Grunsky inequalities

If

 p(w)=\sum_{n=1}^N n^{-1} \lambda_n \Phi_n(w),

then

p(g(z))=\left( \sum_{n=1}^N n^{-1}\lambda_n z^n\right) %2B\left(\sum_{m=1}^\infty \sum_{n=1}^N \lambda_n c_{nm}z^{-m}\right).

Applying Milin's area theorem,

 \sum_{m=1}^\infty m\left|\sum_{n=1}^N  c_{mn}\lambda_n\right|^2 \le \sum_{n=1}^N {1\over n}|\lambda_n|^2.

(Equality holds here if and only if the complement of the image of g has Lebesgue measure zero.)

So a fortiori

 \sum_{m=1}^N m\left|\sum_{n=1}^N  c_{mn}\lambda_n\right|^2 \le \sum_{n=1}^N {1\over n} |\lambda_n|^2.

Hence the symmetric matrix

a_{mn}=(mn)^{1/2} c_{mn},

regarded as an operator on CN with its standard inner product, satisfies

\|Ax\|\le \|x\|.

So by the Cauchy-Schwarz inequality

|(Ax,y)|\le \|x\|\cdot\|y\|.

With  x_n=n^{-1/2}\lambda_n=\overline{y_n}, this gives the Grunsky inequality:

 \left |\sum_{m=1}^N\sum_{n=1}^N  c_{mn} \lambda_m\lambda_n\right|^2 \le \sum_{n=1}^N {1\over n} |\lambda_n|^2,

Criterion for univalence

Let g(z) be a holomorphic function on z > 1 with

 g(z) = z %2B b_0 %2B b_1 z^{-1} %2B b_2 z^{-2} %2B \cdots

Then g is univalent if and only if the Grunsky coefficients of g satisfy the Grunsky inequalities for all N.

In fact the conditions have already been shown to be necessary. To see sufficiency, note that

 \log {g(z) -g(\zeta)\over z -\zeta}=-\sum_{m,n\ge 1}  c_{mn}z^{-m}\zeta^{-n}

makes sense when |z| and |ζ| are large and hence the coefficients cmn are defined. If the Grunsky inequalities are satisfied then it is easy to see that the |cmn| are uniformly bounded and hence the expansion on the left hand side converges for |z|>1 and |ζ|>1. Exponentiating both sides, this implies that g is univalent.

Pairs of univalent functions

Let F(z) and g(ζ) be univalent holomorphic functions on |z| < 1 and |ζ| > 1, such that their images are disjoint in C. Suppose that these functions are normalized so that

g(\zeta)=\zeta%2Ba_0 %2B b_1\zeta^{-1} %2B b_2\zeta^{-2} %2B \cdots

and

 \displaystyle{F(z)=af(z)}

with a ≠ 0 and

f(z) = z %2B a_2z^2 %2B a_3z^3 %2B \cdots

The Grunsky matrix (cmn) of this pair of functions is defined for all non-zero m and n by the formulas:

 \log {g(\zeta)-g(\eta)\over \zeta -\eta} = -\sum_{m,n\ge 1} c_{mn} \zeta^{-m} \eta^{-n}
\log{g(\zeta) -  f(z)\over \zeta} -\log {g(\zeta)\over \zeta}= -\sum_{m,n\ge 1} c_{-m,n} z^m\zeta^{-n}
 \log {f(z)-f(w)\over z -w} -\log{f(z)\over z} -\log{f(w)\over w} =-\sum_{m,n\ge 1} c_{-m,-n} z^m w^n

with cm, –n = cn,m for m, n ≥ 1, so that cmn is a symmetric matrix.

In 1972 the American mathematician James Hummel extended the Grunsky inequalities to this matrix, proving that for any sequence of complex numbers λ±1, ..., λ±N

 |\sum_{n,m\ne 0} c_{mn}\lambda_m\lambda_n| \le \sum_{n\ne 0} {1\over |n|} |\lambda_n|^2.

The proof proceeds by computing the area of the image of the complement of the images of |z| < r < 1 under F and |ζ| > R >1 under g under a suitable Laurent polynomial h(w).

Let Φn and Φ-n denote the Faber polynomials of g and f(z^{-1})^{-1} and set

 h(w) =\sum_{n\ge 1} {\lambda_n\over n} \Phi_n(w) %2B \sum_{n\ge 1} {\lambda_{-n}\over n} \Phi_{-n}({a\over w}).

Then for |z| < 1

 h(F(z))= \sum_{n\ge 1} {\lambda_{-n}\over n}z^{-n} %2B\alpha%2B \sum_{n\ge 1} \alpha_n z^n,\,\,\, \alpha_n=\sum_m c_{-n,m}\lambda_m.

and for |ζ| > 1

 h(g(\zeta))=\sum_{n\ge 1} {\lambda_n\over n} \zeta^n%2B\beta %2B\sum_{n\ge 1} \beta_n \zeta^{-n},\,\,\,\beta_n=\sum_m c_{nm}\lambda_m.

The area equals

 \int |h^\prime(z)|^2 \, dx dy= {1\over 2 i}\int_{C_1} \overline{h}(z) h^\prime(z)\, dz -{1\over 2 i}\int_{C_2} \overline{h}(z) h^\prime(z)\, dz,

where C1 is the image of the circle |ζ| = R under g and C2 is the image of the circle |z| = r under F.

Hence

{1 \over \pi}\iint |h^\prime|^2 \, dx dy = \left[\sum_{n\ge 1} {1\over n} |\lambda_{-n}|^2 - \sum_{n\ge 1}  |\alpha_n|^2 r^{2n}\right] %2B\left[\sum_{n\ge 1}{1\over n}|\lambda_n|^2 - \sum_{n\ge 1}  |\beta_n|^2 R^{-2n}\right].

Since the area is positive, the right hand side must also be positive. Letting r increase to 1 and R decrease to 1, it follows that

\sum_{m\ne 0} |m| \left|\sum_{n\ne 0} c_{mn}\lambda_n\right|^2 \le \sum_{m\ne 0} {1\over |m|} |\lambda_m|^2

with equality if and only if the complement of the images has Lebesgue measure zero.

As in the case of a single function g, this implies the required inequality.

Unitarity

The matrix

 a_{mn}=\sqrt{|mn|} \cdot c_{mn}

of a single function g or a pair of functions F, g is unitary if and only if the complement of the image of g or the union of the images of F and g has Lebesgue measure zero. So, roughly speaking, in the case of one function the image is a slit region in the complex plane; and in the case of two functions the two regions are separated by a closed Jordan curve.

In fact the infinite matrix A acting on the Hilbert space of square summable sequences satisfies

\displaystyle{A^*A=I,}

But if J denotes complex conjugation of a sequence, then

\displaystyle{ JAJ=A^*, \,\, JA^*J=A}

since A is symmetric. Hence

 \displaystyle{AA^*=JA^*AJ=I}

so that A is unitary.

Equivalent forms of Grunsky inequalities

Goluzin inequalities

If g(z) is a normalized univalent function in |z| > 1, z1, ..., zN are distinct points with |zn| > 1 and α1, ..., αN are complex numbers, the Goluzin inequalities, proved in 1947 by the Russian mathematician Gennadi Mikhailovich Goluzin (1906-1853), state that

 \left|\sum_{m=1}^N \sum_{n=1}^N \alpha_m\alpha_n \log{g(z_m) -g(z_n)\over z_m -z_n}\right|^2 \le \sum_{m=1}^N\sum_{n=1}^N \alpha_m\overline{\alpha_n} \log{1\over 1-(z_m\overline{z_n})^{-1}}.

To deduce them from the Grunsky inequalities, let

\displaystyle{\lambda_k=\sum_{n=1}^N \alpha_n z_n^{-k}.}

for k > 0.

Conversely the Grunsky inequalities follow from the Goluzin inequalities by taking

\displaystyle{\alpha_m={1\over N} \sum_{n=1}^N\lambda_n z_n^m.}

where

\displaystyle{z_n=re^{2\pi i n\over N}}

with r > 1, tending to ∞.

Bergman–Schiffer inequalities

Bergman & Schiffer (1951) gave another derivation of the Grunsky inequalities using reproducing kernels and singular integral operators in geometric function theory; a more recent related approach can be found in Baranov & Hedenhalm (2008). Let f(z) be a normalized univalent function in |z| < 1, let z1, ..., zN be distinct points with |zn| < 1 and let α1, ..., αN be complex numbers.

The Bergman-Schiffer inequalities state that

 \left|\sum_{m=1}^N \sum_{n=1}^N \alpha_m \alpha_n\left[ {f^\prime(z_m)f^\prime(z_n) \over(f(z_m)-f(z_n))^2} - {1\over (z_m-z_n)^{2}}\right]  \right| \le \sum_{m=1}^N \sum_{n=1}^N \alpha_m\overline{\alpha_n} {1\over (1-z_m\overline{z_n})^2}.

To deduce these inequalities from the Grunsky inequalities, set

\displaystyle{\lambda_k=k\sum_{n=1}^N \alpha_n z_n^k.}

for k > 0.

Conversely the Grunsky inequalities follow from the Bergman-Schiffer inequalities by taking

\displaystyle{\alpha_m={1\over N} \sum_{n=1}^N{1\over n} \lambda_n z_n^m.}

where

\displaystyle{z_n=re^{2\pi i n\over N}}

with r < 1, tending to 0.

Applications

The Grunsky inequalities imply many inequalities for univalent functions. They were also used by Schiffer and Charzynski in 1960 to give a completely elementary proof of the Bieberbach conjecture for the fourth coefficient; a far more complicated proof had previously been found by Schiffer and Garabedian in 1955. In 1968 Pedersen and Ozawa independently used the Grunsky inequalities to prove the conjecture for the sixth coefficient.[1][2]

In the proof of Schiffer and Charzynski, if

f(x)=z%2Ba_2 z^2 %2B a_3 z^3 %2B a_4z^4 %2B \cdots

is a normalized univalent function in |z| < 1, then

 g(z) =f(z^2)^{-1/2}=z %2B b_1 z^{-1} %2B b_3 z^{-3} %2B \cdots

is an odd univalent function in |z| > 1.

Combining Gronwall's area theorem for f with the Grunsky inequalities for the first 2 x 2 minor of the Grunsky matrix of g leads to a bound for |a4| in terms of a simple function of a2 and a free complex parameter. The free parameter can be chosen so that the bound becomes a function of half the modulus of a2 and it can then be checked directly that this function is no greater than 4 on the range [0,1].

As Milin showed, the Grunsky inequalities can be exponentiated. The simplest case proceeds by writing

 \log {g(z)-g(\zeta)\over z-\zeta} = - \sum_{n\ge 1} a_n(\zeta^{-1})z^{-n}.

with an(w) holomorphic in |w| < 1.

The Grunsky inequalities, with λn = wn imply that

 \displaystyle{\sum_{n\ge 1} n|a_n(w)|^2 \le - \log (1-|w|^2).}

On the other hand if

  \displaystyle{\sum_{m\ge 0} b_m t^m = \exp \sum_{n\ge 1} a_n t^n }

as formal power series, then the first of the Lebedev-Milin inequalities (1965) states that[3][4]

  \displaystyle{\sum_{n\ge 0} |b_n|^2 \le \exp \sum_{n\ge 1} n |a_n|^2. }

Equivalently the inequality states that if g(z) is a polynomal with g(0) = 0, then

{1\over 2\pi} \int_0^{2\pi} |e^g|^2 \, d\theta \le e^A,

where A is the area of g(D),

To prove the inequality, note that the coefficients are determined by the recursive formula

 \displaystyle{b_n={1\over n}\sum_{m=1}^n ma_m b_{n-m}}

so that by the Cauchy-Schwarz inequality

 \displaystyle{|b_n|^2 \le {1\over n} \sum m^2 |a_m|^2 |b_{n-m}|^2.}

The quantities cn obtained by imposing equality here:

\displaystyle{ c_n= {1\over n} \sum m^2 |a_m|^2 c_{n-m}}

satisfy  |b_n|^2 \le c_n and hence, reversing the steps,

\displaystyle{\sum |b_n|^2 \le  \sum c_n = \exp \sum_{m\ge 1} m|a_m|^2.}

In particular defining bn(w) by the identity

\displaystyle{\sum b_n(\zeta^{-1}) z^{-n} = \exp \sum a_m(\zeta^{-1}) z^{-m} = {g(z)-g(\zeta)\over z-\zeta},}

the following inequality must hold for |w| < 1

\displaystyle{ \sum_{n\ge 0} |b_n(w)|^2 \le (1-|w|^2)^{-1}.}

Beurling transform

The Beurling transform (also called the Beurling-Ahlfors transform and the Hilbert transform in the complex plane) provides one of the most direct methods of proving the Grunsky inequalities, following Bergman & Schiffer (1951) and Baranov & Hedenhalm (2008).

The Beurling transform is defined on L2(C) as the operation of multiplication by z/\overline{z} on Fourier transforms. It thus defines a unitary operator. It can also be defined directly as a principal value integral[5]

 (Th)(w)=\lim_{\varepsilon\rightarrow 0} {1\over \pi}\iint_{|z-w|\ge \varepsilon} {h(z)\over (z-w)^2} \,dx dy.

For any bounded open region Ω in C it defines a bounded operator TΩ from the conjugate of the Bergman space of onto the Bergman space of Ω: a square integrable holomorphic function is extended to 0 off Ω to produce a function in L2(C) to which T is applied and the result restricted to Ω, where it is holomorphic. If f is a holomorphic univalent map from the unit disk D onto Ω then the Bergman space of Ω and its conjugate can be identified with that of D and TΩ becomes the singular integral operator with kernel

 K_f(z,w)={f^\prime(z)f^\prime(w)\over (f(z)-f(w))^2}.

It defines a contraction. On the other hand it can be checked that TD = 0 by computing directly on powers \overline{z}^n using Stokes theorem to transfer the integral to the boundary.

It follows that the operator with kernel

  {f^\prime(z)f^\prime(w)\over (f(z)-f(w))^2} \,-\,{1\over (z-w)^2}={\partial^2\over\partial z \partial w} \log {f(z) -f(w)\over z-w} =-\sum_{m,n\ge 1} mnc_{mn} z^{m-1}w^{n-1}

acts as a contraction on the conjugate of the Bergman space of D. Hence, if

 p(z)=\lambda_1 %2B \lambda_2 \overline{z} %2B \lambda_3 \overline{z}^2 %2B \cdots %2B \lambda_N \overline{z}^{N-1},

then

 \sum_{m=1}^N \left|\sum_{n=1}^N c_{mn}\lambda_n\right|^2=\|(T_f-T_z)p\|^2 =\|T_fp\|^2\le \|p\|^2= \sum_{n=1}^N {1\over n} |\lambda_n|^2.

Grunsky operator and Fredholm eigenvalues

If Ω is a bounded domain in C with smooth boundary, the operator TΩ can be regarded as a bounded antilinear contractive operator on the Bergman space A2(Ω). Its realization on D using a univalent function f mapping D onto Ω shows that TΩ is a Hilbert-Schmidt operator. TΩ is called the Grunsky operator of Ω (or f). It is given by the formula

 (T_\Omega u)(z) = \lim_{\varepsilon\rightarrow 0} {1\over \pi} \iint_{|z-w|\ge \varepsilon} {\overline{u(z)}\over (z-w)^2} \,\, dx dy

for u in the Hilbert space H= A2(Ω).

The antilinear operator T = TΩ satisfies the self-adjointness relation

\displaystyle{(Tu,v)=(Tv,u)}

for u, v in H.

Thus A = T2 is a compact self-adjont linear operator on H with

\displaystyle{(Au,u)=(Tu,Tu)=\|Tu\|^2\ge 0,}

so that A is a positive operator. By the spectral theorem for compact self-adjoint operators, there is an orthonormal basis un of H consisting of eigenvectors of A:

 \displaystyle{Au_n=\mu_n u_n,}

where μn is non-negative by the positivity of A. Hence

 \displaystyle{\mu_n=\lambda_n^2}

with λn ≥ 0. Since T commutes with A, it leaves its eigenspaces invariant. The positivity relation shows that it acts trivially on the zero eigenspace. The other non-zero eigenspaces are all finite dimensional and mutually orthogonal. Thus an orthonormal basis can be chosen on each eigenspace so that:

 \displaystyle{Tu_n=\lambda_n u_n.}

(Note that  T(iu_n)=-\lambda_n iu_n by antilinearity of T.)

The non-zero λn (or sometimes their reciprocals) are called the Fredholm eigenvalues of Ω:

\displaystyle{ 0\le \lambda_n \le \|T\| \le 1.}

If Ω is a bounded domain that is not a disk, Ahlfors showed that

\displaystyle{\|T_\Omega\|< 1.}

The Fredholm determinant for the domain Ω is defined by[6][7]

\Delta_\Omega=\det (I- T_\Omega^2) =\prod (1-\lambda_n^2).

Note that this makes sense because A = T2 is a trace class operator.

Schiffer & Hawley (1962) showed that. if 0 lies in Ω and f fixes 0, then[8][9]

 \Delta_\Omega= -{1\over 12\pi}\left[\|\partial_z \log f^\prime\|^2_D %2B\|\partial_z \log g^\prime\|^2_{D^c}  - 2 \|\partial_z\log f(z)/z \|^2_D  - 2 \|\partial_z\log g(z)/z \|^2_{D^c}\right].

Here the norms are in the Bergman spaces of D and its complement Dc and g is a univalent map from Dc onto Ωc fixing ∞.

A similar formula applies in the case of a pair of univalent functions (see below).

Singular integral operators on a closed curve

Let Ω be a bounded simply connected domain in C with smooth boundary C = ∂Ω. Thus there is a univalent holomorphic map f from the unit disk D onto Ω extending to a smooth map between the boundaries S1 and C.

Notes

  1. ^ Duren 1983, p. 131-133
  2. ^ Koepf 2007
  3. ^ Duren 1983, p. 143-144
  4. ^ Apart from the elementary proof of this result presented here, there are several other analytic proofs in the literature. Nikolski (2002, p. 220), following de Branges, notes that it is a consequence of standard inequalities connected with reproducing kernels. Widom (1988) observed that it was an immediate consequence of Szegő's limit formula (1951). Indeed if f is the real-valued trigonometric polynomial on the circle given as twice the real part of a polynomial g(z) vanishing at 0 on the unit disk, Szegő's limit formula states that the Toeplitz determinants of ef increase to eA where A is the area of g(D). The first determinant is by definition just the constant term in ef = |eg|2.
  5. ^ Ahlfors 1966
  6. ^ Schiffer 1959, p. 261
  7. ^ Schiffer & Hawley 1962, p. 246
  8. ^ Schiffer & Hawley 1962, p. 245-246
  9. ^ Takhtajan & Teo 2006

References