Moment expansion of the linear density‐density response function

We present a low rank moment expansion of the linear density‐density response function. The general interacting (fully nonlocal) density‐density response function is calculated by means of its spectral decomposition via an iterative Lanczos diagonalization technique within linear density functional perturbation theory. We derive a unitary transformation in the space of the eigenfunctions yielding subspaces with well‐defined moments. This transformation generates the irreducible representations of the density‐density response function with respect to rotations within SO(3). This allows to separate the contributions to the electronic response density from different multipole moments of the perturbation. Our representation maximally condenses the physically relevant information of the density‐density response function required for intermolecular interactions, yielding a considerable reduction in dimensionality. We illustrate the performance and accuracy of our scheme by computing the electronic response density of a water molecule to a complex interaction potential. © 2015 Wiley Periodicals, Inc.     

First‐order correlation‐kinetic contribution to Kohn–Sham exchange charge density function in atoms,using quantal density functional theory approach

Using the static exchange‐correlation charge density concept, the total integrated exchange‐charge density function is calculated within the nonrelativistic spin‐restricted exchange‐only (i) optimized effective potential model, and (ii) nonvariational local potential derived from the exchange‐only work potential within the quantal density functional theory, for the ground‐state isoelectronic series: Ga⁺, Zn, Cu^?; In⁺, Cd, Ag^?; and Tl⁺, Hg, Au^?. The difference between the exchange charge density function derived from these potentials is employed to evaluate the first‐order correlation‐kinetic contribution to the integrated exchange charge density. This contribution is found to be important for both the intra‐ and inter‐shell regions. Screening effects on the contribution due to the nd¹⁰ (n = 3–5) subshells are discussed through comparisons with similar calculations on Ca, Sr, and Ba, wherein nd¹⁰ electrons are absent. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Link between the kinetic‐ and exchange‐energy functionals in the generalized gradient approximation

An approximate kinetic‐energy functional of the generalized gradient approximation form was derived following the “conjointness conjecture” of Lee, Lee, and Parr. The functional shares the analytical form of its gradient dependency with the exchange‐energy functionals of Becke and Perdew, Burke, and Ernzerhof. The two free parameters of this functional were determined using the exact values of the kinetic energy of He and Xe atoms. A set of 12 closed‐shell atoms was used to test the accuracy of the proposed functional and more than 30 others taken from the literature. It is shown that the conjointness conjecture leads to a very good class of kinetic‐energy functionals. Moreover, the functional developed in this work is shown to be one of the most accurate despite its simple analytical form. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Hosoya polynomials of armchair open‐ended nanotubes

For a connected graph G we denote by d(G,k) the number of vertex pairs at distance k. The Hosoya polynomial of G is H(G,x) = ∑_k≥0 d(G,k)x^k. In this paper, analytical formulae for calculating the polynomials of armchair open‐ended nanotubes are given. Furthermore, the Wiener index, derived from the first derivative of the Hosoya polynomial in x = 1, and the hyper‐Wiener index, from one‐half of the second derivative of the Hosoya polynomial multiplied by x in x = 1, can be calculated. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Correlation energy,correlated electron density,and exchange‐correlation potential in some spherically confined atoms

We report correlation energies, electron densities, and exchange‐correlation potentials obtained from configuration interaction and density functional calculations on spherically confined He, Be, Be²⁺, and Ne atoms. The variation of the correlation energy with the confinement radius R_c is relatively small for the He, Be²⁺, and Ne systems. Curiously, the Lee–Yang–Parr (LYP) functional works well for weak confinements but fails completely for small R_c. However, in the neutral beryllium atom the CI correlation energy increases markedly with decreasing R_c. This effect is less pronounced at the density‐functional theory level. The LYP functional performs very well for the unconfined Be atom, but fails badly for small R_c. The standard exchange‐correlation potentials exhibit significant deviation from the “exact” potential obtained by inversion of Kohn–Sham equation. The LYP correlation potential behaves erratically at strong confinements. © 2016 Wiley Periodicals, Inc.     

An assessment of pure,hybrid, meta,and hybrid‐meta GGA density functional theory methods for open‐shell systems: The case of the nonheme iron enzyme 8R–LOX

