The eChem project has been previously presented as an interactive platform for quantum and computational chemistry education [J. Chem. Educ. 100, 1664-1671]. However, education is only one side of the eChem project. Another aspect is that it highly accelerates method development by means of code prototyping in notebooks. Complex equations can be understood, and algorithms are examined before the actual software programming step is carried out. Here, the benefits of notebooks for code prototyping are illustrated using the example of vibrational spectroscopy-a type of spectroscopy which involves complex equations with a large number of terms.

The eChem project features an e-book published as a web page (10.30746/978-91-988114-0-7), collecting a repository of Jupyter notebooks developed for the dual purpose of explaining and exploring the theory underlying computational chemistry in a highly interactive manner as well as providing a tutorial-based presentation of the complex workflows needed to simulate embedded molecular systems of real biochemical and/or technical interest. For students ranging from beginners to advanced users, the eChem book is well suited for self-directed learning, but workshops led by experienced instructors and targeting student bodies with specific needs and interests can readily be formed from its components. This has been done by using eChem as the base for a workshop directed toward graduate students learning the theory and practices of quantum chemistry, resulting in very positive assessment of the interactive nature of this framework. The members of the eChem team are engaged in both education and research, and as a mirroring activity, we develop the open-source software upon which this e-book is predominantly based. The overarching vision and goal of our work is to provide a science- and education-enabling software platform for quantum molecular modeling on contemporary and future high-performance computing systems, and to document the resulting development and workflows in the eChem book.

In this article, we analyze the electronic structure modifications of triphenylamine (TPA), a well-known electron donor molecule widely used in photovoltaics and optoelectronics, upon deposition on Au(111) at a monolayer coverage. A detailed study was carried out by synchrotron radiation-based photoelectron spectroscopy, near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, scanning tunneling microscopy (STM), and ab initio calculations. We detect a new feature in the pre-edge energy region of the N K-edge NEXAFS spectrum that extends over 3 eV, which we assign to transitions involving new electronic states. According to our calculations, upon adsorption, a number of new unoccupied electronic states fill the energy region between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of the free TPA molecule and give rise to the new feature the pre-edge region of the NEXAFS spectrum. This finding highlights the occurrence of a considerable modification of the electronic structure of TPA. The appearance of new states in the HOMO-LUMO gap of TPA when adsorbed on Au(111) has crucial implications for the design of molecular nanoelectronic devices based on similar donor systems.

Quantum chemistry is a powerful tool. It is now possible to model complex chemical processes even on a laptop getting insights into matter at its fundamental scale.

But quantum chemistry is also very complex. Answering a chemical question requires selecting parameters among a wide variety of choices. Choosing a model system, an electronic structure method, a basis set, a set of properties, and a wide array of parameters which can affect the results in small but sometimes meaningful way… It can be a very daunting task, even for veterans of the field.

Similarly, for those who wish to get a deeper understanding of a method, going through the pages of equation often riddled with inconsistent notations and formulations is very challenging. And at the end, the link between the equation and the computer implementation found in existing softwares can be vague at best.

We believe that a core issue is that humans are not good at learning in abstract terms. We can get very far with a lecture or a textbook, but we will never build as much intuition about how a clock work as by simply breaking one apart and rebuilding it from scratch.

This is exactly the aim of this page, allowing a hands-on approach to computational chemistry. Together we will dismantle the black box that a computational chemistry code often seems to be, go through all the cogs and gears, and build back together some of the main computational methods of modern computational chemistry. We will do this by presenting the underlying equations, all expressed with consistent notations, as well as by suggesting a simple python implementation, to really display in action how the theory is implemented into a practical tool. Additionally, we will put these methods in context by showing how they can be used to address concrete chemical questions, discussing the strengths and weaknesses of each method and how to best use them to solve practical problems.

