Exploring odd-even effects of simple oligomer-like DRCNnT series: a study based on density functional theory/time-dependent density functional theory calculations

Odd-even effects of short-circuit current density and power conversion efficiency (PCE) are an interesting phenomenon in some organic solar cells. Although some explanations have been given, why they behave in such a way is still an open question. In the present work, we investigate a set of acceptor-donor-acceptor simple oligomer-like small molecules, named the DRCNnT (n = 5-9) series, to give an insight into this phenomenon because the solar cells based on them have high PCE (up to 10.08%) and show strong odd-even effects in experiments. By modeling the DRCNnT series and using density functional theory, we have studied the ground-state electronic structures of the DRCNnT (n = 5-9) series in condensed phase. The calculated results reproduce the experimental trends well. Furthermore, we find that the exciton-binding energies of the DRCNnT series may be one of the key parameters to explain this phenomenon because they also show odd-even effects. In addition, by studying the effects of alkyl branch and terminal group on odd-even effects of dipole moment, we find that eliminating one or two alkyl branches does not break the odd-even effects of dipole moments, but eliminating one or two terminal groups does. Finally, we conclude that removing one alkyl branch close to the terminal group of DRCN5T can decrease highest occupied molecular orbital (HOMO) energy (thus increasing open circuit voltage) and increase dipole moment (thus enhancing charge separation and short-circuit current). This could be a new and simple method to increase the PCE of DRCN5T-based solar cells.     

Effects of Ionic Liquids on the Thermodynamics of Hydrogen Activation by Frustrated Lewis Pairs: A Density Functional Theory Study**

Nowadays, hydrogen activation by frustrated Lewis pairs (FLPs) and their applications are one of the emerging research topics in the field of catalysis. Previous studies have shown that the thermodynamics of this reaction is determined by electronic structures of FLPs and solvents. Herein, we investigated systems consisting of typical FLPs and ionic liquids (ILs), which are well known by their large number of types and excellent solvent effects. The density functional theory (DFT) calculations were performed to study the thermodynamics for H₂ activation by both inter- and intra-molecular FLPs, as well as the individual components. The results show that the computed overall Gibbs free energies in ILs are more negative than that computed in toluene. Through the thermodynamics partitioning, we find that ILs favor the H−H cleavage elemental step over the elemental steps of proton attachment, hydride attachment and zwitterionic stabilization. Moreover, the results show that these effects are strongly dependent on the type of FLPs, where intra-molecular FLPs are more affected compared to the inter-molecular FLPs.     

Diagrammatic valence bond studies on hemocyanin

Summary The Diagrammatic Valence Bond studies on the active sites of hemocyanin, consisting of two Cu(I) ions and an oxygen molecule, are performed to find out the stable geometrical pattern and electronic structure. Different parameters used in this theoretical approach are taken from existing literature on highT _c superconductors. Attempts have been made to find out the differences in electronic structure of [Cu₂O₂]⁺² and [Cu₂O₂N₄]⁺² as it is observed that coordination of nitrogen ligand do affect electronic structure i.e. spin excitation gaps and charge and spin density distribution. A comparison of our results with earlier theoretical results are also presented.     

How Can Rotaxanes Be Modified by Varying Functional Groups at the Axle?—A Combined Theoretical and Experimental Analysis of Thermochemistry and Electronic Effects

We present theoretically as well as experimentally determined thermochemical data of the non-covalent interactions in different axle-substituted pseudorotaxanes. The overall interaction energy lies in the region of 35 kJ mol(-1), independent of the substitution pattern at the axle. Because rearrangement energies of 7 and 3 kJ mol(-1) are required for wheel and axle, respectively, the sum of the net interactions of individual non-covalent bonds must exceed 10 kJ mol(-1) to achieve a successful host-guest interaction. The geometrical analysis shows three hydrogen bonds, and the close inspection of the individual dipole moments as well as the individual hydrogen bonds reveals trends according to the different functional groups at the axle. The individual trends for the different hydrogen bonds almost lead to a cancellation of the substitution effects. From solvent-effect considerations it can be predicted that the pseudorotaxane is stable in CHCl(3) and CH(2)Cl(2), whereas it would dethread in water. Comparing experimentally and theoretically calculated Gibbs free enthalpies, we find reasonable agreement if an exchange reaction of one solvent molecule instead of the direct formation reaction is considered.     

Chlorine‐Enhanced Surface Mobility of Au(100)

