Small molecules (CO2 , iPrNCO, and iPrNCNiPr) activation by the metallomimetics (μ-Hydrido) diborane anion: A DFT investigation on mechanism and chemoselectivity controlling

Main-group metallomimetics provide a new way to replace transition metal complexes to activate inert small molecules under mild conditions. In this work, the activation mechanisms of CO2 , iPrNCO, and iPrNCNiPr by (μ-Hydrido) diborane anion ([1H]- ) have been investigated by density functional theory (DFT) calculations. Two different activation sites, BB versus BH bond of [1H]- , are investigated and compared. The results show that these inert molecules can be activated by [1H]- through cycloadditions under mild conditions. The reactions with iPrNCO and iPrNCNiPr are dynamic and thermodynamic controlling, the obtained products are related not only to the energy barrier but also to the stability of the products. Moreover, the competition for BB/BH bond site activation is directly related to the steric effect of small molecules. CO2 , which is without steric hindrance, can only be activated by the BB bond, whereas iPrNCNiPr can only be activated by the BH bond due to the large steric effect. The medium iPrNCO can be activated not only by the BB bond but also by the BH bond. Our study provides theoretical explanations for the reaction activity and chemoselectivity controlling of the title reaction, and displays the potential applications for compounds containing boron-boron bonds and inert small molecule activation.

When, Where and Why Boron Prefers Boron to Nitrogen

Boron is the archetypal Lewis acid, and therefore it is only natural that it prefers to bind nitrogen, its usual Lewis base counterpart. To challenge this assumption, we present a computationally designed bicyclopentane molecule akin to [1.1.1]propellane, but with pyramidal B and N inner atoms bonded by an "inverted" dative bond. Unexpectedly, the dimer of this system prefers to interact via an atypical boron-boron bond over the supposedly obvious boron-nitrogen bond. A molecular orbital analysis shows that the boron in this peculiar entity acts both as an electron donor and an electron acceptor, making the dimerization an amphoteric-amphoteric interaction process.

Ab initio Modeling of Hydrogen Bonding of Remdesivir and Adenosine with Uridine

Remdesivir (RDV) emerged as an effective drug against the SARS-CoV-2 virus pandemic. One of the crucial steps in the mechanism of action of RDV is its incorporation into the growing RNA strand. RDV, an adenosine analogue, forms Watson-Crick (WC) type hydrogen bonds with uridine in the complementary strand and the strength of this interaction will control efficacy of RDV. While there is a plethora of structural and energetic information available about WC H-bonds in natural base pairs, the interaction of RDV with uridine has not been studied yet at the atomic level. In this article, we aim to bridge this gap, to understand RDV and its hydrogen bonding interactions, by employing density functional theory (DFT) at the M06-2X/cc-pVDZ level. The interaction energy, QTAIM analysis, NBO and SAPT2 are performed for RDV, adenosine, and their complex with uridine to gain insights into the nature of hydrogen bonding. The computations show that RDV has similar geometry, energetic, molecular orbitals, and aromaticity as adenosine, suggesting that RDV is an effective adenosine analogue. The important geometrical parameters, such as bond distances and red-shift in the stretching vibrational modes of adenosine, RDV and uridine identify two WC-type H-bonds. The relative strength of these two H-bonds is computed using QTAIM parameters and the computed hydrogen bond energy. Finally, the SAPT2 study is performed at the minima and at non-equilibrium base pair distances to understand the dominant intermolecular physical force. This study, based on a thorough analysis of a variety of computations, suggests that both adenosine and RDV have similar structure, energetic, and hydrogen bonding behaviour.

Cross-Electrophile C-PIII Coupling of Chlorophosphines with Organic Halides: Photoinduced PIII and Aminoalkyl Radical Generation Enabled by Pnictogen Bonding

