A new ab initio intermolecular potential energy surface and predicted rotational spectra of the Ne-H2S complex

We report a new three-dimensional ab initio intermolecular potential energy surface for the Ne-H(2)S complex with H(2)S monomer fixed at its experimental average structure. Using the supermolecular approach, the intermolecular potential energies were evaluated at CCSD(T) (coupled cluster with single and double and perturbative triple excitations) level with large basis sets including bond functions. The full counterpoise procedure was employed to correct the basis set superposition error. The planar T-shaped global minimum is located at the intermolecular distance of 3.51 ? with a well depth of 71.57 cm(-1). An additional planar local minimum was found to be separated from the global minimum with an energy barrier of 23.11 cm(-1). In addition, two first-order and one second-order saddle points were also located. The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to evaluate the rovibrational energy levels for eight isotopic species of the Ne-H(2)S complexes. The rotational transition frequencies for the eight isotopomers were also determined for the ground and first vibrational excited states, which are all in very good agreement with the available experimental values.     

Ab initio potential energy surface and rovibrational spectrum of Ar-HCCCN

We report an ab initio intermolecular potential energy surface of the Ar-HCCCN complex using a supermolecular method. The calculations were performed using the fourth-order M?ller-Plesset theory with the full counterpoise correction for the basis set superposition error and a large basis set including bond functions. The complex was found to have a planar T-shaped structure minimum and a linear minimum with the Ar atom facing the H atom. The T-shaped minimum is the global minimum with the well depth of 236.81 cm(-1). A potential barrier separating the two minima is located at R=5.57 A and theta=20.39 degrees with the height of 151.59 cm(-1). The two-dimensional discrete variable representation was employed to calculate the rovibrational energy levels for Ar-HCCCN. The rovibrational spectra including intensities for the ground state and the first excited intermolecular vibrational state are also presented. The results show that the spectra are mostly b-type (Delta K(a)=+/-1) transitions with weak a-type (Delta K(a)=0) transitions in structure, which are in good agreement with the recent experimental results [A. Huckauf, W. Jager, P. Botschwina, and R. Oswald, J. Chem. Phys. 119, 7749 (2003)].     

Ab initio and analytic intermolecular potentials for Ar-CH(3)OH

Ab initio calculations at the CCSD(T)/aug-cc-pVTZ level of theory were used to characterize the Ar-CH(3)OH intermolecular potential energy surface (PES). Potential energy curves were calculated for four different Ar + CH(3)OH orientations and used to derive an analytic function for the intermolecular PES. A sum of Ar-C, Ar-O, Ar-H(C), and Ar-H(O) two-body potentials gives an excellent fit to these potential energy curves up to 100 kcal mol(-1), and adding an additional r(-n) term to the Buckingham two-body potential results in only a minor improvement in the fit. Three Ar-CH(3)OH van der Waals minima were found from the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ calculations. The structure of the global minimum is in overall good agreement with experiment (X.-C. Tan, L. Sun and R. L. Kuczkowski, J. Mol. Spectrosc., 1995, 171, 248). It is T-shaped with the hydroxyl H-atom syn with respect to Ar. Extrapolated to the complete basis set (CBS) limit, the global minimum has a well depth of 0.72 kcal mol(-1) with basis set superposition error (BSSE) correction. The aug-cc-pVTZ basis set gives a well depth only 0.10 kcal mol(-1) smaller than this value. The well depths of the other two minima are within 0.16 kcal mol(-1) of the global minimum. The analytic Ar-CH(3)OH intermolecular potential also identifies these three minima as the only van der Waals minima and the structures predicted by the analytic potential are similar to the ab initio structures. The analytic potential identifies the same global minimum and the predicted well depths for the minima are within 0.05 kcal mol(-1) of the ab initio values. Combining this Ar-CH(3)OH intermolecular potential with a potential for a OH-terminated alkylthiolate self-assembled monolayer surface (i.e., HO-SAM) provides a potential to model Ar + HO-SAM collisions.     

Ab initio potential energy surface and bound states for the Kr-OCS complex

