Accuracy of intermolecular interaction energies, particularly those of hetero-atom containing molecules obtained by DFT calculations with Grimme's D2, D3 and D3BJ dispersion corrections

Insights into the Adsorption Mechanism of Chlorpyrifos on Activated Carbon Derived from Prickly Pear Seeds Waste: An Experimental and DFT Modeling Study

The removal of chlorpyrifos (CPF) from water was achieved using activated carbon (AC) derived from prickly pear seeds (PPS) wastes, developed through chemical activation with phosphoric acid. Several physico-chemical characterization methods were employed. The determination of surface functions using the Boehm assay indicated that the processed AC predominantly possesses acidic functions. The results obtained from the Boehm assay were corroborated by the pH value of the point of zero charge (pHpzc), which was equal to 2.5. Specific area calculation by the BET (Brunauer Emmett Teller) method revealed a large specific area (SBET) of 1077.66 m2 g⁻1. Adsorption experiments of CPF on AC demonstrated that the pseudo-second order (PSO) model and the Freundlich model were the most suitable for kinetic and isothermal modeling, respectively. The maximum CPF adsorption capacity of the PPS AC was found to be approximately 35 mg g⁻1. A theoretical study employing the density functional theory (DFT) was conducted using the B3LYP/6-311G (d, p) method. The most reliable adsorption energy (Eads) and Gibbs free energy (ΔGads) values between CPF and the functional groups on the AC surface were calculated. Results indicated a strong interaction between the lactone group of AC and CPF (ΔGads = -7.15 kcal mol⁻1, ΔEads = -21.55 kcal mol⁻1) and the hydroxyl group (ΔGads = -6.61 kcal mol⁻1, ΔEads = -20.66 kcal mol⁻1). This study demonstrates that activated carbon possesses significant adsorption power, making it highly effective for depolluting water contaminated by pesticides. The application of the theoretical DFT method enhances the understanding of the adsorption phenomenon of CPF on AC.

Stability of n-alkanes and n-perfluoroalkanes against horizontal displacement on a graphite surface

The stability of adsorbed molecules on surfaces is fundamental and important for various applications, such as coating, lubrication, friction, and self-assembled structure formation. In this study, we investigated the structures and interaction energies (Eint) of propane, n-pentane, n-heptane, perfluoropropane, n-perfluoropentane, and n-perfluoroheptane adsorbed on the surface of C96H24 (a model surface of graphite). The changes in Eint (ΔEint = Eint - Eint(0)) associated with the horizontal displacement from the stable position were calculated using dispersion-corrected density functional theory (DFT; B3LYP-D3), where Eint(0) is the Eint at the stable position. The maximum value of ΔEint (ΔEint(max)) associated with the horizontal displacement increased as the chain length increased. The ΔEint(max) for the three n-alkanes were 1.10, 1.82, and 2.35 kcal mol-1, respectively. The values for n-perfluoroalkanes were 0.57, 0.83, and 1.04 kcal mol-1, respectively. The ΔEint(max) values for the n-alkanes were significantly larger than those for the corresponding n-perfluoroalkanes. The Eint(max) value per carbon atom of the n-alkanes (ca. 0.30 kcal mol-1) is approximately 2.5 times as large as that of n-perfluoroalkanes (ca. 0.12 kcal mol-1). The ΔEint associated with the horizontal displacement of propane and perfluoropropane on circumcoronene (C54H18) obtained by the B3LYP-D3 calculations are close to those obtained by the second order Møller-Plesset (MP2) and dispersion-corrected double hybrid DFT calculations, suggesting the sufficient accuracy of the ΔEint obtained by the B3LYP-D3. Thus, our quantitative analysis revealed the higher stability of n-alkanes against horizontal displacement on a graphite surface than that of n-perfluoroalkanes.

New routes towards azomethine ylide generation from prolines to synthesize diverse N-heterocycles: a DFT supported endo-selective mechanism

Azomethine ylides are generated using either organocatalysts or metal catalysts via a ballet of decarboxylative C-N coupling choreographed by prolines. These strategies enable diastereoselective [3 + 2] cycloaddition, C-C coupling, and ring annulation, providing sustainable routes. The synthesized pyrrolizines and other heterocycles have potential applications in the development of crucial biomolecules and pharmaceuticals. The endoselectivity of the azomethine ylide is realized and supported through DFT calculations.

A dual-cutoff machine-learned potential for condensed organic systems obtained via uncertainty-guided active learning

Machine-learned potentials (MLPs) trained on ab initio data combine the computational efficiency of classical interatomic potentials with the accuracy and generality of the first-principles method used in the creation of the respective training set. In this work, we implement and train a MLP to obtain an accurate description of the potential energy surface and property predictions for organic compounds, as both single molecules and in the condensed phase. We devise a dual descriptor, based on the atomic cluster expansion (ACE), that couples an information-rich short-range description with a coarser long-range description that captures weak intermolecular interactions. We employ uncertainty-guided active learning for the training set generation, creating a dataset that is comparatively small for the breadth of application and consists of alcohols, alkanes, and an adipate. Utilizing that MLP, we calculate densities of those systems of varying chain lengths as a function of temperature, obtaining a discrepancy of less than 4% compared with experiment. Vibrational frequencies calculated with the MLP have a root mean square error of less than 1 THz compared to DFT. The heat capacities of condensed systems are within 11% of experimental findings, which is strong evidence that the dual descriptor provides an accurate framework for the prediction of both short-range intramolecular and long-range intermolecular interactions.

Structural and Electronic Factors Controlling the Efficiency and Rate of Intersystem Crossing to the Triplet State in Thiophene Polycyclic Derivatives

