Fast and Accurate Quantum Crystallography: From Small to Large, from Light to Heavy

TAAM refinement on high-resolution experimental and simulated 3D ED/MicroED data for organic molecules

3D electron diffraction (3D ED), or microcrystal electron diffraction (MicroED), has become an alternative technique for determining the high-resolution crystal structures of compounds from sub-micron-sized crystals. Here, we considered L-alanine, α-glycine and urea, which are known to form good-quality crystals, and collected high-resolution 3D ED data on our in-house TEM instrument. In this study, we present a comparison of independent atom model (IAM) and transferable aspherical atom model (TAAM) kinematical refinement against experimental and simulated data. TAAM refinement on both experimental and simulated data clearly improves the model fitting statistics (R factors and residual electrostatic potential) compared to IAM refinement. This shows that TAAM better represents the experimental electrostatic potential of organic crystals than IAM. Furthermore, we compared the geometrical parameters and atomic displacement parameters (ADPs) resulting from the experimental refinements with the simulated refinements, with the periodic density functional theory (DFT) calculations and with published X-ray and neutron crystal structures. The TAAM refinements on the 3D ED data did not improve the accuracy of the bond lengths between the non-H atoms. The experimental 3D ED data provided more accurate H-atom positions than the IAM refinements on the X-ray diffraction data. The IAM refinements against 3D ED data had a tendency to lead to slightly longer X-H bond lengths than TAAM, but the difference was statistically insignificant. Atomic displacement parameters were too large by tens of percent for L-alanine and α-glycine. Most probably, other unmodelled effects were causing this behaviour, such as radiation damage or dynamical scattering.

Transferable Hirshfeld atom model for rapid evaluation of aspherical atomic form factors

Form factors based on aspherical models of atomic electron density have brought great improvement in the accuracies of hydrogen atom parameters derived from X-ray crystal structure refinement. Today, two main groups of such models are available, the banks of transferable atomic densities parametrized using the Hansen-Coppens multipole model which allows for rapid evaluation of atomic form factors and Hirshfeld atom refinement (HAR)-related methods which are usually more accurate but also slower. In this work, a model that combines the ideas utilized in the two approaches is tested. It uses atomic electron densities based on Hirshfeld partitions of electron densities, which are precalculated and stored in a databank. This model was also applied during the refinement of the structures of five small molecules. A comparison of the resulting hydrogen atom parameters with those derived from neutron diffraction data indicates that they are more accurate than those obtained with the Hansen-Coppens based databank, and only slightly less accurate than those obtained with a version of HAR that neglects the crystal environment. The advantage of using HAR becomes more noticeable when the effects of the environment are included. To speed up calculations, atomic densities were represented by multipole expansion with spherical harmonics up to l = 7, which used numerical radial functions (a different approach to that applied in the Hansen-Coppens model). Calculations of atomic form factors for the small protein crambin (at 0.73 Å resolution) took only 68 s using 12 CPU cores.

Refinement of X-ray and electron diffraction crystal structures using analytical Fourier transforms of Slater-type atomic wavefunctions in Olex2

An implementation of Slater-type spherical scattering factors for X-ray and electron diffraction for elements in the range Z = 1-103 is presented within the software Olex2. Both high- and low-angle Fourier behaviour of atomic electron density and electrostatic potential can thus be addressed, in contrast to the limited flexibility of the four Gaussian plus constant descriptions which are currently the most widely used method for calculating atomic scattering factors during refinement. The implementation presented here accommodates the increasing complexity of the electronic structure of heavier elements by using complete atomic wavefunctions without any interpolation between precalculated tables or intermediate fitting functions. Atomic wavefunctions for singly charged ions are implemented and made accessible, and these show drastic changes in electron diffraction scattering factors compared with the neutral atom. A comparison between the two different spherical models of neutral atoms is presented as an example for four different kinds of X-ray and two electron diffraction structures, and comparisons of refinement results using the existing diffraction data are discussed. A systematic but slight improvement in R values and residual densities can be observed when using the new scattering factors, and this is discussed relative to effects on the atomic displacement parameters and atomic positions, which are prominent near the heavier elements in a structure.

