Ultrafast structural dynamics with different spatial and temporal scales were investigated during photodissociation of carbon monoxide (CO) from iron(II)-heme in bovine myoglobin during the first 3 ps following laser excitation. We used simultaneous X-ray transient absorption (XTA) spectroscopy and X-ray transient solution scattering (XSS) at an X-ray free electron laser source with a time resolution of 80 fs. Kinetic traces at different characteristic X-ray energies were collected to give a global picture of the multistep pathway in the photodissociation of CO from heme. In order to extract the reaction coordinates along different directions of the CO departure, XTA data were collected with parallel and perpendicular relative polarizations of the laser pump and X-ray probe pulse to isolate the contributions of electronic spin state transition, bond breaking, and heme macrocycle nuclear relaxation. The time evolution of the iron K-edge X-ray absorption near edge structure (XANES) features along the two major photochemical reaction co-ordinates, i.e., the iron(II)-CO bond elongation and the heme macrocycle doming relaxation were modeled by time-dependent density functional theory calculations. Combined results from the experiments and computations reveal insight into interplays between the nuclear and electronic structural dynamics along the CO photodissociation trajectory. Time-resolved small-angle X-ray scattering data during the same process are also simultaneously collected, which show that the local CO dissociation causes a protein quake propagating on different spatial and temporal scales. These studies are important for understanding gas transport and protein deligation processes and shed light on the interplay of active site conformational changes and large-scale protein reorganization.

The Chronus Quantum (ChronusQ) software package is an open source (under the GNU General Public License v2) software infrastructure which targets the solution of challenging problems that arise in ab initio electronic structure theory. Special emphasis is placed on the consistent treatment of time dependence and spin in the electronic wave function, as well as the inclusion of relativistic effects in said treatments. In addition, ChronusQ provides support for the inclusion of uniform finite magnetic fields as external perturbations through the use of gauge‐including atomic orbitals. ChronusQ is a parallel electronic structure code written in modern C++ which utilizes both message passing implementation and shared memory (OpenMP) parallelism. In addition to the examination of the current state of code base itself, a discussion regarding ongoing developments and developer contributions will also be provided.

In the United States the national dialogue around climate change remains enormously polarized, which has stifled support to combat it. This problem persists largely because of the spread of misinformation about the issue. "There is no scientific consensus on climate change." "It still gets cold in the winter, so the planet can't be warming." "Environmental regulations kill jobs and hurt the economy." These statements and others aim to mislead the public and cloud understanding of anthropogenic climate change. This in turn stifles progress toward implementing policies to combat climate change. Misinformation strategies are often funded and implemented by corporations that would be impacted by these new policies. This chapter debunks these myths while also exploring the motivations behind the spreading of misinformation. Strategies to combat this misinformation and enact change are discussed throughout.

The Chronus Quantum (ChronusQ) software package is an open source (under the GNU General Public License v2) software infrastructure which targets the solution of challenging problems that arise in ab initio electronic structure theory. Special emphasis is placed on the consistent treatment of time dependence and spin in the electronic wave function, as well as the inclusion of relativistic effects in said treatments. In addition, ChronusQ provides support for the inclusion of uniform finite magnetic fields as external perturbations through the use of gauge-including atomic orbitals (GIAO). ChronusQ is a parallel electronic structure code written in modern C++ which utilizes both message passing (MPI) and shared memory (OpenMP) parallelism. In addition to the examination of the current state of code base itself, a discussion regarding ongoing developments and developer contributions will also be provided. Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720 Department of Chemistry, University of Washington, Seattle, WA, 98195, United States Corresponding Author: dbwy@lbl.gov Dipartimento di Scienze Chimiche, Universit di Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL 32306, United States Department of Chemistry, Yale University, 225 Prospect Street, New Haven, CT, 06520, United States Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, United States Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, United States Corresponding Author: xsli@uw.edu 1

Dynamic electric properties are most commonly determined by applying linear and nonlinear response theory. This is often a sequential process as each order of response depends on the solution for the previous lower order. Response theory is a perturbative approach and is not directly amenable to modeling time-resolved spectroscopies or experiments involving exotic pulse shapes. Nonperturbative interaction between a system and an electric field can be modeled explicitly in time. This makes it possible to more easily resolve higher-order properties and highly nonlinear processes. Time-dependent configuration interaction has asserted itself as a powerful tool for accurately modeling electronic dynamics. We have implemented time-dependent configuration interaction using the graphical unitary group approach in order to study the dynamics of open-shell systems while retaining spin as a good quantum number. This approach has been used to resolve linear and nonlinear electric properties of molecular systems. Important considerations when modeling dynamic electric properties in the time-domain are presented as well as comparisons to properties of broken symmetry solutions.

