Dual-comb photoacoustic spectroscopy

Dual-comb photoacoustic spectroscopy

Play all audios:

Loading...

ABSTRACT Spectrally resolved photoacoustic imaging is promising for label-free imaging in optically scattering materials. However, this technique often requires acquisition of a separate


image at each wavelength of interest. This reduces imaging speeds and causes errors if the sample changes in time between images acquired at different wavelengths. We demonstrate a solution


to this problem by using dual-comb spectroscopy for photoacoustic measurements. This approach enables a photoacoustic measurement at thousands of wavelengths simultaneously. In this


technique, two optical-frequency combs are interfered on a sample and the resulting pressure wave is measured with an ultrasound transducer. This acoustic signal is processed in the


frequency-domain to obtain an optical absorption spectrum. For a proof-of-concept demonstration, we measure photoacoustic signals from polymer films. The absorption spectra obtained from


these measurements agree with those measured using a spectrophotometer. Improving the signal-to-noise ratio of the dual-comb photoacoustic spectrometer could enable high-speed spectrally


resolved photoacoustic imaging. SIMILAR CONTENT BEING VIEWED BY OTHERS MID-INFRARED CROSS-COMB SPECTROSCOPY Article Open access 24 February 2023 PHOTO-ACOUSTIC DUAL-FREQUENCY COMB


SPECTROSCOPY Article Open access 20 August 2020 DUAL-COMB PHOTOTHERMAL SPECTROSCOPY Article Open access 21 April 2022 INTRODUCTION In vivo and in vitro medical imaging and spectroscopy face


the fundamental challenge of strong optical scattering in biological tissues. This challenge has led researchers to use the photoacoustic (PA) effect, where localized light absorption


generates acoustic waves that can be detected from deep within tissue1,2,3,4. PA measurements have an advantage over pure optical methods because light only has to travel through the sample


to the absorber to create an acoustic wave, but the light does not need to travel back out of the sample to a detector. Other advantages of PA measurements are that scattered optical


illumination can still induce a PA signal, and the generated acoustic waves are only weakly scattered by tissue5. In addition, by measuring the PA response versus excitation wavelength, it


is possible to identify and quantify substances based on their unique absorption features. For example, researchers have used this method, called PA spectroscopy, to quantify glucose levels


in vivo6,7,8. A more advanced technique for medical diagnostics and functional imaging, called multispectral optoacoustic tomography (MSOT) or photoacoustic tomography (PAT), has also been


demonstrated. This technique provides images of multiple exogenous and endogenous species at multiple wavelengths and uses spectral unmixing algorithms to determine the spatial distribution


of each species2,9,10. This technique has been used to differentiate adipose tissue from cholesterol11,12, measure intramuscular fat distributions13, map tumor size and shape14,15, identify


white matter loss in spinal cord injuries16, quantify blood oxygen saturation17,18, and even to construct a label-free molecular map of a mouse fibroblast cell19. One challenge of PA


spectroscopy, or PAS, measurements is that they require a broadband optical excitation source with a spectrally resolved detection or illumination method. Traditionally, PA researchers have


used optical parametric oscillators or other tunable lasers, which require sequential image acquisition over narrow optical-frequency bands. Although it is possible in simple ratiometric PAS


to use only a couple of optical-frequency bands, measuring in more bands allows the identification of additional species and improves detection sensitivity20. However, sequentially scanning


the optical frequency can be time consuming and can lead to potential issues with sample drift or damage between images at different optical frequencies9,21. Broadband thermal sources using


a Fourier transform spectrometer (FTS)22,23,24,25 for optical-frequency resolution have been used to overcome sequential scanning, but the low optical power and slow scanning delay arm have


limited the application of these systems. Alternatively, broadband supercontinuum coherent laser sources may provide higher power than a thermal source, but still require either sequential


frequency scanning or physical scanning of an FTS delay arm and can suffer from high laser intensity noise. Optical frequency combs26,27,28 are an alternative light source for broadband PAS


and multispectral PAT because they simultaneously generate thousands of discrete optical-frequency bands. These frequency bands are called comb teeth because, like the teeth of a comb, they


are evenly spaced and very narrow, with absolute linewidths below 120 kHz29. The comb teeth of a frequency comb are temporally and spatially coherent and propagate in a single spatial mode


with high brightness. Recent work has shown that frequency combs are a promising light source for PAS. Researchers have demonstrated highly sensitive and precise measurements of methane


absorption spectra with frequency-comb-based PAS30,31. These demonstrations relied on a FTS for spectral selectivity and thus required physical scanning of the optical pathlength. It is


possible, however, to eliminate the need for the FTS altogether by using a second comb like a local oscillator. By exploiting the coherence of optical-frequency combs, one can capture an


entire absorption spectrum on a single photodetector without wavelength or delay arm scanning, as is done in dual-comb spectroscopy (DCS)32. DCS provides broad spectral coverage, high


sensitivity, and high spectral resolution for identification of multiple chemical species33,34,35,36,37. It also can provide rapid spectral measurements for monitoring dynamic biological


processes38, as well as for coherent Raman imaging39,40. These advantages should be applicable to PAS using a technique analogous to DCS. This method of dual-comb PAS (DCPAS), which was


suggested in ref. 31, can be used across broad spectral windows without sequential wavelength scanning and can be implemented with a compact, stable, and robust dual-comb spectrometer. In


this article, we report an experimental demonstration of DCPAS. We use two self-referenced fiber-laser frequency combs to measure the PA spectra of thick polymer films over an ≈15 THz


spectral bandwidth centered at 173 THz. With coherent averaging we obtain spectrally resolved PA responses of polydimethylsiloxane (PDMS) and paraffin samples immersed in water. These


experiments reveal that dual-comb optical excitation can generate PA signals and that these signals can be used to detect and identify organic materials. RESULTS DCPAS CONCEPT Figure 1


describes the concept of DCPAS. The DCPAS system is based on laser frequency combs whose frequency-domain output is tens of thousands of evenly spaced narrow lines or comb “teeth”. These


teeth are spaced at the frequency-comb’s pulse repetition rate, _f_rep, ((i) in Fig. 1)27. In the dual-comb approach, we use two such frequency combs where the repetition rate of comb 2,


_f_rep,2, is slightly larger than the repetition rate of comb 1, _f_rep,1, such that _f_rep,2 − _f_rep,1 = Δ_f_rep. As a result, the frequencies of subsequent pairs of comb teeth are offset


by increasing increments of Δ_f_rep, as illustrated in Fig. 1. The optical intensity of the combined combs is given by \(I_0 = \left| {E_1 + E_2} \right|^2\), where _E_1 is the electric


field of the light from comb 1 and _E_2 is the electric field of the light from comb 2. This intensity contains a multi-heterodyne term corresponding to the beat notes of each, increasingly


offset, pair of comb teeth. Assuming the zeroth tooth of comb 1, at frequency _v_1,0, overlaps with the zeroth tooth of comb 2, at frequency _v_2,0 = _v_1,0, then the beat notes occur at