The first ab initio potential energy surface of the Kr-OCS complex is developed using the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)]. The mixed basis sets, aug-cc-pVTZ for the O, C, and S atom, and aug-cc-pVQZ-PP for the Kr atom, with an additional (3s3p2d1f) set of midbond functions are used. A potential model is represented by an analytical function whose parameters are fitted numerically to the single point energies computed at 228 configurations. The potential has a T-shaped global minimum and a local linear minimum. The global minimum occurs at R = 7.146 a(0), θ = 105.0° with energy of -270.73 cm(-1). Bound state energies up to J = 9 are calculated for three isotopomers (82)Kr-OCS, (84)Kr-OCS, and (86)Kr-OCS. Analysis of the vibrational wavefunctions and energies suggests the complex can exist in two isomeric forms: T-shaped and quasi-linear. The calculated transition frequencies and spectroscopic constants of the three isotopomers are in good agreement with the experimental values.     

Spectroscopy and potential energy surface of the H2-CO2 van der Waals complex: experimental and theoretical studies

A 4-D ab initio potential energy surface is calculated for the intermolecular interaction of hydrogen and carbon dioxide, using the CCSD(T) method with a large basis set. The surface has a global minimum with a well depth of 212 cm(-1) and an intermolecular distance of 2.98 A for a planar configuration with both the O-C-O and H-H axes perpendicular to the intermolecular axis. Bound state calculations are performed for the H(2)-CO(2) van der Waals complex with H(2) in both the para and ortho spin states, and the binding energy of paraH(2)-CO(2)(50.4 cm(-1)) is found to be significantly less than that of orthoH(2)-CO(2)(71.7 cm(-1)). The surface supports 7 bound intermolecular vibrational states for paraH(2)-CO(2) and 19 for orthoH(2)-CO(2), and the lower rotational levels with J< or = 4 follow an asymmetric rotor pattern. The calculated infrared spectrum of paraH(2)-CO(2) agrees well with experiment. For orthoH(2)-CO(2), the ground state rotational levels allowed by symmetry are found to have (K(a), K(c))=(even, odd) or (odd, even). This somewhat unexpected fact enables the previously observed experimental spectrum to be assigned for the first time, in good agreement with theory, and indicates that the orientation of hydrogen is perpendicular to the intermolecular axis in the ground state of the orthoH(2)-CO(2) complex.     

Density functional theory study of nine-atom germanium clusters: effect of electron count on cluster geometry

Density functional theory (DFT) at the hybrid B3LYP level has been applied to the germanium clusters Ge(9)(z) clusters (z = -6, -4, -3, -2, 0, +2, and +4) starting from three different initial configurations. Double-zeta quality LANL2DZ basis functions extended by adding one set of polarization (d) and one set of diffuse (p) functions were used. The global minimum for Ge(9)(2)(-) is the tricapped trigonal prism expected by Wade's rules for a 2n + 2 skeletal electron structure. An elongated tricapped trigonal prism is the global minimum for Ge(9)(4)(-) similar to the experimentally found structure for the isoelectronic Bi(9)(5+). However, the capped square antiprism predicted by Wade's rules for a 2n + 4 skeletal electron structure is only 0.21 kcal/mol above this global minimum indicating that these two nine-vertex polyhedra have very similar energies in this system. Tricapped trigonal prismatic structures are found for both singlet and triplet Ge(9)(6)(-), with the latter being lower in energy by 3.66 kcal/mol and far less distorted. The global minimum for the hypoelectronic Ge(9) is a bicapped pentagonal bipyramid. However, a second structure for Ge(9) only 4.54 kcal/mol above this global minimum is the C(2)(v)() flattened tricapped trigonal prism structure found experimentally for the isoelectronic Tl(9)(9)(-). For the even more hypoelectronic Ge(9)(2+), the lowest energy structure consists of an octahedron fused to two trigonal bipyramids. For Ge(9)(4+), the global minimum is an oblate (squashed) pentagonal bipyramid with two pendant Ge vertices.     

A new ab initio potential energy surface for the NeCO complex with the vibrational coordinate dependence

A new high quality three-dimensional potential energy surface for the Ne-CO van der Waals complex is developed using the CCSD(T) method and avqz∕avqz+33221 basis set. The ab initio calculation is performed in a total of 1365 configurations with supermolecule method. There is a single global minimum located in a nearly T-shaped geometry. The global minimum energy is -49.4090 cm(-1) at R(e)=6.40a(0) and θ(e)=82.5(°) for V(00). Using the three-dimensional potential energy surface, we have calculated bound rovibrational energy levels up to J = 10 including the Coriolis coupling terms. Compared with the experimental transition frequencies, the theoretical results are in good agreement with the experimental results.     

