The electron affinities of the silicon fluorides SiFn (n=1–5)

Probing the Effects of Size and Charge on the Monohydration and Dihydration of SiF5- and SiF62- via Comparisons with BF4- and PF6. J Phys Chem A 2024;128:5637-5645. [PMID: 38976798 DOI: 10.1021/acs.jpca.4c03430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

This study systematically examines the interactions of the trigonal bipyramidal silicon pentafluoride and octahedral silicon hexafluoride anions with either one or two water molecules, (SiF5-(H2O)n and SiF62-(H2O)n, respectively, where n = 1, 2). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using the CCSD(T) ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for Si, O, and F; denoted as haTZ). Two monohydrate and six dihydrate minima have been identified for the SiF5-(H2O)n systems, whereas one monohydrate and five dihydrate minima have been identified for the SiF62-(H2O)n systems. Both monohydrated anions have a minimum in which the water molecule adopts a symmetric double ionic hydrogen bond (DIHB) motif with C2v symmetry. However, a second unique monohydrate minimum has been identified for SiF5- in which the water molecule adopts an asymmetric DIHB motif along the edge of the trigonal bipyramidal anion between one axial and one equatorial F atom. This Cs structure is more than 2 kcal mol-1 lower in energy than the C2v local minimum at the CCSD(T)/haTZ level of theory. While the interactions between the solvent and ionic solute are quite strong for the monohydrated anions (electronic dissociation energies of ≈12 and ≈24 kcal mol-1 for the SiF5-(H2O)1 and SiF62-(H2O)1 global minima, respectively), these values are nearly perfectly doubled for the dihydrates, with the lowest-energy SiF5-(H2O)2 and SiF62-(H2O)2 minima exhibiting dissociation energies of ≈24 and ≈47 kcal mol-1, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -100 cm-1 for the SiF5-(H2O)2 minimum and -300 cm-1 for the SiF62-(H2O)2 minimum relative to an isolated H2O molecule at the CCSD(T)/haTZ level of theory. This work also corrects the OH stretching frequency shifts for two dihydrate minima of PF6- that were previously erroneously reported ( J. Phys. Chem. A 2020, 124, 8744-8752, DOI: 10.1021/acs.jpca.0c06466).

Strongly Bound Polynuclear Anions Comprising Scandium Fluoride Building Blocks

The stability of polynuclear anions composed of ScF3 building blocks was studied by using ab initio and density functional theory electronic structure methods and flexible basis sets. Thorough exploration of ground state potential energy surfaces of (Sc2F7)-, (Sc3F10)-, and (Sc4F13)- anions which may be viewed as comprising ScF3 fragments and the additional fluorine atom led to determining the isomeric structures thereof. It was found that the most stable isomers which are predicted to dominate at room temperature correspond to the compact structures enabling the formation of a large number of Sc-F-Sc bridging linkages rather than to the chain-like structures. The vertical electron detachment energies of the (ScnF3n+1)- anions were found to be very large (spanning the 10.85-12.29 eV range) and increasing with the increasing number of scandium atoms (n) and thus the ScF3 building blocks involved in the structure. Thermodynamic stability of (ScnF3n+1)- anions (i.e., their susceptibility to fragmentation) was also verified and discussed.

Reproducibility of density functional approximations: How new functionals should be reported

Density functional theory is the workhorse of chemistry and materials science, and novel density functional approximations are published every year. To become available in program packages, the novel density functional approximations (DFAs) need to be (re)implemented. However, according to our experience as developers of Libxc [Lehtola et al., SoftwareX 7, 1 (2018)], a constant problem in this task is verification due to the lack of reliable reference data. As we discuss in this work, this lack has led to several non-equivalent implementations of functionals such as Becke-Perdew 1986, Perdew-Wang 1991, Perdew-Burke-Ernzerhof, and Becke's three-parameter hybrid functional with Lee-Yang-Parr correlation across various program packages, yielding different total energies. Through careful verification, we have also found many issues with incorrect functional forms in recent DFAs. The goal of this work is to ensure the reproducibility of DFAs. DFAs must be verifiable in order to prevent the reappearance of the above-mentioned errors and incompatibilities. A common framework for verification and testing is, therefore, needed. We suggest several ways in which reference energies can be produced with free and open source software, either with non-self-consistent calculations with tabulated atomic densities or via self-consistent calculations with various program packages. The employed numerical parameters-especially the quadrature grid-need to be converged to guarantee a ≲0.1 μEh precision in the total energy, which is nowadays routinely achievable in fully numerical calculations. Moreover, as such sub-μEh level agreement can only be achieved when fully equivalent implementations of the DFA are used, the source code of the reference implementation should also be made available in any publication describing a new DFA.