radio frequencies, _f__i_, given by $$\begin{array}{*{20}{c}} {f_1} & = & {\nu _{2,1} - \nu _{1,1}} & = & {\left[ {\nu _{2,0} + 1 \times f_{{\mathrm{rep}},2}} \right] - [\nu


_{1,0} + 1 \times f_{{\mathrm{rep}},1}]} & = & {1 \times \Delta f_{{\mathrm{rep}}}} \\ {f_2} & = & {\nu _{2,2} - \nu _{1,2}} & = & {\left[ {\nu _{2,0} + 2 \times


f_{{\mathrm{rep}},2}} \right] - \left[ {\nu _{1,0} + 2 \times f_{{\mathrm{rep}},1}} \right]} & = & {2 \times \Delta f_{{\mathrm{rep}}}} \\ {f_3} & = & {\nu _{2,3} - \nu


_{1,3}} & = & {\left[ {\nu _{2,0} + 3 \times f_{{\mathrm{rep}},2}} \right] - \left[ {\nu _{1,0} + 3 \times f_{{\mathrm{rep}},1}} \right]} & = & {3 \times \Delta


f_{{\mathrm{rep}}}} \\ \vdots & = & \vdots & = & \vdots & = & \vdots \\ {f_n} & = & {\nu _{2,n} - \nu _{1,n}} & = & {\left[ {\nu _{2,0} + n \times


f_{{\mathrm{rep}},2}} \right] - \left[ {\nu _{1,0} + n \times f_{{\mathrm{rep}},1}} \right]} & = & {n \times \Delta f_{{\mathrm{rep}}}} \end{array}$$ (1) As the optical frequencies


_v__i,j_ are typically 10 to 12 orders of magnitude larger than Δ_f_rep, _v_1,_i_ is approximately the same as _v_2,_i_, and it is quite accurate to say that light at the optical-frequency


_v__i_ = _v_1,_i_ ≈ _v_2,_i_ is amplitude modulated at the radio frequency _f__i_. When the sample absorbs light from a pair of comb teeth at optical-frequency _v__i_, it produces a PA


pressure wave with an amplitude, _p_, proportional to the amount of light absorbed by the sample ((iii) in Fig. 1)41,42,43. This proportionality is given by \(p \propto I_0\left( {\nu _i}


\right) \times \left[ {1 - \exp \left( { - \mu _{\mathrm{A}}\left( {\nu _i} \right) \times L} \right)} \right] \times \beta\), where _I_0(_v__i_) is the intensity of the dual-comb excitation


at optical-frequency _v__i_, _μ_A(_v__i_) is the sample’s base-e absorption coefficient at optical-frequency _v__i_, _L_ is the sample thickness that contributes to the PA signal, and _β_


is a proportionality constant that depends on sample material and geometry. The amplitude of this pressure wave is modulated at the multi-heterodyne frequency _f__i_, which can be several


orders of magnitude less than the comb repetition frequencies, _f_rep. When the sample absorbs light from many pairs of comb teeth, at optical frequencies _v_0,_v_1._v_2,…,_v__n_, it


produces a PA pressure wave with frequency components centered at the corresponding modulation frequencies _f_0,_f_1,_f_2,…,_f__n_. As there is a one-to-one correlation between the frequency


_v__i_ of the absorbed light and the resulting acoustic modulation frequency, _f__i_, one can infer the sample’s optical absorption spectrum by analyzing the frequency-domain PA signal. As


every comb tooth is simultaneously incident on the sample and every PA frequency is simultaneously recorded, this technique allows parallel acquisition of thousands of spectral elements and


multiplexed measurement of the sample’s absorption spectrum. It is also useful to consider the time-domain picture. In the time-domain, the DCPAS signal is a series of interferograms, just


as in DCS. Each interferogram reflects the time-dependent optical intensity caused by the interference between successive pulses from the two offset frequency combs. In DCPAS, this


time-dependent optical intensity leads to a time-dependent amplitude modulation of the PA pressure wave, which we detect via an ultrasound transducer. The pressure waves are detected by the


transducer to generate a corresponding time-domain voltage interferogram. We then convert these interferograms to frequency-domain spectra by a Fourier transform to generate the


multi-heterodyne acoustic signal discussed above ((iv) in Fig. 1). This acoustic signal is normalized by the optical excitation spectrum recorded by a photodetector ((v) in Fig. 1). Then the


acoustic frequencies are mapped to optical frequencies using the known parameters of the combs32. This yields the sample’s optical absorption spectrum ((vi) in Fig. 1). One advantage of


DCPAS is that the bandwidth and center frequency of the generated acoustic signal, corresponding to the amplitude modulation on the optical signal, can be adjusted independently of the


optical bandwidth and center frequency. This adjustment is made in a controlled fashion by tuning the frequency locking of the two combs to adjust _v_2,0, _v_1,0, and Δ_f_rep, and therefore,


the multi-heterodyne frequencies _f__i_. This allows any optical spectrum to generate a PA signal that can be detected by an ultrasound transducer despite the transducer’s limited


bandwidth. In this demonstration, we use lasers with ≈160 MHz repetition rates and set the difference in repetition rates (tooth spacings), Δ_f_rep = _f_rep,2−_f_rep,1, to be 66.81 Hz. This


allows us to map the roughly 15 THz optical bandwidth (full width at –10 dB) into an acoustic bandwidth of (15 THz) × (Δ_f_rep/_f_rep) ≈6.3 MHz32 while centering this modulation signal at


the peak of the transducer responsivity, 7.5 MHz. In order to maintain this mapping, as well as the optical coherence required for long-term averaging of the signals, both combs are fully


stabilized in the manner described in ref. 44. EXPERIMENTAL SETUP Figure 2 shows the experimental setup for DCPAS. We use two self-referenced erbium fiber frequency combs at 160 MHz


repetition rate. The output of each comb is amplified and filtered to generate light in a 15-THz-wide band centered at ≈173 THz, selected to coincide with low water absorption and the first


overtone absorption of C–H bonds. This wavelength region has previously been used to image arterial plaques (fatty tissue) in the presence of blood, to distinguish protein and fat, and to


image white matter in a spinal cord4,13. A total of about 24 mW of comb light is incident on the samples. The PA signals generated by the samples then propagate through a water bath to a


polytetrafluoride ultrasound transducer with a frequency response that is peaked at 7.5 MHz with an ≈7.5 MHz bandwidth (see Supplementary Fig. 3). The PA and photodetected signals are


coherently summed in a data acquisition system previously used for mid-infrared DCS34,45. The methods section provides more detail of the experimental setup. DCPAS EXPERIMENTAL RESULTS


Figure 3 shows the DCPAS measurement results from a sample consisting of vertically aligned carbon nanotubes (VACNTs) embedded in PDMS (polydimethylsiloxane). This sample was chosen because


VACNTs have strong optical absorption at ≈170 THz46 and therefore generate a strong PA response when illuminated. This provides a clear demonstration that dual-comb optical excitation (Fig. 


3a) induces a PA response (Fig. 3b). Furthermore, the shape of the PA response reveals that the PA pressure wave is an interferogram that corresponds to the optical interferogram incident on


