Metalloid–Organic Intermolecular Complexes with Charge State-Controlled Conformations

Shape alterations of molecular systems, induced by their (electric) charging/discharging, could facilitate useful electronic and/or mechanical functions in molecular-scale devices and machines. The present study reports structures, stabilities, charge distributions, and IR spectra for a group of complexes of a main-group metalloid (boron) atom with hydrocarbon molecules. The considered systems include the smallest species demonstrating the basic principle of operation, as well as their size-extended analogues, generalizing it to larger counterparts based on such units. The system geometries vary considerably between neutral and ionic counterparts and exhibit two–three typical conformations related to twisting by up to about 90 degrees. The predicted structures correlate with specific infrared spectra, which can enable their experimental identification and transformation tracking. The above-mentioned characteristics suggest the potential utility of such systems for intermolecular switches, with the possible spectral monitoring of their functioning.


Introduction
Metal-organic complexes are characteristic representatives of and are at the very heart of inorganic chemistry.They constitute an enormous variety of systems with numerous applications, such as novel materials with unique optical and magnetic properties, catalyses, and ways of utilizing energy (e.g., sun radiation), etc. Transition metals dominate in this area, while the present work contributes to the increasing volume of studies involving main-group metals.
Specifically, here, the focus is on boron complexes with the smallest and larger hydrocarbon molecules.Boron chemistry is a broad and versatile field, in particular due to B standing next to C in the periodic table; thus, organoboron species have properties that match and enrich those of carbon-based species.In particular, boron-carbon based compounds form a large family of systems, with these atoms being bonded to each other.Further, multiple chemical reactions deal with the formation, breaking, and/or transformation of such bonds.In particular, relevant reviews [1,2] (see also the references therein) describe the role of boron in organic synthesis, along with an explosive growth of its applications in materials and pharmaceutical science [3,4].For instance, some protease inhibitors containing boron are approved medications for dealing with viral infections such as HIV/AIDS, hepatitis C, and COVID-19.
One relatively recent example includes the minimal compounds of B and ethylene (C 2 H 4 ) and their shape alterations depending on the system charge state, i.e., upon ionization and electron attachment [5].Namely, for the B-(C 2 H 4 ) 2 ternary complex, the molecules attached on opposite sides of the atom would significantly change their relative orientation within the range from parallel (for cation) to intermediate (for neutral species) to perpendicular (for anion), with the pairs of B-C bonds twisting accordingly.So, the system behaves like a molecular structural switch controlled by electric charge, with hypothetical practical uses in molecular electronics and machinery, both being areas experiencing a hyperactive rise [6][7][8][9] (see the references therein as well).
The aim of the present work was to expand the above study and check the properties and behaviour of analogous systems with acetylene components (having triply bonded dicarbon, unlike in ethylene) as well as those of extended systems with larger molecules, such as those with H atoms replaced by CH 3 groups while keeping the same central structural unit (CC-B-CC).Specifically, here, the butyne (C 4 H 6 ) and butene (C 4 H 8 ) species are considered.The related goal was to test the generality of the charge-induced geometry variations for boron-hydrocarbon complexes with similar bond patterns.The present work confirms such common features of the small and larger systems and compares them in terms of stability, charge distributions, and IR spectra.The latter properties can also serve as a means for the experimental verification of the system formation and transformations.
This study could thus be considered as related to merging the fields of boron chemistry and molecular devices.An earlier directly relevant review [10] addressed, both experimentally and theoretically, reactions of B with neutral hydrocarbons (including methane, ethylene, acetylene, and others), leading to the formation of organoboron compounds.The present study extends the scope of the subject to other molecules, ternary systems, and ionic derivatives.

Results and Discussion
First, the smallest relevant systems are considered, with hydrogen-terminated (doubleand triple-bonded) C=C and C≡C units.Then, their methylated structural extensions with longer carbon skeletons are included.For each group, neutral and ionic species are investigated in terms of a set of parameters, as described below.

Minimal Systems
The smallest considered systems include binary and ternary complexes of B atoms with one or two acetylene and ethylene molecules.