The performance of a range density functional theory functionals combined in a quantum mechanical (QM)/molecular mechanical (MM) approach was investigated in their ability to reliably provide geometries, electronic distributions, and relative energies of a multicentered open‐shell mechanistic intermediate in the mechanism 8R–Lipoxygenase. With the use of large QM/MM active site chemical models, the smallest average differences in geometries between the catalytically relevant quartet and sextet complexes were obtained with the B3LYP^* functional. Moreover, in the case of the relative energies between ⁴II and ⁶II , the use of the B3LYP^* functional provided a difference of 0.0 kcal mol^–1. However, B3LYP^± and B3LYP also predicted differences in energies of less than 1 kcal mol^–1. In the case of describing the electronic distribution (i.e., spin density), the B3LYP^*, B3LYP, or M06‐L functionals appeared to be the most suitable. Overall, the results obtained suggest that for systems with multiple centers having unpaired electrons, the B3LYP^* appears most well rounded to provide reliable geometries, electronic structures, and relative energies. © 2012 Wiley Periodicals, Inc.     

Combined experimental and computational modeling studies on 4‐[(2‐hydroxy‐3‐methylbenzylidene) amino]‐1,5‐dimethyl‐2‐phenyl‐1,2‐dihydro‐3H‐pyrazol‐3‐one

The Schiff base compound, 4‐[(2‐hydroxy‐3‐methylbenzylidene)amino]‐1,5‐dimethyl‐2‐phenyl‐1,2‐dihydro‐3H‐pyrazol‐3‐one, has been synthesized and characterized by IR, UV–vis, and X‐ray single‐crystal determination. Molecular geometry from X‐ray experiment of the title compound in the ground state have been compared using the density functional method (B3LYP) with 6‐31G(d,p) basis set. Calculated results show that density functional theory (DFT) can well reproduce the structure of the title compound. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6‐31G(d,p) basis set by applying the Onsager and the polarizable continuum model (PCM). The results obtained with these methods reveal that the PCM method provided more stable structure than Onsager's method. By using TD‐DFT method, electronic absorption spectra of the title compound have been predicted and a good agreement with the TD‐DFT method and the experimental one is determined. The predicted nonlinear optical properties of the title compound are much greater than ones of urea. In addition, DFT calculations of the title compound, molecular electrostatic potential and NBO analysis were performed at B3LYP/6‐31G(d,p) level of theory. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

A solvation‐free‐energy functional: A reference‐modified density functional formulation

The three‐dimensional reference interaction site model (3D‐RISM) theory, which is one of the most applicable integral equation theories for molecular liquids, overestimates the absolute values of solvation‐free‐energy (SFE) for large solute molecules in water. To improve the free‐energy density functional for the SFE of solute molecules, we propose a reference‐modified density functional theory (RMDFT) that is a general theoretical approach to construct the free‐energy density functional systematically. In the RMDFT formulation, hard‐sphere (HS) fluids are introduced as the reference system instead of an ideal polyatomic molecular gas, which has been regarded as the appropriate reference system of the interaction‐site‐model density functional theory for polyatomic molecular fluids. We show that using RMDFT with a reference HS system can significantly improve the absolute values of the SFE for a set of neutral amino acid side‐chain analogues as well as for 504 small organic molecules. © 2015 Wiley Periodicals, Inc.     

Domain‐Averaged Exchange‐Correlation Energies as a Physical Underpinning for Chemical Graphs

A novel solution to the problem of assigning a molecular graph to a collection of nuclei (i.e. how to draw a molecular structure) is presented. Molecules are universally understood as a set of nuclei linked by bonds, but establishing which nuclei are bonded and which are not is still an empirical matter. Our approach borrows techniques from quantum chemical topology, which showed for the first time the construction of chemical graphs from wave functions, shifting the focus on energetics. This new focus resolves issues surrounding previous topological analyses, in which domain‐averaged exchange‐correlation energies (V_xc), quantities defined in real space between each possible atom pair, hold the key. Exponential decay of V_xc in non‐metallic systems as the intercenter distance increases guarantees a well‐defined hierarchy for all possible V_xc values in a molecule. Herein, we show that extracting the set of atom pairs that display the largest V_xc values in the hierarchy is equivalent to retrieving the molecular graph itself. Notably, domain‐averaged exchange‐correlation energies are transferable, and they can be used to calculate bond strengths. Fine‐grained details resulted to be related to simple stereoelectronic effects. These ideas are demonstrated in a set of simple pilot molecules.     

