Far infrared spectroscopy of water at different temperatures: GHz to THz dielectric spectroscopy of water

A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields

Simulating the dielectric spectra of solvents requires the nuanced definition of inter- and intra-molecular forces. Non-polarizable force fields, while thoroughly benchmarked for dielectric applications, do not capture all the spectral features of solvents, such as water. Conversely, polarizable force fields have been largely untested in the context of dielectric spectroscopy but include charge and dipole fluctuations that contribute to intermolecular interactions. We benchmark non-polarizable force fields and the polarizable force fields AMOEBA03 and HIPPO for liquid water and find that the polarizable force fields can capture all the experimentally observed spectral features with varying degrees of accuracy. However, the non-polarizable force fields miss at least one peak. To diagnose this deficiency, we decompose the liquid water spectra from polarizable force fields at multiple temperatures into static and induced dipole contributions and find that the peak originates from induced dipole contributions. Broadening our inquiry to other solvents parameterized with the AMOEBA09 force field, we demonstrate good agreement between the experimental and simulated dielectric spectra of methanol and formamide. To produce these spectra, we develop a new computational approach to calculate the dielectric spectrum via the fluctuation dissipation theorem. This method minimizes the error in both the low and high frequency portions of the spectrum, improving the overall accuracy of the simulated spectrum and broadening the computed frequency range.

Measurement of Water Uptake and States in Nafion Membranes Using Humidity-Controlled Terahertz Time-Domain Spectroscopy

Perfluorinated sulfonic acid ionomers are well known for their unique water uptake properties and chemical/mechanical stability. Understanding their performance-stability trade-offs is key to realizing membranes with optimal properties. Terahertz time-domain spectroscopy has been demonstrated to resolve water states inside industrially relevant membranes, producing qualitatively agreeable results to conventional gravimetric analysis and prior demonstrations. Using the proposed humidity-controlled terahertz time-domain spectroscopy, here we quantify this detailed water information inside commercially available Nafion membranes at various humidities for direct comparison against literature values from dynamic vapor sorption, differential scanning calorimetry, and Fourier transform infrared spectroscopy on selected samples. Using this technique therefore opens up opportunities for rapid future parameter space investigation for membrane optimization.

Terahertz spectroscopy of thick and diluted water solutions

While bright terahertz sources are used to perform nonlinear experiments, they can be advantageous for high-precision linear measurements of opaque samples. By placing the sample away from the focus, nonlinearities can be suppressed, and sizeable amounts of transmitted radiation detected. Here, this approach is demonstrated for a 0.5 mm thick layer of liquid water in a static sample holder. Variations of the index of refraction as small as (7 ± 2) · 10-4 were detected at 0.58 THz for an aqueous salt solution containing ten millimoles of sodium chloride. To my knowledge, this precision is unprecedented in time-domain spectroscopy studies of diluted aqueous systems or other optically thick and opaque materials.

Organohydrogel-based transparent terahertz absorber via ionic conduction loss

The fast-growing terahertz technologies require high-performance terahertz absorber for suppressing electromagnetic interference. Since the dissipation mechanism in terahertz band usually focuses on electronic conduction loss, almost all terahertz absorbers are constructed with electronically conducting materials being opaque, which limits their applications in scenarios requiring high visible transmittance. Here, we demonstrate a transparent terahertz absorber based on permittivity-gradient elastomer-encapsulated-organohydrogel. Our organohydrogel-based terahertz absorber exhibits a high absorbing property (average reflection loss of 49.03 dB) in 0.5-4.5 THz band with a thin thickness of 700 μm and a high average visible transmittance of 85.51%. The terahertz absorbing mechanism mainly derives from the ionic conduction loss of the polar liquid in organohydrogel. Besides, the hydrophobic and adhesive elastomer coating endows this terahertz absorber high absorbing stability and interfacial adhesivity. This work paves a viable way to designing transparent terahertz absorbers.

Photothermal effects of terahertz-band and optical electromagnetic radiation on human tissues

