Density Functional Theory Studies on Boron Nitride and Silicon Carbide Nanoclusters Functionalized with Amino Acids for Organophosphorus Pesticide Adsorption

: This study compares the properties of B 12 N 12 and Si 12 C 12 nanoclusters functionalized with tyrosine in the adsorption of organophosphorus pesticides, focusing on adsorption energy and electronic stability. The results indicate that B 12 N 12 /tyrosine exhibits more negative adsorption energies than Si 12 C 12 /tyrosine, suggesting stronger interactions and higher adsorption stability. Additionally, B 12 N 12 demonstrates higher ionization energy and chemical hardness, enhancing its electronic stability during the adsorption process. In contrast, Si 12 C 12 has higher electrophilicity and maximum electron transfer capacity, leading to greater variability in adsorption energy and more ﬂ exible electronic structure adjustments. These ﬁ ndings suggest that B 12 N 12 nanoclusters have greater potential and application value as adsorption materials, particularly when modi ﬁ ed with tyrosine. B 12 N 12 /tyrosine demonstrates higher stability and predictability in pesticide adsorption, making it more suitable for related applications.


Introduction
Nanozymes exhibit excellent stability and enhanced catalytic activity, providing a solid foundation for their applications in environmental pollutant monitoring and remediation [1].The use of nanomaterials as novel carriers for enzyme delivery and the regulation of enzyme activity has sparked significant interest in the field of nanobiotechnology for biomedical applications [2].Currently, nanostructure-based drug delivery systems represent one of the most compelling areas of research.Among these, boron nitride nanoclusters show promise as drug carriers for targeted delivery systems.
Previous studies have investigated the interaction mechanisms between the drug Anagrelide (AG) and B12N12, as well as Al-and Ga-doped B12N12 nanocages in both gaseous and aqueous media using DFT with the B3LYP/6-31G(d,p) method [3].Additionally, a comprehensive computational study was conducted using density functional theory (DFT) to investigate the adsorption of the drug molecule penicillamine (PCA) on Al-and Ga-doped B12N12 nanoclusters in aqueous and chloroform environments.The PCA molecule interacts effectively with Al-and Ga-doped clusters through its four nucleophilic sites: amine, carbonyl, hydroxyl, and thiol groups.The computational results indicate that the most stable adsorption complex is achieved when the PCA molecule is adsorbed through its amine group.The adsorption of the PCA molecule reduces the HOMO-LUMO gap and overall hardness of the doped clusters, confirming an increase in the reactivity of the considered clusters for drug delivery purposes [4].
Based on Celaya et al.'s research findings, boron nitride and some of its mixed nanoclusters have the potential to serve as nanocarriers for the targeted delivery of the drug melphalan [5].Using density functional theory (DFT) calculations, the interactions between Cu/Zn-doped boron nitride nanocages and the anticancer drug mercaptopurine (MP) were investigated, leading to the design of a novel drug carrier.The adsorption of MP on Cu/Zn-doped boron nitride nanocages was found to be appropriate [6].Hachem et al. employed density functional theory (DFT) calculations using the PBE1 functional to investigate the structural, electronic, and spectral properties of B16N16, B15GeN16, and B15SiN16 nanocages modified with sulfasalazine (SSZ).When SSZ is adsorbed onto the B16N16, B15GeN16, and B15SiN16 nanocages through its pyridine ring, the system attains its most stable state.Furthermore, the adsorption of SSZ via the SO2 group on the B16N16, B16GeN16, and B16SiN16 nanocages results in lower binding energy and an increased dipole moment, both of which can enhance the sensitivity of the B16GeN16 nanocage to the drug [7].Recently, boron nitride (BN) nanoparticles have been identified as promising candidates for drug delivery systems.Cao et al.'s studies utilized density functional theory (DFT) and time-dependent density functional theory (TDDFT) models with B3LYP and B3PW91 functionals to investigate the loading of curcumin (CUR) in various potential states onto platinum-functionalized B12N12 nanocages, aiming to enhance its solubility and stability.The strong binding of CUR to Pt-B12N12 can be attributed to the substantial charge transfer from CUR to the cage, resulting in significant changes in the dipole moment and energy gap [8].
In Hasan et al.'s studies, the adsorption mechanisms of the Favipiravir (FPV) molecule on the pristine, Zn-functionalized, and Ni-functionalized B12N12 (BN, Zn f-BN, and Ni f-BN) nanocages' exterior surfaces were investigated using DFT/QTAIM methods and the B3LYP/6-31G(d,p) method.Adsorption energy data indicate that the functionalized BN adsorbents can very effectively adsorb the FPV drug compared to the pristine adsorbents.The reduction in the HOMO-LUMO gap by up to 67.79% suggests that the drug can be promisingly detected through the generated electrical signal in the case of f-BN nanocages [9].Hossain et al.'s research using the DFT B3LYP/6-31G (d, p) level of theory revealed that the Ni-B12N12 nanocage is the most prominent nanostructure for forming drug delivery complexes to transport chlorocholine (CM) drugs.This nanocage exhibits the most stable adsorption of CM drugs in both gas and aqueous media.The interaction between the CM drug and the nanostructure was also confirmed through frontier molecular orbital and QTAIM analyses [10].
Soltani et al.'s studies reported the adsorption of neutral and zwitterionic forms of serine on pristine and Pt-modified B12N12 fullerenes using density functional theory (DFT) and time-dependent density functional theory (TD-DFT) calculations.The binding energies of serine on the fullerene surface were investigated through its hydroxyl (-OH), carboxyl (-COOH), and amine (-NH2) functional groups.According to the analysis using the M06-2X functional, the binding energy of the zwitterionic form of serine on the B12N12 fullerene is less stable than that of the neutral form.The results indicate that the most stable chemisorption state of serine occurs through its amino group interacting with the Pt-modified B12N12 fullerene.Additionally, the conductivity of the B12N12 and Pt-modified B12N12 fullerenes is influenced by band gap changes when serine is adsorbed on their external surfaces [11].
In Feng et al.'s studies, computational simulation methods were employed to investigate the removal of organic pollutants from wastewater.Researchers selected two types of nanoporous boron nitride (BN) materials, B16N16 and B12N12, as adsorbents.They conducted a comprehensive evaluation of pure boron nitride adsorbents and glycine-functionalized boron nitride adsorbents (B16N16(Glycine)2/B12N12(Glycine)2) to understand the mechanisms and adsorption capacities for removing two different emerging pollutants.Particularly, the study examined the adsorption capacity and potential changes in spectral and electronic properties (including frontier molecular orbitals, energy gap (ΔEGAP), chemical softness (σ), and hardness (η)) of the pure and functionalized boron nitride adsorbents before and after the adsorption process.The results indicated that the functionalization of nanospheres of boron nitride with glycine enhances the affinity for pollutant adsorption [12].
BN-based materials exhibit satisfactory adsorption capacities for inorganic pollutants, such as heavy metal ions, and organic pollutants, including dyes and pharmaceutical molecules.The interaction mechanisms between pollutants and BN-based materials primarily involve surface complexation, π-π stacking, and electrostatic interactions [13].Onsori and Alipour's research utilizing density functional theory (DFT) calculations explored the reactivity and electronic sensitivity of synthesized B12N12 nanoclusters towards the anticancer drug cisplatin (CP).Upon adsorption of the CP drug, the conduction level of the BN nanoclusters stabilizes significantly, and the valence state shifts to higher energy levels.As a result, the HOMO-LUMO gap decreases notably.Consequently, the BN nanoclusters transform into semiconductors with higher conductivity following the adsorption process.The increase in conductivity can generate electrical signals that facilitate the detection of CP drugs.Furthermore, UV-vis calculations indicate the emergence of a strong peak in the visible region after CP drug adsorption, which aids in drug detection [14].Kamali et al.'s research aimed to investigate the adsorption of the drug metformin (MF) on fullerene-like nanoclusters of boron nitride (B24N24), aluminum nitride (Al24N24), aluminum phosphide (Al24P24), and boron phosphide (B24P24) in both gas and solvent (water) phases using the B3PW91/6-311G(d,p) theoretical level.A comparison of the HOMO-LUMO energy gaps in the pure nanoclusters indicated that the electronic properties of B24N24 nanoclusters can be enhanced by doping with aluminum and phosphorus atoms, compared to the other nanoclusters studied.A thermodynamic analysis revealed that the interaction between MF and the nanoclusters is exothermic and spontaneous [15].
Zhiani employed quantum mechanical methods to investigate the adsorption and binding characteristics of five different amino acids-namely alanine (Ala), arginine (Arg), asparagine (Asn), histidine (His), and cysteine (Cys)-on the surfaces of graphene (Gra) and boron nitride (BN) nanosheets from a molecular perspective.Density functional theory (DFT) and DFT-D3 calculations were used to study the electronic properties and dispersion interactions of the amino acid/adsorbent complexes.The calculations indicated that the polarity of BN nanosheets provides a stronger affinity for the amino acids [16].In Ghasemi et al.'s studies, density functional theory (DFT) calculations were conducted using B3LYP-D and PW91-D functionals to examine the adsorption behavior and detection of metformin on the outer surface of both pure and nitrogen-doped boron nitride (BN) fullerenes.The results indicated that the NH group of the metformin single bond can chemisorb onto the boron atoms of B12N12 and B16N16 fullerenes.The GeB11N12 biosensor demonstrates significant potential for the determination of metformin in environmental systems [17].
Wu et al.'s theoretical studies have demonstrated that incorporating carbon into (BN)12 fullerene makes the hydrogenation reaction on carbon-doped B11N12C clusters thermodynamically favorable and kinetically feasible under ambient conditions.Without the use of metal catalysts, carbon atoms can serve as active centers for dissociating H2 molecules and provide free hydrogen atoms for further hydrogenation on B11N12C fullerene, thereby reducing material costs in practical hydrogen storage applications [18].Esrafili's studies utilizing density functional theory (DFT) investigated the adsorption and catalytic decomposition of N2O molecules on boron nitride nanocages resembling fullerenes (B12N12).The research revealed that the electron-donating properties of the nanocage play a crucial role in the adsorption and activation of N2O.By incorporating carbon atoms into the B12N12 cluster, the results demonstrated that B11N12C or B12N11C exhibits stronger adsorption of N2O compared to the original B12N12 [19].Yoosefian et al. calculated the drugloading capacity of carboxylated carbon nanotubes functionalized with benzylhydrazine for use as nanocarriers for levodopa.In this context, it is crucial to evaluate all adsorption characteristics of the most stable conformer of benzylhydrazine molecules on carboxylated carbon nanotubes.To determine the minimum energy conformer of benzylhydrazine, first-principles quantum mechanical calculations were conducted on the molecular structure, and a conformational analysis of 512 possible isomers was performed.A novel and easily prepared formulation of benzylhydrazine/carboxylated carbon nanotube conjugates was developed, demonstrating high drug-loading efficiency for levodopa, intended for the treatment of Parkinson's disease [20].
This study aims to investigate the stability and electronic structural properties of B12N12 and Si12C12 nanoclusters functionalized with amino acids during the adsorption of organophosphorus (OP) pesticides.By analyzing adsorption energy, ionization energy, chemical hardness, electronegativity, and other electronic structural parameters, we seek to elucidate the superior performance of B12N12 nanoclusters with tyrosine in the adsorption process and compare it with the stability and electronic structural adaptability of Si12C12 nanoclusters with tyrosine.Specifically, this research employs systematic density functional theory (DFT) calculations to explore the energy changes and electronic characteristics of these nanoclusters when adsorbing OP pesticides, thereby determining the potential application value of B12N12 nanoclusters functionalized with amino acids in practical scenarios.
The generalized gradient approximation (GGA) using the Perdew-Burke-Ernzerhof (PBE) functional [23] was selected to compute the exchange-correlation functional of the Hamiltonian operator, utilizing a double numerical plus D-functions (DND) basis set.For adsorption complexes, the basis set incorporated the relativistic effective core potential (ECP) with small-core pseudopotentials.Furthermore, the conductor-like screening model (COSMO) [24,25] was employed to account for solvent effects under aqueous conditions, with the dielectric constant ε set to 78.54, corresponding to water.The self-consistent field (SCF) convergence criteria were established at 1.0 × 10 −5 Ha for energy and 1.0 × 10 −3 Ha/Å for forces.A smearing value of 0.005 Ha was applied for electronic occupation to facilitate convergence.The total energy for each optimized structure was calculated, followed by a frontier molecular orbital (FMO) analysis to determine the HOMO and LUMO energies.The ionization potential and electron affinity calculated using the HOMO and LUMO with the PBE functional are only qualitatively accurate.Further refinement and additional computational methods may be necessary for quantitatively precise results.
Conceptual density functional theory (CDFT) [26,27]  The AutoDock Vina Extended SAMSON Extension was used to dock the organophosphorus pesticides with the AChE enzyme.In this docking procedure, the search domain was based on the entire size of the AChE enzyme.The grid's center was set to ensure the grid box covered the entire enzyme, with default values used for other parameters.From the first 200 docking poses, the best candidate was selected according to the standard scoring function (Figure 1a,b) [28].The molecular modeling involved combining Tyr124 and Tyr337 from the crystallographic structure of AChE with eight different OPs.After docking AChE with the OPs, the structures of Tyr124_OP and Tyr337_OP, which form hydrogen bonds with OP, were extracted and used as the starting structures for modification with nanoclusters, followed by structural optimization (Figures 2 and 3).The distance between the hydrogen atom on the amino group of tyrosine (Tyr337) and the nitrogen atom of the B12N12 nanocluster ranged from 2.8 to 3.5 Å. Functionalizing B12N12 (boron nitride fullerene) with the amino acid tyrosine is particularly effective due to several chemical characteristics.B12N12's high surface reactivity and the electron deficiency of its boron atoms make it highly receptive to nucleophilic attack by tyrosine's phenolic hydroxyl group.The ability of tyrosine to form hydrogen bonds with the nitrogen atoms in B12N12 and engage in π-π interactions with the boron nitride rings further stabilizes the functionalization.Additionally, the polar nature of both tyrosine and B12N12 allows for favorable dipole-dipole interactions, while their size compatibility ensures effective bonding.The total energy of the optimized structures of the B12N12 and Si12C12 nanoclusters functionalized with Tyr124_OP and Tyr337_OP, after subtracting the energy of the nanocluster_Tyr124 (or Tyr337) and OP, represents the reaction energy of the biomimetic nanozyme adsorption of OP at the hydrogen bonding sites of AChE and OP.This adsorption reaction is described as follows: nanocluster_Tyr124 (or Tyr337) + OP → nanocluster_Tyr124 (or Tyr337)_OP The energy change in the adsorption reaction, ΔE, can be expressed as follows: ΔE = E nanocluster_Tyr124 (or Tyr337)_OP − (E nanocluster_Tyr124 (or Tyr337) + E OP)

