Imaging Momentum Orbital Densities of Conformationally Versatile Molecules: A Benchmark Theoretical Study of the Molecular and Electronic Structures of Dimethoxymethane

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Influence of Membrane Phase on the Optical Properties of DPH

The fluorescent molecule diphenylhexatriene (DPH) has been often used in combination with fluorescence anisotropy measurements, yet little is known regarding the non-linear optical properties. In the current work, we focus on them and extend the application to fluorescence, while paying attention to the conformational versatility of DPH when it is embedded in different membrane phases. Extensive hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the influence of the phase- and temperature-dependent lipid environment on the probe. Already, the transition dipole moments and one-photon absorption spectra obtained in the liquid ordered mixture of sphingomyelin (SM)-cholesterol (Chol) (2:1) differ largely from the ones calculated in the liquid disordered DOPC and solid gel DPPC membranes. Throughout the work, the molecular conformation in SM:Chol is found to differ from the other environments. The two-photon absorption spectra and the ones obtained by hyper-Rayleigh scattering depend strongly on the environment. Finally, a stringent comparison of the fluorescence anisotropy decay and the fluorescence lifetime confirm the use of DPH to gain information upon the surrounding lipids and lipid phases. DPH might thus open the possibility to detect and analyze different biological environments based on its absorption and emission properties.

Tuning the HOMO-LUMO Energy Gap of Small Diamondoids Using Inverse Molecular Design

Functionalized diamondoids show great potential as building blocks for various new optoelectronic applications. However, until now, only simple mono and double substitutions were investigated. In this work, we considered up to 10 and 6 sites for functionalization of the two smallest diamondoids, adamantane and diamantane, respectively, in search for diamondoid derivatives with a minimal and maximal HOMO-LUMO energy gap. To this end, the energy gap was optimized systematically using an inverse molecular design methodology based on the best-first search algorithm combined with a Monte Carlo component to escape local optima. Adamantane derivatives were found with HOMO-LUMO gaps ranging from 2.42 to 10.63 eV, with 9.45 eV being the energy gap of pure adamantane. For diamantane, similar values were obtained. The structures with the lowest HOMO-LUMO gaps showed apparent push-pull character. The push character is mainly formed by sulfur or nitrogen dopants and thiol groups, whereas the pull character is predominantly determined by the presence of electron-withdrawing nitro or carbonyl groups assisted by amino and hydroxyl groups via the formation of intramolecular hydrogen bonds. In contrast, maximal HOMO-LUMO gaps were obtained by introducing numerous electronegative groups.

Electron Momentum Spectroscopy Investigation of Molecular Conformations of Ethanol Considering Vibrational Effects

The interpretation of experimental electron momentum distributions (EMDs) of ethanol, one of the simplest molecules having conformers, has confused researchers for years. High-level calculations of Dyson orbital EMDs by thermally averaging the gauche and trans conformers as well as molecular dynamical simulations failed to quantitatively reproduce the experiments for some of the outer valence orbitals. In this work, the valence shell electron binding energy spectrum and EMDs of ethanol are revisited by the high-sensitivity electron momentum spectrometer employing symmetric noncoplanar geometry at an incident energy of 1200 eV plus binding energy, together with a detailed analysis of the influence of vibrational motions on the EMDs for the two conformers employing a harmonic analytical quantum mechanical (HAQM) approach by taking into account all of the vibrational modes. The significant discrepancies between theories and experiments in previous works have now been interpreted quantitatively, indicating that the vibrational effect plays a significant role in reproducing the experimental results, not only through the low-frequency OH and CH₃ torsion modes but also through other high-frequency ones. Rational explanation of experimental momentum profiles provides solid evidence that the trans conformer is slightly more stable than the gauche conformer, in accordance with thermodynamic predictions and other experiments. The case of ethanol demonstrates the significance of considering vibrational effects when performing a conformational study on flexible molecules using electron momentum spectroscopy.

Electron momentum spectroscopy of dimethyl ether taking account of nuclear dynamics in the electronic ground state

