Time-dependent density functional theory with the generalized restricted-unrestricted approach

Exploration of charge transfer analysis and photovoltaics properties of A-D-A type non-fullerene phenazine based molecules to enhance the organic solar cell properties

To intensify the photovoltaic properties of organic solar cells, density functional theory (DFT) based computational techniques were implemented on six non-fullerene A-D-A type small molecules (N1-N6) modified from reference molecule (R) which consists of phenazine fused with 1,4- Dimethyl-4H-3,7-dithia-4-aza- cyclopenta [α] pentalene on both sides with one of its phenyl rings acting as the central donor unit, further attached with 2-(5,6-Difluoro-2-methylene-3-oxo-indan-1-ylidene)-malononitrile acceptor groups at terminal sites. All proposed compounds have a phenazine base modified with a variety of substituents at the terminals. Transition density matrix, density of states, frontier molecular orbitals, intramolecular charge transfer abilities and optoelectronic properties of these compounds were investigated using B3LYP/6-31G (d, p) and B3LYP/6-31G++ (d,p) level of theory. All six designed compounds exhibited a bathochromic sift in their λmax as compared to the R molecule. All designed molecules also have reduced band gap and smaller excitation energy than R. Among all, N6 exhibited highest λmax and lowest bandgap as compared to reference molecule indicating its promising photovoltaic properties. Decreased hole and electron reorganization energy in several of the suggested compounds is indicative of greater charge mobility in them. PTB7-Th donor was employed to calculate open circuit voltage of all investigated molecules. N1-N5 molecules had improved optoelectronic properties, significant probable power conversion efficiency as evident from their absorption aspects, high values of Voc, and fill factor, compared to R molecule. Designed A-D-A type NF based molecules make OSCs ideal for use in wearable devices, building-integrated photovoltaics and smart fabrics.

Efficient side-chain engineering of thieno-imidazole salt-based molecule to boost the optoelectronic attributes of organic solar cells: A DFT approach

This study focused on modeling and density functional theory (DFT) analysis of reference (AI1) and designed structures (AI11-AI15), based on the thieno-imidazole core, in order to create profitable candidates for solar cells. All the optoelectronic properties of the molecular geometries were computed using DFT and time dependent-DFT approaches. The influence of terminal acceptors on the bandgaps, absorption, hole and electron mobilities, charge transfer capabilities, fill factor, dipole moment, etc. Of the recently designed structures (AI11-AI15), as well as reference (AI1), were evaluated. Optoelectronics and chemical parameters of newly architecture geometries were shown to be superior to the cited molecule. The FMOs and DOS graphs also demonstrated that the linked acceptors remarkably improved the dispersion of charge density in the geometries under study, particularly in AI11 and AI14. Calculated values of binding energy and chemical potential confirmed the thermal stability of the molecules. All the derived geometries surpassed the AI1 (Reference) molecule in terms of maximum absorbance ranging from 492 to 532 nm (in chlorobenzene solvent) and a narrower bandgap ranging from 1.76 to 1.99eV. AI15 had the lowest exciton dissociation energy of 0.22eV as well as lowest electrons and hole dissociation energies, while AI11 and AI14 showed highest VOC, fill factor, power conversion efficiency (PCE), IP and EA (owing to presence of strong electron pulling cyano (CN) moieties at their acceptor portions and extended conjugation) than all the examined molecules, implying that they could be used to build elite solar cells with enhanced photovoltaic attributes.

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

The Dalton quantum chemistry program system

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Role of zero-point vibrational corrections to carbon hyperfine coupling constants in organic π radicals

By analyzing a set of organic π radicals, we demonstrate that zero-point vibrational corrections give significant contributions to carbon hyperfine coupling constants, in one case even inducing a sign reversal for the coupling constant. We discuss the implications of these findings for the computational analysis of electron paramagnetic spectra based on hyperfine coupling constants evaluated at the equilibrium geometry of radicals. In particular, we note that a dynamical description that involves the nuclear motion is in many cases necessary in order to achieve a semi-quantitatively predictive theory for carbon hyperfine coupling constants. In addition, we discuss the implications of the strong dependence of the carbon hyperfine coupling constants on the zero-point vibrational corrections for the selection of exchange-correlation functionals in density functional theory studies of these constants.

Spin-adapted open-shell random phase approximation and time-dependent density functional theory. I. Theory

