Magnetic susceptibility of insulators from first principles

Comparing GIPAW with numerically exact chemical shieldings: The role of two-center contributions to the induced current

In this study, we benchmark density functional theory gauge-including projector-augmented-wave (GIPAW) chemical shieldings against molecular shieldings for which basis set completeness has been achieved [Jensen et al., Phys. Chem. Chem. Phys. 18, 21145 (2016)]. We demonstrate the importance of two-center corrections for GIPAW hydrogen shieldings. For the other nuclei studied, standard GIPAW is sufficiently accurate. We find that GIPAW can be pushed to closely approach the basis set limit. The only source of small inaccuracies lies in the contribution to the shielding that is caused by surface currents, which we estimate comparing GIPAW susceptibilities to converged molecular magnetizabilities.

First-Principles Calculation of the Optical Rotatory Power of Periodic Systems: Application on α-Quartz, Tartaric Acid Crystal, and Chiral (n,m)-Carbon Nanotubes

The self-consistent coupled-perturbed (SC-CP) method in the CRYSTAL program has been adapted to obtain electromagnetic optical rotation properties of chiral periodic systems based on the calculation of the magnetic moment induced by the electric field. Toward that end, an expression for the magnetic transition moment is developed, which involves an appropriate electronic angular momentum operator. This operator is forced to be hermitian so that the chiroptical properties are real. In our formulation, the trace of the optical rotatory power matrix is gauge-origin-invariant as long as the electric dipole transition matrix elements are obtained using the velocity (rather than position) operator. On the other hand, the component along the optic axis is invariant in general for uniaxial and biaxial crystals. Under the same conditions, these properties also do not depend on the so-called missing integers that occur in the treatment of the electric dipole moment of quasi-one-dimensional periodic systems or the analogue of missing integers for the case of higher dimensionality. Tests on a model H₂O₂ polymer confirm the formalism and, as desired, show that the calculated properties are independent of the size and definition of the unit cell. In addition, an empirical relation to a finite oligomer gauge-including atomic orbital (GIAO) calculation is found. Applications, with comparison to experiment, are carried for α-quartz, tartaric acid crystal, and carbon nanotubes. Future developments of this initial approach to chiroptical properties in the solid state are noted.

First-principles calculation of 11B solid-state NMR parameters of boron-rich compounds I: the rhombohedral boron modifications and B12X2 (X = P, As, O)

In the present study, solid-state nuclear magnetic resonance (NMR) spectra under magic angle spinning conditions of the rhombohedral structures α-B and B12P2 are reported together with the corresponding parameter sets from first principles calculations on α-B B12X2 (X = P, As, O). With the combination of density functional theory (DFT) and the gauge-including projector-augmented wave (GIPAW) approach as the theoretical tools at hand the computed 11B parameters lead to unambiguous explanation of the measurements. Thereby, we overcome common obstacles of processing recorded NMR spectra of solid-state compounds with several crystallographic positions, in particular non-trivial signal assignments and parameter determination due to peak overlap or even unexpected intensity/area ratios. In fact, we find very good agreement between the theoretical results and measured spectra without applying fitting procedures. Using the Perdew-Burke-Ernzerhof (PBE) functional, the results of the common construction types for pseudopotentials and referencing methods for the chemical shift determination are compared. Suggestions and conclusions from experimental 11B NMR studies on parameters according to the icosahedral positions are critically discussed, for instance the early suspected correlation to chemical shifts is not confirmed. Regarding the electric field gradient (EFG) a detailed explanation for obtaining small deviations amongst all investigated structures of the icosahedral polar sites compared to the equatorial sites is given. Our results show an important link between the exohedral bonding situation of compounds with icosahedral structure elements and the main axis of the EFG and therefore, also measurable quadrupole coupling constants if certain geometrical conditions are fulfilled. Finally, this work also contributes to establishing the number of unique sites measured by solid-state NMR methods within the modification of β-B.

Kagome van-der-Waals Pd3P2S8 with flat band

With the advanced investigations into low-dimensional systems, it has become essential to find materials having interesting lattices that can be exfoliated down to monolayer. One particular important structure is a kagome lattice with its potentially diverse and vibrant physics. We report a van-der-Waals kagome lattice material, Pd₃P₂S_8, with several unique properties such as an intriguing flat band. The flat band is shown to arise from a possible compact-localized state of all five 4d orbitals of Pd. The diamagnetic susceptibility is precisely measured to support the calculated susceptibility obtained from the band structure. We further demonstrate that Pd₃P₂S₈ can be exfoliated down to monolayer, which ultimately will allow the possible control of the localized states in this two-dimensional kagome lattice using the electric field gating.

Numerical study of disorder on the orbital magnetization in two dimensions