Exploiting the full potential of cryo-EM maps

The Coulomb potential maps generated by electron microscopy (EM) experiments contain not only information about the position but also about the charge state of the atom. This feature of EM maps allows the identification of specific ions and the protonation state of amino acid side chains in the sample. Here, we summarize qualitative observations of charges in EM maps, discuss the difficulties in interpreting the charge in Coulomb potential maps with respect to distinguishing it from radiation damage, and outline considerations to implement the correct charge in fitting algorithms.

Enhancing hydrogen positions in X-ray structures of transition metal hydride complexes with dynamic quantum crystallography

Hirshfeld atom refinement (HAR) is a method which enables the user to obtain more accurate positions of hydrogen atoms bonded to light chemical elements using X-ray data. When data quality permits, this method can be extended to hydrogen-bonded transition metals (TMs), as in hydride complexes. However, addressing hydrogen thermal motions with HAR, particularly in TM hydrides, presents a challenge. At the same time, proper description of thermal vibrations can be vital for determining hydrogen positions correctly. In this study, we employ tools such as SHADE3 and Normal Mode Refinement (NoMoRe) to estimate anisotropic displacement parameters (ADPs) for hydrogen atoms during HAR and IAM refinements performed for seven structures of TM (Fe, Ni, Cr, Nb, Rh and Os) and metalloid (Sb) hydride complexes for which both the neutron and the X-ray structures have been determined. A direct comparison between neutron and HAR/SHADE3/NoMoRe ADPs reveals that the similarity between neutron hydrogen ADPs and those estimated with NoMoRe or SHADE3 is significantly higher than when hydrogen ADPs are refined with HAR. Regarding TM-H bond lengths, traditional HAR exhibits a slight advantage over the other methods. However, combining NoMoRe/SHADE3 with HAR results in a minor decrease in agreement with neutron TM-H bond lengths. For the Cr complex, for which high-resolution X-ray data were collected, an investigation of resolution-related effects was possible.

Hirshfeld Atom Refinement of Metal-Organic Complexes: Treatment of Hydrogen Atoms Bonded to Transition Metals

Hydrogen positions in hydrides play a key role in hydrogen storage materials and high-temperature superconductors. Our recently published study of five crystal structures of transition-metal-bound hydride complexes showed that using aspherical atomic scattering factors for Hirshfeld atom refinement (HAR) resulted in a systematic elongation of metal-hydrogen bonds compared to using spherical scattering factors with the Independent Atom Model (IAM). Even though only standard-resolution X-ray data was used, for the highest-quality data, we obtained excellent agreement between the X-ray and the neutron-derived bond lengths. We present an extended version of this study including 10 crystal structures of metal-organic complexes containing hydrogen atoms bonded to transition-metal atoms for which both X-ray and neutron data are available. The neutron structures were used as a benchmark, and the X-ray structures were refined by applying Hirshfeld atom refinement using various basis sets and DFT functionals in order to investigate the influence of the technical aspects on the length of metal-hydrogen bonds. The result of including relativistic effects in the Hamiltonian and using a cluster of multipoles simulating interactions with a crystal environment during wave function calculations was examined. The effect of the data quality on the final result was also evaluated. The study confirms that a high quality of experimental data is the key factor allowing us to obtain significant improvement in transition metal (TM)-hydrogen bond lengths from HAR in comparison with the IAM. Individual adjustments and better choices of the basis set can improve hydrogen positions. Average differences between TM-H bond lengths obtained with various DFT functionals upon including relativistic effects or between double-ζ and triple-ζ basis sets were not statistically significant. However, if all bonds formed by H atoms were considered, significant differences caused by different refinement strategies were observed. Finally, we examined the refinement of atomic thermal motions. Anisotropic refinement of hydrogen thermal motions with HAR was feasible only in some cases, and isotropically refined hydrogen thermal motions were in similar agreement with neutron values whether obtained with HAR or with the IAM.

Extracting Quantitative Information at Quantum Mechanical Level from Noncovalent Interaction Index Analyses

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.

