Original Research ARTICLE
Discrete-event simulation of uncertainty in single-neutron experiments
- 1Department of Applied Physics, Zernike Institute for Advanced Materials, University of Groningen, Groningen, Netherlands
- 2Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, Jülich, Germany
- 3Department of Physics, RWTH Aachen University, Aachen, Germany
A discrete-event simulation approach which provides a cause-and-effect description of many experiments with photons and neutrons exhibiting interference and entanglement is applied to a recent single-neutron experiment that tests (generalizations of) Heisenberg's uncertainty relation. The event-based simulation algorithm reproduces the results of the quantum theoretical description of the experiment but does not require the knowledge of the solution of a wave equation nor does it rely on concepts of quantum theory. In particular, the data satisfies uncertainty relations derived in the context of quantum theory.
Quantum theory has proven extraordinary powerful for describing a vast number of laboratory experiments. The mathematical framework of quantum theory allows for a straightforward (at least in principle) calculation of numbers which can be compared with experimental data as long as these numbers refer to statistical averages of measured quantities, such as for example atomic spectra, the specific heat and magnetic susceptibility of solids. However, as soon as an experiment is able to record the individual clicks of a detector which contribute to the statistical average, a fundamental problem appears. Although quantum theory provides a recipe to compute the frequencies for observing events, it does not account for the observation of the individual events themselves [1–4]. Prime examples are the single-electron two-slit experiment , single-neutron interferometry experiments  and optics experiments in which the click of a detector is assumed to correspond to the arrival of a single photon .
From the viewpoint of quantum theory, the central issue is how it can be that experiments yield definite answers. On the other hand, it is our brain which decides, based on what it perceives through our senses and cognitive capabilities, what a definite answer is and what it is not. According to Bohr  “Physics is to be regarded not so much as the study of something a priori given, but rather as the development of methods of ordering and surveying human experience. In this respect our task must be to account for such experience in a manner independent of individual subjective judgment and therefore objective in the sense that it can be unambiguously communicated in ordinary human language.” This quote may be read as a suggestion to construct a description in terms of events, some of which are directly related to human experience, and the cause-and-effect relations among them. Such an event-based description obviously yields definite answers and if it reproduces the statistical results of experiments, it also provides a description on a level to which quantum theory has no access.
For many interference and entanglement phenomena observed in optics and neutron experiments, such an event-based description has already been constructed, see Michielsen et al. , De Raedt et al. , De Raedt et al.  for recent reviews. The event-based simulation models reproduce the statistical distributions of quantum theory without solving a wave equation but by modeling physical phenomena as a chronological sequence of events. Hereby events can be actions of an experimenter, particle emissions by a source, signal generations by a detector, interactions of a particle with a material and so on [9–11].
The basic premise of our event-based simulation approach is that current scientific knowledge derives from the discrete events which are observed in laboratory experiments and from relations between those events. Hence, the event-based simulation approach is concerned with how we can model these experimental observations but not with what “really” happens in Nature. This underlying premise strongly differs from the assumption that the observed events are signatures of an underlying objective reality which is mathematical in nature but is in line with Bohr's viewpoint expressed in the above quote.
The general idea of the event-based simulation method is that simple rules define discrete-event processes which may lead to the behavior that is observed in experiments. The basic strategy in designing these rules is to carefully examine the experimental procedure and to devise rules such that they produce the same kind of data as those recorded in experiment, while avoiding the trap of simulating thought experiments that are difficult to realize in the laboratory. Evidently, mainly because of the lack of knowledge, the rules are not unique. Hence, it makes sense to use the simplest rules one could think of until a new experiment indicates that the rules should be modified. The method may be considered as entirely “classical” since it only uses concepts which are directly related to our perception of the macroscopic world but the rules themselves are not necessarily those of classical Newtonian dynamics.
The event-based approach has successfully been used for discrete-event simulations of quantum optics experiments such as the single beam splitter and Mach-Zehnder interferometer experiments, Wheeler's delayed choice experiments, a quantum eraser experiment, two-beam single-photon interference experiments and the single-photon interference experiment with a Fresnel biprism, Hanbury Brown-Twiss experiments, Einstein-Podolsky-Rosen-Bohm (EPRB) experiments, and of conventional optics problems such as the propagation of electromagnetic plane waves through homogeneous thin films and stratified media, see Michielsen et al. , De Raedt et al.  and references therein. For applications to single-neutron interferometry experiments see De Raedt et al. [10, 11]. The same methodology has also been employed to perform discrete-event simulations of quantum cryptography protocols  and universal quantum computation .
In this paper, we extend this list by demonstrating that the same approach provides an event-by-event description of recent neutron experiments [14, 15] devised to test (generalizations of) Heisenberg's uncertainty principle. It is shown that the event-by-event simulation generates data which complies with the quantum theoretical description of this experiment. Therefore, these data also satisfy the inequalities which, in quantum theory, express (generalizations of) Heisenberg's uncertainty principle. However, as the event-by-event simulation does not resort to concepts of quantum theory these findings indicate that there is little intrinsically “quantum mechanical” to these inequalities, in concert with the idea that quantum theory can be cast into a “classical” statistical theory[16–28].
2. Experiment and Quantum Theoretical Description
A block diagram of the neutron experiment designed to test uncertainty relations [14, 15] is shown in Figure 1. We now describe this experiment in operational terms and as we go along, we also give the quantum theoretical description in terms of spin 1/2 particles such as neutrons.
Figure 1. Functional block diagram of the neutron experiment [14, 15] to testuncertainty relations, see also Figure 2 in Erhart et al. , Sulyok et al. . Monochromatic neutrons enter from the left. SA1, SA2, SA3: spin analyzers (see text); SF1, SF2, SF3, SF4: spin flippers (see text). The positions of SF1 and the pair (SF2,SF3) are variable. The preparation stage 1 manipulates the magnetic moments of the neutrons such that, depending on the orientation of SA1 and the position of SF1, they leave stage 1 with their moments directed along the z-axis or along a fixed direction in the x − y plane. In the quantum theoretical description (see text), stages 2, 3 perform the successive measurement of the eigenvalues of the spin operators σϕ = σx cos ϕ + σy sin ϕ and σy, respectively. These seven devices, in combination with the static magnetic field B0z and the variable positions of SF1 and (SF2,SF3), are sufficient to determine the expectation values that appear in Heisenberg-like uncertainty relations (see text).
Conceptually, the neutron experiment [14, 15] exploits two different physical phenomena: the motion of a magnetic moment in a static magnetic field and a spin-analyzer that performs a Stern-Gerlach-like selection of the neutrons based on the direction of their magnetic moments.
A magnetic moment S in an external, static magnetic field Be experiences a rotation about the direction of the unit vector e. The unitary transformation that corresponds to such a rotation is given by
where γ is the gyromagnetic ration of the particle, t is the time that the particle interacts with the magnetic field, the variable φ = γtB is introduced as a shorthand for the angle of rotation, and σ = (σx, σy, σz) denote the Pauli-spin matrices.
The spin analyzer acts as a projector. It is a straightforward (see pages 172 and 250 in Ballentine, see 3) to show that within quantum theory, an ideal spin analyzer directed along the unit vector n is represented by the projection operator
Using Equations (1), (2), it is straightforward to construct the quantum theoretical description of each of the three stages in the experimental setup.
Stage 1. The purpose of the first spin analyzer (SA1) is to prepare neutrons with their magnetic moments in the direction of the static magnetic field B0z. Then, the particle travels for some time in a region where the field B0z is present but as its magnetic moment is aligned along z, the magnetic moment does not rotate.
As will become clear later, to test Ozawa's inequality [14, 15, 29], it is necessary to be able to prepare initial states in which the magnetic moment lies in the x − y plane. In the experiment, this is accomplished by putting in place, the spin flipper SF1. The spin flipper (SF1), in essence a static magnetic field aligned along the x-direction, rotates the magnetization about the x-axis by an amount proportional to the magnetic field. For simplicity, it is assumed that this rotation changes the direction of the magnetic moment from z to y [14, 15]. The position of SF1, relative to the direction of flight of the neutrons, is variable. By moving SF1 one can change the time that the particles perform rotations about the z-axis, hence one can control the direction of the magnetic moment in the x − y plane as the particle leaves stage 1.
Quantum theoretically, the action of the components of stage 1 is described by the product of unitary matrices
where θ0 = π/2 or θ0 = 0 if SF1 is in place or not and θ1 is the variable (through the variable position of SF1) rotation angle, the value of which will be fixed later. Obviously, in the case that SF1 is not present, because the incoming neutrons have their moments aligned along the z-direction, these moments do not perform rotations at all.
Stage 2. This stage consists of a pair of spin flippers (SF2,SF3) and a spin analyzer (SA2). The position of (SF2,SF3), relative to the direction of flight of the neutrons, is variable whereas the position of SA2 is fixed. The action of the components of stage 2 is described by the product of matrices
where, as a consequence of the variable position of (SF2,SF3), the rotation angles θ2, θ3, and θ4 change with the position of (SF2,SF3). The value of variable S1 = ±1 labels one of the two possible alignments of the spin polarizer along z. Note that because of the projection M(S1, z), the matrix T2 is not unitary.
Stage 3. The final stage consists of a spin flipper SF4 and a spin analyzer SA3. The time evolution of the magnetic moment as it traverses this stage is given by
where θ5 is a fixed rotation angle. The value of variable S2 = ±1 labels one of the two possible alignments of the spin polarizer along z. The matrix T3 is not unitary.
According to the postulates of quantum theory, the probability to detect a neutron leaving stage 3 is given by Ballentine 
where it is assumed that the detector simply counts every impinging neutron (which in view of the very high detection efficiency is a very good approximation, see Rauch and Werner ). In Equation (6), the initial state of the S = 1/2 quantum system is represented by the density matrix
The real-valued vector a (ǁaǁ ≤ 1) completely determines the initial state of the quantum system. Using Equations (3–7) we find
independent of θ3 and θ5. Recall that by construction of the experimental setup θ2 + θ4 is fixed. We can make contact to the expressions used in the analysis of the neutron experiment [14, 15], by substituting θ1 = 0, θ2 = ϕ + π/2 and θ4 = −ϕ − π/2 where ϕ is called the detuning angle [14, 15]. We obtain
As explained in detail in subsection 2.1, the experimental setup can be interpreted as performing successive measurements of the operators σϕ = σx cos ϕ + σy sin ϕ and σy, their eigenvalues being S1 and S2, respectively. Note that these two operators do not commute unless cos ϕ = 0 and that the observed eigenvalues S1 and S2 of these two operators are correlated, as is evident from the contribution S1S2 sin ϕ in Equation (9).
From Equation (9), the expectation values of the various spin operators follow immediately. Specifically, we have
and as σ2ϕ = σ2y = 1, the variances 〈σ2ϕ〉a − 〈σϕ〉2a and 〈σ2y〉a − 〈σy〉2a are completely determined by Equation (10).
2.1. Filtering-Type Measurements of One Spin-1/2 Particle
The neutron experiment [14, 15] can be viewed as a particular realization of a filtering-type experiment [3, 30]. The layout of such an experiment is shown in Figure 2. According to this diagram, the experiment consists of performing successive Stern-Gerlach-type measurements on one spin-1/2 particle at a time. Conceptually, assuming a stationary particle source, the neutron experiment [14, 15] and the filtering-type experiment shown in Figure 2 are identical, see also Figure 1 in Erhart et al. , Sulyok et al. . In practice, the difference between the filtering-type experiment and the neutron experiment [14, 15] is that in the latter four experiments (labeled by the variables S1 = ±1 and S2 = ±1), are required for each of the two opposite orientations of the two spin analyzers SA2 and SA3 whereas the setup depicted in Figure 2 directly yields the results of the four separate runs.
Figure 2. Diagram of the filtering-type experiment of which the neutron experiment [14, 15] is a special case. Spin-1/2 particles pass through a Stern-Gerlach magnet M0 that projects the spin onto either the b direction or the -b direction. In case of the former (latter) projection, the particle is directed to the Stern-Gerlach magnet M1 (M2). M1 and M2 are assumed to be identical and project the spin onto either the c direction or the -c direction. A “click” of one of the four detectors D+1,1, D−1,1, D+1,2, and D−1,2 signals the arrival of a particle.
The pictorial description of the filtering experiment goes as follows. A particle enters the Stern-Gerlach magnet M0, with its magnetic field along direction b. M0 “sends" the particle either to Stern-Gerlach magnet M1 or M2. The magnets M1 and M2, identical and both with their magnetic field along direction c, redirect the particle once more and finally, the particle is registered by one of the four detectors D+1,1, D−1,1, D+1,2, and D−1,2. This scenario is repeated until the statistical fluctuations of the counts of the four detectors are considered to be sufficiently small.
We label the particles by a subscript α. After the αth particle leaves M1 or M2, it will hit one (but only one) of the four detectors. We assume ideal experiments, that is at any time one and only one out of four detectors fires. We write x(i,j)α = 1 with i = ±1 and j = 1, 2 if the αth particle was detected by detector Di,j and x(i,j)α = 0 otherwise.
We define two new dichotomic variables by
Note that for each incoming particle, only one of the detectors clicks hence only one of the x(i,j)α's is nonzero or, equivalently x(+1, 1)α + x(−1, 1)α + x(+1, 2)α + x(−1, 2)α = 1.
In the quantum theoretical description of the filtering experiment, if S1,α = ±1, the spin has been projected on the ±b direction. Likewise, if S2,α = ±1, the spin has been projected on the ±c direction. In other words, S1,α and S2,α are the eigenvalues of the spin operator projected on the directions b and c, respectively. Then, according to quantum theory, the probability to observe a pair of eigenvalues (S1, S2) is given by Ballentine , De Raedt et al. .
where the state ρ is given by Equation (7) and the M's denote projection operators. It is easy to see that Equation (9) is a particular case of Equation (12).
Note that [M(S1, b), M(S2, c)] ≠ 0 unless b = ±c, [ρ, M(S1, b)] ≠ 0 unless a = ±b, and that [ρ, M(S2, c)] ≠ 0 unless a = ±c. Thus, for virtually all cases of interest, none of the operators in Equation (12) commute, yet quantum theory yields the probability P(S1, S2|a) for all cases. Clearly, the statement that one can determine the eigenvalues of two non-commuting operators in one experiment contradicts the conventional teaching that non-commuting operators cannot be diagonalized simultaneously and therefore cannot be measured simultaneously. The reason for this apparent contradiction is the hidden assumption that diagonalization and the act of measurement in a laboratory (i.e., a click of the detector) are equivalent in some sense. The filtering-type experiment is a clear example which shows that they are not: according to quantum theory the eigenvalues S1 and S2 of the operators σ·b and σ·c, respectively can always be measured simultaneously even though these operators cannot always be diagonalized simultaneously.
3. Event-By-Event Simulation
A minimal, discrete-event simulation model of single-neutron experiments requires a specification of the information carried by the neutrons, of the algorithm that simulates the source and the devices used in the experimental setup (see Figure 1), and of the procedure to analyze the data.
• Messenger: A neutron is regarded as a messenger carrying a message. In principle, there is a lot of freedom to specify the content of the message, the only criterion being that in the end, the simulation should reproduce the results of real laboratory experiments. We adopt Occam's razor as a guiding principle to determine which kind of data the messenger should carry with it, that is we use the minimal amount of data.
The pictorial description that will be used in the following should not be taken literally: it is only meant to help visualize, in terms of concepts familiar from macroscopic physics, the minimal amount of data the messenger should carry.
Picturing the neutron as a tiny classical magnet we can use the spherical coordinates θ and φ to specify the direction of its magnetic moment
relative to the fixed frame of reference defined by the static magnetic field B0z. The messenger should also be aware of the time it takes to move from one point in space to another. The time of flight and the direction of the magnetic moment are conveniently encoded in a message of the type [10, 11]
where ψ(i) = νt + δi, for i = 1, 2 and φ = δ1 − δ2 = ψ(1) − ψ(2). Within the present model, the state of the neutron, that is the message, is completely described by the angles ψ(1), ψ(2) and θ and by rules (to be specified), by which these angles change as the neutron travels through the network of devices. This model suffices to reproduce the results of single-neutron interference and entanglement experiments and of their idealized quantum theoretical descriptions [10, 11].
In specifying the message Equation (14), we exploited the isomorphism between the algebra of Pauli matrices and rotations in three-dimensional space, not because the former connects to quantum mechanics but only because we find this representation most convenient for our simulation work [9–11]. The direction of the magnetic moment follows from Equation (14) through
A messenger with message u at time t and position r that travels with velocity v, along the direction q during a time interval t′ − t, changes its message according to ψ(i)← ψ(i) + ν(t′ − t) for i = 1, 2, where ν is an angular frequency which is characteristic for a neutron that moves with a fixed velocity v. In a monochromatic beam of neutrons, all neutrons have the same value of ν .
In the presence of a magnetic field B = (Bx, By, Bz), the magnetic moment rotates about the direction of B according to the classical equation of motion
Hence, as the messenger passes a region in which a magnetic field is present, the message u changes into the message
where g denotes the neutron g-factor, μN the nuclear magneton, T the time during which the messenger experiences the magnetic field.
In the event-based simulation of the experiment shown in Figure 1, the time-of-flight T determines the angle of rotation of the magnetic moment through Equation (17) and can, so to speak, be eliminated by expressing all operations in terms of rotation angles. However, this simplification is no longer possible in the event-based simulation of single-neutron interference and entanglement experiments [10, 11].
• Source: When the source creates a messenger, its message needs to be initialized. This means that three angles ψ(1), ψ(2) and θ need to be specified. In practice, instead of implementing stage 1 it is more efficient to prepare the messengers such that the corresponding magnetic moments are along a specified, fixed direction. For instance, to mimic fully coherent, spin-polarized neutrons that enter stage 2 with their spin along the x-axis, the source would create messengers with θ = π/2, and without loss of generality, ψ(1) = ψ(2) = 0.
• Spin-flipper: The spin-flipper rotates the magnetic moment by an angle π/2 about the x axis.
• Spin analyzer: This device shares with the magnet used in Stern-Gerlach experiments the property that it diverts incoming particles in directions which depend on the orientation of their magnetic moments relative to the magnetic field inside the device. For appropriate choices of the experimental parameters, the Stern-Gerlach magnet splits the incoming beam of particles into spatially separated beams. If there are two well-separated beams, the action of the device is to align the magnetic moments of the incoming particles either parallel or anti-parallel to the direction of the field.
Ignoring all the details of interaction of the magnetic moments with the Stern-Gerlach magnet, the operation of separating the incoming beam into spatially separated beams is captured by the very simple probabilistic model defined by
where x = − 1, 1 labels the two distinct spatial directions, 0 < r < 1 is a uniform pseudo-random number, Θ(x) is the unit step function and, as before, S = ±1 labels the orientation of the spin analyzer. For each incoming messenger, a new pseudo-random number is being generated. Recall that |mz| ≤ 1 (see Equation (13)) hence the first term of Equation (18) is a number between zero and two. If we would set mz = 〈σ · n〉, the model defined by Equation (18) would generate minus- and plus ones according to the projection operator Equation (2). Applied to the neutron experiments [14, 15], the function of the spin analyzer is to pass particles with say, spin-up, only. In the simulation model this translates to letting the messenger pass if x = 1 and destroy the messenger if x = − 1.
It is of interest to explore the possibility whether different models for the spin analyzer can yield averages over many events which cannot be distinguished from the results obtained by employing the probabilistic model Equation (18). As the extreme opposite to the probabilistic model, we consider a deterministic learning machine (DLM) defined by the update rule
where x = − 1, 1 labels the two distinct spatial directions and − 1 ≤ u ≤ 1 encodes the internal state of the DLM (the equivalent of the seed of the pseudo-random number generator). The parameter 0 ≤ γ < 1 controls the pace at which the DLM learns the value mzS. The properties of the time series of x's generated by the rules Equation (19) have been scrutinized in great detail elsewhere, see Michielsen et al.  and references therein. Suffice it to say here that for many events, the average of the x's converges to mzS and that the x's are highly correlated in time.
Obviously, the DLM-based model is extremely simple and fully deterministic. It may easily be rejected as a viable candidate model by comparing the correlations in experimentally observed time series with those generated by Equation (19). However, if the experiment only provides data about average quantities, there is no way to rule out the DLM model. Unfortunately, the neutron experiments [14, 15] do not provide the data necessary to reject the DLM model, simply because the spin analyzer passes particles with say, spin-up, only. Hence there is no way to compute time correlations. Although we certainly do not want to suggest that the spin analyzers used in the experiments behave in the extreme deterministic manner as described by Equation (19), it is of interest to test whether such a simple deterministic model can reproduce the averages computed from quantum theory and also obeys the same uncertainty relations as the genuine quantum theoretical model.
• Detector: As a messenger enters the detector, the detection count is increased by one and the messenger is destroyed. The detector counts all incoming messengers. Hence, we assume that the detector has a detection efficiency of 100%. This is a good model of real neutron detectors which can have a detection efficiency of 99% or more .
• Simulation procedure and data analysis: First, we establish the correspondence between the initial message uinitial and the description in terms of the density matrix Equation (7). To this end, we remove all devices from stage 1 and 2 and simply count the number of messages that pass SA3 with S2 = 1, for instance. It directly follows from Equation (18) that the relative frequency of counts is given by mz, the projection of the message onto the z-axis. In other words, we would infer from the data that in a quantum theoretical description the z-component of the density matrix az is equal to mz. By performing rotations of the original message it follows by the same argument that a = minitial.
For each pair of settings (S1, S2) of the spin analyzers (SA2,SA3) and each position of the pair of spin flippers (SF2,SF3) represented by a rotation of ϕ about the z-axis (see Section 2), the source sends N messengers through the network of devices shown in Figure 1. The source only creates a new messenger if (i) the previous messenger has been processed by the detector or (ii) the messenger was destroyed by one of the spin analyzers. In other words, direct communication between messengers is excluded. As a device in the network receives a messenger, it processes the message according to the rules specified earlier and sends the messengers with the new message to the next device in the network. If the device is a spin analyzer, it may happen that the messenger is destroyed. The detector counts all messengers that pass SA3 and destroys these messengers.
For a sequence of N messengers all carrying the same initial message a = minitial, this procedure yields a count N(S1, S2|a) (recall that ϕ is fixed during each sequence of N events). Repeating the procedure for the four pairs of settings yields the relative frequencies
Note that the numerator in Equation (20) is not necessarily equal to N because messengers may be destroyed when they enter a spin analyzer. From Equation (20) we compute
• Validation. The event-based model reproduces the results of the quantum theoretical description if, within the usual statistical fluctuations, we find that F(S1, S2|a) ≈ P(S1, S2|a) with P(S1, S2|a) given by Equation (9). This correspondence is most easily established by noting that for fixed ϕ and a, the three expectations Equations (21–23) completely determine Equation (20) and that, likewise, the quantum theoretical distribution Equation (9) is completely determined by the expectations Equations (21–23) with F(S1, S2|a) replaced by P(S1, S2|a). In other words, for the event-based model to reproduce the results of the quantum theoretical description of the neutron experiment [14, 15], it is necessary and sufficient that the simulation results for Equations (21–23) are in agreement with the quantum theoretical results (see Equation (10) 〈S1〉 = ax cos ϕ + ay sin ϕ, 〈S2〉 = sin ϕ 〈S1〉, and 〈S1 S2〉 = sin ϕ). In Figures 3, 4 we present representative results of event-based simulations of the neutron experiment [14, 15], showing that the simulation indeed reproduces the predictions of the quantum theoretical description of the neutron experiment [14, 15]. Comparing Figures 3, 4 we can only conclude that it does not matter whether the model for the spin analyzers uses pseudo-random numbers Equation (18) or is DLM-based Equation (19).
Figure 3. Results for the expectations 〈S1〉 (red solid circles), 〈S2〉 (green solid squares), and 〈S1S2〉 (blue open circles) as obtained by the event-by-event simulation of the neutron experiment [14, 15], using the model Equation (18) for the spin analyzers. The solid lines represent the corresponding quantum theoretical prediction as obtained from Equation (9). Top: incoming particles have magnetization (1, 0, 0). Bottom: incoming particles have magnetization (0, 1, 0). For each pair of settings (S1, S2) of the spin analyzers (SA2,SA3) and each position of the pair of spin flippers (SF2,SF3) represented by a rotation of ϕ about the z-axis, referred to as detuning angle in Erhart et al. , Sulyok et al. , the simulation consists of sending N = 10000 messengers (“neutrons”) into stage 2.
Figure 4. Same as Figure 3 except that the probabilistic model for the spin-analyzer, as specified by Equation (18), has been replaced by a deterministic learning machine with an update rule given by Equation (19). For each pair of settings (S1, S2) of the spin analyzers (SA2,SA3) and each position of the pair of spin flippers (SF2,SF3) represented by a rotation of ϕ about the z-axis, referred to as detuning angle in Erhart et al. 14 (14], Sulyok et al.15 (15], the simulation consists of sending N = 10000 messengers (“neutrons”) into stage 2. The parameter γ = 0.999.
Summarizing: the event-based simulation model of the neutron experiment [14, 15] presented in this section does not rely, in any sense, on concepts of quantum theory yet it reproduces all features of the quantum theoretical description of the experiment. Although the event-based model is classical in nature, it is not classical in the sense that it cannot be described by classical Hamiltonian dynamics.
4. Uncertainty Relations: Theory
The neutron experiment [14, 15] was conceived to test an error-disturbance uncertainty relation proposed by Ozawa . By introducing particular definitions of the measurement error ϵ(A) of an operator A and the disturbance η(B) of an operator B, Ozawa showed that
and MA and MB represent the operators of different measuring devices (implying [MA, MB] = 0) that allow us to read of the value of the measurement of A and B, respectively. Thereby it is assumed that 〈MA〉 = 〈A〉 and 〈MB〉 = 〈B〉, that is that the measurements of A and B are unbiased, implying that ϵ(A) = Δ(MA − A) and η(B) = Δ(MB − B).
In Ozawa's model of the measurement process, the state of the system + measurement devices is represented by a direct product of the wavefunction of the system and the wavefunction of the measurement devices . The operators A and B refer to the dynamical variables of the system while the MA and MB refer to the dynamical variables of two different measurement devices. Furthermore, it is assumed that both the system and the measuring devices (probes) are described by quantum theory, i.e., the time evolution of the whole system is unitary [29, 32]. Although this basic premise is at odds with the fact that experiments yield definite answers [1–3, 33], within the realm of the quantum theoretical model, it “defines” the measurement process, see Allahverdyan et al.  for an extensive review.
Assuming that [C, D] = 0, wehave
or, taking expectation values,
Taking the absolute value of both sides of Equation (31) and using the triangle inequality we find
The derivation of Equation (33) only makes use of the triangle inequality, the notion of a non-negative inner product on a vector space, the Cauchy-Schwarz inequality and the assumption that [C, D] = 0. Therefore Equation (33) is “universally valid” [29, 32, 34] whenever [C, D] = 0.
Inequality Equation (24) directly follows from Equation (33) by substituting C = MA, D = MB and by using the assumption of unbiasedness, meaning 〈MA − A〉 = 0 and 〈MB − B〉 = 0. With the restriction imposed by the assumptions of unbiasedness and [MA, MB] = 0, it is also “universally valid”.
Unlike Equation (24), inequality Equation (34) lacks a mathematical rigorous basis and therefore it is not a surprise that it can be violated . Indeed, the data recorded in the neutron experiment clearly violate Equation (34) [14, 15]. In general, in mathematical probability theory as well as quantum theory, inequalities such as the Craméer-Rao bound , the Robertson inequality , Equations (24, 33) are mathematical identities which result from applications of the Cauchy-Schwarz inequality. Being mathematical identities within the realm of standard arithmetic, they are void of any physical meaning and cannot be violated. Therefore, if an experiment indicates that such an identity (i.e., inequality) might be violated, this can only imply that there is an ambiguity (error) in the mapping between the variables used in the theoretical model and those assigned to the experimental observations [30, 40]. Any other conclusion that is drawn from such a violation cannot be justified on logical/mathematical grounds.
Following Erhart et al. , Sulyok et al. , we assume that the state of the system is represented by the density matrix ρ = |z〉 〈z|, that is the magnetic moment of the neutrons are assumed to be aligned along the z-direction. With A = σx and B = σy we have Erhart et al. , Sulyok et al. 
Note the absence of ħ in Equation (36), in agreement with work that shows that ħ may be eliminated from the basic equations of (low-energy) physics by a re-definition of the units of mass, time, etc. [41, 42].
Conceptually, the application of Equation (24) to the neutron experiment [14, 15] is not as straightforward as it may seem. In a strict sense, in the neutron experiment [14, 15], there are no measurements of the kind envisaged in Ozawa's measurement model. This is most obvious from the quantum theoretical description of the experiment given in Section 2: for fixed S1 and S2, the relative frequency of detector counts is given by Equation (9), and “noise” caused by “probes” does not enter the description. Indeed, from the expressions of ϵ2(A) and η2(B) in terms of spin operators, see Equation (35), it is immediately clear that in order to determine ϵ2(A) and η2(B), there is no need to actually measure a dynamical variable. Moreover, in the laboratory experiment, the values of S1 and S2 are not actually measured but, as they represent the orientation of the spin analyzers SA2 and SA3, are kept fixed for a certain period of time. Unlike in the thought experiment for which Equation (24) was derived, the outcome of an experimental run is not the set of pairs (S1, S2) but rather the number of counts for this particular set of settings. Nevertheless, with some clever manipulations [14, 15], it is possible to express the unit operators that appear in Equation (35) in terms of dynamical variables, the expectations of which can be extracted from the data of single-neutron experiments.
where we used 〈±z|MA|±z〉 = 〈±z|MB|±z〉 = 0 and P(S1, S2|a) is given by Equation (9).
The expressions Equation (37) are remarkable: they show that ϵ2(A) and η2(B) have to be obtained from two incompatible experiments, namely with initial magnetic moments along x and y, respectively. From the point of view of probability theory, this immediately raises the question why, in this particular case, it is possible to derive mathematically meaningful results that involve two different conditional probability distributions with incompatible conditions. As first pointed out by Boole  and generalized by Vorob'ev , this is possible if and only if there exists a “master” probability distribution for the union of all the incompatible conditions. For instance, in two- and three-slit experiments [3, 44–46] such a master probability distribution does not exist by construction of the experiment. Another prominent example is the violation of one or more Bell inequalities which is known to be mathematically equivalent to the statement that a master probability distribution for the relevant combination of experiments does not exist [30, 40, 47, 48]. However, in contrast to these two examples, in the case of the neutron experiment, one can devise a realizable experiment that simultaneously yields all the averages that can be obtained from two experiments (one with a = x and another one with a = y) of the kind shown in Figures 1 or 2.
Our proof is based on the extension of the filtering-type experiment shown in Figure 2 to three dichotomic variables. Imagine that instead of placing detectors in the output beams that emerge from magnets M1 and M2, we place four identical magnets with their magnetic fields in the direction d and count the particles in each of the eight beams. A calculation, similar to the one that lead to Equation (12), yields 
for the probability to observe the given triple (S1, S2, S3). Choosing a = x, c = y, and b = d = x cos ϕ + y sin ϕ it follows immediately that 〈S1〉, 〈S2〉, and 〈S1S2〉 obtained from Equation (12) with a = x agree with the same averages computed from Equation (38). Likewise, 〈S1〉, 〈S2〉, and 〈S1S2〉 obtained from Equation (12) with a = y coincide with 〈S2S3〉, 〈S1S3〉 and 〈S1S2〉 computed from Equation (38).
5. Uncertainty Relations: Event-Based Simulation
In the neutron experiment [14, 15] and therefore also in our event-based simulation, the numerical values of ϵ(A) and η(B) are obtained by counting detection events. Let N(S1, S2|a) denote the count for the case in which the direction of the magnetic moment of the incoming neutrons (after stage 1) is a and the analyzers SA2 and SA3 are along the directions S1 and S2, respectively. Then, we have
As shown in Erhart et al. , Sulyok et al. , the neutron counts observed in the single-neutron experiment yield numerical values of ϵ(A)η(B) + ϵ(A)σ(B) + σ(A)η(B) which are in excellent agreement with the quantum theoretical prediction .
We have already demonstrated that the “classical” event-based simulation model produces results for the averages which, within the statistical errors, cannot be distinguished from those predicted by quantum theory. Therefore, it is to be expected that the data generated by the event-by-event simulation also satisfies the universally valid error-disturbance uncertainty relation Equation (24) and as shown in Figure 5, this is indeed the case. As expected, the data produced by the event-based simulation also violate Equation (34), independent of whether we use pseudo-random numbers (see Equation 18) or the DLM rule (see Equation 19) to model the operation of the spin analyzer.
Figure 5. Uncertainties ϵ(A)η(B) + ϵ(A)σ(B) + σ(A)η(B) (circles) and ϵ(A)η(B) (squares) as obtained from the event-by-event simulations (solid symbols) of the single-neutron experiment [14, 15], the experiment itself ((open symbols), data kindly provided to us by G. Sulyok and Y. Hasagawa) and quantum theory. The solid horizontal line represents the lowerbound in Equation (24). The other solid lines represent the corresponding quantum theoretical prediction as obtained from Equation (35). It is clear that the naive application of the Heisenberg uncertainty relation, ϵ(A)η(B) ≥ 1 is at odds with the prediction of quantum theory, the event-based simulation data and with the experimental data [14, 15]. On the other hand, the experimental data (open symbols, 14, 15) and the results of the event-based simulation (solid symbols) comply with inequality Equation (24). Top: Simulation performed using the probabilistic model for the spin analyzer, see Equation (18). Bottom: Simulation performed using the deterministic learning machine model for the spin analyzer, see Equation (19), with γ = 0.999. In both cases, the number of input events per detuning angle is N = 10000.
Finally, for the sake of completeness, we show that the event-by-event simulation produces data which complies with the standard Heisenberg-Robertson uncertainty relation Δ(σx)Δ(σy) ≥ |〈σz〉|. Without loss of generality, the state of the spin-1/2 particle may be represented by the density matrix Equation (7), also if it is interacting with other degrees of freedom and the inequality Δ(σx)Δ(σy) ≥ |〈σz〉| reads
or, using Equation (7),
The last inequality also trivially follows from the constraint a2x + a2y + a2z ≤ 1. As in the case of Equation (36), there is no ħ in Equation (40), in agreement with the idea that ħ may be eliminated by re-defining the units of mass, time, etc. [41, 42].
The simulation procedure that we use is as follows.
1. Loop over the values (minitial)z = az = −1,… 1 in small steps, e.g., in steps of 0.05.
2. Generate a uniform pseudo-random number 0 < r < 1. Compute and . This step yields a direction of the magnetization minitial which is chosen randomly in the x − y plane.
3. Generate N messengers with message minitial and send them through a spin analyzer aligned along the x-direction. Count the messengers that pass the spin analyzer. Repeat this procedure for spin analyzers aligned along the −x, ± y and ± z-direction. Processing the N messengers yields the counts N(x|a), N(−x|a), etc.
4. Compute the averages 〈σx〉 ≈ (N(x|a) − N(−x|a))/(N(x|a) + N(−x|a)), etc.
5. Go to step 1. as long as az ≤ 1.
6. Plot the results for (1 − 〈σx〉2)(1 − 〈σy〉2) and 〈σz〉2 as a function of az.
The results of the event-based simulation are shown in Figure 6. Within the usual statistical errors, the classical, statistical model produces data which comply with the Heisenberg-Robertson uncertainty relation Equation (40).
Figure 6. Event-based simulation results for (1 − 〈σx〉2)(1 − 〈σy〉2) (blue circles) and 〈σz〉2 (red squares) for different values of − 1 ≤ az ≤ 1. The solid black line represents the quantum theoretical lowerbound a2z. For each value of az, the initial direction of the magnetic moments is (ax, ay, az) where (ax, ay) is a point on the circle with radius chosen using a uniform pseudo-random number. The fluctuations of the data (1 − 〈σx〉2)(1 − 〈σy〉2) reflect the fact that the initial states with different values of az are uncorrelated. For each value of az, N = 10000 messengers were created. The dotted line is a guide to the eyes only.
We have shown that a genuine classical event-based model can produce events such that their statistics satisfies the (generalized) Heisenberg-Robertson uncertainty relation which, according to present teaching, is a manifestation of truly quantum mechanical behavior.
One might be tempted to argue that in the event-based model, the direction of magnetic moment is known exactly and can therefore not be subject to uncertainty. However, this argument is incorrect in that it ignores the fact that the model of the spin analyzers generates [through the use of pseudo-random numbers, see Equation (18) or the update rule Equation (19)] a distribution of outcome frequencies. In fact, as is well-known, the variance of any statistical experiment (including those that are interpreted in terms of quantum theory) satisfies the Cramér-Rao bound, a lower bound on the variance of estimators of a parameter of the probability distribution in terms of the Fisher information . The Cramér-Rao bound contains, as a special case, Robertson's inequality Δ(x)Δ(p) ≥ ħ/2 [16, 19, 23–25, 49]. The observation that a classical statistical model produces data that complies with “quantum theoretical” uncertainty relations is a manifestation of this general mathematical result. The uncertainty relations provide bounds on the statistical uncertainties in the data and, as shown by our event-based simulation of the neutron experiment [14, 15], are not necessarily a signature of quantum physics, conjugate variables, etc.
As mentioned in the introduction, the event-based approach has successfully been applied to a large variety of single-photon and single-neutron experiments that involve interference and entanglement. In the present paper, we have shown that, without any modification, the same simulation approach can also mimic, event-by-event, an experiment that probes “quantum uncertainty.” As none of these demonstrations rely on concepts of quantum theory and as it is unlikely that the success of all these demonstrations is accidental, one may wonder what it is that makes a system genuine “quantum.”
In essence, in our work we adopt Bohr's point of view that “There is no quantum world. There is only an abstract physical description” (reported by Petersen , for a discussion see Plotnitsky ) and that “The physical content of quantum mechanics is exhausted by its power to formulate statistical laws” . Or, to say it differently, quantum theory describes our knowledge of the atomic phenomena rather than the atomic phenomena themselves . In other words, our viewpoint is that quantum theory captures, and does so extremely well, the inferences that we, humans, make on the basis of experimental data . However it does not describe cause-and-effect processes. Quantum theory predicts the probabilities that events occur, but it cannot answer the question “Why are there events?” , very much as Euclidean geometry cannot answer the question “What is a point?.” On a basic level, it is our perceptual and cognitive system that defines, registers and processes events. Events and the rules that create new events are the key elements of the event-based approach. There is no underlying theory that is supposed to give rise to events and everything follows by inference on the basis of the generated data, very much like in real experiments.
The implication of the work presented in our paper is that the beautiful single-neutron experiments [14, 15] can be explained in terms of cause-and-effect processes in an event-by-event manner, without reference to quantum theory and on a level of detail about which quantum theory has nothing to say. Furthermore, our work suggests that the relevance of “quantum theoretical” uncertainty relations to real experiments needs to be reconsidered.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Koen De Raedt, Karl Hess, Thomas Lippert and Seiji Miyashita for many stimulating discussions.
14. Erhart J, Sponar S, Sulyok G, Badurek G, Ozawa M, Hasegawa Y. Experimental demonstration of a universally valid error-distrurbance uncertainty relation in spin measurements. Nat Phys. (2012) 8:185–9. doi: 10.1038/nphys2194
15. Sulyok G, Sponar S, Erhart J, Badurek G, Ozawa M, Hasegawa Y. Violation of Heisenberg's error-disturbance uncertainty relation in neutron-spin measurements. Phys Rev A. (2013) 88:022110. doi: 10.1103/PhysRevA.88.022110
22. Khrennikov AY, Nilsson B, Nordebo S. Classical signal model reproducing quantum probabilities for single and coincidence detections. J Phys Conf. (2011) 361:012030. doi: 10.1088/1742-6596/361/1/012030
31. Kroupa G, Bruckner G, Bolik O, Zawisky M, Hainbuchner M, Badurek G, et al. Basic features of the upgraded S18 neutron interferometer set-up at ILL. Nucl Instrum Methods Phys Res A (2000) 440:604–8. doi: 10.1016/S0168-9002(99)01049-9
33. Leggett AJ. Quantum mechanics at the macroscopic level. In: de Boer J, Dal E, Ulfbeck O, editors. The Lessons of Quantum Theory: Niels Bohr Centenary Symposium. (Amsterdam: Elsevier) 1987. p. 35–58.
53. De Raedt H, Katsnelson MI, Michielsen K. Quantum theory as the most robust description of reproducible experiments. Quantum Phys. (2013). Available online at: http://arxiv.org/abs/1303.4574
Keywords: discrete event simulation, neutron experiments, quantum mechanics, uncertainty relations, foundations of quantum mechanics
Citation: De Raedt H and Michielsen K (2014) Discrete-event simulation of uncertainty in single-neutron experiments. Front. Physics 2:14. doi: 10.3389/fphy.2014.00014
Received: 27 November 2013; Accepted: 14 February 2014;
Published online: 07 March 2014.
Edited by:Sauro Succi, National Research Council of Italy, Italy
Reviewed by:Marc Thilo Figge, Leibniz-Institute for Natural Product Research and Infection Biology-Hans-Knoell-Institute, Germany
Miller Mendoza Jimenez, ETH Zuerich, Switzerland
Copyright © 2014 De Raedt and Michielsen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hans De Raedt, Zernike Institute for Advanced Materials,University of Groningen, Nijenborgh 4, NL-9747AG Groningen, Netherlands e-mail: email@example.com