Density functional theory of polymer-polymer phase separation behavior

Polyatomic density functional theory is applied to a binary polymer blend. The polymer reference interaction site model (PRISM) liquid state theory provides the homogeneous state correlation functions necessary for the application of density functional theory. An effective chi parameter can be recognized from the density functional expression; however, the phase separation criteria does not depend solely upon the chi parameter, rather it depends upon various combinations of the species-dependent direct correlation functions of the blend. The Flory-Huggins chi parameter along with the associated phase diagram is obtained when the monomer volumes of the blend species are equal and for a range of monomer-monomer attractive interactions. Calculations are performed both with and without the assumption of incompressibility. The density functional theory along with the PRISM determined “input” predict that an isotopic polymer blend shows an upper critical solution temperature (UCST) phenomena. © 1995 John Wiley & Sons, Inc.     

量子相空间中微扰理论

在Torres-Vega和Frederick所提出来的量子相空间理论(简作TF量子相空间理论)的框架下,研究了量子相空间表象下的微扰理论,得出了量子相空间理论框架下的非简并微扰理论(perturbation theory),得到了体系存在微扰情况下的能量和波函数,并且以一维电解质在外加电场中的极化率为例,在量子相空间表象下,对量子体系的状态进行了分析.以期探索出将量子相空间理论和各种近似方法结合起来处理真实体系的问题.     

Density functional theory study of the structure and energetics of negatively charged oligopyrroles

First‐principles calculations are used to investigate the electronic properties of negatively charged n‐pyrrole oligomers with n = 2–18. Chains of neutral oligomers are bent, whereas the negatively charged oligomers become almost planar due to accumulation of negative charge at the end monomers. Isomers of short oligomers (n < 6) display negative electron affinity although the corresponding anions are energetically stable. For longer oligomers with n ≥ 6, the electron affinity is small and positive, slowly increasing with oligopyrrole length. Doping of 12‐pyrrole with lithium atoms shows that negative oxidation states are possible due to electron transfer from dopant to oligomer at locations close to dopant. These 12‐pyrrole regions support extra negative charge and exhibit a local structural change from benzenoid to quinoid structure in the C? C backbone conjugation. Comparison between neutral and doped polypyrrole (PPy) indicates that doped polymers displays a substantial depletion of the band gap energy and the appearance of dopant‐based bands in the gap for a 50% per monomer doping level. It is predicted that Li‐doped PPy is not metallic. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first-principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.     

Gas-phase acidity of D-glucose. A density functional theory study

The gas-phase acidity of D-glucopyranose was studied by means of B3LYP calculations combined with 6-31G(d,p) or 6-31+G(d,p) standard basis sets. For each anomer, deprotonation of the various primary and secondary hydroxyl groups was considered. As in solution, the anomeric hydroxyl is found to be the most acidic for both anomers, but only when the 6-31+G(d,p) basis set is used for geometry optimization. Deprotonation of the anomeric hydroxyl induces an important C(1)--O endocyclic bond elongation and subsequently promotes an energetically favored ring-opening process as attested by the very small calculated activation barriers. The results also suggest that interconversion between the various deprotonated alpha- and beta-anomers may easily occur under slightly energetic conditions. B3LYP/6-311+G(2df,2p) calculations led to the an absolute gas-phase acidity of deltaacidGo(298)(alpha-D-glucose) = 1398 kJ mol(-1). This estimate matches well the only experimental value available to date. Finally, this study again confirms that the use of diffuse functions on heavy atoms is necessary to describe anionic systems properly and to achieve good relative and absolute gas-phase acidities.     

Density functional theory study on electron and hole transport properties of organic pentacene derivatives with electron-withdrawing substituent

Attaching electron-withdrawing substituent to organic conjugated molecules is considered as an effective method to produce n-type and ambipolar transport materials. In this work, we use density functional theory calculations to investigate the electron and hole transport properties of pentacene (PENT) derivatives after substituent and simulate the angular resolution anisotropic mobility for both electron and hole transport. Our results show that adding electron-withdrawing substituents can lower the energy level of lowest unoccupied molecular orbital (LUMO) and increase electron affinity, which are beneficial to the electron injection and ambient stability of the material. Also the LUMO electronic couplings for electron transport in these pentacene derivatives can achieve up to a hundred meV which promises good electron transport mobility, although adding electron-withdrawing groups will introduce the increase of electron transfer reorganization energy. The final results of our angular resolution anisotropic mobility simulations show that the electron mobility of these pentacene derivatives can get to several cm(2) V(-1) s(-1), but it is important to control the orientation of the organic material relative to the device channel to obtain the highest electron mobility. Our investigation provide detailed information to assist in the design of n-type and ambipolar organic electronic materials with high mobility performance.     

