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

Base pairs with 5-chloroorotic acid and comparison with the natural nucleobase. Structural and spectroscopic study, and three suggested antiviral modified nucleosides

A structural and spectroscopic study of 5-chloroorotic acid (5-ClOA) biomolecule was carried out by IR and FT-Raman and the results obtained were compared to those achieved in 5-fluoroorotic acid and 5-aminoorotic acid compounds. The structures of all possible tautomeric forms were determined using DFT and MP2 methods. To know the tautomer form present in the solid state, the crystal unit cell was optimized through dimer and tetramer forms in several tautomeric forms. The keto form was confirmed through an accurate assignment of all the bands. For this purpose, an additional improvement in the theoretical spectra was carried out using linear scaling equations (LSE) and polynomic equations (PSE) deduced from uracil molecule. Base pairs with uracil, thymine and cytosine nucleobases were optimized and compared to the natural Watson-Crick (WC) pairs. The counterpoise (CP) corrected interaction energies of the base pairs were also calculated. Three nucleosides were optimized based on 5-ClOA as nucleobase, and their corresponding WC pairs with adenosine. These modified nucleosides were inserted in DNA:DNA and RNA:RNA microhelices, which were optimized. The position of the -COOH group in the uracil ring of these microhelices interrupts the DNA/RNA helix formation. Because of the special characteristic of these molecules they can be used as antiviral drugs.Communicated by Ramaswamy H. Sarma.

Styrylbenzoquinoline dyads as a new type of fluorescing photochromes operating via [2 + 2] photocycloaddition mechanism: Optimization of the structure

We synthesized and studied a novel bichromophoric dyad in which bridging methylene groups link two styrylbenzo[f]quinoline (SBQ) photochromes to a salicylic acid residue. The dyad was designed for use as a fluorescent P-type photochrome acting via a [2 + 2] photocycloaddition (PCA) reaction. Compared to previously studied dyads, a change in the attachment handle and shortening of the bridging groups resulted in simultaneous rise of the quantum yields of both fluorescence and PCA. Under light irradiation, two competitive reversible reactions occurred in the dyad. The first is photoisomerization between the trans- and cis- isomers of the SBQ moieties. The second is PCA. The latter process was predominant and resulted in the formation of the cyclobutane ring bearing two benzo[f]quinoline (BQ) groups. In the ground S0 state, NMR data and DFT calculations indicated the formation of folded dyad conformers whose structure is pre-organized for PCA due to π-stacking interactions of two SBQ moieties. In the excited dyad, steady-state and time-resolved nanosecond fluorescence spectroscopy revealed the formation of an excimer, which was assumed to be a precursor of cyclobutane. Due to the fluorescence properties of SBQ and BQ, both dyad and cyclobutane fluoresce and can serve as a color-correlated multicolor fluorescence photoswitch. A simple approach is proposed for predicting the relationship between the spectral properties of the dyad and cyclobutane, which are the open and closed isomers of a new type of photochromes. The approach uses the dependence of the position of the maximum of the absorption band of an aromatic compound on the size of the π-system, as well as the fact that the sizes of the π-systems of the dyad and cyclobutane are related by a simple relation.

Quantum-centric high performance computing for quantum chemistry

High performance computing (HPC) is renowned for its capacity to tackle complex problems. Meanwhile, quantum computing (QC) provides a potential way to accurately and efficiently solve quantum chemistry problems. The emerging field of quantum-centric high performance computing (QCHPC), which merges these two powerful technologies, is anticipated to enhance computational capabilities for solving challenging problems in quantum chemistry. The implementation of QCHPC for quantum chemistry requires interdisciplinary research and collaboration across multiple fields, including quantum chemistry, quantum physics, computer science and so on. This perspective provides an introduction to the quantum algorithms that are suitable for deployment in QCHPC, focusing on conceptual insights rather than technical details. Parallel strategies to implement these algorithms on quantum-centric supercomputers are discussed. We also summarize high performance quantum emulating simulators, which are considered a viable tool to explore QCHPC. We conclude with challenges and outlooks in this field.

Theoretical Study of Cyanidin-Resveratrol Copigmentation by the Functional Density Theory

Anthocyanins are colored water-soluble plant pigments. Upon consumption, anthocyanins are quickly absorbed and can penetrate the blood-brain barrier (BBB). Research based on population studies suggests that including anthocyanin-rich sources in the diet lowers the risk of neurodegenerative diseases. The copigmentation caused by copigments is considered an effective way to stabilize anthocyanins against adverse environmental conditions. This is attributed to the covalent and noncovalent interactions between colored forms of anthocyanins (flavylium ions and quinoidal bases) and colorless or pale-yellow organic molecules (copigments). The present work carried out a theoretical study of the copigmentation process between cyanidin and resveratrol (CINRES). We used three levels of density functional theory: M06-2x/6-31g+(d,p) (d3bj); ωB97X-D/6-31+(d,p); APFD/6-31+(d,p), implemented in the Gaussian16W package. In a vacuum, the CINRES was found at a copigmentation distance of 3.54 Å between cyanidin and resveratrol. In water, a binding free energy ∆G was calculated, rendering -3.31, -1.68, and -6.91 kcal/mol, at M06-2x/6-31g+(d,p) (d3bj), ωB97X-D/6-31+(d,p), and APFD/6-31+(d,p) levels of theory, respectively. A time-dependent density functional theory (TD-DFT) was used to calculate the UV spectra of the complexes and then compared to its parent molecules, resulting in a lower energy gap at forming complexes. Excited states' properties were analyzed with the ωB97X-D functional. Finally, Shannon aromaticity indices were calculated and isosurfaces of non-covalent interactions were evaluated.

Spodium bonding with noble gas atoms

The nature of the bonding between a neutral group 12 member (Zn3, Cd3 and Hg3) ring and a noble gas atom was explored using quantum chemical simulations. Natural bond orbital, quantum theory of atoms in molecules, symmetry-adapted perturbation theory, and molecular electrostatic potential surface analysis were also used to investigate the type of interaction between the noble gas atom and the metal rings (Zn3, Cd3 and Hg3). The Zn3, Cd3 and Hg3 rings are bonded to the noble gas through non-covalent interactions, which was revealed by the non-covalent interaction index. Additionally, energy decomposition analysis reveals that dispersion energy is the key factor in stabilizing these systems.

Uracil-water interaction revisited - in search of single H-bonded secondary minima

