Ultra-fast switches are essential devices for basic research and technological development. Here, we demonstrate that aqueous solutions of sodium iodide can be used for this purpose. When pumped with an intense optical pulse at 400 nm, these water-based liquids display large and fast responses in the terahertz range, around 1 THz. In a 9M NaI solution at a temperature comprised between 10 and 50 °C, the relative variation of the terahertz peak transmission drops by 20% at the pump–probe overlap and recovers with a fast time constant of ∼70 fs. As the optical properties of the liquid vary on a timescale shorter than the terahertz cycle, it is possible to tailor the shape of the transmitted terahertz fields. In this way, we demonstrate the frequency upshifting of terahertz radiation from about 1 to 3 THz and beyond with an efficiency of 4%.

At the enormous pressures and temperatures of 300 GPa and 7000 K, molecular dynamics calculations predict that liquid water should be metallic and contribute to the magnetic fields of Uranus and Neptune.1 More recently, experiments demonstrated that thin layers of liquid water could form a metallic solution for a few seconds, via charge transfer from a sodium potassium alloy.2 However, liquid water can also display a metallic-like transient response for an ultra-short amount of time following photoionization. When water molecules in the liquid phase are photoionized via multi-photon optical absorption, the initially excited electron wavefunction can have a radius of3–5 ∼40 Å. The electron is thus delocalized over 9000 water molecules. This “conduction band electron” is excited when the energy of the optical excitation exceeds the “bandgap of liquid water,” located somewhere between6 10 and 12 eV. The wavefunction of the initially delocalized electron shrinks to less than7,8 3 Å within ∼0.2 ps, resulting in the familiar absorption of the relaxed, solvated aqueous electron peaking in the near infrared.9 

When electromagnetic radiation propagates from one medium to another with different optical properties, the Fresnel coefficients determine the properties of the reflected and refracted waves at the interface. While the frequency of the electromagnetic radiation is unchanged, the wavevectors of the reflected and transmitted beams change. By analogy with a spatial boundary, reflection and refraction are also possible when the properties of a medium vary abruptly in time, that is, when the variation is much quicker than the duration of a period of the electromagnetic wave. This temporal boundary10–14 can modify the frequency of the electromagnetic wave rather than its wavevector. Following Morgenthaler’s seminal work,10 the frequency up-shifting and down-shifting of electromagnetic radiation by a temporal boundary were demonstrated in gas plasmas,15,16 semiconductors,17–25 transparent conducting oxides,26–30 and metamaterials.31–33 

The Hamm group demonstrated3 that the photoionization of liquid water can be probed in the terahertz (THz) range. A broad, transient absorption feature centered at 1.5 THz decayed within 200 fs and was associated with the delocalized electrons. Previous optical-pump and optical-probe experiments demonstrated that the dissolution in water of iodide anions (I) increases the number of photo-generated electrons.34 Here, we show that a 9M water solution of sodium iodide displays a significant (−20%) and short-lived (∼70 fs) optical-pump terahertz-probe (OPTP) signal at a temperature between 10 and 50 °C. As the optical properties of the liquid are modulated on a timescale that is shorter than the duration of the THz optical cycle, we demonstrate that the THz fields transmitted by a liquid solution are upshifted in frequency from ∼1 to ∼3 THz with an efficiency of 4%.

OPTP measurements, sometimes called time-resolved terahertz spectroscopy (TRTS), were performed on a series of aqueous solutions containing 1, 3, 5, 7, 9, and 11M of sodium iodide (NaI). The samples were prepared by dissolving sodium iodide salt (99.5% purity, supplied by S3 Chemicals) in ultra-pure milli-Q water with a conductivity of 0.055 µS/cm. The initial solution with a concentration of 1 mol/l (or M) was prepared by dissolving 15 g (0.1 mol) in 100 ml of water. After this initial 1M NaI liquid sample was measured, we obtained the 3M solution by dissolving an additional 30 g of NaI in the same solution. This process was repeated until a total of 165 g of NaI was dissolved in 100 ml of pure water, which is close to the solubility limit at room temperature (∼11M). The equilibrium terahertz properties of these solutions were reported previously.35 The experiments were performed in a nitrogen atmosphere with less than 10% humidity and at three temperatures: 10, 30, and 50 °C. A recirculating chiller stabilized the temperature of the reservoir of the liquid within ±0.05 °C. In order to avoid unwanted pump–probe signals from the windows of a static cell or sample holder,36 the measurements were performed on a recirculating liquid jet. We measured a flowing water jet without windows. A water jet is a freestanding sheet of water that is generated by a nozzle; please see Ref. 37 for further details.

