Wiberg Bond Indices Research Articles

Cation-lone pair interaction in Alkali/alkaline earth metal ion-heavier borazine analogue complexes.

The present study is the first report on the formation of alkali/alkaline earth metal ion-heavier borazine analogue complexes via cation-lone pair interaction. Density functional calculations are performed in scrutinizing the complex formation between alkali (Li+, Na+, K+)/alkaline earth (Be2+, Mg2+, Ca2+) metal ions and heavier borazine analogues (HBA) viz. B3P3H6, Al3N3H6, Al3P3H6, Al3As3H6, and Ga3P3-H6. The complexes are found to be stable in gas phase with stabilization energies within the range 26.40 to 324.74 kcalmol-1. The stability can be attributed to the polarizing power of the involved metal ions. Presence of solvent phase exerted notable impact on the stability of the complexes; stability is reduced significantly with the increase in solvent polarity. The process of complexation is exothermic and spontaneous. QTAIM analysis indicated the presence of both ionic and covalent interaction between HBAs and metal ions. HOMO energy, Wiberg bond index, NCI-isosurface and RDG plot analysis revealed the major role of cation-lone pair interaction in the complexation process.

Exploring the activation potential of heme for 2,4-dichlorophenol, 2,4,6-trichlorophenol, and pentachlorophenol

The presence of chlorophenols in water poses a significant threat to human health and the environment. In response to this issue, a study was undertaken to evaluate the catalytic capabilities of chlorinated Heme towards common chlorophenols present in water, such as 2,4-dichlorophenol, 2,4,6-trichlorophenol, and pentachlorophenol. The study employed the B3LYP method, a sophisticated computational technique within density functional theory, to investigate the molecular interactions and transformations involved. It scrutinized structural parameters, Wiberg Bond Indices, which offer insights into the strength and nature of chemical bonds, along with spectroscopic data including infrared vibrational spectra, ultraviolet-visible absorption spectra, and molecular fluorescence spectra. Furthermore, the research analyzed molecular binding energies and orbital energy levels before and after the formation of complexes between Heme and the targeted chlorophenols. The findings indicate that Heme displays a notable activation characteristic towards these chlorophenols. This suggests that Heme could act as an effective catalyst in the degradation of chlorophenols in water, presenting a novel approach to water purification. The theoretical insights derived from this study are invaluable, potentially guiding the development of more efficient catalytic systems for treating chlorophenol-contaminated water, thereby reducing the environmental and health risks associated with these hazardous compounds.

Cooperativity and halonium transfer in the ternary NCI···CH3I···-CN halogen-bonded complex: An ab initio gas phase study.

The strength and nature of the two halogen bonds in the NCI···CH3I···-CN halogen-bonded ternary complex are studied in the gas phase via ab initio calculations. Different indicators of halogen bond strength were employed to examine the interactions including geometries, complexation energies, Natural Bond Order (NBO) Wiberg bond indices, and Atoms in Molecules (AIM)-based charge density topological properties. The results show that the halogen bond is strong and partly covalent in nature when CH3I donates the halogen bond, but weak and noncovalent in nature when CH3I accepts the halogen bond. Significant halogen bond cooperativity emerges in the ternary complex relative to the corresponding heterodimer complexes, NCI···CH3I and CH3I···-CN, respectively. For example, the CCSD(T) complexation energy of the ternary complex (-18.27kcal/mol) is about twice the sum of the complexation energies of the component dimers (-9.54kcal/mol). The halonium transfer reaction that converts the ternary complex into an equivalent one was also investigated. The electronic barrier for the halonium transfer was calculated to be 6.70kcal/mol at the CCSD(T) level. Although the MP2 level underestimates and the MP3 overestimates the barrier, their calculated MP2.5 average barrier (6.44kcal/mol) is close to that of the more robust CCSD(T) level. Insights on the halonium ion transfer reaction was obtained by examining the reaction energy and force profiles along the intrinsic reaction coordinate, IRC. The corresponding evolution of other properties such as bond lengths, Wiberg bond indices, and Mulliken charges provides specific insight on the extent of structural rearrangements and electronic redistribution throughout the entire IRC space. The MP2 method was used for geometry optimizations. Energy calculations were performed using the CCSD(T) method. The aug-cc-pVTZ basis set was employed for all atoms other than iodine for which the aug-cc-pVTZ-PP basis set was used instead.

