Reorientational Correlation Times Research Articles

Reorientation of DABCO (1,4-diazabicyclo[2.2.2]octane) in hydrogen-bonded 1:2 molecular complexes DABCO-2(SBA) with substituted benzoic acids SBA

The molecular complexes of DABCO-2(SBA) with substituted benzoic acids SBA, which have the hydrogen-bonded trimeric unit, were prepared to investigate the DABCO reorientation in the trimeric unit. By the analysis of the temperature dependence of 1H NMR spin-lattice relaxation time, the correlation time of the DABCO reorientation was determined as a function of the temperature. It is suggested that a change of potential energy of the reorientation is caused with increasing rate of reorientation resulting in a structure change of the trimeric unit to a centrosymmetric structure. The correlation between the activation energy of the DABCO reorientation and the bond nature of the hydrogen bond involved in the trimeric unit is also discussed.

Temperature-Dependent Dielectric Relaxation in Ionic Acetamide Deep Eutectics: Partial Viscosity Decoupling and Explanations from the Simulated Single-Particle Reorientation Dynamics and Hydrogen-Bond Fluctuations.

We report here temperature-dependent (293 ≤ T (K) ≤ 336) dielectric relaxation (DR) measurements of (acetamide + LiBr/NO3-/ClO4-) deep eutectic solvents (DESs) in the frequency window of 0.2 ≤ ν (GHz) ≤ 50 and explore, via molecular dynamics simulations, the relative roles for the collective single-particle reorientational relaxations and the H-bond dynamics of acetamide in the measured DR response. In addition, DR measurements of neat molten acetamide were performed. Recorded DR spectra of these DESs require multi-Debye fits and produce well-separated DR time scales that are spread over several picoseconds to ∼1 ns. Simulations suggest DR time scales derive contributions from both the collective reorientational () relaxation and structural H-bond (CHB(t)) dynamics of acetamide. A good correlation between the measured and simulated activation energies further reveals a strong connection between the measured DR and the simulated and CHB(t). Average DR times exhibit a strong fractional viscosity dependence, suggesting substantial microheterogeneity in these media. Simulations of and CHB(t) reveal strong stretched exponential relaxations with a stretching exponent, 0.4 ≤ β ≤ 0.7. The ratio between the average reorientational correlation times of first and second ranks, , deviates appreciably from Debye's law for homogeneous media. Importantly, a pronounced translation-rotation decoupling between the simulated reorientation and center-of-mass diffusion times was observed.

The Water Polymorphism and the Liquid–Liquid Transition from Transport Data

NMR spectroscopic literature data are used, in a wide temperature-pressure range (180–350 K and 0.1–400 MPa), to study the water polymorphism and the validity of the liquid–liquid transition (LLT) hypothesis. We have considered the self-diffusion coefficient DS and the reorientational correlation time τθ (obtained from spin-lattice T1 relaxation times), measured, respectively, in bulk and emulsion liquid water from the stable to well inside the metastable supercooled region. As an effect of the hydrogen bond (HB) networking, the isobars of both these transport functions evolve with T by changing by several orders of magnitude, whereas their pressure dependence become more and more pronounced at lower temperatures. Both these transport functions were then studied according to the Adam–Gibbs model, typical of glass forming liquids, obtaining the water configurational entropy and the corresponding specific heat contribution. The comparison of the evaluated CP,conf isobars with the experimentally measured water specific heat reveals the full consistency of this analysis. In particular, the observed CP,conf maxima and its diverging behaviors clearly reveals the presence of the LLT and with a reasonable approximation the liquid–liquid critical point (LLCP) locus in the phase diagram.

Heterogeneous Orientational Relaxations and Translation-Rotation Decoupling in (Choline Chloride + Urea) Deep Eutectic Solvents: Investigation through Molecular Dynamics Simulations and Dielectric Relaxation Measurements.