Benchmarking Magnetizabilities with Recent Density Functionals

We have assessed the accuracy of the magnetic properties of a set of 51 density functional approximations, including both recently published and already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled-cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, are obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also, some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP, and PBE0 perform rather well. In contrast, unsatisfactory performance is generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current (GIMIC) method. Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of magnetizabilities.

Oxidizing CO2 with superhalogens

The Sb₃F₁₆ species was found to be capable of ionizing the CO₂ molecule.

Stability of Superhalogen Anions in the Aqueous Phase

The issue of stability of superhalogen anions in an aqueous solution is investigated on the basis of theoretical calculations carried out at the CCSD(T)/6-311++G(d,p)//MP2/6-311++G(d,p) level for two representative negatively charged systems (NaF₂^- and AlF₄^-) whose fragmentation products differ in polarity. The presence of a water solvent is simulated independently by employing the polarized continuum solvation model and by involving eight H₂O molecules explicitly to allow interactions at the molecular level. The best estimates of the Gibbs free energies characterizing the AlF₄^- and NaF₂^- fragmentation reactions in a water solvent are evaluated as equal to 33-34 and 12-14 kcal/mol, respectively (assuming the F^- and AlF₃/NaF products) or 14-15 and 26-28 kcal/mol, respectively (assuming the HF and AlF₃OH^-/NaFOH^- products). The corresponding fragmentation routes are suggested to be nonoperative at T = 298.15 K. The conclusion concerning the thermodynamic stability of the AlF₄^- and NaF₂^- superhalogen anions in the aqueous phase is formulated and discussed.

Density functional theory and Raman spectroscopy applied to structure and vibrational mode analysis of 1,1',3,3'-tetraethyl-5,5',6,6'-tetrachloro- benzimidazolocarbocyanine iodide and its aggregate

