Experimental and theoretical studies of malonaldehyde: vibrational analysis of a strongly intramolecularly hydrogen bonded compound

MC-QTAIM analysis reveals an exotic bond in coherently quantum superposed malonaldehyde

The proton between the two oxygen atoms of the malonaldehyde molecule experiences an effective double-well potential in which the proton's wavefunction is delocalized between the two wells. Herein we employ a state-of-the-art multi-component quantum theory of atoms in molecules partitioning scheme to obtain the molecular structure, i.e. atoms in molecules and bonding network, from the superposed ab initio wavefunctions of malonaldehyde. In contrast to the familiar clamped-proton portrayal of malonaldehyde, in which the proton forms a hydrogen basin, for the superposed states the hydrogen basin disappears and two novel hybrid oxygen-hydrogen basins appear instead, with an even distribution of the proton population between the two basins. The interaction between the hybrid basins is stabilizing thanks to an unprecedented mechanism. This involves the stabilizing classical Coulomb interaction of the one-proton density in one of the basins with one-electron density in the other basin. This stabilizing mechanism yields a bond foreign to the known bonding modes in chemistry.

Finding the direct energy-structure correlations in intramolecular aromaticity assisted hydrogen bonding (AAHB)

A predictive model for intramolecular hydrogen bond energy (E_HB) calculation of polyaromatic ortho-hydroxyaldehydes based on a set of small, functionalized hydrocarbons is developed. The complete data set of 18 compounds was used for this study. The model is based on one of four optional categories of molecular descriptors: geometric, spectroscopic, bond order and topological indices. The model of Wiberg bond indices (WBIs) as descriptors of the CC involved bond based on stepwise regression has acceptable prediction abilities for 14 structures of ortho-hydroxyformylobenzo[a]pyrene derivatives already at the semiempirical level. The presented correlation enables a significantly more rapid and quantitative description of the hydrogen bonding strength than the much more time-consuming MTA method. Thus, WBIs are shown to provide a reliable means for fast prescreening of the energy of chelate hydrogen bonds potentially for any polyaromatic derivatives.

Controlled intramolecular H-transfer in malonaldehyde in the electronic ground state mediated through the conical intersection of 1nπ* and 1ππ* excited electronic states

We report photo-isomerization of malonaldehyde in its electronic ground state (S0), mediated by coupled 1nπ*(S1)-1ππ*(S2) excited electronic states, accomplished with the aid of optimally designed ultraviolet (UV)-laser pulses. In particular, control of H-transfer from a configuration predominantly located in the left well (say, reactant) to that in the right well (say, product) of the electronic ground S0 potential energy surface is achieved by a pump-dump mechanism including the nonadiabatic interactions between the excited S1 and S2 states. An interplay between the nonadiabatic coupling due to the conical intersection of the S1 and S2 states and the laser-molecule interaction is found to be imprinted in the time-dependent electronic population. The latter is also examined by employing optimal fields of varying intensities and frequencies of the UV laser pulses. For the purpose of the present study, we constructed a three-state and two-mode coupled diabatic Hamiltonian with the help of adiabatic electronic energies and transition dipole moments calculated by ab initio quantum chemistry methods. The electronic diabatic model is developed using the calculated adiabatic energies of the two excited electronic states (S1 and S2) in order to carry out the dynamics study. The optimal fields for achieving the controlled isomerization are designed within the framework of optimal control theory employing the optimization technique of a multitarget functional using the genetic algorithm. The laser-driven dynamics of the system is treated by numerically solving the time-dependent Schrödinger equation within the dipole approximation. A time-averaged yield of the target product of ∼40% is achieved in the present treatment of dynamics with optimal laser pulses.

Kinetic isotope effects and how to describe them

We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.

Intramolecular hydrogen bonding in malonaldehyde and its radical analogues

High level Brueckner doubles with triples correction method-based ab initio calculations have been used to investigate the nature of intramolecular hydrogen bonding and intramolecular hydrogen atom transfer in cis-malonaldehyde (MA) and its radical analogues. The radicals considered here are the ones that correspond to the homolytic cleavage of C-H bonds in cis-MA. The results suggest that cis-MA and its radical analogues, cis-MA_RS, and cis-MA_RA, both exist in planar geometry. The calculated intramolecular O-H⋯O=C bond in cis-MA is shorter than that in the radical analogues. The intramolecular hydrogen bond in cis-MA is stronger than in its radicals by at least 3.0 kcal/mol. The stability of a cis-malonaldehyde radical correlates with the extent of electron spin delocalization; cis-MA_RA, in which the radical spin is more delocalized, is the most stable MA radical, whereas cis-MA_RS, in which the radical spin is strongly localized, is the least stable radical. The natural bond orbital analysis indicates that the intramolecular hydrogen bonding (O⋯H⋯O) in cis-malonaldehyde radicals is stabilized by the interaction between the lone pair orbitals of donor oxygen and the σ^* orbital of acceptor O-H bond (n → σ^*_OH). The calculated barriers indicate that the intramolecular proton transfer in cis-MA involves 2.2 kcal/mol lower barrier than that in cis-MA_RS.

