Stabilization and structure calculations for noncovalent interactions in extended molecular systems based on wave function and density functional theories

Prototypical π-π dimers re-examined by means of high-level CCSDT(Q) composite ab initio methods

The benzene-ethene and parallel-displaced (PD) benzene-benzene dimers are the most fundamental systems involving π-π stacking interactions. Several high-level ab initio investigations calculated the binding energies of these dimers using the coupled-cluster with singles, doubles, and quasi-perturbative triple excitations [CCSD(T)] method at the complete basis set [CBS] limit using various approaches such as reduced virtual orbital spaces and/or MP2-based basis set corrections. Here, we obtain CCSDT(Q) binding energies using a Weizmann-3-type approach. In particular, we extrapolate the self-consistent field (SCF), CCSD, and (T) components using large heavy-atom augmented Gaussian basis sets [namely, SCF/jul-cc-pV{5,6}Z, CCSD/jul-cc-pV{Q,5}Z, and (T)/jul-cc-pV{T,Q}Z]. We consider post-CCSD(T) contributions up to CCSDT(Q), inner-shell, scalar-relativistic, and Born-Oppenheimer corrections. Overall, our best relativistic, all-electron CCSDT(Q) binding energies are ∆E_e,all,rel = 1.234 (benzene-ethene) and 2.550 (benzene-benzene PD), ∆H₀ = 0.949 (benzene-ethene) and 2.310 (benzene-benzene PD), and ∆H₂₉₈ = 0.130 (benzene-ethene) and 1.461 (benzene-benzene PD) kcal mol^-1. Important conclusions are reached regarding the basis set convergence of the SCF, CCSD, (T), and post-CCSD(T) components. Explicitly correlated calculations are used as a sanity check on the conventional binding energies. Overall, post-CCSD(T) contributions are destabilizing by 0.028 (benzene-ethene) and 0.058 (benzene-benzene) kcal mol^-1, and thus, they cannot be neglected if sub-chemical accuracy is sought (i.e., errors below 0.1 kcal mol^-1). CCSD(T)/aug-cc-pwCVTZ core-valence corrections increase the binding energies by 0.018 (benzene-ethene) and 0.027 (benzene-benzene PD) kcal mol^-1. Scalar-relativistic and diagonal Born-Oppenheimer corrections are negligibly small. We use our best CCSDT(Q) binding energies to evaluate the performance of MP2-based, CCSD-based, and lower-cost composite ab initio procedures for obtaining these challenging π-π stacking binding energies.

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Efficient Density-Fitted Explicitly Correlated Dispersion and Exchange Dispersion Energies

The leading-order dispersion and exchange-dispersion terms in symmetry-adapted perturbation theory (SAPT), E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾, suffer from slow convergence to the complete basis set limit. To alleviate this problem, explicitly correlated variants of these corrections, E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12, have been proposed recently. However, the original formalism (M., Kodrycka , J. Chem. Theory Comput. 2019, 15, 5965-5986), while highly successful in terms of improving convergence, was not competitive to conventional orbital-based SAPT in terms of computational efficiency due to the need to manipulate several kinds of two-electron integrals. In this work, we eliminate this need by decomposing all types of two-electron integrals using robust density fitting. We demonstrate that the error of the density fitting approximation is negligible when standard auxiliary bases such as aug-cc-pVXZ/MP2FIT are employed. The new implementation allowed us to study all complexes in the A24 database in basis sets up to aug-cc-pV5Z, and the E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 values exhibit vastly improved basis set convergence over their conventional counterparts. The well-converged E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 numbers can be substituted for conventional E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾ ones in a calculation of the total SAPT interaction energy at any level (SAPT0, SAPT2+3, ...). We show that the addition of F12 terms does not improve the accuracy of low-level SAPT treatments. However, when the theory errors are minimized in high-level SAPT approaches such as SAPT2+3(CCD)δMP2, the reduction of basis set incompleteness errors thanks to the F12 treatment substantially improves the accuracy of small-basis calculations.

The origin of the regiospecificity of acrolein dimerization

Acrolein dimerization is a intriguing case since the reaction does not occur to form the electronically preferred regioisomeric adduct. Various explanations have been suggested to rationalize this experimental regioselectivity, however, none of these arguments had been convincing enough. In this work, the hetero Diels-Alder acrolein dimerization was theoretically investigated using DFT and MP2 methods. The influence of nucleophilic/electrophilic interactions and non-covalent interactions (NCI) in the regiospecificity of the reaction were analyzed. Our results show that the NCI at the transition state are the key factor controlling the regiospecificity in this reaction. Besides, we found that the choice of calculation method can have an effect on the prediction of the mechanism in the reaction, as all DFT methods forecast a one-step hetero Diels-Alder acrolein dimerization, while MP2 predicts a stepwise description for the lower energy reaction channel.

