Beam splitter

A beam splitter or beamsplitter is an optical device that splits a beam of light into a transmitted and a reflected beam. It is a crucial part of many optical experimental and measurement systems, such as interferometers, also finding widespread application in fibre optic telecommunications.

Schematic illustration of a beam splitter cube.
1 - Incident light
2 - 50% transmitted light
3 - 50% reflected light
In practice, the reflective layer absorbs some light.

Beam splitters

Designs

In its most common form, a cube, a beam splitter is made from two triangular glass prisms which are glued together at their base using polyester, epoxy, or urethane-based adhesives. (Before these synthetic resins, natural ones were used, e.g. Canada balsam.) The thickness of the resin layer is adjusted such that (for a certain wavelength) half of the light incident through one "port" (i.e., face of the cube) is reflected and the other half is transmitted due to FTIR (Frustrated Total Internal Reflection). Polarizing beam splitters, such as the Wollaston prism, use birefringent materials to split light into two beams of orthogonal polarization states.

 

Aluminium-coated beam splitter.

Another design is the use of a half-silvered mirror. This is composed of an optical substrate, which is often a sheet of glass or plastic, with a partially transparent thin coating of metal. The thin coating can be aluminium deposited from aluminium vapor using a physical vapor deposition method. The thickness of the deposit is controlled so that part (typically half) of the light, which is incident at a 45-degree angle and not absorbed by the coating or substrate material, is transmitted and the remainder is reflected. A very thin half-silvered mirror used in photography is often called a pellicle mirror. To reduce loss of light due to absorption by the reflective coating, so-called "Swiss-cheese" beam-splitter mirrors have been used. Originally, these were sheets of highly polished metal perforated with holes to obtain the desired ratio of reflection to transmission. Later, metal was sputtered onto glass so as to form a discontinuous coating, or small areas of a continuous coating were removed by chemical or mechanical action to produce a very literally "half-silvered" surface.

Instead of a metallic coating, a dichroic optical coating may be used. Depending on its characteristics, the ratio of reflection to transmission will vary as a function of the wavelength of the incident light. Dichroic mirrors are used in some ellipsoidal reflector spotlights to split off unwanted infrared (heat) radiation, and as output couplers in laser construction.

A third version of the beam splitter is a dichroic mirrored prism assembly which uses dichroic optical coatings to divide an incoming light beam into a number of spectrally distinct output beams. Such a device was used in three-pickup-tube color television cameras and the three-strip Technicolor movie camera. It is currently used in modern three-CCD cameras. An optically similar system is used in reverse as a beam-combiner in three-LCD projectors, in which light from three separate monochrome LCD displays is combined into a single full-color image for projection.

Beam splitters with single-mode^{[clarification needed]} fiber for PON networks use the single-mode behavior to split the beam.^{[citation needed]} The splitter is done by physically splicing two fibers "together" as an X.

Arrangements of mirrors or prisms used as camera attachments to photograph stereoscopic image pairs with one lens and one exposure are sometimes called "beam splitters", but that is a misnomer, as they are effectively a pair of periscopes redirecting rays of light which are already non-coincident. In some very uncommon attachments for stereoscopic photography, mirrors or prism blocks similar to beam splitters perform the opposite function, superimposing views of the subject from two different perspectives through color filters to allow the direct production of an anaglyph 3D image, or through rapidly alternating shutters to record sequential field 3D video.

Phase shift

 

Phase shift through a beam splitter with a dielectric coating.

Beam splitters are sometimes used to recombine beams of light, as in a Mach–Zehnder interferometer. In this case there are two incoming beams, and potentially two outgoing beams. But the amplitudes of the two outgoing beams are the sums of the (complex) amplitudes calculated from each of the incoming beams, and it may result that one of the two outgoing beams has amplitude zero. In order for energy to be conserved (see next section), there must be a phase shift in at least one of the outgoing beams. For example (see red arrows in picture on the right), if a polarized light wave in air hits a dielectric surface such as glass, and the electric field of the light wave is in the plane of the surface, then the reflected wave will have a phase shift of π, while the transmitted wave will not have a phase shift; the blue arrow does not pick up a phase-shift, because it is reflected from a medium with a lower refractive index. The behavior is dictated by the Fresnel equations.^[1] This does not apply to partial reflection by conductive (metallic) coatings, where other phase shifts occur in all paths (reflected and transmitted). In any case, the details of the phase shifts depend on the type and geometry of the beam splitter.

Classical lossless beam splitter