The modern theory of orbital magnetization (OM) was developed by using Wannier function method, which has a formalism similar with the Berry phase. In this manuscript, we perform a numerical study on the fate of the OM under disorder, by using this method on the Haldane model in two dimensions, which can be tuned between a normal insulator or a Chern insulator at half filling. The effects of increasing disorder on OM for both cases are simulated. Energy renormalization shifts are observed in the weak disorder regime and topologically trivial case, which was predicted by a self-consistentT-matrix approximation. Besides this, two other phenomena can be seen. One is the localization trend of the band orbital magnetization. The other is the remarkable contribution from topological chiral states arising from nonzero Chern number or large value of integrated Berry curvature. If the fermi energy is fixed at the gap center of the clean system, there is an enhancement of |M| at the intermediate disorder, for both cases of normal and Chern insulators, which can be attributed to the disorder induced topological metal state before localization.

Calculation of Diamagnetic Susceptibility Tensors of Organic Crystals: From Coronene to Pharmaceutical Polymorphs

Understanding why crystallization in strong magnetic fields can lead to new polymorphs requires methods to calculate the diamagnetic response of organic molecular crystals. We develop the calculation of the macroscopic diamagnetic susceptibility tensor, χ^cryst, for organic molecular crystals using periodic density functional methods. The crystal magnetic susceptibility tensor, χ^cryst, for all experimentally known polymorphs, and its molecular counterpart, χ^mol, are calculated for flexible pharmaceuticals such as carbamazepine, flufenamic acid, and chalcones, and rigid molecules, such as benzene, pyridine, acridine, anthracene, and coronene, whose molecular magnetic properties have been traditionally studied. A tensor addition method is developed to approximate the crystal diamagnetic susceptibility tensor, χ^cryst, from the molecular one, χ^mol, giving good agreement with those calculated directly using the more costly periodic density functional method for χ^cryst. The response of pharmaceutical molecules and crystals to magnetic fields, as embodied by χ^cryst, is largely determined by the packing in the crystal, as well as the molecular conformation. The anisotropy of χ^cryst can vary considerably between polymorphs though the isotropic terms are fairly constant. The implications for developing a computational method for predicting whether crystallization in a magnetic field could produce a novel or different polymorph are discussed.

Electronic orbital response of regular extended and infinite periodic systems to magnetic fields. I. Theoretical foundations for static case

A theoretical treatment for the orbital response of an infinite, periodic system to a static, homogeneous, magnetic field is presented. It is assumed that the system of interest has an energy gap separating occupied and unoccupied orbitals and a zero Chern number. In contrast to earlier studies, we do not utilize a perturbation expansion, although we do assume the field is sufficiently weak that the occurrence of Landau levels can be ignored. The theory is developed by analyzing results for large, finite systems and also by comparing with the analogous treatment of an electrostatic field. The resulting many-electron Hamilton operator is forced to be hermitian, but hermiticity is not preserved, in general, for the subsequently derived single-particle operators that determine the electronic orbitals. However, we demonstrate that when focusing on the canonical solutions to the single-particle equations, hermiticity is preserved. The issue of gauge-origin dependence of approximate solutions is addressed. Our approach is compared with several previously proposed treatments, whereby limitations in some of the latter are identified.

The Calculation of NMR Chemical Shifts in Periodic Systems Based on Gauge Including Atomic Orbitals and Density Functional Theory

We present here a method that can calculate NMR shielding tensors from first principles for systems with translational invariance. Our approach is based on Kohn-Sham density functional theory and gauge-including atomic orbitals. Our scheme determines the shielding tensor as the second derivative of the total electronic energy with respect to an external magnetic field and a nuclear magnetic moment. The induced current density due to a periodic perturbation from nuclear magnetic moments is obtained through numerical differentiation, whereas the influence of the responding perturbation in terms of the external magnetic field is evaluated analytically. The method is implemented into the periodic program BAND. It employs a Bloch basis set made up of Slater-type or numeric atomic orbitals and represents the Kohn-Sham potential fully without the use of effective core potentials. Results from calculations of NMR shielding constants based on the present approach are presented for isolated molecules as well as systems with one-, two- and three-dimensional periodicity. The reported values are compared to experiment and results from calculations on cluster models.

The PAW/GIPAW approach for computing NMR parameters: a new dimension added to NMR study of solids