The field of wireless communication has witnessed tremendous advancements in the past few decades, leading to more pervasive and ubiquitous networks. Human bodies are continually exposed to electromagnetic radiation, but typically this does not impact the body as the radiation is non-ionizing and the waves carry low power. However, with progress in the sixth generation (6G) of wireless networks and the adoption of the spectrum above 100 GHz in the next few years, higher power radiation is needed to cover larger areas, exposing humans to stronger and more prolonged radiation. Also, water has a high absorption coefficient at these frequencies and could lead to thermal effects on the skin. Hence, there is a need to study the radiation effects on human tissues, specifically the photothermal effects. In this paper, we present a custom-built, multi-physics model to investigate electromagnetic wave propagation in human tissue and study its subsequent photothermal effects. The proposed finite-element model consists of two segments-the first one estimates the intensity distribution along the beam path, while the second calculates the increase in temperature due to the wave distribution inside the tissue. We determine the intensity variation in the tissue using the radiative transfer equation and compare the results with Monte Carlo analysis and existing analytical models. The intensity information is then utilized to predict the rise in temperature with a bio-heat transfer module, powered by Pennes' bioheat equation. The model is parametric, and we perform a systematic photothermal analysis to recognize the crucial variables responsible for the temperature growth inside the tissue, particularly for terahertz and near-infrared optical frequencies. Our numerical model can serve as a benchmark for studying the high-frequency radiation effects on complex heterogeneous media such as human tissue.

Terahertz Kerr Effect of Liquids

In recent years, tremendous advancements have been made in various technologies such as far-infrared, low-frequency Raman, and two-dimensional (2D) Raman terahertz (THz) spectroscopies. A coherent method has emerged from numerous experimental and theoretical investigations of molecular dynamics in liquids by comparing linear and non-linear spectroscopic techniques. Intermolecular hydrogen bond vibration, molecular reorientation motion, and interaction between molecule/ionic solute and hydrogen bonds have been demonstrated to occur in the THz region, which are closely related to their physical/chemical properties and structural dynamics. However, precise probing of various modes of motion is difficult because of the complexity of the collective and cooperative motion of molecules and spectral overlap of related modes. With the development of THz science and technology, current state-of-the-art THz sources can generate pulsed electric fields with peak intensities of the order of microvolts per centimeter (MV/cm). Such strong fields enable the use of THz waves as the light source for non-linear polarization of the medium and in turn leads to the development of the emerging THz Kerr effect (TKE) technique. Many low-frequency molecular motions, such as the collective directional motion of molecules and cooperative motion under the constraint of weak intermolecular interactions, are resonantly excited by an intense THz electric field. Thus, the TKE technique provides an interesting prospect for investigating low-frequency dynamics of different media. In view of this, this paper first summarizes the research work on TKE spectroscopy by taking a solid material without low-frequency molecular motions as an example. Starting from the principle of TKE technology and its application in investigating the properties of solid matter, we have explored the low-frequency molecular dynamics of liquid water and aqueous solutions using TKE. Liquid water is a core of life and possesses many extraordinary physical and biochemical properties. The hydrogen bond network plays a crucial role in these properties and is the main reason for its various kinetic and thermodynamic properties, which differ from those of other liquids. However, the structure of the hydrogen bond network between water and solutes is not well known. Therefore, evaluating the hydrogen bond-related kinetic properties of liquid water is important.

Highly sensitive detection of broadband terahertz waves using aqueous salt solutions

Water-based coherent detection of broadband terahertz (THz) wave has been recently proposed with superior performances, which can alleviate the limited detection bandwidth and high probe laser energy requirement in the solid- and air-based detection schemes, respectively. Here, we demonstrate that the water-based detection method can be extended to the aqueous salt solutions and the sensitivity can be significantly enhanced. The THz coherent detection signal intensity scales linearly with the third-order nonlinear susceptibility χ⁽³⁾ or quadratically with the linear refractive index η₀ of the aqueous salt solutions, while the incoherent detection signal intensity scales quadratically with χ⁽³⁾ or quartically with η₀, proving the underlying mechanism is the four-wave mixing. Both the coherent and incoherent detection signal intensities appear positive correlation with the solution concentration. These results imply that the liquid-based THz detection scheme could provide a new technique to measure χ⁽³⁾ and further investigate the physicochemical properties in the THz band for various liquids.

Attenuated Total Reflection at THz Wavelengths: Prospective Use of Total Internal Reflection and Polariscopy