NMR and IR Investigations of Strong Intramolecular Hydrogen Bonds

For the purpose of this review, strong hydrogen bonds have been defined on the basis of experimental data, such as OH stretching wavenumbers, νOH, and OH chemical shifts, δOH (in the latter case, after correction for ring current effects). Limits for O-H···Y systems are taken as 2800 > νOH > 1800 cm-1, and 19 ppm > δOH > 15 ppm. Recent results as well as an account of theoretical advances are presented for a series of important classes of compounds such as β-diketone enols, β-thioxoketone enols, Mannich bases, proton sponges, quinoline N-oxides and diacid anions. The O···O distance has long been used as a parameter for hydrogen bond strength in O-H···O systems. On a broad scale, a correlation between OH stretching wavenumbers and O···O distances is observed, as demonstrated experimentally as well as theoretically, but for substituted β-diketone enols this correlation is relatively weak.

Theoretical study of the C-H/O-H stretching vibrations in malonaldehyde

IR and Raman spectra of the malonaldehyde molecule and its deuterated analogues were calculated in the B3LYP/cc-pVQZ approximation. Anharmonicity effects were taken into account both in the context of a standard model of the second order perturbation theory and by constructing the potential energy surfaces (PES) with a limited number of dimensions using the Cartesian coordinates of the hydroxyl hydrogen atom and the stretching coordinates of С-Н, C-D, O-H, and O-D bonds. It was shown that in each of the two equivalent forms of the molecule, besides the global minimum, an additional local minimum at the PES is formed with the energy more than 3,000 cm(-1) higher than the energy in the global minimum. Calculations carried out by constructing the 2D and 3D PESs indicate a high anharmonicity level and multiple manifestations of the stretching О-Н vibrations, despite the fact that the model used does not take into account the splitting of the ground-state and excited vibrational energy levels. In particular, the vibration with the frequency 3,258 cm(-1) may be associated with proton transfer to the region of a local minimum of energy. Comparing the results obtained with the experimental data presented in the literature allowed us to propose a new variant of bands assignments in IR and Raman spectra of the molecule in the spectral region 2,500-3,500 cm(-1).

Estimation of the intramolecular O-H···O═C hydrogen bond energy via the molecular tailoring approach. Part I: aliphatic structures

A simple and universal method for the estimation of the intramolecular hydrogen bond (HB) energy (E(HB)) in hydroxycarbonyl aliphatic compounds is proposed by the application of the molecular tailoring approach (MTA) based on calculations at the second-order Møller-Plesset MP2 level. The calculation of EHB can be realized by the one optimization and three single point calculations of the energy for each compound with carbonyl and hydroxyl groups involved in HB. The intramolecular hydrogen bond energies estimated for 153 structures (of 102 compounds) ranged from 1.4 to 13.7 kcal/mol for systems without resonance-assisted hydrogen bonding (RAHB). To verify the method, we show the correlations of the energy (E(HB)) in six-, seven-, and eight-membered HB rings in the optimized multifunctional molecules with the usual geometry descriptors of hydrogen bonds. Moreover, topological parameters from the atoms in molecules (AIM) theory and the calculated infrared and proton NMR spectra are correlated. The effects of conjugation and π-electron delocalization, bifurcation, and cooperativity are discussed, along with the correlation between the strength and geometrical parameters of H bonding.

Vibrational conical intersections as a mechanism of ultrafast vibrational relaxation

Presenting true crossings of adiabatic potential energy surfaces, conical intersections are a paradigm of ultrafast and efficient electronic relaxation dynamics. The same mechanism is shown to apply also for vibrational conical intersections, which may occur when two high-frequency modes (such as OH stretch vibrations) are coupled to low-frequency modes (such as hydrogen bonding modes). By derivation of a model Hamiltonian and its parametrization for a concrete example, malonaldehyde, the conditions that such conical intersections occur are identified and the consequences for the vibrational dynamics and spectra are demonstrated.

A generalized reactive force field for nonlinear hydrogen bonds: hydrogen dynamics and transfer in malonaldehyde