Monohydrated uracil (UW) complexes are stabilized by both O⋯HO and NH⋯O hydrogen bonds (H-bonds), simultaneously participating in forming three stable cyclic structures. The role and contribution of these individual H-bonds (O⋯HO and NH⋯O) to the stability of the three UW complexes are still not understood, because of the technical problems in obtaining their optimized structures by standard geometry optimization. The present study explores a non-standard approach to identify three single H-bonded local minima structures without imaginary frequency using DFT (M06-2X, B3LYP and B3LYP-D3), MP2 and CCSD(T) theories and Dunning's correlation-consistent aug-cc-pVTZ basis set, in both vacuum and aqueous media (CPCM method). The results reveal that these new structures are very shallow secondary minima between two deep wells or next to a deep well of primary minima (double H-bonded structures) in the potential energy surface. The H-bond energy of these single H-bonded complexes is found to be less sensitive to a wide range (about 15-20 degrees) of O⋯HO and NH⋯O angles, and the linearity is preferred in the stable three single H-bonded structures. The technical method used to locate such a shallow minimum is described in detail and may be useful for identifying local minima in other cases where consecutive multiple H-bonded structures are global minima. Energy decomposition (using symmetry adapted perturbation theory (SAPT)) of interaction energy, electron redistribution, and relevant vibrational modes are discussed.

Experimental and Theoretical Survey of Intramolecular Spodium Bonds/σ/π-Holes and Noncovalent Interactions in Trinuclear Zn(II)-Salen Type Complex with OCN- Ions: A Holistic View in Crystal Engineering

In this work, one new centrosymmetric trinuclear Zn(II) complex 1, [{(OCN)Zn(L)}2Zn], using a salen-type ligand (H2L) in the presence of OCN- was synthesized and characterized via elemental, spectral, SEM-EDX, and single-crystal X-ray diffraction (SCXRD) study. In 1, SCXRD reveals two different stereochemical environments of zinc metal ions; one terminal Zn(II) center adopts square pyramidal geometries with the Addison parameter (τ) 0.095, and the central Zn(II) is tetracoordinated tetrahedral geometry. This article provides evidence of the significance and presence of spodium bonds (SpBs) in solid-state crystal structures involving a pseudotetrahedral environment of the central Zn-atom. X-ray structures reveal intramolecular Zn···O SpBs caused by the methoxy (-OCH3) substituent O-atom adjacent to the coordinated phenoxy O-atom. These noncovalent interactions have been thoroughly studied using density functional theory calculations at the RI-BP86[2]-D3[3]/def2-TZVP level of theory that characterizes the nature of SpBs, including the Baders quantum theory of atoms-in-molecules "QTAIM", molecular electrostatic potential (MEP) surface, and noncovalent index plot (NCI). In addition, a unique complex-isomer-based theoretical model has been vividly employed to estimate the SpBs energy in the complex. Natural bond orbital (NBO) analysis also tries to establish the differentiation between σ-hole and π-hole SpBs' natures more authentically. The complex energy frameworks were used to investigate noncovalent interactions. To better understand the different intermolecular interactions, we conducted a Hirshfeld surface, which revealed N···H (15.4%) and O···H (9.1%) contacts and Zn···O (5.1%) (SpBs).

Cooperativity between H-bonds and tetrel bonds. Transformation of a noncovalent C⋯N tetrel bond to a covalent bond

The dimers and trimers formed by imidazole (IM) and F2TO (T = C, Si, Ge) are studied by ab initio calculations. IM can engage in either a NH⋯O H-bond with F2TO or a T⋯N tetrel bond (TB) with the π-hole above the T atom. The latter is a true noncovalent TB for T = C but is a much shorter and stronger covalent bond with F2SiO or F2GeO. When a second IM is added, the cooperativity emerging from its H-bond with the first IM makes it a stronger nucleophile, leading to two minima with F2CO. The first structure contains a long noncovalent C⋯O TB and there is a much shorter covalent bond in the other, with a small energy barrier separating them. The same sort of double minimum occurs when the two IM units are situated parallel to one another in a stacked geometry.

Carbon dioxide adsorption to UiO-66: theoretical analysis of binding energy and NMR properties

UiO-66 is one of the most valuable metal-organic frameworks because of its excellent adsorption capability for gas molecules and its high stability towards water. Herein we investigated adsorption of carbon dioxide (CO2), acetone, and methanol to infinite UiO-66 using DFT calculations on an infinite system under periodic-boundary conditions and post-Hartree-Fock (SCS-MP2 and MP2.5) calculations on cluster models. Three to four molecules are adsorbed at each of four μ-OH groups bridging three Zr atoms in one unit cell (named Site I). Six molecules are adsorbed around three pillar ligands, where the molecule is loosely surrounded by three terephthalate ligands (named Site II). Also, six molecules are adsorbed around the pillar ligand in a different manner from that at Site II, where the molecule is surrounded by three terephthalate ligands (named Site III). Totally fifteen to sixteen CO2 molecules are adsorbed into one unit cell of UiO-66. The binding energy (BE) decreases in the order Site I > Site III > Site II for all three molecules studied here and in the order acetone > methanol ≫ CO2 in the three adsorption sites. At the site I, the protonic H atom of the μ-OH group interacts strongly with the negatively charged O atom of CO2, acetone and methanol, which is the origin of the largest BE value at this site. Although the DFT calculations present these decreasing orders of BE values correctly, the correction by post-Hartree-Fock calculations is not negligibly small and must be added for obtaining better BE values. We explored NMR spectra of UiO-66 with adsorbed CO2 molecules and found that the isotropic shielding constants of the 1H atom significantly differ among no CO2, one CO2 (at Sites I, II, or III), and fifteen CO2 adsorption cases (Sites I to III) but the isotropic 17O and 13C shielding constants change moderately by adsorption of fifteen CO2 molecules. Thus, 1H NMR measurement is a useful experiment for investigating CO2 adsorption.

Structural Study of a La(III) Complex of a 1,2,3-Triazole Ligand with Antioxidant Activity