Pnictogen bonding (PnB) has gained recognition as an appealing strategy for constructing novel architectures and unlocking new properties. Within the synthetic community, the development of a straightforward and much simpler protocol for cross-electrophile C-PIII coupling remains an ongoing challenge with organic halides. In this study, we present a simple strategy for photoinduced PnB-enabled cross-electrophile C-PIII couplings using readily available chlorophosphines and organic halides via merging single electron transfer (SET) and halogen atom transfer (XAT) processes. In this photomediated transformation, the PnB formed between chlorophosphines and alkyl amines facilitates the photogeneration of PIII radicals and α-aminoalkyl radicals through SET. Subsequently, the resulting α-aminoalkyl radicals activate C-X bonds via XAT, leading to the formation of carbon radicals. This methodology offers operational simplicity and compatibility with both aliphatic and aromatic chlorophosphines and organic halides.

A review on medicinally important heterocyclic compounds and importance of biophysical approach of underlying the insight mechanism in biological environment

Heterocyclic derivatives have more interesting biological properties which hold a remarkable place in pharmaceutical industries due to their unique physiochemical properties and ease of adaption in various biological environments. Of many, the above-said derivatives have been recently examined for their promising action against a few malignancies. Specifically, anti-cancer research has benefited from these derivatives' natural flexibility and dynamic core scaffold. In any case, concerning some other promising anti-cancer drugs, heterocyclic derivative doesn't come without deficiencies. To be a successful drug candidate it should poses Absorption, Distribution, Metabolism and Eliminations (ADME) parameter, and must also have good binding interaction towards carrier protein as well as DNA and less in toxic nature, economically feasible. In this review, we described the overview of biologically important heterocyclic derivatives and their main application in medicine. Further, we focus types of biophysical techniques to understand the binding interaction mechanism.Communicated by Ramaswamy H. Sarma.

A fast approximate extension of the interacting quantum atoms energy decomposition to excited states

An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact V_xc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.

Alkynyl Sulfonium Salts Can Be Employed as Chalcogen-Bonding Catalysts and Generate Alkynyl Radicals under Blue-Light Irradiation

Chalcogen bonding (ChB) has emerged as a promising tool in organic synthesis. However, compared with the well-developed selenium- and tellurium-based salt catalysts, the ChB catalysis of sulfonium salts is still unknown. Here, we report a new type of alkynyl-sulfonium salt ChB catalysis for various ionic transformations, including transfer hydrogenation, bromination, bromolactonization, dimerization of 1,1-diphenylethylene, nitro-Michael addition reaction and Ritter reaction. More importantly, the photocapability of ChB was first demonstrated to generate alkynyl radicals for the synthesis of a variety of chalcogenoacetylenes. Mechanistic studies shed light on the mechanism of the photoinduced reactions and confirmed the involvement of alkynyl radicals which are difficult to generate otherwise.

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol^-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.

To be or not to be? that is the entropic, enthalpic, and molecular interaction dilemma in the formation of (water)20 clusters and methane clathrate

A detailed analysis under a comprehensive set of theoretical and computational tools of the thermodynamical factors and of the intermolecular interactions behind the stabilization of a well known set of (water)20 cavities and of the methane clathrate is offered in this work. Beyond the available reports of experimental characterization at extreme conditions of most of the systems studied here, all clusters should be amenable to experimental detection at 1 atm and moderate temperatures since 280 K marks the boundary at which, ignoring reaction paths, formation of all clusters is no longer spontaneous from the 20H2O → (H2O)20 and CH4 + 20H2O → CH4@512 processes. As a function of temperature, a complex interplay leading to the free energy of formation occurs between the destabilizing entropic contributions, mostly due to cluster vibrations, and the stabilizing enthalpic contributions, due to intermolecular interactions and the PV term, is best illustrated by the highly symmetric 512 cage consistently showing signs of stronger intermolecular bonding despite having smaller binding energy than the other clusters. A fluxional wall of attractive non-covalent interactions, arising because of the cumulative effect of a large number of tiny individual charge transfers to the interstitial region, plays a pivotal role stabilizing the CH4@512 clathrate.

Electrostatics Explains the Reverse Lewis Acidity of BH3 and Boron Trihalides: Infrared Intensities and a Relative Energy Gradient (REG) Analysis of IQA Energies

