Research ArticleQuantum Mechanics

In search of multipath interference using large molecules

See allHide authors and affiliations

Science Advances  11 Aug 2017:
Vol. 3, no. 8, e1602478
DOI: 10.1126/sciadv.1602478


The superposition principle is fundamental to the quantum description of both light and matter. Recently, a number of experiments have sought to directly test this principle using coherent light, single photons, and nuclear spin states. We extend these experiments to massive particles for the first time. We compare the interference patterns arising from a beam of large dye molecules diffracting at single, double, and triple slit material masks to place limits on any high-order, or multipath, contributions. We observe an upper bound of less than one particle in a hundred deviating from the expectations of quantum mechanics over a broad range of transverse momenta and de Broglie wavelength.


Many quantum technologies rely on the ability to prepare and manipulate the wave function ψ of a quantum object. For most applications, the advantage of quantum systems over their classical counterparts stems from the superposition principle—quantum states can be coherently split, simultaneously exploring larger regions of phase-space as they evolve than allowed classically. Typically, the information of interest is encoded in the relative phase difference between superpositions of quantum degrees of freedom, such as energy, spin, momentum, or position. In practice, this phase difference is measured by recombining the superposition states to form an interferometer. Examples include atomic clocks (1), quantum processing devices (2, 3), magnetic sensors (4), and matter-wave interferometers using atoms (5, 6) and macromolecules (7).

The interpretation of the wave function continues to intrigue both theorists (811) and experimentalists (12, 13). Although it is not directly observable, it can be inferred from measurements of the probability density P using Born’s ruleEmbedded Image(1)

As long as the underlying quantum dynamics are linear, Eq. 1 can be derived from decision theory (14), entanglement (15), or measurements confined to closed Hilbert spaces (16). However, nonlinear extensions to quantum mechanics for massive objects (1723), which are stimulated by the apparent absence of quantum phenomena in macroscopic systems, can, in principle, result in modifications to Born’s rule (24). This is somewhat akin to the breakdown of the superposition principle in nonlinear optics. Although some nonlinear extensions to quantum mechanics (25, 26) have been ruled out (27) or constrained by experiment (28), recent advances (29) in probing quantum phenomena in systems of increasing mass (3032) or larger state separation (33) motivate further experimental exploration.

Sorkin (34) showed that the quadratic form of Born’s rule implies that no matter how many paths are accessible to a quantum object in configuration space, the probability of a particular outcome can always be decomposed into a series of terms containing at most two of the allowed paths. This insight has guided a number of recent experiments searching for multipath interference using coherent light and single photons (3538), light carrying orbital angular momentum (39), and molecular spins in nuclear magnetic resonance experiments (40). Here, we compare the interference patterns created by free-flying macromolecules behind single, double, and triple slits, providing an explicit test of multipath interference for the center-of-mass wave function of massive quanta.


Figure 1A shows the experimental layout. A vacuum window coated with a thin, homogeneous layer of phthalocyanine (PcH2) provides a source of molecules with a mass M = 515 atomic mass unit (amu). A tightly focused laser beam desorbs them, producing a molecular beam along z with a mean velocity v = 240 m/s and a full width at half maximum Δv = 130 m/s. This corresponds to de Broglie wavelengths in the range of about λdB = 2.5 to 5.0 pm. After propagating a distance L1 = 1.55 m under high vacuum, the molecules impinge on a 25-nm--thick amorphous carbon mask containing a vertical array of one single, two double, and triple slits. The slits have an average transverse opening width of a = 80 nm. The double slits have periodicities d = 100 and 200 nm, and the triple slit has a periodicity of d = 100 nm. Figure 2 shows a transmission electron microscopy (TEM) image of a part of the triple slit (see Materials and Methods and the Supplementary Materials). The slits are arranged vertically above one another on the same substrate and are exposed to the molecular beam simultaneously.

Fig. 1 Experimental setup.

(A) Focused laser source produces a thermal beam of PcH2 molecules, which diffracts at a vertical array of single, double, and triple slits, which are aligned to the local gravitation field g, before landing on a thin quartz detection screen. The deposited molecules are observed using high-resolution fluorescence imaging. (B) Schematic of the triple slit. The openings (black) have a transverse width a = 80 nm, and their centers are separated by a distance d = 100 nm.

Fig. 2 Diffraction mask.

TEM image of a part of the M4 submask of the triple slit. The openings in the mask (white) are a = 80 nm wide and have a period of d = 100 nm.