The 1,2,3-triazole derivative 2-(4-chlorophenyl)-5-(pyrrolidin-1-yl)-2H-1,2,3-triazole-4-carboxylic acid with potential anticancer activity was used as a ligand in complex formation with the lanthanum(III) ion. The molecular structure and vibrational spectra of the complex were optimized at three DFT levels, and the scaled IR and Raman spectra were compared to the experimental ones. Several scaling procedures were used. Through a detailed analysis, the structure predicted for the newly synthetized La(III) complex was confirmed by the good accordance of the calculated/experimental IR and Raman spectra. The best DFT method appeared to be M06-2X with the Lanl2mb basis set, followed closely by Lanl2dz. The effect of the lanthanide atom on the molecular structure and atomic charge distribution of the triazole ring was evaluated. The potential free radical scavenging activity of both the ligand and the complex was investigated in several radical-generating model systems. The potential mechanisms of antioxidant action (hydrogen atom transfer (HAT) and single-electron transfer (SET)) were elucidated.

Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions

Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.

Graphene as Nanocarrier for Gold(I)-Monocarbene Complexes: Strength and Nature of Physisorption

Gold(I) metal complexes are finding increasing applications as therapeutic agents against a variety of diseases. As their potential use as effective metallodrugs is continuously confirmed, the issue of their administration, distribution and delivery to desired biological targets emerges. Graphene and its derivatives possess attractive properties in terms of high affinity and low toxicity, suggesting that they can efficaciously be used as drug nanocarriers. In the present study, we computationally address the adsorption of a gold(I) N-heterocyclic monocarbene, namely, IMeAuCl (where IMe = 1,3-dimethylimidazol-2-ylidene), on graphene. The Au(I) N-heterocyclic carbene family has indeed shown promising anticancer activity and the N-heterocyclic ring could easily interact with planar graphene nanostructures. By means of high-level electronic structure approaches, we investigated the strength and nature of the involved interaction using small graphene prototypes, which allow us to benchmark the best-performing DFT functionals as well as assess the role of the different contributions to total interaction energies. Moreover, realistic adsorption enthalpies and free energy values are obtained by exploiting the optimal DFT method to describe the drug adsorption on larger graphene models. Such values (ΔH_ads = -18.4 kcal/mol and ΔG_ads= -7.20 kcal/mol for the largest C₁₅₀H₃₀ model) indicate a very favorable adsorption, mainly arising from the dispersion component of the interaction, with the electrostatic attraction also playing a non-negligible role.

Tensor Hypercontraction Form of the Perturbative Triples Energy in Coupled-Cluster Theory

We present the working equations for a reduced-scaling method of evaluating the perturbative triples (T) energy in coupled-cluster theory, through the tensor hypercontraction (THC) of the triples amplitudes (t_ijk^abc). Through our method, we can reduce the scaling of the (T) energy from the traditional O(N7) to a more modest O(N5). We also discuss implementation details to aid future research, development, and software realization of this method. Additionally, we show that this method yields submillihartree (mEh) differences from CCSD(T) when evaluating absolute energies and sub-0.1 kcal/mol energy differences when evaluating relative energies. Finally, we demonstrate that this method converges to the true CCSD(T) energy through the systematic increasing of the rank or eigenvalue tolerance of the orthogonal projector, as well as exhibiting sublinear to linear error growth with respect to system size.

Accurate Prediction of Three-Body Intermolecular Interactions via Electron Deformation Density-Based Machine Learning

This work extends the electron deformation density-based descriptor, originally developed in the electron deformation density-based interaction energy machine learning (EDDIE-ML) algorithm to predict dimer interaction energies, to the prediction of three-body interactions in trimers. Using a sequential learning process to select the training data, the resulting Gaussian process regression (GPR) model predicts the three-body interaction energy within 0.2 kcal mol^-1 of the SRS-MP2/cc-pVTZ reference values for the 3B69 and S22-3 trimer data sets. A hybrid kernel function is introduced, which combines contributions from the average and individual atomic environments, allowing the total trimer interaction energy to be predicted in addition to the three-body contribution using the same descriptor. To extend the range and diversity of trimer interaction energies available in the literature, a new data set based on a protein-ligand crystal structure is introduced, consisting of 509 structures of a central ligand with two protein fragments. Benchmark calculations are provided for the new data set, which contains significantly larger molecular interactions than current databases in the literature in addition to charged fragments. Compared to density funtional theory (DFT)- and wavefunction-based methods for calculating the three-body interaction energy, our model makes predictions in a significantly shorter time frame by reducing the number of required SCF calculations from 7 to 4 performed at the PBE0 level of theory, showcasing the utility and efficiency of our Δ-ML method particularly when applied to larger systems.

Beryllium bonding with noble gas atoms

Quantum chemical calculations were carried out to investigate the nature of the bonding between a neutral Be₃ ring and noble gas atom. Electronic structure calculation for these complexes was carried out at different computational levels in association with natural bond orbital, quantum theory of atoms in molecules, electron localization function, symmetry adapted perturbation theory, and molecular electrostatic potential surface analysis of Be₃ complexes. The Be atoms in the Be₃ moiety are chemically bonded to one another, with the BeBe bond dissociation energy being ~125 kJ mol^-1 . The Be₃ ring interacts with the noble gases through non-covalent interactions. The binding energies of the noble gas atoms with the Be₃ ring increases with increase in their atomic number. The non-covalent interaction index, density overlap region indicator and independent gradient model analyses reveal the presence of non-covalent inter-fragment interactions in the complexes. Energy decomposition analysis reveals that dispersion plays the major role towards stabilizing these systems.

How proton transfer affects the helical parameters in DNA:DNA microhelices

Proton transfer reactions are a widespread phenomenon in many areas of the life sciences and it is one of the origins of the spontaneous point mutations during DNA replication. Because of its importance, many studies have been reported on these reactions. However, the present work is the first one focused on the structural geometrical changes by double proton transfer (DPT). Thus, different Watson-Crick (WC) pairs were optimized first in a simple model with one nucleoside base pair, and in a microhelix form with three nucleoside base pairs. The canonical and few tautomeric forms were considered in DNA:DNA microhelices with A-type and B-type helical forms. The stability of these structures and how the DPT process affects the main geometrical parameters was analyzed, in particular the deformation of the helical parameters. The M06-2X DFT method was used for this purpose. The purine/pyrimidine ring in the keto form appears easier to be deformed than when it is in the enol form. The weaker WC base pair formed with mixed microhelices than with nucleobases alone and the significant deformation of the helical and backbone parameters with the DPT appears to complicate this process in microhelices.Communicated by Ramaswamy H. Sarma.

A Cost Effective Scheme for the Highly Accurate Description of Intermolecular Binding in Large Complexes

