Diving into chemical bonding: An immersive analysis of the electron charge rearrangement through virtual reality

Local Charge-Displacement Analysis: Targeting Local Charge-Flows in Complex Intermolecular Interactions

Charge-displacement (CD) analysis has recently proven to be a simple and powerful scheme for quantitatively analyzing the profile of the charge redistribution occurring upon intermolecular interactions along a given interaction axis. However, when two molecular fragments bind through complex interactions involving multiple concurrent charge flows, ordinary CD analysis is capable of providing only an averaged picture of the related charge-flow profiles and no detailed information on each of them. In this article, we combine CD analysis with a Hirshfeld partitioning of the molecular charge redistribution for a local analysis on focused portions of the molecule, allowing for a detailed characterization of one charge flow at a time. The resulting scheme - the local charge-displacement (LCD) analysis - is tested on the intriguing case of the dimethyl sulfide (DMS)-sulfur dioxide (SO2) complex, characterized by concurrent charge flows relating to a sulfur-sulfur homo-chalcogen interaction and a pair of hydrogen bonds. The LCD scheme is then applied to the analysis of multiple hydrogen bonding in the acetic acid dimer, of base-pairing interactions in DNA, and of ambifunctional hydrogen bonding in the ammonia-pyridine complex.

Charge-Flow Profiles along Curvilinear Paths: A Flexible Scheme for the Analysis of Charge Displacement upon Intermolecular Interactions

The Charge-Displacement (CD) analysis has proven to be a powerful tool for a quantitative characterization of the electron-density flow occurring upon chemical bonding along a suitably chosen interaction axis. In several classes of interesting intermolecular interactions, however, an interaction axis cannot be straightforwardly defined, and the CD analysis loses consistency and usefulness. In this article, we propose a general, flexible reformulation of the CD analysis capable of providing a quantitative view of the charge displacement along custom curvilinear paths. The new scheme naturally reduces to ordinary CD analysis if the path is chosen to be a straight line. An implementation based on a discrete sampling of the electron densities and a Voronoi space partitioning is described and shown in action on two test cases of a metal-carbonyl and a pyridine-ammonia complex.

Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy

The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.

Combining Virtual Reality Visualization with Ensemble Molecular Dynamics to Study Complex Protein Conformational Changes

Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics, and their biological function. Currently, it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome "collective variable" enhanced sampling protocols. Here, we describe a framework that combines ensemble MD simulations and virtual reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables a priori. We further show that eMD-VR generated pathways can be combined with Markov state models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts.

Disentanglement of orthogonal hydrogen and halogen bonds via natural orbital for chemical valence: A charge displacement analysis

As known, the electron density of covalently bound halogen atoms is anisotropically distributed, making them potentially able to establish many weak interactions, acting at the same time as halogen bond donors and hydrogen bond acceptors. Indeed, there are many examples in which the halogen and hydrogen bond coexist in the same structure and, if a correct bond analysis is required, their separation is mandatory. Here, the advantages and limitations of coupling the charge displacement analysis with natural orbital for chemical valence method (NOCV-CD) to separately analyze orthogonal weak interactions are shown, for both symmetric and asymmetric adducts. The methodology gives optimal results with intermolecular adducts but, in the presence of an organometallic complex, also intramolecular interactions can be correctly analyzed. Beyond the methodological aspects, it is shown that correctly separate and quantify the interactions can give interesting chemical insights about the systems.

A General User-Friendly Tool for Kinetic Calculations of Multi-Step Reactions within the Virtual Multifrequency Spectrometer Project

We discuss the implementation of a computer program for accurate calculation of the kinetics of chemical reactions integrated in the user-friendly, multi-purpose Virtual Multifrequency Spectrometer tool. The program is based on the ab initio modeling of the involved molecular species, the adoption of transition-state theory for each elementary step of the reaction, and the use of a master-equation approach accounting for the complete reaction scheme. Some features of the software are illustrated through examples including the interconversion reaction of hydroxyacetone and 2-hydroxypropanal and the production of HCN and HNC from vinyl cyanide.

Chemical promenades: Exploring potential-energy surfaces with immersive virtual reality

The virtual-reality framework AVATAR (Advanced Virtual Approach to Topological Analysis of Reactivity) for the immersive exploration of potential-energy landscapes is presented. AVATAR is based on modern consumer-grade virtual-reality technology and builds on two key concepts: (a) the reduction of the dimensionality of the potential-energy surface to two process-tailored, physically meaningful generalized coordinates, and (b) the analogy between the evolution of a chemical process and a pathway through valleys (potential wells) and mountain passes (saddle points) of the associated potential energy landscape. Examples including the discovery of competitive reaction paths in simple A + BC collisional systems and the interconversion between conformers in ring-puckering motions of flexible rings highlight the innovation potential that augmented and virtual reality convey for teaching, training, and supporting research in chemistry.

Theoretical study of glycine amino acid adsorption on graphene oxide

The non-dissociative and dissociative adsorptions of zwitterionic Gly on graphene oxide (GO) was studied in the framework of DFT using a cluster model approach. In this work, the interaction with an epoxy group of GO basal plane was mainly considered. As a comparison, the non-dissociative and dissociative adsorptions of neutral Gly were also taken into account. The non-dissociative adsorption modes for zwitterionic and neutral Gly conformers show binding energies of 12.2 and 14.4 kcal mol^-1, respectively. These molecules are thought to remain over the GO surface due to attractive noncovalent interactions. Two dissociative adsorption modes, for Z-Gly and N-Gly, show smaller binding energies of 7.2 and 8.4 kcal mol^-1, where the deprotonated species links strongly through a C-O or C-N covalent bond to the GO surface. The results obtained in the present theoretical approach to the glycine/graphene oxide system support the fact that glycine can be attached to epoxy groups of graphene oxide basal planes in addition to the anchoring on edge oxidation groups. In summary, we conclude that glycine can be used as a reducing agent as well as a functionalizer of GO sheets.

