Peter Nielaba

Structure and conductance histogram of atomic-sized Au contacts

Many experiments have shown that the conductance histograms of metallic atomic-sized contacts exhibit a peak structure, which is characteristic of the corresponding material. The origin of these peaks still remains as an open problem. In order to shed some light on this issue, we present a theoretical analysis of the conductance histograms of Au atomic contacts. We have combined classical molecular dynamics simulations of the breaking of nanocontacts with conductance calculations based on a tight-binding model. This combination gives us access to crucial information such as contact geometries, forces, minimum cross section, total conductance, and transmission coefficients of the individual conduction channels. The ensemble of our results suggests that the low-temperature Au conductance histograms are a consequence of a subtle interplay between mechanical and electrical properties of these nanocontacts. At variance with other suggestions in the literature, our results indicate that the peaks in the Au conductance histograms are not a simple consequence of conductance quantization or the existence of exceptionally stable radii. We show that the main peak in the histogram close to one quantum of conductance is due to the formation of single-atom contacts and chains of gold atoms. Moreover, we present a detailed comparison with experimental results on Au atomic contacts where the individual channel transmissions have been determined.

A current-driven single-atom memory

The possibility of fabricating electronic devices with functional building blocks of atomic size is a major driving force of nanotechnology. The key elements in electronic circuits are switches, usually realized by transistors, which can be configured to perform memory operations. Electronic switches have been miniaturized all the way down to the atomic scale. However, at such scales, three-terminal devices are technically challenging to implement. Here we show that a metallic atomic-scale contact can be operated as a reliable and fatigue-resistant two-terminal switch. We apply a careful electromigration protocol to toggle the conductance of an aluminium atomic contact between two well-defined values in the range of a few conductance quanta. Using the nonlinearities of the current-voltage characteristics caused by superconductivity in combination with molecular dynamics and quantum transport calculations, we provide evidence that the switching process is caused by the reversible rearrangement of single atoms. Owing to its hysteretic behaviour with two distinct states, this two-terminal switch can be used as a non-volatile information storage element.

Theoretical analysis of the conductance histograms and structural properties of Ag, Pt, and Ni nanocontacts

Conductance histograms are a valuable tool to study the intrinsic conduction properties of metallic atomicsized contacts. These histograms show a peak structure, which is characteristic of the type of metal under investigation. Despite the enormous progress in the understanding of the electronic transport in metallic nanowires, the origin of this peak structure is still a basic open problem. In the present work we tackle this issue, extending our theoretical analysis of Au conductance histograms Dreher et al., Phys. Rev. B 72, 075435 2005 to different types of metals, namely, Ag, Pt and ferromagnetic Ni. We combine classical molecular dynamics simulations of the breaking of nanocontacts with conductance calculations based on a tight-binding model. This combination gives us access to crucial information such as contact geometries, strain forces, minimum cross-sections, the conductance, transmissions of the individual conduction channels, and, in the case of Ni, the spin polarization of the current. We shall also briefly discuss investigations of Al atomic-sized contacts. From our analysis we conclude that the differences in the histograms of these metals are due to i the very different electronic structures, which means different atomic orbitals contributing to the transport and ii the different mechanical properties, which in a case such as Pt lead to the formation of special structures, namely, monoatomic chains. Of particular interest are results for Ni that indicate the absence of any conductance quantization, and show how the current polarization evolves including large fluctuations from negative values in thick contacts to even positive values in the tunneling regime after rupture of the contact. Finally, we also present a detailed analysis of the breaking forces of these metallic contacts, which are compared to the forces predicted from bulk considerations.

The liquid-solid transition of hard discs: first-order transition or Kosterlitz-Thouless-Halperin-Nelson-Young scenario?

We consider the question of whether a two-dimensional hard-disc fluid has a first-order transition from the liquid state to the solid state as in the three-dimensional melting-crystallization transition or whether one has two subsequent continuous transitions, from the liquid to the hexatic phase and then to the solid phase, as proposed by Kosterlitz, Thouless, Halperin, Nelson and Young (KTHNY). Monte Carlo (MC) simulations of the fluid that study the growth of the bond orientational correlation length, and of the crystal are discussed. The emphasis is on a recent consistency test of the KTHNY renormalization group (RG) scenario, where MC simulations are used to estimate the bare elastic constants and dislocation fugacities in the solid, as a function of density, which then are used as starting values for the RG flow. This approach was validated earlier for the XY model as well.

Elastic Properties of 2D Colloidal Crystals from Video Microscopy

