Design rules for the large two-photon absorption diketopyrrolopyrrole-based quadrupolar symmetrical chromophores

Two-photon absorption of a series of symmetrical diketopyrrolopyrrole (DPP) derivatives is studied by means of density functional theory applied to second-order response function. Several important issues in modeling are highlighted which must be addressed for a reliable reproduction of the experimental results. A comparison of the theoretical results with the experimental ones indicates that the computed two-photon absorption (TPA) cross sections show a good agreement for the three-state model in case of symmetrical quadrupolar DPP derivatives with donor ending groups. Although the results of TPA cross section observed for the systems with strong electron accepting ending groups seem to be overestimated for the three-state model. At the same time using an algorithm for the direct calculation of the third-order processes by the solving dynamic coupled-perturbed Hartree–Fock equation shows relatively good qualitative results for both donor and acceptor ending groups; however, it was found to be less accurate quantitatively for the already published experimental results. Thus, to obtain a valid overview, supposedly, both methods have to be considered. In general, in the present paper, a strategy toward molecular tailoring of large TPA cross-sectional materials is described, which can be extremely useful tool to predict materials’ properties prior to synthesis and measurements.     

Magnetic circular dichroism in real-time time-dependent density functional theory

We apply the adiabatic time-dependent density functional theory to magnetic circular dichroism (MCD) spectra using the real-space, real-time computational method. The standard formulas for the MCD response and its A and B terms are derived from the observables in the time-dependent wave function. We find real-time method is well suited for calculating the overall spectrum, particularly at higher excitation energies where individual excited states are numerous and overlapping. The MCD sum rules are derived and intepreted in the real-time formalism; we find that they are very useful for normalization purposes and assessing the accuracy of the theory. The method is applied to MCD spectrum of C(60) using the adiabatic energy functional from the local density approximation. The theory correctly predicts the signs of the A and B terms for the lowest allowed excitations. However, the magnitudes of the terms only show qualitative agreement with experiment.     

A perspective on nonresonant and resonant electronic response theory for time-dependent molecular properties

The development of electronic response theory in quantum chemistry has been reviewed, starting from the early 1970's and reaching the current state-of-the-art. The general theory has been applied to the calculation of a large number of spectroscopic parameters over the years, and it has been implemented for the majority of standard electronic structure methods. Two formulations of response theory, the Ehrenfest expectation value and the quasi-energy derivative formulation, have turned into leading alternatives for the derivation of computationally tractable expressions of response functions, and they are here reviewed with an attempt to, as far as possible, leave out technical details. A set of four steps are identified as common in derivations of response functions, and the two formulations are compared along this series of steps. Particular emphasis is given to the situation when the oscillation of the weak external electromagnetic field is in resonance with a transition frequency of the system. The formation of physically sound response functions in resonance regions of the spectrum is discussed in light of the causality condition and the Kramers-Kronig relations, and it is achieved in wave function theory by means of the introduction of relaxation parameters in a manner that mimics what one sees in density matrix theory. As a working example, equations are illustrated by their application to a two-state model for para-nitroaniline including the ground and the lowest charge-transfer state in the electric dipole approximation.     

噻吩基卟啉衍生物的单光子和双光子吸收性质的理论研究

本文理论上研究了两个系列的噻吩基卟啉衍生物,这种衍生物在可见光区具有大的双光子吸收截面。用密度泛函理论和ZINDO-SOS方法,计算了分子的几何构型、电子结构,单光子和双光子吸收性质。结果显示噻吩单元的数目影响分子的单光子和双光子吸收性质。具有两个或三个噻吩基团的噻吩基卟啉衍生物在较大范围内具有可用于实际应用中的双光子吸收响应,这一性质有利于这类分子在光限幅中的应用。插入乙炔基有利于扩大共轭范围,增加分子的双光子吸收截面。同时,乙炔基团的加入导致了单光子和双光子波长的红移。从高透明性和相对大的非线性光学响应考虑,噻吩基卟啉衍生物是一类有应用前景的双光子吸收材料。     

Two-photon absorption of hydrogen-bonded octupolar molecule clusters