Why Are Some Pnictogen(III) Pincer Complexes Planar and Others Pyramidal?

Geometrically-constrained pnictogen pincer complexes have emerged in recent years as platforms for unique stoichiometric and catalytic chemical transformations. These complexes feature dynamic conformations ranging from fully planar at the pnictogen centre to distorted-pyramidal geometries, as well as variation between phases. Although the valued reactivity of pnictogen pincer complexes is ascribed to their geometries, there is no unified model to explain the observed conformational outcomes across different ligands and pnictogen centres. Here we propose such a model through computational analysis of more than 1300 structures across 64 complexes (16 ligands and 4 heavy pnictogens), explaining the experimental observations and making new predictions. By looking at signatures of bond stability (bond lengths, Wiberg bond indices) and delocalization (NPA charges, Hirshfeld charges), our framework posits a pnictogen-based σ-bonding effect that favours pyramidalization and exists in competition with a ligand-based π-bonding effect that favours planarity. Variations in structure as a function of pnictogen identity, ligand tethering, electronics, and aromaticity can be reconciled with reference to a balance between these two opposing forces. Careful consideration of the σ/π-bonding effects may aid in the rational design of future pnictogen pincer complexes with predictable geometries and reactivities.

Dative Bonding in Quasimetallatranes Containing Group 15 Donors (Y = N, P, and As) and Group 14 Acceptors (M = Si, Ge, Sn, and Pb).

Metallatranes and their analogous fused ring [3.3.0] bicyclic compounds, quasimetallatranes, have emerged as fascinating molecular systems with intriguing structural, bonding, and conformational properties. We present a comprehensive investigation aimed at unraveling the nature of dative bonding and exploring the conformational flexibility of these compounds. We extensively characterize the dative bond between the metal center and the electron pair donor, using a range of modeling techniques. Our analyses involve structural optimizations, molecular orbital examinations, and covalency ratio calculations, which provide a thorough understanding of the bonding interactions responsible for the stability of these systems. The results confirmed the presence of dative bonds, supported by the close proximity between the metal and the electron-donating group, and the observation of overlapping electron density. Our studies reveal a correlation between the size of the electron-donor and the coordinating metal atom, and the strength of the dative interaction, as indicated by the bond length and the Wiberg bond indices. This bond strength, in turn, influences the conformational preferences adopted by these compounds. This investigation sheds light on the fundamental aspects of the fused ring [3.3.0] bicyclic quasimetallatrane compounds and offers valuable insights into their unique properties.

Study of heavy metals adsorption using a silicate-based material: Experiments and theoretical insights

Heavy metal toxicity in water is a serious problem with harmful effects on human health and the ecosystem. This research studied a silicate-based material as an adsorbent for removing four heavy metals from aqueous solutions. The target metal pollutants selected include manganese (Mn2+), copper (Cu2+), cobalt (Co2+), and zinc (Zn2+). First, theoretical tools including potential energy surface analysis, Natural Population Analysis, AIM, Wiberg Bond Index, QTAIM, and topological methods offer profound insights into the nature of interactions present in the Mg2O8Si3M (M = Mn2+, Cu2+, Co2+, Zn2+) clusters. Second, the synthesis and characterization of eco-friendly hydrated amorphous magnesium silicate (MgOSiO2nH2O) was developed. Last, a simple kinetic adsorption test was applied to assess the material selectivity towards heavy metals and support theoretical results. The kinetic adsorption study was analyzed through the pseudo-first and second-order kinetics, Elovich, and the intraparticle diffusion models. The theoretical analysis of the adsorption energies indicates that the adsorption of four metal ions on the Mg2O8Si3 surface is energetically favorable in all cases. The material displayed the following adsorption sequence: Cu2+ (59 mg g-1) > Zn2+(25 mg g-1) ≈ Co2+ (23 mg g-1) > Mn2+ (15 mg g-1). This knowledge can then be used to design and optimize low-cost silicate-based materials for effective heavy metal removal, contributing to efforts to address environmental pollution and protect public health.

