Noncovalent Interactions in Density Functional Theory

Evaluating Noncovalent Interactions in Halogenated Molecules with Double-Hybrid Functionals and a Dedicated Small Basis Set

We present here an extension of our recently developed PBE-QIDH/DH-SVPD basis set to halogen atoms, with the aim of obtaining, for weakly interacting halogenated molecules, interaction energies close to those provided by a large basis set (def2-TZVPP) coupled to empirical dispersion potential. The core of our approach is the split-valence basis set, DH-SVPD, that has been developed for F, Cl, Br, and I atoms using a self-consistent formula, containing only energy terms computed for dimers and the corresponding monomers at the same level of theory. The basis set developed considering four systems, one for each halogen atoms, has been then tested on the X40, X4 × 10 benchmarks as well as on other two, less standard, data sets. Finally, a large system (380 atoms) has been also considered as a "crash" test. Our results show that the simple and nonempirical PBE-QIDH/DH-SVPD approach is able to provide accurate results for interaction energies of all the considered systems and can thus be considered as a cheaper alternative to DH functionals paired with empirical dispersion corrections and a large basis set of triple-ζ quality.

Exploring the Reactivity of Donor-Acceptor Systems through a Combined Conceptual and Constrained DFT Approach

In the context of the conceptual density functional theory (cDFT) and based on the computational efficiency of the constrained DFT (CDFT), we demonstrate that chemical reactivity can be governed by the difference between the local interacting chemical potentials of the reactants (referred as Edual), in agreement with Sanderson's equalization principle. In a proof-of-concept study, we investigated illustrative examples involving typical non-covalent donor-acceptor systems and reactive systems are provided. For the selected systems, our approach reveals significant mimicking between Edual and the DFT-computed intermolecular interaction energy profiles. We further evaluate the influence of the Coulomb and exchange-correlation contributions in Edual. These latter results suggest that numerous potential energy surfaces of clusters can be explored using a Sanderson-like model only based on classical interactions between molecular orbitals domains. To conclude, this study achieved a deeper understanding of the principles of cDFT and assessed, in a wider context, its efficiency in predicting the chemical reactivity.

Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes

Molecular mechanics (MM) simulations have the potential to provide detailed insights into the mechanisms of enzymes that utilize nucleotides as cofactors. In most cases, the activities of these enzymes also require the binding of divalent cations to catalytic sites. However, modeling divalent cations in MM simulations has been challenging. The inclusion of explicit polarization was considered promising, but despite improvements over nonpolarizable force fields and despite the inclusion of "Nonbonded-fix (NB-fix)" corrections, errors in interaction energies of divalent cations with proteins remain large. Importantly, the application of these models fails to reproduce the experimental structural data on Mg2+·Protein·ATP complexes. Focusing on these complexes, here we provide a systematic assessment of the polarizable AMOEBA model and recommend critical changes that substantially improve its predictive performance. Our key results are as follows. We first show that our recent revision of the AMOEBA protein model (AMOEBABIO18-HFC), which contains high field corrections (HFCs) to induced dipoles, dramatically improves Mg2+-protein interaction energies, reducing the mean absolute error (MAE) from 17 to 10 kcal/mol. This further supports the general applicability of AMOEBABIO18-HFC. The inclusion of many-body NB-fix corrections further reduces MAE to 6 kcal/mol, which amounts to less than 2% error. The errors are estimated with respect to vdW-inclusive density functional theory that we benchmark against CCSD(T) calculations and experiments. We also present a new model of ATP with revised polarization parameters to better capture its high field response, as well as new vdW and dihedral parameters. The ATP model accurately predicts experimental Mg2+-ATP binding free energy in the aqueous phase and provides new insights into how Mg2+ associates with ATP. Finally, we show that molecular dynamics (MD) simulations of Mg2+·Kinase·ATP complexes carried out with these improvements lead to a better agreement in global and local catalytic site structures between MD and X-ray crystallography.

