Interactions of B-substituted Borthiin molecules with diborane: A quantum chemistry investigation

Abstract

In this quantum chemistry study, the B₂H₆…B₃S₃H₃ interaction was illustrated at M062X/6-311G(d,p) level of theory. Substitution of hydrogens of borthiin with F, Cl, CN and Me substituents was clarified. Energy decomposition analysis (EDA) was considered to exploring of the nature of the B₂H₆…B₃S₃X₃ interactions. The interacting sites of the two fragments were reflected with molecular electrostatic potential (MEP) surfaces. Substituent effect on the polarity, electronic spatial extent and HOMO-LUMO gap of the B₂H₆…B₃S₃X₃ complexes were stated. Quantum theory of atoms in molecules (QTAIM) analysis was considered to determination of bond critical points (BCP) between B₂H₆ and B₃S₃X₃ molecules. The electron transfer between two molecules was studied with charge decomposition analysis (CDA). Besides, interactions between B₂H₆ and B₃S₃X₃ molecules were exposed with interacting quantum atoms (IQA) approach.

Keywords

Borthiin diborane substitution effect QTAIM IQA

Introduction

Inorganic aromatic molecules have been interesting molecules in the inorganic chemistry. Therefore, attentions of inorganic chemists exited to preparation and characterization of inorganic aromatic rings. Borthiin (B₃S₃H₃) molecule is one of this family compounds. Preparations and characterizations of this molecule have been excited much attentions. In a study synthesis of B₃S₃(NMe₂)₃ molecule has been published.¹ Quantum chemistry tools have been illustrated structure and properties of borthiin molecule.^2–6

Energetics, structural, and spectral parameters have been considered to illustration of hydrogen bonding (H-bonding) interaction.^7–16 A “Hydrogen Bridge” has been used to discrete all of the non-bonded interactions in which hydrogen atom bridge between any two atoms.¹⁷ Numerous boron compounds contain bridge hydrogen atoms. For model, boron hydrides expose unique bonding characters and have involved many researches due to their exclusive chemical properties.^18–22 Diborane (B₂H₆) molecule has been selected as a typical compound for the consideration of an interaction between n-B₁₈H₂₂ and benzene in the molecular crystal.²³ In this paper, the H_bridge atom of n-B₁₈H₂₂ straight interacts with the center of the C₆H₆ cycle.

The diborane … π-systems interaction has been motivated various investigations.^24,25 For example, substitution effects in the interaction of diborane and borazine have been illustrated.²⁶ In other study, a considered hypothetical reaction has been explored electronic structure of the H_bridge bond in diborane.²⁷ A computational research on substituent effect of terminal hydrogens of B₂H₆ with F atoms on the Borazanaphthalene…B₂H₆ interaction has been published.²⁸

In this quantum chemistry research, we attracted to understanding the interaction of B₂H₆ and borthiin. The substitution effect of hydrogens of borthiin on the structural, and electronic properties of theses systems was explored.

Computational methods

Optimization of the studied molecules was provided with hybrid functional of Truhlar and Zhao (M06-2X)functional,²⁹ 6-311G(d,p) basis set^30–33 and the software package of Gaussian 09.³⁴ Computations of harmonic vibrational frequencies were indicated no imaginary frequency for optimized molecules.

Energy decomposition analysis (EDA) was worked to understanding of interaction between the thiophene and nanoclusters with Multiwfn 3.8 software package.^35,36 The interaction energy (ΔE_int) between the B₂H₆ and B₃S₃X₃ molecules was valued as:

Δ E_{int} = Δ E_{polar} + Δ E_{els} + Δ E_{Ex}

where ΔE_polar, ΔE_els and ΔE_Ex are the electron density polarization (the induction term), electrostatic interaction and the exchange repulsion terms, respectively. Steric term (ΔE_steric) is defined as combinations of electrostatic interaction and exchange repulsion terms.

Molecular electrostatic potential (MEP) visualized by the GaussView-5.0 program. Atomic charges Atomic dipole corrected Hirshfeld atomic charge (ADCH) analysis,³⁷ and charge decomposition analysis (CDA)³⁸ were accomplished using Multiwfn 3.8 software package.^35,39

