Importance of London dispersion effects for the packing of molecular crystals: a case study for intramolecular stacking in a bis-thiophene derivative

Modeling Intermolecular Interactions with Exchange-Hole Dipole Moment Dispersion Corrections to Neural Network Potentials

Neural network potentials (NNPs) are an innovative approach for calculating the potential energy and forces of a chemical system. In principle, these methods are capable of modeling large systems with an accuracy approaching that of a high-level ab initio calculation, but with a much smaller computational cost. Due to their training to density-functional theory (DFT) data and neglect of long-range interactions, some classes of NNPs require an additional term to include London dispersion physics. In this Perspective, we discuss the requirements for a dispersion model for use with an NNP, focusing on the MLXDM (Machine Learned eXchange-Hole Dipole Moment) model developed by our groups. This model is based on the DFT-based XDM dispersion correction, which calculates interatomic dispersion coefficients in terms of atomic moments and polarizabilities, both of which can be approximated effectively using neural networks.

Comparison of Density-Functional Theory Dispersion Corrections for the DES15K Database

While density-functional theory (DFT) remains one of the most widely used tools in computational chemistry, most functionals fail to properly account for the effects of London dispersion. Hence, there are many popular post-self-consistent methods to add a dispersion correction to the DFT energy. Until now, the most popular methods have never been compared on equal footing due to not being implemented in the same electronic structure packages. In this work, we performed a large-scale benchmarking study, directly comparing the accuracy of the exchange-hole dipole moment (XDM), D3BJ, D4, TS, MBD, and MBD-NL dispersion models when applied to the recent DES15K database of nearly 15,000 molecular complexes at both expanded and compressed geometries. Our study showed similarly good performance for all dispersion methods (except TS) when applied to neutral complexes. However, they all performed worse for ionic complexes, particularly those involving dications of alkaline earth metals, due to systematic overbinding by the base PBE0 density functional. Investigation of the largest outliers also revealed that only the MBD and MBD-NL methods demonstrate surprising errors for complexes involving alkali metal cations at compressed geometries where they tended to significantly overbind. As we would expect minimal dispersion binding for such complexes, we further investigated the origins of these errors for the potential energy curve of a model cation-π complex. Overall, there is little choice between the XDM, D3BJ, D4, MBD, and MBD-NL dispersion methods for most systems. However, the MBD-based methods are not recommended for complexes involving organic species and alkali or alkaline earth metal cations, for example when modeling Li+ intercalation into graphite.

Tandem High-Pressure Crystallography-Optical Spectroscopy Unpacks Noncovalent Interactions of Piezochromic Fluorescent Molecular Rotors

To develop luminescent molecular materials with predictable and stimuli-responsive emission, it is necessary to correlate changes in their geometries, packing structures, and noncovalent interactions with the associated changes in their optical properties. Here, we demonstrate that high-pressure single-crystal X-ray diffraction can be combined with high-pressure UV-visible absorption and fluorescence emission spectroscopies to elucidate how subtle changes in structure influence optical outputs. A piezochromic aggregation-induced emitter, sym-heptaphenylcycloheptatriene (Ph7C7H), displays bathochromic shifts in its absorption and emission spectra at high pressure. Parallel X-ray measurements identify the pressure-induced changes in specific phenyl-phenyl interactions responsible for the piezochromism. Pairs of phenyl rings from neighboring molecules approach the geometry of a stable benzene dimer, while conformational changes alter intramolecular phenyl-phenyl interactions correlated with a relaxed excited state. This tandem crystallographic and spectroscopic analysis provides insights into how subtle structural changes relate to the photophysical properties of Ph7C7H and could be applied to a library of similar compounds to provide general structure-property relationships in fluorescent organic molecules with rotor-like geometries.

Identification of biomolecule-based electronic materials from a first-principles study of aliphatic amino acids

Biomolecule-based electronic materials can enable health innovations by virtue of their intrinsic bioactivity and physical properties. However, the ultra-wide bandgap and limited piezoelectric properties of most biomaterials prevent them from reaching their full potential. Herein, the electronic structures and electromechanical properties of aliphatic amino acid crystals are investigated based on density functional theory. L-Met is found to be a wide bandgap p-type semiconductor, and the much-reduced bandgap of 2.88 eV is ascribed to the sulphur atoms in L-Met. L-Leu has a shear piezoelectric voltage constant of 2.706 V mN^-1 that is over an order of magnitude higher than that of lead zirconate titanate, and good toughness and ductility are also revealed in L-Leu from mechanical property investigations. This study illustrates a computational approach to find smart and multifunctional biomaterials and inspire their growth and applications.