For beam splitters with two incoming beams, using a classical, lossless beam splitter with electric fields E_a and E_b each incident at one of the inputs, the two output fields E_c and E_d are linearly related to the inputs through

${\displaystyle \mathbf {E} _{\text{out}}={\begin{bmatrix}E_{c}\\E_{d}\end{bmatrix}}={\begin{bmatrix}r_{ac}&t_{bc}\\t_{ad}&r_{bd}\end{bmatrix}}{\begin{bmatrix}E_{a}\\E_{b}\end{bmatrix}}=\tau \mathbf {E} _{\text{in}},}$ 

where the 2×2 element ${\displaystyle \tau }$  is the beam-splitter transfer matrix and r and t are the reflectance and transmittance along a particular path through the beam splitter, that path being indicated by the subscripts. (The values depend on the polarization of the light.)

If the beam splitter removes no energy from the light beams, the total output energy can be equated with the total input energy, reading

${\displaystyle |E_{c}|^{2}+|E_{d}|^{2}=|E_{a}|^{2}+|E_{b}|^{2}.}$ 

Inserting the results from the transfer equation above with ${\displaystyle E_{b}=0}$  produces

${\displaystyle |r_{ac}|^{2}+|t_{ad}|^{2}=1,}$ 

and similarly for then ${\displaystyle E_{a}=0}$ 

${\displaystyle |r_{bd}|^{2}+|t_{bc}|^{2}=1.}$ 

When both ${\displaystyle E_{a}}$  and ${\displaystyle E_{b}}$  are non-zero, and using these two results we obtain

${\displaystyle r_{ac}t_{bc}^{\ast }+t_{ad}r_{bd}^{\ast }=0,}$ 

where "${\displaystyle ^{\ast }}$ " indicates the complex conjugate. It is now easy to show that ${\displaystyle \tau ^{\dagger }\tau =\mathbf {I} }$  where ${\displaystyle \mathbf {I} }$  is the identity, i.e. the beam-splitter transfer matrix is a unitary matrix.

Expanding, it can be written each r and t as a complex number having an amplitude and phase factor; for instance, ${\displaystyle r_{ac}=|r_{ac}|e^{i\phi _{ac}}}$ . The phase factor accounts for possible shifts in phase of a beam as it reflects or transmits at that surface. Then is obtained

${\displaystyle |r_{ac}||t_{bc}|e^{i(\phi _{ac}-\phi _{bc})}+|t_{ad}||r_{bd}|e^{i(\phi _{ad}-\phi _{bd})}=0.}$ 

Further simplifying, the relationship becomes

${\displaystyle {\frac {|r_{ac}|}{|t_{ad}|}}=-{\frac {|r_{bd}|}{|t_{bc}|}}e^{i(\phi _{ad}-\phi _{bd}+\phi _{bc}-\phi _{ac})}}$ 

which is true when ${\displaystyle \phi _{ad}-\phi _{bd}+\phi _{bc}-\phi _{ac}=\pi }$  and the exponential term reduces to -1. Applying this new condition and squaring both sides, it becomes

${\displaystyle {\frac {1-|t_{ad}|^{2}}{|t_{ad}|^{2}}}={\frac {1-|t_{bc}|^{2}}{|t_{bc}|^{2}}},}$ 

where substitutions of the form ${\displaystyle |r_{ac}|^{2}=1-|t_{ad}|^{2}}$  were made. This leads to the result

${\displaystyle |t_{ad}|=|t_{bc}|\equiv T,}$ 

and similarly,

${\displaystyle |r_{ac}|=|r_{bd}|\equiv R.}$ 

It follows that ${\displaystyle R^{2}+T^{2}=1}$ .

Having determined the constraints describing a lossless beam splitter, the initial expression can be rewritten as

${\displaystyle {\begin{bmatrix}E_{c}\\E_{d}\end{bmatrix}}={\begin{bmatrix}Re^{i\phi _{ac}}&Te^{i\phi _{bc}}\\Te^{i\phi _{ad}}&Re^{i\phi _{bd}}\end{bmatrix}}{\begin{bmatrix}E_{a}\\E_{b}\end{bmatrix}}.}$ ^[2]

Applying different values for the amplitudes and phases can account for many different forms of the beam splitter that can be seen widely used.

