Proton transfer in malonaldehyde: From reaction path to Schrödinger’s Cat

Investigating proton tunneling dynamics in the time-dependent Schrödinger equation

Understanding the temporal evolution of the wave function in the proton tunneling reactions allows us to make theoretical predictions on the possible femtosecond spectroscopy patterns. However, the analytical solution of the time-dependent Schrödinger equation of complex molecules is challenging and requires a high computational cost. In this study, we solve the time-dependent Schrödinger equation using the Fourier grid Hamiltonian method, highlighting its simplicity of calculation, even for multidimensional tunneling reactions. The obtained model is applied to studying malonaldehyde's two-dimensional intramolecular proton tunneling, comparing the results with those obtained using other computational methods.

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.

Selective photoisomerisation of 2-chloromalonaldehyde

Isomerization of 2-chloromalonaldehyde (2-ClMA) is explored giving access to new experimental data on this derivative of malonaldehyde, not yet studied much. Experiments were performed isolating 2-ClMA in argon, neon, and para-hydrogen matrices. UV irradiation of the matrix samples induced isomerization to three open enolic forms including two previously observed along with the closed enolic form after deposition. IR spectra of these specific conformers were recorded, and a clear assignment of the observed bands was obtained with the assistance of theoretical calculations. UV spectra of the samples were measured, showing a blue shift of the π^* ← π absorption with the opening of the internal hydrogen bond of the most stable enol form. Specific sequences of UV irradiation at different wavelengths allowed us to obtain samples containing only one enol conformer. The formation of conformers is discussed. The observed selectivity of the process among the enol forms is analyzed.

Asymmetric substitution changes the UV-induced nonradiative decay pathway and the spectra behaviors of β-diketones

Asymmetric substitution has not been termed as an essential factor in studying photo-induced ultrafast dynamics of molecular system. Asymmetric 4-hydroxybut-3-en-2-one (HEO), together with symmetric malonaldehyde (MA) and acetylacetone (AA), have been provided as target sample to study the nonradiative decay (ND) processes of β-diketones. An effective ND pathway of the three molecules is presented that their excited second (S₂) states transfer to first (S₁) state by nonadiabatic surface hopping, and then transfer to triplet (T₁) state by crossing minimum energy crossing point (MECP), after which decay to ground (S₀) state through MECP. More importantly, the asymmetric substitution of HEO induces the proton transfer in the S₁ state and generates a proton-transferred conformer with lowest energy, which does not occur for MA and AA. This change exploits a new ND pathway that the S₁ state decays to the proton transferred T₁ state and then undergoes reverse proton transfer to S₀ state through the MECPs between the three states. The two pathways of HEO give detailed energy and geometric information on surface hopping of S₂/S₁ and MECPs of S₁/T₁/S₀, and interpret the reason of the ND pathway while not spectra emission. This result is significantly different from the previous reported ND pathway of photoisomerization or conical intersection between different states. This work shows that asymmetric substitution changes the molecular structure and then changes their spectra behaviors.

Monitoring transient events in infrared spectra using local mode analysis

Time-resolved Fourier transformed infrared (FTIR) spectroscopy of chemical reactions is highly sensitive to minimal spatiotemporal changes. Structural features are decoded and represented in a comprehensible manner by combining FTIR spectroscopy with biomolecular simulations. Local mode analysis (LMA) is a tool to connect molecular motion based on a quantum mechanics simulation with infrared (IR) spectral features and vice versa. Here, we present the python-based software tool of LMA and demonstrate the novel feature of LMA to extract transient structural details and identify the related IR spectra at the case example of malonaldehyde (MA). Deuterated MA exists in two almost equally populated tautomeric states separated by a low barrier for proton transfer so IR spectra represent a mixture of both states. By state-dependent LMA, we obtain pure spectra for each tautomeric state occurring within the quantum mechanics trajectory. By time-resolved LMA, we obtain a clear view of the transition between states in the spectrum. Through local mode decomposition and the band-pass filter, marker bands for each state are identified. Thus, LMA is beneficial to analyze the experimental spectra based on a mixture of states by determining the individual contributions to the spectrum and motion of each state.

2-Chloromalonaldehyde, a model system of resonance-assisted hydrogen bonding: vibrational investigation