The reaction path for the formation of BX₃–NH₃ (X = H, F, Cl, Br) complexes was divided into two processes: (i) rehybridization of the acid while adopting a pyramidal geometry, and (ii) the complex formation from the pyramidal geometries of the acid and base. The interacting quantum atom (IQA) method was used to investigate the Lewis acidity trend of these compounds. This topological analysis suggests that the boron–halogen bond exhibits a considerable degree of ionicity. A relative energy gradient (REG) analysis on IQA energies indicates that the acid–base complex formation is highly dependent on electrostatic energy. With increasing halogen electronegativity, a higher degree of ionicity of the B–X is observed, causing an increase in the absolute value of X and B charges. This increases not only the attractive electrostatic energy between the acid and base but also enhances the repulsive energy. The latter is the main factor behind the acidity trend exhibited by trihalides. Changes in geometry are relevant only for complexes where BH₃ acts as an acid, where lower steric hindrance facilitates the adoption of the pyramidal geometry observed in the complex. The CCTDP analysis shows that infrared intensities of BX₃–NH₃ are determined mostly by the atomic charges and not by the charge transfer or polarization. The opposite is observed in covalent analogues.

Thermodynamics and Intermolecular Interactions during the Insertion of Anionic Naproxen into Model Cell Membranes

The insertion process of Naproxen into model dimyristoylphosphatidylcholine (DMPC) membranes is studied by resorting to state-of-the-art classical and quantum mechanical atomistic computational approaches. Molecular dynamics simulations indicate that anionic Naproxen finds an equilibrium position right at the polar/nonpolar interphase when the process takes place in aqueous environments. With respect to the reference aqueous phase, the insertion process faces a small energy barrier of ≈5 kJ mol^-1 and yields a net stabilization of also ≈5 kJ mol^-1. Entropy changes along the insertion path, mainly due to a growing number of realizable microstates because of structural reorganization, are the main factors driving the insertion. An attractive fluxional wall of noncovalent interactions is characterized by all-quantum descriptors of chemical bonding (natural bond orbitals, quantum theory of atoms in molecules, noncovalent interaction, density differences, and natural charges). This attractive wall originates in the accumulation of tiny transfers of electron densities to the interstitial region between the fragments from a multitude of individual intermolecular contacts stabilizing the tertiary drug/water/membrane system.

Metal–ligand interactions in complexes of cyclen-based ligands with Bi and Ac

The structural and bonding properties of Bi and Ac complexes with cyclen-based chelating ligands have been studied using relativistic DFT calculations in conjunction with TZ2P all-electron basis sets. Besides the parent cyclen ligand, the study has covered its extensions with pyridine-type (Lpy), carboxylate (DOTA, DOTPA), picolinate (MeDO2PA) and phosphonate (DOTMP) pendant arms. The effect of the cyclen ring size has been probed by increasing it from [12]aneN4 to [16]aneN4. Additional extensions in the DOTA complexes included the H2O ligand at the 9th coordination site as well as the p-SCN-Bn substituent (a popular linker to the targeting vector). The study focuses on the complex stability, the nature of bonding and the differences between Ac and Bi in the complexes. The metal–ligand interactions have been analysed by the Extended Transition State method combined with Natural Orbitals of Chemical Valence theory and Quantum Theory of Atoms in Molecules models.

A Highly Luminescent Nitrogen-Doped Nanographene as an Acid- and Metal-Sensitive Fluorophore for Optical Imaging

