Toward a Nanophotonic Nose: A Compressive Sensing-Enhanced

Differences between “traditional” compressive sensing and the approach taken here. Example numerical code for compressive sensing reconstruction (...
0 downloads 0 Views 1MB Size
Subscriber access provided by TULANE UNIVERSITY

Letter

Towards a nanophotonic nose: a compressive sensingenhanced, optoelectronic mid-infrared spectrometer Benjamin Cerjan, and Naomi J. Halas ACS Photonics, Just Accepted Manuscript • DOI: 10.1021/acsphotonics.8b01503 • Publication Date (Web): 19 Dec 2018 Downloaded from http://pubs.acs.org on December 20, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Towards a nanophotonic nose: a compressive sensing-enhanced, optoelectronic mid-infrared spectrometer Benjamin Cerjan1 and Naomi J. Halas1,2,3,4 1

Department of Physics and Astronomy, 2Department of Electrical and Computer Engineering, 3 Department of Chemistry, and 4Laboratory for Nanophotonics, Rice University, 6100 Main St, Houston, Texas 77005, United States

Abstract Infrared (IR) spectroscopy has been a central tool for chemical analysis for decades, useful in a wide range of fields for the detection and quantification of molecules based on their unique vibrational resonances. Conventional IR spectroscopy relies on bulky, dispersive optics, however, making portability and miniaturization a substantial challenge. Here we demonstrate a micronscale IR spectrometer where spectrally selective detection is performed optoelectronically, based on the wavelength-dependent mid-IR photocurrent responses of an array of Al grating-based detectors fabricated on a doped Si substrate. Compressive sensing techniques extend our resolution, enabling spectral features to be identified with a remarkably small number of detectors. This work demonstrates a CMOS-compatible, readily scalable approach for the fabrication of compact, room-temperature IR spectrometers capable of use in fieldable applications. Keywords: Spectroscopy, infrared, compressive sensing, aluminum, super-resolution, miniature

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The field of miniature integrated spectroscopic devices is an extremely active research area due to the continual drive to reduce size, decrease cost, and increase portability while maintaining functionality. A wide range of novel devices has been investigated in this pursuit including colloidal quantum dots,1,2 pyroelectric devices,3 metal-insulator-semiconductor detectors,4 plasmonic lab-on-a-chip devices,5 graphene photodetectors,6 pixel arrays for chemical imaging,7 plasmonic filter arrays,8 and a variety of other designs.9–15 These types of devices have demonstrated that some, if not all, functional components of a spectrometer can be effectively miniaturized across a broad spectral range, spanning the ultraviolet to the infrared (IR). The midIR “chemical fingerprinting” region, from 2-12 µm, is of particular interest due to the presence of molecular vibrational resonances. These resonances provide signatures for the detection, identification, and quantification of almost any molecular structure. A detector capable of direct, wavelength-dependent detection in the mid-IR, with room-temperature operation, would be highly useful in a wide range of applications including environmental,16 geological,17 and food safety18 monitoring. In particular, there has been substantial effort invested to design a “nanophotonic nose” – an optical device with direct electrical readouts capable of performing gas-phase molecular detection and identification.19–23 In this work we demonstrate an optoelectronic room temperature micron-scale mid-IR spectrometer fabricated using CMOS-compatible materials. Leveraging a recently developed approach to mid-IR sensing,24 we show that by multiplexing the response from as few as 20 detector elements we are capable of accurately reproducing test input spectra with a resolution of 40 nm. Previous work using nanoscale filter elements in the mid-IR have used both a direct peakenvelope7 or a least-squares8 method for reconstructing an unknown input spectrum. By using techniques from the field of compressive sensing25–28 we are able to enhance our resolvable

ACS Paragon Plus Environment

Page 2 of 23

Page 3 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

resolution well below the Shannon-Nyquist limit even with a small number of detectors.29 We also use our spectrometer for contact-free molecular detection of a polymer film placed in the excitation path of the spectrometer, demonstrating a proof-of-principle device capable of performing as a “nanophotonic nose”. Finally, we use simulations to estimate the optimal number of gratings for a complete vibrational spectral reconstruction using this type of optoelectronic spectrometer.

Results Device fabrication and photodetection Our spectrometer, schematically shown in Figure 1, utilizes plasmonic (high purity) aluminum gratings on highly p-doped silicon substrates. The photodetection mechanism used here has been previously reported,24 but, in brief, the detector response is based on the change in device resistance due to free-carrier absorption in the doped silicon substrate. By applying a small voltage and measuring the change in current, a photoresponse due to incident light can be detected. The aluminum gratings here serve a dual purpose: first as plasmonic gratings that filter which wavelengths of light are transmitted through to the substrate, and second, as electrical contacts for wire-bonding to an external chip carrier. This simple combination allows for a spectrometer composed of many individual gratings to be treated in an ensemble.

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The gratings were fabricated using standard electron-beam lithography and clean room procedures (full details in Methods). The gratings in this work range in pitch from 1250-3600 nm with a gap of 20% of the total pitch (e.g. for a 1250 nm pitch there is a 250 nm gap) and a thickness of 300 nm. Following fabrication, they were wire-bonded to a chip carrier to ensure a consistent

Figure 1: Layout and schematic of the optoelectronic spectrometer used in this work. (a) Rendering of the gratings showing that all gratings are oriented in the same direction and how a measurement is taken with the input light focused onto each grating individually (red beam). (b) Schematic representation of the grating array showing the polarization of the incident light and how the gratings change in pitch when going around the array. (c) Optical image of the grating array used in this work (gratings range from 1250–3600 nm in pitch). (d) SEM image of the grating array. (e) Representative responsivity curves for four gratings with pitch 2000 nm (red), 2375 nm (blue), 2625 nm (purple), and 3000 nm (green), as indicated in (b). The curves are normalized relative to the largest recorded responsivity for any of the gratings.

electrical connection through the device. The unique responsivity of each grating as a function of wavelength (e.g. Figure 1(e)) was characterized consecutively by measuring the current through

ACS Paragon Plus Environment

Page 4 of 23

Page 5 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

each grating with a lock-in amplifier while the grating was illuminated with a focused tunable midIR laser polarized perpendicular to the grating (for all responsivity curves see Figure S1). As each grating has a different pitch, they all possess different responsivity curves as a function of wavelength. Once the response curve as a function of wavelength for each grating was measured, we found we were then able to create a proof-of-principle mid-IR spectrometer using as few as 20 gratings. Spectral reconstruction Reconstructing an unknown input signal directed onto a detector element (or elements) requires that the response from the detector be correlated to a particular input intensity at a given wavelength. This requires either a single detector element coupled with a mechanism for wavelength selection, such as a grating or prism, or a large number of detector elements each with

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

different wavelength-dependent response characteristics. The approach taken here is to use the readout from a number of detector elements, each with a unique response as a function of wavelength. To use a collection of detector elements as a spectrometer, the mathematical problem posed is the following: how do we extract an unknown input signal, 𝑥⃑, given the measured response of each device as a result of this input signal, 𝑦⃑, and the previously measured responsivity curves stored in the matrix D? This process is schematically presented in Figure 2. First, an unknown input spectrum is directed onto a single grating and the current through that grating is recorded.

Figure 2: Reconstruction process for an input signal. An unknown input spectrum, 𝑥⃑, is directed on to each grating. Each grating has a unique responsivity function contained in the matrix 𝑫, and generates a current yi in response to the incident light. These currents are then used to estimate the original input spectrum.

This process is then repeated for all of the gratings (20 in this work). Finally, all of the current measurements are combined to reconstruct the original input spectrum. Mathematically, this corresponds to the equation

ACS Paragon Plus Environment

Page 6 of 23

Page 7 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

𝐃𝑥⃑ − 𝑦⃑ = 0 ,

(1)

where 𝑥⃑ is a representation of a continuous input spectrum 𝑥(𝜆) sampled at m points and 𝑦⃑ is our set of m current measurements, making both of size m × 1. D, of size m × m, contains the measured responsivity curves and so, ideally, this equation can be thought of as integrating our input signal times the response function of each of the gratings at a discrete number of wavelengths. To find

−𝑦=0 ,

(1), ⇒ 𝑥⃑ = 𝐃−1 𝑦⃑ ,

assuming that D is invertible. Physically, the condition that D is invertible corresponds to the responsivity functions of the different detector elements being distinct from one another. This is simply a mathematical restatement of what was postulated above, that the detector elements must all have unique response functions in order to be able to reconstruct arbitrary input signals. In practice, since the noise of an experimental measurement typically prevents an exact solution, the most common numerical approach to solving this problem is through non-negative least-squares regression,30 𝑚𝑖𝑛 || 𝐃𝑥⃑ − 𝑦⃑||22 . 𝑥𝑖 ≥0

(2)

In this numerical approach we are looking for the vector, 𝑥⃑, which minimizes the squared value of the L2-norm, || 𝐃𝑥⃑ − 𝑦⃑||2 = ‖𝑧⃑‖2 ≝ √(𝑧1 2 + 𝑧2 2 + ⋯ ) 𝑧1 𝑧⃑ = ( ⋮ ) , 𝑧𝑚

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

subject to the constraint that every element of 𝑥⃑ is greater than or equal to 0, which holds here since the input spectrum intensity is intrinsically positive or 0. Although least-squares minimization is computationally straightforward and efficient, this technique is limited in spectral resolution by the length of 𝑦⃑, that is, the number of unique current measurements. To keep the solution tractable, our final spectrum, 𝑥⃑, can only have as many data points as there are current measurements in 𝑦⃑. If we tried to increase the number of points in our reconstructed spectrum beyond the number of current measurements, we would end up with an underdetermined system that is unsolvable (we would have a system of m equations and n unknowns, with 𝑛 > 𝑚). This means that 𝑥⃑ and 𝑦⃑ are both of length m, and D is m × m, where m is the number of grating elements, since each grating produces a single current measurement. In our system, which spans the spectral range 2200-3800 nm and has m = 20, this corresponds to a spectral window

final reconstructed resolution of 80 nm (number of points) as displayed in Figure 3. If we try to distinguish different input spectra with peak positions closer together than this limit, we are unable to tell them apart using this method of spectral reconstruction, as shown in Figure 3(d). Thus, the most obvious method to improve the spectral resolution is to simply increase the number of

ACS Paragon Plus Environment

Page 8 of 23

Page 9 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

gratings (increasing the length of 𝑦⃑). The improvement in resolution is linear with the increase in number of gratings, and simulations demonstrating this improvement are presented in the next section. However, another mathematical approach may be employed to obtain an increase in spectral resolution such that we are able to achieve sub-Nyquist resolution without increasing the number of gratings.

Figure 3: Spectral reconstruction using both L1 and L2 methods with comparisons to ground truth. Input signals centered at 2506 nm and 3358 nm (a), 3105 nm (b), and 3145 nm (c). While L2 methods are unable to distinguish the peaks in (b) and (c), L1 reconstruction methods are easily able to tell them apart as seen in (d).

In the previous analysis, we ignored the fact that our responsivity matrix, D, can be sampled with higher resolution than our final reconstruction is capable of producing and that this extra information is lost when we restrict D to a size of m × m. We can re-cast our system of equations to be an underdetermined set by noting that D is in fact of size m × n, where 𝑛 > 𝑚

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

corresponds to the resolution with which we sample the responsivity curves that compose D. Not only can we measure D at a higher resolution, but we can reasonably assume on physical grounds that the responsivity curves are smooth and then interpolate between points to further increase the effective resolution of D. In a standard least-squares regression, such as the one above, this extra information is wasted since it must be discarded when performing the fitting to avoid the accumulation of errors.31 However, following the work of Oliver et al.,25 an adjustment to our reconstruction algorithm allows us to make use of this additional data to substantially increase our reconstructed spectral accuracy. By incorporating a few assumptions about the nature of our input signal, we are able to use this additional information to improve the accuracy of our spectral reconstruction. The improved reconstruction method is based on a technique from compressive sensing for solving underdetermined problems by incorporating the L1-norm, ‖𝑧⃑‖1 ≝ |𝑧1 | + |𝑧2 | + ⋯, as opposed to only L2 (or least-squares) norm minimization in Equation 𝑚𝑖𝑛 || 𝐃𝑥⃑ − 𝑦⃑||22 . 𝑥𝑖 ≥0

(2). By incorporating the L1-norm in our numerical algorithm we can exploit the extra information contained in D rather than discarding it without making our system unsolvable. The inclusion of the L1-norm is analogous to the following pair of assumptions: 1) that there exists a set of basis curves that can accurately reproduce our input spectrum and 2) that in this new basis of curves only a relative few coefficients will be nonzero when describing our input spectrum. Together, these assumptions are referred to as the assumption of sparsity (meaning most elements are zero) in a domain (for a rigorous mathematical formulation see Ref. 32). As a point