Therapeutic Potential of B12N12-X (X = Au, Os, and Pt) Nanostructured as Effective Fluorouracil (5Fu) Drug Delivery Materials

In view of the research-substantiated comparative efficiency of nontoxic and bioavailable nanomaterials synergic with human systems for drug delivery, this work was aimed at studying the comparative efficiency of transition metal (Au, Os, and Pt)-decorated B₁₂N₁₂ nanocages in the adsorption of fluorouracil (5Fu), an antimetabolite-classed anticarcinogen administered for cancers of the breast, colon, rectum, and cervix. Three different metal-decorated nanocages interacted with 5Fu drug at the oxygen (O) and fluorine (F) sites, resulting in six adsorbent-adsorbate systems whose reactivity and sensitivity were investigated using density functional theory computation at the B3LYP/def2TZVP level of theory with special emphasis on the structural geometry, electronic, and topology analysis as well as the thermodynamic properties of the systems. While the electronic studies predicted Os@F as having the lowest and most favorable E_gp and E_ad of 1.3306 eV and -11.9 kcal/mol, respectively, the thermodynamic evaluation showed Pt@F to have the most favorable thermal energy (E), heat capacity (C_p), and entropy (ΔS) values as well as negative ΔH and ΔG while the adsorption studies showed that the greatest degree of chemisorption with E_ad magnitude of -204.5023 kcal/mol was observed in energies ranging from -12.0 to 138.4 kcal/mol with Os@F and Au@F at the lower and upper borders. The quantum theory of atoms in molecules results show that the six systems had noncovalent interactions as well as a certain degree of partial covalency but none showed covalent interaction while the noncovalent interaction analysis corroborated this by showing that the six systems had favorable interactions, though of varying degrees, with very little trace of steric hindrance or electrostatic interactions. Overall, the study showed that notwithstanding the good performance of the six adsorbent systems considered, the Pt@F and Os@F showed the most favorable potential for the delivery of 5Fu.

Dissecting Bonding Interactions in Cysteine Dimers

Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.

Microsolvation versus Encapsulation in Mono, Di, and Trivalent Cations

The effects of the formal charge in the stability and bonding of water cavities when solvating a cation are studied here using [X(H₂ O)₂₀ ]^q+ clusters starting with the well known 5¹² isomer of (water)₂₀ , placing a single mono, di, or trivalent X^q+ cation at the interior, and then optimizing and characterizing the resulting clusters. Highly correlated interaction and deformation energies are calculated using the CCSD(T)-DLPNO formalism. Bonding interactions are characterized using the tools provided by the quantum theory of atoms in molecules, natural bond orbitals, and non-covalent surfaces. Our results indicate that water to water hydrogen bonds are sensibly strengthened resulting in strong cooperative effects, which amount to ≈ 2 ${ \approx 2}$  kcal/mol per hydrogen bond in the bare cavity and to larger values for the systems including the cations. Approximate encapsulation, that is, surrounding the cation by a network of hydrogen bonds akin to the well known methane clathrate seems to be preferred by cations with smaller charge densities while microsolvation, that is, cluster structures having explicit X⋯O contacts seem to be preferred by cations with larger charge densities which severely deform the cavity.

Rotamers in Crystal Structures of Xylitol, D-Arabitol and L-Arabitol

Rotamers are stereoisomers produced by rotation (twisting) about σ bonds and are often rapidly interconverting at room temperature. Xylitol-massively produced sweetener-(2R,3r,4S)-pentane-1,2,3,4,5-pentol) forms rotamers from the linear conformer by rotation of a xylitol fragment around the C2-C3 bond (rotamer 1) or the C3-C4 bond (rotamer 2). The rotamers form two distinguishable structures. Small differences in geometry of rotamers of the main carbon chain were confirmed by theoretical calculations; however, they were beyond the capabilities of the X-ray powder diffraction technique due to the almost identical unit cell parameters. In the case of rotamers of similar compounds, the rotations occurred mostly within hydroxyl groups likewise rotations in L-arabitol and D-arabitol, which are discussed in this work. Our results, supported by theoretical calculations, showed that energetic differences are slightly higher for rotamers with rotations within hydroxyl groups instead of a carbon chain.

