No Arabic abstract
Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Greens function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW+f$_{xc}$). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that the different MD approaches lead to up to 1 eV change in the average band gap, thus, we focused on the band gap dependence with the geometrical properties of the system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance, and HOH angle, as well as the intermolecular hydrogen bond O$cdotcdotcdot$O distance, and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O$cdotcdotcdot$O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows a good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, a broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.
Koopmans-compliant (KC) functionals have been shown to provide accurate spectral properties through a generalized condition of piece-wise linearity of the total energy as a function of the fractional addition/removal of an electron to/from any orbital. We analyze the performance of different KC functionals on the GW100 test-set, comparing the ionization potentials (as opposite of the energy of the highest occupied orbital) of these 100 molecules to those obtained from CCSD(T) total energy differences, and experimental results, finding excellent agreement with a mean absolute error of 0.20 eV for the KIPZ functional, that is state-of-the-art for both DFT-based calculations and many-body perturbation theory. We highlight similarities and differences between KC functionals and other electronic-structure approaches, such as dielectric-dependent hybrid functionals and G$_0$W$_0$, both from a theoretical and from a practical point of view, arguing that Koopmans-compliant potentials can be considered as a local and orbital-dependent counterpart to the electronic GW self-energy, albeit already including approximate vertex corrections.
The performance of density functional theory depends largely on the approximation applied for the exchange functional. We propose here a novel screened exchange potential for semiconductors, with parameters based on the physical properties of the underlying microscopic screening and obeying the requirements for proper asymptotic behavior. We demonstrate that this functional is Koopmans-compliant and reproduces a wide range of band gaps. We also show, that the only tunable parameter of the functional can be kept constant upon changing the cation or the anion isovalently, making the approach suitable for treating alloys.
We introduce a spectral density functional theory which can be used to compute energetics and spectra of real strongly--correlated materials using methods, algorithms and computer programs of the electronic structure theory of solids. The approach considers the total free energy of a system as a functional of a local electronic Green function which is probed in the region of interest. Since we have a variety of notions of locality in our formulation, our method is manifestly basis--set dependent. However, it produces the exact total energy and local excitational spectrum provided that the exact functional is extremized. The self--energy of the theory appears as an auxiliary mass operator similar to the introduction of the ground--state Kohn--Sham potential in density functional theory. It is automatically short--ranged in the same region of Hilbert space which defines the local Green function. We exploit this property to find good approximations to the functional. For example, if electronic self--energy is known to be local in some portion of Hilbert space, a good approximation to the functional is provided by the corresponding local dynamical mean--field theory. A simplified implementation of the theory is described based on the linear muffin--tin orbital method widely used in electronic strucure calculations. We demonstrate the power of the approach on the long--standing problem of the anomalous volume expansion of metallic plutonium.
We consider a minimization scheme based on the Householder transport operator for the Grassman manifold, where a point on the manifold is represented by a m x n matrix with orthonormal columns. In particular, we consider the case where m >> n and present a method with asymptotic complexity mn^2. To avoid explicit parametrization of the manifold we use Householder transforms to move on the manifold, and present a formulation for simultaneous Householder reflections for S-orthonormal columns. We compare a quasi-Newton and nonlinear conjugate gradient implementation adapted to the manifold with a projected nonlinear conjugate gradient method, and demonstrate that the convergence rate is significantly improved if the manifold is taken into account when designing the optimization procedure.
We consider a direct optimization approach for ensemble density functional theory electronic structure calculations. The update operator for the electronic orbitals takes the structure of the Stiefel manifold into account and we present an optimization scheme for the occupation numbers that ensures that the constraints remain satisfied. We also compare sequential and simultaneous quasi-Newton and nonlinear conjugate gradient optimization procedures, and demonstrate that simultaneous optimization of the electronic orbitals and occupation numbers improve performance compared to the sequential approach.