Ab-Initio Calculation of Molecular Aggregation Effects: a Coumarin-343 Case Study

We present time-dependent density functional theory (TDDFT) calculations for single and dimerized Coumarin-343 molecules in order to investigate the quantum mechanical effects of chromophore aggregation in extended systems designed to function as a new generation of sensors and light-harvesting devices. Using the single-chromophore results, we describe the construction of effective Hamiltonians to predict the excitonic properties of aggregate systems. We compare the electronic coupling properties predicted by such effective Hamiltonians to those obtained from TDDFT calculations of dimers, and to the coupling predicted by the transition density cube (TDC) method. We determine the accuracy of the dipole-dipole approximation and TDC with respect to the separation distance and orientation of the dimers. In particular, we investigate the effects of including Coulomb coupling terms ignored in the typical tight-binding effective Hamiltonian. We also examine effects of orbital relaxation which cannot be captured by either of these models

Global hybrids from the semiclassical atom theory satisfying the local density linear response

We propose global hybrid approximations of the exchange-correlation (XC) energy functional which reproduce well the modified fourth-order gradient expansion of the exchange energy in the semiclassical limit of many-electron neutral atoms and recover the full local density approximation (LDA) linear response. These XC functionals represent the hybrid versions of the APBE functional [Phys. Rev. Lett. 106, 186406, (2011)] yet employing an additional correlation functional which uses the localization concept of the correlation energy density to improve the compatibility with the Hartree-Fock exchange as well as the coupling-constant-resolved XC potential energy. Broad energetical and structural testings, including thermochemistry and geometry, transition metal complexes, non-covalent interactions, gold clusters and small gold-molecule interfaces, as well as an analysis of the hybrid parameters, show that our construction is quite robust. In particular, our testing shows that the resulting hybrid, including 20\% of Hartree-Fock exchange and named hAPBE, performs remarkably well for a broad palette of systems and properties, being generally better than popular hybrids (PBE0 and B3LYP). Semi-empirical dispersion corrections are also provided.Comment: 12 pages, 4 figure

Unraveling the performance of dispersion-corrected functionals for the accurate description of weakly bound natural polyphenols

Long-range non-covalent interactions play a key role in the chemistry of natural polyphenols. We have previously proposed a description of supramolecular polyphenol complexes by the B3P86 density functional coupled with some corrections for dispersion. We couple here the B3P86 functional with the D3 correction for dispersion, assessing systematically the accuracy of the new B3P86-D3 model using for that the well-known S66, HB23, NCCE31, and S12L datasets for non-covalent interactions. Furthermore, the association energies of these complexes were carefully compared to those obtained by other dispersion-corrected functionals, such as B(3)LYP-D3, BP86-D3 or B3P86-NL. Finally, this set of models were also applied to a database composed of seven non-covalent polyphenol complexes of the most interest.FDM acknowledges financial support from the Swedish Research Council (Grant No. 621-2014-4646) and SNIC (Swedish National Infrastructure for Computing) for providing computer resources. The work in Limoges (IB and PT) is supported by the “Conseil Régional du Limousin”. PT gratefully acknowledges the support by the Operational Program Research and Development Fund (project CZ.1.05/2.1.00/03.0058 of the Ministry of Education, Youth and Sports of the Czech Republic). IB gratefully acknowledges financial support from “Association Djerbienne en France”

Software for the frontiers of quantum chemistry:An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design

Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited

We revisit our investigation of the diffusion Monte Carlo (DMC) simulation of p-DIB molecular crystal polymorphism. [J. Phys. Chem. Lett. 2010, 1, 1789-1794] We perform, for the first time, a rigorous study of finite-size effects and choice of nodal surface on the prediction of polymorph stability in molecular crystals using fixed-node DMC. Our calculations are the largest which are currently feasible using the resources of the K computer and provide insights into the formidable challenge of predicting such properties from first principles. In particular, we show that finite-size effects can influence the trial nodal surface of a small (1×1×1) simulation cell considerably. We therefore repeated our DMC simulations with a 1×3×3 simulation cell, which is the largest such calculation to date. We used a DFT nodal surface generated with the PBE functional and we accumulated statistical samples with ∼6.4×105 core-hours for each polymorph. Our final results predict a polymorph stability consistent with experiment, but indicate that results in our previous paper were somewhat fortuitous. We analyze the finite-size errors using model periodic Coulomb (MPC) interactions and kinetic energy corrections, according to the CCMH scheme of Chiesa, Ceperley, Martin, and Holzmann. We investigate the dependence of the finite-size errors on different aspect ratios of the simulation cell (k-mesh convergence) in order to understand how to choose an appropriate ratio for the DMC calculations. Even in the most expensive simulations currently possible, we show that the finite size errors in the DMC total energies are far larger than the energy difference between the two polymorphs, although error cancellation means that the polymorph prediction is accurate. Finally, we found that the T-move scheme is essential for these massive DMC simulations in order to circumvent population explosions and large time-step biases.Chemistry and Chemical Biolog

Implementation and optimization of DFT-D/COSab with respect to basis set and functional: Application to polar processes of furfural derivatives in solution

The implementation, optimization, and performance of DFT-D, including the effects of solvation, has been tested on applications of polar processes in solution, where dispersion and hydrogen bonding is known to be involved. Solvent effects are included using our ab initio continuum solvation strategy, COSab, a conductor-like continuum solvation model, modified for ab initio in the quantum chemistry program GAMESS. Structure and properties are investigated across various functionals to evaluate their ability to properly model dispersion and solvation effects. The commonly used S22 set with accurate interaction energies of organic complexes has been used for parametrization studies of dispersion parameters and relevant solvation parameters. Dunning’s correlation consistent basis sets, cc-pVnZ (n = D, T), are used in the optimization, together with the Grimme B97-D exchange-correlation functional. Both water (ε = 78.4) and ether (ε = 4.33) environments are considered. Optimized semiempirical dispersion correction parameters and solvent extent radii are proposed for several functionals. We find that special parametrization of the semiempirical dispersion correction when used together in the DFT-D/COSab approach is not necessary. The global performance is quite acceptable in terms of chemical accuracy and suggests that this approach is a reliable as well as economical method for evaluation of solvent effects in systems with dispersive interactions. The resulting theory is applied to a group of push−pull pyrrole systems to illustrate the effects of donor/acceptor and solvation on their conformational and energetic properties

Implementation and optimization of DFT-D with respect to basis set and functional for Study of Dispersion Interactions in Nanoscale Aromatic Hydrocarbons

The implementation, optimization, and performance of various DFT-D schemes have been tested on models for polar−π interactions between arenes spaced at van der Waals distances and on a series of functionalized corannulene derivatives and complexes. For DFT-D schemes involving a semiempirical correction, optimized parameters are proposed for several basis sets. Performance of the different DFT-D strategies is compared, where functionals include some of the most recently proposed, B97D, B2PLYP, BMK, and M06−2X functionals, together with several other well-known functionals. Semiempircally corrected dispersion functionals hold some promise as useful and affordable methods for studies involving large polynuclear aromatic molecules and molecules on metal surfaces

Ab initio quantum chemical computations of substituent effects on triaziridine strain energy and heat of formation

A computational investigation is carried out on the parent triaziridine and as a function of N-substituents. Assessment of heat of formation, ring strain energy, barriers to inversion of nitrogen, and NMR criteria leads to understanding of issues related to vicinal lone pair repulsion and aromatic stabilization. Results lead to the proposal of a potentially flat structure with a aromatic-like triaziridine system, N3(BH2)3