RESUMEN

The use of noncovalent interactions to control the regioselectivity of transition metal-catalyzed C-H functionalization of arenes has received significant attention in recent years. Herein, we present a mechanistic study based on Density Functional Theory (DFT) of palladium(II)-catalyzed meta-selective C-H olefination employing a noncovalent directing template. We analyze the key steps of the mechanism and discuss the origins of reaction selectivity. The role of the directing template was elucidated, demonstrating its essential function in lowering reaction barriers and controlling selectivity. Our results reveal a competition in activation between ortho- and meta-C-H bonds. Contrary to the previous proposal in the literature, hydrogen bonds between the N-H bonds of the urea moiety and the carbonyl oxygen of the substrate predominantly favor ortho-selectivity over meta-selectivity. DFT results, alongside Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interaction Index analysis, suggest that secondary interactions between the R group linked to the urea moiety and the catalyst exert a more pronounced influence compared to the aforementioned hydrogen bonds, directing the selectivity towards the meta C-H bond.

RESUMEN

CONTEXT: For the first time, the use of monocyclic rings C18 and B9N9 as sensors for the sensing of carbazole-based anti-cancer drugs, such as tetrahydrocarbazole (THC), mukonal (MKN), murrayanine (MRY), and ellipticine (EPT), is described using DFT simulations and computational characterization. The geometries, electronic properties, stability studies, sensitivity, and adsorption capabilities of C18 and B9N9 counterparts towards the selected compounds confirm that the analytes interact through active cavities of the C18 and B9N9 rings of the complexes. METHODS: Based on the interaction energies, the sensitivity of surfaces towards EPT, MKN, MRY, and THC analytes is observed. The interaction energy of EPT@B9N9, MKN@B9N9, MRY@B9N9, and THC@B9N9 complexes are observed - 20.40, - 19.49, - 20.07, and - 18.27 kcal/mol respectively which is more exothermic than EPT@C18, MKN@C18, MRY@C18, and THC@C18 complexes are - 16.37, - 13.97, - 13.96, and - 11.39 kcal/mol respectively. According to findings from the quantum theory of atoms in molecules (QTAIM) and the reduced density gradient (RDG), dispersion forces play a significant role in maintaining the stability of these complexes. The electronic properties including FMOs, density of states (DOS), natural bond orbitals (NBO), charge transfer, and absorption studies are carried out. In comparison of B9N9 and C18, the analyte recovery time for C18 is much shorter (9.91 × 10-11 for THC@C18) than that for B9N9 shorter recovery time value of 3.75 × 10-9 for EPT@B9N9. These results suggest that our reported sensors B9N9 and C18 make it faster to detect adsorbed molecules at room temperature. The sensor response is more prominent in B9N9 due to its fine energy gap and high adsorption energy. Consequently, it is possible to think of these monocyclic systems as a potential material for sensor applications.

Asunto(s)

RESUMEN

The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme's dispersion with Becke-Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50-800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.

RESUMEN

Pyridoxal-S-methyl-isothiosemicarbazone (PLITSC) is a member of an important group of ligands characterized by different complexation modes to various transition metals. In this contribution, a new complex containing two differently protonated PLITSC ligands ([Fe(PLITSC-H)(PLITSC)]SO4)∙2.5H2O was obtained. The crystal structure was solved by the X-ray analysis and used further for the optimization at B3LYP/6-311++G(d,p)(H,C,N,O,S)/def2-TZVP(Fe) level of theory. Changes in the interaction strength and bond distance due to protonation were observed upon examination by the Quantum Theory of Atoms in Molecules. The protein binding affinity of [Fe(PLITSC-H)(PLITSC)]SO4 towards transport proteins (Bovine Serum Albumin (BSA) and Human Serum Albumin (HSA)) was investigated by the spectrofluorimetric titration and molecular docking. The interactions with the active pocket containing fluorescent amino acids were examined in detail, which explained the fluorescence quenching. The interactions between complex and DNA were followed by the ethidium-bromide displacement titration and molecular docking. The binding along the minor groove was the dominant process involving complex in the proximity of DNA.

