Tuning nonlinear optical and optoelectronic properties of vinyl coupled triazene chromophores: a density functional theory and time-dependent density functional theory investigation

Triazenes are a unique class of polyazo compounds containing three consecutive nitrogen atoms in an acyclic arrangement and are promising NLO candidates. In the present work, a series of 15 donor-π-acceptor type vinyl coupled triazene derivatives (VCTDs) with different acceptors (-NO(2), -CN, and -COOH) have been designed, and their structure, nonlinear response, and optoelectronic properties have been studied using density functional theory and time-dependent density functional theory methods. B3LYP/6-311g(d,p) optimized geometries of the designed candidates show delocalization from the acceptor to donor through a π-bridge. Molecular orbital composition analysis reveals that HOMO is stabilized by the π-bridge, whereas acceptors play a major role in the stabilization of LUMO. Among the three acceptors, nitro derivatives are found to be efficient NLO candidates, and tri- and di-substituted cyano and carboxylic acid derivatives also show reasonably good NLO response. The effect of solvation on these properties has been studied using PCM calculations. From TDDFT calculations, the computed absorption spectra of these candidates lie in the range of 350-480 nm in the gas phase and have positive solvatochromism. The ground-state stabilization interactions are accounted from NBO calculations. In an effort to substantiate the thermal stability of the designed candidates, computations have been done to identify the weak interactions in the systems through NCI and AIM analysis. In summary, 10 out of 15 designed candidates are found to have excellent NLO and optoelectronic properties.     

Molecular properties via a subsystem density functional theory formulation: a common framework for electronic embedding

In this article, we present a consistent derivation of a density functional theory (DFT) based embedding method which encompasses wave-function theory-in-DFT (WFT-in-DFT) and the DFT-based subsystem formulation of response theory (DFT-in-DFT) by Neugebauer [J. Neugebauer, J. Chem. Phys. 131, 084104 (2009)] as special cases. This formulation, which is based on the time-averaged quasi-energy formalism, makes use of the variation Lagrangian techniques to allow the use of non-variational (in particular: coupled cluster) wave-function-based methods. We show how, in the time-independent limit, we naturally obtain expressions for the ground-state DFT-in-DFT and WFT-in-DFT embedding via a local potential. We furthermore provide working equations for the special case in which coupled cluster theory is used to obtain the density and excitation energies of the active subsystem. A sample application is given to demonstrate the method.     

Theoretical investigation on electronic structure and stability of some inverted cucurbiturils by density functional theory

Structure and stability of diastereoisomers of cucurbit[n]urils (CB[n = 5–10]), the inverted CB[n]s, were investigated by density functional theory (DFT) computations. All the inverted CB[n]s were less stable than their normal CB[n]s and the mono-inverted ones with one inverted glycoluril unit in their structures were more stable than their doubly-inverted isomers. Relative change in dipole moments and molecular electrostatic potentials (MEP) were discussed with the deformation in geometric structure and charge distribution of the normal and inverted CB[n]s.     

Influence of the central metal on the structural,electronic, and spectroscopic properties of naphthalocyanine complexes: density functional theory calculations

A series of metal naphthalocyanine complexes (M = TiO²⁺, Fe²⁺, Co²⁺, Ni²⁺, Cu²⁺, Zn²⁺, Ru²⁺) have been investigated using density functional theory (DFT) and time-dependent DFT methods in vacuo and in the solvent dimethylsulfoxide in order to evaluate the influence of the different metal atoms on the geometries and optical properties of their complexes. The optimized geometries for the complexes without an axial ligand exhibit planar conformations. Most of the absorption bands of the metal complexes are blue-shifted compared to those of the metal-free naphthalocyanine, both in vacuo and in the solvent. The various transition metals could gradually tune the electronic and spectroscopic properties of their naphthalocyanine complexes, which may provide valuable information for tuning the properties of naphthalocyanine complexes for various applications.     

Crystal structure and density functional theory study on structural properties and energies of a isonicotinohydrazide compound

An X-ray and a theoretical study of the structure of the isoniazid derivative N'-(4-dimethylaminobenzylidene)-isonicotinohydrazide monohydrate (1) are reported. In this work, we will report a combined experimental and theoretical study on the molecular structure, vibrational spectra and energies of N'-(4-dimethylaminobenzylidene)-isonicotinohydrazide monohydrate. The calculated parameters are in good agreement with the corresponding X-ray diffraction values. The FTIR spectrum in the range of 400-4000 cm-1 of N'-(4-dimethylaminobenzylidene)-isonicotinohydrazide monohydrate has been recorded. The molecular geometry and vibrational frequencies and energies in the ground state are calculated by using the DFT (B3LYP, PBE1PBE) methods with 6-311G** basis sets. The calculated HOMO and LUMO energies also confirm that charge transfer occurs within the molecule. The geometries and normal modes of vibrations obtained from B3LYP/PBE1PBE/6-311G** calculations are in good agreement with the experimentally observed data.     