Structures and Stabilities
A boron atom attaches to an acetylene molecule sideways (Figure 1), forming two B-C bonds at the expense of a (broken) C-C π-bond and apparently altering the carbon-atom hybridization from sp to (approximately) sp 2 , hence strongly bending C 2 H 2 away from B. This binary complex is considerably stable by a few electronvolts (Table 1).The attachment of a second C 2 H 2 unit symmetrically on the opposite side of B proceeds similarly and results in a structure with the CC units being considerably twisted relative to one another (see Table 1).As a preliminary interpretation, the latter feature could lower the mutual repulsion of the negatively charged CC units.The ternary B-(C 2 H 2 ) 2 complex shows a notable anticooperativity in stability, with the second added molecule increasing the binding by only a third (about 1.1 eV) relative to B-C 2 H 2 .Accordingly, the B-C bonds stretch, the C-C bonds shrink, and the molecules partially unbend (Table 1).Such nonadditivity is consistent with the binding ability of (monovalent) B, now shared between two attached molecules, and with their mutual electrostatic repulsion.
The above system geometries resemble those of the previously studied ethylene-based counterparts B-C 2 H 4 and B-(C 2 H 4 ) 2 [5].The former is significantly less bound (by about 1 eV) relative to B-C 2 H 2 , while the latter is bound equally compared to B-(C 2 H 2 ) 2 (Table 1).Accordingly, the B-C distances are longer, and the interaction nonadditivity is lower than for the acetylene-based species.The D e values for the binary species slightly exceed those obtained earlier, 3.06 eV for B-C 2 H 2 [11] and 2.02 eV for B-C 2 H 4 [12].
Unlike the above complexes with two identical molecules, the C 2 H 2 -B-C 2 H 4 one is predicted to be maximally twisted (Figure 1), so the attached molecules are perpendicular to each other.This is likely due to weaker C-C interactions between the CC units of different lengths.Still, the hetero-molecular complex is slightly more bound than either homo-molecular counterpart (       The above system geometries resemble those of the previously studied ethylenebased counterparts B-C2H4 and B-(C2H4)2 [5].The former is significantly less bound (by about 1 eV) relative to B-C2H2, while the latter is bound equally compared to B-(C2H2)2 (Table 1).Accordingly, the B-C distances are longer, and the interaction nonadditivity is lower than for the acetylene-based species.The De values for the binary species slightly exceed those obtained earlier, 3.06 eV for B-C2H2 [11] and 2.02 eV for B-C2H4 [12].
Unlike the above complexes with two identical molecules, the C2H2-B-C2H4 one is predicted to be maximally twisted (Figure 1), so the attached molecules are perpendicular to each other.This is likely due to weaker C-C interactions between the CC units of different lengths.Still, the hetero-molecular complex is slightly more bound than either homomolecular counterpart (Table 1).This stronger stabilization appears to be due to the increased B-C2H2 interaction, as follows from the B-C bonds being shorter, C-C bond longer, and CCH bending larger than in B-(C2H2)2, with opposite relations for the B-C2H4 fragment relative to that in B-(C2H4)2.Indeed, attaching C2H2 to B-C2H4 adds about 2.3 eV of binding as compared to about 1.4 eV for attaching C2H4 to B-C2H2.The corresponding ions of B-(C 2 H 2 ) 2 have analogous overall structures, with B sandwiched in between two molecules which, however, alter their mutual orientation.In the anionic and cationic complexes, the molecules are, respectively, perpendicular and parallel to one another (Figure 1).This could be qualitatively rationalized in terms of the molecules charged more and less negatively, respectively; accordingly, their mutual repulsion is stronger and weaker, leading to their increased and decreased relative twisting.
The stability increases from neutral B-(C 2 H 2 ) 2 to its anion to its cation.In particular, this is a general trend for all the above ternary complexes (Table 1).Curiously, for B-C 2 H 2 , the anion behaves similarly, while the cation is less stable than both the anion and neutral counterparts.As a result, the cationic complexes even show a cooperative nonadditivity of binding.The B-C distance increases for the anion and decreases for the cation in both binary and ternary complexes.The C-C distance in the CC units shows opposite variations for B-C 2 H 2 but, interestingly, similar variations for B-(C 2 H 2 ) 2 .
The ions of B-(C 2 H 4 ) 2 [5] repeat the shape changes exhibited by B-(C 2 H 2 ) 2 .Here, another correlation could be suggested to interpret the shape of the anion, namely B-(C 2 H 4 ) 2 − being isoelectronic to spiropentane (with central C instead of B) which is known to have perpendicular C 2 H 4 units [13].Curiously, an analogous consideration of anionic B-(C 2 H 2 ) 2 seems to fail, since isoelectronic spiropentadiene differs in shape (with parallel C 2 H 2 units).The stabilization relative to neutral B-(C 2 H 4 ) 2 is stronger for the anion and weaker for the cation compared to the B-(C 2 H 2 ) 2 case, so the two ternary ions are nearly equally stable (Table 1).Additionally, now it is the ternary anion which is bound cooperatively.The B-C bonds shrink for the anion and stretch for the cation, opposite to the B-(C 2 H 2 ) 2 case.The C-C bonds change in counterphase and are thus similar to those in B-(C 2 H 2 ) 2 .The B-C 2 H 4 complex is also stabilized by anionization, while ionization reshapes it via the insertion of B inside the stretched CC unit (making a linear CBC core).The latter system is therefore not considered further here.For B-C 2 H 4 , anionization increases the B-C bond length and decreases the C-C bond length, similar to the variations in B-C 2 H 2 but oppositely to those in B-(C 2 H 4 ) 2 .
The C 2 H 2 -B-C 2 H 4 hybrid, upon ionization, transforms similar to the homo-molecular counterparts, with the molecules becoming parallel to one another, while anionization leaves the shape of the complex (with mutually perpendicular molecules) unchanged (Figure 1).The system stability rises from neutral to anion to cation, thus reproducing the variations in B-(C 2 H 2 ) 2 , even though the increments are somewhat smaller.The related B-C and C-C distance alterations in the C 2 H 2 and C 2 H 4 components follow those in B-(C 2 H 2 ) 2 and B-(C 2 H 4 ) 2 , respectively, except for the C-C distance in C 2 H 2 , decreasing from the neutral species to the anion.
For the binary species, the BSSE correction reduces the dissociation energies by about 0.5 ± 0.1 eV, being maximal for the anions.The values increase to about 0.6-0.7 eV for the ternary systems.

Charge Distributions
The B atom is predicted to be positively charged in all neutral systems studied (Table 2) within the range of 0.45-0.70e.The charge on B is smallest for the acetylene-based and largest for the ethylene-based species and slightly reduces from the binary to ternary complexes.The latter feature can be related directly to the competition between the two attached molecules (pulling the charge in opposite directions) and inversely to the B-C separations (longer in the ternary systems) determining the charge transfer distance.
In the respective anions and cations, the B atom becomes, respectively, less and more positive, up to near-zero or even negative charges in the binary anionic complexes and up to a full unit charge in cationic B-C 2 H 2 .Curiously, B-(C 2 H 2 ) 2 shows a counter-intuitive exception here, with B being less positive for the cation than the anion.

Ionization Energies and Electron Affinities
The energy differences between the cation and neutral system and between the latter and anion conventionally define the ionization energy (IE) and electron affinity (EA), respectively.Their calculated values for the present systems are presented in Table 3.Such values for the B atom are reported as well and can be seen to match well with the experimental data, offering further validation of the employed theoretical approach.The IE values are lower than IE(B) for all systems except B-C 2 H 2 .This could be readily related to the D e values being larger for the cations than for the corresponding neutrals, while this relation is opposite only for B-C 2 H 2 (Table 1).The IEs are lower for the ternary systems compared to the binary systems.Opposite to the situation for the binary species, the acetylene-based ternary system has lower IE compared to the ethylene-based one, while the value for the hetero-molecular case is intermediate between the two homo-molecular counterparts.
In turn, the EA values for all systems are considerably higher than EA(B), consistent with the higher stabilities of the anions relative to the neutrals.The ternary systems are about triply more electrophilic than the binary ones.The ethylene-based species have higher EAs compared to the acetylene-based ones, while the hetero-molecular system has a lower EA compared to either homo-molecular analogue.

Simulated IR Spectra
Upon the formation of a binary complex with B, the intense high-frequency (C-H stretch related) IR band of C 2 H 2 is apparently cancelled (Figure 2), while its low-frequency (bending) doublet is expanded by a group of blue-shifted (by up to 500 cm −1 ) lines.The most intense one of the latter (near 1200 cm −1 ) corresponds to a B-C 2 H 2 axial stretch.For C 2 H 4 , the dominating low-frequency (C-C-H bending) band splits into anti-and symmetric modes and red-shifts upon such a complexation.
The addition of a second acetylene or second ethylene unit leads to the appearance of a dominant mid-range (near 1600 cm −1 ) band associated with an antisymmetric C 2 H 2/4 -B-C 2 H 2/4 axial stretch, while the lower-frequency part of the spectrum is intensified and red-shifted in the former but suppressed in the latter case.The hetero-molecular ternary complex in effect mainly preserves the most intense band of B-C 2 H 2 (near 1200 cm −1 ) and even somewhat recovers its (red-shifted) C-H stretch component.
The anionization of B-(C 2 H 2 ) 2 produces a very intense band near 1300 cm −1 due to (red-shifted) antisymmetric axial stretching and strong high-frequency (antisymmetric C-H stretching) bands (Figure 3).The latter also dominate for cationic B-(C 2 H 2 ) 2 , which preserves the low-frequency group of lines but not the band dominating in the neutral counterpart.The predominant band of B-(C 2 H 4 ) 2 is also removed in both ions, and their spectra exhibit intense high-and low-frequency bands for the anion and cation, respectively.For the cation, the most intense band, near 800 cm −1 , corresponds to the (strongly redshifted) antisymmetric axial stretch.
C2H4, the dominating low-frequency (C-C-H bending) band splits into anti-and symmetric modes and red-shifts upon such a complexation.
The addition of a second acetylene or second ethylene unit leads to the appearance of a dominant mid-range (near 1600 cm −1 ) band associated with an antisymmetric C2H2/4-B-C2H2/4 axial stretch, while the lower-frequency part of the spectrum is intensified and red-shifted in the former but suppressed in the latter case.The hetero-molecular ternary complex in effect mainly preserves the most intense band of B-C2H2 (near 1200 cm −1 ) and even somewhat recovers its (red-shifted) C-H stretch component.The anionization of B-(C2H2)2 produces a very intense band near 1300 cm −1 due to (red-shifted) antisymmetric axial stretching and strong high-frequency (antisymmetric C-H stretching) bands (Figure 3).The latter also dominate for cationic B-(C2H2)2, which preserves the low-frequency group of lines but not the band dominating in the neutral counterpart.The predominant band of B-(C2H4)2 is also removed in both ions, and their spectra exhibit intense high-and low-frequency bands for the anion and cation, respectively.For the cation, the most intense band, near 800 cm −1 , corresponds to the (strongly red-shifted) antisymmetric axial stretch.

Extended Systems
The following systems are structural derivatives of the above minimal species, with two hydrogen atoms replaced by methyl radicals.

Structures and Stabilities
Both B-C4H6 and B-(C4H6)2 systems geometrically resemble the acetylene-based counterparts (Figure 4), including a considerable twist (by about same angle) between the mid-

Extended Systems
The following systems are structural derivatives of the above minimal species, with two hydrogen atoms replaced by methyl radicals.

Structures and Stabilities
Both B-C 4 H 6 and B-(C 4 H 6 ) 2 systems geometrically resemble the acetylene-based counterparts (Figure 4), including a considerable twist (by about same angle) between the middle CC units in the ternary system, as expected due to the identical B-molecule bonding scheme.Accordingly, the binding energies and the B-C and (middle hereafter) C-C bond lengths are almost unchanged relative to the smaller counterparts (Table 4).In particular, this is also consistent with a similar invariability of stability from boronacetylene to boron-methylacetylene when a single H is replaced by CH 3 , as predicted earlier [14].The dimethylacetylene (coinciding with the presently studied C 4 H 6 ) has also been previously studied in terms of making compounds with B [15], including its structure (shown in Figure 4), though with no stability characterization.
The corresponding ions are also structurally similar to the smaller analogues, with the molecules being perpendicular and parallel in the anionic and cationic ternary complexes, respectively.The binding energies increase from neutral to anionic to cationic complexes for both one and two attached molecules, stronger for the latter case and leading to almost additive stabilization for the ternary cation.This generally repeats the case of neutral counterparts except for the opposite relation between the relative stabilities of the binary ionic complexes.The B-C and C-C bonds are slightly longer in the ternary compared to the respective binary systems and vary with the charge state similarly for both cases.
Similarly, B-C 4 H 8 and B-(C 4 H 8 ) 2 also have geometries analogous to those of the ethylene-based species.Here, the ternary system exhibits two different conformers, with two attached molecules oriented in phase or antiphase (Figure 5), whose cases are labelled cis and trans, respectively.The middle CC units are twisted at an angle which is close to that in B-(C 4 H 6 ) 2 for the cis but larger for the trans conformer (Table 4).Visually, the overall CCCC skeletons look more aligned for the trans case.The cis and trans conformers are nearly equally stable to dissociation, while both these and the binary complex are slightly less stable (by about 0.3-0.4eV) than the respective ethylene-based counterparts.While B-C 4 H 8 is considerably less stable (by about 1.3 eV) than B-C 4 H 6 , the difference reduces to about 0.4 eV for the corresponding ternary species, similar to the smaller systems.The B-C and C-C bonds are again slightly stretched and shrunk, respectively, in the ternary relative to binary systems.The trans case, however, shows some asymmetry, with one molecule being slightly farther away from B and, accordingly, shorter than the other.
The respective ions of the cis conformer reproduce the relative orientations of the middle CC units of the molecules found for the smaller (ethylene-based) counterpartsperpendicular for the anion and parallel for the cation-while for the trans case, the angles are slightly different (Table 4).Accordingly, the overall orientations of the CCCC skeletons of the ions are also perpendicular and parallel for the cis but opposite for the trans cases (Figure 5).Unlike for the smaller (ethylene-based) counterparts, the trans anion remains almost unchanged from the corresponding neutral species.Similar to the C 4 H 6 -based species, the binding energies increase from the neutral system to anion to cation more strongly for the ternary complexes (and equally for the cis and trans cases).The ternary anion exhibits a cooperative binding, similar to the ethylene-based analogue.While the cation of B-C 4 H 8 features B inserted into the middle CC unit, similar to the corresponding smaller (ethylene-based) system, being therefore also dropped from further consideration.In anionic B-(C 4 H 8 ) 2 , the attached molecules somewhat converge towards B and stretch, oppositely to the variations in anionic B-C 4 H 8 , and then evolve reversely in the ternary cation.This is similar to the binary C 4 H 6 -based complexes but different from the ternary ones.These alterations resemble those for the smaller analogues.
In the hetero-molecular C 4 H 6 -B-C 4 H 8 system, the middle CC units are near-perpendicular (Figure 6), as in the smaller counterpart (with acetylene and ethylene).The complex is more stable than its homo-molecular analogues (Table 4), also similar to the situation for the respective smaller species.The C 4 H 6 unit is closer than The corresponding ions are also structurally similar to the smaller analogues, with the molecules being perpendicular and parallel in the anionic and cationic ternary complexes, respectively.The binding energies increase from neutral to anionic to cationic complexes for both one and two attached molecules, stronger for the latter case and leading to almost additive stabilization for the ternary cation.This generally repeats the case of neutral counterparts except for the opposite relation between the relative stabilities of the binary ionic complexes.The B-C and C-C bonds are slightly longer in the ternary compared to the respective binary systems and vary with the charge state similarly for both cases.In the hetero-molecular C4H6-B-C4H8 system, the CC units are near-perpendicular (Figure 6), as in the smaller counterpart (with acetylene and ethylene).The complex is more stable than its homo-molecular analogues (Table 4), also similar to the situation for the respective smaller species.The C4H6 unit is closer than C4H8 to B, and for the former, the B-C and C-C bonds are, respectively, shorter and longer than in B-(C4H6)2, while for the latter, they are intermediate between those in cis and trans B-(C4H8)2.Both B-C and C-C bonds are also marginally longer than in the smaller counterpart (C2H2-B-C2H4).The anionic derivative of C4H6-B-C4H8 approximately preserves the near-perpendicular orientation of the middle CC units, while they become near-parallel in the cation (Figure 6).Similar to the homo-molecular analogues, the stabilities increase upon anionization and more so upon ionization (Table 4).For both ions, the binding energies are intermediate between those for B-(C4H6)2 and B-(C4H8)2, unlike for the neutral systems.The C4H6 unit slightly withdraws from B in the anion and returns even closer in the cation, while C4H8 behaves in counterphase while remaining farther from B.
The BSSE energies are about 0.4 eV for all binary systems, comparable to those for the minimal counterparts.The values double for the ternary systems to about 0.7-0.8eV.

Charge Distributions
The extended systems show general similarities in the distributions of charges as compared to the smaller analogues (Table 2).In particular, the B atom is positive in all systems except B-C4H6 (where it is slightly negative) and B-C4H8 (where it is essentially neutral).The largest charge on B is found in cationic B-(C4H8)2 (both conformers) and B- The anionic derivative of C 4 H 6 -B-C 4 H 8 approximately preserves the near-perpendicular orientation of the middle CC units, while they become near-parallel in the cation (Figure 6).Similar to the homo-molecular analogues, the stabilities increase upon anionization and more so upon ionization (Table 4).For both ions, the binding energies are intermediate between those for B-(C 4 H 6 ) 2 and B-(C 4 H 8 ) 2 , unlike for the neutral systems.The C 4 H 6 unit slightly withdraws from B in the anion and returns even closer in the cation, while C 4 H 8 behaves in counterphase while remaining farther from B.
The BSSE energies are about 0.4 eV for all binary systems, comparable to those for the minimal counterparts.The values double for the ternary systems to about 0.7-0.8eV.

Charge Distributions
The extended systems show general similarities in the distributions of charges as compared to the smaller analogues (Table 2).In particular, the B atom is positive in all systems except B-C 4 H 6 (where it is slightly negative) and B-C 4 H 8 (where it is essentially neutral).The largest charge on B is found in cationic B-(C 4 H 8 ) 2 (both conformers) and B-C 4 H 6 , although the smaller (ethylene-and acetylene-based) counterparts show inverse relative B charges between the binary and ternary species.The extended systems consistently exhibit higher B charges, with the only exception being the B-C 4 H 6 vs. B-C 2 H 2 cations.In all extended ternary complexes containing C 4 H 8 , the charge on B expectedly decreases from neutral to anion and then increases in the cation, while the variation is opposite for the B-(C 4 H 6 ) 2 case, similar to the charge evolutions in the smaller analogues.

Ionization Energies and Electron Affinities
The IE and EA values for the extended systems are lower than for the corresponding smaller counterparts, except for the opposite relation between B-C 4 H 8 and B-C 2 H 4 (Table 3).The relations between the values for the extended systems reproduce those for their smaller analogues.

Simulated IR Spectra
The calculated IR intensity spectra for C 4 H 6 (butyne) and C 4 H 8 (butene) exhibit a few distinct bands, including the most intense ones near 3000 cm −1 and duller ones near 1500 and (for butene) 1000 cm −1 .The and relative intensities of these bands match well with the available experimental data [13].
The attachment of B weakens the high-frequency (C-H stretching) band for C 4 H 6 while adding weaker satellites to it for C 4 H 8 (Figure 7).The new bands appear, including the most intense one near 1200 cm −1 for the butyne complex (corresponding to its axial stretch) and the bands between 500 and 1000 cm −1 for the butene complex, with the latter effectively red-shifting the original (C-C-H bending) band of the free molecule near 1000 cm −1 .
In the ternary systems, the high-frequency band recovers its intensity in the di-butyne complex, and new intense bands (stronger and slightly blue-shifted for the cis di-butene complex) appear near 2000 cm −1 due to the antisymmetric-stretching axial vibrations.In addition, the di-butyne complex acquires a group of lines between 500 and 1000 cm −1 , while a similar part of the spectrum diminishes in the cis di-butene complex.The latter is also valid for its trans counterpart, which, however, lacks an intense mid-range band, perhaps due to the asymmetry, which allows one to distinguish the two conformers.The IR spectrum of the hetero-molecular complex is dominated by an intense (asymmetric-stretch) band near 1300 cm −1 , which differentiates this system from both homo-molecular counterparts, and is contributed by the high-frequency bands.For clarity, the dominating band is associated with the system's axial oscillations between the following formal structures: Molecules 2024, 29, 1635 13 of 17 The anionization of both homo-molecular ternary complexes considerably focuses their IR spectra, which are dominated by narrow high-frequency (C-H stretching) bands (Figure 8).The di-butyne system also shows a comparably intense band near 1300 cm −1 , associated with the (red-shifted) antisymmetric-stretch mode.The cations, instead, preserve and intensify the group of bands in the lower frequency range between 500 and 1500 The anionization of both homo-molecular ternary complexes considerably focuses their IR spectra, which are by narrow high-frequency (C-H stretching) bands (Figure 8).The di-butyne system also shows a comparably intense band near 1300 cm −1 , associated with the (red-shifted) antisymmetric-stretch mode.The cations, instead, preserve and intensify the group of bands in the lower frequency range between 500 and 1500 cm −1 , including the antisymmetric-stretching band near 1000 cm −1 , the most intense (and still further red-shifted) band for the di-butyne case.The di-butene systems add another intense band near 800 cm −1 (due to an antisymmetric stretch as well, but with no accompanying CH 3 group rocking), which is most intense for the cis conformer.For the trans di-butene conformer, the spectra of the ions are very similar to those for the cis counterpart (especially for the anions), except for the interchanged relative intensities of the two intense (low-frequency) bands for the cations.

Computational Methods
The present computational study addresses a series of polyatomic intermolecular systems with components interacting both covalently and (especially for extended structures) noncovalently.Accordingly, a consistent MP2 (second-order Moller-Plessett perturbation theory) level approach with aug-cc-pVTZ basis sets [16] for all atoms has been used as a reasonable compromise between the computing time and accuracy (including the reliable handling of relevant dispersion interactions and considerable charge transfer).In particular, spin-unrestricted Hartree-Fock calculations have been used for open-shell systems.Basis set superposition error (BSSE) corrections have been applied via a conventional counterpoise technique [17].In particular, previous calculations of identical accuracy for (composition-analogous) main-group-metal/benzene complexes [18] showed a very close match to corresponding results produced at a (more accurate) CCSD-T level.In addition, the most relevant (B-C) interactions are obviously involved in the BC diatom, which has been characterized and tested thoroughly against experimental data previously [5].The above theoretical framework has been employed by means of the NWChem ab initio suite of programs [19].
All system geometries have been fully optimized, and energy minima have been verified via all-real vibrational frequencies.Encountered transition states have been reoptimized using corresponding eigenvectors.For extended systems, optimizations have been carried out with the (reduced) aug-cc-pVDZ basis sets, followed by single-point energy calculations with the above larger basis sets.
Related IR intensity spectra have also been evaluated using NWChem.In particular, test calculations for constituent small molecules have exhibited a good match to available experimental data (with intense bands near 700 and 3400 cm −1 for C2H2 and intense bands near 1000,

Computational Methods
The present computational study addresses a series of polyatomic intermolecular systems with components interacting both covalently and (especially for extended structures) noncovalently.Accordingly, a consistent MP2 (second-order Moller-Plessett perturbation theory) level approach with aug-cc-pVTZ basis sets [16] for all atoms has been used as a reasonable compromise between the computing time and accuracy (including the reliable handling of relevant dispersion interactions and considerable charge transfer).In particular, spin-unrestricted Hartree-Fock calculations have been used for open-shell systems.Basis set superposition error (BSSE) corrections have been applied via a conventional counterpoise technique [17].In particular, previous calculations of identical accuracy for (composition-analogous) main-group-metal/benzene complexes [18] showed a very close match to corresponding results produced at a (more accurate) CCSD-T level.In addition, the most relevant (B-C) interactions are obviously involved in the BC diatom, which has been characterized and tested thoroughly against experimental data previously [5].The above theoretical framework has been employed by means of the ab initio suite of programs [19].
All system geometries have been fully optimized, and energy minima have been verified via all-real vibrational frequencies.Encountered transition states have been reoptimized using corresponding eigenvectors.For extended systems, optimizations have been carried out with the (reduced) aug-cc-pVDZ basis sets, followed by single-point energy calculations with the above larger basis sets.
Related IR intensity spectra have also been evaluated using NWChem.In particular, test calculations for constituent small molecules have exhibited a good match to available experimental data (with intense bands near 700 and 3400 cm −1 for C 2 H 2 and intense bands near 1000, 1500, and 3200-3300 cm −1 for C 2 H 4 ).Further tests for extended molecules have also compared well to relevant experiments, as discussed in Section 2.
Charges on atoms have been obtained using the NPA (natural population analysis) approach [20], which has been implemented in the Janpa software (http://janpa.sourceforge.net)[21].

Conclusions
A series of complexes of a B atom with one and two hydrocarbon molecules were investigated at a consistent theoretical level (MP2/aug-cc-pVTZ + BSSE) in terms of structures, stabilities, charge distributions, and IR spectra.These properties were also predicted for the anionic and cationic derivatives of the systems, with attention being paid to their reshaping upon anionization/ionization.The molecules studied included minimal dicarbon species such as C 2 H 2 and C 2 H 4 , as well as their tetracarbon extensions (with methyl groups replacing two hydrogens) C 4 H 6 and C 4 H 8 .Both homo-and hetero-molecular ternary complexes were considered, with two conformers for the C 4 H 8 case.
The B atom can bind sideways to the CC unit (middle one for the larger species), with the second molecule attaching on the other side (symmetrically in the homo-molecular case).The related re-hybridization of the C atoms bends the molecules.The extended systems are generally similar to the corresponding minimal ones, although with some deviations.In most ternary complexes, the binding is anti-cooperative, with the second molecule being bound less strongly, except for some ions.The significant binding of about 2-4 eV per molecule suggests a ready formation of such systems, at least as direct kinetic products, even if in competition with an isomerization of the binary species, e.g., with H transfer to B [10,22], apparently over significant energy barriers.
In the neutral ternary systems, the molecules are oriented at an intermediate angle (mainly around 60 • ) to one another, except for in the hetero-molecular species, where the molecules are mutually perpendicular.Electron attachment turns the attached molecules into a perpendicular relative orientation (or keeps it), while ionization makes them parallel, which applies to the middle CC units in the tetracarbon molecules case.The same relative orientations transfer to the entire C 4 H 6 components, while the orientations can alter (up to inversion) for the C 4 H 8 ones due the non-flatness of the carbon skeleton in these complexes.
The above transformations are readily reversible for the neutral-anion pair, with the neutral-species geometry recovering upon electron detachment from the corresponding anion.The predicted EA values are moderate (about 1-2 eV), suggesting the charge/discharge to be feasible, and the anions are somewhat more stable to dissociation than the corresponding neutrals, which supports their lifetime.The latter feature is even more favourable for the cations, while the much larger IE values (about 5-6 eV for the extended ternaries) could make such charge/discharge more problematic.Additionally, "returning" an electron to cation normally produces an excited-state species (up to Rydberg one) rather than a ground-state species, and there may well be considerable differences between the geometries of these two states.The relaxation from the excited to the ground state may be complicated, even though one could imagine this to possibly occur, in particular, via a light emission, which would add an interesting, potentially usable (optical effect) behaviour to such a hypothetical switch.Or, perhaps a partial and/or gradual and thus more reversible charge-transfer might be a more realistic option here, maybe via an interaction with a electronegative entity.At the very least, applications could be limited to neutral-anion switching despite the reduced range of shape alteration.
The simulated IR spectra vary significantly from the free molecules to the binary and ternary complexes (even between conformers) to their ions and can be used for the experimental identification of the species formation and their transformations with the charge state.The variations include the suppression or recovery of the high-(near and above 3000 cm −1 ) and low-frequency bands (around 1000 cm −1 ) characteristic of the free molecules, their spectral shifts and formation of closely located satellites, and the appearance (for homo-molecular ternary systems) or diminishing (in their ions) of new bands in the mid-range (near 2000 cm −1 ).In particular, the IR spectra intensify for ternary homo-molecular complexes and their ions (especially anions).
The above ternary complexes thus appear to be prototype intermolecular systems, by themselves or as units within larger structures or frameworks, capable of offering charge-controlled shape alterations potentially usable for molecular electronics and devices/machines.In addition, their related switching among different geometries (within 90 • for relative orientations of the molecular components) could possibly be tracked in terms of rather sensitive IR spectra.

Table
).This stronger stabilization appears to be due to the increased B-C 2 H 2 interaction, as follows from the B-C bonds being shorter, C-C bond longer, and CCH bending larger than in B-(C 2 H 2 ) 2 , with opposite relations for the B-C 2 H 4 fragment relative to that in B-(C 2 H 4 ) 2 .Indeed, attaching C 2 H 2 to B-C 2 H 4 adds about 2.3 eV of binding as compared to about 1.4 eV for attaching C 2 H 4 to B-C 2 H 2 .

Table 1 .
Equilibrium parameters (dissociation energies, distances, angles) of the minimal systems.

Table 4 .
Equilibrium parameters (dissociation energies, distances, angles) of the extended systems.

Table 4 .
Equilibrium parameters (dissociation energies, distances, angles) of the extended systems.

Charge State D e /eV R e (B-C)/Å R e (C-C) 1 φ e 2 / • θ e (CCCC) 3
unit (bound to B).2CCC bond angle in C 4 H 6 ; CCCC axial-folding angle in C 4 H 8 . 3twisting angle between two opposite middle C-C units.