the sample. The correlation between the optical excitation and the PA response can be probed further by examining the frequency-domain representation of these signals (Fig. 3c). To obtain


the frequency-domain signals, the time-domain signals are apodized with a 4.8 μs window, zero-padded outside of the apodization window, and Fourier transformed. The Fourier transform yields


the spectrally resolved optical excitation and PA response at acoustic frequencies corresponding to the beat frequencies between the teeth of the two frequency combs (top axis of Fig. 3c).


To convert these acoustic or ultrasound frequencies, _f_US, to the optical frequencies, _v_optical, of the comb teeth that generated the signal, we linearly scale the frequency according to:


$$\nu _{{\mathrm{optical}}} = \nu _0 + \frac{{f_{{\mathrm{rep}},1} + f_{{\mathrm{rep}},2}}}{{2 \times \left( {f_{{\mathrm{rep}},2} - f_{{\mathrm{rep}},1}} \right)}} \times


f_{{\mathrm{US}}},$$ (2) where _f_rep,1 is the repetition rate of comb 1, _f_rep,2 is the repetition rate of comb 2, and _v_0 is a constant offset calculated from the comb repetition rates


and carrier-envelope offset frequencies. With the 4.8 μs apodization window used here, the spectral resolution of the measurement is ≈500 GHz or ≈5 nm (due to the zero-padding, the spectrum


is effectively smoothed to 500 GHz). The spectrally resolved PA response shows how a sample’s absorption coefficient depends on optical frequency. At frequencies where the sample absorbs


strongly, there is a strong PA response; whereas, at frequencies where the sample does not absorb, there is no PA response. Therefore, if a sample absorbs uniformly at all frequencies


contained in the optical excitation, the spectrally resolved PA response should have the same shape as the optical excitation spectrum. The spectrally resolved PA response of the VACNT


sample (Fig. 3c, red line) has a frequency dependence that is similar to the optical excitation spectrum (Fig. 3c, black line). This similarity shows that the PA response is dependent on the


optical excitation spectrum. The relatively weaker amplitude of the PA response for optical frequencies \(\nu _{{\mathrm{optical}}} \, \lesssim \, 166\,{\mathrm{THz}}\) and \(\nu


_{{\mathrm{optical}}} \, \gtrsim \, 176\,{\mathrm{THz}}\) is due to a roll-off in the transducer’s responsivity away from its peak responsivity at _f_US = 7.5 MHz (which corresponds to


_v_optical = 171.2 THz). There are small (around 20%) variations in the PA response relative to the optical excitation spectrum for optical frequencies between 166 THz and 176 THz whose


origins are unknown, but could be caused by sample inhomogeneities. Despite these variations, the VACNT target was useful for characterizing the DCPAS system and for demonstrating that


dual-comb excitation can induce a PA interferogram. However, the VACNT’s absorption spectrum is not representative of a biological material because it lacks a characteristic C-H overtone in


the 175-THz band. To observe the C-H overtone absorption, we measured two different polymer thin films: polydimethylsiloxane (PDMS) and paraffin, both of which are much more weakly absorbing


in this spectral region than the VACNT sample. The PA response of PDMS and paraffin samples are shown in Fig. 4a, b, respectively, for several different averaging times. As the polymers


absorb more weakly than the VACNTs and are less efficient at generating acoustic waves, they produce a DCPAS signal that is 50× and 100× weaker than the VACNTs, as can be seen from comparing


Fig. 3c to Fig. 4a, b, respectively. The spectral features of the PDMS response are somewhat distinguishable over the noise with as little as 4 min of signal averaging and are easily


distinguishable with 20 min of averaging. On the other hand, the spectral features of the paraffin response are readily apparent even at 4 min of averaging due to the approximately two times


stronger signal from the paraffin film. Figure 4c, d compare these PA responses to the samples’ optical absorption spectra measured with a spectrophotometer. In order to make this


comparison, we need to account for several effects inherent to the DCPAS measurement. As mentioned in the introduction, the amplitude of the PA pressure wave is proportional to the sample’s


optical absorbance, the intensity of the optical excitation at the sample, and other things such as acoustic impedance mismatch and the sample’s acoustic response. The measured PA response


is proportional to the pressure wave’s amplitude multiplied by the transfer function of the transducer47,48. For weakly absorbing samples like the paraffin and PDMS in these experiments, the


exponential in Beer’s law can be linearized, and we can write the measured PA response as $$V_{{\mathrm{PA}}}\left( {\nu _{{\mathrm{optical}}}} \right) \approx \beta \times \mu


_{\mathrm{A}}\left( {\nu _{{\mathrm{optical}}}} \right) \times I_0\left( {\nu _{{\mathrm{optical}}}} \right) \times H_{{\mathrm{transducer}}}\left( {\nu _{{\mathrm{optical}}}} \right),$$ (3)


where _V_PA is the measured PA response as a function of optical frequency, _v_optical, _β_ is a proportionality factor to convert from sample absorption to PA pressure, which may be


material and geometry dependent, _μ_A(_v_optical) is the sample’s optical absorption coefficient at optical frequency _v_optical, _I_0(_v_optical) is the intensity of the optical excitation


at optical frequency _v_optical, and _H_transducer (_v_optical) is the ultrasound transducer’s frequency-dependent responsivity after scaling to optical frequencies according to Eq. (2). An


apparent optical-frequency dependence of _β_ could arise from acoustic-frequency-dependent pressure generation or attenuation; however, here we assume that _β_ is constant over the acoustic


frequencies used for the amplitude modulation in these experiments. _H_transducer is obtained from calibration curves provided by the manufacturer (Supplementary Fig. 3), and


_I_0(_v_optical) is measured via the photodetected DCS spectrum. Thus, the normalized PA response of a sample is $$\widehat {V_{{\mathrm{PA}}}}\left( {\nu _{{\mathrm{optical}}}} \right) =


\frac{{V_{{\mathrm{PA}}}\left( {\nu _{{\mathrm{optical}}}} \right)}}{{I_0\left( {\nu _{{\mathrm{optical}}}} \right) \times H_{{\mathrm{transducer}}}\left( {\nu _{{\mathrm{optical}}}}


\right)}} = \beta \times \mu _{\mathrm{A}}\left( {\nu _{{\mathrm{optical}}}} \right).$$ (4) Figure 4c, d show the result of normalizing the polymers’ PA responses according to Eq. (4). For


comparison, these figures also show the optical absorbance spectra, ∝_μ_A(_v_optical), of PDMS and paraffin measured in a transmission-mode spectrophotometer. The normalized PA responses


reflect the absorbance spectra from the spectrophotometer measurements. The DCPAS measurement reveals prominent PDMS absorbance peaks at 171.6 THz, 176.0 THz, 177.4 THz. The location of


these peaks are in good agreement with the spectrophotometer absorbance measurement and previously published results49,50. The normalized PA response of paraffin contains well-defined


absorbance peaks at 169.8 THz and 173.1 THz. Although we could not find literature showing paraffin absorbance spectra in the 165 THz to 180 THz window, the absorbance peaks from the DCPAS


