A derivation of the vibrational Hamiltonian in generalized (internal) coordinates within a grid representation using the Fourier Grid Hamiltonian (FGH) method is presented. The objective of the theory is to make possible the solution of vibrational problems in two or more dimensions in terms of internal coordinates. These coordinates are often the ones of choice when the vibrations of interest are localized or when only selected coordinates are considered in a larger system. As in the case of the original FGH method, the matrix elements are easy to evaluate in a fast and robust manner. The method is tested on two different molecular systems, FHFand picolinic acid N-oxide, both containing strong hydrogen bonds. The illustrative problems are twodimensional and are highly anharmonic. The importance of both the coupling terms, as well as the variable reduced masses, required by the formulation of the Hamiltonian in generalized coordinates, are examined, and the formalism is shown to be robust in that identical results are obtained using different sets of internal coordinates applied to the same physical problem. Good agreement between calculated and observed vibrational frequencies is also obtained.

Organic materials are receiving an increasing amount of attention as electrode materials for future post lithium-ion batteries due to their versatility and sustainability. However, their electrochemical reaction mechanism has seldom been investigated. This is a direct consequence of a lack of straightforward and broadly available analytical techniques. Herein, a straightforward in operando attenuated total reflectance infrared spectroscopy method is developed that allows visualization of changes of all infrared active bands that occur as a consequence of reduction/oxidation processes. In operando infrared spectroscopy is applied to the analysis of three different organic polymer materials in lithium batteries. Moreover, this in operando method is further extended to investigation of redox reaction mechanism of poly(anthraquinonyl sulfide) in a magnesium battery, where a reduction of carbonyl bond is demonstrated as a mechanism of electrochemical activity. Conclusions done by the in operando results are complemented by synthesis of model compound and density functional theory calculation of infrared spectra.

Infrared, Raman and INS spectra of picolinic acid N-oxide (PANO) were recorded and examined for the location of the hydronic modes, particularly O-H stretching and COH bending. PANO is representative of strong chelate hydrogen bonds (H-bonds) with its short O...O distance (2.425 A). H-bonding is possibly well-characterized by diffraction, NMR and NQR data and calculated potential energy functions. The analysis of the spectra is assisted by DFT frequency calculations both in the gas phase and in the solid state. The Car-Parrinello quantum mechanical solid-state method is also used for the proton dynamics simulation; it shows the hydron to be located about 99% of time in the energy minimum near the carboxylic oxygen; jumps to the N-O acceptor are rare. The infrared spectrum excels by an extended absorption (Zundel's continuum) interrupted by numerous Evans transmissions. The model proton potential functions on which the theories of continuum formation are based do not correspond to the experimental and computed characteristics of the hydrogen bond in PANO, therefore a novel approach has been developed; it is based on crystal dynamics driven hydronium potential fluctuation. The envelope of one hundred 0 --> 1 OH stretching transitions generated by molecular dynamics simulation exhibits a maximum at 1400 cm-1 and a minor hump at approximately 1600 cm-1. These positions square well with ones predicted for the COH bending and OH stretching frequencies derived from various one- and two-dimensional model potentials. The coincidences with experimental features have to be considered with caution because the CPMD transition envelope is based solely on the OH stretching coordinate while the observed infrared bands correspond to heavily mixed modes as was previously shown by the normal coordinate analysis of the IR spectrum of argon matrix isolated PANO, the present CPMD frequency calculation and the empirical analysis of spectra. The experimental infrared spectra show some unusual characteristics such as large temperature effects on the intensity of some bands, thus presenting a challenge for theoretical band shape treatments. Our calculations clearly show that the present system is characterized by an asymmetric single well potential with no large amplitudes in the hydronium motion, which extends the existence of Zundel-type spectra beyond the established set of hydrogen bonds with large hydronic vibrational amplitudes.

No abstract

2D free energy surfaces V = V(rOH, rO...O) for the intramolecular H-bond in the title compound were calculated by the DFT method and used in the calculation of primary and secondary chemical shifts of the compound dissolved in chloroform and acetonitrile. Solvent effects were accounted for by the SCRF/PCM method. The corresponding two-dimensional chemical shift surfaces with included solvent reaction field were obtained using the Continuous Set of Gauge Transformations approach at the B3LYP/6-311+G(2d,2p) level of theory. The chemical shifts were estimated as quantum averages along the two internal coordinates in the hydrogen bond and along several vibrational levels according to the Boltzmann distribution at room temperature. Fairly good agreement between the experimental and calculated isotope effects was obtained. 1D and 2D NMR spectra of solutions of picolinic acid N-oxide and its deuterated analogue were recorded and assigned.

The I335Y point mutation effect on the kinetics of phenylethylamine decomposition catalyzed by monoamine oxidase A was elucidated by means of molecular simulation. The established empirical valence bond methodology was used in conjunction with the free energy perturbation sampling technique and a classical force field representing the state of reactants and products. The methodology allows for the simulation of chemical reactions, in the present case the breaking of the α-C-H bond in a phenylethylamine substrate and the subsequent hydrogen transfer to the flavin cofactor, resulting in the formation of the N-H bond on flavin. The empirical parameters were calibrated against the experimental data for the simulated reaction in a wild type protein and then used for the calculation of the reaction free energy profile in the I335Y mutant. In very good agreement with the measured kinetic data, mutation increases the free energy barrier for the rate limiting step by slightly more than 1 kcal mol(-1) and consequently decreases the rate constant by about an order of magnitude. The magnitude of the computed effect slightly varies with simulation settings, but always remains in reasonable agreement with the experiment. Analysis of trajectories reveals a major change in the interaction between phenyl rings of the substrate and the neighboring Phe352 residue upon the I335Y mutation due to the increased local polarity, leading to an attenuated quadrupole interaction between the rings and destabilization of the transition state. Additionally, the increased local polarity in the mutant allows for a larger number of water molecules to be present near the active site, effectively shielding the catalytic effect of the enzyme and contributing to the increased barrier.

Polarized Raman spectra of single crystals of the α-polymorphs of protonated and deuterated oxalic acid dihydrate were recorded. The interpretation of the spectra is assisted by periodic DFT calculations using the CRYSTAL06 program and by comparison with the infrared spectra of the polycrystalline material. The agreement between the calculated and observed band wavenumbers is fair in the case of low-anharmonicity modes, but marked differences appear for the stretching modes that are strongly anharmonic. A very broad feature, extending between ∼2000 and 1200 cm −1 , is attributed to OH stretching. Notable is the topping of this feature by distinct bands that can be attributed to C O stretching, H 2 O scissoring and COH bending coupled to C-O stretching. The assignments are supported by isotope effects. However, deuteration does not notably affect the wavenumber limits of the broad OH stretching band, which suggests that the potential governing the proton dynamics is of the asymmetric double-minimum type with a very low barrier. The calculated normal coordinates show a strong participation of the bending modes of water molecules in almost all internal acid motions, as well as in the external phonons.

No abstract