Catalytic Activity of an Asymmetric Uranyl‐Salophen on 2,4‐Dichlorophenol, 2,4,6‐Trichlorophenol, and Pentachlorophenol

ABSTRACTThis study employs density functional theory to investigate whether asymmetric uranyl‐Salophen catalyzes the activation of 2,4‐dichlorophenol, 2,4,6‐trichlorophenol, and pentachlorophenol. Changes in various important parameters such as bond lengths, Wiberg bond indices, infrared vibrational absorption spectra, natural charge distributions, ultraviolet–visible absorption spectra, and 13C NMR chemical shifts of the aromatic rings of chlorophenols before and after forming complexes with asymmetric uranyl‐Salophen are analyzed. This means that asymmetric uranyl‐Salophen indeed activates 2,4‐dichlorophenol, 2,4, 6‐trichlorophenol, and pentachlorophenol, significantly lengthening the aromatic ring bonds of the three chlorophenols and increasing the electrophilicity of the carbon atoms in the aromatic ring. This suggests that the conjugated system of the aromatic ring is also partially disrupted, indicating that these three chlorophenols are indeed catalytically activated. Furthermore, this implies that this asymmetric uranyl‐Salophen may be used to degrade chlorophenols in water.

Predicting Small Molecule Activations Including Dinitrogen Based on an Inorganic Benzene B4N2 Framework.

Although main group species have emerged in the field of dinitrogen activation in recent years, the reported examples are particularly rare in comparison with transition metal complexes due to their significant challenges. Herein, we demonstrate a [4 + 2] cycloaddition reaction of N2 (with an activation energy as low as 12.5 kcal mol-1) initiated by an inorganic benzene via density functional theory calculations. Such N2 activation is supported by the elongated nitrogen-nitrogen bond distance (dNN), decreased vibration frequency (νNN), and weakened Wiberg bond index (WBINN). Subsequently, the "push-pull" electronic effect, formed by introducing a Lewis acid, HB(C6F5)2, facilitates the generation of thermodynamically more stable products. In addition, this inorganic benzene could also be used to activate a series of small molecules, including carbon dioxide, acetylene, ethylene, and acetonitrile with reaction barriers ranging from 4.7 to 11.6 kcal mol-1. Our findings provide an alternative approach to N2 activation and functionalization, theoretically validating the feasibility of the dual Lewis acid strategy for dinitrogen activation.

β-d-Glucopyranose – Silver+ (1:1) complex assisted aromaticity involving electron deficient X3 (X=Be, Mg) core

Seeing the frontier molecular orbitals (FMO) picture of β-d-glucopyranose–silver+ (1:1) complex ([Ag(C6H12O6)]+) reported in a recent article (Mondal, 2024) [53], we used the FMO based approach to obtain inverted sandwich like organometallic complex comprising X3 (X=Be, Mg) core, in silico. Thermodynamic feasibility of obtaining [(C6H12O6)Ag-X3-Ag(C6H12O6)] (X=Be, Mg) complexes are explored. Conceptual density functional theory (CDFT) based reactivity descriptors guide to ascertain the stability order of such complexes. By the help of natural bond orbital analysis, Wiberg bond indices (WBI), bond critical point analysis, electron density descriptors, and gradient isosurface between selected contacts, the nature of interaction between the X3 core and Ag metals are studied. FMO of the complexes further help in ascertaining the nature of binding. Energy decomposition analysis marks the interaction between X3 and Ag metals in [(C6H12O6)Ag-X3-Ag(C6H12O6)] (X=Be, Mg) as partially covalent. Negative values of nucleus independent chemical shift (NICS) at the triangular ring (X3) centres [NICS(0)] indicate aromatic property.

