An open‐source framework for analyzing N‐electron dynamics. I. Multideterminantal wave functions

The aim of the present contribution is to provide a framework for analyzing and visualizing the correlated many‐electron dynamics of molecular systems, where an explicitly time‐dependent electronic wave packet is represented as a linear combination of N‐electron wave functions. The central quantity of interest is the electronic flux density, which contains all information about the transient electronic density, the associated phase, and their temporal evolution. It is computed from the associated one‐electron operator by reducing the multideterminantal, many‐electron wave packet using the Slater‐Condon rules. Here, we introduce a general tool for post‐processing multideterminant configuration‐interaction wave functions obtained at various levels of theory. It is tailored to extract directly the data from the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions. The procedure is implemented in the open‐source Python program det CI@ORBKIT, which shares and builds on the modular design of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). The new procedure is applied to ultrafast charge migration processes in different molecular systems, demonstrating its broad applicability. Convergence of the N‐electron dynamics with respect to the electronic structure theory level and basis set size is investigated. This provides an assessment of the robustness of qualitative and quantitative statements that can be made concerning dynamical features observed in charge migration simulations. © 2017 Wiley Periodicals, Inc.     

Theoretical Shaping of Femtosecond Laser Pulses for Molecular Photodissociation with Control Techniques Based on Ehrenfest′s Dynamics and Time‐Dependent Density Functional Theory

The combination of nonadiabatic Ehrenfest‐path molecular dynamics (EMD) based on time‐dependent density functional theory (TDDFT) and quantum optimal control formalism (QOCT) was used to optimize the shape of ultra‐short laser pulses to achieve photodissociation of a hydrogen molecule and the trihydrogen cation H₃⁺. This work completes a previous one [A. Castro, ChemPhysChem, 2013 , 14, 1488–1495], in which the same objective was achieved by demonstrating the combination of QOCT and TDDFT for many‐electron systems on static nuclear potentials. The optimization model, therefore, did not include the nuclear movement and the obtained dissociation mechanism could only be sequential: fast laser‐assisted electronic excitation to nonbonding states (during which the nuclei are considered to be static), followed by field‐free dissociation. Here, in contrast, the optimization was performed with the QOCT constructed on top of the full dynamic model comprised of both electrons and nuclei, as described within EMD based on TDDFT. This is the first numerical demonstration of an optimal control formalism for a hybrid quantum–classical model, that is, a molecular dynamics method.     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

Sub‐cycle Dynamics of High Harmonic Generation of Ne Atoms Excited by Attosecond Pulses and Driven by Near‐infrared Laser Fields: A Self‐Interaction‐Free TDDFT Approach

In the framework of the self‐interaction‐free time‐dependent density‐functional theory (TDDFT), we have performed three‐dimensional ab initio calculations of Ne atoms in near‐infrared (NIR) laser fields subject to excitation by a single extreme ultraviolet (XUV) attosecond pulse (SAP). The TDDFT equations are solved accurately and efficiently by means of the time‐dependent generalized pseudo spectral (TDGPS) method. We have explored the transient dynamical behavior of the sub‐cycle high harmonic generation (HHG) for transitions from the excited states to the ground state and found oscillation structures with respect to the time delay between the SAP and NIR fields. We investigate the harmonic emission spectrum from singly excited state 2p3s, 2p4s, 2p3d, 2p5s, 2p4d and 2p6s, 2p5d and the virtual states 2p3p‐, 2p4f‐ and 2p4p+ as the function of time delay. We explore the sub‐cycle Stark shift phenomenon in NIR fields and its influence on the photon emission process. Our analysis reveals several novel features of the sub‐cycle transient HHG dynamics and spectra, the quantum interference pattern between different multiphoton excitation pathways, and we identify the mechanisms responsible for the observed peak splitting in the photon emission spectra.     

Electronic structures and spectra of porphyrin with fused benzoheterocycles: DFT and TDDFT‐PCM investigations