Adsorption Energies
Table 1 shows the OP pesticide adsorption energies of the B12N12 and Si12C12 nanoclusters functionalized with Tyr124 and Tyr337.The adsorption energies for the B12N12 nanoclusters with Tyr337 are generally more negative, indicating higher adsorption stability, whereas the adsorption energy variation for Si12C12 is larger, indicating less consistency.These results suggest that B12N12 nanoclusters have greater potential and application value as adsorption materials, particularly when functionalized with Tyr337, due to their higher adsorption stability and consistency.Most of the organophosphorus (OP) pesticides have negative adsorption energies with B12N12_Tyr337, indicating an exothermic and stable adsorption process.For example, the OP pesticide Dia has an adsorption energy of −34.01 eV, Isa −34.63 eV, and Tri −34.70 eV.One pesticide, in particular, shows notably large negative adsorption energy with B12N12_Tyr337: the OP pesticide Sul has an adsorption energy of −34.93 eV, indicating very high stability.The energy variation for different pesticides adsorbed with B12N12_Tyr337 is minimal, indicating consistent adsorption capability across different OP pesticides.The adsorption energy varies from approximately −26.89 eV (OP pesticide Nal) to −34.93 eV (OP pesticide Sul), suggesting stable and consistent adsorption processes.
In contrast, for the Si12C12 nanoclusters functionalized with Tyr337, the adsorption energy varies widely, with some positive and some negative values, indicating less consistency during the adsorption process.For example, the OP pesticide Par has a positive adsorption energy of 34.13 eV, indicating an endothermic and unstable adsorption process, while the OP pesticide Dia has a negative adsorption energy of −11.05 eV.The OP pesticide Sul has an adsorption energy of −8.79 eV, which is more stable than Par but less stable than Tri, which has an adsorption energy of −9.59 eV.The pesticide with the maximum stability on Si12C12_Tyr337 is Chl, with an adsorption energy of −48.39 eV.Therefore, the adsorption energies of the OP pesticides on the B12N12 nanoclusters with Tyr337 are generally more negative than those of the Si12C12 nanoclusters with Tyr337, indicating higher stability.Most of the OP pesticides exhibit negative adsorption energies on the B12N12 nanoclusters with Tyr337, while several have positive energies on the Si12C12 nanoclusters with Tyr337.

