A Signal Processing Model of Quantum Mechanics
Abstract
This paper develops a deterministic model of quantum mechanics as an accumulationandthreshold process. The model arises from an analogy with signal processing in wireless communications. Complex wavefunctions are interpreted as expressing the amplitude and phase information of a modulated carrier wave. Particle transmission events are modeled as the outcome of a process of signal accumulation that occurs in an extra (nonspacetime) dimension.
Besides giving a natural interpretation of the wavefunction and the Born rule, the model accommodates the collapse of the wave packet and other quantum paradoxes such as EPR and the AhanorovBohm effect. The model also gives a new perspective on the ‘relational’ nature of quantum mechanics: that is, whether the wave function of a physical system is “real” or simply reflects the observer’s partial knowledge of the system. We simulate the model for a 2slit experiment, and indicate possible deviations of the model’s predictions from conventional quantum mechanics. We also indicate how the theory may be extended to a field theory.
Keywords:
quantum mechanics Born rule signal processing threshold process quantum paradoxes∎
1 Wavefunction analogy in wireless communications
Several physical systems are characterized by a process of accumulation (of energy, charge, etc.), which leads to an activation event once the accumulation attains a certain threshold. Examples of such systems include lightning and nerve impulse transmission. In many cases the accumulation process is described in terms of a continuous field, while attaining the threshold triggers a discrete event. This simultaneous presence of discrete and continuous aspects is reminiscent of quantum mechanics.
Signal acquisition in wireless digital communications also follows this same general pattern. Consider a mobile receiver moving randomly within a region in which a modulated carrier wave is broadcast. The carrier wave is modulated both in amplitude and phase. In order to detect the broadcasted signal, the receiver accumulates its received signal until a detection threshold is reached. We shall construct a mathematical model of a system, in which that the location where detection occurs obeys a probability distribution reminiscent of the quantum wavefunction. In our model, the wireless signal has the following characteristics:

The carrier frequency is , so that the signal has the general mathematical form . Such a signal is commonly represented by its “complex amplitude” .

The transmitted signal (at the transmitter) has constant complex amplitude over time intervals of length , where ( is called the “chip width” in digital communications ((1))). The probability distribution of complex amplitudes is Gaussian, so that real and imaginary parts are independent, identically distributed (i.i.d) standard normal random variables with mean 0 and variance 1.