The details of the experimental OPTP setup are given elsewhere.38 In short, the broadband terahertz (THz) probe pulses are generated by two-color plasma filamentation in nitrogen gas. A first, collecting off-axis parabolic mirror (OAPM) collimates the THz radiation. The residual optical radiation is filtered out by a hole in this collecting OAPM and by a high-resistivity silicon slab (Eksma optics) placed after the first OAPM. In order to highlight the effects of pulse shaping and frequency upshifting, paper filters were added to reduce the THz bandwidth. A second OAPM focuses the THz radiation into the sample position, where we placed the flowing, windowless liquid jet sample. A third OAPM collimates the THz beam transmitted by the sample, and a fourth OAPM focuses the THz radiation into the detection crystal. The optical pump pulses at ∼400 nm propagate collinearly with the THz beam via a hole in the middle of the second OAPM. The residual pump energy is filtered out by a hole in the third OAPM and by polyethylene filters placed between the third and fourth OAPM. Examples of the scheme of an optical setup dedicated to non-linear terahertz spectroscopy can be found in Ref. 39.

The THz pulses are detected via electro-optical sampling in a 0.1 mm thick gallium phosphide (GaP) crystal with a temporal resolution of 25 fs. The delay of the sampling beam is the electro-optical sampling delay, tEOS. In order to avoid spurious contributions by THz reflections, Fourier transformation (FT) is performed after zero padding the THz fields collected between tEOS = −1 ps and tEOS = +1 ps to 256 points and with a Hanning window function. This corresponds to a frequency resolution of ∼0.16 THz. The magnitude of the FT of a THz field corresponds to its spectrum. The optical pump pulses are ∼50 fs long and centered at 400 nm, as generated by second harmonic generation. In the following, we will refer to the temporal delay of the pump beam as the pump–probe delay, tPP. The pump–probe delay was varied in steps of 50 fs. The peak intensity of a single optical pump pulse amounted to ∼1.2 TW/cm2. We optimized the laser compression to obtain the fastest and largest transient signals from the liquid samples. The spot sizes of the optical and THz beams were measured with a camera (400 µm). As a wavelength of 400 µm corresponds to 0.75 THz, the optical beam does not excite the probe uniformly for frequencies lower than 0.75 THz. For these reasons, the THz spectra shown here are reliable only above 0.75 THz.

Figure 1 displays the THz fields [Fig. 1(a)] and spectra [Fig. 1(b)] transmitted by an exemplary liquid sample at equilibrium, i.e., when the optical pump pulses are blocked. The reference THz probe field [green in Fig. 1(a)] is obtained by turning off the liquid jet, i.e., by measuring the transmission of the empty path. The duration of the almost single-cycle THz fields is 1 ps. The THz probe field transmitted by the liquid sample is shown in black in Fig. 1(a). The optical properties of the aqueous sample are linked to the delayed arrival of the THz pulse transmitted by the sample and to the reduced amplitude of the transmitted THz field. In the following, the maximum value of the THz field transmitted will be referred to as the “THz peak,” which is marked with E0 for the liquid sample in Fig. 1(a). The THz spectra transmitted by reference (green curve) and sample (black curve) are shown in Fig. 1(b) with an identical scaling factor, i.e., the absolute quantities can be compared. The maximum magnitude of the FT spectrum transmitted by the liquid will be referred to as M0 in the following. All the THz probe spectra are centered around 1 THz, which is typical for intense sources, such as table-top-based tilted-front optical rectification in lithium niobate40–42 and accelerator-based superradiant TELBE43 and TeraFERMI.44,45 As indicated by the gray area in Fig. 1(b), the noise level of the spectra is lower than 1%.

FIG. 1.

Terahertz (THz) source and linear optical properties of a 9M NaI aqueous solution at the temperature of 50 °C. (a) The THz field transmitted by the reference, a ∼17 µm thick layer of nitrogen gas (empty path), is shown in green. The THz field transmitted by the free-flowing liquid jet is shown in black. The maximum THz field transmitted by the liquid is marked with E0. The linear optical properties of the liquid are evidenced by the temporal delay of the THz peak and its lower amplitude. The temporal resolution is 25 fs. (b) The corresponding magnitudes of the Fourier transformations (FTs) are shown in black for the reference and in green for the sample. The noise level is indicated with the gray shaded area and is below 1%. The FT magnitudes are shown with identical scaling. The spectrum peaks at around 1 THz. The frequency resolution is ∼0.16 THz.