Results of quantum theory of atoms in molecules (QTAIM) analysis and interacting quantum atoms (IQA) were provided with the AIMAll software package.⁴⁰ The “Promega” basin integration method using fine interatomic surface mesh was considered.

The IQA partitioning pattern employs the QTAIM description of atomic basins to estimate intra- and interatomic energy terms. Sum of these values is expressed as the total energy of the system:

E_{I Q A}^{t o t a l} = \sum_{A = 1}^{N} E_{I Q A} (A)

where N and A are the total number of atoms in the system and an atom, respectively.

Each atomic energy can be extended as a sum of intra- and interatomic contributions:

E_{I Q A} (A) = E_{i n t r a} (A) + \frac{1}{2} \sum_{B \neq A}^{N} V_{i n t e r} (A, B)

with V_inter (A,B) is reflected as the potential energy between atoms A and B.

The intra-atomic term, E_intra (A), is the sum of the (intra-atomic) kinetic energy, T(A), and electron-electron, V_ee (A), and electron-nucleus V_en (A) potential energies.

The interatomic contributions are stated by:

V_{i n t e r} (A, B) = V_{n n} (A, B) + V_{e n} (A, B) + V_{n e} (A, B) + V_{e e} (A, B)

where, e and n depict electrons and nuclei, respectively. The electron-electron potential energy comprises two contributions: the classical Coulombic energy and the exchange-correlation components.

V_{e e} (A, B) = V_{c o u l} (A, B) + V_{X C} (A, B)

All Coulombic terms can be gathered together into the classical potential term V_cl(A,B).

Then, the interatomic term is stated as:

V_{i n t e r} (A, B) = V_{c l} (A, B) + V_{X C} (A, B)

Results and discussion

Figure 1 presents the optimized geometries of the B₂H₆…B₃S₃X₃ (X = H, F, Cl, CN, Me) systems at the M062X/6-311G(d,p) level of theory.

Figure 1.

The optimized geometries and direction of dipole moment vector of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at the M062X/6-311G(d,p) level of theory.

Energy decomposition analysis (EDA)

Energy decomposition analysis (EDA) is useful tool to clarification the B₂H₆ …B₃S₃X₃ (X = H, F, Cl, CN, Me) interactions. The calculated interaction energy data (ΔE_int) are documented in Table 1. It can be inferred stronger interaction arises with replacing of hydrogen atoms with substituents. The interaction strength decreases as Me > F > Cl > CN > H. Figure 2 indicates diagram of −ΔE_int variations of the B₃S₃X₃…B₂H₆ systems.

Figure 2.

Diagram of −ΔE_int variations of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at the M062X/6-311G(d,p) level of theory.

Table 1.

EDA results of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at M062X/6-311G(d,p) level of theory (in kcal/mol).

Molecule	ΔE_int	ΔE_polar	ΔE_steric
B₃S₃H₃…B₂H₆	−3.93	−4.38	0.45
B₃S₃F₃…B₂H₆	−4.61	−6.42	1.81
B₃S₃Cl₃…B₂H₆	−4.45	−5.63	1.17
B₃S₃(CN)₃…B₂H₆	−3.97	−4.92	0.95
B₃S₃Me₃…B₂H₆	−4.72	−4.34	−0.38

The major reason for this trend is attributed to in the nature of the interaction. B₂H₆ and B₃S₃X₃ molecules are classic Lewis acid and Lewis base, respectively. The strength of the Lewis acid-base interaction is straight related to electron donating ability of the base. Me group includes electron-donating ability. It makes the B₃S₃X₃ ring more electron-rich and highly attractive to the electron-deficient B₂H₆.

On the other hand, halogens are highly electronegative and would seem to be electron-withdrawing. But, they also include lone pairs of electrons that can be donated into the ring through a resonance effect. For fluorine and chlorine, this lone pair donation (resonance) is noteworthy and partially frustrates their strong inductive pull. This results in a net electron density that is moderately high, allowing for a decent interaction with B₂H₆.