Molecular dynamics simulations and dielectric relaxation (DR) measurements in the frequency window, 0.2 ≤ ν/GHz ≤ 50, have been performed to explore the heterogeneous reorientation dynamics in [f choline chloride + (1 - f) urea] deep eutectic solvents (DESs) at f = 0.33 and 0.40 in the temperature range 293 ≤ T/K ≤ 333. The solution viscosity varies by more than an order of magnitude. DR measurements in these DESs reveal multiple relaxation timescales-τ1 ∼ 500 ps, τ2 ∼ 100 ps, τ3 ∼ 30 ps, and τ4 ∼ 5 ps. Simulated rank-dependent collective single-particle reorientational (Cl(t), l being the rank) and structural H-bond [CHB(t)] relaxations can explain the microscopic origin of all these DR timescales. The average DR times, ⟨τDR⟩, exhibit a pronounced fractional viscosity dependence, ⟨τDR⟩ ∝ (η/T)p, with p = 0.1. This experimental evidence of pronounced heterogeneous reorientation dynamics in these DESs is supported by a strong translation-rotation decoupling and a significant deviation of the average reorientational correlation times (⟨τl⟩) from Debye's l(l + 1) law. The simulated ratios between the average rotation and translation timescales for both urea and choline correctly reduce to the appropriate hydrodynamic limit at high temperatures. The stretched exponential relaxations of the simulated self-dynamic structure factors and the non-Gaussian single-particle displacement distributions further support strong temporal heterogeneity in these DESs. Dynamic susceptibilities from the simulated four-point correlations exhibit long correlated timescales. Moreover, simulated activation energies estimated from the temperature-dependent C1(t) decays and the translational diffusion coefficients from the velocity autocorrelation functions agree favorably with those from the corresponding DR and the pulsed field gradient nuclear magnetic resonance measurements.

Symmetrized Normal Mode Analysis of the Spin–Lattice Relaxation in Solid Calcium Borohydride

The symmetrized normal mode analysis presented in a previous paper was used to analyze the 11B spin–lattice relaxation of solid calcium borohydride. From the analytical results, the quadrupolar and the dipolar spectral densities and reorientational correlation times of various relaxation interactions were evaluated. Based on the spectral densities, the quadrupolar contribution to the spin–lattice relaxation of 11B nucleus was calculated. As expected from the QCC value of 11B nucleus in solid calcium borohydride, the ratio of the quadrupolar contribution to the dipole–dipole contribution in the 11B spin–lattice relaxation was proved to be higher than 30% at 100°C. The symmetrized normal mode could be a useful method for the interpretation of AX4 spin system containing a quadrupolar nucleus, and thus the determination of the motional parameters of BH4 unit in solid calcium borohydride.

Supercooled nano-droplets of water confined in hydrophobic rubber.

Hydrophobic elastomers are capable of absorbing a small amount of water that forms droplets around hydrophilic sites. These systems allow the study of confinement effects by a hydrophobic environment on the dynamics and thermodynamic behaviour of water molecules. The freezing-melting properties and the dynamics of water inside nano-droplets in butyl rubber are affected, as revealed by differential scanning calorimetry (DSC) and deuterium nuclear magnetic resonance (2H-NMR). Upon cooling down, all water crystalizes with a bimodal droplet population (da = 3.4 nm and db = 4.4 nm) in a temperature range associated with the droplet size distribution. However, the melting temperature is not shifted according to the Gibbs-Thomson equation. The relative decrease of the 2H-NMR longitudinal magnetization is not a single exponential and, by inverse Laplace transformation, it was deduced to be bimodal in agreement with the DSC measurements (T1,a ∼ 10 ms and T1,b ∼ 200 ms). The deduced correlation time of molecular reorientation is longer than that of bulk water and the behaviour with temperature follows the Vogel-Fulcher-Tammann (VFT) equations with a changing fragility as the droplet size is reduced when reducing hydration.

