Anomalous temperature dependence of nuclear quadrupole interactions in strongly hydrogen-bonded systems from first principles

Looking into the dynamics of molecular crystals of ibuprofen and terephthalic acid using 17 O and 2 H nuclear magnetic resonance analyses

Oxygen-17 and deuterium are two quadrupolar nuclei that are of interest for studying the structure and dynamics of materials by solid-state nuclear magnetic resonance (NMR). Here, 17 O and 2 H NMR analyses of crystalline ibuprofen and terephthalic acid are reported. First, improved 17 O-labelling protocols of these molecules are described using mechanochemistry. Then, dynamics occurring around the carboxylic groups of ibuprofen are studied considering variable temperature 17 O and 2 H NMR data, as well as computational modelling (including molecular dynamics simulations). More specifically, motions related to the concerted double proton jump and the 180° flip of the H-bonded (-COOH)2 unit in the crystal structure were looked into, and it was found that the merging of the C=O and C-OH 17 O resonances at high temperatures cannot be explained by the sole presence of one of these motions. Lastly, preliminary experiments were performed with a 2 H-17 O diplexer connected to the probe. Such configurations can allow, among others, 2 H and 17 O NMR spectra to be recorded at different temperatures without needing to tune or to change probe configurations. Overall, this work offers a few leads which could be of use in future studies of other materials using 17 O and 2 H NMR.

Importance of Nuclear Quantum Effects for NMR Crystallography

The resolving power of solid-state nuclear magnetic resonance (NMR) crystallography depends heavily on the accuracy of computational predictions of NMR chemical shieldings of candidate structures, which are usually taken to be local minima in the potential energy. To test the limits of this approximation, we systematically study the importance of finite-temperature and quantum nuclear fluctuations for ¹H, ¹³C, and ¹⁵N shieldings in polymorphs of three paradigmatic molecular crystals: benzene, glycine, and succinic acid. The effect of quantum fluctuations is comparable to the typical errors of shielding predictions for static nuclei with respect to experiments, and their inclusion improves the agreement with measurements, translating to more reliable assignment of the NMR spectra to the correct candidate structure. The use of integrated machine-learning models, trained on first-principles energies and shieldings, renders rigorous sampling of nuclear fluctuations affordable, setting a new standard for the calculations underlying NMR structure determinations.

How does the NMR thermometer work? Application of combined quantum molecular dynamics and GIPAW calculations into the study of lead nitrate

Lead nitrate is an inorganic salt, commonly used for the accurate temperature determination in the solid state NMR spectroscopy, due to the strong temperature dependence of the ²⁰⁷ Pb chemical shift. As the reason for this phenomenon remained unknown, the main purpose of this study was to explain this temperature dependence at the molecular level. To achieve this, combined CASTEP geometry optimization, quantum molecular dynamics at chosen temperatures and GIPAW NMR computations were performed. Due to the previous literature reports on inaccuracy in the calculation of ²⁰⁷ Pb NMR parameters using GIPAW, a large emphasis was put on the optimization of computational method. The application of quantum molecular dynamics provided the simulation of the temperature-dependent vibrational motions and enabled to accurately compute the changes in the value of Pb δ_iso resulting from them. © 2018 Wiley Periodicals, Inc.

Electron-phonon coupling from finite differences

The interaction between electrons and phonons underlies multiple phenomena in physics, chemistry, and materials science. Examples include superconductivity, electronic transport, and the temperature dependence of optical spectra. A first-principles description of electron-phonon coupling enables the study of the above phenomena with accuracy and material specificity, which can be used to understand experiments and to predict novel effects and functionality. In this topical review, we describe the first-principles calculation of electron-phonon coupling from finite differences. The finite differences approach provides several advantages compared to alternative methods, in particular (i) any underlying electronic structure method can be used, and (ii) terms beyond the lowest order in the electron-phonon interaction can be readily incorporated. But these advantages are associated with a large computational cost that has until recently prevented the widespread adoption of this method. We describe some recent advances, including nondiagonal supercells and thermal lines, that resolve these difficulties, and make the calculation of electron-phonon coupling from finite differences a powerful tool. We review multiple applications of the calculation of electron-phonon coupling from finite differences, including the temperature dependence of optical spectra, superconductivity, charge transport, and the role of defects in semiconductors. These examples illustrate the advantages of finite differences, with cases where semilocal density functional theory is not appropriate for the calculation of electron-phonon coupling and many-body methods such as the GW approximation are required, as well as examples in which higher-order terms in the electron-phonon interaction are essential for an accurate description of the relevant phenomena. We expect that the finite difference approach will play a central role in future studies of the electron-phonon interaction.