ACS Paragon Plus Environment

Page 11 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

of reference, for the spectrometer presented here the number on nonzero coefficients was ~3 out of several hundred, though the technique would still work even for more complicated input signals. In essence, these conditions are a statement about the nature of a signal being measured and how it can be represented. As an example, consider a sine wave of frequency ω and unit amplitude. According to the Shannon-Nyquist theorem, we need, at minimum, 2ω samples in order to uniquely determine this signal.29 With fewer samples than this we run the risk of aliasing, or having different sine waves appear the same due to our sampling rate. An example of how this could happen is presented in Figure S2. Alternatively, we can express this signal in the domain consisting of sine waves with varying frequency (ω, 2ω, 3ω…). In this basis, our reconstruction only needs a single point corresponding to a sine wave of frequency ω, while all other values will be 0. Thus, our signal can be described using a sparse representation in this new basis, since it will only have a single nonzero term. This is the central concept of compressed sensing: carefully selecting the domain in which an input signal is represented can allow for the amount of data needed to describe it to be dramatically reduced, or compressed. We note however that the approach taken in this work has several differences from “traditional” compressed sensing, as discussed in the Supplemental Information. In this work our input is a Gaussian spectrum from a laser source with a known width. If we choose our basis curves to be a set of Gaussian curves, with approximately the correct width, spanning the mid-IR, there is (theoretically) only a single coefficient needed to uniquely determine our input signal (i.e. the coefficient corresponding to the center wavelength of our laser). Now, instead of searching for one coefficient value per wavelength in our reconstruction, we are only looking for the magnitude of a single coefficient and are assuming that the rest are zero (or

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 23

negligible). This dramatically reduces the solution space possible in our system of equations and turns our previously intractable problem into one that is eminently solvable numerically. In general, the steps for performing this type of reconstruction are as follows: first, define a basis set of curves that can be summed with varying magnitude to form any potential input signal of interest – this is called our dictionary, 𝚿 (see Figure S3). Then, assume that our input signal, 𝑠⃑, is sparse in the basis defined by the dictionary. Now, we are able to apply a second condition on the solution to our minimization problem – that the solution should be sparse in this new basis. Without this additional constraint there would be no way to select a single solution from the many possible solutions to such an underdetermined problem. Mathematically, these steps correspond to: 1) defining an n × n dictionary matrix, 𝚿, consisting of a set of basis curves; 2) expressing our reconstructed input signal 𝑥⃑, as a new vector, 𝑠⃑, of length n multiplied by this dictionary set, 𝚿𝑠⃑ = 𝑥⃑ ; 3) retaining the high-resolution information about the response of each detector in D, allowing it to be of size m × n; and 4) incorporating the new L1 constraint into our fitting algorithm, 𝑚𝑖𝑛 || 𝐃𝚿𝑠⃑ − 𝑦⃑||22 + 𝜆‖𝑠⃑‖1 . 𝑠𝑖 ≥0