X-ray absorption spectroscopy is a powerful technique to probe local electronic and nuclear structure. There has been extensive theoretical work modeling K-edge spectra from first principles. However, modeling L-edge spectra directly with density functional theory poses a unique challenge requiring further study. Spin orbit coupling must be included in the model, and a noncollinear density functional theory is required. Using the real-time exact two-component method, we are able to variationally include one-electron spin orbit coupling terms when calculating the absorption spectrum. The abilities of different basis sets and density functionals to model spectra for both closed- and open-shell systems are investigated using SiCl4 and three transition metal complexes, TiCl4, CrO2Cl2, and [FeCl6](3-). Although we are working in the real-time framework, individual molecular orbital transitions can still be recovered by projecting the density onto the ground state molecular orbital space and separating contributions to the time evolving dipole moment.

Projected Hartree-Fock (PHF) theory can restore important symmetries to broken symmetry wave functions. Variation after projection (VAP) implementations make it possible to deliberately break and then restore a given symmetry by directly minimizing the projected energy expression. This technique can be applied to any symmetry that can be broken from relaxing constraints on single Slater determinant wave functions. For instance, generalized Hartree-Fock (GHF) wave functions are eigenfunctions of neither Ŝz nor S2. By relaxing these constraints, the wave function can explore a larger variational space and can reach lower energies than more constrained HF solutions. We have implemented spin-projected GHF (SGHF), which retains many of the advantages of breaking symmetry while also being a spin eigenfunction, with some notable improvements over previous implementations. Our new algorithm involves the formation of new intermediate matrices not previously discussed in the literature. Discretization of the necessary integration over the rotation group SO(3) is also accomplished much more efficiently using Lebedev grids. A novel scheme to incrementally build rotated Fock matrices is also introduced and compared with more standard approaches.

Real‐time time‐dependent electronic structure theory is one of the most promising methods for investigating time‐dependent molecular responses and electronic dynamics. Since its first modern use in the 1990s, it has been used to study a wide variety of spectroscopic properties and electronic responses to intense external electromagnetic fields, complex environments, and open quantum systems. It has also been used to study molecular conductance, excited state dynamics, ionization, and nonlinear optical effects. Real‐time techniques describe non‐perturbative responses of molecules, allowing for studies that go above and beyond the more traditional energy‐ or frequency‐domain‐based response theories. Recent progress in signal analysis, accurate treatment of environmental responses, relativistic Hamiltonians, and even quantized electromagnetic fields have opened up new avenues of research in time‐dependent molecular response. After discussing the history of real‐time methods, we explore some of the necessary mathematical theory behind the methods, and then survey a wide (yet incomplete) variety of applications for real‐time methods. We then present some brief remarks on the future of real‐time time‐dependent electronic structure theory. WIREs Comput Mol Sci 2018, 8:e1341. doi: 10.1002/wcms.1341 This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods

Materials and molecular systems exhibiting long-lived electronic coherence can facilitate coherent transport, opening the door to efficient charge and energy transport beyond traditional methods. Recently, signatures of a possible coherent, recurrent electronic motion were identified in femtosecond pump-probe spectroscopy experiments on a binuclear platinum complex, where a persistent periodic beating in the transient absorption signal's anisotropy was observed. In this study, we investigate the excitonic dynamics that underlie the suspected electronic coherence for a series of binuclear platinum complexes exhibiting a range of interplatinum distances. Results suggest that the long-lived coherence can only result when competitive electronic couplings are in balance. At longer Pt-Pt distances, the electronic couplings between the two halves of the binuclear system weaken, and exciton localization and recombination is favored on short time scales. For short Pt-Pt distances, electronic couplings between the states in the coherent superposition are stronger than the coupling with other excitonic states, leading to long-lived coherence.

Photoexcited Nickel(II) tetramesitylporphyrin (NiTMP), like many open-shell metalloporphyrins, relaxes rapidly through multiple electronic states following an initial porphyrin-based excitation, some involving metal centered electronic configuration changes that could be harnessed catalytically before excited state relaxation. While a NiTMP excited state present at 100 ps was previously identified by X-ray transient absorption (XTA) spectroscopy at a synchrotron source as a relaxed (d,d) state, the lowest energy excited state (J. Am. Chem. Soc., 2007, 129, 9616 and Chem. Sci., 2010, 1, 642), structural dynamics before thermalization were not resolved due to the ∼100 ps duration of the available X-ray probe pulse. Using the femtosecond (fs) X-ray pulses of the Linac Coherent Light Source (LCLS), the Ni center electronic configuration from the initial excited state to the relaxed (d,d) state has been obtained via ultrafast Ni K-edge XANES (X-ray absorption near edge structure) on a time scale from hundreds of femtoseconds to 100 ps. This enabled the identification of a short-lived Ni(I) species aided by time-dependent density functional theory (TDDFT) methods. Computed electronic and nuclear structure for critical excited electronic states in the relaxation pathway characterize the dependence of the complex's geometry on the electron occupation of the 3d orbitals. Calculated XANES transitions for these excited states assign a short-lived transient signal to the spectroscopic signature of the Ni(I) species, resulting from intramolecular charge transfer on a time scale that has eluded previous synchrotron studies. These combined results enable us to examine the excited state structural dynamics of NiTMP prior to thermal relaxation and to capture intermediates of potential photocatalytic significance.

