Ab initio prediction of Boron compounds arising from Borozene: Structural and electronic properties

Structure and electronic properties of two unusual boron clusters obtained by fusion of borozene rings has been studied by means of first principles calculations, based on the generalized-gradient approximation of the density functional theory, and the semiempirical tight-binding method was used for the transport calculations. The role of disorder has also been considered with single vacancies and substitutional atoms. Results show that the pure boron clusters are topologically planar and characterized by (3c-2e) bonds, which can explain, together with the aromaticity (estimated by means of NICS), the remarkable cohesive energy values obtained. Such feature makes these systems competitive with the most stable boron clusters to date. On the contrary, the introduction of impurities compromises stability and planarity in both cases. The energy gap values indicate that these clusters possess a semiconducting character, while when the larger system is considered, zero-values of the density of states are found exclusively within the HOMO-LUMO gap. Electron transport calculations within the Landauer formalism confirm these indications, showing semiconductor-like low bias differential conductance for these stuctures. Differences and similarities with Carbon clusters are highlighted in the discussion.Comment: 10 pages, 2 tables, 5 figure

Temperature Sensitive Nanocapsule of Complex Structural Form for Methane Storage

The processes of methane adsorption, storage and desorption by the nanocapsule are investigated with molecular-dynamic modeling method. The specific nanocapsule shape defines its functioning uniqueness: methane is adsorbed under 40 MPa and at normal temperature with further blocking of methane molecules the K@C601+ endohedral complex in the nanocapsule by external electric field, the storage is performed under normal external conditions, and methane desorption is performed at 350 K. The methane content in the nanocapsule during storage reaches 11.09 mass%. The nanocapsule consists of tree parts: storage chamber, junction and blocking chamber. The storage chamber comprises the nanotube (20,20). The blocking chamber is a short nanotube (20,20) with three holes. The junction consists of the nanotube (10,10) and nanotube (8,8); moreover, the nanotube (8,8) is connected with the storage chamber and nanotube (10,10) with the blocking chamber. The blocking chamber is opened and closed by the transfer of the K@C601+ endohedral complex under electrostatic field action

Influence of a knot on the strength of a polymer strand

Many experiments have been done to determine the relative strength of different knots, and these show that the break in a knotted rope almost invariably occurs at a point just outside the `entrance' to the knot. The influence of knots on the properties of polymers has become of great interest, in part because of their effect on mechanical properties. Knot theory applied to the topology of macromolecules indicates that the simple trefoil or `overhand' knot is likely to be present with high probability in any long polymer strand. Fragments of DNA have been observed to contain such knots in experiments and computer simulations. Here we use {\it ab initio} computational methods to investigate the effect of a trefoil knot on the breaking strength of a polymer strand. We find that the knot weakens the strand significantly, and that, like a knotted rope, it breaks under tension at the entrance to the knot.Comment: 3 pages, 4 figure

Combinatorial–computational–chemoinformatics (C3) approach to finding and analyzing low-energy tautomers

Finding the most stable tautomer or a set of low-energy tautomers of molecules is critical in many aspects of molecular modelling or virtual screening experiments. Enumeration of low-energy tautomers of neutral molecules in the gas-phase or typical solvents can be performed by applying available organic chemistry knowledge. This kind of enumeration is implemented in a number of software packages and it is relatively reliable. However, in esoteric cases such as charged molecules in uncommon, non-aqueous solvents there is simply not enough available knowledge to make reliable predictions of low energy tautomers. Over the last few years we have been developing an approach to address the latter problem and we successfully applied it to discover the most stable anionic tautomers of nucleic acid bases that might be involved in the process of DNA damage by low-energy electrons and in charge transfer through DNA. The approach involves three steps: (1) combinatorial generation of a library of tautomers, (2) energy-based screening of the library using electronic structure methods, and (3) analysis of the information generated in step (2). In steps 1–3 we employ combinatorial, computational and chemoinformatics techniques, respectively. Therefore, this hybrid approach is named “Combinatorial*Computational*Chemoinformatics”, or just abbreviated as C3 (or C-cube) approach. This article summarizes our developments and most interesting methodological aspects of the C3 approach. It can serve as an example how to identify the most stable tautomers of molecular systems for which common chemical knowledge had not been sufficient to make definite predictions

Quantum mechanical polar surface area

A correlation has been established between the absorbed fraction of training-set molecules after oral administration in humans and the Quantum Mechanical Polar Surface Area (QMPSA). This correlation holds for the QMPSA calculated with structures where carboxyl groups are deprotonated. The correlation of the absorbed fraction and the QMPSA calculated on the neutral gas phase optimized structures is much less pronounced. This suggests that the absorption process is mainly determined by polar interactions of the drug molecules in water solution. Rules are given to derive the optimal polar/apolar ranges of the electrostatic potential

A vibrational circular dichroism implementation within a Slater-type-orbital based density functional framework and its application to hexa- and hepta-helicenes

We describe the implementation of the rotational strengths for vibrational circular dichroism (VCD) in the Slater-type orbital based Amsterdam Density Functional (ADF) package. We show that our implementation, which makes use of analytical derivative techniques and London atomic orbitals, yields origin independent rotational strengths. The basis set dependence in the particular case of Slater-type basis functions is also discussed. It turns out that the triple zeta STO basis sets with one set of polarization functions (TZP) are adequate for VCD calculations. The origin- dependence of the atomic axial tensors is checked by a distributed origin gauge implementation. The distributed and common origin gauge implementations yield virtually identical atomic axial tensors with the Slater-type basis sets employed here, proving that our implementation yields origin independent rotational strengths. We verify the implementation for a set of benchmark molecules, for which the dependence of the VCD spectra on the particular choice of the exchange–correlation functional is studied. The pure functionals BP86 and OLYP show a particularly good performance. Then, we apply this approach to study the VCD spectra of hexa- and hepta- helicenes. In particular we focus on relationships between the sign of the rotational strengths of the two helicenes

Effect of Peierls transition in armchair carbon nanotube on dynamical behaviour of encapsulated fullerene

The changes of dynamical behaviour of a single fullerene molecule inside an armchair carbon nanotube caused by the structural Peierls transition in the nanotube are considered. The structures of the smallest C20 and Fe@C20 fullerenes are computed using the spin-polarized density functional theory. Significant changes of the barriers for motion along the nanotube axis and rotation of these fullerenes inside the (8,8) nanotube are found at the Peierls transition. It is shown that the coefficients of translational and rotational diffusions of these fullerenes inside the nanotube change by several orders of magnitude. The possibility of inverse orientational melting, i.e. with a decrease of temperature, for the systems under consideration is predicted.Comment: 9 pages, 6 figures, 1 tabl

Exploitation of symmetry in periodic Self-Consistent-Field ab initio calculations: application to large three-dimensional compounds

Symmetry can dramatically reduce the computational cost (running time and memory allocation) of Self-Consistent-Field ab initio calculations for crystalline systems. Crucial for running time is use of symmetry in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the SACO (Symmetry Adapted Crystalline Orbital) basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. We here illustrate the effectiveness of this scheme, following recent advancements in the CRYSTAL code, concerning memory allocation and direct basis set transformation. Quantitative examples are given for large unit cell systems, such as zeolites (all-silica faujasite and silicalite MFI) and garnets (pyrope). It is shown that the full SCF of 3D systems containing up to 576 atoms and 11136 Atomic Orbitals in the cell can be run with a hybrid functional on a single core PC with 500 MB RAM in about 8 h. © 2014 Science China Press and Springer-Verlag Berlin Heidelberg