Aspherical atom refinements on X-ray data of diverse structures including disordered and covalent organic framework systems: a time-accuracy trade-off

Aspherical atom refinement is the key to achieving accurate structure models, displacement parameters, hydrogen-bond lengths and analysis of weak interactions, amongst other examples. There are various quantum crystallographic methods to perform aspherical atom refinement, including Hirshfeld atom refinement (HAR) and transferable aspherical atom model (TAAM) refinement. Both HAR and TAAM have their limitations and advantages, the former being more accurate and the latter being faster. With the advent of non-spherical atoms in Olex2 (NoSpherA2), it is now possible to overcome some limitations, like treating disorder, twinning and network structures, in aspherical refinements using HAR, TAAM or both together. TAAM refinement in NoSpherA2 showed significant improvement in refinement statistics compared with independent atom model (IAM) refinements on a diverse set of X-ray diffraction data. The sensitivity of TAAM towards poor data quality and disorder was observed in terms of higher refinement statistics for such structures. A comparison of IAM with TAAM and HAR in NoSpherA2 indicated that the time taken by TAAM refinements was of the same order of magnitude as that taken by IAM, while in HAR the time taken using a minimal basis set was 50 times higher than for IAM and rapidly increased with increasing size of the basis sets used. The displacement parameters for hydrogen and non-hydrogen atoms were very similar in both HAR and TAAM refinements. The hydrogen-bond lengths were slightly closer to neutron reference values in the case of HAR with higher basis sets than in TAAM. To benefit from the advantages of each method, a new hybrid refinement approach has been introduced, allowing a combination of IAM, HAR and TAAM in one structure refinement. Refinement of coordination complexes involving metal-organic compounds and network structures such as covalent organic frameworks and metal-organic frameworks is now possible in a hybrid mode such as IAM-TAAM or HAR-TAAM, where the metal atoms are treated via either the IAM or HAR method and the organic part via TAAM, thus reducing the computational costs without compromising the accuracy. Formal charges on the metal and ligand can also be introduced in hybrid-mode refinement.

Accurate crystal structure of ice VI from X-ray diffraction with Hirshfeld atom refinement

Water is an essential chemical compound for living organisms, and twenty of its different crystal solid forms (ices) are known. Still, there are many fundamental problems with these structures such as establishing the correct positions and thermal motions of hydrogen atoms. The list of ice structures is not yet complete as DFT calculations have suggested the existence of additional and - to date - unknown phases. In many ice structures, neither neutron diffraction nor DFT calculations nor X-ray diffraction methods can easily solve the problem of hydrogen atom disorder or accurately determine their anisotropic displacement parameters (ADPs). Here, accurate crystal structures of H2O, D2O and mixed (50%H2O/50%D2O) ice VI obtained by Hirshfeld atom refinement (HAR) of high-pressure single-crystal synchrotron and laboratory X-ray diffraction data are presented. It was possible to obtain O-H/D bond lengths and ADPs for disordered hydrogen atoms which are in good agreement with the corresponding single-crystal neutron diffraction data. These results show that HAR combined with X-ray diffraction can compete with neutron diffraction in detailed studies of polymorphic forms of ice and crystals of other hydrogen-rich compounds. As neutron diffraction is relatively expensive, requires larger crystals which can be difficult to obtain and access to neutron facilities is restricted, cheaper and more accessible X-ray measurements combined with HAR can facilitate the verification of the existing ice polymorphs and the quest for new ones.

Influence of modelling disorder on Hirshfeld atom refinement results of an organo-gold(I) compound

Details of the validation of disorder modelling with Hirshfeld atom refinement (HAR) for a previously investigated organo-gold(I) compound are presented here. The impact of refining disorder on HAR results is discussed using an analysis of the differences of dynamic structure factors. These dynamic structure factor differences are calculated from thermally smeared quantum mechanical electron densities based on wavefunctions that include or exclude electron correlation and relativistic effects. When disorder is modelled, the electron densities stem from a weighted superposition of two (or more) different conformers. Here this is shown to impact the relative importance of electron correlation and relativistic effect estimates expressed by the structure factor magnitudes. The role of disorder modelling is also compared with the effect of the treatment of hydrogen anisotropic displacement parameter (ADP) values and atomic anharmonicity of the gold atom. The analysis of ADP values of gold and disordered carbon atoms showed that the effect of disorder significantly altered carbon ADP values and did not influence those of the gold atom.