An open‐source framework for analyzing N‐electron dynamics. II. Hybrid density functional theory/configuration interaction methodology

In this contribution, we extend our framework for analyzing and visualizing correlated many‐electron dynamics to non‐variational, highly scalable electronic structure method. Specifically, an explicitly time‐dependent electronic wave packet is written as a linear combination of N‐electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time‐dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open‐source Python program det CI@ORBKIT, which extends the capabilities of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one‐electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser‐driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher‐level method provided a judicious choice of functional is made. Broadband excitation of a medium‐sized organic chromophore further demonstrates the scalability of the method. In addition, the time‐dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance. © 2017 Wiley Periodicals, Inc.     

Semiconductor nanocrystals possessing broadly size‐ and facet‐dependent optical properties

Cu₂O cubes, octahedra, and rhombic dodecahedra have been shown to exhibit continuous light absorption and emission band shifts with increasing particle sizes from 10 nm to sub‐microcrystals. They also possess clear facet‐dependent optical properties. Ag₃PO₄, Ag₂O, SrTiO₃, and CeO₂ crystals show similar optical size and facet effects. Thus, spectral shifts over a broad size range far beyond the quantum‐size regime should be generally observable in many semiconductor materials. Facet‐dependent optical properties of a semiconductor can be understood to arise from the presence of an ultrathin surface layer with subtle bond and orbital level variations for different crystal faces. Although these optical features seem unexpected, they should be the general behaviors of semiconductor crystals. As more examples of these optical effects are available, we will find that these intrinsic properties of semiconductors have been ignored in the past. Furthermore, if valence and conduction band positions are broadly tunable by particle size, the knowledge should have tremendous impacts on the applications of semiconductors, where band energies are important to efficient interfacial charge transfer.     

Analytical nuclear excited‐state gradients for the Tamm–Dancoff approximation using uncoupled frozen‐density embedding

We report the derivation and implementation of analytical nuclear gradients for excited states using time‐dependent density functional theory using the Tamm–Dancoff approximation combined with uncoupled frozen‐density embedding using density fitting. Explicit equations are presented and discussed. The implementation is able to treat singlet as well as triplet states and functionals using the local density approximation, the generalized gradient approximation, combinations with Hartree–Fock exchange (hybrids), and range‐separated functionals such as CAM‐B3LYP. The new method is benchmarked against supermolecule calculations in two case studies: The solvatochromic shift of the (vertical) fluorescence energy of 4‐aminophthalimide on solvation, and the first local excitation of the benzonitrile dimer. Whereas for the 4‐aminophthalimide–water complex deviations of about 0.2 eV are obtained to supermolecular calculations, for the benzonitrile dimer the maximum error for adiabatic excitation energies is below 0.01 eV due to a weak coupling of the subsystems. © 2017 Wiley Periodicals, Inc.     

A first‐principles study of diatomic NiAl: Ground state,structure, and spectroscopic constants

A computational study of diatomic NiAl is reported. Molecular properties evaluated include the equilibrium bond length (r_e), equilibrium stretching frequency (ω_e), doublet‐quartet energy splitting, and nickel‐aluminum bond strength. Several interesting conclusions have resulted from this research. First, convergence in calculated properties is smoother with recently reported correlation consistent basis sets than earlier basis sets for Ni and Al. Second, with the exception of bond strength, basis set limit properties extrapolated using correlation basis sets are in agreement with reported data. Third, this research suggests that caution may be needed with regard to the use of DFT for developing interatomic potentials for larger scale simulations. For example, B97‐1 showed better agreement with reported r_e for ²NiAl than B3LYP. However, the situation was reversed for the calculation of ω_e. With respect to bond strength, the situation is unclear due to the scatter among experiment and calculations. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Evaluation of the van Deemter equation in terms of open‐ended flow to chromatography