Dibenzo[hi,st]ovalene (DBOV) has excellent photophysical properties, including strong fluorescence and high ambient stability. Moreover, the optical blinking properties of DBOV have enabled optical super-resolution single-molecule localization microscopy with an imaging resolution beyond the diffraction limit. Various organic and inorganic fluorescent probes have been developed for super-resolution imaging, but those sensitive to pH and/or metal ions have remained elusive. Here, we report a diaza-derivative of DBOV (N-DBOV), synthesized in eight steps with a total yield of 15%. Nitrogen (N)-bearing zigzag edges were formed through oxidative cyclization of amino groups in the last step. UV–vis and fluorescence spectroscopy of N-DBOV revealed its promising optical properties comparable to those of the parent DBOV, while cyclic voltammetry and density functional theory calculations highlighted its lower orbital energy levels and potential n-type semiconductor character. Notably, in contrast to that of the parent DBOV, the strong luminescence of N-DBOV is dependent on pH and the presence of heavy metal ions, indicating the potential of N-DBOV in sensing applications. N-DBOV also exhibited pH-responsive blinking, which enables pH-sensitive super-resolution imaging. Therefore, N-DBOV appears to be a highly promising candidate for fluorescence sensing in biology and environmental analytics.

The competition and cooperativity of hydrogen/halogen bond and π-hole bond involving the heteronuclear ethylene analogues

The noncovalent interactions involving heteronuclear ethylene analogues H₂ CEH₂ (E = Si, Ge and Sn) have been studied by the Møller-Plesset perturbation theory to investigate the competition and cooperativity between the hydrogen/halogen bond and π-hole bond. H₂ CEH₂ has a dual role of being a Lewis base and acid with the region of π-electron accumulation above the carbon atom and the region of π-electron depletion (π-hole) above the E atom to participate in the NCX···CE (X = H and Cl) hydrogen/halogen bond and CE···NCY (Y = H, Cl, Li and Na) π-hole bond, respectively. When HCN/ClCN interacts with H₂ CEH₂ by two sites, the strength of hydrogen bond/halogen bond is stronger than that of π-hole bond. The π-hole bond becomes obviously stronger when the metal substituent of YCN (Y = Li and Na) interacting with H₂ CEH₂ , showing the character of partial covalent, its strength is much greater than that of hydrogen/halogen bond. In the ternary complexes, both hydrogen/halogen bond and π-hole bond are simultaneously strengthened compared to those in the binary complexes, especially in the systems containing alkali metal.

Multipolar electrostatics for hairpin and pseudoknots in RNA: Improving the accuracy of force field potential energy function

Molecular dynamics (MD) simulations that rely on force field methods has been widely used to explore the structure and function of RNAs. However, the current commonly used force fields are limited by the electrostatic description offered by atomic charge, dipole and at most quadrupole moments, failing to capture the anisotropic picture of electronic features. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. A multipolar electrostatic model based on high rank multipole moments is described in this work, which allows us to combine polarizability and anisotropy of electron density. RNA secondary structure was taken as a research system, and its substructures including stem, loops (hairpin loop, bulge loop, internal loop, and multi-branch loop), and pseudoknots (H-type and K-type) were investigated, respectively, as well as the hairpin. First, the atom-atom electrostatic properties derived from one chain of a duplex RNA 2MVY in our previous work (Ref. 58) were measured by the pilot RNA systems of hairpin, hairpin loop, stem, and H-type pseudoknot, respectively. The prediction results were not satisfactory. Consequently, to obtain a general set of electrostatic parameters for RNA force fields, the convergence behavior of the atom-atom electrostatic interactions in the pilot RNA systems was explored using high rank atomic multipole moments. The pilot RNA systems were cut into four types of different-sized molecular fragments, and the single nucleotide fragment and nucleotide-paired fragment proved to be the most reasonable systems for base-unpairing regions and base-pairing regions to investigate the convergence behavior of all types of atom-atom electrostatic interactions, respectively. Transferability of the electrostatic properties drawn from the pilot RNA systems to the corresponding test systems was also investigated. Furthermore, the convergence behavior of atomic electrostatic interactions in other substructures including bulge loop, internal loop, multi-branch loop, and K-type pseudoknot was expected to be modeled via the hairpin.

Modeling Hydrogen Release from Water with Borane and Alane Catalysts: A Unified Reaction Valley Approach

