Electronic properties of interfaces and defects from many-body perturbation theory: Recent developments and applications

We review some recent developments in many body perturbation theory (MBPT) calculations that have enabled the study of interfaces and defects. Starting from the theoretical basis of MBPT, Hedin's equations are presented, leading to the CW and CWI' approximations. We introduce the perturbative approach, that is the one most commonly used for obtaining quasiparticle (QP) energies. The practical strategy presented for dealing with the frequency dependence of the self energy operator is based on either plasmon-pole models (PPM) or the contour deformation technique, with the latter being more accurate. We also discuss the extrapolar method for reducing the number of unoccupied states which need to be included explicity in the calculations. The use of the PAW method in the framework of MBPT is also described. Finally, results which have been obtained using, MBPT for band offsets a interfaces and for defects presented, with companies on the main difficulties and cancels.

Schematic representation of the QP corrections (marked with ) to the band edges (E and E-v) and a defect level (F) for a Si/SiO2 interface (Si and O atoms are represented in blue and red, respectively, in the ball and stick model) with an oxygen vacancy leading to a Si-Si bond (the Si atoms involved in this bond are colored light blue).

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

Density functionals from many-body perturbation theory: The band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby , Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid. (c) 2006 American Institute of Physics.

R-matrix Floquet theory of multiphoton processes:XIII. Resonances in H multiphoton detachment

The results of calculations investigating the effects of autodetaching resonances on the multiphoton detachment spectra of H are presented. The R-matrix Floquet method is used, in which the coupling of the ion with the laser field is described non-perturbatively. The laser field is fixed at an intensity of 10 W cm, while frequency ranges are chosen such that the lowest autodetaching states of the ion are excited through a two- or three-photon transition from the ground state. Detachment rates are compared, where possible, to previous results obtained using perturbation theory. An illustration of how non-lowest-order processes, involving autodetaching states, can lead to light-induced continuum structures is also presented. Finally, it is demonstrated that by using a frequency connecting the 1s and 2s states, the probability of exciting the residual hydrogen atom is significantly enhanced.

Partial and total multiphoton detachment rates for H- in a perturbative B-spline-based approach

An ab initio approach has been applied to study multiphoton detachment rates for the negative hydrogen ion in the lowest nonvanishing order of perturbation theory. The approach is based on the use of B splines allowing an accurate treatment of the electronic repulsion. Total detachment rates have been determined for two- to six-photon processes as well as partial rates for detachment into the different final symmetries. It is shown that B-spline expansions can yield accurate continuum and bound-state wave functions in a very simple manner. The calculated total rates for two- and three-photon detachment are in good agreement with other perturbative calculations. For more than three-photon detachment little information has been available before now. While the total cross sections show little structure, a fair amount of structure is predicted in the partial cross sections. In the two-photon process, it is shown that the detached electrons mainly have s character. For four- and six-photon processes, the contribution from the d channel is the most important. For three- and five-photon processes p electrons dominate the electron emission spectrum. Detachment rates for s and p electrons show minima as a function of photon energy. © 1994 The American Physical Society.

Effect of dipole polarizability on positron binding by strongly polar molecules

A model for positron binding to polar molecules is considered by combining the dipole potential outside the molecule with a strongly repulsive core of a given radius. Using existing experimental data on binding energies leads to unphysically small core radii for all of the molecules studied. This suggests that electron–positron correlations neglected in the simple model play a large role in determining the binding energy. We account for these by including the polarization potential via perturbation theory and non-perturbatively. The perturbative model makes reliable predictions of binding energies for a range of polar organic molecules and hydrogen cyanide. The model also agrees with the linear dependence of the binding energies on the polarizability inferred from the experimental data (Danielson et al 2009 J. Phys. B: At. Mol. Opt. Phys. 42 235203). The effective core radii, however, remain unphysically small for most molecules. Treating molecular polarization non-perturbatively leads to physically meaningful core radii for all of the molecules studied and enables even more accurate predictions of binding energies to be made for nearly all of the molecules considered.

Fundamentals of steady state, non-Newtonian rimming flow

Rimming flow on the inner surface of a horizontal rotating cylinder is investigated. Using a scale analysis, a theoretical description is obtained for steady-state non-Newtonian flow. Simple lubrication theory is applied since the Reynolds number is small and the liquid film is thin. Since the Deborah number is very small the flow is viscometric. The shear-thinning number, which characterizes the shear-thinning effect, may be small or large. A general constitutive law for this kind of flow requires only a single function relating shear stress and shear rate that corresponds to a generalized Newtonian liquid. For this case the run-off condition for rimming flow is derived. Provided the run-off condition is satisfied, the existence of a continuous steady-state solution is proved. The rheological models, which show Newtonian behavior at low shear rates with transition to power-law shear thinning at moderate shear rates, are considered. Numerical results are carried out for the Carreau and Ellis models, which exhibit Newtonian behavior near the free surface and power-law behavior near the wall of the rotating cylinder.

Analytical calculation of far infrared spectra of ice in terms of a molecular model