HSAB Parameters
The ionization energy (I) of B12N12_Tyr337 (5.395 eV) is slightly higher than that of Si12C12_Tyr337 (5.208 eV), suggesting that B12N12_Tyr337 has greater electronic stability (Table 2).This aligns with the more stable adsorption energies observed for B12N12_Tyr337.Additionally, the chemical hardness of B12N12_Tyr337 (1.874 eV) is higher than that of Si12C12_Tyr337 (1.064 eV), indicating a greater resistance to electron transfer.Correspondingly, B12N12_Tyr337 exhibits more negative adsorption energies when adsorbing pesticides, reflecting its more stable electronic structure.The higher electrophilicity index (ω) of Si12C12_Tyr337 (8.070 eV) compared to B12N12_Tyr337 (3.310 eV) implies that Si12C12_Tyr337 more readily accepts electrons during adsorption.Consequently, Si12C12_Tyr337 shows greater variability in adsorption energy.The maximum amount of electronic charge transfer (ΔNmax) indicates the maximum amount of electronic charge a molecule can accept during a reaction.Si12C12_Tyr337's ΔNmax (3.895) is significantly higher than that of B12N12_Tyr337 (1.880), suggesting that Si12C12_Tyr337 can accommodate more electrons (Table 2).In summary, the B12N12_Tyr337 nanoclusters, due to their higher ionization energy and chemical hardness, demonstrate more stable adsorption energies, indicating greater electronic stability during adsorption.In contrast, Si12C12_Tyr337, with higher electrophilicity and maximum electron transfer, shows more variability in adsorption energy, indicating greater flexibility in adjusting its electronic structure during adsorption.a The variables include GAP (eV), the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO); I (eV), ionization energy required to remove an electron from a molecule; A (eV), electron affinity indicating the energy change when an electron is added; χ (eV), electronegativity representing the ability to attract electrons; η (eV), chemical hardness measuring resistance to electron distribution changes; S (eV −1 ), chemical softness as the inverse of hardness; ω (eV), electrophilicity index indicating the tendency to accept electrons; and ΔNmax, the maximum amount of electronic charge transfer.
The reactivity parameters of the pesticides in Table 2 provide insights into their chemical behaviors and interactions with nanoclusters.The energy gap (GAP) indicates the stability and reactivity of the molecules; a smaller GAP, like that of parathion (Par) at 2.417 eV, suggests higher reactivity compared to a larger GAP, such as trichlorfon (Tri) at 4.783 eV, indicating lower reactivity.Ionization energy (I) reflects the energy required to remove an electron, with higher values like naled (Nal) at 6.957 eV indicating more stability.Electron affinity (A) shows the energy change when an electron is added; parathion (Par) has a high electron affinity of 3.757 eV, suggesting it can easily gain electrons.Electron affinity is the energy change when an electron is added to a neutral atom or molecule in the gas phase to form a negative ion.A positive EA value indicates an exothermic process where energy is released, meaning the negative ion is more stable than the neutral atom.Conversely, a negative EA value signifies an endothermic process where energy is required, suggesting the neutral atom does not favor electron addition and the negative ion is less stable.Most atoms have positive EAs, indicating they release energy when gaining an electron, leading to a more stable state.However, some atoms, particularly noble gases, have negative EAs due to their stable electron configurations, requiring energy input to add an electron [29].
Electronegativity (χ) reflects the ability to attract electrons, with parathion (Par) having a high value of 4.966 eV, indicating strong electron attraction.Chemical hardness (η) measures the resistance to electron distribution changes, where higher values like sulfotep (Sul) at 2.695 eV suggest greater resistance.Chemical softness (S) is the inverse of hardness; higher softness values like Si12C12_Tyr337 at 0.470 eV −1 indicate greater reactivity.The electrophilicity index (ω) shows the tendency to accept electrons; parathion (Par) has a high electrophilicity of 10.201 eV, suggesting strong electron acceptance capability.The ΔNmax represents the maximum electronic charge a system can accept, with parathion (Par) having the highest value at 4.109, indicating a high capacity to gain electrons.Comparing these parameters helps in understanding the relative stability and reactivity of the pesticides, where parathion (Par) stands out for its high reactivity and electrophilicity, while trichlorfon (Tri) and sulfotep (Sul) exhibit higher stability.
A study explored the adsorption of glycine amino acid and its zwitterionic form onto hexagonal boron nitride (h-BN) and silicon carbide (h-SiC) sheets using density functional theory (DFT) calculations.The results showed that glycine and its conformers tend to chemisorb onto the h-SiC surface more strongly compared to h-BN, indicating that h-SiC has higher affinity due to stronger interaction energies [30].Another study examined the adsorption of various organic molecules on boron nitride nanosheets, showing that functionalization with specific molecules can significantly enhance binding energies.This finding supports the hypothesis that modifying nanoclusters with amino acids can improve their interaction with target molecules, such as organophosphorus pesticides [31].Research into the electronic properties of boron nitride and its adsorption sites indicated that modifications on the nanosheet surfaces, such as with sulfur compounds, can enhance adsorption through non-covalent interactions.These modifications affect the electronic affinity and reactivity of the nanosheets, which is crucial for applications in pesticide binding and detection [30,31].