The influence of core-hole delocalization for x-ray photoelectron, x-ray absorption, and x-ray emission spectrum calculations is investigated in detail using approaches including response theory, transition-potential methods, and ground state schemes. The question of a localized/delocalized vacancy is relevant for systems with symmetrically equivalent atoms, as well as near-degeneracies that can distribute the core orbitals over several atoms. We show that the issues relating to core-hole delocalization are present for calculations considering explicit core-hole states, e.g., when using a core-excited or core-ionized reference state or for fractional occupation numbers. As electron correlation eventually alleviates the issues, but even when using coupled-cluster single-double and perturbative triple, there is a notable discrepancy between core-ionization energies obtained with localized and delocalized core-holes (0.5 eV for the carbon K-edge). Within density functional theory, the discrepancy correlates with the exchange interaction involving the core orbitals of the same spin symmetry as the delocalized core-hole. The use of a localized core-hole allows for a reasonably good inclusion of relaxation at a lower level of theory, whereas the proper symmetry solution involving a delocalized core-hole requires higher levels of theory to account for the correlation effects involved in orbital relaxation. For linear response methods, we further show that if x-ray absorption spectra are modeled by considering symmetry-unique sets of atoms, care has to be taken such that there are no delocalizations of the core orbitals, which would otherwise introduce shifts in absolute energies and relative features.

This joint experimental-theoretical spectroscopy study of the fullerene derivative PC60BM ([6,6]-phenyl-C60-butyric acid methyl ester) aims to improve the understanding of the effect of photooxidation on its electronic structure. We have studied spin-coated thin films of PC60BM by X-ray Photoelectron Spectroscopy (XPS), Near-edge X-ray Absorption Fine Structure (NEXAFS) spectroscopy, and Fourier Transform Infrared Spectroscopy (FTIR), before and after intentional exposure to simulated sunlight in air for different lengths of time. The π* resonance in the C1s NEXAFS spectrum was found to be a very sensitive probe for the early changes to the fullerene cage, while FTIR spectra, in combination with O1s NEXAFS spectra, enabled the identification of the oxidation products. The changes observed in the spectra obtained by these complementary methods were compared with the corresponding Density Functional Theory (DFT) calculated single-molecule spectra of a large set of in silico generated oxidation products of PC60BM where oxygen atoms were attached to the C60 cage. This comparison confirms that photooxidation of PC60BM disrupts the conjugation of the fullerene cage by a transition from sp2 to sp3-hybridized carbon and causes the formation of several oxidation products, earlier proposed for C60. The agreement between experimental and calculated IR spectra suggests moreover the presence of dicarbonyl and anhydride structures on the fullerene cage, in combination with cage opening at the adsorption site. By including PC60BM with physisorbed O2 molecules on the cage in our theoretical description in order to model oxygen diffused through the film, the experimental O1s XPS and O1s NEXAFS spectra could be reproduced.

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s -> pi* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.

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.

A theoretical and experimental study of the gas phase and liquid acetic acid based on resonant inelastic x-ray scattering (RIXS) spectroscopy is presented. We combine and compare different levels of theory for an isolated molecule for a comprehensive analysis, including electronic and vibrational degrees of freedom. The excitation energy scan over the oxygen K-edge absorption reveals nuclear dynamic effects in the core-excited and final electronic states. The theoretical simulations for the monomer and two different forms of the dimer are compared against high-resolution experimental data for pure liquid acetic acid. We show that the theoretical model based on a dimer describes the hydrogen bond formation in the liquid phase well and that this bond formation sufficiently alters the RIXS spectra, allowing us to trace these effects directly from the experiment. Multimode vibrational dynamics is accounted for in our simulations by using a hybrid time-dependent stationary approach for the quantum nuclear wave packet simulations, showing the important role it plays in RIXS.

Quenching of vibrational excitations in resonant inelastic X-ray scattering (RIXS) spectra of liquid acetic acid is observed. At the oxygen core resonance associated with localized excitations at the O-H bond, the spectra lack the typical progression of vibrational excitations observed in RIXS spectra of comparable systems. We interpret this phenomenon as due to strong rehybridization of the unoccupied molecular orbitals as a result of hydrogen bonding, which however cannot be observed in x-ray absorption but only by means of RIXS. This allows us to address the molecular structure of the liquid, and to determine a lower limit for the average molecular chain length.