X-ray constrained wavefunctions based on Hirshfeld atoms. I. Method and review

The X-ray constrained wavefunction (XCW) procedure for obtaining an experimentally reconstructed wavefunction from X-ray diffraction data is reviewed. The two-center probability distribution model used to perform nuclear-position averaging in the original paper [Grimwood & Jayatilaka (2001). Acta Cryst. A57, 87-100] is carefully distinguished from the newer one-center probability distribution model. In the one-center model, Hirshfeld atoms are used, and the Hirshfeld atom based X-ray constrained wavefunction (HA-XCW) procedure is described for the first time, as well as its efficient implementation. In this context, the definition of the related X-ray wavefunction refinement (XWR) method is refined. The key halting problem for the XCW method - the procedure by which one determines when overfitting has occurred - is named and work on it reviewed.

Fragmentation and transferability in Hirshfeld atom refinement

Hirshfeld atom refinement (HAR) is one of the most effective methods for obtaining accurate structural parameters for hydrogen atoms from X-ray diffraction data. Unfortunately, it is also relatively computationally expensive, especially for larger molecules due to wavefunction calculations. Here, a fragmentation approach has been tested as a remedy for this problem. It gives an order of magnitude improvement in computation time for larger organic systems and is a few times faster for metal-organic systems at the cost of only minor differences in the calculated structural parameters when compared with the original HAR calculations. Fragmentation was also applied to polymeric and disordered systems where it provides a natural solution to problems that arise when HAR is applied. The concept of fragmentation is closely related to the transferable aspherical atom model (TAAM) and allows insight into possible ways to improve TAAM. Hybrid approaches combining fragmentation with the transfer of atomic densities between chemically similar atoms have been tested. An efficient handling of intermolecular interactions was also introduced for calculations involving fragmentation. When applied in fragHAR (a fragmentation approach for polypeptides) as a replacement for the original approach, it allowed for more efficient calculations. All of the calculations were performed with a locally modified version of Olex2 combined with a development version of discamb2tsc and ORCA. Care was taken to efficiently use the power of multicore processors by simple implementation of load-balancing, which was found to be very important for lowering computational time.

Hirshfeld atom refinement based on projector augmented wave densities with periodic boundary conditions

Hirshfeld atom refinement (HAR) is an X-ray diffraction refinement method that, in numerous publications, has been shown to give H-atom bond lengths in close agreement with neutron diffraction derived values. Presented here is a first evaluation of an approach using densities derived from projector augmented wave (PAW) densities with three-dimensional periodic boundary conditions for HAR. The results show an improvement over refinements that neglect the crystal environment or treat it classically, while being on a par with non-periodic approximations for treating the solid-state environment quantum mechanically. A suite of functionals were evaluated for this purpose, showing that the SCAN and revSCAN functionals are most suited to these types of calculation.

Vanishing of the atomic form factor derivatives in non-spherical structural refinement - a key approximation scrutinized in the case of Hirshfeld atom refinement

When calculating derivatives of structure factors, there is one particular term (the derivatives of the atomic form factors) that will always be zero in the case of tabulated spherical atomic form factors. What happens if the form factors are non-spherical? The assumption that this particular term is very close to zero is generally made in non-spherical refinements (for example, implementations of Hirshfeld atom refinement or transferable aspherical atom models), unless the form factors are refinable parameters (for example multipole modelling). To evaluate this general approximation for one specific method, a numerical differentiation was implemented within the NoSpherA2 framework to calculate the derivatives of the structure factors in a Hirshfeld atom refinement directly as accurately as possible, thus bypassing the approximation altogether. Comparing wR₂ factors and atomic parameters, along with their uncertainties from the approximate and numerically differentiating refinements, it turns out that the impact of this approximation on the final crystallographic model is indeed negligible.