Linear Correlation
The linear correlation between the HSAB parameters and ΔE for nanoclusters functionalized with Tyr124 is low, indicating that nanoclusters functionalized with Tyr337 exhibit relatively more stable structures when adsorbing OP pesticides (Figure 4).Compared to Tyr124 functionalization, the adsorption energies of nanoclusters functionalized with Tyr337 are generally more predictable, particularly for B12N12 nanoclusters.This suggests that Tyr337 may provide a more stable structural environment for adsorbing different OP pesticides.The energy gap (GAP) represents the energy difference between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO).A smaller GAP indicates higher reactivity.Generally, a smaller GAP may lead to higher adsorption energy (more negative values) due to easier electron transfer.The results show that molecules with smaller GAP values (e.g., Par and Dia) have higher adsorption energies on B12N12_Tyr337.Ionization energy (I) indicates the minimum energy required to remove an electron from a molecule.Lower ionization energy means the molecule loses electrons more easily.The OP pesticides with lower ionization energies (e.g., Sul) have higher adsorption energies on B12N12_Tyr337, indicating a negative correlation between ionization energy and adsorption energy.Electron affinity (A) indicates the energy change when a molecule accepts an electron.Higher electron affinity means the molecule accepts electrons more easily.The OP pesticides with higher electron affinities (e.g., Par) have higher adsorption energies on B12N12_Tyr337, suggesting a positive correlation between electron affinity and adsorption energy.Electronegativity (χ) measures the ability of a molecule to attract electrons.Higher electronegativity means the molecule attracts electrons more easily.The OP pesticides with higher electronegativities (e.g., Par) have higher adsorption energies on B12N12_Tyr337, indicating a positive correlation between electronegativity and adsorption energy.Chemical hardness (η) measures the resistance of a molecule to electron transfer.Higher chemical hardness means electron transfer is more difficult.The OP pesticides with lower chemical hardness (e.g., Par) have higher adsorption energies on B12N12_Tyr337, indicating a negative correlation between chemical hardness and adsorption energy.Softness (S), the inverse of chemical hardness, indicates the ease of electron transfer.The OP pesticides with higher softness (e.g., Par) have higher adsorption energies on B12N12_Tyr337, suggesting a positive correlation between softness and adsorption energy.The electrophilicity index (ω) is a comprehensive measure of a molecule's ability to accept electrons.The OP pesticides with higher electrophilicity indices (e.g., Par) have higher adsorption energies on B12N12_Tyr337, indicating a positive correlation between the electrophilicity index and adsorption energy.The OP pesticides with higher the maximum amount of electronic charge transfer (ΔNmax) values (e.g., Par) have higher adsorption energies on B12N12_Tyr337, suggesting a positive correlation between ΔNmax and adsorption energy.
Consequently, the B12N12 and Si12C12 nanoclusters with Tyr337 exhibit a linear relationship between adsorption energies and several HSAB parameters (e.g., GAP, I, A, χ, η, S, ω, and ΔNmax) when adsorbing the OP pesticides.These parameters collectively influence the stability and energy changes during the adsorption process.Notably, for the B12N12 nanoclusters with Tyr337, the adsorption energies show higher stability, consistent with the trends in HSAB parameters.These results also suggest that B12N12 nanoclusters have greater potential and application value as adsorption materials.

