Source author record

P. Süle

P. Süle appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
3topics
3close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

9 published item(s)

preprint2015arXiv

Time-lapsed graphene moiré superlattice on Cu(111)

The detailed study of the graphene (gr) moiré superlattices emerging due to the mismatch between the substrate's and gr-overlayer crystal lattices is inevitable because of its high technological relevance. However, little is known about the dynamics of moiré superstructures on gr. Here, we report the first classical molecular dynamics simulation (CMD) of the moiré superlattice of graphene on Cu(111) using a new parameterized Tersoff-potential for the graphene/Cu(111) interface fitted in this paper to nonlocal van der Waals density functional theory (DFT) calculations. The interfacial force field with time-lapsed CMD provides superlattices in good quantitative agreement with the available experimental results. The long range coincidence supercells of $2 \times 2$ and $3 \times 3$ with nonequivalent moiré hills have also been identified and analyzed. The moiré superlattice exhibits a pattern which is dynamical rather than statically pinned to the support and can be observed mostly via time lapsing. The instantaneous snapshots of the periodic moiré pattern already at low temperature are weakly disordered lacking the apparent sharpness of the time averaged pattern and scanning tunneling microscopy images. This suggests the existence of competing orders between a static (1st order) and a dynamical (2nd order) moiré superstructures.The revealed random height fluctuations may limit the important electronic properties of supported graphene such as the mobility of charge carriers.

preprint2014arXiv

Anomalous nanoscale diffusion in Pt/Ti: superdiffusive intermixing

Probing the anomalous nanoscale intermixing using molecular dynamics (MD) simulations in Pt/Ti bilayer we characterize the superdiffusive nature of interfacial atomic transport. In particular, the low-energy ($0.5$ keV) ion-sputtering induced transient enhanced intermixing has been studied by MD simulations. {\em Ab initio} density functional calculations have been used to check and reparametrize the employed heteronuclear interatomic potential. We find a robust intermixing in Pt/Ti driven by nanoscale mass-anisotropy . The sum of the square of atomic displacements $\langle R^2 \rangle$ asymptotically scales nonlinearly ($\langle R^2 \rangle \propto t^2$), where $t$ is the time of ion-sputtering, respectively which is the fingerprint of superdiffusive features. This anomalous behavior explains the high diffusity tail in the concentration profile obtained by Auger electron spectroscopy depth profiling (AES-DP) analysis in Pt/Ti bilayer (reported in ref.: P. Süle, et al., J. Appl. Phys., 101, 043502 (2007)). In Ti/Pt bilayer a linear time scaling of $\langle R^2 \rangle \propto t$ has been found at the Ti/Pt interface indicating the suppression of superdiffusive features. These findings are inconsistent with the standard ion-mixing models. Instead a simple accelerative effect of the downward fluxes of energetic particles on the unidirectional fluxes of preferential intermixing of Pt atoms seems to explain the enhancement of interfacial broadening in Pt/Ti. Contrary to this in Ti/Pt the fluxes of recoils are in counterflow with intermixing Pt atoms and hence slows down the nanoscale mass-effect driven ballistic preferential mobility of Pt atoms.

preprint2014arXiv

Bond order potential for FeSi

A new parameter set has been derived for FeSi using the Albe-Erhart-type bond order potential (BOP) and the PONTIFIX code for fitting the parameters on a large training set of various polymorphs. Ab initio calculations are also carried out to study the relative stability of various polymorphs and to use the obtained data in the training set. The original BOP formalism was unable to account for the correct energetic relationship between the B20 ($ε$-FeSi) and B2 (CsCl) phases and notoriusly slightly favors incorrectly the B2 polymorph. In order to correct this improper behavior the BOP potential has been extended by a Columbic term (BOP+C) in order to account for the partial ionic character of FeSi. Using this potential we are able to account for the correct phase order between the most stable B20 and B2 (CsCl) polymorphs when the net atomic charges are properly set. Although this brings in a new somewhat uncertain parameter (the net charges) one can adjust properly the BOP+C potential for specific problems. To demonstrate this we study under high pressure the B2 phase which becomes more stable vs. B20 as it is found experimentally and expected to be taken place in the Earth mantle. The obtained BOP has also been tested for the metallic and semiconducting disilicides ($α$-FeSi$_2$ and $β$-FeSi$_2$) and for the Si/$β$-FeSi$_2$ heterostructure. The obtained BOP, as many other BOP, overestimates the melting point ($T_m$) of the B20 phase by $\sim 1000$ K if the parameters in the radial part of the potential were obtained according to the Pauling relation (regular way).

preprint2014arXiv

Fe impurity induced ion-nanopatterning: atomistic simulations using a new force field for FeSi

