RESUMEN
This work aims at studying why azodicarbonamide (ADCA), a formally apolar compound with good hydrogen bond (HB) acceptors, is soluble only in polar aprotic solvents like dimethyl sulfoxide (DMSO) but not in water. Solubility measurements, as well as quantum mechanical and classical molecular dynamics simulations, were employed to tackle the problem. We found that in the liquid phase a polar conformer of ADCA (µ = 8.7 D), unreported to date, is favoured under the enthalpic drive provided by a highly polar solvent. At the same time, the very high hydrogen bond propensity of water with itself prevents this solvent from providing an effective hydrogen bond-mediated solvation. Solvents bearing good HB acceptors, while lacking strong HB donors, contribute to further stabilizing solute-solvent adducts through weak and fluxional HBs that involve the amide groups of ADCA. Implications for the solubility of ADCA down to µM concentrations were evaluated, also with the aid of classical simulations of solution nanodroplets.
RESUMEN
A one-pot protocol, consisting of a Pd-catalysed carboamination reaction, followed by N-deprotection and oxidative denitrogenation, has been developed for the synthesis of diversely substituted spirocyclopropyloxindoles, in yields up to 73% and with diastereoselectivity close to 1 : 1. Readily accessible starting materials, mild reaction conditions, an easy to operate one-pot procedure and good functional group tolerance make this transformation a versatile tool for the synthesis of substituted spirocyclopropyloxindoles. This protocol successfully works on the gram-scale and allows access to both diastereoisomers separately. A plausible mechanism was proposed, and a series of post-transformations were performed on the obtained products, showing their remarkable synthetic versatility.
RESUMEN
Classical molecular dynamics simulations are carried out to investigate the aggregation of supercooled benzoic acid in confined spaces. Nanocavities, nanotubes and nanolayers are defined by restricting the periodicity of the simulation to zero, one or two dimensions, with boundaries set by adjustable, general, and computationally cheap van der Waals barriers. The effect of different confinement geometries is explored. It is found that the confinement impacts the liquid collective dynamics, strengthening the correlations that affect the motion of distant molecules. Overall, confinement determines up to a tenfold increase of the viscosity of the liquid and strongly slows down the rotational correlation times. Aggregation mediated by interactions with the walls and partial polarization of the liquid are observed. Additionally, transitions to high-density liquid states occur when stiffer barriers are used. In general, a reduced accessible amount of phase space fosters the struggle for a closer packing to relieve unfavorable atom-atom contacts, while maximizing the attractive ones. In benzoic acid, this implies that the hydrogen bond network is organized more efficiently in high density states.
RESUMEN
The noncovalent interaction (NCI) index is nowadays a well-known strategy to detect NCIs in molecular systems. Even though it initially provided only qualitative descriptions, the technique has been recently extended to also extract quantitative information. To accomplish this task, integrals of powers of the electron distribution were considered, with the requirement that the overall electron density can be clearly decomposed as sum of distinct fragment contributions to enable the definition of the (noncovalent) integration region. So far, this was done by only exploiting approximate promolecular electron densities, which are given by the sum of spherically averaged atomic electron distributions and thus represent too crude approximations. Therefore, to obtain more quantum mechanically (QM) rigorous results from NCI index analyses, in this work, we propose to use electron densities obtained through the transfer of extremely localized molecular orbitals (ELMOs) or through the recently developed QM/ELMO embedding technique. Although still approximate, the electron distributions resulting from the abovementioned methods are fully QM and, above all, are again partitionable into subunit contributions, which makes them completely suitable for the NCI integral approach. Therefore, we benchmarked the integrals resulting from NCI index analyses (both those based on the promolecular densities and those based on ELMO electron distributions) against interaction energies computed at a high quantum chemical level (in particular, at the coupled cluster level). The performed test calculations have indicated that the NCI integrals based on ELMO electron densities outperform the promolecular ones. Furthermore, it was observed that the novel quantitative NCI-(QM/)ELMO approach can be also profitably exploited both to characterize and evaluate the strength of specific interactions between ligand subunits and protein residues in protein-ligand complexes and to follow the evolution of NCIs along trajectories of molecular dynamics simulations. Although further methodological improvements are still possible, the new quantitative ELMO-based technique could be already exploited in situations in which fast and reliable assessments of NCIs are crucial, such as in computational high-throughput screenings for drug discovery.
RESUMEN
In a quite recent study [Genoni et al. (2017). IUCrJ, 4, 136-146], it was observed that the X-ray restrained wavefunction (XRW) approach allows a more efficient and larger capture of electron correlation effects on the electron density if high-angle reflections are not considered in the calculations. This is due to the occurrence of two concomitant effects when one uses theoretical X-ray diffraction data corresponding to a single-molecule electron density in a large unit cell: (i) the high-angle reflections are generally much more numerous than the low- and medium-angle ones, and (ii) they are already very well described at unrestrained level. Nevertheless, since high-angle data also contain important information that should not be disregarded, it is not advisable to neglect them completely. For this reason, based on the results of the previous investigation, this work introduces a weighting scheme for XRW calculations to up-weight the contribution of low- and medium-angle reflections, and, at the same time, to reasonably down-weight the importance of the high-angle data. The proposed strategy was tested through XRW computations with both theoretical and experimental structure-factor amplitudes. The tests have shown that the new weighting scheme works optimally if it is applied with theoretically generated X-ray diffraction data, while it is not advantageous when traditional experimental X-ray diffraction data (even of very high resolution) are employed. This also led to the conclusion that the use of a specific external parameter λJ for each resolution range might not be a suitable strategy to adopt in XRW calculations exploiting experimental X-ray data as restraints.
RESUMEN
Despite great advances in X-ray absorption spectroscopy for the investigation of small molecule electronic structure, the application to biosystems of experimental techniques developed within this research field remains a challenge. To partially circumvent the problem, users resort to theoretical methods to interpret or predict the X-ray absorption spectra of large molecules. To accomplish this task, only low-cost computational strategies can be exploited. For this reason, some of them are single Slater determinant wavefunction approaches coupled with multiscale embedding techniques designed to treat large systems of biological interest. Therefore, in this work, we propose to apply the recently developed IMOM/ELMO embedding method to the determination of core-ionized states. The IMOM/ELMO technique resulted from the combination of the single Slater determinant Δself-consistent-field-initial maximum overlap approach (ΔSCF-IMOM) with the QM/ELMO (quantum mechanics/extremely localized molecular orbital) embedding strategy, a method where only the chemically relevant region of the examined system is treated at fully quantum chemical level, while the rest is described through transferred and frozen extremely localized molecular orbitals (ELMOs). The IMOM/ELMO technique was initially validated by computing core-ionization energies for small molecules, and it was afterwards exploited to study larger biosystems. The obtained results are in line with those reported in previous studies that applied alternative ΔSCF approaches. This makes us envisage a possible future application of the proposed method to the interpretation of X-ray absorption spectra of large molecules.
Asunto(s)
Teoría Cuántica , Modelos Moleculares , Fenómenos Químicos , Espectroscopía de Absorción de Rayos XRESUMEN
Quantum chemistry offers a large variety of methods to treat excited states. Many of them are based on a multireference wave function ansatz and are therefore characterized by an intrinsic complexity and high computational costs. To overcome these drawbacks and also some limitations of simpler single-reference approaches (e.g., configuration interaction singles and time-dependent density functional theory), the single-determinant Δself-consistent field-initial maximum overlap method (ΔSCF-IMOM) has been proposed. This strategy substitutes the aufbau principle with a criterion that occupies molecular orbitals at successive SCF iterations on the basis of their maximum overlap with a proper set of guess orbitals for the target excited state. In this way, it prevents the SCF process to collapse to the ground state wave function and provides excited state single Slater determinant solutions to the SCF equations. Here, we propose to extend the applicability of the IMOM to the treatment of localized excited states of large systems. To accomplish this task, we coupled it with the QM/ELMO (quantum mechanics/extremely localized molecular orbitals) strategy, a quantum mechanical embedding method in which the most chemically relevant part of the system is treated with traditional quantum chemical approaches, while the rest is described by extremely localized molecular orbitals transferred from recently constructed libraries or proper model molecules. After presenting the theoretical foundations of the new IMOM/ELMO technique, in this paper, we will show and discuss the results of preliminary test calculations carried out on both model systems (i.e., decanoic acid, decene, decapentaene, and solvated acrolein) and a system of biological interest (flavin mononucleotide in the flavodoxin protein). We observed that, for localized excited states, the new IMOM/ELMO method provides reliable results, and it reproduces the outcomes of fully IMOM calculations within the chemical accuracy threshold (i.e., 0.043 eV) by including only a limited number of atoms in the QM region. Furthermore, the first application of our embedding technique to a larger biological system gave completely plausible results in line with those obtained through more traditional quantum mechanical methods, thus opening the possibility of using the new approach in future investigations of photobiology problems.
RESUMEN
Quantum mechanics/molecular mechanics (QM/MM) calculations are widely used embedding techniques to computationally investigate enzyme reactions. In most QM/MM computations, the quantum mechanical region is treated through density functional theory (DFT), which offers the best compromise between chemical accuracy and computational cost. Nevertheless, to obtain more accurate results, one should resort to wave function-based methods, which however lead to a much larger computational cost already for relatively small QM subsystems. To overcome this drawback, we propose the coupling of our QM/ELMO (quantum mechanics/extremely localized molecular orbital) approach with molecular mechanics, thus introducing the three-layer QM/ELMO/MM technique. The QM/ELMO strategy is an embedding method in which the chemically relevant part of the system is treated at the quantum mechanical level, while the rest is described through frozen ELMOs. Since the QM/ELMO method reproduces results of fully QM computations within chemical accuracy and with a much lower computational effort, it can be considered a suitable strategy to extend the range of applicability and accuracy of the QM/MM scheme. In this paper, other than briefly presenting the theoretical bases of the QM/ELMO/MM technique, we will also discuss its validation on the well-tested deprotonation of acetyl coenzyme A by aspartate in citrate synthase.
RESUMEN
The development of computationally advantageous methods for the study of large systems is a long-standing research topic in theoretical chemistry. Among these techniques, a prominent place is certainly occupied by the multiscale embedding strategies, from the well-known QM/MM (quantum mechanics/molecular mechanics) methods to the latest and promising fully quantum mechanical approaches. In this Feature Article, we will briefly review the recently proposed QM/ELMO (quantum mechanics/extremely localized molecular orbital) scheme, namely a new multiscale embedding strategy in which the most chemically relevant region of the investigated system is treated at fully quantum chemical level, while the remaining part (namely, the environment) is described by means of transferred extremely localized molecular orbitals that remain frozen throughout the computation. Other than highlighting the theoretical bases, here we will also review the main results obtained through all the currently available variants of the novel method. In particular, we will show how the QM/ELMO embedding scheme has been successfully exploited to perform both ground and excited state calculations, reproducing the results of corresponding fully quantum mechanical computations but with a much lower computational cost. A first application to crystallography will be also discussed, and we will describe how the QM/ELMO approach has been recently coupled with the Hirshfeld atom refinement technique to accurately determine the positions of hydrogen atoms from X-ray diffraction data. Given the reliability and quality of the obtained results, future applications of the current versions of the QM/ELMO embedding strategy to different types of chemical problems are to be expected in the near future. Moreover, further algorithmic improvements and methodological developments are also envisaged, such as the development of a polarizable QM/ELMO scheme accounting for the effects of the QM region on the ELMO subsystem or the use of the new embedding approach in the context of quantum crystallography to perform unprecedented accurate refinements of macromolecular crystal structures.
RESUMEN
The positions of hydrogen atoms in molecules are fundamental in many aspects of chemistry. Nevertheless, most molecular structures are obtained from refinements of X-ray data exploiting the independent atom model (IAM), which uses spherical atomic densities and provides bond lengths involving hydrogen atoms that are too short compared to the neutron reference values. To overcome the IAM shortcomings, the wave function-based Hirshfeld atom refinement (HAR) method has been recently proposed, emerging as a promising strategy able to give element-hydrogen bond distances in excellent agreement with the neutron ones in terms of accuracy and precision. In this Letter, we propose a significant improvement of HAR based on the idea of describing the crystal environment explicitly in the underlying wave function calculation through a quantum mechanical embedding strategy that exploits extremely localized molecular orbitals. Test-bed refinements on a crystal structure characterized by strong intermolecular interactions are also discussed.
RESUMEN
The QM/ELMO (quantum mechanics/extremely localized molecular orbital) method is a recently developed embedding technique in which the most important region of the system under examination is treated at fully quantum mechanical level, while the rest is described by means of transferred and frozen extremely localized molecular orbitals. In this paper, we propose the first application of the QM/ELMO approach to the investigation of excited states, and, in particular, we present the coupling of the QM/ELMO philosophy with Time-Dependent Density Functional Theory (TDDFT) and Equation-of-Motion Coupled Cluster with single and double substitutions (EOM-CCSD). The proposed TDDFT/ELMO and EOM-CCSD/ELMO strategies underwent a series of preliminary tests that were already considered for the validation of other embedding methods for excited states. The obtained results showed that the novel techniques allow the accurate description of localized excitations in large systems by only including a relatively small number of atoms in the region treated at fully quantum chemical level. Furthermore, for TDDFT/ELMO, it was also observed that (i) the method enables to avoid the presence of artificial low-lying charge-transfer states that may affect traditional TDDFT calculations, even using functionals that do not take into account long-range corrections, and (ii) the novel approach can be also successfully exploited to investigate local electronic transitions in quite large systems (e.g., reduced model of the Green Fluorescent Protein), and the accuracy of the results can be improved by including a sufficient number of chemically crucial fragments/residues in the quantum mechanical region. Finally, concerning EOM-CCSD/ELMO, it was also seen that, despite the quite crude approximation of an embedding potential given by frozen extremely localized molecular orbitals, the new strategy is able to satisfactorily account for the effects of the environment. This work paves the way to further extensions of the QM/ELMO philosophy for the study of local excitations in extended systems, suggesting the coupling of the QM/ELMO approach with other quantum chemical strategies for excited states, from the simplest ΔSCF techniques to the most advanced and computationally expensive multireferences methods.
RESUMEN
Despite its role in spin density functional theory and it being the basic observable for describing and understanding magnetic phenomena, few studies have appeared on the electron spin density subtleties thus far. A systematic full topological analysis of this function is lacking, seemingly in contrast to the blossoming in the last 20 years of many studies on the topological features of other scalar fields of chemical interest. We aim to fill this gap by unveiling the kind of information hidden in the spin density distribution that only its topology can disclose. The significance of the spin density critical points, the 18 different ways in which they can be realized and the peculiar topological constraints on their number and kind, arising from the presence of positive and negative spin density regions, is addressed. The notion of molecular spin graphs, spin maxima (minima) joining paths, spin basins and of their valence is introduced. We show that two kinds of structures are associated with a spin-polarized molecule: the usual one, defined through the electron density gradient, and the magnetic structure, defined through the spin density gradient and composed in general by at least two independent spin graphs, related to spin density maxima and minima. Several descriptors, such as the spin polarization index, are introduced to characterize the properties of spin density critical points and basins. The study on the general features of the spin density topology is followed by the specific example of the water molecule in the 3B1 triplet state, using spin density distributions of increasing accuracy.
Asunto(s)
Teoría Funcional de la Densidad , Electrones , Modelos Químicos , Algoritmos , Modelos Teóricos , Agua/químicaRESUMEN
Embedding strategies currently provide the best compromise between accuracy and computational cost in modeling chemical properties and processes of large and complex systems. In this framework, different methods have been proposed all over the years, from the very popular QM/MM approaches to the more recent and very promising density matrix and density functional embedding techniques. Here, we present a further development of the quantum mechanics/extremely localized molecular orbital technique (QM/ELMO) method, a recently proposed multiscale embedding strategy in which the chemically active region of the investigated system is treated at a fully quantum mechanical level, while the rest is described by frozen extremely localized molecular orbitals previously transferred from proper libraries or tailor-made model molecules. In particular, in this work we discuss and assess in detail the extension of the QM/ELMO approach to density functional theory and post-Hartree-Fock techniques by evaluating its performances when it is used to describe chemical reactions, bond dissociations, and intermolecular interactions. The preliminary test calculations have shown that, in the investigated cases, the new embedding strategy enables the results of the corresponding fully quantum mechanical computations to be reproduced within chemical accuracy in almost all the cases but with a significantly reduced computational cost, especially when correlated post-Hartree-Fock strategies are used to describe the quantum mechanical subsystem. In light of the obtained results, we already envisage the future application of the new correlated QM/ELMO techniques to the investigation of more challenging problems, such as the modeling of enzyme catalysis, the study of excited states of biomolecules, and the refinement of macromolecular X-ray crystal structures.
RESUMEN
Reported here is an entirely new application of experimental electron density (EED) in the study of magnetic anisotropy of single-molecule magnets (SMMs). Among those SMMs based on one single transition metal, tetrahedral CoII-complexes are prominent, and their large zero-field splitting arises exclusively from coupling between the d x 2 - y 2 and dxy orbitals. Using very low temperature single-crystal synchrotron X-ray diffraction data, an accurate electron density (ED) was obtained for a prototypical SMM, and the experimental d-orbital populations were used to quantify the dxy-d x 2 - y 2 coupling, which simultaneously provides the composition of the ground-state Kramers doublet wave function. Based on this experimentally determined wave function, an energy barrier for magnetic relaxation in the range 193-268â cm-1 was calculated, and is in full accordance with the previously published value of 230â cm-1 obtained from near-infrared spectroscopy. These results provide the first clear and direct link between ED and molecular magnetic properties.
RESUMEN
One of the well-established methods of modern quantum crystallography is undoubtedly the X-ray constrained wavefunction (XCW) approach, a technique that enables the determination of wavefunctions which not only minimize the energy of the system under examination, but also reproduce experimental X-ray diffraction data within the limit of the experimental errors. Initially proposed in the framework of the Hartree-Fock method, the strategy has been gradually extended to other techniques of quantum chemistry, but always remaining limited to a single-determinant ansatz for the wavefunction to extract. This limitation has been recently overcome through the development of the novel X-ray constrained spin-coupled (XCSC) approach [Genoni et al. (2018). Chem. Eur. J. 24, 15507-15511] which merges the XCW philosophy with the traditional spin-coupled strategy of valence bond theory. The main advantage of this new technique is the possibility of extracting traditional chemical descriptors (e.g. resonance structure weights) compatible with the experimental diffraction measurements, without the need to introduce information a priori or perform analyses a posteriori. This paper provides a detailed theoretical derivation of the fundamental equations at the basis of the XCSC method and also introduces a further advancement of its original version, mainly consisting in the use of molecular orbitals resulting from XCW calculations at the Hartree-Fock level to describe the inactive electrons in the XCSC computations. Furthermore, extensive test calculations, which have been performed by exploiting high-resolution X-ray diffraction data for salicylic acid and by adopting different basis sets, are presented and discussed. The computational tests have shown that the new technique does not suffer from particular convergence problems. Moreover, all the XCSC calculations provided resonance structure weights, spin-coupled orbitals and global electron densities slightly different from those resulting from the corresponding unconstrained computations. These discrepancies can be ascribed to the capability of the novel strategy to capture the information intrinsically contained in the experimental data used as external constraints.
RESUMEN
The development of methods for the quantum mechanical study of macromolecules has always been an important challenge in theoretical chemistry. Nowadays, the techniques proposed in this context can be used to investigate very large systems and can be subdivided into two main categories: fragmentation and embedding strategies. In this paper, by modifying and improving the local self-consistent field approach originally proposed for quantum mechanics/molecular mechanics techniques, we introduce the new multiscale embedding quantum mechanics/extremely localized molecular orbital (QM/ELMO) method. The new strategy enables treatment of chemically relevant regions of large biological molecules through usual methods of quantum chemistry while describing the remaining parts of the systems by means of frozen extremely localized molecular orbitals transferred from properly constructed libraries. Test calculations have shown the correct functioning and the high reliability of the new approach, thus anticipating its possible applications to different fields of physical chemistry, such as rational drug design and structural refinements of proteins.
RESUMEN
The electron density distribution of the complex monoanion Cu(CF3)4 in (Bu4N)[Cu(CF3)4] has been studied by high-resolution X-ray single-crystal diffraction and augmented with theoretical calculations. The study finds that the central copper bears an atomic charge of close to +1, while the occupancy of its d x2- y2 orbital is only 1.26. Using topological analysis combined with theoretical calculations, the depopulation of d x2- y2 is shown to be due to significant covalency in the Cu-C bonds. The combination of the monovalent picture and the covalency is interpreted as a confirmation of an inverted ligand field.
RESUMEN
The source function (SF) decomposes the electron density at any point into contributions from all other points in the molecule, complex, or crystal. The SF "illuminates" those regions in a molecule that most contribute to the electron density at a point of reference. When this point of reference is the bond critical point (BCP), a commonly used surrogate of chemical bonding, then the SF analysis at an atomic resolution within the framework of Bader's Quantum Theory of Atoms in Molecules returns the contribution of each atom in the system to the electron density at that BCP. The SF is used to locate the important regions that control the hydrogen bonds in both Watson-Crick (WC) DNA dimers (adenine:thymine (AT) and guanine:cytosine (GC)) which are studied in their neutral and their singly ionized (radical cationic and anionic) ground states. The atomic contributions to the electron density at the BCPs of the hydrogen bonds in the two dimers are found to be delocalized to various extents. Surprisingly, gaining or loosing an electron has similar net effects on some hydrogen bonds concealing subtle compensations traced to atomic sources contributions. Coarser levels of resolutions (groups, rings, and/or monomers-in-dimers) reveal that distant groups and rings often have non-negligible effects especially on the weaker hydrogen bonds such as the third weak CHâ â â O hydrogen bond in AT. Interestingly, neither the purine nor the pyrimidine in the neutral or ionized forms dominate any given hydrogen bond despite that the former has more atoms that can act as source or sink for the density at its BCP. © 2018 Wiley Periodicals, Inc.
Asunto(s)
Emparejamiento Base , ADN/química , Electrones , Teoría Cuántica , Enlace de HidrógenoRESUMEN
Magnetic properties of open-shell systems depend on their unpaired electron density distribution. Accurate spin density (SD) is difficult to retrieve, both from polarized neutron diffraction (PND) data and from quantum approaches, and its interpretation is not trivial. The Source Function is a useful tool to interpret SD distributions and their accuracy. It is here applied to analyze and compare the theoretical SD in a weakly ferromagnetically coupled end-end azido dicopper complex with that in a strongly-coupled end-on complex. The Source Function enables to highlight the origin of the SD differences between the two dicopper complexes and among adopted computational approaches (CASSCF, DFT, UHF). Further insight is provided by partial Source Function SD reconstructions using given subsets of atoms. DFT methods exaggerate electron sharing between copper and the ligands, causing spin delocalization toward them and overestimating metal-ligand spin polarization, while underestimating CASSCF spin information transmission between atoms. CAS(10,10) SD is closer to the PND SD than other adopted methods © 2018 Wiley Periodicals, Inc.
RESUMEN
Two conformational crystal polymorphs of 3-diethylamino-4-(4-methoxyphenyl)-1,1-dioxo-4H-1λ6,2-thiazete-4-carbonitrile (DTC) have been analyzed in the 100 K-room temperature range by single crystal X-ray diffraction and high quality DFT calculations. DTC has strongly polar as well as aliphatic substituents but no hydrogen bonding groups, and thus qualifies as a test molecule for the relative importance of electrostatic vs. dispersion-repulsion terms. The two polymorphs have the same P21/n space group and differ by a flipping of the -OCH3 group, the two conformations being almost equi-energetic and separated by a low barrier. The system is monotropic in the observed temperature range with nearly identical thermal expansion coefficients and energy-temperature slopes, one phase consistently predicted to be more stable in agreement with the relative ease of appearance. Energy decompositions show that the electrostatic term is dominant and stabilizes with decreasing temperature. Dispersion and repulsion show the expected behavior, the former becoming more stabilizing at lower temperature in contrast with increasing repulsion at higher density. Absolute values and trends are very similar in the two phases, explaining the small total energy difference. Geometrical analyses of intermolecular contacts using fingerprint plots, as well as the study of molecular dipole moments as a function of T in the framework of the Quantum Theory of Atoms in Molecules, reveal more details of phase stability.