Step-by-Step Replacement of Cyano Groups by Tricyanovinyls—The Influence on the Acidity

Acid-base properties are the simplest expression of compounds’ coordinating ability. In the present work, we studied in silico how the gas-phase Brønsted acidity (GA) of several polycyano-substituted compounds change when cyano (CN) groups are replaced by 1,2,2-tricyanovinyl (TCNV) groups in (iso)cyanic acid, dicyanoamine, cyanoform, and hydrogen tetracyanoborate. Different tautomers and conformers/isomers are included in this study. Gas-phase acidity values are compared with the acidities of various acids, including percyanated protonated monocarba-closo-dodecaborate (carborane acid) and dodecaborate, as well as hydrogen cyanide and 1,2,2-tricyanoethene. An estimation of acetonitrile (MeCN), dimethylsufoxide (DMSO), and 1,2-dichloroethane (DCE) acidities is presented using the COSMO-RS method and correlation analysis. The strongest acid with four TCNV groups shows remarkable acidic properties.


Introduction
The gas-phase acidity (GA) of a neutral acid AH and proton affinity (PA) of the corresponding anionic base A -refers to the following equation: PA of an anionic base A − is defined as the enthalpy change (∆H acid ) for the gas-phase reaction 1 between a proton and the anion to produce the neutral conjugate acid of the anion.GA of an acid HA is the Gibbs free energy change (∆G acid ) on deprotonation of the acid according to Equation (1).GA values include the entropy factor.GA is directly measurable and comparable with experimental values, and therefore they will be used in discussions in this work.Reaction 1 shows that the smaller the PA or GA value, the more easily the proton is detached from the acid AH; therefore, the stronger the acid is.Acids with GA values around 400-350 kcal/mol can be considered weak acids, acids having GA values lower than 300 kcal/mol are superacids, while acids with a GA value lower than 250 kcal/mol are referred to as hyperacids.
Acid dissociation constant K a (usually expressed as its negative logarithm pK a ) is an equilibrium constant that shows how strong an acid AH is according to Equations ( 2) and (3).S refers to the solvent , A − refers to an acid anion, and a to an activity.The stronger the acid, the lower the pK a value.The pK a value depends strongly on the solvent [1,2].According to the definition, a superacid is an acid that is stronger than sulfuric acid in the respective solvent/medium.pK a values of superacids are lower than 9 in MeCN, 42 in DCE [3], and around 1.4 pK a units in DMSO [4].
Cyanocarbon chemistry gained its first insights quite some time ago with the seminal works by Webster, Middleton, and Engelhardt introducing different polycyanated compounds such as polycyano cyclopentadienides [5], pentacyanopropenide [6], and other compounds with tricyanovinyl substituents [7].Already back then, the authors realized that conjugate acids of anionic cyanocarbon compounds are very strong acids.
Although cyanocarbon acids are not applicable as pure acids, the anions still find plenty of recognition.Computational acidity studies of cyanocarbons are not rare, and have been conducted mainly by Vianello and Maksić [8][9][10][11][12], Leito and Koppel [13,14], and others [15].The high acidifying effect of the cyano (CN) group has been used to probe the boundaries of the conventional acidity of superacids.It has been shown that the acidity of different hydrocarbons can be increased by more than 100 kcal/mol by substituting them extensively with CN groups [8].The strongest organic cyanocarbon acids reach the gas-phase acidity of 230 kcal/mol.
Substitution with multiple cyano groups has such an acidifying effect because the anion formed after deprotonation is highly stabilized by efficient charge delocalization.The negative charge is attracted by the CN groups both via strong inductive and mesomeric (resonance) effects.Moreover, compared to -NO 2 or -SO 2 CF 3 groups, the CN group is very compact, i.e., it has low steric demand.Because the CN group is stick-like, it cannot be "twisted out of plane" of an aromatic system, which is a typical cause of weakening of the acidifying effect of the nitro group, for example.Thus, it is possible to put numerous CN groups into one molecule without creating extensive steric repulsion.For example, pentacyanophenol (pK a (MeCN) = 4.2) is one of the strongest acids among phenols, leaving 2,4,6-(SO 2 CF 3 ) 3 -phenol behind by 1 pK a unit and 2,4,6-trinitrophenol by 7 pK a units in acetonitrile (MeCN).Inserting five-SO 2 CF 3 or nitro groups into the benzene moiety proves impossible.
However, the CN group also has some disadvantages when it comes to creating strong acids.From the acidity and coordination perspective, the lone pair of nitrogen attracts all nucleophiles, including the proton.It has been shown computationally that the monocarba-closo-dodecaborane acid (CB 11 H 12 H) becomes weaker when one of the hydrogens is substituted with a CN group [16].However, replacing all 12 positions in the carborane moiety with the CN groups leads to an enormous acidity increase by 50 kcal/mol.Pentacyanocyclopentadienide C 5 (CN) 5 − has been a notorious anion that cannot be protonated even with perchloric acid [5].Later, it was shown by C. A. Reed that pentacyanocyclopendadiene HC 5 (CN) 5 has a polymeric structure where anions are bridged with protons [17].The black-brown solid that was described to form is common when trying to prepare the pure acid from polycyano anions.
The 1,2,2-tricyanovinyl group (TCNV) accommodates three CN groups.From the point of view of the original acidity center, all CN groups are further away; therefore, inductively, they have a lower influence on the acidity.However, the resonance possibility creates additional acidity centers with the adjacent CN groups.The negative charge of the anion has the possibility to further delocalize, stabilizing the anionic and giving the possibility to create stronger acid compared to the system with only CN groups.
Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.
Table 1.Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Initial Compound
Step 1 Step 2 Step 3 Step 4 Isocyanic acid 1H, dicyanamine 3H, cyanoform 6H [18], and tetracyanoborate 10 are all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.

