Abstract: A multi-species Monte Carlo (MC) model, combined with an analytical surface model, has been developed in order to investigate the general plasma processes occurring during the sputter deposition of complex oxide films in a dual-magnetron sputter deposition system. The important plasma species, such as electrons, Ar+ ions, fast Ar atoms and sputtered metal atoms (i.e. Mg and Al atoms) are described with the so-called multi-species MC model, whereas the deposition of MgxAlyOz films is treated by an analytical surface model. Targetsubstrate distances for both magnetrons in the dual-magnetron setup are varied for the purpose of growing stoichiometric complex oxide thin films. The metal atoms are sputtered from pure metallic targets, whereas the oxygen flux is only directed toward the substrate and is high enough to obtain fully oxidized thin films but low enough to avoid target poisoning. The calculations correspond to typical experimental conditions applied to grow these complex oxide films. In this paper, some calculation results are shown, such as the densities of various plasma species, their fluxes toward the targets and substrate, the deposition rates, as well as the film stoichiometry. Moreover, some results of the combined model are compared with experimental observations. Note that this is the first complete model, which can be applied for large and complicated magnetron reactor geometries, such as dual-magnetron configurations. With this model, we are able to describe all important plasma species as well as the deposition process. It can also be used to predict film stoichiometries of complex oxide films on the substrate.

Abstract: In this paper, an overview is given of modeling investigations carried out in our research group for a better understanding of plasmas used for medical, environmental and nano applications. The focus is both on modeling the plasma chemistry and the plasma-surface interactions. The plasma chemistry provides the densities and fluxes of the important plasma species. This information can be used as input when modeling the plasma-surface interactions. The combination of plasma simulations and plasma – surface interaction simulations provides a more comprehensive understanding of the underlying processes for these applications.

Abstract: The interaction of a wave packet (and in particular the wave front) with a mass barrier is investigated in one dimension. We discuss the main features of the wave packet that are inherent to two-dimensional wave packets, such as compression during reflection, penetration in the case when the energy is lower than the height of the barrier, waving tails, precursors, and the retardation of the reflected and penetrated wave packets. These features depend on the wave-packet envelope function which we demonstrate by considering the case of a rectangular wave packet with sharp front and trailing edges and a smooth Gaussian wave packet. The method of Fourier integral for obtaining the nonstationary solutions is used.

Abstract: The electrical transport properties of a MoGe thin film with a honeycomb array of nanoscale holes are investigated. The critical current of the system shows nonmatching anomalies as a function of applied magnetic field, enabling us to distinguish between multiquanta vortices trapped in the holes and interstitial vortices located between the holes. The number of vortices trapped in each hole is found to be larger than the saturation number predicted for an isolated hole and shows a nonlinear field dependence, leading to the caging effect as predicted from the Ginzburg-Landau (GL) theory. Our experimental results are supplemented by numerical simulations based on the GL theory.

Abstract: We theoretically find that the resistance of a superconducting microbridge or nanowire decreases while the retrapping current I(r) for the transition to the superconducting state increases when one suppresses the magnitude of the order parameter vertical bar Delta vertical bar in the attached superconducting leads. This effect is a consequence of the increased energy interval for diffusion of the “hot” nonequilibrium quasiparticles (induced by the oscillations of vertical bar Delta vertical bar in the center of the microbridge) to the leads. The effect is absent in short microbridges (with length less than the coherence length) and it is relatively weak in long microbridges (with length larger than the inelastic relaxation length of the nonequilibrium distribution function). A nonmonotonous dependence of I(r) on the length of the microbridge is predicted. Our results are important for the explanation of the enhancement of the critical current and the appearance of negative magnetoresistance observed in many recent experiments on superconducting microbridges or nanowires.

Abstract: Iron carbides play a crucial role in steel manufacturing and processing and to a large extent determine the physical properties of steel products. The modified embedded atom method (MEAM) in combination with Lee's Fe-C potential is a good candidate for molecular dynamics simulations on larger Fe-C systems. Here, we investigate the stability and crystal structures of pure iron and binary iron carbides using MEAM and compare them with the experimental data and quantum-mechanical density functional theory calculations. The analysis shows that the Fe-C potential gives reasonable results for the relative stability of iron and iron carbides. The performance of MEAM for the prediction of the potential energy and the calculated lattice parameters at elevated temperature for pure iron phases and cementite are investigated as well. The conclusion is that Lee's MEAM Fe-C potential provides a promising basis for further molecular dynamics simulations of Fe-C alloys and steels at lower temperatures (up to 800 K).