Motivated by experimental studies of two‐dimensional Ostwald ripening on Au(100) electrodes in chlorine‐containing electrolytes, we have studied diffusion processes using density functional theory. We find that chlorine has a propensity to temporary form AuCl complexes, which diffuse significantly faster than gold adatoms. With and without chlorine, the lowest activation energy is found for the exchange mechanism. Chlorine furthermore reduces the activation energy for the detachment from kink sites. Kinetic Monte Carlo simulations were performed on the basis of extensive density functional theory calculations. The island‐decay rate obtained from these Monte Carlo simulations, as well as the decay rate obtained from the theoretical activation energies and frequency factors when inserted into analytical solutions for Ostwald ripening, are in agreement with experimental island‐decay rates in chlorine‐containing electrolytes.     

Systematic studies on the computation of nuclear magnetic resonance shielding constants and chemical shifts: the density functional models

We present a systematic density functional investigation on the prediction of the 13C, 15N, 17O, and 19F NMR properties of 23 molecules with 21 density functionals. Extensive comparisons are made for both 13C magnetic shieldings and chemical shifts with respect to the gas phase experimental data and the best CCSD(T) results. We find that the OPBE and OPW91 exchange-correlation functionals perform significantly better than some popular functionals such as B3LYP and PBE1PBE, even surpassing, in many cases, the standard wavefunction-based method MP2. Further analysis has been performed to explore the individual role played by various exchange and correlation functionals. We find that the B88 and PBE exchange functionals have a too strong tendency of deshielding, leading to too deshielded magnetic shielding constants; whereas the OPTX exchange functional performs remarkably well. We claim that the main source of error arises from the exchange functional, but correlation functional also makes important contribution. We find that the correlation functionals may be grouped into two classes. class A, such as LYP and B98, leads to deshielded NMR values, deteriorating the overall performance; whereas class B, such as PW91 and PBE, generally increases the absolute shieldings, which complements the exchange functionals, leading to improved results in the calculation of NMR data.     

Correlation of bond metallicity measures to electronegativity for binary oxides

We have investigated alkali, alkaline‐earth, and rutile binary oxides within density functional theory (DFT) and Bader's atoms‐in‐molecules theory, focusing on properties of bond and ring critical points, and their relations to band gap and Pauling electronegativity. We find linear relations of kinetic energy density, electron density, and the gap divided by kinetic energy density at the bond critical points to the difference of Pauling electronegativities of the cation and oxygen anion. At the ring critical points of rutile compounds, we also find that some bond metallicity measures are linearly related to the difference of electronegativities. This study extends our knowledge about the relations between bond critical points, band gap, and electronegativity, but also shows for the first time a quantitative relation between quantities at the ring critical points and global properties of the compounds.     

Systematic investigation on the geometric dependence of the calculated nuclear magnetic shielding constants

We perform a systematic examination on the dependence of the calculated nuclear magnetic shielding constants on the chosen geometry for a selective set of density functional methods of B3LYP, PBE0, and OPBE. We find that the OPBE exchange-correlation functional performs remarkably well when either the optimized geometries or the experimental geometries are used. The popular B3LYP and PBE0 functionals have a clear tendency of deshielding, giving shieldings that are usually too low and shifts that are usually too high, at the experimental geometries. Combined with the Hartree-Fock geometries, however, much improved magnetic constants are obtained for B3LYP and PBE0, due to the compensation effect from the systematic underestimation of bond lengths by the Hartree-Fock method.     

Investigation of challenging spin systems using Monte Carlo configuration interaction and the density matrix renormalization group

We investigate if a range of challenging spin systems can be described sufficiently well using Monte Carlo configuration interaction (MCCI) and the density matrix renormalization group (DMRG) in a way that heads toward a more “black box” approach. Experimental results and other computational methods are used for comparison. The gap between the lowest doublet and quartet state of methylidyne (CH) is first considered. We then look at a range of first‐row transition metal monocarbonyls: MCO when M is titanium, vanadium, chromium, or manganese. For these MCO systems we also employ partially spin restricted open‐shell coupled‐cluster (RCCSD). We finally investigate the high‐spin low‐lying states of the iron dimer, its cation and its anion. The multireference character of these molecules is also considered. We find that these systems can be computationally challenging with close low‐lying states and often multireference character. For this more straightforward application and for the basis sets considered, we generally find qualitative agreement between DMRG and MCCI. © 2017 Wiley Periodicals, Inc.     

One‐Dimensional Arsenic Allotropes: Polymerization of Yellow Arsenic Inside Single‐Wall Carbon Nanotubes