Initial Compound
Step 1 Step 2 Step 3 Step 4 Isocyanic acid 1H, dicyanamine 3H, cyanoform 6H [18], and tetracyanoborate 10 are all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.

Initial Compound
Step 1 Step 2 Step 3 Step 4 Isocyanic acid 1H, dicyanamine 3H, cyanoform 6H [18], and tetracyanoborate 10 are all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.

Initial Compound
Step 1 Step 2 Step 3 Step 4 Isocyanic acid 1H, dicyanamine 3H, cyanoform 6H [18], and tetracyanoborate 10 are all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.

Initial Compound
Step 1 Step 2 Step 3 Step 4 Isocyanic acid 1H, dicyanamine 3H, cyanoform 6H [18], and tetracyanoborate 10 are all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.

Initial Compound
Step 1 Step 2 Step 3 Step 4 all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.
Table 1.Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Initial Compound
Step 1 Step 2 Step 3 Step 4 all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.
Table 1.Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Initial Compound
Step 1 Step 2 Step 3 Step 4 all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.
Table 1.Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Initial Compound
Step 1 Step 2 Step 3 Step 4 all well-known compounds and have been chosen as the base structures of current work.Tetracyanoborate 10 is a weakly coordinating anion [19,20], a component for ionic liquids [21][22][23], and its acid form is stable in the diluted solution [3].The compounds substituted by TCNV groups could be prepared from the readily available tetracyanoethylene (TCNE).Compound 2 can be prepared from TCNE by hydrolysis [6], 4 with ammonium acetate [24], and 5 in liquid ammonia [6].Compound 7 has been prepared by combining TCNE and malononitrile [6].Synthesis methods for 8 are available by Webster [7].Compound 9 is not yet prepared, but the synthesis could be envisaged according to published procedures [25].Borates with TCNV groups 11- 14 are not yet prepared.

