Nonlocality in Bell’s Theorem, in Bohm’s Theory, and in Many Interacting Worlds Theorising

“Locality” is a fraught word, even within the restricted context of Bell’s theorem. As one of us has argued elsewhere, that is partly because Bell himself used the word with different meanings at different stages in his career. The original, weaker, meaning for locality was in his 1964 theorem: that the choice of setting by one party could never affect the outcome of a measurement performed by a distant second party. The epitome of a quantum theory violating this weak notion of locality (and hence exhibiting a strong form of nonlocality) is Bohmian mechanics. Recently, a new approach to quantum mechanics, inspired by Bohmian mechanics, has been proposed: Many Interacting Worlds. While it is conceptually clear how the interaction between worlds can enable this strong nonlocality, technical problems in the theory have thus far prevented a proof by simulation. Here we report significant progress in tackling one of the most basic difficulties that needs to be overcome: correctly modelling wavefunctions with nodes.


Introduction
This paper is based loosely on the talk given by its final author at the 2017 Symposium on Emergent Quantum Mechanics (EmQM17). That talk took as its principal inspiration the following questions (from the list provided by the organisers to give focus to the meeting): • Is the universe local or nonlocal?
• What is the future of scientific explanation? Is scientific metaphysics, e.g., the notions of reality, causality, or physical influence, obsolete in mathematical accounts of the quantum world?
However, in view of the conference also being the David Bohm's Centennial Symposium, the speaker also briefly discussed some work relating to this question: • What is David Bohm's legacy for the future of quantum physics? This paper also concerns these two broad topics, but in the opposite ratio to that in the delivered talk. Section 2 of this paper addresses just the question "Is the universe local or nonlocal?" by examining the historical meaning of the term "local". (The second question listed above, which was addressed in the talk, and in References [1][2][3], is not addressed here). That section expresses the views of the invited contributor, not necessarily those of all the authors. Sections 3-6 concern Many Interacting Worlds (MIW), an approach to quantum mechanics inspired by Bohm's original hidden variables theory [4,5], and which thus we hope will be part of the answer to the last question above.
Having introduced the MIW approach in Section 3, we turn in Section 4 to how it can give rise to nonlocality; that is, how it can offer a positive answer to this (here abridged) focus question: • Are nonlocal connections-e.g., "action-at-a-distance"-fundamental elements in a radically new conception of reality?
However, as Section 4 notes, there are still major difficulties faced in demonstrating these connections in MIW simulations. Section 5 introduces a new MIW technique-a higher-order interaction potential-to address one of these difficulties, namely, wavefunction nodes. Section 6 applies this idea to progress the MIW simulation of stable excited state nodes of a quantum particle in one dimension. An alternative method to address the node problem is presented and simulated in Appendix A. Section 7 concludes with a discussion of open challenges for MIW theorising.