The hierarchy of ab initio and DFT methods for describing an intramolecular non-covalent SiN contact in the silicon compounds using electron diffraction geometries

In the series of silatranes XSi(OCH2CH2)3N, 1 (X = Me, 1a; H, 1b; F, 1c) with the known gas electron diffraction (GED) structures, the problematic geometry of 1-methylsilatrane 1a has been revised. In particular, the new value of the SiN distance (dSiN) in 1a turned out to be ∼0.06 Å longer than the generally accepted one. This dSiN resolves the long-standing contradiction between the data of the structural and spectral experiments regarding the sensitivity of 1 to the medium effect. We also performed the ab initio and DFT study of the combined series of silatranes 1a-c, silylalkylamines H3Si(CH2)3NMe2 (2a) and F3SiCH2NMe2 (2b), silylhydrazines F3SiN(Me)NMe2 (2c) and F3SiN(SiMe3)NMe2 (2d), and silyloxyamines ClH2SiONMe2 (2e,f), (F3C)F2SiONMe2 (2g,h) and F3SiONMe2 (2i), in which the GED dSiN values are in a wide range of 2-3 Å. None of the involved quantum chemical methods has succeeded in reproducing all the experimental gas-phase dSiN values in 1a-c, 2a-i with an acceptable accuracy (0.01-0.03 Å). The problems of the used methods, primarily CCSD with the Pople basis sets, are caused by four molecules with the geminal SiNN and SiON fragments (2d,f-i) and dSiN < 2.3 Å. A reasonable hierarchy of computationally accessible theory levels for studying the physicochemical manifestation of the non-covalent intramolecular SiN interactions can be constructed only at dSiN > 2.3 Å: MP2 < PBE0 ∼ B3PW91 ∼ SCS-MP2 < CCSD < CCSD(T).

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Recent Advances in the Fabrication and Application of Graphene Microfluidic Sensors

This review reports the progress of the recent development of graphene-based microfluidic sensors. The introduction of microfluidics technology provides an important possibility for the advance of graphene biosensor devices for a broad series of applications including clinical diagnosis, biological detection, health, and environment monitoring. Compared with traditional (optical, electrochemical, and biological) sensing systems, the combination of graphene and microfluidics produces many advantages, such as achieving miniaturization, decreasing the response time and consumption of chemicals, improving the reproducibility and sensitivity of devices. This article reviews the latest research progress of graphene microfluidic sensors in the fields of electrochemistry, optics, and biology. Here, the latest development trends of graphene-based microfluidic sensors as a new generation of detection tools in material preparation, device assembly, and chip materials are summarized. Special emphasis is placed on the working principles and applications of graphene-based microfluidic biosensors, especially in the detection of nucleic acid molecules, protein molecules, and bacterial cells. This article also discusses the challenges and prospects of graphene microfluidic biosensors.

Effects of structural variations on π-dimer formation: long-distance multicenter bonding of cation-radicals of tetrathiafulvalene analogues

The multicenter (pancake) bonding between cation-radicals of tetramethyltetraselenafulvalene, TMTSF⁺˙, tetramethyltetrathiafulvalene, TMTTF⁺˙, and bis(ethylenedithio)-tetrathiafulvalene, ET⁺,˙ was compared to that of tetrathiafulvalene, TTF⁺˙. To minimize counter-ion effects, the cation-radical salts with weakly coordinating anions (WCA), tetrakis(3,5-trifluoromethylphenyl)borate, dodecamethylcarborane and hexabromocarborane were prepared. Solid-state (X-ray and EPR) measurements revealed diamagnetic π-dimers in the TMTSF and ET salts and the separate monomers in the TTF salts with all WCAs, while TMTTF existed as a dimer in one and a monomer in two salts. The variable-temperature UV-Vis studies of these salts in solution showed that the thermodynamics of formation of the π-bonded dimers of TMTTF⁺˙ was close to that of TTF⁺˙, while TMTSF⁺˙ and ET⁺˙ showed a higher propensity for π-dimerization. These data indicated that the replacement of sulfur with heavier selenium or insertion of ethylenedithia-substituents into the TTF core increases the π-dimers' stability. Yet, computational analysis indicated that the weakly covalent component of π-bonding decreases in the order TTF > TMTTF > TMTSF > ET. The higher stability of the π-dimers of TMTSF⁺˙ and ET⁺˙ cation-radicals was related to a decrease of the electrostatic repulsion between cationic counter-parts and an increase of dispersion components in these associations.