The unified reaction valley approach combined with the local vibrational mode and ring puckering analysis is applied to investigate the hydrogen evolution from water in the presence of small hydrides such as BH₃, metal hydrides as AlH₃, and their derivatives. We studied a series of reactions involving BH₃, AlH₃, B₂H₆, Al₂H₆, and AlH₃BH₃ with one- and two-water molecules, considering multiple reaction paths. In addition, the influence of the aqueous medium was examined. A general reaction mechanism was identified for most of the reactions. Those that deviate could be associated with unusually high reaction barriers with no hydrogen release. The charge transfer along the reaction path suggests that a viable hydrogen release is achieved when the catalyst adopts the role of a charge donor during the chemical processes. The puckering analysis showed that twistboat and boat forms are the predominant configurations in the case of an intermediate six-membered ring formation, which influences the activation barrier. The local mode analysis was used as a tool to detect the H-H bond formation as well as to probe catalyst regenerability. Based on the correlation between the activation energy and the change in the charge separation for cleaving O-H and B(Al)-H bonds, two promising subsets of reactions could be identified along with prescriptions for lowering the reaction barrier individually with electron-donating/withdrawing substituents.

Bonding and metastability for Group 12 dications

Electronic structure and bonding properties of the Group 12 dications M₂ ²⁺ (M = Zn, Cd, Hg) are investigated and electron density-derived quantities are used to characterize the metastability of these species. Of particular interest are the complementary descriptions afforded by the Laplacian of the electron density ∇² ρ(r) and the one-electron Bohm quantum potential (Q =  ∇ 2 ρ r / 2 ρ r ) along the bond path. Further, properties derived from the pair density including the localization-delocalization matrices (LDMs) and the interacting quantum atoms (IQA) energies are analyzed within the framework of the quantum theory of atoms in molecules (QTAIM). From the crossing points of the singlet (ground) and triplet (excited) potential energy curves, the barriers for dissociation (BFD) are estimated to be 25.2 kcal/mol (1.09 eV) for Zn₂ ²⁺ , 22.8 kcal/mol (0.99 eV) for Cd₂ ²⁺ , and 26.4 kcal/mol (1.14 eV) for Hg₂ ²⁺ . For comparison and benchmarking purposes, the case of N₂ ²⁺ is considered as a texbook example of metastability. At the equilibrium geometries, LDMs, which are used here as an electronic fingerprinting tool, discriminate and group together Group 12 M₂ ²⁺ from its isoelectronic Group 11 M₂ . While "classical" bonding indices are inconclusive in establishing regions of metastability in the bonding, it is shown that the one-electron Bohm quantum potential is promising in this regard.

Cooperativity in a cycloalkane-1,2/1,3-polyol corona: Topological hydrogen bonding in 1,2-diol motifs

A corona, consisting of 18 carbon atoms bearing 12 hydroxy groups in a continuous hydrogen-bonded chain, is built up by alternating degenerate conformations of alternating alkane-1,2-diol and 1,3-diol motifs. Geometries, proton nuclear magnetic resonance shifts and interaction energies for the dodecahydroxycyclo-octadecane and selected fragments are determined by density functional calculations at the B3LYP/6-311+G(d,p) level. Cooperative effects of O-H⋯O-H bonding are evident from the simple juxtaposition of these two motifs with a common OH group in butane-1,2,4-triol conformers. Bracketing a 1,2-diol motif with two 1,3-diol motifs in hexane-1,3,4,6-tetrol leads to a structure in which the 1,2-diol motif displays a bond critical point for hydrogen bonding. This is associated with enhancement of the shift of the hydrogen-bonded OH proton and of the corresponding H⋯O interaction energy. The full corona has a complete outer ring of O-H⋯O-H bond paths, and an inner ring of bond paths, due to C-H⋯H-C hydrogen-hydrogen bonding, which result in a central ring critical point. The topological O-H⋯O-H hydrogen bond, never seen in simple alkane-1,2-diols, is associated with cooperative enhancement of the H⋯O interaction energy, but this is not a necessary condition for a bond path: values for topological C-H⋯H-C hydrogen-hydrogen bonds can be as low as -0.4 kcal mol^-1 .

Carboxymethyl cellulose improved adsorption capacity of polypyrrole/CMC composite nanoparticles for removal of reactive dyes: Experimental optimization and DFT calculation

