•A long-standing solvation dynamics issue is addressed via novel theoretical methods.•Molecular water solvent response to solute charge extinction and creation is studied.•The quantitatively unexplained response asymmetry leads to nonlinear response.•Asymmetric response origin is found via energy flows and molecular configurations.•Water molecule translations, and not the faster water rotations, are responsible.Early molecular dynamics simulations discovered an important asymmetry in the speed of water solvation dynamics for charge extinction and charge creation for an immersed solute, a feature representing a first demonstration of the breakdown of linear response theory. The molecular level mechanism of this asymmetry is examined here via a novel energy flux theoretical approach coupled to geometric probes. The results identify the effect as arising from the translational motions of the solute-hydrating water molecules rather than their rotational/librational motions, even though the latter are more rapid and dominate the energy flow.Download high-res image (105KB)Download full-size image

In previous installments it has been shown how a detailed analysis of energy fluxes induced by electronic excitation of a solute can provide a quantitative understanding of the dominant molecular energy flow channels characterizing solvation—and in particular, hydration— relaxation dynamics. Here this work and power approach is complemented with a detailed characterization of the changes induced by such energy fluxes. We first examine the water solvent’s spatial and orientational distributions and the assorted energy fluxes in the various hydration shells of the solute to provide a molecular picture of the relaxation. The latter analysis is also used to address the issue of a possible “inverse snowball” effect, an ansatz concerning the time scales of the different hydration shells to reach equilibrium. We then establish a link between the instantaneous torque, exerted on the water solvent neighbors’ principal rotational axes immediately after excitation and the final energy transferred into those librational motions, which are the dominant short-time energy receptor.

Carbonic, lactic, and pyruvic acids have been generated in aqueous solution by the transient protonation of their corresponding conjugate bases by a tailor-made photoacid, the 6-hydroxy-1-sulfonate pyrene sodium salt molecule. A particular goal is to establish the pKa of carbonic acid H2CO3. The on-contact proton transfer (PT) reaction rate from the optically excited photoacid to the carboxylic bases was derived, with unprecedented precision, from time-correlated single-photon-counting measurements of the fluorescence lifetime of the photoacid in the presence of the proton acceptors. The time-dependent diffusion-assisted PT rate was analyzed using the Szabo–Collins–Kimball equation with a radiation boundary condition. The on-contact PT rates were found to follow the acidity order of the carboxylic acids: the stronger was the acid, the slower was the PT reaction to its conjugate base. The pKa of carbonic acid was found to be 3.49 ± 0.05 using both the Marcus and Kiefer–Hynes free energy correlations. This establishes H2CO3 as being 0.37 pKa units stronger and about 1 pKa unit weaker, respectively, than the physiologically important lactic and pyruvic acids. The considerable acid strength of intact carbonic acid indicates that it is an important protonation agent under physiological conditions.

The protonation of methylamine base CH3NH2 by carbonic acid H2CO3 within a hydrogen (H)-bonded complex in aqueous solution was studied via Car–Parrinello dynamics in the preceding paper (Daschakraborty, S.; Kiefer, P. M.; Miller, Y.; Motro, Y.; Pines, D.; Pines, E.; Hynes, J. T. J. Phys. Chem. B 2016, DOI: 10.1021/acs.jpcb.5b12742). Here some important further details of the reaction path are presented, with specific emphasis on the water solvent’s role. The overall reaction is barrierless and very rapid, on an ∼100 fs time scale, with the proton transfer (PT) event itself being very sudden (<10 fs). This transfer is preceded by the acid–base H-bond’s compression, while the water solvent changes little until the actual PT occurrence; this results from the very strong driving force for the reaction, as indicated by the very favorable acid-protonated base ΔpKa difference. Further solvent rearrangement follows immediately the sudden PT’s production of an incipient contact ion pair, stabilizing it by establishment of equilibrium solvation. The solvent water’s short time scale ∼120 fs response to the incipient ion pair formation is primarily associated with librational modes and H-bond compression of water molecules around the carboxylate anion and the protonated base. This is consistent with this stabilization involving significant increase in H-bonding of hydration shell waters to the negatively charged carboxylate group oxygens’ (especially the former H2CO3 donor oxygen) and the nitrogen of the positively charged protonated base’s NH3+.