The spin-adaptation of single-reference quantum chemical methods for excited states of open-shell systems has been nontrivial. The primary reason is that the configuration space, generated by a truncated rank of excitations from only one component of a reference multiplet, is spin-incomplete. Those "missing" configurations are of higher ranks and can, in principle, be recaptured by a particular class of excitation operators. However, the resulting formalisms are then quite involved and there are situations [e.g., time-dependent density functional theory (TD-DFT) under the adiabatic approximation] that prevent one from doing so. To solve this issue, we propose here a tensor-coupling scheme that invokes all the components of a reference multiplet (i.e., a tensor reference) rather than increases the excitation ranks. A minimal spin-adapted n-tuply excited configuration space can readily be constructed by tensor products between the n-tuple tensor excitation operators and the chosen tensor reference. Further combined with the tensor equation-of-motion formalism, very compact expressions for excitation energies can be obtained. As a first application of this general idea, a spin-adapted open-shell random phase approximation is first developed. The so-called "translation rule" is then adopted to formulate a spin-adapted, restricted open-shell Kohn-Sham (ROKS)-based TD-DFT (ROKS-TD-DFT). Here, a particular symmetry structure has to be imposed on the exchange-correlation kernel. While the standard ROKS-TD-DFT can access only excited states due to singlet-coupled single excitations, i.e., only some of the singly excited states of the same spin (S(i)) as the reference, the new scheme can capture all the excited states of spin S(i)-1, S(i), or S(i)+1 due to both singlet- and triplet-coupled single excitations. The actual implementation and computation are very much like the (spin-contaminated) unrestricted Kohn-Sham-based TD-DFT. It is also shown that spin-contaminated spin-flip configuration interaction approaches can easily be spin-adapted via the tensor-coupling scheme.

Zero-point vibrational corrections to isotropic hyperfine coupling constants in polyatomic molecules

The present work addresses isotropic hyperfine coupling constants in polyatomic systems with a particular emphasis on a largely neglected, but a posteriori significant, effect, namely zero-point vibrational corrections. Using the density functional restricted-unrestricted approach, the zero-point vibrational corrections are evaluated for the allyl radical and four of its derivatives. In addition for establishing the numerical size of the zero-point vibrational corrections to the isotropic hyperfine coupling constants, we present simple guidelines useful for identifying hydrogens for which such corrections are significant. Based on our findings, we critically re-examine the computational procedures used for the determination of hyperfine coupling constants in general as well as the practice of using experimental hyperfine coupling constants as reference data when benchmarking and optimizing exchange-correlation functionals and basis sets for such calculations.

Density functional restricted-unrestricted approach for nonlinear properties: application to electron paramagnetic resonance parameters of square planar copper complexes

The density functional restricted-unrestricted approach for treatments of spin polarization effects in molecular properties using spin restricted Kohn-Sham theory has been extended from linear to nonlinear properties. It is shown that the spin polarization contribution to a nonlinear property has the form of a quadratic response function that includes the zero-order Kohn-Sham operator, in analogy to the lower order case where the spin polarization correction to an expectation value has the form of a linear response function. The developed approach is used to formulate new schemes for computation of electronic g-tensors and hyperfine coupling constants, which include spin polarization effects within the framework of spin restricted Kohn-Sham theory. The proposed computational schemes are in the present work employed to study the spin polarization effects on electron paramagnetic resonance spin Hamiltonian parameters of square planar copper complexes. The obtained results indicate that spin polarization gives rise to sizable contributions to the hyperfine coupling tensor of copper in all investigated complexes, while the electronic g-tensors of these complexes are only marginally affected by spin polarization and other factors, such as choice of exchange-correlation functional or molecular structures, will have more pronounced impact on the accuracy of the results.

Time-dependent density functional theory for nonlinear properties of open-shell systems

This paper presents response theory based on a spin-restricted Kohn-Sham formalism for computation of time-dependent and time-independent nonlinear properties of molecules with a high spin ground state. The developed approach is capable to handle arbitrary perturbations and constitutes an efficient procedure for evaluation of electric, magnetic, and mixed properties. Apart from presenting the derivation of the proposed approach, we show results from illustrating calculations of static and dynamic hyperpolarizabilities of small Si(3n+1)H(6n+3) (n=0,1,2) clusters which mimic Si(111) surfaces with dangling bond defects. The results indicate that the first hyperpolarizability tensor components of Si(3n+1)H(6n+3) have an ordering compatible with the measurements of second harmonic generation in SiO2/Si(111) interfaces and, therefore, support the hypothesis that silicon surface defects with dangling bonds are responsible for this phenomenon. The results exhibit a strong dependence on the quality of basis set and exchange-correlation functional, showing that an appropriate set of diffuse functions is required for reliable predictions of the first hyperpolarizability of open-shell compounds.