Using molecular dynamics (MD) simulations, the spectroscopy and dynamics of malonaldehyde is investigated. To this end, the recently proposed molecular mechanics with proton transfer (MMPT) potential is generalized to nonlinear hydrogen bonds. The calculated properties for malonaldehyde in both gas and condensed phases, including equilibrium geometries, infrared spectra, tunneling splittings, and hydrogen transfer rates, compare well with previous experimental and computational works. In particular, by using a harmonic bath averaged (HBA) Hamiltonian, which is based on a reaction path Hamiltonian, it is possible to estimate the tunneling splitting in an efficient manner. It is found that a zero point corrected barrier of 6.7 kcal/mol and effective masses of 1.234 (i.e., 23.4% larger than the mass of a physical H-atom) and 1.117 (for the physical D-atom) are consistent with the measured splittings of 21.6 and 2.9 cm(-1), respectively. The HBA Hamiltonian also yields a pair of hydrogen transfer fundamentals at 1573 and 1267 cm(-1), similar to results obtained with a reaction surface Hamiltonian on a MP2/6-31G(**) potential energy surface. This amounts to a substantial redshift of more than 1000 cm(-1) which can be rationalized by comparison with weakly (HCO(+): rare gas) and strongly (H(2)O-H(+)-OH(2)) proton-bound systems. Hydrogen transfer rates in vacuum and water were determined from the validated MMPT potential and it is found that the solvent enhances the rate by a factor of 5 at 300 K. The rates of 2.4/ns and 10/ns are commensurate with previous density functional tight binding ab initio MD studies.

The Intramolecular Hydrogen-Bond in Malonaldehyde as Seen by Infrared Spectroscopy. A Four-Dimensional Model Study

The infrared spectrum of the O–H–O fragment of malonaldehyde is studied using a four-dimensional model. This comprises the OH stretching and the two OH bending vibrations as well an O–O ring deformation mode under the assumption of overall Cs symmetry. The full anharmonic potential energy and dipole moment surfaces are calculated using density functional theory and the respective vibrational eigenvalue problem is solved by an iterative Lanczos method. Fundamental and combination transitions are discussed for the normal species and the symmetrically deuterated isotopomer. Special emphasis is paid to the OH/OD stretching region which reveals the signatures of strong mode mixing what renders a simple assignment in terms of fundamental transitions difficult. In addition results for hot transitions are presented which show a rather different OH/OD band due to the topology of the potential energy surface. The influence of H atom tunneling on the spectrum is briefly addressed employing an alternative three-dimensional model which takes into account the in-plane H atom motion as well as the O–O distance.

Experimental matrix isolation study and quantum-mechanics-based normal-coordinate analysis of the anharmonic infrared spectrum of picolinic acid N-oxide

This work is, according to our knowledge, the first experimental matrix isolation study of a molecular system with a very short and strong intramolecular OH...O hydrogen bond. It also includes a satisfying interpretation of its entire infrared spectrum. The interpretation relies on the calculation at the DFT/B3LYP/6-31G(d,p) level of the harmonic spectrum and of the anharmonic relaxed potential energy for the stretching motion of the hydrogen-bonded proton, used with our recently modified quantum-mechanics-based normal-coordinate analysis. An important observation about the anharmonic spectrum obtained from this procedure is that the OH stretch coordinate contributes to several normal modes, mixing extensively with other in-plane internal coordinates, in particular OH-bending and C=O-stretching. The two intense normal modes with the largest contributions from the OH-stretching coordinate to the potential energy distribution and to the intensity are located near 1700 and 1500 cm(-1). A calculated anharmonic spectrum obtained from this procedure agrees with the experimental spectrum (frequencies and intensity distribution), within the limits of the estimated uncertainties for the calculation and experiment, allowing the interpretation of the latter. The agreement for the frequencies is about 1-3%. The anharmonic spectrum calculated using the anharmonic keyword in Gaussian 03w is not in satisfactory agreement with experiment insofar as the OH-stretching mode is concerned.

A reaction surface Hamiltonian study of malonaldehyde

We report calculations using a reaction surface Hamiltonian for which the vibrations of a molecule are represented by 3N-8 normal coordinates, Q, and two large amplitude motions, s(1) and s(2). The exact form of the kinetic energy operator is derived in these coordinates. The potential surface is first represented as a quadratic in Q, the coefficients of which depend upon the values of s(1),s(2) and then extended to include up to Q(6) diagonal anharmonic terms. The vibrational energy levels are evaluated by solving the variational secular equations, using a basis of products of Hermite polynomials and appropriate functions of s(1),s(2). Our selected example is malonaldehyde (N=9) and we choose as surface parameters two OH distances of the migrating H in the internal hydrogen transfer. The reaction surface Hamiltonian is ideally suited to the study of the kind of tunneling dynamics present in malonaldehyde. Our results are in good agreement with previous calculations of the zero point tunneling splitting and in general agreement with observed data. Interpretation of our two-dimensional reaction surface states suggests that the OH stretching fundamental is incorrectly assigned in the infrared spectrum. This mode appears at a much lower frequency in our calculations due to substantial transition state character.

