The fundamental properties of laser-induced plasma in liquid water, such as the ultrafast electron migration and solvation, have not yet been clarified. We use 1650-nm femtosecond laser pulses to induce the plasma in a stable free-flowing water film under the strong field ionization mechanism. Moreover, we adopt intense terahertz (THz) pulses to probe the ultrafast temporal evolution of quasifree electrons of the laser-induced plasma in water on the subpicosecond scale. For the first time, the THz wave absorption signal with a unique two-step decay characteristic in time domain is demonstrated, indicating the significance of electron solvation in water. We employ the Drude model combined with the multilevel intermediate model and particle-in-a-box model to simulate and analyze the key information of quasifree electrons, such as the frequency-domain absorption characteristics and solvation ratio. In particular, we observe that the solvation capacity of liquid water decreases with the increase of pumping energy. Up to ∼50 % of quasifree electrons cannot be captured by traps associated with the bound states as the pumping energy increases to 90 μJ / pulse. The ultrafast electron evolution in liquid water revealed by the optical-pump/THz-probe experiment provides further insights into the formation and evolution mechanisms of liquid plasma.
Liquid water is the most common liquid on Earth, and its basic physical and chemical properties are of great significance.1,2 In particular, electron migration and solvation associated with photoionization in liquid water is a unique and complex process, which has attracted many researchers to explore various segments of the process. For example, Laenen et al.3 used midinfrared femtosecond probe pulses, and Birkedal et al.4 employed near-infrared/visible femtosecond probe pulses to observe the latter stage of the electron solvation process. Recently, Savolainen et al.5 studied the initial stage of this process by adopting terahertz (THz) probe pulses and revealed the initial delocalization length of the excited excess electrons in water. These previous studies constitute a motion picture of an excess electron, which first delocalizes in the conduction band and then is solvated within about 1 to 2 ps.4–13 This process is considered to include many transient and highly active intermediates, resulting in a large blueshift in the absorption spectrum, from the THz or gigahertz region to the visible range. Nonetheless, previous studies mainly investigated the behavior of low-density quasifree electrons excited by small pump power density, which is far from enough for further understanding the transient evolution of photoexcited plasma in liquid water.
Laser-induced plasma has attracted widespread attention due to its diverse physical phenomena. For example, plasma excited in the gaseous environment has been widely investigated and applied in many fields, such as femtosecond pulse compression,14–16 broadband THz wave emission,17–20 and self-healing of light spot.21 However, compared with the numerous studies of plasma in the gaseous environment, the properties of photoexcited plasma in the liquid environment have not been fully revealed, especially in water with large molecular density and strong molecular interactions. Theoretically, laser-induced plasma in liquid water will be accompanied by more complex and stronger nonlinear effects than those in gas, since water has a higher nonlinear coefficient,22 a lower excitation threshold, and a higher electron density.23,24 In the latest research, the femtosecond laser-induced plasma in water has been used to emit strong THz waves.25–28 In these studies, the physical mechanism of plasma formation in liquid water is considered as an extension of the theory of gas plasma. However, the effect of water molecules on plasma is not fully unraveled. In addition, to the best of our knowledge, there are no reports on the electron migration and solvation process related to strong laser excited plasma in water, especially the initial stage of the interaction between the laser and water on the subpicosecond timescale.
In this paper, we use intense THz pulses (peak electric field strength of ) and a stable free-flowing water film to probe the absorption characteristics of the plasma in liquid water. The femtosecond laser pulses at 1650 nm are used to excite the plasma in liquid water under the strong field mechanism. We observe that the absorption characteristics of the THz wave by plasma show a clear two-step decay process in the time domain, which can be explained by the partial quasifree electron solvation. We connect this THz wave absorption characteristic by plasma with the change of quasifree electron density through the Drude model and combine the multilevel intermediate model with the particle-in-a-box model to analyze the temporal evolution process of quasifree electrons. The absorption changes of THz waves and proportions of solvated electrons in water under various pumping energies are also revealed.
Experimental Results and Analysis
The schematic diagram of the optical-pump THz probe spectroscopy system is shown in Fig. 1(a). The intense THz pulses are generated by the organic 4-N,N-dimethylamino-4′-N′-methyl-stilbazoliumtosylate (DAST) crystal under excitation with the femtosecond laser pulses (0.4 mJ per pulse, repetition rate of 1 kHz, and center wavelength of 1550 nm) from an optical parametric amplifier (OPA, Spectra Physics). To reach the diffraction limited focusing and allow for high field (peak electric field strength of ), we use the THz bullet lambda cubic wave front correction scheme.29 The idler-wave output from the OPA is used as pump laser for photoionization (repetition rate of 1 kHz and center wavelength of 1650 nm) and focused onto the water film to excite the plasma by a lens with a focal length of 4 in. placed in front of parabolic mirror 1 (PM1). The linearly polarized THz pulses reflected from PM1 and linearly polarized optical pump pulses are transmitted collinearly through the liquid water. In order to efficiently record the evolution process on the subpicosecond time scale, we adjust the relative time delay between the THz probe and 1650 nm pump pulse with an accuracy of 10 fs. The residual NIR is filtered out using a set of low-pass filters (QMC Instruments Ltd.) with cut-off frequencies of 18 THz. The THz energy is controlled by two THz polarizers. The THz pulses passing through the water film are collected by two off-axis parabolic mirrors and collinearly focused with the probe pulses onto the GaP crystal for electro-optical sampling. Since the measured broadband THz pulse will be additionally distorted by the frequency-dependent factors in the electro-optical detection crystal (such as phase matching, reflection, dispersive propagation, and absorption), we use the full complex response function of the GaP detector to reconstruct all of the spectral information collected in this work.30–32 A wire-guided, gravity-driven apparatus is employed in this work, which can produce optically stable, flowing thin liquid film and avoid interference caused by a traditional cuvette.33,34 Two aluminum wires with diameters of are separated by to generate a stable free-flowing water film. The thickness of the water film is by controlling the water flow rate of . An optical second-harmonic intensity autocorrelator is used to measure and calibrate the thickness.27 Within the time interval of 1 ms between two adjacent pulses, the distance that the gravity-driven water film falls is about , which is much larger than the pump spot size of . Therefore, the adopted free-flowing water film in our study can effectively avoid the artifacts from cuvette windows and the heat accumulation effect in the excitation volume. All of the measurements were conducted at the room temperature of . For all of the measurements, the system was purged with dry nitrogen gas to improve the signal-to-noise ratio.
Figure 1(b) shows the frequency-domain amplitudes of the THz electric field pulses propagating through the water film and nitrogen gas with the optical pump beam blocked, and the black curves in Fig. 2(a) (and its inset) show the corresponding time-domain waveforms. Compared with the spectrum in nitrogen gas, the liquid water film has a very strong absorption effect for the high-frequency components of THz pulses. The THz spectrum transmitted through water film extends from to with the peak at , and there is almost a flat phase for the entire measured spectrum. We fix the THz scanning time delay () at the peak of the THz pulse and scan the pump–probe time delay () to obtain the transient temporal evolution of the THz wave absorption change () induced by the formed plasma in water on the subpicosecond time scale, as shown in Fig. 2(b). In comparison, the temporal evolution of THz wave absorption characteristics induced by the formed plasma in nitrogen gas (without the water film) is also presented in Fig. 2(b). The temporal characteristics and corresponding mechanism in the gas phase have been widely investigated.35,36 The results of nitrogen gas in our experiment are consistent with previous studies. The experimental evidence in Fig. 2 also shows that in the presence of photoexcited plasma, the THz wave is significantly absorbed in water with a peak modulation more than five times larger than that in nitrogen gas under the same optical pump power. The corresponding time-domain THz pulses under the maximal absorption by plasma in water and nitrogen gas are shown in Fig. 2(a) and its inset, respectively. Since the absorption of THz waves by plasma is isotropic, the change of polarization angles between the THz pulse and pump pulse has no effect on our results. Moreover, we find that the THz wave absorption by plasma in water rises rapidly within a few hundred femtoseconds and then exhibits a unique two-step decay characteristic, which is very different from the monotonic decay characteristic of nitrogen gas. The first step of decay for water is an ultrafast process within the subpicosecond scale, and the following is a relatively slow process with a decay time exceeding over hundreds of picoseconds [as shown in the inset of Fig. 2(b)]. We use a double exponential decay model11 to describe the above two-step decay process (, where , are derived from measured results of solvation rate , which is closely related to the pump energy; here, , for the pump energy of ). The fitted results, depicted by the orange line in Fig. 2(b) (and its inset), match the experimental data well, and the obtained time constants are and , respectively.
To explain the unique absorption characteristics in liquid water, particles produced or modified by photoionization, including ions, excited neutral molecules, and electrons in bound and quasifree states, need to be considered. The THz response may be influenced by the movement of electrons or the disturbance of the nuclear spectrum of the disturbed water network. Related to the charged nuclear response, for example, the change in the molar extinction coefficient of water caused by NaCl dissolution is .37 For the plasma density with the order of in this work, the produced charged nuclear response contributes to the THz wave absorption change, which is several orders of magnitude smaller than the measured peak modulation magnitude. In addition, the temperature increase due to photoionization in water is also a factor to be considered. In our experiment, the adopted stable free-flowing water film is constantly refreshed, which significantly limits the thermal accumulation. Considering the factors such as the heat capacity and absorption coefficient of liquid water, we estimate that each optical pulse () increases the temperature of the excited volume of water film by no more than 0.7°C, and the corresponding absorption change is .38 However, the above reasons cannot explain the large THz wave absorption change with maximal value of nearly 50% in this work. Therefore, we believe that the measured modulation mainly originates from the strong interaction between low effective mass electrons and THz pulses.39 In liquid water, the electrons in the plasma induced by the ultrashort laser pulses exist in two states, which are solvated electrons and quasifree electrons. The solvated electrons are considered to be trapped in the bound state and do not contribute to the conductivity. They are generated from ionization within about 1 to 2 ps and have a long lifetime (∼hundreds of picoseconds) with the characteristic of an optical resonance centered at 720 nm (bandwidth of about 330 nm).5,13,23 We do not rule out the weak contribution of solvated electrons to THz wave absorption, but this obviously cannot explain the dramatic THz response in the initial stage from no pump beam to the maximal absorption on the subpicosecond scale in this work. Therefore, we mainly attribute the rapid rise of the modulation signal [as shown in Fig. 2(b)] associated with the THz wave absorption characteristics to the contribution of quasielectrons in the initial stage of plasma formation.
The decay process for quasifree electrons of plasma in liquid water includes geminate recombination, diffusion, and solvation processes. The solvation process dominates on the subpicosecond scale, and the decay time is at least three orders of magnitude faster than other processes.13,40,41 Therefore, we believe that the first step of decay in Fig. 2(b) is mainly caused by the solvation of partial quasifree electrons. This also explains the monotonic decay characteristic in nitrogen gas due to the absence of the solvation process. The second-step of decay is considered to be the result of the quasifree electrons being recaptured by the traps associated with the bound states during the demise of hydrated electrons, which mainly reflects the lifetime of hydrated electrons. In general, the formation of hydrated electrons from photoexcited plasma occurs on a time scale of 1 to 2 ps.5,13,23 In this work, the modulation signal rapidly decreases with the time constant of in the first-step of decay. This can be explained by the multilevel intermediate model related to the solvation process.10 Specifically, some quasifree electrons transition to the first intermediates with weak THz wave absorption capacity on a time scale of . Here, we use the Drude model to evaluate the scattering process for the quasifree electron transitioning into the first intermediate and analyze the THz wave absorption characteristics based on the quasifree electron density. First, Eq. (1) is adopted to establish the relationship between the alternation of the composite conductivity  and the THz wave absorption change:42 which is extremely weak for liquid water with a strong absorption coefficient in the THz range. Then, the relationship between the induced change in composite conductivity and the quasifree electron density is established, as shown in Eq. (2):
The experimental frequency-domain characteristics of THz wave absorption  under maximal modulation by plasma in liquid water are shown in Fig. 3. The simulated results based on the above Drude model are also shown as the orange dotted line. An obvious absorption peak near is demonstrated. The simulation results are basically consistent with the experimental data. Note that the measured absorption is slightly lower than the simulated value in the low-frequency range. This is mainly because the THz pulse tail carrying low-frequency components always lags behind the main pulse () to pass through the plasma, and then the quasifree electron density has been partially attenuated. In addition, the absorption dip at in the measurement is mainly derived from the characteristics of the DAST source itself. According to the Drude model, the lifetime of quasifree electrons is calculated to be , which is approximately in agreement with the experimental measurement of . In addition, the calculated maximum value of the quasifree electron density of the plasma excited in liquid water is with the adopted pump energy of .
Furthermore, we measure the transient evolution processes of THz wave absorption by plasma in water with the excitation of different pumping energies and calculate the corresponding modulation degree at 1.5 THz under the maximal absorption . The results are shown in Fig. 4. We find that the laser-excited peak quasifree electron density has a linear relationship with the pumping energy in the range from 40 to , which can be explained by the mechanism of the strong field ionization.43 In our experiment, the peak quasifree electron density tends to saturate when the pumping energy is higher than .
As the pumping energy increases, the plateau feature after the first step of decay shows an upward trend, as shown in the inset of Fig. 4(a). This is due to the decrease in the proportion of quasifree electrons that are solvated with the increase of pumping energy. For the excited low-density quasifree electrons (such as the condition of the pumping energy ), the ratio of the quasifree electron density () to the solvated electron density () at the equilibrium state depends on the relative mobilities, i.e., . and represent the mobility of solvated electrons and quasifree electrons, respectively. The electron mobility is related to the lifetime of the electron (). Therefore, for quasifree electrons with a lifetime of and solvated electrons with a lifetime of , more than 99.8% of the quasifree electrons will be solvated at the equilibrium state. We believe that the low-density quasifree electrons will be completely solvated after . However, the energy-dependent experimental data of Fig. 4(a) show that the modulation degree in the plateau region is about 0.8% at the pumping energy of , which comes from thermal effects. According to Ref. 38, the modulation degree of 0.8% is caused by a temperature difference of , which is basically consistent with the temperature estimation in our present setup. In addition, the temperature response will not be significant as the pumping energy increases. As shown in Fig. 4(a), as the pumping energy increases from 40 to , the modulation degree at the plateau stage increases significantly from to . However, even the thermal effect induced by () only causes a modulation degree of 2% in THz absorption. Therefore, considering the increase of the quasifree electron density in the plasma and the factor that the plasma will destroy the network structure of liquid water, we attribute the increase of the modulation degree at the plateau stage with increasing pump pulse energy [the inset of Fig. 4(a)] mainly to the saturation of the traps associated with the bound states. When the traps are saturated, more quasifree electrons cannot be relaxed through the solvation process. Excluding the thermal effect, we calculate the ratio of quasifree electrons solvated with different pumping energies [i.e., solvation rate , as shown in Fig. 4(b)]. Specifically, the low-density quasifree electrons are almost all captured by the traps related to the bound state with the weak pump pulse excitation (). However, the traps gradually tend to be saturated, and the solvation ratio decreases as the pump pulse energy exceeds the threshold ( in this work). When the quasifree electron density in plasma is saturated, the ratio of quasifree electrons solvated drops to with pumping energy higher than .
The size of the electron in the plasma after photoionization is also important information and can be regarded as the radius of gyration of the spin density. The particle-in-a-box model is adopted in our work, which describes the solvation of electrons as a continuously changing process.5,44 At the initial stage of plasma formation, electrons usually have a large radius of gyration, which is basically equal to their initial average ejection length. The electrons at this time can be regarded as quasifree electrons in an unbound state. Then, the water molecules solvate the electrons and form a cage in a bound state. The electron’s gyration radius gradually decreases until the quasifree electrons are solvated within a time scale of 1 to 2 ps. The relationship between the resonance absorption frequency () of electrons and the radius of gyration () can be expressed as5,44 Based on Eq. (3), we calculate the initial radius of gyration of quasifree electrons in liquid water at an absorption peak of 1.5 THz as in our experiment, which is in agreement with the previously published results.5 As the radius of gyration continuously reduces by the restraint of the water molecule cage, the resonant absorption frequency continues to blueshift and leaves the THz observation window after . Eventually, the resonant absorption frequency converges to the near visible spectrum at around 720 nm with the corresponding radius of gyration of about 2.5 Å.5
THz time-domain spectroscopy can provide instantaneous information about the frequency-dependent conductivity of the medium excited by optical pump pulses. In this work, we explore the optical-pump/THz-probe experiment to investigate the transient dynamic evolution of quasifree electrons of plasma in water on the subpicosecond time scale. The time-resolved transient THz wave absorption measurements prove that the quasifree electrons in the laser-excited plasma in water will be solvated rapidly. The measured temporal transmittance data contain valuable information, which can be supported by a combination of the Drude model, the multilevel intermediate model, and the particle-in-a-box model. We conclude that the initial delocalization radius of gyration of quasifree electrons is , the scattering time is , and the residence time of electrons in the bound state is up to . The peak quasifree electron density increases linearly, while the pumping energy increases from 40 to and tends to saturate with a value of . In addition, we observe that as the quasifree electron density increases, the traps related to the bound states appear to saturate, resulting in a large number of quasifree electrons not being completely solvated. When the pump energy increases to , the proportion of solvated electrons gradually decreases to about . This work can provide valuable insights on the fundamental aspects of the charge transport process in water and lay the foundation for further understanding the physicochemical properties and transient evolution of the femtosecond laser pulse excited plasma in water. We believe that the proposed method can also be applied to study the physical mechanisms of plasma deformation for the ionic aqueous solution and other liquids.
This work was supported by the Beijing Natural Science Foundation (Grant No. JQ18015), the National Natural Science Foundation of China (Grant Nos. 61935001, 61905271, 12074272), the Guangdong Basic and Applied Basic Research Foundation (Grant No. 2020A1515011083), and the Beijing Advanced Innovation Center for Imaging Theory and Technology. The authors declare no conflicts of interest.
Yong Tan received his BS degree from Chang’an University, Xi’an, China, in 2016. He is currently pursuing his PhD in instrument science and technology at Beijing Institute of Technology, Beijing, China. His research interests include instrument science and technology and terahertz emission. His research primarily concerns generation and detection of terahertz waves.
Hang Zhao received her BS degree from Shenyang University of technology, Shenyang, China, in 2014. She received her MS degree from Capital Normal University, Beijing, China, in 2017. She is currently pursuing her PhD in instrument science and technology at Beijing Institute of Technology, Beijing, China. Her research primarily concerns generation and detection of terahertz waves.
Rui Zhang received his BEng degree in electronic information engineering from Xidian University, Xi’an, China, in 2010. He received his PhD in mechanics from Peking University, Beijing, China, in 2016. His research interests include terahertz spectroscopy and imaging, biomedical application of terahertz techniques, and terahertz wave generation. He is now an associate professor at Shenzhen Institute of Advanced Technology.
Yuejin Zhao received his PhD in optical engineering from Beijing Institute of Technology, Beijing, China, in 1990. He is currently a professor with the School of Optoelectronics, Beijing Institute of Technology. He has disciplinary responsibility as a team leader of “Instrument Science and Technology” with Beijing Institute of Technology. His current research interests include terahertz imaging technology based on MEMS-infrared imaging technology, space-optical technology, and intelligent photoelectric instrument development.
Cunlin Zhang received his PhD in photoelectric engineering and optical engineering from Beijing Institute of Technology, Beijing, China, in 1992. He is currently a professor with the Capital Normal University and also the director of Beijing Key Laboratory of Terahertz Spectrum and Imaging. He has disciplinary responsibility as a team leader of “Terahertz Spectrum and Imaging” with Beijing Institute of Technology. His current research interests include terahertz and infrared spectrum and imaging.
Xi-Cheng Zhang is the Parker Givens Chair and former director of the Institute of Optics, University of Rochester (UR), New York, USA. With his BS degree in 1982 from Peking University, he received his MS and PhD degrees in physics from Brown University, Rhode Island, in 1983 and 1986, respectively. He has been recognized with many honors and awards, including the Humboldt Prize, from the Alexander von Humboldt Foundation (2018); IRMMW-THz Kenneth Button Prize (2014); OSA William F. Meggers Award (2012); IEEE Photonics Society William Streifer Scientific Achievement Award (2011); and the Rensselaer William H. Wiley 1866 Award (2009), among others.
Liangliang Zhang received her PhD in instrument science and technology from Beijing Institute of Technology, Beijing, China, in 2008. She is currently a professor with the Department of Physics, Capital Normal University. Her research interests include a wide variety of topics in the area of terahertz emission and terahertz detection.