The shape effect and its consequences for polar surfaces and for heterogeneous catalysis.

In this paper we develop the shape effect, which is relevant for crystalline materials whose size is larger than that of the thermodynamic limit. According to this effect the electronic properties of one surface of a crystal depend upon all of its surfaces, i.e. on the overall shape. At first, qualitative mathematical arguments are presented for the existence of this effect based on the conditions for the stability of polar surfaces. Our treatment explains why such surfaces are observed even though earlier theory indicated that they should not exist. Then, models are developed from which it is found computationally that changing the shape of a polar crystal can substantially alter the magnitude of its surface charges. Apart from surface charges, it follows that the crystal shape will also significantly affect bulk properties, most notably polarization and piezoelectric responses. Additional model calculations show a strong shape effect on the activation energy for heterogeneous catalysis primarily through local surface charges rather than a non-local/long range electrostatic potential.

Properties of Naked Silver Clusters with Up to 100 Atoms as Found with Embedded-Atom and Density-Functional Calculations.

The structural and energetic properties of small silver clusters Agn with n = 2-100 atoms are reported. For n = 2-100 the embedded atom model for the calculation of the total energy of a given structure in combination with the basin-hopping search strategy for an unbiased structure optimization has been used to identify the energies and structures of the three energetically lowest-lying isomers. These optimized structures for n = 2-11 were subsequently studied further through density-functional-theory calculations. These calculations provide additional information on the electronic properties of the clusters that is lacking in the embedded-atom calculations. Thereby, also quantities related to the catalytic performance of the clusters are studied. The calculated properties in comparison to other available theoretical and experimental data show a good agreement. Previously unidentified magic (i.e., particularly stable) clusters have been found for n>80. In order to obtain a more detailed understanding of the structural properties of the clusters, various descriptors are used. Thereby, the silver clusters are compared to other noble metals and show some similarities to both copper and nickel systems, and also growth patterns have been identified. All vibrational frequencies of all the clusters have been calculated for the first time, and here we focus on the highest and lowest frequencies. Structural effects on the calculated frequencies were considered.

Mechanistic insights into aerobic oxidative cleavage of glycol catalyzed by an Anderson-type polyoxometalate [IMo6O24]5.

A computational investigation of the aerobic oxidative C-C bond cleavage reaction of glycol catalyzed by an Anderson-type heteropolyanion HPA [IMo6O24]5- in the presence of acetonitrile as solvent has been performed at the WB97XD/6-31G(d,p)/lanl2dz level. Two reaction pathways have been identified. The catalytic cycle of each pathway consists of three steps: oxidation cleavage of a glycol molecule by the HPA, oxidation of the HPA by one dioxygen molecule, and, finally, oxidation of a second glycol and regeneration of the catalyst. These reaction pathways have been thoroughly investigated in terms of energetic, natural bond orbital (NBO), natural charges, and geometrical parameters. It is found that (i) even though the top oxygen atoms of the Anderson heteropolyanion are not the most negatively charged ones, they are more likely to react with the diol hydroxyl groups, (ii) a direct relationship between the presence of the iodine ion I(VII) and the studied oxidation reaction could not be identified, and (iii) in terms of energy, the transfer of the two hydrogen atoms is the most energetic step.

DFT evaluation of structural, electronic and variation properties for complex carbohydrates with biological interest.

The synthetic bicyclic bis(hemiacetals) compounds 1,5-pyranose-9,7-pyranoses, with a structural analogy to the bicyclic monosaccharide Bradyrhizose, have been described here based on a theoretical approach, using DFT calculations with the B3LYP functional combined with the 6-311 + G(d,p) basis set. First, we have performed a geometrical and electronic properties description of (1 R,9S), (1S,9S) and (1S,9R)-1,5-pyranose-9,7-pyranoses. Results analysis indicated that, slight differences in the three-dimensional orientations of their atoms lead to an enormous difference in chemical reactivity. Consequently, (1S,9S) and (1S,9R) isomers are predicted to be the most resembling the natural bradyrhizose in structural features. To enhance the performance of these two isomers, a set of modifications through functional groups attached to the reactive sites were determined by local reactivity descriptors. Subsequently, in order to get more information on the obtained derivatives for both isomers, HOMO, LUMO, Egap and four electronic parameters were calculated and compared. The substituted systems show a good performance in chemical reactivity than the unmodified parent compounds.Communicated by Ramaswamy H. Sarma.

