Rotational motion in liquid water is anisotropic: a nuclear magnetic resonance and molecular dynamics simulation study

Experimental NMR measurements of the deuterium and (17)O T(1) relaxation times in deuterium-enriched liquid water have been performed from 275 to 350 K. These relaxation times can yield rotational correlation times of appropriate molecule-fixed unit vectors if the quadrupole coupling constants and asymmetry parameters are known. We determine the latter from ab initio studies of water clusters and experimental chemical shift measurements. We find that the rotational correlation time for the OD bond vector in D(2)(16)O varies from 5.8 ps at 275 K to 0.86 ps at 350 K, and that the rotational correlation time for the out-of-plane vector of dilute D(2)(17)O in D(2)(16)O varies from 4.4 ps at 275 K to 0.64 ps at 350 K. These results indicate that the rotational motion of water is anisotropic. Molecular dynamics simulations of liquid water are in good agreement with these experiments at the higher temperatures, but the simulation results are considerably faster than experiment at the lower temperatures.     

Dynamics of dilute water in carbon tetrachloride

A dilute solution of water in a hydrophobic solvent, such as carbon tetrachloride (CCl4), presents an opportunity to study the rotational properties of water without the complicating effects of hydrogen bonds. We report here the results of theoretical, experimental, and semiempirical studies of a 0.03 mole percent solution of water in CCl4. It is shown that for this solution there are negligible water-water interactions or water-CCl4 interactions; theoretical and experimental values for proton NMR chemical shifts (deltaH) are used to confirm the minimal interactions between water and the CCl4. Calculated ab initio values and semiempirical values for oxygen-17 and deuterium quadrupole coupling constants (chi) of water/CCl4 clusters are reported. Experimental values for the 17O, 2H, and 1H NMR spin-lattice relaxation times, T1, of 0.03 mole percent water in dilute CCl4 solution at 291 K are 94+/-3 ms, 7.0+/-0.2 s, and 12.6+/-0.4 s, respectively. These T1 values for bulk water are also referenced. "Experimental" values for the quadrupole coupling constants and relaxation times are used to obtain accurate, experimental values for the rotational correlation times for two orthogonal vectors in the water molecule. The average correlation time, tauc, for the position vector of 17O (orthogonal to the plane of the molecule) in monomer water, H2(17)O, is 91 fs. The average value for the deuterium correlation time for the deuterium vector in 2H2O is 104 fs; this vector is along the OD bond. These values indicate that the motion of monomer water in CCl4 is anisotropic. At 291 K, the oxygen rotational correlation time in bulk 2H2(17)O is 2.4 ps, the deuterium rotational correlation time in the same molecule is 3.25 ps. (Ropp, J.; Lawrence, C.; Farrar, T. C.; Skinner, J. L. J. Am. Chem. Soc. 2001, 123, 8047.) These values are a factor of about 20 longer than the tauc value for dilute monomer water in CCl4.     

Machine learning and semi-empirical calculations: a synergistic approach to rapid,accurate, and mechanism-based reaction barrier prediction

Modern QM modelling methods, such as DFT, have provided detailed mechanistic insights into countless reactions. However, their computational cost inhibits their ability to rapidly screen large numbers of substrates and catalysts in reaction discovery. For a C–C bond forming nitro-Michael addition, we introduce a synergistic semi-empirical quantum mechanical (SQM) and machine learning (ML) approach that allows the prediction of DFT-quality reaction barriers in minutes, even on a standard laptop using widely available modelling software. Mean absolute errors (MAEs) are obtained that are below the accepted chemical accuracy threshold of 1 kcal mol⁻¹ and substantially better than SQM methods without ML correction (5.71 kcal mol⁻¹). Predictive power is shown to hold when the ML models are applied to an unseen set of compounds from the toxicology literature. Mechanistic insight is also achieved via the generation of full SQM transition state (TS) structures which are found to be very good approximations for the DFT-level geometries, revealing important steric interactions in some TSs. This combination of speed, accuracy, and mechanistic insight is unprecedented; current ML barrier models compromise on at least one of these important criteria.

A synergistic approach that combines machine learning with semi-empirical methods enables the fast and accurate prediction of DFT-quality reaction barriers, with mechanistic insights available from semi-empirical transition state geometries.     

Measuring the stiffness of concrete shear walls during dynamic tests

Results from a series of experiments designed to measure the stiffness of low-aspect-ratio, reinforced-concrete shear walls subjected to simulated seismic inputs on a shake table are reported. The geometry of the test structures allows them to be modeled as single-degree-of-freedom systems. Forces were estimated from accelerometer measurements on masses attached to the structures. Dynamic relative-displacement measurements were obtained from groups of strain gages wired in series to act as one continuous gage. Because this method measures relative displacements, potential sources of error associated with unspecified base motion are avoided. tiffness values determined from the relative displacement measurements were compared with stiffness values determined indirectly from frequency-response functions. Measured, accelerometer data were used to calculate the frequency-response functions. The stiffness values determined from the relative-displacement measurements gave results similar to those given by mechanics-of-materials beam theory that accounts for shear deformation. The stiffness values determined from the frequency-response functions were considerably less than those from the theory.     