FIG. 1.

Terahertz (THz) source and linear optical properties of a 9M NaI aqueous solution at the temperature of 50 °C. (a) The THz field transmitted by the reference, a ∼17 µm thick layer of nitrogen gas (empty path), is shown in green. The THz field transmitted by the free-flowing liquid jet is shown in black. The maximum THz field transmitted by the liquid is marked with E0. The linear optical properties of the liquid are evidenced by the temporal delay of the THz peak and its lower amplitude. The temporal resolution is 25 fs. (b) The corresponding magnitudes of the Fourier transformations (FTs) are shown in black for the reference and in green for the sample. The noise level is indicated with the gray shaded area and is below 1%. The FT magnitudes are shown with identical scaling. The spectrum peaks at around 1 THz. The frequency resolution is ∼0.16 THz.

Close modal

The difference between the arrival times of the THz peak transmitted by the empty path and the THz peak transmitted by the sample varied between about 50 and 80 fs for all the water samples studied here. By assuming, for simplicity, an index of refraction equal to 2, this temporal delay corresponds to a thickness of the liquid jet between 15 and 24 µm for all the samples. The ratio between the amplitude of the THz peak transmitted by the liquid sample and the THz peak transmitted by the empty path varied between about 60% and 65% in all the investigated solutions. Our setup is designed to perform high-resolution transient and relative measurements, but is not equally suited to quantify the absolute optical properties of the samples in the THz range because the thickness of the liquid jet is affected by the temperature and salt concentration. The equilibrium optical properties of these liquids can be obtained with THz time-domain spectroscopy and a static cell, as reported earlier.35 

The overall pump-induced response can be quantified by performing a “pump scan” or “fixed-gate scan.” In this experiment, the electro-optical sampling beam is kept at the fixed time delay corresponding to the maximum field transmission (E0), while the pump–probe delay (tPP) is varied. The relative variation of the THz peak field (ΔE/E0 = E0,pumpON/E0,pumpOFF − 1) is measured as a function of the pump–probe delay. The result is shown in Fig. 2(a) for the 9M NaI sample at three temperatures. The transient change of the THz peak drops by 20% at zero pump–probe delay (tPP = 0 ps) and depends weakly on the temperature. To our knowledge, this is the largest OPTP signal ever reported on a liquid sample. After the maximum transient response is reached, the relative variation of the THz transmission decays. This decay can be fit to a bi-exponential function, which is shown with the black curve in Fig. 2(a). The fastest component has a decay constant of τ1 = (70 ± 30) fs and accounts for most of the size of the transient signal. The transient THz signal and its decay constants are similar in all the investigated samples, but the amplitude varies with the salt concentration. As shown in Fig. 2(b), the size of the transient THz signal is roughly proportional to the concentration of sodium iodide in water, thus proving the crucial role of this salt in the processes measured here. A maximum signal of −24% is found for the 11M NaI solution at 50 °C. In this case, however, we have data only at 30 and 50 °C because the salt came out of solution at 10 °C. For these reasons, we highlight in Figs. 2(a) and 3 the results of the 9M NaI water sample, which is stable between 10 and 50 °C. As a further control, we note that similar measurements on pure water revealed much smaller pump–probe signals, about ΔE/E0 ∼ −10−4 at the pump–probe overlap, tPP = 0. The experimental results on pure water will be included in a separate publication.

FIG. 2.

Optical-pump terahertz-probe (OPTP) transient response measured at the terahertz (THz) peak. (a) The curves displays the relative variation of the peak THz field transmitted by the 9M NaI solution, ΔE/E0, as a function of pump–probe delay, tPP. We repeated the measurements by keeping the liquid sample at three temperatures: 10 °C (blue), 30 °C (red), and 50 °C (orange). The dynamics at pump–probe delays larger than tPP > 0 ps can be fit to the same bi-exponential function for all the studied temperatures (black curve). The temporal resolution is 50 fs. (b) We performed measurements as a function of the concentration of sodium iodide: 1, 3, 5, 7, 9, and 11M. The pump–probe signal at pulse overlap increases with the solute concentration. The data are the average of 10× independent measurements at each temperature. The vertical error bars correspond to the standard error of the mean. When not visible, the error bars are smaller than the thickness of the curves.

FIG. 2.