Density functional theory study of tautomerization of 2‐aminothiazole in the gas phase and in solution

The amino/imino tautomeric equilibrium in the isolated, mono‐, di‐, and trihydrate forms and dimer of 2‐aminothiazole, and the effects of hydration or self‐assistance on the transition state structures corresponding to proton transfer from the amino to imino form, have been investigated by the B3LYP method in conjunction with 6‐31+G(d,p) and 6‐311+G(3df,2p) basis sets in the gas phase and in solution. The amino form has been found to be the predominant tautomer. The tautomeric barrier heights for water‐ and self‐assisted tautomerization reactions are significantly lower than that from the amino to imino form by the intramolecular proton transfer, showing the catalytic effect of water molecule(s) and the important role of 2‐aminothiazole itself for intermolecular proton transfer. Comparison between the tautomeric barriers demonstrates that the self‐association tautomerization through the dimerization is the most favorable pathway. Bulk solvent effects have been taken into account using the polarizable continuum model (PCM) of water and CCl₄. The polar medium is favorable for the population of the imino form. The amino/imino equilibrium is also analyzed using the aromaticity index nucleus‐independent chemical shift (NICS); the NICS values for the amino form (about ?10 ppm) are more negative than the imino species (about ?8 ppm), showing that the amino form is more stable. The time‐dependent density functional theory (TDDFT) calculations of electronic absorption spectra suggest that the λ_max of dimer is 255 nm. The oscillator strength of the imino forms is less than the amino form, and increases with the polarity of the solvents. All calculations for the tautomerization of 2‐aminothiazole are in reasonable line with the available experiments. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Analytical gradients

A full implementation of analytical energy gradients for molecular and periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian‐type orbitals as basis functions. Its key component is a combination of density fitting (DF) approximation and continuous fast multipole method (CFMM) that allows for an efficient calculation of the Coulomb energy gradient. For exchange‐correlation part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097) is extended to energy gradients. Computational efficiency and asymptotic O(N) scaling behavior of the implementation is demonstrated for various molecular and periodic model systems, with the largest unit cell of hematite containing 640 atoms and 19,072 basis functions. The overall computational effort of energy gradient is comparable to that of the Kohn–Sham matrix formation. © 2016 Wiley Periodicals, Inc.     

磷化镍晶体结构的密度泛函理论研究

运用密度泛函理论系统地研究了已知的磷化镍晶体如Ni3P、 Ni12P5、 Ni2P、 Ni5P4、 NiP、 NiP2 和 NiP3等的结构、成键以及相关的热力学稳定性,对于结构简易化合物（Ni2P和NiP3）的弹性行为进行了预测。这些数据有助于更好的理解磷化镍的催化活性。     

黄嘌呤及其互变异构体的密度泛函理论研究

嘌呤碱及其衍生物在生物系统中起重要作用。对人具有兴奋和利尿作用的茶碱和咖啡碱就是黄嘌呤的甲基衍生物。黄嘌呤存在多种互变异构体,从理论计算的角度研究这些互变异构体的几何结构、电子结构及相对稳定性,进一步研究溶剂对其结构和性质的影响是有意义旧。密度泛函理论方法既考虑了电子相关,又较其它CI（组态相互作用）或MPn（n级微扰）方法节省机时。其计算结果较好,被广泛应用于研究各种化合物。本文采用密度泛函B3LYP／6—311G方法对14种黄嘌呤可能的互变异构体（见图1）,分别在气相和水相中进行几何构型全自由度优化和能量计算,讨论了异构体的相对稳定性,水的溶剂化作用对异构体的能量、几何构型、电荷分布和偶极矩的影响,探讨了溶剂极性对异构体的能量和偶极矩的影响。     

Density functional theory rationalization of the substituent effects in trifluoromethyl-pyridinol derivatives

