16 December 2016 Dark GPC: extended nodal beam areas from binary-only phase
Author Affiliations +
Optical Engineering, 55(12), 125102 (2016). doi:10.1117/1.OE.55.12.125102
We show a simplified method of generating extended regions of destructive interference with near arbitrary shapes using the generalized phase contrast (GPC) method. For Gaussian input beams, GPC typically results in a intensified user-defined input mask shape against a dark background. In this work, we investigate conditions wherein GPC’s synthetic reference wave destructively interferes with what is typically the foreground pattern. Using alternate conditions for the input phase mask, the locations of light and darkness are interchanged with respect to typical GPC output mappings. We show experimentally how “dark GPC” allows the dark regions to be easily reshaped using a binary-only phase mask encoded on a spatial light modulator. Similar to standard GPC, the method does not require complex calculations or the fabrication of complex gray-level phase elements. The simplified approach and flexibility in the output shapes make dark GPC attractive for applications such as optical trapping of low-index particles or superresolution microscopy like stimulated emission depletion.
Villangca, Bañas, Palima, and Glückstad: Dark GPC: extended nodal beam areas from binary-only phase



Structured light distributions have been used in applications such as optical trapping and manipulation,12.3 optical sorting,4,5 advanced microscopy, and selective uncaging or excitation6 in neurophotonics or optogenetics.7 In some applications, however, it is desirable to divert light from a particular region, and “light shaping” techniques are used to create well-defined patterns of darkness, on top of a bright background. For example, bounded regions of darkness are used in defining potential distributions for atom trapping, optical trapping of low-index media,89.10 and in stimulated emission depletion microscopy.11

Another important imaging application is coronography in astronomy. Modern uses divert light from an overpowering bright source (e.g., star) to increase the visibility of some weaker signals (e.g., an orbiting planet). Early coronagraphs use an amplitude mask and were originally designed by Lyot in the 1930s. However, Lyot’s amplitude masking has inherent limitations, such as completely blocking light within its physical extent and thereby inherently blocking objects of interest.12 A more light-efficient alternative is to implement phase-only coronography where undesired light can be canceled by destructive interference, as is the case for so-called common path interferometers. The phase-only approach typically uses a disk-shaped phase mask filter in the focal plane of a telescope. The disk-shaped phase mask filter is positioned can introduce a π-phase shift to a small portion of the Airy disk pattern.13 The unperturbed and phase-shifted portions of the incident light can then destructively interfere at the exit pupil of the coronagraphic setup. Instead of a single-phase shifting region, an improvement in the design uses a quadrant arrangement of π-phase shifting regions, which is more robust compared with the disk phase mask.14 A vortex spatial filter has been later introduced by Swartzlander.15 A subwavelength surface-relief grating coined the annular groove phase mask by Mawet et al.12 was also used to create a spiral phase of topological charge lp=2. Mawet et al. have also analytically shown the formation of a nodal area.

In the aforementioned examples, the geometry of the resulting nodal area is circular. Currently, techniques for creating structured darkness typically use optical vortices created by static phase plates16 or computer-generated holograms.17 Laguerre–Gaussian beams or other similar beams containing singularity have limited control over the lateral shape of the dark region. Having control over the darkness shapes while maintaining sharp transitions from dark to bright brings more possibilities and flexibility. For example, engineering such dark regions to shapes other than circular can lead to specifically designed dipole traps.18 Potential uses in laser-materials processing can include shielding or cutting out user-defined shapes.

Recently, a 4f approach using optimized phase elements at the Fourier plane has been demonstrated.19 Although successful in creating so-called “dark nodal areas,” the complexity of designing and fabricating such vortex-like phase elements can be prohibitive for practical applications. In principle, any light shaping technique, even photon-inefficient amplitude modulation, can be used to form voids on top of a background of light. Furthermore, one could always mathematically and philosophically argue that any darkness is a result of destructive interference. We set our work apart by having a clearer physical description in which waves are actually destructively interfering. Further, we deliberately tweak the superposed waves, so they are out-of-phase within desired patterns that can take up a wide variety of shapes.