measurement closely match the peaks from our spectrophotometer measurement, and the peak at 173.1 THz agrees well with the location of a reported absorbance peak51. The results shown in Fig.


 4 prove that DCPAS is sensitive to the absorption spectra of polymer samples. However, the signal-to-noise ratio (SNR) of this proof-of-concept DCPAS system is low for weakly absorbing


features such as the C-H overtone bands in these polymer films. Thus, long averaging times are currently required to obtain spectrally resolved responses like those in Fig. 4. Figure 5a


shows the time-domain SNR versus averaging time for both VACNTs and the polymer samples. (The time-domain SNR is defined in the Methods section and Supplementary Note 1.) In all cases the


SNR is proportional to the square root of averaging time, as expected for white noise. The VACNT sample requires 4 orders of magnitude less averaging time than the polymer samples to reach


the same SNR because the VACNTs produce an ≈100× stronger PA signal than the polymers. Figure 5b shows that the SNR is linearly proportional to the power of the optical excitation. This


proportionality indicates that the SNR of the system is limited by detector noise and not subject to saturation effects even for the strongly absorbing VACNT sample. This lack of saturation


suggests that higher SNR measurements of organic films are feasible with a higher-power dual-comb source. DISCUSSION The reported results highlight the potential of DCPAS for label-free


imaging without spectral scanning. As shown in Fig. 5, a high SNR ( >300) can be reached through coherent averaging. However, the averaging times required for high spectral SNR are likely


too long to be useful, especially for medical imaging applications, and it is useful to consider potential routes to reduce the acquisition time. First, the signal strength could be


improved by operating at a wavelength range with larger molecular absorption cross-sections. For example, the PA signal from the VACNT sample has a much higher SNR than that of the PDMS


sample because the VACNTs absorb light much more strongly than the PDMS in this wavelength region. In this proof-of-principle experiment, we focused on the optical window at ≈173 THz because


it was easily accessible with our laser system, contains first overtone absorption by C–H bonds, and has lower water absorption than adjacent frequency bands13. However, mid-infrared


frequency combs34,35,38,52 could measure much stronger fundamental C-H rovibrational transitions. These stronger absorption features (up to several orders of magnitude) could provide a


proportional increase in the PA signal strength and a corresponding quadratic decrease in the required averaging time. However, strong water absorption would limit mid-infrared DCPAS to


<100 μm penetration depths in biological tissue because the 1/e penetration depth is only 32 μm for excitation at _v_optical = 60 THz53. Nevertheless, mid-infrared DCPAS could be useful


for measuring thin samples in the same way that pulsed PAS can measure thin samples at mid-infrared wavelengths8. Moving to visible/near-infrared wavelengths would enable measurements of


oxygenated and deoxygenated hemoglobin. This visible/near-infrared hemoglobin absorption is about ten times stronger than lipid absorption at around 173 THz without any competing absorption


from water4. Additionally, other exogenous agents and dyes with strong absorption features could be measured with visible/near-infrared wavelengths9. The acquisition time could also be


reduced by use of a higher comb repetition frequency. In DCPAS, the fundamental sample point spacing of the measured spectrum is set by the comb repetition rate, _f_rep. In this


demonstration, the spectral sample spacing is _f_rep = 160 MHz, much finer than any absorption feature of interest. For this reason, we apply a 4.8-μs apodization window before the Fourier


transformation. This apodization improves the spectral SNR by a factor of \(\sqrt {500\,{\mathrm{GHz/}}160\,{\mathrm{MHz}}} \approx 56\) and smooths the frequency-domain response to a


spectral resolution of ≈500 GHz. This 500-GHz resolution is acceptable for the envisioned applications because the absorption features in tissue samples are typically 100 s of GHz to 1000 s


of GHz wide. For instance, the full-widths at half-maximum of lipid absorption spectra at 172 THz are ≈6000 GHz, and the absorption peaks of plaque lipids and adipose lipids are separated by


about 1200 GHz12. However, rather than smoothing the frequency-domain response via apodization, it would be better if the fundamental sample spacing, _f_rep, matched the desired resolution.


In other words, there is an SNR penalty because the fundamental spectral sample spacing is not matched to the desired resolution. This SNR penalty can be understood by recognizing that a


wider tooth-spacing, with a fixed total power, results in a greater power per comb tooth. Taking this into account, one finds that the required averaging time should scale linearly with


1/_f_rep for a fixed output power and fixed SNR up until the repetition rate matches the desired resolution54. It is possible to generate frequency combs with 10 GHz repetition rates from a


mode locked laser55,56, quantum cascade laser-based combs38, or the use of electro-optic frequency combs57,58, while microresonator-based frequency combs could provide 500 GHz repetition


rates59,60. A 500 GHz repetition rate would result in a ≈3000-fold reduction in averaging time, or, for a given averaging time a \(\sqrt {3000} = 55 \times\) improvement in SNR. However,


there are challenges to increasing the repetition rate while simultaneously maintaining the average comb output power and spectral coverage. Nonetheless, recent reports of dark-pulse


microresonators and dissipative Kerr soliton combs have shown >200 GHz repetition rate and on-chip powers exceeding 10 mW61,62. With further improvements in efficiency and output


coupling, these comb sources could be well-suited for DCPAS. If we assume a 500-GHz repetition rate, the same excitation power used here, and a tenfold stronger absorbance by illuminating in


the visible or mid-infrared, the acquisition time is reduced by a factor of 300,000 so that a 2-h acquisition time drops to 24 milliseconds. Depending on the application, this could be


further reduced by increasing the incident optical power. As shown in Fig. 5b, a linear increase in optical power leads to a linear improvement in SNR and to a corresponding quadratic


decrease in the required averaging time. In some cases, additional considerations will limit the incident power for a given application. For example maximum permissible exposure limits for


in vivo samples63 and potential thermal damage for in vitro samples must be avoided. In addition to increasing the incident optical power, the 24-millisecond acquisition time could also be


decreased by decreasing the measurement noise. As the PA noise floor is limited by the transducer noise, improved transducer readout electronics or lower-noise transducing elements would


directly improve the SNR43 and lead to a quadratic reduction in acquisition time. In this initial demonstration of DCPAS, we did not choose an optimal transducer, but rather the 7.5 MHz


transducer used is similar to transducers typically employed for conventional PAS experiments. The ability to adjust the multi-heterodyne frequencies of the optical excitation in DCPAS would


allow the use of transducers with center frequencies and bandwidths both ranging from ≈1 MHz to ≈100 MHz. The arguments mentioned above suggest that DCPAS could be useful for applications


such as spectrally resolved medical imaging. In addition to considering these arguments, we also evaluated the potential of DCPAS by directly comparing DCPAS with a PAS system using a


tunable pulsed optical parametric oscillator (OPO). As shown in Supplementary Fig. 1 and Supplementary Fig. 2, the current DCPAS implementation achieves a similar SNR to that obtained with


pulsed PAS with 20 nJ pulse energies after accounting for wavelength multiplexing. Existing fiber-coupled OPO-based PAS systems can provide up to 20× more pulse energy and a correspondingly