Asunto(s)

RESUMEN

Chalcogen bond (ChB) catalysis, as a new type in the field of non-covalent bond catalysis, has become a hot research topic in the field of organocatalysis in recent years. In the present work, we investigated the catalytic performance of a series of hypervalent ChB catalysis based on the intramolecular Aza-Michael reaction of aminochalcone. The reaction includes the carbon-nitrogen bond coupling step (key step) and the proton transfer step. The catalytic performance of mono-dentate pentafluorophenyl chalcogen bond donor ChB1 was comparable to that of bis-dentate chalcogen bond donor ChB4, and stronger than that of mono-dentate chalcogen bond donors ChB2 and ChB3. The formation of the chalcogen bond between the catalyst and the carbonyl oxygen atom of the reactant, causing the charge rearrangement of the reactant and C(1) charge of the -C-Ph group to become more positive, thereby the ChB catalysis promoted the nucleophile reaction. The electron density of the chalcogen bond of the pre-complex, the most positive electrostatic potentials of the catalyst, and the NPA charge of the key atom are proportional to the Gibbs energy barrier of the C-N bond coupling process, which provides an idea to predict the catalytic activity of the ChB catalysis.

RESUMEN

In this proof-of-concept paper, we show how exchange-correlation effects can be simply recovered for interatomic energies within the interacting quantum atoms decomposition when local, gradient generalized, or meta-gradient generalized approximations are used in density functional theory (DFT) calculations. We also demonstrate how inhomogeneity and non-local effects can be introduced even from a pure local scheme, without resorting to any orbital information. Finally, we provide numerical evidence on a database of selected energetic molecules that this decomposition scheme can be efficiently used to build accurate models for the prediction of molecular energies from an initial "cheap" DFT calculation.

RESUMEN

A microwave-assisted synthesis of 7-amino-1,2,4-triazolo[1,5-a][1,3,5]triazine-2-propanamides was developed using a three-component, catalyst-free reaction of cyanamide and trimethyl orthoformate with 3-(5-amino-1H-1,2,4-triazol-3-yl)propanamides (3). The reaction tolerated structurally diverse substrates and proceeded chemo- and regio-selectively, affording the target compounds in high purity in 5-10 minutes. The convenient chromatography-free isolation and purification of the products add practicality to this method. The structural features of the prepared compounds were investigated using dynamic NMR spectroscopy, X-ray crystallography and computational chemistry calculations. X-ray crystallography performed on a representative compound, 3-(7-amino-1,2,4-triazolo[1,5-a][1,3,5]triazin-2-yl)-N-(4-benzyl)propanamide (4 l), showed the overall molecular conformation to adopt the shape of the letter C. Notable localisation of π-electron density is found within the 1,2,4-triazolo[1,5-a][1,3,5]triazine system; a relatively short C-NH2 bond is consistent with restricted rotation about this bond. This study also presents a detailed analysis of the molecular interactions in 4 l using DFT and QTAIM methods with a focus on the hydrogen-bonding and π-stacking interactions that influence the molecular packing of 4 l. The findings reveal the significant roles of N-H⋅O, N-H⋅N and C-H⋅N interactions, along with electrostatically enhanced π⋅π contacts. A broad screening for insecticidal, fungicidal and herbicidal properties identified several compounds with potent herbicidal activity against Matricaria inodora.

RESUMEN