A straightforward approach in creating a dark region is amplitude masking. This approach can be extended to project arbitrary and dynamic patterns using a spatial light modulator (SLM) such as a digital micromirror device (DMD).20,21 The fast switching of DMD elements enables rapid adaptive beam shaping and can project grayscale images as well. Phase-only masking is an alternative way to beam shaping and can achieve much higher photon efficiencies for light projection. Phase-only masking is performed by modifying an incoming light wavefront, so constructive interference occurs at desired target locations at the output. To create a dark region by phase-only masking, a portion of the incident beam is made to be out-of-phase to form destructive interference when recombined at the output projection plane. The creation of a dark core by destructive interference is fairly easy to grasp physically if we think of it as sum of areas. If for example π-phase-shifted regions are regarded as negative areas and unmodulated regions as positive areas, then they will result in a net sum of zero after performing an optical Fourier transformation. A phase filter with π-phase-shifted radially opposite regions added at the focal plane will, however, result in a destructive interference within the geometric area of the pupil after taking the Fourier transform. This has been observed in the four-quadrant, vortex, and annular groove phase filters mentioned earlier. The disk phase mask is different as it makes the “opposite” areas cancel not straightforward. Our method similarly uses a common path interferometer setup, where we show that we can switch operation between constructive to destructive interference by placing appropriate binary masks and filters. We present calculations based on on-axis amplitude matching to obtain parameters to aid in the production of the masks and filters and to adapt to arbitrary patterns.

The generalized phase contrast (GPC) method is normally used for efficient phase-only shaping of light into speckle-free contiguous optical distributions. The method is useful for applications such as static beam shaping, optical manipulation, or excitation in two-photon optogenetics. It operates by synthesizing a reference wave that interferes with an imaged phase-modulated input light. The input (0π) phase modulation can have near-arbitrary or user-defined two-dimensional distributions and directly governs where constructive and destructive interferences happen at the output. Although seldom emphasized, standard “bright GPC” works by creating an extended region of destructive interference that is negative to the shaped foreground. Furthermore, one can “poke” dark holes on top of a bright GPC pattern, offering greater freedom on the location of the dark regions as used, for example, to manipulate microparticles with low refractive indices.9,22 In this work, however, we maximize the area of destructive interference by finding alternative modalities and conditions that make the synthetic reference wave (SRW) out-of-phase with the imaged foreground pattern. A practical advantage of having a larger darkness region is the utilization of more pixels in the SLM, which leads to more detailed or well-defined darkness regions. The higher resolution may make the dark patterns practical for laser materials processing, for example.

A standard GPC light shaping setup derives from a 4f imaging configuration wherein a laser-illuminated phase-only mask is transformed into a corresponding shaped optical intensity. GPC maps the phase input into an intensity output using a binary phase contrast filter (PCF) situated at the Fourier plane that is responsible for forming the SRW. The SRW, formed from the phase-shifted low-frequency components of the input, interferes with the 4f-imaged input phase mask, hence increasing the intensity within the mask’s shaped foreground region while forming darkness outside this region. The interfering fields can be seen clearly in Fig. 1. Mathematically, the GPC output can be written as


In this manner, it is easy to see that constructive or destructive interference will depend on the balance of the amplitudes between the two terms in Eq. (1). The amplitude of the SRW is dependent on the size of the phase mask and PCF and implies that an optimal pairing can be found. To optimize the light throughput (efficiency) in the shaped region, we have previously derived the optimally matching sizes of the phase mask and PCF given the laser source’s Gaussian half-waist, w0, wavelength, λ, and Fourier lens’s focal length, f.6 In practical usage, the PCF size is held fixed while the input phase mask is swappable or programmable through a dynamically reconfigurable SLM. Our theoretical analysis has previously shown that the optimal PCF phase-shifting radius, Δrf, should be 1.1081 times the Fourier Gaussian half-waist, wf. A matching phase mask size has also been identified for different shapes and has been demonstrated experimentally to yield efficiently shaped beams.23

Fig. 1

A GPC configuration generating an SRW that destructively interferes with imaged input foreground pattern.


For a GPC setup, generating darkness means finding configurations for which the SRW destructively interferes with the shaped foreground region. Since the PCF is usually fixed in our applications, we look for phase mask distributions, φ(x,y), that will cause Eq. (1) to nullify at the phase-encoded shaped region. A quick look on the effect of the size of the phase mask distribution for a circular pattern of radius ΔR=ζw0 is shown in Fig. 2. At some radius corresponding to ζ1.02, destructive interference occurs, forming a nodal area.

Fig. 2

