Stabilization of zwitterions in solution: phosphinic and phosphonic acid GABA analogues

The influence of treatment of electron correlation, size of basis set and choice of solvation model on the predicted stabilization of zwitterionic phosphinic and phosphonic acid gamma-aminobutyric acid (GABA) analogues is investigated using ab initio molecular orbital methods and density functional theory. Density functional theory with the B3LYP functional and a composite basis set composed of the 6-31+G(2df) basis for phosphorus atoms and the 6-31+G(d,p) basis set for all other atoms is found to give an acceptable tradeoff between accuracy and computational expense. Either directly optimizing zwitterionic conformers within the conductor-like screening solvation model (COSMO) or pre-optimizing as dihydrates and then applying the COSMO solvation model give an acceptable treatment of solvation in terms of determining stable solvated structures, although directly optimizing within COSMO is simpler and less computationally expensive. With this protocol, cis-constrained phosphinic and phosphonic acid GABA analogues, which exhibit lower affinities for GABAC receptors, are found to possess only folded, intramolecularly hydrogen bonded low energy conformers. Trans-constrained analogues, on the other hand, exhibit higher affinities for GABAC receptors and are found to exist only as partially folded stable conformers. Conformationally flexible analogues can attain folded, partially folded or fully extended conformations and also have high biological activity. These results imply that the partially folded conformation is likely to be the most biologically active.     

Prediction of iron K-edge absorption spectra using time-dependent density functional theory

Iron K-edge X-ray absorption pre-edge features have been calculated using a time-dependent density functional approach. The influence of functional, solvation, and relativistic effects on the calculated energies and intensities has been examined by correlation of the calculated parameters to experimental data on a series of 10 iron model complexes, which span a range of high-spin and low-spin ferrous and ferric complexes in O(h) to T(d) geometries. Both quadrupole and dipole contributions to the spectra have been calculated. We find that good agreement between theory and experiment is obtained by using the BP86 functional with the CP(PPP) basis set on the Fe and TZVP one of the remaining atoms. Inclusion of solvation yields a small improvement in the calculated energies. However, the inclusion of scalar relativistic effects did not yield any improved correlation with experiment. The use of these methods to uniquely assign individual spectral transitions and to examine experimental contributions to backbonding is discussed.     

态定(线性响应)-极化连续模型/含时密度泛函方法研究一种有机发光材料的吸收和发射光谱

3-(二氰亚甲基)-5,5-二甲基-1-(3-[9-(2-乙基-己基)-咔唑基]-乙烯基)环己烷(DCDHCC)是一种用于光电器件中的有机染料,它具有良好的发光特性.我们使用含时密度泛函方法(TD-PBE0,TD-BMK和TD-M06)以及极化连续模型(PCM)计算了该材料在溶剂中的吸收和发射特性.计算中使用了线性响应(LR)、态定(SS)两种溶剂模型和6-31G(d)、6-31+G(d,p)两种基组.计算了DCDHCC在苯、四氢呋喃和丙酮溶剂中的吸收和发射光谱,并与实验观测进行了比较.结果表明:对于吸收光谱的计算,杂化函数的影响大于基组和溶剂模型,在三种函数中BMK更适于研究DCDHCC的吸收光谱;而对于发射光谱,基组的影响最大,基组通过影响激发态构型从而影响发射光谱,对于激发态构型的优化需要使用6-31+G(d,p)基组.我们希望这些研究能对今后设计类似的发光分子有帮助.     

The molecular, electronic structures and vibrational spectra of metal-free, N,N'-dideuterio and magnesium tetra-2,3-pyridino-porphyrazines: Density functional calculations

A theoretical investigation of the fully optimized geometries and electronic structures of the metal-free (TPdPzH(2)), N,N'-dideuterio (TPdPzD(2)), and magnesium (TPdPzMg) tetra-2,3-pyridino-porphyrazine has been conducted based on density functional theory. The optimized geometries at density functional theory level for these compounds are reported here for the first time. A comparison between the different molecules for the geometry, molecular orbital, and atomic charge is made. The substituent effect of the N atoms on the molecular structures of these compounds is discussed. The IR and Raman spectra for these three compounds have also been calculated at density functional B3LYP level using the 6-31G(d) basis set. Detailed assignments of the NH, NM, and pyridine ring vibrational bands in the IR and Raman spectra have been made based on assistance of animated pictures. The simulated IR spectra of TPdPzH(2) are compared with the experimental absorption spectra, and very good consistency has been found. The isotope effect on the IR and Raman spectra is also discussed.     