The transfer matrix appears to have 6 amplitude and phase parameters, but it also has 2 constraints: ${\displaystyle R^{2}+T^{2}=1}$  and ${\displaystyle \phi _{ad}-\phi _{bd}+\phi _{bc}-\phi _{ac}=\pi }$ . To include the constraints and simplify to 4 independent parameters, we may write^[3] ${\displaystyle \phi _{ad}=\phi _{0}+\phi _{T},\phi _{bc}=\phi _{0}-\phi _{T},\phi _{ac}=\phi _{0}+\phi _{R}}$  (and from the constraint ${\displaystyle \phi _{bd}=\phi _{0}-\phi _{R}-\pi }$ ), so that

${\displaystyle {\begin{aligned}\phi _{T}&={\tfrac {1}{2}}\left(\phi _{ad}-\phi _{bc}\right)\\\phi _{R}&={\tfrac {1}{2}}\left(\phi _{ac}-\phi _{bd}+\pi \right)\\\phi _{0}&={\tfrac {1}{2}}\left(\phi _{ad}+\phi _{bc}\right)\end{aligned}}}$ 

where ${\displaystyle 2\phi _{T}}$  is the phase difference between the transmitted beams and similarly for ${\displaystyle 2\phi _{R}}$ , and ${\displaystyle \phi _{0}}$  is a global phase. Lastly using the other constraint that ${\displaystyle R^{2}+T^{2}=1}$  we define ${\displaystyle \theta =\arctan(R/T)}$  so that ${\displaystyle T=\cos \theta ,R=\sin \theta }$ , hence

${\displaystyle \tau =e^{i\phi _{0}}{\begin{bmatrix}\sin \theta e^{i\phi _{R}}&\cos \theta e^{-i\phi _{T}}\\\cos \theta e^{i\phi _{T}}&-\sin \theta e^{-i\phi _{R}}\end{bmatrix}}.}$ 

A 50:50 beam splitter is produced when ${\displaystyle \theta =\pi /4}$ . The dielectric beam splitter above, for example, has

${\displaystyle \tau ={\frac {1}{\sqrt {2}}}{\begin{bmatrix}1&1\\1&-1\end{bmatrix}},}$ 

i.e. ${\displaystyle \phi _{T}=\phi _{R}=\phi _{0}=0}$ , while the "symmetric" beam splitter of Loudon ^[2] has

${\displaystyle \tau ={\frac {1}{\sqrt {2}}}{\begin{bmatrix}1&i\\i&1\end{bmatrix}},}$ 

i.e. ${\displaystyle \phi _{T}=0,\phi _{R}=-\pi /2,\phi _{0}=\pi /2}$ .

Use in experiments

Beam splitters have been used in both thought experiments and real-world experiments in the area of quantum theory and relativity theory and other fields of physics. These include:

• The Fizeau experiment of 1851 to measure the speeds of light in water
• The Michelson–Morley experiment of 1887 to measure the effect of the (hypothetical) luminiferous aether on the speed of light
• The Hammar experiment of 1935 to refute Dayton Miller's claim of a positive result from repetitions of the Michelson-Morley experiment
• The Kennedy–Thorndike experiment of 1932 to test the independence of the speed of light and the velocity of the measuring apparatus
• Bell test experiments (from ca. 1972) to demonstrate consequences of quantum entanglement and exclude local hidden-variable theories
• Wheeler's delayed choice experiment of 1978, 1984 etc., to test what makes a photon behave as a wave or a particle and when it happens
• The FELIX experiment (proposed in 2000) to test the Penrose interpretation that quantum superposition depends on spacetime curvature
• The Mach–Zehnder interferometer, used in various experiments, including the Elitzur–Vaidman bomb tester involving interaction-free measurement; and in others in the area of quantum computation

Quantum mechanical description

In quantum mechanics, the electric fields are operators as explained by second quantization and Fock states. Each electrical field operator can further be expressed in terms of modes representing the wave behavior and amplitude operators, which are typically represented by the dimensionless creation and annihilation operators. In this theory, the four ports of the beam splitter are represented by a photon number state ${\displaystyle |n\rangle }$  and the action of a creation operation is ${\displaystyle {\hat {a}}^{\dagger }|n\rangle ={\sqrt {n+1}}|n+1\rangle }$ . The following is a simplified version of Ref.^[3] The relation between the classical field amplitudes ${\displaystyle {E}_{a},{E}_{b},{E}_{c}}$ , and ${\displaystyle {E}_{d}}$  produced by the beam splitter is translated into the same relation of the corresponding quantum creation (or annihilation) operators ${\displaystyle {\hat {a}}_{a}^{\dagger },{\hat {a}}_{b}^{\dagger },{\hat {a}}_{c}^{\dagger }}$ , and ${\displaystyle {\hat {a}}_{d}^{\dagger }}$ , so that