Finally, the CN group is a powerful electron-withdrawing group. It pull out electron density through both a strong inductive effect (due to the electronegative nitrogen) and a strong resonance effect (via its π-bonds). This dual withdrawal leaves the B₃S₃X₃ ring electron-poor, creating it much less possible to interact strongly with B₂H₆.

The negative ΔE_polar values for all B₂H₆…B₃S₃X₃ systems expose that these data stabilize complex systems. The most stability is found in the presence of fluorine substituent.

It can be provided positive ΔE_steric values for these systems, except X = Me. Accordingly, these values destabilize B₂H₆ …B₃S₃X₃ systems. The most instability is seen in the presence of fluorine substituent. Hence, the main driving force of the B₂H₆ …B₃S₃X₃ interactions is polarization or orbital-interaction stabilization.

Polarity

Computed dipole moment values of the B₃S₃H_3−nX _n …B₂H₆ (X = H, F, Cl, CN, Me; n = 0–3) systems are given in Table 2. These values show that replacing of hydrogen atoms of borthiin with F, Cl, Me, and CN substitutions increase polarity. These values reveal increase as H < F < Me < Cl < CN. Orientations of dipole moment vector of the molecules are presented in Figure 1. It can be found, most polarity is seen in the presence of CN group. This group is a long and rigid spacer that maximizes the distance between separated charges, yielding the largest dipole despite minimal electron donation. Cl substituent is large and polarizable, letting a smaller charge to be separated over a longer distance, increasing dipole moment. Me group donates the most density, but the charge is spread over a diffuse volume, limiting the dipole. Finally, F substituent gives small density, so charge separation is modest.

Table 2.

Dipole moment (μ, debye), electronic spatial extent (ESE, au) and atomic dipole corrected Hirshfeld atomic charge (ADCH) of the bridged H and S atom in the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at M062X/6-311G(d,p) level of theory.

Molecule	μ	ESE	Q(H)	Q(S1)	Q(S2)
B₃S₃H₃…B₂H₆	0.17	1647.61	0.0271	−0.2505	−0.2517
B₃S₃F₃…B₂H₆	0.27	2501.87	0.0207	0.0139	0.0241
B₃S₃Cl₃…B₂H₆	0.31	3876.13	0.0175	−0.3583	−0.3726
B₃S₃(CN)₃…B₂H₆	0.56	3808.19	0.0149	0.1820	0.1833
B₃S₃Me₃…B₂H₆	0.20	2797.30	0.0309	−0.3923	−0.0791

Atomic dipole corrected Hirshfeld atomic charge (ADCH)

The calculated atomic dipole corrected Hirshfeld atomic charge (ADCH) on the bridged H and S atoms of B₃S₃H_3−nX _n …B₂H₆ systems are gathered in Table 2. Computed ADCH values of bridged hydrogen of B₂H₆ equals to 0.1167. on the other hand, calculated ADCH values of S atom of B₃S₃X₃ molecules are −0.0446 (H), −0.0695 (F) −0.0338 (Cl), 0.0548 (CN) and −0.0960 (Me). The smaller charges of bridged H atom are found B₃S₃H_3−nX _n …B₂H₆ systems than in isolated B₂H₆ molecule (Figure 3). These charges decrease as X = Me > H > F > Cl > CN. It can be gotten from the outcomes that the ADCH values on the H_bridge atom that is involved in the interaction with the π-cloud reduce with respect to the isolated state in the presence of X = F, Cl and CN. In contrast, the ADCH value on the H_bridge atom enhances with respect to the isolated state in the presence of X = Me.

Figure 3.

The calculated Hirshfeld atomic charges on the bridged H of B₃S₃H_3−nX _n …B₂H₆ in complexed and isolated B₂H₆ molecule in the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at the M062X/6-311G(d,p) level of theory.

Electronic spatial extent (ESE)

ESE is defined as the expectation value of the multiplied of electron density by the distance from a molecule mass center. More it is physical property of the electron density volume round the molecule. As the extent value enhances as electron cloud becomes more defused. Calculated ESE data of the B₃S₃H_3−nX _n …B₂H₆ systems are tabulated in Table 2.