Capabilities of the attenuated total reflection (ATR) at THz wavelengths for increased sub-surface depth characterisation of (bio-)materials are presented. The penetration depth of a THz evanescent wave in biological samples is dependent on the wavelength and temperature and can reach 0.1–0.5 mm depth, due to the strong refractive index change ∼0.4 of the ice-water transition; this is quite significant and important when studying biological samples. Technical challenges are discussed when using ATR for uneven, heterogeneous, high refractive index samples with the possibility of frustrated total internal reflection (a breakdown of the ATR reflection mode into transmission mode). Local field enhancements at the interface are discussed with numerical/analytical examples. Maxwell’s scaling is used to model the behaviour of absorber–scatterer inside the materials at the interface with the ATR prism for realistic complex refractive indices of bio-materials. The modality of ATR with a polarisation analysis is proposed, and its principle is illustrated, opening an invitation for its experimental validation. The sensitivity of the polarised ATR mode to the refractive index between the sample and ATR prism is numerically modelled and experimentally verified for background (air) spectra. The design principles of polarisation active optical elements and spectral filters are outlined. The results and proposed concepts are based on experimental conditions at the THz beamline of the Australian Synchrotron.

Low-frequency vibrations of water molecules in DNA minor groove

Water molecules around the DNA form the hydration shell having different structural and dynamical features in different regions of the double helix. In the DNA minor groove, water molecules are highly ordered and in the case of AT nucleotide sequence, the formation of a hydration spine is observed. In the present research, the vibrations of the hydration spine have been studied to establish the mode of translational vibrations of water molecules in the DNA low-frequency spectra (water-spine vibrations). Using the developed phenomenological model with the parameters determined for different nucleotides of the DNA fragment CGCGAATTCGCG, the frequencies of vibrations of the hydration spine have been obtained within 185 ± 20 cm[Formula: see text] depending on type of nucleotide. The obtained frequencies are in the same region as the translational vibrations of water molecules in the bulk. To select the mode of water-spine vibrations from those modes that are present in the bulk water, the dynamics of DNA with different nucleotide contents has been analyzed, and the possible influence of heavy water has been estimated. The determined features of the mode of water vibrations in the hydration spine of DNA minor groove indicate that this mode may be observed in the experimental spectra.

Ultrafast hydrogen bond dynamics of liquid water revealed by terahertz-induced transient birefringence

The fundamental properties of water molecules, such as their molecular polarizability, have not yet been clarified. The hydrogen bond network is generally considered to play an important role in the thermodynamic properties of water. The terahertz (THz) Kerr effect technique, as a novel tool, is expected to be useful in exploring the low-frequency molecular dynamics of liquid water. Here, we use an intense and ultrabroadband THz pulse (peak electric field strength of 14.9 MV/cm, centre frequency of 3.9 THz, and bandwidth of 1-10 THz) to resonantly excite intermolecular modes of liquid water. Bipolar THz field-induced transient birefringence signals are observed in a free-flowing water film. We propose a hydrogen bond harmonic oscillator model associated with the dielectric susceptibility and combine it with the Lorentz dynamic equation to investigate the intermolecular structure and dynamics of liquid water. We mainly decompose the bipolar signals into a positive signal caused by hydrogen bond stretching vibration and a negative signal caused by hydrogen bond bending vibration, indicating that the polarizability perturbation of water presents competing contributions under bending and stretching conditions. A Kerr coefficient equation related to the intermolecular modes of water is established. The ultrafast intermolecular hydrogen bond dynamics of water revealed by an ultrabroadband THz pump pulse can provide further insights into the transient structure of liquid water corresponding to the pertinent modes.

Polyelectrolytes induce water-water correlations that result in dramatic viscosity changes and nuclear quantum effects

Ions interact with water via short-ranged ion-dipole interactions. Recently, an additional unexpected long-ranged interaction was found: The total electric field of ions influences water-water correlations over tens of hydration shells, leading to the Jones Ray effect, a 0.3% surface tension depression. Here, we report such long-range interactions contributing substantially to both molecular and macroscopic properties. Femtosecond elastic second harmonic scattering (fs-ESHS) shows that long-range electrostatic interactions are remarkably strong in aqueous polyelectrolyte solutions, leading to an increase in water-water correlations. This increase plays a role in the reduced viscosity, which changes more than two orders of magnitude with polyelectrolyte concentration. Using D₂O instead of H₂O shifts both the fs-ESHS and the viscosity curve by a factor of ~10 and reduces the maximum viscosity value by 20 to 300%, depending on the polyelectrolyte. These phenomena cannot be explained using a mean-field approximation of the solvent and point to nuclear quantum effects.