${\displaystyle \left({\begin{matrix}{\hat {a}}_{c}^{\dagger }\\{\hat {a}}_{d}^{\dagger }\end{matrix}}\right)=\tau \left({\begin{matrix}{\hat {a}}_{a}^{\dagger }\\{\hat {a}}_{b}^{\dagger }\end{matrix}}\right)}$ 

where the transfer matrix is given in classical lossless beam splitter section above:

${\displaystyle \tau =\left({\begin{matrix}r_{ac}&t_{bc}\\t_{ad}&r_{bd}\end{matrix}}\right)=e^{i\phi _{0}}\left({\begin{matrix}\sin \theta e^{i\phi _{R}}&\cos \theta e^{-i\phi _{T}}\\\cos \theta e^{i\phi _{T}}&-\sin \theta e^{-i\phi _{R}}\end{matrix}}\right).}$ 

Since ${\displaystyle \tau }$  is unitary, ${\displaystyle \tau ^{-1}=\tau ^{\dagger }}$ , i.e.

${\displaystyle \left({\begin{matrix}{\hat {a}}_{a}^{\dagger }\\{\hat {a}}_{b}^{\dagger }\end{matrix}}\right)=\left({\begin{matrix}r_{ac}^{\ast }&t_{ad}^{\ast }\\t_{bc}^{\ast }&r_{bd}^{\ast }\end{matrix}}\right)\left({\begin{matrix}{\hat {a}}_{c}^{\dagger }\\{\hat {a}}_{d}^{\dagger }\end{matrix}}\right).}$ 

This is equivalent to saying that if we start from the vacuum state ${\displaystyle |00\rangle _{ab}}$  and add a photon in port a to produce

${\displaystyle |\psi _{\text{in}}\rangle ={\hat {a}}_{a}^{\dagger }|00\rangle _{ab}=|10\rangle _{ab},}$ 

then the beam splitter creates a superposition on the outputs of

${\displaystyle |\psi _{\text{out}}\rangle =\left(r_{ac}^{\ast }{\hat {a}}_{c}^{\dagger }+t_{ad}^{\ast }{\hat {a}}_{d}^{\dagger }\right)|00\rangle _{cd}=r_{ac}^{\ast }|10\rangle _{cd}+t_{ad}^{\ast }|01\rangle _{cd}.}$ 

The probabilities for the photon to exit at ports c and d are therefore ${\displaystyle |r_{ac}|^{2}}$  and ${\displaystyle |t_{ad}|^{2}}$ , as might be expected.

Likewise, for any input state ${\displaystyle |nm\rangle _{ab}}$ 

${\displaystyle |\psi _{\text{in}}\rangle =|nm\rangle _{ab}={\frac {1}{\sqrt {n!}}}\left({\hat {a}}_{a}^{\dagger }\right)^{n}{\frac {1}{\sqrt {m!}}}\left({\hat {a}}_{b}^{\dagger }\right)^{m}|00\rangle _{ab}}$ 

and the output is

${\displaystyle |\psi _{\text{out}}\rangle ={\frac {1}{\sqrt {n!}}}\left(r_{ac}^{\ast }{\hat {a}}_{c}^{\dagger }+t_{ad}^{\ast }{\hat {a}}_{d}^{\dagger }\right)^{n}{\frac {1}{\sqrt {m!}}}\left(t_{bc}^{\ast }{\hat {a}}_{c}^{\dagger }+r_{bd}^{\ast }{\hat {a}}_{d}^{\dagger }\right)^{m}|00\rangle _{cd}.}$ 

Using the multi-binomial theorem, this can be written

${\displaystyle {\begin{aligned}|\psi _{\text{out}}\rangle &={\frac {1}{\sqrt {n!m!}}}\sum _{j=0}^{n}\sum _{k=0}^{m}{\binom {n}{j}}\left(r_{ac}^{\ast }{\hat {a}}_{c}^{\dagger }\right)^{j}\left(t_{ad}^{\ast }{\hat {a}}_{d}^{\dagger }\right)^{(n-j)}{\binom {m}{k}}\left(t_{bc}^{\ast }{\hat {a}}_{c}^{\dagger }\right)^{k}\left(r_{bd}^{\ast }{\hat {a}}_{d}^{\dagger }\right)^{(m-k)}|00\rangle _{cd}\\&={\frac {1}{\sqrt {n!m!}}}\sum _{N=0}^{n+m}\sum _{j=0}^{N}{\binom {n}{j}}r_{ac}^{\ast j}t_{ad}^{\ast (n-j)}{\binom {m}{N-j}}t_{bc}^{\ast (N-j)}r_{bd}^{\ast (m-N+j)}\left({\hat {a}}_{c}^{\dagger }\right)^{N}\left({\hat {a}}_{d}^{\dagger }\right)^{M}|00\rangle _{cd},\\&={\frac {1}{\sqrt {n!m!}}}\sum _{N=0}^{n+m}\sum _{j=0}^{N}{\binom {n}{j}}{\binom {m}{N-j}}r_{ac}^{\ast j}t_{ad}^{\ast (n-j)}t_{bc}^{\ast (N-j)}r_{bd}^{\ast (m-N+j)}{\sqrt {N!M!}}\quad |N,M\rangle _{cd},\end{aligned}}}$ 