Protonation by carbonic acid H2CO3 of the strong base methylamine CH3NH2 in a neutral contact pair in aqueous solution is followed via Car–Parrinello molecular dynamics simulations. Proton transfer (PT) occurs to form an aqueous solvent-stabilized contact ion pair within 100 fs, a fast time scale associated with the compression of the acid–base hydrogen-bond (H-bond), a key reaction coordinate. This rapid barrierless PT is consistent with the carbonic acid-protonated base pKa difference that considerably favors the PT, and supports the view of intact carbonic acid as potentially important proton donor in assorted biological and environmental contexts. The charge redistribution within the H-bonded complex during PT supports a Mulliken picture of charge transfer from the nitrogen base to carbonic acid without altering the transferring hydrogen’s charge from approximately midway between that of a hydrogen atom and that of a proton.

Infrared spectroscopy has been used to characterize the solvent effect on the OH stretching vibrations νOH of phenol, 1-naphthol, 2-naphthol, 1-hydroxypyrene, and ethanol. We distinguish the dielectric (nonspecific) effect of the solvent on ΔνOH, the observed red-shifts in νOH, from the much larger red-shift caused by direct hydrogen (H)-bonding interactions with the solvents. To isolate the solvent dielectric constant ε effect on νOH, the OH oscillator was also studied when it is already H-bonded with an invariant oxygen base, dimethyl sulfoxide. We find that ΔνOH depends importantly on ΔPA, the difference between the proton affinities of the conjugate base of the proton donor and the proton acceptor. For a given H-bonded complex, νOH tends to vary inversely with ε, exhibiting different slopes for polar and nonpolar solvents, i.e., solvents comprising molecules with and without a permanent dipole moment, respectively. We use a two-state valence-bond-based theory to analyze our experimental data. This demonstrates that the OH oscillator acquires a more ionic-like character in the vibrational excited state, i.e., charge transfer; this results in a stronger H-bond in a more anharmonic potential for the OH vibration. The theory distinguishes between nonpolar and polar solvents and successfully accounts for the observed 1/ε and ΔPA variations.

Infrared spectroscopy measurements were used to characterize the OH stretching vibrations in a series of similarly structured fluoroethanols, RCH2OH (R = CH3, CH2F, CHF2, CF3), a series which exhibits a systematic increase in the molecule acidity with increasing number of F atoms. This study, which expands our earlier efforts, was carried out in non-hydrogen-bonding solvents comprising molecules with and without a permanent dipole moment, with the former solvents being classified as polar solvents and the latter designated as nonpolar. The hydrogen bond interaction in donor–acceptor complexes formed in solution between the fluorinated ethanol H-donors and the H-acceptor base DMSO was investigated in relation to the solvent dielectric and to the differences ΔPA of the gas phase proton affinities (PAs) of the conjugate base of the fluorinated alcohols and DMSO. We have observed that νOH decreases as the acidity of the alcohol increases (ΔPA decreases) and that νOH varies inversely with ε, exhibiting different slopes for nonpolar and polar solvents. These 1/ε slopes tend to vary linearly with ΔPA, increasing with increasing acidity. These experimental findings, including the ΔPA trends, are described with our recently published two-state Valence Bond-based theory for acid–base H-bonded complexes. Lastly, the correlation of the alcohol’s conjugate base PAs with Taft σ* values of the fluorinated ethyl groups CHnF3–nCH2– provides a connection of the inductive effects for these groups with the acidity parameter ΔPA associated with the H-bonded complexes.