Initial Maximum Overlap Method for Large Systems by the Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Technique

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.

Relativistic Hirshfeld atom refinement of an organo-gold(I) compound

The main goal of this study is the validation of relativistic Hirshfeld atom refinement (HAR) as implemented in Tonto for high-resolution X-ray diffraction datasets of an organo-gold(I) compound. The influence of the relativistic effects on statistical parameters, geometries and electron density properties was analyzed and compared with the influence of electron correlation and anharmonic atomic motions. Recent work in this field has indicated the importance of relativistic effects in the static electron density distribution of organo-mercury compounds. This study confirms that differences in electron density due to relativistic effects are also of significant magnitude for organo-gold compounds. Relativistic effects dominate not only the core region of the gold atom, but also influence the electron density in the valence and bonding region, which has measurable consequences for the HAR refinement model parameters. To study the effects of anharmonic motion on the electron density distribution, dynamic electron density difference maps were constructed. Unlike relativistic and electron correlation effects, the effects of anharmonic nuclear motion are mostly observed in the core area of the gold atom.

Three-Layer Multiscale Approach Based on Extremely Localized Molecular Orbitals to Investigate Enzyme Reactions

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.

lamaGOET: an interface for quantum crystallography

In quantum crystallography, theoretical calculations and crystallographic refinements are closely intertwined. This means that the employed software must be able to perform both quantum-mechanical calculations and crystallographic least-squares refinements. So far, the program Tonto is the only one able to do that. The lamaGOET interface described herein deals with this issue since it interfaces dedicated quantum-chemical software (the widely used Gaussian package and the specialized ELMOdb program) with the refinement capabilities of Tonto. Three different flavours of quantum-crystallographic refinements of the dipetide glycyl-l-threonine dihydrate are presented to showcase the capabilities of lamaGOET: Hirshfeld atom refinement (HAR), HAR-ELMO, namely HAR coupled with extremely localized molecular orbitals, and X-ray constrained wavefunction fitting.

The advanced treatment of hydrogen bonding in quantum crystallography

Although hydrogen bonding is one of the most important motifs in chemistry and biology, H-atom parameters are especially problematic to refine against X-ray diffraction data. New developments in quantum crystallography offer a remedy. This article reports how hydrogen bonds are treated in three different quantum-crystallographic methods: Hirshfeld atom refinement (HAR), HAR coupled to extremely localized molecular orbitals and X-ray wavefunction refinement. Three different compound classes that form strong intra- or intermolecular hydrogen bonds are used as test cases: hydrogen maleates, the tripeptide l-alanyl-glycyl-l-alanine co-crystallized with water, and xylitol. The differences in the quantum-mechanical electron densities underlying all the used methods are analysed, as well as how these differences impact on the refinement results.

Towards accurate and precise positions of hydrogen atoms bonded to heavy metal atoms

A comparison of five X-ray structures of transition-metal-bound hydride complexes, successfully refined using Hirshfeld Atom Refinement (HAR) against low resolution X-ray diffraction data (including the positions and ADPs of all hydrogen atoms), with neutron structures shows that using aspherical atomic scattering factors instead of spherical ones results in systematic elongation of metal-hydrogen bonds, which in the case of the highest-quality data leads to excellent agreement of the X-ray and the neutron-derived bond lengths.

Climbing Jacob's Ladder of Structural Refinement: Introduction of a Localized Molecular Orbital-Based Embedding for Accurate X-ray Determinations of Hydrogen Atom Positions

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.

A Step toward the Quantification of Noncovalent Interactions in Large Biological Systems: The Independent Gradient Model-Extremely Localized Molecular Orbital Approach