Wideband far infrared (FIR) spectra of complex permittivity e(p) of ice are calculated in terms of a simple analytical theory based on the method of dipolar autocorrelation functions. The molecular model represents a revision of the model recently presented for liquid water in Adv. Chem. Phys. 127 (2003) 65. A composite two-fractional model is proposed. The model is characterised by three phenomenological potential wells corresponding to the three FIR bands observed in ice. The first fraction comprises dipoles reorienting in a rather narrow and deep hat-like well; these dipoles generate the librational band centred at the frequency approximate to 880 cm(-1). The second fraction comprises elastically interacting particles; they generate two nearby bands placed around frequency 200 cm(-1). For description of one of these bands the harmonic oscillator (HO) model is used, in which translational oscillations of two charged molecules along the H-bond are considered. The other band is produced by the H-bond stretch, which governs hindered rotation of a rigid dipole. Such a motion and its dielectric response are described in terms of a new cut parabolic (CP) model applicable for any vibration amplitude. The composite hat-HO-CP model results in a smooth epsilon(nu) ice spectrum, which does not resemble the noise-like spectra of ice met in the known literature. The proposed theory satisfactorily agrees with the experimental ice spectrum measured at - 7 degrees C. The calculated longitudinal optic-transverse optic (LO-TO) splitting occurring at approximate to 250 cm(-1) qualitatively agrees with the measured data. (c) 2004 Elsevier B.V. All rights reserved.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Single-ionization of helium at Ti:sapphire wavelengths: rates and scaling laws

We present a numerical and theoretical study of intense-field single-electron ionization of helium at 390 nm and 780 nm. Accurate ionization rates (over an intensity range of (0.175-34) X10^14 W/ cm^2 at 390 nm, and (0.275 - 14.4) X 10^14 W /cm^2 at 780 nm) are obtained from full-dimensionality integrations of the time-dependent helium-laser Schroedinger equation. We show that the power law of lowest order perturbation theory, modified with a ponderomotive-shifted ionization potential, is capable of modelling the ionization rates over an intensity range that extends up to two orders of magnitude higher than that applicable to perturbation theory alone. Writing the modified perturbation theory in terms of scaled wavelength and intensity variables, we obtain to first approximation a single ionization law for both the 390 nm and 780 nm cases. To model the data in the high intensity limit as well as in the low, a new function is introduced for the rate. This function has, in part, a resemblance to that derived from tunnelling theory but, importantly, retains the correct frequency-dependence and scaling behaviour derived from the perturbative-like models at lower intensities. Comparison with the predictions of classical ADK tunnelling theory confirms that ADK performs poorly in the frequency and intensity domain treated here.

Choline saccharinate and choline acesulfamate: Ionic liquids with low toxicities

Choline saccharinate and choline acesulfamate are two examples of hydrophilic ionic liquids, which can be prepared from easily available starting materials (choline chloride and a non-nutritive sweetener). The (eco)toxicity of these ionic liquids in aqueous solution is very low in comparison to other types of ionic liquids. A general method for the synthesis and purification of hydrophilic ionic liquids is presented. The method consists of a silver-free metathesis reaction, followed by purification of the ionic liquid by ion-exchange chromatography. The crystal structures show a marked difference in hydrogen bonding between the two ionic liquids, although the saccharinate and the acesulfamate anions show structural similarities. The optimized structures, the energetics, and the charge distribution of cation-anion pairs in the ionic liquids were studied by density functional theory (DFT) and second-order (Moller-Plesset) perturbation theory calculations. The occupation of the non-Lewis orbitals was considered to obtain a qualitative picture of the Lewis structures. The calculated interaction energies and the dipole moments for the ion pairs in the gas phase were discussed.

Three-photon detachment of electrons from the fluorine negative ion

Absolute three-photon detachment cross sections are calculated for the fluorine negative ion within the lowest-order perturbation theory. The Dyson equation of the atomic many-body theory is used to obtain the ground-state 2p wavefunction with correct asymptotic behaviour, corresponding to the true (experimental) binding energy. We show that in accordance with the adiabatic theory this is crucial for obtaining absolute values of the multiphoton cross sections. Comparisons with other calculations and experimental data are presented.

Phase behavior of mixtures of colloidal platelets and nonadsorbing polymers

The phase behavior of a model system of colloidal platelets and nonadsorbing polymers is studied using computer simulations and perturbation theory. The equation of state for the pure platelet reference system is obtained by Monte Carlo simulations, and the free volume fraction accessible to polymers is measured by a trial insertion method. The free volume fraction is also calculated using scaled particle theory. Subsequently, the phase diagram for platelet-polymer mixtures is calculated. For a platelet aspect ratio L/D=0.1 and a polymer to platelet size ratio d/D>0.2, we observe coexistence between two isotropic phases with different densities. For smaller polymers d/D

Calculation of gamma spectra for positron annihilation on molecules

Calculations of gamma spectra for positron annihilation for a selection of molecules, including methane and its fluoro-substitutes, ethane, propane, butane and benzene are presented. The contribution to the ?-spectra from individual molecular orbitals is obtained from electron momentum distributions calculated using the density functional theory (DFT) based B3LYP/TZVP model. For positrons thermalised to room temperature, the calculation, in its simplest form, effectively treats the positron as a plane wave and gives positron annihilation ?-spectra linewidths that are broader (30-40%) than experiment, although the main chemical trends are reproduced. The main physical reason for this is the neglect of positron repulsion from the nuclei. We show that this effect can be incorporated through momentum-dependent correction factors, determined from positron-atom calculations, e.g., many-body perturbation theory. Inclusion of these factors in the calculation gives linewidths that are in improved agreement with experiment.