The performance of Gaussian basis sets for density functional theory-based calculations of core-electron spectroscopies is assessed. The convergence of core-electron binding energies and core-excitation energies using a range of basis sets including split-valence, correlation-consistent, polarisation-consistent and individual gauge for localised orbitals basis sets is studied. For self-consistent field calculations of core-electron binding energies and core-excitation energies of first-row elements, relatively small basis sets can accurately reproduce the values of much larger basis sets, with the IGLO basis sets performing particularly well. Calculations for the K-edge of second-row elements are more challenging, and of the smaller basis sets, pcSseg-2 has the best performance. For the correlation-consistent basis sets, inclusion of core-valence correlation functions is important, with the cc-pCVTZ basis set giving accurate results. Time-dependent density functional theory-based calculations of core-excitation energies show less sensitivity to the basis set with relatively small basis sets, such as pcSseg-1 or pcSseg-2, reproducing the values for much larger basis sets accurately. In contrast, time-dependent density functional theory calculations of X-ray emission energies are highly dependent on the basis set, but the IGLO-II, IGLO-III and pcSseg-2 basis sets provide a good level of accuracy.

Quantum dynamics simulations are used to simulate the ultrafast X-ray Absorption Near-Edge Structure (XANES) spectra of photoexcited pyrazine including two strongly coupled electronically excited states and four normal mode degrees of freedom.

Resonant inelastic soft X-ray scattering maps for the water molecule are simulated by combining quantum chemical calculations of X-ray spectroscopy with ab initio molecular dynamics. The resonant inelastic scattering intensity is computed using the Kramers-Heisenberg formalism, which accounts for channel interference and polarization anisotropy. Algebraic diagrammatic construction and density functional theory-based approaches for the calculation of the X-ray transition energies and transition dipole moments of the absorption and emission processes are explored. Conformational sampling of both ground and core-excited intermediate states allows the effects of ultrafast dynamics on the computed maps to be studied. Overall, it is shown how resonant inelastic scattering maps can be simulated with a computationally efficient protocol that can be extended to investigate larger systems.

Self-consistent field methods for excited states offer an attractive low-cost route to study not only excitation energies but also properties of excited states. Here, we present the generalization of two selfconsistent field methods, the maximum overlap method (MOM) and the σ-SCF method, to calculate excited states in strong magnetic fields and investigate their stability and accuracy in this context. These methods use different strategies to overcome the well-known variational collapse of energybased optimizations to the lowest solution of a given symmetry. The MOM tackles this problem in the definition of the orbital occupations to constrain the self-consistent field procedure to converge on excited states, while the σ-SCF method is based on the minimization of the variance instead of the energy. To overcome the high computational cost of the variance minimization, we present a new implementation of the σ-SCF method with the resolution of identity approximation, allowing the use of large basis sets, which is an important requirement for calculations in strong magnetic fields. The accuracy of these methods is assessed by comparison with the benchmark literature data for He, H 2 , and CH + . The results reveal severe limitations of the variance-based scheme, which become more acute in large basis sets. In particular, many states are not accessible using variance optimization. Detailed analysis shows that this is a general feature of variance optimization approaches due to the masking of local minima in the optimization. In contrast, the MOM shows promising performance for computing excited states under these conditions, yielding results consistent with available benchmark data for a diverse range of electronic states.

We have examined how the self-interaction error in density-functional theory (DFT) calculations affects energies calculated on large systems (600-1000 atoms) involving several charged groups. We employ 18 different quantum mechanical (QM) methods, including Hartree-Fock, as well as pure, hybrid, and range-separated DFT methods. They are used to calculate reaction and activation energies for three different protein models in vacuum, in a point-charge surrounding, or with a continuum-solvent model. We show that pure DFT functionals give rise to a significant delocalization of the charges in charged groups in the protein, typically by ∼0.1 e, as evidenced from the Mulliken charges. This has a clear effect on how the surroundings affect calculated reaction and activation energies, indicating that these methods should be avoided for DFT calculations on large systems. Fortunately, methods such as CAM-B3LYP, BHLYP, and M06-2X give results that agree within a few kilojoules per mole, especially when the calculations are performed in a point-charge surrounding. Therefore, we recommend these methods to estimate the effect of the surroundings with large QM systems (but other QM methods may be used to study the intrinsic reaction and activation energies).

The nitrogen K-edge resonant inelastic x-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft x-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s→π* and 1s→Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s→π* and π→π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.

Glyoxalase I (GlxI) is a member of the glyoxalase system, which is important in cell detoxification and converts hemithioacetals of methylglyoxal (a cytotoxic byproduct of sugar metabolism that may react with DNA or proteins and introduce nucleic acid strand breaks, elevated mutation frequencies and structural or functional changes of the proteins) and glutathione into Dlactate. GlxI accepts both the S and R enantiomers of hemithioacetal, but converts them to only the S-D enantiomer of lactoylglutathione. Interestingly, the enzyme shows this unusual specificity with a rather symmetric active site (a Zn ion coordinated to two glutamate residues; Glu-99 and Glu-172), making the investigation of its reaction mechanism challenging. Herein, we have performed a series of combined quantum mechanics and molecular mechanics calculations to study the reaction mechanism of GlxI. The substrate can bind to the enzyme in two different modes, depending on the direction of its alcoholic proton (H2; toward Glu-99 or Glu-172). Our results show that the S substrate can react only if H2 is directed toward Glu-99 and the R substrate only if H2 is directed toward Glu-172. In both cases, the reactions lead to the experimentally observed S-D enantiomer of the product. In addition, the results do not show any low-energy paths to the wrong enantiomer of the product from neither the S nor the R substrate. Previous studies have presented several opposing mechanisms for the conversion of R and S enantiomers of the substrate to the correct enantiomer of the product. Our results confirm one of them for the S substrate, but propose a new one for the R substrate.

The charge distribution associated with individual components in functionalised ionic liquids (ILs) can be tuned by careful manipulation of the substituent groups incorporated into the ions.