DFT study of linear and nonlinear optical properties of donor-acceptor substituted stilbenes, azobenzenes and benzilideneanilines

A theoretical analysis of the linear and nonlinear optical properties of six push–pull π-conjugated molecules with stilbene, azobenzene and benzilideneaniline as a backbone is presented. The photophysical properties of the investigated systems were determined by using response functions combined with density functional theory (DFT). Several different exchange-correlation potentials were applied in order to determine parameters describing the one- and two-photon spectra of the studied molecules. In particular, the recently proposed Coulomb-attenuated model (CAM-B3LYP) was used to describe charge-transfer (CT) excited states. In order to compare theoretical predictions with available experimental data, calculations with inclusion of solvent effects were performed. The BLYP and the CAM-B3LYP functionals were found to yield values of two-photon absorption (TPA) probabilities closer to experimental values than the B3LYP functional or the HF wavefunction. Moreover, molecular static hyperpolarisabilities were determined using both DFT and second-order Møller-Plesset perturbation (MP2) theory. Likewise, the CAM-B3LYP functional was found to outperform other applied exchange-correlation potentials in determining first hyperpolarisability (β). Moreover, it was confirmed on a purely theoretical basis that the presence of a –C=C– bridge between the phenyl rings leads to a much larger nonlinear optical response in comparison with a –N=N– bridge.     

Spectroscopic investigation on glutamic acid by Coulomb-attenuating and double hybrid density functional theory methods

This study deals with the identification of glutamic acid by means of quantum chemical approach. FT-IR, FT-Raman and UV–vis spectra were recorded in the region 4000–400, 4000–50 cm^? 1 and 200–600 nm, respectively. CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p) calculations were performed to obtain the optimised molecular structures, vibrational frequencies and corresponding vibrational assignment, thermodynamic properties and natural bonding orbital (NBO) analysis. The results show that the obtained optimised geometric parameters (bond lengths, bond angles and bond dihedrals) and vibrational frequencies were found to be in good agreement with the experimental results. The calculations of the electronic spectra were compared with the experimental ones. Furthermore, highest occupied molecular orbital and lowest unoccupied molecular orbital analyses and UV–vis spectral analysis were also performed to determine the energy band gaps and transition states. NBO analysis, calculated using density functional theory methods (CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p)), was induced to find inter-molecular atoms. ¹³C and ¹H NMR isotropic chemical shifts were calculated and the assignments made were compared with the ChemDraw Ultra values.     

Substituent effects on the properties of fluorene-thieno[3,4-b]pyrazine derivatives for light-emitting applications

Five thieno[3,4-b]pyrazine-based model compounds were studied to explore the effects of the substituent groups (alkyl or aryl) on the structure, atomic charge, optical properties, ionization potential (IP), electron affinity (EA), and reorganization energy. Theoretical calculations were carried out by density functional theory (DFT) using the B3LYP hybrid function combined and CAM-B3LYP with the 6-31G(d) basis set. The lowest-lying absorption and emission spectra of 9,9'-diethylhexylfluorene-alt-5,7-dithien-2-yl-thieno[3,4-b] pyrazine (FDDTTP) with alkyl groups showed a blue-shift, while those of FDDTTP with aryl groups exhibited a red-shift. The results agree well with analytical data from reorganization energies. IPs are brought down by both alkyl and aryl groups. However, EAs are raised only by aryl units. The results indicate that aryl groups are more helpful in forming excitions for FDDTTP molecules. Consequently, FDDTTP with aryl groups are more efficient acceptor segments for designing donor-acceptor copolymers than those with alkyl groups.     

Many-body energies during proton transfer in an aqueous system

The energetics of the mechanism of proton transfer from a hydronium ion to one of the water molecules in its first solvation shell are studied using density functional theory and the Møller–Plesset perturbation (MP2) method. The potential energy surface of the proton transfer mechanism is obtained at the B3LYP and MP2 levels with the 6-311++G** basis set. Many-body analysis is applied to the proton transfer mechanism to obtain the change in relaxation energy, two-body, three-body and four-body energies when proton transfer occurs from the hydronium ion to one of the water molecules in its first solvation shell. It is observed that the binding energy (BE) of the complex decreases during the proton transfer process at both levels of theory. During the proton transfer process, the % contribution of the total two-body energy to the binding energy of the complex increases from 62.9 to 68.09% (39.9 to 45.95%), and that of the total three-body increases from 25.9 to 27.09% (24.16 to 26.17%) at the B3LYP/6-311++G** (MP2/ 6-311++G**) level. There is almost no change in the water–water–water three-body interaction energy during the proton transfer process at both levels of theory. The contribution of the relaxation energy and the total four-body energy to the binding energy of the complex is greater at the MP2 level than at the B3LYP level. Significant differences are found between the relaxation energies, the hydronium–water interaction energies and the four-body interaction energies at the B3LYP and MP2 levels.     

