Conformers of diheteroaryl ketones and thioketones: a quantum chemical study of their properties and fundamental intramolecular energetic effects

Theoretical insights into the spectroscopic properties of ferrocenyl hetaryl ketones

Quantum chemical calculations have been carried out to elucidate the electronic structure as well as to draw structure-property relationships for a series of ferrocenyl hetaryl ketones by means of simulated NMR, IR and UV-vis spectra. In this series, the list of hetaryl groups included furan-2-yl, thiophen-2-yl, selenophen-2-yl, 1H-pyrrol-2-yl and N-methylpyrrol-2-yl. Density functional theory was employed to determine the ground-state properties of the five ketones while their excited-state properties were modeled using a broad range of theoretical methods, namely from time-dependent density functional theory to multiconfigurational and multireference ab initio approaches. The patterns in the ¹³C and ¹⁷O chemical shifts of the carbonyl group were explained by the geometrical twist of hetaryl rings and by the electronic parameters corresponding to π-bonds conjugation and group hardness. Furthermore, the corresponding ¹³C and ¹⁷O shielding constants were analyzed in terms of both their dia/paramagnetic and Lewis/non-Lewis contributions within the framework of natural chemical shielding theory. The pattern in the vibrational frequency of the carbonyl bond was connected with changes in its bond length and bond order. It was established that the electronic absorption spectra of the studied ketones are largely characterized by low-intensity d → π* transitions in the visible region and the dominant high-intensity π → π* transition in the UV region. Finally, the theoretical methods best suited for modeling the excited-state properties of such ketones were designated.

Interacting Quantum Atoms-A Review

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

Structure of Diferrocenyl Thioketone: From Molecule to Crystal

Ferrocenyl-functionalized thioketones have recently been recognized as useful building blocks for sulfur-containing compounds with potential applications in materials chemistry. This work is devoted to a single representative of such thioketones, namely diferrocenyl thioketone (Fc₂CS), whose structure has been determined here for the first time. Both X-ray crystallography and a wide variety of quantum-chemical methods were used to explore the structure of Fc₂CS. In addition to the X-ray structure determination, intermolecular interactions occurring in the crystal structure of Fc₂CS were examined in detail by quantum-chemical methods. These methods were also an invaluable tool in studying the molecular structure of Fc₂CS, from the gas phase to solutions and to its crystal. Intramolecular interactions governing the conformational behavior of an isolated Fc₂CS molecule were deduced from quantum-chemical analyses carried out in orbital space and real space. Our experimental and theoretical results indicate that the main structural features of an isolated Fc₂CS molecule in its lowest-energy geometry are retained both upon solvation and in the crystal. The tilt of ferrocenyl groups is only slightly affected by crystal packing forces that are dominated by dispersion. Nonetheless, a network of intermolecular interactions, such as H···H, C···H and S···H, was detected in the Fc₂CS crystal but each of them is fairly weak.

Does the Intra-Atomic Deformation Energy of Interacting Quantum Atoms Represent Steric Energy?

We show that the mutual, through-space compression of atomic volume experienced by approaching topological atoms causes an exponential increase in the intra-atomic energy of those atoms, regardless of approach orientation. This insight was obtained using the modern energy partitioning method called interacting quantum atoms (IQA). This behaviour is consistent for all atoms except hydrogen, which can behave differently depending on its environment. Whilst all atoms experience charge transfer when they interact, the intra-atomic energy of the hydrogen atom is more vulnerable to these changes than larger atoms. The difference in behaviour is found to be due to hydrogen's lack of a core of electrons, which, in heavier atoms, consistently provide repulsion when compressed. As such, hydrogen atoms do not always provide steric hindrance. In accounting for hydrogen's unusual behaviour and demonstrating the exponential character of the intra-atomic energy in all other atoms, we provide evidence for IQA's intra-atomic energy as a quantitative description of steric energy.

An interacting quantum atom study of model SN 2 reactions (X- ···CH3 X, X = F, Cl, Br, and I)