Structures and excitation energies of Zn–tetraarylporphyrin analogues: A theoretical study

The photophysical properties of β-substituted Zn–tetraarylporphyrin (ZnTAP) analogues used as dyes in dye-sensitized solar cells were studied using density functional theory (DFT). Singlet-excitation energy calculations of ZnTAP analogues were performed using time-dependent DFT with B3LYP, B3PW91, PBE0 exchange–correlation functionals at 6-31G(d) and 6-31+G(d) basis sets using B3LYP/6-31G(d) geometries. The PBE0 functional at 6-31+G(d) basis set provided a better correlation with the experimental data for both B- and Q-bands. The inclusion of solvation effect in the calculations provided a good agreement in terms of B:Q_ave ratio of the oscillator strengths for both analogues with the experimental values. Analogue 2 has a higher and a more balanced charge-carrier transport rates than analogue 1. In general, the addition of an electron-donating group in the meso-substituent (analogue 2) resulted in a narrower band gap, higher oscillator strength, a more red-shifted absorption spectra, and better charge-transfer characteristics than analogue 1.     

Spectroscopic, quantum chemical DFT/HF study and synthesis of [2.2.1] hept-2'-en-2'-amino-N-azatricyclo [3.2.1.0(2,4)] octane

The compound 4-N-bicyclo [2.2.1] hept-2'-en-2'-amino-N-azatricyclo [3.2.1.0(2,4)] octane (2) has been synthesized and characterized by elemental analysis, IR, UV-vis, mass and NMR. Density functional theory (DFT) and Hartree-Fock (HF) calculations have been carried out for the title compound by using the standard 6-31G* basis set. The calculated results show that the predicted geometry can well reproduce the structural parameters. Predicted vibrational frequencies have been assigned and compared with experimental IR spectra and they complement each other. The theoretical electronic absorption spectra have been calculated by using CIS, TD-DFT and ZINDO methods. The (13)C NMR and (1)H NMR of compound (2) have been calculated by means of Becke 3-Lee-Yang-Parr (B3LYP) density functional method with 6-31G* basis set. Comparison between the experimental and the theoretical results indicates that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties. On the basis of vibrational analyses, the thermodynamic properties of the title compound at different temperatures have been calculated.     

MP2, DFT and ab initio calculations on thioxanthone

Density functional theory (DFT), HF and MP2 calculations have been carried out to investigate thioxanthone molecule using the standard 6-31+G(d,p) basis set. The results of MP2 calculations show a butterfly structure for thioxanthone. The calculated results show that the predicted geometry can well reproduce the structural parameters. The predicted vibrational frequencies were assigned and compared with experimental IR spectra. A good harmony between theory and experiment is found. The theoretical electronic absorption spectra have been calculated using CIS method. ¹³C and ¹H NMR of the title compound have been calculated by means of B3LYP density functional method with 6-31+G(d,p) basis set. The comparison of the experimental and the theoretical results indicate that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties.     

Investigation of the structure, the optical properties, and the photophysics of some indolocarbazoles having terminal aromatic rings

Structures, optical properties, and photophysics of ladder indolo[3,2-b]carbazoles substituted symmetrically by phenylene and thiophene rings have been investigated theoretically and experimentally. The ground state optimized structures were obtained using the density functional theory (DFT) as approximated by the B3LYP functional and employing the 6-31G^* basis set. All derivatives were found nonplanar in their electronic ground states. The character and the energy of the singlet–singlet electronic transitions have been investigated by applying the time-dependent density functional theory (TDDFT) to the correspondingly optimized-ground-state geometries. The ab initio restricted configuration interaction (singles) method (RCIS/6-31G^*) was adopted to obtain the first singlet excited-state structures (S₁) of the molecule. TDDFT calculations performed on the S₁ optimized geometries was used to obtain emission energies. UV–vis and fluorescence spectroscopies were analyzed in conjunction with theoretical calculations. The computed excitation and emission energies were found in reasonable agreement with the experimental absorption and fluorescence spectra. Finally, the photophysical behavior of the indolocarbazoles have been studied by means of steady state and time resolved fluorescence. The overall data have allowed the determination of the rate constants for the radiative and nonradiative decay processes. Both theoretical and experimental data show that the replacement of phenylene rings by thiophene units induces a red shift in the absorption and fluorescence spectra. This behavior is interpreted in terms of the electron donor properties of the thiophene ring. On the other hand, the change of the substitutional pattern, from 2,8 to 3,9, causes a significant hypsochromic shift of the absorption and fluorescence bands.     