In 2001, Mauri and Pickard introduced the gauge including projected augmented wave (GIPAW) method that enabled for the first time the calculation of all-electron NMR parameters in solids, i.e. accounting for periodic boundary conditions. The GIPAW method roots in the plane wave pseudopotential formalism of the density functional theory (DFT), and avoids the use of the cluster approximation. This method has undoubtedly revitalized the interest in quantum chemical calculations in the solid-state NMR community. It has quickly evolved and improved so that the calculation of the key components of NMR interactions, namely the shielding and electric field gradient tensors, has now become a routine for most of the common nuclei studied in NMR. Availability of reliable implementations in several software packages (CASTEP, Quantum Espresso, PARATEC) make its usage more and more increasingly popular, maybe indispensable in near future for all material NMR studies. The majority of nuclei of the periodic table have already been investigated by GIPAW, and because of its high accuracy it is quickly becoming an essential tool for interpreting and understanding experimental NMR spectra, providing reliable assignments of the observed resonances to crystallographic sites or enabling a priori prediction of NMR data. The continuous increase of computing power makes ever larger (and thus more realistic) systems amenable to first-principles analysis. In the near future perspectives, as the incorporation of dynamical effects and/or disorder are still at their early developments, these areas will certainly be the prime target.

Electrical polarization and orbital magnetization: the modern theories

Macroscopic polarization P and magnetization M are the most fundamental concepts in any phenomenological description of condensed media. They are intensive vector quantities that intuitively carry the meaning of dipole per unit volume. But for many years both P and the orbital term in M evaded even a precise microscopic definition, and severely challenged quantum-mechanical calculations. If one reasons in terms of a finite sample, the electric (magnetic) dipole is affected in an extensive way by charges (currents) at the sample boundary, due to the presence of the unbounded position operator in the dipole definitions. Therefore P and the orbital term in M--phenomenologically known as bulk properties--apparently behave as surface properties; only spin magnetization is problemless. The field has undergone a genuine revolution since the early 1990s. Contrary to a widespread incorrect belief, P has nothing to do with the periodic charge distribution of the polarized crystal: the former is essentially a property of the phase of the electronic wavefunction, while the latter is a property of its modulus. Analogously, the orbital term in M has nothing to do with the periodic current distribution in the magnetized crystal. The modern theory of polarization, based on a Berry phase, started in the early 1990s and is now implemented in most first-principle electronic structure codes. The analogous theory for orbital magnetization started in 2005 and is partly work in progress. In the electrical case, calculations have concerned various phenomena (ferroelectricity, piezoelectricity, and lattice dynamics) in several materials, and are in spectacular agreement with experiments; they have provided thorough understanding of the behaviour of ferroelectric and piezoelectric materials. In the magnetic case the very first calculations are appearing at the time of writing (2010). Here I review both theories on a uniform ground in a density functional theory (DFT) framework, pointing out analogies and differences. Both theories are deeply rooted in geometrical concepts, elucidated in this work. The main formulae for crystalline systems express P and M in terms of Brillouin-zone integrals, discretized for numerical implementation. I also provide the corresponding formulae for disordered systems in a single k-point supercell framework. In the case of P the single-point formula has been widely used in the Car-Parrinello community to evaluate IR spectra.

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

Quantum theory of orbital magnetization and its generalization to interacting systems

Based on standard perturbation theory, we present a full quantum derivation of the formula for the orbital magnetization in periodic systems. The derivation is generally valid for insulators with or without a Chern number, for metals at zero or finite temperatures, and at weak as well as strong magnetic fields. The formula is shown to be valid in the presence of electron-electron interaction, provided the one-electron energies and wave functions are calculated self-consistently within the framework of the exact current and spin-density functional theory.

NMR Solvent Shifts of Adenine in Aqueous Solution from Hybrid QM/MM Molecular Dynamics Simulations

We present first principles calculations of the NMR solvent shift of adenine in aqueous solution. The calculations are based on snapshots sampled from a molecular dynamics simulation, which were obtained via a hybrid quantum-mechanical/mechanical modeling approach, using an all-atom force field (TIP3P). We find that the solvation via the strongly fluctuating hydrogen bond network of water leads to nontrivial changes in the NMR spectra of the solutes regarding the ordering of the resonance lines. Although there are still sizable deviations from experiment, the overall agreement is satisfactory for the 1H and 15N NMR shifts. Our work is another step toward a realistic first-principles prediction of NMR chemical shifts in complex chemical environments.

Current Densities and Nucleus-Independent Chemical Shift Maps from Reciprocal-Space Density Functional Perturbation Theory Calculations

A method to calculate condensed-matter nucleus independent chemical shift maps (NICS maps) from first principles in the framework of density functional theory is presented. I use a pseudopotential plane-wave approach in which the electronic current density and the NICS map are obtained from an inverse Fourier transformation of the induced magnetic field represented in reciprocal space (G space). Due to its intrinsically periodic description, the method is suitable for isolated molecules (by using a supercell technique) and for condensed-phase systems like solids. The periodic NICS method was applied to hydrogen-bonded calixhydroquinone nanotubes, crystalline graphite, and two carbon nanotube systems.

Orbital magnetization in periodic insulators

