-
Layer group classification of two-dimensional materials
Authors:
Jingheng Fu,
Mikael Kuisma,
Ask Hjorth Larsen,
Kohei Shinohara,
Atsushi Togo,
Kristian S. Thygesen
Abstract:
The symmetry of a crystal structure with a three-dimensional (3D) lattice can be classified by one of the 230 space group types. For some types of crystals, e.g. crystalline films, surfaces, or planar interfaces, it is more appropriate to assume a 2D lattice. With this assumption, the structure can be classified by one of the 80 layer group types. We have implemented an algorithm to determine the…
▽ More
The symmetry of a crystal structure with a three-dimensional (3D) lattice can be classified by one of the 230 space group types. For some types of crystals, e.g. crystalline films, surfaces, or planar interfaces, it is more appropriate to assume a 2D lattice. With this assumption, the structure can be classified by one of the 80 layer group types. We have implemented an algorithm to determine the layer group type of a 3D structure with a 2D lattice, and applied it to more than 15,000 monolayer structures in the Computational 2D Materials (C2DB) database. We compare the classification of monolayers by layer groups and space groups, respectively. The latter is defined as the space group of the 3D bulk structure obtained by repeating the monolayer periodically in the direction perpendicular to the 2D lattice (AA-stacking). By this correspondence, nine pairs of layer group types are mapped to the same space group type due to the inability of the space group to distinguish the in-plane and out-of-plane axes. In total 18% of the monolayers in the C2DB belong to one of these layer group pairs and are thus not properly classified by the space group type.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
Accurate Dielectric Response of Solids: Combining the Bethe-Salpeter Equation with the Random Phase Approximation
Authors:
Amalie H. Søndersted,
Mikael Kuisma,
Jakob K. Svaneborg,
Mark Kamper Svendsen,
Kristian S. Thygesen
Abstract:
The Bethe-Salpeter equation (BSE) can provide an accurate description of low-energy optical spectra of insulating crystals - even when excitonic effects are important. However, due to high computational costs it is only possible to include a few bands in the BSE Hamiltonian. As a consequence, the dielectric screening given by the real part of the dielectric function can be significantly underestim…
▽ More
The Bethe-Salpeter equation (BSE) can provide an accurate description of low-energy optical spectra of insulating crystals - even when excitonic effects are important. However, due to high computational costs it is only possible to include a few bands in the BSE Hamiltonian. As a consequence, the dielectric screening given by the real part of the dielectric function can be significantly underestimated by the BSE. Here we show that universally accurate optical response functions can be obtained by combining a four-point BSE-like equation for the irreducible polarisability with a two-point Dyson equation which includes the higher-lying transitions within the random phase approximation (RPA). The new method is referred to as BSE+. It has a computational cost comparable to the BSE but a much faster convergence with respect to the size of the electron-hole basis. We use the method to calculate refractive indices and electron energy loss spectra for a test set of semiconductors and insulators. In all cases the BSE+ yields excellent agreement with experimental data across a wide frequency range and outperforms both BSE and RPA.
△ Less
Submitted 30 November, 2023;
originally announced November 2023.
-
GPAW: An open Python package for electronic-structure calculations
Authors:
Jens Jørgen Mortensen,
Ask Hjorth Larsen,
Mikael Kuisma,
Aleksei V. Ivanov,
Alireza Taghizadeh,
Andrew Peterson,
Anubhab Haldar,
Asmus Ougaard Dohn,
Christian Schäfer,
Elvar Örn Jónsson,
Eric D. Hermes,
Fredrik Andreas Nilsson,
Georg Kastlunger,
Gianluca Levi,
Hannes Jónsson,
Hannu Häkkinen,
Jakub Fojt,
Jiban Kangsabanik,
Joachim Sødequist,
Jouko Lehtomäki,
Julian Heske,
Jussi Enkovaara,
Kirsten Trøstrup Winther,
Marcin Dulak,
Marko M. Melander
, et al. (22 additional authors not shown)
Abstract:
We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually indepen…
▽ More
We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE) providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation (BSE), variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support of GPU acceleration has been achieved with minor modifications of the GPAW code thanks to the CuPy library. We end the review with an outlook describing some future plans for GPAW.
△ Less
Submitted 16 April, 2024; v1 submitted 23 October, 2023;
originally announced October 2023.
-
Excitonic insulators and superfluidity in 2D bilayers without external fields
Authors:
Fredrik Nilsson,
Mikael Kuisma,
Sahar Pakdel,
Kristian S. Thygesen
Abstract:
We explore a new platform for realising excitonic insulators, namely van der Waals (vdW) bilayers comprising 2D Janus materials. In previous studies, Type-II heterobilayers have been brought to the excitonic insulating regime by tuning the band alignment using external gates. In contrast, the Janus bilayers of the current work represent intrinsic excitonic insulators. We first conduct ab initio ca…
▽ More
We explore a new platform for realising excitonic insulators, namely van der Waals (vdW) bilayers comprising 2D Janus materials. In previous studies, Type-II heterobilayers have been brought to the excitonic insulating regime by tuning the band alignment using external gates. In contrast, the Janus bilayers of the current work represent intrinsic excitonic insulators. We first conduct ab initio calculations to obtain the quasiparticle band structures, screened Coulomb interaction, and interlayer exciton binding energies of the bilayers. These ab initio-derived quantities are then used to construct a BCS-like Hamiltonian of the exciton condensate. By solving the mean-field gap equation, we identify 16 vdW Janus bilayers with insulating ground states and superfluid properties. Our calculations expose a new class of advanced materials that are likely to exhibit novel excitonic phases at low temperatures, and highlight the subtle competition between interlayer hybridization, spin-orbit coupling, and dielectric screening that govern their properties.
△ Less
Submitted 9 August, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Dipolar coupling of nanoparticle-molecule assemblies: An efficient approach for studying strong coupling
Authors:
Jakub Fojt,
Tuomas P. Rossi,
Tomasz J. Antosiewicz,
Mikael Kuisma,
Paul Erhart
Abstract:
Strong light-matter interactions facilitate not only emerging applications in quantum and non-linear optics but also modifications of materials properties. In particular the latter possibility has spurred the development of advanced theoretical techniques that can accurately capture both quantum optical and quantum chemical degrees of freedom. These methods are, however, computationally very deman…
▽ More
Strong light-matter interactions facilitate not only emerging applications in quantum and non-linear optics but also modifications of materials properties. In particular the latter possibility has spurred the development of advanced theoretical techniques that can accurately capture both quantum optical and quantum chemical degrees of freedom. These methods are, however, computationally very demanding, which limits their application range. Here, we demonstrate that the optical spectra of nanoparticle-molecule assemblies, including strong coupling effects, can be predicted with good accuracy using a subsystem approach, in which the response functions of the different units are coupled only at the dipolar level. We demonstrate this approach by comparison with previous time-dependent density functional theory calculations for fully coupled systems of Al nanoparticles and benzene molecules. While the present study only considers few-particle systems, the approach can be readily extended to much larger systems and to include explicit optical-cavity modes.
△ Less
Submitted 13 January, 2021;
originally announced January 2021.
-
Real-Time Time-Dependent Density Functional Theory Implementation of Electronic Circular Dichroism Applied to Nanoscale Metal-Organic Clusters
Authors:
Esko Makkonen,
Tuomas P. Rossi,
Ask Hjorth Larsen,
Olga Lopez-Acevedo,
Patrick Rinke,
Mikael Kuisma,
Xi Chen
Abstract:
Electronic circular dichroism (ECD) is a powerful spectroscopical method for investigating chiral properties at the molecular level. ECD calculations with the commonly used linear-response time-dependent density functional theory (LR-TDDFT) framework can be prohibitively costly for large systems. To alleviate this problem, we present here an ECD implementation for the projector augmented-wave meth…
▽ More
Electronic circular dichroism (ECD) is a powerful spectroscopical method for investigating chiral properties at the molecular level. ECD calculations with the commonly used linear-response time-dependent density functional theory (LR-TDDFT) framework can be prohibitively costly for large systems. To alleviate this problem, we present here an ECD implementation for the projector augmented-wave method in the real-time-propagation TDDFT (RT-TDDFT) framework in the open-source GPAW code. Our implementation supports both local atomic basis set and real-space finite-difference representations of wave functions. We benchmark our implementation against an existing LR-TDDFT implementation in GPAW for small chiral molecules. We then demonstrate the efficiency of our local atomic basis set implementation for a large hybrid nanocluster.
△ Less
Submitted 16 July, 2020;
originally announced July 2020.
-
Hot-Carrier Generation in Plasmonic Nanoparticles: The Importance of Atomic Structure
Authors:
Tuomas P. Rossi,
Paul Erhart,
Mikael Kuisma
Abstract:
Metal nanoparticles are attractive for plasmon-enhanced generation of hot carriers, which may be harnessed in photochemical reactions. In this work, we analyze the coherent femtosecond dynamics of photon absorption, plasmon formation, and subsequent hot-carrier generation through plasmon dephasing using first-principles simulations. We predict the energetic and spatial hot-carrier distributions in…
▽ More
Metal nanoparticles are attractive for plasmon-enhanced generation of hot carriers, which may be harnessed in photochemical reactions. In this work, we analyze the coherent femtosecond dynamics of photon absorption, plasmon formation, and subsequent hot-carrier generation through plasmon dephasing using first-principles simulations. We predict the energetic and spatial hot-carrier distributions in small metal nanoparticles and show that the distribution of hot electrons is very sensitive to the local structure. Our results show that surface sites exhibit enhanced hot-electron generation in comparison to the bulk of the nanoparticle. While the details of the distribution depend on particle size and shape, as a general trend lower-coordinated surface sites such as corners, edges, and {100} facets exhibit a higher proportion of hot electrons than higher-coordinated surface sites such as {111} facets or the core sites. The present results thereby demonstrate how hot carriers could be tailored by careful design of atomic-scale structures in nanoscale systems.
△ Less
Submitted 29 July, 2020; v1 submitted 27 February, 2020;
originally announced February 2020.
-
Space partitioning of exchange-correlation functionals with the projector augmented-wave method
Authors:
H. Levämäki,
M. Kuisma,
K. Kokko
Abstract:
We implement a Becke fuzzy cells type space partitioning scheme for the purposes of exchange-correlation within the GPAW projector augmented-wave method based density functional theory code. Space partitioning is needed in the situation where one needs to treat different parts of a combined system with different exchange-correlation functionals. For example, bulk and surface regions of a system co…
▽ More
We implement a Becke fuzzy cells type space partitioning scheme for the purposes of exchange-correlation within the GPAW projector augmented-wave method based density functional theory code. Space partitioning is needed in the situation where one needs to treat different parts of a combined system with different exchange-correlation functionals. For example, bulk and surface regions of a system could be treated with functionals that are specifically designed to capture the distinct physics of those regions. Here we use the space partitioning scheme to implement the quasi-nonuniform exchange-correlation scheme, which is a useful practical approach for calculating metallic alloys on the generalized gradient approximation level. We also confirm the correctness of our implementation with a set of test calculations.
△ Less
Submitted 15 October, 2018;
originally announced October 2018.
-
Plasmon excitations in mixed metallic nanoarrays
Authors:
Kevin Conley,
Neha Nayyar,
Tuomas P. Rossi,
Mikael Kuisma,
Volodymyr Turkowski,
Martti Puska,
Talat S. Rahman
Abstract:
We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different ato…
▽ More
We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different atomic species and the hybridization between d and sp electrons. As a result it is possible to tune the position of the surface plasmon resonance, split it to several peaks, and eventually achieve broadband absorption of radiation. Mixing the arrays with transition metals can strongly attenuate the plasmonic behaviour. We analyze the origin of these phenomena and show that they arise from rich interactions between single-particle electron-hole and collective electron excitations. The tunability of the plasmonic response of arrays of atomic chains, which can be realized on solid surfaces, opens wide possibilities for their applications. In the present study we obtain guidelines how the desired properties can be achieved.
△ Less
Submitted 30 September, 2018;
originally announced October 2018.
-
Electric field-controlled reversible order-disorder switching of a metal tip surface
Authors:
Ludvig de Knoop,
Mikael Juhani Kuisma,
Joakim Löfgren,
Kristof Lodewijks,
Mattias Thuvander,
Paul Erhart,
Alexandre Dmitriev,
Eva Olsson
Abstract:
While it is well established that elevated temperatures can induce surface roughening of metal surfaces, the effect of a high electric field on the atomic structure at ambient temperature has not been investigated in detail. Here we show with atomic resolution using in situ transmission electron microscopy how intense electric fields induce reversible switching between perfect crystalline and diso…
▽ More
While it is well established that elevated temperatures can induce surface roughening of metal surfaces, the effect of a high electric field on the atomic structure at ambient temperature has not been investigated in detail. Here we show with atomic resolution using in situ transmission electron microscopy how intense electric fields induce reversible switching between perfect crystalline and disordered phases of gold surfaces at room temperature. Ab initio molecular dynamics simulations reveal that the mechanism behind the structural change can be attributed to a vanishing energy cost in forming surface defects in high electric fields. Our results demonstrate how surface processes can be directly controlled at the atomic scale by an externally applied electric field, which promotes an effective decoupling of the topmost surface layers from the underlying bulk. This opens up opportunities for development of active nanodevices in e.g. nanophotonics and field-effect transistor technology as well as fundamental research in materials characterization and of yet unexplored dynamically-controlled low-dimensional phases of matter.
△ Less
Submitted 7 May, 2018;
originally announced May 2018.
-
Interlayer exciton dynamics in van der Waals heterostructures
Authors:
Simon Ovesen,
Samuel Brem,
Christopher Linderälv,
Mikael Kuisma,
Paul Erhart,
Malte Selig,
Ermin Malic
Abstract:
Exciton binding energies of hundreds of meV and strong light absorption in the optical frequency range make transition metal dichalcogenides (TMDs) promising for novel optoelectronic nanodevices. In particular, atomically thin TMDs can be stacked to heterostructures enabling the design of new materials with tailored properties. The strong Coulomb interaction gives rise to interlayer excitons, wher…
▽ More
Exciton binding energies of hundreds of meV and strong light absorption in the optical frequency range make transition metal dichalcogenides (TMDs) promising for novel optoelectronic nanodevices. In particular, atomically thin TMDs can be stacked to heterostructures enabling the design of new materials with tailored properties. The strong Coulomb interaction gives rise to interlayer excitons, where electrons and holes are spatially separated in different layers. In this work, we reveal the microscopic processes behind the formation, thermalization and decay of these fundamentally interesting and technologically relevant interlayer excitonic states. In particular, we present for the exemplary MoSe$_2$-WSe$_2$ heterostructure the interlayer exciton binding energies and wave functions as well as their time- and energy-resolved dynamics. Finally, we predict the dominant contribution of interlayer excitons to the photoluminescence of these materials.
△ Less
Submitted 23 April, 2018;
originally announced April 2018.
-
libvdwxc: A library for exchange-correlation functionals in the vdW-DF family
Authors:
Ask Hjorth Larsen,
Mikael Kuisma,
Joakim Löfgren,
Yann Pouillon,
Paul Erhart,
Per Hyldgaard
Abstract:
We present libvdwxc, a general library for evaluating the energy and potential for the family of vdW-DF exchange--correlation functionals. libvdwxc provides an efficient implementation of the vdW-DF method and can be interfaced with various general-purpose DFT codes. Currently, the GPAW and Octopus codes implement interfaces to libvdwxc. The present implementation emphasizes scalability and parall…
▽ More
We present libvdwxc, a general library for evaluating the energy and potential for the family of vdW-DF exchange--correlation functionals. libvdwxc provides an efficient implementation of the vdW-DF method and can be interfaced with various general-purpose DFT codes. Currently, the GPAW and Octopus codes implement interfaces to libvdwxc. The present implementation emphasizes scalability and parallel performance, and thereby enables \textit{ab initio} calculations of nanometer-scale complexes. The numerical accuracy is benchmarked on the S22 test set whereas parallel performance is benchmarked on ligand-protected gold nanoparticles ($\text{Au}_{144}(\text{SC}_{11}\text{NH}_{25})_{60}$) up to 9696 atoms.
△ Less
Submitted 7 June, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
Kohn-Sham decomposition in real-time time-dependent density-functional theory: An efficient tool for analyzing plasmonic excitations
Authors:
Tuomas P. Rossi,
Mikael Kuisma,
Martti J. Puska,
Risto M. Nieminen,
Paul Erhart
Abstract:
The real-time-propagation formulation of time-dependent density-functional theory (RT-TDDFT) is an efficient method for modeling the optical response of molecules and nanoparticles. Compared to the widely adopted linear-response TDDFT approaches based on, e.g., the Casida equations, RT-TDDFT appears, however, lacking efficient analysis methods. This applies in particular to a decomposition of the…
▽ More
The real-time-propagation formulation of time-dependent density-functional theory (RT-TDDFT) is an efficient method for modeling the optical response of molecules and nanoparticles. Compared to the widely adopted linear-response TDDFT approaches based on, e.g., the Casida equations, RT-TDDFT appears, however, lacking efficient analysis methods. This applies in particular to a decomposition of the response in the basis of the underlying single-electron states. In this work, we overcome this limitation by developing an analysis method for obtaining the Kohn-Sham electron-hole decomposition in RT-TDDFT. We demonstrate the equivalence between the developed method and the Casida approach by a benchmark on small benzene derivatives. Then, we use the method for analyzing the plasmonic response of icosahedral silver nanoparticles up to Ag$_{561}$. Based on the analysis, we conclude that in small nanoparticles individual single-electron transitions can split the plasmon into multiple resonances due to strong single-electron-plasmon coupling whereas in larger nanoparticles a distinct plasmon resonance is formed.
△ Less
Submitted 8 March, 2017;
originally announced March 2017.
-
Localized surface plasmon resonance in silver nanoparticles: Atomistic first-principles time-dependent density-functional theory calculations
Authors:
Mikael Kuisma,
Arto Sakko,
Tuomas P. Rossi,
Ask H. Larsen,
Jussi Enkovaara,
Lauri Lehtovaara,
Tapio T. Rantala
Abstract:
We observe using ab initio methods that localized surface plasmon resonances in icosahedral silver nanoparticles enter the asymptotic region already between diameters of 1 and 2 nm, converging close to the classical quasistatic limit around 3.4 eV. We base the observation on time-dependent density-functional theory simulations of the icosahedral silver clusters Ag$_{55}$ (1.06 nm), Ag$_{147}$ (1.6…
▽ More
We observe using ab initio methods that localized surface plasmon resonances in icosahedral silver nanoparticles enter the asymptotic region already between diameters of 1 and 2 nm, converging close to the classical quasistatic limit around 3.4 eV. We base the observation on time-dependent density-functional theory simulations of the icosahedral silver clusters Ag$_{55}$ (1.06 nm), Ag$_{147}$ (1.60 nm), Ag$_{309}$ (2.14 nm), and Ag$_{561}$ (2.68 nm). The simulation method combines the adiabatic GLLB-SC exchange-correlation functional with real time propagation in an atomic orbital basis set using the projector-augmented wave method. The method has been implemented for the electron structure code GPAW within the scope of this work. We obtain good agreement with experimental data and modeled results, including photoemission and plasmon resonance. Moreover, we can extrapolate the ab initio results to the classical quasistatically modeled icosahedral clusters.
△ Less
Submitted 8 April, 2015; v1 submitted 24 March, 2015;
originally announced March 2015.
-
Kohn-Sham potential with discontinuity for band gap materials
Authors:
M. Kuisma,
J. Ojanen,
J. Enkovaara,
T. T. Rantala
Abstract:
We model a Kohn-Sham potential with a discontinuity at integer particle numbers derived from the GLLB approximation of Gritsenko et al. We evaluate the Kohn-Sham gap and the discontinuity to obtain the quasiparticle gap. This allows us to compare the Kohn-Sham gaps to those obtained by accurate many-body perturbation theory based optimized potential methods. In addition, the resulting quasiparticl…
▽ More
We model a Kohn-Sham potential with a discontinuity at integer particle numbers derived from the GLLB approximation of Gritsenko et al. We evaluate the Kohn-Sham gap and the discontinuity to obtain the quasiparticle gap. This allows us to compare the Kohn-Sham gaps to those obtained by accurate many-body perturbation theory based optimized potential methods. In addition, the resulting quasiparticle band gap is compared to experimental gaps. In the GLLB model potential, the exchange-correlation hole is modeled using a GGA energy density and the response of the hole to density variations is evaluated by using the common-denominator approximation and homogeneous electron gas based assumptions. In our modification, we have chosen the PBEsol potential as the GGA to model the exchange hole, and add a consistent correlation potential. The method is implemented in the GPAW code, which allows efficient parallelization to study large systems. A fair agreement for Kohn-Sham and the quasiparticle band gaps with semiconductors and other band gap materials is obtained with a potential which is as fast as GGA to calculate.
△ Less
Submitted 31 May, 2010; v1 submitted 1 March, 2010;
originally announced March 2010.