Tautomeric equilibria and pi electron delocalization for some monohydroxyarenes--quantum chemical studies

Keto-enol tautomeric interconversions and variations of the pi-electron distribution were studied for 11 isolated monohydroxyarenes at the DFT(B3LYP)/6-311++G(2df,2p) level. For two monohydroxyarenes (phenol and 9-anthrol), the PCM model of solvation (water) was also applied to the DFT geometries. The geometry-based HOMA index was applied to estimate pi-electron delocalization in the keto and enol tautomeric forms. Thermodynamic parameters of tautomeric interconversions (DeltaET, DeltaGT, TDeltaST, pKT) were calculated to estimate relative stabilities of individual tautomers and their percentage contents in the tautomeric mixtures. In almost all cases, the aromatic enol forms are strongly favored. An exception is 9-anthrol, which prefers its keto form. The resonance stabilization of this form comes from the central ring. Generally, aromaticity is the main factor that influences tautomeric equilibria in monohydroxyarenes. Hydration effect is considerably smaller and it does not change the tautomeric preference.

Investigation of the release mechanism of a sparingly water-soluble drug from solid dispersions in hydrophilic carriers based on physical state of drug, particle size distribution and drug–polymer interactions

In the present study the release mechanism of the sparingly water-soluble drug felodipine (FELO) from particulate solid dispersions in PVP or PEG was investigated. FT-IR data indicated that a N-H...O hydrogen bond is formed between FELO and polymers. The drug-polymer interaction was theoretically studied with the density functional theory with the B3LYP exchange correlation function. The interaction energies have been estimated at -31.8 kJ/mol for PVP and -18.8 kJ/mol for PEG. Also, detailed vibrational analysis of the complexes showed that the red shift of the N-H bond stretching in FELO molecule due to H-bonding was higher in the FELO-PVP complex than in the FELO-PEG complex. Both the experimental and theoretical data indicated that a stronger interaction of FELO with PVP than with PEG was developed. The interactions of FELO with the polymer appeared to control the physical state (amorphous or crystalline) and the particle size of FELO in the solid dispersions. In the FELO/PVP dispersions, the drug is found as amorphous nanoparticles whereas in FELO/PEG dispersions the drug is dispersed as crystalline microparticles. The size of drug particles in the dispersion was also influenced by drug proportion, with an increase in drug content of the dispersion resulting in increased drug particle size. The particle size of drug, the proportion of drug in the dispersion and the properties of the polymer (molecular weight) appeared to determine the mechanism of drug release from the solid dispersions, which was drug diffusion (through the polymer layer)-controlled at low drug contents and drug dissolution-controlled at high drug contents. In situ DLS measurements indicate that the large initial particles of FELO/PVP and FELO/PEG solid dispersions with low drug content (10-20 wt%) are very rapidly decreased to smaller particles (including nanoparticles) during dissolution, leading to the observed impressive enhancement of FELO release rate from these dispersions.

UV and IR Photoisomerization of Acetylacetone Trapped in a Nitrogen Matrix

UV- and IR-induced photoisomerization of acetylacetone trapped in a nitrogen matrix at 4.3 K have been carried out using a tunable optical parametric oscillator type laser, or a mercury vapor lamp, coupled with Fourier Transform IR and UV spectroscopies. After deposition, the main form present in the cryogenic matrix is that chelated (enol). Upon UV irradiation, the intramolecular H bond is broken leading to nonchelated isomers among seven possible open forms. These forms have then been irradiated by resonant pi* <-- pi UV irradiation, or by resonant nuOH irradiation. The selective UV irradiation allows us to suggest a first vibrational assignment while the nuOH irradiation leads us to observe interconversions between the nonchelated isomers. In order to support our vibrational assignment, we have carried out theoretical calculations at the B3LYP/6-311++G(2d,2p) level of theory. This study shows that only five isomers are observed among eight postulated.

Vibrational assignment and structure of trifluorobenzoylacetone. A density functional theoretical study

