No Arabic abstract
CeCuAl3 and CeAuAl3, crystallizing in the non-centrosymmetric BaNiSn3 tetragonal structure, are known mainly for their unusual neutron scattering spectra involving additional excitations ascribed to vibron quasi-bound quantum state in CeCuAl3 and anti-crossing of phonon and crystal field excitations in CeAuAl3. In this work, we present results of nuclear magnetic resonance experiments on their lanthanum analogues - LaCuAl3 and LaAuAl3. The character of nuclear magnetic resonance spectra of 139La, 27Al, and 65Cu measured in LaAuAl3 and LaCuAl3 is dominated by electric quadrupole interaction. The spectral parameters acquired from experimental data are confronted with values obtained from the electronic structure calculations. The results show remarkable diffrences for the two compounds. The 139La spectrum in LaAuAl3 can be interpreted by a single spectral component corresponding to uniform environment of La atoms in the crystal structure, whereas for LaCuAl3 the spectrum decomposition yields a wide distribution of spectral parameters, which is not possible to explain by a single La environment, and multiple non-equivalent La positions in the crystal structure are required to interpret the spectrum.
A magnetic skyrmion observed experimentally in chiral magnets is a topologically protected spin texture. For their unique properties, such as high mobility under current drive, skyrmions have huge potential for applications in next-generation spintronic devices. Defects naturally occurring in magnets have profound effects on the static and dynamical properties of skyrmions. In this work, we study the effect of an atomic defect on a skyrmion using the first-principles calculations within the density functional theory, taking MnSi as an example. By substituting one site of Mn or Si with different elements, we can tune the pinning energy. The effects of pinning by an atomic defect can be understood qualitatively within a phenomenological model.
Localized basis sets in the projector augmented wave formalism allow for computationally efficient calculations within density functional theory (DFT). However, achieving high numerical accuracy requires an extensive basis set, which also poses a fundamental problem for the interpretation of the results. We present a way to obtain a reduced basis set of atomic orbitals through the subdiagonalization of each atomic block of the Hamiltonian. The resulting local orbitals (LOs) inherit the information of the local crystal field. In the LO basis, it becomes apparent that the Hamiltonian is nearly block-diagonal, and we demonstrate that it is possible to keep only a subset of relevant LOs which provide an accurate description of the physics around the Fermi level. This reduces to some extent the redundancy of the original basis set, and at the same time it allows one to perform post-processing of DFT calculations, ranging from the interpretation of electron transport to extracting effective tight-binding Hamiltonians, very efficiently and without sacrificing the accuracy of the results.
We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections.
Linear scaling methods for density-functional theory (DFT) simulations are formulated in terms of localised orbitals in real-space, rather than the delocalised eigenstates of conventional approaches. In local-orbital methods, relative to conventional DFT, desirable properties can be lost to some extent, such as the translational invariance of the total energy of a system with respect to small displacements and the smoothness of the potential energy surface. This has repercussions for calculating accurate ionic forces and geometries. In this work we present results from textsc{onetep}, our linear scaling method based on localised orbitals in real-space. The use of psinc functions for the underlying basis set and on-the-fly optimisation of the localised orbitals results in smooth potential energy surfaces that are consistent with ionic forces calculated using the Hellmann-Feynman theorem. This enables accurate geometry optimisation to be performed. Results for surface reconstructions in silicon are presented, along with three example systems demonstrating the performance of a quasi-Newton geometry optimisation algorithm: an organic zwitterion, a point defect in an ionic crystal, and a semiconductor nanostructure.
We propose a method to decompose the total energy of a supercell containing defects into contributions of individual atoms, using the energy density formalism within density functional theory. The spatial energy density is unique up to a gauge transformation, and we show that unique atomic energies can be calculated by integrating over Bader and charge-neutral volumes for each atom. Numerically, we implement the energy density method in the framework of the Vienna ab initio simulation package (VASP) for both norm-conserving and ultrasoft pseudopotentials and the projector augmented wave method, and use a weighted integration algorithm to integrate the volumes. The surface energies and point defect energies can be calculated by integrating the energy density over the surface region and the defect region, respectively. We compute energies for several surfaces and defects: the (110) surface energy of GaAs, the mono-vacancy formation energies of Si, the (100) surface energy of Au, and the interstitial formation energy of O in the hexagonal close-packed Ti crystal. The surface and defect energies calculated using our method agree with size-converged calculations of the difference between the total energies of the system with and without the defect. Moreover, the convergence of the defect energies with size can be found from a single calculation.