The pnictogen nanomaterials, including phosphorene and arsenene, display remarkable electronic and chemical properties. Yet, the structural diversity of these main group elements is still poorly explored. Here we fill single‐wall carbon nanotubes with elemental arsenic from the vapor phase. Using electron microscopy, we find chains of highly reactive As₄ molecules as well as two new one‐dimensional allotropes of arsenic: a single‐stranded zig‐zag chain and a double‐stranded zig‐zag ladder. These linear structures are important intermediates between the gas‐phase clusters of arsenic and the extended sheets of arsenene. Raman spectroscopy indicates weak electronic interaction between the arsenic and the nanotubes which implies that the formation of the new allotropes is driven primarily by the geometry of the confinement. The relative stabilities of the new arsenic structures are estimated computationally. Band‐gap calculations predict that the insulating As₄ chains become semiconducting, once converted to the zig‐zag ladder, and form a fully metallic allotrope of arsenic as the zig‐zag chain.     

Optical modulation of amplified emission in a polyfluorene-diarylethene blend

A novel system for the modulation of amplified emission based on a polyfluorene/diarylethene (namely F8BT/DTP) blend is shown. The high sensitivity of amplified spontaneous emission (ASE) is exploited to achieve efficient emission modulation with a low-intensity control signal. Modulation is then characterized by photoluminescence (PL) lifetime measurements, photocurrent experiments, and density functional theory calculations. This system can also act as a photocurrent switch based on the same principle. This technique may represent a useful tool for fluorescence quenching and sensing as well as find application in organic photonics.     

四唑负离子与脒正离子复合物稳定性的理论研究

采用密度泛函B3LYP/6-31+G(d,p) 方法和自洽反应场极化连续模型(PCM)研究四唑负离子与脒类(乙脒或苄脒)正离子形成的两种形式的复合物(末端和侧端)及甲酸根负离子与脒类(乙脒或苄脒)正离子形成的复合物在气相和二甲亚砜（DMSO）溶剂中的稳定性。在气相中,四唑-乙脒和四唑-苄脒复合物的相互作用能(∆E)末端分别比侧端的大3.56和3.72 kJ/mol,表明末端复合物稍占优势。甲酸与乙脒和苄脒的复合物的相互作用能(∆E)分别比四唑与乙脒和苄脒的复合物的大59.35和58.99 kJ/mol,表明脒与甲酸形成复合物时相互作用更强。溶剂DMSO的作用使得所有复合物的相互作用能变小,但脒与四唑的相互作用仍比脒与甲酸的弱。前者的结合常数与后者的相比只有1/315（乙脒）和1/218（苄脒）,这与实验结果相一致。     

Density Functional Theory Calculations Revealing Metal‐like Band Structures for Ultrathin Germanium (111) and (211) Surface Layers

To find out if germanium possesses facet‐dependent electrical‐conductivity properties, surface‐state density functional theory (DFT) calculations were performed on one to six layers of germanium (100), (110), (111), and (211) planes. Tunable Ge(100) and Ge(110) planes always present the same semiconducting band structure with a band gap of 0.67 eV expected of bulk germanium. In contrast, one, two, four, and five layers of Ge(111) and Ge(211) plane models show metal‐like band structures with continuous density of states (DOS) throughout the entire band. For three and six layers of Ge(111) and Ge(211) plane models, the normal semiconducting band structure was obtained. The plane layers with metal‐like band structures also show Ge−Ge bond‐length deviations and bond distortions, as well as significantly different 4s and 4p frontier‐orbital electron counts and relative percentages integrated over the valence and conduction bands from those of the semiconducting state. These differences should contribute to strikingly dissimilar band structures. The calculation results suggest the observation of facet‐dependent electrical‐conductivity properties of germanium materials; when making transistors from germanium, the facet effects with shrinking dimensions approaching 3 nm may also need to be considered.     

Modified Guanines as Constituents of Smart Ligands for Nucleic Acid Quadruplexes

Repetitive guanine‐rich nucleic acid sequences play a crucial role in maintaining genome stability and the cell life cycle and represent potential targets for regulatory drugs. Recently, it has been demonstrated that guanine‐based ligands with a porphyrin core can be used as markers of G‐quadruplex assemblies in cell tissues. Herein, model systems of guanine‐based ligands are explored by DFT methods. The energies of formation of modified guanine tetrads and those of modified tetrads stacked on the top of natural guanine tetrads have been calculated. The interaction energy has been decomposed into contributions from hydrogen bonding, stacking, and ion coordination and a twist–rise potential energy scan has been performed to find the individual local minima. Energy decomposition analysis reveals the impact of various substituents (F, Cl, Br, I, Me, NMe₂) on individual energy terms. In addition, cooperative reinforcement in forming the modified and stacked tetrads, as well as the frontier orbitals participating in the hydrogen‐bonding framework involving the HOMO–LUMO gap between the occupied σ_HOMO on the proton‐accepting C=O and =N? groups and unoccupied σ_LUMO on the N?H groups, has been studied. The investigated systems are demonstrated to have a potential in ligand development, mainly due to stacking enhancement compared with natural guanine, which is used as a reference.     