Nonlinear optical properties DFT calculations of polyacethylene and copolymers models substituted with aldimines chromophores as side chains.

Aldimine derivatives chromophores grafted on polyacetylene oligomers were first designed to investigate the nonlinear optical (NLO) response of the resulting materials using CAMB3LYP method. The effects of different factors such as the chain length separating the substitutions as well as the configurations and the orientations of these latter, were examined and discussed. In a second part of this paper, NLO responses, in particular the static hyperpolarizabilities of polyvinyl oligomers substituted by aldimine chromophores via carboxylic groups and by pyrrolidone groups were calculated for different configurations. The stability and the hydrogen bonding in each oligomer were also discussed.

Efficient model photosensitizers based on metallocenyl complexes with thiophene-N = N-pyrimidine as π-conjugated bridge and cyanoacrylate as an anchoring group: a density functional theory study.

Eight push-pull systems involving containing four transition metals (iron, ruthenium, cobalt, and nickel), metallocenes as donor groups, cyanoacrylate as electron attractor group, and thiophene-N = N- pyrimidine derivatives as π-conjugated bridges were designed and studied using DFT and TD-DFT methods involving B3LYP and CAM-B3LYP functionals combined with the cc-pVDZ/LANL2DZ basis sets. The main purpose of this work is to determine the effect of metallocene in improving the photosensitization property of such chromophores. This was done by calculating their light-harvesting efficiency LHE as well as other properties employed for DSSC application. The considered dyes were first studied in the gas phase, then in the presence of TiO2 nanoparticles representing the semi-conductor, and finally in the presence of a specific implicit solvent. The presence of iron as metal involved in the metallocene group supplemented by extending the π-conjugated bridge by a cyanovinyl spacer was demonstrated so as to give the most optimal response taking into account the lower cost and toxicity as well as the friendliness to the environment of iron as metal.

Role of the Backbone when Optimizing Functional Groups─A Theoretical Study Based on an Improved Inverse-Design Approach.

We present an improved inverse-design approach for automatically identifying molecular (or other) systems with optimal values for prechosen properties. The new approach uses SMILES (simplified molecular input line entry system) to describe molecular structures efficiently, a genetic algorithm to optimize the molecules automatically, and the DFTB+ (self-consistent charge density functional tight-binding) method to calculate electronic properties. Thereby, almost every class of materials─even macromolecules or monomers─can be studied easily. Without crossover operators but with only mutation operators, the genetic algorithm is more adaptive to SMILES while keeping its efficiency. DFTB+ is more accurate than the DFTB method used in our previous inverse-design approach for the study of excited states and charge transfer processes. The improved approach is applied to optimize benzene, pyridine, pyridazine, pyrimidine, and pyrazine derivatives for seven electronic properties, which all are highly relevant and important for the performance of molecules in solar cells. We found that for some electronic properties, the precise composition and structure of the backbone have remarkable impacts on the value of the electronic properties and/or on the set of functional groups that leads to the best performance. On the contrary, for other properties, these effects are less pronounced. The reasonable optimal functional groups and/or substitution patterns are reported.

Theoretical Study on Non-Linear Optics Properties of Polycyclic Aromatic Hydrocarbons and the Effect of Their Intercalation with Carbon Nanotubes.

Results of a theoretical study devoted to comparing NLO (non-linear optics) responses of derivatives of tetracene, isochrysene, and pyrene are reported. The static hyperpolarizability ß, the dipole moment µ, the HOMO and LUMO orbitals, and their energy gap were calculated using the CAM-B3LYP density functional combined with the cc-pVDZ basis set. The para-disubstituted NO2-tetracene-N(CH3)2 has the highest NLO response, which is related to a large intramolecular charge transfer. Adding vinyl groups to the para-disubstituted NO2-tetracene-N(CH3)2 results in an increase in the NLO responses. We further investigated the effect of the intercalation of various push-pull molecules inside an armchair single-walled carbon nanotube. The intercalation leads to increased NLO responses, something that depends critically on the position of the guest molecule and/or on functionalization of the nanotube by donor and attractor groups.

Kinetics and Thermodynamics of CO Oxidation by (TiO2)6.

