An open-source program named VeloxChem has been developed for the calculation of electronic real and complex linear response functions at the levels of Hartree–Fock and Kohn–Sham density functional theories. With an object-oriented program structure written in a Python/C++ layered fashion, VeloxChem enables time-efficient prototyping of novel scientific approaches without sacrificing computational efficiency, so that molecular systems involving up to and beyond 500 second-row atoms (or some 10,000 contracted and in part diffuse Gaussian basis functions) can be routinely addressed. In addition, VeloxChem is equipped with a polarizable embedding scheme for the treatment of the classical electrostatic interactions with an environment that in turn is modeled by atomic site charges and polarizabilities. The underlying hybrid message passing interface (MPI)/open multiprocessing (OpenMP) parallelization scheme makes VeloxChem suitable for execution in high-performance computing cluster environments, showing even slightly beyond linear scaling for the Fock matrix construction with use of up to 16,384 central processing unit (CPU) cores. An efficient—with respect to convergence rate and overall computational cost—multifrequency/gradient complex linear response equation solver enables calculations not only of conventional spectra, such as visible/ultraviolet/X-ray electronic absorption and circular dichroism spectra, but also time-resolved linear response signals as due to ultra-short weak laser pulses. VeloxChem distributed under the GNU Lesser General Public License version 2.1 (LGPLv2.1) license and made available for download from the homepage https://veloxchem.org. This article is categorized under: Software > Quantum Chemistry Electronic Structure Theory > Density Functional Theory Theoretical and Physical Chemistry > Spectroscopy.

The Gator program has been developed for computational spectroscopy and calculations of molecular properties using real and complex propagators at the correlated level of wave function theory. Currently, the focus lies on methods based on the algebraic diagrammatic construction (ADC) scheme up to the third order of perturbation theory. An auxiliary Fock matrix-driven implementation of the second-order ADC method for excitation energies has been realized with an underlying hybrid MPI/OpenMP parallelization scheme suitable for execution in high-performance computing cluster environments. With a modular and object-oriented program structure written in a Python/C++ layered fashion, Gator additionally enables time-efficient prototyping of novel scientific approaches, as well as interactive notebook-driven training of students in quantum chemistry. This article is categorized under: Computer and Information Science > Computer Algorithms and Programming Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry.

We present the derivation and implementation of an efficient matrix-folded formalism for the evaluation of complex response functions at the level of the second-order algebraic-diagrammatic construction ADC(2) scheme. The adopted hybrid MPI/OpenMP parallelization strategy enables the utilization of vast computational resources without significant performance losses as demonstrated by consistent size-scalings for the computationally most demanding steps obtained from calculations of the UV/vis spectra of a guanine oligomer series ranging up to 1,032 contracted basis functions. Further analysis of the convergence behavior of the involved iterative subspace algorithm revealed the superiority of a frequency-separated treatment of response equations even for a large spectral window including 101 frequencies. We demonstrate the applicability to general quantum mechanical operators by the first reported electric circular dichroism spectrum calculated with the complex polarization propagator approach at the ADC(2) level of theory.

While the anomalous non-additive size-dependencies of static dipole polarizabilities and van der Waals C-6 dispersion coefficients of carbon fullerenes are well established, the widespread reported scalings for the latter (ranging from N-2.2 to N-2.8) call for a comprehensive first-principles investigation. With a highly efficient implementation of the linear complex polarization propagator, we have performed Hartree-Fock and Kohn-Sham density functional theory calculations of the frequency-dependent polarizabilities for fullerenes consisting of up to 540 carbon atoms. Our results for the static polarizabilities and C-6 coefficients show scalings of N-1.2 and N-2.2, respectively, thereby deviating significantly from the previously reported values obtained with the use of semi-classical/empirical methods. Arguably, our reported values are the most accurate to date as they represent the first ab initio or first-principles treatment of fullerenes up to a convincing system size.

We demonstrate that contemporary computational resources allow for accurate theoretical studies of systems matching recent advances in experimental helicene chemistry. Concerned with first-principles calculations of carbohelicenes, our work surpasses CH[12] as the largest system investigated to date and unravels trends in the electronic structure of the low-lying states of the homologous series. Utilizing a highly efficient implementation of linear response algorithms, we present electronic circular dichroism (CD) spectra of carbohelicenes ranging from CH[5] to CH[30] at the level of Kohn-Sham density-functional theory. Our results for a systematic increase in system size show the emergence of new CD bands that subsequently rise to intensities dominating the spectrum. The spectral band progressions exhibit a periodicity directly linked to the number of overlapping layers of conjugation. While our findings rectify the current understanding of the electronic structure of carbohelicenes, they also serve as a general call for caution regarding the extrapolation of trends from small system ranges.