Quantitative atom-atom potentials from rotational tunneling: their extraction and their use

Methyl rotor quantum states and the effect of chemical environment in organic crystals: γ-picoline and toluene

Using a set of first-principles calculations, we have studied the methyl tunnel splitting for molecular crystals of γ-picoline and toluene. The effective rotational potential energy surface of the probe methyl rotor along the tunneling path is evaluated using first-principles electronic structure calculations combined with the nudged elastic band method. The tunnel splitting is calculated by an explicit diagonalization of the one-dimensional time-independent Hamiltonian matrix. The effects of chemical environment and rotor-rotor coupling on the rotational energy barriers were investigated. It is found that more dense packing of the molecules in toluene compared to that in γ-picoline gives rise to a larger rotational barrier which in turn yields a considerably smaller tunnel splitting. Moreover, it turned out that coupled motion of the face-to-face methyl groups in γ-picoline has a significant effect on the reduction of the rotational barrier. Our results are in good agreement with the experimentally observed tunnel splitting.

Molecular flexibility and structural instabilities in crystalline l-methionine

We have investigated the dynamics in polycrystalline samples of l-methionine related to the structural transition at about 307K by incoherent inelastic and quasielastic neutron scattering, X-ray powder diffraction as well as ab-initio calculations. l-Methionine is a sulfur amino acid which can be considered a derivative of alanine with the alanine R-group CH3 exchanged by CH3S(CH2)2. Using X-ray powder diffraction we have observed at ~190K an anomalous drop of the c-lattice parameter and an abrupt change of the β-monoclinic angle that could be correlated to the anomalies observed in previous specific heat measurements. Distinct changes in the quasielastic region of the neutron spectra are interpreted as being due to the onset and slowing-down of reorientational motions of the CH3-S group, are clearly distinguished above 130K in crystalline l-methionine. Large-amplitude motions observed at low frequencies are also activated above 275K, while other well-defined vibrations are damped. The ensemble of our results suggests that the crystalline structure of l-methionine is dynamically highly disordered above 275K, and such disorder can be linked to the flexibility of the molecular thiol-ether group.

Methyl dynamics flattens barrier to proton transfer in crystalline tetraacetylethane

We analyze the interplay between proton transfer in the hydrogen-bond bridge, O···H···O, and lattice dynamics in the model system tetraacetylethane (TAE) (CH(3)CO)(2)CH═CH(COCH(3))(2) using density functional theory. Lattice dynamics calculations and molecular dynamics simulations are validated against neutron scattering data. Hindrance to the cooperative reorientation of neighboring methyl groups at low temperatures gives a preferred O atom for the bridging proton. The amplitude of methyl torsions becomes larger with increasing temperature, so that the free-energy minimum for the proton becomes flat over 0.2 Å. For the isolated molecule, however, we show an almost temperature-independent symmetric double-well potential persists. This difference arises from the much higher barriers to methyl torsion in the crystal that make the region of torsional phase space that is most crucial for symmetrization poorly accessible. Consequently, the proton-transfer potential remains asymmetric though flat at the base, even at room temperature in the solid.

The effect of host relaxation and dynamics on guest molecule dynamics in H2/tetrahydrofuranhydrate

We use ab initio molecular dynamics simulations to obtain classically the effects of H2O cage motions on the potential-energy surface (PES) of encapsulated H2 in the H2/tetrahydrofuran-hydrate system. The significant differences between the PES for the H2 in rigid and flexible cages that we find will influence calculation of the quantum dynamics of the H2. Part of these differences arises from the relaxation of the H2O cage around the classical H2, with a second part arising from the coupling of both translational and rotational motions of H2 with the H20 cage. We find that isotopic substitution of 2H for 1H of the H2O cage affects the coupling, which has implications for experiments that require the use of 2H2O, including inelastic neutron scattering that uses 2H2O cages in order to focus on the H2 guest dynamics. Overall, this work emphasizes the importance of taking into account cage dynamics in any approach used to understand the dynamics of H2 guests in porous framework materials.

X-ray diffraction and inelastic neutron scattering study of 1:1 tetramethylpyrazine chloranilic acid complex: temperature, isotope, and pressure effects

The x-ray diffraction studies of the title complex were carried out at room temperature and 14 K for H/D (in hydrogen bridge) isotopomers. At 82 K a phase transition takes place leading to a doubling of unit cells and alternation of the hydrogen bond lengths linking tetramethylpyrazine (TMP) and chloranilic acid molecules. A marked H/D isotope effect on these lengths was found at room temperature. The elongation is much smaller at 14 K. The infrared isotopic ratio for O-H(D)...N bands equals to 1.33. The four tunnel splittings of methyl librational ground states of the protonated complex required by the structure are determined at a temperature T=4.2 K up to pressures P=4.7 kbars by high resolution neutron spectroscopy. The tunnel mode at 20.6 microeV at ambient pressure shifts smoothly to 12.2 microeV at P=3.4 kbars. This is attributed to an increase of the strength of the rotational potential proportional to r(-5.6). The three other tunnel peaks show no or weak shifts only. The increasing interaction with diminishing intermolecular distances is assumed to be compensated by a charge transfer between the constituents of deltae/e approximately 0.02 kbar(-1). The phase transition observed between 3.4 and 4.7 kbars leads to increased symmetry with only two more intense tunneling bands. In the isotopomer with deuterated hydrogen bonds and P=1 bar all tunnel intensities become equal in consistency with the low temperature crystal structure. The effect of charge transfer is confirmed by a weakening of rotational potentials for those methyl groups whose tunnel splittings were independent of pressure. Density functional theory calculations for the model TMP.(HF)2 complex and fully ionized molecule TMP+ point out that the intramolecular rotational potential of methyl groups is weaker in the charged species. They do not allow for the unequivocal conclusions about the role of the intermolecular charge transfer effect on the torsional frequencies.

