No Arabic abstract
The deviation of the electron density around the nuclei from spherical symmetry determines the electric field gradient (EFG), which can be measured by various types of spectroscopy. Nuclear Quadrupole Resonance (NQR) is particularly sensitive to the EFG. The EFGs, and by implication NQR frequencies, vary dramatically across materials. Consequently, searching for NQR spectral lines in previously uninvestigated materials represents a major challenge. Calculated EFGs can significantly aid at the search inception. To facilitate this task, we have applied high-throughput density functional theory calculations to predict EFGs for 15187 materials in the JARVIS-DFT database. This database, which will include EFG as a standard entry, is continuously increasing. Given the large scope of the database, it is impractical to verify each calculation. However, we assess accuracy by singling out cases for which reliable experimental information is readily available and compare them to the calculations. We further present a statistical analysis of the results. The database and tools associated with our work are made publicly available by JARVIS-DFT ( https://www.ctcms.nist.gov/~knc6/JVASP.html ) and NIST-JARVIS API ( http://jarvis.nist.gov ).
I summarize Density Functional Theory (DFT) in a language familiar to quantum field theorists, and introduce several apparently novel ideas for constructing {it systematic} approximations for the density functional. I also note that, at least within the large $K$ approximation ($K$ is the number of electron spin components), it is easier to compute the quantum effective action of the Coulomb photon field, which is related to the density functional by algebraic manipulations in momentum space.
We analyze recent measurements [R. Blinc, V. V. Laguta, B. Zalar, M. Itoh and H. Krakauer, J. Phys. : Cond. Mat., v.20, 085204 (2008)] of the electric field gradient on the oxygen site in the perovskites SrTiO3 and BaTiO3, which revealed, in agreement with calculations, a large difference in the EFG for these two compounds. In order to analyze the origin of this difference, we have performed density functional electronic structure calculations within the local-orbital scheme FPLO. Our analysis yields the counter-intuitive behavior that the EFG increases upon lattice expansion. Applying the standard model for perovskites, the effective two-level p-d Hamiltonian, can not explain the observed behavior. In order to describe the EFG dependence correctly, a model beyond this usually sufficient p-d Hamiltonian is needed. We demonstrate that the counter-intuitive increase of the EFG upon lattice expansion can be explained by a s-p-d model, containing the contribution of the oxygen 2s states to the crystal field on the Ti site. The proposed model extension is of general relevance for all related transition metal oxides with similar crystal structure.
We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $alpha$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $alpha$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant).
The methods of density-functional perturbation theory may be used to calculate various physical response properties of insulating crystals including elastic, dielectric, Born charge, and piezoelectric tensors. These and other important tensors may be defined as second derivatives of the total energy with respect to atomic-displacement, electric-field, or strain perturbations, or as mixed derivatives with respect to two of these perturbations. The resulting tensor quantities tend to be coupled in complex ways in polar crystals, giving rise to a variety of variant definitions. For example, it is generally necessary to distinguish between elastic tensors defined under different electrostatic boundary conditions, and between dielectric tensors defined under different elastic boundary conditions. Here, we describe an approach for computing all of these various response tensors in a unified and systematic fashion. Applications are presented for two materials, wurtzite ZnO and rhombohedral BaTiO3, at zero temperature.
In this paper the relationship between the density functional theory of freezing and phase field modeling is examined. More specifically a connection is made between the correlation functions that enter density functional theory and the free energy functionals used in phase field crystal modeling and standard models of binary alloys (i.e., regular solution model). To demonstrate the properties of the phase field crystal formalism a simple model of binary alloy crystallization is derived and shown to simultaneously model solidification, phase segregation, grain growth, elastic and plastic deformations in anisotropic systems with multiple crystal orientations on diffusive time scales.