Solvent effects on the structures of the hydrated copper dication clusters

Understanding the solvation process requires an understanding of how the solvent affects the molecular cluster configurations. The current study examined the impact of solvents on the structures of the hydrated Cu2+ cation, as well as the temperature dependence, and the energetics of the hydrated Cu2+ cation clusters from sizes n=1 from n=10 (Cu2+(H2O)n=1−10). The potential energy surfaces (PESs) in the solvent phase (represented by a dielectric continuum medium) at the MP2/6-31++G(d,p) level of computation have been comprehensively investigated. We employed the integral equation formalism polarized continuum model (IEF-PCM) to describe the dielectric continuum medium. It is pointed out that changes in PESs of hydrated copper dication isomers and structural morphologies occur in the continuum medium. Thus, the relative energies are smallest at T=0 K and for the cluster populations of Cu2+(H2O)n=1−10 are modified. Thus, in the solvent phase, the competition between the conformers is more severe than in the gas phase. Furthermore, highly coordinated conformers are more favorable, and favored conformers in the gas phase are disfavored in the solvent phase. The structural investigations are validated by the analysis of Wiberg Bond Index (WBI) at 300 K. In the solvent phase, the WBIs of the most stable isomers are the smallest and the WBIs of axial bonds are lower than those of equatorial bonds. Concerning the energy analysis, the binding energies were predicted for enormous sizes (at saturation) using a fit function. Thus, the new values of the binding energies obtained are -5467.8 and ▪ in the gas and solvent phases respectively. All these observed differences are due to the effects of the implicit solvent.

First principles study on interactions in inorganic molecular crystals at zero dimensions

Inorganic molecular crystals (IMCs) are attracting significant scientific interest due to their distinctive anisotropic crystal structure. The molecular units in the form of cages are the fundamental unit of these materials and most of their properties are dependent on the geometry, their arrangement and the interactions present in between these cages. This work is carried out to shed light on the structural arrangement and the interaction properties of three IMCs including antimony oxide (Sb2O3), phosphorous tri-selenium (P4Se3) and phosphorous trioxide (P4O6) in zero dimensions. The reported calculations were carried out using density Functional Theory (DFT) to investigate the interactions between the cages based on ETS-NOCV, NBO, Wiberg bond indices and QTAIM analysis. The results revealed weak VdWs association along the inter-cages and significant covalent bonding within the cage. The investigation shows that Sb2O3 and P4Se3 IMCs have potential applications in high-pressure and stress-bearing materials.

Palladium(II)-catalyzed annulation of N-methoxy amides and arynes: computational mechanistic insights and substituents effects.

The combined use of transition metal-catalyzed C-H activation with aryne annulation reactions has emerged as an important strategy in organic synthesis. In this study, the mechanisms of the palladium(II)-catalyzed annulation reaction of N-methoxy amides and arynes were computationally investigated by density functional theory. The role of methoxy amide as a directing group was elucidated through the calculation of three different pathways for the C-H activation step, showing that the pathway where amide nitrogen acts as a directing group is preferable. At the reductive elimination transition state, an unstable seven-membered ring is formed preventing the lactam formation. A substituent effect study based on an NBO analysis, Hammet, and using a More O'Ferall-Jenks plot indicates that the C-H activation step proceeds via an electrophilic concerted metalation-deprotonation (eCMD) mechanism. The results show that electron-withdrawing groups increase the activation barrier and contribute to an early Pd-C bond formation and a late C-H bond breaking when compared with electron-donating substituents. Our computational results are in agreement with the experimental data provided in the literature. All calculations were performed using Gaussian 16 software. Geometry optimizations, frequency analyses at 393.15 K, and IRC calculations were conducted at the M06L/Def2-SVP level of theory. Corrected electronic energies, NBO charges, and Wiberg bond indexes were computed at the M06L/Def2-TZVP//M06L/Def2-SVP level of theory. Implicit solvent effects were considered in all calculations using the SMD model, with acetonitrile employed as the solvent.