(3)

Here we have introduced λ as a thresholding parameter to select for sparsity in the final solution (it will not be used to represent wavelength at any later point in the text). That is, we are looking for the vector 𝑠⃑ that has the fewest number of nonzero elements possible, while still minimizing our original equation and setting values smaller than λ to 0.

ACS Paragon Plus Environment

Page 13 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

This type of problem, commonly referred to as basis pursuit denoising (BPDN), has been well studied in the convex optimization literature,33 and is known to be effective at reproducing input signals from a relatively small number of measurements in a wide range of applications.34,35 (MATLAB code for numerically solving Equation 𝑚𝑖𝑛 || 𝐃𝚿𝑠⃑ − 𝑦⃑||22 + 𝜆‖𝑠⃑‖1 . 𝑠𝑖 ≥0

(3) is provided in the Supplemental Information). Using this more sophisticated fitting algorithm, the achievable spectral resolution is doubled, and we are able to resolve peaks that were indistinguishable using an L2-norm reconstruction. As can be seen in Figure 3(e) (solid lines), the L1-norm based spectral reconstruction is able to accurately reconstruct two input signals separated by only 40 nm using the same experimental data as the L2-norm reconstruction. This is well below the Shannon-Nyquist limit imposed by the limited number of data points and is a clear demonstration of the achievable improvements when using more advanced reconstructive methods and some general assumptions about the nature of the input signal. We stress that this type of analysis has been demonstrated to be both robust against noise36 and well-suited for a wide variety of input signals,37 not just those that are explicitly members of the dictionary matrix as in this work. In any number of signalacquisition problems the use of such a reconstruction algorithm has the potential for improvements in signal reconstruction accuracy.