Thiophene polycyclic derivatives are widely used in organic light-emitting diodes, photovoltaics, and medicinal chemistry applications. Understanding the electronic and structural factors controlling their intersystem crossing rates is paramount for these applications to be successful. This study investigates the photophysical, electronic structure, and excited state dynamics of 1,2-benzodiphenylene sulfide, benzo[b]naphtho[1,2-d]thiophene, and benzo[b]naphtho[2,3-d]thiophene in polar aprotic and non-polar solvents. Steady-state absorption and emission spectroscopy, femtosecond transient absorption spectroscopy, and DFT and TD-DFT calculations are employed. Low fluorescence quantum yields of 1.2 to 2.7 % are measured in acetonitrile and cyclohexene, evidencing that the primary relaxation pathways in these thiophene derivatives are nonradiative. Linear interpolation of internal coordinates calculations predict that an S-C bond elongation reaction coordinate facilitates the efficient intersystem crossing to the T1 state. Excitation of 1,2-benzodiphenylene sulfide and benzo[b]naphtho[1,2-d]thiophene at 350 nm or benzo[b]naphtho[2,3-d]thiophene at 365 nm, populates the lowest-energy 1ππ* state, which relaxes to the 1ππ* minimum in tens of picoseconds or intersystem crosses to the triplet manifold in ca. 500 ps to 1.1 ns depending on the position at which the benzene rings are added. Excitation at 266 nm does not affect the intersystem crossing rates. Laser photodegradation experiments demonstrate that the thiophene polycyclic derivatives are highly photostable.

Two pairs of new isobenzofuranone enantiomers from a soil-derived fungus Penicillium canescens DWS225

Two pairs of new isobenzofuranone derivative enantiomers, (±)-penicifurans E (1) and (±)-penicifurans F (2), together with four know compounds (3-6) were isolated from the solid fermentation of Penicillium canescens DWS225. The structures of these enantiomers were elucidated by extensive NMR spectroscopic data, and their absolute configurations were assigned by the experimental and calculated ECD data. The neuroprotective effects of all the isolates against oxygen-glucose deprivation/reperfusion injury in pheochromocytoma-12 cells (PC12) were investigated.

Ligand Characterization and DNA Intercalation of Ru(II) Polypyridyl Complexes: A Local Vibrational Mode Study

We investigated in this work ruthenium-ligand bonding across the RuN framework in 12 Ru(II) polypyridyl complexes in the gas phase and solution for both singlet and triplet states, in addition to their affinity for DNA binding through π-π stacking interactions with DNA nucleobases. As a tool to assess the intrinsic strength of the ruthenium-ligand bonds, we determined local vibrational force constants via our local vibrational mode analysis software. We introduced a novel local force constant that directly accounts for the intrinsic strength of the π-π stacking interaction between DNA and the intercalated Ru(II) complex. According to our findings, [Ru(phen)2(dppz)]2+ and [Ru(phen)2(11-CN-dppz)]2+ provide an intriguing trade-off between photoinduced complex excitation and the strength of the subsequent π-π stacking interaction with DNA. [Ru(phen)2(dppz)]2+ displays a small singlet-triplet splitting and a strong π-π stacking interaction in its singlet state, suggesting a favorable photoexcitation but potentially weaker interaction with DNA in the excited state. Conversely, [Ru(phen)2(11-CN-dppz)]2+ exhibits a larger singlet-triplet splitting and a stronger π-π stacking interaction with DNA in its triplet state, indicating a less favorable photoinduced transition but a stronger interaction with DNA postexcitation. We hope our study will inspire future experimental and computational work aimed at the design of novel Ru-polypyridyl drug candidates and that our new quantitative measure of π-π stacking interactions in DNA will find a general application in the field.

Structure and Chiroptical Properties of Anthra[1,2-a]anthracene-1-yl Dimers as New Biaryls

Dimers of anthra[1,2-a]anthracene-1-yl units and its mesityl derivative were synthesized by Ni(0)-mediated coupling of the corresponding chloro derivatives as new biaryls. The X-ray analysis and DFT calculations revealed that two polycyclic aromatic units with nonplanar deformations took a twisted conformation about the single bond as a chiral axis. Enantiomers of the nonsubstituted compound were resolved by chiral HPLC, and the enantiopure samples showed intense Cotton effects at 321 nm in the circular dichroism (CD) spectra and emission bands at 449 nm in the circularly polarized luminescence (CPL) spectra with dissymmetry factor of |glum| 3.6×10-3. The absolute stereochemistry of this biaryl was determined by the theoretical calculation of CD spectrum by the time-dependent DFT method. The barrier to enantiomerization was determined to be 108 kJ mol-1 at 298 K. The dynamic process proceeded via a stepwise mechanism involving the helical inversion of each aromatic unit and the rotation about the biaryl axis as analyzed by the DFT calculations.

New lignans from Phyllanthodendron dunnianum

A new lignan phyllanins A (1) and a lignan phyllanins B (2) for which the absolute configuration was determined for the first time, along with four known lignans (3-6) were isolated from the branch and leaf extracts of Phyllanthodendron dunnianum. Their planar structures were mainly determined by a combination of 1D and 2D NMR, HRESIMS spectral analyses, and the absolute configurations of the compounds 1 and 2 were established by DFT GIAO 13C NMR and electronic circular dichroism (ECD) calculations. In addition, all these six lignans were firstly tested for the antibacterial activities against MRSA, Staphylococcus aureus, Enterococcus faecalis, Pseudomonas aeruginosa and Escherichia coli. Among these compounds, 2 and 5 showed potential antibacterial activities against MRSA and S. aureus with MIC values of 4 and 8 µg/mL, respectively.

Ochrolines A-C, three new indole diketopiperazines from cultures of endophytic fungi Bionectria ochroleuca SLJB-2