Computed ESE values reveal replacing of hydrogen atoms of borthiin with F, Cl, Me, and CN substitutions increase ESE values. The ESE values increase as H < F < Me < CN < Cl. It must be noticed the ESE is dominated by two geometric factors: The intrinsic size of the substituent atom (Van der Waals radius) and the bond length between the ring and the substituent. The trend is a competition between how far the nuclei are spread out (geometry) and how tightly the electrons are held (electronegativity). It can be seen largest ESE in the presence of Cl substituent. It includes 3rd shell electrons (3s, 3p). These electrons are much farther from the nucleus than the 2nd shell electrons of F, C, or N. Besides, Cl substituent contains has lone pairs that are very “soft” and polarizable. Even though the B-Cl bond is strong, the electron cloud surrounding chlorine is obviously huge and diffuse. The slight negative charge that develops on Cl in the complex resides in these large, fluffy orbitals, causing in the largest ESE value. In linear CN group, N atom is small, the triple bond creates a rigid rod. This extends the substituent far away from the central borthiin ring. The π-electrons in the triple bond are also quite mobile. However, CN withdraws electron density, the distance of those electrons from the center of the complex is large, increasing ESE significantly. Me group has more significant volume than F substituent, however Me donates electrons (making the cloud denser), the physical presence of the carbon and three hydrogens pushes the electron cloud outward further than Fluorine's tight cloud. Fluorine electronegative substituent pulls electron density very close to its nucleus. Thus, while the nuclei are slightly further apart, the electrons are held tightly, limiting the ESE.

Frontier orbitals

Plots of frontier orbitals of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems are presented in Figure 4. It can be deduced most contribution is belonged to B₃S₃X₃ fragment.

Figure 4.

Plots of frontier orbitals of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at the M062X/6-311G(d,p) level of theory.

Evaluated energy values of frontier orbitals of these systems are gathered in Table 3. These values show that the stability of HOMO decreases as X = CN > F > Cl > H > Me. On the other hand, the stability of LUMO decreases as X = CN > Cl > H > F > Me.

Table 3.

Frontier orbital energy and global indices of reactivity for the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at M062X/6-311G(d,p) level of theory (in eV).

Molecule	E(HOMO)	E(LUMO)	Gap	η	μ	ω
B₂H₆	−10.56	0.12	10.68	5.34	−5.22	2.55
B₃S₃H₃	−9.63	−1.32	8.32	4.16	−5.48	3.61
B₃S₃F₃	−9.89	−0.75	9.14	4.57	−5.32	3.10
B₃S₃Cl₃	−9.86	−1.67	8.19	4.10	−5.76	4.05
B₃S₃CN₃	−10.75	−3.16	7.58	3.79	−6.96	6.38
B₃S₃Me₃	−9.09	−0.73	8.36	4.18	−4.91	2.88
B₃S₃H₃…B₂H₆	−9.65	−1.42	8.23	4.12	−5.54	3.72
B₃S₃F₃…B₂H₆	−9.88	−0.85	9.03	4.51	−5.37	3.19
B₃S₃Cl₃…B₂H₆	−9.86	−1.74	8.13	4.06	−5.80	4.14
B₃S₃(CN)₃…B₂H₆	−10.71	−3.19	7.52	3.76	−6.95	6.42
B₃S₃Me₃…B₂H₆	−9.17	−0.84	8.34	4.17	−5.00	3.00

Computed HOMO-LUMO gap values of these systems are given in Table 3. These values decreases X = F > H > Me > Cl > CN. This trend in understandable in terms of substituent character. In the presence of F substituent, LUMO energy raises (via σ* mixing or lone-pair antibonding interaction) whereas E(HOMO) slightly lowering. Therefore, largest HOMO-LUMO gap can be observed. In the presence of Me group, E(HOMO) raises and E(LUMO) raises less respect to X = H. So, smaller HOMO-LUMO gap is found than X = H. Cl lowers E(HOMO) moderately and lowers E(LUMO) significantly (via d-orbital acceptance). Thus, smaller HOMO-LUMO can be seen. Finally, CN decreases E(HOMO) and decreases E(LUMO) even more (via π*- acceptance). Hence, smallest HOMO-LUMO gap can be created.

