Abstract

Quantum-mechanical electronic structure computations were employed to explore initial steps towards a comprehensive design of polycarborane architectures through assembly of molecular units. Aspects considered were (i) the striking modification of geometrical parameters through substitution, (ii) endohedral carboranes and proposed ejection mechanisms for energy/ion/atom/energy storage/transport, (iii) the excited state character in single and dimeric molecular units, and (iv) higher architectural constructs. A goal of this work is to find optimal architectures where atom/ion/energy/spin transport within carborane superclusters is feasible in order to modernize and improve future photoenergy processes.

1. Introduction

The number of new borane and carborane molecules and their crystal structures has increased enormously due to the rich variety of synthetic routes discovered in the last decades [1, 2]. We have utilized a bottom-to-top approach and electronic structure computations to design icosahedral borane and carborane molecular architectures by starting with single molecules as well as with dimeric units connected directly or through bridges [316]. Knowledge of the electronic, magnetic, optical, and other properties of such single and dimeric units illuminates the evaluation of the global properties of 1D, 2D, or 3D superclusters and nanostructures [1722]. The use of light to modify properties of nanostructures is at the heart of current cutting-edge technology. Thus, the irradiation of photonic crystals can change their conducting properties [23] and control spin transport. [24, 25] We now use electronic structure computations to explore the first design stages of carborane-based architectures having potential uses in molecular and nano-photochemistry.

2. Computational Methodology

For the computations of the electronic structure in the compounds studied in this work, we employed different methods. As a benchmark for triplet and other excited states, we used the complete active space self-consistent field (CASSCF) [26] and their perturbative methods (CASPT2) [27], as implemented in the MOLCAS [28] program. For the optimization of geometries in the groundstate and triplet states, use was made of the well-known B3LYP hybrid which combines Hartree-Fock and density functional theory [2931] as implemented in the Gaussian03 program [32]. The Becke, three-parameter, Lee-Yang-Parr (B3LYP) hybrid exchange-correlation ( ) functional is constructed as a linear combination of the Hartree-Fock exact exchange functional ( ) and a given number of exchange and correlation explicit density functionals (x and c as subindex correspond to exchange and correlation respectively): = + 0.20 + 0.72 + 0.81 , where LDA and GGA stand for local density approximation and generalized gradient approximation. As regards to the computation of singlet-triplet energy gaps in the biradicals presented in this work, we used the CASSSCF/CASPT2 method as benchmark model—the multireference wave functions being based on an active space of six electrons distributed in six molecular orbitals (6,6)—together with the broken-symmetry (U)B3LYP [3335] approach: in the latter, the singlet-triplet energy gaps are estimated as = , where , , , and stand for energy of the broken-symmetry solution, energy of the triplet state, and expectation value of (square of spin operator) for the broken-symmetry and triplet states, respectively. In all studies we have used the 6-31 and 6-31+ basis sets, which are of double-zeta quality and contain polarization functions for all nonhydrogen atoms (6-31 ) and polarization plus diffuse functions for all nonhydrogen atoms (6-31+ ) respectively. Larger basis sets have been used when necessary for checking the validity and accuracy of the computed properties.

3. Results and Discussion

A goal in the nanoscience of nD-dimensional ( = ) networks made up of single 0D units (molecules/clusters/atoms/ions) is to ascertain the global or local properties of the system as function of the building units. What changes or modifications result from external perturbations (e.g., a electromagnetic field)? If a single unit can be modified, to what extent do the properties of a network built from such units change, if at all? We have addressed such questions by starting with the well-known basic icosahedral ortho, meta, and para carborane as well as borane units [36] to construct D dimensional networks and to explore their properties [18, 19, 3740]. While these compounds have been known since the 1960s, the properties of D dimensional networks made up thereof are less known. This section describes the changes in geometric and electronic properties, , that these systems undergo as a function of substituents , charges , spins , and the nature of the wave functions , and thus = ( ). We also consider endohedral and dimeric structures connected directly or with bridging units.

We have established that striking variations (over 1  !) in cage carbon-carbon distances ( ) occur in 1,2-disubtituted -carboranes and in the corresponding dianions (formed formally by double proton abstraction) [8]. Figure 1 shows distances ranging from 1.6  to 2.5  as a function of charge and substituent; the B3LYP/6-31 model was used in the computations. Here we should emphasize that the feeding of electronic charge into the cage causes a gradual transformation from a closo toward a nido, arachno, or hypho type of geometrical structure [41]. The triplet states of such systems [9] (see Figure 2) also reveal rich structural variety.