Substituent effects on benzdiyne: a density functional theory study

Density functional theory (DFT) studies were performed to investigate the effect of substituents on the properties of benzdiyne derivatives. Twelve substituted benzdiynes-C(6)X(2), where X = F, Cl, Br, Me, CF(3), CN, OH, NO(2), NH(2), OMe, NMe(2), and Ph-were considered along with the unsubstituted 1,4-benzdiyne. The structures, vibrational frequencies, and IR intensities of these benzdiynes were studied with a popular three-parameter hybrid density functional (B3LYP) combined with the split-valence 6-31G(d) basis set and Dunning's correlation-consistent polarized triple-zeta (cc-pVTZ) basis set. The relative stabilities of the substituted benzdiynes were studied with the help of reaction energies of isodesmic reactions, which showed that the electron-withdrawing groups destabilized the benzdiynes more than they did the corresponding benzenes, whereas the electron-donating groups stabilized the benzdiynes more than they did their benzene counterparts. Correlation analyses revealed that field/inductive effects played a more important role than did resonance effects. The changes in atomic charges and spin populations due to the substituents were also studied. The asymmetric nu(Ctbd1;C) stretching modes obtained were close to the 1500-cm(-)(1) mark. Reinvestigation of the experimental results supported these results; a weak IR band at 1486 cm(-)(1) was assigned to this asymmetric stretching mode in C(6)(CF(3))(2) F. Some other benzdiynes also had large IR intensity values for their asymmetric nu(Ctbd1;C) vibrational modes due to the coupling with other vibrational modes. Heats of formation for the substituted benzdiynes were obtained from the reaction energies calculated at the B3LYP/cc-pVTZ level of theory.     

Characterization of the high-pressure structural transition and thermodynamic properties in sodium chloride: a computational investigation on the basis of the density functional theory

Using first-principles calculations, the elastic constants, the thermodynamic properties, and the structural phase transition between the B1 (rocksalt) and the B2 (cesium chloride) phases of NaCl are investigated by means of the pseudopotential plane-waves method. The calculations are performed within the generalized gradient approximation to density functional theory with the Perdew-Burke-Ernzerhof exchange-correlation functional. On the basis of the third-order Birch-Murnaghan equation of states, the transition pressure Pt between the B1 phase and the B2 phase of NaCl is determined. The calculated values are generally speaking in good agreement with experiments and with similar theoretical calculations. From the theoretical calculations, the shear modulus, Young's modulus, rigidity modulus, and Poisson's ratio of NaCl are derived. According to the quasi-harmonic Debye model, we estimated the Debye temperature of NaCl from the average sound velocity. Moreover, the pressure derivatives of elastic constants, partial differentialC11/partial differentialP, partial differentialC12/partial differentialP, partial differentialC44/partial differentialP, partial differentialS11/partial differential P, partial differentialS12/partial differentialP, and partial differentialS44/partial differentialP, for NaCl crystal are investigated for the first time. This is a quantitative theoretical prediction of the elastic and thermodynamic properties of NaCl, and it still awaits experimental confirmation.     

Structural, electronic, and bonding properties of zeolite Sn-beta: a periodic density functional theory study

The structural, electronic, and the bonding properties of the zeolite Sn-beta (Sn-BEA) have been investigated by using the periodic density functional theory. Each of the nine different T-sites in BEA were substituted by Sn atoms and all the nine geometries were completely optimized by using the plane-wave basis set in conjunction with the ultra-soft pseudopotential. On the basis of the structural and the electronic properties, it has been demonstrated that the substitution of Sn atoms in the BEA framework is an endothermic process and hence the incorporation of Sn in the BEA is limited. The lowest unoccupied molecular orbitals (LUMO) energies have been used to characterize the Lewis acidity of each T-site. On the basis of the relative cohesive energy and the LUMO energy, the T2 site is shown to be the most favorable site for the substitution Sn atoms in the BEA framework.     

Differentiability of Lieb functional in electronic density functional theory

