arxiv: v2 [math.pr] 15 Nov 2016

Transcription

1 STEIN S METHOD, MANY INTERACTING WORLDS AND QUANTUM MECHANICS Ian W. McKeague, Erol Peköz, and Yvik Swan Columbia University, Boston University, and Université de Liège arxiv: v2 [math.pr] 5 Nov 206 Abstract Hall, Deckert and Wiseman (204 recently proposed that quantum theory can be understood as the continuum limit of a deterministic theory in which there is a large, but finite, number of classical worlds. A resulting Gaussian limit theorem for particle positions in the ground state, agreeing with quantum theory, was conjectured in Hall, Deckert and Wiseman (204 and proven by McKeague and Levin (206 using Stein s method. In this article we propose new connections between Stein s method and Many Interacting Worlds (MIW theory. In particular, we show that quantum position probability densities for higher energy levels beyond the ground state may arise as distributional fied points in a new generalization of Stein s method. These are then used to obtain a rate of distributional convergence for conjectured particle positions in the first energy level above the ground state to the (two-sided Mawell distribution; new techniques must be developed for this setting where the usual density approach Stein solution (see Chatterjee and Shao (20 has a singularity. Keywords: Interacting particle system, Higher energy levels, Mawell distribution, Stein s method INTRODUCTION Hall, Deckert and Wiseman (204 proposed a many interacting worlds (MIW theory for interpreting quantum mechanics in terms of a large but finite number of classical worlds. In the case of the MIW harmonic oscillator, an energy minimization argument was used to derive a recursion giving the location of the oscillating particle as viewed in each of the worlds. Hall et al. conjectured that the empirical distribution of these locations converges to Gaussian as the total number of worlds N increases. McKeague and Levin (206 recently proved such a result and provided a rate of convergence. More specifically, McKeague and Levin showed that if,... N is a decreasing, zero-mean sequence of real numbers satisfying the recursion relation n+ = n + + n, (. then the empirical distribution of the n tends to standard Gaussian when N. Here n represents the location of the oscillating particle in the nth world, and the Gaussian limit distribution agrees with quantum theory for a particle in the lowest energy (ground state.

2 The hypothesized correspondence with quantum theory suggests that stable configurations should also eist at higher energies in the MIW theory. Moreover, the empirical distributions of these configurations should converge to distributions with densities of the form p k ( = (He k( 2 ϕ(, k! R, (.2 where ϕ( is the standard normal density, He k ( = ( k e 2 /2 dk d k e 2 /2 is the (probabilist s kth Hermite polynomial, and k is a non-negative integer. The ground state discussed above corresponds to k = 0 and has the standard Gaussian limit. However, the question of how to characterize higher energy MIW states corresponding to k is still unresolved as far as we know. The energy minimization approach of Hall, Deckert and Wiseman (204 starts with an analysis of the Hamiltonian for the MIW harmonic oscillator: H 0 (, p = E(p + V ( + U 0 (, where the locations of particles (having unit mass in the N worlds are specified by = (,..., N with > 2 >... > N, and their momenta by p = (p,..., p N. Here E(p = N n= p2 n/2 is the kinetic energy, V ( = N n= 2 n is the potential energy (for the parabolic trap, and U 0 ( = N ( n+ n n n n= is called the interworld potential, where 0 = and N+ =. In the ground state, there is no movement because all the momenta p n have to vanish for the total energy to be minimized. In this case, as mentioned above, Hall, Deckert and Wiseman (204 showed that the particle locations n satisfy (. and McKeague and Levin (206 showed that the empirical distribution tends to a standard Gaussian distribution. Our contribution in the present article is to derive an interworld potential for the second energy state (k = and show that the empirical distribution of the configuration that minimizes the corresponding Hamiltonian has a limit distribution that again agrees with quantum theory. The interworld potential in this case is shown to be U ( = 9 N ( n+ 3 n 3 n 3 n (.3 n n= and the minimizer of the corresponding Hamiltonian H (, p = E(p+V (+U ( is shown to satisfy the recursion ( n 3 n+ = 3 n 3. (.4 Further, we show that if,..., N is a decreasing, zero-mean solution, then the empirical distribution of the n converges to the (two-sided Mawell distribution 2 i= i 2

3 having density p ( = 2 e 2 /2 / 2π. The entire sequence,... N should be viewed as indeed by N, though we suppress notation for this dependence and write,... N instead of,n,..., N,N. We also give a rate of convergence using a new etension of Stein s method. Our approach is generalizable to recursions that converge to the distributions of other higher energy states of the quantum harmonic oscillator, although we do not pursue such etensions here. We initially thought that the MIW interpretation could be based on a universal interworld potential function U 0 that applies to all energy levels, with the densities p k ( then arising as limits of local minima of H 0. However, this idea turned out to be analytically unworkable. Here we propose an alternative approach in terms of adapting the interworld potential to each higher energy level. Minimizing the resulting Hamiltonian is then tractable and the solution can be shown to converge to p k (, at least in the case k =. Hall, Deckert and Wiseman (204 derived their interworld potential U 0 as a discretization of Bohm s quantum potential summed over the particle ensemble, see Bohm (952. The challenge in general is to etend this derivation to higher-energy wave functions in a way that leads to an eplicit recursion minimizing the resulting Hamiltonian, and to show that it agrees with p k ( in the limit. A major contribution here, in addition to providing a rate of convergence, is a general method for finding such interworld potential functions and their associated particle recursions. Stein s method (see Stein (986, Chen, Goldstein and Shao (200 and Ross (20 is a well established technique for obtaining eplicit error bounds for distributional limit theorems. The usual density approach (see Chatterjee and Shao (20 for applying Stein s method to arbitrary random variables does not seem to apply in cases where the density function vanishes at a point that is not at the endpoints of the range of the random variable (here we have p (0 = 0 and the random variable can take both negative and positive values; in this case the solution to the Stein equation will have a singularity and also unbounded derivatives, and this therefore requires the new technique we give here to handle such distributions. It is interesting to note that, while there are plenty of eamples of Stein s method applied to distributions with a density having a zero at the endpoints of the range of the random variable (the gamma and beta distributions, for eample, there have been no eamples that we know of where there are zeros inside the range; the higher energy distributions p k (, for k > 0, appear to be the first such distributions considered. The price one has to pay with our approach for handling these zeros is that more complicated estimates must be made from the couplings. In this case we have an eplicit representation of the recursion, and therefore the coupling, and this leads to the possibility of making these estimates. In Section 2 we generalize the argument of Hall, Deckert and Wiseman (204 to derive the interworld potential, and show how it leads to the solution (.4. In Section 3 we introduce the notion of a generalized zero-bias transformation, and show that the distributional properties of eigenstates of the quantum harmonic oscillator can be characterized in terms of fied points of this transformation. Also, we derive the generalized zero-bias distribution for the empirical distribution of general configurations. Section 4 develops our results based on the new etension of Stein s method to show convergence of the configuration that minimizes the Hamiltonian of the second energy state. 3

4 2 INTERWORLD POTENTIALS FOR HIGHER ENERGY STATES Hall, Deckert and Wiseman (204 introduced their MIW theory from the perspective of the de Broglie Bohm interpretation of quantum mechanics, which is mathematically equivalent to standard quantum theory. They used this approach to construct an ansatz for the conjectured interworld potential U 0 governing the ground state wave function of the quantum harmonic oscillator. In this section we introduce an etended version of this ansatz aimed at providing a MIW characterization of the higher energy eigenstates. Our argument follows along the lines of Section IIIA of Hall, Deckert and Wiseman (204 with the major difference being that we now need to introduce a more general way of approimating the density of particle location for a stationary wave function ψ(, namely for a density of the form p( = ψ( 2 = b(ϕ(, where b( is a non-negative, even, smooth function having finitely many zeros. Here b represents a baseline that varies more rapidly than ϕ(. Let > 2 >... > N. Bohm s quantum potential summed over the ensemble { n } is defined by U ψ ( = N [p ( n /p( n ] 2 (2. n= where we are using dimensionless units. An approimation to p( n based on ignoring ϕ( is given (up to a normalizing constant by where B( = 0 p ( n p( n p( n = b( n B( n B( n+, b(t dt is the cumulative baseline function. This suggests p( [ n p( n ( n n p( n B( n B( n+ B( n B( n ] b( n, where we set B( 0 = and B( N+ =. Our proposed ansatz for the interworld potential is then based on inserting the above epression into (2. to obtain U b ( = N [ ] 2 B( n+ B( n b( n 2. (2.2 B( n B( n n= Note that our earlier assumptions about b imply that B is strictly increasing, so U b is well-defined. In the simplest cases b( = and b( = 2 the above epression for U b agrees with the interworld potentials U 0 and U defined in the Introduction. Specializing to the case b( = 2, the following argument characterizes the minimizer of the Hamiltonian H (i.e., the ground state when the interworld potential is U in terms of a solution to the recursion (.4. In any ground state the particles do not move, so the kinetic energy E vanishes. Then, adapting the argument of Hall, Deckert and Wiseman (204 to apply to H, we have 9(N 2 = 9 [ N 3 n+ 3 n 3 n+ 3 n n= ] 2 4

5 [ N ( = 9 3 n= n+ 3 n [ N ( 9 3 n+ 3 n n= U ( V (, 3 n 3 n 3 n 3 n ] 2 2 n( n 3 N /2 n 2 4 n ] [ N ] ( n 3 N /2 n 2 n= where the first inequality is Cauchy Schwarz. So U 9(N 2 /V, leading to H = U + V 9(N 2 /V + V 6(N with the last inequality being equality for V = 3(N. It follows that H is minimized when V = 3(N, the mean 3 N of {3 n, n =,..., N} vanishes, and [ ] = α n 3 n+ 3 n 3 n 3 n for some constant α. The sum of the right of the above display telescopes, leading to the recursion (.4 by rearranging and noting that α = V/(N = 3. Density Figure : Eample with b( = 2, N = 22, showing the piecewise constant density having mass /(N uniformly distributed over the intervals between successive n compared with the Mawell density, where the breaks in the histogram are the successive n satisfying the recursion (.4. The following lemma provides the basic properties we need to ensure the eistence of a solution of the Mawell recursion (.4 that minimizes the Hamiltonian H, as well as ensuring that the solution is unique. This result is analogous to Lemma of McKeague and Levin (206 concerning solutions of (., but the difference here is that the variance is 3, agreeing with the Mawell distribution (rather than close to standard normal in the case of (.. 5

6 Lemma 2.. Suppose N is even. Every zero-median solution,..., N of (.4 satisfies: (P Zero-mean: N = 0. (P2 Mawell variance: N = 3(N. (P3 Symmetry: n = N+ n for n =,..., N. Further, there eists a unique solution,..., N such that (P and (P4 Strictly decreasing: >... > N hold. This solution has the zero-median property, and thus also satisfies (P2 and (P3. Proof. The proof follows identical steps to the proof of Lemma of McKeague and Levin (206, apart from the variance property (P2, which is proved using (P and (P3 as follows. Denote S n = n i= i for n =,..., N, and set S 0 = 0. Using (.4 we can write N N N 3(N = 3 S n Sn = S n ( 3 n 3 n+ = [(S n + n 3 n S n 3 n+] = n= n= N [S n 3 n S n 3 n+ + 2 n] n= = N S N 3 N, where we used the recursion in the second equality, and the last equality is from a telescoping sum. (P3 implies S N = 0, so S N = / N, and (P2 follows. Although in the sequel we concentrate on the case k = (see Figure, to conclude this section we briefly discuss general densities of the form p k given in (.2. The above argument for b( = 2 can be etended to general U b under the condition that B( is proportional to b(, which is the case when b( is proportional to r for some even non-negative integer r (but not for the square of the kth Hermite polynomial unless k = 0 or. Under this condition, it can be shown that the minimizer of the Hamiltonian based on U b is a symmetric solution of the recursion ( n i B( n+ = B( n. (2.3 b( i We have not been able to show that this recursion minimizes the Hamiltonian for general b, but our numerical results suggest that it is very close if not identical to a minimizer. With k = 2 we have b( = ( 2 2 /2, B( = 5 /0 3 /3 + /2, and the symmetric solution of the resulting recursion produces a remarkably good agreement with p k, see Figure 2. i= n= 6

7 Density Figure 2: Eample with b( = He k ( 2 /k! for k = 2, N = 4, where the breaks in the histogram are the successive n satisfying the recursion (2.3 and the red curve is p k (. 3 GENERALIZED ZERO-BIAS TRANSFORMATIONS Let W be a symmetric random variable and b: R R a non-negative function such that σ 2 = E[W 2 /b(w ] <. Goldstein and Reinert (997 gives a distributional fied point characterization of the Gaussian distribution, which we generalize in the definition below. Definition 3.. If there is a random variable W such that [ ] [ ] σ 2 E f (W = E W f(w b(w b(w for all absolutely continuous functions f : R R such that E W f(w /b(w <, we say that W has the b-generalized-zero-bias distribution of W. Remark 3.2. Goldstein and Reinert (997 study the case b( = and show that W has the same distribution as W if and only if W has a Gaussian distribution. Distributional fied point characterizations for eponential, gamma and other nonnegative distributions and the connection with Stein s method have been studied in Peköz and Röllin (20, Peköz, Röllin and Ross (203, and Peköz, Röllin and Ross (206. Remark 3.3. By a routine etension of the proof of Proposition 2. of Chen, Goldstein and Shao (200, it can be shown that there eists a unique distribution for W, and it is absolutely continuous with density [ ] p ( b(e W b(w W. 7

8 We note in passing that the σ 2 is misplaced in the first display of Chen et al. s proposition, which corresponds to b( =, the usual zero-bias distribution of W. The composition of the b-generalized-zero-bias transformation with the (/b-generalizedzero-bias transformation is the usual zero-bias transformation. Remark 3.4. With ϕ the standard normal density and b a ϕ-integrable function, if W has density p( = b(ϕ(, (3. then its distribution is a fied point for the b-generalized-zero-bias transformation since p t ( = b( p(t dt = b( tϕ(t dt = p(. b(t The following result gives the b-generalized-zero-bias distribution of the uniform distribution on N points. Proposition 3.5. Given an integer N >, let > 2 >... > N be such that b( n > 0 for all n. Let P N be the empirical distribution of the n : P N (A = #{n: n A} N for any Borel set A R. Under the symmetry condition n = N n+ for n =,..., N, the b-generalized-zero-bias distribution P N of P N is defined, and has density [ n ] p i ( b( b( i for n+ < n (n =,..., N, and p ( = 0 if > or N. Proof. Immediate from Remark 3.3. i= Recall the following distances between distribution functions F and G. Kolmogorov distance is The d K (F, G = sup F ( G(, R and the Wasserstein distance is d W (F, G = sup h H R h df R h dg where H = {h: R R Lipschitz with h } and is the supremum norm. Using Proposition.2 in Ross (20, these two metrics are seen to be related by d K (F, G 2C d W (F, G if G has density bounded by C. Restricting attention to the special case b( = 2, we can now state our main result, along with an important corollary. 8

9 Theorem 3.6. Suppose W is constructed on the same probability space as the zeromean random variable W and is distributed according to the 2 -generalized-zero-bias distribution of W. Let M have the two-sided Mawell density 2 e 2 /2 / 2π. Then there eist positive finite constants λ, λ 2, λ 3 and λ 4 such that d W (L (W, L (M λ E W W + λ 2 E [ W W W ] + λ 3 E W W + λ 4E W W. (3.2 Proof. The inequality follows immediately from Theorem 4.4. Finiteness of the constants (along with eplicit upper bounds is detailed in Proposition 4.5. The following corollary gives a rate of convergence of the solution to (.4 to the two-sided Mawell distribution in terms of the Wasserstein distance; we postpone the proof until Section 4.3. Corollary 3.7. Suppose,... N is a monotonic, zero-mean, finite sequence of real numbers satisfying (.4, let P N be the empirical distribution of these values, and let M be as in Theorem 3.6. Then there is a constant c > 0 such that log N d W (P N, L (M c N. 4 THE STEIN EQUATION AND ITS SOLUTIONS 4. General considerations Let X have a density p defined as in (3.. The first step is to identify an appropriate Stein equation and bound its solutions. Let h be such that E[ h(x] = 0. There are many possible starting points. The well known density approach (see Chatterjee and Shao (20 starts with the Stein equation which is easily solved to yield f( = f ( + p ( f( = h( p( h(ub(uϕ(u du. b(ϕ( If b vanishes at a point, as it does in the two-sided Mawell case where b(0 = 0, then this solution f will have a singularity and we cannot carry on with the usual program for applying Stein s method. For these types of distributions we propose a new approach; the price we pay here is the necessity to bound several additional quantities concerning the couplings we obtain. The eplicit nature of the recursion here allows us to compute these quantities. In view of Definition 3. it is natural to consider the Stein equation f (w b(w w f(w b(w = h(w (4. 9

10 which can be solved using the usual normal approimation solution, but the resulting estimates will again rest on properties of f/b which is unbounded in the cases we are interested in. Because of this, we introduce an original route which leads to the correct order bounds we are seeking. First, following Ley, Reinert and Swan (206 we introduce an integral operator associated with the Gaussian density: h T ϕ h(w := ϕ(w w h(uϕ(u du (4.2 (called the inverse Stein operator mapping functions with Gaussian-mean zero and sufficiently well-behaved tails into bounded functions. Net define the Stein kernel of X (or, equivalently, of p by τ X ( = p( up(u du. The Stein kernel satisfies the integration by parts formula E [τ X (Xf (X] = E [Xf(X] (4.3 for all sufficiently regular functions f. Stein kernels were introduced in Stein (986, Cacoullos and Papathanasiou (989, and have proven to be of great use in Gaussian analysis, see, e.g., Nourdin and Peccati (2009 and Chatterjee (2009. Remark 4.. As discussed in the Introduction, our concern in this paper is with symmetric densities of the form p k ( given in (.2. The Stein kernel of p k is τ k ( = and direct integration leads to uhe k (u 2 ϕ(u du/(ϕ(he k ( 2 τ ( = , τ 2 ( = ( 2 2, τ 3 ( = ( (4.4 The net result applies to general densities, but in the following sections we focus on further results for the special case τ (. We aim to assess the proimity between the law of X and some W by estimating the Wasserstein distance between their distributions. Suppose that there eists a W following the b-generalized-zero-bias distribution of W and defined on the same probability space as W. To each integrable test function h we associate the function f := f h, the solution to (4. with h = h Eh(X. This association is unique in the sense that there eists only one absolutely continuous version of f satisfying (4. at all points, see e.g. Chen, Goldstein and Shao (200. We then write (supposing here and in the sequel that σ 2 = [ E [h(w ] E [h(x] = E f (W b(w W f(w ] b(w 0

11 [ = E f (W b(w f ] (W. (4.5 b(w One of the major differences between the present setting and the usual applications of Stein s method is that here we cannot bound the right hand side of (4.5 directly because solutions f of (4. have a singularity at 0. In order to bypass this difficulty it was therefore necessary to introduce some further tweaking of the method which we now detail. Such results should also have an intrinsic interest for users of Stein s method as it bears natural generalizations for density functions of the form p( = b(p 0 ( with p 0 and b well chosen. Proposition 4.2. Let b( be a nonnegative even function with support in (, such that lim ± b(ϕ( = 0. Suppose furthermore that b is absolutely continuous and integrable w.r.t. ϕ with integral b(ϕ( d =. Let X be a random variable with density b(ϕ(. Then τ X ( = + T ϕ b ( b( (4.6 under the convention that the ratio is set to zero at all points such that b( = 0 and T ϕ b ( 0. Further, with h defined by (4., g h ( = is the unique bounded solution of the ODE b(u h(uϕ(u du b(ϕ( + b (uϕ(u du (4.7 τ X (g ( g( = h. (4.8 Proof. Integrating by parts in the definition of the Stein kernel for p = bϕ we get (assuming that lim ± b(ϕ( = 0 + yp(y dy = = b(ϕ( + b(y( ϕ (y dy b (yϕ(ydy so that (4.6 follows by definition (4.2 of the inverse Stein operator. For the second claim note how τ X (g ( g( = (τ X(g(p( p( so that the conclusion (along with unicity follows from the same argument as in Nourdin and Peccati (202, Proposition Intuition (supported e.g. by Stein (986, Lesson VI or the more recent work Döbler (205 encourages us to claim that functions (4.7 will have satisfactory behavior. It is thus natural to seek a connection between equations of the form (4. and (4.8. This we summarize in the net lemma.

12 Lemma 4.3. Let all notations be as above and introduce the function g = g f defined at all through f ( f( = τ X (g ( g(. b( Then E Proof. Since and [ f (W b(w f ] (W = E [W (τ b(w X (W g(w W (τ X (W g(w ] + E [τ X (W g (W τ X (W g (W ]. (4.9 f ( f( b( = (f(ϕ( b(ϕ( τ X (g ( g( = (b(τ X(g(ϕ( b(ϕ( at all for which b( 0, we deduce that f and g are mutually defined by f = (bτ X g. This in turn gives ( f ( b b( = ( b( τ X( + τ X( g( + τ X (g ( =: ψ(g( + τ X (g ( which, combined with ψ( = (τ X ( (that is easily derived using the various definitions involved, leads to the useful identity from which (4.9 is directly derived. f ( b( = (τ X( g( + τ X (g ( (4.0 In view of Lemma 4.3 and identity (4.5 it remains to find bounds on the two terms on the right hand side of ( Approimating the two-sided Mawell distribution Theorem 4.4. Let p( = 2 ϕ(, and take f a solution to the Stein equation f (w/w 2 wf(w/w 2 = h(w, (4. where h is a function having bounded first derivative and zero-mean under p. Set c = h. Then for any coupling of W and W on a joint probability space such that W has the 2 -generalized zero biased distribution for W, [ f E (W (W f ] (W λ 2 (W 2 E W W + λ 2 E [ W W W ] + λ 3 E W W + λ 4E W W (4.2 with λ 6c, λ 2 7c, λ 3 8c and λ 4 22c. (4.3 2

13 Proof. With b( = 2 we have τ X ( = + 2/ 2 and ψ( = 2/, so that (4.9 becomes [ 2 = E W g(w 2 ] [( W g(w + E + 2 ( g (W + 2 ] g (W (W 2 (W [( 2 = 2E W ] [ ] g(w + 2E W W (g(w g(w [ + E [g (W g (W ] + 2E (W 2 g (W ] (W 2 g (W. The first two terms are dealt with easily to get [( 2 E W ] [ g(w + 2E W W (g(w g(w ] [ 2 g E W ] [ ] W + 2 g E W W W. For the last two terms we introduce the function to get on the one hand so that E [g (W g (W ] = E χ( = g (/ [ ] W g (W W g (W W W = E [(W W χ(w ] + E [W (χ(w χ(w ] E [g (W g (W ] χ E [ W W ] + χ E [ W (W W ] and, on the other hand [ E (W 2 g (W ] (W 2 g (W [ = E = E W χ(w ] W χ(w [( W ] χ(w + E W so that [ 2 E (W 2 g (W ] [ (W 2 g (W 2 χ E W ] W [ ] W (χ(w χ(w [ ] + 2 χ E W W W. Combining these different estimates we obtain (4.2, with λ, λ 2, λ 3 and λ 4 epressed in terms of χ, χ, g and g as follows: λ = χ, λ 2 = χ, λ 3 = 2( g + χ and λ 4 = 2( g + χ. The inequalities in (4.3 are proved in the Proposition 4.5 below. 3

14 The net step is to bound χ, χ, g and g in a non trivial way; this we achieve in the net proposition. Proposition 4.5. Let h : R R be absolutely continuous and integrable with respect to p( = 2 ϕ(. Set c = h which we suppose to be finite. Let X p, define { e 2 /2 g 0 ( = y2 (h(y E [h(x] e y2 /2 dy if > 0 e 2 /2 y2 (h(y E [h(x] e y2 /2 dy if 0 (4.4 and set Then g( = g 0( and χ( = g (. (4.5 g 3c, g 4c, χ 6c and χ 7c. Remark 4.6. The function g 0 defined in (4.4 satisfies g 0( g 0 ( = ( 2 h( E[Z 2 h(z] (4.6 with Z ϕ a standard Gaussian random variable. Remark 4.7. The function g defined in (4.5 satisfies (g(τ(p( p( = h( E [h(x] with X p and τ( = + 2/ 2. Proof. In order to simplify future notations we introduce Φ( = ϕ(t dt, Φ( = ϕ(t dt Υ( = /2 e2 t2 e t2 /2 dt and Ῡ( = /2 e2 t2 e t2 /2 dt. Using the identity b a t 2 e t2 /2 dt = ae a2 /2 be b2 /2 + we deduce that Υ( = + e 2 /2 and thus b a e t2 /2 dt, a < b, (4.7 e t2 /2 dt and Ῡ( = + e2 /2 e t2 /2 dt Υ(, Ῡ( + π 2 Υ( at all R and lim Ῡ( = lim =. (4.8 The proof is now broken down into several steps. Step : rewrite the solutions. Following Chen, Goldstein and Shao (200, page 39 we rewrite the test functions in term of their derivatives (still with Z a standard normal random variable h(y E [h(x] = h(y E [ Z 2 h(z ] 4

15 y = = y z 2 (h(y h(z ϕ(z dz ( y z 2 h (t dt ϕ(z dz z y z 2 ( z y h (t dt ϕ(z dz. Changing the order of integration then using (4.7 leads to the rhs becoming y [ t ] [ ] h (t z 2 ϕ(z dz du h (t z 2 ϕ(z dz dt y t = = and thus h (t [ tϕ(t + h (ttϕ(t dt + h(y E [h(x] = y t y ] ϕ(zdz dt h (tφ(t dt h (tφ(tdt We deduce the following useful bound h( E [h(x] y y y [ h (t tϕ(t + h (t Φ(t dt, t ] ϕ(z dz dt h (t Φ(tdt E [Zh (Z]. (4.9 c 2( + / 2π + 2/π c. (4.20 Plugging (4.9 in (4.4 leads to (we restrict the discussion to > 0, the other case following by symmetry g 0 ( = E[Zh (Z]Υ( + e 2 /2 e 2 /2 y y y 2 e y2 /2 h (tφ(t dtdy y 2 e y2 /2 h (t Φ(t dtdy =: I( =: II( =: III( To deal with the quantities II( and III( we again interchange integrations to get ( II( = e 2 /2 y 2 e y2 /2 dy h (tφ(t dt ( + e 2 /2 y 2 e y2 /2 dy h (tφ(t dt and = Υ( III( = e 2 /2 = e 2 /2 t h (tφ(t dt + e 2 /2 ( t y 2 e y2 /2 dy h (t Φ(t dt ( e 2 /2 Υ( e t2 /2 Υ(t 5 e t2 /2 Υ(th (tφ(t dt h (t Φ(t dt

16 = Υ( and thus if 0 we have h (t Φ(t dt e 2 /2 g 0 ( = E[Zh (Z]Υ( + Υ( Υ( h (t Φ(t dt + e 2 /2 h (tφ(t dt By a similar argument we deduce that if < 0 then g 0 ( = E[Zh (Z]Υ( + Υ( Υ( h (tφ(t dt + e 2 /2 h (t Φ(t dt e t2 /2 Υ(th (t Φ(t dt e t2 /2 Υ(th (t dt. (4.2 e t2 /2 Υ(th (t dt. (4.22 Step 2: a bound on g. Supposing h c we can use (4.2 and the first claim in (4.8 to deduce that for 0: g 0 ( ce Z ( + π/2 + c( + π/2 + c( + π/2 Φ(t dt + ce 2 /2 Φ(t dt e t2 /2 (t + π/2 dt. The last two terms decrease strictly to 0 as, with maimum value c/2 and c( + π/2, respectively. The first term is equal to c( 2/π + and the second one is equal to c( + π/2 Φ(t dt = c( + π/2 (Φ( + ϕ( c( 2 + ( π/2 + / 2π + /2. Similar (symmetric bounds hold for 0 and thus, collecting all these estimates, we may conclude: g( = g 0( 3c. ( Step 3: a bound on g. Here we start by rewriting the derivative as g ( = g 0( ( g 0(. ( Using (4.23, the second summand is easily seen to be uniformly bounded (by 3c. We are left with the first summand for which we start by rewriting the numerator, for 0, using (4.2: g 0( = Υ (E[Zh (Z] + Υ ( 6 h (tφ(t dt

17 which leads to Υ ( h (t Φ(t dt + Υ( ( h (Φ( + h ( Φ( + e 2 /2 Υ(th (te t2 /2 dt e 2 /2 Υ(h (e 2 /2 g 0( = Υ (E[Zh (Z] + Υ ( + e 2 /2 h (tφ(t dt Υ ( h (t Φ(t dt Υ(th (te t2 /2 dt. (4.25 Now we can use the fact that Υ ( = e 2 /2 /2 e t2 dt for all 0 as well as ( all the arguments outlined at the previous step to deduce the bound: g 0( 2 c + 2( + / 2π + π 2π 2c( + whence g 0( c c. (4.26 Similar (symmetric arguments hold also for negative and thus g ( 4c. Step 4: a bound on χ( = g (/. Using (4.24 we know that χ( = g 0( ( ( g 0(. ( The second summand in (4.27 is bounded using (4.22 to get For the first summand we use (4.25 to deduce ( g 0( 3c. (4.28 g 0( = Υ ( + Υ ( + e 2 /2 E[Zh (Z] h (tφ(t dt Υ ( Υ(th (te t2 /2 dt. h (t Φ(t dt At this stage it is useful to remark that, for 0, the function Υ (/ is strictly decreasing with maimal value ( ( π/2 and hence g 0 ( c π + 2π c (2 + 3 and thus g 0 ( ( c which, combined with (4.28, leads (after applying the symmetric arguments for 0 to χ( 6c. 7

18 Step 5: a bound on χ. Direct computations using (4.6 χ( = ( g ( ( h( E[Z 2 h(z] and thus χ ( = 23 g 0 ( ( ( 2 h( E[Z 2 h(z] g 0 ( h (. Using the bounds 2 3 /( , 2/( 2 + 2, (2 2 /( and /( as well as (4.23, (4.26 and (4.20 we conclude (after applying the symmetric arguments for 0 χ ( 7c. 4.3 Verifying bounds on epectations In this section we find bounds on the epectations in Theorem 3.6 in order to prove Corollary 3.7. We will make use of the following lemma. Lemma 4.8. If,..., N is the unique strictly decreasing zero-mean solution of (.4, then = O( log N. Proof. To simplify the notation, note that it suffices to consider the rescaled recursion 3 n+ = 3 n Sn, where S n is defined in the proof of Lemma 2.. By epressing 3 as a telescoping sum, m 3 = ( 3 n 3 n++ 3 m = n= m n= m Sn + 3 m n= (n/ + 3 m (+log m+ 3 m, where we have used Euler s approimation to the harmonic sum for the last inequality. By the variance property (P2 (in this rescaled case N = N we have that is bounded away from zero (as a sequence indeed by N and m is bounded, so 2 m/ is bounded. Dividing the above display by, we then obtain = O( log N. Proof of Corollary 3.7. From Proposition 3.5 and the recursion (.4, note that p ( puts mass /(N on each interval between successive n, so it is easy to create a coupling of W P N with W p ( such that W W n n+ when W [ n+, n ]. For a detailed proof of such a coupling, see the construction given in McKeague and Levin (206. From Lemma 4.8 we then have E W W N ( n n+ = 2 N N = O n= 8 ( log N N. (4.29

19 Second, using W = O( log N it follows immediately that ( log N E[ W W W ] = O. N Third, the zero-median property gives 2 3 m = 3 m 3 m+ = S m where m = N/2 +, so m / N. By symmetry E W W = E W W W ( m+, m] + 2 (m/ m = m /m, m n= E W W W ( n+, n]. From Proposition 3.5 note that p ( 2 for ( m+, m ]. Also using the fact that p ( puts mass /(N on this interval, the first term above can be written 6 m ( 3 m(n 0 ( 2 3 d m m (N = O. N The second term is bounded above by the telescoping sum 2 N m n= ( n+ n = 2 ( ( = O, N m N so we have Fourth, E E W W ( = O. N W W ( log N N E W W = O N using W m / N and (4.29. The Corollary now follows from Theorem 3.6. ACKNOWLEDGEMENTS The research of Ian McKeague was partially supported by NSF Grant DMS and NIH Grant 2R0GM The research of Yvik Swan was partially by the Fonds de la Recherche Scientifique - FNRS under Grant no F as well as IAP Research Network P7/06 of the Belgian State (Belgian Science Policy. We also thank the Institute for Mathematical Sciences at National University of Singapore for support during the Workshop on New Directions in Stein s Method (May 8 29, 205 where work on the paper was initiated. 9

20 REFERENCES D. Bohm (952. A suggested interpretation of the quantum theory in terms of hidden variables. I. Phys. Rev. 85, T. Cacoullos and V. Papathanasiou (989. Characterizations of distributions by variance bounds. Statist. Probab. Lett. 7, S. Chatterjee (2009. Fluctuations of eigenvalues and second order Poincaré inequalities. Probab. Theory Related Fields S. Chatterjee and Q.-M. Shao (20. Non-normal approimation by Stein s method of echangeable pairs with application to the Curie Weiss model. Ann. App. Probab. 2, L. Chen, L. Goldstein and Q.-M. Shao (200. Normal Approimation by Stein s Method. Springer Verlag. C. Döbler (205. Stein s method of echangeable pairs for the beta distribution and generalizations. Electron. J. Probab. 20, 34. L. Goldstein and G. Reinert (997. Stein s method and the zero bias transformation with application to simple random sampling. Ann. Appl. Probab. 7, M. J. W. Hall, D. A. Deckert and H. M. Wiseman (204. Quantum phenomena modeled by interactions between many classical worlds. Phys. Rev. X 4, C. Ley, G. Reinert and Y. Swan (206. Stein s method for comparison of univariate distributions. I. W. McKeague and B. Levin (206. Convergence of empirical distributions in an interpretation of quantum mechanics. Ann. Appl. Probab I. Nourdin and G. Peccati (2009. Stein s method on Wiener chaos. Probab. Theory Related Fields I. Nourdin and G. Peccati (202. Normal approimations with Malliavin calculus: from Stein s method to universality. Vol. 92. Cambridge University Press. E. Peköz and A. Röllin (20. New rates for eponential approimation and the theorems of Rényi and Yaglom. Ann. Probab. 39, E. Peköz, A. Röllin and N. Ross (203. Degree asymptotics with rates for preferential attachment random graphs. Ann. Appl. Probab. 23, E. Peköz, A. Röllin and N. Ross (206. Generalized gamma approimation with rates for urns, walks and trees. Ann. Probab., Vol. 44, No. 3, pp N. Ross (20. Fundamentals of Stein s method. Probab. Surv. 8, C. Stein (986. Approimate Computation of Epectations. Institute of Mathematical Statistics Lecture Notes Monograph Series, 7. Institute of Mathematical Statistics, Hayward, CA. 20