Design of Carborane Molecular Architectures via Electronic Structure Computations
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.
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 [3–16]. 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 [17–22]. 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  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)  and their perturbative methods (CASPT2) , as implemented in the MOLCAS  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 [29–31] as implemented in the Gaussian03 program . 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 [33–35] 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  to construct D dimensional networks and to explore their properties [18, 19, 37–40]. 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) . 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 . The triplet states of such systems  (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  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 .
Inserting endohedral atoms/ions inside the carborane cage is another possibility for changing the properties of the single units forming the D networks [10–12, 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  (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, CC, 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) . Figure 8 shows the spin density of the triplet and the “broken-symmetry” singlet biradical states.
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.
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.
T. B. Marder and Z. Lin, Eds., Contemporary Metal Boron Chemistry I: Borylenes, Boryls, Borane Sigma-Complexes, and Borohydrides, Structure and Bonding, Springer, Berlin, Germany, 2008.
C. Viñas, G. Barberà, J. M. Oliva, F. Teixidor, A. J. Welch, and G. M. Rosair, “Are halocarboranes suitable for substitution reactions? The case for 3-I-1,2-closo-: molecular orbital calculations, aryldehalogenation reactions, NMR interpretation of closo-carboranes, and molecular structures of 1-Ph-3-Br-1,2-closo- and 3-Ph-1,2-closo-,” Inorganic Chemistry, vol. 40, no. 26, pp. 6555–6562, 2001.View at: Google Scholar
J. Llop, C. Viñas, J. M. Oliva et al., “C obaltabisdicarbollide anion —as doping agent on intelligent membranes for ion capture,” Journal of Organometallic Chemistry, vol. 657, pp. 232–238, 2002.View at: Google Scholar
J. M. Oliva, F. Teixidor, and C. Viñas, “C-C electronic interaction in 1,2-substituted 1,2- ortho-carboranes: to bond or not to bond?” in Conceptual Aspects of Chemical Bonding, abstracts of papers of the American Chemical Society 226th National Meeting: 664-ORGN part 2, September 2003.View at: Google Scholar
F. Teixidor, G. Barberà, A. Vaca et al., “Are methyl groups electron-donating or electron-withdrawing in boron clusters? Permethylation of o-carborane,” Journal of the American Chemical Society, vol. 127, no. 29, pp. 10158–10159, 2005.View at: Google Scholar
J. M. Oliva, N. L. Allan, P. V. R. Schleyer, C. Viñas, and F. Teixidor, “Strikingly long distances in 1,2-disubstituted ortho-carboranes and their dianions,” Journal of the American Chemical Society, vol. 127, pp. 13538–13547, 2005.View at: Google Scholar
J. M. Oliva and L. Serrano-Andrés, “A computational study of the lowest singlet and triplet states of neutral and dianionic 1,2-substituted icosahedral and octahedral o-carboranes,” Journal of Computational Chemistry, vol. 27, no. 4, pp. 524–535, 2006.View at: Google Scholar
L. Serrano-Andrés and J. M. Oliva, “Photochemical window mechanism for controlled atom release in carborane endohedral boxes: Theoretical evidence,” Chemical Physics Letters, vol. 432, no. 1–3, pp. 235–239, 2006.View at: Google Scholar
J. I. Burgos, L. Serrano-Andrés, J. M. Oliva, and D. J. Klein, “On the effect of radical character, substitution and atom encapsulation on the volume of icosahedral (car)boranes,” Afinidad, vol. 65, no. 533, pp. 32–38, 2008.View at: Google Scholar
L. Serrano-Andrés, D. J. Klein, P. V. R. Schleyer, and J. M. Oliva, “What electronic structures and geometries of carborane mono- and ortho-, meta-, and para-diradicals are preferred?” Journal of Chemical Theory and Computation, vol. 4, no. 8, pp. 1338–1347, 2008.View at: Google Scholar
J. M. Oliva, D. J. Klein, P. V. R. Schleyer, and L. Serrano-Andrés, “Design of carborane molecular architectures with electronic structure computations: from endohedral and polyradical systems to multidimensional networks,” Pure and Applied Chemistry, vol. 81, no. 4, pp. 719–729, 2009.View at: Publisher Site | Google Scholar
O. Crespo, C. Gimeno, A. Laguna, I. Ospino, G. Aullón, and J. M. Oliva, “Organometallic gold complexes of carborane. Theoretical comparative analysis of ortho, meta, and para derivatives and luminescence studies,” Dalton Transactions, vol. 19, pp. 3807–3813, 2009.View at: Publisher Site | Google Scholar
H. Hiura and T. Kanayama, “Carborane superclusters formed by ion–molecule reactions in an ion trap,” Journal of Molecular Structure, vol. 735-736, pp. 367–374, 2005.View at: Google Scholar
W. Jiang, D. E. Harwell, M. D. Mortimer, C. B. Knobler, and M. F. Hawthorne, “Palladium-catalyzed coupling of ethynylated P-carborane derivatives—synthesis and structural characterization of modular ethynylated P-carborane molecules,” Inorganic Chemistry, vol. 35, no. 15, pp. 4355–4359, 1996.View at: Google Scholar
J. Müller, K. Baše, T. F. Magnera, and J. Michl, “Rigid-rod oligo--carboranes for molecular tinkertoys. An inorganic Langmuir-Blodgett film with a functionalized outer surface,” Journal of the American Chemical Society, vol. 114, p. 9721, 1992.View at: Google Scholar
J. N. Hohman, P. Zhang, E. I. Morin et al., “Self-assembly of carboranethiol isomers on Au111: intermolecular interactions determined by molecular dipole orientations,” ACS Nano, vol. 3, no. 3, pp. 527–536, 2009.View at: Google Scholar
S. Khatua, J. M. Guerrero, K. Claytor et al., “Micrometer-scale translation and monitoring of individual nanocars on glass,” ACS Nano, vol. 3, no. 2, pp. 351–356, 2009.View at: Google Scholar
J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the Flow of Light, Princeton University Press, Princeton, NJ, USA, 2nd edition, 2008.
N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and B. J. van Wees, “Electronic spin transport and spin precession in single graphene layers at room temperature,” Nature, vol. 448, pp. 571–574, 2007.View at: Google Scholar
B. O. Roos, “Ab initio methods in quantum chemistry,” in Advances in Chemical Physics, K. P. Lawley, Ed., vol. 69, pp. 399–445, John Wiley & Sons, New York, NY, USA, 1987.View at: Google Scholar
M. Merchán and L. Serrano-Andrés, “Ab initio methods for excited states,” in Computational Photochemistry, M. Olivucci, Ed., vol. 81, pp. 743–754, Elsevier, Amsterdam, The Netherlands, 2005.View at: Google Scholar
K. Andersson, M. Barysz, A. Bernhardsson et al., MOLCAS, Version 5.0., Department of Theoretical Chemistry, Lund, Sweden, 2000.
A. D. Becke, “A new mixing of Hartree-Fock and local density-functional theories,” The Journal of Chemical Physics, vol. 98, no. 2, pp. 1372–1377, 1993.View at: Google Scholar
J. P. Perdew, M. Ernzerhof, and K. Burke, “Rationale for mixing exact exchange with density functional approximations,” Journal of Chemical Physics, vol. 105, p. 9982, 1996.View at: Google Scholar
S. H. Vosko, L. Wilk, and M. Nusair, “Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis,” Canadian Journal of Physics, vol. 58, pp. 1200–1211, 1980.View at: Google Scholar
M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., GAUSSIAN 03, Gaussian, Wallingford, Conn, USA, 2004.
L. Noodleman, C. Y. Peng, D. A. Case, and J. M. Mouesca, “Orbital interactions, electron delocalization, and spin interactions in iron-sulfur clusters,” Chemical Reviews, vol. 144, pp. 199–244, 1995.View at: Google Scholar
A. Ovchinnikov and J. K. Labanowski, “Simple spin correction of unrestricted density-functional calculation,” Physical Review A, vol. 53, no. 6, pp. 3946–3952, 1996.View at: Google Scholar
K. Yamaguchi, T. Kawasaki, Y. Takano, Y. Kigarawa, Y. Yamashita, and H. Fujita, “Electronic excitations and first hyperpolarizability of 2-methyl-4-nitroaniline clusters,” International Journal of Quantum Chemistry, vol. 90, pp. 359–370, 2002.View at: Google Scholar
R. N. Grimes, Carboranes, Academic Press, New York, NY, USA, 1970.
X. Yang, W. Jiang, C. B. Knobler, and M. F. Hawthorne, “Chemical structures and schematic representations of a crown ether,” Journal of the American Chemical Society, vol. 114, p. 9719, 1992.View at: Google Scholar
B. T. King, B. C. Noll, A. J. McKinley, and J. Michl, “Dodecamethylcarba-closo-dodecaboranyl (CB(11)Me(12)(.)), a stable free-radical,” Journal of the American Chemical Society, vol. 118, no. 44, pp. 10902–10903, 1996.View at: Google Scholar
L. Eriksson, K. Vyakaranam, J. Ludvik, and J. Michl, “Synthesis of the isolable biradicals and trans-,” The Journal of Organic Chemistry, vol. 72, no. 7, pp. 2351–2356, 2007.View at: Google Scholar
I. Zharov, A. Saxena, J. Michl, and J. Miller, “C-hydroxy-closo-dodecaboranes,” Inorganic Chemistry, vol. 36, p. 6033, 1997.View at: Google Scholar
J. Michl and V. Bonacic-Koutecky, Electronic Aspects of Organic Photochemistry, Wiley-Interscience, New York, NY, USA, 1990.
E. D. Jemmis and M. M. Balakrishnarajan, “Ab initio predictions on novel stuffed polyhedral boranes,” Journal of the American Chemical Society, vol. 122, no. 30, pp. 7392–7393, 2000.View at: Google Scholar
O. P. Charkin, N. M. Klimenko, D. Moran, A. M. Mebel, D. O. Charkin, and P. V. R. Schleyer, “Theoretical study of complexes of closo-borane, alane, and gallane anions with cations of light metals inside and outside of icosahedral clusters (, Al, and Ga),” Journal of Physical Chemistry A, vol. 106, no. 47, pp. 11594–11602, 2002.View at: Publisher Site | Google Scholar