Unveiling the structural and bonding properties of AuSi2- and AuSi3- clusters: A comprehensive analysis of anion photoelectron spectroscopy and abinitio calculations.

Silicon clusters infused with transition metals, notably gold, exhibit distinct characteristics crucial for advancing microelectronics, catalysts, and energy storage technologies. This investigation delves into the structural and bonding attributes of gold-infused silicon clusters, specifically AuSi2- and AuSi3-. Utilizing anion photoelectron spectroscopy and abinitio computations, we explored the most stable isomers of these clusters. The analysis incorporated Natural Population Analysis, electron localization function, molecular orbital diagrams, adaptive natural density partitioning, and Wiberg bond index for a comprehensive bond assessment. Our discoveries reveal that cyclic configurations with the Au atom atop the Si-Si linkage within the fundamental Si2 and Si3 clusters offer the most energetically favorable structures for AuSi2- and AuSi3- anions, alongside their neutral counterparts. These anions exhibit notable highest occupied molecular orbital-lowest unoccupied molecular orbital gaps and significant σ and π bonding patterns, contributing to their chemical stability. Furthermore, AuSi2- demonstrates π aromaticity, while AuSi3- showcases a distinctive blend of σ antiaromaticity and π aromaticity, crucial for their structural robustness. These revelations expand our comprehension of gold-infused silicon clusters, laying a theoretical groundwork for their potential applications in high-performance solar cells and advanced functional materials.

Size Matters: Computational Insights into the Crowning of Noble Gas Trioxides.

In pursuit of enhancing the stability of the highly explosive and shock-sensitive compound XeO3, we performed quantum chemical calculations to investigate its possible complexation with electron-rich crown ethers, including 9-crown-3, 12-crown-4, 15-crown-5, 18-crown-6, and 21-crown-7, as well as their thio analogues. Furthermore, we expanded our study to other noble gas trioxides (NgO3), namely, KrO3 and ArO3. The basis set superposition error (BSSE) corrected interaction energies for these adducts range from -13.0 kcal/mol to -48.2 kcal/mol, which is notably high for σ-hole-mediated noncovalent interactions. The formation of these adducts was observed to be more favorable with the increase in the ring size of the crowns and less favorable while going from XeO3 to ArO3. A comprehensive analysis by various computational tools such as the mapping of the electrostatic potential (ESP), Wiberg bond indices (WBIs), Bader's theory of atoms-in-molecules (AIM), natural bond orbital (NBO) analysis, noncovalent interaction (NCI) plots, and energy decomposition analysis (EDA) revealed that the C-H···O interactions, as well as dispersion interactions, play a pivotal role in stabilizing adducts involving larger crowns. A noteworthy outcome of our study is the revelation of a coordination number of 9 for xenon in the complex formed between XeO3 and the thio analogue of 18-crown-6, which is higher than the largest number reported to date.

Exploring Unconventional σ-Hole Interactions: Computational Insights into the Interaction of XeO3 with Non-Aromatic Coordinating Solvents.