Simulating Periodic Systems on a Quantum Computer Using Molecular Orbitals

The variational quantum eigensolver (VQE) is one of the most appealing quantum algorithms to simulate electronic structure properties of molecules on near-term noisy intermediate-scale quantum devices. In this work, we generalize the VQE algorithm for simulating periodic systems. However, the numerical study of a one-dimensional (1D) infinite hydrogen chain using existing VQE algorithms shows a remarkable deviation of the ground-state energy with respect to the exact full configuration interaction (FCI) result. Here, we present two schemes to improve the accuracy of quantum simulations for periodic systems. The first one is a modified VQE algorithm, which introduces a unitary transformation of Hartree-Fock orbitals to avoid the complex wave function. The second one is combining VQE with the quantum subspace expansion approach (VQE/QSE). Numerical benchmark calculations demonstrate that both of the two schemes provide an accurate description of the potential energy curve of the 1D hydrogen chain. In addition, excited states computed with the VQE/QSE approach also agree very well with FCI results.

Tetrel Bonding and Other Non-Covalent Interactions Assisted Supramolecular Aggregation in a New Pb(II) Complex of an Isonicotinohydrazide

A new supramolecular Pb(II) complex [PbL(NO₂)]_n was synthesized from Pb(NO₃)₂, N’-(1-(pyridin-2-yl)ethylidene)isonicotinohydrazide (HL) and NaNO₂. [PbL(NO₂)]_n is constructed from discrete [PbL(NO₂)] units with an almost ideal N₂O₃ square pyramidal coordination environment around Pb(II). The ligand L⁻ is coordinated through the 2-pyridyl N-atom, one aza N-atom, and the carbonyl O-atom. The nitrite ligand binds in a κ²-O,O coordination mode through both O-atoms. The Pb(II) center exhibits a hemidirected coordination geometry with a pronounced coordination gap, which allows a close approach of two additional N-atoms arising from the N=C(O) N-atom of an adjacent molecule and from the 4-pyridyl N-atom from the another adjacent molecule, yielding a N₄O₃ coordination, constructed from two Pb–N and three Pb–O covalent bonds, and two Pb⋯N tetrel bonds. Dimeric units in the structure of [PbL(NO₂)]_n are formed by the Pb⋯N=C(O) tetrel bonds and intermolecular electrostatically enforced π⁺⋯π⁻ stacking interactions between the 2- and 4-pyridyl rings and further stabilized by C–H⋯π intermolecular interactions, formed by one of the methyl H-atoms and the 4-pyridyl ring. These dimers are embedded in a 2D network representing a simplified uninodal 3-connected fes (Shubnikov plane net) topology defined by the point symbol (4∙8²). The Hirshfeld surface analysis of [PbL(NO₂)] revealed that the intermolecular H⋯X (X = H, C, N, O) contacts occupy an overwhelming majority of the molecular surface of the [PbL(NO₂)] coordination unit. Furthermore, the structure is characterized by intermolecular C⋯C and C⋯N interactions, corresponding to the intermolecular π⋯π stacking interactions. Notably, intermolecular Pb⋯N and, most interestingly, Pb⋯H interactions are remarkable contributors to the molecular surface of [PbL(NO₂)]. While the former contacts are due to the Pb⋯N tetrel bonds, the latter contacts are mainly due to the interaction with the methyl H-atoms in the π⋯π stacked [PbL(NO₂)] molecules. Molecular electrostatic potential (MEP) surface calculations showed marked electrostatic contributions to both the Pb⋯N tetrel bonds and the dimer forming π⁺⋯π⁻ stacking interactions. Quantum theory of atoms in molecules (QTAIM) analyses underlined the tetrel bonding character of the Pb⋯N interactions. The manifold non-covalent interactions found in this supramolecular assembly are the result of the proper combination of the polyfunctional multidentate pyridine-hydrazide ligand and the small nitrito auxiliary ligand.

Isotope Effects Induced by Molecular Compression

