skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Molecular origins of bulk viscosity in liquid water
The rapid equilibrium fluctuations of water molecules are intimately connected to the rheological response; molecular motions resetting the local structure and stresses seen as flow and volume changes. In the case of water or hydrogen bonding liquids generally, the relationship is a non-trivial consideration due to strong directional interactions complicating theoretical models and necessitating clear observation of the timescale and nautre of the associated equilibrium motions. Recent work has illustrated a coincidence of timescales for short range sub-picosecond motions and the implied timescale for the shear viscosity response in liquid water. Here, neutron and light scattering methods are used to experimentally illustrate the timescale of bulk viscosity and provide a description of the associated molecular relaxation. Brillouin scattering has been used to establish the timescale of bulk viscosity; and borrowing the Maxwell approach, the ratio of the bulk viscosity, ζ , to the bulk modulus, K , yields a relaxation time, τ B , which emerges on the order of 1–2 ps in the 280 K to 303 K temperature range. Inelastic neutron scattering is subsequently used to describe the motions of water and heavy water at the molecular scale, providing both coherent and incoherent scattering data. A rotational (alternatively described as localized) motion of water protons on the 1–2 ps timescale is apparent in the incoherent scattering spectra of water, while the coherent spectra from D 2 O on the length scale of the first sharp diffraction peak, describing the microscopic density fluctuations of water, confirms the relaxation of water structure at a comparable timescale of 1–2 ps. The coincidence of these three timescales provides a mechanistic description of the bulk viscous response, with the local structure resetting due to rotational/localized motions on the order of 1–2 ps, approximately three times slower than the relaxations associated with shear viscosity. In this way we show that the shear viscous response is most closely associated with changes in water network connectivity, while the bulk viscous response is associated with local density fluctuations.  more » « less
Award ID(s):
1836556
PAR ID:
10300504
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
17
ISSN:
1463-9076
Page Range / eLocation ID:
9494 to 9502
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The effect of the network-to-molecular structural transformation with increasing phosphorus content in P x Se 100− x (30 ≤ x ≤ 67) supercooled liquids on their shear-mechanical response is investigated using oscillatory shear rheometry. While network liquids with 30 ≤ x ≤ 40 are characterized by shear relaxation via a network bond scission/renewal process, a Maxwell scaling of the storage (G′) and loss (G″) shear moduli, and a frequency-independent viscosity at low frequencies, a new relaxation process emerges in liquids with intermediate compositions (45 ≤ x ≤ 50). This process is attributed to an interconversion between network and molecular structural moieties. Predominantly molecular liquids with x ≥ 63, on the other hand, are characterized by a departure from Maxwell behavior as the storage modulus shows a linear frequency scaling G′(ω) ∼ ω over nearly the entire frequency range below the G′–G″ crossover and a nearly constant ratio of G″/G′ in the terminal region. Moreover, the dynamic viscosity of these rather fragile molecular liquids shows significant enhancement over that of network liquids at frequencies below the dynamical onset and does not reach a frequency-independent regime even at frequencies that are four orders of magnitude lower than that of the onset. Such power-law relaxation behavior of the molecular liquids is ascribed to an extremely broad distribution of relaxation timescales with the coexistence of rapid rotational motion of individual molecules and cooperative dynamics of transient molecular clusters, with the latter being significantly slower than the shear relaxation timescale. 
    more » « less
  2. Methylation occurs in a myriad of systems with protective and regulatory functions. 8-methoxypyrene-1,3,6-trisulfonate (MPTS), a methoxy derivative of a photoacid, serves as a model system to study effects of methylation on the excited state potential energy landscape. A suite of spectroscopic techniques including transient absorption, wavelength-tunable femtosecond stimulated Raman spectroscopy (FSRS), and fluorescence quantum yield measurements via steady-state electronic spectroscopy reveal the energy dissipation pathways of MPTS following photoexcitation. Various solvents enable a systematic characterization of the H-bonding interaction, viscosity, and dynamic solvation that influence the ensuing relaxation pathways. The formation of a charge-transfer state out of the Franck–Condon region occurs on the femtosecond-to-picosecond solvation timescale before encountering a rotational barrier. The rotational relaxation correlates with the H-bond donating strength of solvent, while the rotational time constant lengthens as solvent viscosity increases. Time-resolved excited-state FSRS, aided by quantum calculations, provides crucial structural dynamics knowledge and reveals the sulfonate groups playing a dominant role during solvation. Several prominent vibrational motions of the pyrene ring backbone help maneuver the population toward the more fluorescent state. These ultrafast correlated electronic and nuclear motions ultimately govern the fate of the photoexcited chromophore in solution. Overall, MPTS in water displays the highest probability to fluoresce, while the aprotic and more viscous dimethyl sulfoxide enhances the nonradiative pathways. These mechanistic insights may apply robustly to other photoexcited chromophores that do not undergo excited-state proton transfer or remain trapped in a broad electronic state and also provide design principles to control molecular optical responses with site-specific atomic substitution. 
    more » « less
  3. It is generally accepted that the limit on the stable rotation of neutron stars is set by gravitational-radiation reaction (GRR) driven instabilities, which cause the stars to emit gravitational waves that carry angular momentum away from them. The instability modes are moderated by the shear viscosity and the bulk viscosity of neutron star matter. Among the GRR instabilities, the f-mode instability plays a historically predominant role. In this work, we determine the instability periods of this mode for three different relativistic models for the nuclear equation of state (EoS) named DD2, ACB4, and GM1L. The ACB4 model for the EoS accounts for a strong first-order phase transition that predicts a new branch of compact objects known as mass-twin stars. DD2 and GM1L are relativistic mean field (RMF) models that describe the meson-baryon coupling constants to be dependent on the local baryon number density. Our results show that the f-mode instability associated with m=2 sets the limit of stable rotation for cold neutron stars (T≲1010 K) with masses between 1M⊙ and 2M⊙. This mode is excited at rotation periods between 1 and 1.4 ms (∼20% to ∼40% higher than the Kepler periods of these stars). For cold hypothetical mass-twin compact stars with masses between 1.96M⊙ and 2.10M⊙, the m=2 instability sets in at rotational stellar periods between 0.8 and 1 millisecond (i.e., ∼25% to ∼30% above the Kepler period). 
    more » « less
  4. Abstract Thermodynamical and dynamical aspects of the climate system response to an-thropogenic forcing are often considered in two distinct frameworks: The former in the context of the forcing-feedback framework; the latter in the context of eddy-mean flow feedbacks and large-scale thermodynamic constraints. Here we use experiments with the dynamical core of a general circulation model (GCM) to provide insights into the relationships between these two frameworks. We first demonstrate that the climate feedbacks and climate sensitivity in a dynamical core model are determined by its prescribed thermal relaxation timescales. We then perform two experiments: One that explores the relationships between the thermal relaxation timescale and the climatological circulation; and a second that explores the relationships between the thermal relaxation timescale and the circulation response to a global warming-like forcing perturbation. The results indicate that shorter relaxation timescales (i.e., lower climate sensitivities in the context of a dynamical core model) are associated with 1) a more vigorous large-scale circulation, including increased thermal diffusivity and stronger and more poleward storm tracks and eddy-driven jets and 2) a weaker poleward displacement of the storm tracks and eddy-driven jets in response to the global warming-like forcing perturbation. Interestingly, the circulation response to the forcing perturbation effectively disappears when the thermal relaxation timescales are spatially uniform, suggesting that the circulation response to homogeneous forcing requires spatial inhomogeneities in the local feedback parameter. Implications for anticipating the circulation response to global warming and thermodynamic constraints on the circulation are discussed. 
    more » « less
  5. The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly-stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine the coil radius of gyration, Rg, within PECs using small angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 0C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η0, versus chain length N, was η0 ~ N3.1, in agreement with “sticky reptation” theory. The lifetime and activation energy, Ep, of a pair between polyanion and polycation repeat units, Pol+Pol-, were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2Ep, showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN- acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol+Pol- pairing efficiently, which decreases effective crosslinking and decreases viscosity. An equation was derived that quantitatively predicts this effect. 
    more » « less