Researcher profile

Stefanos Papanikolaou

Stefanos Papanikolaou contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
27works
0followers
18topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

27 published item(s)

preprint2022arXiv

Yielding in multi-component metallic glasses: Universal signatures of elastic modulus heterogeneities

Sheared multi-component bulk metallic glasses are characterized by both chemical and structural disorder that define their properties. We investigate the behavior of the local, microstructural elastic modulus across the plastic yielding transition in six Ni-based multi-component glasses, that are characterized by compositional features commonly associated with solid solution formability. We find that elastic modulus fluctuations display consistent percolation characteristics pointing towards universal behavior across chemical compositions and overall yielding sharpness characteristics. Elastic heterogeneity grows upon shearing via the percolation of elastically soft clusters within an otherwise rigid amorphous matrix, confirming prior investigations in granular media and colloidal glasses. We find clear signatures of percolation transition with spanning clusters that are universally characterized by scale-free characteristics and critical scaling exponents. The spatial correlation length and mean cluster size tend to diverge prior to yielding, with associated critical exponents that exhibit fairly weak dependence on compositional variations as well as macroscopic stress-strain curve details.

preprint2021arXiv

Shear banding instability in \remove{high entropy} multi-component metallic glasses: Interplay of composition and short-range order}

The shear-banding instability in quasi-statically driven bulk metallic glasses emerges from collective dynamics, mediated by shear transformation zones and associated non-local elastic interactions. It is also phenomenologically known that sharp structural features of shear bands are typically correlated to the sharpness of the plastic yielding transition, being predominant in commonly studied alloys composed of multiple different elements, that have very different atomic radii. However, in the opposite limit \remove{of high-entropy multicomponent alloys,} where elements' radii are relatively similar, plastic yielding of bulk metallic glasses is highly dependent on compositional and ordering features. In particular, a known mechanism at play involves the formation of short-range order dominated by icosahedra-based clusters. Here, we report on atomistic simulations of multi-component metallic glasses with different chemical compositions showing that the degree of strain localization is largely controlled by the interplay between composition-driven icosahedra-ordering and collectively-driven shear transformation zones. By altering compositions, strain localization ranges from diffuse homogenized patterns to singular crack-like features. We quantify the dynamical yielding transition by measuring the atoms' susceptibility to plastic rearrangements, strongly correlated to the local atomic structure. We find that the abundance of short-range ordering of icosahedra within rearranging zones increases glassy materials' capacity to delocalize strain. The kind of plastic yielding can be often qualitatively inferred by the commonly used compositional descriptor that characterizes element associations, the misfit parameter $δ_a$, and also by uncommon ones, such as shear-band width and shear-band dynamics' correlation parameters.

preprint2020arXiv

Inference and Prediction of Nanoindentation Response in FCC Crystals: Methods and Discrete Dislocation Simulation Examples

Nanoindentation, a common technique for probing the mechanical properties of crystalline materials, exhibits both surface and bulk-dominated responses that are linked to the parent crystal's elasticity and plasticity. For FCC crystals, the nanoindentation response is primarily controlled by the indented-grain's crystalline orientation and dislocation density. However, the indentation properties, such as hardness, have ambiguous significance at the nanoscale due to tip-specific size effects. Through extensive discrete dislocation dynamics simulations of flat-punch nanoindentation on FCC single crystals, we show that a one-to-one correspondence between an indented location and the pre-existing dislocation density is possible with the help of statistical, machine-learning algorithms. We cluster and classify various experimentally plausible ensembles with varying dislocation density and/or crystalline orientation, using nanoindentation force-depth curves or post-indent surface displacement images. We use this classification to statistically predict the nanoindentation response of a given location, through the average response of classified load-displacement curves. In addition, we demonstrate that our simulation results qualitatively agree with experimental data on common FCC single crystals. Ultimately, we propose a novel pathway to uncover the depth dependence of local surface material properties for more complete, non-destructive, microstructural characterization, through nanoindentation.

preprint2019arXiv

The $Λ$-invariant and topological pathways to influence sub-micron strength and crystal plasticity