Plate theory and adsorption theory are the main tools available for understanding chromatographic experiments. Both theories predict a Gaussian distribution of solute molecules within the tubular system. However, Gaussian concentration distributions are observed predominantly at slow linear flow rates, while asymmetric concentration distributions are observed at the linear flow rates most used in chromatography. Allegedly, this asymmetry originates from an inhomogeneous distribution of grain sizes in the column and column overload. However, it is an experimental fact that the distribution of chemicals within an injected volume of solute changes as a function of time, while the response is measured simultaneously. Accordingly, the obtained signal cannot be compared to the theory before some type of time‐based deconvolution of the data has been performed. Adjustments to high‐performance liquid chromatography data were thus proposed through empirical equations that describe the relevant time values, peak height, peak area, and parameters of the van Deemter equation. It was proposed that the transfer of solute from the front to the rear part of the pulse during laminar open‐ended flow occurs at rate that depends on the linear flow rate, and to a lesser extent, on properties of the response function.     

TD‐DFT benchmarks: A review

Time‐Dependent Density Functional Theory (TD‐DFT) has become the most widely‐used theoretical approach to simulate the optical properties of both organic and inorganic molecules. In this contribution, we review TD‐DFT benchmarks that have been performed during the last decade. The aim is often to pinpoint the most accurate or adequate exchange‐correlation functional(s). We present both the different strategies used to assess the functionals and the main results obtained in terms of accuracy. In particular, we discuss both vertical and adiabatic benchmarks and comparisons with both experimental and theoretical reference transition energies. More specific benchmarks (oscillator strengths, excited‐state geometries, dipole moments, vibronic shapes, etc.) are summarized as well. © 2013 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.     

Accelerating Kohn‐Sham response theory using density fitting and the auxiliary‐density‐matrix method

An extension of the formulation of the atomic‐orbital‐based response theory of Larsen et al., JCP 113, 8909 (2000) is presented. This new framework has been implemented in LSDalton and allows for the use of Kohn‐Sham density‐functional theory with approximate treatment of the Coulomb and Exchange contributions to the response equations via the popular resolution‐of‐the‐identity approximation as well as the auxiliary‐density matrix method (ADMM). We present benchmark calculations of ground‐state energies as well as the linear and quadratic response properties: vertical excitation energies, polarizabilities, and hyperpolarizabilities. The quality of these approximations in a range of basis sets is assessed against reference calculations in a large aug‐pcseg‐4 basis. Our results confirm that density fitting of the Coulomb contribution can be used without hesitation for all the studied properties. The ADMM treatment of exchange is shown to yield high accuracy for ground‐state and excitation energies, whereas for polarizabilities and hyperpolarizabilities the performance gain comes at a cost of accuracy. Excitation energies of a tetrameric model consisting of units of the P700 special pigment of photosystem I have been studied to demonstrate the applicability of the code for a large system.     

Orbital functional theory of linear response and excitation

Orbital functional theory (OFT) is based on a rule that determines a single‐determinant reference state Φ for any exact N‐electron eigenstate Ψ. An OFT model postulates an explicit correlation energy functional E_c of occupied orbital functions {?_i} and occupation numbers {n_i}. The orbital Euler–Lagrange equations are analogous to Kohn–Sham equations, but do not in general contain local potential functions. Time‐dependent Hartree–Fock theory is generalized in OFT to a formally exact linear response theory that includes electronic correlation. In the exchange‐only limit, the theory reduces to the random‐phase approximation of many‐body theory. The formalism determines excitation energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Methoxy and Methyl Group Rotation: Solid‐State NMR 1H Spin‐Lattice Relaxation,Electronic Structure Calculations,X‐ray Diffractometry,and Scanning Electron Microscopy

We report solid‐state ¹H nuclear magnetic resonance (NMR) spin‐lattice relaxation experiments, X‐ray diffractometry, field‐emission scanning electron microscopy, and both single‐molecule and cluster ab initio electronic structure calculations on 1‐methoxyphenanthrene ( 1 ) and 3‐methoxyphenanthrene ( 2 ) to investigate the rotation of the methoxy groups and their constituent methyl groups. The electronic structure calculations and the ¹H NMR relaxation measurements can be used together to determine barriers for the rotation of a methoxy group and its constituent methyl group and to develop models for the two coupled motions.