Optical-pump terahertz-probe (OPTP) transient response measured at the terahertz (THz) peak. (a) The curves displays the relative variation of the peak THz field transmitted by the 9M NaI solution, ΔE/E0, as a function of pump–probe delay, tPP. We repeated the measurements by keeping the liquid sample at three temperatures: 10 °C (blue), 30 °C (red), and 50 °C (orange). The dynamics at pump–probe delays larger than tPP > 0 ps can be fit to the same bi-exponential function for all the studied temperatures (black curve). The temporal resolution is 50 fs. (b) We performed measurements as a function of the concentration of sodium iodide: 1, 3, 5, 7, 9, and 11M. The pump–probe signal at pulse overlap increases with the solute concentration. The data are the average of 10× independent measurements at each temperature. The vertical error bars correspond to the standard error of the mean. When not visible, the error bars are smaller than the thickness of the curves.

Close modal
FIG. 3.

Full terahertz (THz) field detected via electro-optical sampling (tEOS), with (red curves) and without (black curves) optical pump set to zero electro-optical delay, tPP = 0 ps. (a) The THz field transmitted when the pump is blocked is shown with the black curve. When the optical pump pulse is turned on at tPP = 0 (black arrow), we detect the red THz field instead. (b) Magnitude of the Fourier transformation (FT) of the THz fields transmitted when the pump is off (black spectrum) and when the pump is on (red spectrum). The FT magnitudes are displayed with an identical scaling factor. From the red and black curves in panel (a), it is possible to estimate the relative variation of the THz field obtained by turning on the optical pump pulses, which is reported in panel (c). We estimate the relative change of the THz spectrum (d) from the spectra in panel (b). The efficiency of the frequency upshifting is estimated to the ratio between the purple shaded area in panel (d) and the gray shaded area in panel (b). The temporal resolution is 25 fs in panels (a) and (c), and the frequency resolution is ∼0.16 THz in panels (b) and (d). As demonstrated in Fig. 2, the transient THz signals are independent of the liquid temperature between 10 and 50 °C, within the experimental error. Thus, the data in all panels of Fig. 3 are averaged over 9× independent measurements: 3× were performed at the temperature of 10 °C, 3× at 30 °C, and 3× at 50 °C. The vertical error bars correspond to the standard error of the mean. When not visible, the error bars are smaller than the thickness of the curves.

FIG. 3.

Full terahertz (THz) field detected via electro-optical sampling (tEOS), with (red curves) and without (black curves) optical pump set to zero electro-optical delay, tPP = 0 ps. (a) The THz field transmitted when the pump is blocked is shown with the black curve. When the optical pump pulse is turned on at tPP = 0 (black arrow), we detect the red THz field instead. (b) Magnitude of the Fourier transformation (FT) of the THz fields transmitted when the pump is off (black spectrum) and when the pump is on (red spectrum). The FT magnitudes are displayed with an identical scaling factor. From the red and black curves in panel (a), it is possible to estimate the relative variation of the THz field obtained by turning on the optical pump pulses, which is reported in panel (c). We estimate the relative change of the THz spectrum (d) from the spectra in panel (b). The efficiency of the frequency upshifting is estimated to the ratio between the purple shaded area in panel (d) and the gray shaded area in panel (b). The temporal resolution is 25 fs in panels (a) and (c), and the frequency resolution is ∼0.16 THz in panels (b) and (d). As demonstrated in Fig. 2, the transient THz signals are independent of the liquid temperature between 10 and 50 °C, within the experimental error. Thus, the data in all panels of Fig. 3 are averaged over 9× independent measurements: 3× were performed at the temperature of 10 °C, 3× at 30 °C, and 3× at 50 °C. The vertical error bars correspond to the standard error of the mean. When not visible, the error bars are smaller than the thickness of the curves.

Close modal

The results in Fig. 2 prove that the optical-pump excitation strongly perturbs the THz properties of NaI water solutions. The pump-induced effect is the largest over a timescale τ1 ∼ 70 fs, which is much shorter than the duration of an optical cycle of the THz field (∼1 ps, see Fig. 1). As this liquid sample could slice or clip the THz field with sub-cycle resolution, we can use it as a temporal boundary for THz radiation.