A solid understanding of the Lieb functional F_L is important because of its centrality in the foundations of electronic density functional theory. A basic question is whether directional derivatives of F_L at an ensemble‐V‐representable density are given by (minus) the potential. A widely accepted purported proof that F_L is Gâteaux differentiable at EV‐representable densities would say, “yes.” But that proof is fallacious, as shown here. F_L is not Gâteaux differentiable in the normal sense, nor is it continuous. By means of a constructive approach, however, we are able to show that the derivative of F_L at an EV‐representable density ρ₀ in the direction of ρ₁ is given by the potential if ρ₀ and ρ₁ are everywhere strictly greater than zero, and they and the ground state wave function have square integrable derivatives through second order. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Cyclic and ladder hydrogen bonded cyanamide oligomers: a density functional theory and many-body analysis approach

This work reports hydrogen bonding interaction in cyclic and ladder oligomers using density functional theory method. Many-body analysis technique has been used to study the nature of interactions between different molecules and their contribution to the binding energy of a respective hydrogen bonded oligomers. Hydrogen bonds in cyclic trimer to pentamer are stronger than those in corresponding ladder structures. Cyanamide monomer shows the lowest energy at B3LYP/aug-cc-pvdz level among different methods used here with the same basis set. The geometrical parameters for cyanamide monomer obtained at B3LYP/aug-cc-pvdz level are in excellent agreement with the experimental determinations. Cyclic structures are more stable than the ladder. In cyclic oligomers not only total two-body energies, but higher body energies also contribute significantly to the binding energy of a respective complex whereas in ladder, only total two-body energies contribute significantly and higher-body energies are almost negligible for cyanamide trimer to pentamer.     

Small fullerenes with BN belts: a density functional theory investigation

Density functional calculations have been carried out on a series of BCN hybrid fullerenes with certain substitution patterns in comparison with their parent compounds Cn (n = 30, 32, 36, 38, 40, 44, 48, 50, 52). The substitutional structures, energy gaps between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, ionization potentials, electron affinities, as well as molecular electrostatic potentials have been systematically investigated. The following important points of BCN hybrid fullerenes are stressed: The present studied fullerenes, comprising tubular "belt" and polar "cap", could be divided into three types of structure; each has different indexes of tubular structure and terminal caps. The properties of BCN hybrid fullerenes depend on the type of "tubular belt + polar cap" structures, especially, the HOMO and LUMO characters and MEPs of BCN fullerene are strongly governed by their structure types.     

Geometries, stabilities, and electronic properties of small anion Mg-doped gold clusters: a density functional theory study

First-principle density functional theory is used for studying the anion gold clusters doped with magnesium atom. By performing geometry optimizations, the equilibrium geometries, relative stabilities, and electronic and magnetic properties of [Au(n)Mg]? (n = 1-8) clusters have been investigated systematically in comparison with pure gold clusters. The results show that doping with a single Mg atom dramatically affects the geometries of the ground-state Au(n+1)? clusters for n = 2-7. Here, the relative stabilities are investigated in terms of the calculated fragmentation energies, second-order difference of energies, and highest occupied?lowest unoccupied molecular orbital energy gaps, manifesting that the ground-state [Au(n)Mg]? and Au(n+1)? clusters with odd-number gold atoms have a higher relative stability. In particular, it should be noted that the [Au?Mg]? cluster has the most enhanced chemical stability. The natural population analysis reveals that the charges in [Au(n)Mg]? (n = 2-8) clusters transfer from the Mg atom to the Au frames. In addition, the total magnetic moments of [Au(n)Mg]? clusters exhibit an odd-even oscillation as a function of cluster size, and the magnetic effects mainly come from the Au atoms.     

Excitonic effects in a time-dependent density functional theory

Excited state properties of one-dimensional molecular materials are dominated by many-body interactions resulting in strongly bound confined excitons. These effects cannot be neglected or treated as a small perturbation and should be appropriately accounted for by electronic structure methodologies. We use adiabatic time-dependent density functional theory to investigate the electronic structure of one-dimensional organic semiconductors, conjugated polymers. Various commonly used functionals are applied to calculate the lowest singlet and triplet state energies and oscillator strengths of the poly(phenylenevinylene) and ladder-type (poly)(para-phenylene) oligomers. Local density approximations and gradient-corrected functionals cannot describe bound excitonic states due to lack of an effective attractive Coulomb interaction between photoexcited electrons and holes. In contrast, hybrid density functionals, which include long-range nonlocal and nonadiabatic corrections in a form of a fraction of Hartree-Fock exchange, are able to reproduce the excitonic effects. The resulting finite exciton sizes are strongly dependent on the amount of the orbital exchange included in the functional.     

Growth pattern and electronic properties of acetonitrile clusters: a density functional study