Density functional theory (DFT) and time‐dependent DFT (TDDFT) are applied to study seven asymmetric π‐conjugated porphyrins with extended benzoheterocycles: quinoline, indole, benzoimidazole, benzothiazole, benzooxazole, 2,1,3‐benzothiadiazole, and 2,1,3‐benzoxadiazole. The solvation effects on the excitation energies for these porphyrin derivatives in chloroform are taken into account by using the continuum model (C‐PCM) combined with TDDFT, and this method makes a closer agreement with the experimental values, especially for the B‐bands of these objects. Great efforts have been made on investigating the influences of the fused aromatic units of the porphyrins on the absorption properties as these can be particularly important for many applications. Benzoheterocycle introduction and solvent effects have been systemically investigated, and close agreement is obtained between calculated and measured UV–vis spectra. These theoretical data could shed light on future synthetic chemistry. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

On low‐lying excited states of extended nanographenes

Low‐lying excited states of planarly extended nanographenes are investigated using the long‐range corrected (LC) density functional theory (DFT) and the spin‐flip (SF) time‐dependent density functional theory (TDDFT) by exploring the long‐range exchange and double‐excitation correlation effects on the excitation energies, band gaps, and exciton binding energies. Optimizing the geometries of the nanographenes indicates that the long‐range exchange interaction significantly improves the C C bond lengths and amplify their bond length alternations with overall shortening the bond lengths. The calculated TDDFT excitation energies show that long‐range exchange interaction is crucial to provide accurate excitation energies of small nanographenes and dominate the exciton binding energies in the excited states of nanographenes. It is, however, also found that the present long‐range correction may cause the overestimation of the excitation energy for the infinitely wide graphene due to the discrepancy between the calculated band gaps and vertical ionization potential (IP) minus electron affinity (EA) values. Contrasting to the long‐range exchange effects, the SF‐TDDFT calculations show that the double‐excitation correlation effects are negligible in the low‐lying excitations of nanographenes, although this effect is large in the lowest excitation of benzene molecule. It is, therefore, concluded that long‐range exchange interactions should be incorporated in TDDFT calculations to quantitatively investigate the excited states of graphenes, although TDDFT using a present LC functional may provide a considerable excitation energy for the infinitely wide graphene mainly due to the discrepancy between the calculated band gaps and IP–EA values. © 2017 Wiley Periodicals, Inc.     

Star-shaped Organic Molecules That Comprise a 1,3,5-Trisubs-tituted Benzene Core and Three Oligoaryleneethynylene Arms as Light-emitting Materials

Star‐shaped rigid molecules that comprise a 1,3,5‐trisubstitued benzene core and three oligoaryleneethynylene arms have great potential application in organic light‐emitting devices (OLEDs). Their optical and electronic properties are tuned by the star‐shaped molecular size. To reveal the relationship between the properties and structures, we perform a systemic investigation for these organic molecules. The ground and excited state molecules are studied using density functional theory (DFT), the ab initio HF, and the single excitation configuration interaction (CIS), respectively. And the electronic absorption and emission spectra are investigated with time‐dependent density functional theory (TDDFT) and Zerner's intermediate neglect of differential overlap (ZINDO) methods. The results show that the HOMOs, LUMOs, energy gaps, ionization potentials (IP), electron affinities (EA), absorption and emission spectra are controlled by the star‐shaped molecular size, which favor the hole and electron injection into OLEDs. With increasing the molecular conjugated length, the absorption and emission spectra exhibit red shifts to some extent and are in good agreement with the experimental ones. Also, the calculated emission spectra range from 330 to 440 nm. All the calculated show that the star‐shaped molecules are promising as blue light emitting materials     

Computation of magnetic circular dichroism by sum‐over‐states summations

Magnetic circular dichroism (MCD) spectroscopy has been established as a convenient method to study electronic structure, in particular for small symmetric organic molecules. Newer applications on more complex systems are additionally stimulated by the latest availability of precise quantum‐chemical techniques for the spectral simulations. In this work, a sum over states (SOS) summation is reexamined as an alternative to the derivative techniques for the MCD modeling. Unlike in previous works, the excited electronic states are calculated by the time‐dependent density functional theory (TDDFT). A gradient formulation of the MCD intensities is also proposed, less dependent on the origin choice than the standard expressions. The dependencies of the results on the basis set, number of electronic states, and coordinate origin are tested on model examples, including large symmetric molecules with degenerate electronic states. The results suggest that the SOS/TDDFT approach is a viable and accurate technique for spectral simulation. It may even considerably reduce the computational time, if compared with the traditional MCD computational procedures based on the response theory. © 2013 Wiley Periodicals, Inc.     