The chelated enol isomer of 2-chloromalonaldehyde (2-ClMA) is experimentally characterized for the first time by IR and Raman spectroscopies. The spectra are obtained by trapping the molecule in cryogenic matrices and analyzed with the assistance of theoretical calculations. Experiments were performed in argon, neon and para-hydrogen matrices. The results highlight puzzling matrix effects, beyond site effects, which are interpreted as due to a tunneling splitting of the vibrational levels related to the proton transfer along the internal hydrogen bond (IHB). 2-ClMA is thus one of the very few molecules in which the H tunneling has been observed in cryogenic matrices. The comparison with its parent molecule (malonaldehyde) shows experimentally and theoretically the weakening of the IHB upon chlorination, with a reduced cooperative effect in the resonance assisted hydrogen bond. In addition, the Cl substitution induces an important stabilization of two open enol conformers. These two open forms appear in the spectra of as-deposited samples, meaning that, in contrast with other well-studied molecules of the same family (β-dialdehydes and β-diketones), they are present in the gas phase at room temperature.

Muon-Substituted Malonaldehyde: Transforming a Transition State into a Stable Structure by Isotope Substitution

Isotope substitutions are usually conceived to play a marginal role on the structure and bonding pattern of molecules. However, a recent study [Angew. Chem. Int. Ed. 2014, 53, 13706-13709; Angew. Chem. 2014, 126, 13925-13929] further demonstrates that upon replacing a proton with a positively charged muon, as the lightest radioisotope of hydrogen, radical changes in the nature of the structure and bonding of certain species may take place. The present report is a primary attempt to introduce another example of structural transformation on the basis of the malonaldehyde system. Accordingly, upon replacing the proton between the two oxygen atoms of malonaldehyde with the positively charged muon a serious structural transformation is observed. By using the ab initio nuclear-electronic orbital non-Born-Oppenheimer procedure, the nuclear configuration of the muon-substituted species is derived. The resulting nuclear configuration is much more similar to the transition state of the proton transfer in malonaldehyde rather than to the stable configuration of malonaldehyde. The comparison of the "atoms in molecules" (AIM) structure of the muon-substituted malonaldehyde and the AIM structure of the stable and the transition-state configurations of malonaldehyde also unequivocally demonstrates substantial similarities of the muon-substituted malonaldehyde to the transition state.

New QM/MM implementation of the DFTB3 method in the gromacs package

The approximate density-functional tight-binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle-mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free-energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron-induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization.

Symmetry breaking in the axial symmetrical configurations of enolic propanedial, propanedithial, and propanediselenal: pseudo Jahn–Teller effect versus the resonance-assisted hydrogen bond theory

The origin of the symmetry breaking in the axial symmetrical configurations of enolic propanedial (1), propanedithial (2), and propanediselenal (3) have been investigated by means of time-dependence density functional theory and natural bond orbital interpretations. The results obtained at the quantum chemistry composite (G2MP2, CBS-QB3), ab initio molecular orbital (MP2/6-311++G**), and hybrid density functional theory (B3LYP/6-311++G**) levels of theory showed that the hydrogen-centered synchronous axial symmetrical (C2v) configurations of compounds 1–3 possessing the maximum π-electron delocalization within the M1=C2–C3=C4–M5–H6 keto-enol groups are less stable than their corresponding plane symmetrical (Cs) forms. Importantly, the symmetry breaking in the C2v configurations of the enol forms of compounds 1–3 to their corresponding plane symmetrical Cs configurations is due to the pseudo Jahn–Teller effect (PJTE) by mixing the ground A1 and excited B2 electronic states resulting in a PJT (A1 + B2) ⊗ b2 problem. We may expect that by the decrease of the energy gaps between reference states in the C2v forms that are involved in the PJTE decrease from compound 1 to compound 3, the PJT stabilization energy (PJTSE) may increase but the results obtained showed that the corresponding PJTSEs decrease. This fact can be justified by the increase of the electron delocalizations from the nonbonding orbitals of the C=M moieties to the antibonding orbitals of the H–M bonds, which leads to an increase of the π-electron delocalization within the M1=C2–C3=C4–M5–H6 keto-enol groups. In confrontation between the impacts of the resonance-assisted hydrogen bond and PJTE in the structural and configurational properties of compounds 1–3, PJTE has an overwhelming contribution and causes the symmetry breaking of the C2v configurations to their corresponding Cs forms. The correlations between the structural parameters, synchronicity indices, natural charges, PJTSEs, electron delocalizations, and the hardness of compounds 1–3 have been investigated.