A theory is presented for the proton stretch vibrational frequency νAH for hydrogen (H−) bonded complexes of the acid dissociation type, that is, AH···B ⇔ A–···HB+(but without complete proton transfer), in both polar and nonpolar solvents, with special attention given to the variation of νAH with the solvent’s dielectric constant ε. The theory involves a valence bond (VB) model for the complex’s electronic structure, quantization of the complex’s proton and H-bond motions, and a solvent coordinate accounting for nonequilibrium solvation. A general prediction is that νAH decreases with increasing ε largely due to increased solvent stabilization of the ionic VB structure A–···HB+ relative to the neutral VB structure AH···B. Theoretical νAH versus 1/ε slope expressions are derived; these differ for polar and nonpolar solvents and allow analysis of the solvent dependence of νAH. The theory predicts that both polar and nonpolar slopes are determined by (i) a structure factor reflecting the complex’s size/geometry, (ii) the complex’s dipole moment in the ground vibrational state, and (iii) the dipole moment change in the transition, which especially reflects charge transfer and the solution phase proton potential shapes. The experimental proton frequency solvent dependence for several OH···O H-bonded complexes is successfully accounted for and analyzed with the theory.

We report a study of DNA deformations using a coarse-grained mechanical model and quantitatively interpret the allosteric effects in protein–DNA binding affinity. A recent single-molecule study (Kim et al. Science 2013, 339, 816) showed that when a DNA molecule is deformed by specific binding of a protein, the binding affinity of a second protein separated from the first protein is altered. Experimental observations together with molecular dynamics simulations suggested that the origin of the DNA allostery is related to the observed deformation of DNA’s structure, in particular, the major groove width. To unveil and quantify the underlying mechanism for the observed major groove deformation behavior related to the DNA allostery, here we provide a simple but effective analytical model where DNA deformations upon protein binding are analyzed and spatial correlations of local deformations along the DNA are examined. The deformation of the DNA base orientations, which directly affect the major groove width, is found in both an analytical derivation and coarse-grained Monte Carlo simulations. This deformation oscillates with a period of 10 base pairs with an amplitude decaying exponentially from the binding site with a decay length lD ≈10 base pairs as a result of the balance between two competing terms in DNA base-stacking energy. This length scale is in agreement with that reported from the single-molecule experiment. Our model can be reduced to the worm-like chain form at length scales larger than lP but is able to explain DNA’s mechanical properties on shorter length scales, in particular, the DNA allostery of protein–DNA interactions.

The ultrafast librational (hindered rotational) relaxation of a rotationally excited H2O molecule in pure liquid water is investigated by means of classical nonequilibrium molecular dynamics simulations and a power and work analysis. This analysis allows the mechanism of the energy transfer from the excited H2O to its water neighbors, which occurs on a sub-100 fs time scale, to be followed in molecular detail, i.e., to determine which water molecules receive the energy and in which degrees of freedom. It is found that the dominant energy flow is to the four hydrogen-bonded water partners in the first hydration shell, dominated by those partners’ rotational motion, in a fairly symmetric fashion over the hydration shell. The minority component of the energy transfer, to these neighboring waters’ translational motion, exhibits an asymmetry in energy reception between hydrogen-bond-donating and -accepting water molecules. The variation of the energy flow characteristics with rotational axis, initial rotational energy excitation magnitude, method of excitation, and temperature is discussed. Finally, the relation of the nonequilibrium results to equilibrium time correlations is investigated.

The first, key step of water oxidation catalysis by the ruthenium blue dimer transition metal complex has been studied via density functional methods and with extensive explicit solvation, starting from the oxidized catalytically active form of the dimer. This step is the rate-limiting O–O single bond formation. This reaction is found to involve several proton transfers through a proton relay chain, synergetically coupled to electron flow through the μ-oxo bridge of the dimer. The barrier for the O–O formation step is found to arise primarily from the surrounding aqueous solvent, suggesting that it might be substantially lowered in suitable environments. Some remarks are given concerning the following, penultimate step prior to the formation of dioxygen and the stable form of the dimer, in which it is suggested that another proton relay chain is at play.

Here we briefly review our theoretical work on nitric acid dissociation at an aqueous interface, an issue important in a wide range of atmospheric contexts, employing both quantum chemical methods and Car-Parrinello molecular dynamics simulations. A first aspect concerns whether or not this dissociation, via proton transfer to a neighboring water molecule, occurs as a function of the acid’s location atop and at various depths below the surface. Further aspects concern the molecular level details of the mechanism, i.e. which motions are important, of the first proton transfer to form a contact ion pair, and of the subsequent proton transfer to produce a solvent-separated ion pair from that contact ion pair.