Optical absorption and emission properties of fluoranthene, benzo[k]fluoranthene, and their derivatives. A DFT study

Fluoranthene and benzo[k]fluoranthene-based oligoarenes are good candidates for organic light-emitting diodes (OLEDs). In this work, the electronic structure and optical properties of fluoranthene, benzo[k]fluoranthene, and their derivatives have been studied using quantum chemical methods. The ground-state structures were optimized using the density functional theory (DFT) methods. The lowest singlet excited state was optimized using time-dependent density functional theory (TD-B3LYP) and configuration interaction singles (CIS) methods. On the basis of ground- and excited-state geometries, the absorption and emission spectra have been calculated using the TD-DFT method with a variety of exchange-correlation functionals. All the calculations were carried out in chloroform medium. The results show that the absorption and emission spectra calculated using the B3LYP functional is in good agreement with the available experimental results. Unlikely, the meta hybrid functionals such as M06HF and M062X underestimate the absorption and emission spectra of all the studied molecules. The calculated absorption and emission wavelength are more or less basis set independent. It has been observed that the substitution of an aromatic ring significantly alters the absorption and emission spectra.     

Profiling energetics and spectroscopic signatures in prototropic tautomers of asymmetric phthalocyanine analogues

Density functional theory (DFT) and time-dependent density functional theory (TDDFT) were used to explain discrepancies in UV-vis and MCD spectra of the metal-free tribenzo[b,g,l]thiopheno[3,4-q]porphyrazine (1), substituted tribenzo[b,g,l]porphyrazine (2), and 2,3-bis(methylcarboxyl)phthalocyanine (3). On the basis of gas-phase and polarized continuum solvation model (PCM) DFT and TDDFT calculations, it was suggested that both NH tautomers contribute to the spectroscopic signature of 1, whereas the Q-band region of 2 and 3 is dominated by a single NH tautomer. For all tested compounds, it was found that the combination of the BP86 exchange-correlation functional, 6-31G(d) basis set, and TDDFT-PCM approach provides the best accuracy in energies of the Q(x)- and Q(y)-bands of the individual NH tautomers as well as correctly describes their relative energy differences, which are important in understanding of experimental spectroscopy of the target systems.     

Solvatochromic and thermochromic shifts of electronic spectra of polar solute molecules in a mixture of polar and nonpolar solvent; the role of solvent-solvent interactions

A theoretical model is proposed to describe the influence of the concentration of a polar solvent and the temperature of a solution on the electronic spectra of a polar solute in a binary solvent mixture. It is shown that the interaction between molecules of the polar solvent in the first solvation shell makes the significant contribution to the formation of absorption and fluorescence bands of the solute. An experimental study of solvatochromic and thermochromic shifts of steady-state fluorescence spectra of 3-amino-N-methylphthalimide in decalin--propanol mixture for different values of propanol mole fraction is carried out. Good qualitative agreement between the experimental data and calculation results is observed.     

From Oligo(Phenyleneethynylene) Monomers to Supramolecular Helices: The Role of Intermolecular Interactions in Aggregation

Supramolecular helices that arise from the self-assembly of small organic molecules via non-covalent interactions play an important role in the structure and properties of the corresponding materials. Here we study the supramolecular helical aggregation of oligo(phenyleneethynylene) monomers from a theoretical point of view, always guiding the studies with experimentally available data. In this way, by systematically increasing the number of monomer units, optimized n-mer geometries are obtained along with the corresponding absorption and circular dichroism spectra. For the geometry optimizations we use density functional theory together with the B3LYP-D3 functional and the 6–31G** basis set. For obtaining the spectra we resort to time-dependent density functional theory using the CAM-B3LYP functional and the 3–21G basis set. These combinations of density functional and basis set were selected after systematic convergence studies. The theoretical results are analyzed and compared to the experimentally available spectra, observing a good agreement.     

Time dependent density functional theory calculations for electronic circular dichroism spectra and optical rotations of conformationally flexible chiral donor-acceptor dyad