Molecular level insights into the mechanism and thermodynamics of CO oxidation by a (TiO2)6 cluster have been obtained through density functional calculations. Thereby, in this study, as an example, two different structural isomers of (TiO2)6 are considered with the purpose of understanding the interplay between local structure and activity for the CO oxidation reaction. Active sites in the two isomeric forms were identified on the basis of global and local reactivity descriptors. For the oxidation of CO to CO2, the study considered both sequential and simultaneous adsorption of CO and O2 on (TiO2)6 cluster through the ER and LH mechanisms, respectively. Three different pathways were obtained for CO oxidation by (TiO2)6 cluster, and the mechanistic route of each pathway were identified by locating the transition-state and intermediate structures. The effect of temperature on the rate of the reaction was investigated within the harmonic approximation. The structure-dependent activity of the cluster was rationalized through reactivity descriptors and analysis of the frontier orbitals.

Optimizing small conjugated molecules for solar-cell applications using an inverse-design method.

Small organic conjugated molecules are key elements for low-cost photovoltaic devices. One example is cyanopyridone molecules. By modifying these molecules, for instance through optimally chosen functional groups attached to the backbone, their properties can be improved. However, the very large number of possible modifications makes it difficult to identify the best performing molecules. In the present work, we have used a computational inverse-design approach (PooMa) to identify the positions and types of functional groups attached to a modified cyanopyridone that lead to the best performance in solar-energy harvesting. A QSPR model based on five electronic descriptors has been used to determine the properties of solar cells. Our approach uses a genetic algorithm to search the chemical space containing 184 (104,976) substituted cyanopyridone systems and predicts out of those the best 20 molecules with optimal performance efficiencies (PCE). PooMa uses the Density-Functional Tight-Binding (DFTB) method for calculating the electronic properties. DFTB is a fast method with acceptable accuracy and, therefore, can be used on a normal desktop without expensive hard- or software. In order to get further information about our suggested systems, a DFT method and its derivative TD-DFT are applied.

Competition between tubular, planar and cage geometries: a complete picture of structural evolution of Bn (n = 31-50) clusters.

Stimulated by the early theoretical prediction of B80 fullerene and the experimental finding of the B40 cage, the structures of medium-sized boron clusters have attracted intensive research interest during the last decade, but a complete picture of their size-dependent structural evolution remains a puzzle. Using a genetic algorithm combined with density-functional theory calculations, we have performed a systematic global search for the low-lying structures of neutral Bn clusters with n = 31-50. Diverse structural patterns, including tubular, quasi-planar, cage, core-shell, and bilayer, are demonstrated for the ground-state Bn clusters; for certain cluster sizes, unprecedented geometries are predicted for the first time. Their stabilities at finite temperatures are evaluated, and the competition mechanism between various patterns is elucidated. Chemical bonding analysis reveals that the availability of localized σ bonds and delocalized π bonds in the Bn clusters play a key role in their structural stability. Our results provide important insights into the bonding pattern and growth behavior of medium-sized boron clusters, which lay the foundation for experimental design and synthesis of boron nanostructures.

Theoretical study of the mechanism behind the site- and enantio-selectivity of C-H functionalization catalysed by chiral dirhodium catalyst.

The C-H functionalization is very important for the synthesis of pharmaceuticals and complex natural products. Rhodium carbenoids, obtained when a dirhodium(ii) catalyst containing a crown formed by chiral ligands reacts with diazo compounds with both an electron donating group and an electron withdrawing group, play an important part in controlling site- and enantio-selectivity for functionalization of non-activated C-H bonds. It has earlier been demonstrated that the tertiary C-H bond is more favored to be functionalized inside the crown of the dirhodium catalyst with S-configuration ligands compared with the secondary and primary C-H bonds although the latter possess weaker steric effects. We argue that the higher site- and enantio-selectivity for some types of C-H bond functionalization can be related to intermolecular hydrogen bonding, steric hindrance, and weak interactions when the dirhodium catalyst is interacting with the chiral ligands.

Ultranarrow heterojunctions of armchair-graphene nanoribbons as resonant-tunnelling devices.