In small volumes, sample dimensions are known to strongly influence mechanical behavior, especially strength and crystal plasticity. This correlation fades away at the so-called mesoscale, loosely defined at several micrometers in both experiments and simulations. However, this picture depends on the entanglement of the initial defect configuration. In this paper, we study the effect of sample dimensions with a full control on dislocation topology, through the use of a novel observable for dislocation ensembles, the $Λ$-invariant, that depends only on mutual dislocation linking: It is built on the natural vortex character of dislocations and it has a continuum/discrete correspondence that may assist multiscale modeling descriptions. We investigate arbitrarily complex initial dislocation microstructures in sub-micron-sized pillars, using three-dimensional discrete dislocation dynamics simulations for finite volumes. We demonstrate how to engineer nanoscale dislocation ensembles that appear virtually independent from sample dimensions, either by biased-random dislocation loop deposition or by sequential mechanical loads of compression and torsion.

preprint2019arXiv

The Emergence and Role of Dipolar Dislocation Patterns in Discrete and Continuum Formulations of Plasticity

The plasticity transition at the yield strength of a crystal typically signifies the tendency of dislocation defects towards relatively unrestricted motion. For an isolated dislocation the motion is in the slip plane with velocity proportional to the shear stress, while due to the long range interaction dislocation ensembles move towards satisfying emergent collective elastoplastic modes. Such collective motions have been discussed in terms of the elusively defined backstress. In this paper, we develop a two-dimensional stochastic continuum dislocation dynamics theory that clarifies the role of backstress and demonstrates a precise agreement with the collective behavior of its discrete counterpart, as a function of applied load and with only three essential free parameters. The main ingredients of the continuum theory is the evolution equations of statistically stored and geometrically necessary dislocation densities, which are driven by the long-range internal stress, a stochastic flow stress term and, finally, two local diffusion-like terms. The agreement is shown primarily in terms of the patterning characteristics that include the formation of dipolar dislocation walls.

preprint2016arXiv

Obstacles and sources in dislocation dynamics: Strengthening and statistics of abrupt plastic events in nanopillar compression

Mechanical deformation of nanopillars displays features that are distinctly different from the bulk behavior of single crystals: Yield strength increases with decreasing size and plastic deformation comes together with strain bursts or/and stress drops (depending on loading conditions) with a very strong sensitivity of the stochasticity character on material preparation and conditions. The character of the phenomenon is standing as a paradox: While these bursts resemble the universal, widely independent of material conditions, noise heard in bulk crystals using acoustic emission (AE) techniques, they strongly emerge primarily with decreasing size and increasing strength in nanopillars. In this paper, we present a realistic but minimal discrete dislocation plasticity model for the elasto-plastic deformation of nanopillars that is consistent with the main experimental observations of nano pillar compression experiments and provides a clear insight to this paradox. With increasing sample size, the model naturally transitions between the typical progressive behavior of nanopillars to a behavior that resembles evidence for bulk mesoscale plasticity. The combination of consistent strengthening, large flow stress fluctuations and critical avalanches is only observed in the {\it depinning regime} where obstacles are much stronger than dislocation sources; in contrast, when dislocation source strength becomes comparable to obstacle barriers, then yield strength size effects are absent but plasticity avalanche dynamics is strongly universal, across sample width and aspect-ratio scales. Finally, we elucidate the mechanism that leads to the connection between depinning and size effects in our model dislocation dynamics. In this way, our model builds a way towards unifying statistical aspects of mechanical deformation across scales.

preprint2015arXiv

"Irregularization" of Systems of Conservation Laws

We explore new ways of regulating defect behavior in systems of conservation laws. Contrary to usual regularization schemes (such as a vanishing viscosity limit), which attempt to control defects by making them smoother, our schemes result in defects which are \textit{more singular}, and we thus refer to such schemes as "irregularizations". In particular, we seek to produce \textit{delta shock} defects which satisfy a condition of \textit{stationarity}. We are motivated to pursue such exotic defects by a physical example arising from dislocation dynamics in materials physics, which we describe.

preprint2014arXiv

Ising nematic fluid phase of hard-core dimers on the square lattice