In a “probe scan” or “fixed-pump scan” experiment, the THz field is fully recorded by delaying the electro-optical sampling beam (tEOS). When the optical pump is blocked, we obtain the THz transmission of the liquid at equilibrium, which is shown with the black trace in Fig. 3(a). When the optical pump is turned on at tPP = 0 [black arrow in Fig. 3(a)], the transmitted THz field is perturbed—see the red trace in Fig. 3(a). It is evident from the raw data that the pump pulse induces a strong asymmetry in the THz pulse, which implies that the pump interaction changes the frequency components of the THz field. The THz spectra obtained by FT are shown in Fig. 3(b). As expected, the spectra are different for pump on [red in Fig. 3(b)] and pump off [black in Fig. 3(b)]. Please note that the THz spectra are shown with the same scaling factor and can be quantitatively compared. When the optical pump pulses perturb the liquid, the transmitted THz magnitude is weaker at ∼2 THz and stronger above ∼3 THz. The purple curve in Fig. 3(c) is equal to the difference between the red and black curves in Fig. 3(a), thus showing the relative change of the THz field transmission. Similarly, Fig. 3(d) shows the relative change of the magnitude of the FT, which highlights the redistribution of THz energy induced by the optical pump pulses. The THz magnitude decreases by about −10% between 1 and 2 THz and increases by about +4% at 3 THz, +1% at 4 THz, and +0.5% at 5 THz. The efficiency of the upshifting from a frequency range at around 1 THz to frequencies close to 3 THz amounts to 4%. We calculate this efficiency from the area of the increased magnitude of the FT between 2.5 and 5.5 THz when the pump pulse is on [purple shaded area in Fig. 3(d)] and divide it by the area of the FT magnitude between 0 and 4 THz when the pump is off [gray shaded area in Fig. 3(b)].

Two main processes could contribute to the experimental observations reported here: absorption by delocalized electrons and the impulsive response of the photo-excited solute. As suggested by the increasing value for increasing NaI content [Fig. 2(b)], the transient THz signal is related to the photoionization of I, which results in the generation of a precursor to the solvated aqueous electron3 and the neutral iodine atom,34 I. The precursor to the solvated aqueous electron was described before as a large, delocalized electron, with a radius of up to 40 Å that decays within ∼200 fs. According to a particle in a box model,3 a delocalized electron with a radius of 40 Å absorbs THz radiation at around 1.5 THz. Thus, the generation of short-lived delocalized electrons is expected to reduce the THz transmission within 200 fs, which is consistent with the experimental observations reported here. A second process involving the THz field emitted by a photo-excited solute in water was described by Ahmed et al.46 They performed OPTP experiments on a free-flowing liquid jet of water containing a photo-absorbing dye (Coumarin). The pump pulse resonantly excites a molecular state of the solute while THz radiation probes the transient response. They found a fast and negative signal most prominent at pump–probe overlap, ΔE/E0 ∼ −6 · 10−3 at tPP = 0. They also performed non-equilibrium molecular dynamics simulations and explained this negative signal as the THz field emitted by the solute and its surrounding water upon photoexcitation of the solute molecule. By extension, also in our case, the photoionization of I results in a charge redistribution in liquid water, for which the emission of a THz field with the opposite phase is expected, thus resulting in a decreased transmission of the THz signal. Further theory and experimental works are required to pinpoint the role of the THz emission triggered by the photoionization of water.

In summary, this Letter demonstrates that water solutions containing sodium iodide display very large and very fast optical-pump terahertz-probe signals at the temperature of 10, 30, and 50 °C. The relative variation of the transmitted THz field decreases by up to −20% at the pump–probe overlap and recovers with a fast time constant of ∼70 fs. Moreover, the spectral content of the THz radiation was manipulated by recording the THz fields at a temporally fixed pump excitation. The sub-cycle slicing of the THz field redistributes its spectral component, resulting in the frequency upshifting to ∼3 THz with an efficiency of 4% in the 9M NaI solution. The results presented here could have significant consequences. In the shorter term, we propose that solutions of sodium iodide should be used as reference samples in the alignment and optimization of OPTP or TRTS experiments on liquids. In addition, this technique could be used to tailor THz fields at low frequency.40–45 In the long term, especially if combined with suitable window materials, this work could open the door to “water electronics.” The advantage could be the fast decay time of the THz photo-switch, τ1 ∼ 70 fs. For example, to our knowledge, the fastest decay time of an OPTP modulation in a typical semiconductor is larger than47–49 >150 fs. The liquid switches reported here could be used at a faster rate or, equivalently, have a larger bandwidth (1/70 fs−1 > 1/150 fs−1).

We acknowledge support from the Cluster of Excellence RESOLV (Grant No. EXC 2033 – 390677874) funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), by the ERC Advanced Grant No. 695437 (THz Calorimetry) and by the DFG Project No. 509442914. We acknowledge support by the Open Access Publication Funds of the Ruhr-Universität Bochum. We are particularly grateful to Thorsten Ockelmann, Sashary Ramos and Martina Havenith. This work was dedicated to my late father.

The authors have no conflicts to disclose.