Bimetallic CpMM'Nacnac molecules with group 2 and 12 metals (M=Be, Mg, Ca, Zn, Cd, Hg) that contain novel metal-metal bonding have been investigated in a theoretical study of their molecular and electronic structure, thermodynamic stability, and metal-metal bonding. In all cases the metal-metal bonds are characterized as electron-sharing covalent single bonds from natural bond orbital (NBO) and energy-decomposition analysis with natural orbitals of chemical valence (EDA-NOCV) analysis. The sum of [MM'] charges is relatively constant, with all complexes exhibiting a [MM']2+ core. Quantum theory of atoms in molecules (QTAIM) analysis indicates the presence of non-nuclear attractors (NNA) in the metal-metal bonds of the BeBe, MgMg, and CaCa complexes. There is substantial electron density (0.75-1.33 e) associated with the NNAs, which indicates that these metal-metal bonds, while classified as covalent electron-sharing bonds, retain significant metallic character that can be associated with reducing reactivity of the complex. The predicted stability of these complexes, combined with their novel covalent metal-metal bonding and potential as reducing agents, make them appealing targets for the synthesis of new metal-metal bonds.

RESUMEN

This study delves into the nature of individual hydrogen bonds and the relationship between metal cations and hydrogen bonding in the Watson-Crick guanine-cytosine (GC) base pair and its alkali and alkaline earth cation-containing complexes (Mn+-GC). The findings reveal how metal cations affect the nature and strength of individual hydrogen bonds. The study employs interacting quantum atoms (IQA) analysis to comprehensively understand three individual hydrogen bonds within the GC base pair and its cationic derivatives. These analyses unveil the nature and strength of hydrogen bonds and serve as a valuable reference for exploring the impact of cations (and other factors) on each hydrogen bond. All the H ⋯ D interactions (H is hydrogen and D is oxygen or nitrogen) in the GC base pair are primarily electrostatic in nature, with the charge transfer component playing a substantial role. Introducing a metal cation perturbs all H ⋯ D interatomic interactions in the system, weakening the nearest hydrogen bond to the cation (indicated by a) and reinforcing the other (b and c) interactions. Notably, the interaction a, the strongest H ⋯ D interaction in the GC base pair, becomes the weakest in the Mn+-GC complexes. A broader perspective on the stability of GC and Mn+-GC complexes is provided through interacting quantum fragments (IQF) analysis. This approach considers all pairwise interactions between fragments and intra-fragment components, offering a complete view of the factors that stabilize and destabilize GC and Mn+-GC complexes. The IQF analysis underscores the importance of electron sharing, with the dominant contribution arising from the inter-fragment exchange-correlation term, in shaping and sustaining GC and Mn+-GC complexes. From this point of view, alkaline and alkaline earth cations have distinct effects, with alkaline cations generally weakening inter-fragment interactions and alkaline earth cations strengthening them. In addition, IQA and IQF calculations demonstrate that the hydration of cations led to small changes in the hydrogen bonding network. Finally, the IQA interatomic energies associated with the hydrogen bonds and also inter-fragment interaction energies provide robust indicators for characterizing hydrogen bonds and complex stability, showing a strong correlation with total interaction energies.

Asunto(s)

RESUMEN

Observational data show complex organic molecules in the interstellar medium (ISM). Hydrogenation of small unsaturated carbon double bond could be one way for molecular complexification. It is important to understand how such reactivity occurs in the very cold and low-pressure ISM. Yet, there is water ice in the ISM, either as grain or as mantle around grains. Therefore, the addition of atomic hydrogen on double-bonded carbon in a series of seven molecules have been studied and it was found that water catalyzes this reaction. The origin of the catalysis is a weak charge transfer between the π MO of the unsaturated molecule and H atom, allowing a stabilizing interaction with H2O. This mechanism is rationalized using the non-covalent interaction and the quantum theory of atoms in molecules approaches.

RESUMEN

CONTEXT: This research aims to offer a deeper understanding of the bonding interactions between M-Se and M-CO and how these interactions change across the group 6 transition metal series: [Se2M3(CO)10]2- (M = Cr, Mo, W). It also seeks to explore the impact of carbonyl groups on M-M interactions within the clusters. Seven criteria, which are based on QTAIM properties, have been considered and compared with the corresponding criteria in other transition metal clusters. The results confirm that no such bond critical points or bond baths occur between transition metals, which instead have 5c-7e bonding interactions delocalized over their five-membered M3(µ-Se)2 ring, as evidenced by the non-negligible nonbonding delocalization indices. The topological properties of three bond clusters, Cr-Se, Mo-Se, and W-Se, resemble those of "intermediate closed shell characters," which combine covalent and electrostatic properties. Source function calculations indicated that the bonded Se atom contributed the most to each Cr-Se and Mo-Se bcp. The OCO atoms and nonbonded Se atoms also contributed to some extent. However, metal atoms act as sinks rather than as sources of electron density. In contrast, the majority of the metal atoms, both bonded and nonbonded, contribute to Cr-W bcps. Analysis of the delocalization indices δ(M…O) in the three clusters indicates that CO significantly contributes to Cr π-back donation in cluster 1. In contrast, no π-back donation occurs from CO to Mo or W in clusters 2 or 3, respectively. METHODS: The B3P86 hybrid functional was used for computations in the Gaussian 09 software. The LanL2DZ basis set was employed for Cr, Mo, and W, while the 6-31G (d, p) basis set was used for C, O, and Se atoms. We performed QTAIM analysis using the AIM2000 and Multiwfn packages, incorporating B3P86/WTBS for Cr, Mo, and W atoms. The 6-311++G(3df,3pd) basis set was used for C, O, and Se atoms. Additionally, we utilized the ELF and SF.

RESUMEN

N-(Benzothiazole-2-yl)pyrrolamide DNA gyrase inhibitors with benzyl or phenethyl substituents attached to position 3 of the benzothiazole ring or to the carboxamide nitrogen atom were prepared and studied for their inhibition of Escherichia coli DNA gyrase by supercoiling assay. Compared to inhibitors bearing the substituents at position 4 of the benzothiazole ring, the inhibition was attenuated by moving the substituent to position 3 and further to the carboxamide nitrogen atom. A co-crystal structure of (Z)-3-benzyl-2-((4,5-dibromo-1H-pyrrole-2-carbonyl)imino)-2,3-dihydrobenzo[d]-thiazole-6-carboxylic acid (I) in complex with E. coli GyrB24 (ATPase subdomain) was solved, revealing the binding mode of this type of inhibitor to the ATP-binding pocket of the E. coli GyrB subunit. The key binding interactions were identified and their contribution to binding was rationalised by quantum theory of atoms in molecules (QTAIM) analysis. Our study shows that the benzyl or phenethyl substituents bound to the benzothiazole core interact with the lipophilic floor of the active site, which consists mainly of residues Gly101, Gly102, Lys103 and Ser108. Compounds with substituents at position 3 of the benzothiazole core were up to two orders of magnitude more effective than compounds with substituents at the carboxamide nitrogen. In addition, the 6-oxalylamino compounds were more potent inhibitors of E. coli DNA gyrase than the corresponding 6-acetamido analogues.

Asunto(s)

RESUMEN

This article explores the possible presence of a pentacle valence bond structure in C 5 ${_5 }$ cyclic molecules. At this end, we have used quantum chemistry tools to elucidate the possible arrangement and the nature of chemical bonds within linear, cyclic, and three-dimensional structures only formed by five carbon atoms. While the linear structure is clearly the most stable one, local minima were obtained for both bi- and three-dimensional structures. Using the localization-delocalization matrices approach, we characterize both the minimum linear structure and the cyclic ones. Interestingly, the linear structure is a combination of ionic and covalent bonds, albeit the four distances are almost identical, when using Density Functional Theory. For cyclic C 5 ${_5 }$ , the pentacle bonding arrangement emerges as a significant Lewis structure, indicative of an unusual formal configuration characterized by five intersecting C-C bonds. Our calculations show that this pentacle arrangement in cyclic C 5 ${_5 }$ scheme is also present in the more known cyclo-pentadienyl molecule.

RESUMEN

The synthesis of two new hexahydroisoquinoline-4-carbonitrile derivatives (3a and 3b) is reported along with spectroscopic data and their crystal structures. In compound 3a, the intramolecular O-H···O hydrogen bond constraints the acetyl and hydroxyl groups to be syn. In the crystal, inversion dimers are generated by C-H···O hydrogen bonds and are connected into layers parallel to (10-1) by additional C-H···O hydrogen bonds. The layers are stacked with Cl···S contacts 0.17 Å less than the sum of the respective van der Waals radii. The conformation of the compound 3b is partially determined by the intramolecular O-H···O hydrogen bond. A puckering analysis of the tetrahydroisoquinoline unit was performed. In the crystal, O-H···O and C-H···O hydrogen bonds together with C-H···π(ring) interactions form layers parallel to (01-1) which pack with normal van der Waals interactions. To understand the binding efficiency and stability of the title molecules, molecular docking, and 100 ns dynamic simulation analyses were performed with CDK5A1. To rationalize their structure-activity relationship(s), a DFT study at the B3LYP/6-311++G** theoretical level was also done. The 3D Hirshfled surfaces were also taken to investigate the crystal packings of both compounds. In addition, their ADMET properties were explored.Communicated by Ramaswamy H. Sarma.

Asunto(s)

RESUMEN

Almost all phase-change memory materials (PCM) contain chalcogen atoms, and their chemical bonds have been denoted both as 'electron-deficient' [sometimes referred to as 'metavalent'] and 'electron-rich' ['hypervalent', multicentre]. The latter involve lone-pair electrons. We have performed calculations that can discriminate unambiguously between these two classes of bond and have shown that PCM have electron-rich, 3c-4e ('hypervalent') bonds. Plots of charge transferred between (ET) and shared with (ES) neighbouring atoms cannot on their own distinguish between 'metavalent' and 'hypervalent' bonds, both of which involve single-electron bonds. PCM do not exhibit 'metavalent' bonding and are not electron-deficient; the bonding is electron-rich of the 'hypervalent' or multicentre type.

RESUMEN

The current study comprehensively investigates the adsorption behavior of chromium (Cr(III)) in wastewater using Algerian kaolinite clay. The structural and textural properties of the kaolinite clay are extensively characterized through a range of analytical methods, including XRD, FTIR, SEM-EDS, XPS, laser granulometry, N2 adsorption isotherm, and TGA-DTA. The point of zero charge and zeta potential are also assessed. Chromium adsorption reached equilibrium within five minutes, achieving a maximum removal rate of 99% at pH 5. Adsorption equilibrium is modeled using the Langmuir, Freundlich, Temkin, Elovich, and Dubinin-Radushkevitch equations, with the Langmuir isotherm accurately describing the adsorption process and yielding a maximum adsorption capacity of 8.422 mg/g for Cr(III). Thermodynamic parameters suggest the spontaneous and endothermic nature of Cr(III) sorption, with an activation energy of 26.665 kJ/mol, indicating the importance of diffusion in the sorption process. Furthermore, advanced DFT computations, including COSMO-RS, molecular orbitals, IGM, RDG, and QTAIM analyses, are conducted to elucidate the nature of adsorption, revealing strong binding interactions between Cr(III) ions and the kaolinite surface. The integration of theoretical and experimental data not only enhances the understanding of Cr(III) removal using kaolinite but also demonstrates the effectiveness of this clay adsorbent for wastewater treatment. Furthermore, this study highlights the synergistic application of empirical research and computational modeling in elucidating complex adsorption processes.

RESUMEN

This study reveals a new non-covalent interaction called a π-hole halogen bond, which is directional and potentially non-linear compared to its sister analog (σ-hole halogen bond). A π-hole is shown here to be observed on the surface of halogen in halogenated molecules, which can be tempered to display the aptness to form a π-hole halogen bond with a series of electron density-rich sites (Lewis bases) hosted individually by 32 other partner molecules. The [MP2/aug-cc-pVTZ] level characteristics of the π-hole halogen bonds in 33 binary complexes obtained from the charge density approaches (quantum theory of intramolecular atoms, molecular electrostatic surface potential, independent gradient model (IGM-δginter)), intermolecular geometries and energies, and second-order hyperconjugative charge transfer analyses are discussed, which are similar to other non-covalent interactions. That a π-hole can be observed on halogen in halogenated molecules is substantiated by experimentally reported crystals documented in the Cambridge Crystal Structure Database. The importance of the π-hole halogen bond in the design and growth of chemical systems in synthetic chemistry, crystallography, and crystal engineering is yet to be fully explicated.

Asunto(s)

RESUMEN

Even though molecules are fundamentally quantum entities, the concept of a molecule retains certain classical attributes concerning its constituents. This includes the empirical separability of a molecule into its three-dimensional, rigid structure in Euclidean space, a framework often obtained through experimental methods like X-Ray crystallography. In this work, we delve into the mathematical implications of partitioning a molecule into its constituent parts using the widely recognized Atoms-In-Molecules (AIM) schemes, aiming to establish their validity within the framework of Information Theory concepts. We have uncovered information-theoretical justifications for employing some of the most prevalent AIM schemes in the field of Chemistry, including Hirshfeld (stockholder partitioning), Bader's (topological dissection), and the quantum approach (Hilbert's space definition). In the first approach we have applied the generalized principle of minimum relative entropy derived from the Sharma-Mittal two-parameter functional, avoiding the need for an arbitrary selection of reference promolecular atoms. Within the ambit of topological-information partitioning, we have demonstrated that the Fisher information of Bader's atoms conform to a comprehensive theory based on the Principle of Extreme Physical Information avoiding the need of employing the Schwinger's principle, which has been proven to be problematic. For the quantum approach we have presented information-theoretic justifications for conducting Löwdin symmetric transformations on the density matrix to form atomic Hilbert spaces generating orthonormal atomic orbitals with maximum occupancy for a given wavefunction.