Molecular electrostatic potential (MEP)

Now, we used MEP surfaces for exploring of EDA results. The MEP surfaces of the B₂H₆ and B₃S₃H₃ molecules can be help to exploring the binding sites in a rational approach for the B₃S₃H₃…B₂H₆ complex. MEP surfaces of these molecules are shown in Figure 5. The MEP maps obviously exemplify negative areas S atom of B₃S₃H₃ and a positive area around bridged hydrogen of B₂H₆ molecule. Therefore, interaction between two fragments is occurred between S atom of B₃S₃H₃ and bridged hydrogen of B₂H₆ molecule.

Figure 5.

The MEP surfaces of the B₂H₆ and B₃S₃H₃ molecules.

Charge decomposition analysis (CDA)

Now, we illustrated charge transfer between two fragments of B₃S₃X₃…B₂H₆ systems using of decomposition analysis (CDA) theory. Table 4 lists the CDA results of these systems. The CDA outcomes express larger electrons donated from B₃S₃X₃ fragment to the diborane molecule in the presence of X = H, F, Me. This variation can be attributed to better electron donating ability of these substituents. In contrast, greater electrons donated from diborane to B₃S₃X₃ fragment in presence of X = Cl, CN. This variation is compatible with to better electron withdrawing ability of these substituents.

Table 4.

Charge decomposition analysis (CDA) results of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at M062X/6-311G(d,p) level of theory.

Molecule	d	b	d–b	r
B₃S₃H₃…B₂H₆	0.004319	0.019351	−0.015031	−0.131100
B₃S₃F₃…B₂H₆	0.011732	0.028434	−0.016702	−0.138058
B₃S₃Cl₃…B₂H₆	0.005192	−0.006369	0.011561	−0.136704
B₃S₃(CN)₃…B₂H₆	0.020377	−0.003036	0.023413	−0.131560
B₃S₃Me₃…B₂H₆	0.003316	0.012329	−0.009013	−0.143154

d, b and r terms are the number of electrons donated from fragment B₂H₆ to fragment B₃S₃X₃, the number of electrons back donated from fragment B₃S₃X₃ to fragment B₂H₆, the number of electrons involved in repulsive polarization, respectively.

A minor (but non-negligible) value of r term for these interaction exposing a noticeable electron repulsion effect between the fragments.

QTAIM

Now, we provided worth idea about the bond critical points (BCP) of S…H of B₂H₆…B₃S₃X₃ molecules with QTAIM computations. Two bond critical points are accepted at S1…H and S2…H bonds. Computed ρ_BCP(S…H) values are recorded in Table 5. It can be observed the ρ_BCP(S…H) values decrease as X = Me > H > Cl > F > CN. This trend explains with variation of S…H distances. S…H distances decrease with increasing of ρ_BCP(S…H) values.

Table 5.

S…H bond distance (in Å) and QTAIM results of the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at M062X/6-311G(d,p) level of theory.

Molecule	r(S…H)	ρ (e.Å⁻³)	∇²ρ (e.Å⁻⁵)	V (e².Å⁻⁴)	G (e².Å⁻⁴)	H (e².Å⁻⁴)	\|V(r)\|/G(r)
B₃S₃H₃…B₂H₆	3.088	0.0064	0.0189	−0.0029	0.0038	0.0009	0.7528
B₃S₃F₃…B₂H₆	3.112	0.0062	0.0187	−0.0027	0.0037	0.0010	0.7375
B₃S₃Cl₃…B₂H₆	3.106	0.0062	0.0188	−0.0028	0.0038	0.0010	0.7438
B₃S₃(CN)₃…B₂H₆	3.163	0.0056	0.0174	−0.0025	0.0034	0.0009	0.7209
B₃S₃Me₃…B₂H₆	3.048	0.0069	0.0204	−0.0032	0.0041	0.0010	0.7662