On the potential application of DFT methods in predicting the interaction-induced electric properties of molecular complexes. Molecular H-bonded chains as a case of study

A detailed analysis of the selected DFT functionals for the calculations of interaction-induced dipole moment, polarizability and first-order hyperpolarizability has been carried out. The hydrogen-bonded model chains consisting of HF, H(2)CO and H(3)N molecules have been chosen as a case study. The calculations of the components of the static electric properties using the diffuse Dunning's basis set (aug-cc-pVDZ) have been performed employing different types of density functionals (B3LYP, LC-BLYP, PBE0, M06-2X and CAM-B3LYP). Obtained results have been compared with those gained at the CCSD(T) level of theory. The counterpoise correction scheme, namely site-site function counterpoise, has been applied in order to eliminate basis set superposition error. The performed tests allow to conclude that the DFT functionals can provide a useful tool for prediction of the interaction-induced electric properties, however a caution has to be urged to their decomposition to the two- and many-body terms.     

The electrostatic embedding contribution to DFT calculations of ligand-amino acid residues interaction

In this work, we demonstrate that the inclusion of long-range interactions has a significant impact on the estimation of ligand–protein binding energies. Within the scope of the electrostatically embedded adaptation of the molecular fragmentation with conjugated caps (EE-AMFCC) scheme, we unveil the role played by long-range contributions in distinct levels of quantum mechanical calculations. As a prototypical system, we consider ibuprofen coupled to the human serum albumin. In particular, we show that some relevant ligand–residue interaction energies can only be accurately captured in density functional theory (DFT) approaches when the electrostatic background is properly represented by an explicit point charge distribution.

Open image in new window

Graphical Abstract (left) The binding site FA3/FA4 of HSA containing the attached IBU. (right) Absolute value of difference between the biding energies calculated including the electrostatic embedding and the energies calculated without the electrostatic embedding using the HF, B3LYP, CAM-B3LYP, and MP2 methodologies

     

“Russian doll” complexes of [n]cycloparaphenylenes: a theoretical study

The formation of "Russian doll" complexes consisting of [n]cycloparaphenylenes was predicted using quantum chemistry tools. The electronic structures of multiple inclusion complexes containing up to four macrocycles were explored at the M06-2X/6-31G* level of theory. The binding energy between the macrocycles increases from the center to the periphery of the complex and can be >60?kcal?mol(-1) for macrocycles containing 14 and 19 repeating units. It has been demonstrated that additional electrostatic interactions originating from the asymmetric electron density distribution observed when comparing the concave and convex macrocycle sides are responsible for the high binding energies in these Russian doll complexes. Oxidation or reduction of the Russian doll complexes creates polarons that are delocalized across the complexes. In the case of polaron cations, most of the polarons are localized at the macrocycle with the smallest ionization potential; for polaron anions, the negative charge is localized across the outer rings of the complex. Because anion polarons are more delocalized than cation polarons, the relaxation energies of the polaron anions were found to be smaller than those of the polaron cations.     

Excitons in Metal‐Halide Perovskites

The unprecedented increase of the power conversion efficiency of metal‐halide perovskite solar cells has significantly outpaced the understanding of their fundamental properties. One of the biggest puzzles of perovskites has been the exciton binding energy, which has proved to be difficult to determine experimentally. Many contradictory reports can be found in the literature with values of the exciton binding energy from a few meV to a few tens of meV. In this review the results of the last few years of intense investigation of the exciton physic in perovskite materials are summarized. In particular a critical overview of the different experimental approaches used to determine exciton binding energy is provided. The problem of exciton binding energy in the context of the polar nature of perovskite crystals and related polaron effects which have been neglected to date in most of work is discussed. It is shown that polaron effects can reconcile at least some of the experimental observations and controversy present in the literature. Finally, the current status of the exciton fine structure in perovskite materials is summarized. The peculiar carrier–phonon coupling can help to understand the intriguing efficiency of light emission from metal‐halide perovskites.     