Molecular structure and vibrational frequencies of 4,4,4-trifluoro-1-phenyl-1,3-butanedione, known as trifluorobenzoylacetone (TFBA), have been investigated by means of density functional theory (DFT) calculations. The results were compared with those of benzoylacetone (BA), acetylacetone (AA), and trifluoroacetylacetone (TFAA). Comparing the calculated and experimental band frequencies and intensities suggests coexisting of both stable cis-enol conformers in comparable proportions in the sample. The energy difference between the two stable chelated enol forms is negligible, 0.96 kcal/mol, calculated at B3LYP/6-311++G** level of theory. The molecular stability and the hydrogen bond strength were investigated by applying the natural bond orbital (NBO) theory and geometry calculations. The theoretical calculations and spectroscopic results indicate that the hydrogen bond strength of TFBA is between those of TFAA and AA, considerably weaker than that of BA.

Vibrational assignment and structure of dibenzoylmethane. A density functional theoretical study

Molecular structure and vibrational frequencies of 1,3-diphenyl-1,3-propanedione, known as dibenzoylmethane (DBM), have been investigated by means of density functional theory (DFT) calculations. The results were compared with those of benzoylacetone (BA) and acetylacetone (AA), the parent molecule. IR and Raman spectra of DBM and its deuterated analogue were clearly assigned. The calculated hydrogen bond energy of DBM is 16.15 kcal/mol, calculated at B3LYP/6-311++G** level of theory, which is 0.28 kcal/mol more than that of AA. This result is in agreement with the vibrational and NMR spectroscopy results. The molecular stability and the hydrogen bond strength were investigated by applying the Natural Bond Orbital analysis (NBO) and geometry calculations. The theoretical calculations indicate that the hydrogen bond in DBM is relatively stronger than that in BA and AA.

Structure and vibrational assignment of the enol form of 1,1,1-trifluoro-2,4-pentanedione

Molecular structure of 1,1,1-trifluoro-pentane-2,4-dione, known as trifluoro-acetylacetone (TFAA), has been investigated by means of Density Functional Theory (DFT) calculations and the results were compared with those of acetylacetone (AA) and hexafluoro-acetylacetone (HFAA). The harmonic vibrational frequencies of both stable cis-enol forms were calculated at B3LYP level of theory using 6-31G** and 6-311++G** basis sets. We also calculated the anharmonic frequencies at B3LYP/6-31G** level of theory for both stable cis-enol isomers. The calculated frequencies, Raman and IR intensities, and depolarization ratios were compared with the experimental results. The energy difference between the two stable cis-enol forms, calculated at B3LYP/6-311++G**, is only 5.89 kJ/mol. The observed vibrational frequencies and Raman and IR intensities are in excellent agreement with the corresponding values calculated for the most stable conformation, 2TFAA. According to the theoretical calculations, the hydrogen bond strength for the most stable conformer is 57 kJ/mol, about 9.5kJ/mol less than that of AA and about 14.5 kJ/mol more than that of HFAA. These hydrogen bond strengths are consistent with the frequency shifts for OH/OD stretching and OH/OD out-of-plane bending modes upon substitution of CH(3) groups with CF(3) groups. By comparing the vibrational spectra of both theoretical and experimental data, it was concluded that 2TFAA is the dominant isomer.

Experimental and Theoretical UV Characterizations of Acetylacetone and Its Isomers

Cryogenic matrix isolation experiments have allowed the measurement of the UV absorption spectra of the high-energy non-chelated isomers of acetylacetone, these isomers being produced by UV irradiation of the stable chelated form. Their identification has been done by coupling selective UV-induced isomerization, infrared spectroscopy, and harmonic vibrational frequency calculations using density functional theory. The relative energies of the chelated and non-chelated forms of acetylacetone in the S0 state have been obtained using density functional theory and coupled-cluster methods. For each isomer of acetylacetone, we have calculated the UV transition energies and dipole oscillator strengths using the excited-state coupled-cluster methods, including EOMCCSD (equation-of-motion coupled-cluster method with singles and doubles) and CR-EOMCCSD(T) (the completely renormalized EOMCC approach with singles, doubles, and non-iterative triples). For dipole-allowed transition energies, there is a very good agreement between experiment and theory. In particular, the CR-EOMCCSD(T) approach explains the blue shift in the electronic spectrum due to the formation of the non-chelated species after the UV irradiation of the chelated form of acetylacetone. Both experiment and CR-EOMCCSD(T) theory identify two among the seven non-chelated forms to be characterized by red-shifted UV transitions relative to the remaining five non-chelated isomers.

Proton tunneling estimates for malonaldehyde vibrations from supersonic jet and matrix quenching experiments

FTIR tunneling splittings for a range of fundamental excitations of malonaldehyde are determined from adiabatic cooling and symmetry breaking experiments in three different molecular states (isolated, coated by Ar layers and embedded in bulk Ar matrices), showing that there is room for improvement in available theoretical models.