Modeling Optical Spectra of Large Organic Systems Using Real-Time Propagation of Semiempirical Effective Hamiltonians

Super-resolution techniques to simulate electronic spectra of large molecular systems

An accurate treatment of electronic spectra in large systems with a technique such as time-dependent density functional theory is computationally challenging. Due to the Nyquist sampling theorem, direct real-time simulations must be prohibitively long to achieve suitably sharp resolution in frequency space. Super-resolution techniques such as compressed sensing and MUSIC assume only a small number of excitations contribute to the spectrum, which fails in large molecular systems where the number of excitations is typically very large. We present an approach that combines exact short-time dynamics with approximate frequency space methods to capture large narrow features embedded in a dense manifold of smaller nearby peaks. We show that our approach can accurately capture narrow features and a broad quasi-continuum of states simultaneously, even when the features overlap in frequency. Our approach is able to reduce the required simulation time to achieve reasonable accuracy by a factor of 20-40 with respect to standard Fourier analysis and shows promise for accurately predicting the whole spectrum of large molecules and materials.

Modeling of High-Harmonic Generation in the C60 Fullerene Using Ab Initio, DFT-Based, and Semiempirical Methods

We report calculations of the high-harmonic generation spectra of the C60 fullerene molecule carried out by employing a diverse set of real-time time-dependent quantum chemical methods. All methodologies involve expanding the propagated electronic wave function in bases consisting of the ground and singly excited time-independent eigenstates obtained through the solution of the corresponding linear-response equations. We identify the correlation and exchange effect in the spectra by comparing the results from methods relying on the Hartree-Fock reference determinant with those obtained using approaches based on the density functional theory with different exchange-correlation functionals. The effect of the full random-phase approximation treatment of the excited electronic states is also analyzed and compared with the configuration interaction singles and the Tamm-Dancoff approximation. We also showcase the fact that the real-time extension of the semiempirical method INDO/S can be effectively applied for an approximate description of laser-driven dynamics in large systems.

ATR-far-ultraviolet spectroscopy: a challenge to new σ chemistry

This review reports the recent progress on ATR-far ultraviolet (FUV) spectroscopy in the condensed phase. ATR-FUV spectroscopy for liquids and solids enables one to explore various topics in physical chemistry, analytical chemistry, nanoscience and technology, materials science, electrochemistry, and organic chemistry. In this review, we put particular emphasis on the three major topics: (1) studies on electronic transitions and structures of various molecules, which one cannot investigate via ordinary UV spectroscopy. The combined use of ATR-FUV spectroscopy and quantum chemical calculations allows for the investigation of various electronic transitions, including σ, n-Rydberg transitions. ATR-FUV spectroscopy may open a new avenue for σ-chemistry. (2) ATR-FUV spectroscopy enables one to measure the first electronic transition of water at approximately 160 nm without peak saturation. Using this band, one can study the electronic structure of water, aqueous solutions, and adsorbed water. (3) ATR-FUV spectroscopy has its own advantages of the ATR method as a surface analysis method. ATR-FUV spectroscopy is a powerful technique for exploring a variety of top surface phenomena (∼50 nm) in adsorbed water, polymers, graphene, organic materials, ionic liquids, and so on. This review briefly describes the principles, characteristics, and instrumentation of ATR-FUV spectroscopy. Next, a detailed description about quantum chemical calculation methods for FUV and UV regions is given. The recent application of ATR-FUV-UV spectroscopy studies on electronic transitions from σ orbitals in various saturated molecules is introduced first, followed by a discussion on the applications of ATR-FUV spectroscopy to studies on water, aqueous solutions, and adsorbed water. Applications of ATR-FUV spectroscopy in the analysis of other materials such as polymers, ionic liquids, inorganic semiconductors, graphene, and carbon nanocomposites are elucidated. In addition, ATR-FUV-UV-vis spectroscopy focusing on electrochemical interfaces is outlined. Finally, FUV-UV-surface plasmon resonance studies are discussed.

Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group

In this work, we simulate the electron dynamics in molecular systems with the time-dependent density matrix renormalization group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact nonrelativistic electronic Hamiltonian. We show that by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.

A new release of MOPAC incorporating the INDO/S semiempirical model with CI excited states

The semiempirical INDO/S Hamiltonian is incorporated into a new release of MOPAC2016. The MOPAC2016 software package has long been at the forefront of semiempirical quantum chemical methods (SEQMs) for small molecules, proteins, and solids and until this release has included only NDDO-type SEQMs. The new code enables the calculation of excited states using the INDO/S Hamiltonian combined with a configuration interaction (CI) approach using single excitations (CIS), single and double excitations (CISD), or multiple reference determinants (MRCI) where reference determinants are generated using a complete active space (CAS) approach. The capacity to perform excited-state calculations beyond the CIS level makes INDO/CI one of the few low-cost computational methods capable of accurately modeling states with substantial double-excitation character. Solvent corrections to the ground-state and excited-state energies can be computed using the COSMO implicit solvent model, incorporating state-specific corrections to the excited states based on the solvent refractive index. This code produces physically reasonable electronic structures, absorption spectra, and solvatochromic shifts at low computational costs for systems up to hundreds of atoms, and for both organic molecules and metal clusters.

NWChem: Past, present, and future

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.

Optical excitations of chlorophyll a and b monomers and dimers

A necessary first step in the development of technologies such as artificial photosynthesis is understanding the photoexcitation process within the basic building blocks of naturally occurring light harvesting complexes (LHCs). The most important of these building blocks in biological LHCs such as LHC II from green plants are the chlorophyll a (Chl a) and chlorophyll b (Chl b) chromophores dispersed throughout the protein matrix. However, efforts to describe such systems are still hampered by the lack of computationally efficient and accurate methods that are able to describe optical absorption in large biomolecules. In this work, we employ a highly efficient linear combination of atomic orbitals (LCAOs) to represent the Kohn-Sham (KS) wave functions at the density functional theory (DFT) level and perform time-dependent density functional theory (TDDFT) calculations in either the reciprocal space and frequency domain (LCAO-TDDFT-k-ω) or real space and time domain (LCAO-TDDFT-r-t) of the optical absorption spectra of Chl a and b monomers and dimers. We find that our LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations reproduce results obtained with a plane-wave (PW) representation of the KS wave functions (PW-TDDFT-k-ω) but with a significant reduction in computational effort. Moreover, by applying the Gritsenko, van Leeuwen, van Lenthe, and Baerends solid and correlation derivative discontinuity correction Δ_x to the KS eigenenergies, with both LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t methods, we are able to semiquantitatively reproduce the experimentally measured photoinduced dissociation results. This work opens the path to first principles calculations of optical excitations in macromolecular systems.