Hydration Structure and Dynamics of the Favipiravir Antiviral Drug: A Molecular Modelling Approach

Abstract The hydration structure of the Favipiravir antiviral drug, at infinite dilution in water, was investigated by employing a systematic molecular modelling approach. An effective interaction potential model was employed for Favipiravir, using the intramolecular geometry and charge distribution from quantum chemical calculations performed in the present treatment and adopting well-established Lennard-Jones parameters. The hydration structure and related dynamics were further investigated by means of classical molecular dynamics simulations. These calculations have revealed the existence of different types of hydrogen bonds between Favipiravir and the surrounding water molecules, with continuous lifetimes in the sub picosecond range and intermittent lifetimes in the range of 0.8–5.4 ps. The self-diffusion coefficient of Favipiravir at 298.15 K was found to be three times lower than the value obtained for water in solution, while comparable to the values measured for other common painkillers, anti-inflammatory drugs, antibiotics and corticosteroids for asthma treatment. It was also revealed that the rotational motions of Favipiravir are more retarded in comparison with water and this is reflected in the calculated reorientational correlation times of specific intramolecular vectors. The results obtained could be useful for further pharmacokinetic and computer-aided docking studies to evaluate the efficiency of this antiviral drug.

Dynamic Properties of Water Confined in Graphene-Based Membrane: A Classical Molecular Dynamics Simulation Study.

We performed molecular dynamics simulations of water molecules inside a hydrophobic membrane composed of stacked graphene sheets. By decreasing the density of water molecules inside the membrane, we observed that water molecules form a droplet through a hydrogen bond with each other in the hydrophobic environment that stacked graphene sheets create. We found that the water droplet translates as a whole body rather than a dissipate. The translational diffusion coefficient along the graphene surface increases as the number of water molecules in the droplet decreases, because the bigger water droplet has a stronger van der Waals interaction with the graphene surface that hampers the translational motion. We also observed a longer hydrogen bond lifetime as the density of water decreased, because the hydrophobic environment limits the libration motion of the water molecules. We also calculated the reorientational correlation time of the water molecules, and we found that the rotational motion of confined water inside the membrane is anisotropic and the reorientational correlation time of confined water is slower than that of bulk water. In addition, we employed steered molecular dynamics simulations for guiding the target molecule, and measured the free energy profile of water and ion penetration through the interstice between graphene sheets. The free energy profile of penetration revealed that the optimum interlayer distance for desalination is ~10 Å, where the minimum distance for water penetration is 7 Å. With a 7 Å interlayer distance between the graphene sheets, water molecules are stabilized inside the interlayer space because of the van der Waals interaction with the graphene sheets where sodium and chloride ions suffer from a 3–8 kcal/mol energy barrier for penetration. We believe that our simulation results would be a significant contribution for designing a new graphene-based membrane for desalination.

Translational and Rotational Diffusion in Liquid Water at Very High Pressure: A Simulation Study.

Translational and rotational diffusion coefficients of liquid water have been computed from molecular dynamics simulation with a recent polarizable potential at 298, 400, and 550 K at very high pressure. At 298 K, the model reproduces the initial increase and the occurrence of a maximum for the translational and rotational diffusion coefficients when the pressure increases. At 400 and 550 K, translational and rotational diffusion coefficients are found to monotonically decrease when pressure increases in the gigapascal range, with the translational coefficient decreasing faster than the rotational one. At 400 K, such an evolution of the rotational diffusion coefficient contrasts with quasielastic neutron scattering results predicting a near independence of the rotational diffusion with a pressure increase above ≃0.5 GPa. An interpretation is proposed to explain this discrepancy. The pressure dependence of the translation-rotation coupling is analyzed. The anisotropy of rotational diffusion is investigated by computing the rotational diffusion tensor in a molecular system of axes and the reorientational correlation times of rank 1 and rank 2 of the inertia axes and of the OH bond vector. Deviation of the simulation data with respect to the predictions of the isotropic Debye model of rotational diffusion are quantified and can be used to estimate experimental rotational diffusion coefficients from experimental reorientational correlation times.