The independent gradient model (IGM) is a recent electron density-based computational method that enables to detect and quantify covalent and noncovalent interactions. When applied to large systems, the original version of the technique still relies on promolecular electron densities given by the sum of spherically averaged atomic electron distributions, which leads to approximate evaluations of the inter- and intramolecular interactions occurring in systems of biological interest. To overcome this drawback and perform IGM analyses based on quantum mechanically rigorous electron densities also for macromolecular systems, we coupled the IGM approach with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) that allow fast and reliable reconstructions of polypeptide and protein electron densities. The validation tests performed on small polypeptides and peptide dimers have shown that the novel IGM-ELMO strategy provides results that are systematically closer to the fully quantum mechanical ones and outperforms the IGM method based on the crude promolecular approximation, but still keeping a quite low computational cost. The results of the test calculations carried out on proteins have also confirmed the trends observed for the IGM analyses conducted on small systems. This makes us envisage the future application of the novel IGM-ELMO approach to unravel complicated noncovalent interaction networks (e.g., in protein-protein contacts) or to rationally design new drugs through molecular docking calculations and virtual high-throughput screenings.

HgH2 meets relativistic quantum crystallography. How to teach relativity to a non-relativistic wavefunction

The capability of X-ray constrained wavefunction (XCW) fitting to introduce relativistic effects into a non-relativistic wavefunction is tested. It is quantified how much of the reference relativistic effects can be absorbed in the non-relativistic XCW calculation when fitted against relativistic structure factors of a model HgH₂ molecule. Scaling of the structure-factor sets to improve the agreement statistics is found to introduce a significant systematic error into the XCW fitting of relativistic effects.

Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Strategy for Excited States: Coupling to Time-Dependent Density Functional Theory and Equation-of-Motion Coupled Cluster

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.

Accurate crystal structures and chemical properties from NoSpherA2

The relationship between the structure and the properties of a drug or material is a key concept of chemistry. Knowledge of the three-dimensional structure is considered to be of such importance that almost every report of a new chemical compound is accompanied by an X-ray crystal structure - at least since the 1970s when diffraction equipment became widely available. Crystallographic software of that time was restricted to very limited computing power, and therefore drastic simplifications had to be made. It is these simplifications that make the determination of the correct structure, especially when it comes to hydrogen atoms, virtually impossible. We have devised a robust and fast system where modern chemical structure models replace the old assumptions, leading to correct structures from the model refinement against standard in-house diffraction data using no more than widely available software and desktop computing power. We call this system NoSpherA2 (Non-Spherical Atoms in Olex2). We explain the theoretical background of this technique and demonstrate the far-reaching effects that the improved structure quality that is now routinely available can have on the interpretation of chemical problems exemplified by five selected examples.

Hirshfeld atom like refinement with alternative electron density partitions

Hirshfeld atom refinement is one of the most successful methods for the accurate determination of structural parameters for hydrogen atoms from X-ray diffraction data. This work introduces a generalization of the method [generalized atom refinement (GAR)], consisting of the application of various methods of partitioning electron density into atomic contributions. These were tested on three organic structures using the following partitions: Hirshfeld, iterative Hirshfeld, iterative stockholder, minimal basis iterative stockholder and Becke. The effects of partition choice were also compared with those caused by other factors such as quantum chemical methodology, basis set, representation of the crystal field and a combination of these factors. The differences between the partitions were small in terms of R factor (e.g. much smaller than for refinements with different quantum chemistry methods, i.e. Hartree-Fock and coupled cluster) and therefore no single partition was clearly the best in terms of experimental data reconstruction. In the case of structural parameters the differences between the partitions are comparable to those related to the choice of other factors. We have observed the systematic effects of the partition choice on bond lengths and ADP values of polar hydrogen atoms. The bond lengths were also systematically influenced by the choice of electron density calculation methodology. This suggests that GAR-derived structural parameters could be systematically improved by selecting an optimal combination of the partition and quantum chemistry method. The results of the refinements were compared with those of neutron diffraction experiments. This allowed a selection of the most promising partition methods for further optimization of GAR settings, namely the Hirshfeld, iterative stockholder and minimal basis iterative stockholder.

Quantum Crystallography in the Last Decade: Developments and Outlooks

