Theoretical study on the structure and UV–visible spectrum of pyridine–chloranil charge-transfer complex

The ground-state structure of the charge-transfer complex formed by pyridine (Py) as electron donor and chloranil (CA) as acceptor has been studied by full geometry optimization at the MP2 and DFT levels of theory. Binding energies were calculated and counterpoise corrections were used to correct the BSSE. Both MP2 and DFT indicate that the pyridine binds with chloranil to form an inclined T-shape structure, with the pyridine plane perpendicular to the chloranil. The CP and ZPE corrected binding energies were calculated to be 14.21 kJ/mol by PBEPBE/6-31G(d) and 23.21 kJ/mol by MP2/6-31G(d). The charge distribution of the ground state Py–CA complex was evaluated with the natural population analysis, showing a net charge transfer from Py to CA. Analysis of the frontier molecular orbitals reveals a σ–π interaction between CA and Py, and the binding is reinforced by the attraction of the O₇ atom of CA with the H₂₃ atom of Py. TD-DFT calculations have been performed to analyze the UV–visible spectrum of Py–CA complex, revealing both the charge transfer transitions and the weak symmetry-relieved chloranil π–π^* transition in the UV–visible region.     

Ab initio and DFT studies of hydrogen bond interactions in difluoroacetic acid dimer

Density functional theory (DFT) and ab initio calculations were performed for difluoroacetic acid (DFA). Eight theoretically possible conformers were considered and by using conformational analysis only three stable conformers were found. The hydrogen bonding interaction of DFA complex has been investigated using DFT and ab initio methods for cis conformers. Stabilization energies of dimers including basis set superposition error and ZPE were found in the range 8.89–13.08 kcal mol⁻¹. It was found that EFC dimer is slightly more stable. Red shift of O–H bond in the range −226.3 to 505.7 cm⁻¹ predicted for dimers. The natural bond orbital analysis was applied to characterize nature of the interaction.     

Vibrational spectroscopic and density functional theory studies of chloranil-imidazole interaction

The present work reports vibrational spectra and density functional theory calculations for chloranil, imidazole and their complexes. The experimentally observed infrared and Raman bands have been assigned with the help of calculated vibrational frequencies and potential energy distribution analysis. Some bands of chloranil and imidazole have been found to shift on the complex formation due to partial electronic charge transfer from imidazole to chloranil. The charge transfer between these molecules is also corroborated by the electronic absorption spectroscopy and calculations. The theoretical values of the interaction energy of various possible chloranil-imidazole interactions suggest that the two molecules interact preferably via N and H atoms of imidazole and CO group of chloranil with their molecular planes almost perpendicular to each other.     

Comparative study of imidazole hydration: Ab initio and electrostatic calculations vs. Cambridge structural database analysis

We studied geometries and energies of complexes between water and neutral or protonated imidazole by ab initio molecular orbital calculations using the 4-31G basis set with and without the counterpoise correction. Positions of hydration sites and relative binding energies could be also estimated by using the electrostatic field map of imidazole as calculated by our bond increment method. The reliability of the calculations is confirmed by comparing the geometries of the imidazole-water complex to the experimental ones from the Cambridge Structural Database. These were obtained by X-ray diffraction studies on crystals with water bound to a molecule containing the imidazole fragment.     

Asymmetry in platinum acetylide complexes: confinement of the triplet exciton to the lowest energy ligand

To determine structure-optical property relationships in asymmetric platinum acetylide complexes, we synthesized the compounds trans-Pt(PBu3)2(C[triple bond]CC6H5)(C[triple bond]C-C6H4-C[triple bond]CC6H5) (PE1-2), trans-Pt(PBu3)2(C[triple bond]CC6H5)(C[triple bond]C-C6H4-C[triple bond]C-C6H4-C[triple bond]CC6H5) (PE1-3) and trans-Pt(PBu3)2(C[triple bond]C-C6H4-C[triple bond]CC6H5)(C[triple bond]C-C6H4-C[triple bond]C-C6H4-C[triple bond]CC6H5) (PE2-3) that have different ligands on either side of the platinum and compared their spectroscopic properties to the symmetrical compounds PE1, PE2 and PE3. We measured ground state absorption, fluorescence, phosphorescence and triplet state absorption spectra and performed density functional theory (DFT) calculations of frontier orbitals, lowest lying singlet states, triplet state geometries and energies. The absorption and emission spectra give evidence the singlet exciton is delocalized across the central platinum atom. The phosphorescence from the asymmetric complexes comes from the largest ligand. Time-dependent (TD) DFT calculations show the S1 state has mostly highest occupied molecular orbital (HOMO) --> lowest unoccupied molecular orbital (LUMO) character, with the LUMO delocalized over the chromophore. In the asymmetric chromophores, the LUMO resides on the larger ligand, suggesting the S1 state has interligand charge transfer character. The triplet state geometries obtained from the DFT calculations show distortion on the lowest energy ligand, whereas the other ligand has the ground state geometry. The calculated trend in the triplet state energies agrees very well with the experimental trend. Calculations of triplet state spin density also show the triplet exciton is confined to one ligand. In the asymmetric complexes the spin density is confined to the largest ligand. The results show Kasha's rule applies to these complexes, where the triplet exciton moves to the lowest energy ligand.     