We report a systematic theoretical study on the growth pattern and electronic properties of acetonitrile clusters [(CH(3)CN)(n) (n = 1, 9, 12)] using density functional approach at the B3LYP6-31++G(d,p) level. Although we have considered a large number of configurations for each cluster, the stability of the lowest energy isomer was verified from the Hessian calculation. It is found that the lowest energy isomer of the dimer adopts an antiparallel configuration. For trimer and tetramer, cyclic ring structures were found to be favored over the dipole stabilized structure. In general, it is found that the intermolecular CH...N interactions play a significant role in the stabilization of the cyclic layered geometry of acetonitrile clusters. A critical comparison between trimer and tetramer clusters suggests that the three member cyclic ring is more stable than four member rings. The growth motif for larger clusters (n = 5-9, 12) follows a layered pattern consisting of three or four membered rings, which, in fact, is used as the building block. Based on the stability analysis, it is found that clusters with an even number of molecular entities are more stable than the odd clusters, except trimer and nonamer. The exceptional stability of these two clusters is attributed to the formation of trimembered cyclic rings, which have been found to form the building blocks for larger clusters.     

Octacoordinate carbons encaged inside carborane clusters: a density functional theory investigation

This work focuses on the computational design and characterization of a novel series of endohedral carborane clusters containing octacoordinate carbon centers. The structural and bonding features and the thermodynamic and kinetic stabilities are discussed extensively based on density functional theory calculations. These nonclassical carboranes are fascinating in structure not only for the octacoordinate carbon center but also for the surrounding carbon and boron ligands with inverted bonding configuration. These endohedral carboranes are higher in energy than the corresponding exohedral isomers due to the high strain in the system. A new stability rule based on the donor-acceptor model is proposed to predict the stability ordering for these carborane isomers. In addition, some of these octacoordinate carboranes might have relatively high kinetic stabilities, which is rather hopeful for the experimental syntheses.     

Structure, spectroscopy, and reactivity properties of porphyrin pincers: a conceptual density functional theory and time-dependent density functional theory study

Porphyrin and pincer complexes are both important categories of compounds in biological and catalytic systems. The idea to combine them is computationally investigated in this work. By employment of density functional theory (DFT), conceptual DFT, and time-dependent DFT approaches, structure, spectroscopy, and reactivity properties of porphyrin pincers are systematically studied for a selection of divalent metal ions. We found that the porphyrin pincers are structurally and spectroscopically different from their precursors and are more reactive in electrophilic and nucleophilic reactions. A few quantitative linear/exponential relationships have been discovered between bonding interactions, charge distributions, and DFT chemical reactivity indices. These results are implicative in chemical modification of hemoproteins and understanding chemical reactivity in heme-containing and other biologically important complexes and cofactors.     

Calculation of the electronic and optical properties of indium tin oxide by density functional theory

Density functional theory (DFT) was used to calculate the bulk electronic and optical properties of indium tin oxide (ITO). The ITO model was constructed replacing indium atoms with tin atoms in the cubic unit cell of indium oxide. To allow more possibilities for tin atom substitution than afforded by the forty-atom primitive cell of indium oxide all eighty atoms of the unit cell were included in the stoichiometry (In_32−xSn_xO₄₈) using periodic boundary conditions. A number of properties of ITO were calculated including the optical band gap, charge carrier density and plasma frequency. The dependence of the electronic and optical properties of ITO on a variety of parameters such as the tin content, cubic lattice parameter and the distance between adjacent tin atoms was investigated. The electronic and optical properties agreed well with experimental data and allowed insight into the origin of the electronic and optical properties of ITO.     

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.     

Determination of structures, stabilities, and electronic properties for bimetallic cesium-doped gold clusters: a density functional theory study

The equilibrium geometric structures, stabilities, and electronic properties of bimetallic Au(n)Cs (n = 1-10) and pure gold Au(n) (n ≤ 11) clusters have been systematically investigated by using density functional theory with meta-generalized gradient approximation. The optimized geometries show that one Au atom capped on Au(n-1)Cs structures and Cs atom capped Au(n) structures for different sized Au(n)Cs (n = 1-10) clusters are two dominant growth patterns. Theoretical calculated results indicate that the most stable isomers have three-dimensional structures at n = 4 and 6-10. Averaged atomic binding energies, fragmentation energies, and second-order difference of energies exhibit a pronounced even-odd alternations phenomenon. The same even-odd alternations are found in the highest occupied-lowest unoccupied molecular orbital gaps, vertical ionization potential, vertical electron affinity, and hardnesses. In addition, it is found that the charge in corresponding Au(n)Cs clusters transfers from the Cs atom to the Au(n) host in the range of 0.851-1.036 electrons.