Charge-transfer octupolar molecules can form clusters in solution through intermolecular hydrogen bonds. In the present work we explore the role of such clustering on two-photon absorption (TPA) spectra assuming 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as a model system. Using density functional quadratic response theory we examine different cluster structures of TATB dimers, trimers, and tetramers taken from snapshots of molecular dynamics simulations. In comparison with the TPA spectrum of a monomer, significant red shifts of charge-transfer states are predicted for all chosen clusters, which mainly is the result of the distortion of the structures induced by the aggregation. The TPA spectra for dimers and trimers show strong conformation dependence, whereas they turn out to be more stable for tetramers. Enhancements of TPA absorption have also been found for clusters containing less distorted molecules connected by hydrogen bonds.     

Two photon absorption of extended substituted phenylenevinylene oligomers: A TDDFT study

Time dependent density functional theory has been applied to analyze the effect of substituent groups on one-photon (OPA) and two-photon (TPA) absorption properties of extended amine-terminated phenylenevinylene oligomers. All investigated molecules are characterized by increased TPA activity. Calculated TPA cross-sections for these compounds do not show TPA enhancements observed in experimental measurements. Stabilization of the TPA active excited state in these chromophores leads to the small (0.2 eV) energy separation between the OPA and TPA allowed excited states, which agrees well with experimental data and may lead to the enhancement of the TPA response due to the strong vibronic couplings.     

Effects of conjugation in length and dimension on two-photon properties of fluorene-based chromophores

We report calculated two-photon (TPA) absorption spectra based upon results obtained from quadratic response time-dependent density functional theory for fluorene-based donor-π-acceptor molecules. Coulomb attenuated functionals with a long-range exchange contribution are applied to predict TPA excitation energies and cross-sections. Observed spectra are explained, and the effects of conjugation and multibranching on the TPA spectra are discussed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Contribution of the Mark S. Gordon 65th Birthday Festschrift Issue.     

Calculations of nonlinear response properties using the intermediate state representation and the algebraic-diagrammatic construction polarization propagator approach: two-photon absorption spectra

An earlier proposed approach to molecular response functions based on the intermediate state representation (ISR) of polarization propagator and algebraic-diagrammatic construction (ADC) approximations is for the first time employed for calculations of nonlinear response properties. The two-photon absorption (TPA) spectra are considered. The hierarchy of the first- and second-order ADC∕ISR computational schemes, ADC(1), ADC(2), ADC(2)-x, and ADC(3/2), is tested in applications to H(2)O, HF, and C(2)H(4) (ethylene). The calculated TPA spectra are compared with the results of coupled cluster (CC) models and time-dependent density-functional theory (TDDFT) calculations, using the results of the CC3 model as benchmarks. As a more realistic example, the TPA spectrum of C(8)H(10) (octatetraene) is calculated using the ADC(2)-x and ADC(2) methods. The results are compared with the results of TDDFT method and earlier calculations, as well as to the available experimental data. A prominent feature of octatetraene and other polyene molecules is the existence of low-lying excited states with increased double excitation character. We demonstrate that the two-photon absorption involving such states can be adequately studied using the ADC(2)-x scheme, explicitly accounting for interaction of doubly excited configurations. Observed peaks in the experimental TPA spectrum of octatetraene are assigned based on our calculations.     

Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Gyo?rffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.     

Calculation of two-photon absorption spectra of donor-pi-acceptor compounds in solution using quadratic response time-dependent density functional theory

Linear and quadratic response time-dependent density functional theories have been applied to calculate the photophysical properties of donor-pi-acceptor molecules which are known to have large nonlinear absorption. The linear absorption and two-photon absorption spectra predicted using hybrid functionals, including the Coulomb-attenuated model, with continuum solvation models are reported and compared to experiment and to previous theoretical predictions. While the quadratic response with these functionals overestimated the TPA cross sections relative to experiment when a Gaussian linewidth function was used, a fairly good agreement was obtained when a Lorentzian linewidth function was applied. In addition, the comparison of the TPA cross sections calculated by the sum over states with those calculated by the two-state approximation indicates the importance of the higher energy states in TPA, particularly in nondegenerate experiments.     

Auxiliary density perturbation theory

A new approach, named auxiliary density perturbation theory, for the calculation of second energy derivatives is presented. It is based on auxiliary density functional theory in which the Coulomb and exchange-correlation potentials are expressed by auxiliary function densities. Different to conventional coupled perturbed Kohn-Sham equations the perturbed density matrix is obtained noniteratively by solving an inhomogeneous equation system with the dimension of the auxiliary function set used to expand the auxiliary function density. A prototype implementation for the analytic calculation of molecular polarizabilities is presented. It is shown that the polarizabilities obtained with the newly developed auxiliary density perturbation approach match quantitative with the ones from standard density functional theory if augmented auxiliary function sets are used. The computational advantages of auxiliary density perturbation theory are discussed, too.     