An artificial sodium-selective subnanochannel

Single-ion selectivity with high precision has long been pursued for fundamental bioinspired engineering and applications such as in ion separation and energy conversion. However, it remains a challenge to develop artificial ion channels to achieve single-ion selectivity comparable to their biological analogs, especially for high Na⁺/K⁺ selectivity. Here, we report an artificial sodium channel by subnanoconfinement of 4'-aminobenzo-15-crown-5 ethers (15C5s) into ~6-Å-sized metal-organic framework subnanochannel (MOFSNC). The resulting 15C5-MOFSNC shows an unprecedented Na⁺/K⁺ selectivity of tens to 10² and Na⁺/Li⁺ selectivity of 10³ under multicomponent permeation conditions, comparable to biological sodium channels. A co-ion-responsive single-file transport mechanism in 15C-MOFSNC is proposed for the preferential transport of Na⁺ over K⁺ due to the synergetic effects of size exclusion, charge selectivity, local hydrophobicity, and preferential binding with functional groups. This study provides an alternative strategy for developing potential single-ion selective channels and membranes for many applications.

Syntheses, Structural Characterization, and Cytotoxicity Assessment of Novel Mn(II) and Zn(II) Complexes of Aroyl-Hydrazone Schiff Base Ligand

This work describes the syntheses, structural characterization, and biological profile of Mn(II)- and Zn(II)-based complexes 1 and 2 derived from the aroyl-hydrazone Schiff base ligand (L1). The synthesized compounds were thoroughly characterized by elemental analysis, Fourier transform infrared spectroscopy (FTIR), UV-vis, electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), and single-crystal X-ray diffraction (s-XRD). Density functional theory (DFT) studies of complexes 1 and 2 were performed to ascertain the structural and electronic properties. Hirshfeld surface analysis was used to investigate different intermolecular interactions that define the stability of crystal lattice structures. To ascertain the therapeutic potential of complexes 1 and 2, in vitro interaction studies were carried out with ct-DNA and bovine serum albumin (BSA) using analytical and multispectroscopic techniques, and the results showed more avid binding of complex 2 than complex 1 and L1. The antioxidant potential of complexes 1 and 2 was examined against the 2,2-diphenyl picrylhydrazyl (DPPH) free radical, which revealed better antioxidant ability of the Mn(II) complex. Moreover, the antibacterial activity of synthesized complexes 1 and 2 was tested against Gram-positive and Gram-negative bacteria in which complex 2 demonstrated more effective bactericidal activity than L1 and complex 1 toward Gram-positive bacteria. Furthermore, the in vitro cytotoxicity assessment of L1 and complexes 1 and 2 was carried out against MDA-MB-231 (triple negative breast cancer) and A549 (lung) cancer cell lines. The cytotoxic results revealed that the polymeric Zn(II) complex exhibited better and selective cytotoxicity against the A549 cancer cell line as was evidenced by its low IC₅₀ value.

On the use of a volume constraint to account for thermal expansion effects on the low-frequency vibrations of molecular crystals

A volume-constraint method is presented as a means to capture the influence of thermal expansion on the low-frequency vibrations in molecular crystals. In particular, the room-temperature terahertz absorption spectra of L-tartaric acid, α-lactose monohydrate, and α-para-aminobenzoic acid (PABA) have been simulated using dispersion-corrected, solid-state density functional theory (DFT-D). By comparing the normal modes obtained with a unit cell optimised without constraints to those obtained with a unit cell optimised while constrained to keep its experimental volume, wholesale improvements to the resultant spectrum is achieved when using the constrained geometry by inhibiting cell contraction. These improvements are demonstrated over a range of popular density functionals and basis sets up to triple-zeta complexity. A correlation method is then presented as a means to quantitatively compare the vibrational pattern of normal modes obtained from both unit cells. This analysis reveals that thermal expansion can effect the character and relative frequency of normal modes, with the choice of geometry ultimately affecting the assignment of the experimental absorptions. The sensibility of using the experimental volume as an approximation is then discussed, where it is speculated that large basis sets or hybrid functionals are necessary to ensure that the thermal expansion effect is not overestimated. The low-frequency absorption spectrum of PABA is then fully characterised using the PBE-D3BJ/6-311G(2d,2p) method.

Mixed Ni(II) and Co(II) complexes of nalidixic acid drug: Synthesis, characterization, DNA/BSA binding profile and in vitro cytotoxic evaluation against MDA-MB-231 and HepG2 cancer cell lines

