Researcher profile

Hannes Jónsson

Hannes Jónsson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

Indirect Mechanism of Au adatom Diffusion on the Si(100) Surface

Calculations of the diffusion of a Au adatom on the dimer reconstructed Si(100)-2x1 surface reveal an interesting mechanism that differs significantly from a direct path between optimal binding sites, which are located in between dimer rows. Instead, the active diffusion mechanism involves promotion of the adatom to higher energy sites on top of a dimer row and then fast migration along the row, visiting ca. a hundred sites at room temperature, before falling back down into an optimal binding site. This top-of-row mechanism becomes more important the lower the temperature is. The calculations are carried out by finding minimum energy paths on the energy surface obtained from density functional theory within the PBEsol functional approximation followed by kinetic Monte Carlo simulations of the diffusion over a range of temperature from 200 K to 900 K. While the activation energy for the direct diffusion mechanism is calculated to be 0.84 eV, the effective activation energy for the indirect mechanism is on average 0.56 eV.

preprint2022arXiv

Reassignment of magic numbers for icosahedral Au clusters: 310, 564, 928 and 1426

Icosahedral Au clusters with three and four shells of atoms are found to deviate significantly from the commonly assumed Mackay structures. By introducing additional atoms in the surface shell and creating a vacancy in the center of the cluster, the calculated energy per atom can be lowered significantly, according to several different descriptions of the interatomic interaction. Analogous icosahedral structures with five and six shells of atoms are generated using the same structural motifs and are similarly found to be more stable than Mackay icosahedra. The lowest energy per atom is obtained with clusters containing 310, 564, 928 and 1426 atoms, as compared with the commonly assumed magic numbers of 309, 561, 923 and 1415. Some of the vertices in the optimized clusters have a hexagonal ring of atoms, rather than a pentagon, with the vertex atom missing. An inner shell atom in some cases moves outwards by more than an Ångström into the surface shell at the vertex site. This feature, as well as the wide distribution of nearest-neighbor distances in the surface layer, can strongly influence the catalytic properties of icosahedral clusters. The structural optimization is initially carried out using the GOUST method with atomic forces estimated with the EMT empirical potential function, but the atomic coordinates are then refined by minimization using electron density functional theory (DFT) or Gaussian approximation potential (GAP). A single energy barrier is found to separate the Mackay icosahedron from a lower energy structure where a string of atoms moves outwards in a concerted manner from the center so as to create a central vacancy while placing an additional atom in the surface shell.

preprint2022arXiv

Simulations of the Electrochemical Oxidation of Pt Nanoparticles of Various Shapes

The activity and stability of a platinum nanoparticle (NP) is not only affected by its size but additionally depends on its shape. To this end, simulations can identify structure-property relationships to make a priori decisions on the most promising structures. While activity is routinely probed by electronic structure calculations on simplified surface models, modeling the stability of NP model systems in electrochemical reactions is challenging due to the long timescale of relevant processes such as oxidation beyond the point of reversibility. In this work, a routine for simulating electrocatalyst stability is presented. The procedure is referred to as GREG after its main ingredients - a grand-canonical simulation approach using reactive force fields to model electrochemical reactions as a function of the galvanic cell potential. The GREG routine is applied to study the oxidation of 3 nm octahedral, cubic, dodecahedral, cuboctahedral, spherical, and tetrahexahedral platinum NPs. The oxidation process is analyzed using adsorption isobars as well as interaction energy heat maps that provide the basis for constructing electrochemical phase diagrams. Onset potentials for surface oxidation increase in the sequence cube ~= dodecahedron <= octahedron <= tetrahexahdron < sphere < cuboctahedron, establishing a relationship between oxidation behavior and surface facet structure. The electrochemical results are rationalized using structural and electronic analysis.

preprint2022arXiv

Transferable Potential Function for Flexible H$_2$O Molecules Based on the Single Center Multipole Expansion

A potential function is presented for describing a system of flexible H$_2$O molecules based on the single center multipole expansion (SCME) of the electrostatic interaction. The model, referred to as SCME/f, includes the variation of the molecular quadrupole moment as well as the dipole moment with changes in bond length and angle so as to reproduce results of high level electronic structure calculations. The multipole expansion also includes fixed octupole and hexadecapole moments, as well as anisotropic dipole-dipole, dipole-quadrupole and quadrupole-quadrupole polarizability tensors. The model contains five adjustable parameters related to the repulsive interaction and damping functions in the electrostatic and dispersion interactions. Their values are adjusted to reproduce the lowest energy isomers of small clusters, (H$_2$O)$_n$ with $n=2-6$, as well as measured properties of the ice Ih crystal. Subsequent calculations of the energy difference between the various isomer configurations of the clusters show that SCME/f gives good agreement with results of electronic structure calculations and represents a significant improvement over the previously presented rigid SCME potential function. Analysis of the vibrational frequencies of the clusters and structural properties of ice Ih crystal show the importance of accurately describing the variation of the quadrupole moment with molecular structure.

preprint2021arXiv

Assessment of the Accuracy of Density Functionals for Calculating Oxygen Reduction Reaction on Nitrogen Doped Graphene