higher SNR. However, given the scaling arguments above, the use of higher repetition rate combs could improve DCPAS SNR by a factor of 55×. Therefore, an improved DCPAS implementation would


be competitive with PAS systems that use tunable pulsed lasers while enabling many spectral elements to be acquired simultaneously without the extra time and potential systematics introduced


by spectral scanning. Finally, DCPAS would also benefit from a more accurate normalization procedure to remove systematic wavelength-dependent variations in the PA response. Our


normalization procedure assumes that optical scattering, sample heating, pressure wave generation, and pressure wave propagation are independent of the multi-heterodyne frequencies in the


optical intensity. Any frequency-dependent effects in these processes would skew the inferred absorption spectrum because of the coupling between ultrasonic and optical frequencies


delineated by Eq. (2). One could address this problem by developing a model for the frequency dependence of PA signal generation under dual-comb excitation. This frequency dependence would


be included in the _β_ proportionality factor in Eq. (4). Alternatively, one could use a reference sample with a known (and ideally flat) absorption spectrum to characterize the frequency


response of the DCPAS system23. In applications where it is not necessary to determine quantitative absorption spectra, these issues might be completely circumvented by applying principal


component analysis to detect specific substances in DCPAS measurements8. Another limitation of the normalization procedure used here is that the optical excitation spectrum is only measured


before the light hits the sample. This could be improved by using optical inversion calculations to obtain the optical spectrum as a function of depth. In that case, one could use a


depth-resolved version of Eq. (4) to normalize the PA response64. Overall, the results demonstrated here show that dual-comb excitation can generate PA signals. These PA signals, in turn,


reflect the absorption spectra of sample materials and can be used to detect and identify different organic solids. Furthermore, we propose several strategies to increase DCPAS acquisition


speeds and to develop more accurate normalization procedures. With these improvements, DCPAS could become a useful new technology for high-speed, label-free spectroscopy and possibly


imaging. The authors would like to acknowledge a recent preprint demonstrating DCPAS for detection of gaseous acetylene65. This shows yet another potential application of DCPAS. METHODS


SAMPLES We used three different samples for DCPAS in this work. One sample consists of a ≈1 mm × 2 mm film of vertically aligned carbon nanotubes (VACNTs) embedded in an ≈240 μm thick film


of poly(dimethylsiloxane) (PDMS). The VACNTs are multi-walled CNTs with an estimated tube diameter of ≈10 nm to ≈20 nm based on measurements of similarly prepared samples. The thickness of


the VACNT layer varies from ≈15 μm to ≈40 μm (see Supplementary Fig. 4). The PDMS film is ≈2.5 mm wide and ≈5 mm long, with several millimeter-wide regions that contain no VACNTs (see


Supplementary Fig. 5). One of these regions is used as a second PA sample that consists of only PDMS. The third sample is a piece of paraffin-based paper (Parafilm “M”, ≈120 μm thick) cut


into a ≈5 mm × 5 mm square (the use of trade names is necessary to specify the experimental results and cannot imply endorsement by the National Institute of Standards and Technology.). Each


of these samples is laminated onto the top of a N-BK7 wedge by gently applying pressure with a pair of tweezers. The wedge is sealed against an opening in the bottom of a water bath so that


the sample is immersed in the water (location (v) in Fig. 2) and oriented toward the transducer. The distance from the sample to the transducer is equal to the transducer focal length of


approximately 2 cm. We also measured the optical absorption spectra of the PDMS and paraffin in transmission-mode with a NIR spectrophotometer (Perkin-Elmer Lambda 1050 equipped with a


3-detector module [The use of trade names is necessary to specify the experimental results and cannot imply endorsement by the National Institute of Standards and Technology]) for comparison


with the spectra obtained in the DCPAS system. The same PDMS sample was used for DCPAS and spectrophotometer measurements. The spectrophotometer illuminated a portion of the PDMS sample


that contained no VACNTs; however, this illumination was not in the exact same location as the dual-comb excitation. The paraffin sample used for DCPAS could not be used for


spectrophotometer measurements because it displayed too much optical scattering for transmission measurements. For transmission measurements, we reduced the optical scattering by placing the


paraffin between two glass microscope slides, heating the glass to 70° C, and applying pressure. Spectrophotometer measurements of this sample were made after allowing it to cool to room


temperature. DCPAS SYSTEM DESIGN Our DCPAS system uses two self-referenced fiber frequency combs with repetition rates, _f_rep, of ≈160 MHz, described in more detail in refs. 44 and 66. Each


comb emits ≈40 mW of light centered at 192 THz (1560 nm) with a 3-dB bandwidth of ≈45 nm. This output is filtered with a 1550 nm ± 5 nm bandpass filter and then amplified to ≈400 mW with a


triple-pumped erbium-doped fiber amplifier. After amplification, the light from each comb is spectrally broadened in ≈3 cm of highly nonlinear fiber with a dispersion of ≈2.2 ps nm−1 km−1.


This yields ≈60 mW of light from each comb spanning 158 THz to 187 THz (location i in Fig. 2). This 158 THz to 187 THz light from each comb is combined on a beam-splitter, transmitted


through a 187 THz low-pass filter and a linear polarizer, and coupled into single-mode PM1550 fiber (location ii in Fig. 2) for transport to the PA setup. The light is then launched into


free space by a _f_ = 50.8 mm reflective collimator. About 88% of the light is transmitted through a CaF2 wedge and ≈6% is reflected by both the front and back faces. The back reflection is


coupled into single-mode fiber and detected by an extended InGaAs photodiode (location (iii) in Fig. 2) to record the optical excitation spectrum. The front reflection is coupled into a


single-mode fiber delay line, for ≈100 μs of delay. Then it is detected by a standard InGaAs photodiode to provide a trigger for the FPGA data acquisition system. The ≈25 mW of light


transmitted through the CaF2 wedge is focused by a _f_ = 45 mm achromatic lens to a spot size of 4 μm on the sample to generate the PA signal. Absorption by the sample generates pressure


waves that are detected by a polytetrafluoride ultrasound transducer with a frequency response that is peaked at 7.5 MHz with a ≈7.5 MHz bandwidth (Supplementary Note 2) and a focal length


of 2 cm (Olympus Panametrics A320S [the use of trade names is necessary to specify the experimental results and cannot imply endorsement by the National Institute of Standards and


Technology]). This transducer is mounted on a 3-D translation stage and positioned to maximize the PA signal, which occurs when the transducer focus overlaps the optical focus of the


incident comb light. In addition, the sample is translated vertically to place it at the common focal plane of the comb light and transducer. As there is 2 cm (transducer focal length) of


water between the sample and the transducer and because the absorption coefficient of water is _μ_A > 6.5 cm−1 for the optical frequencies in this experiment, no light reaches the face of


the transducer. We have also verified that when the sample is translated out of the common focus of the illumination and transducer, no signal is measured. After the ultrasound transducer,


the PA signal is amplified by ≈80 dB and filtered with a 15 MHz low-pass filter before it is input to the data acquisition system, as shown in Fig. 2. The optical excitation signal is


