The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures.

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches.

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0-0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang-Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang-Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0-0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

Structure Prediction of Self-Assembled Dye Aggregates from Cryogenic Transmission Electron Microscopy, Molecular Mechanics, and Theory of Optical Spectra.

Cryogenic transmission electron microscopy (cryo-TEM) studies suggest that TTBC molecules self-assemble in aqueous solution to form single-walled tubes with a diameter of about 35 Å. In order to reveal the arrangement and mutual orientations of the individual molecules in the tube, we combine information from crystal structure data of this dye with a calculation of linear absorbance and linear dichroism spectra and molecular dynamics simulations. We start with wrapping crystal planes in different directions to obtain tubes of suitable diameter. This set of tube models is evaluated by comparing the resulting optical spectra with experimental data. The tubes that can explain the spectra are investigated further by molecular dynamics simulations, including explicit solvent molecules. From the trajectories of the most stable tube models, the short-range ordering of the dye molecules is extracted and the optimization of the structure is iteratively completed. The final structural model is a tube of rings with 6-fold rotational symmetry, where neighboring rings are rotated by 30° and the transition dipole moments of the chromophores form an angle of 74° with respect to the symmetry axis of the tube. This model is in agreement with cryo-TEM images and can explain the optical spectra, consisting of a sharp red-shifted J-band that is polarized parallel to to the symmetry axis of the tube and a broad blue-shifted H-band polarized perpendicular to this axis. The general structure of the homogeneous spectrum of this hybrid HJ-aggregate is described by an analytical model that explains the difference in redistribution of oscillator strength inside the vibrational manifolds of the J- and H-bands and the relative intensities and excitation energies of those bands. In addition to the particular system investigated here, the present methodology can be expected to aid the structure prediction for a wide range of self-assembled dye aggregates.