Calculated ∇²ρ_BCP(S…H) values are recorded in Table 5. These positive values are well-matched with the S…H closed-shell interactions.

Calculated density of total electron energy (H, S…H) values are recorded in Table 5. The positive H values are reflected as a van der Waals interaction indicator.

Evaluated |V(r)|/G(r) ratio values are summarized in Table 5. |V(r)|/G(r) < 1 values are characteristic of van der Waals interaction.⁴¹

Interacting quantum atoms (IQA) approach

Interacting quantum atoms (IQA) method is valuable tool to exploring of metal-ligand interactions mechanisms in the B₂H₆…B₃S₃X₃ systems was explored with the atomic energy partitioning method, namely.^42–45 The overall interaction energy and electrostatic part of the interaction energy for the S-H bond in B₂H₆…B₃S₃X₃ complexes are presented in Table 6.

Table 6.

Decomposition of two-bodied interaction energies within the IQA framework for the B₃S₃X₃…B₂H₆ (X = H, F, Cl, CN, Me) systems at M062X/6-311G(d,p) level of theory (in a.u).

Molecule	$V_{n e}^{S H}$	$V_{e n}^{S H}$	$V_{e e}^{S H}$	$V_{n n}^{S H}$	$V_{C}^{S H}$	$V_{c l}^{S H}$	$V_{X C}^{S H}$	$E_{i n t}^{S H}$	$V_{X C}^{S H} / E_{i n t}^{S H}$
B₃S₃H₃…B₂H₆	−4.35	−2.94	4.66	2.74	4.66	0.12	−0.01	0.11	0.05
B₃S₃F₃…B₂H₆	−4.33	−2.91	4.63	2.72	4.63	0.11	−0.01	0.11	0.05
B₃S₃Cl₃…B₂H₆	−4.33	−2.92	4.63	2.73	4.63	0.11	−0.01	0.10	0.05
B₃S₃(CN)₃…B₂H₆	−4.26	−2.86	4.54	2.68	4.54	0.10	0.00	0.10	0.05
B₃S₃Me₃…B₂H₆	−4.40	−2.98	4.71	2.78	4.72	0.12	−0.01	0.11	0.05

$E_{i n t}^{S H}$ values

Positive $E_{int}^{SH}$ data express _that S…H interactions are repulsive in the B₂H₆…B₃S₃X₃ systems.

$V_{c l}^{S H}$ values

Positive $V_{cl}^{SH}$ values indicate the S…H interactions for B₂H₆…B₃S₃X₃ complexes are described by extremely great electrostatic repulsion. These values are all lesser in absolute value than the primary exchange interaction energy.

Exchange-interaction ratio (EIR)

The exchange-interaction ratio (EIR) is expressed as:

EIR = \frac{V_{XC}^{SH}}{E_{int}^{SH}}

Because $V_{XC}^{SH}$ is continually negative, the sign of EIR can be considered to recognize an interatomic interaction.

For 0 < EIR ≪ 1 (when $| V_{XC}^{AB} | ≪ | E_{int}^{AB} |$ ) an interaction is subjected by electrostatic terms (the lesser EIR the greater degree of iconic character is detected) however a value greater than 1 (when $| V_{XC}^{SH} | ≫ | E_{int}^{SH} |$ ) points at a mainly classical covalent bond.

Computed EIR values of S…H bonds in the B₃S₃X₃…B₂H₆ systems indicate degree of ionic character.

Conclusion

Quantum chemistry exploration of replacement effects of hydrogens of borthiin with X = F, Cl, Me, CN on the interaction with B₂H₆ at M062X/6-311G(d,p) level of theory exposed more significance interaction energetically, in the presence of methyl group. largest electronic spatial extent value was found in the presence of Cl substituent. The most polarity of these complexes was belonged to CN group. Molecular orbital analysis revealed decreasing of HOMO-LUMO gap values as: X = F > H > Me > Cl > CN. The MEP maps obviously exemplified negative areas S atom of B₃S₃H₃ and a positive area around bridged hydrogen of B₂H₆ molecule. The CDA outcomes expressed larger electrons donated from B₃S₃X₃ fragment to the diborane molecule in the presence of X = H, F, Me. In contrast, greater electrons donated from diborane to B₃S₃X₃ fragment in presence of X = Cl, CN. Computed EIR values of S…H bonds in the B₃S₃X₃…B₂H₆ systems indicated degree of ionic character.