A systematic investigation is performed on the electronic transport properties of armchair-graphene nanoribbon (AGNR) heterojunctions using spin-polarized density functional theory calculations in combination with the non-equilibrium Green's function formalism. 9-AGNR and 5-AGNR structures are used to form a single-well configuration by sandwiching a 5-AGNR between two 9-AGNRs. At the same time, these 9-AGNRs are matched at the left and right to electrodes, 9 and 5 being the number of carbon dimers as width. This heterojunction mimics an electronic device with two potential barriers (9-AGNR) and one quantum well (5-AGNR) where quasi-bound states are confined. First, we study the ground state properties, and then we calculate the electron transport properties of this device as a function of the well width. We show the presence of electronic tunnelling resonances between the barriers by delocalized electron density inside the well's structure. This is corroborated by transmission curves, localized densities of states (LDOS), current-vs.-bias voltage results, and the trend of the resonances as a function of the well width. This work shows that carbon AGNRs may be used as resonant-tunnelling devices for applications in nanoelectronics.

Density Functional Theory Descriptors for Ionic Liquids and the Introduction of a Coulomb Correction.

As a result of continuing ionic liquid research, it becomes clearer that charge transfer in ionic liquids has a physical reality. In a recent publication, we demonstrated the utility of simple density functional theory descriptors to estimate charge transfer for a large number of ion combinations, which is possible because the ions are treated separately. A major disadvantage found was that the charge transfer was systematically overestimated. In this work, we introduce a correction to account for the losses in Coulomb attraction when charge is transferred from the anion to the cation. We find that accounting for these losses is important to describe charge transfer in ionic liquids appropriately. The advantage that the calculations can be performed separately on the individual, isolated ions is maintained. The corrected as well as the uncorrected charge transfer have been calculated for over 4000 cation-anion combinations at the R(O)B3LYP/6-311+G(2d,p)//RB3LYP/6-31+G(d,p) level of theory. With the correction, the absolute values for the charge transfer are no longer unrealistically high and agree well with other charge transfer estimates from the literature. In general, the cumulative nature of the Haven ratio is now correctly mirrored in the relationship between the corrected theoretical charge transfer and the experimental estimate from the Nernst-Einstein relation. Earlier findings on the similarities between ether-functionalized and nonfunctionalized ionic liquids are confirmed. However, we also observe inconsistencies when using the experimental charge transfer estimates together with the ionicity interpretation of the Haven ratio. These can be interpreted as a hint toward the latter premise being wrong.

Temperature and isomeric effects in nanoclusters.

The canonical thermodynamics of clusters is determined quantum mechanically in the general case of multiple minima of the potential energy surface (PES) using the superposition approximation. As an illustration of the consequences of our approach, we study in detail the thermodynamic properties of CuN clusters with N from 2 to 150 as a function of cluster size, temperature, and number of isomers. Thereby, for instance, solid-solid transition temperatures for several cluster sizes are determined. We show that the cluster vibrations have a strong impact on the stability of the clusters and that this can be observed not only at medium and high temperatures but also at low temperatures and even at T = 0 K. Thus, including zero-point corrections can change the relative energetic ordering of different isomers. This isomeric effect results in an oscillatory dependence of the heat capacity on cluster size at moderate and high temperatures. Moreover, for the identification of magic clusters at non-vanishing temperature, the Helmholtz free energy is analyzed as a function of cluster size and temperature within a one-, two-, and three-minima model of the PES. Thereby, we demonstrate that the concept of magic clusters is strongly temperature dependent so that in several cases clusters change from being magic or non-magic at T = 0 K to be the opposite at non-vanishing temperature. We emphasize that all these effects are not specific for copper clusters alone but can also be observed in other metal or semiconductor nanoclusters.

Application of an inverse-design method to optimizing porphyrins in dye-sensitized solar cells.

Dye-sensitized solar cells (DSSCs) have attracted much interest during the past few decades. However, it is still a tremendous challenge to identify organic molecules that give an optimal power conversion efficiency (PCE). Here, we apply our recently developed, inverse-design method for this issue with the special aim of identifying porphyrins with promisingly high PCE. It turns out that the calculations lead to the prediction of 15 new molecules with optimal performances and for which none so far has been studied. These porphyrin derivatives will in the near future be synthesized and subsequently tested experimentally. Our inverse-design approach, PooMa, is based on the strategy of providing suggestions for molecular systems with optimal properties. PooMa has been developed as a tool that requires minimal resources and, therefore, builds on various approximate methods. It uses genetic algorithm to screen thousands (or often more) of molecules. For each molecule, the density-functional tight-binding (DFTB) method is used for calculating the electronic properties. In the present work, five different electronic properties are determined, all of which are related to optical performance. Subsequently, a quantitative structure-property relationship (QSPR) model is constructed that can predict the PCE through those five electronic properties. Finally, we benchmark our results through more accurate DFT calculations that give further information on the predicted optimal molecules.