There has been a growing interest in quantitative predictions of the intermolecular binding energy of large complexes. One of the most important quantum chemical techniques capable of such predictions is the domain-based local pair natural orbital (DLPNO) scheme for the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)], whose results are extrapolated to the complete basis set (CBS) limit. Here, the DLPNO-based focal-point method is devised with the aim of obtaining CBS-extrapolated values that are very close to their canonical CCSD(T)/CBS counterparts, and thus may serve for routinely checking a performance of less expensive computational methods, for example, those based on the density-functional theory (DFT). The efficacy of this method is demonstrated for several sets of noncovalent complexes with varying amounts of the electrostatics, induction, and dispersion contributions to binding (as revealed by accurate DFT-based symmetry-adapted perturbation theory (SAPT) calculations). It is shown that when applied to dimeric models of poly(3-hydroxybutyrate) chains in its two polymorphic forms, the DLPNO-CCSD(T) and DFT-SAPT computational schemes agree to within about 2 kJ/mol of an absolute value of the interaction energy. These computational schemes thus should be useful for a reliable description of factors leading to the enthalpic stabilization of extended systems.

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys-Bernardi function counterpoise (CP) procedure "overcorrects" supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ∼150 atoms and ligand-protein models with ∼300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. This is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP "overcorrects" for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

Environmental Control of the Magnetic Behavior of Transition Metal Complexes: Density Functional Theory Study of Zeolite Y Embedded Complexes [M(bpy)3]2+@Y (M = Fe2+, Co2+)