First Experimental Quantitative Charge Density Studies of Advanced Intermediate of Vitamin D Analogues

Vitamins D are a group of fat-soluble secosteroids which play a regulatory role in the functioning of most cells. Rational design of new vitamin D analogs, of increased therapeutic potency and lowered calcemic side effects, requires high-resolution initial structures and a deep understanding of interactions with the molecular targets. In this paper, using quantum crystallography, we present the first determination of the experimental quantitative charge density of an advanced intermediate of vitamin D analogues as well as a reconstruction of the theoretical electron density of final vitamin D analogues. Application of these methods allows for topological and electrostatic interaction energy analysis. We showed that the A-ring chair conformation has a significant influence on the topological properties of vitamin D compounds. Moreover, the interactions between the CD-ring and side-chain additionally stabilize the crystal structure. These results are supported by our theoretical calculations and previous biological studies.

Insight into the Mechanisms of Low Coverage Adsorption of N-Alcohols on Single Walled Carbon Nanohorn

We report for the first time the chromatographic study of n-alcohols (from methanol to butanol) adsorption on single walled carbon nanohorn (SWCNH). Using measured temperature dependence of adsorption isotherms (373-433 K) the isosteric adsorption enthalpy is calculated and compared with the data reported for a graphite surface. It is concluded that a graphite surface is more homogeneous, and the enthalpy of adsorption on SWCNHs at zero coverage correlates well with molecular diameter and polarizability, suggesting leading role of dispersive interactions, i.e., no heteroatoms presence in the walls of SWCNH structures. Next using modern DFT approach we calculate the energy of n-alcohols interactions with a graphene sheet and with a single nanocone finally proposing a more realistic-double nanocone model. Obtained results suggest alcohols entrapping between SWCNH with OH groups located toward nanocones ends, leading to the conclusions about very promising future applications of SWCNHs in catalytic reactions with participation of n-alcohols.

Alkylammonium Cation Affinities of Nitrogenated Organobases: The Roles of Hydrogen Bonding and Proton Transfer

Alkylammonium cation affinities of 64 nitrogen-containing organobases, as well as the respective proton transfer processes from the alkylammonium cations to the base, have been computed in the gas phase by using DFT methods. The guanidine bases show the highest proton transfer values (191.9-233 kJ mol^-1 ) whereas the cis-2,2'-biimidazole presents the largest affinity towards the alkylammonium cations (>200 kJ mol^-1 ) values. The resulting data have been compared with the experimentally reported proton affinities of the studied nitrogen-containing organobases revealing that the propensity of an organobase for the proton transfer process increases linearly with its proton affinity. This work can provide a tool for designing senors for bioactive compounds containing amino groups that are protonated at physiological pH.

Computational Tools to Rationalize and Predict the Self-Assembly Behavior of Supramolecular Gels

Supramolecular gels form a class of soft materials that has been heavily explored by the chemical community in the past 20 years. While a multitude of experimental techniques has demonstrated its usefulness when characterizing these materials, the potential value of computational techniques has received much less attention. This review aims to provide a complete overview of studies that employ computational tools to obtain a better fundamental understanding of the self-assembly behavior of supramolecular gels or to accelerate their development by means of prediction. As such, we hope to stimulate researchers to consider using computational tools when investigating these intriguing materials. In the concluding remarks, we address future challenges faced by the field and formulate our vision on how computational methods could help overcoming them.

Cause, Regulation and Utilization of Dye Aggregation in Dye-Sensitized Solar Cells