The influence of α-substitution in the structure, bonding and thermochemical properties of trifluoromethyl-pyridinol derivatives and their protonated counterparts has been studied by means of density functional theory. The geometries of the neutral and protonated species were optimized at the B3-LYP/6-311G(d,p) level of theory. Final energies were obtained through single point B3-LYP/6-311+G(3df,2p) calculations.The relative orientation of the different substituents within the heterocycle ring favours the formation of unexpected intramolecular hydrogen bonds (IHB), which have been characterized by means of the Atoms in Molecules theory of Bader. Although weak, these IHB are of great importance for understanding the gas phase structure and the thermodynamical properties of these compounds. Surprisingly, most of the substituted investigated pyridinols present proton affinities below or close to that calculated for the unsubstituted pyridine molecule. Only pyridinols bearing strong σ or π donor activating groups show proton affinities greater than that of pyridine.     

Exchange and correlation energy in density functional theory

Several different versions of density functional theory (DFT) that satisfy Hohenberg–Kohn theorems are characterized by different definitions of a reference or model state determined by an N‐electron ground state. A common formalism is developed in which exact Kohn–Sham equations are derived for standard Kohn–Sham theory, for reference‐state density functional theory, and for unrestricted Hartree–Fock (UHF) theory considered as an exactly soluble model Hohenberg–Kohn theory. A natural definition of exchange and correlation energy functionals is shown to be valid for all such theories. An easily computed necessary condition for the locality of exchange and correlation potentials is derived. While it is shown that in the UHF model of DFT the optimized effective potential (OEP) exchange satisfies this condition by construction, the derivation shows that this condition is not, in general, sufficient to define an exact local exchange potential. It serves as a test to eliminate proposed local potentials that are not exact for ground states. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 521–525, 2000     

Gas‐phase pyrolysis mechanisms of 3‐anilino‐1‐propanol: Density functional theory study

The gas‐phase pyrolytic decomposition mechanisms of 3‐anilino‐1‐propanol with the products of aniline, ethylene, and formaldehyde or N‐methyl aniline and aldehyde were studied by density functional theory. The geometries of the reactant, transition states, and intermediates were optimized at the B3LYP/6‐31G (d, p) level. Vibration analysis was carried out to confirm the transition state structures, and the intrinsic reaction coordinate method was performed to search the minimum energy path. Four possible reaction channels are shown, including two concerted reactions of direct pyrolytic decomposition and two indirect channels in which the reactant first becomes a ring‐like intermediate, followed by concerted pyrogenation. One of the concerted reactions in the direct pyrolytic decomposition has the lowest activation barrier among all the four channels, and so, it occurs more often than others. The results appear to be consistent with the experimental outcomes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

白藜芦醇清除自由基活性的密度泛函理论研究