In order to control the explosiveness and shock sensitivity of XeO3 , we have investigated its plausible interaction with various non-aromatic coordinating solvents, serving as potential Lewis base donors, through density functional theory (DFT) calculations. Out of twenty six such solvents, the top ten were thus identified and then thoroughly examined by employing various computational tools such as the mapping of the electrostatic potential surface (MESP), Wiberg bond indices (WBIs), non-covalent interaction (NCI) plots, Bader's theory of atoms-in-molecules (AIM), natural bond orbital (NBO) analysis, and the energy decomposition analysis (EDA). The amphoteric nature of XeO3 was also explored by investigating the extent of back donation from the lone pair of Xe to the antibonding orbital of the donating atom/group of the solvent molecules. The C-H…O interactions were also found to be a contributing factor in the stabilization of these adducts. Although these aerogen-bonding interactions were found to be predominantly electrostatic, significant contributions from the orbital contributions, as well as dispersion interactions, were observed. The top three non-aromatic solvents (among the twenty six studied) which form the strongest adducts with XeO3 are proposed to be hexamethylphosphoramide (HMPA), N,N'-dimethylpropyleneurea (DMPU) and tetramethylethylenediamine (TMEDA).

Near equivalence of polarizability and bond order flux metrics for describing covalent bond rearrangements.

Identification of the breaking point for the chemical bond is essential for our understanding of chemical reactivity. The current consensus is that a point of maximal electron delocalization along the bonding axis separates the different bonding regimes of reactants and products. This maximum transition point has been investigated previously through the total position spread and the bond-parallel components of the static polarizability tensor for describing covalent bond breaking. In this paper, we report that the first-order change of the Wiberg and Mayer bond index with respect to the reaction coordinate, the bond flux, is similarly maximized and is nearly equivalent with the bond breaking points determined by the bond-parallel polarizability. We investigate the similarites and differences between the two bonding metrics for breaking the nitrogen triple bond, twisting around the ethene double bond, and a set of prototypical reactions in the hydrogen combustion reaction network. The Wiberg-Mayer bond flux provides a simpler approach to calculating the point of bond dissociation and formation and can yield greater chemical insight through bond specific information for certain reactions where multiple bond changes are operative.

A theoretical study of the oxidation of benzene by manganese oxide clusters: formation of quinone intermediates.

Manganese oxides (MnxOy) have been widely applied in various chemical industrial processes owing to their long lifetime, low cost and high abundance. They have been used as co-reactants for the elimination of volatile organic compounds (VOCs); however, their oxidation mechanism is not clearly established. In this theoretical study, interaction capacities between benzene (C6H6) and MnxOy clusters, which were modeled with MnO2 and Mn2O3 molecules, were investigated by quantum chemical computations using density functional theory (DFT) with the PBE-D3 functional. The interaction capacity between C6H6 and MnxOy was evaluated, and the probing of the initial stage of the C6H6 oxidation at a molecular level offers an in-depth oxidation reaction path. Interaction energies computed in several spin states, along with the analysis of the electron distribution using the quantum theory of atoms in molecules, natural bond orbital and Wiberg bond index techniques as well as local softness values and MO energies of fragments, point out that the interaction between C6H6 and Mn2O3 is stronger than that with MnO2, amounting to -43 and -35 kcal mol-1, respectively, and the metal atom is identified as the primary active site. During the oxide cluster-assisted oxidation, benzene firstly undergoes an oxidation reaction by active oxygen to generate intermediates such as hydroquinone and benzoquinone. The pathway involving p-benzoquinone as the product (noted as PR1) is the most energetically favored one through a transition structure lying at 19 kcal mol-1, below the energy reference of the reactants, leading to an energy barrier significantly lower than that of 36 kcal mol-1 found for the gas phase oxidation reaction with molecular oxygen without the assistance of the oxide clusters. Potential energy profiles illustrating the reaction paths and molecular mechanisms were described in detail.

Chemiluminescence of a Firefly Luciferin Analogue Reveals that Formation of the Key Intermediate Responsible for Excited State Generation Occurs on a Fully Concerted Step.