The Changes in HSAB Parameters after Adsorption
The GAP values decrease after adsorption (Table 3), indicating significant changes in the electronic structure during adsorption, leading to a reduced energy gap.The GAP of B12N12_Tyr337 decreases from 3.747 eV to 1.726 eV (B12N12_Tyr337_Par), and the GAP of Si12C12_Tyr337 decreases from 2.128 eV to 1.741 eV (Si12C12_Tyr337_Par).
The ΔGAP (eV) values for the B12N12_Tyr337 and Si12C12_Tyr337 nanoclusters are presented in Figure 5.The results reveal significant differences in the adsorption characteristics and stability between the two types of nanoclusters.The ΔGAP values for the B12N12_Tyr337 nanoclusters range from 0.002 eV to 2.021 eV.The highest ΔGAP is observed for B12N12_Tyr337_Par with a value of 2.021 eV, indicating a substantial energy change upon modification.Conversely, the lowest ΔGAP is noted for B12N12_Tyr337_Isa at 0.002 eV, suggesting minimal energy change.Overall, the average ΔGAP for the B12N12_Tyr337 nanoclusters is relatively high, indicating a stable interaction when they are modified with Tyr337.For the Si12C12_Tyr337 nanoclusters, the ΔGAP values vary from −0.027 eV to 0.607 eV.The maximum ΔGAP is observed for Si12C12_Tyr337_Sul at 0.607 eV, signifying a moderate energy change.The minimum ΔGAP is recorded for Si12C12_Tyr337_Dia at −0.027 eV, indicating a slight decrease in energy.The average ΔGAP values for the Si12C12_Tyr337 nanoclusters are lower than those of B12N12_Tyr337, implying less stability.
The ΔGAP values for the B12N12_Tyr337 nanoclusters are generally higher and more consistent compared to those of the Si12C12_Tyr337 nanoclusters.This indicates that the B12N12_Tyr337 nanoclusters exhibit more significant energy changes, associated with higher stability and stronger interactions with adsorbed species.In contrast, the Si12C12_Tyr337 nanoclusters show a broader variation and lower ΔGAP values, suggesting less stable interactions.The higher stability of the B12N12_Tyr337 nanoclusters, as evidenced by their ΔGAP values, suggests they have greater potential as adsorption materials.Their consistent and substantial energy changes upon modification highlight their suitability for applications requiring stable and robust adsorption properties.The variability and lower stability of the Si12C12_Tyr337 nanoclusters indicate they may be less effective for such purposes.The analysis demonstrates that B12N12 nanoclusters, particularly when modified with Tyr337, show greater potential and application value as adsorption materials due to their higher stability and stronger interactions.This makes them more promising candidates for practical applications where stable adsorption is crucial compared to Si12C12_Tyr337 nanoclusters.
The ionization energy (I) changes slightly after adsorption.The ionization energy of B12N12_Tyr337 increases from 5.395 eV to 5.449 eV (B12N12_Tyr337_Par), and that of Si12C12_Tyr337 increases from 5.208 eV to 5.530 eV (Si12C12_Tyr337_Par).This suggests an enhancement in the stability of the electronic structure during adsorption.The electron affinity (A) also changes significantly after adsorption.The electron affinity of B12N12_Tyr337 increases from 1.648 eV to 3.723 eV (B12N12_Tyr337_Par), and that of Si12C12_Tyr337 increases from 3.080 eV to 3.789 eV (Si12C12_Tyr337_Par).This indicates that the system becomes more electron-accepting after adsorption.
The adsorption process also results in an increased electronegativity and electrophilicity index and decreased chemical hardness, indicating that the electronic structure of the nanoclusters becomes more flexible and reactive during adsorption.Moreover, B12N12 exhibits greater chemical hardness and stability after adsorption, while Si12C12 shows higher electrophilicity and electron transfer capacity.The electronegativity (χ) increases after adsorption, indicating a stronger attraction for electrons.The electronegativity of B12N12_Tyr337 increases from 3.522 eV to 4.586 eV (B12N12_Tyr337_Par), and that of Si12C12_Tyr increases from 4.144 eV to 4.660 eV (Si12C12_Tyr337_Par).The chemical hardness (η) generally decreases after adsorption, making the system more prone to electron transfer.The chemical hardness of B12N12_Tyr337 decreases from 1.874 eV to 0.863 eV (B12N12_Tyr337_Par), and that of Si12C12_Tyr337 decreases from 1.064 eV to 0.871 eV (Si12C12_Tyr337_Par).The softness (S) increases, indicating a greater ability to respond to external electrons.The softness of B12N12_Tyr337 increases from 0.267 eV⁻¹ to 0.579 eV⁻¹ (B12N12_Tyr337_Par), and that of Si12C12_Tyr337 increases from 0.470 eV⁻¹ to 0.574 eV⁻¹ (Si12C12_Tyr337_Par).The electrophilicity index (ω) significantly increases after adsorption, indicating a greater tendency to accept electrons.The electrophilicity index of B12N12_Tyr337 increases from 3.310 eV to 12.185 eV (B12N12_Tyr337_Par), and that of Si12C12_Tyr337 increases from 8.070 eV to 12.470 eV (Si12C12_Tyr337_Par).The maximum amount of electronic charge transfer (ΔNmax) also increases, indicating that the system can accept more electrons.The ΔNmax of B12N12_Tyr337 increases from 1.880 to 5.314 (B12N12_Tyr337_Par), and that of Si12C12_Tyr337 increases from 3.895 to 5.353 (Si12C12_Tyr337_Par).