Experimental studies of the oxygen reduction reaction (ORR) at nitrogen doped graphene electrodes have reported a remarkably low overpotential, on the order of 0.5 V, similar to Pt based electrodes. Theoretical calculations using density functional theory have lent support for this claim. However, other measurements have indicated that transition metal impurities are actually responsible for the ORR activity, thereby raising questions about the reliability of both the experiments and the calculations. In order to assess the accuracy of the theoretical calculations, various generalized gradient approximation (GGA), meta-GGA and hybrid functionals are employed here and calibrated against high-level wave function based coupled cluster calculations (CCSD(T)) of the overpotential as well as self-interaction corrected density functional calculations and published quantum Monte Carlo calculations of O adatom binding to graphene. The PBE0 and HSE06 hybrid functionals are found to give more accurate results than the GGA and meta-GGA functionals, as would be expected, and for low dopant concentration, 3.1%, the overpotential is calculated to be 1.0 V. The GGA and meta-GGA functionals give a lower estimate by as much as 0.4 V. When the dopant concentration is doubled, the overpotential calculated with hybrid functionals drops, while it increases in GGA functional calculations. The opposite trends result from different potential determining steps, the *OOH species being of central importance in the hybrid functional calculations while the reduction of *O determines the overpotential obtained in GGA and meta-GGA calculations. The results presented here are mainly based on calculations of periodic representations of the system, but a comparison is also made with molecular flake models which are found to give erratic results.

preprint2021arXiv

Elastic Collision Based Dynamic Partitioning Scheme for Hybrid Simulations

The scattering-adapted flexible inner region ensemble separator (SAFIRES) is a partitioning scheme designed to divide a simulation cell into two regions to be treated with different computational methodologies. SAFIRES prevents particles from crossing between regions and resolves boundary events through elastic collisions of the particles mediated by the boundary, conserving energy and momenta. A multiple-time-step propagation algorithm is introduced where the time step is scaled automatically to identify the moment a collision occurs. If the length of the time step is kept constant, the new propagator reduces to a regular algorithm for Langevin dynamics, and to the velocity Verlet algorithm for classical dynamics if the friction coefficient is set to zero. SAFIRES constitutes the exact limit of the premise behind boundary-based methods such as FIRES, BEST, and BCC which take advantage of the indistinguishability of molecules on opposite sides of the separator. It gives correct average ensemble statistics despite the introduction of an ensemble separator. SAFIRES is tested in simulations where the molecules on the two sides are treated in the same way, for a Lennard-Jones (LJ) liquid and a LJ liquid in contact with a surface, as well as for liquid modelling simulations using the TIP4P force field. Simulations using SAFIRES are shown to reproduce the unconstrained reference simulations without significant deviations.

preprint2020arXiv

Magnetic Skyrmion Annihilation by Quantum Mechanical Tunneling

Magnetic skyrmions are nano-scale magnetic states that could be used in various spintronics devices. A central issue is the mechanism and rate of various possible annihilation processes and the lifetime of metastable skyrmions. While most studies have focused on classical over-the-barrier mechanism for annihilation, it is also possible that quantum mechanical tunneling through the energy barrier takes place. Calculations of the lifetime of magnetic skyrmions in a two-dimensional lattice are presented and the rate of tunneling compared with the classical annihilation rate. A remarkably strong variation in the crossover temperature and the lifetime of the skyrmion is found as a function of the values of parameters in the extended Heisenberg Hamiltonian, i.e. the out-of-plane anisotropy, Dzyaloshinskii-Moriya interaction (DMI) and applied magnetic field. Materials parameters and conditions are identified where the onset of tunneling could be observed on a laboratory time scale. In particular, it is predicted that skyrmion tunneling could be observed in the PdFe/Ir(111) system when an external magnetic field on the order of 6 T is applied.

preprint2020arXiv

Skyrmions in antiferromagnets: thermal stability and the effect of external field and impurities

Calculations of skyrmions in antiferromagnets (AFMs) are presented, and their properties compared with skyrmions in corresponding ferromagnets (FMs). The rates of skyrmion collapse and escape through the boundary of a track, as well as the binding to and collapse at a non-magnetic impurity, are calculated as a function of applied magnetic field. The activation energy for skyrmion annihilation is the same in AFMs and corresponding FMs in the absence of an applied magnetic field. The pre-exponential factor in the Arrhenius rate law is, however, different because skyrmion dynamics is different in the two systems. An applied magnetic field has opposite effects on skyrmions in the two types of materials. In AFMs the rate of collapse of skyrmions as well as the rate of escape through the edge of a magnetic strip decreases slightly with increasing field, while these rates increase strongly for a skyrmion in the corresponding FMs when the field is directed antiparallel to the magnetization in the center of the skyrmion. A non-magnetic impurity is less likely to trap a skyrmion in AFMs especially in the presence of a magnetic field. This, together with the established fact that a spin polarized current moves skyrmions in AFMs in the direction of the current, while in FMs skyrmions move at an angle to the current, demonstrates that skyrmions in AFMs have several advantageous properties over skyrmions in FMs for memory and spintronic devices.