A Theoretical Study of the Interaction Between Cytosine and BX<Subscript>3</Subscript> (X = F,Cl) Systems

The B3LYP method of DFT and HF theories of ab initio with 6-311+G^** basis sets were used to predict the geometries of the cytosine-BX₃ (X,=F, Cl) complex systems. Four conformers were obtained with no imaginary frequencies, respectively. The binding energies, enthalpies and Gibbs energies of cytosine-BX₃ have been obtained. The analyses of the combinations between cytosine and BX₃ using the natural bond orbital (NBO) method and thermodynamics indicate that the complexes (a) and (e), which depend on the proton affinities of the oxygen on the cytosine and boron in BX₃, are the most stable ones with their combination energies of −234.21 and −228.23 kJ.mol⁻¹ (B3LYP method, BSSE corrected). Based on the calculation results, a reasonable method was employed to calculate the change in the enthalpies and Gibbs energies to form eight complexes in the gaseous state at 298.15K and 101.325 kPa. It can be shown that the conformers (a) and (e) are the most stable and form readily.     

Influence of the d orbital occupation on the nature and strength of copper cation-pi interactions: threshold collision-induced dissociation and theoretical studies

Threshold collision-induced dissociation techniques are employed to determine the bond dissociation energies of a wide variety of copper cation-pi complexes, Cu(+)(pi-ligand), where pi-ligand = benzene, flurobenzene, chlorobenzene, bromobenzene, iodobenzene, phenol, toluene, anisole, pyrrole, N-methylpyrrole, indole, naphthalene, aniline, N-methylaniline, and N,N-dimethylaniline. The primary and lowest energy dissociation pathway corresponds to the endothermic loss of the intact neutral pi-ligand for all complexes except those to N-methylpyrrole, indole, aniline, N-methylaniline, and N,N-dimethylaniline. In the latter complexes, the primary dissociation pathway corresponds to loss of the intact ligand accompanied by charge transfer, thereby producing a neutral copper atom and ionized pi-ligand. Fragmentation of the pi-ligands is also observed at elevated energies in several cases. Theoretical calculations at the B3LYP/6-311G(d,p) level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Multiple low-energy conformers are found for all of the copper cation-pi complexes. Theoretical bond dissociation energies are determined from single point energy calculations at the B3LYP/6-311+G(3df,2p) level of theory using the B3LYP/6-311G(d,p) optimized geometries. The agreement between theory and experiment is very good for most complexes. The nature and strength of the binding in these copper cation-pi complexes are studied and compared with the corresponding cation-pi complexes to Na(+). Natural bond orbital analyses are carried out to examine the influence of the d orbital occupation on copper cation-pi interactions.     

胸腺嘧啶-BH3构型与性质的DFT和MP2研究

Geometries and combination energies are predicated at B3LYP / 6-31G（d）and MP2 / 6-31G（d）level for thymine-BH3 complexes and 5 geometries have been obtained. Then single point energy calculations using larger basis sets（6-311 + G（2df）and aug-cc-pVDZ）and vibrational analysis and natural bond orbital analysis are carried out on the 5 optimized conformers. The outcome indicates that the conformers with the boron atom combined with O directly are relatively stable ones,（a）and（b）,with the combination energies of 90. 4 and 88. 0 kJ / mol （B3LYP / 6-31G（d）,BSSE corrected）. The fact is that the nitrogen atom offers electron to the empty atomic orbital of boron which produces the conformers（c）and（d）. Only one conformer is found which is formed because two carbon atoms offer π electron to the empty orbital of boron. The charge transference exists in all the conformers. The combination energies have a good line relation with their charge transference. The calculated results show that when the complex forms their IR spectrum moved to the red side and the frequency shifts are relative to the stabilities of the complexes.     

Ionization Energies and Dyson Orbials of Allyl Alcohol and Allyl Mercaptan Conformers

