Improved quasiparticle wave functions and mean field for G(0)W(0) calculations: Initialization with the COHSEX operator

The GW approximation to the electron self-energy has become a standard method for ab initio calculation of excited-state properties of condensed-matter systems. In many calculations, the G W self-energy operator, E, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0 W0 scheme. However, there are known situations in which this diagonal Go Wo approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc-COHSEX-PG W, involves construction of an improved mean field using the static limit of GW, known as COHSEX (Coulomb hole and screened exchange), which is significantly simpler to treat than GW W. In this scheme, frequency-dependent self energy E(N), is constructed and taken to be diagonal in the COHSEX orbitals after the system is solved self-consistently within this formalism. The second method is called off diagonal-COHSEX G W (od-COHSEX-PG W). In this method, one does not self-consistently change the mean-field starting point but diagonalizes the COHSEX Hamiltonian within the Kohn-Sham basis to obtain quasiparticle wave functions and uses the resulting orbitals to construct the G W E in the diagonal form. We apply both methods to a molecular system, silane, and to two bulk systems, Si and Ge under pressure. For silane, both methods give good quasiparticle wave functions and energies. Both methods give good band gaps for bulk silicon and maintain good agreement with experiment. Further, the sc-COHSEX-PGW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

On the positivity of the Jansen-He? operator for arbitrary mass

Iantchenko, A.; Jakuba?a-Amundsen, D.H., (2003) 'On the positivity of the Jansen-He? operator for arbitrary mass', Annales of the Institute Henri Poincar? 4 pp.1083-1099 RAE2008

Low rank subspace clustering (LRSC)

We consider the problem of fitting a union of subspaces to a collection of data points drawn from one or more subspaces and corrupted by noise and/or gross errors. We pose this problem as a non-convex optimization problem, where the goal is to decompose the corrupted data matrix as the sum of a clean and self-expressive dictionary plus a matrix of noise and/or gross errors. By self-expressive we mean a dictionary whose atoms can be expressed as linear combinations of themselves with low-rank coefficients. In the case of noisy data, our key contribution is to show that this non-convex matrix decomposition problem can be solved in closed form from the SVD of the noisy data matrix. The solution involves a novel polynomial thresholding operator on the singular values of the data matrix, which requires minimal shrinkage. For one subspace, a particular case of our framework leads to classical PCA, which requires no shrinkage. For multiple subspaces, the low-rank coefficients obtained by our framework can be used to construct a data affinity matrix from which the clustering of the data according to the subspaces can be obtained by spectral clustering. In the case of data corrupted by gross errors, we solve the problem using an alternating minimization approach, which combines our polynomial thresholding operator with the more traditional shrinkage-thresholding operator. Experiments on motion segmentation and face clustering show that our framework performs on par with state-of-the-art techniques at a reduced computational cost.

On the similarity of Sturm-Liouville operators with non-Hermitian boundary conditions to self-adjoint and normal operators

We consider one-dimensional Schrödinger-type operators in a bounded interval with non-self-adjoint Robin-type boundary conditions. It is well known that such operators are generically conjugate to normal operators via a similarity transformation. Motivated by recent interests in quasi-Hermitian Hamiltonians in quantum mechanics, we study properties of the transformations and similar operators in detail. In the case of parity and time reversal boundary conditions, we establish closed integral-type formulae for the similarity transformations, derive a non-local self-adjoint operator similar to the Schrödinger operator and also find the associated “charge conjugation” operator, which plays the role of fundamental symmetry in a Krein-space reformulation of the problem.

Third order trace formula

In J. Funct. Anal. 257 (2009) 1092-1132, Dykema and Skripka showed the existence of higher order spectral shift functions when the unperturbed self-adjoint operator is bounded and the perturbation is Hilbert-Schmidt. In this article, we give a different proof for the existence of spectral shift function for the third order when the unperturbed operator is self-adjoint (bounded or unbounded, but bounded below).

Generalized rayleigh principle and its applications