We present a model of classical hard-core dimers on the square lattice that contains an Ising nematic phase in its phase diagram. We consider a model with an attractive interaction for parallel dimers on a given plaquette of the square lattice and an attractive interaction for neighboring parallel dimers on the same row ({\it viz} column) of the lattice. By extensive Monte carlo simulations we find that with a finite density of holes the phase diagram has, with rising temperatures, a columnar crystalline phase, an Ising nematic liquid phase and a disordered fluid phase, separated by Ising continuous phase transitions. We present strong evidence for the Ising universality class of both transitions. The Ising nematic phase can be interpreted as either an intermediate classical thermodynamic phase (possibly of a strongly correlated antiferromagnet) or as a phase of a 2D quantum dimer model using the Rokhsar-Kivelson construction of exactly solvable quantum Hamiltonians.

preprint2014arXiv

Shearing a glass and the role of pin-delay in models of interface depinning

When a disordered solid is sheared, yielding is followed by the onset of intermittent response that is characterized by slip in local regions usually labeled shear-transformation zones (STZ). Such intermittent response resembles the behavior of earthquakes or contact depinning, where a well-defined landscape of pinning disorder prohibits the deformation of an elastic medium. Nevertheless, a disordered solid is evidently different in that pinning barriers of particles are due to neighbors that are also subject to motion. Microscopic yielding leads to destruction of the local microstructure and local heating. It is natural to assume that locally a liquid emerges for a finite timescale before cooling down to a transformed configuration. For including this characteristic transient in glass depinning models, we propose a general mechanism that involves a "pin-delay" time Tpd, during which each region that slipped evolves as a fluid. The new timescale can be as small as a single avalanche time-step. We demonstrate that the inclusion of this mechanism causes a drift of the critical exponents towards higher values for the slip sizes $τ$, until a transition to permanent shear-banding behavior happens causing almost oscillatory, stick-slip response. Moreover, it leads to a proliferation of large events that are highly inhomogeneous and resemble sharp slip band formation. Our model appears to be qualitatively consistent with recent experiments and simulations of disordered solids under shear.

preprint2014arXiv

The statistics of frictional families

We develop a theoretical description for mechanically stable frictional packings in terms of the difference between the total number of contacts required for isostatic packings of frictionless disks and the number of contacts in frictional packings, $m=N_c^0-N_c$. The saddle order $m$ represents the number of unconstrained degrees of freedom that a static packing would possess if friction were removed. Using a novel numerical method that allows us to enumerate disk packings for each $m$, we show that the probability to obtain a packing with saddle order $m$ at a given static friction coefficient $μ$, $P_m(μ)$, can be expressed as a power-series in $μ$. Using this form for $P_m(μ)$, we quantitatively describe the dependence of the average contact number on friction coefficient for static disk packings obtained from direct simulations of the Cundall-Strack model for all $μ$ and $N$.

preprint2013arXiv

Classical Topological Order in Abelian and Non-Abelian Generalized Height Models

We present Monte Carlo simulations on a new class of lattice models in which the degrees of freedom are elements of an abelian or non-abelian finite symmetry group G, placed on directed edges of a two-dimensional lattice. The plaquette group product is constrained to be the group identity. In contrast to discrete gauge models (but similar to past work on height models) only elements of symmetry-related subsets S of G are allowed on edges. These models have topological sectors labeled by group products along topologically non-trivial loops. Measurement of relative sector probabilities and the distribution of distance between defect pairs are done to characterize the types of order (topological or quasi-LRO) exhibited by these models. We present particular models in which fully local non-abelian constraints lead to global topological liquid properties.

preprint2013arXiv

Computational studies of the glass-forming ability of model bulk metallic glasses