Twelve conformations of a chiral donor-acceptor (charge-transfer) dyad and six conformations of its dimer complex were structurally optimized by using the Kohn-Sham density functional theory (BLYP/TZV2P) incorporating a recently developed empirical correction scheme that uses C6/R6 potentials for van der Waals interactions (DFT-D). Subsequent time-dependent DFT calculations with BH-LYP and B3-LYP functionals (with triple-zeta basis set) were performed to obtain theoretical circular dichroism (CD) spectra. The experimental CD spectra obtained independently were properly reproduced by averaging the calculated spectra of individual conformers according to a Boltzmann population derived from single-point SCS-MP2 energies. The optical rotations of the monomer were also calculated by using the same functionals with an aug-cc-pVDZ basis set. Dielectric continuum solvation models (COSMO) applied to correct the relative energies from the isolated molecule calculations resulted in conformer distributions that piled the same or even poorer level of agreement with the experimental CD spectrum. Our results clearly show the advantage of the DFT-D method for the geometry optimization of large systems with donor-acceptor interactions and the TD-DFT/BH-LYP calculations for reproducing the experimental CD spectra. As compared with the calculated optical rotations, the wealthy information embedded in the experimental/calculated CD spectra is requisite for the configurational and/or conformational analyses of relatively large and flexible chiral organic molecules in solution.     

Direct determination of absolute configuration of methyl-substituted phenyloxiranes: combined experimental and theoretical approach

Three possible methyl-substituted phenyloxiranes have been synthesized in enantioenriched form (89-99% enantiomeric excess (ee)), and their vibrational absorption (VA) and vibrational circular dichroism (VCD) spectra have been recorded. The experimental spectra are compared to theoretical spectra obtained from quantum mechanical calculations (density functional theory with the B3LYP hybrid exchange correlation functional with 6-31++G*, aug-cc-pVDZ, or aug-cc-pVTZ basis set) and related to the physical structure of the compounds. The absolute configuration could be established directly in each case by comparing experimental and theoretical spectra. In addition, we have been able to document the changes that occur both in structures and in the VA and VCD spectra due to substituent effects on the oxirane ring.     

Spectral shifts of the n → π* and π → π* transitions of uracil based on a modified form of solvent reorganization energy

According to our recent studies on the nonequilibrium solvation, the solvent reorganization energy is found to be the cost of maintaining the residual polarization P', which equilibrates with the extra electric field E(ex). On the basis of this solvent reorganization energy and the well-established equilibrium solvation energy, a novel and reasonable expression for the spectral shift of the electronic absorption spectra is proposed in this work. Furthermore, the two lowest transitions of uracil in aqueous solution are investigated as test cases with the TDDFT/6-311++G** method. The obtained spectral shift is 0.48 eV for n → π* transition and -0.14 eV for π → π* transition, agreeing well with available experimental results. The contributions to the shift are discussed and the electrostatic plus polarization components are found to be crucial for the electronic absorption spectra of uracil in aqueous solution.     

The n → π* transition in nitroso compounds: A TD-DFT study

The visible spectra of a series of various nitroso dyes have been evaluated by using time-dependent density functional theory approach and explicitly taking into account bulk solvent effects. Using the PBE0 functional with the 6-311++G(3d,3p) atomic basis set, the agreement between theoretical and experimental n → π* transition energies is excellent for all nitroso species. The mean absolute error is limited to 0.05 eV and the maximal deviation is 0.21 eV for a set of 22 cases. As a blind test, we have predicted the absorption spectra for unknown species, in order to assess the impact of the atom bound to the NO group.     

Experimental and theoretical study on the absorption and fluorescence properties of substituted aryl hydrazones of 1,8-naphthalimide

Absorption and fluorescence spectra in acetonitrile for a series of substituted aryl hydrazones of N-hexyl-1,8-naphthalimide are studied with the aim of potential application of the compounds for enzyme activity localization. The influence of the substituents on the spectral characteristics has been evaluated. The absorption and fluorescence energies of substituted aryl-1,8-naphthalimide hydrazones have been calculated with the PCM TDDFT formalism. The M06 and PBE0 functionals, combined with the 6-31+G(d) atomic basis set, have been found to accurately model the excited state properties of the present set of solvated fluorophores. Absorption and fluorescence spectral characteristics have been rationalized in terms of experimental and theoretical electronic indices in order to assess their predictive abilities for application in designing analogues with good emitting properties. An excellent linear dependence is established between the experimental fluorescence and Hammett σ(p)(+) substituent constants and on the other hand σ(p)(+) constants correlate with the theoretically calculated values for the electrostatic potential at nuclei (EPN). A model for predicting the fluorescence properties of substituted hydrazones by means of EPN is drawn, including the polysubstituted derivatives, where Hammett constants are not applicable.     