In this work, herein we report the synthesis, structural characterization and in vitro cytotoxic evaluation of two mixed Co(II)/Ni(II)-nalidixic acid-bipyridyl complexes (1 and 2). The structural analysis of metal complexes 1 and 2 was carried out by analytical and multispectroscopic techniques (FT-IR, UV-vis, EPR, sXRD). The crystallographic details of complexes 1 and 2 revealed a monoclinic crystal system with P21/c space group. DFT studies of complexes were performed to get electronic structure and localization of HOMO and LUMO electron densities. Hirshfeld surface analysis of metal complexes 1 and 2 was employed to understand the various intermolecular interactions (C-H···O, N-H···H and O-H···O) that define the stability of crystal lattice structures. The comparative interaction studies of complex 1 and complex 2 with DNA/BSA were performed by diverse multispectroscopic and analytical techniques to evaluate their chemotherapeutic potential. The magnitude of the DNA binding propensity and binding mode was verified by calculating K_b, K and K_sv values. Higher binding affinity was observed in case of complex 2via intercalative mode. Furthermore, the cytotoxic assessment of complexes 1 and 2 was examined against MDA-MB-231 (triple negative human breast cancer cell line) and HepG2 (liver carcinoma cell line) employing MTT assay which revealed remarkably effecient and specific cytotoxic activity of complex 2.

Capture of Gaseous Iodine in Isoreticular Zirconium-Based UiO-n Metal-Organic Frameworks: Influence of Amino Functionalization, DFT Calculations, Raman and EPR Spectroscopic Investigation

A series of Zr-based UiO-n MOF materials (n=66, 67, 68) have been studied for iodine capture. Gaseous iodine adsorption was collected kinetically from a home-made set-up allowing the continuous measurement of iodine content trapped within UiO-n compounds, with organic functionalities (-H, -CH₃ , -Cl, -Br, -(OH)₂ , -NO₂ , -NH₂ , (-NH₂ )₂ , -CH₂ NH₂ ) by in-situ UV-Vis spectroscopy. This study emphasizes the role of the amino groups attached to the aromatic rings of the ligands connecting the {Zr₆ O₄ (OH)₄ } brick. In particular, the preferential interaction of iodine with lone-pair groups, such as amino functions, has been experimentally observed and is also based on DFT calculations. Indeed, higher iodine contents were systematically measured for amino-functionalized UiO-66 or UiO-67, compared to the pristine material (up to 1211 mg/g for UiO-67-(NH₂ )₂ ). However, DFT calculations revealed the highest computed interaction energies for alkylamine groups (-CH₂ NH₂ ) in UiO-67 (-128.5 kJ/mol for the octahedral cavity), and pointed out the influence of this specific functionality compared with that of an aromatic amine. The encapsulation of iodine within the pore system of UiO-n materials and their amino-derivatives has been analyzed by UV-Vis and Raman spectroscopy. We showed that a systematic conversion of molecular iodine (I₂ ) species into anionic I^- ones, stabilized as I^- ⋅⋅⋅I₂ or I₃ ^- complexes within the MOF cavities, occurs when I₂ @UiO-n samples are left in ambient light.

Designing a Solution-Stable Distannene: The Decisive Role of London Dispersion Effects in the Structure and Properties of {Sn(C6H2-2,4,6-Cy3)2}2 (Cy = Cyclohexyl)

The reaction of 1 equiv of the dimeric lithium salt of a new London dispersion effect donor ligand {Li(C₆H₂-2,4,6-Cy₃)·OEt₂}₂ (Cy = cyclohexyl) with SnCl₂ afforded the distannene {Sn(C₆H₂-2,4,6-Cy₃)₂}₂ (1). The distannene remains dimeric in solution, as indicated by its room-temperature ¹¹⁹Sn NMR signal (δ = 361.3 ppm) and its electronic spectrum, which is invariant over the temperature range of -10 to 100 °C. The formation of the distannene, which has a short Sn-Sn distance of 2.7005(7) Å and greatly enhanced stability in solution compared to that of other distannenes, is due to increased interligand London dispersion (LD) attraction arising from multiple close approaches of ligand C-H moieties across the Sn-Sn bond. DFT-D4 calculations revealed a dispersion stabilization of dimer 1 of 38 kcal mol^-1 and a dimerization free energy of ΔG_dimer = -6 kcal mol^-1. In contrast, the reaction of 2 equiv of the similarly shaped but less bulky, less hydrogen-rich Li(C₆H₂-2,4,6-Ph₃)·(OEt₂)₂ with SnCl₂ yielded the monomeric stannylene Sn(C₆H₂-2,4,6-Ph₃)₂ (2), which is unstable in solution at ambient temperature.