Molecular polarizability anisotropy of liquid water revealed by terahertz-induced transient orientation

Reaction pathways of biochemical processes are influenced by the dissipative electrostatic interaction of the reagents with solvent water molecules. The simulation of these interactions requires a parametrization of the permanent and induced dipole moments. However, the underlying molecular polarizability of water and its dependence on ions are partially unknown. Here, we apply intense terahertz pulses to liquid water, whose oscillations match the timescale of orientational relaxation. Using a combination of terahertz pump / optical probe experiments, molecular dynamics simulations, and a Langevin dynamics model, we demonstrate a transient orientation of their dipole moments, not possible by optical excitation. The resulting birefringence reveals that the polarizability of water is lower along its dipole moment than the average value perpendicular to it. This anisotropy, also observed in heavy water and alcohols, increases with the concentration of sodium iodide dissolved in water. Our results enable a more accurate parametrization and a benchmarking of existing and future water models.

The intermolecular dynamics of liquid water impact most biological processes. Here, the authors use intense terahertz electromagnetic pulses to generate a transient, out-of-equilibrium state of the water network to show that the molecules become oriented and probe the polarizability of this anisotropic state.

Electromagnetic-radiation absorption by water

Why does a microwave oven work? How does biological tissue absorb electromagnetic radiation? Astonishingly, we do not have a definite answer to these simple questions because the microscopic processes governing the absorption of electromagnetic waves by water are largely unclarified. This absorption can be quantified by dielectric loss spectra, which reveal a huge peak at a frequency of the exciting electric field of about 20 GHz and a gradual tailing off toward higher frequencies. The microscopic interpretation of such spectra is highly controversial and various superpositions of relaxation and resonance processes ascribed to single-molecule or molecule-cluster motions have been proposed for their analysis. By combining dielectric, microwave, THz, and far-infrared spectroscopy, here we provide nearly continuous temperature-dependent broadband spectra of water. Moreover, we find that corresponding spectra for aqueous solutions reveal the same features as pure water. However, in contrast to the latter, crystallization in these solutions can be avoided by supercooling. As different spectral contributions tend to disentangle at low temperatures, this enables us to deconvolute them when approaching the glass transition under cooling. We find that the overall spectral development, including the 20 GHz feature (employed for microwave heating), closely resembles the behavior known for common supercooled liquids. Thus water's absorption of electromagnetic waves at room temperature is not unusual but very similar to that of glass-forming liquids at elevated temperatures, deep in the low-viscosity liquid regime, and should be interpreted along similar lines.

Far infrared spectroscopy of hydrogen bonding collective motions in complex molecular systems

Far infrared spectroscopy as a tool for the study of inter and intramolecular interactions in complex molecular structures.

The origin of the Debye relaxation in liquid water and fitting the high frequency excess response

Debye relaxation in the microwave/THz region of the dielectric spectra of water is related to propagation of defects.

New terahertz dielectric spectroscopy for the study of aqueous solutions

We present the development of a high precision, tunable far-infrared (terahertz) frequency-domain dielectric spectrometer for studying the dynamics of biomolecules in aqueous solutions in the gigahertz-to-terahertz frequency. As an important benchmark system, we report on the measurements of the absorption and refractive index for liquid water in the frequency range from 5 GHz to 1.12 THz (0.17-37.36 cm(-1) or 0.268-60 mm). The system provides a coherent radiation source with power up to 20 mW in the gigahertz-to-terahertz region. The dynamic range of our instrument reaches 10(12) and the system achieves a spectral resolution of less than 100 Hz. The temperature of samples can be controlled precisely with error bars of ±0.02 °C from 0 °C to 90 °C. Given these attributes, our spectrometer provides unique capabilities for the accurate measurement of even very strongly absorbing materials such as aqueous solutions.

Collective hydration dynamics of guanidinium chloride solutions and its possible role in protein denaturation: a terahertz spectroscopic study