Chemical bonding in cuprous complexes with simple nitriles: octet rule and resonance concepts versus quantitative charge-redistribution analysis

Resonance structures for six cuprous complexes with simple nitriles are interpreted by means of a quantitative analysis of charge redistribution upon copper-nitrile bonding.

Influence of halogen bonding on gold(i)-ligand bond components and DFT characterization of a gold-iodine halogen bond

A gold(i) complex bearing nitrogen acyclic carbene (NAC) and selenourea (SeU) has been used to verify whether the second-sphere SeI halogen bond (XB) is able to modify the Dewar-Chatt-Duncanson components of the Au-C and Au-Se bonds. The chosen system was found to be thermically unstable but it allowed an in-depth theoretical study by means of Energy Decomposition Analysis, Natural Bond Orbital and Natural Orbitals for Chemical Valence methods, coupled with Charge Displacement analysis. Indeed, in the presence of iodoperfluoroalkanes as XB donors, iodine interacts with the lone pair of the coordinated selenium, enhancing the Au ← C σ donation and depressing the Au → C π back-donation, as demonstrated also by the increase of the rotational barrier of the C-N bond of the NAC (see G. Ciancaleoni and others, Chem. - Eur. J., 2015, 21, 2467). On the other hand, in the presence of N-iodosuccinimide (NIS), the gold directly establishes a XB with the iodine by using its d lone pairs. This AuI XB is favored by the low steric hindrance of the ligands coordinated to the gold and the presence of the amino protons of SeU, which establish additional hydrogen bonds with the NIS. Also in this case, the effect is to increase the σ acidity and decrease the π basicity of the metal.

Potential-Energy Surfaces for Ring-Puckering Motions of Flexible Cyclic Molecules through Cremer-Pople Coordinates: Computation, Analysis, and Fitting

Ring-puckering motion in 12 flexible cyclic molecules is investigated by calculation and analysis of two-dimensional potential-energy surfaces (PESs) using the so-called ring-puckering coordinates proposed by Cremer and Pople. The PESs are calculated by means of density-functional theory using a B2PLYP-D3BJ exchange-correlation functional with a maug-cc-pVTZ basis set, and results are compared to the available experimental and theoretical data. Special care is devoted to the aspect of symmetry in such two-dimensional PESs, which are here reported for the first time also for molecules whose planar form has symmetry lower than D_{5 h} or C_{2 v}. The issue of PES fitting and that of solving the nuclear dynamics using ring-puckering coordinates are also addressed. Analytical formulations of the computed PESs using suitable functional forms with a limited set of parameters are provided.

Study of Commodity VR for Computational Material Sciences

Recent advancements in virtual reality (VR) devices and software environments make it possible to easily incorporate this technology for many applications, including computational materials science. For studying three-dimensional (3D) structure models and related chemical information, we focused on using a commodity VR device (VIVE) and an authoring tool (Unity). To visualize 3D chemical structures, disturbances like judder due to dropped frames should be eliminated from the VR experience to improve simulations. We propose a simple evaluation method that is straightforward for the nonexpert or novice VR user. We examine the major visualization representations including ball, ball and stick, and isosurface systems. For systematic benchmark measurements, a pendulum from the VR device was used to generate periodic oscillatory motion during measurements of a time series in frames per second (fps). For VIVE with a refresh rate of 90 Hz, judder occurred when less than 90 fps. We demonstrated the system size limitations for the results of molecular dynamics simulations of phase separation of ABA block copolymers and experimental observations of filler morphologies in rubber.

Interplay between Gold(I)-Ligand Bond Components and Hydrogen Bonding: A Combined Experimental/Computational Study

The influence of weak interactions on the donation/back-donation bond components in the complex [(NHC)Au(SeU)]⁺ (NHC = N-heterocyclic carbene; SeU = selenourea) has been studied by coupling experimental and theoretical techniques. In particular, NMR ¹H and pulsed-field gradient spin-echo titrations allowed us to characterize the hydrogen bond (HB) between the -NH₂ moieties of SeU and the anions PF₆ ^- and ClO₄ ^-, whereas ⁷⁷Se NMR spectroscopy allowed us to characterize the Au-Se bond. Theoretically, the Au-Se and Au-C orbital interactions have been decomposed using the natural orbital for the chemical valence framework and the bond components quantified through the charge displacement analysis. This methodology provides the quantification of the Dewar-Chatt-Duncanson (DCD) components for the Au-C and Au-Se bonds in the absence and presence of the second-sphere HB. The results presented here show that the anion has a dual mode action: it modifies the conformation of the cation by ion pairing (and this already influences the DCD components) and it induces new polarization effects that depend on the relative anion/cation relative orientation. The perchlorate polarizes SeU, enhancing the Se → Au σ donation and the Au → C back-donation and depressing the C → Au σ donation. On the contrary, the hexafluorophosphate depresses both the Se → Au and C → Au σ donations.

Ferrocenes with simple chiral substituents: an in-depth theoretical and experimental VCD and ECD study

Ferrocenes bearing chiral pendants are investigated through VCD and ECD. The VCD spectra are best interpreted by GVPT2-anharmonic DFT calculations. Diagnostic bands related to the absolute configuration of the title compounds are found in both kinds of spectra.

Vibrational circular dichroism under the quantum magnifying glass: from the electronic flow to the spectroscopic observable

A chemically intuitive method to analyse and interpret vibrational circular dichroism spectra based on the vibrational transition current density.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?