Temperature dependence of X-ray absorption and nuclear magnetic resonance spectra: probing quantum vibrations of light elements in oxides

Probing the quantum thermal fluctuations of nuclei in light-element oxides using XANES and NMR spectroscopies.

Combining solid-state NMR spectroscopy with first-principles calculations - a guide to NMR crystallography

Recent advances in the application of first-principles calculations of NMR parameters to periodic systems have resulted in widespread interest in their use to support experimental measurement. Such calculations often play an important role in the emerging field of "NMR crystallography", where NMR spectroscopy is combined with techniques such as diffraction, to aid structure determination. Here, we discuss the current state-of-the-art for combining experiment and calculation in NMR spectroscopy, considering the basic theory behind the computational approaches and their practical application. We consider the issues associated with geometry optimisation and how the effects of temperature may be included in the calculation. The automated prediction of structural candidates and the treatment of disordered and dynamic solids are discussed. Finally, we consider the areas where further development is needed in this field and its potential future impact.

Temperature Dependence of NMR Parameters Calculated from Path Integral Molecular Dynamics Simulations

The influence of temperature on NMR chemical shifts and quadrupolar couplings in model molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory (DFT) calculations of shielding and electric field gradient (EFG) tensors. An approach based on convoluting calculated shielding or EFG tensor components with probability distributions of selected bond distances and valence angles obtained from DFT-PIMD simulations at several temperatures is used to calculate the temperature effects. The probability distributions obtained from the quantum PIMD simulations, which includes nuclear quantum effects, are significantly broader and less temperature dependent than those obtained with conventional DFT molecular dynamics or with 1D scans through the potential energy surface. Predicted NMR observables for the model systems were in excellent agreement with experimental data.

Kinetics of Coupled Double Proton and Deuteron Transfer in Hydrogen-Bonded Ribbons of Crystalline Pyrazole-4-carboxylic Acid

The proton tautomerism of pyrazole-4-carboxylic acid (PCA) has been studied by a combination of 15N CPMAS and 2H NMR spectroscopy and relaxometry. Down to 250 K, PCA forms a hydrogen bonded ribbon where adjacent carboxylic and pyrazole groups are linked by an OH···N and an O···HN hydrogen bond, forming either the tautomeric state A or B. Down to about 250 K, the tautomerism is fast on the NMR timescale and degenerate, corresponding to a phase exhibiting dynamic proton disorder. At lower temperatures, a transition to an ordered phase is observed with localized protons, assigned to an all-syn conformation adopting the sequence of tautomeric states ..ABABA.. The longitudinal 15N relaxation times T 1 of PCA-15N2 have been measured at 9.12 MHz (2.1 T). Because of the low field, a chemical shift anisotropy mechanism could be neglected, and the data were analyzed in terms of a dipolar 1H-15N relaxation mechanism, yielding the rate constants k HH. The rate constants k HD and k DD were obtained from the measurement and analysis of the 2H T 1 values of PCA-15N1-d0.9 and PCA-15N1-d0.1 measured at 46.03 MHz. Within the margin of error, no kinetic isotope effects could be detected, in contrast to previous results reported for the very fast tautomerism of solid benzoic acid dimers and the much slower tautomerism of solid 3,5-diphenyl-4-brompyrazole (DPBrP) dimers. The Arrhenius curves of all three systems were simulated using the Bell–Limbach tunneling model. Evidence for a major heavy atom motion for the tautomerism of PCA is obtained, associated with small angle reorientation of PCA molecules around the molecular axis. The observed proton order-disorder transition and the mechanism of the observed rate process are discussed in terms of a coupling of adjacent tautomeric states.

Vibrational averaging of the chemical shift in crystalline α-glycine

Averaging of the chemical shift over the molecular motion improves the simulated data and provides additional information about the temperature dependence and system dynamics. However, crystal modeling is difficult due to the limited precision of the plane-wave density functional theory (DFT) methods and approximate vibrational schemes. On the glycine example, we investigate how the averaging can be achieved within the periodic boundary conditions at the DFT level. The nuclear motion is modeled with the vibrational configuration interaction, with other simplified quantum anharmonic schemes, and the classical Born-Oppenheimer molecular dynamics (BOMD). The results confirm a large vibrational contribution to the isotropic shielding values. Both the first and second derivatives of the shielding were found important for the quantum averaging. The first derivatives influence the shielding mostly due to the anharmonic character of the CH and NH stretching modes, whereas second derivatives produce most vibrational corrections associated with the lower-frequency vibrational modes. Temperature excitations of the lowest-frequency vibrational states and the expansion of the crystal cell both determine the temperature dependence of nuclear magnetic resonance parameters. The vibrational quantum approach as well as classical BOMD schemes provided temperature dependencies of the chemical shifts that are consistent with the previous experimental data.