In this experiment, the transverse coherence length of the molecular beam at the mask 2λdBL1x is limited by the size of the molecular source [Δx = 1.6(1) μm] and ranges from 5 to 10 μm. This exceeds the maximum transverse separation between slits in each mask by more than an order of magnitude. Behind the mask, the molecules propagate a further distance L2 = 0.58 m before they are adsorbed on a quartz plate where they are imaged using laser-induced fluorescence.

Consider the three idealized slits illustrated in Fig. 1B. If the spatial coherence of a molecular beam encountering these apertures extends over all three slits, the final wave function is described by the coherent sum of the individual wave functions transiting each slitEmbedded Image(2)This results in a probability density PABC = |ψABC| = |ψA + ψB + ψC|2, which can be written in terms of one- and two-path interference contributions only asEmbedded Image(3)Here, Embedded Image represents the probability density expected, according to Born’s rule, when all three paths are available. PAB is the measured probability density when only paths A and B are open, PA applies when only path A is open, and so on. By comparing measured three-path interference patterns with a combination of one- and two-path patterns, one can therefore provide an upper bound to the contribution of multipath interference using Sorkin’s criterion (34) through a nonzero measurement of the dimensionless parameterEmbedded Image(4)Here, x is the transverse position on the detection screen and λdB is the de Broglie wavelength of the interfering quanta.

Within an experimental run, we deposit about 2 × 104 molecules in an area of the detection screen measuring 80 × 280 μm2. This number density is sufficiently dilute, even in the densest regions, that Nγ, the number of fluorescence photons detected, remains proportional to the number of molecules in a given region, N = αNγ.

According to Eq. 4, seven different interference patterns are required to place bounds on ε(x, λdB). However, in this experiment, the interference patterns PAB and PBC are expected to differ only by a translation of 100 nm along x. This is two orders of magnitude smaller than the spatial separation of the diffraction fringes and below the limit of our spatial resolution. We can therefore simplify Eq. 4 asEmbedded Image(5)

Here, N is the number of molecules, subscripts label the number of paths available, and superscripts denote the respective slit width and separation. This simplification is valid because Embedded Image, Embedded Image, Embedded Image and Embedded Image. This allows only four submasks to be used: M1 with a single slit, M2 with two slits separated by d, M3 with two slits separated by 2d, and M4 with three slits separated by d.

Figure 3 (A to D) shows typical fluorescence images for molecules diffracted at submasks M1 to M4 in a region of 80 × 280 μm2. This region was chosen to be slightly larger than the single-slit diffraction pattern, which sets the natural length scale for these experiments. In all images, we observe a broad transverse distribution of molecules, consistent with a single-slit diffraction envelope. However, where more than one path is available, additional fringes, consistent with λdB = xd/L2 , are observed. Here, m is an integer labeling the diffraction order, and we have made the small-angle approximation sin(θ) ≈ θ ≈ x/L2. Slower molecules have more time to fall as they travel from the source to the screen, resulting in larger separations between diffraction orders than those of faster particles. The molecular velocity, and therefore the de Broglie wavelength, is encoded in the vertical position of molecules on the detection screen. Hence, we are able to postselect molecules with different de Broglie wavelengths after each experiment. Figure 3E shows the corresponding value of ε(x, λdB) for the interference patterns shown in Fig. 3 (A to D) after subtraction in accordance with Eq. 5. We see that the deviation from zero is typically at the one-molecule level or fewer.

Fig. 3 Molecule diffraction patterns.

Fluorescence images of molecules after diffraction at (A) a single slit, (B) a double slit with period 2d, (C) a double slit with period d, and (D) a triple slit with period d. Slower molecules fall further under gravity in the time it takes them to reach the screen. This results in a larger separation between diffraction orders further down the screen. All images are aligned vertically with respect to gravity. In (D), we have highlighted the center of the m = 0 fringe. (E) Sorkin parameter ε(x, λdB) for the interference patterns (A to D) calculated in accordance with Eq. 5.

