Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree–Fock, Kohn–Sham, multiconfigurational self-consistent-field, Møller–Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

We provide a new and definitive spectral assignment for the absorption, emission, high-resolution fluorescence excitation, linear dichroism, and/or magnetic circular dichroism spectra of 32 chlorophyllides in various environments. This encompases all data used to justify previous assignments and provides a simple interpretation of unexplained complex decoherence phenomena associated with Qx → Qy relaxation. Whilst most chlorophylls conform to the Gouterman model and display two independent transitions Qx (S2) and Qy (S1), strong vibronic coupling inseparably mixes these states in chlorophyll-a. This spreads x-polarized absorption intensity over the entire Q-band system to influence all exciton-transport, relaxation and coherence properties of chlorophyll-based photosystems. The fraction of the total absorption intensity attributed to Qx ranges between 7% and 33%, depending on chlorophyllide and coordination, and is between 10% and 25% for chlorophyll-a. CAM-B3LYP density-functional-theory calculations of the band origins, relative intensities, vibrational Huang-Rhys factors, and vibronic coupling strengths fully support this new assignment.

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Using the long-range correction (LC) density functional theory (DFT) scheme introduced by Iikura et al. [J. Chem. Phys. 115, 3540 (2001)] and the Coulomb-attenuating model (CAM-B3LYP) of Yanai et al. [Chem. Phys. Lett. 393, 51 (2004)], we have calculated a series of properties that are known to be poorly reproduced by standard functionals: Bond length alternation of pi-conjugated polymers, polarizabilities of delocalized chains, and electronic spectra of extended dyes. For each of these properties, we present cases in which traditional hybrid functionals do provide accurate results and cases in which they fail to reproduce the correct trends. The quality of the results is assessed with regard to experimental values and/or data arising from electron-correlated wave function approaches. It turns out that (i) both LC-DFT and CAM-B3LYP provide an accurate bond length alternation for polyacetylene and polymethineimine, although for the latter they decrease slightly too rapidly with chain length. (ii) The LC generalized gradient approximation and MP2 polarizabilities of long polyphosphazene and polymethineimine oligomers agree almost perfectly. In the same way, CAM-B3LYP corrects the major part of the B3LYP faults. (iii) LC and CAM techniques do not help in correcting the nonrealistic evolution with chain length of the absorption wavelengths of cyanine derivatives. In addition, though both schemes significantly overestimate the ground to excited state transition energy of substituted anthraquinone dyes, they provide a more consistent picture once a statistical treatment is performed than do traditional hybrid functionals.

While density functional theory (DFT) has been proven to be extremely useful for the prediction of thermodynamic and spectroscopic properties of molecules, to date most functionals used in common implementations of DFT display a systematic failure to predict the properties of charge-transfer processes. While this is explicitly manifest in Rydberg transitions of atoms and molecules and in molecular charge-transfer spectroscopy, it also becomes critical for systems containing extended conjugation such as polyenes and other conducting polymers, porphyrins, chlorophylls, etc. A new density functional, a Coulomb-attenuated hybrid exchange-correlation functional (CAM-B3LYP), has recently been developed specifically to overcome these limitations, and it has been shown to properly predict molecular charge-transfer spectra. Here, we demonstrate that it predicts qualitatively reasonable spectra for porphyrin, some oligoporphyrins, and chlorophyll. However, alternate density functionals developed to overcome the same limitations such as current-density functional theory are shown, in their present implementation, to remain inadequate. The CAM-B3LYP results are shown to be in excellent agreement with complete-active-space plus second-order Møller-Plesset perturbation theory and symmetry-adapted cluster configuration interaction calculations: These depict the N and higher bands of porphyrins and chlorophylls as being charge-transfer bands associated with localization of molecular orbitals on individual pyrrole rings. The validity of the basic Gouterman model for the spectra of porphyrins and chlorophylls is confirmed, rejecting modern suggestions that non-Gouterman transitions lie close in energy to the Q-bands of chlorophylls. As porphyrins and chlorophylls provide useful paradigms for problems involving extended conjugation, the results obtained suggest that many significant areas of nanotechnology and biotechnology may now be realistically treated by cost-effective density-functional-based computational methods.

Coupled cluster singles and doubles linear response (CCLR) calculations have been carried out for excitation energies and dipole transition strengths for the lowest excitations in LiH, CH+, and C4 and the results compared with the results from a CI-like approach to equation of motion coupled cluster (EOMCC) . The transition strengths are similar in the two approaches for single molecule calculations on small systems. However, the CCLR approach gives size-intensive dipole transition strengths, while the EOMCC formalism does not. Thus, EOMCC calculations can give unphysically dipole transition strengths, e.g., in EOMCC calculations on a sequence of noninteracting LiH systems we obtained a negative dipole strength for the lowest totally symmetric dipole allowed transition for 19 or more noninteracting LiH systems. The CCLR approach is shown to be a very attractive "black box" approach for the calculation of transition moments.

Articles you may be interested inPolarized atomic orbitals for self-consistent field electronic structure calculations Nuclear shielding calculations are presented for multiconfigurational self-consistent field wave functions using London atomic orbitals (gauge invariant atomic orbitals). Calculations of nuclear shieldings for eight molecules (H 2 0, H 2 S, CH 4 , N 2 , CO, HF, F 2 , and S02) are presented and compared to corresponding individual gauges for localIzed orbitals (IGLO) results. The London results show better basis set convergence than IGLO, especially for heavier atoms. It is shown that the choice of active space is crucial for determination of accurate nuclear shielding constants.8178

Defect states in 2D materials present many possible uses but both experimental and computational characterization of their spectroscopic properties is difficult. We provide and compare results from 13 DFT and ab initio computational methods for up to 25 excited states of a paradigm system, the V N C B defect in hexagonal boron nitride (h-BN). Studied include include: (i) potentially catastrophic effects for computational methods arising from the multi-reference nature of the closed-shell and open-shell states of the defect, which intrinsically involves broken chemical bonds, (ii) differing results from DFT and time-dependent DFT (TDDFT) calculations, (iii) comparison of cluster models to periodic-slab models of the defect, (iv) the starkly differing effects of nuclear relaxation on the various electronic states as broken bonds try to heal that control the widths of photoabsorption and photoemission spectra, (v) the effect of zero-point energy and entropy on free-energy differences, (vi) defect-localized and conduction/valence band transition natures, and (vii) strategies needed to ensure that the lowest-energy state of a defect can be computationally identified. Averaged state-energy differences of 0.3 eV are found between CCSD(T) and MRCI energies, with thermal effects on free energies sometimes also being of this order. However, DFT-based methods can perform very poorly. Simple generalized-gradient functionals like PBE fail at the most basic level and should never be applied to defect states. Hybrid functionals like HSE06 work very well for excitations within the triplet manifold of the defect, with an accuracy equivalent to or perhaps exceeding the accuracy of the ab initio methods used. However, HSE06 underestimates triplet state energies by on average 0.7 eV compared to closed-shell singlet states, while open-shell singlet states are predicted to be too low in energy by 1.0 eV. This leads to miss-assignment of the ground state of the V N C B defect. Long-range corrected functionals like CAM-B3LYP are shown to work much better and to represent the current entry level for DFT calculations on defects. As significant differences between cluster and periodic-slab models are also found, the widespread implementation of such functionals in periodic codes is in urgent need.