In this paper, we mainly deal with cigenvalue problems of non-self-adjoint operator. To begin with, the generalized Rayleigh variational principle, the idea of which was due to Morse and Feshbach, is examined in detail and proved more strictly in mathematics. Then, other three equivalent formulations of it are presented. While applying them to approximate calculation we find the condition under which the above variational method can be identified as the same with Galerkin's one. After that we illustrate the generalized variational principle by considering the hydrodynamic stability of plane Poiseuille flow and Bénard convection. Finally, the Rayleigh quotient method is extended to the cases of non-self-adjoint matrix in order to determine its strong eigenvalne in linear algebra.

Model-based compressive Harmonic-Aware Matching Pursuit: An evaluation

This paper addresses devising a reliable model-based Harmonic-Aware Matching Pursuit (HAMP) for reconstructing sparse harmonic signals from their compressed samples. The performance guarantees of HAMP are provided; they illustrate that the introduced HAMP requires less data measurements and has lower computational cost compared with other greedy techniques. The complexity of formulating a structured sparse approximation algorithm is highlighted and the inapplicability of the conventional thresholding operator to the harmonic signal model is demonstrated. The harmonic sequential deletion algorithm is subsequently proposed and other sparse approximation methods are evaluated. The superior performance of HAMP is depicted in the presented experiments. © 2013 IEEE.

Decay spectrum and decay subspace of normal operators

Let A be a self-adjoint operator on a Hilbert space. It is well known that A admits a unique decomposition into a direct sum of three self-adjoint operators A(p), A(ac) and A(sc) such that there exists an orthonormal basis of eigenvectors for the operator A(p) the operator A(ac) has purely absolutely continuous spectrum and the operator A(sc) has purely singular continuous spectrum. We show the existence of a natural further decomposition of the singular continuous component A c into a direct sum of two self-adjoint operators A(sc)(D) and A(sc)(ND). The corresponding subspaces and spectra are called decaying and purely non-decaying singular subspaces and spectra. Similar decompositions are also shown for unitary operators and for general normal operators.

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).

Nonlinear optics from an ab initio approach by means of the dynamical Berry phase: Application to second- and third-harmonic generation in semiconductors

We present an ab initio real-time-based computational approach to study nonlinear optical properties in condensed matter systems that is especially suitable for crystalline solids and periodic nanostructures. The equations of motion and the coupling of the electrons with the external electric field are derived from the Berry-phase formulation of the dynamical polarization [Souza et al., Phys. Rev. B 69, 085106 (2004)]. Many-body effects are introduced by adding single-particle operators to the independent-particle Hamiltonian. We add a Hartree operator to account for crystal local effects and a scissor operator to correct the independent particle band structure for quasiparticle effects. We also discuss the possibility of accurately treating excitonic effects by adding a screened Hartree-Fock self-energy operator. The approach is validated by calculating the second-harmonic generation of SiC and AlAs bulk semiconductors: an excellent agreement is obtained with existing ab initio calculations from response theory in frequency domain [Luppi et al., Phys. Rev. B 82, 235201 (2010)]. We finally show applications to the second-harmonic generation of CdTe and the third-harmonic generation of Si. 

On the group-velocity property for wave-activity conservation laws

The density and the flux of wave-activity conservation laws are generally required to satisfy the group-velocity property: under the WKB approximation (i.e., for nearly monochromatic small-amplitude waves in a slowly varying medium), the flux divided by the density equals the group velocity. It is shown that this property is automatically satisfied if, under the WKB approximation, the only source of rapid variations in the density and the flux lies in the wave phase. A particular form of the density, based on a self-adjoint operator, is proposed as a systematic choice for a density verifying this condition.

RIEMANN ZETA ZEROS AND PRIME NUMBER SPECTRA IN QUANTUM FIELD THEORY

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Domínios de potências fracionárias de operadores matriciais segundo Lasiecka-Triggiani

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Root system of singular perturbations of the harmonic oscillator type operators

We analyze perturbations of the harmonic oscillator type operators in a Hilbert space H, i.e. of the self-adjoint operator with simple positive eigenvalues μ k satisfying μ k+1 − μ k ≥ Δ > 0. Perturbations are considered in the sense of quadratic forms. Under a local subordination assumption, the eigenvalues of the perturbed operator become eventually simple and the root system contains a Riesz basis.