Working in the Wannier representation, we derive an expression for the orbital magnetization of a periodic insulator. The magnetization is shown to be comprised of two contributions, an obvious one associated with the internal circulation of bulklike Wannier functions in the interior, and an unexpected one arising from net currents carried by Wannier functions near the surface. Each contribution can be expressed as a bulk property in terms of Bloch functions in a gauge-invariant way. Our expression is verified by comparing numerical tight-binding calculations for finite and periodic samples.

Orbital Magnetization in Extended Systems

While the orbital magnetic dipole moment of any finite sample is well-defined, it becomes ill-defined in the thermodynamic limit as a result of the unboundedness of the position operator. Effects due to surface currents and to bulk magnetization are not easily disentangled. The corresponding electrical problem, where surface charges and bulk polarization appear as entangled, was solved about a decade ago by the modern theory of polarization, based on a Berry phase. We follow a similar path here, making progress toward a bulk expression for the orbital magnetization in an insulator represented by a lattice-periodic Hamiltonian with broken time-reversal symmetry. We therefore limit ourselves to the case where the macroscopic (i.e. cell-averaged) magnetic field vanishes. We derive an expression for the contribution to the magnetization arising from the circulating currents internal to the bulk Wannier functions, and then transform to obtain a Brillouin zone integral involving the occupied Bloch orbitals. A version suitable for practical implementation in discretized reciprocal space is also derived, and the gauge invariance of both versions is explicitly shown. However, tests on a tight-binding model indicate the presence of additional edge currents, and it remains to be determined whether these can be related to the bulk band structure.

First-Principles Calculation of 17O and 25Mg NMR Shieldings in MgO at Finite Temperature: Rovibrational Effect in Solids

The temperature dependence of (17)O and (25)Mg NMR chemical shifts in solid MgO have been calculated using a first-principles approach. Density functional theory, pseudopotentials, a plane-wave basis set, and periodic boundary conditions were used both to describe the motion of the nuclei and to compute the NMR chemical shifts. The chemical shifts were obtained using the gauge-including projector augmented wave method. In a crystalline solid, the temperature dependence is due to both (i) the variation of the averaged equilibrium structure and (ii) the fluctuation of the atoms around this structure. In MgO, the equilibrium structure at each temperature is uniquely defined by the cubic lattice parameters, which we take from experiment. We evaluate the effect of the fluctuations within a quasiharmonic approximation. In particular, the dynamical matrix, defining the harmonic Hamiltonian, has been computed for each equilibrium volume. This harmonic Hamiltonian was used to generate nuclear configurations that obey quantum statistical mechanics. The chemical shifts were averaged over these nuclear configurations. The results reproduce the previously published experimental NMR data measured on MgO between room temperature and 1000 degrees C. It is shown that the chemical shift behavior with temperature cannot be explained by thermal expansion alone. Vibrational corrections due to the fluctuations of atoms around their equilibrium position are crucial to reproduce the experimental results.

Ab initio calculations in a uniform magnetic field using periodic supercells

We present a formulation of ab initio electronic structure calculations in a finite magnetic field, which retains the simplicity and efficiency of techniques widely used in first principles molecular dynamics simulations, based on plane-wave basis sets and Fourier transforms. In addition we discuss results obtained with this method for the energy spectrum of interacting electrons in quantum wells, and for the electronic properties of dense fluid deuterium in a uniform magnetic field.

Nonlocal pseudopotentials and magnetic fields

We show how to describe the coupling of electrons to nonuniform magnetic fields in the framework of the widely used norm-conserving pseudopotential approximation for electronic structure calculations. Our derivation applies to magnetic fields that are smooth on the scale of the core region. The method is validated by application to the calculation of the magnetic susceptibility of molecules within density functional theory (DFT) in the local density approximation. Our results are compared with high-quality all-electron DFT results obtained using Gaussian basis sets and another recently proposed pseudopotential formalism.

Coupling of nonlocal potentials to electromagnetic fields

Nonlocal Hamiltonians are used widely in first-principles quantum calculations; the nonlocality stems from eliminating undesired degrees of freedom, e.g., core electrons. To date, attempts to couple nonlocal systems to external electromagnetic (EM) fields have been heuristic or limited to weak or long wavelength fields. Using Feynman path integrals, we derive an exact, closed-form coupling of arbitrary EM fields to nonlocal systems. Our results justify and clarify the couplings used to date and are essential for systematic computation of linear and especially nonlinear responses.

Atomic structure of icosahedral B4C boron carbide from a first principles analysis of NMR spectra

Density functional theory is demonstrated to reproduce the 13C and 11B NMR chemical shifts of icosahedral boron carbides with sufficient accuracy to extract previously unresolved structural information from experimental NMR spectra. B4C can be viewed as an arrangement of 3-atom linear chains and 12-atom icosahedra. According to our results, all the chains have a CBC structure. Most of the icosahedra have a B11C structure with the C atom placed in a polar site, and a few percent have a B (12) structure or a B10C2 structure with the two C atoms placed in two antipodal polar sites.