To compare results from different experimental runs as well as our values to earlier measurements with photons (35, 37, 40), we rescale ε(x, λdB) by the total number of molecules detected for each de Broglie wavelength, NtotaldB), within the triple-slit interference patternEmbedded Image(6)This ratio compares the number of molecules that deviate from the predictions of quantum mechanics to the total number of molecules detected for a given de Broglie wavelength and therefore provides a measure of any potential violation that is independent of the particle flux and integration time. Note that the precise definition of κ varies throughout the existing literature. In many cases, it is the intensity of particles detected at the center of the interference pattern that is used. However, in far-field molecule interferometry, the interaction energy between the interfering particles and the diffractive element can alter the relative intensity of each diffraction order, resulting in species- and diffraction-dependent values of κ. In contrast, the total number of molecules detected remains independent of the diffraction process at low energies.

Figure 4A shows the normalized Sorkin parameter κ(x, λdB) averaged over five experimental runs, each of which takes about 8 hours to accumulate. We see that any residual structure remaining in ε(x, λdB) from Fig. 3E is greatly reduced by the combined action of averaging and normalizing to the total number of molecules detected. To produce this figure, we have combined the four interference patterns of Fig. 3 (A to D) in accordance with Eqs. 4 to 6. The absence of pronounced peaks qualitatively constrains multipath contributions to the interference patterns generated in our molecule interferometer.

Fig. 4 Bounding multipath interference.

(A) Normalized Sorkin parameter κ(x, λdB). (B) κ(0, λdB), the normalized Sorkin parameter as a function of de Broglie wavelength at the center of the m = 0 interference fringe. (C) κ(x, 3.5 pm), the normalized Sorkin parameter as a function of transverse position on the detection screen for a mean de Broglie wavelength of 3.5 pm. Each figure shows the average over five different experimental runs. In (B) and (C), black lines show the mean, whereas shaded gray areas enclose the 1σ SE. The spatial resolution in all figures is limited by the point-spread function of the imaging system, which is 1.6 μm, approximately four pixels in the plane of the molecular detection screen.

In a single measurement, we obtain interference patterns for particles with de Broglie wavelengths spanning more than an entire octave as well as a broad range of transverse positions. More details relating to the conversion of vertical position on the charge-coupled device (CCD) camera to de Broglie wavelength can be found in Materials and Methods and the Supplementary Materials. To quantitatively constrain multipath interference, we focus on κ recorded at the center of the zeroth diffraction fringe and a transverse slice with a de Broglie wavelength of 3.5 pm. Figure 4 (B and C) shows κ(0, λdB) and κ(x, 3.5 pm) averaged over five experimental runs. The black line depicts the mean value, with the shaded region illustrating the accompanying 1σ SE. We see no significant deviation of κ from zero beyond our statistical uncertainty across the entire range of de Broglie wavelengths or transverse momenta probed.


Our molecular experiment provides a direct constraint on contributions to multipath interference using massive particles by comparing the interference patterns resulting from single-, double-, and triple-slit experiments. Matter waves have a different dispersion relation than light, which has enabled us to place bounds on multipath interference across a range of de Broglie wavelengths spanning more than an octave, in addition to a broad range of transverse momenta, in a single experimental run.

We observe that, on average, less than one molecule in a hundred deviates from the coherent, Schrödinger description of quantum mechanics when compared with the corresponding local number of molecules of the same de Broglie wavelength in a triple-slit experiment. This results in an upper limit of |κ|2 ≤ 10−2 across de Broglie wavelengths ranging from 2.5 to 5.0 pm and transverse positions ranging over 80 μm centered on the zeroth interference fringe.

Unlike photons, which are detected destructively, here, the interfering quanta are durable molecules that stick to a screen where they can be detected in a nondestructive manner, even days after deposition (41). Whether the wavelike description of quantum objects extends to higher masses (29, 42) remains an open question that motivates further experiments using larger, more massive, and complex particles. However, the long integration time required to produce a single-molecule interference pattern limits our statistics, which are the dominant uncertainty in this experiment. Further progress hinges on the development of new sources for slow, cold, and collimated macromolecular, cluster, or nanoparticle beams.

Direct tests of multipath interference using beams of molecules are best implemented in a single diffraction grating scheme, such as the one used in our present work. State-of-the-art material diffraction gratings have already reached their natural size limit, with masks as thin as a single atomic layer (43); further advances may therefore require optical gratings or new methods of producing effective slits for matter waves.

High-order interference terms, resulting in finite values of κ, are predicted by quantum mechanics. These additional terms can be understood using the Feynman path integral approach, wherein nonlinear paths also add up constructively to contribute to the final interference pattern (44, 45). For a triple-slit experiment, such as the one presented here, the maximum value of this multipath contribution scales as Embedded Image. The relatively high mass of our interfering quanta, which are derived from a thermal source, results in a short de Broglie wavelength. This suppresses the multipath contribution to the interference patterns well below our current precision. However, matter-wave interferometry using cold, low-mass atoms, such as helium or lithium, can achieve very long de Broglie wavelengths and may be able to reach the necessary precision.