Footnotes

ORCID iD

Reza Ghiasi

Ethical approval

Not applicable.

Consent to participate

All the co-authors consent to participate.

Consent to publish

All the co-authors consent to publish

Author's contribution

Reza Ghiasi: Supervision, Conceptualization, Data curation, Investigation, Formal analysis

Nahid Shajari: Methodology, Data curation, Writing, Review and Editing

Funding

The authors received no financial support for the research, authorship, and/or publication of this article.

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Data availability

Data will be made available on request.

References

Forstner

Mutterties

. Boron-sulfur ring compounds. Inorg Chem 1965; 5: 164.

Ghiasi

Monajjemi

. Theoretical study of borthiin and its derivatives: structure and aromaticity. J Sulfur Chem 2007; 28: 505–511.

Ghiasi

. Theoretical study of the properties of fluoroborthiin and fluoroboroxine. Journal of Molecular Structure: THEOCHEM 2008; 853: 77–81.

Ghiasi

Moghimi

. Theoretical study of the interactions between Borthiin and fluorinated Borthiins with difluorine. Phosphorus, Sulfur Silicon Relat Elem 2010; 185: 1964–1971.

Ghiasi

. Molecular interaction of H2 and H2O with borthiin: a theoretical study. Russ Chem Bull 2012; 61: 248–252.

Ghiasi

Tale

Daneshdoost

, et al. Nature of the interactions between B12N12 and inorganic aromatic molecules (B3N3H6, B3P3H6, B3O3H3, B3S3H3). Russ J Inorg Chem 2023; 68: 753–760.

Pimentel

Mc Clellan

. The hydrogen bond. San Francisco: W. H. Freeman, 1960.

Joesten

Schaad

. Hydrogen bonding. New York: Marcel Dekker, 1974.

Jeffrey

. An introduction to hydrogen bonding. New York: Oxford University Press, 1997.

10.

Scheiner

. Hydrogen bonding. A theoretical perspective. New York: Oxford University Press, 1997.

11.

Steiner

. The hydrogen bond in the solid state. Angew Chem, Int Ed Engl 2002; 41: 48–76.

12.

Grabowski

. Understanding. hydrogen bonds. London: Royal Society of Chemistry, 2020.

13.

Grabowski

. Hydrogen bonding - new insights. Amsterdam: Springer, 2006.

14.

Novoa

. Intermolecular Interactions in Crystals: Fundamentals of Crystal Engineering, 2017.

15.

Jabłoński

. Hydrogen bonds. Molecules 2023; 28: 1616.

16.

Scheiner

. Chemical Bond and Intermolecular Interactions. Open Educational Resources Universitas Airlangga, 2025.

17.

Desiraju

. Hydrogen bridges in crystal engineering: interactions without borders. Acc Chem Res 2002; 35: 565–573.

18.

Eisch

. Boron-carbon multiple bonds. Adv Organomet Chem 1996; 39: 355–391.

19.

Mesbah

Prasang

Hofmann

, et al. Dianions of tetraboranes(4): puckered aromatic four-membered rings and their reactions with conservation of aromaticity. Angew Chem, Int Ed 2003; 42: 1717–1719.

20.

Ziegler

Rauhut

. Accurate vibrational configuration interaction calculations on diborane and its isotopologues. J Phys Chem A 2019; 123: 3367–3373.

21.

Silva

Moreira

Costa

RFd

, et al. Elastic electron scattering by diborane(6) and diborane(4) molecules. Phys Chem A 2022; 126: 6710–6718.

22.

Chou

S-L

J-I

Peng

Y-C

, et al. Identification of cyc-B₃H₃ with three bridging B–H–B bonds in a six-membered ring. ACS Omega 2017; 2: 529–535.