Ammonium ions in alkali metal halide crystals: tunnelling and spin relaxation

We have used low energy inelastic neutron scattering spectroscopy to examine the tunnelling spectroscopy of the ammonium ion in the (NH4)0.02Rb(x)K(0.98-x)I system. The concentration of different species were varied as x increased, this was followed systematically and the first consistent assignment scheme for these features is given. Differences were also found for the relaxation rate of the spin temperature inversions that could be generated in these species. At a critical concentration--about x = 0.04 mole fraction--the relaxation rates of the species changed dramatically.

Low frequency internal modes of 1,2,4,5-tetramethylbenzene, tetramethylpyrazine and tetramethyl-1,4-benzoquinone INS, Raman, infrared and theoretical DFT studies

The results of inelastic neutron scattering (INS), Raman and infrared (IR) studies on 1,2,4,5-tetramethylbenzene (durene), tetramethylpyrazine (TMP) and tetramethyl-1,4-benzoquinone (TMBQ) in the solid state are reported. The observed frequencies are analyzed on the basis of DFT calculations. The low frequency region, below 400 cm(-1), related to the torsional and bending out-of-plane vibrations of the CH(3) groups, is of particular interest. The detailed analysis is possible due to the simulation of the INS spectra by using the auntie-CLIMAX program. It is shown that the observed low frequency INS bands are dramatically shifted, compared to the calculated ones, towards higher frequencies. Although one cannot exclude deficiencies of theoretical methods as applied to low frequency modes, it seems more probable the interpretation based on an existence of non-conventional CH(...)pi, CH(...)N, CH(...)O hydrogen bonds formed by the methyl groups in crystalline phases.

Elastic, quasielastic, and inelastic neutron-scattering studies on the charge-transfer hexamethylbenzene-tetracyanoquinodimethane complex

The 1:1 hexamethylbenzene (HMB)-tetracyanoquinodimethane (TCNQ) complex shows a first-order phase transition at 230/218 K (heating/cooling) with no change of the space group. The neutron-diffraction studies reveal that this transition is related to a freezing of the rotation of methyl groups. The results for 100 K enabled precise determination of configuration of HMB.TCNQ complexes. The planes of HMB and TCNQ molecules from small angle (6 degrees) so that the dicyanomethylene group approaches the HMB molecule to a distance of 3.34 angstroms. The conformation of methyl groups was exactly determined. The quasielastic neutron-scattering spectra can be interpreted in terms of 120 degrees jumps with different activation barrier in low- and high-temperature phases, equal to 3.7 and 1.8 kJ/mol, respectively. These values are lower than that for neat HMB (6 kJ/mol). The conclusion can be drawn that the methyl groups can reorient more freely in the complex. This conclusion is in agreement with the results of inelastic neutron-scattering studies of low-frequency modes assigned to torsional vibrations of methyl groups. These frequencies are lower than those for neat HMB. The analyzed increase of frequencies of these modes as compared with free molecules can be interpreted as due to formation of unconventional C-H...Y hydrogen bonds which are more pronounced in crystals of neat HMB than in those of HMB.TCNQ. The low-frequency librational modes can be treated as a sensitive measure of unconventional hydrogen bonds formed by the CH3 groups.

Rotational tunneling of methyl groups in low temperature phases of mesitylene: potentials and structural implications

Mesitylene can be stabilized at He temperature in three solid phases of so far unknown crystal structures. Rotational tunneling of methyl groups is based on rotational potentials and used to characterize structural aspects. In phase III found after the first fast cooling of the sample three nonequivalent methyl rotors with splittings of 2.7, 4.1 and 16.3 microeV are observed. Three other unresolved bands are identified by their librational modes. In the second phase II the metastability is emphasized by tunneling energies still changing at temperatures T< or = 12 K. Above this temperature tunneling bands at 6.6, 12.5, 15.0 and 18.3 microeV evolve in the manner characteristic of coupling to phonons. In the equilibrium phase I a single tunnel splitting of 10.2 microeV represents all methyl groups. A unit cell containing a single molecule at a site of threefold symmetry explains quantitatively this spectrum. Phases II and III most likely contain two nonequivalent molecules in the unit cell with no local symmetry in phase II and a mirror plane in phase III. The good moderator properties for neutrons are most likely not connected to the low energy tunneling bands but to a dense vibrational phonon density of states.