The chemiluminescence (CL) reaction of eight different 2-(4-hydroxyphenyl)-4,5-dihydrothiazole-4-carboxylate esters with an organic superbase and oxygen was investigated through a kinetic and computational study. These esters are all analogues to the luciferin substrate involved in efficient firefly bioluminescence. The kinetic data obtained from CL emission and light absorption assays were used in the context of linear free energy relationships (LFER); we obtained the Hammett reaction constant ρ = +1.62 ± 0.09 and the Brønsted constant βlg = -0.39 ± 0.04. These observations from LFER, together with activation parameters obtained from Arrhenius plots, suggest that the formation of the high-energy intermediate (HEI) 1,2-dioxetanone occurs via a concerted mechanism during the rate-determining step of the reaction. Calculations performed using density functional theory support a late transition state for HEI formation within the reaction mechanism pathway, which was described considering geometric parameters, Wiberg bond indices from natural bond order analysis, and the atomic charges derived from the electrostatic potential.

Theoretical insights into the extraction behaviors of neptunium(VI) and plutonium(IV) based on the structure of organophosphorus ligands

Organophosphorus ligands such as TBP, TiAP, and DMHMP have exhibited excellent performance in recovering actinides from spent fuel. In this work, the molecular geometries and properties of TBP, TiAP and DMHMP were investigated using density functional theory calculations. Furthermore, the extraction mechanism of ligands for actinides (Np(Ⅵ), Pu(Ⅳ)) was further elucidated by simulating the microstructures and extraction reactions of metal–ligand complexes. The results demonstrate that the complexation ability of the three ligands on actinide cations (NpO22+ and Pu4+) follows the order of DMHMP > TiAP > TBP. The electrostatic potential (ESP) analysis indicates that the nucleophilic ability of DMHMP is stronger than that of the other two ligands. The frontier molecular orbital analysis of the three ligands represents that DMHMP has the highest HOMO energy, suggesting that it has the strongest electron-donating capability and is more likely to bond with metal ions. The values of Wiberg bond indices (WBI) suggest that the MO bonds in DMHMP complexes have more covalency. According to the QTAIM analysis, the interactions between actinide cations and the ligands are predominantly ionic in nature. The molecular orbital analysis of the complexes shows that the M(NO3)n·2DMHMP (MNpO22+ and Pu4+) complexes are more stable, which is supported by thermodynamic energy analysis. This work has clarified the complexing properties of actinide cations with three ligands, shedding light on the extraction mechanisms of organophosphorus ligands for actinide cations. It is anticipated to lay the theoretical foundation for the efficient recovery of critical actinide elements in spent fuel reprocessing, which will also provide innovative approaches for the design and development of related separation processes.

Utility of (MgO)12 nanocage as a chemical sensor for recognition of amphetamine drug: A computational inspection

DFT calculations on sensor-drug interactions are necessary for understanding binding mechanisms, predicting sensor performance, evaluating stability and reactivity, and rational design of sensor materials. We scrutinized the adsorption of amphetamine (AFE) on the pure magnesium oxide nano-cage (MgONC) by applying density functional theory. All geometries and single point energy computations were optimized at M06–2X/6–311 G (d, p). Furthermore, we performed an analysis of the natural bond orbital (NBO) and evaluated the values of partial natural charges. Additionally, we investigated donor-acceptor (D-A) interactions and examined the Wiberg bond index (WBI) in greater depth. The MgONC was capable of adsorbing AFE with greater strength with the energy of adsorption (Eads) of −48.19 kcal/mol (for stable configurations). Moreover, the NBO method demonstrated more effective D-A interactions between AFE and the MgONC. Based on the computations, for the most stable configuration, there was a substantial alteration in the HOMO-LUMO gap of the MgONC following the drug adsorption, thus increasing the electrical conductance (EC) of the MgONC. The sensing mechanism is related to the gap difference, which depends on the change in the EC. We adopted the conventional transition state theory for the prediction of recovery time. The computations indicated that the MgONC+ AFE configuration had a short recovery time for the desorption of AFE. Finally, based on our findings, we could conclude that the MgONC is an appropriate choice for the improvement of effective AFE sensors. DFT study of drug sensors will focus on enhancing sensitivity, selectivity, and stability while exploring novel materials and optimizing performance through theoretical simulations and analysis.