The conformers of allyl alcohol and allyl mercaptan were studied with B3LYP/aug-cc-pVTZ method. Their relative energies were calculated at MP3, MP4(SDQ), and CCSD(T) levels. The most stable conformers for these two molecules are Gauche-gauche' (Gg'). The theo-retical photoelectron spectra simulated with the calculated ionization energies demonstrate that there are at least four conformers in allyl alcohol and four conformers in allyl mercaptan in the gas-phase experiments. The Dyson orbitals of the highest occupied molecular orbital (HOMO) and the next HOMO (HOMO-1) of allyl mercaptan Gg' conformer show strongly mixing n_S and π_C=C characteristics, which may be due to the resonance and inductive effects between π_C=C and n_S in HOMO-1 and HOMO.     

A computational characterization of boron-oxygen multiple bonding in HN=CH-CH=CH-NH-BO

Structures, relative energies, and bonding characteristics for various conformers of 3-imino-N-(oxoboryl)prop-1-en-1-amine, HN=CH-CH=CH-NH-BO, and the corresponding borocycle (-HN=CH-CH=CH-NH-B-)O are discussed using results from second-order M?ller-Plesset (MP2) perturbation theory with the Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. These MP2 results are compared to those from computationally efficient density functional theory (DFT) calculations using the LDA, PBE, TPSS, BLYP, B3LYP, BVP86, OLYP, O3LYP, and PBE1PBE functionals in conjunction with the economical Pople-type 6-311++G(d,p) basis set to evaluate the suitability of these DFT/6-311++G(d,p) levels for use with larger boron-containing systems. The effects of an aqueous environment were incorporated into the calculations using COSMO methodology. The calculated boron-oxygen bond lengths, orbital compositions, and bond orders in all the (acyclic) HN=CH-CH=CH-NH-BO conformers were consistent with the presence of a boron-oxygen triple bond, similar to that found in H-BO and H2N-BO. The (-HN=CH-CH=CH-NH-B-)O borocycle is predicted to be planar (C2v symmetry), and it is approximately 30 kcal/mol lower in energy than any of the (acyclic) HN=CH-CH=CH-NH-BO conformers; the boron-oxygen bond in this borocycle has significant double bond character, a bonding scheme for which there has been only one experimental structure reported in the literature (Vidovic, D. ; et al. J. Am. Chem. Soc. 2005, 127, 4566- 4569).     

The weak hydrogen bond in the fluorobenzene-ammonia van der Waals complex: Insights into the effects of electron withdrawing substituents on pi versus in-plane bonding

The fluorobenzene-ammonia van der Waals complex has been studied using a combination of two-color resonance enhanced multiphoton ionization (REMPI) spectroscopy, counterpoise corrected RICC2 ab initio molecular orbital calculations, and multidimensional Franck-Condon analysis. The experimental REMPI spectrum is characterized by a dominant, blueshifted band origin, and weak activity in intermolecular vibrational modes. RICC2 geometry optimizations and numerical vibrational frequency calculations of the neutral ground and first excited states have been performed on a number of different structural isomers of the complex using basis sets ranging from augmented double-zeta to quadruple-zeta level. Ground state basis set superposition error corrected zero-point binding energies show the in-plane sigma complex, forming a pseudo-six-membered ring connecting the fluorine atom and ortho-hydrogen, to be consistently the most stable of all six conformations considered, at all levels of theory. Comparison of computed zero-point excitation energies for the most stable pi and sigma conformers with fluorobenzene show that the sigma complex is the only conformer predicted to exhibit a spectral blueshift upon electronic excitation. The computed neutral ground and first excited state geometries and frequencies were used to perform multidimensional Franck-Condon simulations of the S(1)-S(0) vibronic spectrum for each of the most stable conformers. These simulations yielded null spectra for transitions involving the most stable of the pi complexes, pi(bridge); a spectrum rich in strong intermolecular vibrational structure for the second of the pi complexes, in complete contrast to the experimental spectrum; and for the sigma complex, a spectrum exhibiting weak intermolecular activity in line with that observed experimentally. This last simulation allowed an almost complete vibrational assignment of the intermolecular structure in the REMPI spectrum. The agreement between computational results and experiment overwhelmingly favors assignment of the spectrum to the in-plane sigma complex.     

Influence of geometry relaxation on the energies of the S1 and S2 states of violaxanthin, zeaxanthin, and lutein