Translational and rotational diffusion of rod shaped molecules by molecular dynamics simulations

The results of molecular dynamics simulations of the dynamical evolution of assemblies of linear rigid rods of variable aspect ratio, a, and number density, ρ, in the isotropic phase are reported. The rods consist of m equally spaced sites interacting with the Weeks-Chandler-Andersen repulsive pair potential, where 2 < m < 16. With increasing m, features specific to long rods, such as anisotropic self-diffusion, become apparent. There is also an increasing separation between the characteristic relaxation times of the torque, angular velocity, and reorientational time correlation functions with increasing density. The latter is exponential at high densities even for dimers. The isotropic translational diffusion coefficient, Di, and rotational diffusion coefficient, Dr, are reported as a function of m and ρ or volume fraction, ξ. The mDi data scale with ξ throughout much of the simulated range, while the rotational diffusion coefficients scale approximately as m3Dr against ρ at low densities but as ∼m6Dr at high ρ, consistent with theories of colloidal and noncolloidal rod-containing liquids. The crossover density between the two regimes is parameterized in analytic form. The probability distribution functions for displacements and angular jumps in a given time show evidence of non-Gaussian behavior with increasing density. The shear viscosity and Di scale approximately as m and m-1, respectively, in the semidilute regime, which is consistent with a Stokes-Einstein-like relationship. At high concentrations, a frustrated or glassy structure formed in which the rods were randomly oriented.

Solution dynamics of the natural bioactive molecule capsaicin: a relaxation study

Nuclear magnetic resonance spectroscopy is a straightforward technique for studying molecular dynamics that range in timescale from picosecond (motions faster than molecular reorientation) to those that occur in real-time. This approach is important to highlight the behavior of bioactive molecules in solution, and to acquire information about action mechanisms and potential pharmacological effects. Proton and carbon-13 spin–lattice relaxation experiments were performed to calculate the reorientational correlation time for protonated carbons. Capsaicin showed complex dynamical properties and the results revealed two regions with different dynamical properties: the aliphatic region with fast reorientation motions and the aromatic region with slow motions.

Simultaneous determination of deuteron quadrupole coupling constants and rotational correlation times: the model case of hydrogen bonded ionic liquids.

We show that deuteron quadrupole coupling constants (DQCCs), and reorientational correlation times of molecular bonds N-D that are involved in hydrogen bonding, can be determined from NMR T1 relaxation time experiments simultaneously. For this purpose, we used trialkylammonium-based protic ionic liquids (PILs) as model compounds. They exhibit high viscosities and wide liquid ranges that allow measurements far beyond the extreme narrowing region (ω0τc ≪ 1). The T1 minima already occur at temperatures significantly above room temperature. We obtain reasonable DQCCs for the liquid phase if anisotropic motion is considered. The DQCCs are very small due to attractive Coulomb interaction between the cation and anion, which is further enhanced by hydrogen bonding. The DQCCs strongly depend on the interaction strength of the anion but are independent of the alkyl chain length of the trialkyl ammonium cations pointing to the exclusive cation-anion interaction along the hydrogen bond.

Rotational and translational dynamics and their relation to hydrogen bond lifetimes in an ionic liquid by means of NMR relaxation time experiments and molecular dynamics simulation.