Molecular detection To demonstrate the viability of our microscale IR spectrometer for molecular detection, we placed a quartz window coated with a thin layer (~2-3 µm) of poly(methyl methacrylate) (PMMA) in the optical excitation path of the spectrometer. Similar to the spectral measurements above, we

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

directed a fixed wavelength (centered at 3387 nm to overlap the CH2 and CH3 stretching region) from a tunable mid-IR laser onto each grating, sequentially, and measured the response with and without the PMMA layer for 60 seconds each (full details in Methods). After reconstructing the resulting spectra, we were then able to detect the change in optical intensity due to the vibrational absorption of the PMMA, as compared against an uncoated quartz window, Figure 4(a,b) (solidsymbol lines). As can be seen, there is a clear drop in measured intensity due to the presence of the thin molecular layer. To verify that the measured change in intensity we observed was due to

Figure 4: Molecular detection using an integrated optoelectronic spectrometer. (a) Experimental (solidsymbol) and simulated (dot-dashed) spectral reconstruction with and without PMMA in the beam path (molecular structure inset). (b) Difference between blank and PMMA-coated quartz windows and measured absorption of PMMA (orange, right axis). (c) Simulated spectral reconstruction (black) using 20 (i), 50 (ii), 100 (iii), 150 (iv), and 500 (v) gratings (spectra offset for clarity) along with the simulated input signal (red dashed line, (vi)). Experimental absorption data (orange) from PMMA on a quartz window (right axis).

