Bonding analysis of water clusters using quasi-atomic orbitals

The quasi-atomic orbital (QUAO) bonding analysis introduced by Ruedenberg and co-workers is used to develop an understanding of the hydrogen bonds in small water clusters, from the dimer through the hexamer (bag, boat, book, cyclic, prism and cage conformers). Using kinetic bond orders as a metric, it is demonstrated that as the number of waters in simple cyclic clusters increases, the hydrogen bonds strengthen, from the dimer through the cyclic hexamer. However, for the more complex hexamer isomers, the strength of the hydrogen bonds varies, depending on whether the cluster contains double acceptors and/or double donors. The QUAO analysis also reveals the three-center bonding nature of hydrogen bonds in water clusters.

3-D bond-paths of QTAIM and the stress tensor in neutral lithium clusters, Lim (m = 2-5), presented on the Ehrenfest force molecular graph

In this investigation we set out to understand the origins of non-nuclear attractors (NNAs) found for neutral lithium clusters Lim (m = 2-5) on the QTAIM molecular graph but not on the Ehrenfest force F(r) molecular graph. Therefore, we pursued the stress tensor σ(r) without using the dependency on the QTAIM partitioning, since previously σ(r) was only calculated within the QTAIM partitioning, to see if any indication of NNA character can be determined. Because the stress tensor σ(r) lacks an associated scalar- or vector-field as is the case for QTAIM and the Ehrenfest F(r) partitioning schemes respectively, a stress tensor σ(r) partitioning scheme cannot be constructed. Therefore, to overcome this difficulty we use next generation QTAIM, constructed from the most preferred directions of electronic charge density accumulation, to calculate the stress tensor σ(r) 3-D bond-paths on the Ehrenfest force F(r) molecular graph. Using next generation 3-D bond-paths within the Ehrenfest force F(r) partitioning, we can classify the degree of NNA character in the absence of NNAs. A much higher degree of NNA character is found to be present for the stress tensor σ(r) 3-D bond-paths than for the corresponding QTAIM or Ehrenfest force F(r) 3-D bond-paths. The stabilizing effect of the NNA is demonstrated by undertaking Li2 bond-path compression and stretching distortions sufficient to cause the annihilation of the NNA. The compression and stretching distortions also lead to a large increase in the 3-D bond-path asymmetry and persistent bond-path torsion respectively.

Non-nuclear attractors in small charged lithium clusters, Limq (m = 2-5, q = ±1), with QTAIM and the Ehrenfest force partitioning

In this investigation we explore the function and existence of the non-nuclear attractor (NNA) for a series of small charged lithium clusters Limq (m = 2-5, q = ±1) using QTAIM and the Ehrenfest force F(r) partitioning schemes. The NNAs were found to be present in all of the Limq (m = 2-5, q = ±1) clusters for QTAIM, in contrast none were found for F(r). We discovered that the anionic and cationic lithium dimers are limiting cases for minimal and maximal impact of the NNA related to the relative sparseness of total charge density ρ(r) distributions respectively. Evidence is found that the NNA in the anionic dimer is in the process of being annihilated by two neighboring BCPs. We provide a measure of the size of the NNA and find for Limq (m = 2-5, q = ±1) that larger NNAs correlate with increased Li-Li separations. The NNA was determined to be a persistent feature by varying the Li separations for the cationic and anionic dimers. Very large Li separations failed to induce an NNA in the F(r) anionic dimer and therefore we conclude that F(r) is unable to detect NNAs. The metallicity ξ(rb) was also used to measure the sparseness of the distribution of ρ(r) and significant metallic character, on the basis of ξ(rb) > 1, was present for QTAIM but not for F(r), providing further evidence that F(r) cannot detect NNAs. Advantages of the use of Ehrenfest force F(r) partitioning scheme are discussed that include the design of nano-devices through tuning of the Ehrenfest potential VF(b) by the application of external forces such as a constant electric or strain field.

Exploration of the forbidden regions of the Ramachandran plot (ϕ-ψ) with QTAIM

A new QTAIM interpretation of the Ramachandran plot is formulated from the most and least facile eigenvectors of the second-derivative matrix of the electron density with a set of 29 magainin-2 peptide conformers. The presence of QTAIM eigenvectors associated with the most and least preferred directions of electronic charge density explained the role of hydrogen bonding, HH contacts and the glycine amino acid monomer in peptide folding. The highest degree of occupation of the QTAIM interpreted Ramachandran plot was found for the glycine amino acid monomer compared with the remaining backbone peptide bonds. The mobility of the QTAIM eigenvectors of the glycine amino acid monomer was higher than for the other amino acids and was comparable to that of the hydrogen bonding, explaining the flexibility of the magainin-2 backbone. We experimented with a variety of hybrid QTAIM-Ramachandran plots to highlight and explain why the glycine amino acid monomer largely occupies the 'forbidden' region on the Ramachandran plot. In addition, the new hybrid QTAIM-Ramachandran plots contained recognizable regions that can be associated with concepts familiar from the conventional Ramachandran plot whilst retaining the character of the QTAIM most and least preferred regions.

Fe3+chelating quinoline–hydrazone hybrids with proven cytotoxicity, leishmanicidal, and trypanocidal activities

Neglected tropical diseases cause great concern in developing countries where there are millions of reported infected humans. Our calculations support a direct relationship between biological activity and the Fe³⁺chelating ability of the shown set of quinoline–hydrazone hybrids.

Why do antifreeze proteins require a solenoid?