Is the Universe Local or Nonlocal?
The answer to this EmQM17 focus question hinges, of course, on what one means by the term "local" (assuming that "nonlocal" is simply its complement). Perhaps surprisingly, it seems [2,6] that the word was not used in the context of interpreting EPR quantum correlations prior to the 1964 paper of Bell [7]. In that paper, Bell proved his 1964 Bell's theorem; to quote [7]: In a theory in which parameters are added to quantum mechanics to determine the results of individual measurements, without changing the statistical predictions, there must be a mechanism whereby the setting of one measuring device can influence the reading of another instrument, however remote.
In other words, some quantum phenomena are incompatible with the joint assumption of predetermination (or causality; Bell used both terms) and locality (or separability; Bell used both terms). In the above quote, it is the negation of locality that is characterised, in a way consistent with Bell's earlier definition of locality; to quote [7]: It is the requirement of locality, or more precisely that the result of a measurement on one system be unaffected by operations on a distant system with which it has interacted in the past, that creates the essential difficulty.
Thus, Bell intended to be (somewhat) precise about what he meant by locality. Unfortunately he did not give a general mathematical definition, nor did he define terms like "unaffected" or (in the first quote and elsewhere) "influence".
In the theorem he proves that the role "locality" plays is the following. The assumption of predetermination means that an arbitrarily long time before any measurements are performed, there existed in the world a collection of hidden variables, λ, that, together with the future measurement settings, determines the future outcomes. Adding the "vital assumption [of locality]" implies that the outcome A of one party (Alice, say) cannot depend on the settings b of a distant party (Bob, say), but only by her own setting a. That is, in a theory θ with local predetermination of outcomes, there exists a function, A θ (a, λ) such that A = A θ (a, λ), and likewise for Bob. From this, Bell was able to derive his famous theorem, that there exist sets of measurement on entangled quantum systems whose results cannot be explained by any such model. (Note that Bell made an implicit assumption in 1964, related to freedom of choice, that P(λ|a, b) = P(λ). Here we follow Bell, as such an assumption is certainly necessary [8].) In the above formulation, "locality" has a precise meaning only in the context that predetermination has already been assumed. This is fine for the purpose to which Bell puts it. If one wanted to broaden the definition so that it applied independently of the assumption of predetermination then, it has been argued [2,6,9,10], the natural reading of Bell's verbal definition would be as follows. An arbitrary theory contains initial variables λ, which may or may not be hidden, and which may or may not be sufficient to determine all outcomes, and is described (in the limited context we are considering) by theoretical probabilities P θ (A, B|a, b, λ). The theory is local if there exists a function P θ (A|a, λ) such that and likewise for Bob. Others have argued [11][12][13][14], to the contrary, that this definition does not fit with the other use Bell makes of the concept of locality in his 1964 paper, which is the first paragraph of his Section II. There Bell implies that, according to the EPR argument [15], locality plus perfect correlations of outcomes implies predetermination. This shows, if nothing else, that Bell did think, in 1964, that locality was an assumption that had meaning prior to the assumption of predetermination. However, the obvious meaning, Equation (1), does not work in the EPR argument. Even one of Bell's most ardent admirers [13] was forced to admit this [14]: It is simply not clear how to translate Bell's words here (about locality) into a sharp mathematical statement in terms of which the EPR argument might be rigorously rehearsed. . . . [I]t must be admitted that Bell's recapitulation of the EPR argument in this paragraph leaves something to be desired.
Regardless, is there really a problem here? Bell does not say that he believes that locality plus perfect correlations implies predetermination. He merely says, at the beginning of that paragraph, ". . . the EPR argument is as follows" (my emphasis). In the preceding (opening) paragraph he is even weaker: "The paradox of Einstein, Podolsky and Rosen was advanced as an argument . . . " (my emphasis again). In a follow up paper in 1971 [16], Bell is weaker still. He gives three motivations for the assumption of predetermination (in a section entitled, unambiguously, "Motivations"). An EPR-style argument is the third motivation; he does not claim any logical deduction from locality but merely appeals to the intuition of the reader for the reasonableness of predetermination by hidden variables (see Reference [6] for details). The 1969 paper of CHSH [17] seems actively skeptical of the EPR argument, saying it "led them to infer that quantum mechanics is not a complete theory" (my emphasis). It is thus clear that, at this time, the EPR argument was far from being regarded as a rigorous proof for the necessity of hidden variables. (For a discussion of the extent to which the EPR argument is such a proof, see Reference [18].) The situation with regard to the EPR argument changed rapidly, at least in Bell's mind, after he formulated, in 1976 [19], a concept that does allow one to infer predetermination from perfect correlations. This was the concept of "local causality", stated most succinctly in a later paper [20] A consequence . . . of "local causality" [is] the outcomes [in the two labs] having no dependence on one another nor on the settings of the remote [measurement], but only on the local [measurement settings] and on the past causes.
In the situation considered above, a theory is "locally causal" only if there exists a function and likewise for Bob. Note the appearance of B as a conditional variable on the left hand side, which distinguishes this concept from "locality" in Equation (1 Bell clearly (and, I think [2,3], rightly) thought local causality to be a more natural concept than locality as per Equation (1), as he never used the latter concept again. Regrettably, however, he did not abandon the word "locality". Rather, beginning even in 1976 [19], he sometimes used "local" as short-hand for "locally causal", and, a few years later, was apparently convinced that local causality was the concept that he (and EPR) had always used [21] (for details, see [2]). However, at least in his final word on the subject [20], Bell showed his preference unequivocally for the terminology "local causality" over "locality".
Thus we may return to the primary question posed above-Is the universe local or nonlocal? If by "locality" one means "locally causal", the concept Bell promoted for most of his career in quantum foundations [22], then the answer (barring more exotic possibilities such as "superdeterminism" [20], retrocausality [23], and the subjectivity of macroreality [24]) is that the universe is nonlocal; it violates local causality. To avoid confusion, we might agree to say that the universe is Bell-nonlocal [2]. If, on the other hand, one adopts the definition of "local" indicated by Bell's 1964 paper and commonly used in text books [25,26], then the answer is that we cannot say whether the universe is local or nonlocal. Operational quantum mechanics satisfies this weaker sense of locality, simply because it does not feature signalling faster than light, and denies the need for giving any account for quantum correlations beyond an operational one. We can only say that the universe is nonlocal, in this strict sense, if we make some other assumptions about its nature, such as determinism.

David Bohm's Legacy: Permission to Theorise Radically New Conceptions of Reality
The most famous quantum theory which does make the assumption of determinism is of course David Bohm's [4,5]. Indeed, Bohm's theory was an inspiration to Bell who summarised the result of his 1964 paper in the introduction-cum-abstract as [A] hidden variable interpretation of elementary quantum theory [4,5] has been explicitly constructed. That particular interpretation has . . . a grossly nonlocal structure. This is characteristic, according to the result to be proved here, of any such theory which reproduces exactly the quantum mechanical predictions.
Why Bell considered Bohm's interpretation to be "grossly nonlocal", rather than nonlocal simpliciter, is unclear. Perhaps it was because the theory is nonlocal even in situations where there is an obvious local hidden variable theory, as in the EPR-correlations [15], or the EPR-Bohm correlations [27].
Unlike operational quantum mechanics, Bohm's theory is a precise and universal physical theory. Restricting to the case of interacting nonrelativistic scalar particles for simplicity of discussion, it takes the universe to be described by a universal wavefunction Ψ(q), obeying Schrödinger's equation, where q is the vectorised list of the coordinates of all the particles. However, it also postulates a single point in configuration space, x, which encodes the real positions of all these particles. This "marvellous point" [28] or "world-particle" [29] has a deterministic equation of motionẋ = v Ψ (x) guided vicinally by Ψ(q). (Note that "vicinal" is a synonym of "local" in the latter's quotidian sense, introduced here to avoid any possible confusion with "local" in the technical sense defined in Section 2.) Here, "vicinal guiding" means that the world-particle's velocity v Ψ (x) depends on Ψ(q) and finitely many derivatives, evaluated at q = x. However, vicinal in configuration space is not vicinal in 3D space-the positions of Bohmian particles in one region of 3D space can affect the motion of an arbitrarily distant particle if entanglement is present. Since it is the position of Bohmian particles that encodes what an experimenter decides to measure, this gives rise to the 'gross' nonlocality of Bohmian mechanics which Bell noted in 1964.
Bohm's original proposal [4,5] actually used a second-order dynamical equationẍ = a vicinal (x) + a Ψ (x). The 3D-vicinal acceleration a vicinal (x) is given by Newton's laws, involving inter-particle potentials which drop off with 3D-distance. The nonlocal effects in Bohm's theory arise from a separate, 3D-nonvicinal, quantum acceleration a Ψ (x). Bohm's publication of an explicitly nonlocal theory seems to have made it acceptable for other physicists to publish realist approaches to quantum mechanics in direct opposition to the Copenhagen interpretation. This included: de Broglie in 1956 [30], reviving his unpublished idea from 30 years earlier which prefigured much of Bohm's work and used the first-order dynamics described earlier; and Everett in 1957 [31], introducing the relative state interpretation, more popularly known as the many worlds interpretation.
Taking inspiration from both Bohm and Everett, two of us plus Deckert introduced in 2014 [29] what we called the Many Interacting Worlds (MIW) approach to quantum mechanics. We suggested that it might be possible to reproduce quantum phenomena without a universal wavefunction Ψ(q) (except to define initial conditions). In its place we postulated an enormous, but countable, ensemble X = x j : j of points x j in configuration space (similar ideas have been proposed earlier by a number of authors [32][33][34][35], but they considered a continuum of worlds, which, in our view, leads to some of the same conceptual issues that Everett's interpretation faces-see also [36]). Each point is a world-particle, just as Bohmian mechanics postulates, and the dynamics is intended to reproduce a deterministic Bohmian trajectory for each world-particle. However, the trajectory of a world-particle is guided not by a wavefunction, but by the locations of nearby world-particles in configuration space: where N j ⊂ X is the set of world-particles in the vicinity (somehow defined) of world-particle x j , and a MIW is some fixed (j-independent) function describing the interaction of many worlds.
In the MIW approach, probabilities arise only because observers are ignorant of which world x j they actually occupy, and so assign an equal weighting to all worlds compatible with the macroscopic state of affairs they perceive, in accordance with Laplace's principle. In a typical experiment, where the outcome is indeterminate in operational quantum mechanics, the final configurations of the worlds in the MIW approach can be grouped into different subsets, X o ⊂ X of world-particles, still extremely large in number, based on shared macroscopic properties corresponding to the different possible outcomes o. In Everett's approaches, these groups correspond to branches of the universal wavefunction, and one has to argue that the square modulus of the coefficient of each branch somehow manifests as the correct probability for the experimenter. In the MIW approach, the operational quantum probabilities will be equal to the number of worlds in each group divided by the number of worlds at the start.

General Considerations
Just as Bohmian dynamics is vicinal in configuration space but gives rise to nonlocality in 3D space, so too is the type of MIW dynamics just described. Indeed, it must be nonlocal because, like Bohmian mechanics, it is deterministic. However, we can be more specific about how this nonlocality, in Bell's 1964 sense, arises in the MIW approach.
Consider a branch Ψ E (in the Everettian sense) of the universal wavefunction Ψ. In the MIW approach, such a branch will correspond to a subset X E ⊂ X of world-particles, still extremely large in number, that are close together (on a macroscopic scale) in configuration space and mostly far (on a macroscopic scale) from the rest of X. Thus, to a good approximation, this set of world-particles will evolve autonomously, on the time scale of interest. Now consider a particular world-particle, x j ∈ X E , and a particular part of the vector x j , comprising a lower-dimensional vector b j , that contains the variables that encode the macroscopic fact of the decision by an experimenter, Bob, about what experiment to perform. That is, the decision is the same in every world under consideration here. However, there is no reason in the theory for b j to be correlated with the other variables in x j , until the experiment is actually performed. Once it is performed, the value of b j will have a direct (second-order in time) effect on the vicinal (in the 3D sense) variables in x j , say B j . Then, in the MIW approach, the change in B j will cause a change (again second-order in time) in x k ∈ N j (note that N j is only a tiny subset of X E ). Now say that in Ψ E , there is entanglement between the system corresponding to B j , and a distant system corresponding to some other observables A j . Then the change in x k will be a change not just in B k but also in A k . This is because the MIW acceleration function a MIW pays no heed to vicinality in the 3D sense-it cares only about the fact that the worlds x k ∈ N j are in the vicinity of the world x j in configuration space. Finally, the same interaction will cause a (second-order in time) change in A j , the distant observables in the original world j, since for at least some of the x k ∈ N j we will have x j ∈ N k .
Thus it is clear in principle how the MIW approach can give rise to nonlocality in Bell's 1964 sense, with Bob's choice b j causing a change in the variables A j , variables which correspond to the outcome observed by Alice in her distant, perhaps even space-like separated, lab, all in the same world j. It might seem that Bob could use this to signal to Alice faster than light. However, that is not so, for essentially the same reason that signal-locality is respected in Bohmian mechanics. Alice and Bob are ignorant of which world j they inhabit. They know only the macroscopic configuration X E in which their world must lie. A successful MIW theory would thus exhibit nonlocality at the hidden level of an individual world, but not at the observable level, averaging over all the worlds in X E .

Simulations
Can we explicitly simulate this in our MIW approach, showing, for example, the violation of a Bell inequality, or even just the EPR paradox, without signalling? Here we run up against the limitations in the development of our approach.
In our 2014 paper [29], we introduced a toy MIW theory to describe a very simple universe, comprising a single particle in 1D. We showed analytically that, in the limit of a large number of worlds, it gave the correct ground state distribution for a harmonic oscillator, and the correct description for the first and second moments of free particle evolution. We showed numerically that it could reproduce, at least qualitatively, the double-slit interference phenomenon. Furthermore, we argued that it can plausibly reproduce other generic quantum phenomena such as barrier tunnelling and reflection.
To demonstrate nonlocal correlations between measurements on distant quantum systems, we obviously require more than one particle. One might think one would require particles to model the measuring apparatuses as well as having at least two entangled particles. However, since the position of a particle is reified in the MIW approach, similarly to Bohm's theory [4,5], a particle can always represent its own measured position (this differs from its momentum, as any measurement of momentum will rely on the action of the quantum potential in Bohm's theory, or the interworld interaction potential in our approach). Moreover, if we consider a particle on a spring, with an externally controllable spring constant k, then either its initial position or its initial momentum can be encoded in position at a nominated time. This is because such evolution turns momentum information into position information at time T/4 = (π/2) √ m/k, and back into position information at time T/2. (Alternatively, the same can be achieved even by free evolution with an externally controlled mass m, an assumption which does not violate any of the MIW framework. By reducing the mass, the time taken for the momentum to be encoded into the position can be made arbitrarily small; by increasing the mass, the position remains as it was initially to an arbitrarily good approximation. Thus, either position or momentum could be measured at a nominated time.) In this way, nonlocal correlations can be explored with a world of just two 1D particles. This scenario could be used not only to demonstrate the EPR paradox, but also Bell nonlocality, as proven by Bell himself in 1986 [37].
Recently, we have generalised MIW for a 2D particle, which is equivalent to two 1D particles, with some numerical success [38]. However, that study was restricted to finding ground states-stationary states with no nodes in the wavefunction. The same paper did consider finding the first excited states for a 1D particle, but only for symmetric potentials for which the node at x = 0 can be put in by hand. This restriction was necessary because, it is now clear, the toy MIW model of Reference [29] cannot reproduce stationary excited states at all. More generally, it cannot be expected to quantitatively reproduce dynamical quantum evolution in which nodes appear and disappear. In the proposal of Bell [37], to violate a Bell inequality by free evolution of two entangled particles, the required initial wavefunction has two nodes and can be expected to develop more during evolution. Thus, to study nonlocality in the MIW approach we certainly need to go beyond the toy model of Reference [29]. Progress in this direction is the topic of the next two sections.

MIW Beyond the Toy Model
The MIW toy model, for a nonrelativistic 1D quantum particle of mass m moving in a potential V(q), has N worlds each containing a 1D particle. Denoting the positions and momenta of these world-particles by X = (x 1 , x 2 , . . . , x N ) and P = (p 1 , p 2 , . . . , p N ), their motion is described in the toy model by a Hamiltonian of the form [29] H(X, P) = The first term is just the sum of N classical Hamiltonians, one for each world. It is the second term, a potential energy responsible for interactions between the worlds, which accounts for quantum phenomena: For evaluation purposes we formally define x 0 := −∞ and x N+1 = ∞, and the ordering x 1 < x 2 < · · · < x N has been assumed (this ordering is preserved under time-evolution, due to the repulsive nature of the potential). The motion of each world-particle is determined by the usual Hamiltonian equations of motion, i.e., Note that the interworld potential U toy (X) vanishes in the classical limith = 0, and also for the case N = 1. In either of these cases, each world evolves independently, according to Newton's laws. More generally, however, forh = 0 and N > 1 the interworld potential U toy (X) leads to forces on each world that act to reproduce quantum phenomena such as Ehrenfest's theorem, spreading of wave packets, tunneling through a barrier, and interference effects [29]. For the case of a 1D oscillator, V(q) = 1 2 mω 2 q 2 , it has further been shown, in the limit N → ∞, that the average energy per world of the MIW ground state converges to the quantum groundstate energy 1 2h ω [29], and that the corresponding stationary distribution of worlds samples the usual quantum Gaussian probability distribution [29,39].
The form of U toy (X) in Equation (4) above is a sum of three-body terms, leading to a force on the n-th world that depends on the positions of the two neighbouring worlds on either side. However, while this form is sufficient to reproduce the quantum phenomena noted above, we have found that it is too simple to model the behaviour of quantum wave functions with nodes. Hence, as noted in the previous section, we must turn to more complex forms of the interworld potential, involving interaction between greater numbers of neighbouring worlds.
Fortunately, there is a great deal of freedom in choosing this potential in the MIW approach. It is possible this freedom could be curtailed via suitable physically-motivated axioms (such as, for example, requiring Ehrenfest's theorem to hold). Here, however, we take a nonaxiomatic approach, to show how the interworld potential of the toy model can be straightforwardly generalised to allow direct interactions between an arbitrary number of worlds, in a manner corresponding to greater accuracy in approximation of the Bohmian acceleration. This greater accuracy supports a corresponding expectation of being able to successfully model wave function nodes.
In particular, it was shown in section II.D of Reference [29] that a suitable general form for the interworld potential in the 1D case is obtained by replacing U toy (X) in the Hamiltonian (3) by where P n and P n are approximations to a probability density P(q) at q = x n , and where P(q) corresponds to some smoothing of the empirical distribution of the world positions {x 1 , x 2 , . . . , x N }. Equation (6) may be regarded as an approximation of the quantum potential in Bohm's theory [4,5] (for which P(q) = |Ψ(q)| 2 ), and hence is expected to reproduce quantum evolution more and more closely as this approximation is improved. The toy model potential in Equation (4) is the simplest such approximation [29]. To go to higher order approximations, we tried two different methods. The first we call the rational smoothing method and the second the equivariance method. In rational smoothing we take a systematic approach to approximate P(q) and its derivative by ratios of polynomials, where the order of these polynomials determines the number of neighbours that each world directly interacts with. This approach is developed in the next subsection, followed by an example. It will be applied to the description of nodes in Section 6.2. rational smoothing is less computational resource-intensive compared to the equivariance Method. The latter is discussed in Appendix A.

Constructing Generalised Interworld Potentials
If P(q) is some smooth probability density as above, which approximates the distribution of world positions, {x 1 , x 2 , . . . , x N }, then the cumulative distribution C(x) := x −∞ dq P(q) must satisfy (at least approximately) Thus, the area under P(q) is divided into N neighbourhoods, each containing one world and having area 1/N. It will be assumed for simplicity that P(q) is nonzero almost everywhere, implying that C(x) is strictly monotonic and hence invertible.
To systematically approximate the quantum force at x n to a given accuracy, we first define the inverse of the cumulative distribution by y(u) := C −1 (u). It follows immediately from Equation (7) that Further, differentiating u = C(y) with respect to y and using d/dy ≡ (1/y )(d/du), gives P(y) = 1/y , P (y) = −y /(y ) 3 , and hence that P (y) P(y) = − y (y ) 2 .
Our aim is now to approximate this last expression by a function of the world positions, and then substitute this approximation into the right hand side of Equation (6) to obtain a corresponding form for the interworld potential U(X).
In particular, expanding y(u n+c ) in a Taylor series about y(u n ) gives the approximation n denotes the lth derivative of y(u) at u = u n . The first L derivatives of y can therefore be approximated to this accuracy by choosing a set of coefficients {α cl } such that ∑ c α cl c l = l! δ ll , l, l = 1, 2, . . . , L.
It is shown how to construct suitable {α cl } below. Equations (10) and (11) immediately yield the approximations accurate to O(1/N L ), for l = 1, 2, . . . , L. Finally, substitution into Equation (9) and hence into Equation (6) gives the corresponding interworld potential, Note that this interworld potential is translation invariant, and scales as 1/λ 2 under x n → λx n . It follows that analogues of Ehrenfest's theorem and wavepacket spreading hold for all such potentials [29].
An advantage of Hamiltonian based methods is that we can use a symplectic numerical integrator to make sure that the total energy is conserved [40]. The quantum force for the mth world-particle can be calculated from Equations (5) and (13) as:

Examples
To obtain an explicit example of the interworld potential in Equation (13), let M denote an L × C matrix with coefficients M lc = c l , and A denote the C × L matrix with coefficients α cl , where c ranges over some set of C integers. Equation (11) can then be written in the matrix form The existence of a solution requires C ≥ L. Further, since the values of α c1 and α c2 are required in Equation (13) It follows that the simplest interworld potential constructed in this way corresponds to C = L = 2. Labelling the two values of c by c = ±1 corresponds, via Equation (12), to approximating the derivatives of y(x n ) via the values of the nearest neighbours x n±1 . Ordering the values of c as −1, 1, the corresponding coefficients {α cl } in Equation (13) are then given by Surprisingly, this is not actually equivalent to the toy model of Reference [29], even though it involves the same number of neighbours in the potential.
The simplest higher-order interworld potential corresponds to C = L = 3, but odd values of C necessarily introduce an unphysical left-right asymmetry, with each world being coupled to different numbers of neighbouring worlds on either side via U α,L,n (X). To preserve symmetry we therefore next consider the case C = L = 4. Labelling the four values of c by c = ±1, ±2, the derivatives of y(x n ) in Equation (12) can be calculated with values of the nearest and next nearest neighbours x n±1 and x n±2 . If we order the values of c as −2, −1, 1, 2, the coefficients {α cl } in Equation (13) are then given by

Numerical Results
To test our higher-order methods against the original toy model, we apply them to the problem of finding the ground and the first excited state of a harmonic oscillator. It is the latter, containing a node, where the toy model fails and it is necessary to use a higher-order interworld potential. In this section we use dimensionless configuration coordinates and dimensionless times T := ωt/2π, (19) where ω is the harmonic oscillator frequency. Figure 1b shows the result of applying the toy model in Equations (3) and (4) to the ground state of a harmonic oscillator, corresponding to V(q) = 1 2 mω 2 x 2 . In this test we distributed 50 world-particles with x n determined by inverting Equation (7) for the groundstate probability density P (0) t (q) (Figure 1a). We then evolved these under the single-particle harmonic oscillator potential 1 2 mω 2 x 2 n and the quantum interworld potential. We see that as expected for a stationary quantum potential, the classical and quantum forces cancel each other and the world-particles stay stationary. The slight oscillations for the world-particles near the boundary is due to differences between the toy-model and Bohmian potentials in areas with high curvature and low sampling. These differences imply an exact stationary state for the MIW potential that has slightly different values for the x n for any finite value of N [29].

Toy Model
However, if we distribute world-particles based on the probability density of the first excited state of a harmonic oscillator (Figure 2a) and apply the same model, the world-particles will not stay stationary as expected (Figure 2b). We also implemented the nearest-neighbour interworld potential defined by Equation (17), and found very similar behaviour to the original toy model. That is, it fails to support stationary configurations corresponding to excited energy eigenstates. The unsuccessful simulation results in Figure 3 show that the problem of the nodes cannot be fixed by increasing the sampling. In the case of a node, because the probability density is zero, there is always a second order curve between the two world-particles adjacent to the node that can never be correctly estimated by toy model approximation. This leads to a poor approximation of the Bohmian force for the world-particles near the node and instead of staying stationary the world-particles move towards the gap and fill the gap.
The general problem with areas of low probability density, and in particular in the region of the node of the quantum state at x = 0 in Figure 2a, is that sampling is low and nearest-neighbour approximations are not valid for calculating the quantum potential. If the probability is low but non-zero, theoretically, we can increase the number of worlds until we reach a good sampling in those areas. To test the simulation with higher sampling, we tried 5000 worlds. To reduce the simulation time, we only focused on the 20 world-particles around the node (Figure 3a). To apply the correct boundary condition for truncated area, we kept five world-particles near each boundary artificially fixed.

Higher Order Potential
To test a higher order approximation, we applied the rational smoothing model in Equation (13) for the case L = 4 (equivalent to a 5-world interaction), to the first excited state of a harmonic oscillator with 5000 worlds. This corresponds to the example in Equation (17).
In the first test we kept all the worlds artificially stationary except for the two worlds around the node. Figure 4 shows that these two middle worlds stay stationary which means that the quantum potential for those is well approximated. To apply the boundary condition, the rest of the world-particles were kept fixed. For the evolution, rational smoothing with L = 4 (equivalent to 5-world approximation) is used. Time step is 10 −9 . The two world-particles stay stationary as expected.
In the second test we tried the same scenario but this time simulated 10 worlds around the node and kept the rest stationary. Figure 5 shows that the simulation was successful in the sense that world-particles stayed close to the starting position and did not move towards the gap i.e., the low density region, corresponding to a node of the associated excited state, is preserved by the interaction, and the evolution is approximately stationary. The obvious difference compared to Figure 4 is the oscillations around the starting point. The oscillations appear chaotic, as might be expected for a system of nonlinearly coupled harmonic oscillators. These oscillations might also be due to the time step being too long. We were not able to test this conjecture because the test in Figure 5 took a few days on our desktop computer and, due to time constraints, we could not run it with smaller time steps.
We repeated the same test with 7-worlds approximation (L = 6). Figure 6 shows that this higher order approximation decreases oscillation compared to the results in Figure 5. Hence, the simulation of a wave function with nodes, which failed for the toy model (Figures 2 and 3), is seen to become more and more accurately modelled in the MIW approach as the number of directly interacting worlds is increased (Figures 5 and 6). Similar convergence might also be possible with the equivariance Method presented in Appendix A. However, the results of the simplest test, with a potential involving five interacting worlds, were not as positive as those of the rational smoothing explored in this section. They are also more numerically intensive, so we did not pursue it further.

Conclusions
Bell's theorem of 1964 showed that any deterministic interpretation of quantum mechanics must be nonlocal (Section 2). Bohm's theory of 1952 is the example par excellence of a nonlocal, deterministic theory (Section 3). The theoretical approach we introduced in 2014-Many Interacting Worlds-is also deterministic and, if it is to succeed in replicating Bohmian mechanics (and thus all quantum phenomena), must be nonlocal (Section 4). It is conceptually clear how an interworld potential can lead to nonlocality in Bell's 1964 sense, and, perhaps surprisingly, this could in principle be simulated with a universe of only two particles, in 1D (one spatial dimension), if an externally controlled spatially localized potential is allowed, as discussed in Section 4.
Unfortunately, this proof-of-principle simulation of Bell-nonlocality cannot be done using the toy model for an interworld potential introduced in Reference [29]. The reason is that it cannot deal with wavefunctions having nodes, and nodes are certainly necessary for modelling the entangled states and measurements necessary for violating a Bell inequality [37]. When one prepares a distribution of worlds, in one dimension, corresponding to a stable excited quantum state, the dynamics of the toy model causes the gap between worlds where the node should collapse, and no stationary configuration is reached. Dealing with nodes is a problem in many quantum simulation methods based on Bohmian mechanics [41]. Nodes should not be a problem for interpretations involving a continuum of worlds [32][33][34][35], as they are formulated to be exactly equivalent to quantum mechanics. However, as remarked in Section 3, our view is that these interpretations do not solve the conceptual problems of the Everettian many-worlds interpretation.
Here we showed that this problem of nodes in our discrete MIW approach may not be fundamental, but rather may be a result of using a too simple form for the inter-world interaction potential. By using a higher-order approximation to define our interworld potential from Bohm's quantum potential, we were able to show that a gap in the configuration of worlds, corresponding to the node of the first excited harmonic oscillator energy eigenstate, can remain open for at least a full harmonic oscillator period, and perhaps indefinitely. The world configurations were not exactly stationary, but rather had high frequency irregular oscillations. However, by increasing the order of the approximation, the size of the oscillations could be reduced.
Our simulations considered only the dynamics near the node, with more distant worlds artificially held fixed. Whether the node would remain stable if all worlds were allowed to evolve according to the MIW dynamics is an open question. In addition, our simulations were restricted to one particle in 1D. The MIW approach has been successfully used to simulate one particle in 2D (or, equivalently, two particles in 1D) [38,42] and 3D [42], but only to find the ground state configuration. Combining these research directions to be able to simulate stable excited states for two particles in 1D, and beyond, is a challenge for future work. Finally, to realise a simulation of the Bell experiment described in Section 4 would require simulating not just stationary nodes but also dynamical nodes, that may appear and then disappear in an instant, which are notoriously difficult to deal with in Bohmian-inspired numerical approaches [41]. Thus, much work remains to be done, but the positive results reported here are encouraging.

Abbreviations
The following abbreviations are used in this manuscript: MIW Many Interacting Worlds EPR Einstein, Podolsky, and Rosen 1D One (spatial) dimension 2D Two (spatial) dimensions 3D Three (spatial) dimensions

Appendix A. A Equivariance Method
Here we introduce an alternative higher order method to rational smoothing, motivated by the equivariance property of Bohmian mechanics [11]. In particular, for a given configuration of worlds, {x 1 , x 2 , . . . , x n }, we construct a smooth polynomial probability density, P n (x), in the region of each world x n , and use this to approximate the Bohmian potential. The coefficients of the polynomial are determined by requiring that area under P n (x) between x n and x n+1 is equal to the constant value 1/N for each world, analogously to Equation (7). The accuracy of this equal-probability or "equivariance" method will increase with the degree of the polynomial.
Here we will illustrate the equivariance method for the case where P n (x) is third-order, corresponding to direct interactions between sets of five adjacent worlds. Thus, P n (x) = a n + b n x + c n x 2 + d n x 3 . (A1) The equivariance method then requires that x n+1 x n P n (x)dx = const. = N −1 .
To evaluate the quantum force, we use an interworld interaction potential of the form of Equation (6) and P n as in Equation (A1) (with a n , b n , c n and d n taken from (A10)) to give U =h 2 8m N ∑ n=1 b n + 2c n x n + 3d n x 2 n a n + b n x n + c n x 2 n + d n x 3 n 2 . (A11) This may be compared with the corresponding interaction potential U α,L in Equation (13) obtained via rational smoothing. The quantum force on each particle is evaluated similarly to Equation (14). We used Matlab to evaluate the analytical derivatives of these terms for the simulation. The resulting equation is too lengthy to include here (c. 360,000 characters for N = 5). Figure A1 shows the result of applying the equivariance method to the first neighbours of the node in the first excited state of a harmonic oscillator. It shows that these world-particles do not move into the gap from their original position, but rather undergo some oscillatory motion. This is poorer behaviour than the corresponding rational smoothing simulation shown in Figure 4. The simulations were also slower because of the complexity of the analytical form of the force law, mentioned above. For these reasons we have not pursued this method further. Figure A1. The first excited state of a harmonic oscillator is simulated using 5-world approximation in equivariance Mmethod. Five-thousand worlds are used and only the two world-particles next to the node are simulated. The rest are kept stationary, similar to Figure 4 for the rational smoothing case.