Using the supramolecular approach developed for the study of the guest-host interactions in the zeolite Y encapsulated [Fe(bpy)₃]²⁺ compound: [Fe(bpy)₃]²⁺@Y (bpy = 2,2'-bipyridine) [Vargas et al., J. Chem. Theory Comput. 2009, 5, 97-115], we apply density functional theory (DFT) to the study of the influence of zeolite Y encapsulation on the structural and energetic properties of [Co(bpy)₃]²⁺ in the low-spin (LS) and high-spin (HS) states, while revisiting [Fe(bpy)₃]²⁺@Y. Although the accurate prediction of the HS-LS energy difference ΔE_HL^el remains challenging for current DFT methods, they give accurate estimates of its variation Δ(ΔE_HL^el) in a series of complexes of a given transition metal ion. Therefore, denoting [M(bpy)₃]²⁺@Y_SM as the supramolecular model of the inclusion compounds, the values of ΔE_HL^el for the bpy complexes in the gas phase and in the supercage of zeolite Y were determined by combining the DFT estimates of Δ(ΔE_HL^el) in the series {[M(NCH)₆]²⁺, [M(bpy)₃]²⁺, and [M(bpy)₃]²⁺@Y_SM}, with accurate CCSD(T) estimates of ΔE_HL^el in the benchmark complexes [M(NCH)₆]²⁺ (M = Fe, Co) [Lawson Daku et al., J. Chem. Theory Comput., 2012, 8, 4216-4231]. Generalized gradient approximations as well as global and range-separated hybrids were employed. In order to better account for the key role of dispersion, they were also augmented with the semiempirical D2, D3BJ, and D3BJM dispersion corrections when available. The use of the D3BJ and D3BJM corrections led to similar results, and this is only with the use of the D2 scheme that (i) the free and encapsulated [Fe(bpy)₃]²⁺ are correctly predicted as LS species and that (ii) the encapsulation of both complexes translates into a destabilization of their HS state with respect to their LS state. The increase of the HS-LS energy difference is smaller for [Co(bpy)₃]²⁺ than [Fe(bpy)₃]²⁺ because the HS-LS molecular volume difference ΔV_HL in [Co(bpy)₃]²⁺ is ∼50% smaller than in [Fe(bpy)₃]²⁺. Periodic DFT calculations performed on crystalline [M(bpy)₃]²⁺@Y show that the employed [M(bpy)₃]²⁺@Y_SM supramolecular model allows the influence of encapsulation on the geometry and the spin-state energetics of [M(bpy)₃]²⁺ (M = Fe, Co) to be quantitatively captured.

Probing and Manipulating Noncovalent Interactions in Functional Polymeric Systems

Noncovalent interactions, which usually feature tunable strength, reversibility, and environmental adaptability, have been recognized as driving forces in a variety of biological and chemical processes, contributing to the recognition between molecules, the formation of molecule clusters, and the establishment of complex structures of macromolecules. The marriage of noncovalent interactions and conventional covalent polymers offers the systems novel mechanical, physicochemical, and biological properties, which are highly dependent on the binding mechanisms of the noncovalent interactions that can be illuminated via quantification. This review systematically discusses the nanomechanical characterization of typical noncovalent interactions in polymeric systems, mainly through direct force measurements at microscopic, nanoscopic, and molecular levels, which provide quantitative information (e.g., ranges, strengths, and dynamics) on the binding behaviors. The fundamental understandings of intermolecular and interfacial interactions are then correlated to the macroscopic performances of a series of noncovalently bonded polymers, whose functions (e.g., stimuli-responsiveness, self-healing capacity, universal adhesiveness) can be customized through the manipulation of the noncovalent interactions, providing insights into the rational design of advanced materials with applications in biomedical, energy, environmental, and other engineering fields.

The hitchhiker's guide to dynamic ion-solvent clustering: applications in differential ion mobility spectrometry

This article highlights the fundamentals of ion-solvent clustering processes that are pertinent to understanding an ion's behaviour during differential mobility spectrometry (DMS) experiments. We contrast DMS with static-field ion mobility, where separation is affected by mobility differences under the high-field and low-field conditions of an asymmetric oscillating electric field. Although commonly used in mass spectrometric (MS) workflows to enhance signal-to-noise ratios and remove isobaric contaminants, the chemistry and physics that underpins the phenomenon of differential mobility has yet to be fully fleshed out. Moreover, we are just now making progress towards understanding how the DMS separation waveform creates a dynamic clustering environment when the carrier gas is seeded with the vapour of a volatile solvent molecule (e.g., methanol). Interestingly, one can correlate the dynamic clustering behaviour observed in DMS experiments with gas-phase and solution-phase molecular properties such as hydrophobicity, acidity, and solubility. However, to create a generalized, global model for property determination using DMS data one must employ machine learning. In this article, we provide a first-principles description of differential ion mobility in a dynamic clustering environment. We then discuss the correlation between dynamic clustering propensity and analyte physicochemical properties and demonstrate that analytes exhibiting similar ion-solvent interactions (e.g., charge-dipole) follow well-defined trends with respect to DMS clustering behaviour. Finally, we describe how supervised machine learning can be used to create predictive models of molecular properties using DMS data. We additionally highlight open questions in the field and provide our perspective on future directions that can be explored.

The tetrel bonding role in supramolecular aggregation of lead(II) acetate and a thiosemicarbazide derivative

A new Pb^II coordination complex [PbL(OAc)], which was readily synthesized from a mixture of Pb(OAc)₂·3H₂O and 1-(pyridin-2-yl)benzylidene-4-phenylthiosemicarbazide (HL) is reported. The crystal structure analysis of [PbL(OAc)] showed that the Pb^II cation is N,N',S-chelated by the tridentate pincer-type ligand L and by the oxygen atoms of the acetate anion. In addition, the metal centre forms Pb...O and Pb...S tetrel bonds with an adjacent complex molecule, yielding a 1D zigzag polymeric chain, which is reinforced by N-H...O hydrogen bonds and π...π interactions. These chains are interlinked by C-H...py non-covalent interactions, realized between one of the acetate hydrogen atoms and the pyridine rings. According to the Hirshfeld surface analysis, the crystal packing is mainly characterized by intermolecular H...H, H...C and H...O contacts, followed by H...N, H...S, C...C, C...N, Pb...H, Pb...O and Pb...S contacts. The FTIR and ¹H NMR spectra of [PbL(OAc)] testify to the deprotonation of the parent ligand HL, while the acetate ligand exhibits an anisobidentate coordination mode as established by means of single-crystal X-ray diffraction and FTIR spectroscopy. Lastly, theoretical calculations at the PBE0-D3/def2-TZVP level of theory have been used to analyze and characterize the Pb...O and Pb...S tetrel bonds observed in the crystal of [PbL(OAc)], using a combination of QTAIM (Quantum Theory of Atoms in Molecules) and NCIPlot (Non-Covalent Interaction Plot) computational tools.

Performance of Dispersion-Inclusive Density Functional Theory Methods for Energetic Materials

Molecular crystals of energetic materials (EMs) are denser than typical molecular crystals and are characterized by distinct intermolecular interactions between nitrogen-containing moieties. To assess the performance of dispersion-inclusive density functional theory (DFT) methods, we have compiled a data set of experimental sublimation enthalpies of 31 energetic materials. We evaluate the performance of three methods: the semilocal Perdew-Burke-Ernzerhof (PBE) functional coupled with the pairwise Tkatchenko-Scheffler (TS) dispersion correction, PBE with the many-body dispersion (MBD) method, and the PBE-based hybrid functional (PBE0) with MBD. Zero-point energy contributions and thermal effects are described using the quasi-harmonic approximation (QHA), including explicit treatment of thermal expansion, which we find to be non-negligible for EMs. The lattice energies obtained with PBE0+MBD are the closest to experimental sublimation enthalpies with a mean absolute error of 9.89 kJ/mol. However, the state-of-the-art treatment of vibrational and thermal contributions makes the agreement with experiment worse. Pressure-volume curves are also examined for six representative materials. For pressure-volume curves, all three methods provide reasonable agreement with experimental data with mean absolute relative errors of 3% or less. Most of the intermolecular interactions typical of EMs, namely nitro-amine, nitro-nitro, and nitro-hydrogen interactions, are more sensitive to the choice of the dispersion method than to the choice of the exchange-correlation functional. The exception is π-π stacking interactions, which are also very sensitive to the choice of the functional. Overall, we find that PBE+TS, PBE+MBD, and PBE0+MBD do not perform as well for energetic materials as previously reported for other classes of molecular crystals. This highlights the importance of testing dispersion-inclusive DFT methods for diverse classes of materials and the need for further method development.

Analytic Gradients for the Long-Range-Corrected XYG3 Type of Doubly Hybrid Density Functionals: Theory, Implementation, and Assessment

An analytic gradient algorithm for the long-range-corrected (lrc-) XYG3 type of doubly hybrid functionals has been derived and implemented. Due to the introduction of a long-range second-order perturbation (lrPT2) as the correction, the Z-vector equations need to be modified, where the construction of a unique total Lagrangian of the method is required. Geometry optimizations using lrc-XYG3 and lrc-XYGJ-OS have then been applied to the A21 data set that consists of 21 noncovalently interacting systems with CCSD(T) reference data for structures, in order to testify the idea of introducing lrPT2 for a better description of intermolecular interactions. While lrc-XYG3 was observed to offer a better description in dispersion-dominant systems as compared to the XYG3 functional, lrc-XYGJ-OS stands out with an overall better balanced performance in reproducing intermolecular geometries.

Weakly Bound Complex Formation between HCN and CH3Cl: A Matrix-Isolation and Computational Study

The matrix-isolated infrared spectrum of a hydrogen cyanide-methyl chloride complex was investigated in a solid argon matrix. HCN and CH₃Cl were co-condensed onto a substrate held at 10 K with an excess of argon gas, and the infrared spectrum was measured using Fourier-transform infrared spectroscopy. Quantum chemical geometry optimization, harmonic frequency, and natural bonding orbital calculations indicate stabilized hydrogen- and halogen-bonded structures. The two resulting weakly bound complexes are both composed of one CH₃Cl molecule bound to a (HCN)₃ subunit, where the three HCN molecules are bound head-to-tail in a ring formation. Our study suggests that─in the presence of CH₃Cl─the formation of (HCN)₃ is promoted through complexation. Since HCN aggregates are an important precursor to prebiotic monomers (amino acids and nucleobases) and other life-bearing polymers, this study has astrophysical implications toward the search for life in space.

Generation of Basis Sets for Accurate Molecular Calculations: Application to Helium Atom and Dimer

A new approach for basis set generation is reported and tested in helium atom and dimer. The basis sets thus computed, named sigma, range from DZ to 5Z and consist of the same composition as Dunning basis sets but with a different treatment of contractions. The performance of the sigma sets is analyzed for energy and other properties of He atom and He dimer, and the results are compared with those obtained with Dunning and ANO basis sets. The sigma basis sets and their extended versions up to triple augmented provide better energy values than Dunning basis sets of the same composition, and similar values to those attained with the currently available ANO. Extrapolation to complete basis set of correlation energy is compared between the sigma basis sets and those of Dunning, showing the better performance of the former in this respect.

Quantum algorithms for electronic structures: basis sets and boundary conditions

The advantages of quantum computers are believed to significantly change the research paradigm of chemical and materials sciences, where computational characterization and theoretical design play an increasingly important role. It is especially desirable to solve the electronic structure problem, a central problem in chemistry and materials science, efficiently and accurately with well-designed quantum algorithms. Various quantum electronic-structure algorithms have been proposed in the literature. In this article, we briefly review recent progress in this direction with a special emphasis on the basis sets and boundary conditions. Compared to classical electronic structure calculations, there are new considerations in choosing a basis set in quantum algorithms. For example, the effect of the basis set on the circuit complexity is very important in quantum algorithm design. Electronic structure calculations should be performed with an appropriate boundary condition. Simply using a wave function ansatz designed for molecular systems in a material system with a periodic boundary condition may lead to significant errors. Artificial boundary conditions can be used to partition a large system into smaller fragments to save quantum resources. The basis sets and boundary conditions are expected to play a crucial role in electronic structure calculations on future quantum computers, especially for realistic systems.

The interplay of intra- and intermolecular errors in modeling conformational polymorphs

Conformational polymorphs of organic molecular crystals represent a challenging test for quantum chemistry because they require careful balancing of the intra- and intermolecular interactions. This study examines 54 molecular conformations from 20 sets of conformational polymorphs, along with the relative lattice energies and 173 dimer interactions taken from six of the polymorph sets. These systems are studied with a variety of van der Waals-inclusive density functionals theory models; dispersion-corrected spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2D); and domain local pair natural orbital coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. We investigate how delocalization error in conventional density functionals impacts monomer conformational energies, systematic errors in the intermolecular interactions, and the nature of error cancellation that occurs in the overall crystal. The density functionals B86bPBE-XDM, PBE-D4, PBE-MBD, PBE0-D4, and PBE0-MBD are found to exhibit sizable one-body and two-body errors vs DLPNO-CCSD(T) benchmarks, and the level of success in predicting the relative polymorph energies relies heavily on error cancellation between different types of intermolecular interactions or between intra- and intermolecular interactions. The SCS-MP2D and, to a lesser extent, ωB97M-V models exhibit smaller errors and rely less on error cancellation. Implications for crystal structure prediction of flexible compounds are discussed. Finally, the one-body and two-body DLPNO-CCSD(T) energies taken from these conformational polymorphs establish the CP1b and CP2b benchmark datasets that could be useful for testing quantum chemistry models in challenging real-world systems with complex interplay between intra- and intermolecular interactions, a number of which are significantly impacted by delocalization error.