High-Resolution Solid-State NMR Studies of Poly(vinyl phosphonic acid) Proton-Conducting Polymer: Molecular Structure and Proton Dynamics

The structure and the local proton mobility of poly(vinyl phosphonic acid) were studied by solid-state NMR under fast magic-angle spinning. At elevated temperatures, the signature of the hydrogen-bonded P-OH protons is observed in 1H magic-angle spinning (MAS) NMR as a single resonance at 10.5 ppm. Both 1H double-quantum NMR and variable-temperature experiments demonstrate that P-OH protons are mobile and thus able to contribute to proton conductivity. Below room temperature, two different types of hydrogen-bonded P-OH resonances are observed at 10.5 and 15 ppm, and 1H double-quantum NMR demonstrates that these protons are immobile on the NMR time scale. By means of first-principles calculations of a model polymer, we have assigned the additional hydrogen-bonded species at lower temperatures to phosphonic acid anhydride and charged anhydride. Also, in the 31P MAS NMR spectrum, two distinct resonances appear, arising from "normal" phosphonic acid and phosphonic acid anhydride. 31P double-quantum NMR experiments reveal that there is no phase segregation between normal and phosphonic acid anhydride and the condensation reaction occurs randomly throughout the system. The formation of acid anhydride leads to a decrease in proton conductivity through two mechanisms, (1) decrease in the number of charge carriers and (2) blockage of charge transport pathways through immobilization of charge carriers together with a hindered reorientation of the anhydride group. Our results provide strong evidence for these mechanisms by demonstrating that the conductivity is greatly influenced by the presence of phosphonic acid anhydride.

First-Principles Molecular Dynamics Evaluation of Thermal Effects on the NMR1JLi,C Spin–Spin Coupling

Car-Parrinello (CP) molecular dynamics were applied to sample conformations of various models of organolithium aggregates which are chosen to estimate (1)J(Li,C) NMR coupling constants. The results show that the deviations from the values computed using static (optimized) geometries are small provided no large-amplitude motions occur within the timescale of the simulations. In the case of the vinyllithium dimer, for which rotation of the vinyl chain is observed, this approach allows analysis of the various contributions to the experimentally measured constants. For the trisolvated methyllithium monomer, partial decoordination of solvating dimethyl ether is observed and results in a significant shift of (1)J(Li,C). All these results highlight that a varied physicochemical machinery is hidden behind general empirical formulas, such as the Bauer-Winchester-Schleyer rule used experimentally.

First Principles and Experimental 1H NMR Signatures of Solvated Ions: The Case of HCl(aq)

A combined experimental and ab initio study is presented of the 1H NMR chemical shift distribution of aqueous hydrogen chloride solution as a function of acid concentration, based on Car-Parrinello molecular dynamics simulations and fully periodic NMR chemical-shift calculations. The agreement of computed and experimental spectra is very good. From first-principles calculations, we can show that the individual contributions of Eigen and Zundel ions, regular water molecules, and the chlorine solvation shell to the NMR line are very distinct and almost independent of the acid concentration. From the computed instantaneous NMR distributions, it is further possible to characterize the average variation in hydrogen-bond strength of the different complexes.

Bulk Chemical Shifts in Hydrogen-Bonded Systems from First-Principles Calculations and Solid-State-NMR

We present an analysis of bulk (1)H NMR chemical shifts for a series of biochemically relevant molecular crystals in analogy to the well-known solvent NMR chemical shifts. The term bulk shifts denotes the change in NMR frequency of a gas-phase molecule when it undergoes crystallization. We compute NMR parameters from first-principles electronic structure calculations under full periodic boundary conditions and for isolated molecules and compare them to the corresponding experimental fast magic-angle spinning solid-state NMR spectra. The agreement between computed and experimental lines is generally very good. The main phenomena responsible for bulk shifts are packing effects (hydrogen bonding and pi-stacking) in the condensed phase. By using these NMR bulk shifts in well-ordered crystalline model systems composed of biologically relevant molecules, we can understand the individual spectroscopic signatures of packing effects. These local structural driving forces, hydrogen bonding, pi-stacking, and related phenomena, stand as a model for the forces that govern the assembly of much more complex supramolecular aggregates. We show to which accuracy condensed-phase ab initio calculations can predict structure and structure-property relationships for noncovalent interactions in complex supramolecular systems.