Bulk metallic glasses (BMGs) are produced by rapidly thermally quenching supercooled liquid metal alloys below the glass transition temperature at rates much faster than the critical cooling rate R_c below which crystallization occurs. The glass-forming ability of BMGs increases with decreasing R_c, and thus good glass-formers possess small values of R_c. We perform molecular dynamics simulations of binary Lennard-Jones (LJ) mixtures to quantify how key parameters, such as the stoichiometry, particle size difference, attraction strength, and heat of mixing, influence the glass-formability of model BMGs. For binary LJ mixtures, we find that the best glass-forming mixtures possess atomic size ratios (small to large) less than 0.92 and stoichiometries near 50:50 by number. In addition, weaker attractive interactions between the smaller atoms facilitate glass formation, whereas negative heats of mixing (in the experimentally relevant regime) do not change R_c significantly. These studies represent a first step in the development of computational methods for quantitatively predicting glass-formability.

preprint2013arXiv

Isostaticity at Frictional Jamming

Amorphous packings of frictionless, spherical particles are isostatic at jamming onset, with the number of constraints (contacts) equal to the number of degrees of freedom. Their structural and mechanical properties are controlled by the interparticle contact network. In contrast, amorphous packings of frictional particles are typically hyperstatic at jamming onset. We perform extensive numerical simulations in two dimensions of the geometrical asperity (GA) model for static friction, to further investigate the role of isostaticity. In the GA model, interparticle forces are obtained by summing up purely repulsive central forces between periodically spaced circular asperities on contacting grains. We compare the packing fraction, contact number, mobilization distribution, and vibrational density of states using the GA model to those generated using the Cundall-Strack (CS) approach. We find that static packings of frictional disks obtained from the GA model are mechanically stable and isostatic when we consider interactions between asperities on contacting particles. The crossover in the structural and mechanical properties of static packings from frictionless to frictional behavior as a function of the static friction coefficient coincides with a change in the type of interparticle contacts and the disappearance of a peak in the density of vibrational modes for the GA model. These results emphasize that mesoscale features of the model for static friction play an important role in determining the properties of granular packings.

preprint2013arXiv

Quasi-periodic events in crystal plasticity and the self-organized avalanche oscillator

When external stresses in a system - physical, social or virtual - are relieved through impulsive events, it is natural to focus on the attributes of these avalanches. However, during the quiescent periods in between, stresses may be relieved through competing processes, such as slowly flowing water between earthquakes or thermally activated dislocation flow between plastic bursts. Such unassuming, smooth responses can have dramatic effects on the avalanche properties. Here we report a thorough experimental investigation of slowly compressed Ni microcrystals, covering three orders of magnitude in nominal strain rate, that exhibits unconventional quasi-periodic avalanche bursts and higher critical exponents as the strain rate is decreased. Our analytic and computational study naturally extends dislocation avalanche modeling to incorporate dislocation relaxations and reveals the emergence of the self-organized avalanche oscillator, a novel critical state exhibiting oscillatory approaches toward a depinning critical point. We demonstrate that the predictions of our theory are faithfully exhibited in our experiments.

preprint2013arXiv

Scaling theory of continuum dislocation dynamics in three dimensions: Self-organized fractal pattern formation

We focus on mesoscopic dislocation patterning via a continuum dislocation dynamics theory (CDD) in three dimensions (3D). We study three distinct physically motivated dynamics which consistently lead to fractal formation in 3D with rather similar morphologies, and therefore we suggest that this is a general feature of the 3D collective behavior of geometrically necessary dislocation (GND) ensembles. The striking self-similar features are measured in terms of correlation functions of physical observables, such as the GND density, the plastic distortion, and the crystalline orientation. Remarkably, all these correlation functions exhibit spatial power-law behaviors, sharing a single underlying universal critical exponent for each type of dynamics.

preprint2013arXiv

Universal properties of magnetization dynamics in polycrystalline ferromagnetic films

We investigate the scaling behavior in the statistical properties of Barkhausen noise in ferromagnetic films. We apply the statistical treatment usually employed for bulk materials in experimental Barkhausen noise time series measured with the traditional inductive technique in polycrystalline ferromagnetic films having different thickness from 100 up to 1000 nm, and investigate the scaling exponents. Based on this procedure, we can group the samples in a single universality class, characterized by exponents τ\sim 1.5, α\sim 2.0, and 1/σνz \sim \vartheta \sim 2.0. We interpret these results in terms of theoretical models and provide experimental evidence that a well-known mean-field model for the dynamics of a ferromagnetic domain wall in three-dimensional ferromagnets can be extended for films. We identify that the films present an universal three-dimensional magnetization dynamics, governed by long-range dipolar interactions, even at the smallest thicknesses, indicating that the two-dimensional magnetic behavior commonly verified for films cannot be generalized for all thickness ranges.