The remarkable ability of guanidinium chloride (GdmCl) to denature proteins is a well studied yet controversial phenomenon; the exact molecular mechanism is still debatable, especially the role of hydration dynamics, which has been paid less attention. In the present contribution, we have addressed the issue of whether the collective hydrogen bond dynamics of water gets perturbed in the presence of GdmCl and its possible impact on the denaturation of a globular protein human serum albumin (HSA), using terahertz (THz) time domain spectroscopy (TTDS) in the frequency range of 0.3-2.0 THz. The collective hydrogen bond dynamics is determined by fitting the obtained complex dielectric response in a multiple Debye relaxation model. To compare the results, the studies were extended to two more salts: tetramethylguanidinium chloride (TMGdmCl) and sodium chloride (NaCl). It was concluded that the change in hydration dynamics plays a definite role in the protein denaturation process.

What is the primary mover of water dynamics?

We present evidence that the microscopic origin of both the excess wing and the main relaxation process of pure water is the same.

Intermolecular vibrational study in liquid water and ice by using far infrared spectroscopy with synchrotron radiation of MIRRORCLE 20

Far infrared absorption measurements for distilled water and ice Ih were performed in the frequency range from 20 cm(-1) to 1000 cm(-1) with Fourier Transform Infrared Spectrometer (FTIR) utilizing SR of a portable synchrotron. Four vibrational bands were separated from measured spectra in liquid water. We found that a peak frequency of 40±1 cm(-1) did not depend on the temperature in a range between 10.0°C and 70.0°C; however, any low energy excitation modes were not observed in the ice spectrum. It is concluded that this band is caused by collective vibrations specific to the dynamical structure in liquid water.

Modeling terahertz heating effects on water

We apply Kirchhoff's heat equation to model the influence of a CW terahertz beam on a sample of water, which is assumed to be static. We develop a generalized model, which easily can be applied to other liquids and solids by changing the material constants. If the terahertz light source is focused down to a spot with a diameter of 0.5 mm, we find that the steady-state temperature increase per milliwatt of transmitted power is 1.8?C/mW. A quantum cascade laser can produce a CW beam in the order of several milliwatts and this motivates the need to estimate the effect of beam power on the sample temperature. For THz time domain systems, we indicate how to use our model as a worst-case approximation based on the beam average power. It turns out that THz pulses created from photoconductive antennas give a negligible increase in temperature. As biotissue contains a high water content, this leads to a discussion of worst-case predictions for THz heating of the human body in order to motivate future detailed study. An open source Matlab implementation of our model is freely available for use at www.eleceng.adelaide.edu.au/thz.

Simultaneous composition and thickness measurement of paper using terahertz time-domain spectroscopy

We present a noncontact method for quantitative composition and thickness monitoring of flat sheet products using terahertz time-domain spectroscopy. We apply the method to obtain simultaneous measurement of thickness and moisture content of paper sheets. The paper is modeled as an effective medium of water mixed with fibers, and model parameters are estimated from fits to the measured transmission amplitude. We demonstrate the method on two different paper samples and obtain uncertainties that are comparable with existing sensor technology. Monte Carlo simulations indicate that these uncertainties can be reduced further by at least an order of magnitude.

Characterization of the TIP4P-Ew water model: vapor pressure and boiling point

The liquid-vapor-phase equilibrium properties of the previously developed TIP4P-Ew water model have been studied using thermodynamic integration free-energy simulation techniques in the temperature range of 274-400 K. We stress that free-energy results from simulations need to be corrected in order to be compared to the experiment. This is due to the fact that the thermodynamic end states accessible through simulations correspond to fictitious substances (classical rigid liquids and classical rigid ideal gases) while experiments operate on real substances (liquids and real gases, with quantum effects). After applying analytical corrections the vapor pressure curve obtained from simulated free-energy changes is in excellent agreement with the experimental vapor pressure curve. The boiling point of TIP4P-Ew water under ambient pressure is found to be at 370.3+/-1.9 K, about 7 K higher than the boiling point of TIP4P water (363.7+/-5.1 K; from simulations that employ finite range treatment of electrostatic and Lennard-Jones interactions). This is in contrast to the approximately +15 K by which the temperature of the density maximum and the melting temperature of TIP4P-Ew are shifted relative to TIP4P, indicating that the temperature range over which the liquid phase of TIP4P-Ew is stable is narrower than that of TIP4P and resembles more that of real water. The quality of the vapor pressure results highlights the success of TIP4P-Ew in describing the energetic and entropic aspects of intermolecular interactions in liquid water.