In this study, polypyrrole/carboxymethyl cellulose nanocomposite particles (PPy/CMC NPs) were synthesized and applied for removal of reactive red 56 (RR56)and reactive blue 160 (RB160) as highly toxic dyes. The amount of CMC was found significantly effective on the surface adsorption efficiency. Different optimization methods including the genetic programming, response surface methodology, and artificial neural network (ANN) were used to optimize the effect of different parameters including pH, adsorption time, initial dye concentration and adsorbent dose. The maximum adsorption of RR56 and RB160 were found under the following optimum conditions: pH of 4 and 5, adsorption time of 55 min and 52 min for RR56 and RB160, respectively, initial dye concentration of 100 mg/L and adsorbent dose of 0.09 g for both dyes. were obtained for RR56 and RB160, respectively. Also, the results indicated that ANN method could predict the experimental adsorption data with higher accuracy than other methods. The analysis of ANN results indicated that the adsorbent dose is the main factor in RR56 removal, followed by time, pH and initial concentration, respectively. However, initial concentration mostly determines the RB160 removal process. The isotherm data for both dyes followed the Langmuir isotherm model with a maximum adsorption capacity of 104.9 mg/g and 120.7 mg/g for RR56 and RB160, respectively. In addition, thermodynamic studies indicated the endothermic adsorption process for both studied dyes. Moreover, DFT calculations were carried out to obtain more insight into the interactions between the dyes and adsorbent. The results showed that the hydrogen bondings and Van der Waals interactions are dominant forces between the two studied dyes and PPy/CMC composite. Furthermore, the interaction energies calculated by DFT confirmed the experimental adsorption data, where PPy/CMC resulted in higher removal of both dyes compared to PPy. The developed nanocomposite showed considerable reusability up to 3 cylces of the batch adsorption process.

Cooperativity in alkane-1,2- and 1,3-polyols: NMR, QTAIM, and IQA study of O─H… OH and C─H… OH bonding interactions

Proton nuclear magnetic resonance chemical shifts and atom-atom interaction energies for alkanepolyols with 1,2-diol and 1,3-diol repeat units, and for their 1:1 pyridine complexes, are computed by density functional theory calculations. In the 1,3-polyols, based on a tG'Gg' repeat unit, the only important intramolecular hydrogen bonding interactions are O─H^… OH. By quantum theory of atoms in molecules analysis of the electron density, unstable bond and ring critical points are found for such interactions in 1,2-polyols with tG'g repeat units, from butane-1,2,3,4-tetrol onwards and in their pyridine complexes from propane-1,2,3-triol onwards. Several features (OH proton shifts and charges, and interaction energies computed by the interacting quantum atoms approach) are used to monitor the dependence of cooperativity on chain length: This is much less regular in 1,2-polyols than in 1,3-polyols and by most criteria has a higher damping factor. Well defined C─H^… OH interactions are found in butane-1,2,3,4-tetrol and higher members of the 1,2-polyol series, as well as in their pyridine complexes: There is no evidence for cooperativity with O─H^… OH bonding. For the 1,2-polyols, there is a tenuous empirical relationship between the existence of a bond critical point for O─H^… OH hydrogen bonding and the interaction energies of competing exchange channels, but the primary/secondary ratio is always less than unity.

On the Interaction between Superatom Al12Be and DNA Nucleobases/Base Pairs: Bonding Nature and Potential Applications in O2 Activation and CO Oxidation

The interaction between quasi-chalcogen superatom Al₁₂Be and DNA nucleobases/base pairs has been explored by searching for the most stable Al₁₂Be-X (X = DNA bases and base pairs) complexes. Our results reveal that Al₁₂Be prefers to combine with guanine by two Al-O and Al-N bonds rather than the other DNA bases, no matter in free state or base pair. The formed Al-N and Al-O bonds between Al₁₂Be and DNA bases proved to be strong polar covalent bonds by the Wiberg bond index, nature bond orbitals, atoms in molecules theory, localized molecular orbitals, and electron localization functions analyses. More importantly, it is found that the formed global minimum of Al₁₂Be-G has the ability to activate an oxygen molecule into a peroxide dianion ¹O₂ ^2-, which can further catalyze the CO oxidation via the Eley-Rideal mechanism with a small energy barrier of 7.78 kcal/mol. We hope that this study could not only provide an in-depth understanding on the intermolecular interaction between metallic superatoms and DNA at the molecular level but also attract more interest in designing and synthesizing superatom-based heterogeneous catalysts with DNA/nucleobases as basic building blocks.

