QM/MM study of the absorption spectra of DsRed.M1 chromophores

We report geometries and vertical excitation energies for the red and green chromophores of the DsRed.M1 protein in the gas phase and in the solvated protein environment. Geometries are optimized using density functional theory (DFT, B3LYP functional) for the isolated chromophores and combined quantum mechanical/molecular mechanical (QM/MM) methods for the protein (B3LYP/MM). Vertical excitation energies are computed using DFT/MRCI, OM2/MRCI, and TDDFT as QM methods. In the case of the red chromophore, there is a general blue shift in the excitation energies when going from the isolated chromophore to the protein, which is caused both by structural changes and by electrostatic interactions with the environment. For the lowest ππ* transition, these two factors contribute to a similar extent to the overall DFT/MRCI shift of 0.4 eV. An enlargement of the QM region to include active‐site residues does not change the DFT/MRCI excitation energies much. The DFT/MRCI results are closest to experiment for both chromophores. OM2/MRCI and TDDFT overestimate the first vertical excitation energy by 0.3–0.5 and 0.2–0.4 eV, respectively, relative to the experimental or DFT/MRCI values. The experimental gap of 0.35 eV between the lowest ππ* excitation energies of the red (cis‐acylimine) and green (trans‐peptide) forms is well reproduced by DFT/MRCI and TDDFT (0.32 and 0.37 eV, respectively). A histogram spectrum for an equal mixture of the two forms, generated by OM2/MRCI calculations on 450 snapshots along molecular dynamics trajectories, matches the experimental spectrum quite well, with a gap of 0.23 eV and an overall blue shift of about 0.3 eV. DFT/MRCI appears as an attractive choice for calculating excitation energies in fluorescent proteins, without the shortcomings of TDDFT and computationally more affordable than CASSCF‐based approaches. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

DFT/TDDFT Studies on the Electronic Structures and Spectral Properties of Carbazole‐based Blue Light‐emitting Dendrimers

The purpose of this paper is to provide an in‐depth investigation of the electronic and optical properties of two series of carbazole‐based blue light‐emitting dendrimers, including 1 – 6 six oligomers. These materials show great potential for application in organic light‐emitting diodes as efficient blue‐light and red‐light emitting materials due to the tuning of the optical and electronic properties by the use of different electron donors (D) and electron acceptors (A). The geometric and electronic structures of these compounds in the ground state are calculated using density functional theory (DFT) and the ab initio HF, whereas the lowest singlet excited states were optimized by ab initio single excitation configuration interaction (CIS). All DFT calculations are performed using the B3LYP functional on 6‐31G* basis set. The outcomes show that the highest occupied molecular orbitals (HOMOs), lowest occupied molecular orbitals (LUMOs), energies gaps, ionization potentials, electron affinities and reorganization energies of each molecular are affected by different D and A moieties and different substitute positions.     

A semiempirical long-range corrected exchange correlation functional including a short-range Gaussian attenuation (LCgau-B97)

We applied an improved long‐range correction scheme including a short‐range Gaussian attenuation (LCgau) to the Becke97 (B97) exchange correlation functional. In the optimization of LCgau‐B97 functional, the linear parameters are determined by least squares fitting. Optimizing μ parameter (0.2) that controls long‐range portion of Hartree‐Fock (HF) exchange to excitation energies of large molecules (Chai and Head‐Gordon, J Chem Phys 2008, 128, 084106) and additional short‐range Gaussian parameters (a = 0.15 and k = 0.9) that controls HF exchange inclusion ranging from short‐range to mid‐range (0.5–3 Å) to ground state properties achieved high performances of LCgau‐B97 simultaneously on both ground state and excited state properties, which is better than other tested semiempirical density functional theory (DFT) functionals, such as ωB97, ωB97X, BMK, and M0x‐family. We also found that while a small μ value (～0.2) in LC‐DFT is appropriate to the local excitation and intramolecular charge‐transfer excitation energies, a larger μ value (0.42) is desirable in the Rydberg excitation‐energy calculations. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Ethylene coordination,insertion, and chain transfer at a cationic aluminum center: A comparative study with Ab Initio correlated level and density functional methods