absorption by the molecules in the beam path and not some other effect, we also performed simulations of the spectral reconstruction process (dash-dot lines) accounting for the separately measured absorbance of the PMMA (see Figure S4). We hypothesize that the discrepancy seen between the simulated and experimental curves (Figure 4(b)) is due to experimental noise

ACS Paragon Plus Environment

Page 14 of 23

Page 15 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

producing a larger shift in the measured spectrum than was theoretically predicted. However, the peak heights match exceptionally well, indicating that the reduction in amplitude is due to the presence of the PMMA layer. Briefly, these simulations were performed as follows (see Supplemental Information for sample code). The predicted output of the laser at a given wavelength was modulated by either the measured absorbance due to a blank quartz window or a window with a layer of PMMA. These simulated spectra were then used as inputs to the measured responsivity curves to generate a series of ideal current values for the gratings. These idealized currents were then fed into the spectral reconstruction algorithm to generate simulated reconstructed spectra. Finally, these simulated reconstructions were scaled by a constant factor to correlate their intensity with that of the experimental reconstructions. The relative peak heights match very well between the simulated and experimental spectra, demonstrating that the spectrometer is accurately measuring the reduction in light intensity due to the molecular absorption. This is particularly evident in the difference between the reconstructed spectra with and without the molecular layer, Figure 4(b), where we see a prominent feature that lines up very well with the molecular absorption of PMMA (at ~3350 nm). Based on the observed noise in the measurements, we estimate that the minimum detectable layer would be a factor of ~4 smaller than what was measured here, enabling the detection of a molecular layer approximately 0.6 µm thick without increasing the integration time of the measurements. To demonstrate the achievable resolution with this type of spectrometer, we performed simulations with increasing numbers of gratings until we were able to completely reproduce the vibrational structure of PMMA, Figure 4(c), with a blackbody source. These simulations were performed by generating random responsivity curves, using those curves to generate current

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

values, and then finally using those currents to do spectral reconstruction. As can be seen, by increasing the number of gratings from 20 to 500 we are able to fully capture the vibrational fingerprint of PMMA.

Discussion In this work we have demonstrated a fully integrated spectrometer capable of reproducing spectra with resolution below the Shannon-Nyquist limit using algorithms adopted from the convex optimization literature. By using such algorithms many of the miniature spectrometers in the literature could substantially increase their achievable resolution. Additionally, we showed that this spectrometer can be used directly for molecular detection without any modification. Finally, we performed simulations to demonstrate that an optoelectronic spectrometer of this type is capable of directly measuring the vibrational modes of a sample. Since this spectrometer is CMOS compatible, straightforward to fabricate, inexpensive, and operates at room temperature, we believe that fully-integrated spectrometers such as this one provide an approach towards developing a compact “nanophotonic nose” capable of gas, liquid, or solid-phase molecular sensing in a miniature package with a direct electrical readout. This type of spectrometer shows promise in an expanding range of applications where compact, fieldable IR spectrometers are of increasing need, including biosensing, environmental monitoring, and chemical imaging.

Materials and Methods Device fabrication

ACS Paragon Plus Environment

Page 16 of 23

Page 17 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

