The ONETEP linear-scaling density functional theory program.

The Journal of chemical physics 152:17 (2020) 174111-174111

Authors:

Joseph CA Prentice, Jolyon Aarons, James C Womack, Alice EA Allen, Lampros Andrinopoulos, Lucian Anton, Robert A Bell, Arihant Bhandari, Gabriel A Bramley, Robert J Charlton, Rebecca J Clements, Daniel J Cole, Gabriel Constantinescu, Fabiano Corsetti, Simon M-M Dubois, Kevin KB Duff, Jos茅 Mar铆a Escart铆n, Andrea Greco, Quintin Hill, Louis P Lee, Edward Linscott, David D O'Regan, Maximillian JS Phipps, Laura E Ratcliff, 脕lvaro Ruiz Serrano, Edward W Tait, Gilberto Teobaldi, Valerio Vitale, Nelson Yeung, Tim J Zuehlsdorff, Jacek Dziedzic, Peter D Haynes, Nicholas DM Hine, Arash A Mostofi, Mike C Payne, Chris-Kriton Skylaris

Abstract:

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.

Competing pairing interactions responsible for the large upper critical field in a stoichiometric iron-based superconductor CaKFe4As4

Physical Review B American Physical Society 101 (2020) 134502

Authors:

Matthew Bristow, William Knafo, Pascal Reiss, William Meier, Paul C Canfield, Stephen J Blundell, Amalia Coldea

Abstract:

The upper critical field of multiband superconductors is an important quantity that can reveal details about the nature of the superconducting pairing. Here we experimentally map out the complete upper-critical-field phase diagram of a stoichiometric superconductor, CaKFe4As4, up to 90 T for different orientations of the magnetic field and at temperatures down to 4.2K. The upper critical fields are extremely large, reaching values close to ∼3 Tc at the lowest temperature, and the anisotropy decreases dramatically with temperature, leading to essentially isotropic superconductivity at 4.2K. We find that the temperature dependence of the upper critical field can be well described by a two-band model in the clean limit with band-coupling parameters favoring intraband over interband interactions. The large Pauli paramagnetic effects together with the presence of the shallow bands is consistent with the stabilization of an FFLO state at low temperatures in this clean superconductor.

Anomalous high-magnetic field electronic state of the nematic superconductors FeSe1-xSx

Phys. Rev. Research 2, 013309 (2020) (2020)

Authors:

M Bristow, P Reiss, AA Haghighirad, Z Zajicek, SHIV Singh, T Wolf, D Graf, W Knafo, A McCollam, AMALIA Coldea

Abstract:

Understanding superconductivity requires detailed knowledge of the normal electronic state from which it emerges. A nematic electronic state that breaks the rotational symmetry of the lattice can potentially promote unique scattering relevant for superconductivity. Here, we investigate the normal transport of superconducting FeSe$_{1-x}$S$_x$ across a nematic phase transition using high magnetic fields up to 69 T to establish the temperature and field-dependencies. We find that the nematic state is an anomalous non-Fermi liquid, dominated by a linear resistivity at low temperatures that can transform into a Fermi liquid, depending on the composition $x$ and the impurity level. Near the nematic end point, we find an extended temperature regime with $T^{1.5}$ resistivity. The transverse magnetoresistance inside the nematic phase has as a $H^{1.55}$ dependence over a large magnetic field range and it displays an unusual peak at low temperatures inside the nematic phase. Our study reveals anomalous transport inside the nematic phase, driven by the subtle interplay between the changes in the electronic structure of a multi-band system and the unusual scattering processes affected by large magnetic fields and disorder

Anomalous high-magnetic field electronic state of the nematic superconductors FeSe 1 鈭 x S x

Physical Review Research American Physical Society 2:1 (2020) 013309

Authors:

M Bristow, P Reiss, AA Haghighirad, Z Zajicek, SHIV Singh, T Wolf, D Graf, W Knafo, A McCollam, Amalia Coldea

Abstract:

Understanding superconductivity requires detailed knowledge of the normal electronic state from which it emerges. A nematic electronic state that breaks the rotational symmetry of the lattice can potentially promote unique scattering relevant for superconductivity. Here, we investigate the normal transport of superconducting FeSe 1 鈭 x S x across a nematic phase transition using high-magnetic fields up to 69 T to establish the temperature and field dependencies. We find that the nematic state is dominated by a linear resistivity at low temperatures that evolves towards Fermi-liquid behavior, depending on the composition x and the impurity level. Near the nematic end point, we find an extended temperature regime with 鈭 T 1.5 resistivity, different from the behavior found near an antiferromagnetic critical point. The variation of the resistivity exponent with temperature reflects the importance of the nematoelastic coupling that can also suppress divergent critical fluctuations at the nematic end point. The transverse magnetoresistance inside the nematic phase has a 鈭 H 1.55 dependence over a large magnetic field range and it displays an unusual peak at low temperatures inside the nematic phase. Our study reveals anomalous transport inside the nematic phase, influenced by both changes in the electronic structure and the scattering with the lattice and spin fluctuations.

Competing pairing interactions responsible for the large upper critical field in a stoichiometric iron-based superconductor, CaKFe$_4$As$_4$

(2020)

Authors:

Matthew Bristow, William Knafo, Pascal Reiss, William Meier, Paul C Canfield, Stephen J Blundell, Amalia I Coldea