23.

Hamilton

EJM

Kultyshev

, et al. A stacking interaction between a bridging hydrogen atom and aromatic π density in the n-B₁₈H₂₂–benzene system. Chem-Eur J 2006; 12: 2571–2578.

24.

Min

Shore

, et al. Nature of “hydrogen bond” in the diborane−benzene complex: covalent, electrostatic, or dispersive?. Inorg Chem 2007; 46: 3956.

25.

Tian

H-B

Bai

, et al.

Unique interactions between diborane and π orbitals: blue- or red-shifted hydrogen bonding?

J Phys Chem A 2008; 112: 8121–8128.

26.

Ravinder

Subramanian

. Substitution effects of diborane on the interaction with borazine (inorganic benzene). J Phys Chem A 2010; 114: 565–557.

27.

Liao

. Interpreting the electronic structure of the hydrogen-bridge bond in B2H6 through a hypothetical reaction. Struct Chem 2012; 23: 525–527.

28.

Ghiasi

Shajari

. Replacement effects of terminal hydrogens of diborane with fluorine on the interaction with borazanaphthalene. Results in Chemistry 2026; 27: 103225.

29.

Zhao

Truhla

. Comparative DFT study of van der waals complexes: rare-gas dimers, alkaline-earth dimers, zinc dimer, and zinc-rare-gas dimers. J Phys Chem 2006; 110: 5121–5129.

30.

Hay

. Gaussian Basis sets for molecular calculations. The representation of 3d orbitals in transition-metal atoms. J Chem Phys 1977; 66: 4377–4384.

31.

Krishnan

Binkley

Seeger

, et al. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J Chem Phys 1980; 72: 650–654.

32.

McLean

Chandler

. Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z=11–18. J Chem Phys 1980; 72: 5639–5648.

33.

Wachters

AJH

. Gaussian Basis set for molecular wavefunctions containing third-row atoms. J Chem Phys 1970; 52: 1033–1036.

34.

Frisch

Trucks

Schlegel

, et al. Wallingford, CT: Gaussian Inc., 2013.

35.

Chen

. Multiwfn: a multifunctional wavefunction analyzer. J Comp Chem 2012; 33: 580–592.

36.

. A comprehensive electron wavefunction analysis toolbox for chemists, multiwfn. J Chem Phys 2024; 161: 082503.

37.

Chen

. Atomic dipole moment corrected Hirshfeld population method. J Theor Comput Chem 2012; 11: 63.

38.

Xiao

. Generalized charge decomposition analysis (GCDA) method. J Adv Phys Chem 2015; 4: 111–124.

39.

Chen

. Quantitative analysis of molecular surface based on improved marching tetrahedra algorithm. J Mol Graph Model 2012; 38: 314–323.

40.

Keith

. TK Gristmill Software, AIMAll, Version 19.10.12. Overland Park, KS, USA, 2019.

41.

Cukrowski

Lange

JHd

MM . Physical nature of interactions in Zn^II complexes with 2,2′-bipyridyl: quantum theory of atoms in molecules (QTAIM), interacting quantum atoms (IQA), noncovalent interactions (NCI), and extended transition state coupled with natural orbitals for chemical valence (ETS-NOCV) comparative studies. J Phys Chem A, 2014; 118: 623–637.

42.

Blanco

Pandás

Mckee

. Bay-type H···H “bonding” in cis-2-butene and related species: QTAIM versus NBO description. J Comput Chem 2014; 35: 1499–1508.

43.

Pendas

Blanco

Francisco

. Chemical fragments in real space: definitions, properties, and energetic decompositions. J Comput Chem 2007; 28: 161–184.

44.

Blanco

Pendas

Francisco

. Interacting quantum atoms: A correlated energy decomposition scheme based on the quantum theory of atoms in molecules. J Chem Theory Comput 2005; 1: 1096–1109.

45.

Francisco

Pendas

Blanco

MAA

. A molecular energy decomposition scheme for atoms in molecules. J Chem Theory Comput 2006; 2: 90–102.