preprint2011arXiv

Avalanche Spatial Structure and Multivariable Scaling Functions; Sizes, Heights, Widths, and Views through Windows

We introduce a systematic method for extracting multivariable universal scaling functions and critical exponents from data. We exemplify our insights by analyzing simulations of avalanches in an interface using simulations from a driven quenched Kardar-Parisi-Zhang (qKPZ) equation. We fully characterize the spatial structure of these avalanches- we report universal scaling functions for size, height and width distributions, and also local front heights. Furthermore, we resolve a problem that arises in many imaging experiments of crackling noise and avalanche dynamics, where the observed distributions are strongly distorted by a limited field of view. Through artificially windowed data, we show these distributions and their multivariable scaling functions may be written in terms of two control parameters, the window size and the characteristic length scale of the dynamics. For the entire system and the windowed distributions we develop accurate parameterizations for the universal scaling functions, including corrections to scaling and systematic error bars, facilitated by a novel software environment SloppyScaling.

preprint2011arXiv

Dielectric breakdown and avalanches at non-equilibrium metal-insulator transitions

Motivated by recent experiments on the finite temperature Mott transition in VO2 films, we propose a classical coarse-grained dielectric breakdown model where each degree of freedom represents a nanograin which transitions from insulator to metal with increasing temperature and voltage at random thresholds due to quenched disorder. We describe the properties of the resulting non-equilibrium metal-insulator transition and explain the universal characteristics of the resistance jump distribution. We predict that by tuning voltage, another critical point is approached, which separates a phase of "bolt"-like avalanches from percolation-like ones.

preprint2011arXiv

Straining the Identity of Majorana Fermions

We propose an experimental setup of an interferometer for the observation of neutral Majorana fermions on topological insulator - superconductor - ferromagnet junctions. We show that the extended lattice defects naturally present in materials, dislocations, induce spin currents on the edges while keeping the bulk time-reversal symmetry intact. We propose a simple two terminal conductance measurement in an interferometer formed by two edge point contacts, which reveals the nature of Majorana states through the effect of dislocations. The zero temperature magneto-conductance changes from even oscillations with period phi/2 (phi is the flux quantum hc/e) to odd oscillations with period phi, when non-trivial dislocations are present and the Majorana states are sufficiently strongly coupled. Additionally, the conductance acquires a notable asymmetry as a function of the incident electron energy, due to the topological influence of the dislocations, while resonances appear at the coupling energy of Majorana states.

preprint2010arXiv

A minimal model of plasma membrane heterogeneity requires coupling cortical actin to criticality

We present a minimal model of plasma membrane heterogeneity that combines criticality with connectivity to cortical cytoskeleton. Our model is motivated by recent observations of micron-sized critical fluctuations in plasma membrane vesicles that are detached from their cortical cytoskeleton. We incorporate criticality using a conserved order parameter Ising model coupled to a simple actin cytoskeleton interacting through point-like pinning sites. Using this minimal model, we recapitulate several experimental observations of plasma membrane raft heterogeneity. Small (r~20nm) and dynamic fluctuations at physiological temperatures arise from criticality. Including connectivity to cortical cytoskeleton disrupts large fluctuations, prevents macroscopic phase separation at low temperatures (T<=22°C), and provides a template for long lived fluctuations at physiological temperature (T=37°C). Cytoskeleton-stabilized fluctuations produce significant barriers to the diffusion of some membrane components in a manner that is weakly dependent on the number of pinning sites and strongly dependent on criticality. More generally, we demonstrate that critical fluctuations provide a physical mechanism to organize and spatially segregate membrane components by providing channels for interaction over large distances.

preprint2010arXiv

A Striped Holographic Superconductor