The performance of correlated ab initio methods and DFT methods was compared for the propagation and chain transfer steps of ethylene polymerization by a model aluminum–amidinate system, [{HC(NH)₂}AlCH₂CH₃]⁺. All methods agree that the main chain transfer mechanism is β‐hydrogen transfer to the monomer (BHT), and that this is substantially easier than propagation; implications for the real Jordan system are discussed briefly. Counterpoise corrections are necessary to obtain reasonable olefin complexation energies. Activation energies are consistently lower at DFT (BP86, B3LYP) than at ab initio levels [MP2, MP3, MP4, CI, CCSD(T)]; the differences are particularly large (16 kcal/mol) for the BHT reaction. This is suggested to be related to the known problem of DFT in describing hydrogen bridged systems. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 398–410, 2000     

Excited state properties,fluorescence energies,and lifetimes of a poly(fluorene‐phenylene), based on TD‐DFT investigation

The structural and electronic properties of fluorene‐phenylene copolymer (FP)_n, n = 1–4 were studied by means of quantum chemical calculations based on density functional theory (DFT) and time dependent density functional theory (TD‐DFT) using B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest singlet excited state. It was found that (FP)_n is nonplanar in its ground state while the electronic excitations lead to planarity in its S₁ state. Absorption and fluorescence energies were calculated using TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods. Vertical excitation energies and fluorescence energies were obtained by extrapolating these values to infinite chain length, resulting in extrapolated values for vertical excitation energy of 2.89 and 2.87 eV, respectively. The S₁ ← S₀ electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is distinguishing in terms of oscillator strength. Fluorescence energies of (FP)_n calculated from TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods are 2.27 and 2.26 eV, respectively. Radiative lifetimes are predicted to be 0.55 and 0.51 ns for TD‐B3LYP/SVP and TD‐B3LYP/SVP+ calculations, respectively. These fundamental information are valuable data in designing and making of promising materials for LED materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

The sensing mechanism of fluoride anion for a novel molecule D2: Desilylation and intramolecular charge transfer

In this work, the sensing mechanism of a new fluoride chemosensor 12‐([tert‐butyldiphenylsilyl]oxy)‐8a,13a‐dihydro‐7H‐benzo[de]benzo[4,5]imidazo[2,1‐a]‐isoquinolin‐7‐one (abbreviated as D2) is investigated using density functional theory (DFT) and time‐dependent DFT (TDDFT) methods. The theoretical electronic spectra (vertical excitation energies and fluorescence peak) reproduced previous experimental results (D. Li et al., Spectrochim. Acta A Mol. Biomol. Spectrosc. 2017 , 185, 173), which confirms the rationality of the theoretical level used in this work. The constructed potential energy curve of the desilylation process suggests that the low barrier could be responsible for the rapid response to fluoride anions. Analyses of the binding energies show that only fluoride anion can be detected by D2 chemosensor in dimethylsulfoxide (DMSO). In view of the excitation process, the strong intramolecular charge transfer (ICT) process of the S₀ → S₁ transition explains the red shift of the absorption peak of the D2 sensor with the addition of fluoride anions. This work not only presents a straightforward sensing mechanism of sensing of the fluoride anion by the D2 chemosensor but should also play an important role in the synthesis and design of fluorescent sensors in future.     

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.     

Synthesis,Structures, and Properties of Highly Strained Cyclophenylene–Ethynylenes with Axial and Helical Chirality

Single and double cyclophenylene–ethynylenes (CPEs) with axial and helical chirality have been synthesized by the Sonogashira cross‐coupling of di‐ and tetraethynyl biphenyls with a U‐shaped prearomatic diiodoparaphenylene followed by reductive aromatization. X‐ray crystallographic analyses and DFT calculations revealed that the CPEs possess highly twisted bent structures. Bend angles on the edge of the paraphenylene units were close to the value of [5]cycloparaphenylene (CPP)—the smallest CPP to date. The double and single CPEs possessed stable chirality despite flexible biphenyl structures because of the high strain in the diethynyl–paraphenylene moiety. In both the single and double CPEs, orbital interactions along the biphenyl axis were observed by DFT calculations in LUMO and LUMO+2 of the single CPE and LUMO+1 of the double CPE, which likely cause lowering of these orbital energies. Concerning chiroptical properties: boosting of the g_abs value was observed in the biphenyl‐based double CPE, as well as the binaphthyl‐based single CPE, compared to the biphenyl‐based single CPE.     

A theoretical study of bond energies in model Si–H–Cl molecules using density functional approaches for representing Si surface chemistry