We report a concerted theoretical and experimental effort to determine the reorientational dynamics as well as hydrogen bond lifetimes for the doubly ionic hydrogen bond +OH⋯O- in the ionic liquid (2-hydroxyethyl)trimethylammonium bis(trifluoromethylsulfonyl)imide [Ch][NTf2] by using a combination of NMR relaxation time experiments, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. Due to fast proton exchange, the determination of rotational correlation times is challenging. For molecular liquids, 17O-enhanced proton relaxation time experiments have been used to determine the rotational correlation times for the OH vectors in water or alcohols. As an alternative to those expensive isotopic substitution experiments, we employed a recently introduced approach which is providing access to the rotational dynamics from a single NMR deuteron quadrupolar relaxation time experiment. Here, the deuteron quadrupole coupling constants (DQCCs) are obtained from a relation between the DQCC and the δ 1H proton chemical shifts determined from a set of DFT calculated clusters in combination with experimentally determined proton chemical shifts. The NMR-obtained rotational correlation times were compared to those obtained from MD simulations and then related to viscosities for testing the applicability of popular hydrodynamic models. In addition, hydrogen bond lifetimes were derived, using hydrogen bond population correlation functions computed from MD simulations. Here, two different time domains were observed: The short-time contributions to the hydrogen lifetimes and the reorientational correlation times have roughly the same size and are located in the picosecond range, whereas the long-time contributions decay with relaxation times in the nanosecond regime and are related to rather slow diffusion processes. The computed average hydrogen bond lifetime is dominated by the long-time process, highlighting the importance and longevity of hydrogen-bonded ion pairs in these ionic liquids.

Orientational dynamics in a room temperature ionic liquid: Are angular jumps predominant?

Reorientational dynamics of the constituent ions in a room temperature ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]), are explored via molecular dynamics simulations, and several features of orientation dynamics are summarized. The anion, [PF6]-, not only exhibits a higher propensity to orientation jumps than the cation, [BMIM]+ but also accesses a wider jump angle distribution and larger peak-angle. Jump and waiting time distributions for both the ions depict power-law dependences, suggesting temporally heterogeneous dynamics for the medium. This heterogeneity feature is further highlighted by the finding that the simulated first rank (ℓ = 1) and second rank (ℓ = 2) average reorientational correlation times reflect a severe break-down of Debye's ℓ(ℓ + 1) law for orientational diffusion in an isotropic homogeneous medium. Simulated average H-bond lifetime resides between the mean orientation jump and waiting times, while the structural H-bond relaxation suggests, as in normal liquids, a pronounced presence of translational motion of the partnering ions. Average simulated jump trajectories reveal a strong rotation-translation coupling and indicate relatively larger changes in spatial and angular arrangements for the anion during an orientation jump. In fact, a closer inspection of all these results points toward more heterogeneous dynamics for [PF6]- than [BMIM]+. This is a new observation and may simply be linked to the ion-size. However, such a generalization warrants further study.

Revisiting imidazolium based ionic liquids: Effect of the conformation bias of the [NTf2] anion studied by molecular dynamics simulations

We study ionic liquids composed of 1-alkyl-3-methylimidazolium cations and bis(trifluoromethyl-sulfonyl)imide anions ([C n MIm][NTf2]) with varying chain-length n = 2, 4, 6, 8 by using molecular dynamics simulations. We show that a reparametrization of the dihedral potentials as well as charges of the [NTf2] anion leads to an improvement of the force field model introduced by Köddermann, Paschek, and Ludwig [ChemPhysChem 8, 2464 (2007)] (KPL-force field). A crucial advantage of the new parameter set is that the minimum energy conformations of the anion (trans and gauche), as deduced from ab initio calculations and Raman experiments, are now both well represented by our model. In addition, the results for [C n MIm][NTf2] show that this modification leads to an even better agreement between experiment and molecular dynamics simulation as demonstrated for densities, diffusion coefficients, vaporization enthalpies, reorientational correlation times, and viscosities. Even though we focused on a better representation of the anion conformation, also the alkyl chain-length dependence of the cation behaves closer to the experiment. We strongly encourage to use the new NGOLP (Neumann, Golub, Odebrecht, Ludwig, Paschek) force field for the [NTf2] anion instead of the earlier KPL parameter set for computer simulations aiming to describe the thermodynamics, dynamics, and also structure of imidazolium-based ionic liquids.