Partitioning a Molecule into the Atomic Basins and the Resultant Atomic Charges from Quantum Chemical Topology Analysis of the Kohn-Sham Potential

Quantum chemical topology (QCT) solidifies the chemical basic concepts demonstrating how a molecular system is intrinsically partitioned into its components and what the interaction lines between them are. Here, QCT analysis using a Kohn-Sham one-electron potential (KSpot) in KS equation as a scalar function is initiated and explored, showing KSpot and its resultant electron force lines have novel spatial features which reveal that an atom in a molecule is a spatial basin governed by its nucleus as a 3D-attractor that terminates all the electron force lines defined by the negative gradient of KSpot and that a chemical bond line is just a minimum path of KSpot for the electron motion. Particularly, the atomic charges from this KSpot QCT analysis are moderate and good, having much lower dependence on basis sets chosen for computation. This may provide a platform for the study of molecular structures and properties, intra- and intermolecular electrostatic interaction, energy decomposition, and construction of force field.

A Cornucopia of Iridium Nitrogen Compounds Produced from Laser-Ablated Iridium Atoms and Dinitrogen

The reaction of laser-ablated iridium atoms with dinitrogen molecules and nitrogen atoms yield several neutral and ionic iridium dinitrogen complexes such as Ir(N₂ ), Ir(N₂ )⁺ , Ir(N₂ )₂ , Ir(N₂ )₂ ^- , IrNNIr, as well as the nitrido complexes IrN, Ir(N)₂ and IrIrN. These reaction products were deposited in solid neon, argon and nitrogen matrices and characterized by their infrared spectra. Assignments of vibrational bands are supported by ab initio and first principle calculations as well as ^14/15 N isotope substitution experiments. The structural and electronic properties of the new dinitrogen and nitrido iridium complexes are discussed. While the formation of the elusive dinitrido complex Ir(N)₂ was observed in a subsequent reaction of IrN with N atoms within the cryogenic solid matrices, the threefold coordinated iridium trinitride Ir(N)₃ could not be observed so far.

Efficient Computation of the Nonlinear Schrödinger Equation with Time-Dependent Coefficients

Motivated by the limited work performed on the development of computational techniques for solving the nonlinear Schrödinger equation with time-dependent coefficients, we develop a modified Runge–Kutta pair with improved periodicity and stability characteristics. Additionally, we develop a modified step size control algorithm, which increases the efficiency of our pair and all other pairs included in the numerical experiments. The numerical results on the nonlinear Schrödinger equation with a periodic solution verified the superiority of the new algorithm in terms of efficiency. The new method also presents a good behaviour of the maximum absolute error and the global norm in time, even after a high number of oscillations.

Theoretical Investigation of Cyano-Chalcogen Dimers and Their Importance in Molecular Recognition

In this manuscript the different noncovalent interactions established between (HYCN)₂ dimers (Y=S, Se and Te) have been studied at the MP2 and CCSD(T) level of theory. Several homodimers have been taken into account, highlighting the capacity of these compounds to act both as electron donor and acceptor. The main properties studied were geometries, binding energy (E_b ), and molecular electrostatic potential (MEP). Given the wide application of chalcogen bonds, and more specifically of cyano-chalcogen moieties in molecular recognition, natural bond orbital (NBO), "atoms-in-molecules" (AIM), and electron density shift (EDS) analysis were also used to analyse the different noncovalent interactions upon complexation. The presence of hydrogen, chalcogen and dipole-dipole interactions was confirmed and their implications on molecular recognition were analysed.