Elastic constants of two-dimensional (2D) colloidal crystals are determined by measuring strain fluctuations induced by Brownian motion of particles. Paramagnetic colloids confined to an air-water interface of a pendant drop are crystallized under the action of a magnetic field, which is applied perpendicular to the 2D layer. Using video microscopy and digital image processing we measure fluctuations of the microscopic strain obtained from random displacements of the colloidal particles from their mean (reference) positions. From these we calculate system-size dependent elastic constants, which are extrapolated using finite-size scaling to obtain their values in the thermodynamic limit. The data are found to agree rather well with zero-temperature calculations.

Protostellar collapse and fragmentation using an MHD gadget

Although the influence of magnetic fields is regarded as vital in the star formation process, only a few magnetohydrodynamics (MHD) simulations have been performed on this subject within the smoothed particle hydrodynamics method. This is largely due to the unsatisfactory treatment of non-vanishing divergence of the magnetic field. Recently smoothed particle magnetohydrodynamics (SPMHD) simulations based on Euler potentials have proven to be successful in treating MHD collapse and fragmentation problems, however these methods are known to have some intrinsic difficulties. We have performed SPMHD simulations based on a traditional approach evolving the magnetic field itself using the induction equation. To account for the numerical divergence, we have chosen an approach that subtracts the effects of numerical divergence from the force equation, and additionally we employ artificial magnetic dissipation as a regularization scheme. We apply this realization of SPMHD to a widely known setup, a variation of the ‘Boss and Bodenheimer standard isothermal test case’, to study the impact of the magnetic fields on collapse and fragmentation. In our simulations, we concentrate on setups, where the initial magnetic field is parallel to the rotation axis. We examine different field strengths and compare our results to other findings reported in the literature. We are able to confirm specific results found elsewhere, namely the delayed onset of star formation for strong fields, accompanied by the tendency to form only single stars. We also find that the ‘magnetic cushioning effect’, where the magnetic field is wound up to form a ‘cushion’ between the binary, aids binary fragmentation in a case where previously only formation of a single protostar was expected.

Simulation of the shape memory effect in a NiTi nano model system

Abstract The shape memory behavior of a NiTi nanoparticle is analyzed by molecular dynamics simulations. After a detailed description of the equilibrium structures of the used model potential, the multi variant martensitic ground state, which depends on the geometry of the particle, is discussed. Tensile load is applied, changing the variant configuration to a single domain state with a remanent strain after unloading. Heating the particle leads to a shape memory effect without a phase transition to the austenite, but by variant reorientation and twin boundary formation at a certain temperature. These processes are described by stress–strain and strain–temperature curves, together with a visualization of the microstructure of the nanoparticle. Results are presented for five different Ni concentrations in the vicinity of 50%, showing for example, that small deviations from this ideal composition can influence the critical temperature for shape recovery significantly.

Layer reduction in driven 2D-colloidal systems through microchannels

The transport behavior of a system of gravitationally driven superparamagnetic colloidal particles is investigated. The motion of the particles through a narrow channel is governed by magnetic dipole interactions, and a layered structure forms parallel to the walls. The arrangement of the particles is perturbed by diffusion and the motion induced by gravity leading to a density gradient along the channel. Our main result is the reduction of the number of layers. Experiments and Brownian dynamics simulations show that this occurs due to the density gradient along the channel.

Diffusional Relaxation in Dimer Deposition

In deposition of dimers on a 1D lattice substrate, we find by analytical arguments, supported by numerical Monte Carlo simulations, that the effect of added diffusional relaxation is to allow the full, saturation coverage, 100%, for large times. This limiting coverage is approached according to the ~ 1/√t power law preceded, for fast diffusion, by the mean-field crossover regime with the intermediate ~ 1/t behavior.

Molecular dynamics study of the thermopower of Ag, Au, and Pt nanocontacts

Using molecular dynamics simulations of many junction stretching processes combined with tight-bindingbased electronic structure and transport calculations, we analyze the thermopower of silver (Ag), gold (Au), and platinum (Pt) atomic contacts. In all cases we observe that the thermopower vanishes on average within the standard deviation and that its fluctuations increase for a decreasing minimum cross section of the junctions. However, we find a suppression of the fluctuations of the thermopower for the s-valent metals Ag and Au, when the conductance originates from a single, perfectly transmitting channel. Essential features of the experimental results for Au, Ag, and copper (Cu) of Ludoph and van Ruitenbeek [Phys. Rev. B 59, 12290 (1999)], as yet unaddressed by atomistic studies, can hence be explained by considering the atomic and electronic structure at the disordered narrowest constriction of the contacts. For the multivalent metal Pt our calculations predict the fluctuations of the thermopower to be larger by one order of magnitude as compared to Ag and Au, and suppressions of the fluctuations as a function of the conductance are absent. Main features of our results are explained in terms of an extended single-level model.