Real-space, real-time calculation of dynamic hyperpolarizabilities

The finite-difference method to calculate hyperpolarizabilities is generalized for dynamical case. The calculation of the dynamical hyperpolarizabilities from non-perturbative, explicitly time-dependent single particle states obtained in the framework of the time-dependent density functional theory, is implemented in real space and real time. The optical response functions up to the third order are extracted in frequency domain. The present approach is free of deficiencies associated with atom centered basis sets and allows treatment of large molecules. The calculated results are in good agreement with experiments and with other theoretical calculations for various test cases.     

Turbo charging time-dependent density-functional theory with Lanczos chains

We introduce a new implementation of time-dependent density-functional theory which allows the entire spectrum of a molecule or extended system to be computed with a numerical effort comparable to that of a single standard ground-state calculation. This method is particularly well suited for large systems and/or large basis sets, such as plane waves or real-space grids. By using a superoperator formulation of linearized time-dependent density-functional theory, we first represent the dynamical polarizability of an interacting-electron system as an off-diagonal matrix element of the resolvent of the Liouvillian superoperator. One-electron operators and density matrices are treated using a representation borrowed from time-independent density-functional perturbation theory, which permits us to avoid the calculation of unoccupied Kohn-Sham orbitals. The resolvent of the Liouvillian is evaluated through a newly developed algorithm based on the nonsymmetric Lanczos method. Each step of the Lanczos recursion essentially requires twice as many operations as a single step of the iterative diagonalization of the unperturbed Kohn-Sham Hamiltonian. Suitable extrapolation of the Lanczos coefficients allows for a dramatic reduction of the number of Lanczos steps necessary to obtain well converged spectra, bringing such number down to hundreds (or a few thousands, at worst) in typical plane-wave pseudopotential applications. The resulting numerical workload is only a few times larger than that needed by a ground-state Kohn-Sham calculation for a same system. Our method is demonstrated with the calculation of the spectra of benzene, C(60) fullerene, and of chlorophyll a.     

Theoretical investigation of two-photon absorption allowed excited states in symmetrically substituted diacetylenes by ab initio molecular-orbital method

Symmetrically substituted diacetylene compounds, which shows large two-photon absorption (TPA) cross sections, have been theoretically investigated by the ab initio molecular-orbital method employing several theoretical models including the configuration interaction with single excitation (CIS), random phase approximation (RPA), and time-dependent density-functional theory (TDDFT) methods. The calculated excited energies are overestimated by CIS or RPA, whereas underestimated by TDDFT with the B3LYP parametrization for both one-photon absorption (OPA) and TPA allowed states. The lowest OPA state is well described by the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) transition. On the other hand, lower TPA allowed states can be represented as the superposition of the HOMO-LUMO+1 and HOMO-1-LUMO transitions, giving rise to two TPA allowed states. The absorption intensity for the lower TPA state of the diacetylenes molecules is discussed in terms of the alternancy symmetry and its breaking. The symmetry property is differently manifested for neutral and dicationic diacetylenes. Introduction of charges breaks the alternancy symmetry, which gives rise to an increase in the TPA cross sections at the lower frequency. The upper TPA state is calculated to show huge TPA cross sections, which reproduces the enhancement of the TPA cross section experimentally observed for one of the diacetylenes at the higher-frequency region. The enhancement is discussed employing an index defined as the ratio of the transition polarizability and its static limit, which represents the degree of influence of one-photon resonance on the TPA intensity. The huge TPA cross sections are found to be due to a near-resonance effect. The present theoretical calculation approves the previously proposed assumption based on the four-state (dual three-state) model, which consists of the ground, one OPA allowed, and two TPA allowed states.     

A new method for direct calculation of total energy of protein

A new scheme is developed for efficient quantum mechanical calculation of total energy of protein based on a recently developed MFCC (molecular fractionation with conjugate caps) approach. In this scheme, the linear-scaling MFCC method is first applied to calculate total electron density of protein. The computed electron density is then employed for direct numerical integration in density functional theory (DFT) to yield total energy of protein, with the kinetic energy obtained by a proposed ansatz. Numerical studies are carried out to calculate torsional energies of two polypeptides using this approach and the energies are shown to be in good agreement with the corresponding full system DFT calculation.     

A heuristic model of damped quantum rotation effects in nuclear magnetic resonance spectra