Charge-dependent cavity radii for an accurate dielectric continuum model of solvation with emphasis on ions: aqueous solutes with oxo, hydroxo, amino, methyl, chloro, bromo, and fluoro functionalities

Dielectric continuum solvation models are widely used because they are a computationally efficacious way to simulate equilibrium properties of solutes. With advances that allow for molecular-shaped cavities, they have reached a high level of accuracy, in particular for neutral solutes. However, benchmark tests show that existing schemes for defining cavities are unable to consistently predict accurately the effects of solvation on ions, especially anions. This work involves the further development of a protocol put forth earlier for defining the cavities of aqueous solutes, with resulting advances that are most striking for anions. Molecular cavities are defined as interlocked spheres around atoms or groups of atoms in the solute, but the sphere radii are determined by simple empirically based expressions involving the effective atomic charges of the solute atoms (derived from molecular electrostatic potential) and base radii. Both of these terms are optimized for the different types of atoms or functional groups in a training set of neutral and charged solutes. Parameters in these expressions for radii were fitted by minimizing residuals between calculated and measured standard free energies of solvation (DeltaG(s)*), weighted by the uncertainty in the measured value. The calculations were performed using density functional theory with the B3LYP functional and the 6-311+G** basis set and the COnductor-like Screening MOdel (COSMO). The optimized radii definitions reproduce DeltaG(s)* of neutral solutes and singly charged ions in the training set to within experimental uncertainty and, more importantly, accurately predict DeltaG(s)* of compounds outside the training set, in particular anions (J. Phys. Chem. A 2003, 107, 5778). Inherent to this approach, the cavity definitions reflect the strength of specific solute-water interactions. We surmise that this feature underlies the success of the model, referred to as the CD-COSMO model for Charge-Dependent (also Camaioni-Dupuis) COSMO model. These findings offer encouragement that we can keep extending this scheme to other functional groups and obtain better accuracy in using continuum solvation models to predict equilibrium properties of aqueous ionic solutes. The approach is illustrated for a number of test cases, including the determination of acidities of an amine base, a study of the tautomerization equilibrium of a zwitterionic molecule (glycine), and calculating solvation energies of transition states toward a full characterization of reaction pathways in aqueous phase, here in S(N)2 exchange reactions. The calculated reaction barriers in aqueous solution are in excellent agreement with experimental values.     

Ab initio procedure for aqueous-phase pKa calculation: the acidity of nitrous acid

We present an ab initio procedure for accurately calculating aqueous-phase pKa values and apply it to study the acidity of nitrous acid (HNO2, or HONO). The aqueous-phase pK(a) of nitrous acid was obtained from calculated gas-phase acidities and solvation free energies via a thermodynamic cycle and the solvation model chemistry of Barone et al. (J. Chem. Phys. 1997, 107, 3210). Solvation free energies were calculated at the HF/6-31G(d) level using the dielectric-polarizable continuum and the integral equation formalism-polarizable continuum solvent models (D-PCM and IEF-PCM, respectively), with the D-PCM model yielding the most accurate pKa values. For HF free energies of solvation, significant improvements in accuracy could be made by moving to the larger 6-311++G(3df,3pd) and aug-cc-pVQZ basis sets. Solvation free energies were also calculated using the density functional theory (DFT) methods B3LYP, TPSS, PBE0, B1B95, VSXC, B98 and O3LYP, with the most accurate methods being TPSS and VSXC, which provided average errors of less than 0.11 pKa units. Solvation free energies calculated with the different DFT methods were relatively insensitive to the basis set used. Our theoretical calculations are compared with experimental results obtained using stopped flow spectrophotometry. The pKa of nitrous acid was measured as 3.16 at 25 degrees C, and the enthalpy and entropy of nitrous acid dissociation were calculated from measurements as 6.7 kJ mol(-1) and -38.4 J mol(-1) K(-1), respectively, between 25 and 45 degrees C. The UV/visible absorption spectra of the nitrite ion and nitrous acid were also examined, and molar extinction coefficients were obtained for each.     

Successful application of TD-DFT in transient absorption spectra assignment

[Reaction: see text] Despite the wide application of time-resolved spectra technology in investigating various chemical and biological processes, it is still a great challenge to give an unambiguous spectra assignment for observed transient species, e.g., the curcumin-derived radicals. Here we report that time-dependent density functional theory with the B3LYP functional and an appropriate basis set can give reliable absorption spectra for carbon- and oxygen-centered radicals derived from curcumin, which indicates the potential of theoretical methods in helping assign transient spectra.