Search

Nuclear magnetic resonance (NMR) is a powerful spectroscopic technique that is sensitive to the local atomic structure of matter. Computational predictions of NMR parameters can help to interpret experimental data and validate structural models, and machine learning (ML) has emerged as an efficient route to making such predictions. Here, we systematically study graph-neural-network approaches to representing and learning tensor quantities for solid-state NMR -- specifically, the anisotropic magnetic shielding and the electric field gradient. We assess how the numerical accuracy of different ML models translates into prediction quality for experimentally relevant NMR properties: chemical shifts, quadrupolar coupling constants, tensor orientations, and even static 1D spectra. We apply these ML models to a structurally diverse dataset of amorphous SiO$_2$ configurations, spanning a wide range of density and local order, to larger configurations beyond the reach of traditional first-principles methods, and to the dynamics of the $\alpha\unicode{x2013}\beta$ inversion in cristobalite. Our work marks a step toward streamlining ML-driven NMR predictions for both static and dynamic behavior of complex materials, and toward bridging the gap between first-principles modeling and real-world experimental data., Comment: 13 pages, 7 figures

Over the last two decades, following the early developments on maximally localized Wannier functions, an ecosystem of electronic-structure simulation techniques and software packages leveraging the Wannier representation has flourished. This environment includes codes to obtain Wannier functions and interfaces with first-principles simulation software, as well as an increasing number of related post-processing packages. Wannier functions can be obtained for isolated or extended systems (both crystalline and disordered), and can be used to understand chemical bonding, to characterize electric polarization, magnetization, and topology, or as an optimal basis set, providing very accurate interpolations in reciprocal space or large-scale Hamiltonians in real space. In this review, we summarize the current landscape of techniques, materials properties and simulation codes based on Wannier functions that have been made accessible to the research community, and that are now well integrated into what we term a \emph{Wannier function software ecosystem}. First, we introduce the theory and practicalities of Wannier functions, starting from their broad domains of applicability to advanced minimization methods using alternative approaches beyond maximal localization. Then we define the concept of a Wannier ecosystem and its interactions and interoperability with many quantum simulations engines and post-processing packages. We focus on some of the key properties and capabilities that are empowered by such ecosystem\textemdash from band interpolations and large-scale simulations to electronic transport, Berryology, topology, electron-phonon couplings, dynamical mean-field theory, embedding, and Koopmans functionals\textemdash concluding with the current status of interoperability and automation. [...], Comment: Review article, revised version, 60 pages, 12 figures, references in alphabetical order. To appear in Rev. Mod. Phys. https://journals.aps.org/rmp/accepted/e3072EabZ471ee0bd0ff81e5ef7e01214bbaff57a

Structure factors obtained from diffraction experiments are one of the most important quantities for characterizing the electronic and structural properties of materials. Methods for calculating this quantity from plane-wave density functional theory (DFT) codes are typically prohibitively expensive to perform, requiring the electron density to be constructed and evaluated on dense real-space grids. Making use of the projector functions found in both the Vanderbilt ultrasoft pseudopotential and projector augmented wave methods, we implement an approach to calculate structure factors which avoids the use of a dense grid by separating the rapidly changing contributions to the electron density and treating them on logarithmic radial grids. This approach is successfully validated against structure factors obtained from all-electron DFT and experiments for three prototype systems, allowing structure factors to be obtained at all-electron accuracy at a fraction of the cost of previous approaches for plane-wave DFT.