The damped quantum rotation (DQR) theory describes temperature effects in NMR spectra of hindered molecular rotators composed of identical atoms arranged in regular N-gons. In the standard approach, the relevant coherent dynamics are described quantum mechanically and the stochastic, thermally activated motions classically. The DQR theory is consistent. In place of random jumps over one, two, etc., maxima of the hindering potential, here one has damping processes of certain long-lived coherences between spin-space correlated eigenstates of the rotator. The damping-rate constants outnumber the classical jump-rate constants. The jump picture is recovered when the former cluster appropriately around only as many values as the number of the latter. The DQR theory was confirmed experimentally for hindered methyl groups in solids and even in liquids above 170 K. In this paper it is shown that for three-, four-, and sixfold rotators, the Liouville space equations of NMR line shapes, derived previously with the use of the quantum mechanical reduced density matrix approach, can be be given a heuristic justification. It is based on an equation of motion for the effective spin density matrix, where the relevant spin hamiltonian contains randomly fluctuating terms. The occurrence of the latter can be rationalized in terms of fluctuations of the tunneling splittings between the torsional sublevels of the rotator, including momentary liftings of the Kramers degeneracies. The question whether such degeneracy liftings are physical or virtual is discussed. The random terms in the effective hamiltonian can be Monte Carlo modeled as piecewise constant in time, which affords the stochastic equation of motion to be solved numerically in the Hilbert spin space. For sixfold rotators, this way of calculating the spectra can be useful in the instances where the Liouville space formalism of the original DQR theory is numerically unstable.     

Optical and Electrical Properties of Triphenylamine Derivatives for Dye-sensitized Solar Cells and Designing of Novel Molecule

With density functional theory(DFT) method, the optimization of molecular configurations and the calculation of frontier molecular orbitals were achieved for triphenylamine(TPA)-based dye-sensitized solar cell materials at the B3LYP/6-31G(d, p) level. Time-dependent density functional theory(TD-DFT) was applied to calculating the probability of the transition from the ground state to the excited state. And UV-Vis absorption spectra were derived with Franck-Condon approximation. The conjugation length, substitution groups and spatial effects show a slight influence on the dihedral angle of the TPA group. The increase of conjugation length may cause a smaller energy gap as well as a higher highest occupied molecular orbital(HOMO) and a lower lowest unoccupied molecular orbital (LUMO). The introduction of methoxyl group and TPA group could lower the energy gap while the HOMO and LUMO were elevated in energy.     

Average excitation energies from time-dependent density functional response theory

The authors present an occupation number averaging scheme for time-dependent density functional response theory (TD-DFRT) in frequency domain. The known problem that TD-DFRT within the local (spin) density approximation (LDA/LSDA) inaccurately predicts Rydberg and charge-transfer excitation energies has been reexamined from the methodology of linear response, without explicit correction of the exchange-correlation potential. The working equations of TD-DFRT are adapted to treat arbitrary difference of orbital occupation numbers, using the nonsymmetric matrix form of Casida's formulation of TD-DFRT [M. E. Casida, in Recent Advances in Density Functional Methods, edited by D. P. Chong (World Scientific, Singapore, 1995), Pt. I, p. 155]. The authors' scheme is applied to typical closed-shell and open-shell molecular systems by examining the dependence of excitation energies on the fraction of excited electron. Good performance of this modified linear response scheme is shown, and is consistent with the authors' previous examination by the real-time propagation approach, suggesting that the calculation of average excitation energies might be one of the ways to better decode excitation energies from LDA/LSDA. Different techniques for treating singlet, triplet, and doublet states are discussed.     

Numerical solution of Dalgarno-Lewis equations by a mapped Fourier grid method

Inhomogeneous radial differential equations emerging in applications of standard perturbation theory are numerically solved by a novel approach making use of Fourier grid methods in conjunction with a simple mapping scheme. The proposed algorithm is applied along the lines of the Dalgarno-Lewis method [Proc. R. Soc. London, Ser. A 223, 70 (1955)] to the calculation of the static dipole polarizabilities and hyperpolarizabilities of 1s, 2s, and 2p states of hydrogen atom and their frequency dependent dynamic dipole polarizabilities. The high efficiency and accuracy of the algorithm are demonstrated for the above test cases, where exact values are available. Then, the frequency dependent dipole polarizability of the ground state of lithium atom is computed by a variationally stable method combining an effective local potential approach with a second-order energy correction. The obtained results are in perfect agreement with other elaborate theoretical approaches.