The Power of Depth for Feedforward Neural Networks

arXiv:1512.03965v4 [cs.LG] 9 May 2016

Ronen Eldan
Weizmann Institute of Science
ronen.eldan@weizmann.ac.il

Ohad Shamir
Weizmann Institute of Science
ohad.shamir@weizmann.ac.il

Abstract
We show that there is a simple (approximately radial) function on Rd , expressible by a small 3-layer
feedforward neural networks, which cannot be approximated by any 2-layer network, to more than a
certain constant accuracy, unless its width is exponential in the dimension. The result holds for virtually
all known activation functions, including rectified linear units, sigmoids and thresholds, and formally
demonstrates that depth – even if increased by 1 – can be exponentially more valuable than width for
standard feedforward neural networks. Moreover, compared to related results in the context of Boolean
functions, our result requires fewer assumptions, and the proof techniques and construction are very
different.

1

Introduction and Main Result

Learning via multi-layered artificial neural networks, a.k.a. deep learning, has seen a dramatic resurgence of
popularity over the past few years, leading to impressive performance gains on difficult learning problems,
in fields such as computer vision and speech recognition. Despite their practical success, our theoretical
understanding of their properties is still partial at best.
In this paper, we consider the question of the expressive power of neural networks of bounded size. The
boundedness assumption is important here: It is well-known that sufficiently large depth-2 neural networks,
using reasonable activation functions, can approximate any continuous function on a bounded domain (Cybenko [1989], Hornik et al. [1989], Funahashi [1989], Barron [1994]). However, the required size of such
networks can be exponential in the dimension, which renders them impractical as well as highly prone to
overfitting. From a learning perspective, both theoretically and in practice, our main interest is in neural
networks whose size is bounded.
For a network of bounded size, a basic architectural question is how to trade off between its width and
depth: Should we use networks that are narrow and deep (many layers, with a small number of neurons per
layer), or shallow and wide? Is the “deep” in “deep learning” really important? Or perhaps we can always
content ourselves with shallow (e.g. depth-2) neural networks?
Overwhelming empirical evidence as well as intuition indicates that having depth in the neural network
is indeed important: Such networks tend to result in complex predictors which seem hard to capture using
shallow architectures, and often lead to better practical performance. However, for the types of networks
used in practice, there are surprisingly few formal results (see related work below for more details).
In this work, we consider fully connected feedforward neural networks, using a linear output neuron
and some non-linear activation function on the other neurons, such as the commonly-used rectified linear
unit (ReLU, σ(z) = max{z, 0}), as well as the sigmoid (σ(z) = (1 + exp(−z))−1 ) and the threshold
(σ(z) = 1 {z ≥ 0}). Informally speaking, we consider the following question: What functions on Rd
1

expressible by a network with `-layers and w neurons per layer, that cannot be well-approximated by any
network with < ` layers, even if the number of neurons is allowed to be much larger than w?
More specifically, we consider the simplest possible case, namely the difficulty of approximating functions computable by 3-layer networks using 2-layer networks, when the networks are feedforward and fully
connected. Following a standard convention, we define a 2-layer network of width w on inputs in Rd as
x 7→

w
X

vi σ (hwi , xi + bi )

(1)

i=1

where σ : R → R is the activation function, and vi , bi ∈ R, wi ∈ Rd , i = 1, . . . , w are parameters of the
network. This corresponds to a set of w neurons
P computing x 7→ σ(hwi ,1xi + bi ) in the first layer, whose
output is fed to a linear output neuron x 7→ w
i=1 vi xi in the second layer . Similarly, a 3-layer network of
width w is defined as


w
w
X
X
vi,j σ (hwi,j , xi + bi,j ) + ci  ,
(2)
ui σ 
i=1

j=1

where ui , ci , vi,j , bi,j ∈ R, wi,j ∈ Rd , i, j = 1, . . . , w are parameters of the network. Namely, the outputs
of the neurons in the first layer are fed to neurons in the second layer, and their outputs in turn are fed to a
linear output neuron in the third layer.
Clearly, to prove something on the separation between 2-layer and 3-layer networks, we need to make
some assumption on the activation function σ(·) (for example, if σ(·) is the identity, then both 2-layer and
3-layer networks compute linear functions, hence there is no difference in their expressive power). All we
will essentially require is that σ(·) is universal, in the sense that a sufficiently large 2-layer network can
approximate any univariate Lipschitz function which is non-constant on a bounded domain. More formally,
we use the following assumption:
Assumption 1. Given the activation function σ, there is a constant cσ ≥ 1 (depending only on σ) such that
the following holds: For any L-Lipschitz function f : R → R which is constant outside a bounded interval
RL
[−R, R], and for any δ, there exist scalars a, {αi , βi , γi }w
i=1 , where w ≤ cσ δ , such that the function
h(x) = a +

w
X

αi · σ(βi x − γi )

i=1

satisfies
sup |f (x) − h(x)| ≤ δ.
x∈R

This assumption is satisfied by the standard activation functions we are familiar with. First of all, we
provide in Appendix A a constructive proof for the ReLU function. For the threshold, sigmoid, and more
general sigmoidal functions (e.g. monotonic functions which satisfy limz→∞ σ(z) = a, limz→−∞ σ(z) = b
for some a 6= b in R), the proof idea is similar, and implied by the proof of Theorem 1 of Debao [1993]2 .
1
Note that sometimes one also adds a constant bias parameter b to the output neuron, but this can be easily simulated by a
“constant” neuron i in the first layer where wi = 0 and vi , bi are chosen appropriately. Also, sometimes the output neuron is
defined to have a non-linearity as well, but we stick to linear output neurons, which is a very common and reasonable assumption
for networks computing real-valued predictions.
2
Essentially, a single neuron with such a sigmoidal activation can express a (possibly approximate) single-step function, a
combination of w such neurons can express a function with w such steps, and any L-Lipschitz function which is constant outside
[−R, R] can be approximated to accuracy δ with a function involving O(RL/δ) steps.

2

Finally, one can weaken the assumed bound on w to any poly(R, L, 1/δ), at the cost of a worse polynomial
dependence on the dimension d in Thm. 1 part 1 below (see Subsection 4.4 for details).
In addition, for technical reasons, we will require the following mild growth and measurability conditions, which are satisfied by virtually all activation functions in the literature, including the examples
discussed earlier:
Assumption 2. The activation function σ is (Lebesgue) measurable and satisfies
|σ(x)| ≤ C(1 + |x|α )
for all x ∈ R and for some constants C, α > 0.
Our main result is the following theorem, which implies that there are 3-layer networks of width polynomial in the dimension d, which cannot be arbitrarily well approximated by 2-layer networks, unless their
width is exponential in d:
Theorem 1. Suppose the activation function σ(·) satisfies assumption 1 with constant cσ , as well as assumption 2. Then there exist universal constants c, C > 0 such that the following holds: For every dimension
d > C, there is a probability measure µ on Rd and a function g : Rd → R with the following properties:
√
1. g is bounded in [−2, +2], supported on {x : kxk ≤ C d}, and expressible by a 3-layer network of
width Ccσ d19/4 .
2. Every function f , expressed by a 2-layer network of width at most cecd , satisfies
Ex∼µ (f (x) − g(x))2 ≥ c.
The proof is sketched in Sec. 2, and is formally presented in Sec. 4. Roughly speaking, g approximates
a certain radial function g̃, depending only on the norm of the input. With 3 layers, approximating radial
functions (including g̃) to arbitrary accuracy is straightforward, by first approximating the squared norm
function, and then approximating the univariate function acting on the norm. However, performing this
approximation with only 2 layers is much more difficult, and the proof shows that exponentially many
neurons are required to approximate g̃ to more than constant accuracy. We conjecture (but do not prove) that
a much wider family of radial functions also satisfy this property.
We make the following additional remarks about the theorem:
Remark 1 (Activation function). The theorem places no constraints on the activation function σ(·) beyond
assumptions 1 and 2. In fact, the inapproximability result for the function g̃ holds even if the activation functions are different across the first layer neurons, and even if they are chosen adaptively (possibly depending
on g̃), as long as they satisfy assumption 2.
Remark 2 (Constraints on the parameters). The theorem places no constraints whatsoever on the parameters of the 2-layer networks, and they can take any values in R. This is in contrast to related depth separation
results in the context of threshold circuits, which do require the size of the parameters to be constrained (see
discussion of related work below).
Remark 3 (Properties of g). At least for specific activation functions such as the ReLU, sigmoid, and
threshold, the proof construction implies that g is poly(d)-Lipschitz, and the 3-layer network expressing it
has parameters bounded by poly(d).
3

Related Work
On a qualitative level, the question we are considering is similar to the question of Boolean circuit lower
bounds in computational complexity: In both cases, we consider functions which can be represented as a
combination of simple computational units (Boolean gates in computational complexity; neurons in neural
networks), and ask how large or how deep this representation needs to be, in order to compute or approximate
some given function. For Boolean circuits, there is a relatively rich literature and some strong lower bounds.
A recent example is the paper Rossman et al. [2015], which shows for any d ≥ 2 an explicit depth d,
linear-sized circuit on {0, 1}n , which cannot be non-trivially approximated by depth d − 1 circuits of size
polynomial in n. That being said, it is well-known that the type of computation performed by each unit in
the circuit can crucially affect the hardness results, and lower bounds for Boolean circuits do not readily
translate to neural networks of the type used in practice, which are real-valued and express continuous
functions. For example, a classical result on Boolean circuits states that the parity function over {0, 1}d
cannot be computed by constant-depth Boolean circuits whose size is polynomial in d (see for instance
Håstad [1986]). Nevertheless, the parity function can in fact be easily computed by a simple 2-layer, O(d)width real-valued neural network with most reasonable activation functions3 .
A model closer to ours is a threshold circuit, which is a neural network where all neurons (including
the output neuron) has a threshold activation function, and the input is from the Boolean cube (see Parberry [1994] for a survey). For threshold circuits, the main known result in our context is that computing
inner products mod 2 over d-dimensional Boolean vectors cannot be done with a 2-layer network with
poly(d)-sized parameters and poly(d) width, but can be done with a small 3-layer network (Hajnal et al.
[1993]). Note that unlike neural networks in practice, the result in Hajnal et al. [1993] is specific to the
non-continuous threshold activation function, and considers hardness of exact representation of a function
by 2-layer circuits, rather than merely approximating it. Following the initial publication of our paper, we
were informed (Martens [2015]) that the proof technique, together with techniques in the papers (Maass
et al. [1994], Martens et al. [2013])), can possibly be used to show that inner product mod 2 is also hard
to approximate, using 2-layer neural networks with continuous activation functions, as long as the network
parameters are constrained to be polynomial in d, and that the activation function satisfies certain regularity
conditions4 . Even so, our result does not pose any constraints on the parameters, nor regularity conditions
beyond assumptions 1,2. Moreover, we introduce a new proof technique which is very different, and demonstrate hardness of approximating not the Boolean inner-product-mod-2 function, but rather functions in Rd
with a simple geometric structure (namely, radial functions).
Moving to networks with real-valued outputs, one related field is arithmetic circuit complexity (see
Shpilka and Yehudayoff [2010] for a survey), but the focus there is on computing polynomials, which can
be thought of as neural networks where each neuron computes a linear combination or a product of its inputs.
Again, this is different than most standard neural networks used in practice, and the results and techniques
do not readily translate.
Recently, several works in the machine learning community attempted to address questions similar to
the one we consider here. Pascanu et al. [2013], Montufar et al. [2014] consider the number of linear regions
3
See Rumelhart et al. [1986], Figure 6, where reportedly the structure was even found automatically by back-propagation.
For a threshold activation function σ(z) = 1 {z ≥ 0} and input x = (x1 , . . . , xd ) ∈ {0, 1}d , the network is given by
Pd+1
Pd
i+1
1
1
x 7→
σ
i=1 (−1)
j=1 xj − i + 2 . In fact, we only need σ to satisfy σ(z) = 1 for z ≥ 2 and σ(z) = 0 for