As an important member of third generation solar cell, dye-sensitized solar cells (DSSCs) have the advantages of being low cost, having an easy fabrication process, utilizing rich raw materials and a high-power conversion efficiency (PCE), prompting nearly three decades as a research hotspot. Recently, increasing the photoelectric conversion efficiency of DSSCs has proven troublesome. Sensitizers, as the most important part, are no longer limited to molecular engineering, and the regulation of dye aggregation has become a widely held concern, especially in liquid DSSCs. This review first presents the operational mechanism of liquid and solid-state dye-sensitized solar cells, including the influencing factors of various parameters on device efficiency. Secondly, the mechanism of dye aggregation was explained by molecular exciton theory, and the influence of various factors on dye aggregation was summarized. We focused on a review of several methods for regulating dye aggregation in liquid and solid-state dye-sensitized solar cells, and the advantages and disadvantages of these methods were analyzed. In addition, the important application of quantum computational chemistry in the study of dye aggregation was introduced. Finally, an outlook was proposed that utilizing the advantages of dye aggregation by combining molecular engineering with dye aggregation regulation is a research direction to improve the performance of liquid DSSCs in the future. For solid-state dye-sensitized solar cells (ssDSSCs), the effects of solid electrolytes also need to be taken into account.

Application of XDM to ionic solids: The importance of dispersion for bulk moduli and crystal geometries

Dispersion corrections are essential in the description of intermolecular interactions; however, dispersion-corrected functionals must also be transferrable to hard solids. The exchange-hole dipole moment (XDM) model has demonstrated excellent performance for non-covalent interactions. In this article, we examine its ability to describe the relative stability, geometry, and compressibility of simple ionic solids. For the specific cases of the cesium halides, XDM-corrected functionals correctly predict the energy ranking of the B1 and B2 forms, and a dispersion contribution is required to obtain this result. Furthermore, for the lattice constants of the 20 alkali halides, the performance of XDM-corrected functionals is excellent, provided that the base functional's exchange enhancement factor properly captures non-bonded repulsion. The mean absolute errors in lattice constants obtained with B86bPBE-XDM and B86bPBE-25X-XDM are 0.060 Å and 0.039 Å, respectively, suggesting that delocalization error also plays a minor role in these systems. Finally, we considered the calculation of bulk moduli for alkali halides and alkaline-earth oxides. Previous claims in the literature that simple generalized gradient approximations, such as PBE, can reliably predict experimental bulk moduli have benefited from large error cancellations between neglecting both dispersion and vibrational effects. If vibrational effects are taken into account, dispersion-corrected functionals are quite accurate (4 GPa-5 GPa average error), again, if non-bonded repulsion is correctly represented. Careful comparisons of the calculated bulk moduli with experimental data are needed to avoid systematic biases and misleading conclusions.

Improved Basis-Set Incompleteness Potentials for Accurate Density-Functional Theory Calculations in Large Systems

The accurate calculation of chemical properties using density-functional theory (DFT) requires the use of a nearly complete basis set. In chemical systems involving hundreds to thousands of atoms, the cost of the calculations place practical limitations on the number of basis functions that can be used. Therefore, in most practical applications of DFT to large systems, there exists a basis-set incompleteness error (BSIE). In this article, we present the next iteration of the basis-set incompleteness potentials (BSIPs), one-electron potentials designed to correct for basis-set incompleteness error. The ultimate goal associated with the development of BSIPs is to allow the calculation of molecular properties using DFT with near-complete-basis-set results at a computational cost that is similar to a small basis set calculation. In this work, we develop BSIPs for 10 atoms in the first and second rows (H, B-F, Si-Cl) and 15 common basis sets of the Pople, Dunning, Karlsruhe, and Huzinaga types. Our new BSIPs are constructed to minimize BSIE in the calculation of reaction energies, barrier heights, noncovalent binding energies, and intermolecular distances. The BSIPs were obtained using a training set of 15 944 data points. The fitting approach employed a regularized linear least-squares method with variable selection (the LASSO method), which results in a much better fit to the training data than our previous BSIPs while, at the same time, reducing the computational cost of BSIP development. The proposed BSIPs are tested on various benchmark sets and demonstrate excellent performance in practice. Our new BSIPs are also transferable; i.e., they can be used to correct BSIE in calculations that employ density functionals other than the one used in the BSIP development (B3LYP). Finally, BSIPs can be used in any quantum chemistry program that have implemented effective-core potentials without changes to the software.