We have measured electronic and Raman scattering spectra of 1,1',3,3'-tetraethyl-5,5',6,6'-tetrachloro-benzimidazolocarbocyanine iodide (TTBC) in various environments, and we have calculated the ground state geometric and spectroscopic properties of the TTBC cation in the gas and solution phases (e.g., bond distances, bond angles, charge distributions, and Raman vibrational frequencies) using density functional theory. Our structure calculations have shown that the ground state equilibrium structure of a cis-conformer lies ∼200 cm(-1) above that of a trans-conformer and both conformers have C(2) symmetry. Calculated electronic transitions indicate that the difference between the first transitions of the two conformers is about 130 cm(-1). Raman spectral assignments of monomeric- and aggregated-TTBC cations have been aided by density functional calculations at the same level of the theory. Vibrational mode analyses of the calculated Raman spectra reveal that the observed Raman bands above 700 cm(-1) are mainly associated with the in-plane deformation of the benzimidazolo moieties, while bands below 700 cm(-1) are associated with out-of-plane deformations of the benzimidazolo moieties. We have also found that for the nonresonance excited experimental Raman spectrum of aggregated-TTBC cation, the Raman bands in the higher-frequency region are enhanced compared with those in the nonresonance spectrum of the monomeric cation. For the experimental Raman spectrum of the aggregate under resonance excitation, however, we find new Raman features below 600 cm(-1), in addition to a significantly enhanced Raman peak at 671 cm(-1) that are associated with out-of-plane distortions. Also, time-dependent density functional theory calculations suggest that the experimentally observed electronic transition at ∼515 nm (i.e., 2.41 eV) in the absorption spectrum of the monomeric-TTBC cation predominantly results from the π → π∗ transition. Calculations are further interpreted as indicating that the observed shoulder in the absorption spectrum of TTBC in methanol at 494 nm (i.e., 2.51 eV) likely results from the ν(") = 0 → ν' = 1 transition and is not due to another electronic transition of the trans-conformer-despite the fact that measured and calculated NMR results (not provided here) support the prospect that the shoulder might be attributable to the 0-0 band of the cis-conformer.

Electron binding capabilities of some silylenes having small singlet-triplet splittings or triplet ground states

Several silyl and alkaline metal substituted silylenes have been investigated using the CAS-ACPF method in conjunction with the aug-cc-pVTZ basis sets. Silylsilylene and disilylsilylene are found to have singlet ground states with DeltaEST(-) values of 0.676 and 0.319 eV, respectively. The adiabatic ground state electron affinities are found to be 1.572 and 2.361 eV for HSiSiH(3) and Si(SiH(3))(2). respectively. Both silylenes possesses a stable 2A1 excited negative ion state, with respective adiabatic EA values of 0.037 and 1.000 eV. In contrast, all silylenes with at least one akaline metal substituent exhibit triplet neutral ground states. The metalated silylenes HSiLi, HSiNa, LiSiSiH(3), NaSiLi, SiLi(2), and SiNa(2) have adiabatic ground state EAs somewhat below 1 eV, but each of these negatively charged system possesses up to three bound excited negative ion states, some of which are dipole-bound states.

Global minimum structure searches via particle swarm optimization

Novel implementation of the evolutionary approach known as particle swarm optimization (PSO) capable of finding the global minimum of the potential energy surface of atomic assemblies is reported. This is the first time the PSO technique has been used to perform global optimization of minimum structure search for chemical systems. Significant improvements have been introduced to the original PSO algorithm to increase its efficiency and reliability and adapt it to chemical systems. The developed software has successfully found the lowest-energy structures of the LJ(26) Lennard-Jones cluster, anionic silicon hydride Si(2)H(5) (-), and triply hydrated hydroxide ion OH(-) (H(2)O)(3). It requires relatively small population sizes and demonstrates fast convergence. Efficiency of PSO has been compared with simulated annealing, and the gradient embedded genetic algorithm.

Induced hyper-valence of carbon in metal:fluorocarbon complexes

Charge-transfer in combination with geometric features induces a pentavalent state of the carbon atom in metastable super-fluorinated alkali-metal compounds, M:C2F7 (M = K to Cs).

Metastable Intermolecular Charge-Transfer Complexes with a Pentavalent Carbon Atom

A novel type of metastable complexes of metal (such as Na) atoms with a super-fluorinated carbon molecule is investigated, with a carbon atom exhibiting a unique, pentavalent character. It is induced by the charge transfer from the alkali metal component, followed by a geometric compression of the ion pair system. Analysis of the electron-density distribution confirms the real chemical nature of the extra C-F bond. Structure and stability of the system are characterized ab initio, and a spectrum of electronic perturbations is considered. The ways of forming the systems are discussed, and the spectroscopic parameters are predicted, facilitating their detection in experiments.

Singlet-triplet splittings and ground- and excited-state electron affinities of selected cyanosilylenes, XSiCN (X = H, F, Cl, CH3, SiH3, CN)

Several cyanosilylenes, XSiCN, (X = H, F, Cl, CH3, SiH3, CN) have been investigated using the RHF-ACPF and CAS(2,2)-ACPF methods in conjunction with the aug-cc-pVTZ basis sets. All silylenes are found to have singlet ground states. The ground-state electron affinities are found to be rather high, i.e., 1.832, 1.497, 1.896, 1.492, 2.235, and 2.631 eV for HSiCN, FSiCN, ClSiCN, H3CSiCN, H3SiSiCN, and Si(CN)2, respectively. The existence of bound excited negative ion states has been discovered for the first time within these silylenes. All these bound excited anion states belong to the totally symmetric irreducible representations and can be characterized as dipole-bound negative ion states. All triplet excited states have even larger dipole moments than the singlet states and are, therefore, "dressed" by dipole-bound negative ion states, which correspond to Feshbach resonances.

Calculation of Electron Detachment Energies for Water Cluster Anions: An Appraisal of Electronic Structure Methods, with Application to (H2O)20- and (H2O)24-

We present benchmark calculations of vertical electron detachment energies (VDEs) for various conformers of (H2O)n-, using both wave function and density functional methods, in sequences of increasingly diffuse Gaussian basis sets. For small clusters (n < or = 6), a systematic examination of VDE convergence reveals that it is possible to converge this quantity to within approximately 0.01 eV of the complete-basis limit, using a highly diffuse but otherwise economical Pople-style basis set of double-zeta quality, with 28 atom-centered basis functions per water molecule. Floating-center basis functions can be useful but are not required to obtain accurate VDEs. Second-order Møller-Plesset perturbation (MP2) theory suffices to obtain VDEs that are within 0.05 eV of the results from both experiment and coupled-cluster theory, and which always err toward underbinding the extra electron. In contrast to these consistent predictions, VDEs calculated using density functional theory (DFT) vary widely, according to the fraction of Hartree-Fock exchange in a given functional. Common functionals such as BLYP and B3LYP overestimate the VDE by 0.2-0.5 eV, whereas a variant of Becke's "half and half" functional is much closer to coupled-cluster predictions. Exploratory calculations for (H2O)20- and (H2O)24- cast considerable doubt on earlier calculations that were used to assign the photoelectron spectra of these species to particular cluster isomers.

Structures and electron affinities of the di-arsenic fluorides As2Fn/As2F n− (n= 1-8)

Developments in the preparation of new materials for microelectronics are focusing new attention on molecular systems incorporating several arsenic atoms. A systematic investigation of the As2Fn/As2Fn- systems was carried out using Density Functional Theory methods and a DZP++ quality basis set. Global and low-lying local geometric minima and relative energies are discussed and compared. The three types of neutral-anion separations reported in this work are: the adiabatic electron affinity (EAad), the vertical electron affinity (EAvert), and the vertical detachment energy (VDE). Harmonic vibrational frequencies pertaining to the global minimum for each compound are reported. From the first four studied species (As2Fn, n=1-4), all neutral molecules and their anions are shown to be stable with respect to As-As bond breaking. The neutral As2F molecule and its anion are predicted to have Cs symmetry. We find the trans F-As-As-F isomer of C2h symmetry and a pyramidalized vinylidene-like As-As-F2- isomer of Cs symmetry to be the global minima for the As2F2 and As2F2- species, respectively. The lowest lying minima of As2F3 and As2F3- are vinyl radical-like structures F-As-As-F2 of Cs symmetry. The neutral As2F4 global minimum is a trans-bent (like Si2H4) F2-As-As-F2 isomer of C2 symmetry, while its anion is predicted to have an unusual fluorine-bridged (C(1)) structure. The global minima of the neutral As2Fn species, n=5-8, are weakly bound complexes, held together by dipole-dipole interactions. All such structures have the AsFm-AsFn form, where (m,n) is (2,3) for As2F5, (3,3) for As2F6, (4,3) for As2F7), and (5,3) for As2F8. For As2F8 the beautiful pentavalent F4As-AsF4 structure (analogous to the stable AsF5 molecule) lies about 30 kcal/mol above the AsF3 . . . AsF5 complex. The stability of AsF(5) depends crucially on the strong As-F bonds, and replacing one of these with an As-As bond (in F4As-AsF4) has a very negative impact on the molecule's stability. The anions As2Fn-, n=5-8, are shown to be stable with respect to the As-As bond breaking, and we predict that all of them have fluorine-bridged or fluorine-linked structures. The zero-point vibrational energy corrected adiabatic electron affinities are predicted to be 2.28 eV (As2F), 1.95 eV (As2F2), 2.39 eV (As2F3), 1.71 eV (As2F4), 2.72 eV (As2F5), 1.79 eV (As2F6), 5.26 eV (As2F7), and 3.40 eV (As2F8) from the BHLYP method. Vertical detachment energies are rather large, especially for species with fluorine-bridged global minima, having values up to 6.45 eV (As2F7, BHLYP).