z ≤ − 21 , so the construction easily generalizes to other activation functions (such as a ReLU or a sigmoid), possibly by using a
small linear combination of them to represent such a σ.
4
See remark 20 in Martens et al. [2013]. These conditions are needed for constructions relying on distributions over a finite set
(such as the Boolean hypercube). However, since we consider continuous distributions on Rd , we do not require such conditions.

4

which can be expressed by ReLU networks of a given width and size, and Bianchini and Scarselli [2014]
consider the topological complexity (via Betti numbers) of networks with certain activation functions, as a
function of the depth. Although these can be seen as measures of the function’s complexity, such results
do not translate directly to a lower bound on the approximation error, as in Thm. 1. Delalleau and Bengio
[2011], Martens and Medabalimi [2014] and Cohen et al. [2015] show strong approximation hardness results
for certain neural network architectures (such as polynomials or representing a certain tensor structure),
which are however fundamentally different than the standard neural networks considered here.
Quite recently, Telgarsky [2015] gave a simple and elegant construction showing that for any k, there are
k-layer, O(1) wide ReLU networks on one-dimensional data, which can express a sawtooth function on [0, 1]
which oscillates O(2k ) times, and moreover, such a rapidly oscillating function cannot be approximated by
poly(k)-wide ReLU networks with o(k/ log(k)) depth. This also implies regimes with exponential separation, e.g. that there are k 2 -depth networks, which any approximating k-depth network requires Ω(exp(k))
width. These results demonstrate the value of depth for arbitrarily deep, standard ReLU networks, for a
single dimension and using functions which have an exponentially large Lipschitz parameter. In this work,
we use different techniques, to show exponential separation results for general activation functions, even if
the number of layers changes by just 1 (from two to three layers), and using functions in Rd whose Lipschitz
parameter is polynomial in d.

2

Proof Sketch

In a nutshell, the 3-layer network we construct approximates a radial function with bounded support (i.e.
one which depends on the input x only via its Euclidean norm kxk, and is 0 for any x whose norm is larger
than some threshold). With 3 layers, approximating radial functions is rather straightforward: First, using
assumption 1, we can construct a linear combination of neurons expressing the univariate mapping z 7→ z 2
arbitrarily well in any bounded domain. Therefore, by adding these combinations together, P
one for each
coordinate, we can have our network first compute (approximately) the mapping x 7→ kxk2 = i x2i inside
any bounded domain, and then use the next layer to compute some univariate function of kxk2 , resulting
in an approximately radial function. With only 2 layers, it is less clear how to approximate such radial
functions. Indeed, our proof essentially indicates that approximating radial functions with 2 layers can
require exponentially large width.
To formalize this, note that if our probability measure µ has a well-behaved density function which can
be written as ϕ2 (x) for some function ϕ, then the approximation guarantee in the theorem, Eµ (f (x) −
g(x))2 , can be equivalently written as
Z
Z
(f (x) − g(x))2 ϕ2 (x)dx =
(f (x)ϕ(x) − g(x)ϕ(x))2 dx = kf ϕ − gϕk2L2 .
(3)
In particular, we will consider a density function which equals ϕ2 (x), where ϕ is the inverse Fourier transform of the indicator 1 {x ∈ B}, B being the origin-centered unit-volume Euclidean ball (the reason for this
choice will become evident later). Before continuing, we note that a formula for ϕ can be given explicitly
(see Lemma 2), and an illustration of it in d = 2 dimensions is provided in Figure 2. Also, it is easily
verified thatR ϕ2 (x) is indeed
It is clearly non-negative, and by isometry of the Fourier
R a density function:
R
transform, ϕ2 (x)dx = ϕ̂2 (x)dx = 1 {x ∈ B}2 dx, which equals 1 since B is a unit-volume ball.
Our goal now is to lower bound the right hand side of Eq. (3). To continue, we find it convenient to
consider the Fourier transforms fc
ϕ, gc
ϕ of the functions f ϕ, gϕ, rather than the functions themselves. Since

5

Figure 1: The left figure represents ϕ(x) in d = 2 dimensions. The right figure represents a cropped and
re-scaled version, to better show the oscillations of ϕ beyond the big origin-centered bump. The density of
the probability measure µ is defined as ϕ2 (·)
the Fourier transform is isometric, the above equals
kfc
ϕ − gc
ϕk2L2 .
Luckily, the Fourier transform of functions expressible by a 2-layer network has a very particular form.
Specifically, consider any function of the form
f (x) =

k
X

fi (hvi , xi),

i=1

where fi : R → R (such as 2-layer networks as defined earlier). Note that f may not be square-integrable, so
formally speaking it does not have a Fourier transform in the standard sense of a function on Rd . However,
assuming |fi (x)| grows at most polynomially as x → ∞ or x → −∞, it does have a Fourier transform
in the more general sense of a tempered distribution (we referSthe reader to the proof for a more formal
discussion). This distribution can be shown to be supported on i span{vi }: namely, a finite collection of
lines5 . The convolution-multiplication principle implies that fc
ϕ equals fˆ ? ϕ̂, or the convolution of fˆ with
S
ˆ
the indicator of a unit-volume ball B. Since f is supported on i span{vi }, it follows that
Supp(fc
ϕ) ⊆ T :=