Dynamics study of the reaction OH- + C2H2-->C2H- + H2O with crossed beams and density-functional theory calculations

The proton transfer reaction between OH- and C2H2, the sole reactive process observed over the collision energy range from 0.37 to 1.40 eV, has been studied using the crossed beam technique and density-functional theory (DFT) calculations. The center of mass flux distributions of the product C2H- ions at three different energies are highly asymmetric, characteristic of a direct process occurring on a time scale much less than a rotational period of any transient intermediate. The maxima in the flux distributions correspond to product velocities and directions close to those of the precursor acetylene reactants. The reaction quantitatively transforms the entire exothermicity into internal excitation of the products, consistent with an energy release motif in which the proton is transferred early, in a configuration in which the forming bond is extended. This picture is supported by DFT calculations showing that the first electrostatically bound intermediate on the reaction pathway is the productlike C2H- H2O species. Most of the incremental translational energy in the two higher collision energy experiments appears in product translational energy, and provides an example of induced repulsive energy release characteristic of the heavy+light-heavy mass combination.     

Reaction dynamics of OH+(3Sigma-)+C2H2 studied with crossed beams and density functional theory calculations

The reactions between OH+(3Sigma-) and C2H2 have been studied using crossed ion and molecular beams and density functional theory calculations. Both charge transfer and proton transfer channels are observed. Products formed by carbon-carbon bond cleavage analogous to those formed in the isoelectronic O(3P)+C2H2 reaction, e.g., 3CH2 + HCO+, are not observed. The center of mass flux distributions of both product ions at three different energies are highly asymmetric, with maxima close to the velocity and direction of the precursor acetylene beam, characteristic of direct reactions. The internal energy distributions of the charge transfer products are independent of collision energy and are peaked at the reaction exothermicity, inconsistent with either the existence of favorable Franck-Condon factors or energy resonance. In proton transfer, almost the entire reaction exothermicity is transformed into product internal excitation, consistent with mixed energy release in which the proton is transferred with both the breaking and forming bonds extended. Most of the incremental translational energy in the two higher-energy experiments appears in product translational energy, providing an example of induced repulsive energy release.     

Dynamic nuclear polarization at 9T using a novel 250GHz gyrotron microwave source

In this communication, we report enhancements of nuclear spin polarization by dynamic nuclear polarization (DNP) in static and spinning solids at a magnetic field strength of 9T (250 GHz for g=2 electrons, 380 MHz for 1H). In these experiments, 1H enhancements of up to 170+/-50 have been observed in 1-13C-glycine dispersed in a 60:40 glycerol/water matrix at temperatures of 20K; in addition, we have observed significant enhancements in 15N spectra of unoriented pf1-bacteriophage. Finally, enhancements of approximately 17 have been obtained in two-dimensional 13C-13C chemical shift correlation spectra of the amino acid U-13C, 15N-proline during magic angle spinning (MAS), demonstrating the stability of the DNP experiment for sustained acquisition and for quantitative experiments incorporating dipolar recoupling. In all cases, we have exploited the thermal mixing DNP mechanism with the nitroxide radical 4-amino-TEMPO as the paramagnetic dopant. These are the highest frequency DNP experiments performed to date and indicate that significant signal enhancements can be realized using the thermal mixing mechanism even at elevated magnetic fields. In large measure, this is due to the high microwave power output of the 250 GHz gyrotron oscillator used in these experiments.     

Surface tailoring for controlled protein adsorption: effect of topography at the nanometer scale and chemistry

Protein adsorption behavior is at the heart of many of today's research fields including biotechnology and materials science. With understanding of protein-surface interactions, control over the conformation and orientation of immobilized species may ultimately allow tailor-made surfaces to be generated. In this contribution protein-surface interactions have been examined with particular focus on surface curvature with and without surface chemistry effects. Silica spheres with diameters in the range 15-165 nm with both hydrophilic and hydrophobic surface chemistries have been used as model substrates. Two proteins differing in size and shape, bovine serum albumin (BSA) and bovine fibrinogen (Fg), have been used in model studies of protein binding with detailed secondary structure analysis being performed using infrared spectroscopy (IR) on surface-bound proteins. Although trends in binding affinity and saturation values were similar for both proteins, albumin is increasingly less ordered on larger substrates, while fibrinogen, in contrast, loses secondary structure to a greater extent when adsorbing onto particles with high surface curvature. These effects are compounded by surface chemistry, with both proteins becoming more denatured on hydrophobic surfaces. Both surface chemistry and topography play key roles in determining the structure of the bound proteins. A model of the binding characteristics of these two proteins onto surfaces having differing curvature and chemistry is presented. We propose that properties of an adsorbed protein layer may be guided through careful consideration of surface structure, allowing the fabrication of materials/surface coatings with tailored bioactivity.