Density Functional Theory Descriptors for Ionic Liquids and the Charge-Transfer Interpretation of the Haven Ratio.

One of the few properties common to all ionic liquids is their inherent electrical conductivity, which makes them promising candidates for advanced electrochemical applications. A central finding in this respect is that the measured conductivity is almost always lower than the one obtained from the Nernst-Einstein relation. There has been much dispute about whether correlated motion, charge transfer, or some sort of aggregation is the reason for this difference. In this work, we apply density functional theory-based descriptors to estimate the charge transfer in ionic liquids, which allow predictions for a large number of systems with minimal effort. The theoretical charge transfer was obtained from vertical ionization potentials and electron affinities at the RB3LYP/6-311+G(2d,p)//RB3LYP/6-31+G(d,p) level of theory. To be able to compare and classify the values obtained with this approach, another measure for charge transfer, available directly from the Nernst-Einstein relation, is introduced. The two quantities show significant correlation for some subsets of ionic liquids for which a sufficient amount of information is available. Additionally, the purely theoretical charge-transfer values allow for identifying interesting systems that should be the subject of further investigation.

A structural DFT study of MM, GG, MG, and GM alginic acid disaccharides and reactivity of the MG metallic complexes.

The density functional theory method using the B3LYP/6-31G(d,p) level of theory was used to perform isoenergetic maps in order to determine the lower energy conformers of four disaccharides constituting alginic acids, which are based on ß-D-mannuronic (M) and α-L-guluronic acid (G), called MM, GG, MG, and GM. The preferred structures are combined to monovalent (Li+, Na+, and K+) cations and further fully optimized, and an isoenergetic map corresponding to the complex (MG2-, 2Na+) was performed. Then, the reactivity of MG complexes with mono- and bivalent cations was studied using the global nucleophilic index. The position selectivity was also predicted using the local nucleophilic indices. It was demonstrated that experimental trends of relative reactivity and regioselectivity of the complexes are correctly predicted using these empirical indices of reactivity. Graphical abstract MM, GG, MG, and GM alginic acid disaccharides and reactivity of the MG metallic complexes.

Structure and stability of propellane-like E[Formula: see text].

Propellane-like structures are an important building block for group-14 compounds. Therefore, in the present work, we study theoretically various such structures of the form E[Formula: see text] starting with the pure molecules for which E = E' = E″ = C, Si, Ge, or Sn. Subsequently, we study the systems with E ≠ E' = E″ and finally, we consider some selected cases molecules for which E, E', and E″ all are different. Special emphasis is put on identifying structural trends for the molecules with at least two different group-14 elements. The resulting scheme is, we believe, generally valid for group-14 based systems with more than one type of group-14 elements.

From properties to materials: An efficient and simple approach.

We present an inverse-design method, the poor man's materials optimization, that is designed to identify materials within a very large class with optimized values for a pre-chosen property. The method combines an efficient genetic-algorithm-based optimization, an automatic approach for generating modified molecules, a simple approach for calculating the property of interest, and a mathematical formulation of the quantity whose value shall be optimized. In order to illustrate the performance of our approach, we study the properties of organic molecules related to those used in dye-sensitized solar cells, whereby we, for the sake of proof of principle, consider benzene as a simple test system. Using a genetic algorithm, the substituents attached to the organic backbone are varied and the best performing molecules are identified. We consider several properties to describe the performance of organic molecules, including the HOMO-LUMO gap, the sunlight absorption, the spatial distance of the orbitals, and the reorganisation energy. The results show that our method is able to identify a large number of good candidate structures within a short time. In some cases, chemical/physical intuition can be used to rationalize the substitution pattern of the best structures, although this is not always possible. The present investigations provide a solid foundation for dealing with more complex and technically relevant systems such as porphyrins. Furthermore, our "properties first, materials second" approach is not limited to solar-energy harvesting but can be applied to many other fields, as briefly is discussed in the paper.