where ${\displaystyle M=n+m-N}$  and the ${\displaystyle {\tbinom {n}{j}}}$  is a binomial coefficient and it is to be understood that the coefficient is zero if ${\displaystyle j\notin \{0,n\}}$  etc.

The transmission/reflection coefficient factor in the last equation may be written in terms of the reduced parameters that ensure unitarity:

${\displaystyle r_{ac}^{\ast j}t_{ad}^{\ast (n-j)}t_{bc}^{\ast (N-j)}r_{bd}^{\ast (m-N+j)}=(-1)^{j}\tan ^{2j}\theta (-\tan \theta )^{m-N}\cos ^{n+m}\theta \exp -i\left[(n+m)(\phi _{0}+\phi _{T})-m(\phi _{R}+\phi _{T})+N(\phi _{R}-\phi _{T})\right].}$ 

where it can be seen that if the beam splitter is 50:50 then ${\displaystyle \tan \theta =1}$  and the only factor that depends on j is the ${\displaystyle (-1)^{j}}$  term. This factor causes interesting interference cancellations. For example, if ${\displaystyle n=m}$  and the beam splitter is 50:50, then

${\displaystyle {\begin{aligned}\left({\hat {a}}_{a}^{\dagger }\right)^{n}\left({\hat {a}}_{b}^{\dagger }\right)^{m}&\to \left[{\hat {a}}_{a}^{\dagger }{\hat {a}}_{b}^{\dagger }\right]^{n}\\&=\left[\left(r_{ac}^{\ast }{\hat {a}}_{c}^{\dagger }+t_{ad}^{\ast }{\hat {a}}_{d}^{\dagger }\right)\left(t_{bc}^{\ast }{\hat {a}}_{c}^{\dagger }+r_{bd}^{\ast }{\hat {a}}_{d}^{\dagger }\right)\right]^{n}\\&=\left[{\frac {e^{-i\phi _{0}}}{\sqrt {2}}}\right]^{2n}\left[\left(e^{-i\phi _{R}}{\hat {a}}_{c}^{\dagger }+e^{-i\phi _{T}}{\hat {a}}_{d}^{\dagger }\right)\left(e^{i\phi _{T}}{\hat {a}}_{c}^{\dagger }-e^{i\phi _{R}}{\hat {a}}_{d}^{\dagger }\right)\right]^{n}\\&={\frac {e^{-2in\phi _{0}}}{2^{n}}}\left[e^{i(\phi _{T}-\phi _{R})}\left({\hat {a}}_{c}^{\dagger }\right)^{2}+e^{-i(\phi _{T}-\phi _{R})}\left({\hat {a}}_{d}^{\dagger }\right)^{2}\right]^{n}\end{aligned}}}$ 

where the ${\displaystyle {\hat {a}}_{c}^{\dagger }{\hat {a}}_{d}^{\dagger }}$  term has cancelled. Therefore the output states always have even numbers of photons in each arm. A famous example of this is the Hong–Ou–Mandel effect, in which the input has ${\displaystyle n=m=1}$ , the output is always ${\displaystyle |20\rangle _{cd}}$  or ${\displaystyle |02\rangle _{cd}}$ , i.e. the probability of output with a photon in each mode (a coincidence event) is zero. Note that this is true for all types of 50:50 beam splitter irrespective of the details of the phases, and the photons need only be indistinguishable. This contrasts with the classical result, in which equal output in both arms for equal inputs on a 50:50 beam splitter does appear for specific beam splitter phases (e.g. a symmetric beam splitter ${\displaystyle \phi _{0}=\phi _{T}=0,\phi _{R}=\pi /2}$ ), and for other phases where the output goes to one arm (e.g. the dielectric beam splitter ${\displaystyle \phi _{0}=\phi _{T}=\phi _{R}=0}$ ) the output is always in the same arm, not random in either arm as is the case here. From the correspondence principle we might expect the quantum results to tend to the classical one in the limits of large n, but the appearance of large numbers of indistinguishable photons at the input is a non-classical state that does not correspond to a classical field pattern, which instead produces a statistical mixture of different ${\displaystyle |n,m\rangle }$  known as Poissonian light.