Compounds under Study
Four base structures-(iso)cyanic acid 1H, dicyanamine 3H, cyanoform 6H, and hydrogen tetracyanoborate 10H-are denoted as initial compounds (either in acid or anionic form).The CN groups of these compounds were stepwise replaced by TCNV groups (steps 1 to 4), as shown in Table 1.
Table 1.Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Initial Compound
Step 1 Step 2 Step 3 Step 4 In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.
In the case of each tautomer, different conformers and isomers (diastereomers, Figure 2) were included in the calculations.The TCNV group is always planar in the case of anions and tautomers T1, T2, and T4-T5.For anion and tautomers T1, T2, and T6, the free rotation about the C-X bond can be envisaged.If the torsion angle between C-C-X-C is 0 • (marked with red color in Figure 2), then the TCNV group is in syn conformation; when the torsion angle is 180 • , the anti conformer exists.In the case of tautomers T3, T4, and T5, there is a formal double bond between C and X atoms; free rotation is impossible, and therefore cis and trans isomers are distinguished.In the case of tautomer T3, where the proton is attached to the terminal carbon atom, the carbon atom has sp 3 character, and the TCNV group is not planar.In there, the conformers can be distinguished as syn(H) (torsion angle between H-C-C-X atoms is 0 • ), anti(H) (180 • ), and eclipsed(H) (≈90 • ), as shown in Figure 2.This type of distinction between different diastereomers has also been used when X = B. Distinguishing between diastereomers and tautomers is only needed for computations to see which form is more or less stable.In reality, a mixture of tautomers and different conformers may exist.This fact is taken into account when solution pK a values are calculated.
In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.Figure 1.Schematic presentation of tautomers calculated in the current work.T denotes a tautomer which is referred to everywhere in the text.In case X = B, the tautomers are zwitterions rather than neutral compounds.
In the case of each tautomer, different conformers and isomers (diastereomers, Figure 2) were included in the calculations.The TCNV group is always planar in the case of anions and tautomers T1, T2, and T4-T5.For anion and tautomers T1, T2, and T6, the free rotation about the C-X bond can be envisaged.If the torsion angle between C-C-X-C is 0° (marked with red color in Figure 2), then the TCNV group is in syn conformation; when the torsion angle is 180°, the anti conformer exists.In the case of tautomers T3, T4, and T5, there is a formal double bond between C and X atoms; free rotation is impossible, and therefore cis and trans isomers are distinguished.In the case of tautomer T3, where the proton is attached to the terminal carbon atom, the carbon atom has sp 3  TCNV group is not planar.In there, the conformers can be distinguished as syn(H) (torsion angle between H-C-C-X atoms is 0°), anti(H) (180°), and eclipsed(H) (≈90°), as shown in Figure 2.This type of distinction between different diastereomers has also been used when X = B. Distinguishing between diastereomers and tautomers is only needed for computations to see which form is more or less stable.In reality, a mixture of tautomers and different conformers may exist.This fact is taken into account when solution pKa values are calculated.

Gas-Phase Acidities
The final results and the experimental or theoretical literature data are presented in Table 3.The whole set of calculations with the E, H, and G values and COSMO-RS pKa values (see below) are available in Supplementary Materials and the data repository [26].
According to the workflow, the DFT BP86/def-TZVP method was used to carry out gas-phase calculations (i.e., geometry optimization, and frequency calculations) for all the tautomers and diastereomers.Altogether, approximately 400 gas-phase calculations were carried out, and more than 300 GA values were obtained (Figure 3).GA values were always calculated for the structurally most similar acid-base pairs.This means, for example, for the syn-syn configuration acid, the syn-syn structured anion is used.Current GA values do not contain information from different tautomers or conformers; they have been calculated only using the most stable structurally similar forms.

anti syn trans cis syn(H) anti(H) eclipsed(H)
Anion, T1, T2, T6: T3, T4, T5: T3:  In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.In each case, different tautomers were included in the calculations (Figure 1).If compounds bearing only CN groups are considered, tautomers T1 and T2 exist.T2, T4, and T5 are imino taoutomers.Tautomers T1-T5 have resonance structures written in neutral form, and T6 can be written in zwitterionic form.In the case of central atom B, the tautomers do not correspond to the structural formula in Figure 1 due to the peculiarity of boron bonding; they correspond more to zwitterionic resonance structures rather than neutral forms.
Altogether, 14 title compounds were studied.Additionally, to compare the acidities, five related compounds were added to the study (Table 2).We look at each anion or acid from the perspective of the central atom O, N, C, and B, sometimes in the text also referred to as X.CN and TCNV groups are always looked at as substituents, even if protonated, without further specification of their names.