Asymptotic Pairwise Dispersion Corrections Can Describe Layered Materials Accurately

A recent study by Tawfik et al. [ Phys. Rev. Mater. 2018, 2, 034005] found that few density functionals, none of which are asymptotic pairwise dispersion methods, describe the geometry and binding of layered materials accurately. Here, we show that the exchange-hole dipole moment (XDM) dispersion model attains excellent results for graphite, hexagonal BN, and transition-metal dichalcogenides. Contrary to what has been argued, successful modeling of layered materials does not necessitate meta-GGA exchange, nonlocal correlation functionals, or the inclusion of three-body dispersion terms. Rather, a GGA functional, combined with a simple asymptotic pairwise dispersion correction, can be reliably used, provided that it properly accounts for the geometric dependence of the dispersion coefficients. The overwhelming contribution to the variation of the pairwise dispersion coefficients comes from the immediate vicinity of an atom and is already present for single layers. Longer-range and interlayer effects are examined in detail for graphite.

Modified Density Functional Dispersion Correction for Inorganic Layered MFX Compounds (M = Ca, Sr, Ba, Pb and X = Cl, Br, I)

MFX (M = Ca, Ba, Sr, Pb and X = Cl, Br, I) compounds have received considerable attention due to their technological application as X-ray detectors, pressure sensors, and optical data storage materials, when doped with rare-earth ions. MFX compounds belong to the class of layered materials with a tetragonal Matlockite crystal structure, characterized by weakly stacked double-halide layers along the crystallographic c-axis. These layers predominantly determine phase transitions, elastic, and mechanical properties. However, the correct description of the lattice parameter c is a challenge for most standard DFT functionals, which tend to overestimate the lattice parameter c. Because of the weak interactions between the halide layers, dispersion-corrected functionals seem to be a better choice. We investigated 11 different inorganic layered MFX compounds for which experimental data are available, with standard and dispersion-corrected functionals to assess their performance in reproducing the lattice parameter c, structural, and vibrational properties of the MFX compounds. Our results revealed that these functionals do not describe the weak interactions between the halide layers in a balanced way. Therefore, we modified Grimme's popular DFT-D2 dispersion correction scheme in two different ways by (i) replacing the dispersion coefficients and van der Waals radii with those of noble gas atoms or (ii) increasing the van der Waals radii of the MFX atoms up to 40%. Comparison with the available experimental data revealed that the latter approach applied to the PBE (Perdew-Burke-Ernzerhof)-D2 functional with 30% increased van der Waals radii, which we coined PBE-D2* (S_rvdw 1.30) is best suited to fine-tune the description of the weak interlayer interactions in MFX compounds, thus significantly improving the description of their structural, vibrational, and mechanical properties. Work is in progress applying this new, computationally inexpensive scheme to other inorganic layered compounds and periodic systems with weakly stacked layers.

Influence of the dipolar interactions on the relative stability in spin crossover systems

Molecules exhibiting a spin-crossover transition have been proposed for a number of applications such as molecular switches, spintronic tunable interfaces, and single molecule gates. Both the rational design of new spin-crossover systems and the improvement of the properties of the already existing ones require a theoretical understanding of the relative energy of the high (HS) and low spin state (LS) molecules in the solid-state. This has proved to be very challenging so far. Here, we shed some light on the importance of considering the symmetry and the geometry of the crystallographic cell to correctly evaluate the influence of the dipolar interactions on the relative energies of the molecular complex in both different spin states. Moreover, in the case of Fe(SCN)₂ (phen)₂ dipolar interactions are found to play an important role for the stabilization of the LS complex. © 2016 Wiley Periodicals, Inc.