k
[

(span{vi } + B) .

i=1

In words, the support of fc
ϕ is contained in a union of tubes of bounded radius passing through the origin.
This is the key property of 2-layer networks we will use to derive our main theorem. Note that it holds
5
Roughly speaking, this is because each function x 7→ fi (hvi , xi) is constant in any direction perpendicular to vi , hence do
not have non-zero Fourier components in those directions. In one dimension, this can be seen by the fact that the Fourier transform
of the constant 0 function is the Dirac delta function, which equals 0 everywhere except at the origin.

6

regardless of the exact shape of the fi functions, and hence our proof will also hold if the activations in the
network are different across the first layer neurons, or even if they are chosen in some adaptive manner.
To establish our theorem, we will find a function g expressible by a 3-layer network, such that gc
ϕ has
a constant distance (in L2 space) from any function supported on T (a union of k tubes as above). Here is
where high dimensionality plays a crucial role: Unless k is exponentially large in the dimension, the domain
T is very sparse when one considers large distances from the origin, in the sense that
V old−1 (T ∩ rSd−1 )
. ke−d
V old−1 (rSd−1 )
(where Sd−1 is the d-dimensional unit Euclidean sphere, and V old−1 is the d − 1-dimensional Hausdorff
measure) whenever r is large enough with respect to the radius of B. Therefore, we need to find a function
g so that gc
ϕ has a lot of mass far away from the origin, which will ensure that kfc
ϕ − gc
ϕk2L2 will be large.
Specifically, we wish to find a function g so that gϕ is radial (hence gc
ϕ is also radial, so having large mass
in any direction implies large mass in all directions), and has a significant high-frequency component, which
implies that its Fourier transform has a significant portion of its mass outside of the ball rB.
The construction and analysis of this function constitutes the technical bulk of the proof. The main
difficulty in this step is that even if the Fourier transform ĝ of g has some of its L2 mass on high frequencies,
it is not clear that this will also be true for gc
ϕ = g ? 1 {B} (note that while convolving with a Euclidean ball
increases the average distance from the origin in the L1 sense, it doesn’t necessarily do the same in the L2
sense).
We overcome this difficulty by considering a random superposition of indicators of thin shells: Specifically, we consider the function
N
X
g̃(x) =
i gi (x),
(4)
i=1

where i ∈ {−1, +1}, N = poly(d), and
√ gi (x) = 1 {kxk ∈ ∆i }, where ∆i are disjoint intervals of
width O(1/N ) on values in the range Θ( d). Note that strictly speaking, we cannot take our hard-toapproximate function g to equal g̃, since g̃ is discontinuous and therefore cannot be expressed by a 3-layer
neural network with continuous activations functions. However, since our probability distribution ϕ2 can
be shown to have bounded density in the support of Eq. (4), we can use a 3-layer network toP
approximate
2
such a function arbitrarily well with respect to the distribution ϕ (for example, by computing i i gi (x) as
above, with each hard indicator function gi replaced by a Lipschitz function, which differs from gi on a set
with arbitrarily small probability mass). Letting the function g be such a good approximation, we get that
if no 2-layer network can approximate the function g̃ in Eq. (4), then it also cannot approximate its 3-layer
approximation g.
Let us now explain why the function defined in Eq. (4) gives us what we need. For large N , each gi is
supported on a thin Euclidean
shell, hence gi ϕ is approximately the same as ci gi for someP
constant ci . As
PN
c̃
a result, g̃(x)ϕ(x) ≈ i=1 i ci gi (x), so its Fourier transform (by linearity) is gϕ(w) ≈ N
i=1 i ci gˆi (w).
Since gi is a simple indicator function, its Fourier transform gˆi (w) is not too difficult to compute explicitly,
and involves an appropriate Bessel function which turns out to have a sufficiently large mass sufficiently far
away from the origin.
Knowing that each summand gi has a relatively large mass on high frequencies, our only remaining
objective is to find a choice for the signs i so that the entire sum will have the same property. This is attained
by a random choice of signs: it is an easy observation that given an orthogonal projection P in a P
Hilbert space

H, and any sequence of vectors v1 , ..., vN ∈ H such that |P vi | ≥ δ|vi |, one has that E |P i i vi |2 ≥
7

P
δ 2 i |vi |2 when the signs i are independent Bernoulli ±1 variables. Using this observation with P being
the projection onto the subspace spanned by functions supported on high frequencies and with the functions
ĝi , it follows that there is at least one choice of the i ’s so that a sufficiently large portion of g̃’s mass is on
high frequencies.

3

Preliminaries

We begin by defining some of the standard notation we shall use. We let N and R denote the natural and real
numbers, respectively. Bold-faced letters denote vectors in d-dimensional Euclidean space Rd , and plainfaced letters to denote either scalars or functions (distinguishing
R 2 between them should be clear from context).
L2 denotes the space of squared integrable functions ( x f (x)dx
< ∞, where the integration is over Rd ),
R
and L1 denotes the space of absolutely integrable functions ( x |f (x)|dx < ∞). k · k denotes
R the Euclidean
=
f (x)g(x)dx),
f,
g,
we
have
hf,
gi
norm, h·, ·iL2 denotes inner product in L2 space (for functions
L2
R
2
2
k·kL2 denotes the standard norm in L2 space ( kf kL2 = x f (x) dx), and k·kL2 (µ) denotes the L2 space
R
norm weighted by a probability measure µ (namely kf k2L2 (µ) = f (x)2 dµ(x)). Given two functions f, g,
we let f g be shorthand for the function x 7→ f (x) · g(x), and f + g be shorthand for x 7→ f (x) + g(x).
Given two sets A, B in Rd , we let A + B = {a + b : a ∈ A, b ∈ B} and AC = {a ∈ Rd : a ∈
/ A}
Fourier Transform. For a function f : R → R, our convention for the Fourier transform is
Z
ˆ
f (w) =
exp (−2πixw) f (x)dx
R

whenever the integral is well defined. This is generalized for f : Rd → R by
Z
fˆ(w) =
exp (−2πihx, wi) f (x)dx.

(5)

Rd

Radial Functions. A radial function f : Rd 7→ R is such that f (x) = f (x0 ) for any x, x0 such that
kxk = kx0 k. When dealing with radial functions, which are invariant to rotations, we will somewhat abuse
notation and interchangeably use vector arguments x to denote the value of the function at x, and scalar
arguments r to denote the value of the same function for any vector of norm r. Thus, for a radial function
f : Rd → R, f (r) equals f (x) for any x such that kxk = r.
Euclidean Spheres and Balls. Let Sd−1 be the unit Euclidean sphere in Rd , Bd be the d-dimensional
unit Euclidean ball, and let Rd be the radius so that Rd Bd has volume one. By standard results, we have the
following useful lemma:
q
√
√
1/d
Lemma 1. Rd = π1 Γ d2 + 1
, which is always between 15 d and 12 d.
Bessel Functions. Let Jν : R 7→ R denote the Bessel function
kind, of order ν. The Bessel
P∞of the first
2m+ν
(−1)m
function has a few equivalent definitions, for example Jν (x) = m=0 m!Γ(m+ν+1) x2
where Γ(·) is
the Gamma function. Although itqdoes nothave a closed form,
J
(x)
has
an
oscillating
shape,
which for
ν


2
asymptotically large x behaves as πx
cos − (2ν+1)π
+ x . Figure 3 illustrates the function for ν = 20. In
4
appendix C, we provide additional results and approximations for the Bessel function, which are necessary
for our proofs.

8

0.25
0.2
0.15
0.1
0.05
0
−0.05
−0.1
−0.15
−0.2

0

50

100

150

Figure 2: Bessel function of the first kind, J20 (·)

4

Proof of Thm. 1

In this section, we provide the proof of Thm. 1. Note that some technical proofs, as well as some important
technical lemmas on the structure of Bessel functions, are deferred to the appendix.

4.1

Constructions

As discussed in Sec. 2, our theorem rests on constructing a distribution µ and an appropriate function g,
which is easy to approximate (w.r.t. µ) by small 3-layer networks, but difficult to approximate using 2-layer
networks. Thus, we begin by formally defining g, µ that we will use.
dµ
First, µ will be defined as the measure whose density is dx
= ϕ2 (x), where ϕ(x) is the Fourier transform
of the indicator
of a unit-volume
Euclidean ball 1 {w ∈ Rd Bd }. Note that since the Fourier transform is an
R
R
2
isometry, Rd ϕ(x) dx = Rd 1 {w ∈ Rd Bd }2 dw = 1, hence µ is indeed a probability measure. The form
of ϕ is expressed by the following lemma:
Lemma 2. Let ϕ(x) be the Fourier transform of 1 {w ∈ Rd Bd }. Then


Rd d/2
ϕ(x) =
Jd/2 (2πRd kxk).
kxk
The proof appears in Appendix B.1.
To define our hard-to-approximate function, we introduce some notation. Let α ≥ 1 and γ be some
large numerical constants to be determined later, and set N = γd2 , assumed to be an integer (essentially,
we need α, γ to be sufficiently large so that all the lemmas we construct below would hold). Consider the
intervals





√
√
i
i−1
α d, 1+
α d , i = 1, 2, . . . , N.
∆i =
1+
N
N
We split the intervals to “good” and “bad” intervals using the following definition:
Definition 1. ∆i is a good interval (or equivalently, i is good) if for any x ∈ ∆i
2
Jd/2
(2πRd x) ≥

Otherwise, we say that ∆i is a bad interval.
9

1
.
80πRd x

For any i, define
(
1 {x ∈ ∆i }
gi (x) =
0

i good
i bad

(6)

By definition of a “good” interval and Lemma 2, we see that gi is defined to be non-zero, when the value of
ϕ on the corresponding interval ∆i is sufficiently bounded away from 0, a fact which will be convenient for
us later on.
Our proof will revolve around the L2 function
g̃(x) =

N
X

i gi (x),

i=1

which as explained in Sec. 2, will be shown to be easy to approximate arbitrarily well with a 3-layer network,
but hard to approximate with a 2-layer network.

4.2

Key Lemmas

In this subsection, we collect several key technical lemmas on gi and ϕ, which are crucial for the main proof.
The proofs of all the lemmas can be found in Appendix B.
The following lemma ensures that ϕ(x) is sufficiently close to being a constant on any good interval:
Lemma 3. If d ≥ 2, α ≥ c and N ≥ cα3/2 d2 (for some sufficiently large universal constant c), then inside
any good interval ∆i , ϕ(x) has the same sign, and
supx∈∆i |ϕ(x)|
≤ 1 + d−1/2 .
inf x∈∆i |ϕ(x)|
The following lemma ensures that the Fourier transform ĝi of gi has a sufficiently large part of its L2
mass far away from the origin:
Lemma 4. Suppose N ≥ 100αd3/2 . Then for any i,
Z
Z
1
2
gˆi (w)dw ≥
gˆi 2 (w)dw,
2 Rd
(2Rd Bd )C
where gˆi is the Fourier transform of gi .
The following lemma ensures that gc
i ϕ also has sufficiently large L2 mass far away from the origin:
Lemma 5. Suppose that α ≥ C, N ≥ Cα3/2 d2 and d > C, where C > 0 is a universal constant. Then for
any i,
Z
Z
1
2
[
((gi ϕ)(w)) dw ≥
(ϕ(x)gi (x))2 dx.
4
C
d
(2Rd Bd )
R
The following lemma ensures that a linear combination of the gi ’s has at least a constant L2 (ϕ2 ) probability mass.
Lemma 6. Suppose that α ≥ c and N ≥ c(αd)3/2 for some sufficiently large universal constant c, then for
every choice of i ∈ {−1, +1}, i = 1, . . . , N , one has
!2
Z X
N
0.003
.
i gi (x) ϕ2 (x)dx ≥
α
i=1

10

P
can be approxFinally, the following lemma guarantees that the non-Lipschitz function N
i=1 i gi (x)
PN
2
imated by a Lipschitz function (w.r.t. the density ϕ ). This will be used to show that i=1 i gi (x) can
indeed be approximated by a 3-layer network.
Lemma 7. Suppose that d ≥
√2. For√any choice of i ∈ {−1, +1}, i = 1, . . . , N , there exists an N -Lipschitz
function f , supported on [α d, 2α d] and with range in [−1, +1], which satisfies
Z
f (x) −

N
X

!2
ϕ2 (x)dx ≤

i gi (x)

i=1

4.3

3
√ .
d

α2

Inapproximability of the Function g̃ with 2-Layer Networks

The goal of this section is to prove the following proposition.
Proposition 1. Fix a dimension d, suppose that d > C, α > C and N ≥ Cα3/2 d2 and let k be an integer
satisfying
k ≤ cecd
(7)
with c, C > 0 being P
universal constants. There exists a choice of i ∈ {−1, +1}, i = 1, . . . , N , such that
d
the function g̃(x) = N
i=1 i gi (kxk) has the following property. Let f : R → R be of the form
f (x) =

k
X

fi (hx, vi i)

(8)

i=1

for vi ∈ Sd−1 , where fi : R → R are measurable functions satisfying
|fi (x)| ≤ C 0 (1 + |x|κ )
for constants C 0 , κ > 0. Then one has
kf − g̃kL2 (µ) ≥ δ/α
where δ > 0 is a universal constant.
The proof of this proposition requires a few intermediate steps. In the remainder of the section, we will
assume that N, d, α are chosen to be large enough to satisfy the assumptions of Lemma 6 and Lemma 5. In
other words we assume that d > C and N ≥ Cα3/2 d2 for a suitable universal constant C > 0. We begin
with the following:
Lemma 8. Suppose that d, N are as above. There exists a choice of i ∈ {−1, 1}, i = 1, . . . , N such that
2
\ !
X

i gi ϕ (w) dw ≥ c

Z
(2Rd Bd )C

i

for a universal constant c > 0.

11

Proof. Suppose that each i is chosen independently and uniformly at random from {−1, +1}. It suffices to
show that



2
!
Z
\
X



E
i gi ϕ (w) dw ≥ c
(2Rd Bd )C

i

for some universal constant c > 0, sincethat would ensure there exist some choice of 1 , . . . , N satisfying
the lemma statement. Define h(w) = 1 w ∈ (2Rd Bd )C and consider the operator
c
P (g) = ĝh.
This is equivalent to removing low-frequency components from g (in the Fourier domain), and therefore is
an orthogonal projection. According to Lemma 5 and isometry of the Fourier transform, we have
kP (gi ϕ)k2L2 ≥

1
kgi k2L2 (µ)
4

(9)

for every good i. Moreover, an application of Lemma 6, and the fact that hgi , gj iL2 = 0 for any i 6= j (as
gi , gj have disjoint supports) tells us that
N
X
i=1

kgi k2L2 (µ)

=

N
X

2

≥c

gi

i=1

for a universal constant c > 0. We finally get,



2

!
\
Z
Z
N
X



E
i gi ϕ (w) dw = E 
(2Rd Bd )C

(10)

L2 (µ)

N
X

Rd

i=1

=E

N
X

!
i P (gi ϕ) (x)

=


dx

i=1
2

i P (gi ϕ)

i=1
N
X

!2

L2

E[i j ]hP (gi ϕ), P (gj ϕ)iL2

i,j=1

=

N
X

kP (gi ϕ)k2L2

i=1
Eq. (9)

≥

N
1 X
kgi k2L2 (µ)
4
i,j=1

Eq. (10)

≥

c/4.

for a universal constant c > 0.
Claim 2. Let f be a function such that f ϕ ∈ L2 , and is of the form in Eq. (8). Suppose that the functions fi
are measurable functions satisfying
|fi (x)| ≤ C(1 + |x|α )
(11)
12

for constants C, α > 0. Then,
Supp(fc
ϕ) ⊂

k
[

(Span{vi } + Rd Bd )

(12)

i=1

Proof. Informally, the proof is based onP
the convolution-multiplication and linearity principles of the Fourier
transform, which imply that if f =
i fi , where each fi as well as ϕ have a Fourier transform, then
P
P ˆ
fc
ϕ = i fd
ϕ
=
f
?
ϕ̂.
Roughly
speaking,
in our case each fˆi (x) = fˆi (hx, vi i) (as a function in Rd )
i
i i
is shown to be supported on Span{vj }, so its convolution with ϕ̂ (which is an indicator for the ball Rd Bd )
must be supported on Span{vi } + Rd Bd . Summing over i gives the stated result.
Unfortunately, this simple analysis is not formally true, since we are not guaranteed that fi has a Fourier
transform as a function in L2 (this corresponds to situations where the integral in the definition of the Fourier
transform in Eq. (5) does not converge). However, at least for functions satisfying the claim’s conditions,
the Fourier transform still exists as a generalized function, or tempered distribution, over Rd , and using this
object we can attain the desired result.
We now turn to provide the formal proof and constructions, starting with a description of tempered
distributions and their relevant properties (see [Hunter and Nachtergaele, 2001, Chapter 11] for a more
complete survey). To start, let S denote the space of Schwartz functions 6 on Rd . A tempered distribution
µ in our context is a continuous linear operator from S to R (this can also be viewed as an element in the
dual space S ∗ ). In particular, any measurable function h : Rd 7→ R, which satisfies a polynomial growth
condition similar to Eq. (11), can be viewed as a tempered distribution defined as
Z
ψ 7→ hh, ψi :=
h(x)ψ(x)dx,
Rd

where ψ ∈ S. Note that the growth condition ensures that the integral above is well-defined. The Fourier
transform ĥ of a tempered distribution h is also a tempered distribution, and defined as
hĥ, ψi := hh, ψ̂i,
where ψ̂ is the Fourier transform of ψ. It can be shown that this directly generalizes the standard notion of
Fourier transforms of functions. Finally, we say that a tempered distribution h is supported on some subset
of Rd , if hh, ψi = 0 for any function ψ ∈ S which vanishes on that subset.
With these preliminaries out of the way, we turn to the setting considered in the claim. Let fˆi denote
the Fourier transform of fi (possibly as a tempered distribution, as described above), the existence of which
is guaranteed by the fact that fi is measurable and by Eq. (11). We also define, for ψ : Rd → R and
1 ≤ i ≤ N , a corresponding function ψi : R → R by
ψi (x) = ψ(xvi ),
and define the tempered distribution µi (over Schwartz functions in Rd ) as
hµi , ψi := hfi , ψ̂i i,
6

This corresponds 
to infinitely-differentiable functions ψ : Rd 7→R with super-polynomially decaying values and derivatives,
α

α1

α2
1
d ∂
in the sense that supx xα
1 · x2 · · · xd ∂xα1 ·
1

∂ α2
α
∂x2 2

···

∂ αd
α
∂xd d

f (x)

13

< ∞ for all indices α1 , . . . , αd .

which is indeed an element of S ∗ by the linearity of the Fourier transform, by the continuity of ψ → ψi (x)
with respect to the topology of S and by the dominated convergence theorem. Finally, define
f˜i (x) = fi (hx, vi i).
Using the fact that7
Z
ĝ(x)dx = g(0)

(13)

Rd

for any g ∈ S, recalling that vi has unit norm, and letting vi⊥ denote the subspace of vectors orthogonal to
vi , we have the following for any ψ ∈ S:
Z
Z Z
˜
˜
hfi , ψ̂i =
fi (x)ψ̂(x)dx =
f˜i (x + vi y)ψ̂(x + vi y)dxdy
(14)
Rd

fi (y)

=
R

Z

vi⊥

Z

=

fi (y)
R

Z

vi⊥

fi (y)
R

Z

ψ̂(x + vi y)dxdy
Z

vi⊥

ψ(w) exp(−2πihw, x + vi yi)dwdxdy
Rd

Z Z

Z

=

vi⊥

R

vi⊥

ψ(w1 + vi w2 ) exp(−2πihw1 + vi w2 , x + vi yi)dw1 dw2 dxdy

Z

=
(13)

R

Z

Z

Z

Z

fi (y)

exp(−2πiw2 y)

R

R

Z

Z

vi⊥

vi⊥

ψ(w1 + vi w2 ) exp(−2πihw1 , xi)dw1 dxdw2 dy

=
fi (y) exp(−2πiw2 y)ψ(vi w2 )dw2 dy
R
R
Z
Z
=
fi (y)ψ̂(vi y)dy =
fi (y)ψ̂i (y)dy = hfi , ψ̂i i = hµi , ψi,
R

R

where the use of Fubini’s theorem is justified by the fact that ψ ∈ S.
We now use the convolution-multiplication theorem (see e.g., [Hunter and Nachtergaele, 2001, Theorem
11.35]) according to which if f, g ∈ L1 then
f[
? g = fˆĝ.

(15)

Using this, we have the following for every ψ ∈ S:
d̃
hfi ϕ, ψi = hf˜i ϕ, ψ̂i = hf˜i , ϕψ̂i
(15)

= hf˜i , ϕ̂[
? ψi

(14)

= hµi , ϕ̂ ? ψi = hµi , 1 {Rd Bd } ? ψi.

d̃
Based on this equation, we now claim that hfi ϕ, ψi = 0 for any ψ ∈ S supported on the complement of
d̃
Span{vi } + Rd Bd . This would imply that the tempered distribution fi ϕ is supported in Span{vi } + Rd Bd ,
7


R
RR
R
R
ĝ(x)dx =
g(x) exp(−2πihx, wi)dwdx =
g(x) exp(−2πihx, wi) · 1dw dx =
This is because
g(x)δ(x)dx = g(0), where δ(·) is the Dirac delta function, which is the Fourier transform of the constant 1 function. See
also [Hunter and Nachtergaele, 2001, Chapter 11, Example 11.31].
R

14

S
and therefore fc
ϕ is supported in ki=1 (Span{vi } + Rd Bd ) (by linearity of the Fourier transform and the
P
fact that f = ki=1 f˜i ). Since the Fourier transform of f ϕ as a tempered distribution coincides with the
standard one (as we assume f ϕ ∈ L2 ), the result follows.
It remains to prove that fd
i ϕ(ψ) = 0 for any ψ ∈ S supported on the complement of Span{vi } + Rd Bd .
For such ψ, by definition of a convolution, 1 {Rd Bd } ? ψ is supported on the complement of Span{vi }.
However, µi is supported on Span{vi } (since if ψ vanishes on vi , then ψi is the zero function, hence ψ̂i is
also the zero function, and hµi , ψi = hfi , ψ̂i i = 0). Therefore, hµi , 1 {Rd Bd } ? ψi = 0, which by the last
d̃
displayed equation, implies that hfi ϕ, ψi = 0 as required.
Lemma 9. Let q, w be two functions of unit norm in L2 . Suppose that q satisfies
Supp(q) ⊂

k
[

(Span{vj } + Rd Bd )

(16)

j=1

R
for some k ∈ N. Moreover, suppose that w is radial and that 2Rd Bd w(x)2 dx ≤ 1 − δ for some δ ∈ [0, 1].
Then
hq, wiL2 ≤ 1 − δ/2 + k exp(−cd)
where c > 0 is a universal constant.
Proof. Define A = (2Rd Bd )C and denote
T =

k
[

(Span{vj } + Rd Bd )

j=1

so that T contains the support of q(·). For each r > 0, define
h(r) =

V old−1 (rSd−1 ∩ T )
,
V old−1 (rSd−1 )

where Sd−1 is the Euclidean sphere in Rd . Since T is star shaped, the function h(r) is non-increasing. We
claim that there exists a universal constant c > 0 such that
h(2Rd ) ≤ k exp(−cd).

(17)

Indeed, fix v ∈ Sd−1 and consider the tube T0 = Span{v} + Rd Bd . Let z be a uniform random point in
2Rd Sd−1 . We have by a standard calculation (See e.g., [Sodin, 2007, Section 2])
Pr(z ∈ T0 ) = Pr(kzk2 − hz, vi2 ≤ Rd2 )
= Pr(4Rd2 − hz, vi2 ≤ Rd2 ) = Pr(|hz, vi| ≥
R1
2 (d−3)/2
√
3/2 (1 − t )
= R1
dt ≤ exp(−cd).
2 (d−3)/2
0 (1 − t )
Using a union bound and the definition of h, equation Eq. (17) follows.
Next, define
R
kxkSd−1 q(y)dHd−1 (y)
q̃(x) =
V old−1 (kxkSd−1 )
15

√

3Rd )

to be the averaging of q(·) with respect to rotations (in the above formula Hd−1 denotes the d−1 dimensional
Hausdorff measure, i.e. the standard measure in d − 1 dimensions). We have the following: Since w(·) is
radial and has unit L2 norm, and we assume q(·) is supported on T , we have
Z
Z
(1)
w(x)q̃(x)dx
w(x)q(x)dx =
A

A
(2)

≤ kwkL2 kq̃1 {A}kL2
sZ
∞
(3)
=
q̃(r)2 V old−1 (rSd−1 )dr
2Rd

sZ

2
Z
1
=
V old−1
q(y)dHd−1 (y) dr
V old−1 (rSd−1 ) rSd−1 ∩T
2Rd
sZ

2
Z
∞
1
2
d−1
=
h(r) V old−1 (rS )
q(y)dHd−1 (y) dr
V old−1 (rSd−1 ∩ T ) rSd−1 ∩T
2Rd
sZ


Z
∞
(4)
1
2
2
d−1
q(y) dHd−1 (y) dr
≤
h(r) V old−1 (rS )
V old−1 (rSd−1 ∩ T ) rSd−1 ∩T
2Rd
sZ
Z
∞
=
h(r)
q 2 (y)dHd−1 (y)dr
∞

2Rd
(5)

≤



(rSd−1 )

rSd−1

(6)
p
h(2Rd ) kq1 {A}kL2 ≤ k exp(−cd/2).

In the above, (1) follows from w(·) being radial; (2) follows from Cauchy-Schwartz; (3) follows from w(·)
having unit L2 norm; (4) follows from the fact that the term being squared is the expectation of q(y) where
y is uniformly distributed in rSd−1 ∩ T , and we have (Ey [q(y)])2 ≤ Ey [q 2 (y)] by Jensen’s inequality; (5)
follows from h(·) being non-increasing; and (6) follows from Eq. (17) and the fact that q(·) has unit L2
norm.
As a result of these calculations, we have
Z
Z
hq, wiL2 =
w(x)q(x)dx +
w(x)q(x)dx
A
AC
sZ
 C
= k exp(−cd/2) + 1 ·
≤ k exp(−cd/2) + kqkL2 w1 A
w2 (x)dx
L2
≤ k exp(−cd/2) +

√

AC

1 − δ.

where we used the assumption that q(·) is unit norm and that
√
Since 1 − δ ≤ 1 − 12 δ for any δ ∈ [0, 1], the result follows.

R

AC

w2 (x)dx =

R

(2Rd Bd )C

w2 (x)dx ≤ 1 − δ.

Proof of Proposition 1. Define
g̃(x) =

X

i gi (|x|)

i

where (i ) are the signs provided by Lemma 8. According to Lemma 6, we have
kg̃kL2 (µ) ≥ c1 /α
16

(18)

for a universal constant c1 > 0. Note that the function g̃ is bounded between −1 and 1. Define the function
c̃
gϕ
. By construction (hence according to Lemma 8) we have
w = kg̃ϕk
L2

R

Z

w(x)2 dx = 1 −

2Rd Bd

kg̃ϕk2L2

2Rd Bd

R
≤1−

c̃ 2 dx
gϕ(x)

2Rd Bd

c̃ 2 dx
gϕ(x)

kϕk2L2

≤ 1 − c2 ,

for a universal constant c2 > 0, where in the first inequality we used the fact that |g̃(x)| ≤ 1 for all x.
Next, define the function q =

fc
ϕ
kf ϕkL

, where f is an arbitrary function having the form in Eq. (8).
2

Thanks to the assumptions on the functions fi , we may invoke8 Claim 2, by which we observe that the
functions w, q satisfy the assumptions of Lemma 9. Thus, as a consequence of this lemma we obtain that
hq, wiL2 ≤ 1 − c2 /2 + k exp(−c3 d)

(19)

for a universal constant c3 > 0. Next, we claim that since kqkL2 = kwkL2 = 1, we have for every scalars
β1 , β2 > 0 that
β2
kβ1 q − β2 wkL2 ≥
kq − wkL2 .
(20)
2
Indeed, we may clearly multiply both β1 and β2 by the same constant affecting the correctness of the
formula, thus we may assume that β2 = 1. It thus amounts to show that for two unit vectors v, u in a Hilbert
space, one has that minβ>0 kβv − uk2 ≥ 14 kv − uk2 . We have

min kβv − uk2 = min β 2 kvk2 − 2βhv, ui + kuk2
β
β

= min β 2 − 2βhv, ui + 1
β

1
= 1 − hv, ui2 = kv − uk2
2
which in particular implies formula Eq. (20).
Combining the above, and using the fact that q, w have unit L2 norm, we finally get


c̃
kf − g̃kL2 (µ) = kf ϕ − g̃ϕkL2 = fc
ϕ − gϕ
= kf ϕkL2 q(·) − kg̃ϕkL2 w(·)
L2

Eq. (20)

≥
Eq. (18)

≥
Eq. (19)

≥

L2

1
1
kq − wkL2 kg̃ϕkL2 = kq − wkL2 kg̃kL2 (µ)
2
2
c1
1p
2(1 − hq, wiL2 )
2
α
√
c1 c2
c1 p
2 max(c2 /2 − k exp(−c3 d), 0) ≥
2α
4α

where in the last inequality, we use the assumption in Eq. (7), choosing c = min{c2 /4, c3 }. The proof is
complete.
8

Claim 2 also requires that f ϕ is an L2 function, but we can assume this without loss of generality: Otherwise, kf ϕkL2 = ∞,
and since g̃ϕ is an L2 function with kg̃ϕkL2 < ∞, we would have kf − g̃k2L2 (µ) = kf ϕ − g̃ϕk2L2 = ∞, in which case the
proposition we wish to prove is trivially true.

17

4.4

Approximability of the Function g̃ with 3-Layer Networks

The next ingredient
PN missing for our proof is the construction of a 3-layer function which approximates the
function g̃ = i=1 i gi .
Proposition 2. There is a universal constant C > 0 such that the following holds. Let δ ∈ (0, 1). Suppose
that d ≥ C and that the functions gi are constructed as in Eq. (6). For any choice of i ∈ {−1, +1},
i = 1, . . . , N , there exists a function g expressible by a 3-layer network of width at most 8cδσ α3/2 N d11/4 +1,
and with range in [−2, +2], such that
g(x) −

N
X

√

≤

i gi (kxk)

i=1

L2 (µ)

3
+ δ.
αd1/4

The proof of this proposition relies on assumption 1, which ensures that we can approximate univariate
functions using our activation function. As discussed before Thm. 1, one can also plug in weaker versions
of the assumption (i.e. worse polynomial dependence of the width w on R, L, 1/δ), and get versions of
proposition 2 where the width guarantee has worse polynomial dependence on the parameters N, α, d, δ.
This would lead to versions of the Thm. 1 with somewhat worse constants and polynomial dependence on
the dimension d.
For this proposition, we need a simple intermediate result, in which an approximation for radially symmetric Lipschitz functions in Rd , using assumption 1, is constructed.
Lemma 10. Suppose the activation function σ satisfies assumption 1. Let f be an L-Lipschitz function
supported on [r, R], where r ≥ 1. Then for any δ > 0, there exists a function g expressible by a 3-layer
2 2
network of width at most 2cσ√d rδR L + 1, such that
sup g(x) − f (kxk) < δ.
x∈Rd

Proof. Define the 2R-Lipschitz function
l(x) = min{x2 , R2 },
which is constant outside [−R, R], as well as the function
`(x) =

d
X

l(xi ) =

i=1

d
X

min{x2i , R2 }

i=1

P
on Rd . Applying assumption 1 on l(·), we can obtain a function ˜l(x) having the form a+ w
i=1 αi σ(βi x−γi )
so that
√
rδ
sup ˜l(x) − l(x) ≤
,
dL
x∈R
and where the width parameter w is at most

2cσ√dR2 L
.
rδ

˜
`(x)
=

Consequently, the function

d
X
i=1

18

˜l(xi )

can be expressed in the form a +

Pw

i=1 αi σ(βi x

− γi ) where w ≤

2cσ√
d2 R2 L
,
rδ

and it holds that

√
˜ − `(x) ≤
sup `(x)
x∈Rd

Next, define

rδ
.
L

(21)

( √
f ( x) x ≥ 0
s(x) =
.
0
x<0

L
-Lipschitz and
Since f is L-Lipschitz and supported on {x : r ≤ x ≤ R}, it follows that s is 2√
r
2
2
supported
Pw on the interval [−R , R ]. Invoking assumption 1 again, we can construct a function s̃(x) =
a + i=1 αi σ(βi x − γi ) satisfying
sup |s̃(x) − s(x)| < δ/2,
(22)
x∈R

where w ≤

R2 L

cσ√

rδ

.

˜ which by definition of s̃, `,
˜ has the form
Now, let us consider the composition g = s̃ ◦ `,


w
w
X
X
a+
ui σ 
vi,j σ (hwi,j , xi + bi,j ) + ci 
i=1

(23)

j=1

for appropriate scalars a, ui , ci , vi,j , bi,j and vectors wi,j , and where w is at most


2cσ d2 R2 L cσ R2 L
2cσ d2 R2 L
√
√
max
, √
=
.
rδ
rδ
rδ
Eq. (23) is exactly a 3-layer network (compare to Eq. (2)), except that there is an additional constant term a.
a
However, by increasing w by 1, we can simulate a by an additional neuron x 7→ σ(σ(0)+z)
·σ(σ(h0, xi)+z),
where z is some scalar such that σ(σ(0) + z) 6= 0 (note that if there is no such z, then σ is the zero function,
and therefore cannot satisfy assumption 1). So, we can write the function g as a 3-layer network (as defined
in Eq. (2)), of width at most
2cσ d2 R2 L
√
+ 1.
rδ
All the remains now is to prove that supx∈Rd |g(x) − f (kxk)| ≤ δ. To do so, we note that for any x ∈ Rd ,
we have
˜
|g(x) − f (kxk)| = s̃(`(x))
− f (kxk)
˜
˜
˜ − s(`(x)) + |s(`(x)) − f (kxk)|
≤ s̃(`(x))
− s(`(x)
+ s(`(x)
p
˜
˜
˜ − s(`(x)) + f ( `(x)) − f (kxk) .
= s̃(`(x))
− s(`(x)
+ s(`(x)
Let us consider each of the three absolute values:
• The first absolute value term is at most δ/2 by Eq. (22).
• The second absolute value term, since s is
most δ/2 by Eq. (21).

L
√
-Lipschitz,
2 r

19

is at most

L ˜
√
|`(x)
2 r

− `(x)|, which is at

• As to the third absolute value term, if kxk2 ≤ R2 , then `(x) = kxk2 and the term is zero. If
kxk2 > R2 , then it is p
easy to verify that `(x) ≥ R2 , and since f is continuous and supported on
[r, R], it follows that f ( `(x) = f (kxk) = 0, and again, we get a zero.
Summing the above, we get that |g(x) − f (kxk)| ≤

δ
2

+

δ
2

+ 0 = δ as required.

We are now ready to prove Proposition 2, which is essentially a combination of Lemmas 7 and 10.
Proof of Proposition 2. First, invoke Lemma 7 to obtain an N -Lipschitz function h with range in [−1, +1]
which satisfies
v
!2
uZ
√
N
N
u
X
X
3
t
2
˜
f (x) −
h(x) −
i gi (x) ϕ (x)dx ≤
.
(24)
i gi (x)
=
1/4
αd
Rd
i=1
i=1
L2 (µ)

√
√
Next, we use Lemma 10 with R = 2α d, r = α d, L = N to construct a function g expressible by a
3-layer network of width at most 8cδσ α3/2 N d11/4 + 1, satisfying supx∈Rd |g(x) − h(x)| ≤ δ. This implies
that kg − hkL2 (µ) ≤ δ, and moreover, that the range of g is in [−1 − δ, 1 + δ] ⊆ [−2, +2] (since we assume
δ < 1). Combining with Eq. (24) and using the triangle inequality finishes the proof.

4.5

Finishing the Proof

We are finally ready to prove our main theorem.
Proof of Theorem 1. The proof is a straightforward combination of propositions 1 and 2 (whose conditions
can be verified to follow immediately from the assumptions used in the theorem). We first choose α = C
and N = dCα3/2 d2 e with the constant C taken from the statement of Proposition 1. By invoking this
proposition we obtain signs i ∈ {−1, 1} and a universal constant δ1 > 0 for which any function f expressed
by a bounded-size 2-layer network satisfies
kg̃ − f kL2 (µ) ≥ δ1 ,

(25)

P
where g̃(x) = N
i=1 i gi (kxk). Next, we use Proposition 2 with δ = min{δ1 /2, 1} to approximate g̃ by a
function g expressible by a 3-layer network of width at most
16cσ 3/2 11/4
16cσ 3/2 5/2 2 11/4
α Nd
+1 =
C dC d ed
+ 1 ≤ C 0 cσ d19/4 ,
δ
δ
(where C 0 is a universal constant depending on the universal constants C, δ1 ), so that
kg̃ − gkL2 (µ) ≤ δ ≤ δ1 /2.

(26)

Combining Eq. (25) and Eq. (26) with the triangle inequality, we have that kf − gkL2 (µ) ≥ δ1 /2 for any
2-layer function f . The proof is complete.
Acknowledgements
OS is supported in part by an FP7 Marie Curie CIG grant, the Intel ICRI-CI Institute, and Israel Science
Foundation grant 425/13. We thank James Martens and the anonymous COLT 2016 reviewers for several
helpful comments.
20

References
Andrew R. Barron. Approximation and estimation bounds for artificial neural networks. Machine Learning,
14(1):115–133, 1994.
M. Bianchini and F. Scarselli. On the complexity of shallow and deep neural network classifiers. In ESANN,
2014.
N. Cohen, O. Sharir, and A. Shashua. On the expressive power of deep learning: A tensor analysis. arXiv
preprint arXiv:1509.05009, 2015.
G. Cybenko. Approximation by superpositions of a sigmoidal function. Mathematics of control, signals and
systems, 2(4):303–314, 1989.
C. Debao. Degree of approximation by superpositions of a sigmoidal function. Approximation Theory and
its Applications, 9(3):17–28, 1993.
O. Delalleau and Y. Bengio. Shallow vs. deep sum-product networks. In NIPS, pages 666–674, 2011.
DLMF. NIST Digital Library of Mathematical Functions. http://dlmf.nist.gov/, Release 1.0.10 of 2015-0807, 2015. URL http://dlmf.nist.gov/.
Ken-Ichi Funahashi. On the approximate realization of continuous mappings by neural networks. Neural
networks, 2(3):183–192, 1989.
Loukas Grafakos and Gerald Teschl. On fourier transforms of radial functions and distributions. Journal of
Fourier Analysis and Applications, 19(1):167–179, 2013.
András Hajnal, Wolfgang Maass, Pavel Pudlák, Márió Szegedy, and György Turán. Threshold circuits of
bounded depth. Journal of Computer and System Sciences, 46(2):129–154, 1993.
J. Håstad. Almost optimal lower bounds for small depth circuits. In Proceedings of the eighteenth annual
ACM symposium on Theory of computing, pages 6–20. ACM, 1986.
K. Hornik, M. Stinchcombe, and H. White. Multilayer feedforward networks are universal approximators.
Neural networks, 2(5):359–366, 1989.
John K. Hunter and Bruno Nachtergaele. Applied analysis. World Scientific Publishing Co., Inc., River
Edge, NJ, 2001.
Ilia Krasikov. Approximations for the bessel and airy functions with an explicit error term. LMS Journal of
Computation and Mathematics, 17(01):209–225, 2014.
W. Maass, G. Schnitger, and E. Sontag. A comparison of the computational power of sigmoid and boolean
threshold circuits. In V. P. Roychowdhury, K. Y. Siu, and A. Orlitsky, editors, Theoretical Advances in
Neural Computation and Learning, pages 127–151. Springer, 1994.
J. Martens and V. Medabalimi. On the expressive efficiency of sum product networks. arXiv preprint
arXiv:1411.7717, 2014.
James Martens. Private Communication, 2015.
21

James Martens, Arkadev Chattopadhya, Toni Pitassi, and Richard Zemel. On the representational efficiency
of restricted boltzmann machines. In NIPS, pages 2877–2885, 2013.
G. F Montufar, R. Pascanu, K. Cho, and Y. Bengio. On the number of linear regions of deep neural networks.
In Advances in Neural Information Processing Systems, pages 2924–2932, 2014.
I. Parberry. Circuit complexity and neural networks. MIT press, 1994.
R. Pascanu, G. Montufar, and Y. Bengio. On the number of inference regions of deep feed forward networks
with piece-wise linear activations. arXiv preprint arXiv, 1312, 2013.
B. Rossman, R. Servedio, and L.-Y. Tan. An average-case depth hierarchy theorem for boolean circuits. In
FOCS, 2015.
D. E. Rumelhart, G. E. Hinton, and R. J. Williams. Learning internal representations by error propagation.
In D. E. Rumelhart and J.L. McClelland, editors, Parallel distributed Processing, volume 1. MIT Press,
1986.
A. Shpilka and A. Yehudayoff. Arithmetic circuits: A survey of recent results and open questions. Foundations and Trends R in Theoretical Computer Science, 5(3–4):207–388, 2010.
Sasha Sodin. Tail-Sensitive Gaussian Asymptotics for Marginals of Concentrated Measures in High Dimension. In Vitali D. Milman and Gideon Schechtman, editors, Geometric Aspects of Functional Analysis,
volume 1910 of Lecture Notes in Mathematics, pages 271–295. Springer Berlin Heidelberg, 2007.
M. Telgarsky. Representation benefits of deep feedforward networks. arXiv preprint arXiv:1509.08101,
2015.

22

A

Approximation Properties of the ReLU Activation Function

In this appendix, we prove that the ReLU activation function satisfies assumption 1, and also prove bounds
on the Lipschitz parameter of the approximation and the size of the required parameters. Specifically, we
have the following lemma:
Lemma 11. Let σ(z) = max{0, z} be the ReLU activation function, and fix L, δ, R > 0. Let f : R → R
RL
which is constant outside an interval [−R, R]. There exist scalars a, {αi , βi }w
i=1 , where w ≤ 3 δ , such that
the function
w
X
h(x) = a +
αi σ(x − βi )
(27)
i=1

is L-Lipschitz and satisfies
sup f (x) − h(x) ≤ δ.

(28)

x∈R

Moreover, one has |αi | ≤ 2L and w ≤ 3 RL
δ .
Proof. If one has 2RL < δ, then the results holds trivially because we may take the function h to be the 0
function (with width parameter w = 0). Otherwise, we must have R ≥ δ/2L, so by increasing the value of
R by a factor of at most 2, we may assume without loss of generality that there exists an integer m such that
R = mδ/L.
Let h be the unique piecewise linear function which coalesces with f on points of the form δ/Li,
i ∈ Z ∩ [−m, m], is linear in the intervals (wδ/L, (w + 1)δ/L) and is constant outside [−R, R]. Since
f is L-Lipschitz, equation Eq. (28) holds true. It thus suffices to express h as a function having the form
Eq. (27). Let βi = iδ/L, choose a = h(−R) and set
αi = h0 (βi +

δ
2L )

− h0 (βi −

δ
2L ),

− m ≤ i ≤ m.

Then clearly equation Eq. (27) holds true. Moreover, we have |αi | ≤ 2L, which completes the proof.

B
B.1

Technical Proofs
Proof of Lemma 2

By definition of the Fourier transform,
Z

exp(−2πix> w)dw.

ϕ(x) =
w:kwk≤Rd

Since ϕ(x) is radial (hence rotationally invariant), let us assume without loss of generality that it equals re1 ,
where r = kxk and e1 is the first standard basis vector. This means that the integral becomes
!
Z
Z
Z
Rd

exp(−2πirw1 )dw =
w:kwk≤Rd

exp(−2πirw1 )
w1 =−Rd

w2 ...wd :

23

Pd

j=2

wj2 ≤Rd2 −w12

dw2 . . . dwd

dw1 .

The expression inside the parenthesis is simply the volume of a ball of radius Rd2 − w12
Vd−1 be the volume of a unit ball in Rd−1 , this equals
Z

Rd

1/2

in Rd−1 . Letting



d−1
exp(−2πirw1 ) Vd−1 (Rd2 − w12 ) 2 dw1 .

w1 =−Rd

Performing the variable change z = arccos(w1 /Rd ) (which implies that as w1 goes from −Rd to Rd , z
goes from π to 0, and also Rd cos(z) = w1 and −Rd sin(z)dz = dw1 ), we can rewrite the integral above as
Z π
 d−1
Vd−1
Rd2 − Rd2 cos2 (z) 2 exp(−2πirRd cos(z))Rd sin(z)dz
z=0
Z π
d
sind (z) exp (−2πirRd cos(z)) dz.
= Vd−1 Rd
z=0

Since we know that this integral must be real-valued (since we’re computing the Fourier transform ϕ(x),
which is real-valued and even), we can ignore the imaginary components, so the above reduces to
Z π
d
sind (z) cos (2πrRd cos(z)) dz.
(29)
Vd−1 Rd
z=0

By a standard formula for Bessel functions (see Equation 10.9.4. in DLMF), we have
Jd/2 (x) =

(x/2)d/2

π 1/2 Γ d+1
2

Z

π

sind (z) cos(x cos(z))dz,

0

which by substituting x = 2πrRd and changing sides, implies that
Z
0

π


π 1/2 Γ d+1
2
sin (z) cos(2πrRd cos(z))dz =
Jd/2 (2πrRd ).
(πrRd )d/2
d

Plugging this back into Eq. (29), we get the expression
d/2 π
Vd−1 Rd

Plugging in the explicit formula Vd−1 =

d+1
2
(πr)d/2

1/2 Γ

π (d−1)/2
,
Γ( d+1
2 )



Rd
r


Jd/2 (2πrRd ).

this simplifies to

d/2
Jd/2 (2πRd r).

Recalling that this equals ϕ(x) where kxk = r, the result follows.

B.2

Proof of Lemma 3

By Lemma 2,

ϕ(x) =

Rd
x

d/2
Jd/2 (2πRd x).

24

Moreover,
√ using the definition of a good interval, and the fact that the maximal value in any interval is at
most 2α d, we have
1
1
(30)
|Jd/2 (2πRd x)| ≥ √
≥p
√ .
80πRd x
160πRd α d
√
Since x (in any interval) is at least α d, then Jd/2 (2πRd x) is 2πRd -Lipschitz in x by Lemma 12. Since
√
α d
N ,

Eq. (30) implies that Jd/2 (2πRd x) (and hence ϕ(x)) does not change

√
√ 3/2
√
signs in the interval, provided that N > 2 160 παRd d
. Recalling that Rd ≤ 12 d, this is indeed
satisfied by the lemma’s conditions.
Turning to the second part of the lemma, assuming ϕ(x) is positive without loss of generality, and using
the Lipschitz property of Jd/2 (·) and Eq. (30), we have
the width of each interval only

supx∈∆i ϕ(x)
≤
inf x∈∆i ϕ(x)

supx∈∆i



inf x∈∆i



Rd
x

d/2

supx∈∆i Jd/2 (2πRd x)
d/2 · inf
x∈∆i Jd/2 (2πRd x)
Rd
x
√

dα d
inf x∈∆i Jd/2 (2πRd x) + 2πRN
≤
·
inf x∈∆i Jd/2 (2πRd x)
√ !d/2
!
√ q
√
inf x∈∆i x + αN d
2πRd α d
1+
≤
80πRd α d
inf x∈∆i x
N
!
√ !d/2
√
√
α d
2 80(παRd d)3/2
√
≤ 1+
1+
N
Nα d
!
√

d/2
1
2 80(παd/2)3/2
≤ 1+
1+
,
N
N



supx∈∆i x
inf x∈∆i x

d/2

which is less than 1 + d−1/2 provided that N ≥ cα3/2 d2 for some universal constant c.

B.3

Proof of Lemma 4

The result is trivially true for a bad interval i (where gi is the 0 function, hence both sides of the inequality
in the lemma statement are 0), so we will focus on the case that i is a good interval.
√
√
For simplicity, let us denote the interval ∆i as [`, ` + δ], where δ = N1 and ` is between α d and 2α d.
1
Therefore, the conditions in the lemma imply that δ ≤ 50d`
. Also, we drop the i subscript and refer to gi as
g.
Since, g is a radial function, its Fourier transform is also radial, and is given by
Z

∞

ĝ(w) = ĝ(kwk) = 2π


g(s)

s=0

s
kwk

d/2−1
Jd/2−1 (2πskwk)s ds,

(see for instance Grafakos and Teschl [2013], section 2, and references therein). Using this formula, and
switching to polar coordinates (letting Ad denote the surface area of a unit sphere in Rd ), we have the

25

following:
Z

2Rd

Z

2

ĝ (w)dw =
2Rd Bd
Z 2Rd

=

 Z
Ad rd−1 2π
Z

2Rd

g(s)

r=0

r

2
Jd/2−1 (2πsr)s ds

dr

s=0

2Rd

s

r=0

2

`+δ

Z
r

= 4π Ad

 s d/2−1

2
g(s)sd/2 Jd/2−1 (2πsr) ds dr

Z
r

Z

2

r=0
∞

s=0
∞

r=0

= 4π 2 Ad

Ad rd−1 ĝ 2 (r)dr

d/2

Jd/2−1 (2πsr) ds

dr.

(31)

s=`

By Lemma 12, |Jd/2−1 (x)| ≤ 1, hence Eq. (31) can be upper bounded by
2

2Rd

Z

Z

d/2

r

4π Ad

s

r=0

2

`+δ

ds

dr ≤ 4π 2 Ad

2Rd


2
r δ(` + δ)d/2 dr

r=0

s=`

≤ 4π 2 Ad δ 2 (` + δ)d

Z

Z

2Rd

r dr = 8π 2 Ad δ 2 (` + δ)d Rd2 .

r=0

Overall, we showed that

Z

ĝ 2 (w)dw ≤ 8π 2 Rd2 Ad δ 2 (` + δ)d .

(32)

2Rd Bd

R
Let us now turn to consider Rĝ 2 (w)dw, where the integration is over all of w ∈ Rd . By isometry of the
Fourier transform, this equals g 2 (x)dx, so
Z

Z

ĝ 2 (w)dw =

Z

g 2 (x)dx =

Rd

∞

Ad rd−1 g 2 (r)dr =

r=0

Z

`+δ

Ad rd−1 dr ≥ Ad δ`d−1 .

r=`

Combining this with Eq. (32), we get that
R


2
8π 2 Rd2 Ad δ 2 (` + δ)d
δ d
2Rd Bd ĝ (w)dw
2
2
R
= 8π Rd `δ 1 +
.
≤
2
`
Ad δ`d−1
Rd ĝ (w)dw
Since we assume δ ≤
above is at most

1
50d` ,

and it holds that 1 +


1 d
50d

≤ exp(1/50) and Rd ≤

1
2

√

d by Lemma 1, the



1 d
1
1
2π d`δ 1 +
≤ 2π 2 d`δ exp(1/50) ≤ 2π 2 exp(1/50) < .
50d
50
2
2

R

Overall, we showed that
R

2
2Rd Bd ĝ (w)dw
R
ĝ 2 (w)dw
Rd

2
(2Rd Bd )C ĝ (w)dw
R
2
Rd ĝ (w)dw

≤ 21 , and therefore
R

=

Rd

R
ĝ 2 (w)dw − 2Rd Bd ĝ 2 (w)dw
1
1
R
≥ 1− =
2
2
2
Rd ĝ (w)dw

as required.

26

B.4

Proof of Lemma 5

The result is trivially true for a bad interval i (where gi is the 0 function, hence both sides of the inequality
in the lemma statement are 0), so we will focus on the case that i is a good interval.
Define a = supx∈∆i ϕ(x). Using Lemma 3, we have that ϕ(x) does not change signs in the interval ∆i .
Suppose without loss of generality that it is positive. Moreover, by the same lemma we have that
|ϕ(x) − a| ≤ d−1/2 a, ∀x ∈ ∆i
Consequently, we have that
Z

Z

(((ϕ\
− a)gi )(w))2 dw ≤

(((ϕ\
− a)gi )(w))2 dw

(33)

Rd

(2Rd Bd )C

Z

((ϕ − a)gi (x))2 dx
Z
−1
(agi (x))2 dx.
≤d

=

Rd

Rd

Next, by choosing the constant C to be large enough, we may apply Lemma 4, which yields that
Z
Z
1
2
[
((agi )(w)) dw ≥
(agi (x))2 dx.
2
C
d
(2Rd Bd )
R

(34)

By the triangle inequality, we have that for two vectors u, v in a normed space, one has kvk2 ≥ kuk2 −
2kvkkv − uk. This teaches us that
Z
Z
2
2
[
[
((gi ϕ)(w)) dw ≥
((ag
i )(w)) dw
(2Rd Bd )C

(2R Bd )C

sZd

2
[
((ag
i )(w)) dw

−2

sZ

(2Rd Bd )C
Eq. (33)

Z

2
−1/2
[
((ag
i )(w)) dw − 2d

≥
(2Rd Bd
Eq. (34)

≥

B.5

(((ϕ\
− a)gi )(w))2 dw

(2Rd Bd )C

Z

(agi (x))2 dx

Rd

)C

1
(1 − 4d−1/2 )
2

Z

1
(agi (x)) dx ≥
4
Rd
2

Z

(ϕ(x)gi (x))2 dx.

Rd

Proof of Lemma 6

Since the gi for different i have disjoint supports (up to measure-zero sets), the integral in the lemma equals
Z X
N
i=1

2

2

(i gi (x)) ϕ (x)dx =

Z X
N

gi2 (x)ϕ2 (x)dx

Z
=

ϕ2 (x)dx,

x:kxk∈good ∆i

i=1

where we used the definition of gi . Switching to polar coordinates (letting Ad be the surface area of the unit
sphere in Rd ), and using the definition of ϕ from Lemma 2, this equals
Z
Z
Rdd 2
Jd/2 (2πRd r)dr
Ad
rd−1 ϕ2 (r)dr = Ad
r∈good ∆i
r∈good ∆i r
27

Recalling that Ad =

dπ d/2
Γ( d2 +1)

and that Rdd = π −d/2 Γ

d
2


+ 1 by Lemma 1, this equals

2 (2πR r)
Jd/2
d

Z
d

r

r∈good ∆i

dr.

(35)

√
√
We now claim that for any r ∈ [α d, 2α d] (that is, in any interval),
2
(2πRd r) ≥
Jd/2

1
40πRd r

=⇒

r ∈ good ∆i ,

(36)

which would imply that we can lower bound Eq. (35) by
Z
d

√
2α d

√
α d


2 (2πR r) 
Jd/2
d
1
2
dr.
1 Jd/2 (2πRd r) ≥
r
40πRd r

(37)

To see √
why Eq. (36) holds, consider an r which satisfies the left hand side. The width of its interval is at
α d
most N , and by Lemma 12, Jd/2 (2πRd r) is at most 2πRd -Lipschitz in r. Therefore, for any other r0 in
the same interval as r, it holds that
√
r
1
2πRd α d
0
Jd/2 (2πRd r ) ≥
−
,
40πRd r
N
q
1
by the condition on N in the lemma statement, and the facts
which can be verified to be at least 80πR
dr
√
√
1
1
2
that r ≤ 2α d, Rd ≤ 2 d. As a result, Jd/2 (2πRd r0 ) ≥ 80πR
for any r0 in the same interval as r, which
dr
implies that r is in a good interval.
We now continue by taking Eq. (37), and performing the variable change x = 2πRd r, leading to
Z
d

√
4πRd α d

√
2πRd α d


2 (x) 
Jd/2
1
2
1 Jd/2
(x) ≥
dx.
x
20x

√
Applying Lemma 15 with β = 2πRd α/ d (which by Lemma 1, is between 2πα/5 and πα, hence satisfies
the conditions of Lemma 15 if α is large enough), this is at least
d

0.005
0.003
0.005
≥
≥
,
βd
2πα/5
α

from which the lemma follows.

B.6

Proof of Lemma 7

For any i, define
(
max{1, N dist(x, ∆C
i )}
ǧi (x) =
0

i good
i bad

where dist(x, ∆C
i ) is the distance of x from the boundaries of ∆i . Note that for bad i, this is the same as
gi (x), whereas for good i, it is an N -Lipschitz approximation of gi (x).

28

P
Let f (x) = N
i=1 ǧ(x), and note that since the support of ǧi are disjoint, f is also N Lipschitz. With
this definition, the integral in the lemma becomes
N
X

Z

!2
i (ǧi (x) − gi (x))

ϕ2 (x)dx.

i=1

Since the support of ǧi (x) − gi (x) is disjoint for different i, this equals
Z X
N

2

2

(ǧi (x) − gi (x)) ϕ (x)dx =

i=1

N Z
X

(ǧi (x) − gi (x))2 ϕ2 (x)dx.

i=1

Switching to polar coordinates (using Ad to denote the surface area of the unit sphere in Rd ), and using the
definition of ϕ from Lemma 2, this equals
N Z
X
i=1

∞

Ad r

d−1

2 2

(ǧi (r) − gi (r)) ϕ (r)dr =

0

N Z
X
i=1

0

∞

Ad

Rdd
2
(ǧi (r) − gi (r))2 Jd/2
(2πRd r)dr.
r

Using the definition of Rd from Lemma 1, and the fact that Ad =
N Z
X
i=1

∞

0

dπ d/2
,
Γ( d2 +1)

this equals

d
2
(ǧi (r) − gi (r))2 Jd/2
(2πRd r)dr.
r

Now, note that by definition of gˇi , gi , their difference |ǧi (r) − gi (r)| can be non-zero (and at √
most 1)√only
for r belonging to two sub-intervals of width N1 within the interval ∆i (which itself lies in [α d, 2α d]).
√
2 (2πR r)
Moreover, for such r (which is certainly at least α d), we can use Lemma 14 to upper bound Jd/2
d
1.3
by αd . Overall, we can upper bound the sum of integrals above by
N
X
i=1

C

3
d
2 1.3
√ ·
√ .
·
<
2
α d N αd
α d

Technical Results On Bessel functions

Lemma 12. For any ν ≥ 0 and x, |Jν (x)| ≤ 1. Moreover, for any ν ≥ 1 and x ≥ 3ν, Jν (x) is 1-Lipschitz
in x.
Proof. The bound on the magnitude follows from equation 10.14.1 in DLMF.
The derivative of Jν (x) w.r.t. x is given by −Jν+1 (x) + (ν/x)Jν (x) (see equation 10.6.1 in DLMF).
Since |Jν+1 (x)| and |Jν (x)|, for ν ≥ 1, are at most √12 (see equation 10.14.1 in DLMF), we have that the

magnitude of the derivative is at most √12 1 + xν ≤ √12 1 + 13 < 1.
To prove the lemmas below, we will need the following explicit approximation result for the Bessel
function Jd/2 (x), which is an immediate corollary of Theorem 5 in Krasikov [2014], plus some straight√
forward approximations (using the facts that for any z ∈ (0, 0.5], we have 1 − z 2 ≥ 1 − 0.3z and
0 ≤ z arcsin(z) ≤ 0.6z):
29

Lemma 13 (Krasikov [2014]). If d ≥ 2 and x ≥ d, then
s


2
(d + 1)π
+ fd,x x
Jd/2 (x) −
cos −
πcd,x x
4
where

r

d2 − 1
, fd,x = cd,x +
1−
4x2

cd,x =

√

≤ x−3/2 ,
√

d2 − 1
arcsin
2x

d2 − 1
2x

!
.

Moreover, assuming x ≥ d,
1 ≥ cd,x ≥ 1 −
and
1.3 ≥ 1 +
Lemma 14. If d ≥ 2 and r ≥

0.15 d
≥ 0.85
x

0.3 d
0.15 d
≥ fd,x ≥ 1 −
≥ 0.85
x
x

√
d, then

1.3
2
Jd/2
(2πRd r) ≤ √ .
r d
√
√
Proof. Using Lemma 13 (which is justified since r ≥ d and Rd ≥ 15 d by Lemma 1), the fact that cos is
at most 1, and the assumption d ≥ 2,
r
2
Jd/2 (2πRd r) ≤
+ (2πRd r)−3/2
π · 0.85 · 2πRd r
!
r
1
2
1
=√
+
0.85π 2πRd r
2πRd r
s
!
r
5
2
5
√
+ √ √
≤
0.85π 2π d d
2π dr
s
!
r
5
2
5
√
≤
+
0.85π 4π
2π dr
Overall, we have that
2
Jd/2
(2πRd r)

5
√
≤
2πr d

r

2
5
+
0.85π 4π

!2
≤

1.3
√ .
r d

Lemma 15. For any β ≥ 1, d ≥ 2 such that βd ≥ 127, it holds that
Z

2βd

βd

2 (x)
Jd/2

x


·1

2
Jd/2
(x)

30

1
≥
20x


dx ≥

0.005
.
βd

Proof. For any a, b ≥ 0, we have a · 1 {a ≥ b} ≥ a − b. Therefore,


Z 2βd
1
1
2
2
(x) ≥
(x) · 1 Jd/2
dx
· Jd/2
x
20x
βd


Z 2βd
1
1
2
≥
dx
· Jd/2 (x) −
x
20x
βd
Z 2βd
Z 2βd
1
1
1 2
Jd/2 (x)dx −
dx
=
x
20 βd x2
βd
Z 2βd
1 2
1
=
Jd/2 (x)dx −
.
x
40βd
βd
We now wish to use Lemma 13 and plug in the approximation for Jd/2 (x). To do so, let a = Jd/2 (x), let b
be its approximation from Lemma 13, and let  = x−3/2 the bound on the approximation from the lemma.
Therefore, we have |a − b| ≤ . This implies
a2 ≥ b2 − (2|b| + ),

(38)

which follows from
b2 − a2 = (b + a)(b − a) ≤ (|b| + |a|)|b − a| ≤ (|b| + |b| + ) = (2|b| + ).
Eq. (38) can be further simplified, since by definition of b and Lemma 13,
s
r
1
2
2
≤
≤√ .
|b| ≤
πcd,x x
π · 0.85 · x
x
Plugging this back into Eq. (38), plugging in the definition of a, b, and recalling that cd,x ≤ 1 and x ≥ d ≥ 2,
we get that

 

(d + 1)π
2
2
2
−3/2
2
cos −
+ fd,x x − √ + x
x−3/2
Jd/2 (x) ≥
πcd,x x
4
x


2
(d + 1)π
2
≥
cos −
+ fd,x x − 3x−2 .
πx
4
Therefore,
Z

2βd

1 2
J (x)dx
x d/2

βd

≥
=

2
π

Z

2
π

Z

2βd

βd
2βd
βd



Z 2βd
1
(d + 1)π
2
cos
−
+
f
x
dx
−
3
x−3 dx
d,x
x2
4
βd


1
(d + 1)π
9
+ fd,x x dx − 2 2 .
cos2 −
x2
4
8β d

To compute the integral above, we will perform a variable change, but first lower bound the integral in a
more convenient form. A straightforward calculation (manually or using a symbolic computation toolbox)
reveals that
r
∂
d2 − 1
(fd,x x) = 1 −
,
∂x
4x2
31

which according to Lemma 13, equals cd,x , which is at most 1. Using this and the fact that fd,x ≥ 0.85 by
the same lemma ,


Z 2βd
(d + 1)π
1
2
cos −
+ fd,x x dx
x2
4
βd



Z 2βd
1
(d + 1)π
∂
2
≥
cos
−
+
f
x
(f
x)
dx
d,x
d,x
x2
4
∂x
βd



Z 2βd
0.852
(d + 1)π
∂
2
≥
cos −
+ fd,x x
(fd,x x) dx
(fd,x x)2
4
∂x
βd

Using the variable change z = fd,x x, and the fact that 1.3 ≥ fd,x ≥ 0.85, the above equals
2

Z

fd,2βd 2βd

0.85

fd,βd βd





Z 1.7βd
1
(d + 1)π
1
(d + 1)π
2
2
2
cos −
+ z dz ≥ 0.85
cos −
+ z dz
2
z2
4
4
1.3βd z






(d+1)π
∂
1
z
We now perform integration by parts. Note that cos2 − (d+1)π
+
z
=
+
sin
−
+
2z
,
4
∂z 2
4
2
and sin is always bounded by 1, hence


Z 1.7βd
1
(d + 1)π
2
cos −
+ z dz
2
4
1.3βd z




(d+1)π
z
1
Z 1.7βd z + 1 sin − (d+1)π + 2z
+
sin
−
+
2z
1.7βd
2
4
2
2
4
2
=
dz
+2
2
3
z
z
1.3βd
1.3βd



(d+1)π

Z 1.7βd 
sin
−
+
2z
1.7βd
2
1
1
1


≥
+
+
−
dz
2z
4z 2
z 2 2z 3
1.3βd
1.3βd



(d+1)π


+
2z
sin
−
1.7βd
1.7βd
2
1
1
1


=
+
+ − + 2
2
2z
4z
z 4z
1.3βd
1.3βd




1 + sin − (d+1)π
+ 2z
1.7βd
2
1

= − +
2
2z
4z
1.3βd



(d+1)π


1
+
sin
−
+
2z
1.7βd
1.7βd
2
1

+
= −
2
2z
4z
1.3βd
1.3βd


0.09
1+1
≥
−
βd
4(1.3βd)2


1
1
=
0.09 −
.
βd
3.38βd

32

Concatenating all the lower bounds we attained so far, we showed that


Z 2βd
1
1
2
2
dx
· J (x) · 1 Jd/2 (x) ≥
x d/2
9x
βd


1
2
1
9
2 1
≥−
+ 0.85
0.09 −
−
40βd 8β 2 d2 π
βd
3.38βd


1.27
1
0.015 −
.
≥
βd
βd
If βd ≥ 127, this is at least

0.005
βd ,

from which the lemma follows.

33