Gas-Phase Acidities
The final results and the experimental or theoretical literature data are presented in Table 3.The whole set of calculations with the E, H, and G values and COSMO-RS pK a values (see below) are available in Supplementary Materials and the data repository [26].
According to the workflow, the DFT BP86/def-TZVP method was used to carry out gas-phase calculations (i.e., geometry optimization, and frequency calculations) for all the tautomers and diastereomers.Altogether, approximately 400 gas-phase calculations were carried out, and more than 300 GA values were obtained (Figure 3).GA values were always calculated for the structurally most similar acid-base pairs.This means, for example, for the syn-syn configuration acid, the syn-syn structured anion is used.Current GA values do not contain information from different tautomers or conformers; they have been calculated only using the most stable structurally similar forms.
After the most stable forms were received, other calculation methods were employed to evaluate the suitability of the def-TZVP method (Table 4).These different methods were used only to calculate GA values for the most stable acid-base pairs resulting from the BP86/def-TZVP method.It was decided that GA values obtained from the B3LYP/6-311 + G(d,p) method were to be used as selected GA values (Table 3 and Figure 4) to make the results comparable with Maksic, Vianello, and Koppel's work (e.g., [9,27]).This calculation method has a sufficiently high precision to obtain data for larger molecules and is generally accepted as a precise method for H and G values [28].The G4(MP2) method [29], known as the method obtaining exact energies and which matches generally very well with the experimental data, was also tested among other accurate methods, such as W1RO and CBS-APNO; however, all these methods are highly resource-consuming and it was not possible to get results for all the compounds.The only exception is hydrogen borates, where the initial three replacements increase acidity a humble 3 to 5 kcal/mol, but the last replacement increase it almost to 17 kcal/mol.This phenomenon comes from the fact that the protonation site changes from a relatively stable T2 tautomer to T6.
Divalent oxygen substituted with the CN group exhibits two tautomers: OH-acid (cyanic acid, 1H-T1) and 29.8 kcal/mol more stable NH-acid (isocyanic acid 1H-T2).Replacement of the CN group with the TCNV group in isocyanic acid 1H leads to an enormous acidity increase by 52.6 kcal/mol, reaching the GA value of 282.1 kcal/mol for 2H-T1.Such an impressive acidity increase comes from the fact that the size of the anion as well as the charge delocalization possibility increases considerably.Small [CNO] − anion has a limited ability to delocalize the negative charge, and therefore, the acid form does not deprotonate easily.When the TCNV group replaces the CN group, the negative charge is stabilized by extensive delocalization and induction by CN groups.Figure 5 also illustrates that CN groups alone have a much higher negative charge concentration on the nitrogen (red color).In other cases, the GA value decreases from the initial to the final compound (i.e., fully substituted by TCNV groups) by around 35-30 kcal/mol.The more CN groups the initial acid contains, the lower the acidity increase is.This indicates that only a certain number of CN groups is necessary for a stable anion and strong acid.All the GA values from all the tested calculation methods stay within 11 kcal/mol (Table 4).Usually, they remain between 3 and 5 kcal/mol.The B3LYP/6-311 + G(d,p) calculations result in the highest acidities (i.e., the lowest GA values); however, in most cases, they match pretty well with G4(MP2) GA values.The most significant differences are for the molecules that are extensively substituted with TCNV groups, where perhaps the electron diffusion and its precision in the calculation may play a crucial role.BP86/def-TVZP values, in fact, describe all the title compounds sufficiently well.
In Figure 3, all the GA values calculated in this work are depicted as circles, colorcoded based on different tautomers.It gives an excellent overview of how, moving down to compounds with more TCNV groups, the GA values of different tautomers come closer to each other.In the case of the first five compounds, groups of tautomers lay relatively far from each other.In the case of C-centred molecules, the different GA values of different tautomers are closer or even mixed.For acids 6H, 7H, and 8H, tautomer T2 is the most stable form.For the acid 9H, the tautomers T1, T4, and T5 have very similar GA values, while the GA values for tautomer T1 vary widely, ranging from 242 to 251 kcal/mol, depending on the conformation.In all cases, the T6 tautomer is one of the least stable acids except for hydrogen borates: the T6 tautomer will become more and more stable compared to T4 and T5 tautomers, and when a single CN group is missing, the T6 tautomer is the most stable acid form.
In general, when the anion contains a single CN group (T2), then it is the most basic site of the anion.When the single CN group is missing, then in three cases out of four (central atom O, N, and C), the T1 tautomer is the most stable neutral form.In the case of hydrogen borates, T1 cannot be the most stable form; instead, tautomer T6 is the most stable acid form.
Figure 4 presents the overview of the change in GA values during the replacement of CN groups by TCNV groups of initial compounds (B3LYP/6-311 + G(d,p) basis set).It is seen how every additional replacement increases acidity less than the first replacement.The only exception is hydrogen borates, where the initial three replacements increase acidity a humble 3 to 5 kcal/mol, but the last replacement increase it almost to 17 kcal/mol.This phenomenon comes from the fact that the protonation site changes from a relatively stable T2 tautomer to T6.
Divalent oxygen substituted with the CN group exhibits two tautomers: OH-acid (cyanic acid, 1H-T1) and 29.8 kcal/mol more stable NH-acid (isocyanic acid 1H-T2).Replacement of the CN group with the TCNV group in isocyanic acid 1H leads to an enormous acidity increase by 52.6 kcal/mol, reaching the GA value of 282.1 kcal/mol for 2H-T1.Such an impressive acidity increase comes from the fact that the size of the anion as well as the charge delocalization possibility increases considerably.Small [CNO] − anion has a limited ability to delocalize the negative charge, and therefore, the acid form does not deprotonate easily.When the TCNV group replaces the CN group, the negative charge is stabilized by extensive delocalization and induction by CN groups.Figure 5 also illustrates that CN groups alone have a much higher negative charge concentration on the nitrogen (red color).In other cases, the GA value decreases from the initial to the final compound (i.e., fully substituted by TCNV groups) by around 35-30 kcal/mol.The more CN groups the initial acid contains, the lower the acidity increase is.This indicates that only a certain number of CN groups is necessary for a stable anion and strong acid.