The detectors were fabricated on boron-doped silicon substrates with resistivity of 0.01-0.02 Ω-cm. PMMA A4 (Microchem) was used as an electron-beam resist, followed by development in 1:3 methyl isobutyl keytone (MIBK)/ isopropyl alcohol (IPA) (Microchem). After development, the oxide layer on the substrate was etched in a 10:1 solution of buffered oxide etchant (NH4:HF, Microchem) followed by a 45 second rinse in deionized water to grow a uniform thin oxide layer. The sample was then dried using N 2 gas and immediately placed into the evaporation chamber. 300 nm of pure aluminum (99.99%, Ted Pella) was deposited at 5x10 -7 Torr and a rate of ~5 Å/sec as determined by a quartz crystal microbalance. Lift-off was accomplished by soaking the sample in room-temperature acetone overnight, followed by a brief sonication and IPA rinse. A large (~2 mm radius) ground connection was made to the silicon substrate using indium solder after scratching through the oxide layer with a diamond-tipped scribe. Every grating on the chip, as well as the ground connection, was then wire-bonded with 1 mil gold wire to an external chip carrier (Kyocera) to facilitate consistent repeated measurements. IR measurements of molecules FTIR transmission spectra were taken on a Hyperion 80v spectrometer (Bruker) using a silicon carbide globar source; the background for all measurements was an empty beam path. The molecular samples were prepared by spin-coating PMMA A4 on a quartz window (Ted Pella) before heating at 180° C for approximately 5 minutes to cure and dry.

Electro-optical characterization An optical parametric oscillator (OPO, Coherent) pumped by a Ti:Sapphire laser (Coherent) was used to generate mid-IR light for this study. The OPO idler outputs mid-IR light

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

from roughly 2-4 m with a pulse frequency of 80 MHz, a pulse width of ~200fs, and an average power of ~5 mW at the device. The beam was chopped to allow photocurrent signals to be recorded using a lock-in amplifier (Signal Recovery) in “current mode” which applies a small voltage (~0.15 mV) in opposition to the open-circuit potential of the device. The laser was split using a CaF2 50/50 beamsplitter (Thorlabs) so that its power could be monitored in parallel with the current through the devices. The laser was then focused onto the grating portion of each device individually using a reflective objective (Thorlabs, 15x, 0.30 NA). For responsivity curves, the laser scanned its spectral range in 25 nm intervals while both laser power measurements and current values were recorded. For spectral measurements, the idler was set to a target wavelength and then current and laser power values for each grating were recorded at a constant wavelength for 60 seconds at a rate of 1 Hz. For molecular detection measurements, the beam was additionally passed through either a blank quartz window, or a quartz window coated with a target molecule as described above. For the molecular detection measurements, the laser was centered at 3387 nm to overlap with the CH stretching region of the PMMA.

Spectral reconstruction L2-norm spectral reconstructions were performed using straightforward linear programming least-squares minimization. L1-norm spectral reconstructions were calculated using a gradient-descent method with a fixed step size based on the Lipschitz condition to ensure convergence. Pseudocode for both fitting algorithms and all simulations is supplied in the Supplemental Information. Simulated reconstructions for molecular detection were based on the measured IR spectra of bare quartz and PMMA-coated quartz windows. To determine the number

ACS Paragon Plus Environment

Page 18 of 23

Page 19 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

of gratings necessary for direct vibrational spectral readout, random strength and width responsivity curves with peak positions spanning the range 2.2-3.8 µm were used to approximate physical gratings. These responsivity curves were then used to generate idealized currents for a given input spectrum. Finally, these responsivity curves and currents were used to produce spectral reconstructions just as the experimentally measured values were.

Acknowledgements The authors would like to thank Jesse Adams and Vivek Boominathan for their assistance with the compressive sensing algorithm and Sam Gottheim for enlightening conversation. This work was supported by National Science Foundation grant ECCS-1610229 and Defense Threat Reduction Agency grant HDTRA1-16-0042.

Supporting Information – Responsivity and SEM characterization for all gratings. Introductory compressive sensing images. IR spectra for uncoated- and PMMA-coated quartz windows. Picture of experimental setup. Differences between “traditional” compressive sensing and the approach taken here. Example numerical code for compressive sensing reconstruction.

References:

(1)

Bao, J.; Bawendi, M. G. A Colloidal Quantum Dot Spectrometer. Nature 2015, 523 (7558), 67–70.

(2)

Keuleyan, S.; Lhuillier, E.; Brajuskovic, V.; Guyot-Sionnest, P. Mid-Infrared HgTe Colloidal Quantum Dot Photodetectors. Nat. Photonics 2011, 5 (8), 489–493.

(3)