Abstract: The ground-state properties of a paramagnetic Mott insulator at half-filling are investigated in the presence of an external electric field using the inhomogeneous Gutzwiller approximation for a single-band Hubbard model in a slab geometry. We find that the metal-insulator transition is shifted toward higher Hubbard repulsions by applying an electric field perpendicular to the slab. The main reason is the accumulation of charges near the surface. The spatial distribution of site-dependent quasiparticle weight shows that it is maximal in a few layers beneath the surface, while the central sites where the field is screened have a very low quasiparticle weight. Our results show that above a critical-field value, states near the surface will be metallic, while the bulk quasiparticle weight is extremely suppressed but never vanishing, even for large Hubbard repulsions above the bulk zero-field critical value. Below the critical-field value, our results hint toward an insulating state in which the electric field is totally screened and the slab is again at half-filling.

Abstract: We study nonlinear flux dynamics in a current-carrying type-I superconductor. The stray magnetic field of the current induces the intermediate state, where nucleation of flux domains is discretized to a single fluxoid at a time, while their final shape (tubular or laminar), size, and nucleation rate depend on applied current and edge conditions. The current induces opposite flux domains on opposite sides of the sample, and subsequently drives them to annihilation-which is also discretized, as a sequence of vortex-antivortex pairs. The discretization of both nucleation and annihilation leaves measurable traces in the voltage across the sample and in locally probed magnetization. The reported dynamic phenomena thus provide an unambiguous proof of a flux quantum being the smallest building block of the intermediate state in type-I superconductors.

Abstract: We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain epsilon(zz), which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K point. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

Abstract: The coupling of angstrom-sized electron probes with spin-polarized electronic transitions shows that the inelastically scattered probe electron is in a mixed state containing electron vortices with nonzero orbital angular momentum. These electrons create an asymmetric intensity distribution in energy filtered diffraction patterns, giving access to maps of the magnetic moments with atomic resolution. A feasibility experiment shows evidence of the predicted effect. Potential applications are column-by-column maps of magnetic ordering, and the creation of angstrom-sized free electrons with orbital angular momentum by inelastic scattering in a thin ferromagnetic foil.

Abstract: Superconductivity in nanowires made of weak coupling superconductor materials is investigated using a Green function approach. We show that these are multigap systems in which the ratio Delta(T)/k(B)T(c) is to a large extent similar to what is observed in some high-T-c two-gap systems, such as MgB2 and some of the Fe-based superconductors. On the other hand, because of confinement, the superfluid density has a temperature behavior of the form n(s)(T) = 1 – (T/T-c)(3) near T-c, thus deviating from the BCS behavior for bulk superconductors.

Abstract: The effective Hamiltonian for an electron in a quasi-one-dimensional quantum ring in the presence of spin-orbit interactions is derived. We demonstrate that, when both coupling types are simultaneously present, the effective Hamiltonian derived by the lowest-radial-state approximation produces energy spectra and charge densities which deviate strongly from the exact ones. For equal Rashba and Dresselhaus coupling constants the lowest-radial-state approximation opens artifactal avoided crossings in the energy spectra and deforms the circular symmetry of the confined charge densities. In this case, there does not exist a ring thin enough to justify the restriction to the lowest radially quantized energy state. We derive the effective Hamiltonian accounting for both the lowest and the first excited radial states, and show that the inclusion of the latter restores the correct features of the exact solution. Relation of this result to the states of a quantum wire is also discussed.

Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: We utilize an exact variational numerical procedure to calculate the ground-state properties of a polaron in the presence of Rashba and linear Dresselhaus spin-orbit coupling. We find that when the linear Dresselhaus spin-orbit coupling approaches the Rashba spin-orbit coupling, the Van Hove singularity in the density of states will be shifted away from the bottom of the band and finally disappear when the two spin-orbit couplings are tuned to be equal. The effective mass will be suppressed; the trend will become more significant for low phonon frequency. The presence of two dominant spin-orbit couplings will make it possible to tune the effective mass with more varied observables.