The ability of silylenes to bind excess electrons: electron affinities of SiX(2), and SiXY species (X,Y=H,CH(3),SiH(3),F,Cl,Br)

Recently, Ishida and co-workers have isolated silylene radical anions via the one-electron reduction of isolable cyclic dialkylsilylenes, discovering these corresponding radical anions to be relatively stable at low temperatures. Herein we report theoretical predictions of the adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies of a series of methyl, silyl, and halosubstituted silylene compounds. This research utilizes the carefully calibrated DZP++ basis with the combination of the popular nonhybrid and hybrid DFT functionals, BLYP, B3LYP, and BHHLYP. The level of theory employed and the ensemble of species under study confirm the ability of silylenes to bind excess electrons with Si(SiH(3))(2) being the most effective, having a predicted AEA of 1.95 eV. While it is known that methyl substituents have a diminishing effect on the computed electron affinities (EAs), it is shown that fluorine shows an analogous negative effect. Similarly, previous suggestions that Si(CH(3))(2) will not bind an electron appear incorrect, with EA[Si(CH(3))(2)] predicted here to be 0.46 eV.

Cyclic perfluorocarbon radicals and anions having high global warming potentials (GWPs): structures, electron affinities, and vibrational frequencies

Adiabatic electron affinities, optimized molecular geometries, and IR-active vibrational frequencies have been predicted for small cyclic hydrocarbon radicals C(n)H(2)(n)(-)(1) (n = 3-6) and their perfluoro counterparts C(n)F(2)(n)(-)(1) (n = 3-6). Total energies and optimized geometries of the radicals and corresponding anions have been obtained using carefully calibrated (Chem. Rev. 2002, 102, 231) density functional methods, namely, the B3LYP, BLYP, and BP86 functionals in conjunction with the DZP++ basis set. The predicted electron affinities show that only the cyclopropyl radical tends to bind electrons among the hydrocarbon radicals studied. The trend for the perfluorocarbon (PFC) radicals is quite different. The electron affinities increase with expanding ring size until n = 5 and then slightly decrease at n = 6. Predicted electron affinities of the hydrocarbon radicals using the B3LYP hybrid functional are 0.24 eV (C(3)H(5)/C(3)H(5)(-)), -0.19 eV (C(4)H(7)/C(4)H(7)(-)), -0.15 eV (C(5)H(9)/C(5)H(9)(-)), and -0.11 eV (C(6)H(11)/C(6)H(11)(-)). Analogous electron affinities of the perflurocarbon radicals are 2.81 eV (C(3)F(5)/C(3)F(5)(-)), 3.18 eV (C(4)F(7)/C(4)F(7)(-)), 3.34 eV (C(5)F(9)/C(5)F(9)(-)), and 3.21 eV (C(6)F(11)/C(6)F(11)(-)).