The ratio of field amplitude to transmitted signal amplitude (denoted by depends on the field location , but is independent of time. For mathematical simplicity, we assume that the ratio assumes one of a finite set of complex values ; and that within the (finite) region of interest, the sets all have equal area (see Figure 1).
The receiver has the following characteristics:

The receiver consists of an oscillating circuit with natural frequency , which is driven by the signal field at the receiver’s current location.

The receiver moves slowly enough so that its field amplitude does not change significantly over time intervals of length , where is an integer .

The receiver moves in such a way that its position uniformly samples the entire region of interest (for instance, by random walk).

Our mathematical proof (see Appendix) requires that the receiver’s fields over the time intervals and are statistically independent whenever . Strictly speaking, a receiver moving under random walk will not satisfy this condition: instead, the receiver would have to make uniformlydistributed random jumps at times A rigorous treatment with randomwalk motion would require a more careful analysis.

The receiver detects the signal when the receiver’s amplitude exceeds a fixed threshold .
Note this simple model does not include any effects from polarization, propagation delay or Doppler phase shifting. With the above assumptions, the field at receiver position r can be expressed (with the aid of complex amplitudes) as follows:
(1) 
where

= (field amplitude at ) / (signal amplitude at transmitter)

may be written as , where are i.i.d. standard normal random variables.

denotes the “ceiling” function, i.e. next largest integer greater than .
We now suppose that the trajectory of the receiver is given by the function . It follows that the equation for the amplitude of the driven oscillator is:
(2) 
This equation may be expressed as the real part of the complex equation
(3) 
The solution of (3) which satisfies is
(4) 
According to our assumptions, the factor in the second integrand oscillates rapidly compared to the rest of the integrand, which causes the second integral to be negligible compared to the first. The model assumptions imply that can be treated as constant over time intervals of length . Using the notation , we have:
(5) 
The oscillation at time has complex amplitude , where
(6) 
According to the model assumptions, each is one of the values . Define the random variable to be the index corresponding to random variable : that is
(7) 
Define as the time index at which first passes a given threshold :
(8) 
Our goal is to evaluate the probability distribution of corresponding to the first passing of the threshold:
(9) 
This corresponds to the probability distribution of the value of field at the location of detection. It is a wellknown fact in signal processing that the rate of accumulation of a random signal is proportional to the signal power, which is in turn proportional to the squared signal amplitude ((1)). It stands to reason that given a signal that assumes different power levels at different times but with equal probabilities, the chance of the accumulated signal passing a fixed threshold while at a certain power level should be proportional to that power level. This is in fact the case; in the Appendix we prove that
(10) 
which are exactly the Born probabilities for the spatial wavefunction .
2 Single quantum detection event model
In this section, we present a model (based on the model in the previous section) that explains quantum detection probabilities. The model is discretized for conceptual clarity and computational tractability; it is fairly straightforward to see how the model could be taken to a continuous limit.
We emphasize that the probability distribution in the previous model arose from the outcome of a process. The process involved sampling the entire region of potential detection before the actual detection was made. This representative sampling was necessary in order for the field strengths to translate into relative probabilities. We want similar characteristics for the quantum process.
In our previous model the process variable is time; this was appropriate because we were only concerned about the spatial position of the receiver at the moment of detection. However, in quantum mechanics, we are concerned about the location of detection events within spacetime. It is impossible to have a process that unfolds in time that at the same time samples all spacetime locations before determining the detection location. For this reason, it is necessary to introduce a new process variable, so that the process of signal accumulation takes place in a nonobservable dimension which we will call the dimension.
Our wireless communications model had a physical receiver which moves within the state space of possible detection locations. Quantum detection (say of a particle on a screen) does not appear to have any corresponding receiver. We therefore introduce the notion of a detectron, which plays the same role as the receiver in our previous model. The accumulation takes place as the detectron jumps around and uniformly samples the set of all potential detection locations. This “jumping around” takes place in the dimension; for fixed , the detectron’s spacetime location is fixed. We emphasize that the detectron is a mathematical construct, and should not be considered as a physical particle; we will say more about the physical nature of detectrons in Section 5.
We also postulate a carrier wave that oscillates as a function of (not as a function of time) having the mathematical form . The frequency is unknown, and does not correspond to any measurable quantity in spacetime. The signal has the following characteristics:

The signal has constant complex amplitude over intervals of length , where . The distribution of complex amplitudes is meanzero Gaussian, with i.i.d. standard normal real and imaginary parts;

The signal is multiplied by a complex field amplitude which is independent of . For mathematical simplicity, we assume that the amplitude assumes one of a finite set of complex values ; and that within the spacetime confines of the detector, the sets all have equal 4volume.
The detectron has the following characteristics:

Associated with the detectron is an oscillator (which varies sinusoidally with ) with natural frequency , which is driven by the signal field at the detectron’s current spacetime location;

The detectron moves in spacetime (as a function of ) slowly enough so that its field amplitude does not change significantly over aintervals of length , where is an integer ;

The detectron moves in such a way that it uniformly samples the spacetime extent of the detector.

The detectron becomes a detection when its oscillator’s amplitude exceeds a fixed threshold .
We can apply this model to the twoslit setup shown in Figure 2. The detectron moves within the spacetime confines of the detection screen. The complex field amplitude corresponds to the conventional Schrödinger wavefunction at the screen, which in the ray approximation is given by:
(11) 
where , are the (observable) wave number and frequency, and , .
We simulated this system using MATLAB, with the parameters shown in Table 1 . We only considered a single time slice, and restricted to the portion of the screen. The locations were discretized into 100 bins; the detectron jumped uniformly randomly from bin to bin every iteration steps. At each iteration, the signal was incremented by , where are i.i.d. complex random variables with standard normal real and imaginary parts. Each time the detection threshold was reached, a detection was logged and the simulation was restarted. Altogether 100,000 detections were logged. Figure 3 shows the detection probability distribution obtained in the simulation. Agreement is very close with the theoretical result , with given by (11). Note that in the simulation, distribution peaks are slightly lower than theoretical values. If quantum probabilities are indeed the result of such an accumulation process, it is possible that measured probability values may be lower than the conventional quantum prediction. Unfortunately, it is not possible to predict the extent of the lowering from our model, because it depends on details of the accumulation process that are not accessible to measurement.
Parameter symbol  Parameter signifigance (distances in wavelengths)  Value 

1/2 the distance between slits  5  
Distance from source to slit screen  
Distance from slit screen to detection screen  
Number from slit screen to detection screen  400  
Number of detections  
Screen halfwidth  
Number of bins (discretization)  100  
Detection threshold  500 
3 Wavefunction formation via accumulation
The preceding section describes how to obtain quantumlike detection probabilities given that a certain “broadcast field” is present. However, it provides no mechanism for the creation of the broadcast field itself. In this section, we show how the broadcast field can be modeled as the result of a process of accumulation that parallels the signal accumulation described above. The wellknown pathintegral expression for the propagator is given by ((2)):
(12) 
Here is the action, and the notation denotes an equallyweighted summation over all possible paths from to . This integral may be seen as the outcome of an accumulation process. We may envision a succession of carrierwave blips, where each blip corresponds to a single path and makes a differential contribution to the field which is proportional to (as shown in Figure 4). Recall that corresponds to the amplitude and phase of a modulated carrier wave; thus expresses the influence of a source at on the amplitude and phase of the wave at when a blip passes between them.
These blips may be associated with detectrons as follows. If is a possible event detection location, then a path that passes through can be identified with a detectron location of . We postulate that simultaneously with causing an incremental change in the field , the blip also increments the overall complex detection signal amplitude by , as described in the communication model in the previous section. In this way, the blips perform a dual mathematical function: they both build up the field, and furnish the uniform random sampling of possible detection sites that is required to obtain Bornrule probabilities (as was shown in the previous section). The process of detection signal buildup and detection is shown in Figure 5.
To see how this works in practice, we focus specifically on the twoslit experiment shown in Figure 2 in two space dimensions, using nonrelativistic electrons of fixed energy as particles. We consider the spatial distribution of detections at time ; due to invariance, this distribution will be independent of . The detection screen corresponds to spatial locations . We assume the source is configured so that , where and denote the spatial locations of the two slits ( and , respectively). We also assume that is large enough so that the signal accumulation process has negligible effect on the size of . We replace with negative (since only negative times contribute to detection at ) and obtain an expression the Schrödinger kernel for paths that exit through the upper slit :
(13) 
where is a constant of proportionality. Similarly, the kernel for paths that exit through the lower slit is . The theoretical expression for the wavefunction is
(14) 
which evaluates to
(15) 
where
(16) 
It was computationally intractable for us to simulate the path integrals that give rise to the Schröedinger kernels themselves. Instead, we assumed the kernels and simulated the cumulative effect of different source points. We use the following algorithm to accumulate both the fields at different screen locations and the overall detection signal:
Initialize: for all detection locations
While
If divides
Change current detectron location ;
Choose random time (uniformly distributed);
End If
;
;
Signal=Signal + ;
End While
Record detection at location
Parameter symbol  Description  Value (in mks units) 

1/2 the distance between slits  m  
Distance from slits to detection screen  m  
Number of iterations between detectron jumps  25  
Number of detections  3600  
Minimum value for random time  m sec  
Maximum value for random time  sec  
Electron frequency  Hz  
Screen halfwidth  m  
Number of bins (discretization)  43  
Detection threshold 
Table 2 shows the parameters of the simulation. Figure 6 shows the results for 3600 detections. The figure shows counts from 1/2 of the detection screen (), which was discretized into 43 bins (with the edge of bin 43 at the center of the screen). The simulation relative counts per bin and the accumulated field are plotted, as well as theoretical detection probabilities from (15). In general, the theoretical curve lies within error bars of the simulation counts; the simulation counts are consistently slightly higher than theoretical probabilities near the interference pattern nulls.
4 Multiple quantum events
In the preceding sections we have described how accumulation of the signal associated with the moving detectron eventually leads to detection when the accumulated signal passes a threshold. This description treats a single quantum event in isolation without considering interactions between quantum events. Accordingly, we now postulate that every spacetime event corresponds to a detectron, and all detectron signals at a given instant are multiplied before accumulation. The situation thus remains as shown in Figure 5, except that the jumps in detection signal magnitude are caused by products of detectron signals. In fact, the effects of this multiplication are already included in the model in Section 2 via the random Gaussian factors which multiply the field at the detectron: these fluctuations in the signal amplitude correspond to the variation of other detectron signals.
Along with the change in detection signal, at each instant the complex amplitude of the field at each detectron location is incremented as follows. Let be the (timeordered) sequence detectron spacetime locations for the preuniverse at . Then the complex wavefunction amplitude at is incremented by , where is a path from to . Given that the paths and are statistically independent for , this gives rise to the expression (12) for the propagator.
In summary, each preuniverse shown in Figure 5 gives a single contribution to the overall detection signal, that accumulates as the preuniverses unfold in the dimension. Once the detection signal attains a threshold for a particular value of , the spacetime universe that we experience is actualized as a “snapshot” at that instant.
5 Extension to a field theory
In the previous sections, we have presented detectrons as moving placeholders for possible detection sites. This picture presumes a division of the universe between fields and detectors. This division of course is not realistic, for the detectors themselves are represented by fields in their own right. We may eliminate this dichotomy by identifying detectrons as field configurations (including the transmitted particle field and the field representing the atoms of the detection screen) that produce contributions to the detection signal. The transition to a field theory may then be accomplished by replacing paths with field configurations. In other words, the preuniverse at each instant has a particular field configuration that produces an overall contribution to the detection signal. As increases, the field configurations vary and their contributions to the detection signal are accumulated. At the instant where the threshold is attained, we arrive at the field configuration corresponding to the observable universe.
6 Explanations of quantum paradoxes
Collapse of the wave packet: In our model, the wave function is seen as an “actual” field that develops via a process of accumulation. The field is not merely a representation of the observer’s partial knowledge; but on the other hand, the field is not directly observable via physical events in spacetime. The apparent “collapse” of the wave packet is due to the fact that the physical universe is only a single “slice” of the entire process.
EPR and Bell’s inequality: Nonlocal effects pose no problem for this model, for the model itself is inherently nonlocal. An EPR experiment where the two spin detectors are aligned parallel will always detect antialigned particles because both detections correspond to the same blip and are thus perfectly correlated. If the detectors are not parallel, the detection probabilities are still determined according to the quantum expression appropriate for that configuration.
AharanovBohm Effect: The AharanovBohm effect shows that fields that are localized in a region where a particle can never be detected can still have an effect on the motion of the particle. This poses no difficulty to our model, because in our model the socalled particle is not an object that travels through spacetime but rather a correlated series of detection events.
Identical particle statistics: More work is needed to introduce spin into the model. However, in light of our framework it is not surprising that “particles” obey special statistics, because “particles” are not separate objects at all. What we call a “particle” is simply a series of correlated events.
7 Comparison With Other Interpretations of Quantum Mechanics
We briefly compare our interpretation with other alternative interpretations of quantum mechanics.
Everett’s “Manyworlds” interpretation(3) requires exponential plethorization of spacetimes. Our model, which embeds spacetime within one additional dimension, possesses a much simpler state space.
Bohm’s quantum mechanics(4) posits that particles such as electrons are able to track along with pilot waves. This appears to imply that these particles have some sort of inner structure. In our model particles are not “objects” at all, so no such complications appear.
Cramer’s transactional quantum mechanics(5) interprets as a wave traveling backwards in time, but gives no explanation why should be interpreted as a probability. Furthermore, transactional quantum mechanics is not very clear about the order in which “transactions” are determined. In our model, all transactions are determined “simultaneously” (at , and a single accumulation process is used to determine all interaction events.
We also remark that none of these alternative models explains why the wavefunction is complex, nor why the squared amplitude is interpreted as a probability.
8 Possible experimental verification
If true, then our model indicates that the usual formula for a quantum wavefunction is a statistical approximation, and small deviations from the probabilities predicted by the wave quation should be expected. In particular, in our simulations we consistently found that detection rates near theoretical wavefunction nulls were higher that the conventional quantum predictions. Unfortunately, the size of these effects would depend on aspects of the process that cannot be directly measured.
9 Conclusions
Our model presents a radically different picture of reality. Traveling “particles” are replaced with series of detection events; as a visual analogy, imagine a series of fireflies in a line that flash successively, giving the impression that a single firefly is moving along the line. Our model replaces temporal causality with atemporal causality; past, present, future are actualized together as the result of a process that occurs in a different dimension. Apparent temporal causality is due to correlation and not causation. The wavefunction is given a physical interpretation as a dynamical field; and the Born rule based on the wavefunction is derived as the natural result of a thresholding process involving this field. The model includes possible differences from conventional quantum mechanics. A lowering of peak probabilities in quantum interference patterns compared to the conventional quantummechanical prediction is a possible result of the model.
10 Appendix: Mathematical derivation of the Born rule
In this section, we prove the Born probability rule,
(17) 
for the wireless communication scenario described in Section 1. We will use the notation and definitions of that section.
In order to investigate the dependence of on , for each fixed we will investigate the event
(18) 
conditioned on fixed sequences of initial s, corresponding to the events
(19) 
We shall show that
(20) 
where is independent of . The events for fixed partition the sample space and . Furthermore, the events partition the event , so we obtain
(21)  
which is the desired result.
We prove (20) as follows. Conditioned on event , we have
(22) 
were the have i.i.d. standard normal real and imaginary parts (we write this as: ). It follows that is a random walk in the complex plane with independent (but not identically distributed) steps. We also have
(23) 
We shall assume that , so is much greater than any individual term in . It follows that the distribution of for all sample paths can be approximated as a standard Brownian motion , where the time variable is given by
(24) 
The sample paths comprised in the event correspond in the Brownian motion picture to sample paths for which for all . For these sample paths, the distribution of corresponds to the position probability density for a standard Brownian motion with absorbing barrier at .
Now there is a close connection between Brownian motion and the heat equation as follows. Let be the probability density at time of a Brownian motion with absorbing barrier at . Then can be found by solving the heat equation with corresponding boundary and initial conditions, which in this case are:

Boundary conditions: for ;

Initial conditions: ,
where (…) is the Dirac delta function. We do not need the complete solution for (which can be expressed in terms of the Bessel functions , but we will make use of the following properties:
 a)

is radial, so we may write as
 b)

is for and ;
 c)