Zero-point energies (ZPEs) of hydroxyl ion and hydrogen and water molecules, free and compressed in C₆₀ cages, are computed; the excess energy acquired by molecules under compression is in the range 2-3 kcal/mol and depends on the isotopes. The differences in ZPE of compressed isotopic molecules strongly exceed those of the free molecules, resulting in the large deuterium and tritium isotope effects. These effects induced by compression are suggested as a probe for testing molecular compression of enzymatic sites; they may be important for understanding enormously large isotope effects observed in some enzymatic reactions, where they are attributed to the tunneling.

Conformational preferences of cationic β-peptide in water studied by CCSD(T), MP2, and DFT methods

The conformational preferences of the cationic nylon-3 βNM [(3R,4)-diaminobutanoic acid, dAba] dipeptide in water were explored as the first step to understand the mode of action of polymers of βNM against phylogenetically diverse and intrinsically drug-resistant pathogenic fungi. The CCSD(T), MP2, M06-2X, ωB97X-D, B2PLYP-D3BJ, and DSD-PBEP86-D3BJ levels of theory with various basis sets were assessed for relative energies of the 45 local minima of the cationic Ac-dAba-NHMe located at the SMD M06-2X/6-31+G(d) level of theory in water against the benchmark CCSD(T)/CBS-limit energies in water. The best performance was obtained at the double-hybrid DSD-PBEP86-D3BJ/def2-QZVP level of theory with RMSD = 0.12 kcal/mol in water. The M06-2X/def2-QZVP level of theory predicted reasonably the conformational preference with RMSD = 0.38 kcal/mol in water and may be an alternative level of theory with marginal deviations for the calculation of conformational energies of relatively longer cationic peptides in water. In particular, the H₁₄–helical structures appeared to be the most feasible conformations for the cationic Ac-dAba-NHMe populated at 48–64% by relative free energies in water. The hexamer built from the H₁₄–structure of the cationic Ac-dAba-NHMe adopted a left-handed 3₁₄-helix, which has a slightly narrower radius and a longer rise than the regular 3₁₄-helix of β-peptides. Hence, the 3₁₄-helices of oligomers or polymers of the cationic dAba residues are expected to be the active conformation to exhibit the ability to bridge between charged lipid head groups that might cause a local depression or invagination of the membrane of fungi.

Molecular Interpretation of Pharmaceuticals’ Adsorption on Carbon Nanomaterials: Theory Meets Experiments

The ability of carbon-based nanomaterials (CNM) to interact with a variety of pharmaceutical drugs can be exploited in many applications. In particular, they have been studied both as carriers for in vivo drug delivery and as sorbents for the treatment of water polluted by pharmaceuticals. In recent years, the large number of experimental studies was also assisted by computational work as a tool to provide understanding at molecular level of structural and thermodynamic aspects of adsorption processes. Quantum mechanical methods, especially based on density functional theory (DFT) and classical molecular dynamics (MD) simulations were mainly applied to study adsorption/release of various drugs. This review aims to compare results obtained by theory and experiments, focusing on the adsorption of three classes of compounds: (i) simple organic model molecules; (ii) antimicrobials; (iii) cytostatics. Generally, a good agreement between experimental data (e.g. energies of adsorption, spectroscopic properties, adsorption isotherms, type of interactions, emerged from this review) and theoretical results can be reached, provided that a selection of the correct level of theory is performed. Computational studies are shown to be a valuable tool for investigating such systems and ultimately provide useful insights to guide CNMs materials development and design.

Benchmark Experimental Gas-Phase Intermolecular Dissociation Energies by the SEP-R2PI Method

The gas-phase ground-state dissociation energy D0( S0) of an isolated and cold bimolecular complex is a fundamental measure of the intermolecular interaction strength between its constituents. Accurate D0 values are important for the understanding of intermolecular bonding, for benchmarking high-level theoretical calculations, and for the parameterization of dispersion-corrected density functionals or force-field models that are used in fields ranging from crystallography to biochemistry. We review experimental measurements of the gas-phase D0( S0) and D0( S1) values of 55 different M⋅S complexes, where M is a (hetero)aromatic molecule and S is a closed-shell solvent atom or molecule. The experiments employ the triply resonant SEP-R2PI laser method, which involves M-centered ( S0 → S1) electronic excitation, followed by S1 → S0 stimulated emission spanning a range of S0 state vibrational levels. At sufficiently high vibrational energy, vibrational predissociation of the M⋅S complex occurs. A total of 49 dissociation energies were bracketed to within ≤1.0 kJ/mol, providing a large experimental database of accurate noncovalent interactions.

Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods

Molecular crystal structure prediction is increasingly being applied to study the solid form landscapes of larger, more flexible pharmaceutical molecules. Despite many successes in crystal structure prediction, van der Waals-inclusive density functional theory (DFT) methods exhibit serious failures predicting the polymorph stabilities for a number of systems exhibiting conformational polymorphism, where changes in intramolecular conformation lead to different intermolecular crystal packings. Here, the stabilities of the conformational polymorphs of o-acetamidobenzamide, ROY, and oxalyl dihydrazide are examined in detail. DFT functionals that have previously been very successful in crystal structure prediction perform poorly in all three systems, due primarily to the poor intramolecular conformational energies, but also due to the intermolecular description in oxalyl dihydrazide. In all three cases, a fragment-based dispersion-corrected second-order Møller-Plesset perturbation theory (MP2D) treatment of the crystals overcomes these difficulties and predicts conformational polymorph stabilities in good agreement with experiment. These results highlight the need for methods which go beyond current-generation DFT functionals to make crystal polymorph stability predictions truly reliable.

Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy

The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

POM-based dyes featuring rigidified bithiophene π linkers: potential high-efficiency dyes for dye-sensitized solar cells

A series of POM-based dyes with a triphenylamine electron donor group, cyanoacrylic acid electron acceptor group and different π linkers of thiophene derivatives were systematically investigated to analyze the influence of a rigidified bithiophene with fastening atoms (C and N) on the performance of dye-sensitized solar cells (DSSCs) based on density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations.

Inclusion of Van der Waals Interactions in DFT using Wannier Functions without empirical parameters

We describe a method for including van der Waals (vdW) interactions in Density Functional Theory (DFT) using the Maximally-Localized Wannier functions (MLWFs), which is free from empirical parameters. With respect to the previous DFT/vdW-WF2 version, in the present DFT/vdW-WF2-x approach, the empirical, short-range, damping function is replaced by an estimate of the Pauli exchange repulsion, also obtained by the MLWFs properties. Applications to systems contained in the popular S22 molecular database and to the case of adsorption of Ar on graphite, and Xe and water on graphene, indicate that the new method, besides being more physically founded, also leads to a systematic improvement in the description of systems where vdW interactions play a significant role.

Supramolecular lead(ii) architectures engineered by tetrel bonds

The structures, including tetrel bonding, of Pb^II coordination compounds assembled from N′-(pyridin-2-ylmethylene)picolinohydrazide, N′-(pyridin-2-ylmethylene)nicotinohydrazide and N′-(1-(pyridin-2-yl)ethylidene)isonicotinohydrazide ligands are discussed.

A new spodium bond driven coordination polymer constructed from mercury(ii) azide and 1,2-bis(pyridin-2-ylmethylene)hydrazine

In this manuscript, the synthesis and X-ray characterization of a new spodium bond driven coordination polymer constructed from mercury(ii) azide and 1,2-bis(pyridin-2-ylmethylene)hydrazine are reported.

On the importance of π-hole spodium bonding in tricoordinated HgII complexes

The synthesis and X-ray characterization of two new tri-coordinated Hg complexes where the planar Hg atom participates in π-hole spodium bonding.

Atlas of putative minima and low-lying energy networks of water clusters n = 3-25

We report a database consisting of the putative minima and ∼3.2 × 10⁶ local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as low-lying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation (MP2)/aug-cc-pVTZ level of theory. In total, we have identified 3 138 303 networks corresponding to local minima of the clusters n = 3-25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1-3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but not limited to neural networks and machine and deep learning) to understand the properties of water, nature's most important substance.

Mechanically rigid supramolecular assemblies formed from an Fmoc-guanine conjugated peptide nucleic acid

The variety and complexity of DNA-based structures make them attractive candidates for nanotechnology, yet insufficient stability and mechanical rigidity, compared to polyamide-based molecules, limit their application. Here, we combine the advantages of polyamide materials and the structural patterns inspired by nucleic-acids to generate a mechanically rigid fluorenylmethyloxycarbonyl (Fmoc)-guanine peptide nucleic acid (PNA) conjugate with diverse morphology and photoluminescent properties. The assembly possesses a unique atomic structure, with each guanine head of one molecule hydrogen bonded to the Fmoc carbonyl tail of another molecule, generating a non-planar cyclic quartet arrangement. This structure exhibits an average stiffness of 69.6 ± 6.8 N m-1 and Young's modulus of 17.8 ± 2.5 GPa, higher than any previously reported nucleic acid derived structure. This data suggests that the unique cation-free "basket" formed by the Fmoc-G-PNA conjugate can serve as an attractive component for the design of new materials based on PNA self-assembly for nanotechnology applications.