Proteins whose presence prevents water from freezing in living organisms at temperatures below 0 °C are referred to as antifreeze proteins. This group includes molecules of varying size (from 30 to over 300 aa) and variable secondary/supersecondary conformation. Some of these proteins also contain peculiar structural motifs called solenoids. We have applied the fuzzy oil drop model in the analysis of four categories of antifreeze proteins: 1 - very small proteins, i.e. helical peptides (below 40 aa); 2 - small globular proteins (40-100 aa); 3 - large globular proteins (>100 aa) and 4 - proteins containing solenoids. The FOD model suggests a mechanism by which antifreeze proteins prevent freezing. In accordance with this theory, the presence of the protein itself produces an ordering of water molecules which counteracts the formation of ice crystals. This conclusion is supported by analysis of the ordering of hydrophobic and hydrophilic residues in antifreeze proteins, revealing significant variability - from perfect adherence to the fuzzy oil drop model through structures which lack a clearly defined hydrophobic core, all the way to linear arrangement of alternating local minima and maxima propagating along the principal axis of the solenoid (much like in amyloids). The presented model - alternative with respect to the ice docking model - explains the antifreeze properties of compounds such as saccharides and fatty acids. The fuzzy oil drop model also enables differentiation between amyloids and antifreeze proteins.

QTAIM and Stress Tensor Characterization of Intramolecular Interactions Along Dynamics Trajectories of a Light-Driven Rotary Molecular Motor

A quantum theory of atoms in molecules (QTAIM) and stress tensor analysis was applied to analyze intramolecular interactions influencing the photoisomerization dynamics of a light-driven rotary molecular motor. For selected nonadiabatic molecular dynamics trajectories characterized by markedly different S₁ state lifetimes, the electron densities were obtained using the ensemble density functional theory method. The analysis revealed that torsional motion of the molecular motor blades from the Franck-Condon point to the S₁ energy minimum and the S₁/S₀ conical intersection is controlled by two factors: greater numbers of intramolecular bonds before the hop-time and unusually strongly coupled bonds between the atoms of the rotor and the stator blades. This results in the effective stalling of the progress along the torsional path for an extended period of time. This finding suggests a possibility of chemical tuning of the speed of photoisomerization of molecular motors and related molecular switches by reshaping their molecular backbones to decrease or increase the degree of coupling and numbers of intramolecular bond critical points as revealed by the QTAIM/stress tensor analysis of the electron density. Additionally, the stress tensor scalar and vector analysis was found to provide new methods to follow the trajectories, and from this, new insight was gained into the behavior of the S₁ state in the vicinity of the conical intersection.

A QTAIM and stress tensor investigation of the torsion path of a light-driven fluorene molecular rotary motor

The utility of the QTAIM/stress tensor analysis method for characterizing the photoisomerization of light driven molecular rotary machines is investigated on the example of the torsion path in fluorene molecular motor. The scalar and vector descriptors of QTAIM/stress tensor reveal additional information on the bonding interactions between the rotating units of the motor, which cannot be obtained from the analysis of the ground and excited state potential energy surfaces. The topological features of the fluorene motor molecular graph display that, upon the photoexcitation a certain increase in the torsional stiffness of the rotating bond can be attributed to the increasing topological stability of the rotor carbon atom attached to the rotation axle. The established variations in the torsional stiffness of the rotating bond may cause transfer of certain fraction of the torsional energy to other internal degrees of freedom, such as the pyramidalization distortion. © 2016 Wiley Periodicals, Inc.

Microsolvation of methylmercury: structures, energies, bonding and NMR constants (199Hg,13C and17O)

Hartree–Fock (HF) and second order perturbation theory (MP2) calculations within the scalar and full relativistic frames were carried out in order to determine the equilibrium geometries and interaction energies between cationic methylmercury (CH₃Hg⁺) and up to three water molecules.

Unusual bonding modes of perfluorobenzene in its polymeric (dimeric, trimeric and tetrameric) forms: entirely negative fluorine interacting cooperatively with entirely negative fluorine

The Laplacian of the charge density plot for the tetramer of perfluorobenzene, displaying the attraction between entirely negative fluorine atoms, as well as the windmill type intermolecular topology between those atoms.

Hybrid QTAIM and electrostatic potential-based quantum topology phase diagrams for water clusters

The topological diversity of sets of isomers of water clusters (W = H₂O)_n, 7 ≤ n ≤ 10, is analyzed employing the scalar fields of total electronic charge density ρ(r) and the molecular electrostatic potential (MESP).

Advances in understanding of chemical bonding: inputs from experimental and theoretical charge density analysis

The development of charge density analysis has undergone a major renaissance in the last two decades. In recent years, the characterization of bonding features associated with atoms in molecules and in crystals has been explored using high-resolution X-ray diffraction data (laboratory or synchrotron) complemented by high level ab initio theoretical calculations. The extraction of one electron topological properties, namely, electrostatic charges, dipole moment and higher moments, electrostatic potential, electric field gradients, in addition to evaluation of the local kinetic and potential energy densities, have contributed toward an understanding of the electron density distributions in molecular solids. New topological descriptors, namely, the source function (SF) and electron localization function (ELF) provide additional information as regards characterization of the topology of the electron density. In addition, delocalization indices have also been developed to account for bonding features pertinent to M-M bonds. The evaluation of these properties have contributed significantly toward the understanding of intra- and intermolecular bonding features in organic, inorganic, and biomolecules in the crystalline phase, with concomitant applications in the understanding of chemical reactivity and material/biological properties. In recent years, the focus has strongly shifted toward the understanding of structure-property relationships in organometallic complexes containing labile M-C bonds in the crystal structure with subsequent implications in catalysis. This perspective aims to highlight the major developments in electron density measurements in the past few years and provides pointers directed toward the potential use of this technique in future applications for an improved understanding of chemical bonding in systems that have been unexplored.