Theoretical investigation of the He-I2(E3Πg) ion-pair state: ab initio intermolecular potential and vibrational levels

We present a theoretical study on the potential energy surface and vibrational bound states of the E electronic excited state of the HeI(2) van der Waals system. The interaction energies are computed using accurate ab initio methods and large basis sets. Relativistic small-core effective core potentials in conjunction with a quintuple-zeta quality basis set are employed for the heavy iodine atoms in multireference configuration interaction calculations for the (3)A' and (3)A" states. For the representation of the potential energy surface we used a general interpolation technique for constructing potential surfaces from ab initio data based on the reproducing kernel Hilbert space method. The surface presents global and local minima for T-shaped configurations with well-depths of 33.2 and 4.6 cm(-1), respectively. Vibrational energies and states are computed through variational quantum mechanical calculations. We found that the binding energy of the HeI(2)(E) T-shaped isomer is 16.85 cm(-1), in excellent agreement with recent experimental measurements. In lieu of more experimental data we also report our predictions on higher vibrational levels and we analyze the influence of the underlying surface on them. This is the first attempt to represent the potential surface of such a highly excited electronic state of a van der Waals complex, and it demonstrates the capability of the ab initio technology to provide accurate results for carrying out reliable studies to model experimental data.     

A new potential energy surface and predicted infrared spectra of He-CO2: dependence on the antisymmetric stretch of CO2

A new potential energy surface involving the antisymmetric Q(3) normal mode of CO(2) for the He-CO(2) van der Waals complex is constructed at the coupled-cluster singles and doubles with noniterative inclusion of connected triple [CCSD(T)] level with augmented correlation-consistent quadruple-zeta (aug-cc-pVQZ) basis set plus bond functions. Two vibrationally adiabatic potentials with CO(2) at both the ground and the first excited vibrational states are generated from the integration of the three-dimensional potential over the Q(3) coordinate. The potential has a T-shaped global minimum and two equivalent linear local minima. The bound rovibrational energy levels are obtained using the radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm. The observed band origin shift of the complex (0.0946 cm(-1)) is successfully reproduced by our calculation (0.1034 cm(-1)). The infrared spectra of the complex are also predicted. The fundamental band is in excellent agreement with the experiment. Most of the transitions corresponding to the observed hot band [M. J. Weida et al., J. Chem. Phys. 101, 8351 (1994)] are assigned reasonably.     

Ab initio study of the intermolecular potential energy surface in the ion-induced-dipole hydrogen-bonded O2(-)(X2Πg)-H2(X1Σg(+)) complex

This work presents the first investigation on the intermolecular potential energy surface of the ground electronic state of the O2(-)(2Πg)-H2(1Σg(+)) complex. High level correlated ab initio calculations were carried out using the Hartree-Fock spin-unrestricted coupled cluster singles and doubles including perturbative triples correction [RHF-UCCSD(T)]/aug-cc-pVXZ levels of calculations, where XZ = DZ, TZ, QZ, and 5Z. Results of full geometry optimization and the intermolecular potential energy surface (IPES) calculations show four equivalent minimum energy structures of L-shaped geometry with Cs symmetry at equilibrium along the 2A″ surface of the complex. For these equilibrium minimum energy structures, the most accurate value for the dissociation energy (De) was calculated as 1407.7 cm(-1), which was obtained by extrapolating the counterpoise (CP) corrected De values to the complete basis set (CBS) limit. This global minimum energy structure is stabilized by an ion-induced-dipole hydrogen bond. Detailed investigations of the IPES show that the collinear structure is unstable, while the C2v geometries present saddle points along the 2A″ surface. The barrier height between the two equivalent structures that differs in whether the hydrogen-bonded hydrogen atom is above or below the axis that connects centers of masses of the H2 and O2(-) moieties within the complex was calculated as 70 cm(-1). This suggests that the complex exhibits large amplitude motion. The barrier height to rotation of the H2 moiety by 180° within the complex is 1020 cm(-1). Anharmonic oscillator calculations predicted a strong H-H stretch fundamental transition at 3807 cm(-1). Results of the current work are expected to stimulate further theoretical and experimental investigations on the nature of intermolecular interactions in complexes that contain the superoxide radical and various closed-shell molecules that model atmospheric and biological molecules. These studies are fundamental to understanding the role of the O2(-) anion in chemistry in the atmosphere and in biological systems.     