Odd carbon long linear chains HC2n+1H (n = 4-11): properties of the neutrals and radical anions

The optimized geometries, adiabatic electron affinities, vertical electron affinities, vertical electron detachment energies (for the anions), and IR-active vibrational frequencies have been predicted for the long linear carbon chains HC(2)(n)()(+1)H (n = 4-11). The B3LYP density functional combined with DZP and TZ2P basis sets was used in this theoretical study. These methods have been extensively calibrated versus experiment for the prediction of electron affinities (Chem. Rev. 2002, 102, 231). The computed physical properties are discussed and compared with the even carbon chains HC(2)(n)()H. The predicted electron affinities form a remarkably regular sequence: 2.12 eV (HC(9)H), 2.42 eV (HC(11)H), 2.66 eV (HC(13)H), 2.85 eV (HC(15)H), 3.01 eV (HC(17)H), 3.14 eV (HC(19)H), 3.25 eV (HC(21)H), and 3.35 eV (HC(23)H). These electron affinities are as much as 0.4 eV higher than those for analogous even carbon chains. The predicted structures display an intermediate cumulene-polyacetylene type of bonding, with the inner carbons appearing cumulenic and the outer carbons polyacetylenic.

Molecules for materials: germanium hydride neutrals and anions. Molecular structures, electron affinities, and thermochemistry of GeHn/GeHn- (n = 0-4) and Ge2Hn/Ge2Hn(-) (n = 0-6)