A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach

A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol^-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.

Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory

This work studies the underlying nature of H-bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry-adapted perturbation theory (SAPT). Quantum theory of atoms-in-molecules (QTAIM) descriptors and other HB indices were calculated based on high-quality density functional theory wave functions. We propose a new and rigorous classification of H-bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as "very weak" HBs with a BE ≥ -2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as "weak-to-medium" HBs with a BE varying between -2.5 and -14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into "medium" HBs with a BE in the range of -11.0 to -15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into "strong" HBs whose BE is more negative than -15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function-based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; ρ_BCP ) of HBs. The fitted equation for neutral complexes is BE/kcal/mol =  - 223.08 × ρ_BCP /a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE/kcal/mol =  - 332.34 × ρ_BCP /a. u. - 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. © 2019 Wiley Periodicals, Inc.

Electronic structure benchmark calculations of inorganic and biochemical carboxylation reactions

Carboxylation reactions represent a very special class of chemical reactions that is characterized by the presence of a carbon dioxide (CO₂ ) molecule as reactive species within its global chemical equation. These reactions work as fundamental gear to accomplish the CO₂ fixation and thus to build up more complex molecules through different technological and biochemical processes. In this context, a correct description of the CO₂ electronic structure turns out to be crucial to study the chemical and electronic properties associated with this kind of reactions. Here, a systematic study of CO₂ electronic structure and its contribution to different carboxylation reaction electronic energies has been carried out by means of several high-level ab initio post-Hartree Fock (post-HF) and density functional theory (DFT) calculations for a set of biochemistry and inorganic systems. We have found that for a correct description of the CO₂ electronic correlation energy it is necessary to include post-CCSD(T) contributions (beyond the gold standard). These high-order excitations are required to properly describe the interactions of the four π-electrons associated with the two degenerated π-molecular orbitals of the CO₂ molecule. Likewise, our results show that in some reactions it is possible to obtain accurate reaction electronic energy values with computationally less demanding methods when the error in the electronic correlation energy compensates between reactants and products. Furthermore, the provided post-HF reference values allowed to validating different DFT exchange-correlation functionals combined with different basis sets for chemical reactions that are relevant in biochemical CO₂ fixing enzymes. © 2019 Wiley Periodicals, Inc.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Prediction of self-assembly of adenosine analogues in solution: a computational approach validated by isothermal titration calorimetry

The recent discovery of the role of adenosine-analogues as neuroprotectants and cognitive enhancers has sparked interest in these molecules as new therapeutic drugs. Understanding the behavior of these molecules in solution and predicting their ability to self-assemble will accelerate new discoveries. We propose a computational approach based on density functional theory, a polarizable continuum solvation description of the aqueous environment, and an efficient search procedure to probe the potential energy surface, to determine the structure and thermodynamic stability of molecular clusters of adenosine analogues in solution, using caffeine as a model. The method was validated as a tool for the prediction of the impact of small structural variations on self-assembly using paraxanthine. The computational results were supported by isothermal titration calorimetry experiments. The thermodynamic parameters enabled the quantification of the actual percentage of dimer present in solution as a function of concentration. The data suggest that both caffeine and paraxanthine are present at concentrations comparable to the ones found in biological samples.

A Simple Model for Halogen Bond Interaction Energies

Halogen bonds are prevalent in many areas of chemistry, physics, and biology. We present a statistical model for the interaction energies of halogen-bonded systems at equilibrium based on high-accuracy ab initio benchmark calculations for a range of complexes. Remarkably, the resulting model requires only two fitted parameters, X and B—one for each molecule—and optionally the equilibrium separation, R e , between them, taking the simple form E = X B / R e n . For n = 4 , it gives negligible root-mean-squared deviations of 0.14 and 0.28 kcal mol - 1 over separate fitting and validation data sets of 60 and 74 systems, respectively. The simple model is shown to outperform some of the best density functionals for non-covalent interactions, once parameters are available, at essentially zero computational cost. Additionally, we demonstrate how it can be transferred to completely new, much larger complexes and still achieve accuracy within 0.5 kcal mol - 1 . Using a principal component analysis and symmetry-adapted perturbation theory, we further show how the model can be used to predict the physical nature of a halogen bond, providing an efficient way to gain insight into the behavior of halogen-bonded systems. This means that the model can be used to highlight cases where induction or dispersion significantly affect the underlying nature of the interaction.