The quantum chemical topology method has been used to analyze the energetic profiles in the X-  + CH3 X → XCH3  + X- SN 2 reactions, with X = F, Cl, Br, and I. The evolution of the electron density properties at the BCPs along the reaction coordinate has been analysed. The interacting quantum atoms (IQA) method has been used to evaluate the intra-atomic and interatomic energy variations along the reaction path. The different energetic terms have been examined by the relative energy gradient method and the ANANKE program, which enables automatic and unbiased IQA analysis. Four of the six most important IQA energy contributions were needed to reproduce the reaction barrier common to all reactions. The four reactions considered share many common characteristics but when X = F a number of particularities occur. © 2017 Wiley Periodicals, Inc.

Extension of the interacting quantum atoms (IQA) approach to B3LYP level density functional theory (DFT)

An interaction between two atoms, bonded or non-bonded, consists of interatomic contributions: electrostatic energy, exchange energy and electronic correlation energy. Together with the intra-atomic energy of an atom, these contributions are the basic components of the Interacting Quantum Atom (IQA) energy decomposition scheme. Here, we investigate IQA's proper use in conjunction with an explicit implementation of the B3LYP functional. The recovery of the total molecular energy from the IQA components is emphasised, for the first time. A systematic study of three model systems of biological relevance, N-methylacetamide (NMA), the doubly capped tripeptide GlyGlyGly and an alloxan dimer, shows the stabilization effect of B3LYP on most of the interatomic exchange energies (V) compared to their Hartree-Fock values. Diagrams of exchange energies versus interatomic distance show the clustering of interactions, one cluster for each 1,n (n = 1 to 6 where the atoms are separated by n - 1 bonds). The positioning of some V values outside their expected cluster marks interesting interactions.

Heteroatom and solvent effects on molecular properties of formaldehyde and thioformaldehyde symmetrically disubstituted with heterocyclic groups C4H3Y (where Y = O-Po)

In this work several molecular properties of symmetrically disubstituted formaldehyde and thioformaldehyde have been studied using a quantum chemistry approach based on density functional theory. Five-membered heteroaromatic rings containing a single group 16 heteroatom were taken into account as the substituents (i.e., furan-2-yl, thiophen-2-yl, selenophen-2-yl, tellurophen-2-yl, and the experimentally as yet unknown polonophen-2-yl). For the resulting ten formaldehyde and thioformaldehyde derivatives, the geometry, energetics, frontier molecular orbitals, dipole moment and polarizability of their molecules were examined in order to establish the effect of ring heteroatom on these properties. Furthermore, these properties were also determined for the molecules in three solvents of low polarity (benzene, chloroform, and dichloromethane) in order to expand our study to include solvent effects. The dipole moment and polarizability of the investigated molecules show regular variations when the ring heteroatom descends through group 16 and the solvent polarity grows. The heteroatom and/or solvent effects on the part of the studied properties are, however, more complex. An attempt is made to rationalize the observed variations in the molecular properties. The conformational behavior of the investigated molecules was also explored and the conformationally weighted values of dipole moment and polarizability are presented.

Graphical abstract

The ANANKE relative energy gradient (REG) method to automate IQA analysis over configurational change

Much chemical insight ultimately comes down to finding out which fragment of a total system behaves like the total system, in terms of an energy profile. A simple example is that of the water dimer, where this system is regarded as held together by a hydrogen bond. The hydrogen bond consists of two atoms (H···O), which energetically behave similarly to the total system (H₂O)₂. However, from a quantum mechanical point of view, each atom in the total system interacts with any other atom. Thus, the view that the hydrogen bond by itself governs the energetic stability of the water dimer needs rigorous justification. In this work, we propose a method that provides such a justification, in general, but only illustrated on the water dimer here. This method is based on the topological energy partitioning method called interacting quantum atoms (IQA). The method is implemented in the program ANANKE, which calculates correlations between the energy profile of the total system and those of subsystems (or fragments). ANANKE acts on the IQA energy contributions obtained for a sequence of full-system geometries controlled by a coordinate of interest (e.g. the O···H distance in the water dimer). Although applied only for the water dimer in this work, the method is general and able to explain the gauche effect, the torsional barrier in biphenyl, the arrow-pushing scheme of an enzymatic reaction (peptide hydrolysis in the HIV-1 Protease active site), and halogen-alkane nucleophilic substitution (S_N2) reactions. Those applications will appear elsewhere as separate and elaborated case studies; here we focus on the details of the ANANKE method and its justification, using the water dimer as a concrete case.