Solution Phase Acidities
Computational solution phase pK a values in acetonitrile (MeCN) and dimethylsulfoxide (DMSO) were obtained using the COSMO-RS approach [37][38][39][40].COSMO-RS is generally unsuitable to get absolute pK a values; therefore, COSMO-RS pK a values must be corrected using experimental pK a values for similar acids [41].For this reason, 12 auxiliary acids, whose pK a values in MeCN are known, were used in the correlation to obtain corrected pK a (MeCN) values for the title compounds (Table 5 and Figure 6).Most of the known pK a (MeCN) values were, in fact, not directly measured in MeCN but were obtained from the correlation between pK a (DCE) and pK a (MeCN) values [3].It means that the literature pK a (MeCN) values are already obtained from the correlation.All 12 auxiliary acids were 2-substituted-1,1,3,3-tetracyanopropenes (2-X-TCNP); T4 tautomers of these acids were the most stable forms and only those tautomers were used to calculate COSMO-RS pK a values.In the case of 3,4-(MeO) 2 -C 6 H 4 -and CN-CH 2 -substituted TCNP-s, two different conformers were used for COSMO-RS pK a values.For the title compounds, different tautomers and conformers were used to calculate COSMO-RS pK a values.Table 5.A list of auxiliary acids used to obtain corrected pK a (MeCN) values, 2H and 10H, is excluded from the correlation.X-TCNP = 2-X-1,1,3,3-tetracyanopropene [3].In addition to 12 auxiliary acids, 5 title compounds, whose literature pK a (MeCN) values are also known, were also included in the correlation to correct COSMO-RS values.As seen, two of those compounds, [O-TCNV]H 2H and [(CN) 4 B]H 10H, deviate significantly from the correlation line (Figure 6).2H is an OH acid, which does not suit the correlation of NH-acids, and hydrogen borates, even though protonated on the CN group, can deviate because of the peculiarity of the C-N-H fragment (see below).pK a (MeCN) LITERATURE = 1.003 (0.052) • pK a (MeCN) COSMO-RS + 2.14 (0.11) r 2 = 0.966; S = 0.416; n = 15 (4)