In this review article, we report on the recent progresses in the field of quantum crystallography that has witnessed a massive increase of production coupled with a broadening of the scope in the last decade. It is shown that the early thoughts about extracting quantum mechanical information from crystallographic experiments are becoming reality, although a century after prediction. While in the past the focus was mainly on electron density and related quantities, the attention is now shifting toward determination of wavefunction from experiments, which enables an exhaustive determination of the quantum mechanical functions and properties of a system. Nonetheless, methods based on electron density modelling have evolved and are nowadays able to reconstruct tiny polarizations of core electrons, coupling charge and spin models, or determining the quantum behaviour at extreme conditions. Far from being routine, these experimental and computational results should be regarded with special attention by scientists for the wealth of information on a system that they actually contain.

TAAM: a reliable and user friendly tool for hydrogen-atom location using routine X-ray diffraction data

Hydrogen is present in almost all of the molecules in living things. It is very reactive and forms bonds with most of the elements, terminating their valences and enhancing their chemistry. X-ray diffraction is the most common method for structure determination. It depends on scattering of X-rays from electron density, which means the single electron of hydrogen is difficult to detect. Generally, neutron diffraction data are used to determine the accurate position of hydrogen atoms. However, the requirement for good quality single crystals, costly maintenance and the limited number of neutron diffraction facilities means that these kind of results are rarely available. Here it is shown that the use of Transferable Aspherical Atom Model (TAAM) instead of Independent Atom Model (IAM) in routine structure refinement with X-ray data is another possible solution which largely improves the precision and accuracy of X-H bond lengths and makes them comparable to averaged neutron bond lengths. TAAM, built from a pseudoatom databank, was used to determine the X-H bond lengths on 75 data sets for organic molecule crystals. TAAM parametrizations available in the modified University of Buffalo Databank (UBDB) of pseudoatoms applied through the DiSCaMB software library were used. The averaged bond lengths determined by TAAM refinements with X-ray diffraction data of atomic resolution (d_min ≤ 0.83 Å) showed very good agreement with neutron data, mostly within one single sample standard deviation, much like Hirshfeld atom refinement (HAR). Atomic displacements for both hydrogen and non-hydrogen atoms obtained from the refinements systematically differed from IAM results. Overall TAAM gave better fits to experimental data of standard resolution compared to IAM. The research was accompanied with development of software aimed at providing user-friendly tools to use aspherical atom models in refinement of organic molecules at speeds comparable to routine refinements based on spherical atom model.

Localized Molecular Orbital-Based Embedding Scheme for Correlated Methods

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.

fragHAR: towards ab initio quantum-crystallographic X-ray structure refinement for polypeptides and proteins

The first ab initio aspherical structure refinement against experimental X-ray structure factors for polypeptides and proteins using a fragmentation approach to break up the protein into residues and solvent, thereby speeding up quantum-crystallographic Hirshfeld atom refinement (HAR) calculations, is described. It it found that the geometric and atomic displacement parameters from the new fragHAR method are essentially unchanged from a HAR on the complete unfragmented system when tested on dipeptides, tripeptides and hexapeptides. The largest changes are for the parameters describing H atoms involved in hydrogen-bond interactions, but it is shown that these discrepancies can be removed by including the interacting fragments as a single larger fragment in the fragmentation scheme. Significant speed-ups are observed for the larger systems. Using this approach, it is possible to perform a highly parallelized HAR in reasonable times for large systems. The method has been implemented in the TONTO software.

On the use of the Obara–Saika recurrence relations for the calculation of structure factors in quantum crystallography

Modern methods of quantum crystallography are techniques firmly rooted in quantum chemistry and, as in many quantum chemical strategies, electron densities are expressed as two-centre expansions that involve basis functions centred on atomic nuclei. Therefore, the computation of the necessary structure factors requires the evaluation of Fourier transform integrals of basis function products. Since these functions are usually Cartesian Gaussians, in this communication it is shown that the Fourier integrals can be efficiently calculated by exploiting an extension of the Obara–Saika recurrence formulas, which are successfully used by quantum chemists in the computation of molecular integrals. Implementation and future perspectives of the technique are also discussed.