Benchmarking approximate density functional theory for s/d excitation energies in 3d transition metal cations

Holthausen has recently provided a comprehensive study of density functional theory for calculating the s/d excitation energies of the 3d transition metal cations. This study did not include the effects of scalar relativistic effects, and we show here that the inclusion of scalar relativistic effects significantly alters the conclusions of the study. We find, contrary to the previous study, that local functionals are more accurate for the excitation energies of 3d transition method cations than hybrid functionals. The most accurate functionals, of the 38 tested, are SLYP, PBE, BP86, PBELYP, and PW91.     

Excited-state properties and environmental effects for protonated schiff bases: a theoretical study.

Complete active space self-consistent field (CASSCF), multireference configuration interaction (MRCI), density functional theory (DFT), time dependent DFT (TDDFT) and the singles and doubles coupled-cluster (CC2) methodologies have been used to study the ground state and excited states of protonated and neutral Schiff bases (PSB and SB) as models for the retinal chromophore. Systems with two to four conjugated double bonds are investigated. Geometry relaxation effects are studied in the excited pipi* state using the aforementioned methods. Taking the MRCI results as reference we find that CASSCF results are quite reliable even though overshooting of geometry changes is observed. TDDFT does not reproduce bond alternation well in the pipi* state. CC2 takes an intermediate position. Environmental effects due to solvent or protein surroundings have been studied in the excited states of the PSBs and SBs using a water molecule and solvated formate as model cases. Particular emphasis is given to the proton transfer process from the PSB to its solvent partner in the excited state. It is found that its feasibility is significantly enhanced in the excited state as compared to the ground state, which means that a proton transfer could be initiated already at an early step in the photodynamics of PSBs.     

Parameterization of peptide 13C carbonyl chemical shielding anisotropy in molecular dynamics simulations.

NMR chemical shielding anisotropy (CSA) relaxation is an important tool in the study of dynamical processes in proteins and nucleic acids in solution. Herein, we investigate how dynamical variations in local geometry affect the chemical shielding anisotropy relaxation of the carbonyl carbon nucleus, using the following protocol: 1) Using density functional theory, the carbonyl (13)C' CSA is computed for 103 conformations of the model peptide group N-methylacetamide (NMA). 2) The variations in computed (13)C' CSA parameters are fitted against quadratic hypersurfaces containing cross terms between the variables. 3) The predictive quality of the CSA hypersurfaces is validated by comparing the predicted and de novo calculated (13)C' CSAs for 20 molecular dynamics snapshots. 4) The CSA fluctuations and their autocorrelation and cross correlation functions due to bond-length and bond-angle distortions are predicted for a chemistry Harvard molecular mechanics (CHARMM) molecular dynamics trajectory of Ca(2+)-saturated calmodulin and GB3 from the hypersurfaces, as well as for a molecular dynamics (MD) simulation of an NMA trimer using a quantum mechanically correct forcefield. We find that the fluctuations can be represented by a 0.93 scaling factor of the CSA tensor for both R(1) and R(2) relaxations for residues in helix, coil, and sheet alike. This result is important, as it establishes that (13)C' relaxation is a valid tool for measurement of interesting dynamical events in proteins.     

Experimental Electron Densities and Chemical Bonding

Some recent results of charge density analysis by X-ray and neutron diffraction are discussed. Problems that have been studied in a number of laboratories include the nature of single, double, and triple bonds, lone-pair hybridization, bonding in some metals, alloys, and organometallic compounds, and the derivation of physical properties from X-ray diffraction densities. At the present stage of development of methods studies of series of related compounds are feasible and expected to find widespread application.     

Oxidation Effects on the Electronic Properties of Hydrogenated Silicon Clusters

The geometric and electronic properties of some hydrogenated silicon clusters in the presence of oxygen on the surface have been investigated.The density functional theory with generalized gradient approximation functional was applied in our calculations.By calculating the total energy,the double bond Si=O is shown to be more stable than the bridge bond Si-O-Si for large size oxidized clusters.The results of Mulliken population analysis indicate that a so-called passivation effect is enhanced by oxidization effects.From the energy band structures and density of states,we find that some localized states are induced by the p-orbital of O atom mainly and reduce the energy gaps substantially.