Compound
Correlation 4 contains only CN acids, whose protonation site is in the CN group.It is not the most proper way to correct pK a (MeCN) values for OH (2H) and CH (9H) acids, and also, perhaps hydrogen borates (10H-14H) would need a different correlation; however, the scarcity of data does not allow creating an individual specific correlation equation for every group of compounds.Therefore, Equation ( 4) was used to obtain corrected computational pK a (MeCN) values for all the studied compounds.Most of the values are consistent with the literature data.It is assumed that even if the absolute corrected pK a (MeCN) values deviate from the actual pK a values, then relative pK a values for hydrogen borates, as well as cyclic carboranes and boranes, should be consistent with each other.It is also possible that the literature pK a (MeCN) value of 10H is deviating since it is already in the literature obtained from the correlation.pK a (DMSO) values in the current work are obtained from pK a (MeCN) values by averaging two different values calculated from previously published pK a (MeCN) vs pK a (DMSO) correlation equations [34,42].COSMO-RS pK a (DMSO) values are also available; the correlation between acquired pK a (DMSO) values and COSMO-RS pK a (DMSO) values is very good (see SM), and it could be used to correct COSMO-RS pK a values.pK a values for 1,2-dichloroethane (DCE) were also obtained from the pK a (MeCN) values using a previously published correlation [3,36].
It is seen from Table 3 that our method strongly downgrades the acidity of HCN in MeCN (5 units) and DMSO (4 units) compared to the literature values.The literature pK a value in MeCN is obtained computationally, and the pK a value in DMSO is directly measured.The pK a value of HCN depends very strongly on the solvation because its conjugate base is a very small anion.Computational methods that do not take solvation specifically into account, as the method used in our work, can give considerably wrong values.The pK a values of HCN should be under deeper investigation in the future.

Structural Details
In Figure 7, the structures of the initial anions and their conjugate acids are presented.Anion 1 is stick-like, 3 is bent, 6 is planar, and 10 is a tetrahedron.The most stable tautomer in the case of all initial acids is T2, i.e., the proton attached to the nitrogen on the CN group.H-N-C is bent except for the hydrogen borates, where the angle is almost always close to 180 • .
The geometry of other anions corresponds to their parent initial anion.TCNV groups are always planar, and the angles between substituents are close to the C-X-C angle shown in Figure 7. Angles between TCNV groups are generally larger than the angles between CN groups.
It is interesting to see how the most stable forms of anions that contain two TCNV groups differ from each other (Figure 8).The more single CN groups are bound directly to the central atom, the more the TCNV groups are compressed.The anti-anti conformer of anion 5 is 4.2 kcal/mol (BP86/def-TZVP results) more stable than the syn-syn conformer.Anti-syn 8 is 0.6 and 0.7 kcal/mol more stable than syn-syn and anti-anti, respectively.Synsyn 12 is only 0.1 kcal/mol more stable than anti-syn, but anti-anti conformer already has 0.8 kcal/mol higher energy.This shows that the fewer substituents the anion contains, the more rigid the structure is.The most stable conformers of the corresponding conjugate acids of the anions 5 and 8 are the same.For 12H-T2, the most stable conformer is anti-syn, and other conformers are 1.1 kcal/mol less stable.

Structural Details
In Figure 7, the structures of the initial anions and their conjugate acids are presented.Anion 1 is stick-like, 3 is bent, 6 is planar, and 10 is a tetrahedron.The most stable tautomer in the case of all initial acids is T2, i.e., the proton attached to the nitrogen on the CN group.H-N-C is bent except for the hydrogen borates, where the angle is almost always close to 180°.The geometry of other anions corresponds to their parent initial anion.TCNV groups are always planar, and the angles between substituents are close to the C-X-C angle shown in Figure 7. Angles between TCNV groups are generally larger than the angles between CN groups.
It is interesting to see how the most stable forms of anions that contain two TCNV groups differ from each other (Figure 8).The more single CN groups are bound directly to the central atom, the more the TCNV groups are compressed.The anti-anti conformer of anion 5 is 4.2 kcal/mol (BP86/def-TZVP results) more stable than the syn-syn conformer.Anti-syn 8 is 0.6 and 0.7 kcal/mol more stable than syn-syn and anti-anti, respectively.Synsyn 12 is only 0.1 kcal/mol more stable than anti-syn, but anti-anti conformer already has 0.8 kcal/mol higher energy.This shows that the fewer substituents the anion contains, the more rigid the structure is.The most stable conformers of the corresponding conjugate acids of the anions 5 and 8 are the same.For 12H-T2, the most stable conformer is antisyn, and other conformers are 1.1 kcal/mol less stable.