Time‐dependent density functional calculations on the electronic spectra of the neutral nickel complex [Ni(LISQ)2] (LISQ = 3,5‐di‐tert‐butyl‐o‐diiminobenzosemiquinonate(1−)) and its monoanion and dication

The electronic spectrum of the neutral nickel complex [Ni(L^ISQ)₂] (L^ISQ = 3,5‐di‐tert‐butyl‐o‐diiminobenzosemiquinonate(1^?)) and the spectra of its anion and dication have been calculated by means of time‐dependent density functional theory. The electronic ground state of the neutral complex exhibits an open shell singlet diradical character. The mandatory multireference problem for this electronic ground state has been treated approximately by using the unrestricted and spin symmetry broken Kohn‐Sham Slater determinant as the wave function for the noninteracting reference system in the time‐dependent density functional calculations. A reasonable agreement with observed transition energies and band intensities has been achieved. This holds also for the long wavelength transitions that are shown to be of charge transfer type. The charge distributions in the electronic ground state and the corresponding low lying excited states, however, are rather similar. Thus, the known failure of standard time‐dependent density functional theory to describe improperly long range charge transfer transitions is absent in this work. The applied computational scheme might be adequate for calculating electronic spectra of transition metal complexes with noninnocent ligands. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Structural and solvent effects on the spectroscopic properties of 1, 8‐naphthalimide derivatives: A density functional study

The molecular structures of 1, 8‐naphthalimide derivatives were investigated at density functional theory level within framework of PBE1PBE/6‐31G*. The vertical ionization potential and their delocalization energy of the X‐ray solid structure and gas‐phase optimized structure were explored. The configuration difference between them was attributed to the π‐π interaction of the solid effect, which has negligible effect on their absorption spectra. Solid effect also weakens the intramolecular interaction. Their absorption and luminescent spectra in gas and solvent phase were calculated by time‐dependent density functional theory (TDDFT) and conductor polarizable continuum models (CPCM)‐TDDFT approaches. Obvious red shifts from the solvent effect were found. Substituents on the imides will not improve their spectra properties a lot, whereas substituents on the naphthalene of naphthalimide would modify their properties to emit different spectra. Systematical deviation of vertical excitation energy from absorption and emission spectra, obtained by CPCM‐PBEPBE/6‐31G* and CIS‐CPCM‐PBEPBE/6‐31G* models, were about 0.05 eV and 0.02 eV compared with the experimental values. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

ERKALE—A flexible program package for X‐ray properties of atoms and molecules

ERKALE is a novel software program for computing X‐ray properties, such as ground‐state electron momentum densities, Compton profiles, and core and valence electron excitation spectra of atoms and molecules. The program operates at Hartree–Fock or density‐functional level of theory and supports Gaussian basis sets of arbitrary angular momentum and a wide variety of exchange‐correlation functionals. ERKALE includes modern convergence accelerators such as Broyden and ADIIS and it is suitable for general use, as calculations with thousands of basis functions can routinely be performed on desktop computers. Furthermore, ERKALE is written in an object oriented manner, making the code easy to understand and to extend to new properties while being ideal also for teaching purposes. © 2012 Wiley Periodicals, Inc.     

Analysis of self‐interaction correction for describing core excited states

Core‐excitation energies are calculated by the self‐interaction‐corrected time‐dependent density functional theory (SIC‐TDDFT) and SIC‐delta‐self‐consistent field (SIC‐ΔSCF) methods. For carbon monoxide, SIC‐TDDFT severely overestimates core‐excitation energies, while the SIC‐ΔSCF method using Kohn–Sham density functional theory (KS‐DFT) slightly overestimates. These behaviors are attributed to the fact that the self‐interaction errors in the total and orbital energies considerably differ. We evaluate the difference of the self‐interaction errors for the Slater exchange functional. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Time‐dependent density functional theory study on the electronic excited‐state hydrogen‐bonding dynamics of 4‐aminophthalimide (4AP) in aqueous solution: 4AP and 4AP–(H2O)1,2 clusters