The reliability of density functional theory (DFT) methods for calculating Si(SINGLE BOND)2H, Si(SINGLE BOND)Cl, and Si(SINGLE BOND)Si bond energies is examined in reactions involving molecules and small clusters representing various surface sites appropriate for Si surface chemistry. Results are presented for systematic studies using a valence double-zeta polarization basis for both all-electron calculations and valence–electron calculations employing effective core potentials (ECPs). All-electron DFT results are comparable to much more demanding MP4, G2, and MC–SCF–CI calculations for computed bond energies. Whereas the use of ECPs introduces systematic energy differences of ca. 3–5 kcal/mol compared to AE results, depending on the type of bond involved, the use of ECPs for carrying out calculations on larger clusters is discussed where AE calculations become more computationally demanding. The convergence of Si bond energies as a function of replacing hydrogens with silyl groups is examined. In constructing models to describe etching processes involving Cl species on Si surfaces, the need for incorporating differences in thermochemistries for one-, two-, and three-coordinate Si surface sites is emphasized. Comparisons of semiempirical approaches for thermochemistries of Si-containing species find these methods somewhat less reliable for obtaining reliable bond energies compared to computationally more demanding DFT and ab initio correlated models. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 2075–2085, 1997     

Configuration interaction and density functional study of the influence of lithium cation complexation on vertical and adiabatic excitation energies of enones

The changes in the excited state energies of representative cyclic enones (cyclopentenone and cyclohexenone) induced by lithium ion coordination have been examined using ab initio and DFT methods. Quantitative estimates of the vertical triplet state energies were obtained using configuration interaction calculations at the CIS and CIS(D) levels with the 6‐31+G(d) basis. Inclusion of perturbative doubles corrections has a marked effect on the relative energies of the n–π* and π–π* triplet states. At both CI and CIS(D) levels, lithium complexation is predicted to raise the energy of the n–π* triplet state much more than the π–π* triplet. The trends obtained at the CIS(D) level are reproduced using B3LYP/6‐31+G(d) calculations. Adiabatic excitation energies were also computed by carrying out geometry optimization of the triplet states at the B3LYP level. While the separation between the geometry optimized n–π* and π–π* triplet states is very small for the parent enones, the π–π* triplet is clearly favored in the lithium complexes. These results suggest the possibility of reversing the reactive photoexcited state in enones through cation complexation. The conclusions provide a rationale for interesting variations in product distributions observed for enones in cation exchanged zeolites. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1598–1604, 2001     

Impact of replacement of the central benzene ring in anthracene by a heterocyclic ring on electronic excitations and reorganization energies in anthratetrathiophene molecules

The impact of changing the central benzene ring on the electronic excitations and reorganization energies (λ) of the anthratetrathiophene (ATT) molecules is studied by density functional theory (DFT) and time‐dependent DFT (TD‐DFT) quantum chemical calculations. The effect of changing the position of the sulfur atom at the periphery of anthracene on the optical and charge transfer properties is also studied. The calculated results suggest that the HOMO, LUMO, HOMO–LUMO energy gap, ionization potential (IP), electron affinity (EA), hole extraction potential (HEP), electron extraction potential (EEP), and reorganization energies (λ) are affected by replacing the central ring with different heterocyclic rings and the position of the sulfur atom. In addition, all molecules show good hole‐ and electron‐transport properties. This work may be helpful for future design and preparation of high‐performance charge‐transport materials.     

Design of low band gap polymers employing density functional theory—hybrid functionals ameliorate band gap problem

Band gaps in solids and excitation energies in finite systems are underestimated significantly if estimated from differences between eigenvalues obtained within the local spin density approximation (LSDA). In this article we present results on 20 small- and medium-sized π-systems which show that HOMO–LUMO energy differences obtained with the B3LYP, B3P86, and B3PW91 functionals are in good agreement with vertical excitation energies from UV-absorption spectra. The improvement is a result of the use of the exact Hartree–Fock exchange with hybrid methods. Negative HOMO energies and negative LUMO energies do not provide good estimates for IPs and EAs. In contrast to Hartree–Fock theory, where IPs are approximated well and EAs are given poorly, DFT hybrid methods underestimate IPs and EAs by about the same amount. LSDA yields reasonable EAs but poor IPs. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1943–1953, 1997     

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     