采用密度泛函理论(DFT)方法,从静态与动态两大方面分析了白藜芦醇分子酚羟基在不同溶剂中清除自由基活性的能力大小.通过白藜芦醇的结构参数、前线轨道理论、3种抽氢反应机制等方面分析了分子活性位与其性质的关系.探讨了白藜芦醇分子不同位置酚羟基清除·OH和·OOH的抗氧化机理,得到了该分子与·OOH发生抽氢反应时的过渡态结构.计算结果表明,在任何溶剂中白藜芦醇分子C(4')位上酚羟基活性最高,发生抽氢反应时反应热最小,是高活性位点.     

First‐order correlation‐kinetic contribution to Kohn–Sham exchange charge density function in atoms,using quantal density functional theory approach

Using the static exchange‐correlation charge density concept, the total integrated exchange‐charge density function is calculated within the nonrelativistic spin‐restricted exchange‐only (i) optimized effective potential model, and (ii) nonvariational local potential derived from the exchange‐only work potential within the quantal density functional theory, for the ground‐state isoelectronic series: Ga⁺, Zn, Cu^?; In⁺, Cd, Ag^?; and Tl⁺, Hg, Au^?. The difference between the exchange charge density function derived from these potentials is employed to evaluate the first‐order correlation‐kinetic contribution to the integrated exchange charge density. This contribution is found to be important for both the intra‐ and inter‐shell regions. Screening effects on the contribution due to the nd¹⁰ (n = 3–5) subshells are discussed through comparisons with similar calculations on Ca, Sr, and Ba, wherein nd¹⁰ electrons are absent. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Dimetallic sulfide endohedral metallofullerene Sc2S@C76: Density functional theory characterization

In terms of density functional theory combined with statistic mechanics computations, we investigated a dimetallic sulfide endohedral fullerene Sc₂S@C₇₆ which has been synthesized without any characterization in experiments. Our theoretical study reveals that Sc₂S@T_d(19151)‐C₇₆ which satisfies the isolated‐pentagon rule (IPR) possesses the lowest energy, followed by three non‐IPR structures (Sc₂S@C_2v(19138)‐C₇₆, Sc₂S@C_s (17490)‐C₇₆, and Sc₂S@C₁(17459)‐C₇₆). To clarify the relative stabilities of those isomers at high temperatures, enthalpy–entropy interplay has been taken into consideration. Calculation results indicate that three species Sc₂S@T_d(19151)‐C₇₆, Sc₂S@C_2v(19138)‐C₇₆, and Sc₂S@C₁(17459)‐C₇₆ have noticeable molar fractions at the fullerene‐formation temperature region (500–3000K), and the Sc₂S@C₁(17459)‐C₇₆ with one pentagon pair becomes the most predominant isomer above 1800 K, suggesting that the unexpected non‐IPR structure is thermodynamically favorable at elevated temperatures. In addition, the structural characteristics, electron features, UV‐vis‐NIR adsorptions, and ¹³C NMR spectra of those three stable structures are introduced to assist experimental identification and characterization in future. © 2014 Wiley Periodicals, Inc.     

Density functional calculations on triply excited states of lithium isoelectronic sequence

Triply excited states of many-electron atomic systems are characterized by the presence of strong electron correlation, closeness to more than one threshold, and degeneracy with many continua; therefore, they offer unusual challenges to theoretical methodologies. In the present article, we computed with reasonable accuracy all the n=2 intrashell triply excited states (2s²2p ²P; 2s2p^{2 2}D, ⁴P, ²P, ²S; and 2p^{3 2}D, ²P, ⁴S) of three-electron atomic systems (Z=2, 3, 4, 6, 8, 10) by using a density functional formalism developed recently in our laboratory, based on the nonvariational Harbola–Sahni exchange potential in conjunction with a parametrized local Wigner and Lee–Yang–Parr correlation potentials. Nonrelativistic energies and densities are obtained by solving a Kohn–Sham-type differential equation. The calculated results are compared with available experimental and other theoretical data. The 2p³(⁴S)→1s2p²(⁴P) transition wavelength for the isoelectronic series is also computed. The overall good agreement of our results with the literature data indicates the reliability of the present density functional methodology for excited states of many-electron systems. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65 : 317–332, 1997     

Orbital functional theory of linear response and excitation

Orbital functional theory (OFT) is based on a rule that determines a single‐determinant reference state Φ for any exact N‐electron eigenstate Ψ. An OFT model postulates an explicit correlation energy functional E_c of occupied orbital functions {?_i} and occupation numbers {n_i}. The orbital Euler–Lagrange equations are analogous to Kohn–Sham equations, but do not in general contain local potential functions. Time‐dependent Hartree–Fock theory is generalized in OFT to a formally exact linear response theory that includes electronic correlation. In the exchange‐only limit, the theory reduces to the random‐phase approximation of many‐body theory. The formalism determines excitation energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Density functional theory study of p‐chloroaniline adsorption on Pd surfaces and clusters

The adsorption mode of aromatic molecules on transition metal surfaces plays a key role in their catalytic transformation. In this study, by means of density functional theory calculations, we systematically investigate the adsorption of p‐chloroaniline on a series of Pd surfaces, including stepped surfaces, flat surfaces, and clusters. The adsorption energies of p‐chloroaniline on these substrates [Pd(221), Pd(211), Pd(111), Pd(100), Pd₁₃‐icosahedral, Pd₁₃‐cubo‐octahedron, Pd₅₅] are ?1.90, ?2.13, ?1.70, ?2.11, ?2.53, ?2.65, ?2.23 eV, respectively. Benzene ring is adsorpted on catalyst rather than amine group in p‐chloroaniline molecular. A very good linear relationship is further found between the adsorption energies of p‐chloroaniline and the d‐band center of both Pd surfaces and clusters. The lower of d‐band center of Pd models, the stronger adsorption of p‐chloroaniline on catalysts. In addition, the frontier molecular orbital and density of states analysis explain the adsorption energy sequence: cluster Pd₁₃ > stepped Pd(221) surface > flat Pd(111) surface. © 2014 Wiley Periodicals, Inc.