Figure 2 shows the optimized geometries of the triplet states of the same o-carboranes and the dianions considered in Figure 1 as singlet ground states; Figure 3 displays the corresponding spin density. The triplet states display interesting distortions and geometrical arrangements: for instance, a hexagonal window appears in the triplet state geometry of -carborane (Figure 2(a)) and the distance has increased by 1  ! The amino groups in triplet 1,2- -1,2- are essentially planar (Figure 2(b)) and the distance has increased by 0.7  . After double proton abstraction, the distances of the triplet dianions (Figures 2(c) and 2(d)) remain similar for o-carborane but increase up to 2.85 Å in the 1,2- -1,2- triplet state geometry (Figure 2(d)).

The spin density (Figure 3) in these clusters is concentrated mainly around the exo groups NH2 an and on the cage carbon atoms. For -carborane (Figure 2(a)), the spin density is spread around the “open” frontal hexagonal window (Figure 3(a)).

With respect to electronic structure properties, results of a comprehensive study [9] of the adiabatic singlet-triplet energy gaps— = ( )—in a series of substituted -carboranes are plotted in Figure 4. The adiabatic singlet-triplet energy gap is computed as the energy difference between the triplet and the singlet state optimized geometries: = . For the results displayed in Figure 4, computations at the (U)B3LYP/6-31 level of theory were used, the latter calibrated with the more accurate CASPT2 method. As some of the gaps are within the visible region of the electromagnetic spectrum, the corresponding systems should display photochemical absorption/emission properties, if intersystem crossing processes are plausible [42].

Inserting endohedral atoms/ions inside the carborane cage is another possibility for changing the properties of the single units forming the D networks [1012, 43, 44]. Although endohedral carboranes have not been realized experimentally, recent computational studies suggest than an insertion/ejection mechanism should be plausible using photon energies around 6.4 eV for Li+ and [10] (see Figure 5).

The implementation of ejection mechanisms similar to that shown in Figure 5, in an D-dimensional carborane network, opens up many opportunities for energy/atom/ion transport via photochemical molecular switches in molecular and nano-photochemistry.

The possibilities of using endohedral boranes and carboranes can be exploited further within nanoscience applications. For instance, the dramatic changes in the dipole moments of the dimeric units displayed in Figure 6, as a function of the position of the endohedral atom, call for potential applications such as molecular condensers or largely orientable molecular devices in electric fields.

These descriptions of (car)boranes from electronic structure computations illustrate the possibilities of changing their properties by the effects of exo substitution—specifically in -carboranes—or by endo insertion/ejection mechanisms. If such endohedral compounds can be prepared, they surely will respond to external perturbations, such as electromagnetic radiation. We are now exploring the consequences of extending systems based on the units considered above in different dimensions. What properties of finite (infinite) aggregates will result? The answer is not straightforward; several theoretical methodologies must be applied.

Biradical systems with nearly degenerate singlet-triplet energy gaps, such as those in Figure 7, are examples of how the properties of single units can change with spin as the parameter. The dimeric system built by joining two radicals by an acetylene bridge, C C , illustrates this principle as well. High-level quantum chemical computations showed that this biradical has a singlet ground state with a triplet state very slightly higher in energy ( 0.005 eV, on the order of at room temperature) [15]. Figure 8 shows the spin density of the triplet and the “broken-symmetry” singlet biradical states.

4. Conclusions

The properties of borane and carborane clusters can be varied by changing (i) the exo substituents, (ii) the wave function (e.g., singlet versus triplet states), and (iii) the endohedral atoms or ions inserted into monomeric or dimeric units, resulting in very large dipole moments. We also have predicted a theoretical mechanism for atom/ion ejection via a photochemical switch. Our electronic structure computations have established patterns for building higher architectural constructs starting from modified carborane building units. Our goal, to control atom, ion, spin, or energy transport inside borane and carborane networks in different dimensions, aims at the very heart of nanoscience and nanotechnology. The elaboration of higher architectures depends on experimental breakthroughs in assembling molecules, like those begun in the 1990s with 1D connections forming linear rods [18, 19]. Higher-dimensional constructs, for example, those in Figure 9, can include endohedral atoms inside the cages. The geometries can be modified by exo substituents. Abstracting hydrogen atoms also would change the properties as the spins of the resulting unpaired electrons might be delocalized in different directions. These possibilities, with the as yet unknown constructs in Figure 8, are research goals and starting points for the elaboration of finite or infinite networks capable of interacting with electromagnetic fields.

Acknowledgments

This research has been supported by the Spanish MEC/FEDER, projects MAT2006-13646-C03-02, CTQ2007-61260, and CSD2007-0010 Consolider-Ingenio in Molecular Nanoscience, and in the USA by the National Science Foundation Grant CHE-0716718, Georgia, and the Welch Foundation of Houston, Texas, through Grant BD-0894.