Terahertz direct-contact refractometer/spectrometer

A method for terahertz (THz) reflection spectroscopy is presented that utilizes direct contact between the backside of the pulsed emitter and the sample. Changes in the complex reflectivity of the emitter-sample interface are used to measure the spectrum. This method is especially useful for spectroscopy of liquid samples and soft materials. It simplifies the optical apparatus, has good sensitivity, and has capability for use in the construction of miniaturized THz probes. The problem of nonplanar wave fronts inside the emitter is discussed.

Free energy of liquid water from a computer simulation via cell theory

A method to calculate the free energy of water from computer simulation is presented. Based on cell theory, it approximates the potential energy surface sampled in the simulation by an anisotropic six-dimensional harmonic potential to model the three hindered translations and three hindered rotations of a single rigid water molecule. The potential is parametrized from the magnitude of the forces and torques measured in the simulation. The entropy of these six harmonic oscillators is calculated and summed with a conformational term to give the total entropy. Combining this with the simulation enthalpy yields the free energy. The six water models examined are TIP3P, SPC, TIP4P, SPC/E, TIP5P, and TIP4P-Ew. The results reproduce experiment well: free energies for all models are within 1.6 kJ mol(-1) and entropies are within 3.6 J K(-1) mol(-1). Approximately two-thirds of the entropy comes from translation, a third from rotation, and 5% from conformation. Vibrational frequencies match those in the experimental infrared spectrum and assist in their assignment. Intermolecular quantum effects are found to be small, with free energies for the classical oscillator lying 0.5-0.7 kJ mol(-1) higher than in the quantum case. Molecular displacements and vibrational and zero point energies are also calculated. Altogether, these results validate the harmonic oscillator as a quantitative model for the liquid state.

Probing Intermolecular Interactions in Water/Ionic Liquid Mixtures by Far-infrared Spectroscopy

Far-infrared spectra in the range from 600 to 20 cm-1 of two hydrophilic (1-ethyl-3-methylimidazolium tetrafluoroborate and 1-butyl-3-methylimidazolium tetrafluoroborate) and one hydrophobic (1-butyl-3-methylimidazolium hexafluorophosphate) ionic liquids and their mixtures with water at different concentrations are reported. Shifts of the librational water bands depending on the nature of the anion are found to be related to the strength of the interaction between the water molecules and the anions. For both hydrophilic ionic liquids, the librational band is centered around 460 cm-1, whereas for the hydrophobic ionic liquid, it is shifted to 388 cm-1, indicating less hindered rotation of single water molecules. Multivariate curve resolution, paying special attention to the spectral range from 50 to 350 cm-1, was used to investigate the presence of different species with increasing water concentration. For both hydrophilic ionic liquids, a band located at 153 cm-1 was resolved into two different contributions. A small contribution at 202 cm-1 can be attributed to intermolecular interactions between water molecules forming dimers. The major contribution (centered at 148 cm-1) corresponds to water molecules that do not bond to each other via H-bonding. It is therefore assigned to a hindered translation arising from the stretching of the hydrogen bond between BF4- anions and water molecules. Formation of water dimers in the hydrophobic ionic liquid does not occur. Furthermore, the spectral contribution of the stretching of H-bonds between water molecules and PF6- cannot be unambiguously detected, which indicates an extremely weak interaction between water molecules and this anion.

Signatures of the hydrogen bonding in the infrared bands of water

Infrared spectroscopy measurements have been completed over a wide range of frequencies allowing to measure the evolution of both intramolecular and intermolecular vibrational modes in water as a function of temperature. Emphasis is made on the high frequency OH stretching band and the so-called connectivity band that lies in the far infrared region. The substructures of the two infrared bands are analyzed in terms of different levels of connectivity of the water molecules, along the statements of the percolation model. Both band profiles appear to be related to the different degrees of connectivity of water molecules. Comparison of the data with the predictions of the percolation model shows good agreement as for the temperature evolution of liquid water. This work provides additional support to the interpretation of water bands substructures as signatures of its very specific connectivity pattern.