Three new indole diketopiperazines, ochrolines A-C (1-3), along with three known compounds (4-6), were isolated and identified from the EtOAc extract of the solid fermentation of Bionectria ochroleuca SLJB-2. Notably, compound 1 featured a natural rarely-occurring caged skeleton with a 6/5/6/7 heterotetracyclic bridged ring system. The structures including absolute configurations of 1-3 were fully accomplished by extensive spectroscopic analyses, DFT GIAO 13C NMR and electronic circular dichroism (ECD) calculations. The plausible biogenetic pathways of these new indole diketopiperazines were also proposed. Moreover, the cytotoxic activity screening revealed that compound 2 exhibited moderate inhibitory effect against A549 with inhibition rate of 57.44% at the concentration of 50 μM and compound 1 exhibited mild inhibitory activities against A549, Hela and MCF-7.

Ion mobility mass spectrometry uncovers regioselectivity in the carboxylate-assisted C-H activation of palladium N-heterocyclic carbene complexes

Carboxylate-assisted Pd-catalyzed C-H bond activation constitutes a mild and versatile synthetic tool to efficiently and selectively cleave inert C-H bonds. Herein, we demonstrate a simple method to experimentally evaluate both reactivity and selectivity in such systems using mass spectrometry (MS) methods. The N-heterocyclic carbene (NHC) cations [(NHC)PdX]+, bearing as X- ligand bases commonly used to promote the C-H activation (carboxylates and bicarbonate), are generated in the gas-phase by ESI-MS. Their C-H bond activation at the N-bound groups of the NHC is then studied using Collision Induced Dissociation (CID) experiments. Ion Mobility Spectrometry (IM)-MS is exploited to identify a number of regioisomers associated with the distinctive site selective C-H activations. It is demonstrated that such C-H activation concomitant with acetic acid release occurs from a mixture of activated [(NHC-H)Pd(CH3CO2H)]+ and non-activated [(NHC)Pd(CH3CO2)]+ complexes. The identity of the X-type ligands (X = Cl-, carboxylates and bicarbonate) has a significant impact on the regioisomer branching ratio upon CID conditions. IM-MS in conjunction with a DFT mechanistic study is presented for the acetate-assisted C-H activation of the [(NHC)Pd(CH3CO2)]+ cation featuring butyl and aryl as N-donor groups.

Experimental and Theoretical Screening of Core Gold Nanoparticles and Their Binding Mechanism to an Anticancer Drug, 2-Thiouracil

This study demonstrated the capability of two readily available optical spectroscopy tools, namely UV-Vis absorption spectrophotometry and Raman/surface-enhanced Raman spectroscopy, to select in a rapid and noninvasive manner the most homogenous gold nanoparticle (AuNP) models and to identify their chemical binding mechanism to 2-thiouracil (2-TU). 2-TU is an anticancer drug of great promise in the antiproliferative and photothermal therapies of cancer. The citrate-capped AuNPs emerged as the most stable as well as time- and cost-effective AuNP model out of the three widely used colloidal nanocores (citrate-, borohydride-citrate-, and sodium dodecyl sulfate (SDS)-capped AuNPs) that were examined. 2-TU chemically attached to the relatively monodispersed AuNPs via a chemisorption mechanism. The 2-TU-AuNPs complex formed through the covalent bonding of the S atom of 2-TU to the nanosurface in a vertical orientation. The spectroscopic results were then confirmed with the help of density functional theory (DFT) calculations and other physicochemical characterization tools for nanomaterials such as transmission electron microscopy (TEM), dynamic light scattering (DLS), and zeta potential. Overall, the purified 2-TU-AuNPs were found to be spherical, had an average diameter of 25 ± 2 nm, a narrow size distribution (1-30 nm), a sharp localized surface plasmon resonance (LSPR) peak at 525 nm, and a negative surface charge (-14 mV).

A zinc-containing porphyrin aluminum MOF in sorption of diethyl sulfide vapor: mechanistic experimental and computational study

We report a mechanistic study of the interactions in the sorption of volatile organic sulfur compound (VOSC) diethyl sulfide (DES) by zinc porphyrin aluminum MOF (actAl-MOF-TCPPZn) compound 3. First, interactions were studied under dynamic conditions with the vapor of DES in flowing air, using in situ time-dependent ATR-FTIR spectroscopy in a controlled atmosphere with a new facile spectroscopic mini-chamber. The first binding site includes μ(O-H) and COO- groups as detected by characteristic peak shifts. Control experiments with a model compound, which lacks porosity and these groups, show no peak shifts. An additional insight was obtained by DFT computations using small clusters. The kinetics of sorption of DES by compound 3 is of the Langmuir adsorption model and pseudo-first order with rate constant robs = 0.442 ± 0.056 min-1. Sorption of DES under static conditions in saturated vapor results in stoichiometric adsorption complex [Al-MOF-TCPPZn]1(DES)4 characterized by spectroscopic, structural and gravimetric methods; the adsorbed amount is very high (381 mg g-1 sorbent). The repetitive sorption and desorption of DES are conducted, with facile regeneration. Finally, the mechanistic details were determined by Raman and photoluminescence (PL) spectroscopy using a confocal Raman microscope. Photoexcitation of compound 3 at 405 nm into the Soret band of the metalloporphyrin linker shows the characteristic PL peaks of Q-bands: the purely electronic Q(0-0) and first vibronic Q(0-1) bands. Upon interaction with DES, preferential quenching of PL from the Q(0-0) band occurs with a significant increase of the signal of the vibronic Q(0-1) band, reflecting bonding to the metalloporphyrin ring. Compound 3 is of interest to mechanistic studies of VOSCs, their removal from air, and optical chemo-sensing.

Koninginins X-Z, Three New Polyketides from Trichoderma koningiopsis SC-5