Dao, T. D.; Ishii, S.; Yokoyama, T.; Sawada, T.; Sugavaneshwar, R. P.; Chen, K.; Wada, Y.; Nabatame, T.; Nagao, T. Hole Array Perfect Absorbers for Spectrally Selective

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Midwavelength Infrared Pyroelectric Detectors. ACS Photonics 2016, 3 (7), 1271–1278. (4)

Lin, C.-H.; Liu, C. W. Metal-Insulator-Semiconductor Photodetectors. Sensors Journal, IEEE 2010, 10, 8797–8826.

(5)

Schwarz, B.; Reininger, P.; Ristanić, D.; Detz, H.; Andrews, A. M.; Schrenk, W.; Strasser, G. Monolithically Integrated Mid-Infrared Lab-on-a-Chip Using Plasmonics and Quantum Cascade Structures. Nat. Commun. 2014, 5.

(6)

Yao, Y.; Shankar, R.; Rauter, P.; Song, Y.; Kong, J.; Loncar, M.; Capasso, F. HighResponsivity Mid-Infrared Graphene Detectors with Antenna-Enhanced Photocarrier Generation and Collection. Nano Lett. 2014, 14 (7), 3749–3754.

(7)

Tittl, A.; Leitis, A.; Liu, M.; Yesilkoy, F.; Choi, D.-Y.; Neshev, D. N.; Kivshar, Y. S.; Altug, H. Imaging-Based Molecular Barcoding with Pixelated Dielectric Metasurfaces. Science (80-. ). 2018, 360 (6393), 1105–1109.

(8)

Craig, B.; Shrestha, V. R.; Meng, J.; Cadusch, J. J.; Crozier, K. B. Experimental Demonstration of Infrared Spectral Reconstruction Using Plasmonic Metasurfaces. Opt. Lett. 2018, 43 (18), 4481.

(9)

Huang, E.; Ma, Q.; Liu, Z. Etalon Array Reconstructive Spectrometry. Sci. Rep. 2017, 7 (1), 40693.

(10) Kuyken, B.; Ideguchi, T.; Holzner, S.; Yan, M.; Hänsch, T. W.; Van Campenhout, J.; Verheyen, P.; Coen, S.; Leo, F.; Baets, R.; et al. An Octave-Spanning Mid-Infrared Frequency Comb Generated in a Silicon Nanophotonic Wire Waveguide. Nat. Commun. 2015, 6, 1–6. (11) Zobenica, Ž.; Van Der Heijden, R. W.; Petruzzella, M.; Pagliano, F.; Leijssen, R.; Xia, T.; Midolo, L.; Cotrufo, M.; Cho, Y. J.; Van Otten, F. W. M.; et al. Integrated Nano-OptoElectro-Mechanical Sensor for Spectrometry and Nanometrology. Nat. Commun. 2017, 8 (1), 1–8. (12) Yu, Z.; Liu, L.; Wang, K.; Wang, J.; Li, J.; Zhao, L. Ultraviolet Wavelength Identification Using Energy Distribution of Hot Electrons. ACS Omega 2017, 2 (7), 3710–3715. (13) Yang, T.; Peng, J.; Li, X.; Shen, X.; Zhou, X.; Huang, X.; Huang, W.; Ho, H. Compact Broadband Spectrometer Based on Upconversion and Downconversion Luminescence. Opt. Lett. 2017, 42 (21), 4375–4378. (14) Hummon, M. T.; Kang, S.; Bopp, D. G.; Li, Q.; Westly, D. A.; Kim, S.; Fredrick, C.; Diddams, S. A.; Srinivasan, K.; Aksyuk, V.; et al. Photonic Chip for Laser Stabilization to an Atomic Vapor with 10 −11 Instability. Optica 2018, 5 (4), 443. (15) Alavirad, M.; Mousavi, S. S.; Roy, L.; Berini, P. Schottky-Contact Plasmonic Dipole Rectenna Concept for Biosensing. Opt. Express 2013, 21 (4), 4328. (16) Rosenqvist, Å.; Milne, A.; Lucas, R.; Imhoff, M.; Dobson, C. A Review of Remote Sensing Technology in Support of the Kyoto Protocol. Environ. Sci. Policy 2003, 6 (5), 441–455.

ACS Paragon Plus Environment

Page 20 of 23

Page 21 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