A molecular orbital selection approach for fast calculations of specific rotation with density functional theory

In this work, we describe a simple approach to select the most important molecular orbitals (MOs) to compute the optical rotation tensor through linear response (LR) Kohn-Sham density functional theory (KS-DFT). Taking advantage of the iterative nature of the algorithms commonly used to solve the LR equations, we select the MOs with contributions to the guess perturbed density that are larger than a certain threshold and solve the LR equations with the selected MOs only. We propose two criteria for the selection, and two definitions of the selection threshold. We then test the approach with two functionals (B3LYP and CAM-B3LYP) and two basis sets (aug-cc-pVDZ and aug-cc-pVTZ) on a set of 51 organic molecules with specific rotation spanning five orders of magnitude, 10⁰–10⁴ deg (dm⁻¹ (g/mL)⁻¹). We show that this approach indeed can provide very accurate values of specific rotation with estimated speedup that ranges from 2 to 8× with the most conservative selection criterion, and up to 20 to 30× with the intermediate criterion.     

Experimental and theoretical DFT (B3LYP,X3LYP,CAM-B3LYP and M06-2X) study on electronic structure,spectral features,hydrogen bonding and solvent effects of 4-methylthiadiazole-5-carboxylic acid

Detailed structural, electronic and spectroscopic study of 4-methylthiadiazole-5-carboxylic acid, one of the simplest 1,2,3-thiadiazole derivatives has been performed using density functional theory at four different functionals (B3LYP, X3LYP, CAM-B3LYP and M06-2X). The two possible conformers and their dimeric forms have been investigated for the stability and hence for the calculation of molecular properties of the title compound. Vibrational analysis has been performed with the help of experimental FT-IR and FT-Raman spectra. NBO analysis has been performed to estimate the N–H—O=C hydrogen bond strength and to evaluate the intra and inter molecular charge transfer in the system. Intermolecular hydrogen-bond strength has also been computed using Atoms in Molecules (AIM) theory. To visualise spatial domain, key sites of electron transitions and electron density difference between ground as well as excited states, and their 2D and 3D plots have been computed. Solvent effect on the intermolecular hydrogen bonding have also been investigated using solvents of different polarities. Non-linear optical properties, molecular electrostatic potential surface map (MESP), thermodynamic potentials at different temperatures have also been computed and plotted.     

A comparative study of semi-squaraine and squaraine dyes using computational techniques: tuning the charge transfer/biradicaloid character by substitution

Semi-squaraines (SMSQ) are known as donor-acceptor (D-A) type molecules whereas squaraines (SQ), which differs from SMSQ by an extra donor group, are more or less biradicaloids in nature. The effect of the additional donor group in SQ, which changes the nature of the molecule, on geometrical and electronic structure are studied here and compared with the corresponding SMSQ. It is noticed from the geometrical parameters that, a strong resonance exists in SQ whereas disparity in carbon-carbon bond lengths of central C₄ ring is seen in SMSQ dyes. The increasing and decreasing of antibonding interactions between central C₄ ring and side donor groups cause destabilization of HOMO and stabilization of LUMO respectively in case of SQ compared to SMSQ molecules. This leads to decreasing the HOMO-LUMO gap and promotes biradicaloid character of SQ. The absorption maxima obtained by using TD-DFT method with BLYP, B3LYP, BHandHLYP, CAM-B3LYP and M06-2X functionals are not in good agreement with experimental results. On the other hand SAC-CI method gives better results for all the molecules. From this work we can evolve a design principle of these molecules which play a role as sensitizers in dye sensitized solar cells.     

TD-DFT investigation of D–π–A organic dyes with thiophene moieties as π-spacers for use as sensitizers in DSSCs

The geometrical, conformational, and electronic properties of a series of D–π–A metal-free dyes designed for use as sensitizers in DSSCs were studied using DFT and TD-DFT methods. A substituted triphenylamine moiety was used as the donor group and 2-cyanoacrylic acid as the acceptor group in these dyes. They also contained conjugated bridging π-linker groups containing two or more thiophene rings to enhance the intramolecular charge transfer. The B3LYP, M06-HF, ωB97XD and CAM-B3LYP functionals were utilized in combination with the 6-31G(d,p) basis set for the calculations. The dye solvation process was taken into account via the polarizable continuum model. To rationalize the relationships between dye structure and the photochemical properties of the dyes when used as sensitizers in DSSCs, the vertical excitation energies, the light-harvesting efficiencies, the free-energy changes during the process of injecting an electron into the surface of a TiO₂ nanocrystalline semiconductor, and the open-circuit potentials were calculated for all of the dyes in the solvent THF using the above methods. The results of these computations are discussed and compared with the available corresponding experimental data.     