Structure elucidation, in vitro binding studies and ROS-dependent anti-cancer activity of Cu(II) and Zn(II) phthaloylglycinate(phen) complexes against MDA-MB-231 cells

New mononuclear Cu(II) and Zn(II)-based complexes 1 [Cu(L)2(diimine)HOCH3] and 2 [Zn(L)2(diimine)] have been synthesized as anti-cancer chemotherapeutics targeted to tRNA. The structure elucidation of complexes 1 and 2 was carried out by spectroscopic and single X-ray diffraction studies. In vitro interaction studies of complexes 1 and 2 with ct-DNA/tRNA were performed by employing various biophysical techniques to evaluate and predict their interaction behavior and preferential selectivity at biomolecular therapeutic targets. The corroborative results of the interaction studies demonstrated that complexes 1 and 2 exhibited avid binding propensity via intercalative mode of binding toward ct-DNA/tRNA. Electrophoretic assay revealed that the complexes 1 and 2 were able to promote single- and double-strand cleavage of the plasmid DNA at low micromolar concentrations under physiological conditions in the absence of an additional oxidizing or reducing agent. RNA hydrolysis studies revealed that the complexes 1 and 2 could promote tRNA cleavage in a concentration and time-dependent manner. The cytotoxic potential of complexes 1 and 2 was evaluated against the MDA-MB-231 cell line, which showed that the complexes were able to inhibit the cell growth in a dose-dependent manner. The intracellular ROS production and mitochondrial superoxide anion assay revealed that the complexes 1 and 2 induce a dose-dependent activity, suggesting the involvement of ROS-mediated mitochondrial apoptotic pathway leading to cell death.

Novel Polymorph of Favipiravir-An Antiviral Medication

Various solid forms of pharmaceutically important compounds exhibit different physical properties and bioactivity; thus, knowledge of the structural landscape and prediction of spontaneous polymorph transformations for an active pharmaceutical ingredient is of practical value for the pharmaceutical industry. By recrystallization from ethyl acetate, a novel polymorph of 6-fluoro-3-hydroxypyrazine-2-carboxamide (trademark favipiravir, RNA polymerase inhibitor) was obtained and characterized using differential scanning calorimetry (DSC), infra-red spectroscopy and powder X-ray diffraction (XRD) analysis. The favipiravir molecule in two polymorphs realizes similar H-bonding motifs, but the overall H-bonded networks differ. Based on periodic density functional theory calculations, the novel tetragonal polymorph with two interpenetrated H-bonded networks is slightly less stable than the orthorhombic one with the zst topology of the underlying H-bonded net that is in accord with experimentally observed powder XRD patterns of slow conversion of the tetragonal phase to the orthorhombic one. However, topological analysis of net relations revealed that no transformations can be applied to convert H-bonded networks in the experimental unit cells, and DSC data indicate no solid-state reactions at heating.

Density functional theory predictions of the mechanical properties of crystalline materials

The DFT-predicted mechanical properties of crystalline materials are crucial knowledge for their screening, design, and exploitation.

Water soluble ionic Co(II), Cu(II) and Zn(II) diimine-glycinate complexes targeted to tRNA: structural description, in vitro comparative binding, cleavage and cytotoxic studies towards chemoresistant prostate cancer cells

Four new water soluble Co(ii), Cu(ii) and Zn(ii) ionic metal complexes (1-4) [Cu(diimine)(H2O)2(glycinate)]+[glycinate]-, [Co(diimine)(H2O)4]+[glycinate]- and [Zn(diimine) (H2O)4]+[glycinate]-, where diimine = 2,2'-bipyridine (1-3) and 1,10-phenanthroline (4) were synthesized and thoroughly characterized by spectroscopic and single X-ray crystallographic studies. Complex 1 possesses a triclinic crystal system with a penta-coordinated geometry whereas complexes 2-4 crystallized in an isostructural monoclinic system having distorted octahedral geometry. Density functional theory (DFT) studies for complexes 1-4 were performed to correlate their geometrical parameters and to calculate the energy of frontier molecular orbitals. The corroborative results of spectroscopic and voltammetric studies with ct-DNA and tRNA revealed that the complexes bind noncovalently via an electrostatic mode of binding with specificity for tRNA as compared to ct-DNA. Gel electrophoresis experiments revealed that all the complexes unwind the plasmid pBR322 DNA at low micromolar concentrations (2-9 μM) following an oxidative mechanism for Cu(ii) and Co(ii) complexes (1, 2 and 4) whereas the Zn(ii) complex (3) mediates DNA cleavage by the hydrolytic pathway. The tRNA cleavage showed concentration and time dependent activity of the complexes to promote RNA hydrolysis. Furthermore, the BSA binding ability of complexes 1-4 was monitored, which revealed that the complexes could quench the intrinsic fluorescence in a static manner. Complexes 1-4 were found to be non-toxic towards normal prostate epithelial cells, PNT2, but were potent against chemoresistant metastatic prostate cancer cells, Du145, with GI50 values ranging from 12.75-37 μM. Complexes 1 and 2 also showed cytotoxic activity against cancer stem cells having GI50 values of 14.70 and 14.90 μM, respectively. Molecular docking studies were performed with DNA and tRNA which further validated the spectroscopic analysis demonstrating the higher binding affinity of the complexes towards tRNA.