Koninginins X-Z (1-3), three novel polyketides, were isolated from the solid fermentation of the endophytic fungus Trichoderma koningiopsis SC-5. Their structures, including the absolute configurations, were comprehensively characterized by a combination of NMR spectroscopic methods, HRESIMS, 13C NMR, DFT GIAO 13C NMR, and electronic circular dichroism calculations as well as single crystal X-ray diffraction. In addition, all the compounds were evaluated for antifungal activity against Candida albicans.

Hydrogen-Bond-Driven Peptide Nanotube Formation: A DFT Study

DFT calculations were carried out to examine geometries and binding energies of H-bond-driven peptide nanotubes. A bolaamphiphile molecule, consisting of two N-α amido glycylglycine head groups linked by either one CH2 group or seven CH2 groups, is used as a building block for nanotube self-assembly. In addition to hydrogen bonds between adjacent carboxy or amide groups, nanotube formation is also driven by weak C-H· · ·O hydrogen bonds between a methylene group and the carboxy OH group, and between a methylene group and an amide O=C group. The intratubular O-H· · ·O=C hydrogen bonds account for approximately a third of the binding energies. Binding energies calculated with the wB97XD/DGDZVP method show that the hydrocarbon chains play a stabilizing role in nanotube self-assembly. The shortest nanotube has the length of a single monomer and a diameter than increases with the number of monomers. Lengthening of the tubular structure occurs through intertubular O-H· · ·O=C hydrogen bonds. The average intertubular O-H· · ·O=C hydrogen bond binding energy is estimated to change with the size of the nanotubes, decreasing slightly towards some plateau value near 15 kcal/mol according to the wB97XD/DGDZVP method.

Bracelet-like Complexes of Lithium Fluoride with Aromatic Tetraamides, and Their Potential for LiF-Mediated Self-Assembly: A DFT Study

Geometries and binding energies of complexes between a LiF molecule and a model aromatic tetraamide are obtained using various DFT methods. The tetraamide consists of a benzene ring and four amides positioned so that the LiF molecule can bind via Li⋯O=C or N-H⋯F interactions. The complex with both interactions is the most stable one, followed by the complex with only N-H⋯F interactions. Doubling the size of the former resulted in a complex with a LiF dimer sandwiched between the model tetraamides. In turn, doubling the size of the latter resulted in a more stable tetramer with bracelet-like geometry having the two LiF molecules also sandwiched but far apart from each other. Additionally, all methods show that the energy barrier to transition to the more stable tetramer is small. The self-assembly of the bracelet-like complex mediated by the interactions of adjacent LiF molecules is demonstrated by all computational methods employed.

Molecular Simulation of the Interaction of Diclofenac with Halogen (F, Cl, Br)-Encapsulated Ga12As12 Nanoclusters

Diclofenac is one of the most frequently consumed over-the-counter anti-inflammatory agents globally, and several reports have confirmed its global ubiquity in several environmental compartments. Therefore, the need to develop more efficient monitoring/sensing devices with high detection limits is still needed. Herein, quantum mechanical simulations using density functional theory (DFT) computations have been utilized to evaluate the nanosensing efficacy and probe the applicability of Ga₁₂As₁₂ nanostructure and its engineered derivatives (halogen encapsulation F, Br, Cl) as efficient adsorbent/sensor materials for diclofenac. Based on the DFT computations, it was observed that diclofenac preferred to interact with the adsorbent material by assuming a flat orientation on the surface while interacting via its hydrogen atoms with the As atoms at the corner of the GaAs cage forming a polar covalent As-H bond. The adsorption energies were observed to be in the range of -17.26 to -24.79 kcal/mol and therefore suggested favorable adsorption with the surface. Nonetheless, considerable deformation was observed for the Br-encapsulated derivative, and therefore, its adsorption energy was observed to be positive. Additionally, encapsulation of the GaAs nanoclusters with halogens (F and Cl) enhanced the sensing attributes by causing a decrease in the energy gap of the nanocluster. And therefore, this suggests the feasibility of the studied materials as potentiometric sensor materials. These findings could offer some implications for the potential application of GaAs and their halogen-encapsulated derivatives for electronic technological applications.

Liver-cell protective pyridones from the fungi Tolypocladium album dws120

Five previously undescribed pyridone derivatives, tolypyridones I-M, were identified from the solid rice medium fermented by Tolypocladium album dws120, along with two known compounds tolypyridone A (or trichodin A) and pyridoxatin. Their planar structures and partial relative configurations have been determined by careful interpretation of their spectroscopic data. The full assignment of the relative and absolute configurations of tolypyridones I-M was achieved by gauge-independent atomic orbital ¹³C NMR calculation, quantitative nuclear Overhauser effects based interatomic distance calculation, and electronic circular dichroism calculation. In addition, we have fully determined the configuration of tolypyridone A by X-ray diffraction analysis. In bioassay, tolypyridones I was able to restore cell viability and inhibit the release of alanine aminotransferase and aspartate aminotransferase for ethanol-induced LO2 cell, suggesting its potential as a liver protective agent.

New Insights into Adsorption Properties of the Tubular Au26 from AIMD Simulations and Electronic Interactions