This report will describe our recent studies of transition metal complex structural dynamics on the fs and ps time scales using an X-ray free electron laser source, Linac Coherent Light Source (LCLS). Ultrafast XANES spectra at the Ni K-edge of nickel(ii) tetramesitylporphyrin (NiTMP) were measured for optically excited states at a timescale from 100 fs to 50 ps, providing insight into its sub-ps electronic and structural relaxation processes. Importantly, a transient reduced state Ni(i) (π, 3dx2-y2) electronic state is captured through the interpretation of a short-lived excited state absorption on the low-energy shoulder of the edge, which is aided by the computation of X-ray transitions for postulated excited electronic states. The observed and computed inner shell to valence orbital transition energies demonstrate and quantify the influence of the electronic configuration on specific metal orbital energies. A strong influence of the valence orbital occupation on the inner shell orbital energies indicates that one should not use the transition energy from 1s to other orbitals to draw conclusions about the d-orbital energies. For photocatalysis, a transient electronic configuration could influence d-orbital energies up to a few eV and any attempt to steer the reaction pathway should account for this to ensure that external energies can be used optimally in driving desirable processes. NiTMP structural evolution and the influence of the porphyrin macrocycle conformation on relaxation kinetics can be likewise inferred from this study.

The interaction between a quantum mechanical system and plane wave light is usually modeled within the electric dipole approximation. This assumes that the intensity of the incident field is constant over the length of the system and transition probabilities are described in terms of the electric dipole transition moment. For short wavelength spectroscopies, such as X-ray absorption, the electric dipole approximation often breaks down. Higher order multipoles are then included to describe transition probabilities. The square of the magnetic dipole and electric quadrupole are often included, but this results in an origin-dependent expression for the oscillator strength. The oscillator strength can be made origin-independent if all terms through the same order in the wave vector are retained. We will show the consequences and potential pitfalls of using either of these two expressions. It is shown that the origin-dependent expression may violate the Thomas-Reiche-Kuhn sum rule and the origin-independent expression can result in negative transition probabilities.

Single-reference techniques based on coupled-cluster (CC) theory, in the forms of linear response (LR) or equation of motion (EOM), are highly accurate and widely used approaches for modeling valence absorption spectra. Unfortunately, these equations with singles and doubles (LR-CCSD and EOM-CCSD) scale as O(N⁶), which may be prohibitively expensive for the study of high-energy excited states using a conventional eigensolver. In this paper, we present an energy-specific non-Hermitian eigensolver that is able to obtain high-energy excited states (e.g., XAS K-edge spectrum) at low computational cost. In addition, we also introduce an improved trial vector for iteratively solving the EOM-CCSD equation with a focus on high-energy eigenstates. The energy-specific EOM-CCSD approach and its low-scaling alternatives are applied to calculations of carbon, nitrogen, oxygen, and sulfur K-edge excitations. The results are compared to other implementations of CCSD for excited states, energy-specific linear response time-dependent density functional theory (TDDFT), and experimental results with multiple statistical metrics are presented and evaluated.

X-ray absorption spectroscopy (XAS) has become a powerful technique in chemical physics, because of advances in synchrotron technology that have greatly improved its temporal and spectroscopic resolution. Our recent work on energy-specific time-dependent density functional theory (ES-TDDFT) allows for the direct calculation of excitation energies in any region of the absorption spectrum, from UV-vis to X-ray. However, the ability of different density functional theories to model X-ray absorption spectra (XAS) of light elements has not yet been verified for ES-TDDFT. This work is a calibration of the ability of existing DFT kernels and basis sets to reproduce experimental K-edge excitation energies. Results were compared against 30 different transitions from gas-phase experiments. We focus on six commonly used density functionals (BHandHLYP, B3LYP, PBE1PBE, BP86, HSE06, LC-ωPBE) and various triple-ζ basis sets. The effects of core and diffuse functions are also investigated.