Precise knowledge of the excitation energies of the lowest excited states S(1) and S(2) of the carotenoids violaxanthin, lutein, and zeaxanthin is a prerequisite for a fundamental understanding of their role in light harvesting and photoprotection during photosynthesis. By means of density functional theory (DFT) and time-dependent DFT (TDDFT), the electronic and structural properties of the ground and first and second excited states are studied in detail. According to our calculations, all-s-cis-zeaxanthin and s-cis-lutein conformers possess lower total ground-state energies than the corresponding s-trans conformers. Thus, only s-cis isomers are probably physiologically relevant. Furthermore, the influence of geometric relaxation on the energies of the ground state and S(1) and S(2) states has been studied in detail. It is demonstrated that the energies of these states change significantly if the carotenoid adopts the equilibrium geometry of the S(1) state. Considering these energetic effects in the interpretation of S(1) excitation energies obtained from fluorescence and transient absorption spectroscopy shifts the S(1) excitation energies about 0.2 eV to higher energy above the excitation energy of the chlorophyll a.     

C_4H_4Y(Y=O,S,Se)与BX_3(X=H,F,Cl)作用的计算化学研究

采用量子化学的密度泛函B3LYP和二阶微扰MP2(full)方法对C4H4Y(Y=O,S,Se)与BX3(X=H,F,Cl)形成的电子授受型复合物进行了研究,所得18个复合物的构型包括BX3位于C=C双键上方的π-p作用型和B与O,S,Se直接作用的n-p作用型.体系C4H4Y-BH3以n-p作用型较为稳定,体系C4H4Y-BF3,C4H4Y-BCl3的π-p和n-p作用型复合物稳定性相当.对各复合物的几何构型、振动频率和自然键轨道分析表明,复合物的形成过程中均存在几何构型的改变、电荷的转移和振动频率的变化,它们的变化规律与复合物稳定性的变化规律基本一致,即按H,F,Cl的顺序依次降低.     

羰基镍簇Ni(CO)n(n=1-4)的结构和Ni-CO键解离性质的密度泛函理论研究

在密度泛函理论框架下, 应用不同泛函计算了配合物Ni(CO)n(n=1~4)的平衡几何构型和振动频率. 考察了泛函和基组重叠误差对预测Ni—CO键解离能的影响. 计算结果表明, 用杂化泛函能得到与实验一致的优化几何构型和较合理的振动频率. 对Ni(CO)n(n=2~4)体系, 用“纯”泛函, 如BP86和BPW91, 可得到与CCSD(T)更符合、 并与实验值接近的解离能. 当解离产物出现单个金属原子或离子(如金属羰基配合物的完全解离)时, BSSE校正项的计算中应保持金属部分的电子结构一致. 只有考虑配体基组和不考虑配体基组两种情况下金属的电子构型与配合物中金属的构型一致时, 才能得到合理的BSSE校正, 从而预测合理的解离能.     

A joint QM/MD study on α-, β- and γ-cyclodextrins in selective complexation with cathinone

In this article, density functional theory (DFT) calculations and 30 ns molecular dynamic (MD) simulations were performed to investigate the ability of α-, β- and γ-cyclodextrins (CDs) to form selective complexes with cathinone. DFT calculations in the gas phase, water, chloroform and methanol reveal that the solvents, reduce the stability of the complexes. Optimized structures confirm that α-CD cannot encapsulate cathinone, completely, while other CDs showed an opposite behavior. DFT calculations indicate that cathinone has the most stable complex with γ-CD in comparison to the α- and β-CDs. Natural bond orbital and quantum theory of atoms in molecules analyses reveal that the electrostatic interactions between cathinone and CDs are the driving force of the complex formation. MD simulations confirm that different solvents play an important role in the stability of the cathinone complexes and the obtained MD results are in good agreement with the DFT calculations.     

Ab initio study of Eu3+–L (L=H2O,H2S,NH2CH3, S(CH3)2, imidazole) complexes

The ground states and binding energies of Eu³⁺–L (L=H₂O,H₂S,NH₂CH₃,S(CH₃)₂, imidazole) complexes has been determined using ab initio techniques. The binding is mostly electrostatic as expected. The empty f orbital is different for the S compounds, being a π-like orbital, while for the O and N containing ligands it is a σ-like orbital. However, the range in the binding energies for the different f holes is small.     

Excited‐State Hydrogen Bonding Strengthening and Weakening with Intramolecular Charge Transfer in Resorufin‐Water Complexes: A TD‐DFT Study