filtered with a 50 MHz low-pass filter after the photodetector and before the data acquisition system. The PA signal and the optical excitation signal are digitized on separate channels of


the data acquisition system, and the resulting time-domain signals are interferograms, as shown in Fig. 3a, b. For the SNR calculations used in Fig. 5, the time-domain SNR is defined as the


peak height of the interferogram divided by the standard deviation of the noise outside of the interferogram centerburst (see Supplementary Fig. 2 and Supplementary Note 1). As the SNR is


low for the weakly absorbing polymer samples, multiple interferograms are coherently averaged in real-time as is common in DCS34,45. We implement this signal averaging on a data acquisition


system that consists of a multi-channel 14-bit analog-digital converter and a field-programmable gate array (FPGA). The firmware in the FPGA removes residual phase noise from the


interferograms before co-adding them for long-term averaging34,45. To measure slow phase drifts, the phase correction firmware requires a stronger signal than is present in a single PA


interferogram. We resolve this issue by adding the optical DCS voltage to the DCPAS voltage as a delayed trigger signal as shown in Fig. 2. DATA AVAILABILITY Data and processing scripts are


available on request from the authors. REFERENCES * Beard, P. Biomedical photoacoustic imaging. _Interface Focus_ 1, 602–631 (2011). Article  PubMed  PubMed Central  Google Scholar  * Wang,


L. V. & Hu, S. Photoacoustic tomography: in vivo imaging from organelles to organs. _Science_ 335, 1458–1462 (2012). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Drexler,


W. et al. Optical coherence tomography today: speed, contrast, and multimodality. _J. Biomed. Opt._ 19, 071412 (2014). Article  ADS  PubMed  Google Scholar  * Hui, J. et al. Bond-selective


photoacoustic imaging by converting molecular vibration into acoustic waves. _Photoacoustics_ 4, 11–21 (2016). Article  PubMed  PubMed Central  Google Scholar  * Zhang, H. F., Maslov, K.,


Stoica, G. & Wang, L. V. Functional photoacoustic microscopy for high-resolution and noninvasive in vivo imaging. _Nat. Biotechnol._ 24, 848–851 (2006). Article  CAS  PubMed  Google


Scholar  * Dasa, M. K., Markos, C., Janting, J. & Bang, O. Multispectral photoacoustic sensing for accurate glucose monitoring using a supercontinuum laser. _J. Opt. Soc. Am. B_ 36,


A61–A65 (2019). Article  ADS  CAS  Google Scholar  * Ghazaryan, A., Ovsepian, S. V. & Ntziachristos, V. Extended near-infrared optoacoustic spectrometry for sensing physiological


concentrations of glucose. _Front. Endocrinol._ 9, 112 (2018). Article  Google Scholar  * Pleitez, M. A. et al. In vivo Noninvasive monitoring of glucose concentration in human epidermis by


mid-infrared pulsed photoacoustic spectroscopy. _Anal. Chem._ 85, 1013–1020 (2013). Article  CAS  PubMed  Google Scholar  * Taruttis, A. & Ntziachristos, V. Advances in real-time


multispectral optoacoustic imaging and its applications. _Nat. Photonics_ 9, 219–227 (2015). Article  ADS  CAS  Google Scholar  * Wang, L. V. & Yao, J. A practical guide to photoacoustic


tomography in the life sciences. _Nat. Methods_ 13, 627–638 (2016). Article  CAS  PubMed  PubMed Central  Google Scholar  * Dasa, M. K. et al. High-pulse energy supercontinuum laser for


high-resolution spectroscopic photoacoustic imaging of lipids in the 1650-1850 nm region. _Biomed. Opt. Express_ 9, 1762–1770 (2018). Article  CAS  PubMed  PubMed Central  Google Scholar  *


Wu, M., Jansen, K., van der Steen, A. F. W. & van Soest, G. Specific imaging of atherosclerotic plaque lipids with two-wavelength intravascular photoacoustics. _Biomed. Opt. Express_ 6,


3276–3286 (2015). Article  CAS  PubMed  PubMed Central  Google Scholar  * Wang, P., Wang, H.-W., Sturek, M. & Cheng, J.-X. Bond-selective imaging of deep tissue through the optical


window between 1600 and 1850 nm. _J. Biophotonics_ 5, 25–32 (2012). Article  PubMed  Google Scholar  * Hudson, S. V. et al. Targeted noninvasive imaging of EGFR-expressing orthotopic


pancreatic cancer using multispectral optoacoustic tomography. _Cancer Res._ 74, 6271–6279 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * Li, R. et al. Assessing breast


tumor margin by multispectral photoacoustic tomography. _Biomed. Opt. Express_ 6, 1273–1281 (2015). Article  PubMed  PubMed Central  Google Scholar  * Wu, W., Wang, P., Cheng, J.-X. &


Xu, X.-M. Assessment of white matter loss using bond-selective photoacoustic imaging in a rat model of contusive spinal cord injury. _J. Neurotrauma_ 31, 1998–2002 (2014). Article  PubMed 


PubMed Central  Google Scholar  * Wang, X., Xie, X., Ku, G., Wang, L. V. & Stoica, G. Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using


high-resolution photoacoustic tomography. _J. Biomed. Opt._ 11, 024015 (2006). Article  ADS  PubMed  CAS  Google Scholar  * Fainchtein, R., Stoyanov, B. J., Murphy, J. C., Wilson, D. A.


& Hanley, D. F. Local determination of hemoglobin concentration and degree of oxygenation in tissue by pulsed photoacoustic spectroscopy. _Biomed. Optoacoust._ 3916, 19–33 (2000).


Article  ADS  CAS  Google Scholar  * Shi, J. et al. High-resolution, high-contrast mid-infrared imaging of fresh biological samples with ultraviolet-localized photoacoustic microscopy. _Nat.


Photonics_ 13, 609–615 (2019). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Tzoumas, S., Nunes, A., Deliolanis, N. C. & Ntziachristos, V. Effects of multispectral


excitation on the sensitivity of molecular optoacoustic imaging. _J. Biophotonics_ 8, 629–637 (2015). Article  PubMed  Google Scholar  * Taruttis, A., Claussen, J., Razansky, D. &


Ntziachristos, V. Motion clustering for deblurring multispectral optoacoustic tomography images of the mouse heart. _J. Biomed. Opt._ 17, 016009 (2012). Article  ADS  PubMed  Google Scholar


  * Lloyd, L. B., Riseman, S. M., Burnham, R. K., Eyring, E. M. & Farrow, M. M. Fourier transform photoacoustic spectrometer. _Rev. Sci. Instrum._ 51, 1488–1492 (1980). Article  ADS  CAS


  Google Scholar  * Teng, Y. C. & Royce, B. S. H. Quantitative Fourier transform IR photoacoustic spectroscopy of condensed phases. _Appl. Opt._ 21, 77–80 (1982). Article  ADS  CAS 


PubMed  Google Scholar  * Busse, G. & Bullemer, B. Use of the opto-acoustic effect for rapid scan Fourier spectroscopy. _Infrared Phys._ 18, 255–256 (1978). Article  ADS  CAS  Google