Experimental design

The experiment is housed in a 2.14-m-long vacuum chamber with a background pressure below 10−7 mbar (41). At this pressure, collisions with background gas can be safely neglected. Furthermore, collisions predominantly prevent the molecules from reaching the small detection region, whereby they reduce the count rate without affecting the shape of the diffraction patterns. Each experimental run took about 8 hours. Approximately 50 mg of PcH2 (M = 515 amu) was evaporated onto a CF100 vacuum window, resulting in a homogeneous layer, only a few micrometers thick (41). Evaporating molecules from the coated window using a tightly focused 60-mW laser beam with a wavelength of 421 nm produced a dilute molecular beam along z, with, on average, less than one molecule passing through the grating at any given time. Particles leave the surface from a circular spot with a diameter of Δx = 1.6(1) μm, measured using an optical microscope. The beam was replenished in situ by translating the window laterally to expose new regions of the surface to the laser beam. This ensures that the position of the beam relative to the diffraction mask remains fixed. After desorption, the molecules propagate a distance of L1 = 1.55 m before coherently illuminating the diffraction mask consisting of four sets of slits (M1 to M4). The angular spread of the molecular beam is well described by a cosine distribution, resulting in a flux of particles varying by <10−6 across the 0.6-mrad solid angle required to simultaneously illuminate all slits. The mask was positioned relative to the molecular beam with a precision of 25 μm using an xy translation stage, and the slits were orientated with respect to gravity (y axis) to better than a few milliradians.

After passing through, the mask molecules arrive at a quartz plate where they stick. The molecular diffraction patterns were then revealed by illuminating the quartz plate using a 661-nm laser and observing the resulting fluorescence from the deposited molecules. The illumination laser beam has a top-hat intensity distribution produced by reflecting a Gaussian beam from the surface of a spatial light modulator (SLM) before passing through a rotating ground-glass diffuser. Each interference pattern was imaged using a 20-s exposure, during which the SLM cycled through 100 different interference masks, averaging over the laser beam phase fronts, resulting in an average intensity that varies by a few percent over the entire region of interest. Fluorescence photons were collected using a 20× objective (Zeiss Ultrafluar), which focuses them onto an electron-multiplying CCD camera (Andor iXon 885-KCSVP) containing 1003 pixels × 1004 pixels with an area of 8 × 8 μm2 each. The point spread function of the imaging system is about four times larger than the size of a single pixel and limits our spatial resolution. We therefore averaged the two-dimensional (2D) interferograms over blocks of 4 pixels × 4 pixels.

Fluorescence photons were separated from the illumination light using a 711-nm band-pass filter with a bandwidth of 25 nm and an optical density of >7. Special care is taken to avoid saturation and nonlinear effects by having only thin layers of molecules deposited. On average, each molecule contributes 1800 CCD counts (more details are described in the Supplementary Materials).

Diffraction mask

The mask was fabricated into a 25-nm-thick amorphous carbon membrane (Ted Pella). The slits were machined using a focused gallium ion beam (Raith ionLine) with a kinetic energy of 35 kV, a one-dimensional (1D) beam density of 48 × 103 pC/m, a current of 19.2 pA, and a dwell time of 0.72 ms. During the writing process, the ion beam was translated across the substrate in steps of 1 nm. The focused waist of the ion beam is about 60 nm, requiring each individual subslit to be exposed twice to realize a slit width of 80 nm.

The morphology of the mask was investigated using TEM. The mask extends over a height (y axis) of 958 μm and consists of four submasks: M1 containing a single slit, M2 and M3 containing double slits, and M4 containing a triple slit. Each slit is 21 μm long and consists of a vertical array of 19 subslits, separated from one another by 92(5)-μm-long -support bars, which stabilize the mask and prevent deformations of the ultrathin membrane. The submasks were separated from one another by 300 μm to ensure that the diffraction patterns from neighboring slits do not overlap on the detection screen. The average width of the subslits within a given slit varies by less than 1 nm for all slits (SD from the 19 subslits). For submask M1, the opening width is 86(5) nm. For submask M2, the opening widths are 75(5) and 81(5) nm with a separation of 102(5) nm. For submask M3, the opening widths are 79(5) nm and a = 82(5) nm, and their separation is 2d = 195(5) nm. For submask M3, the opening widths are 78(5), 83(5), and 81(5) nm, with separations of 97(5) and 98(5) nm. Here, error bars are the systematic uncertainty of the electron microscope and are consistent with the statistical variations we observed between slits. The 230:1 aspect ratio of the slits and their large opening fraction, a/d = 0.8, helped increase the molecular transmission and therefore signal to noise.