QM/MM calculation of solvent effects on absorption spectra of guanine

Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited‐state calculations were extracted from ground‐state molecular dynamics (MD) simulations using the self‐consistent‐charge density functional tight binding (SCC‐DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited‐state calculations used time‐dependent density functional theory (TDDFT) and the DFT‐based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto‐N7H and keto‐N9H guanine, with particular focus on solvent effects in the low‐energy spectrum of the keto‐N9H tautomer. When compared with the vertical excitation energies of gas‐phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC‐DFTB‐based rather than B3LYP‐based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas‐phase and solution spectra, typically ca. 0.1–0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent‐induced structural changes and from electrostatic solute–solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Quantum chemical comparison of vertical,adiabatic, and 0‐0 excitation energies: The PYP and GFP chromophores

A number of benchmark studies investigating the performance of quantum chemical methods for calculating vertical excitation energies are today available in the literature. However, less established is the variation between methods in their estimates of the differences between vertical, adiabatic, and 0‐0 excitation energies. To this end, such excitation energies are here calculated for the bright S₁ states of the anionic chromophores of the photoactive yellow protein (PYP) and the green fluorescent protein (GFP) in the gas phase using configuration interaction singles, complete active space self‐consistent field, coupled‐cluster singles and doubles, and time‐dependent density functional theory methods. Although the estimates of the excitation energies vary by more than 1 eV between the methods, the differences between the different types of excitation energies are found to be relatively method‐insensitive, varying by ～0.1 eV only for these particular chromophores. Specifically, the adiabatic energies are uniformly 0.10–0.17 (PYP) and 0.06–0.17 eV (GFP) lower than the vertical energies, and the 0‐0 energies are similarly 0.09–0.14 (PYP) and 0.07–0.17 eV (GFP) lower than the adiabatic energies. © 2012 Wiley Periodicals, Inc.     

Combined experimental and computational modeling studies on 3,5‐dimethyl‐pyrazole‐1‐carbodithioic acid benzyl ester

The title compound, 3,5‐Dimethyl‐pyrazole‐1‐carbodithioic acid benzyl ester, has been synthesized and structurally characterized by X‐ray single crystal diffraction, elemental analysis, IR spectra, and UV‐Vis spectrum. The crystal belongs to orthorhombic, space group P212121, with a = 5.3829(15), b = 11.193(3), c = 21.824(6) Å, V = 1315.0(6) ų, and Z = 4. The molecules are connected via intermolecular C–H···N hydrogen bonds into 1D infinite chains. The crystal structure is consolidated by the intramolecular C–H···S hydrogen bonds. Furthermore, Density functional theory (DFT) calculations of the structure, stabilities, orbital energies, composition characteristics of some frontier molecular orbitals and Mulliken charge distributions of the title compound were performed by means of Gaussian 03W package and taking B3LYP/6‐31G(d) basis set. The time‐dependent DFT (TD‐DFT) calculations have been employed to calculate the electronic spectrum of the title compound, and the UV‐Vis spectra has been discussed on this basis. The results show that DFT method at B3LYP/6‐31G(d) level can well reproduce the structure of the title compound. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Intramolecular electronic communication in a dimethylaminoazobenzene–fullerene C60 dyad: An experimental and TD‐DFT study

An electronically push–pull type dimethylaminoazobenzene–fullerene C₆₀ hybrid was designed and synthesized by tailoring N,N‐dimethylaniline as an electron donating auxochrome that intensified charge density on the β‐azonitrogen, and on N‐methylfulleropyrrolidine (NMFP) as an electron acceptor at the 4 and 4′ positions of the azobenzene moiety, respectively. The absorption and charge transfer behavior of the hybrid donor‐bridge‐acceptor dyad were studied experimentally and by performing TD‐DFT calculations. The TD‐DFT predicted charge transfer interactions of the dyad ranging from 747 to 601 nm were experimentally observed in the UV‐vis spectra at 721 nm in toluene and dichloromethane. A 149 mV anodic shift in the first reduction potential of the N?N group of the dyad in comparison with the model aminoazobenzene derivative further supported the phenomenon. Analysis of the charge transfer band through the orbital picture revealed charge displacement from the n_(N?N) (nonbonding) and π _(N?N) type orbitals centered on the donor part to the purely fullerene centered LUMOs and LUMO+n orbitals, delocalized over the entire molecule. The imposed electronic perturbations on the aminoazobenzene moiety upon coupling it with C₆₀ were analyzed by comparing the TD‐DFT predicted and experimentally observed electronic transition energies of the dyad with the model compounds, NMFP and (E)‐N,N‐dimethyl‐4‐(p‐tolyldiazenyl)aniline (AZNME). The n_(N?N) → π*_(N?N) and π_(N?N) → π*_(N?N) transitions of the dyad were bathochromically shifted with a significant charge transfer character. The shifting of π_(N?N) → π*_(N?N) excitation energy closer to the n → π*_(N?N) in comparison with the model aminoazobenzene emphasized the predominant existence of charge separated quinonoid‐like ground state electronic structure. Increasing solvent polarity introduced hyperchromic effect in the π_(N?N) → π*_(N?N) electronic transition at the expense of transitions involved with benzenic states, and the extent of intensity borrowing was quantified adopting the Gaussian deconvolution method. On a comparative scale, the predicted excitation energies were in reasonable agreement with the observed values, demonstrating the efficiency of TD‐DFT in predicting the localized and the charge transfer nature of transitions involved with large electronically asymmetric molecules with HOMO and LUMO centered on different parts of the molecular framework. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A detailed theoretical study on the excited‐state hydrogen‐bonding dynamics and the proton transfer mechanism for a novel white‐light fluorophore

In this work, density functional theory (DFT) and time‐dependent DFT (TDDFT) methods were used to investigate the excited‐state dynamics of the excited‐state hydrogen‐bonding variations and proton transfer mechanism for a novel white‐light fluorophore 2‐(4‐[dimethylamino]phenyl)‐7‐hyroxy‐6‐(3‐phenylpropanoyl)‐4H‐chromen‐4‐one ( 1 ). The methods we adopted could successfully reproduce the experimental electronic spectra, which shows the appropriateness of the theoretical level in this work. Using molecular electrostatic potential (MEP) as well as the reduced density gradient (RDG) versus the product of the sign of the second largest eigenvalue of the electron density Hessian matrix and electron density (sign[λ₂]ρ), we demonstrate that an intramolecular hydrogen bond O1–H2···O3 should be formed spontaneously in the S₀ state. By analyzing the chemical structures, infrared vibrational spectra, and hydrogen‐bonding energies, we confirm that O1–H2·O3 should be strengthened in the S₁ state, which reveals the possibility of an excited‐state intramolecular proton transfer (ESIPT) process. On investigating the excitation process, we find the S₀ → S₁ transition corresponding to the charge transfer, which provides the driving force for ESIPT. By constructing the potential energy curves, we show that the ESIPT reaction results in a dynamic equilibrium in the S₁ state between the forward and backward processes, which facilitates the emission of white light.     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

How dependent are molecular and atomic properties on the electronic structure method? Comparison of Hartree‐Fock,DFT, and MP2 on a biologically relevant set of molecules

This article compares molecular properties and atomic properties defined by the quantum theory of atoms in molecules (QTAIM) obtained from three underlying levels of theory: MP2(full), density functional theory (DFT) (B3LYP), and Hartree‐Fock (H‐F). The same basis set (6‐311++G(d,p)) has been used throughout the study. The calculations and comparisons were applied to a set of 30 small molecules representing common fragments of biological molecules. The molecular properties investigated are the energies and the electrostatic moments (up to and including the quadrupoles), and the atomic properties include electron populations (and atomic charge), atomic dipolar and quadrupolar polarizations, atomic volumes, and corrected and raw atomic energies. The Cartesian distance between dipole vectors and the Frobenius distance between the quadrupole tensors calculated at the three levels of theory provide a measure of their correlation (or lack thereof). With the exception of energies (atomic and molecular), it is found that both DFT and H‐F are in excellent agreement with MP2, especially with regards to the electrostatic mutipoles up to the quadrupoles, but DFT and MP2 agree better in almost all studied properties (with the exception of molecular geometries). QTAIM properties whether obtained from H‐F, DFT(B3LYP), or MP2 calculations when used in the construction of empirical correlations with experiment such as quantitative structure‐activity‐(or property)‐relationships (QSAR/QSPR) are equivalent (because the properties calculated at the three levels are very highly correlated among themselves with r² typically >0.95, and therefore preserving trends). These results suggest that the massive volume of results that were published in the older literature at the H‐F level is valid especially when used to study trends or in QSAR or QSPR studies, and, as long as our test set of molecules is representative, there is no pressing need to re‐evaluate them at other levels of theory except when inadequate basis sets were used by today's standards. Extensive tabulation of molecular and atomic properties at the three theoretical levels is available in the Supporting Information, including optimized geometries, molecular energies, virial ratios, molecular electrostatic moments up to and including hexadecapoles, atomic populations, atomic volumes, atomic electrostatic moments up to and including the quadrupoles, and atomic energies. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010