1H-T2 3H-T2 6H-T2 10H-T2 < H-N-C [°]
126.8 126.6 123.9 179.9 DG [kcal/mol] 29.8 10.1 0.9  Further analysis of the structures and energetics for each compound gro presented in the Supplementary Material.Further analysis of the structures and energetics for each compound gro presented in the Supplementary Material.

Figure 1 .
Figure 1.Schematic presentation of tautomers calculated in the current work.T denotes a tautomer which is referred to everywhere in the text.In case X = B, the tautomers are zwitterions rather than neutral compounds.

Figure 3 .
Figure 3.Each circle indicates calculated GA values with all possible color-coded calculation methods based on different tautomers.

Figure 4
Figure 4 presents the overview of the change in GA values during the replacement of CN groups by TCNV groups of initial compounds (B3LYP/6-311 + G(d,p) basis set).It is seen how every additional replacement increases acidity less than the first replacement.Theonly exception is hydrogen borates, where the initial three replacements increase acidity a humble 3 to 5 kcal/mol, but the last replacement increase it almost to 17 kcal/mol.This phenomenon comes from the fact that the protonation site changes from a relatively stable T2 tautomer to T6.Divalent oxygen substituted with the CN group exhibits two tautomers: OH-acid (cyanic acid, 1H-T1) and 29.8 kcal/mol more stable NH-acid (isocyanic acid 1H-T2).Replacement of the CN group with the TCNV group in isocyanic acid 1H leads to an enormous acidity increase by 52.6 kcal/mol, reaching the GA value of 282.1 kcal/mol for 2H-T1.Such an impressive acidity increase comes from the fact that the size of the anion as well as the charge delocalization possibility increases considerably.Small [CNO] − anion has a limited ability to delocalize the negative charge, and therefore, the acid form does not deprotonate easily.When the TCNV group replaces the CN group, the negative charge is stabilized by extensive delocalization and induction by CN groups.Figure5also illustrates that CN groups alone have a much higher negative charge concentration on the nitrogen (red color).In other cases, the GA value decreases from the initial to the final compound (i.e., fully substituted by TCNV groups) by around 35-30 kcal/mol.The more CN groups the initial acid contains, the lower the acidity increase is.This indicates that only a certain number of CN groups is necessary for a stable anion and strong acid.

Figure 3 .
Figure 3.Each circle indicates calculated GA values with all possible color-coded calculation methods based on different tautomers.Molecules 2023, 28, x FOR PEER REVIEW 9 of 16

Figure 4 .
Figure 4. Overview of the change in GA values (B3LYP/6-311 + G(d,p)) in stepwise replacement of CN groups by TCNV groups of title compounds; additional acids are included.Protonation sites are indicated with *.

Figure 4 .
Figure 4. Overview of the change in GA values (B3LYP/6-311 + G(d,p)) in stepwise replacement of CN groups by TCNV groups of title compounds; additional acids are included.Protonation sites are indicated with *.

Figure 4 .
Figure 4. Overview of the change in GA values (B3LYP/6-311 + G(d,p)) in stepwise replacement of CN groups by TCNV groups of title compounds; additional acids are included.Protonation sites are indicated with *.

Figure 5 .
Figure 5.The surface of the charge density of anions 3, 4, and 5, computed using the COSMO-RS approach.The red area indicates a negative partial charge on the CN groups and the green area indicates close to zero partial charge.As a comparison, the GA values of hydrogen cyanide (HCN) 15H and its altered analog [TCNV]H 16H are included.HCN is a weak acid in the gas phase, having practically the same acidity as isocyanic acid [CNO]H 1H.The most stable T1 tautomer of [TCNV]H is 18.7 kcal/mol stronger acid than HCN.The most stable tautomers are both CH-acids; in the case of 16H, the negative charge can delocalize into only one CN group.The anionic acid B12(CN)12H -17H is as strong as cyanoform.The B12(CN)12H2 18H acid is almost as strong as [B-TCNV4]H, 14H.Carborane acid, 19H, with 12 CN groups, is 4.4 kcal/mol stronger acid than 14H.Tricyanovinyl alcohol [O-TCNV]H 2H and cyanoform [(CN)3C]H 6H are considered strong acids, and they have very similar acidities in the gas phase and the solution.Three