Dynamical transition in molecular glasses and proteins observed by spin relaxation of nitroxide spin probes and labels.

In glassy substances and biological media, dynamical transitions are observed in neutron scattering that manifests itself as deviations of the translational mean-squared displacement, 〈x2〉, of hydrogen atoms from harmonic dynamics. In biological media, the deviation occurs at two temperature intervals, at ∼100-150 K and at ∼170-230 K, and it is attributed to the motion of methyl groups in the former case and to the transition from harmonic to anharmonic or diffusive motions in the latter case. In this work, electron spin echo (ESE) spectroscopy-a pulsed version of electron paramagnetic resonance-is applied to study the spin relaxation of nitroxide spin probes and labels introduced in molecular glass former o-terphenyl and in protein lysozyme. The anisotropic contribution to the rate of the two-pulse ESE decay, ΔW, is induced by spin relaxation appearing because of restricted orientational stochastic molecular motion; it is proportional to 〈α2〉τc, where 〈α2〉 is the mean-squared angle of reorientation of the nitroxide molecule around the equilibrium position and τc is the correlation time of reorientation. The ESE time window allows us to study motions with τc < 10-7 s. For glassy o-terphenyl, the 〈α2〉τc temperature dependence shows a transition near 240 K, which is in agreement with the literature data on 〈x2〉. For spin probes of essentially different size, the obtained data were found to be close, which evidences that motion is cooperative, involving a nanocluster of several neighboring molecules. For the dry lysozyme, the 〈α2〉τc values below 260 K were found to linearly depend on the temperature in the same way as it was observed in neutron scattering for 〈x2〉. As spin relaxation is influenced only by stochastic motion, the harmonic motions seen in ESE must be overdamped. In the hydrated lysozyme, ESE data show transitions near 130 K for all nitroxides, near 160 K for the probe located in the hydration layer, and near 180 K for the label in the protein interior. For this system, the two latter transitions are not observed in neutron scattering. The ESE-detected transitions are suggested to be related with water dynamics in the nearest hydration shell: with water glass transition near 130 K and with the onset of overall water molecular reorientations near 180 K; the disagreement with neutron scattering is ascribed to the larger time window for ESE-detected motions.

Atomistic Insight into Tetraalkylphosphonium Bis(oxalato)borate Ionic Liquid/Water Mixtures. 2. Volumetric and Dynamic Properties.

Atomistic molecular dynamics simulations have been performed to investigate volumetric quantities and dynamic properties of binary trihexyltetradecylphosphonium bis(oxalato)borate ([P6,6,6,14][BOB]) ionic liquid (IL)/water mixtures with different water concentrations. The predicted liquid densities for typical [P6,6,6,14][BOB] IL/water mixtures are consistent with available experimental data with a relative discrepancy of less than 3%. The liquid densities and excess molar volumes of all studied [P6,6,6,14][BOB] IL/water mixtures are characterized by concave and convex features, respectively, within full water concentration range. The dynamic properties of [P6,6,6,14] cations, [BOB] anions, and water molecules are particularly analyzed through calculation of velocity autocorrelation functions, diffusion coefficients, and reorientational autocorrelation functions and correlation times. The translational and reorientational mobilities of three species become faster upon increasing water concentration in [P6,6,6,14][BOB] IL/water mixtures and present complex dynamical characteristics arising from three distinct microscopic diffusion features within the full water concentration range. The obtained striking volumetric quantities and particular dynamic properties are well correlated to microscopic liquid structural organization and distinct local ionic environment of all studied [P6,6,6,14][BOB] IL/water mixtures.

Effect of solvent composition on the structural and dynamical properties of sodium chloride solutions in water-methanol mixtures