A.B. and C.H. contributed equally to this work

Adrian Buchmann: Data curation (equal); Formal analysis (supporting); Investigation (equal); Methodology (supporting); Resources (equal); Validation (equal); Writing – review & editing (supporting). Claudius Hoberg: Data curation (equal); Formal analysis (supporting); Investigation (equal); Methodology (supporting); Resources (equal); Software (equal); Validation (equal). Fabio Novelli: Conceptualization (lead); Data curation (equal); Formal analysis (lead); Funding acquisition (equal); Investigation (equal); Methodology (lead); Project administration (lead); Resources (equal); Software (equal); Supervision (lead); Validation (equal); Visualization (lead); Writing – original draft (lead); Writing – review & editing (lead).

The data that support the findings of this study are available within the article.

1.
C.
Cavazzoni
,
G. L.
Chiarotti
,
S.
Scandolo
,
E.
Tosatti
,
M.
Bernasconi
, and
M.
Parrinello
,
Science
283
,
44
(
1999
).
2.
P. E.
Mason
,
H. C.
Schewe
,
T.
Buttersack
,
V.
Kostal
,
M.
Vitek
,
R. S.
McMullen
,
H.
Ali
,
F.
Trinter
,
C.
Lee
,
D. M.
Neumark
,
S.
Thürmer
,
R.
Seidel
,
B.
Winter
,
S. E.
Bradforth
, and
P.
Jungwirth
,
Nature
595
,
673
(
2021
).
3.
J.
Savolainen
,
F.
Uhlig
,
S.
Ahmed
,
P.
Hamm
, and
P.
Jungwirth
,
Nat. Chem.
6
,
697
(
2014
).
4.
T. W.
Kee
,
D. H.
Son
,
P.
Kambhampati
, and
P. F.
Barbara
,
J. Phys. Chem. A
105
,
8434
(
2001
).
5.
N.
Linz
,
S.
Freidank
,
X.-X.
Liang
,
H.
Vogelmann
,
T.
Trickl
, and
A.
Vogel
,
Phys. Rev. B
91
,
134114
(
2015
).
6.
A. P.
Gaiduk
,
T. A.
Pham
,
M.
Govoni
,
F.
Paesani
, and
G.
Galli
,
Nat. Commun.
9
,
247
(
2018
).
7.
L.
Turi
and
P. J.
Rossky
,
Chem. Rev.
112
,
5641
(
2012
).
8.
J. M.
Herbert
,
Phys. Chem. Chem. Phys.
21
,
20538
(
2019
).
9.
F.-Y.
Jou
and
G. R.
Freeman
,
J. Phys. Chem.
83
,
2383
(
1979
).
10.
F. R.
Morgenthaler
,
IEEE Trans. Microwave Theory Tech.
6
,
167
(
1958
).
11.
S. C.
Wilks
,
J. M.
Dawson
, and
W. B.
Mori
,
Phys. Rev. Lett.
61
,
337
(
1988
).
12.
Y.
Xiao
,
D. N.
Maywar
, and
G. P.
Agrawal
,
Opt. Lett.
39
,
574
(
2014
).
13.
B. W.
Plansinis
,
W. R.
Donaldson
, and
G. P.
Agrawal
,
Phys. Rev. Lett.
115
,
183901
(
2015
).
14.
E.
Galiffi
,
R.
Tirole
,
S.
Yin
,
H.
Li
,
S.
Vezzoli
,
P. A.
Huidobro
,
M. G.
Silveirinha
,
R.
Sapienza
,
A.
Alù
, and
J. B.
Pendry
,
Adv. Photonics
4
,
014002
(
2022
).
15.
R. L.
Savage
,
R. P.
Brogle
,
W. B.
Mori
, and
C.
Joshi
,
IEEE Trans. Plasma Sci.
21
,
5
(
1993
).
16.
N.
Yugami
,
T.
Niiyama
,
T.
Higashiguchi
,
H.
Gao
,
S.
Sasaki
,
H.
Ito
, and
Y.
Nishida
,
Phys. Rev. E
65
,
036505
(
2002
).
17.
D.
Hashimshony
,
A.
Zigler
, and
K.
Papadopoulos
,
Phys. Rev. Lett.
86
,
2806
(
2001
).
18.
A.
Nishida
,
N.
Yugami
,
T.
Higashiguchi
,
T.
Otsuka
,
F.
Suzuki
,
M.
Nakata
,
Y.
Sentoku
, and
R.
Kodama
,
Appl. Phys. Lett.
101
,
161118
(
2012
).
19.
I.
Geltner
,
Y.
Avitzour
, and
S.
Suckewer
,
Appl. Phys. Lett.
81
,
226
(
2002
).
20.
M.
Shalaby
,
M.
Peccianti
,
D. G.
Cooke
,
C. P.
Hauri
, and
R.
Morandotti
,
Appl. Phys. Lett.
106
,
051110
(
2015
).
21.
J.
Bühler
,
J.
Allerbeck
,
G.
Fitzky
,
D.
Brida
, and
A.
Leitenstorfer
,
Optica
5
,
821
(
2018
).
22.
B.
Mayer
,
C.
Schmidt
,
J.
Bühler
,
D. V.
Seletskiy
,
D.
Brida
,
A.
Pashkin
, and
A.
Leitenstorfer
,
New J. Phys.
16
,
063033
(
2014
).
23.
K.
Lee
,
J.
Park
,
S.
Lee
,
S.
Baek
,
J.
Park
,
F.
Rotermund
, and
B.
Min
,
Nanophotonics
11
,
2045
(
2022
).
24.
F.
Miyamaru
,
C.
Mizuo
,
T.
Nakanishi
,
Y.
Nakata
,
K.
Hasebe
,
S.
Nagase
,
Y.
Matsubara
,
Y.
Goto
,
J.
Pérez-Urquizo
,
J.
Madéo
, and
K. M.
Dani
,
Phys. Rev. Lett.
127
,
053902
(
2021
).
25.
A.
Schiff-Kearn
,
L.
Gingras
,
S.
Bernier
,
N.
Chamanara
,
K.
Agarwal
,
J.-M.
Ménard
, and
D. G.
Cooke
,
Commun. Phys.
4
,
1
(
2021
).
26.
J.
Bohn
,
T. S.
Luk
,
C.
Tollerton
,
S. W.
Hutchings
,
I.
Brener
,
S.
Horsley
,
W. L.
Barnes
, and
E.
Hendry
,
Nat. Commun.
12
,
1017
(
2021
).
27.
P.
Guo
,
R. D.
Schaller
,
L. E.
Ocola
,
B. T.
Diroll
,
J. B.
Ketterson
, and
R. P. H.
Chang
,
Nat. Commun.
7
,
12892
(
2016
).
28.
J.
Bohn
,
T. S.
Luk
,
S.
Horsley
, and
E.
Hendry
,
Optica
8
,
1532
(
2021
).
29.
K.
Pang
,
M. Z.
Alam
,
Y.
Zhou
,
C.
Liu
,
O.
Reshef
,
K.
Manukyan
,
M.
Voegtle
,
A.
Pennathur
,
C.
Tseng
,
X.
Su
,
H.
Song
,
Z.
Zhao
,
R.
Zhang
,
H.
Song
,
N.
Hu
,
A.
Almaiman
,
J. M.
Dawlaty
,
R. W.
Boyd
,
M.
Tur
, and
A. E.
Willner
,
Nano Lett.
21
,
5907
(
2021
).
30.
Y.
Zhou
,
M. Z.
Alam
,
M.
Karimi
,
J.
Upham
,
O.
Reshef
,
C.
Liu
,
A. E.
Willner
, and
R. W.
Boyd
,
Nat. Commun.
11
,
2180
(
2020
).
31.
K.
Lee
,
J.
Son
,
J.
Park
,
B.
Kang
,
W.
Jeon
,
F.
Rotermund
, and
B.
Min
,
Nat. Photonics
12
,
765
(
2018
).
32.
K.
Lee
,
J.
Park
,
J.
Son
,
B. J.
Kang
,
W. T.
Kim
,
S. C.
Lee
,
B.
Min
, and
F.
Rotermund
,
Opt. Express
27
,
12762
(
2019
).
33.
J.
Tunesi
,
L.
Peters
,
J. S.
Totero Gongora
,
L.
Olivieri
,
A.
Fratalocchi
,
A.
Pasquazi
, and
M.
Peccianti
,
Phys. Rev. Res.
3
,
1
(
2021
).
34.
J. A.
Kloepfer
,
V. H.
Vilchiz
,
V. A.
Lenchenkov
,
A. C.
Germaine
, and
S. E.
Bradforth
,
J. Chem. Phys.
113
,
6288
(
2000
).
35.
P.
Zalden
,
L.
Song
,
X.
Wu
,
H.
Huang
,
F.
Ahr
,
O. D.
Mücke
,
J.
Reichert
,
M.
Thorwart
,
P. K.
Mishra
,
R.
Welsch
,
R.
Santra
,
F. X.
Kärtner
, and
C.
Bressler
,
Nat. Commun.
9
,
2142
(
2018
).
36.
M.
Sajadi
,
M.
Wolf
, and
T.
Kampfrath
,
Opt. Express
23
,
28985
(
2015
).
37.
A.
Buchmann
,
C.
Hoberg
, and
M.
Havenith
, in (
IEEE
).
38.
C.
Hoberg
,
P.
Balzerowski
, and
M.
Havenith
,
AIP Adv.
9
,
035348
(
2019
).
39.
H. Y.
Hwang
,
S.
Fleischer
,
N. C.
Brandt
,
B. G.
Perkins
,
M.
Liu
,
K.
Fan
,
A.
Sternbach
,
X.
Zhang
,
R. D.
Averitt
, and
K. A.
Nelson
,
J. Mod. Opt.
62
,
1447
(
2015
).
40.
J.
Hebling
,
K.-L.
Yeh
,
M. C.
Hoffmann
,
B.
Bartal
, and
K. A.
Nelson
,
J. Opt. Soc. Am. B
25
,
B6
(
2008
).
41.
F.
Meyer
,
N.
Hekmat
,
T.
Vogel
,
A.
Omar
,
S.
Mansourzadeh
,
F.
Fobbe
,
M.
Hoffmann
,
Y.
Wang
, and
C. J.
Saraceno
,
Opt. Express
27
,
30340
(
2019
).
42.
H.
Hirori
,
A.
Doi
,
F.
Blanchard
, and
K.
Tanaka
,
Appl. Phys. Lett.
98
,
091106
(
2011
).
43.
B.
Green
,
S.
Kovalev
,
V.
Asgekar
,
G.
Geloni
,
U.
Lehnert
,
T.
Golz
,
M.
Kuntzsch
,
C.
Bauer
,
J.
Hauser
,
J.
Voigtlaender
,
B.
Wustmann
,
I.
Koesterke
,
M.
Schwarz
,
M.
Freitag
,
A.
Arnold
,
J.
Teichert
,
M.
Justus
,
W.
Seidel
,
C.
Ilgner
,
N.
Awari
,
D.
Nicoletti
,
S.
Kaiser
,
Y.
Laplace
,
S.
Rajasekaran
,
L.
Zhang
,
S.
Winnerl
,
H.
Schneider
,
G.
Schay
,
I.
Lorincz
,
A. A.
Rauscher
,
I.
Radu
,
S.
Mährlein
,
T. H.
Kim
,
J. S.
Lee
,
T.
Kampfrath
,
S.
Wall
,
J.
Heberle
,
A.
Malnasi-Csizmadia
,
A.
Steiger
,
A. S.
Müller
,
M.
Helm
,
U.
Schramm
,
T.
Cowan
,
P.
Michel
,
A.
Cavalleri
,
A. S.
Fisher
,
N.
Stojanovic
, and
M.
Gensch
,
Sci. Rep.
6
,
22256
(
2016
).
44.
F.
Novelli
,
C. Y.
Ma
,
N.
Adhlakha
,
E. M.
Adams
,
T.
Ockelmann
,
D.
Das Mahanta
,
P.
Di Pietro
,
A.
Perucchi
, and
M.
Havenith
,
Appl. Sci.
10
,
5290
(
2020
).
45.
C.
Millon
,
J.
Schmidt
,
S.
Ramos
,
E. P.
van Dam
,
A.
Buchmann
,
C.
Saraceno
, and
F.
Novelli
,
AIP Advances
12
,
115319
(
2022
).
46.
S.
Ahmed
,
A.
Pasti
,
R. J.
Fernández-Terán
,
G.
Ciardi
,
A.
Shalit
, and
P.
Hamm
,
J. Chem. Phys.
148
,
234505
(
2018
).
47.
S. S.
Prabhu
,
S. E.
Ralph
,
M. R.
Melloch
, and
E. S.
Harmon
,
Appl. Phys. Lett.
70
,
2419
(
1997
).
48.
H.
Němec
,
A.
Pashkin
,
P.
Kužel
,
M.
Khazan
,
S.
Schnüll
, and
I.
Wilke
,
J. Appl. Phys.
90
,
1303
(
2001
).
49.
G.
Segschneider
,
F.
Jacob
,
T.
Löffler
,
H. G.
Roskos
,
S.
Tautz
,
P.
Kiesel
, and
G.
Döhler
,
Phys. Rev. B
65
,
125205
(
2002
).