Relative stability of planar versus double-ring tubular isomers of neutral and anionic boron cluster B20 and B20-

High-level ab initio molecular-orbital methods have been employed to determine the relative stability among four neutral and anionic B20 isomers, particularly the double-ring tubular isomer versus three low-lying planar isomers. Calculations with the fourth-order Moller-Plessset perturbation theory [MP4(SDQ)] and Dunning's correlation consistent polarized valence triple zeta basis set as well as with the coupled-cluster method including single, double, and noniteratively perturbative triple excitations and the 6-311G(d) basis set show that the double-ring tubular isomer is appreciably lower in energy than the three planar isomers and is thus likely the global minimum of neutral B20 cluster. In contrast, calculations with the MP4(SDQ) level of theory and 6-311+G(d) basis set show that the double-ring anion isomer is appreciably higher in energy than two of the three planar isomers. In addition, the temperature effects on the relative stability of both 10B20- and 11B20- anion isomers are examined using the density-functional theory. It is found that the three planar anion isomers become increasingly more stable than the double-ring isomer with increasing the temperature. These results are consistent with the previous conclusion based on a joint experimental/simulated anion photoelectron spectroscopy study [B. Kiran et al., Proc. Natl. Acad. Sci. U.S.A. 102, 961 (2005)], that is, the double-ring anion isomer is notably absent from the experimental spectra. The high stability of the double-ring neutral isomer of B20 can be attributed in part to the strong aromaticity as characterized by its large negative nucleus-independent chemical shift. The high-level ab initio calculations suggest that the planar-to-tubular structural transition starts at B20 for neutral clusters but should occur beyond the size of B20- for the anion clusters.     

Ab initio investigation of the NH(X)-N2 van der Waals complex

The NH-N(2) van der Waals complex has been examined at the CCSD(T) level of theory using aug-cc-pVDZ and aug-cc-pVTZ basis sets. The full basis set superposition error correction was applied. Two minimum energy structures were located for the electronic ground state. The global minimum corresponds to a linear geometry of the complex (NH-N-N), with D(e)=236 cm(-1) and R(c.m.)=4.22 A. The secondary minimum corresponds to a T-shaped geometry of C(2v) symmetry, where the nitrogen atom of the H-N moiety points toward the center of mass of the N(2) unit, aligned with the a-inertial axis of the complex. The binding energy and R(c.m.) value for the secondary minimum were 144 cm(-1) and 3.63 A, respectively. This potential energy surface is consistent with the properties of matrix-isolated NH-N(2), and it is predicted that linear NH-N(2) will be a stable complex in the gas phase at low temperatures.     

Ab initio intermolecular potential energy surface, bound states, and microwave spectra for the van der Waals complex Ne-HCCCN

We report two ab initio intermolecular potential energy surfaces for Ne-HCCCN using a supermolecular method. The calculations were performed at the fourth-order M?ller-Plesset (MP4) and the coupled cluster singles-and-doubles with noniterative inclusion of connected triples [CCSD(T)] levels with the full counterpoise correction for the basis set superposition error and a large basis set including bond functions. The complex was found to have a planar T-shaped structure minimum and a linear minimum with the Ne atom facing the H atom. The two-dimensional discrete variable representation method was employed to calculate the rovibrational bound states. In addition, the microwave spectra including intensities for the ground vibrational state were predicted. The results show that the spectrum is dominated by b-type (DeltaK(a) = +/-1) transitions with very weak a-type (DeltaK(a) = 0) transitions. The calculated results at the CCSD(T) potential are in good agreement with those at MP4 potential.     

Rotational spectra and structure of the Ar2-H2S complex: pulsed nozzle Fourier transform microwave spectroscopic and ab initio studies