Conclusions
This study demonstrates that B12N12_Tyr337 nanoclusters exhibit stronger interactions when adsorbing OP pesticide molecules compared to Si12C12_Tyr337.B12N12_Tyr337 possesses higher ionization energy and chemical hardness, indicating greater electronic stability.The adsorption energies of B12N12_Tyr337 are more negative and stable, suggesting that B12N12_Tyr337 provides a more stable structural environment for OP pesticide adsorption.In contrast, although Si12C12_Tyr337 has a higher electron acceptance capability, its adsorption energies vary more significantly, indicating lower consistency.Therefore, B12N12 nanoclusters modified with amino acid tyrosine have greater potential and value for applications as adsorption and biosensing materials.

Figure 2 .
Figure 2.After docking AChE with the OPs (Par, Cou, Dia, and Nal), the structures of Tyr337_OP were extracted and used as the starting structures for modification with B12N12 nanoclusters.These modified structures were then subjected to DFT geometry optimization.

Figure 3 .
Figure 3.After docking AChE with the OPs (Tri, Sul, Isa, and Chl), the structures of Tyr337_OP were extracted and used as the starting structures for modification with B12N12 nanoclusters.These modified structures were then subjected to DFT geometry optimization.

Table 2 .
The HSAB parameters before adsorption a.