What is “many-body” dispersion and should I worry about it?

“Many-body” dispersion can refer to two distinct phenomena, here termed electronic and atomic many-body effects, both of which cause the dispersion energy to be non-additive.

Soft experimental constraints for soft interactions: a spectroscopic benchmark data set for weak and strong hydrogen bonds

An experimental benchmark data base on rotational constants, vibrational properties and energy differences for weakly and more strongly hydrogen-bonded complexes and their constituents from the spectroscopic literature is assembled. It is characterized in detail and finally contracted to a more compact, discriminatory set (ENCH-51, for Experimental Non-Covalent Harmonic with 51 entries). The meeting points between theory and experiment consist of equilibrium rotational constants and harmonic frequencies and energies, which are back-corrected from experimental observables and are very easily accessible by quantum chemical calculations. The relative performance of B3LYP-D3, PBE0-D3 and M06-2X density functional theory predictions with a quadruple-zeta basis set is used to illustrate systematic errors, error compensation and selective performance for structural, vibrational and energetical observables. The current focus is on perspectives and different benchmarking methodologies, rather than on a specific theoretical method or a specific class of compounds. Extension of the data base in chemical, observable and quantum chemical method space is encouraged.

DL-3-Aminoisobutyric acid: vibrational, NBO and AIM analysis of N-H⋯O bonded-zwitterionic dimer model

A zwitterionic dimer model constructed of inter-molecular -N-H⋯O bonding has been proposed for the solid sample of DL-3-Aminoisobutyric acid consistent with IR absorption and Raman spectral features measured in the 3500-400/50 cm^-1. This zwitterionic dimer model in water as solvent has been computed at B3LYP/6-311++G(d,p) and B3LYP-D3/6-311++G(d,p) levels including Grimme's dispersion correction associated with the -N-H⋯O interaction and SCRF-SMD method. Of the several possible monomer and dimer conformational structures, the most stable dimer constructed of two zwitterion monomer units has produced vibrational modes due to the -NH₃ ⁺ cation and -CO₂‾ anion involved in the -N-H⋯O bonding in fair agreement with the observed broad but composite IR modal features near the 3500-2000 cm^-1. Except for the frequency of asymmetric stretching mode of the -NH₃ ⁺ cation, its symmetric and bending modes agree with the observed values. As for the -CO₂‾ anion, the frequencies of all of its modes are in good agreement with the experiment. Natural bond orbital (NBO), molecular electrostatic potential (MEP), atoms-in-molecules (AIM) and non-covalent interaction (NCI) analyses have been used to understand electronic characterization of the -N-H⋯O bonding.

Molecular Structures Polymorphism the Role of F…F Interactions in Crystal Packing of Fluorinated Tosylates

The peculiarities of interatomic interactions formed by fluorine atoms were studied in four tosylate derivatives p-CH3C6H4OSO2CH2CF2CF3 and p-CH3C6H4OSO2CH2(CF2)nCHF2 (n = 1, 5, 7) using X-ray diffraction and quantum chemical calculations. Compounds p-CH3C6H4OSO2CH2(CF2)nCHF2 (n = 1, 5) were crystallized in several polymorph modifications. Analysis of intermolecular bonding was carried out using QTAIM approach and energy partitioning. All compounds are characterized by crystal packing of similar type and the contribution of intermolecular interactions formed by fluorine atoms to lattice energy is raised along with the increase of their amount. The energy of intra- and intermolecular F…F interactions is varied in range 0.5–13.0 kJ/mol. Total contribution of F…F interactions to lattice energy does not exceed 40%. Crystal structures of studied compounds are stabilized mainly by C-H…O and C-H…F weak hydrogen bonds. The analysis of intermolecular interactions and lattice energies in polymorphs of p-CH3C6H4OSO2CH2(CF2)nCHF2 (n = 1, 5) has shown that most stabilized are characterized by the least contribution of F…F interactions.