Scholar  * Farrow, M. M., Burnham, R. K. & Eyring, E. M. Fourier‐transform photoacoustic spectroscopy. _Appl. Phys. Lett._ 33, 735–737 (1978). Article  ADS  Google Scholar  * Diddams, S.


A. The evolving optical frequency comb [Invited]. _J. Opt. Soc. Am. B_ 27, B51–B62 (2010). Article  CAS  Google Scholar  * Ye, J. & Cundiff, S. T. (eds). _Femtosecond Optical Frequency


Comb: Principle, Operation, and Applications_. (Springer Science & Business Media, Boston, 2005). * Newbury, N. R. Searching for applications with a fine-tooth comb. _Nat. Photonics_ 5,


186–188 (2011). Article  ADS  CAS  Google Scholar  * Truong, G.-W. et al. Accurate frequency referencing for fieldable dual-comb spectroscopy. _Opt. Express_ 24, 30495–30504 (2016). Article


  ADS  CAS  PubMed  Google Scholar  * Karhu, J. et al. Broadband photoacoustic spectroscopy of 14CH4 with a high-power mid-infrared optical frequency comb. _Opt. Lett._ 44, 1142–1145 (2019).


Article  ADS  CAS  PubMed  Google Scholar  * Sadiek, I., Mikkonen, T., Vainio, M., Toivonen, J. & Foltynowicz, A. Optical frequency comb photoacoustic spectroscopy. _Phys. Chem. Chem.


Phys._ 20, 27849–27855 (2018). Article  CAS  PubMed  Google Scholar  * Coddington, I., Newbury, N. & Swann, W. Dual-comb spectroscopy. _Optica_ 3, 414–426 (2016). Article  ADS  Google


Scholar  * Okubo, S. et al. Ultra-broadband dual-comb spectroscopy across 1.0–1.9 µm. _Appl. Phys. Express_ 8, 082402 (2015). Article  ADS  CAS  Google Scholar  * Ycas, G. et al.


High-coherence mid-infrared dual-comb spectroscopy spanning 2.6 to 5.2 μm. _Nat. Photonics_ 12, 202–208 (2018). Article  ADS  CAS  Google Scholar  * Muraviev, A. V., Smolski, V. O., Loparo,


Z. E. & Vodopyanov, K. L. Massively parallel sensing of trace molecules and their isotopologues with broadband subharmonic mid-infrared frequency combs. _Nat. Photonics_ 12, 209–214


(2018). Article  ADS  CAS  Google Scholar  * Waxman, E. M. et al. Intercomparison of open-path trace gas measurements with two dual frequency comb spectrometers. _Atmos. Meas. Tech._ 10,


3295–3311 (2017). Article  ADS  PubMed  PubMed Central  Google Scholar  * Ideguchi, T., Poisson, A., Guelachvili, G., Picqué, N. & Hänsch, T. W. Adaptive real-time dual-comb


spectroscopy. _Nat. Commun._ 5, 3375 (2014). Article  ADS  PubMed  PubMed Central  CAS  Google Scholar  * Klocke, J. L. et al. Single-shot sub-microsecond mid-infrared spectroscopy on


protein reactions with quantum cascade laser frequency combs. _Anal. Chem._ 90, 10494–10500 (2018). Article  CAS  PubMed  Google Scholar  * Ideguchi, T. et al. Coherent Raman spectro-imaging


with laser frequency combs. _Nature_ 502, 355–358 (2013). Article  ADS  CAS  PubMed  Google Scholar  * Ideguchi, T. et al. Microfluidic single-particle chemical analyzer with dual-comb


coherent Raman spectroscopy. _Opt. Lett._ 43, 4057–4060 (2018). Article  ADS  CAS  PubMed  Google Scholar  * Schmid, T. Photoacoustic spectroscopy for process analysis. _Anal. Bioanal.


Chem._ 384, 1071–1086 (2006). Article  CAS  PubMed  Google Scholar  * Kuthirummal, N. Listening to nanomaterials: photoacoustic spectroscopy. _J. Chem. Educ._ 86, 1238–1240 (2009). Article 


CAS  Google Scholar  * Yao, J. & Wang, L. V. Sensitivity of photoacoustic microscopy. _Photoacoustics_ 2, 87–101 (2014). Article  PubMed  PubMed Central  Google Scholar  * Sinclair, L.


C. et al. Invited article: a compact optically coherent fiber frequency comb. _Rev. Sci. Instrum._ 86, 081301 (2015). Article  ADS  CAS  PubMed  Google Scholar  * Roy, J., Deschênes, J.-D.,


Potvin, S. & Genest, J. Continuous real-time correction and averaging for frequency comb interferometry. _Opt. Express_ 20, 21932–21939 (2012). Article  ADS  PubMed  Google Scholar  *


Lehman, J., Yung, C., Tomlin, N., Conklin, D. & Stephens, M. Carbon nanotube-based black coatings. _Appl. Phys. Rev._ 5, 011103 (2018). Article  ADS  CAS  Google Scholar  * Maslov, K.


& Wang, L. V. Photoacoustic imaging of biological tissue with intensity-modulated continuous-wave laser. _J. Biomed. Opt._ 13, 024006 (2008). Article  ADS  PubMed  Google Scholar  *


Sivaramakrishnan, M., Maslov, K., Zhang, H. F., Stoica, G. & Wang, L. V. Limitations of quantitative photoacoustic measurements of blood oxygenation in small vessels. _Phys. Med. Biol._


52, 1349–1361 (2007). Article  PubMed  Google Scholar  * Cai, D. K., Neyer, A., Kuckuk, R. & Heise, H. M. Optical absorption in transparent PDMS materials applied for multimode


waveguides fabrication. _Opt. Mater._ 30, 1157–1161 (2008). Article  ADS  CAS  Google Scholar  * Turek, I., Tarjányi, N., Martinček, I. & Káčik, D. Effect of mechanical stress on optical


properties of polydimethylsiloxane. _Opt. Mater._ 36, 965–970 (2014). Article  ADS  CAS  Google Scholar  * McGoverin, C. M., Lewis, K., Yang, X., Bostrom, M. P. G. & Pleshko, N. The


contribution of bone and cartilage to the near-infrared spectrum of osteochondral tissue. _Appl. Spectrosc._ 68, 1168–1175 (2014). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 


* Schliesser, A., Picqué, N. & Hänsch, T. W. Mid-infrared frequency combs. _Nat. Photonics_ 6, 440–449 (2012). Article  ADS  CAS  Google Scholar  * Robertson, C. W. & Williams, D.


Lambert Absorption coefficients of water in the infrared. _J. Opt. Soc. Am._ 61, 1316–1320 (1971). Article  ADS  CAS  Google Scholar  * Newbury, N. R., Coddington, I. & Swann, W. C.


Sensitivity of coherent dual-comb spectroscopy. _Opt. Express_ 18, 7929–7945 (2010). Article  ADS  CAS  PubMed  Google Scholar  * Krainer, L. et al. Tunable picosecond pulse-generating laser