This paper reports the rotational spectrum and structure of the Ar2-H2S complex and its HDS and D2S isotopomers. The ground state structure has heavy-atom C2v symmetry with the two Ar atoms indistinguishable and H2S freely rotating as evinced by the fact that asymmetric top energy levels with Kp=odd levels are missing. The rotational constants for the parent isotopomer are: A=1733.115(1) MHz, B=1617.6160(5) MHz and C=830.2951(2) MHz. Unlike the Ar-H2S complex, the Ar2-H2S does not show an anomalous isotopic shift in rotational constants on deuterium substitution. However, the intermolecular potential is still quite floppy, leading to very different centrifugal distortion constants for the three isotopomers. The Ar-Ar and Ar-c.m.(H2S) distances are determined to be 3.820 A and 4.105 A, respectively. The A rotational constants for Ar2-H2S/HDS/D2S isotopomers are very close to each other and to the B constant of free Ar2, indicating that H2S does not contribute to the moment of inertia about the a-axis. Ab initio calculations at MP2 level with aug-cc-pVQZ basis set lead to an equilibrium C2v minimum structure with the Ar-Ar line perpendicular to the H-H line and the S away from Ar2. The centrifugal distortion constants, calculated using the ab initio force field, are in reasonable agreement with the experimental values. However, they do not show the variation observed for different isotopmers. The binding energy of Ar2-H2S has been determined to be 507 cm-1(6.0 kJ mol-1) by CBS extrapolation after correcting for basis set superposition error. Potential energy scans point out that the barrier for internal rotation of H2S about its b axis is only 10 cm-1 and it is below the zero point energy (13.5 cm-1) in this torsional degree of freedom. Internal rotation of H2S about its a- and c-axes also have small barriers of about 50 cm-1 only, suggesting that H2S is extremely floppy within the complex.     

A new "spectroscopic" potential energy surface for formaldehyde in its ground electronic state

We report a new "spectroscopic" potential energy surface (PES) of formaldehyde (H(2)(12)C(16)O) in its ground electronic state, obtained by refining an ab initio PES in a least-squares fitting to the experimental spectroscopic data for formaldehyde currently available in the literature. The ab initio PES was computed using the CCSD(T)/aug-cc-pVQZ method at 30 840 geometries that cover the energy range up to 44 000 cm(-1) above equilibrium. Ro-vibrational energies of formaldehyde were determined variationally for this ab initio PES by means of the program TROVE [Theoretical ROtation-Vibration Energies; S. N. Yurchenko, W. Thiel, and P. Jensen, J. Mol. Spectrosc. 245, 126 (2007)]. The parameter values in the analytical representation of the PES were optimized in fittings to 319 ro-vibrational energies with J = 0, 1, 2, and 5. The initial parameter values in the fittings were those of the ab initio PES, the ro-vibrational eigenfunctions obtained from this PES served as a basis set during the fitting process, and constraints were imposed to ensure that the refined PES does not deviate unphysically from the ab initio one in regions of configuration space not sampled by the experimental data. The resulting refined PES, referred to as H(2)CO-2011, reproduces the available experimental J ≤ 5 data with a root-mean-square error of 0.04 cm(-1).     

Fourier transform spectroscopy and direct potential fit of a shelflike state: application to E(4)1Σ(+) KCs