The time‐dependent density functional theory (TDDFT) method has been carried out to investigate the excited‐state hydrogen‐bonding dynamics of 4‐aminophthalimide (4AP) in hydrogen‐donating water solvent. The infrared spectra of the hydrogen‐bonded solute?solvent complexes in electronically excited state have been calculated using the TDDFT method. We have demonstrated that the intermolecular hydrogen bond C? O···H? O and N? H···O? H in the hydrogen‐bonded 4AP?(H₂O)₂ trimer are significantly strengthened in the electronically excited state by theoretically monitoring the changes of the bond lengths of hydrogen bonds and hydrogen‐bonding groups in different electronic states. The hydrogen bonds strengthening in the electronically excited state are confirmed because the calculated stretching vibrational modes of the hydrogen bonding C?O, amino N? H, and H? O groups are markedly red‐shifted upon photoexcitation. The calculated results are consistent with the mechanism of the hydrogen bond strengthening in the electronically excited state, while contrast with mechanism of hydrogen bond cleavage. Furthermore, we believe that the transient hydrogen bond strengthening behavior in electroniclly excited state of chromophores in hydrogen‐donating solvents exists in many other systems in solution. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Time dependent wave packet treatment of 2ΠgN2− and 3Σ−NO− shape resonances using two‐dimensional surfaces for electron‐N2 and NO interactions

The ²Π_gN and ³Σ^?NO^? resonances in electron‐N₂ and NO collisions have been treated using both nuclear and electronic degrees of freedom and a two‐dimensional (2D) time dependent wave packet approach to ascertain the importance of nonlocality in electron–nuclear interaction. The results so obtained are compared with vibrational excitation cross‐sections obtained experimentally and those from other theoretical/numerical approaches using 1D local complex potential, 2D model with a combination of the exterior complex scaling method and a finite‐element implementation of the discrete‐variable representation. The results obtained provide detailed insight into the nuclear dynamics induced by electron–molecule collision and reveal that while for resonant excitation of lower vibrational modes, the nonlocal effect may not be as critical but importance of nonlocal effects may increase with increase in quanta of resonant vibrational excitation. © 2012 Wiley Periodicals, Inc.     

Electronic structure and spectroscopic properties of 6‐aminophenanthridine and its derivatives: Insights from density functional theory

6‐Aminophenanthridine (6AP) and its derivatives show important biological activities as antiprion compounds and inhibitors of the protein folding activity of the ribosome. Both of these activities depend on the RNA binding property of these compounds, which has been recently characterized by fluorescence spectroscopy. Hence, fundamental insights into the photophysical properties of 6AP compounds are highly important to understand their biological activities. In this work, we have calculated electronic structures and optical properties of 6AP and its three derivatives 6AP8CF₃, 6AP8Cl, and 6APi by density functional theory (DFT) and time‐dependent density functional theory (TDDFT). Our calculated spectra show a good agreement with the experimental absorption and fluorescence spectra, and thus, provide deep insights into the optical properties of the compounds. Furthermore, comparing the results obtained with four different hybrid functionals, we demonstrate that the accuracy of the functionals varies in the order B3LYP > PBE0 > M062X > M06HF. © 2015 Wiley Periodicals, Inc.     

Structural,electronic, and optical properties of phenol‐pyridyl boron complexes for light‐emitting diodes

The coordination chemistry of polydentate chelating ligands that contain mixed pyridinephenol donor sets has been a sought‐after target of study and is a possible extension to the chemistry of polypyridines. In this article, seven compounds, which are the four‐coordinate boron complexes containing the mixed phenol‐pyridyl group, have been studied by theoretical calculation. They can function as charge transport materials and emitters, with high efficiency and stability. To reveal the relationship between the structures and properties of these bifunctional or multifunctional electroluminescent materials, the ground and excited state geometries were optimized at the B3LYP/6‐31G(d), HF/6‐31G(d), and CIS/6‐31G(d) levels, respectively. The ionization potentials (IPs) and electron affinities (EAs) were computed. The mobilities of hole and electron in these compounds were studied computationally based on the Marcus electron transfer theory. The lowest excitation energies, and the maximum absorption and emission wavelengths of these compounds were calculated by time‐dependent density functional theory method. As a result of these calculations, the values of HOMO, LUMO, energy gaps, IPs, EAs, and the balance between the hole‐ and electron‐transfer are greatly improved with the substitution of carbazole in compound 6 . The calculated emission spectra of the seven studied molecules can almost cover the full UV‐vis range (from 447.4 to 649.3 nm). Also, the Stokes shifts are unexpectedly large, ranging from 139.4 to 335.1 nm. This will result in the relatively long fluorescence lifetimes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Relationship between orbital energy gaps and excitation energies for long‐chain systems