5 Figure 5 .
Figure 5.The surface of the charge density of anions 3, 4, and 5, computed using the COSMO-RS approach.The red area indicates a negative partial charge on the CN groups and the green area indicates close to zero partial charge.As a comparison, the GA values of hydrogen cyanide (HCN) 15H and its altered analog [TCNV]H 16H are included.HCN is a weak acid in the gas phase, having practically the same acidity as isocyanic acid [CNO]H 1H.The most stable T1 tautomer of [TCNV]H is 18.7 kcal/mol stronger acid than HCN.The most stable tautomers are both CH-acids; in the case of 16H, the negative charge can delocalize into only one CN group.The anionic acid B 12 (CN) 12 H -17H is as strong as cyanoform.The B 12 (CN) 12 H 2 18H acid is almost as strong as [B-TCNV 4 ]H, 14H.Carborane acid, 19H, with 12 CN groups, is 4.4 kcal/mol stronger acid than 14H.Tricyanovinyl alcohol [O-TCNV]H 2H and cyanoform [(CN) 3 C]H 6H are considered strong acids, and they have very similar acidities in the gas phase and the solution.Three CN groups seem to work as efficiently even though the protonation site (O vs. N) and the conformation are different.

Figure 6 .
Figure 6.Blue dots are used to create the correlation between literature and COSMO-RS computational values (Equation (4)) of the compounds listed in Table 5, whereas red points (2H and 10H) have been omitted from the correlation.

Figure 7 .
Figure 7. Geometries of initial anions 1, 3, 6, and 10, and their most stable acid forms.H-N-C angles are included.ΔG values (kcal/mol) show how much the T2 tautomer is more stable relative to the T1 tautomer at BP86/def-TZVP level.Here and in other figures (also in Supplementary Materials) red color corresponds to oxygen, green to carbon, blue to nitrogen and pink to boron atoms.

Figure 7 .Figure 8 .
Figure 7. Geometries of initial anions 1, 3, 6, and 10, and their most stable acid forms.H-N-C angles are included.∆G values (kcal/mol) show how much the T2 tautomer is more stable relative to the T1 tautomer at BP86/def-TZVP level.Here and in other figures (also in Supplementary Materials) red color corresponds to oxygen, green to carbon, blue to nitrogen and pink to boron atoms.Molecules 2023, 28, x FOR PEER REVIEW 1

Figure 9 .
Figure 9. Geometry of anions 9 and 14.The surface of the charge density of anion 14.

Figure 8 .Figure 8 .
Figure 8.The most stable conformers of anions 5, 8, and 12 containing two TCNV groups.In Figure 9, the geometries of the bulkiest anions are shown.The C-C-C angles of the symmetrical anion 9 are all 119.9 • , and the central part containing four carbon atoms is planar; only TCNV groups are bent out of the plane at 39.6 • .Tetrahedral anion 14 has C-B-C angles all between 109.2 and 109.7 • .The surface of the charge density shows that the negative charge concentration is highest on the N atoms of CN groups, but the charge is very effectively delocalized.

Figure 9 .
Figure 9. Geometry of anions 9 and 14.The surface of the charge density of anion 14.

14 Figure 9 .
Figure 9. Geometry of anions 9 and 14.The surface of the charge density of anion 14.

Table 1 .
Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Table 1 .
Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Table 1 .
Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Table 1 .
Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Table 1 .
Schematic presentation of the anions of title compounds studied in the current work.In the text, anions' conjugate acids are denoted as NumberH with tautomer number (T1-T6), if necessary.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 2 .
Other compounds as anions which are used for comparison.

Table 3 .
Theoretical (DFT B3LYP/6-311 + G(d,p)) PA and GA values (kcal/mol).All literature values are shown in brackets, whereas values in parentheses have a computational origin.Corrected COSMO-RS pK a (MeCN) values, pK a (DMSO), and pK a (DCE) values are derived from correlation analysis.

Table 4 .
GA values obtained using different calculation methods.Assigned GA values are in bold.