RESUMEN

A DFT-based computational study is carried out to delve into the interplay between hyperconjugation and rehybridization effects underlying the formation of blue- or red-shifting H-bonds (HBs) in 1 : 1 complexes of cyclic ethers (HB acceptor) of varying ring-size with haloforms, CHF3 and CHCl3 (HB donor). The calculations reveal that with decreasing angular strain (increasing ring-size) of the cyclic ethers, the extent of blue-shift increases for 1 : 1 complexes with CHF3, while a reverse sequence is observed with CHCl3, eventually leading to a red-shifting HB in the oxepane : CHCl3 complex. It is noted that the trend in the shift of C-H stretching fundamental is not mirrored by the C-H bond length or interaction energies for both the systems studied, that is, the low sensitivity of the changes on the strain on the O-atom of HB acceptor (cyclic ethers) is to be emphasized.

RESUMEN

CONTEXT: This study delves into the chemical nuances of thiophenols and their derivatives through a comprehensive computational analysis, moving beyond traditional energetic perspectives such as bond dissociation enthalpy and S-H dissociation dynamics. By employing the overlap model along with its topological descriptors (OP/TOP), quantum theory of atoms in molecules (QTAIM), and local vibrational mode (LVM) theories, the research provides a deeper understanding of the S-H and C-S bonding scenarios in substituted thiophenols. The investigation follows the electron-donating capacity of S-H substituent variation with the nature and positioning of other ring substituents. Energy profile analyses indicate distinct stability differences in the cis and trans conformations of meta- and para-PhSH systems, influenced by the electron-donating strength of these substituents. The study also uncovers significant variations in S-H bond distances and descriptor values, particularly in para-substituted PhSH, reflecting the influence of electron-donating or withdrawing substituents. In contrast, alterations at the meta-position show minimal effects on C-S bond descriptors, while para-substitutions markedly influence C-S bond characteristics, demonstrating a clear correlation with the electron-donating or withdrawing capabilities of the substituents. This research sheds light on the intricate bond dynamics in aromatic systems with diverse substituents, highlighting the complex interaction between electronic effects and molecular conformation. METHODS: The study employs the ω B97X-D/Def2TZVP level of theory for molecular geometries, ensuring accurate characterization of structures as true minima via analytical harmonic frequency determination. The electronic properties of S-H and C-S bonds in variously substituted thiophenols were analyzed using OP/TOP, QTAIM, and LVM methodologies. Computational processes, including conformational scans, geometry optimizations, and vibrational frequency calculations, were conducted using Gaussian 09, with ultra-fine integration grids and tight convergence criteria for the SCF procedure. Bond descriptors were computed utilizing ChemBOS, Multiwfn, and LModeA software, providing a robust and detailed examination of bond properties.