The difference between the excitation energies and corresponding orbital energy gaps, the exciton binding energy, is investigated based on time‐dependent (TD) density functional theory (DFT) for long‐chain systems: all‐trans polyacetylenes and linear oligoacenes. The optimized geometries of these systems indicate that bond length alternations significantly depend on long‐range exchange interactions. In TDDFT formalism, the exciton binding energy comes from the two‐electron interactions between occupied and unoccupied orbitals through the Coulomb‐exchange‐correlation integral kernels. TDDFT calculations show that the exciton binding energy is significant when long‐range exchange interactions are involved. Spin‐flip (SF) TDDFT calculations are then carried out to clarify double‐excitation effects in these excitation energies. The calculated SF‐TDDFT results indicate that double‐excitation effects significantly contribute to the excitations of long‐chain systems. The discrepancies between the vertical ionization potential minus electron affinity (IP–EA) values and the HOMO–LUMO excitation energies are also evaluated for the infinitely long polyacetylene and oligoacene using the least‐square fits to estimate the exciton binding energy of infinitely long systems. It is found that long‐range exchange interactions are required to give the exciton binding energy of the infinitely long systems. Consequently, it is concluded that long‐range exchange interactions neglected in many DFT calculations play a crucial role in the exciton binding energies of long‐chain systems, while double‐excitation correlation effects are also significant to hold the energy balance of the excitations. © 2016 Wiley Periodicals, Inc.     

Electronic structure and molecular orbital study of the first excited state of the high‐efficiency blue OLED material bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex from ab initio and TD‐B3LYP

Bis(2‐methyl‐8‐quinolinolato)aluminum(III) hydroxide complex (AlMq₂OH) is used in organic light‐emitting diodes (OLEDs) as an electron transport material and emitting layer. By means of ab initio Hartree–Fock (HF) and density functional theory (DFT) B3LYP methods, the structure of AlMq₂OH was optimized. The frontier molecular orbital characteristics and energy levels of AlMq₂OH have been analyzed systematically to study the electronic transition mechanism in AlMq₂OH. For comparison and calibration, bis(8‐quinolinolato)aluminum(III) hydroxide complex (Alq₂OH) has also been examined with these methods using the same basis sets. The lowest singlet excited state (S₁) of AlMq₂OH has been studied by the singles configuration interaction (CIS) method and time‐dependent DFT (TD‐DFT) using a hybrid functional, B3‐LYP, and the 6‐31G* basis set. The lowest singlet electronic transition (S₀ → S₁) of AlMq₂OH is π → π* electronic transitions and primarily localized on the different quinolate ligands. The emission of AlMq₂OH is due to the electron transitions from a phenoxide donor to a pyridyl acceptor from another quinolate ligand including C → C and O → N transference. Two possible electron transfer pathways are presented, one by carbon, oxygen, and nitrogen atoms and the other via metal cation Al³⁺. The comparison between the CIS‐optimized excited‐state structure with the HF ground‐state structure indicates that the geometric shift is mainly confined to the one quinolate and these changes can be easily understood in terms of the nodal patterns of the highest occupied and lowest unoccupied molecular orbitals. On the basis of the CIS‐optimized structure of the excited state, TD‐B3‐LYP calculations predict an emission wavelength of 499.78 nm. An absorption wavelength at 380.79 nm on the optimized structure of B3LYP/6‐31G* was predicted. They are comparable to AlMq2OH 485 and 390 nm observed experimentally for photoluminescence and UV‐vis absorption spectra of AlMq₂OH solid thin film on quartz, respectively. Lending theoretical corroboration to recent experimental observations and supposition, the reasons for the blue‐shift of AlMq2OH were revealed. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004