Recently, we revealed the electronic nature of the tubular Au₂₆ based on spherical aromaticity. The peculiar structure of the Au₂₆ could be an ideal catalyst model for studying the adsorptions of the Au nanotubes. However, through Google Scholar, we found that no one has reported connections between the structure and reactivity properties of Au₂₆. Here, three kinds of molecules are selected to study the fundamental adsorption behaviors that occur on the surface of Au₂₆. When one CO molecule is adsorbed on the Au₂₆, the σ-hole adsorption structure is quickly identified as belonging to a ground state energy, and it still maintains integrity at a temperature of 500 K, where σ donations and π-back donations take place; however, two CO molecules make the structure of Au₂₆ appear with distortions or collapse. When one H₂ is adsorbed on the Au₂₆, the H-H bond length is slightly elongated due to charge transfers to the anti-bonding σ* orbital of H₂. The Au₂₆-H₂ can maintain integrity within 100 fs at 300 K and the H₂ molecule starts moving away from the Au₂₆ after 200 fs. Moreover, the Au₂₆ can act as a Lewis base to stabilize the electron-deficient BH₃ molecule, and frontier molecular orbitals overlap between the Au₂₆ and BH₃.

Infinitene as two fused helicoidal trails of fused rings: evaluation of the magnetic behavior of [12]infinitene and anionic species displaying global aromaticity and antiaromaticity

The unique formation of an infinity-shaped carbon backbone made exclusively from fused benzene rings has recently been achieved. The structure of [12]infinitene can be viewed as two fused [6]helicene structures with a central crossover section, depicting a global aromatic behavior along with the overall structure, with deshielding regions along both helicoidal axes. In addition, the ¹³C-NMR characteristics are discussed. The formation of a cumulative region involving the shielding regions from the aromatic rings is depicted along with the overall aesthetically pleasant structural backbone, which is enhanced at the crossover section. For the evaluated dianionic counterpart, the structure shows a deshielding region above the fused-ring trail and a helicoidal shielding region, ascribed to a global antiaromatic counterpart. The aromaticity is recovered and enhanced at the tetranionic state. Thus, the neutral and tetranionic states are able to build up a long-ranged shielding region, given by the global aromatic behavior, with an enhanced shielding region at the center of the crossover section displaying π-π stacked rings.

Ethanol Oxidation Reaction Mechanism on Gold Nanowires from Density Functional Theory

Thin gold nanowires (NWs) are materials that could be used as support in different chemical reactions. Using density functional theory (DFT) it was shown that NWs that form linear atomic chains (LACs) are suitable for stimulating chemical reactions. To this end, the oxidation reaction of ethanol supported on the LACs of Au-NWs was investigated. Two types of LACs were used for the study, one pure and the other with an oxygen impurity. The results showed that the oxygen atom in the LAC fulfills important functions throughout the reaction pathway. Before the chemical reaction, it was observed that the LAC with impurity gains structural stability, that is, the oxygen acts as an anchor for the gold atoms in the LAC. In addition, the LAC was shown to be sensitive to disturbances in its vicinity, which modifies its nucleophilic character. During the chemical reaction, the oxidation of ethanol occurs through two different reaction paths and in two stages, both producing acetaldehyde (CH₃ CHO). The different reaction pathways are a consequence of the presence of oxygen in the LAC (oxygen conditions the formation of reaction intermediates). In addition, the oxygen in the LAC also modifies the kinetic behavior in both reaction stages. It was observed that, by introducing an oxygen impurity in the LAC, the activation energy barriers decrease ∼69 % and ∼97 % in the first and second reaction stages, respectively.

DFT calculations in solution systems: solvation energy, dispersion energy and entropy

DFT calculations of reaction mechanisms in solution have always been a hot topic, especially for transition-metal-catalyzed reactions. The calculation of solvation energy is performed using either the polarizable continuum model (PCM) or the universal solvation model SMD. The PCM calculation is very sensitive to the choice of atomic radii to form a cavity, where the self-consistent isodensity PCM (SCI-PCM) has been recognized as the best choice and our IDSCRF radii can provide a similar cavity. Moving from a gas-phase case to a solution case, dispersion energy and entropy should be carefully treated. The solvent-solute dispersion is also important in solution systems, and it should be calculated together with the solute dispersion. Only half of the solvent-solute dispersion energy from the PCM calculation belongs to the solute molecules to maintain a thermal equilibrium between a solute molecule and its cavity, similar to the treatment of electrostatic energy. Relative solute dispersion energy should also be shared equally with the newly formed cavity. The entropy change from a gas phase to a liquid phase is quite large, but the modern quantum chemistry programs can only calculate the gas-phase translational entropy based on the idea-gas equation. In this review, we will provide an operable method to calculate the solution translational entropy, which has been coded in our THERMO program.

Three-in-one: exploration of co-encapsulation of cabazitaxel, bicalutamide and chlorin e6 in new mixed cyclodextrin-crosslinked polymers

Three-in-one: a single bCyD polymer easily prepared in water is used to co-encapsulate cabazitaxel and bicalutamide with chlorin e6 affording a nanoplatform to implement multimodal cancer therapy.

Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity

This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.

Assessing the Performance of Al12N12 and Al12P12 Nanostructured Materials for Alkali Metal Ion (Li, Na, K) Batteries

This study focused on the potential of aluminum nitride (Al₁₂N₁₂) and aluminum phosphide (Al₁₂P₁₂) nanomaterials as anode electrodes of lithium-ion (Li-ion), sodium-ion (Na-ion), and potassium-ion (K-ion) batteries as investigated via density functional theory (DFT) calculations at PBE0-D3, M062X-D3, and DSDPBEP86 as the reference method. The results show that the Li-ion battery has a higher cell voltage with a binding energy of -1.210 eV and higher reduction potential of -6.791 kcal/mol compared to the sodium and potassium ion batteries with binding energies of -0.749 and -0.935 eV and reduction potentials of -6.414 and -6.513 kcal/mol, respectively, using Al₁₂N₁₂ material. However, in Al₁₂P₁₂, increases in the binding energy and reduction potential were observed in the K-ion battery with values -1.485 eV and -7.535 kcal/mol higher than the Li and Na ion batteries with binding energy and reduction potential -1.483, -1.311 eV and -7.071, -7.184 eV, respectively. Finally, Al₁₂N₁₂ and Al₁₂P₁₂ were both proposed as novel anode electrodes in Li-ion and K-ion batteries with the highest performances.