The GeH(n) (n = 0-4) and Ge(2)H(n) (n = 0-6) systems have been studied systematically by five different density functional methods. The basis sets employed are of double-zeta plus polarization quality with additional s- and p-type diffuse functions, labeled DZP++. For each compound plausible energetically low-lying structures were optimized. The methods used have been calibrated against a comprehensive tabulation of experimental electron affinities (Chemical Reviews 102, 231, 2002). The geometries predicted in this work include yet unknown anionic species, such as Ge(2)H(-), Ge(2)H(2)(-), Ge(2)H(3)(-), Ge(2)H(4)(-), and Ge(2)H(5)(-). In general, the BHLYP method predicts the geometries closest to the few available experimental structures. A number of structures rather different from the analogous well-characterized hydrocarbon radicals and anions are predicted. For example, a vinylidene-like GeGeH(2) (-) structure is the global minimum of Ge(2)H(2) (-). For neutral Ge(2)H(4), a methylcarbene-like HGë-GeH(3) is neally degenerate with the trans-bent H(2)Ge=GeH(2) structure. For the Ge(2)H(4) (-) anion, the methylcarbene-like system is the global minimum. The three different neutral-anion energy differences reported in this research are: the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). For this family of molecules the B3LYP method appears to predict the most reliable electron affinities. The adiabatic electron affinities after the ZPVE correction are predicted to be 2.02 (Ge(2)), 2.05 (Ge(2)H), 1.25 (Ge(2)H(2)), 2.09 (Ge(2)H(3)), 1.71 (Ge(2)H(4)), 2.17 (Ge(2)H(5)), and -0.02 (Ge(2)H(6)) eV. We also reported the dissociation energies for the GeH(n) (n = 1-4) and Ge(2)H(n) (n = 1-6) systems, as well as those for their anionic counterparts. Our theoretical predictions provide strong motivation for the further experimental study of these important germanium hydrides.

What is the nature of polyacetylene neutral and anionic chains HC(2n)H and HC(2n)H(-) (n = 6-12) that have recently been observed?

The optimized geometries, adiabatic electron affinities, and IR-active vibrational frequencies have been predicted for the long linear carbon chains HC(2n)H. The B3LYP density functional combined with the DZP basis set was used in this theoretical study. The computed physical properties are discussed. The predicted electron affinities form a remarkably regular sequence: 1.78 (HC(12)H), 2.08 (HC(14)H), 2.32 (HC(16)H), 2.53 (HC(18)H), 2.69 (HC(20)H), 2.83 (HC(22)H), and 2.95 eV (HC(24)H). The predicted structures display an alternating triple and very short single bond pattern, with the degree of bond alternation significantly less for the radical anions.

Molecules for materials: structures, thermochemistry, and electron affinities of the digermanium fluorides Ge2Fn/Ge2Fn- (n = 1-6): a wealth of unusual structures

A systematic investigation of Ge2Fn/Ge2Fn- systems was carried out with five density functional theory (DFT) methods in conjunction with DZP++ basis sets. For each compound various structures, including minima, transition states and other energetically low lying stationary points, were optimized. The geometries and relative energies are discussed and compared. Adiabatic electron affinities, vertical electron affinities and anion vertical detachment energies are reported. Three types of dissociation energies pertaining to the global minima for each compound are reported. The theoretical predictions are in good agreement with the limited available experimental results. Many unusual structural features are predicted for these systems. Neutral Ge2F is predicted to have a bridged C2v structure, while its anion is very floppy, with the bridged structure very slightly favoured. The Ge2F2 molecule is predicted to have the butterfly structure known from experiment for Si2H2, while the Ge2F3- ion has a trans-bent structure. Ge2F3 is predicted to have an unprecedented FGe-F-GeF structure with no Ge-Ge bond, while its anion has a somewhat more conventional monobridged structure, analogous to that of the nonclassical vinyl cation. Neutral Ge2F4 has a dibridged structure of C2h symmetry, while its anion has a trans-bent structure with a very long Ge-Ge bond. The Ge2F5 molecule is doubly bridged and has no Ge-Ge bond, while the anion is of the type F2Ge-F-GeF2, again with no Ge-Ge bond. Ge2F6 has the anticipated ethane structure, as does its anion, but with a very long Ge-Ge bond. The adiabatic electron affinities (EAad) are predicted to be 2.12 (Ge2F), 2.03 (Ge2F2), 2.02 (Ge2F3), 1.64 (Ge2F4), 4.57 (Ge2F5), and 2.66 eV (Ge2F6), respectively, by the BHLYP method, which is regarded as the best method in the present paper for predicting EAs. Comparisons with the analogous C2Fn and Si2Fn systems reveals some interesting trends and differences. For example, while C2F6 will not capture an electron, Si2F6 is predicted to have small electron affinity (0.73 eV), while that of Ge2F6 is substantial. The same trend to larger EAs on going down the Periodic Table is seen for the X2F5 systems, with 1.77 (C2F5), 2.68 (Si2F5), and 4.57 eV (Ge2F5). However, the EAs do not follow a monotonic trend for the X2F, X2F3 and X2F4 systems with respect to the series X = C, Si, Ge.