Theory and practice of modeling van der Waals interactions in electronic-structure calculations

The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.

Unique Surface Enhanced Raman Scattering Substrate for the Study of Arsenic Speciation and Detection

In this study, a three-dimensional surface enhanced Raman scattering (SERS) substrate comprised of silver coated gold nanorods (Ag/AuNRs) decorated on electrospun polycaprolactone (PCL) fibers has been applied,  for the first time, to quantitative analytical measurements on various arsenic species: p-arsanilic acid ( pAsA), roxarsone (Rox), and arsenate (As^V), with a demonstrated sensitivity below 5 ppb. As^V detection in a solution of common salt ions has been demonstrated, showing the tolerance of the substrate to more complex environments. pAsA adsorption behavior on the substrate surface has been investigated in detail using these unique SERS substrates. Calculations based on density functional theory (DFT) support the spectral observation for pAsA. This substrate also has been shown to serve as a platform for in situ studies of arsenic desorption and reduction. This SERS substrate is potentially an excellent environmental sensor for both fundamental studies and practical applications.

Expanding the Range of Force Fields Available for ONIOM Calculations: The SICTWO Interface

The ONIOM scheme is one of the most popular QM/MM approaches, but its extended application has been so far hindered by the limited availability of force fields in most practical implementations. This paper describes a simple software code to overcome this limitation, and its application to three representative chemical problems. The "Shell Interface for Combining Tinker With ONIOM" (SICTWO) program gives access to all force fields available in the Tinker molecular mechanics program from a Gaussian09 or Gaussian16 calculation. The first application presented is the geometry optimization of dithienobicyclo-[4.4.1]-undeca-3,8-diene-11-one ethylene glycol ketal. A variety of force fields were tested for the MM part, and the molecular structure using the ONIOM(B3LYP:OPLS-AA) method shows good agreement with the experimental results. The second problem is related with enantioselectivity of the vanadium-catalyzed asymmetric oxidation of 1,2-bis( tert-butyl)-disulfide. ONIOM(B3LYP:MM3) results, obtained with SICTWO, are consistent with those of a previous IMOMM(B3LYP:MM3) study. In the third application, we have combined AMOEBA09 polarizable force field with ONIOM to study a water molecule binding on water ice. Calculated ONIOM(ωB97X-D:AMOEBA09) binding energies are consistent with the ωB97X-D results. These studies show SICTWO as an easy-to-use tool that can further expand the use of the multiscale ONIOM method within computational chemistry and computational biology.

The adsorption of alcohols on strained Pt3Ni(111) substrates: a density functional investigation within the D3 van der Waals correction

In this theoretical study, we address the effect of strain and alloying on the adsorption of methanol, ethanol and glycerol on Pt₃Ni(111) surfaces.

A DFT study of RuO4 interactions with porous materials: metal–organic frameworks (MOFs) and zeolites

The potential use of zeolite and MOF materials for the capture of RuO₄ has been investigated for the first time. A hydrated form of HKUST-1 could be a promising sorbent due to its ability to form multiple hydrogen bonds.

The role of the alkali and chalcogen atoms on the stability of the layered chalcogenide [Formula: see text] (A = alkali-metal; M = metal-cations; Q = chalcogen) compounds: a density functional theory investigation within van der Waals corrections

There is a great interest to design two-dimensional (2D) chalcogenide materials, however, our atomistic understanding of the major physical parameters that drive the formation of 2D or three-dimensional (3D) chalcogenides is far from satisfactory, in particular, for complex quaternary systems. To address this problem, we selected a set of quaternary 2D and 3D chalcogenide compounds, namely, [Formula: see text] (A  =  Li, K, Cs; Q  =  S, Se, Te), which were investigated by density functional theory calculations within van der Waals (vdW) corrections. Employing experimental crystal structures and well designed crystal modifications, we found that the average atomic radius of the alkali-metal, A, and chalcogen, Q, species play a crucial role in the stability of the 2D structures. For example, the 2D structures are energetically favored for smaller [Formula: see text] and larger [Formula: see text] average atomic radius, while 3D structures are favored at intermediate average atomic radius. Those results are explained in terms of strain minimization and Coulomb repulsion of the anionic species in the structure framework. Furthermore, the equilibrium lattice parameters are in excellent agreement with experimental results. Thus, the present insights can help in the design of stable quartenary 2D chalcogenide compounds.

Glycerol adsorption on a defected Pt6/Pt(100) substrate: a density functional theory investigation within the D3 van der Waals correction