for all ;
These properties can be mathematically proven, but are also intuitive consequences of the physical interpretation of as an evolving temperature distribution within a disk where the boundary is held at zero temperature. In light of property c), at time can be approximated near the boundary as
(25) 
It follows from our identification of with that
(26) 
Since ,it follows that the terms for . By rotating in the complex plane we have that
(27)  
where .
In the case where and , we have
(28) 
Thus the condition can be replaced to a very close approximation by the condition (see also Figure 7) and
(29)  
Note that for since . Now for sample paths with the event is the complement of the event . It follows in summary that
(30) 
Note that unless ; and since , the approximation (26) holds on this range. Our integral becomes
(31) 
Changing variable to , we have
(32)  
Acknowledgements.
Thanks to Walter Wilcox for many helpful suggestions.Footnotes
 email: thron@ct.tamus.edu
 email: thron@ct.tamus.edu
 The actual distribution arising from a product of random signals will be lognormal rather than normal: we used normal random variables for computational simplicity. The results are not sensitive to the particular shape of the distribution.
References
 J. Proakis, Digital Communications Edition, McGraw Hill (2000).
 H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets Edition, World Scientific, Singapore (2004).
 B. S. DeWitt, R. N Graham, eds, The ManyWorlds Interpretation of Quantum Mechanics, Princeton Series in Physics, Princeton University Press (1973).
 D. Bohm and B. J. Hiley, The Undivided Universe: An Ontological Interpretation of Quantum Theory, Routledge & Kegan Paul, London (1993).
 J.G. Cramer, “An Overview of the Transactional Interpretation of Quantum Mechanics”, International Journal of Theoretical Physics 27, 227 (1988).