Rigorous derivation is given in the Fearn–Loudon 1987 paper^[4] and extended in Ref ^[3] to include statistical mixtures with the density matrix.

Non-symmetric beam-splitter

In general, for a non-symmetric beam-splitter, namely a beam-splitter for which the transmission and reflection coefficients are not equal, one can define an angle ${\displaystyle \theta }$  such that

${\displaystyle {\begin{cases}|R|=\sin(\theta )\\|T|=\cos(\theta )\end{cases}}}$ 

where ${\displaystyle R}$  and ${\displaystyle T}$  are the reflection and transmission coefficients. Then the unitary operation associated with the beam-splitter is then

${\displaystyle {\hat {U}}=e^{i\theta \left({\hat {a}}_{a}^{\dagger }{\hat {a}}_{b}+{\hat {a}}_{a}{\hat {a}}_{b}^{\dagger }\right)}.}$ 

Application for quantum computing

In 2000 Knill, Laflamme and Milburn (KLM protocol) proved that it is possible to create a universal quantum computer solely with beam splitters, phase shifters, photodetectors and single photon sources. The states that form a qubit in this protocol are the one-photon states of two modes, i.e. the states |01⟩ and |10⟩ in the occupation number representation (Fock state) of two modes. Using these resources it is possible to implement any single qubit gate and 2-qubit probabilistic gates. The beam splitter is an essential component in this scheme since it is the only one that creates entanglement between the Fock states.

Similar settings exist for continuous-variable quantum information processing. In fact, it is possible to simulate arbitrary Gaussian (Bogoliubov) transformations of a quantum state of light by means of beam splitters, phase shifters and photodetectors, given two-mode squeezed vacuum states are available as a prior resource only (this setting hence shares certain similarities with a Gaussian counterpart of the KLM protocol).^[5] The building block of this simulation procedure is the fact that a beam splitter is equivalent to a squeezing transformation under partial time reversal.

Diffractive beam splitter

This section is an excerpt from Diffractive beam splitter.[edit]

 

7x7 matrix using green laser and diffractive beam splitter.

The diffractive beam splitter^[6][7]

(also known as multispot beam generator or array beam generator) is a single optical element that divides an input beam into multiple output beams.^[8] Each output beam retains the same optical characteristics as the input beam, such as size, polarization and phase. A diffractive beam splitter can generate either a 1-dimensional beam array (1xN) or a 2-dimensional beam matrix (MxN), depending on the diffractive pattern on the element. The diffractive beam splitter is used with monochromatic light such as a laser beam, and is designed for a specific wavelength and angle of separation between output beams.

See also

• Power dividers and directional couplers

References

 

Wikimedia Commons has media related to Beam splitters.

1. Zetie, K P; Adams, S F; Tocknell, R M, How does a Mach–Zehnder interferometer work? (PDF), retrieved 13 February 2014
2. R. Loudon, The quantum theory of light, third edition, Oxford University Press, New York, NY, 2000.
3. Campos, Richard; Bahaa, Saleh; Malvin, Teich (Aug 1, 1989). "Quantum mechanical lossless beam splitter: SU(2) symmetry and photon statistics". Physical Review A. 40 (3): 1371.
4. Fearn, H.; Loudon, R. (1987). "Quantum theory of the lossless beam splitter". Optics Communications. 64 (6): 485–490. doi:10.1016/0030-4018(87)90275-6.
5. Chakhmakhchyan, Levon; Cerf, Nicolas (2018). "Simulating arbitrary Gaussian circuits with linear optics". Physical Review A. 98: 062314. arXiv:1803.11534. doi:10.1103/PhysRevA.98.062314.
6. Diffraction Gratings and Applications, Loewen, Erwin C. and Popov, Evgeny. Marcel Dekker, Inc. 1997.
7. Digital diffractive optics: an introduction to planar diffractive optics and related technology, Bernard C. Kress, Patrick Meyrueis , 2005.
8. Diffractive Optics – Design, Fabrication and Test, O'Shea, Suleski, Kathman and Prather, 2004. p.83