(17) Schneevoigt, N. J.; van der Linden, S.; Thamm, H. P.; Schrott, L. Detecting Alpine Landforms from Remotely Sensed Imagery. A Pilot Study in the Bavarian Alps. Geomorphology 2008, 93 (1–2), 104–119. (18) Börjesson, T.; Eklöv, T.; Jonsson, A.; Sundgren, H.; Schnürer, J. Electronic Nose for Odor Classification of Grains. Cereal Chem. 1996, 73 (4), 457–461. (19) Shakoor, A.; Cheah, B. C.; Hao, D.; Al-Rawhani, M.; Nagy, B.; Grant, J.; Dale, C.; Keegan, N.; McNeil, C.; Cumming, D. R. S. Plasmonic Sensor Monolithically Integrated with a CMOS Photodiode. ACS Photonics 2016, 3 (10), 1926–1933. (20) Fedyanin, D. Y.; Stebunov, Y. V. All-Nanophotonic NEMS Biosensor on a Chip. Sci. Rep. 2015, 5, 1–12. (21) Zhou, X.; Lee, S.; Xu, Z.; Yoon, J. Recent Progress on the Development of Chemosensors for Gases. Chem. Rev. 2015, 115 (15), 7944–8000. (22) Spinelle, L.; Gerboles, M.; Kok, G.; Persijn, S.; Sauerwald, T. Review of Portable and LowCost Sensors for the Ambient Air Monitoring of Benzene and Other Volatile Organic Compounds. Sensors (Switzerland) 2017, 17 (7). (23) Nagle, H. T.; Gutierrez-Osuna, R.; Schiffman, S. S. The How and Why of Electronic Noses. Spectrum, IEEE 1998, 35 (9), 22–31. (24) Zheng, B.; Zhao, H.; Cerjan, B.; Yazdi, S.; Ringe, E.; Nordlander, P.; Halas, N. J. A RoomTemperature Mid-Infrared Photodetector for on-Chip Molecular Vibrational Spectroscopy. Appl. Phys. Lett. 2018, 113 (10), 101105. (25) Oliver, J.; Lee, W.; Park, S.; Lee, H.-N. Improving Resolution of Miniature Spectrometers by Exploiting Sparse Nature of Signals. Opt. Express 2012, 20 (3), 2613. (26) Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: New York, 2004. (27) Starling, D. J.; Storer, I.; Howland, G. A. Compressive Sensing Spectroscopy with a Single Pixel Camera. Appl. Opt. 2016, 55 (19), 5198. (28) Starling, D. J.; Ranalli, J. Compressive Sensing for Spatial and Spectral Flame Diagnostics. Sci. Rep. 2018, 8 (1), 2556. (29) Shannon, C. E. Communication in the Presence of Noise. Proc. IRE 1949, 37 (1), 10–21. (30) Chen, D.; Plemmons, R. J. Nonnegativity Constraints in Numerical Analysis. In The Birth of Numerical Analysis; WORLD SCIENTIFIC, 2009; pp 109–139. (31) Chang, C.-C.; Lee, H.-N. On the Estimation of Target Spectrum for Filter-Array Based Spectrometers. Opt. Express 2008, 16 (2), 1056–1061. (32) Duarte, M. F.; Eldar, Y. C. Structured Compressed Sensing: From Theory to Applications. IEEE Trans. Signal Process. 2011, 59 (9), 4053–4085.

ACS Paragon Plus Environment

ACS Photonics 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(33) Donoho, D. L. For Most Large Underdetermined Systems of Equations, the Minimal l1Norm near-Solution Approximates the Sparsest near-Solution. Commun. Pure Appl. Math. 2006, 59 (7), 907–934. (34) Donoho, D. L. Compressed Sensing. IEEE Trans. Inf. Theory 2006, 52 (4), 1289–1306. (35) Baraniuk, R. G. Compressive Sensing [Lecture Notes]. IEEE Signal Process. Mag. 2007, 24 (July), 118–121. (36) Figueiredo, M.; Nowak, R.; Wright, S. J. Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems, IEEE J. Sel. Top. Signa 2007, vol (4), 1pp586-597. (37) Zhang, S.; Dong, Y.; Fu, H.; Huang, S. L.; Zhang, L. A Spectral Reconstruction Algorithm of Miniature Spectrometer Based on Sparse Optimization and Dictionary Learning. Sensors 2018, 18 (2).

CONTRIBUTIONS Author contributions: B.C. and N.J.H. designed research; B.C. performed research and analysis; B.C. and N.J.H. wrote the paper. COMPETEING INTERESTS The authors declare no competing financial interests.

ACS Paragon Plus Environment

Page 22 of 23

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Photonics

Towards a nanophotonic nose: a compressive sensing-enhanced, optoelectronic mid-infrared spectrometer

ACS Paragon Plus Environment