Analysis of diffraction patterns

All four submasks were exposed simultaneously to the molecular beam. However, the limited field of view of the imaging system requires each interference pattern to be imaged separately. To determine κ for a given experimental run, we first centered the individual images Embedded Image with respect to each other. The horizontal and vertical axes of each image were arranged according to the center of the m = 0 interference fringe, corresponding to x = 0, and the median velocity, respectively. We can align the transverse position of each interference pattern with a precision better than 1 μm in the plane of the detection screen and a vertical precision corresponding to a de Broglie wavelength of 5 fm. Further details are shown in the Supplementary Materials. The molecular velocity for a given vertical pixel was determined from the separation, δx, between the m = ± 1 diffraction orders of the triple-slit diffraction pattern according to v = 4πℏL2/Mdδx.

The dark counts per pixel N0 were subtracted for each interferogram individually. Within a given experimental run, the fractional change in N0 between different interference patterns is less than 5 × 10−3. Dark counts were determined by measuring the mean number of counts detected in a square region of 50 pixels × 50 pixels where no molecules are expected. The size and region used to estimate N0 affect the measurement of κ(x, λdB) by less than one part in 102 when the region of interest is varied between 20 and 200 pixels. This is two orders of magnitude smaller than the statistical uncertainty. The precise position of the same region of interest relative to the center of the interference pattern has a similarly small effect on the measurement of κ(x, λdB) when the distance is varied between about 700 and 900 pixels.

Surface interactions

Deviations in the slit morphology alter the position-dependent Casimir-Polder interaction between individual molecules and the mask. This has been studied in detail for matter-wave diffraction at material gratings, including the effect of thermally populated rotational and vibrational states on far-field diffraction (46, 47). For each slit, the interference pattern recorded is the ensemble average over the entire phase-space distribution of the molecular beam and the mask. The modification to the interference patterns are therefore negligible for our slits, which are extremely uniform—varying by less than 5 nm from slit to slit. In a previous work by Sinha et al. (35), the seven different interference patterns necessary to measure ε(x, λdB) were achieved by blocking slits using an additional mask. In our apparatus, this is not possible because of the extremely small opening apertures necessary for our matter waves, which are typically 1000 times smaller than their optical counterparts. In our interferometer, the additional blocking mask would need to be in the vicinity of 100 nm of the diffraction mask, where the strong Casimir interaction between the two masks would modify the transmission function of the unblocked beams. At these short distances, fabrication imperfections in the additional blocking mask would also contribute a systematic distortion of the interference patterns, with a similar magnitude to deviations arising from the diffraction masks themselves.


Supplementary material for this article is available at

Diffraction mask

Molecule detection

Interference patterns

fig. S1. TEM images of the four submasks.

fig. S2. Intensity distribution of the 661-nm illumination light.

fig. S3. Surface migration of PcH2 molecules behind a grating with a period of 2 μm.

fig. S4. Molecule number as a function of detector position.

fig. S5. Raw data before dark count subtraction.

fig. S6. Determining molecular velocities from fringe spacing.

fig. S7. Reconstructing molecular velocity distributions.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank B. Dakic and C. Brukner for helpful discussions and T. Susi for help with the recording of TEM images. Funding: We are grateful for financial support from the Austrian Science Fund (FWF) (W1210-3), the European Research Council (320694), and the European Commission (304886). J.P.C. is supported by a Vienna Center for Quantum Science and Technology fellowship. C.B. acknowledges support from the Alexander von Humboldt-Foundation through a Feodor-Lynen fellowship. Author contributions: J.P.C. and M.A. devised the experiment. Y.L. and O.C. prepared the diffraction mask. C.B. and C.K. acquired the data. J.P.C. analyzed the data. All authors discussed the results and contributed to the writing of the paper. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from J.P.C. (joseph.cotter{at}

Stay Connected to Science Advances

Navigate This Article