Synthesis, spectroscopic findings and crystal engineering of Pb(ii)-Salen coordination polymers, and supramolecular architectures engineered by σ-hole/spodium/tetrel bonds: a combined experimental and theoretical investigation

Spontaneous self-assembly is one of the available synthetic routes to achieve structurally versatile and unique crystal complexes with selected metal-ligand combinations in the spirit of pseudohalides. In this endeavour, we designed a novel 1D coordination polymer (CP), [(Cd)(Pb)(L)(η¹-NCS)(η¹-SCN)] _n (1), using a compartmental Salen ligand (H₃L) in the presence of NaSCN. The characterization of the CP was accomplished using several spectroscopic techniques: MALDI-TOF, PXRD, SEM, EDX mapping, and single-crystal X-ray crystallography. The CP crystallizes in the monoclinic space group P2₁/c with Z = 4. SCXRD reveals Cd(ii) and Pb(ii) metal ions fulfilled distorted square pyramidal and hemi-directed coordination spheres. Cd(ii) is placed in the inner N₂O₂ and heavier Pb(ii) in the outer O₄ compartments of the de-protonated form of the ligand [L]^2-. Supramolecular interactions in the intricate crystal structure produced attractive molecular architectures of the compound. The flexible aliphatic -OH pendent group coordinates with the Pb(ii) ions. This unique binding further elevates the supramolecular crystal topographies. The supramolecular interactions were authenticated by Hirshfeld surface analysis (HSA). The observation of the recurring unconventional tetrel bonds was rationalized by DFT calculations and surface plots of molecular electrostatic potential (MEP). In the 1D polymeric chain in the complex, the O-atom of the -OH groups shows a tetrel bonding interaction with the Pb atom. We have found that the combination of QTAIM/NCI and QTAIM/ELF plots helps reveal the nature of these contacts. Moreover, the QTAIM/ELF plot determines the donor-acceptor interaction between the O-atom and the Pb atom, establishing the σ-hole. Agreeably, the σ-hole interaction also helps Pb(ii) serve as a Lewis acid in the complex. Finally, spodium and tetrel bonds are formed, possible thanks to a hemi-directional coordination sphere of the Pb atoms in the polymer described.

Emissive and reactive excimers in a covalently-linked supramolecular multi-chromophoric system with a balanced rigid-flexible structure

A novel multi-chromophoric system, triad, in which two styrylbenzoquinoline (SBQ) photochromes are connected by a balanced rigid-flexible linker comprising 2,3-naphthylene framework (a residue of 3-oxy-2-naphthoic acid) and tetramethylene groups, was designed and synthesized to study an excimer formation in the excited state. The ¹H NMR data testified that triad exists in solution as folded conformers with asymmetric parallel-displaced SBQ units. Under light irradiation, in the triad, competitive photoisomerization and [2 + 2] photocycloaddition reactions were observed, both reactions being reversible. The photocycloaddition resulted in a tetrasubstituted cyclobutane. The red-shifted fluorescence spectrum and the appearance of a long-lived component in the triad fluorescence decay indicated formation of an 'emissive' excimer. The photocycloaddition is assumed to occur in a 'reactive' excimer, in which the ethylene groups of the SBQ photochromes are located at a distance sufficient for the formation of the σ-bonds between them. Quantum-chemical density functional theory (DFT) calculations at M06-2X/6-31G* level predicted the existence of the triad conformers with π-stacking interaction of SBQ photochromes, the structure of which is pre-organized for the excimer formation and photocycloaddition. For the first time, both emissive and reactive excimers were experimentally observed in the multi-chromophoric system with two diarylethylene photochromes undergoing [2 + 2] photocycloaddition.

The effect of midbond functions on interaction energies computed using MP2 and CCSD(T)