We study inhomogeneous solutions of a 3+1-dimensional Einstein-Maxwell-scalar theory. Our results provide a holographic model of superconductivity in the presence of a charge density wave sourced by a modulated chemical potential. We find that below a critical temperature superconducting stripes develop. We show that they are thermodynamically favored over the normal state by computing the grand canonical potential. We investigate the dependence of the critical temperature on the modulation&#39;s wave vector, which characterizes the inhomogeneity. We find that it is qualitatively similar to that expected for a weakly coupled BCS theory, but we point out a quantitative difference. Finally, we use our solutions to compute the conductivity along the direction of the stripes.

preprint2010arXiv

Bending crystals: Emergence of fractal dislocation structures

We provide a minimal continuum model for mesoscale plasticity, explaining the cellular dislocation structures observed in deformed crystals. Our dislocation density tensor evolves from random, smooth initial conditions to form self-similar structures strikingly similar to those seen experimentally - reproducing both the fractal morphologies and some features of the scaling of cell sizes and misorientations analyzed experimentally. Our model provides a framework for understanding emergent dislocation structures on the mesoscale, a bridge across a computationally demanding mesoscale gap in the multiscale modeling program, and a new example of self-similar structure formation in non-equilibrium systems.

preprint2010arXiv

Beyond power laws: Universality in the average avalanche shape

We report the measurement of multivariable scaling functions for the temporal average shape of Barkhausen noise avalanches, and show that they are consistent with the predictions of simple mean-field theories. We bypass the confounding factors of time-retarded interactions (eddy currents) by measuring thin permal- loy films, and bypass thresholding effects and amplifier distortions by applying Wiener deconvolution. We find experimental shapes that are approximately symmetric, and track the evolution of the scaling function. We solve a mean- field theory for the magnetization dynamics and calculate the form of the scaling function in the presence of a demagnetizing field and a finite field ramp-rate, yielding quantitative agreement with the experiment.

preprint2010arXiv

Electronic States of Graphene Grain Boundaries

We introduce a model for amorphous grain boundaries in graphene, and find that stable structures can exist along the boundary that are responsible for local density of states enhancements both at zero and finite (~0.5 eV) energies. Such zero energy peaks in particular were identified in STS measurements [J. Červenka, M. I. Katsnelson, and C. F. J. Flipse, Nature Physics 5, 840 (2009)], but are not present in the simplest pentagon-heptagon dislocation array model [O. V. Yazyev and S. G. Louie, Physical Review B 81, 195420 (2010)]. We consider the low energy continuum theory of arrays of dislocations in graphene and show that it predicts localized zero energy states. Since the continuum theory is based on an idealized lattice scale physics it is a priori not literally applicable. However, we identify stable dislocation cores, different from the pentagon-heptagon pairs, that do carry zero energy states. These might be responsible for the enhanced magnetism seen experimentally at graphite grain boundaries.

preprint2010arXiv

First- order versus unconventional phase transitions in three-dimensional dimer models

We study the phase transition between the Coulomb liquid and the columnar crystal in the 3D classical dimer model, which was found to be continuous in the O(3) universality class. In addition to nearest neighbor interactions which favor parallel dimers, further neighbor interactions are allowed in such a manner that the cubic symmetry of the original system remains intact. We show that the transition in the presence of weak additional, symmetry preserving interactions is first-order. However the universality class of the transition remains continuous when the additional interactions are weakly repulsive. In this way, we verify the existence of a multicritical point near the unperturbed transition and we identify a critical line of unconventional transitions between the Coulomb liquid phase and the $6-$fold columnar phase.

preprint2009arXiv

Nodal-antinodal dichotomy from pairing disorder in d-wave superconductors

We study the basic features of the local density of states (LDOS) observed in STM experiments on high-T$_c$ d-wave superconductors in the context of a minimal model of a d-wave superconductor which has {\it weakly} modulated off-diagonal disorder. We show that the low and high energy features of the LDOS are consistent with the observed experimental patterns and in particular, the anisotropic local domain features at high energies. At low energies, we obtain not only the scattering peaks predicted by the octet model, but also weak features that should be experimentally accessible. Finally, we show that the emerging features of the LDOS lose their correspondence with the features of the imposed disorder, as its complexity increases spatially.