The ion-bombardment induced nanopatterning of Si(001) has been simulated by atomistic simulations with and without Fe impurity. The surface contamination has been simulated by using a new force field developed for FeSi. This is a fitted bond order potential (BOP) given for Si and for Fe by Albe et al. This BOP formula has been optimized simultaneously for FeSi, Si and for Fe. Using this new force field we are able to follow the ion-beam assisted deposition for Si in the presence of Fe contamination in the surface region. As an overall result, we get an unexpectedly rich variety of nanopatterns formed by the reorganization of the crater rims of the individual ion impacts. The previously thought simple atomistically roughened surfaces show unprecedented landscapes and topography with nanoscale features. The characteristic size of the units of the pattern is in the range of a few nms. Typical of the occurred pattern is the network of interconnected elongated adatom islands. We also see the self-organization of this pattern upon ion-bombardments. At $50^{\circ}$ impact angle we get a nanoporous surface (sponge-like) both for Fe-contaminated and Fe-free simulations. At $70^{\circ}$ of impact (grazing angle of incidence) the pattern resembles to that of elongated atomic chains (adatom islands) along the ion impact direction. This latter pattern could be understood as a prepattern state towards rippling. At lower angles ($30^{\circ}$) nanoholes rule the landsdcape. The obtained pattern corresponds to low fluence experiments which are used to consider as simple roughening without showing any sing of patterning.

preprint2014arXiv

The DFT and molecular dynamics multiscale study of the corrugation of graphene on Ru(0001): the unexpected stability of the moire-buckled structure

Results from first principles density functional theory (DFT) calculations and classical molecular dynamics (CMD) simulations are presented on moire-corrugation of graphene (gr). We find that the moire-corrugated graphene could be surprisingly stable against the perfectly flat gr-sheet as pointed out by CMD simulations and DFT calculations. We also show that using the cost-effective CMD approach one can simulate graphene on e.g. Ru(0001) with a correct binding registry and reasonable corrugation and adhesion energy. A new force field has been parameterized for the interface using an angular-dependent Abell-Tersoff potential. The newly parameterized Abell-Tersoff interface potential provides correct moire superstructures in accordance with scanning tunnelling microscopy images and with DFT results. Based on ab initio DFT calculations, we also find that the CMD moire superstructure can be used as a preoptimized structure for DFT calculations and for further geometry optimization. The nearly flat gr (the corrugation $ξ\approx 0.2$ $\hboxÅ$) on Ru(0001) is slightly energetically unfavorable vs. the moire-corrugated gr-system ($ξ\approx 2.0$ $\hboxÅ$) as revealed by van der Waals DFT structural relaxation.

preprint2013arXiv

The classical molecular dynamics simulation of graphene on Ru(0001) using a fitted Tersoff interface potential

The accurate molecular dynamics simulation of weakly bound adhesive complexes, such as supported graphene, is challenging due to the lack of an adequate interface potential. Instead of the widely used Lennard-Jones potential for weak and long-range interactions we use a newly parameterized Tersoff-potential for graphene/Ru(0001) system. The new interfacial force field provides adequate moir$\acute{e}$ superstructures in accordance with scanning tunnelling microscopy images and with DFT results. In particular, the corrugation of $ξ\approx 1.0 \pm 0.2$ $\hboxÅ$ is found which is somewhat smaller than found by DFT approaches ($ξ\approx 1.2$ $\hboxÅ$) and is close to STM measurements ($ξ\approx 0.8 \pm 0.3$ $\hboxÅ$).The new potential could open the way towards large scale simulations of supported graphene with adequate moir$\acute{e}$ supercells in many fields of graphene research. Moreover, the new interface potential might provide a new strategy in general for getting accurate interaction potentials for weakly bound adhesion in large scale systems in which atomic dynamics is inaccessible yet by accurate DFT calculations.

preprint2004arXiv

What is the real driving force of ion beam mixing?

Molecular dynamics simulations have been used to study the driving force of ion irradiation induced interfacial mixing in metal bilayers in which the relative mass of the constituents is considerable. We find no apparent effect of chemical forces, such as heat of mixing or cohesive energy up to 8 keV ion energy, although a considerable number of liquid and high energy particles (hot atoms) persist up to even 20 ps during the thermal spike. This result is in direct conflict with the widely accepted theory of thermal spike mixing (chemical interdiffusion model). The supersaturation of vacancies also occurs and which induces a thermally activated intermixing of the lighter constituent of the bilayer. The delay and the decoupling of the intermixing of the light constituent is explained as a backscattering effect at the interface: the interface acts as a diffusional barrier for high energy particles due to the large difference in atomic masses. The heavier atoms are predominantly ejected to the overlayer at the beginning of the thermal spike while the light atoms are injected to the bulk at the beginning of the cooling period (in Ti/Pt) or during the thermal spike with some time delay (Al/Pt).

preprint2003arXiv

Condensation energy, charge ordering and inter-layer coupling in cuprates

The correlation between the condensation energy and the critical temperature is studied within a charge ordered superlattice bilayer model in which pairing is supported by inter-layer Coulomb energy gain (potential energy driven superconductivity).The 2D pair-condensate can be characterized by a charge ordered state with a "checkerboard" like pattern seen by scanning tunneling microscopy.The drop of the c-axis dielectric screening can be the primary source of the condensation energy at optimal doping.We find that Coulomb energy gain occurs along the c-axis, which is proportional to the measured condensation energy ($U_0$) and to $T_c$: $E_c^{3D} \approx 2 (ξ_{ab}/a_0+1)^2 U_0 \approx k_B T_c$ and is due to inter-layer charge complementarity (charge asymmetry of the boson condensate)where $ξ_{ab}$ is the coherence length of the condensate and $a_0 \approx 3.9 Å$ is the in-plane lattice c onstant.The static c-axis dielectric constant $ε_c$ is calculated for various cuprates and compared with the available experimental data.