In this article we use MP2 and CCSD(T) calculations for the A24 and S66 data sets to explore how midbond functions can be used to generate cost effective counterpoise corrected supramolecular interaction energies of noncovalent complexes. We use the A24 data set to show that the primary role of midbond functions is not to approach the complete basis set limit, but rather to ensure a balanced description of the molecules and the interaction region (unrelated to the basis set superposition error). The need for balance is a consequence of using atom centered basis sets. In the complete basis set limit, the error will disappear, but reaching the complete basis set limit is not feasible beyond small systems. For S66 we investigate the need for increasing the number of midbond centers. Results show that adding a second midbond center increases the accuracy, but the effect is secondary to changing the atom centered basis set. Further, by comparing calculations using the 3s3p2d1f1g midbond set with using aug-cc-pVDZ and aug-cc-pVTZ as midbond sets, we see that the requirements for the midbond set to be effective, is not just that it contains diffuse functions, but also that high angular momentum functions are included. By comparing two approaches for placing midbond centers we show that results are not particularly sensitive to placement as long as the placement is reasonable.

Evidence of cluster formation of croconic acid with Ag, Au and Cu cages, enhancement of electronic properties and Raman activity

Investigation of the adsorption properties croconic acid (CCA) with metal clusters (mC: Ag, Au and Cu) are reported using DFT method. CCA is found to form stable cluster with transition metal clusters of copper, silver and gold. The drug-cluster complexaton energy is slightly more for the copper nanocluster-drug complex. Non-covalent interaction analysis indicated that strong interactions and weak van der Waal interaction is present between drug and metal clusters. Dipole moment of the drug-gold cluster is found to be higher than that of the other systems. SERS studies demonstrates improved Raman signals for multiple wavenumbers of all CCA-metal cluster complexes. Mulliken charge analysis show that all CCA oxygen atom's charge changes due to the interactions with the mCs. Clustering of CCA with metal cages enhances the medicinal properties and the metal nanoclusters will act as a drug carrier of CCA.

Digital-intellectual design of microporous organic polymers

Microporous organic polymers (MOPs) are a new class of microporous materials. Due to their high porosity, large pore volume, and large surface area, MOPs exhibit excellent performance in gas adsorption and storage, membrane separation, ion capture, heterogeneous catalysis, light energy conversion and storage, capacitance, and other fields. However, selecting high-performance materials for specific applications from thousands of candidate MOPs is a key problem. Traditional design strategies for new materials with targeted properties, including trial-and-error and relying on the experiences of domain experts, are time- and cost-consuming. With the rapid development of computation technology and theoretical chemistry, the discovery of new materials is no longer a purely experimental subject. Breaking away from the traditional trial-and-error strategy for materials discovery, materials design is emerging and gaining increasing attention. In addition, the ability to collect "big data" has greatly improved and has further stimulated the development of new methods for materials design and discovery. In this perspective, we examine how data-driven techniques combine artificial intelligence (AI) and human expertise, playing a significant role in the design of MOPs. Such analytics can significantly reduce time-to-insight and accelerate the cost-effective materials discovery, which is the goal for designing future MOPs.

DFT Functionals for Modeling of Polyethylene Chains Cross-Linked by Metal Atoms. DLPNO-CCSD(T) Benchmark Calculations

Density functional theory (DFT) functionals for calculations of binding energies (BEs) of the polyethylene (PE) chains cross-linked by selected metal atoms (M) are benchmarked against DLPNO-CCSD(T) and DLPNO-CCSD(T1) data. PEX-M-PEX complexes as compared with plain parallel PEX···PEX chains with X = 3-9 carbon atoms are model species characterized by a cooperative effect of covalent C-M-C bonds and interchain dispersion interactions. The accuracy of DLPNO-CC methods was assessed by a comparison of BEs with the canonical CCSD(T) results for small PE3-M-PE3 complexes. Functionals for PEX···PEX and closed-shell PEX-M-PEX complexes (M = Be, Mg, Zn) were benchmarked against DLPNO-CCSD(T) BEs; open-shell complexes (M = Li, Ag, Au) were benchmarked against the DLPNO-CCSD(T1) method with iterative triples. Three dispersion corrections were combined with 25 DFT functionals for calculations of BEs with respect to PEX-M and PEX fragments employing def2-TZVPP and def2-QZVPP basis sets. Accuracy to within 5% for the closed-shell PEX-M-PEX complexes was achieved with five functionals. Less accurate are functionals for the open-shell PEX-M-PEX complexes; only two functionals deviate by less than 15% from DLPNO-CCSD(T1). Particularly problematic were PEX-Li-PEX complexes. A reasonable overall performance across all complexes in terms of the mean absolute percentage error is found for the range-separated hybrid functionals ωB97X-D3 and CAM-B3LYP/D3(BJ)-ABC.

Theoretical Investigation of the 4,5-Dibromorodamine Methyl Ester (TH9402) Photosensitizer Used in Photodynamic Therapy: Photophysics, Reactions in the Excited State, and Interactions with DNA

Photosensitizer (PS) molecules play a critical role in photodynamic therapy of cancer and the understanding of the molecular mechanism involved in the photophysics of these compounds, and their reactions in the excited state are, therefore, of great interest for the development of this technique. In this article, the photophysics of the cationic PS 4,5-dibromorodamine methyl ester (TH9402), its electron- and energy-transfer reactions in the excited triplet state, with molecular oxygen, nitric oxide, guanosine-5'-monophosphate (GMP), and guanine, and the interaction with DNA were evaluated. Time-dependent density functional theory calculations at the TPSSh/Def2-TZVP//B3LYP/Def2-TZVP level of theory in water solution reveals that the PS has a bright S₁ state 2.33 eV above the ground state that produces a fluorescent rate constant of 5.40 × 10⁷ s^-1, calculated using Fermi's golden rule within a path integral formalism. Once excited to the bright state, the main intersystem crossing (ISC) channel involves the coupling with the T₂ state just below S₁ (S₁ → T₂ → T₁) with an overall ISC rate constant of 10.1 × 10⁷ s^-1, in good agreement with the experimental data. Excited-state reaction thermodynamics, computed at the M06-2X/Def2-TZVP//B3LYP/Def2-TZVP level of theory in water, showed that from all the excited-state electron-transfer reactions studied, only the transfer from GMP to the PS is thermodynamically favorable, independent of the protonation state of guanosine, which indicates a possible DNA photo-oxidation mechanism for the PS. Triplet-triplet energy-transfer reactions from TH9402 to molecular oxygen, producing reactive singlet oxygen, and to the deprotonated guanosine, producing ³GMP^2-, are also thermodynamically favorable, with ΔG = -2.0 and -24.0 kcal//mol, respectively. However, the energy transfer to the monoprotonated guanosine is not favorable, (ΔG = 36.1), suggesting that in the DNA double-strand environment, this energy-transfer process may not be observed. The results show that the PS can act through electron transfer and triplet-triplet energy-transfer reactions involved in mechanism types I and II in photodynamic therapy. Interactions of TH9402 with the d(AGACGTCT)₂ octanucleotide revealed that the PS can intercalate between the d(GpC)-d(CpG) base pairs in three different orientations and, upon intercalation, the π → π* transition of the PS shows a bathochromic shift up to 90 nm and up to 60% decrease in intensity. Interactions through groove binding showed a smaller bathochromic shift of 52.2 nm and a 56% decrease in intensity of the main transition band.