Efficient Computation of the Interaction Energies of Very Large Non-covalently Bound Complexes

We present a new benchmark set consisting of 16 large non-covalently bound systems (LNCI16) ranging from 380 up to 1988 atoms and featuring diverse interaction motives. Gas-phase interaction energies are calculated with various composite DFT, semi-empirical quantum mechanical (SQM), and force field (FF) methods and are evaluated using accurate DFT reference values. Of the employed QM methods, PBEh-3c proves to be the most robust for large systems with a relative mean absolute deviation (relMAD) of 8.5% with respect to the reference interaction energies. r2SCAN-3c yields an even smaller relMAD, at least for the subset of complexes for which the calculation could be converged, but is less robust for systems with smaller HOMO–LUMO gaps. The inclusion of Fock-exchange is therefore important for the description of very large non-covalent interaction (NCI) complexes in the gas phase. GFN2-xTB was found to be the best performer of the SQM methods with an excellent result of only 11.1% deviation. From the assessed force fields, GFN-FF and GAFF achieve the best accuracy. Considering their low computational costs, both can be recommended for routine calculations of very large NCI complexes, with GFN-FF being clearly superior in terms of general applicability. Hence, GFN-FF may be routinely applied in supramolecular synthesis planning.

1 Introduction

2 The LNCI16 Benchmark Set

3 Computational Details

4 Generation of Reference Values

5 Results and Discussion

6 Conclusions

Five new secondary metabolites from an endophytic fungus Phomopsis sp. SZSJ-7B

Two previously undescribed lactones, phomolides A and B (1 and 2), and three new sesquiterpenoids, phomenes A-C (3-5), together with one known compound, colletotricholide A (6), were isolated from the endophytic fungus Phomopsis sp. SZSJ-7B. Their chemical structures, including the absolute configurations, were comprehensively established by extensive analyses of NMR, high-resolution electrospray ionization mass spectrometry, electronic circular dichroism powered by theoretical calculations, and X-ray diffractions. Moreover, the cytotoxic and antibacterial activities of compounds 1-6 were also evaluated, and the results demonstrated that compound 2 showed significant antibacterial effects towards methicillin-resistant Staphylococcus aureus and S. aureus strains with minimum inhibitory concentration as low as 6.25 μg/ml, which was comparable to that of the clinical drug vancomycin. Moreover, all compounds showed no cytotoxic activity.

Magnetic response properties of carbon nano-onions

The magnetic response of a number of double- and triple-layer carbon nano-onions (CNOs) is analyzed by calculating the magnetically induced current density and the induced magnetic field using the pseudo-π model. Qualitatively the same magnetic response was obtained in calculations at the all-electron level. The calculations show that the CNOs exhibit strong net diatropic (paratropic) ring currents when the external magnetic field points perpendicularly to one of the six-membered (five-membered) rings. They are deshielded inside and shielded outside the CNO; the latter dominates for larger CNOs. The magnetic response originates from a combination of spherical paratropic current densities on the inside of each carbon layer and diatropic current densities on the outside of them. The quantitative differences in the aromaticity of the CNOs as compared to single fullerenes are discussed in terms of ring-current strengths. The magnetic response of some of the CNOs is approximately the sum of the magnetic response of the individual layers, whereas deviations are significant for CNOs containing C₈₀. For the largest CNOs, the deviation from the sum of the fullerene contributions is larger, especially when the external magnetic field is perpendicular to a six-membered ring.

Reassignment of the structures of pestalopyrones A-D

Pestalopyrones A-D are four unusual tricyclic pyrone derivatives with flexible chiral structures, isolated from the endophytic fungus Pestalotiopsis neglecta S3. The full elucidation of their structures was a challenging task, and remained unsolved in the original article. Herein, the relative configurations of pestalopyrones A and pestalopyrones B were unambiguously assigned by detailed analyses on spectroscopic data and GIAO ¹³C NMR calculation method with sorted training sets (STS). The planar structures of pestalopyrones C and pestalopyrones D were revised by reinterpretation of their reported spectroscopic data, and then their relative configurations were deduced by STS GIAO ¹³C NMR calculation and NOE analysis. The absolute configurations of all the mentioned compounds were determined by the comparison of their experimental and calculated ECD curves.

Non-covalent interactions in the glutathione peroxidase active center and their influence on the enzyme activity

In this computational work, the structure of the active center of the enzyme glutathione peroxidase (in three forms -SeH, -SeOH and -Se(O)OH) and the non-covalent interactions in it were investigated using modern quantum chemistry methods. The non-covalent interactions are described in detail. The presence of σ-hole interactions (chalcogen, tetrel and pnictogen bonds) formed mostly by a selenium atom as an electrophile in the glutathione peroxidase active center is confirmed for the first time. It is shown that a number of non-covalent interactions stabilize intermediates along the catalytic cycle and that modelling of the whole enzyme active center is necessary for accurate predictions of thermodynamic parameters, in particular, activation barriers.

Non-intersecting ring currents in [12]infinitene

The aromaticity of the newly synthesized [12]infinitene is addressed via analysis of the magnetically induced current density and the induced magnetic field. Our calculations reveal that [12]infinitene responds to an external magnetic field by creating two current-density pathways that flow diatropically along the edges of the molecule. The current-density pathways do not intersect. The entire structure is completely shielded suggesting that the infinitene molecule is aromatic, contrary to what the Möbius rule for twisted annulene structures predicts. We also studied the dication of [12]infinitene, which sustains two paratropic ring currents flowing along the edges. The space between the stacked rings at the crossing point is shorter for the dication as compared to the neutral molecule. Hence, a strong through-space current density appears at the crossing point of π-π stacked rings.