The geometric structures, infrared spectra and hydrogen bond binding energies of the various hydrogen‐bonded Res^?‐water complexes in states S0 and S1 have been calculated using the density functional theory (DFT) and time‐dependent density functional theory (TD‐DFT) methods, respectively. Based on the changes of the hydrogen bond lengths and binding energies as well as the spectral shifts of the vibrational mode of the hydroxyl groups, it is demonstrated that hydrogen bonds HB‐II, HB‐III and HB‐IV are strengthened while hydrogen bond HB‐I is weakened in the four singly hydrogen‐bonded Res^?‐Water complexes upon photoexcitation. When the four hydrogen bonds are formed simultaneously between one resorufin anion and four water molecules in the Res^?‐4Water complex, all the hydrogen bonds are weakened in both the ground and excited states compared with those in the corresponding singly hydrogen‐bonded Res^?‐Water complexes. Furthermore, in complex Res^?‐4Water, hydrogen bonds HB‐II and HB‐IV are strengthened while hydrogen bonds HB‐I and HB‐III are weakened after the electronic excitation. The hydrogen bond strengthening and weakening in the various hydrogen‐bonded Res^?‐water complexes should be due to the redistribution of the charges among the four heteroatoms (O_1‐3 and N₁) within the resorufin molecule upon the optical excitation.     

Hypothesis of Hemoprotein Sensor Confirmed by Ab initio Quantum-Chemical Method

The nature of the chemical bond of complexes of iron and cobalt porphyrinates with ligands is studied by the quantum-chemical method in the Hartree–Fock self-consistent field approximation using the 3-21G basis set. The addition of oxygen molecule to the MP and MPIm complexes (M = Fe, Co; Im is imidazole) is established to be more favorable than water addition. However, imidazole, which is the second ligand in the MPImO₂ and MPImH₂O complexes (M = Fe, Co), increases the M–O₂ and M–H₂O binding energies for iron, but decreases them for cobalt. The Co atom is bound with the porphyrin ring more strongly than the iron atom due to the larger total overlap of the atomic orbitals. The calculations of the binding energy in the complexes demonstrate similar changes in the structures of the spatial conformation of the deoxy form (FeP + H₂O) of iron porphyrinate and the oxy form (CoP + O₂) of cobalt porphyrinate. This is an argument in favor of the hypothesis of hemoprotein sensor of partial oxygen stress in tissues.     

Theoretical study of quinolines-I(2) intermolecular interaction and implications on dye-sensitized solar cell performance

The monomer and intermolecular charge-transfer complexes of 13 different quinoline derivatives with diiodine were studied using ab initio molecular orbital (MO) and density functional theory (DFT) methods. Calculations revealed that the sigma* orbital of iodine interacts with the nitrogen lone pair in the quinoline ring. The open-circuit photovoltage (V(oc)) values of an Ru(II) complex dye-sensitized nanocrystalline TiO(2) solar cell with an I(-)/I(3) (-) redox electrolyte in acetonitrile using quinoline additives were compared to the computational calculations on the intermolecular interaction between quinolines and I(2). The optimized geometries, frequency analyses, Mulliken population analyses, natural bond orbital (NBO) analyses, and interaction energies indicate that the V(oc) value of the solar cell is higher when quinoline complexes more favorably interact with I(2). Therefore, the interaction between the quinoline additives and iodine redox electrolyte is an important factor for controlling dye-sensitized solar cell performance.     

High-resolution electron spectroscopy, preferential metal-binding sites, and thermochemistry of lithium complexes of polycyclic aromatic hydrocarbons

Polycyclic aromatic hydrocarbons are model systems for studying the mechanisms of lithium storage in carbonaceous materials. In this work, Li complexes of naphthalene, pyrene, perylene, and coronene were synthesized in a supersonic metal-cluster beam source and studied by zero-electron-kinetic-energy (ZEKE) electron spectroscopy and density functional theory calculations. The adiabatic ionization energies of the neutral complexes and frequencies of up to nine vibrational modes in the singly charged cations were determined from the ZEKE spectra. The metal-ligand bond energies of the neutral complexes were obtained from a thermodynamic cycle. Preferred Li∕Li(+) binding sites with the aromatic molecules were determined by comparing the measured spectra with theoretical calculations. Li and Li(+) prefer the ring-over binding to the benzene ring with a higher π-electron content and aromaticity. Although the ionization energies of the Li complexes show no clear correlation with the size of the aromatic molecules, the metal-ligand bond energies increase with the extension of the π-electron network up to perylene, then decrease from perylene to coronene. The trends in the ionization and metal-ligand bond dissociation energies of the complexes are discussed in terms of the orbital energies, local quadrupole moments, and polarizabilities of the free ligands and the charge transfer between the metal atom and aromatic molecules.