We report an ab initio investigation to obtain an improved atomistic understanding of the adsorption properties of glycerol on a defected Pt₆/Pt(100) substrate.

Molecular structure and vibrations of NTCDA monolayers on Ag(111) from density-functional theory and infrared absorption spectroscopy

The structure and vibrational properties of the metal-organic interface of 1,4,5,8-naphthalene-tetracarboxylic dianhydride (NTCDA) on Ag(111) were analysed using Fourier-transform infrared absorption spectroscopy in conjunction with density functional theory calculations including dispersion forces (PBE-D3). Mode assignments and polarizations as well as molecular distortions were determined for four adsorption geometries of NTCDA on top and bridge sites aligned either parallel or perpendicular to the Ag rows and compared to accurate calculations of the free molecule. This enables an in-depth understanding of surface effects on the computed and experimental vibrational spectra of the adsorbed NTCDA molecule. The molecule-substrate interaction comprises two major and equally important contributions: non-directional van der Waals forces between molecule and surface, and covalent bonding of the acyl oxygen atoms with underlying Ag atoms, which is quantified by charge-transfer analysis. Furthermore, adsorption energy calculations showed that the molecular axis of flat-lying NTCDA is oriented preferably in parallel to the Ag rows. The molecule is subject to particular distortions from the planar gas phase structure with covalent bonding leading to downward bending of the acyl oxygen atoms and Pauli repulsion to upward bending of the carbon core. In parallel, strong buckling of the silver surface was identified. As found in previous studies, the lowest unoccupied molecular orbital (LUMO) of the molecule slips below the Fermi level and becomes partially populated upon adsorption. Excitation of totally symmetric vibrational modes then leads to substantial interfacial dynamical charge transfer, which is convincingly reproduced in the calculated IR spectra.

Aliphatic–aromatic stacking interactions in cyclohexane–benzene are stronger than aromatic–aromatic interaction in the benzene dimer

Stacking interactions between cyclohexane and benzene were studied in crystal structures from the Cambridge Structural Database and by ab initio calculations.

Adsorption of water and ethanol on noble and transition-metal substrates: a density functional investigation within van der Waals corrections

We report an extensive density functional theory investigation of water and ethanol adsorption on several Cu-, Pt-, and Au-based substrates including substrates with low-coordinated sites due to intermixing of Pt–Cu and Pt–Au in the topmost surface.

Consistent structures and interactions by density functional theory with small atomic orbital basis sets

A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of "low-cost" electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT methods and reach that of triple-zeta AO basis set second-order perturbation theory (MP2/TZ) level at a tiny fraction of computational effort. Periodic calculations conducted for molecular crystals to test structures (including cell volumes) and sublimation enthalpies indicate very good accuracy competitive to computationally more involved plane-wave based calculations. PBEh-3c can be applied routinely to several hundreds of atoms on a single processor and it is suggested as a robust "high-speed" computational tool in theoretical chemistry and physics.

Conformational and vibrational analyses of meta-tyrosine: An experimental and theoretical study

M-tyrosine is one kind of positional isomer of tyrosine which is widely applied in agrichemical, medicinal chemistry, and food science. However, the structural and vibrational features of m-tyrosine have not been reported or systematically investigated. In this work, potential energy surface (PES) calculations were used for searching and determining the stable zwitterionic conformers of m-tyrosine, and the Raman spectra of m-tyrosine and deuterated m-tyrosine were measured and interpreted based on theoretical computation. For the spectral simulation, several DFT-based quantum chemistry (QC) methods were employed, and the M06-2X functional with SMD solvent model was found to be best in reproducing the Raman spectra and geometrical property. As such, this study has not only presented a detailed study of m-tyrosine's vibrational property which is lack in the literature, but also may shed some light on the optimal choice of QC methods for calculation of conformations and vibrational properties of zwitterionic amino acids.

Conformational studies of poly(vinylidene fluoride), poly(trifluoroethylene) and poly(vinylidene fluoride-co-trifluoroethylene) using density functional theory