Chalcogen···Chalcogen Bonding in Molybdenum Disulfide, Molybdenum Diselenide and Molybdenum Ditelluride Dimers as Prototypes for a Basic Understanding of the Local Interfacial Chemical Bonding Environment in 2D Layered Transition Metal Dichalcogenides

An attempt was made, using computational methods, to understand whether the intermolecular interactions in the dimers of molybdenum dichalcogenides MoCh2 (Ch = chalcogen, element of group 16, especially S, Se and Te) and similar mixed-chalcogenide derivatives resemble the room temperature experimentally observed interactions in the interfacial regions of molybdenites and their other mixed-chalcogen derivatives. To this end, MP2(Full)/def2-TVZPPD level electronic structure calculations on nine dimer systems, including (MoCh2)2 and (MoChCh′2)2 (Ch, Ch′ = S, Se and Te), were carried out not only to demonstrate the energetic stability of these systems in the gas phase, but also to reproduce the intermolecular geometrical properties that resemble the interfacial geometries of 2D layered MoCh2 systems reported in the crystalline phase. Among the six DFT functionals (single and double hybrids) benchmarked against MP2(full), it was found that the double hybrid functional B2PLYPD3 has some ability to reproduce the intermolecular geometries and binding energies. The intermolecular geometries and binding energies of all nine dimers are discussed, together with the charge density topological aspects of the chemical bonding interactions that emerge from the application of the quantum theory of atoms in molecules (QTAIM), the isosurface topology of the reduced density gradient noncovalent index, interaction region indicator and independent gradient model (IGM) approaches. While the electrostatic surface potential model fails to explain the origin of the S···S interaction in the (MoS2)2 dimer, we show that the intermolecular bonding interactions in all nine dimers examined are a result of hyperconjugative charge transfer delocalizations between the lone-pair on (Ch/Ch′) and/or the π-orbitals of a Mo–Ch/Ch′ bond of one monomer and the dπ* anti-bonding orbitals of the same Mo–Ch/Ch′ bond in the second monomer during dimer formation, and vice versa. The HOMO–LUMO gaps calculated with the MN12-L functional were 0.9, 1.0, and 1.1 eV for MoTe2, MoSe2 and MoS2, respectively, which match very well with the solid-state theoretical (SCAN-rVV10)/experimental band gaps of 0.75/0.88, 0.90/1.09 and 0.93/1.23 eV of the corresponding systems, respectively. We observed that the gas phase dimers examined are perhaps prototypical for a basic understanding of the interfacial/inter-layer interactions in molybdenum-based dichalcogenides and their derivatives.

Short Pyridine-Furan Springs Exhibit Bistable Dynamics of Duffing Oscillators

The intensive development of nanodevices acting as two-state systems has motivated the search for nanoscale molecular structures whose dynamics are similar to those of bistable mechanical systems, such as Euler arches and Duffing oscillators. Of particular interest are the molecular structures capable of spontaneous vibrations and stochastic resonance. Recently, oligomeric molecules that were a few nanometers in size and exhibited the bistable dynamics of an Euler arch were identified through molecular dynamics simulations of short fragments of thermo-responsive polymers subject to force loading. In this article, we present molecular dynamics simulations of short pyridine-furan springs a few nanometers in size and demonstrate the bistable dynamics of a Duffing oscillator with thermally-activated spontaneous vibrations and stochastic resonance.

Evaluation of packing single and multiple atoms and molecules in the porous organic cage CC3-R