Maximally-localised Wannier functions (MLWFs) are routinely used to compute from first-principles advanced materials properties that require very dense Brillouin zone integration and to build accurate tight-binding models for scale-bridging simulations. At the same time, high-throughput (HT) computational materials design is an emergent field that promises to accelerate the reliable and cost-effective design and optimisation of new materials with target properties. The use of MLWFs in HT workflows has been hampered by the fact that generating MLWFs automatically and robustly without any user intervention and for arbitrary materials is, in general, very challenging. We address this problem directly by proposing a procedure for automatically generating MLWFs for HT frameworks. Our approach is based on the selected columns of the density matrix method (SCDM) and we present the details of its implementation in an AiiDA workflow. We apply our approach to a dataset of 200 bulk crystalline materials that span a wide structural and chemical space. We assess the quality of our MLWFs in terms of the accuracy of the band-structure interpolation that they provide as compared to the band-structure obtained via full first-principles calculations. Finally, we provide a downloadable virtual machine that can be used to reproduce the results of this paper, including all first-principles and atomistic simulations as well as the computational workflows., Comment: In addition to the main manuscript and supplemental materials, we have added in the Materials Cloud entry a dataset with the Wannierized band structures for all 200 materials (which can be downloaded from \url{https://archive.materialscloud.org/record/file?file_id=22842ba6-5528-48d7-9005-daa8d6a32d9d&record_id=425&filename=Vitale-2020-all-bands.pdf})

Phase imaging in electron microscopy is sensitive to the local potential, including charge redistribution from bonding. We demonstrate that electron ptychography provides the necessary sensitivity to detect this subtle effect by directly imaging the charge redistribution in single layer boron nitride. Residual aberrations can be measured and corrected post-collection, and simultaneous atomic number contrast imaging provides unambiguous sub-lattice identification. Density functional theory calculations confirm the detection of charge redistribution.

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials., Comment: 26 pages, 8 figures

We propose modifications to the functional form of the SCAN density functional to eliminate numerical instabilities. This is necessary to allow reliable, automatic generation of pseudopotentials (including PAW potentials). The regularized SCAN is designed to match the original form very closely, and we show that its performance remains comparable.

We extend the Vanderbilt ultrasoft pseudopotential scheme by adding kinetic energy density terms, in order to use meta-GGA exchange potentials, such as the Becke-Johnson or Tran-Blaha potentials, in the planewave-pseudopotential implementation of Density Functional Theory. Having implemented kinetic energy augmentation and non-linear core correction terms in the CASTEP density functional package, we evaluate the validity of our approach by comparing the calculated electronic structure of isolated atoms and semiconductor crystals to all-electron benchmark calculations. Based on our results, we provide recommendations for the practical use of the Tran-Blaha exchange in planewave-pseudopotential codes.

We present a thorough density functional theory study of the magneto-electric (ME) effect in Cr$_{2}$O$_{3}$. The spin-lattice ME tensor $\alpha$ was determined in the low-field and spin flop (SF) phases, using the method of dynamical magnetic charges, and found to be the sum of three distinct components. Two of them, a large relativistic "cycloidal" term and a small longitudinal term, are independent on the spin orientation. The third, only active in the SF phases is also of relativistic origin and arises from magnetic-field-induced chirality, leading to a non-toroidal ME response., Comment: 5 pages with 2 figures, 2 pages supplementary material

The pressure dependence of the Fermi surface of arsenic is examined using the technique of Wannier interpolation, enabling a dense sampling of the Brillouin zone and the ability to capture fine features within it. Focusing primarily on the A7 to simple-cubic phase transition, we find that this semimetal to metal transition is accompanied by the folding of Fermi surfaces. The pressure dependence of the density of states (DOS) of arsenic indicates that the onset of the Peierls-type cubic to rhombohedral distortion is signified by the appearance of emerging van Hove singularities in the DOS, especially around the Fermi level. As we noted in an earlier study, high levels of convergence are consequently required for an accurate description of this transition., Comment: 19 pages, 21 figures, 16 animated gif files

We present a method for the first-principles calculation of nuclear magnetic resonance (NMR) J-coupling in extended systems using state-of-the-art ultrasoft pseudopotentials and including scalar-relativistic effects. The use of ultrasoft pseudopotentials is allowed by extending the projector augmented wave (PAW) method of Joyce et. al [J. Chem. Phys. 127, 204107 (2007)]. We benchmark it against existing local-orbital quantum chemical calculations and experiments for small molecules containing light elements, with good agreement. Scalar-relativistic effects are included at the zeroth-order regular approximation (ZORA) level of theory and benchmarked against existing local-orbital quantum chemical calculations and experiments for a number of small molecules containing the heavy row six elements W, Pt, Hg, Tl, and Pb, with good agreement. Finally, 1J(P-Ag) and 2J(P-Ag-P) couplings are calculated in some larger molecular crystals and compared against solid-state NMR experiments. Some remarks are also made as to improving the numerical stability of dipole perturbations using PAW., Comment: 18 pages, 7 figures, 10 tables. The following article has been submitted to The Journal of Chemical Physics. After it is published, it will be found at http://jcp.aip.org

The electronic ground state of a periodic system is usually described in terms of extended Bloch orbitals, but an alternative representation in terms of localized "Wannier functions" was introduced by Gregory Wannier in 1937. The connection between the Bloch and Wannier representations is realized by families of transformations in a continuous space of unitary matrices, carrying a large degree of arbitrariness. Since 1997, methods have been developed that allow one to iteratively transform the extended Bloch orbitals of a first-principles calculation into a unique set of maximally localized Wannier functions, accomplishing the solid-state equivalent of constructing localized molecular orbitals, or "Boys orbitals" as previously known from the chemistry literature. These developments are reviewed here, and a survey of the applications of these methods is presented. This latter includes a description of their use in analyzing the nature of chemical bonding, or as a local probe of phenomena related to electric polarization and orbital magnetization. Wannier interpolation schemes are also reviewed, by which quantities computed on a coarse reciprocal-space mesh can be used to interpolate onto much finer meshes at low cost, and applications in which Wannier functions are used as efficient basis functions are discussed. Finally the construction and use of Wannier functions outside the context of electronic-structure theory is presented, for cases that include phonon excitations, photonic crystals, and cold-atom optical lattices., Comment: 62 pages. Accepted for publication in Reviews of Modern Physics

We report on our investigation of the crystal structure of arsenic under compression, focusing primarily on the pressure-induced A7 to simple cubic (sc) phase transition. The two-atom rhombohedral unit cell is subjected to pressures ranging from 0 GPa to 200 GPa; for each given pressure, cell lengths and angles, as well as atomic positions, are allowed to vary until the fully relaxed structure is obtained. We find that the nearest and next-nearest neighbor distances give the clearest indication of the occurrence of a structural phase transition. Calculations are performed using the local density approximation (LDA) and the PBE and PW91 generalized gradient approximations (GGA-PBE and GGA-PW91) for the exchange-correlation functional. The A7 to sc transition is found to occur at 21+/-1 GPa in the LDA, at 28+/-1 GPa in the GGA-PBE and at 29+/-1 GPa in the GGA-PW91; no volume discontinuity is observed across the transition in any of the three cases. We use k-point grids as dense as 66X66X66 to enable us to present reliably converged results for the A7 to sc transition of arsenic., Comment: To be published in Physical Review B; material supplementary to this article is available at arXiv:0810.1693

The material presented here is supplementary to the article entitled: "Density-functional investigation of the rhombohedral to simple cubic phase transition of arsenic"; it deals with the convergence issues involved in studying a semi-metal to metal structural phase transition such as the A7 to sc transition of arsenic., Comment: This material is supplementary to arXiv:0810.1692

A method to calculate NMR J-coupling constants from first principles in extended systems is presented. It is based on density functional theory and is formulated within a planewave-pseudopotential framework. The all-electron properties are recovered using the projector augmented wave approach. The method is validated by comparison with existing quantum chemical calculations of solution-state systems and with experimental data. The approach has been applied to verify measured J-coupling in a silicophosphate structure, Si5O(PO4)6, Comment: 9 pages

While the intrinsic anomalous Hall conductivity is normally written in terms of an integral of the electronic Berry curvature over the occupied portions of the Brillouin zone, Haldane has recently pointed out that this quantity (or more precisely, its ``non-quantized part'') may alternatively be expressed as a Fermi-surface property. Here we present an {\it ab-initio} approach for computing the anomalous Hall conductivity that takes advantage of this observation by converting the integral over the Fermi sea into a more efficient integral on the Fermi surface only. First, a conventional electronic-structure calculation is performed with spin-orbit interaction included. Maximally-localized Wannier functions are then constructed by a post-processing step in order to convert the {\it ab-initio} electronic structure around the Fermi level into a tight-binding-like form. Working in the Wannier representation, the Brillouin zone is sampled on a large number of equally spaced parallel slices oriented normal to the total magnetization. On each slice, we find the intersections of the Fermi-surface sheets with the slice by standard contour methods, organize these into a set of closed loops, and compute the Berry phases of the Bloch states as they are transported around these loops. The anomalous Hall conductivity is proportional to the sum of the Berry phases of all the loops on all the slices. Illustrative calculations are performed for Fe, Co and Ni., Comment: 12 pages, 9 figures

We present an efficient first-principles approach for calculating Fermi surface averages and spectral properties of solids, and use it to compute the low-field Hall coefficient of several cubic metals and the magnetic circular dichroism of iron. The first step is to perform a conventional first-principles calculation and store the low-lying Bloch functions evaluated on a uniform grid of k-points in the Brillouin zone. We then map those states onto a set of maximally-localized Wannier functions, and evaluate the matrix elements of the Hamiltonian and the other needed operators between the Wannier orbitals, thus setting up an ``exact tight-binding model.'' In this compact representation the k-space quantities are evaluated inexpensively using a generalized Slater-Koster interpolation. Because of the strong localization of the Wannier orbitals in real space, the smoothness and accuracy of the k-space interpolation increases rapidly with the number of grid points originally used to construct the Wannier functions. This allows k-space integrals to be performed with ab-initio accuracy at low cost. In the Wannier representation, band gradients, effective masses, and other k-derivatives needed for transport and optical coefficients can be evaluated analytically, producing numerically stable results even at band crossings and near weak avoided crossings., Comment: 12 pages, 7 figures

We compute from first principles the infrared dispersion of the nonlinear susceptibility $\chi^{(2)}$ in zincblende semiconductors. At terahertz frequencies the nonlinear susceptibility depends not only on the purely electronic response $\chi^{(2)}_{\infty}$, but also on three other parameters $C_1$, $C_2$ and $C_3$ describing the contributions from ionic motion. They relate to the TO Raman polarizability, the second-order displacement-induced dielectric polarization, and the third-order lattice potential. Contrary to previous theory, we find that mechanical anharmonicity ($C_3$) dominates over electrical anharmonicity ($C_2$), which is consistent with recent experiments on GaAs. We predict that the sharp minimum in the intensity of second-harmonic generation recently observed for GaAs between $\omega_{\rm TO}/2$ and $\omega_{\rm TO}$ does not occur for several other III-V compounds., Comment: 9 pages, 3 figures; updated bibliography

The intrinsic anomalous Hall effect in ferromagnets depends on subtle spin-orbit-induced effects in the electronic structure, and recent ab-initio studies found that it was necessary to sample the Brillouin zone at millions of k-points to converge the calculation. We present an efficient first-principles approach for computing the anomalous Hall conductivity. We start out by performing a conventional electronic-structure calculation including spin-orbit coupling on a uniform and relatively coarse k-point mesh. From the resulting Bloch states, maximally-localized Wannier functions are constructed which reproduce the ab-initio states up to the Fermi level. The Hamiltonian and position-operator matrix elements, needed to represent the energy bands and Berry curvatures, are then set up between the Wannier orbitals. This completes the first stage of the calculation, whereby the low-energy ab-initio problem is transformed into an effective tight-binding form. The second stage only involves Fourier transforms and unitary transformations of the small matrices set up in the first stage. With these inexpensive operations, the quantities of interest are interpolated onto a dense k-point mesh and used to evaluate the anomalous Hall conductivity as a Brillouin zone integral. The present scheme, which also avoids the cumbersome summation over all unoccupied states in the Kubo formula, is applied to bcc Fe, giving excellent agreement with conventional, less efficient first-principles calculations. Remarkably, we find that more than 99% of the effect can be recovered by keeping a set of terms depending only on the Hamiltonian matrix elements, not on matrix elements of the position operator., Comment: 16 pages, 7 figures