First hyperpolarizability of cyclooctatetraene modulated by alkali and alkaline earth metals

In the present investigation, the first hyperpolarizability of alkali and alkaline earth metal derivatives of cyclooctatetraene (COT) has been calculated using BHHLYP and CAM-B3LYP functional for 6-311++G(d,p), 6-311++G(3df,3pd), and aug-pc 2 basis sets. Introduction of Na/K atoms at the axial position of COT and Li, Na, K/Be, Mg, Ca metal atoms and cyanide groups at the equatorial sites leads to lager enhancement of first hyperpolarizability. The ring charge density can account for the variation of first hyperpolarizability. The two state model has been invoked to explain the variation of first hyperpolarizability.     

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β _μ ) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree–Fock, Møller–Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning’s correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β _μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β _μ (9H)/β _μ (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β _μ values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

Structure and electronic properties of azadirachtin

We performed a combined DFT and Monte Carlo ¹³C NMR chemical-shift study of azadirachtin A, a triterpenoid that acts as a natural insect antifeedant. A conformational search using a Monte Carlo technique based on the RM1 semiempirical method was carried out in order to establish its preferred structure. The B3LYP/6-311++G(d,p), wB97XD/6-311++G(d,p), M06/6-311++G(d,p), M06-2X/6-311++G(d,p), and CAM-B3LYP/6-311++G(d,p) levels of theory were used to predict NMR chemical shifts. A Monte Carlo population-weighted average spectrum was produced based on the predicted Boltzmann contributions. In general, good agreement between experimental and theoretical data was obtained using both methods, and the ¹³C NMR chemical shifts were predicted highly accurately. The geometry was optimized at the semiempirical level and used to calculate the NMR chemical shifts at the DFT level, and these shifts showed only minor deviations from those obtained following structural optimization at the DFT level, and incurred a much lower computational cost. The theoretical ultraviolet spectrum showed a maximum absorption peak that was mainly contributed by the tiglate group.     

Critical test of some computational methods for prediction of NMR <Superscript>1</Superscript>H and <Superscript>13</Superscript>C chemical shifts

Performance of 18 DFT functionals (B1B95, B3LYP, B3PW91, B97D, BHandHLYP, BMK, CAM-B3LYP, HSEh1PBE, M06-L, mPW1PW91, O3LYP, OLYP, OPBE, PBE1PBE, tHCTHhyb, TPSSh, wB97xD, VSXC) in combinations with six basis sets (cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, aug-cc-pVTZ, IGLO-II, and IGLO-III) and three methods for calculating magnetic shieldings (GIAO, CSGT, IGAIM) was tested for predicting ¹H and ¹³C chemical shifts for 25 organic compounds, for altogether 86 H and 88 C atoms. Proton shifts varied between 1.03 ppm to 12.00 ppm and carbon shifts between 7.87 ppm to 209.28 ppm. It was found that the best method for calculating ¹³C shifts is PBE1PBE/aug-cc-pVDZ with CSGT or IGAIM approaches (mae?=?1.66 ppm), for ¹H the best results were obtained with HSEh1PBE, mPW1PW91, PBE1PBE, CAM-B3LYP, and B3PW91 functionals with cc-pVTZ basis set and with CSGT or IGAIM approaches (mae?=?0.28 ppm). We found that often larger basis sets do not give better results for chemical shifts. The best basis sets for calculating ¹H and ¹³C chemical shifts were cc-pVTZ and aug-cc-pVDZ, respectively. CSGT and IGAIM NMR approaches can perform really well and are in most cases better than popular GIAO approach.

Open image in new window

Graphical Abstract Mean absolute errors for ¹H and ¹³C chemical shifts and computational times of neutral toluene molecule with aug-cc-pVDZ basis set and CSGT approach

     

Computational and spectroscopic data correlation study of <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>'-bisarylmalonamides (Part II)

To complement a previous UV study, we present a quantitative evaluation of substituent effects on spectroscopic data (¹H and ¹³C NMR chemical shifts as well as FT-IR absorption frequency) applied to N,N'-bisarylmalonamides, using simple and extended Hammett equations as well as the Swain-Lupton equation. Furthermore, the DFT CAM-B3LYP/6–311+G(d,p) method was applied to study the impact of different solvents on the geometry of the molecules and their spectral data. Additionally, experimental data are correlated with theoretical results; excellent linear dependence was obtained. The overall results presented in this paper show that N,N'-bisarylmalonamides are prominent candidates for model molecules.     

Assessment of the photosensitization properties of cationic porphyrins in interaction with DNA nucleotide pairs

We present a theoretical assessment of the photosensitization properties of meso-mono(N-methylpyridyl) triphenylporphyrin (1, MmPyP⁺), which interacts with DNA nucleotide pairs [adenine (A)-thymine (T); guanine (G)-cytosine (C)] via an external binding mode. The photosensitization properties of the arrangements 1A, 1T, 1G and 1C were investigated. A set of density functionals (B3LYP, PBE0, CAM-B3LYP, M06-2X, B97D) with the 6-31G(d) basis set was used to calculate the electronic absorption spectra in solution (water) following TD-DFT methodology. In all the arrangements, with the exception of 1C, the functional PBE0 produced the lowest deviation of the Soret band (0.1–0.2 eV). Using this functional, we show that the porphyrin–nucleotide interaction is stabilized, as reflected by a larger HOMO–LUMO gap than free porphyrin. A more important effect of the interaction corresponds to the red-shift of the Soret band of MmPyP⁺, which is in agreement with experimental results. This behavior could be explained by the higher symmetry found in arrangements with a lower dipole moment, and by the more symmetrical distribution of electronic density along the molecular orbitals, which provokes electronic transitions of lower energy. The structural model allowed us to show that MmPyP⁺ improves the characteristics as a photosensitizer when it interacts with nucleotide pairs due to the longer wavelength required for the Soret band. Results obtained for porphyrins with larger monocationic substituents (2, MmAP+; 3, MONPP+) do not lead to the same behavior. Although the structural model is insufficient to describe porphyrin photosensitization, it suggests that improvements in this property are produced by the inclusion of a cationic charge in the pyridyl ring and a smaller size of the substituent leading to a better communication in the porphyrin–nucleotide pair.     

Evolution of Hot Polaron States with a Nanosecond Lifetime in a Manganite Perovskite

Understanding and controlling the relaxation process of optically excited charge carriers in solids with strong correlations is of great interest in the quest for new strategies to exploit solar energy. Usually, optically excited electrons in a solid thermalize rapidly on a femtosecond to picosecond timescale due to interactions with other electrons and phonons. New mechanisms to slow down thermalization will thus be of great significance for efficient light energy conversion, e.g., in photovoltaic devices. Ultrafast optical pump–probe experiments in the manganite Pr_0.65Ca_0.35MnO_3, a photovoltaic, thermoelectric, and electrocatalytic material with strong polaronic correlations, reveal an ultraslow recombination dynamics on a nanosecond‐time scale. The nature of long living excitations is further elucidated by photovoltaic measurements, showing the presence of photodiffusion of excited electron–hole polaron pairs. Theoretical considerations suggest that the excited charge carriers are trapped in a hot polaron state. Escape from this state is possible via a slow dipole‐forbidden recombination process or via rare thermal fluctuations toward a conical intersection followed by a radiation‐less decay. The strong correlation between the excited polaron and the octahedral dynamics of its environment appears to be substantial for stabilizing the hot polaron.     

ONIOM DFT/PM3 calculation on the interaction between STI-571 and abelson tyrosine kinase

The ONIOM2 (B3LYP/6–31G (d, p): PM3) and B3LYP/6–31G (d, p) methods were applied to investigate the interaction between STI-571 and abelson tyrosine kinase binding site. The complex of N-[4-methyl-3-(4-pyridin-3-yl-pyrimidin-2-ylamino)- phenyl]-benzamide (part of STI-571) and related 16 amino acid residues were found at B3LYP/6–31G (d, p) level to have hydrogen bonds and π....π stacking interaction, their binding energy via HAF optimization was −20.4 kcal mol⁻¹. The results derived from this study agreed well with the reported observation. Figure Optimized structure of STI-571 and Thr315 in abelson tyrosine kinase based on ONIOM2 method