with repetition rate exceeding 10 GHz. _Electron. Lett._ 38, 225–227 (2002). Article  ADS  CAS  Google Scholar  * Bartels, A., Heinecke, D. & Diddams, S. A. Passively mode-locked 10 GHz


femtosecond Ti:sapphire laser. _Opt. Lett._ 33, 1905–1907 (2008). Article  ADS  CAS  PubMed  Google Scholar  * Torres-Company, V. & Weiner, A. M. Optical frequency comb technology for


ultra-broadband radio-frequency photonics. _Laser Photonics Rev._ 8, 368–393 (2014). Article  ADS  Google Scholar  * Carlson, D. R. et al. Ultrafast electro-optic light with subcycle


control. _Science_ 361, 1358–1363 (2018). Article  ADS  CAS  PubMed  Google Scholar  * Pasquazi, A. et al. Micro-combs: a novel generation of optical sources. _Phys. Rep._ 729, 1–81 (2018).


Article  ADS  MathSciNet  CAS  MATH  Google Scholar  * Kippenberg, T. J., Holzwarth, R. & Diddams, S. A. Microresonator-based optical frequency combs. _Science_ 332, 555–559 (2011).


Article  ADS  CAS  PubMed  Google Scholar  * Xue, X., Wang, P.-H., Xuan, Y., Qi, M. & Weiner, A. M. Microresonator Kerr frequency combs with high conversion efficiency: high-efficiency


microcombs. _Laser Photonics Rev._ 11, 1600276 (2017). Article  ADS  Google Scholar  * Karpov, M., Pfeiffer, M. H. P., Liu, J., Lukashchuk, A. & Kippenberg, T. J. Photonic chip-based


soliton frequency combs covering the biological imaging window. _Nat. Commun._ 9, 1146 (2018). Article  ADS  PubMed  PubMed Central  CAS  Google Scholar  * Laser Institute of America.


_American National Standard for Safe Use of Lasers_. (Laser Institute of America, 2014). * Cox, B., Laufer, J. G., Arridge, S. R. & Beard, P. C. Quantitative spectroscopic photoacoustic


imaging: a review. _J. Biomed. Opt._ 17, 061202 (2012). Article  ADS  PubMed  Google Scholar  * Wildi, T., Voumard, T., Brasch, V., Yilmaz, G., & Herr, T. Photo-acoustic dual-frequency


comb spectroscopy. arXiv:2004.04691 (2020). * Herman, D. et al. Femtosecond timekeeping: slip-free clockwork for optical timescales. _Phys. Rev. Appl_ 9, 044002 (2018). Article  ADS  CAS 


Google Scholar  Download references ACKNOWLEDGEMENTS We would like to acknowledge G. Ycas for fruitful discussion and development of FPGA code. We acknowledge helpful comments from E.M.


Waxman and A.S. Kowligy. J.T.F. and A.M.G. both conducted their research while the authors held NRC Research Associateship awards at NIST. AUTHOR INFORMATION Author notes * Gabriel M.


Colacion Present address: Optical Science and Engineering, University of New Mexico, 1313 Goddard, SE, Albuquerque, NM, 87106, USA * Eli V. Hoenig Present address: Pritzker School of


Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, ERC 387, Chicago, IL, 60637, USA * Edgar F. Perez Present address: Institute for Research in Electronics and Applied


Physics, University of Maryland, 8279 Paint Branch Drive, College Park, MD, 20742-3511, USA AUTHORS AND AFFILIATIONS * National Institute of Standards and Technology, Applied Physics


Division, 325 Broadway, Boulder, CO, 80305, USA Jacob T. Friedlein, Esther Baumann, Kimberly A. Briggman, Gabriel M. Colacion, Fabrizio R. Giorgetta, Aaron M. Goldfain, Daniel I. Herman, Eli


V. Hoenig, Jeeseong Hwang, Nathan R. Newbury, Edgar F. Perez, Christopher S. Yung, Ian Coddington & Kevin C. Cossel * Department of Physics, University of Colorado, Boulder, CO, 80309,


USA Esther Baumann, Gabriel M. Colacion, Fabrizio R. Giorgetta, Daniel I. Herman, Eli V. Hoenig & Edgar F. Perez Authors * Jacob T. Friedlein View author publications You can also search


for this author inPubMed Google Scholar * Esther Baumann View author publications You can also search for this author inPubMed Google Scholar * Kimberly A. Briggman View author publications


You can also search for this author inPubMed Google Scholar * Gabriel M. Colacion View author publications You can also search for this author inPubMed Google Scholar * Fabrizio R.


Giorgetta View author publications You can also search for this author inPubMed Google Scholar * Aaron M. Goldfain View author publications You can also search for this author inPubMed 


Google Scholar * Daniel I. Herman View author publications You can also search for this author inPubMed Google Scholar * Eli V. Hoenig View author publications You can also search for this


author inPubMed Google Scholar * Jeeseong Hwang View author publications You can also search for this author inPubMed Google Scholar * Nathan R. Newbury View author publications You can also


search for this author inPubMed Google Scholar * Edgar F. Perez View author publications You can also search for this author inPubMed Google Scholar * Christopher S. Yung View author


publications You can also search for this author inPubMed Google Scholar * Ian Coddington View author publications You can also search for this author inPubMed Google Scholar * Kevin C.


Cossel View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS J.T.F helped design and build the DCPAS system. He also acquired and analyzed the


data and wrote the manuscript. E.B., G.M.C., D.H., E.V.H., E.F.P., and helped design and build the DCPAS system. E.B., K.A.B., F.R.G., J.H., N.R.N., I.C., and K.C.C. provided guidance for


experimental design and interpretation of data. F.R.G. and K.C.C. helped with data analysis. C.S.Y. and J.H. fabricated the PDMS/VACNT sample. A.M.G., K.A.B., and J.T.F. performed the PAS


measurements with the optical parametric oscillator. All authors approve of the submitted version of the manuscript and agree to be personally accountable for his/her own contributions and


to ensure that questions related to the accuracy or integrity of any part of the work, even ones in which he/she was not personally involved, are appropriately investigated, resolved, and


the resolution documented in the literature. CORRESPONDING AUTHOR Correspondence to Kevin C. Cossel. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests.


ADDITIONAL INFORMATION PEER REVIEW INFORMATION _Nature Communications_ thanks Frans Harren, Markku Vainio and the other, anonymous, reviewer(s) for their contribution to the peer review of


this work. Peer reviewer reports are available. PUBLISHER’S NOTE Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION PEER REVIEW FILE RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International


License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,


provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons


license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by


statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit


http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Friedlein, J.T., Baumann, E., Briggman, K.A. _et al._ Dual-comb photoacoustic


spectroscopy. _Nat Commun_ 11, 3152 (2020). https://doi.org/10.1038/s41467-020-16917-y Download citation * Received: 27 November 2019 * Accepted: 29 May 2020 * Published: 19 June 2020 * DOI:


https://doi.org/10.1038/s41467-020-16917-y SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not


currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative