The stability of cage and ring isomers for carbon and boron nitride clusters

Mechanistic Insights into the Formation of Lithium Fluoride Nanotubes

Born-Oppenheimer molecular dynamics (BOMD) and periodic density functional theory (DFT) calculations have been applied for describing the mechanism of formation of lithium fluoride (LiF) nanotubes with cubic, hexagonal, octagonal, decagonal, dodecagonal, and tetradecagonal cross-sections. It has been shown that high energy structures, such as nanowires, nanorings, nanosheets, and nanopolyhedra are transient species for the formation of stable nanotubes. Unprecedented (LiF)_n clusters (n≤12) were also identified, some of them lying less than 10 kcal mol^-[1] above their respective global minima. Such findings indicate that stochastic synthetic techniques, such as laser ablation and chemical vapor deposition, should be combined with a template-driven procedure in order to generate the nanotubes with adequate efficiency. Apart from the stepwise growth of LiF units, the formation of nanotubes was also studied by rolling up a planar square sheet monolayer, which could be hypothetically produced from the exfoliation of the FCC crystal structure. It was shown that both pathways could lead to the formation of alkali halide nanotubes, a still unprecedented set of one-dimensional materials.

The chemical space of B, N-substituted polycyclic aromatic hydrocarbons: Combinatorial enumeration and high-throughput first-principles modeling

Combinatorial introduction of heteroatoms in the two-dimensional framework of aromatic hydrocarbons opens up possibilities to design compound libraries exhibiting desirable photovoltaic and photochemical properties. Exhaustive enumeration and first-principles characterization of this chemical space provide indispensable insights for rational compound design strategies. Here, for the smallest seventy-seven Kekulean-benzenoid polycyclic systems, we reveal combinatorial substitution of C atom pairs with the isosteric and isoelectronic B, N pairs to result in 7 453 041 547 842 (7.4 tera) unique molecules. We present comprehensive frequency distributions of this chemical space, analyze trends, and discuss a symmetry-controlled selectivity manifestable in synthesis product yield. Furthermore, by performing high-throughput ab initio density functional theory calculations of over thirty-three thousand (33k) representative molecules, we discuss quantitative trends in the structural stability and inter-property relationships across heteroarenes. Our results indicate a significant fraction of the 33k molecules to be electronically active in the 1.5-2.5 eV region, encompassing the most intense region of the solar spectrum, indicating their suitability as potential light-harvesting molecular components in photo-catalyzed solar cells.

Possible mechanism of BN fullerene formation from a boron cluster: Density-functional tight-binding molecular dynamics simulations

We simulate the formation of a BN fullerene from an amorphous B cluster at 2000 K by quantum mechanical molecular dynamics based on the density-functional tight-binding method. We run 30 trajectories 200 ps in length, where N atoms are supplied around the target cluster, which is initially an amorphous B36 cluster. Most of the incident N atoms are promptly incorporated into the target cluster to form B-N-B bridges or NB3 pyramidal local substructures. BN fullerene formation is initiated by alternating BN ring condensation. Spontaneous atomic rearrangement and N2 dissociation lead to the construction of an sp(2) single-shelled structure, during which the BN cluster undergoes a transition from a liquid-like to a solid-like state. Continual atomic rearrangement and sporadic N2 dissociation decrease the number of defective rings in the BN cluster and increase the number of six-membered rings, forming a more regular shell structure. The number of four-membered rings tends to remain constant, and contributes to more ordered isolated-tetragon-rule ring placement.

Cubic C8 : An Observable Allotrope of Carbon?

Ab initio and DFT calculations are used to investigate the structure, electronic properties, spectra and reactivity of cubic C8 , which is predicted to be aromatic according to Hirsch's rule. Although highly strained and with a small amount of diradical character, the carbon cube represents a surprisingly deep minimum and should therefore be observable as an isolated molecule. It is, however, predicted to be very reactive, both with itself and triplet oxygen. Calculated IR, Raman, and UV/Vis spectra are provided to aid identification of cubic C8 should it be synthesized.

Boron-nitride and aluminum-nitride "Pringles" and flapping motion

Motivated by the recent successful synthesis of a new nanocarbon, namely, a warped, double-concave graphene "Pringle" (Nat. Chem., 2013, 5, 739), we investigate properties of warped boron-nitride (BN) and aluminum-nitride (AlN) analogues, i.e., the non-planar B40N40H30 and Al40N40H30 "Pringles" using density functional theory (DFT) calculations. Particular attention is placed on the effect of non-hexagonal rings on the stability and physical properties of BN and AlN Pringles. We find that the warped BN and AlN Pringles with one pentagon and five heptagons are stable without imaginary frequencies. Both the warped B40N40H30 and Al40N40H30 Pringles are expected to be flexible in solution as both can periodically change their shape in a dynamic "flapping" fashion due to their much lower activation barrier of racemization compared to that of the C80H30 counterpart. Since the warped B40N40H30 possesses a smaller HOMO-LUMO gap than the planar B39N39H30, it is expected that incorporating non-hexagonal ring defects by design can be an effective way to modify electronic properties of BN-based nanoplates.

A new generation of B(n)N(n) rings as a supplement to boron nitride tubes and cages

In B(n)N(n) cages or tubes, when the quasi-borazine rings are attached to each other through a pair of common atoms of B and N, the bonding structure is named class A. On the other hand, there are some B(n)N(n) rings including a covalent bond between two atoms of B and N, which are named class B. In all previous studies, both reports of synthesis and theoretical calculation of boron nitride tubes and cages, the quasi-borazine units are attached together like class A. There are no theoretical or experimental reports from class B compounds except for a brief study in our previous works (Struct. Chem. 2012, 23, 551-580; J. Phys. Chem. C 2010, 114, 15315.). In this study, we have used two kinds of boron nitride rings from a twisted BN sheet in the same chirality created by different mechanisms. For (4, 4) chirality, the molecules B(16)N(16) and B(15)N(15) are found to respectively represent class A and B, and for (5, 5) chirality the molecules B(20)N(20) and B(18)N(18) are respectively again of class A and B. The structure of class A rings is similar to boron nitride tubes, but we have shown that it is impossible to produce a macromolecule of class B form as tubes or cages, because there is much more instability and intermolecular tension in macro forms of class B. This is the main reason that the class B molecules are rare and, because of their small size, have not yet been synthesized, although we have some suggestions for the synthesis of these kinds of molecules. The stability and electromagnetic properties with hybrid density functional theory using the EPR-III and EPR-II basis sets for explanation of hyperfine parameters and spin densities, electrical potential, and isotropic Fermi coupling constant of these rings have been studied by the nonbonded interaction models. Normal mode analyses including aromaticity have been investigated by using the nucleus independent chemical shift values at the ring center. Interaction energy and gain in energy aid in describing the stability that is promoted upon gradual binding with molecular hydrogen, and a linear relationship occurred between them.

Advances in computational studies of energy materials

We review recent developments and applications of computational modelling techniques in the field of materials for energy technologies including hydrogen production and storage, energy storage and conversion, and light absorption and emission. In addition, we present new work on an Sn2TiO4 photocatalyst containing an Sn(II) lone pair, new interatomic potential models for SrTiO3 and GaN, an exploration of defects in the kesterite/stannite-structured solar cell absorber Cu2ZnSnS4, and report details of the incorporation of hydrogen into Ag2O and Cu2O. Special attention is paid to the modelling of nanostructured systems, including ceria (CeO2, mixed Ce(x)O(y) and Ce2O3) and group 13 sesquioxides. We consider applications based on both interatomic potential and electronic structure methodologies; and we illustrate the increasingly quantitative and predictive nature of modelling in this field.

Boron nitride cages from B12N12 to B36N36: square-hexagon alternants vs boron nitride tubes

The structures and stabilities of square-hexagon alternant boron nitrides (Bx Nx , x=12-36) vs their tube isomers containing octagons, decagons and dodecagons have been computed at the B3LYP density functional level of theory with the correlation-consistent cc-pVDZ basis set of Dunning. It is found that octagonal B20N20 and B24N24 tube structures are more stable than their square-hexagon alternants by 18.6 and 2.4 kcal mol(-1), respectively, while the square-hexagon alternants of other cages are more stable. Trends in stability as a function of cluster size are discussed.