Different conformations of systems consisting of poly(vinylidene fluoride) (PVDF), poly(trifluoroethylene) (PTrFE) and poly(vinylidene fluoride-co-trifluoroethylene) (P(VDF-TrFE)) were investigated using density functional theory with dispersion correction. It was found that the trans-gauche-trans-gauche' (TGTG') conformation of a single PVDF chain is the lowest energy conformer. Crystals of PVDF were modelled using between two to five chains with up to 12 repeat units in each chain and, in agreement with experiment, structures comprised partly or completely of chains with the TGTG' conformation are more stable than structures built up from chains with all-trans (TTTT) conformation. This indicates that an all-trans segment or chain will not induce the growth of a larger crystal with the same chain conformations. In contrast, the energetically preferred structure of PTrFE chains is an all-trans (TTTT) conformation, and the results indicate that copolymerization of vinylidene fluoride with trifluoroethylene can facilitate the formation of the all-trans PVDF conformations. This is probably due to increased intramolecular repulsion between the fluorine atoms and an increased intermolecular attraction in the crystal structure.

The thermochemistry of london dispersion-driven transition metal reactions: getting the 'right answer for the right reason'

Reliable thermochemical measurements and theoretical predictions for reactions involving large transition metal complexes in which long-range intramolecular London dispersion interactions contribute significantly to their stabilization are still a challenge, particularly for reactions in solution. As an illustrative and chemically important example, two reactions are investigated where a large dipalladium complex is quenched by bulky phosphane ligands (triphenylphosphane and tricyclohexylphosphane). Reaction enthalpies and Gibbs free energies were measured by isotherm titration calorimetry (ITC) and theoretically 'back-corrected' to yield 0 K gas-phase reaction energies (ΔE). It is shown that the Gibbs free solvation energy calculated with continuum models represents the largest source of error in theoretical thermochemistry protocols. The ('back-corrected') experimental reaction energies were used to benchmark (dispersion-corrected) density functional and wave function theory methods. Particularly, we investigated whether the atom-pairwise D3 dispersion correction is also accurate for transition metal chemistry, and how accurately recently developed local coupled-cluster methods describe the important long-range electron correlation contributions. Both, modern dispersion-corrected density functions (e.g., PW6B95-D3(BJ) or B3LYP-NL), as well as the now possible DLPNO-CCSD(T) calculations, are within the 'experimental' gas phase reference value. The remaining uncertainties of 2-3 kcal mol(-1) can be essentially attributed to the solvation models. Hence, the future for accurate theoretical thermochemistry of large transition metal reactions in solution is very promising.

Polymorphism in two biologically active dihydropyrimidinium hydrochloride derivatives: quantitative inputs towards the energetics associated with crystal packing

A new polymorph belonging to the tetrahydropyrimidinium class of compounds, namely 6-(4-chlorophenyl)-5-(methoxycarbonyl)-4-methyl-2-(3-(trifluoromethylthio)phenylamino)-3,6-dihydropyrimidin-1-ium chloride, and a hydrate of 2-(3-bromophenylamino)-6-(4-chlorophenyl)-5-(methoxycarbonyl)-4-methyl-3,6-dihydropyrimidin-1-ium chloride, have been isolated and characterized using single-crystal X-ray diffraction (XRD). A detailed comprehensive analysis of the crystal packing in terms of the associated intermolecular interactions and a quantification of their interaction energies have been performed for both forms of the two different organic salts (A and B) using X-ray crystallography and computational methods such as density functional theory (DFT) quantum mechanical calculations, PIXEL lattice-energy calculations (with decomposition of total lattice energy into the Coulombic, polarization, dispersion and repulsion contribution), the calculation of the Madelung constant (the EUGEN method), Hirshfeld and two-dimensional fingerprint plots. The presence of ionic [N-H](+)···Cl(-) and [C-H](+)···Cl(-) hydrogen bonds mainly stabilizes the crystal packing in both forms A and B, while in the case of B·H2O [N-H](+)···O(water) and O(water)-H···Cl(-) hydrogen bonds along with [N-H](+)···Cl(-) and [C-H](+)···Cl(-) provide stability to the crystal packing. The lattice-energy calculations from both PIXEL and EUGEN methods revealed that in the case of A, form (I) (monoclinic) is more stable whereas for B it is the anhydrous form that is more stable. The analysis of the `Madelung mode' of crystal packing of two forms of A and B and its hydrates suggest that differences exist in the position of the charged ions/atoms in the organic solid state. The R/E (distance-energy) plots for all the crystal structures show that the molecular pairs in their crystal packing are connected with either highly stabilizing (due to the presence of organic R(+) and Cl(-)) or highly destabilizing Coulombic contacts. The difference in crystal packing and associated intermolecular interactions between polymorphs (in the case of A) or the hydrates (in the case of B) have been clearly elucidated by the analysis of Hirshfeld surfaces and two-dimensional fingerprint plots. The relative contributions of the various interactions to the Hirshfeld surface for the cationic (dihydropyrimidinium) part and anionic (chloride ion) part for the two forms of A and B and its hydrate were observed to be different.