Molecular dynamics (MD) simulations have been performed on sodium chloride solutions in water-methanol mixtures under ambient conditions with methanol mole fractions (xm) ranging from 0.09 to 0.87. MD simulations are performed with the flexible SPC water model and the flexible six-site model for methanol. The ions are modeled as charged Lennard–Jones spheres. The structural properties of the solutions are discussed on the basis of radial distribution functions. A chain cluster of methanol is predominant in the mixture where xm exceeds 0.7. The number of methanol molecules are higher than water molecules in the first coordination shell of sodium and chloride ions beyond xm=0.5 and 0.48 respectively. With increasing xm, the diffusion coefficients for water (DH2O) and for methanol (DCH3OH) initially decrease. From xm=0.29 DCH3OH increases while DH2O rises only from about 0.7. The two diffusion coefficients are equal at xm=0.8. The deduced first and second reorientational correlation times for the OH and the dipole moment vectors for water speed up with increasing xm. Through the τ1/τ2 ratio, the reorientational motion of water molecules can be entirely ascribed to a jump. OH methanol vector reorientational occurs by a mechanism similar to that of water.

Kinetics of Rapid Covalent Bond Formation of Aniline with Humic Acid: ESR Investigations with Nitroxide Spin Labels

The bioavailability of many soil contaminants depends on their interaction with the soil organic matter. The paper presents a new approach of using stable paramagnetic spin labels for investigating the kinetics of covalent binding of specific xenobiotic functional groups with humic acids, a major organic matter fraction. Leonardite humic acid (LHA) was incubated with the nitroxide spin labels amino-TEMPO (4-amino-2,2,6,6-Tetramethylpiperidin-1-oxyl) and anilino-NO (2,5,5-Trimethyl-2-(3-aminophenyl)pyrrolidin-1-oxyl), respectively, which contain an aliphatic or aromatic functionality susceptible to interaction with LHA. Electron spin resonance spectra of LHA samples without and with the enzyme laccase were recorded at X-band frequency (9.43 GHz) at room temperature and neutral pH. Binding was detected by a pronounced broadening of the spectral lines after incubation of LHA for both spin labels. The development of a broad signal component in the spectrum of anilino-NO indicated the immobilization due to strong binding of the aniline group. The reorientational correlation time of bound anilino-NO is more than two orders of magnitude greater than that of the free label. The ratio of the amount of bound to the unbound species was used to determine the kinetics of the covalent bond formation. Reaction rate constants of 0.16 and 0.01 min−1 were determined corresponding to half-times of 4.3 and 69.3 min, respectively. Treatment of LHA with laccase enhanced the amount of the reacting anilino-NO species by a factor of 7.6, but left the reaction rate unaltered. Oxidative radical coupling was excluded using the spin trap agent n-tert-butyl-alpha-phenylnitrone.

Deuteron quadrupole coupling constants and reorientational correlation times in protic ionic liquids.

We describe a method for the accurate determination of deuteron quadrupole coupling constants χD for N-D bonds in triethylammonium-based protic ionic liquids (PILs). This approach was first introduced by Wendt and Farrar for O-D bonds in molecular liquids, and is based on the linear relationship between the deuteron quadrupole coupling constants χD, and the proton chemical shifts δ(1)H, as obtained from DFT calculated properties in differently sized clusters of the compounds. Thus the measurement of δ(1)H provides an accurate estimate for χD, which can then be used for deriving reorientational correlation-times τND, by means of NMR deuteron quadrupole relaxation time measurements. The method is applied to pure PILs including differently strong interacting anions. The obtained χD values vary between 152 and 204 kHz, depending on the cation-anion interaction strength, intensified by H-bonding. We find that considering dispersion corrections in the DFT-calculations leads to only slightly decreasing χD values. The determined reorientational correlation times indicate that the extreme narrowing condition is fulfilled for these PILs. The τc values along with the measured viscosities provide an estimate for the volume/size of the clusters present in solution. In addition, the correlation times τc, and the H-bonded aggregates were also characterized by molecular dynamics (MD) simulations.