The absorption of multiple atoms and molecules, including Kr, Xe, CH₄, CO₂, C₂H₂, H₂O, and SF₆, within CC3-R, a Porous Organic Cage (POC), was calculated and analyzed. The CC3-R molecule has one central cavity and four window sites. Most adsorbents were modeled with either one unit in the central cavity, four units in the window sites, or with five units in both sites. For Xe, the most favorable site was the central one. The CO₂ molecule binds about 3 kcal mol^-1 in free energy more strongly than CH₄ in the central cavity of CC3-R at 300 K which may be enough to allow useful discrimination. Four C₂H₂ units and four CO₂ units are calculated to bind similarly inside CC3-R (ΔH(298 K) = -8.6 and -7.7 kcal mol^-1 per unit, respectively). Since H₂O is smaller, more waters can easily fit inside. For twelve water molecules, the binding enthalpy per water is ΔH(298 K) = -16.4 kcal mol^-1. For comparison, the binding enthalpy of (H₂O)₁₂ at the same level of theory (B3LYP/6-31G(d,p)-D3BJ//M06-2X/6-31G(d)) is predicted to be -12.3 kcal mol^-1 per water. Finally, the dimerization of CC3-R and the association of CC3-R with CC3-S was studied as well as 3 to 9 iodine atoms enclosed in CC3-R.

Transient absorption spectroscopy of the electron transfer step in the photochemically activated polymerizations of N-ethylcarbazole and 9-phenylcarbazole

The polymerization of photoexcited N-ethylcarbazole (N-EC) in the presence of an electron acceptor begins with an electron transfer (ET) step to generate a radical cation of N-EC (N-EC˙+). Here, the production of N-EC˙+ is studied on picosecond to nanosecond timescales after N-EC photoexcitation at a wavelength λex = 345 nm using transient electronic and vibrational absorption spectroscopy. The kinetics and mechanisms of ET to diphenyliodonium hexafluorophosphate (Ph2I+PF6-) or para-alkylated variants are examined in dichloromethane (DCM) and acetonitrile (ACN) solutions. The generation of N-EC˙+ is well described by a diffusional kinetic model based on Smoluchowski theory: with Ph2I+PF6-, the derived bimolecular rate coefficient for ET is kET = (1.8 ± 0.5) × 1010 M-1 s-1 in DCM, which is consistent with diffusion-limited kinetics. This ET occurs from the first excited singlet (S1) state of N-EC, in competition with intersystem crossing to populate the triplet (T1) state, from which ET may also arise. A faster component of the ET reaction suggests pre-formation of a ground-state complex between N-EC and the electron acceptor. In ACN, the contribution from pre-reaction complexes is smaller, and the derived ET rate coefficient is kET = (1.0 ± 0.3) × 1010 M-1 s-1. Corresponding measurements for solutions of photoexcited 9-phenylcarbazole (9-PC) and Ph2I+PF6- give kET = (5 ± 1) × 109 M-1 s-1 in DCM. Structural modifications of the electron acceptor to increase its steric bulk reduce the magnitude of kET: methyl and t-butyl additions to the para positions of the phenyl rings (para Me2Ph2I+PF6- and t-butyl-Ph2I+PF6-) respectively give kET = (1.2 ± 0.3) × 1010 M-1 s-1 and kET = (5.4 ± 1.5) × 109 M-1 s-1 for reaction with photoexcited N-EC in DCM. These reductions in kET are attributed to slower rates of diffusion or to steric constraints in the ET reaction.

Azide···Oxygen Interaction: A Crystal Engineering Tool for Conformational Locking

We have designed, synthesized and crystallized 36 compounds, each containing an azide group and an oxygen atom separated by three bonds. Crystal structure analysis revealed that each of these molecules adopts a conformation in which the azide and oxygen groups orient syn to each other with a short O ··· N b  contact. Geometry-optimized structures [using  M06-2X/6-311G(d,p) level of theory ] also showed the syn conformation in all 36 of these cases, suggesting that this not merely a crystal packing effect. Quantum topological analysis using Bader's Atoms in Molecules (AIM) theory revealed bond paths and bond critical points (BCP) in these structures suggesting its nature and energetics to be similar to weak hydrogen bonding. The NCI-RDG plot clearly revealed the attractive interaction consisting of electrostatic or dispersive components in all the 36 systems. NBO analysis suggested a weak orbital-relaxation (charge-transfer) contribution of energy for a few  (sp2) O-donor systems. Natural population analysis (NPA) and molecular electrostatic potential mapping (MESP) of these crystal structures further revealed the existence of favorable azide-oxygen interaction.   A CSD search indicated the frequent and consistent occurrence of this interaction and its role dictating the syn conformation of azide and oxygen in molecules where these groups are separated by 2-4 bonds.

Exceptionally Long Covalent CC Bonds-A Local Vibrational Mode Study

For decades one has strived to synthesize a compound with the longest covalent C−C bond applying predominantly steric hindrance and/or strain to achieve this goal. On the other hand electronic effects have been added to the repertoire, such as realized in the electron deficient ethane radical cation in its D3d form. Recently, negative hyperconjugation effects occurring in diamino-o-carborane analogs such as di-N,N-dimethylamino-o-carborane have been held responsible for their long C−C bonds. In this work we systematically analyzed CC bonding in a diverse set of 53 molecules including clamped bonds, highly sterically strained complexes such as diamondoid dimers, electron deficient species, and di-N,N-dimethylamino-o-carborane to cover the whole spectrum of possibilities for elongating a covalent C−C bond to the limit. As a quantitative intrinsic bond strength measure, we utilized local vibrational CC stretching force constants ka(CC) and related bond strength orders BSO n(CC), computed at the ωB97X-D/aug-cc-pVTZ level of theory. Our systematic study quantifies for the first time that whereas steric hindrance and/or strain definitely elongate a C−C bond, electronic effects can lead to even longer and weaker C−C bonds. Within our set of molecules the electron deficient ethane radical cation, in D3d symmetry, acquires the longest C−C bond with a length of 1.935 Å followed by di-N,N-dimethylamino-o-carborane with a bond length of 1.930 Å. However, the C−C bond in di-N,N-dimethylamino-o-carborane is the weakest with a BSO n value of 0.209 compared to 0.286 for the ethane radical cation; another example that the longer bond is not always the weaker bond. Based on our findings we provide new guidelines for the general characterization of CC bonds based on local vibrational CC stretching force constants and for future design of compounds with long C−C bonds.

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

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

Linker Expansion and Its Impact on Switchability in Pillared-Layer MOFs

Linker elongation is an important method to systematically adjust porosity and pore size in isoreticular MOFs. In flexible structures, this approach opens the possibility for the systematic analysis of the building blocks and their contribution to the overall flexible behavior enabling tuning of the framework responsivity toward molecular stimuli. In this work, we report two new compounds isoreticular to the highly flexible pillared layer structure DUT-8(Ni) ([Ni₂(2,6-ndc)₂(dabco)]_n, 2,6-ndc = 2,6-naphthalenedicarboxylate, dabco = 1,4-diazabicylo[2.2.2]octane). Aromatic linker 2,6-ndc was substituted by longer carboxylic linkers, namely, 4,4'-biphenyldicarboxylate (4,4'-bpdc) and 4,4'-stilbenedicarboxylate (4,4'-sdc), while the dabco pillar was retained. The structural response of the new compounds toward the desolvation and adsorption of various fluids was studied using advanced in situ PXRD techniques, demonstrating distinct differences in the flexible behavior of three compounds and disclosing the impact of linker structure on the framework response. Theoretical calculations provide mechanistic insights and an energetic rationale for the pronounced differences in switchability observed. The energetics of linker bending and linker-linker dispersion interactions govern the phase transitions in investigated MOFs.

Biodegradation of anionic polyacrylamide mediated by laccase and amidase: docking, virtual mutation based on affinity and DFT study

The aim of this work was to document the elucidation of a mechanism as a reference.