The influence of nuclear dynamics in the electronic ground state on the (e,2e) momentum profiles of dimethyl ether has been analyzed using the harmonic analytical quantum mechanical and Born-Oppenheimer molecular dynamics approaches. In spite of fundamental methodological differences, results obtained with both approaches consistently demonstrate that molecular vibrations in the electronic ground state have a most appreciable influence on the momentum profiles associated to the 2b1, 6a1, 4b2, and 1a2 orbitals. Taking this influence into account considerably improves the agreement between theoretical and newly obtained experimental momentum profiles, with improved statistical accuracy. Both approaches point out in particular the most appreciable role which is played by a few specific molecular vibrations of A1, B1, and B2 symmetries, which correspond to C-H stretching and H-C-H bending modes. In line with the Herzberg-Teller principle, the influence of these molecular vibrations on the computed momentum profiles can be unraveled from considerations on the symmetry characteristics of orbitals and their energy spacing.

Theoretical study of molecular vibrations in electron momentum spectroscopy experiments on furan: an analytical versus a molecular dynamical approach

The influence of thermally induced nuclear dynamics (molecular vibrations) in the initial electronic ground state on the valence orbital momentum profiles of furan has been theoretically investigated using two different approaches. The first of these approaches employs the principles of Born-Oppenheimer molecular dynamics, whereas the so-called harmonic analytical quantum mechanical approach resorts to an analytical decomposition of contributions arising from quantized harmonic vibrational eigenstates. In spite of their intrinsic differences, the two approaches enable consistent insights into the electron momentum distributions inferred from new measurements employing electron momentum spectroscopy and an electron impact energy of 1.2 keV. Both approaches point out in particular an appreciable influence of a few specific molecular vibrations of A1 symmetry on the 9a1 momentum profile, which can be unravelled from considerations on the symmetry characteristics of orbitals and their energy spacing.

Ring-puckering effects on electron momentum distributions of valence orbitals of oxetane

The binding energy spectra and electron momentum distributions for the outer-valence molecular orbitals of oxetane have been measured utilizing (e, 2e) electron momentum spectrometer with non-coplanar asymmetric geometry at the impact energy of 2500 eV. The experimental momentum distributions were compared with the density functional theory calculations employing B3LYP hybrid functional with aug-cc-pVTZ basis set. It was found that the calculation at planar geometry (C2v) completely fails to interpret the large "turn-up" at low momentum region in electron momentum distribution of the highest occupied molecular orbital (HOMO) 3b1, while the calculations considering the thermal abundances of planar (C2v) and bent (Cs) conformers or the thermally populated vibrational states of ring-puckering motion have significantly improved the agreement. The results indicate that the ring-puckering motion of oxetane has a strong effect on the electron density distribution of HOMO.

Electron momentum spectroscopy of 1-butene: a theoretical analysis using molecular dynamics and molecular quantum similarity

The results of experimental studies of the valence electronic structure of 1-butene by means of electron momentum spectroscopy (EMS) have been reinterpreted on the basis of molecular dynamical simulations in conjunction with the classical MM3 force field. The computed atomic trajectories demonstrate the importance of thermally induced nuclear dynamics in the electronic neutral ground state, in the form of significant deviations from stationary points on the potential energy surface and considerable variations of the C-C-C-C dihedral angle. These motions are found to have a considerable influence on the computed spectral bands and outer-valence electron momentum distributions. Euclidean distances between spherically averaged electron momentum densities confirm that thermally induced nuclear motions need to be fully taken into account for a consistent interpretation of the results of EMS experiments on conformationally flexible molecules.

Photoelectron and electron momentum spectroscopy of tetrahydrofuran from a molecular dynamical perspective

The results of experimental studies of the valence electronic structure of tetrahydrofuran employing He I photoelectron spectroscopy as well as Electron Momentum Spectroscopy (EMS) have been reinterpreted on the basis of Molecular Dynamical simulations employing the classical MM3 force field and large-scale quantum mechanical simulations employing Born-Oppenheimer Molecular Dynamics in conjunction with the dispersion corrected ωB97XD exchange-correlation functional. Analysis of the produced atomic trajectories demonstrates the importance of thermal deviations from the lowest energy path for pseudorotation, in the form of considerable variations of the ring-puckering amplitude. These deviations are found to have a significant influence on several outer-valence electron momentum distributions, as well as on the He I photoelectron spectrum.