Reconstruction of frozen-core all-electron orbitals from pseudo-orbitals

All-Electron BSE@GW Method for K-Edge Core Electron Excitation Energies

We present an accurate computational approach to calculate absolute K-edge core electron excitation energies as measured by X-ray absorption spectroscopy. Our approach employs an all-electron Bethe-Salpeter equation (BSE) formalism based on GW quasiparticle energies (BSE@GW) using numeric atom-centered orbitals (NAOs). The BSE@GW method has become an increasingly popular method for the computation of neutral valence excitation energies of molecules. However, it was so far not applied to molecular K-edge excitation energies. We discuss the influence of different numerical approximations on the BSE@GW calculation and employ in our final setup (i) exact numeric algorithms for the frequency integration of the GW self-energy, (ii) G₀W₀ and BSE starting points with ∼50% of exact exchange, (iii) the Tamm-Dancoff approximation and (iv) relativistic corrections. We study the basis set dependence and convergence with common Gaussian-type orbital and NAO basis sets. We identify the importance of additional spatially confined basis functions as well as of diffuse augmenting basis functions. The accuracy of our BSE@GW method is assessed for a benchmark set of small organic molecules, previously used for benchmarking the equation-of-motion coupled cluster method [Peng et al., J. Chem. Theory Comput., 2015, 11, 4146], as well as the medium-sized dibenzothiophene (DBT) molecule. Our BSE@GW results for absolute excitation energies are in excellent agreement with the experiment, with a mean average error of only 0.63 eV for the benchmark set and with errors <1 eV for the DBT molecule.

A classical treatment of optical tunneling in plasmonic gaps: extending the quantum corrected model to practical situations

The optical response of plasmonic nanogaps is challenging to address when the separation between the two nanoparticles forming the gap is reduced to a few nanometers or even subnanometer distances. We have compared results of the plasmon response within different levels of approximation, and identified a classical local regime, a nonlocal regime and a quantum regime of interaction. For separations of a few Ångstroms, in the quantum regime, optical tunneling can occur, strongly modifying the optics of the nanogap. We have considered a classical effective model, so called Quantum Corrected Model (QCM), that has been introduced to correctly describe the main features of optical transport in plasmonic nanogaps. The basics of this model are explained in detail, and its implementation is extended to include nonlocal effects and address practical situations involving different materials and temperatures of operation.

Hyperfine interactions at nitrogen interstitial defects in diamond

Diamond has many extreme physical properties and it can be used in a wide range of applications. In particular it is a highly effective particle detection material, where radiation damage is an important consideration. The WAR9 and WAR10 are electron paramagnetic resonance centres seen in irradiated, nitrogen-containing diamond. These S = 1/2 defects have C(2v) and C(1h) symmetry, respectively, and the experimental spectra have been interpreted as arising from nitrogen split-interstitial centres. Based upon the experimental and theoretical understanding of interstitial nitrogen defect structures, the AIMPRO density functional code has been used to assess the assignments for the structures of WAR9 and WAR10. Although the calculated hyperfine interaction tensors are consistent with the measured values for WAR9, the thermal stability renders the assignment problematic. The model for the WAR10 centre yields principal directions of the hyperfine tensor at variance with observation. Alternative models for both centres are discussed in this paper, but no convincing structures have been found.

Inner-shell spectroscopy by the Gaussian and augmented plane wave method

We present an approach for calculating near-edge X-ray absorption spectra at the density functional theory level, which is suited for condensed matter simulations. The method is based on the standard solution of the all-electron KS equations with a modified core-hole potential, which reproduces the relaxation of the orbitals induced by the promotion of the core electron to an unoccupied valence level. The all-electron description of the charge density is based on the Gaussian and augmented plane wave formalism. The reliability of the proposed method is assessed by comparing the computed spectra of some small molecules in the gas phase to the experimental spectra reported in literature. The sensitivity of the computed spectra to the local environment, i.e. the specific bonds formed by the absorbing atom or the presence of hydrogen bonds, open promising perspective for this technique as a predictive tool in the investigation of a more complex system of an unknown structure. The straightforward extension of the method to condensed matter is demonstrated by the calculation of the C K-edge in diamond.

First-principle molecular dynamics with ultrasoft pseudopotentials: parallel implementation and application to extended bioinorganic systems

We present a plane-wave ultrasoft pseudopotential implementation of first-principle molecular dynamics, which is well suited to model large molecular systems containing transition metal centers. We describe an efficient strategy for parallelization that includes special features to deal with the augmented charge in the contest of Vanderbilt's ultrasoft pseudopotentials. We also discuss a simple approach to model molecular systems with a net charge and/or large dipole/quadrupole moments. We present test applications to manganese and iron porphyrins representative of a large class of biologically relevant metalorganic systems. Our results show that accurate density-functional theory calculations on systems with several hundred atoms are feasible with access to moderate computational resources.

Importance of quantum tunneling in vacancy-hydrogen complexes in diamond

Our ab initio calculations of the hyperfine parameters for negatively charged vacancy-hydrogen and nitrogen-vacancy-hydrogen complexes in diamond compare static defect models and models which account for the quantum tunneling behavior of hydrogen. The static models give rise to hyperfine splittings that are inconsistent with the experimental electron paramagnetic resonance data. In contrast, the hyperfine parameters for the quantum dynamical models are in agreement with the experimental observations. We show that the quantum motion of the proton is crucial to the prediction of symmetry and hyperfine constants for two simple defect centers in diamond. Static a priori methods fail for these systems.

X-ray absorption spectra of water within a plane-wave Car-Parrinello molecular dynamics framework

We describe the implementation of a simple technique to simulate core-level spectra within the Car-Parrinello plane-waves molecular dynamics framework. The x-ray absorption (XA) spectra are generated using the transition potential technique with the effect of the core hole included through a specifically developed pseudopotential for the core-excited atom. Despite the lack of 1s core orbitals in the pseudopotential treatment, the required transition moments are accurately calculated without reconstruction of the all-electron orbitals. The method is applied to the oxygen XA spectra of water in its various aggregation states, but it is transferable to any first-row element. The computed spectra are compared favorably with the results from all-electron cluster calculations, as well as with experimental data. The periodicity of the plane-wave technique improves the description of condensed phases. The molecular dynamics simulation enables in principle a proper treatment of thermal effects and dynamical averaging in complex systems.

Calculation of near-edge x-ray-absorption fine structure at finite temperatures: Spectral signatures of hydrogen bond breaking in liquid water

We calculate the near-edge x-ray-absorption fine structure of H(2)O in the gas, hexagonal ice, and liquid phases using heuristic density-functional based methods. We present a detailed comparison of our results with experiment. The differences between the ice and water spectra can be rationalized in terms of the breaking of hydrogen bonds around the absorbing molecule. In particular the increase in the pre-edge absorption feature from ice to water is shown to be due to the breaking of a donor hydrogen bond. We also find that in water approximately 19% of hydrogen bonds are broken.

Concentration of small ring structures in vitreous silica from a first-principles analysis of the Raman spectrum

Using a first-principles approach, we calculate Raman spectra for a model structure of vitreous silica. We develop a perturbational method for calculating the dielectric tensor in an ultrasoft pseudopotential scheme and obtain Raman coupling tensors by finite differences with respect to atomic displacements. For frequencies below 1000 cm(-1), the parallel-polarized Raman spectrum of vitreous silica is dominated by oxygen bending motions, showing a strong sensitivity to the intermediate range structure. By modeling the Raman coupling, we derive estimates for the concentrations of three- and four-membered rings from the experimental intensities of the Raman defect lines.