GPC output central intensity as a function of the relative radius, ζ=ΔR/w0, of a circular phase mask for a fixed PCF (η=Δrf/wf=1.1081). Inset (Video 1): profile of the GPC output at ζ value corresponding to the lowest central intensity with line scan taken horizontally across the center of the output pattern (Video 1, QuickTime Mov, 954 KB [URL: http://dx.doi.org/10.1117/1.OE.55.12.125102.1]).


From an engineering standpoint, it will be helpful to have a working relation between phase mask and PCF given in terms of the parameters η and ζ. In general, an analytic form for the SRW cannot be found easily, especially since the phase mask can take in arbitrary nonanalytic forms. However, for slowly varying phase mask distribution such as a circle, it can be approximated by a sum of Gaussian functions. Consider an input beam such as a laser beam with a Gaussian profile a(r). When the laser beam passes through a phase mask containing the target pattern φ(r), it is modulated as Ein(r)=a(r)exp[iφ(r)]. For a circular phase mask, the modulation takes the form of φ(r)=πcirc(r/ζw0), and the input field becomes Ein(r)=a(r)2circ(r/ζw0)a(r). The first lens takes the Fourier transform of Ein and is multiplied by the transfer function of the PCF. The circular π-phase-shifting region of the PCF has a radius of Δrf corresponding to frequency cutoff Δfr=ηf0, and thus the transfer function is given by H(fr)=12circ(fr/ηf0). The PCF can be seen as a low-pass filter where the low-spatial frequency components are π-phase shifted. The field is once again Fourier transformed by the second lens and the output is given by


The term ESRW(r)=F1{2circ(fr/ηf0)·F[Ein(r)]} represents the SRW. For an input circular phase mask, Ein(r) will have a circ function in the phase term and be multiplied with a Gaussian profile, and the ESRW(r) will involve convolution of Gaussian and jinc functions. To simplify the calculation, the jinc is approximated with a Gaussian function. For circ(fr/ηf0), which transforms into 2πη2f02jinc(2πηf0r), the Gaussian approximation is given by πη2f02exp(π2η2f02r2/d02), where d0=0.37π. The parameter d0 is chosen to match the central part of the jinc. The convolution is now between two Gaussian functions and simply results in a broadened Gaussian function. However, this does not account for the negative values in the original jinc function, and the amplitude should be corrected. The correct central amplitude can be analytically obtained6 as πf0[1exp(η2)]. This correction is done for every instance of Fourier transform of a bounded Gaussian function. Finally, we arrive at the following approximate equation for the GPC output:


where kη2=(η2+d02)/η2, kζ2=(ζ2+d02)/ζ2 and kη/kζ2=(η2/kη2+d02)/η2/kζ2. Assuming the input field has a unit on-axis amplitude, we also require the SRW to have a value equal to one. Thus, the “darkness” condition can be written as ESRW(0)=Ein(0). For the circular π-phase shifted mask that we used, the input amplitude is Ein(0)=1; thus,


To see whether Eq. (4) will lead to the creation of a dark area, we compute for the efficiency, which we define as the power within the area defined by the phase pattern divided by the incident power. Dark nodal areas are therefore those resulting in low efficiency. An efficiency map is shown in Fig. 3, where an overlay of Eq. (4) shows that on-axis amplitude matching is sufficient to identify combinations of parameters that will result in darkness. The parameters that will result in bright GPC can also be seen. The insets in Fig. 3 show the calculated GPC outputs for selected ζη pairs. Interestingly, there are parameters where the on-axis SRW amplitude is zero and the Gaussian input reappears.

Fig. 3

Efficiency map for different η and ζ pairs. The bands show the regions of bright and dark GPCs. An overlay of Eq. (4) represented as contour lines for ESRW(0)=  1 and 0 shows that on-axis amplitude matching can be used to get optimal parameters. Line scans of the intensity profiles are taken horizontally across the center.



Experiments and Results

The dynamic GPC setup used for the dark GPC experimental verifications is described in an earlier work23 and illustrated in Fig. 4. A liquid crystal on silicon SLM (Hamamatsu Photonics) is used to generate arbitrary binary input phase patterns. The SLM has an area of 16×12  mm2 and pixel pitch of 20  μm. For this proof-of-principle demonstration, the SLM has been illuminated with a 532-nm filtered super-continuum laser, polarized and expanded such that 2w0=4  mm (corresponding to 200 SLM pixels). The Fourier lens used has a focal length of 100 mm, and the PCF’s radius, Δrf, is 9.4  μm. To form arbitrary patterns of darkness, binary bitmap images were directly drawn on the SLM window and mapped to 0 and π phase shifts. The images were scaled, so the balance between the imaged input beam and the SRW can be attained. Alternatively, the phase filter can be modified by means of iterative optimization but requires a dynamic phase filter.19 Our approach is much simpler. The experimental results are shown in Fig. 5. The results show that noncircular dark patterns are produced by designing them to have the same zero-order strength as the circular case.

Fig. 4

Schematic illustration of the dynamic GPC light shaper used for generating dark GPC patterns (schematic based on Ref. 23) and the actual experimental setup.


Fig. 5

Dark GPC experimental results showing different darkness shapes.


The sensitivity of dark GPC makes it useful for positioning the PCF in a GPC light shaper using a calibration mask.23 However, this sensitivity also visualizes faint amounts of light, due to the combined phase aberrations in the setup. The presence of extra rings outside the dark regions can be solved depending on the intended applications. In coronography, it is common practice to add a Lyot stop at the exit pupil plane to remove these.



We have shown how to generate extended regions of destructive interference using a GPC light shaping system by finding alternate conditions wherein the SRW is out-of-phase with the corresponding input foreground. Our direct binary mapping approach for generating dark regions is far simpler than previously reported techniques and thus allows us to experimentally demonstrate arbitrarily shaped dark regions. Instead of using multilevel phase elements or iterative optimizations techniques,19 we take advantage of GPC’s use of a simpler and reusable binary PCF and interchangeable binary phase input patterns. As GPC can operate over a broad wavelength range,24 it would be interesting to see whether dark GPC exhibits the same robustness and can be used for multispectral applications found in microscopy.


This work has been financially supported by the Enhanced Spatial Light Control Advanced Optical Fibres (e-space) project by Innovation Fund Denmark and GAP funding from the Copenhagen Cleantech Cluster (CCC). We thank the following for lending us equipment: Hamamatsu Photonics Central Research Laboratory for the spatial light modulator and NKT Photonics for the super-continuum laser. Portions of this work were presented at the Proc. SPIE 9764, Complex Light and Optical Forces X in 2016, entitled “Dark GPC.”



H.-U. Ulriksen et al., “Independent trapping, manipulation and characterization by an all-optical biophotonics workstation,” J. Eur. Opt. Soc. Rapid Publ. 3, 8034 (2008).http://dx.doi.org/10.2971/jeos.2008.08034Google Scholar


J. Glückstad et al., “Wave-guided optical waveguides tracked and coupled using dynamic diffractive optics,” in Workshop on Information Optics (WIO 2014), Neuchâtel, Switzerland (2014).Google Scholar


K. Dholakia and T. Čižmár, “Shaping the future of manipulation,” Nat. Photonics 5(6), 335–342 (2011).http://dx.doi.org/10.1038/nphoton.2011.80Google Scholar


A. Bañas et al., “Cell sorting using efficient light shaping approaches,” Proc. SPIE 9764, 97640F (2016).http://dx.doi.org/10.1117/12.2214165Google Scholar


I. Perch-Nielsen et al., “Parallel particle identification and separation for active optical sorting,” J. Opt. A: Pure Appl. Opt. 11(3), 034013 (2009).http://dx.doi.org/10.1088/1464-4258/11/3/034013Google Scholar


A. Bañas et al., “GPC light shaper for speckle-free one- and two-photon contiguous pattern excitation,” Opt. Express 7102(5), 5299–5310 (2014).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.22.005299Google Scholar


E. Papagiakoumou et al., “Scanless two-photon excitation of channelrhodopsin-2,” Nat. Methods 7(10), 848–854 (2010).1548-7091http://dx.doi.org/10.1038/nmeth.1505Google Scholar


K. T. Gahagan and G. A. Swartzlander, “Optical vortex trapping of particles,” Opt. Lett. 21(11), 827–829 (1996).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.21.000827Google Scholar


P. J. Rodrigo, V. R. Daria and J. Glückstad, “Real-time interactive optical micromanipulation of a mixture of high- and low-index particles,” Opt. Express 12(7), 1417–1425 (2004).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.12.001417Google Scholar


V. R. Daria, M. A. Go and H.-A. Bachor, “Simultaneous transfer of linear and orbital angular momentum to multiple low-index particles,” J. Opt. 13(4), 044004 (2011).http://dx.doi.org/10.1088/2040-8978/13/4/044004Google Scholar


S. W. Hell and J. Wichmann, “Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy,” Opt. Lett. 19(11), 780–782 (1994).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.19.000780Google Scholar


D. Mawet et al., “Annular groove phase mask coronagraph,” Astrophys. J. 633, 1191–1200 (2005).ASJOAB0004-637Xhttp://dx.doi.org/10.1086/462409Google Scholar


F. Roddier and C. Roddier, “Stellar coronagraph with phase mask,” Publ. Astron. Soc. Pac. 109, 815–820 (1997).ASPSFO1050-3390Google Scholar


D. Rouan et al., “The four quadrant phase mask coronagraph. I. Principle,” Publ. Astron. Soc. Pac. 112(777), 1479–1486 (2000).PASPAU0004-6280http://dx.doi.org/10.1086/317707Google Scholar


G. A. Swartzlander, “Peering into darkness with a vortex spatial filter,” Opt. Lett. 26(8), 497–499 (2001).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.26.000497Google Scholar


V. V. Kotlyar et al., “Generation of phase singularity through diffracting a plane or Gaussian beam by a spiral phase plate,” J. Opt. Soc. Am. A. Opt. Image Sci. Vis. 22(5), 849–861 (2005).http://dx.doi.org/10.1364/JOSAA.22.000849Google Scholar


N. R. Heckenberg et al., “Generation of optical phase singularities by computer-generated holograms,” Opt. Lett. 17(3), 221 (1992).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.17.000221Google Scholar


J. G. Lee et al., “Analogs of basic electronic circuit elements in a free-space atom chip,” Sci. Rep. 3, 1034 (2013).SRCEC32045-2322http://dx.doi.org/10.1038/srep01034Google Scholar


G. J. Ruane et al., “Nodal areas in coherent beams,” Optica 2(2), 147–150 (2015).http://dx.doi.org/10.1364/OPTICA.2.000147Google Scholar


J. Liang et al., “High-precision laser beam shaping using a binary-amplitude spatial light modulator,” Appl. Opt. 49(8), 1323–1330 (2010).APOPAI0003-6935http://dx.doi.org/10.1364/AO.49.001323Google Scholar


J. Liang, “Grayscale laser image formation using a programmable binary mask,” Opt. Eng. 51, 108201 (2012).http://dx.doi.org/10.1117/1.OE.51.10.108201Google Scholar


V. R. Daria, P. J. Rodrigo and J. Glückstad, “Dynamic array of dark optical traps,” Appl. Phys. Lett. 84(3), 323 (2004).APPLAB0003-6951http://dx.doi.org/10.1063/1.1642752Google Scholar


A. Bañas et al., “GPC light shaper: static and dynamic experimental demonstrations,” Opt. Express 22(20), 23759–23769 (2014).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.22.023759Google Scholar


O. Kopylov et al., “GPC light shaping a supercontinuum source,” Opt. Express 23(3), 1894–1905 (2015).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.23.00184Google Scholar


Mark Jayson Villangca received his master’s degree in physics from the University of the Philippines, while working on beam shaping using computer-generated holograms. He received his PhD in photonics engineering from the Technical University of Denmark as a member of the Programmable Phase Optics group while working on light-driven microrobots. He also works with generalized phase contrast and digital holography for beam shaping.

Andrew Rafael Bañas received his PhD from DTU Fotonik. His works with the Programmable Phase Optics group have been featured in Optics Express and OPN. These include applying Fourier optics or electrodynamics to maximize experiments. He has also designed and built hardware and software for beam shaping and optical manipulation systems. He is currently pursuing tech-transfer activities that include applications of generalized phase contrast and cell sorting for studying disease.

Darwin Palima received his PhD in physics from the University of the Philippines and moved to Denmark to work as a postdoc. He works as an associate professor at DTU Fotonik, Technical University of Denmark. He coauthored a book on generalized phase contrast and actively publishes in peer-reviewed journals and conference proceedings. He teaches biophotonics and optical engineering while pursuing research interests that include computer-generated holograms, generalized phase contrast, optical trapping and micromanipulation, microscopy, and biophotonics applications.

Jesper Glückstad received his PhD and DSc degrees and is professor at DTU Fotonik and former guest professor in biophotonics at LIT. He established ppo.dk and received the DOPS Award in 2000. Elected “Scientist of the Year” by Ib Henriksen’s Foundation in 2005. He is fellow of OSA and SPIE as the first from Denmark based on his invention of GPC and its derivatives published in a Springer book in 2009. Founder of OptoRobotix.com and GPCphotonics.com.

Mark Jayson Villangca, Andrew R. Bañas, Darwin Palima, Jesper Glückstad, "Dark GPC: extended nodal beam areas from binary-only phase," Optical Engineering 55(12), 125102 (16 December 2016). http://dx.doi.org/10.1117/1.OE.55.12.125102
Submission: Received ; Accepted

Destructive interference

Spiral phase plates

Spatial light modulators

Beam shaping

Phase shift keying

Phase contrast

Binary data

Back to Top