The paper presents high-resolution experimental study and a direct potential construction of a shelflike state E(4)(1)Σ(+) of the KCs molecule converging to K(4(2)S) + Cs(5(2)D) atomic limit; such data are of interest for selecting optical paths for producing and monitoring cold polar diatomics. The collisionally enhanced laser induced fluorescence (LIF) spectra corresponding to both spin-allowed E(4)(1)Σ(+) → X(1)(1)Σ(+) and spin-forbidden E(4)(1)Σ(+) → a(1)(3)Σ(+) transitions of KCs were recorded in visible region by Fourier transform spectrometer with resolution of 0.03 cm(-1). Overall about 1650 rovibronic term values of the E(4)(1)Σ(+) state of (39)K(133)Cs and (41)K(133)Cs isotopologues nonuniformly covering the energy range [16987, 18445] cm(-1) above the minimum of the ground X-state were determined with the uncertainty of 0.01 cm(-1). Experimental data field is limited by vibrational levels v' ∈ [2, 74] with rotational quantum numbers J' ∈ [1, 188]. The closed analytical form for potential energy curve (PEC) based on Chebyshev polynomial expansion (CPE) was implemented to a direct potential fit (DPF) of the experimental term values of the most abundant (39)K(133)Cs isotopologue. Besides analyticity, regularity, correct long-range behavior, and nice convergence properties, the CPE form demonstrated optimal balance on flexibility and constraint for the DPF of a shelflike state aggravated by a limited data set. The mass-invariant properties of the CPE PEC were tested by the prediction of rovibronic term values of the (41)K(133)Cs isotopomer which coincided with their experimental counterparts with standard deviation of 0.0048 cm(-1). The CPE modeling is compared with the highly flexible pointwise inverted perturbation approach model, as well as with conventional Dunham analysis of restricted data set v' ≤ 50. Reliability of the empirical PEC is additionally confirmed by good agreement between the calculated and experimental relative intensity distributions in the long E(v') → X(v") LIF progressions.     

Potential energy surface, van der Waals motions, and vibronic transitions in phenol-argon complex

The structure and intermolecular vibrational energy levels of the phenol-Ar complex are calculated from its potential energy surface. This surface is constructed from a large set of the interaction energy values computed using second-order Moller-Plesset perturbation theory with the augmented correlation consistent polarized valence double-zeta basis set. The global minimum in the potential energy surface corresponds to a cluster structure with Ar located over the geometric center of the phenol ring at a distance of 3.510 A and shifted by 0.1355 A towards oxygen. The calculated dissociation energy of 371 cm(-1) is in accordance with the experiment. Additional local minima higher in energy are with Ar placed in the phenol plane. However, they are too shallow to form the bound states corresponding to planar isomers. The deformation of the potential energy surface shape, created by the interaction of Ar with the phenolic oxygen, is responsible for a pronounced intermode mixing. As a result, a set of hybrid stretching-bending states appears which cannot be described in terms of the standard models. The intermode coupling is reflected in the vibronic structure of the S1-S0 electronic transition. The intensities of the vibronic bands are calculated from the electronic transition dipole moment surfaces determined using the ab initio single-excitation configuration interaction method. They allow us to correct and complete the assignment of the spectra observed in phenol-Ar, as well as in the analogous complexes of phenol with Kr and Xe.     

Nitrous oxide dimer: a new potential energy surface and rovibrational spectrum of the nonpolar isomer

The spectrum of nitrous oxide dimer was investigated by constructing new potential energy surfaces using coupled-cluster theory and solving the rovibrational Schro?dinger equation with a Lanczos algorithm. Two four-dimensional (rigid monomer) global ab initio potential energy surfaces (PESs) were made using an interpolating moving least-squares (IMLS) fitting procedure specialized to describe the interaction of two linear fragments. The first exploratory fit was made from 1646 CCSD(T)/3ZaP energies. Isomeric minima and connecting transition structures were located on the fitted surface, and the energies of those geometries were benchmarked using complete basis set (CBS) extrapolations, counterpoise (CP) corrections, and explicitly correlated (F12b) methods. At the geometries tested, the explicitly correlated F12b method produced energies in close agreement with the estimated CBS limit. A second fit to 1757 data at the CCSD(T)-F12b/VTZ-F12 level was constructed with an estimated fitting error of less than 1.5?cm(-1). The second surface has a global nonpolar O-in minimum, two T-shaped N-in minima, and two polar minima. Barriers between these minima are small and some wave functions have amplitudes in several wells. Low-lying rovibrational wave functions and energy levels up to about 150?cm(-1) were computed on the F12b PES using a discrete variable representation/finite basis representation method. Calculated rotational constants and intermolecular frequencies are in very close agreement with experiment.     

A photoelectron spectroscopy and ab initio study of B21-: negatively charged boron clusters continue to be planar at 21

The structures and chemical bonding of the B(21)(-) cluster have been investigated by a combined photoelectron spectroscopy and ab initio study. The photoelectron spectrum at 193 nm revealed a very high adiabatic electron binding energy of 4.38 eV for B(21)(-) and a congested spectral pattern. Extensive global minimum searches were conducted using two different methods, followed by high-level calculations of the low-lying isomers. The global minimum of B(21)(-) was found to be a quasiplanar structure with the next low-lying planar isomer only 1.9 kcal/mol higher in energy at the CCSD(T)/6-311-G* level of theory. The calculated vertical detachment energies for the two isomers were found to be in good agreement with the experimental spectrum, suggesting that they were both present experimentally and contributed to the observed spectrum. Chemical bonding analyses showed that both isomers consist of a 14-atom periphery, which is bonded by classical two-center two-electron bonds, and seven interior atoms in the planar structures. A localized two-center two-electron bond is found in the interior of the two planar isomers, in addition to delocalized multi-center σ and π bonds. The structures and the delocalized bonding of the two lowest lying isomers of B(21)(-) were found to be similar to those in the two lowest energy isomers in B(19)(-).