Molecular Tailoring Approach for Estimating Individual Intermolecular Interaction Energies in Benzene Clusters

There is no general method available for the estimation of individual intermolecular interaction energies in weakly bound molecular clusters, and such studies are limited only to the dimer. Recently, we proposed a molecular tailoring approach-based method for the estimation of individual O-H···O hydrogen bond energies in water clusters. In the present work, we extend the applicability of this method for estimating the individual intermolecular interaction energies in benzene clusters, which are expected to be small. The basis set superposition error (BSSE)-corrected individual intermolecular interaction energies in linear (LN) benzene clusters, LN-(Bz)_n n = 3-7, were calculated to be in the range from -1.75 to -2.33 kcal/mol with the cooperativity contribution falling between 0.05 and 0.20 kcal/mol, calculated at the MP2.5/aug-cc-pVDZ level of theory. In the case of non-linear (NLN) benzene clusters, NLN-(Bz)_n n = 3-5, the BSSE-corrected individual intermolecular interaction energies exhibit a wider range from -1.16 to -2.55 kcal/mol with cooperativity contribution in the range from 0.02 to -0.61 kcal/mol. The accuracy of these estimated values was validated by adding the sum of interaction energies to the sum of monomer energies. These estimated molecular energies of clusters were compared with their actual calculated values. The small difference (<0.3 kcal/mol) in these two values suggests that our estimated individual intermolecular interaction energies in benzene clusters are quite reliable.

Implicit water model within the Zimm-Bragg approach to analyze experimental data for heat and cold denaturation of proteins

Studies of biopolymer conformations essentially rely on theoretical models that are routinely used to process and analyze experimental data. While modern experiments allow study of single molecules in vivo, corresponding theories date back to the early 1950s and require an essential update to include the recent significant progress in the description of water. The Hamiltonian formulation of the Zimm-Bragg model we propose includes a simplified, yet explicit model of water-polypeptide interactions that transforms into the equivalent implicit description after performing the summation of solvent degrees of freedom in the partition function. Here we show that our model fits very well to the circular dichroism experimental data for both heat and cold denaturation and provides the energies of inter- and intra-molecular H-bonds, unavailable with other processing methods. The revealed delicate balance between these energies determines the conditions for the existence of cold denaturation and thus clarifies its absence in some proteins.

Supramolecular structures of NiII and CuII with the sterically demanding Schiff base dyes driven by cooperative action of preagostic and other non-covalent interactions

This work reports on synthesis and extensive experimental and theoretical investigations on photophysical, structural and thermal properties of the NiII and CuII discrete mononuclear homoleptic complexes [Ni(L I,II)2] and [Cu(L I,II)2] fabricated from the Schiff base dyes o-HOC6H4-CH=N-cyclo-C6H11 (HL I) and o-HOC10H6-CH=N-cyclo-C6H11 (HL II), containing the sterically crowding cyclo-hexyl units. The six-membered metallocycles adopt a clearly defined envelope conformation in [Ni(L II)2], while they are much more planar in the structures of [Ni(L I)2] and [Cu(L I,II)2]. It has been demonstrated by in-depth bonding analyses based on the ETS-NOCV and Interacting Quantum Atoms energy-decomposition schemes that application of the bulky substituents, containing several C-H groups, has led to the formation of a set of classical and unintuitive intra- and inter-molecular interactions. All together they are responsible for the high stability of [Ni(L I,II)2] and [Cu(L I,II)2]. More specifically, London dispersion dominated intramolecular C-H⋯O, C-H⋯N and C-H⋯H-C hydrogen bonds are recognized and, importantly, the attractive, chiefly the Coulomb driven, preagostic (not repulsive anagostic) C-H⋯Ni/Cu interactions have been discovered despite their relatively long distances (∼2.8-3.1 Å). All the complexes are further stabilized by the extremely efficient intermolecular C-H⋯π(benzene) and C-H⋯π(chelate) interactions, where both the charge-delocalization and London dispersion constituents appear to be crucial for the crystal packing of the obtained complexes. All the complexes were found to be photoluminescent in CH2Cl2, with [Cu(L II)2] exhibiting the most pronounced emission - the time-dependent density-functional-theory computations revealed that it is mostly caused by metal-to-ligand charge-transfer transitions.

First-principles calculations of hybrid inorganic-organic interfaces: from state-of-the-art to best practice

The computational characterization of inorganic-organic hybrid interfaces is arguably one of the technically most challenging applications of density functional theory. Due to the fundamentally different electronic properties of the inorganic and the organic components of a hybrid interface, the proper choice of the electronic structure method, of the algorithms to solve these methods, and of the parameters that enter these algorithms is highly non-trivial. In fact, computational choices that work well for one of the components often perform poorly for the other. As a consequence, default settings for one materials class are typically inadequate for the hybrid system, which makes calculations employing such settings inefficient and sometimes even prone to erroneous results. To address this issue, we discuss how to choose appropriate atomistic representations for the system under investigation, we highlight the role of the exchange-correlation functional and the van der Waals correction employed in the calculation and we provide tips and tricks how to efficiently converge the self-consistent field cycle and to obtain accurate geometries. We particularly focus on potentially unexpected pitfalls and the errors they incur. As a summary, we provide a list of best practice rules for interface simulations that should especially serve as a useful starting point for less experienced users and newcomers to the field.

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.