Showing 1–50 of 137 results for author: Alava, M

Atomistic-Level Analysis of Nanoindentation-Induced Plasticity in Arc--Melted NiFeCrCo Alloys: The role of stacking faults

Authors: F. J. Dominguez-Gutierrez, A. Olejarz, M. Landeiro Dos Reis, E. Wyszkowska, D. Kalita, W. Y. Huo, I. Jozwik, L. Kurpaska, S. Papanikolaou, M. J. Alava, K. Muszka

Abstract: Concentrated solid solution alloys (CSAs) have attracted attention for their promising properties; however, current manufacturing methods face challenges in complexity, high costs, and limited scalability, raising concerns about industrial viability. The prevalent technique, arc melting, yields high-purity samples with complex shapes. In this study, we explore nanoindentation tests at room tempera… ▽ More Concentrated solid solution alloys (CSAs) have attracted attention for their promising properties; however, current manufacturing methods face challenges in complexity, high costs, and limited scalability, raising concerns about industrial viability. The prevalent technique, arc melting, yields high-purity samples with complex shapes. In this study, we explore nanoindentation tests at room temperature where arc-melted samples exhibit larger grain sizes, diminishing the effects of grain boundaries on the results. Motivated by these findings, our investigation focuses on the atomistic-level exploration of plasticity mechanisms, specifically dislocation nucleation and propagation during nanoindentation tests. The intricate chemistry of NiFeCrCo CSA influences pile-ups and slip traces, aiming to elucidate plastic deformation by considering both pristine and pre-existing stacking fault tetrahedra. Our analysis scrutinizes dynamic deformation processes, defect nucleation, and evolution, complemented by stress-strain and dislocation densities-strain curves illustrating the hardening mechanism of defective materials. Additionally, we examine surface morphology and plastic deformation through atomic shear strain and displacement mappings. This integrated approach provides insights into the complex interplay between material structure and mechanical behavior, paving the way for an enhanced understanding and potential advancements in CSA applications. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Simulated multi-component CuZr(Al) metallic glasses akin to experiments

Authors: Rene Alvarez-Donado, Silvia Bonfanti, Mikko Alava

Abstract: We study a three-component CuZrAl metallic glass system by means of a combined Monte Carlo and Molecular Dynamics simulations scheme. This hybrid method allows us to generate equilibrated samples at temperatures below the conventional glass transition for the first time, achieving a more stable glassy regime. By using a realistic potential for the interactions of metallic species, we explore the k… ▽ More We study a three-component CuZrAl metallic glass system by means of a combined Monte Carlo and Molecular Dynamics simulations scheme. This hybrid method allows us to generate equilibrated samples at temperatures below the conventional glass transition for the first time, achieving a more stable glassy regime. By using a realistic potential for the interactions of metallic species, we explore the kinetics, thermodynamics, and rheology of a CuZrAl glass, and then compare these findings with those of the ubiquitous CuZr. Remarkably, the resulting sheared glassy configurations show an abrupt stress drop corresponding to the shear band, akin to experimental observations. Our results pave the way for theoretical studies of complex metallic glasses and offer comparisons with experiments. △ Less

Submitted 9 January, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 5 pages, supplementary information added to the end of the manuscript

Plastic deformation mechanisms during nanoindentation of W, Mo, V body-centered cubic single crystals and their corresponding W-Mo, W-V equiatomic random solid solutions

Authors: F. J. Dominguez-Gutierrez, S. Papanikolaou, S. Bonfanti, M. J. Alava

Abstract: Deformation plasticity mechanisms in alloys and compounds may unveil the material capacity towards optimal mechanical properties. We conduct a series of molecular dynamics (MD) simulations to investigate plasticity mechanisms due to nanoindentation in pure tungsten, molybdenum and vanadium body-centered cubic single crystals, as well as the also body-centered cubic, equiatomic, random solid soluti… ▽ More Deformation plasticity mechanisms in alloys and compounds may unveil the material capacity towards optimal mechanical properties. We conduct a series of molecular dynamics (MD) simulations to investigate plasticity mechanisms due to nanoindentation in pure tungsten, molybdenum and vanadium body-centered cubic single crystals, as well as the also body-centered cubic, equiatomic, random solid solutions (RSS) of tungsten--molybdenum and tungsten--vanadium alloys. Our analysis focuses on a thorough, side-by-side comparison of dynamic deformation processes, defect nucleation, and evolution, along with corresponding stress--strain curves. We also check the surface morphology of indented samples through atomic shear strain mapping. As expected, the presence of Mo and V atoms in W matrices introduces lattice strain and distortion, increasing material resistance to deformation and slowing down dislocation mobility of dislocation loops with a Burgers vector of 1/2 $\langle 111 \rangle$. Our side-by-side comparison displays a remarkable suppression of the plastic zone size in equiatomic W--V RSS, but not in equiatomic W--Mo RSS alloys, displaying a clear prediction for optimal hardening response equiatomic W--V RSS alloys. If the small-depth nanoindentation plastic response is indicative of overall mechanical performance, it is possible to conceive a novel MD-based pathway towards material design for mechanical applications in complex, multi-component alloys. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Striation lines in intermittent fatigue crack growth in an Al alloy

Authors: Anniina Kinnunen, Ivan V. Lomakin, Tero Mäkinen, Kim Widell, Juha Koivisto, Mikko J. Alava

Abstract: Fatigue failure of crystalline materials is a difficult problem in science and engineering, and recent results have shown that fatigue crack growth can occur in intermittent jumps which have fat-tailed distributions. As fatigue crack propagation is known to leave markings -- called striations -- on the fracture surface, the distances between these should also have fat-tailed distributions, if the… ▽ More Fatigue failure of crystalline materials is a difficult problem in science and engineering, and recent results have shown that fatigue crack growth can occur in intermittent jumps which have fat-tailed distributions. As fatigue crack propagation is known to leave markings -- called striations -- on the fracture surface, the distances between these should also have fat-tailed distributions, if the crack propagation is intermittent. Here, we combine macroscale crack tip tracking in fatigue crack growth measurements of aluminum 5005 samples with \emph{post-mortem} scanning electron microscopy imaging of the striation lines. We introduce two different methods for extracting the striation line spacing from the images. What we find is a similar distribution of striation spacings as jump sizes using one of our methods, but the average striation spacing does not correlate with the crack growth rate. We conclude that we observe avalanche-like crack propagation, reflected in both the macroscale crack tip tracking as well as the analysis of the fracture surfaces. Our results show that the fracture surfaces can be used to study the intermittency of fatigue crack propagation and in development of crack-resistant materials. The advantages and disadvantages of the two methods introduced are discussed. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Journal ref: Phys. Rev. Materials 7, 053602, 2023

Failure precursors and failure mechanisms in hierarchically patterned paper sheets in tensile and creep loading

Authors: Mahshid Pournajar, Tero Mäkinen, Seyyed Ahmad Hosseini, Paolo Moretti, Mikko Alava, Michael Zaiser

Abstract: Quasi-brittle materials endowed with (statistically) self-similar hierarcical microstructures show distinct failure patterns that deviate from the standard scenario of damage accumulation followed by crack nucleation-and-growth. Here we study the failure of paper sheets with hierarchical slice patterns as well as non-hierarchical and unpatterned reference samples, considering both uncracked sample… ▽ More Quasi-brittle materials endowed with (statistically) self-similar hierarcical microstructures show distinct failure patterns that deviate from the standard scenario of damage accumulation followed by crack nucleation-and-growth. Here we study the failure of paper sheets with hierarchical slice patterns as well as non-hierarchical and unpatterned reference samples, considering both uncracked samples and samples containing a macroscopic crack. Failure is studied under displacement-controlled tensile loading as well as under creep conditions. Acoustic emission records and surface strain patterns are recorded alongside stress-strain and creep curves. The measurements demonstrate that hierarchical patterning efficiently mitigates against strain localization and crack propagation. In tensile loading, this results in a significantly increased residual strength of cracked samples. Under creep conditions, for a given range of lifetimes hierarchically patterned samples are found to sustain larger creep strains at higher stress levels; their creep curves show unusual behavior characterized by multiple creep rate minima due to the repeated arrest of emergent localization bands. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Quasi-localized modes in crystalline high entropy alloys

Authors: Silvia Bonfanti, Roberto Guerra, Rene Alvarez-Donado, Pawel Sobkowicz, Stefano Zapperi, Mikko Alava

Abstract: High Entropy Alloys (HEAs) are designed by mixing multiple metallic species in nearly the same amount to obtain crystalline or amorphous materials with exceptional mechanical properties. Here we use molecular dynamics simulations to investigate the role of positional and compositional disorder in determining the low-frequency vibrational properties of CrMnFeCoNi HEAs. Our results show that the exp… ▽ More High Entropy Alloys (HEAs) are designed by mixing multiple metallic species in nearly the same amount to obtain crystalline or amorphous materials with exceptional mechanical properties. Here we use molecular dynamics simulations to investigate the role of positional and compositional disorder in determining the low-frequency vibrational properties of CrMnFeCoNi HEAs. Our results show that the expected dependence of the density of states on the frequency as $D(ω)\simω^4$ is recovered for amorphous HEAs and is also observed for partially crystallized alloys with deviations that depend on the degree of crystallization. We find that the quasi-localized vibrations are still visible in crystalline HEAs, albeit suppressed compared to the corresponding amorphous alloys. Our work offers a unified perspective to describe HEA mechanical properties in terms of their vibrational density of states. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 8 pages, 7 figures

Prediction of steel nanohardness by using graph neural networks on surface polycrystallinity maps

Authors: Kamran Karimi, Henri Salmenjoki, Katarzyna Mulewska, Lukasz Kurpaska, Anna Kosińska, Mikko Alava, Stefanos Papanikolaou

Abstract: As a bulk mechanical property, nanoscale hardness in polycrystalline metals is strongly dependent on microstructural features that are believed to be heavily influenced from complex features of polycrystallinity -- namely, individual grain orientations and neighboring grain properties. We train a graph neural network (GNN) model, with each grain center location being a graph node, to assess the pr… ▽ More As a bulk mechanical property, nanoscale hardness in polycrystalline metals is strongly dependent on microstructural features that are believed to be heavily influenced from complex features of polycrystallinity -- namely, individual grain orientations and neighboring grain properties. We train a graph neural network (GNN) model, with each grain center location being a graph node, to assess the predictability of micromechanical responses of nano-indented low-carbon 310S stainless steel (alloyed with Ni and Cr) surfaces, solely based on surface polycrystallinity, captured by electron backscatter diffraction maps. The grain size distribution ranges between $1-100~μ$m, with mean grain size at $18~μ$m. The GNN model is trained on a set of nanomechanical load-displacement curves, obtained from nanoindentation tests and is subsequently used to make predictions of nano-hardness at various depths, with sole input being the grain locations and orientations. Model training is based on a sub-standard set of $\sim10^2$ hardness measurements, leading to an overall satisfactory performance. We explore model performance and its dependence on various structural/topological grain-level descriptors, such as the grain size and number of nearest neighbors. Analogous GNN model frameworks may be utilized for quick, inexpensive hardness estimates, for guidance to detailed nanoindentation experiments, akin to cartography tool developments in the world exploration era. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Detection of the onset of yielding and creep failure from digital image correlation

Authors: Tero Mäkinen, Agata Zaborowska, Małgorzata Frelek-Kozak, Iwona Jóźwik, Łukasz Kurpaska, Stefanos Papanikolaou, Mikko J. Alava

Abstract: There are a multitude of applications in which structural materials would be desired to be nondestructively evaluated, while in a component, for plasticity and failure characteristics. In this way, safety and resilience features can be significantly improved. Nevertheless, while failure can be visible through cracks, plasticity is commonly invisible and highly microstructure-dependent. Here, we sh… ▽ More There are a multitude of applications in which structural materials would be desired to be nondestructively evaluated, while in a component, for plasticity and failure characteristics. In this way, safety and resilience features can be significantly improved. Nevertheless, while failure can be visible through cracks, plasticity is commonly invisible and highly microstructure-dependent. Here, we show that an equation-free method based on principal component analysis is capable of detecting yielding and tertiary creep onset, directly from strain fields that are obtained by digital image correlation, applicable on components, continuously and nondestructively. We demonstrate the applicability of the method to yielding of Ni-based Haynes 230 metal alloy polycrystalline samples, which are also characterized through electron microscopy and benchmarked using continuum polycrystalline plasticity modeling. Also, we successfully apply this method to yielding during uniaxial tension of Hastelloy X polycrystalline samples, and also to the onset of tertiary creep in quasibrittle fiber composites under uniaxial tension. We conclude that there are key features in the spatiotemporal fluctuations of local strain fields that can be used to infer mechanical properties. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Journal ref: Physical Review Materials 6, 103601, 2022

Wood compression in four-dimensional in situ tomography

Authors: Tero Mäkinen, Alisa Halonen, Juha Koivisto, Mikko J. Alava

Abstract: Wood deformation, in particular when subject to compression, exhibits scale-free avalanche-like behavior as well as structure-dependent localization of deformation. We have taken three-dimensional (3D) x-ray tomographs during compression with constant stress rate loading. Using digital volume correlation, we obtain the local total strain during the experiment and compare it to the global strain an… ▽ More Wood deformation, in particular when subject to compression, exhibits scale-free avalanche-like behavior as well as structure-dependent localization of deformation. We have taken three-dimensional (3D) x-ray tomographs during compression with constant stress rate loading. Using digital volume correlation, we obtain the local total strain during the experiment and compare it to the global strain and acoustic emission. The wood cells collapse layer by layer throughout the sample starting from the softest parts, i.e., the spring wood. As the damage progresses, more and more of the softwood layers throughout the sample collapse, which indicates damage spreading instead of localization. In 3D, one can see a fat-tailed local strain rate distribution, indicating that inside the softwood layers, the damage occurs in localized spots. The observed log-normal strain distribution is in agreement with this view of the development of independent local collapses or irreversible deformation events. A key feature in the mechanical behavior of wood is then in the complex interaction of localized deformation between or among the annual rings. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Journal ref: Physical Review Materials 6, L070601, 2022

Dislocation plasticity in equiatomic NiCoCr alloys: The effect of short-range order

Authors: Amir H. Naghdi, Kamran Karimi, Axel E. Poisvert, Amin Esfandiarpour, Rene Alvarez, Pawel Sobkowicz, Mikko Alava, Stefanos Papanikolaou

Abstract: Equiatomic NiCoCr solid solutions have been recently shown to display exceptional mechanical properties, that have been suggested to be linked to nanostructural short-range order (SRO) features that may arise from thermal treatments, such as annealing or/and aging. Here, we use hybrid Monte Carlo-Molecular Dynamics simulations to gain further insights of thermal effects on the SRO formation as wel… ▽ More Equiatomic NiCoCr solid solutions have been recently shown to display exceptional mechanical properties, that have been suggested to be linked to nanostructural short-range order (SRO) features that may arise from thermal treatments, such as annealing or/and aging. Here, we use hybrid Monte Carlo-Molecular Dynamics simulations to gain further insights of thermal effects on the SRO formation as well as the edge dislocation plasticity mechanisms of equiatomic NiCoCr face-centered cubic solid solution. For that purpose, we utilize two well-known NiCoCr interatomic potentials, one of which displays well documented SRO, believed to be linked to experimental evidence and labeled as the Li-Sheng-Ma potential, while the other (Farkas-Caro) does not. We use these two potentials to discern short range ordering (from inherent randomness in random solid solutions) and understand how SROs influence dislocation depinning dynamics in various thermal annealing scenarios. In this context, we used robust, scale-dependent metrics to infer a characteristic SRO size in the Li-Sheng-Ma case by probing local concentration fluctuations which otherwise indicate uncorrelated patterns in the Farkas-Caro case in a close agreement with random alloys. Using relevant order parameters, we also report on the drastic increase of chemical ordering within the stacking fault region. More importantly, we find that the Li-Sheng-Ma potential leads to excellent edge dislocation depinning strength with low stacking fault width. Our findings indicate an enhanced roughening mechanism due to the SROs-misfit synergy that leads to significant improvements in dislocation glide resistance. We argue that the improvements in alloy strength have their atomistic origins in the interplay between nanoscopic SROs and atomic-level misfit properties. △ Less

Submitted 22 October, 2022; originally announced October 2022.

Velocity of interfaces with short and long ranged elasticity under sinusoidal creep

Authors: Juha Savolainen, Mikko Alava

Abstract: Plenty of research on elastic interfaces has been done on systems where the interface is pushed with a constant force. We studied the average velocity of an interface under a sinusoidal driving in the creep region, considering both short-range elastic systems, such as magnetic domain walls during a hysteresis loop, and long ranged systems such as fractures. We obtained a modified version of the cr… ▽ More Plenty of research on elastic interfaces has been done on systems where the interface is pushed with a constant force. We studied the average velocity of an interface under a sinusoidal driving in the creep region, considering both short-range elastic systems, such as magnetic domain walls during a hysteresis loop, and long ranged systems such as fractures. We obtained a modified version of the creep velocity with approximate power-law behaviour and a material dependent exponent for short ranged systems and simpler behaviour for long-range elasticity. We discuss whether the model can be applied to fatigue fractures, or if extra physics is needed. △ Less

Submitted 27 March, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 8 pages, 4 figures

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

Authors: Kamran Karimi, Mikko J. Alava, Stefanos Papanikolaou

Abstract: 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 t… ▽ More 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. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Multiscale nanoindentation modeling of concentrated solid solutions: A continuum plasticity model

Authors: K. Frydrych, F. J. Dominguez-Gutierrez, M. J. Alava, S. Papanikolaou

Abstract: Recently developed single-phase concentrated solid-solution alloys (CSAs) contain multiple elemental species in high concentrations with different elements randomly arranged on a crystalline lattice. These chemically disordered materials present excellent physical properties, including high-temperature thermal stability and hardness, with promising applications to industries at extreme operating e… ▽ More Recently developed single-phase concentrated solid-solution alloys (CSAs) contain multiple elemental species in high concentrations with different elements randomly arranged on a crystalline lattice. These chemically disordered materials present excellent physical properties, including high-temperature thermal stability and hardness, with promising applications to industries at extreme operating environments. The aim of this paper is to present a continuum plasticity model accounting for the first time for the behaviour of a equiatomic five-element CSA, that forms a face-centered cubic lattice. The inherent disorder associated with the lattice distortions caused by an almost equiatomic distribution of atoms, is captured by a single parameter α that quantifies the relative importance of an isotropic plastic contribution to the model. This results in multiple plasticity mechanisms that go beyond crystallographic symmetry-based ones, common in the case of conventional single element metals. We perform molecular dynamics simulations of equiatomic CSAs: NiFe, NiFeCr, NiFeCrCo, and Cantor alloys to validate the proposed continuum model which is implemented in the finite element method and applied to model nanoindentation tests for three different crystallographic orientations. We obtain the representative volume element model by tracking the combined model yield surface. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Effect of thresholding on avalanches and their clustering for interfaces with long-range elasticity

Authors: Juha Savolainen, Lasse Laurson, Mikko Alava

Abstract: Avalanches are often defined as signals higher than some detection level in bursty systems. The choice of the detection threshold affects the number of avalanches, but it can also affect their temporal correlations. We simulated the depinning of a long-range elastic interface and applied different thresholds including a zero one on the data to see how the sizes and durations of events change and h… ▽ More Avalanches are often defined as signals higher than some detection level in bursty systems. The choice of the detection threshold affects the number of avalanches, but it can also affect their temporal correlations. We simulated the depinning of a long-range elastic interface and applied different thresholds including a zero one on the data to see how the sizes and durations of events change and how this affects temporal avalanche clustering. Higher thresholds result in steeper size and duration distributions and cause the avalanches to cluster temporally. Using methods from seismology, the frequency of the events in the clusters was found to decrease as a power-law of time, and the size of an event in a cluster was found to help predict how many events it is followed by. The results bring closer theoretical studies of this class of models to real experiments, but also highlight how different phenomena can be obtained from the same set of data. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 8 pages, 8 figures

Journal ref: Phys. Rev. E, 105, 054152, 2022

Atomistic simulations of nanoindentation in single crystalline tungsten: The role of interatomic potentials

Authors: F. J. Dominguez-Gutierrez, P. Grigorev, A. Naghdi, Q. Q. Xu, J. Byggmastar, G. Y. Wei, T. D. Swinburne, S. Papanikolaou, M. J. Alava

Abstract: Computational modeling is usually applied to aid experimental exploration of advanced materials to better understand the fundamental plasticity mechanisms during mechanical testing. In this work, we perform Molecular dynamics (MD) simulations to emulate experimental room temperature spherical-nanoindentation of crystalline W matrices by different interatomic potentials: EAM, modified EAM, and a re… ▽ More Computational modeling is usually applied to aid experimental exploration of advanced materials to better understand the fundamental plasticity mechanisms during mechanical testing. In this work, we perform Molecular dynamics (MD) simulations to emulate experimental room temperature spherical-nanoindentation of crystalline W matrices by different interatomic potentials: EAM, modified EAM, and a recently developed machine learned based tabulated Gaussian approximation potential (tabGAP) for describing the interaction of W-W. Results show similarities between load displacements and stress-strain curves, regardless of the numerical model. However, a discrepancy is observed at early stages of the elastic to plastic deformation transition showing different mechanisms for dislocation nucleation and evolution, that is attributed to the difference of Burgers vector magnitudes, stacking fault and dislocation glide energies. Besides, contact pressure is investigated by considering large indenters sizes that provides a detailed analysis of screw and edge dislocations during loading process. Furthermore, the glide barrier of this kind of dislocations are reported for all the interatomic potentials showing that tabGAP model presents the most accurate results with respect to density functional theory calculations and a good qualitative agreement with reported experimental data △ Less

Submitted 11 December, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Edge dislocations in multi-component solid solution alloys: Beyond traditional elastic depinning

Authors: A. Esfandiarpour, S. Papanikolaou, M. Alava

Abstract: High-entropy alloys (HEA) form solid solutions with large chemical disorder and excellent mechanical properties. We investigate the origin of HEA strengthening in face-centered cubic (FCC) single-phase HEAs through molecular dynamics simulations of dislocations, in particular, the equiatomic $\rm CrCoNi$, $\rm CrMnCoNi$, $\rm CrFeCoNi$, $\rm CrMnFeCoNi$, $\rm FeNi$, and also,… ▽ More High-entropy alloys (HEA) form solid solutions with large chemical disorder and excellent mechanical properties. We investigate the origin of HEA strengthening in face-centered cubic (FCC) single-phase HEAs through molecular dynamics simulations of dislocations, in particular, the equiatomic $\rm CrCoNi$, $\rm CrMnCoNi$, $\rm CrFeCoNi$, $\rm CrMnFeCoNi$, $\rm FeNi$, and also, $\rm Fe_{0.4}Mn_{0.27}Ni_{0.26}Co_{0.05}Cr_{0.02}$, $\rm Fe_{0.7}Ni_{0.11}Cr_{0.19}$. The dislocation correlation length $ξ$, roughness amplitude $R_{a}$, and stacking fault widths $W_{SF}$ are tracked as a function of stress. All alloys are characterized by a well defined depinning stress ($σ_c$) and we find a novel regime where exceptional strength is observed, and a fortuitous combination takes place, of small stacking fault widths and large dislocation roughness $R_{a}$. Thus the depinning of two partials seems analogous to unconventional domain wall depinning in disordered magnetic thin films. This novel regime is identified in specific compositions commonly associated with exceptional mechanical properties ($\rm CrCoNi$, $\rm CrMnCoNi$, $\rm CrFeCoNi$, and $\rm CrMnFeCoNi$). Yield stress from analytical solute-strengthening models underestimates largely the results in these cases. A possible strategy for increasing strength in multi-component single-phase alloys is the combined design of stacking fault width and element-based chemical disorder. △ Less

Submitted 31 October, 2021; originally announced November 2021.

Comments: 5 pages, 5 figures, it was submitted to Physical review letter journal (submission date :10/5/2021)

Direct detection of plasticity onset through total-strain profile evolution

Authors: Stefanos Papanikolaou, Mikko J. Alava

Abstract: Plastic yielding in solids strongly depends on various conditions, such as temperature and loading rate and indeed, sample-dependent knowledge of yield points in structural materials promotes reliability in mechanical behavior. Commonly, yielding is measured through controlled mechanical testing at small or large scales, in ways that either distinguish elastic (stress) from total deformation measu… ▽ More Plastic yielding in solids strongly depends on various conditions, such as temperature and loading rate and indeed, sample-dependent knowledge of yield points in structural materials promotes reliability in mechanical behavior. Commonly, yielding is measured through controlled mechanical testing at small or large scales, in ways that either distinguish elastic (stress) from total deformation measurements, or by identifying plastic slip contributions. In this paper we argue that instead of separate elastic/plastic measurements, yielding can be unraveled through statistical analysis of total strain fluctuations during the evolution sequence of profiles measured in-situ, through digital image correlation. We demonstrate two distinct ways of precisely quantifying yield locations in widely applicable crystal plasticity models, that apply in polycrystalline solids, either by using principal component analysis or discrete wavelet transforms. We test and compare these approaches in synthetic data of polycrystal simulations and a variety of yielding responses, through changes of the applied loading rates and the strain-rate sensitivity exponents. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: 9 pages, 6 figures

Journal ref: Phys. Rev. Materials 5, 083602 (2021)

Avalanche correlations and stress-strain curves in discrete dislocation plasticity

Authors: Henri Salmenjoki, Lasse Laurson, Mikko J. Alava

Abstract: The sequence of deformation bursts during plastic deformation exhibits scale-free features. In addition to the burst or avalanche sizes and the rate of avalanches the process is characterized by correlations in the series which become manifest in the resulting shape of the stress-strain curve. We analyze such features of plastic deformation with 2D and 3D simulations of discrete dislocation dynami… ▽ More The sequence of deformation bursts during plastic deformation exhibits scale-free features. In addition to the burst or avalanche sizes and the rate of avalanches the process is characterized by correlations in the series which become manifest in the resulting shape of the stress-strain curve. We analyze such features of plastic deformation with 2D and 3D simulations of discrete dislocation dynamics models and we show, that only with severe plastic deformation the ensuing memory effects become negligible. The role of past deformation history and dislocation pinning by disorder are studied. In general, the correlations have the effect of reducing the scatter of the individual stress-strain curves around the mean one. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 9 pages, 9 figures. Submitted

Journal ref: Phys. Rev. Materials 5, 073601 (2021)

Nanoindentation of single crystalline Mo: Atomistic defect nucleation and thermomechanical stability

Authors: F. J. Domínguez-Gutiérrez, S. Papanikolaou, A. Esfandiarpour, P. Sobkowicz, M. Alava

Abstract: The mechanical responses of single crystalline Body-Centered Cubic (BCC) metals, such as molybdenum (Mo), outperform other metals at high temperatures, so much so that they are considered as excellent candidates for applications under extreme conditions, such as the divertor of fusion reactors. The excellent thermomechanical stability of molybdenum at high temperatures (400-1000$^{\rm o}$C) has al… ▽ More The mechanical responses of single crystalline Body-Centered Cubic (BCC) metals, such as molybdenum (Mo), outperform other metals at high temperatures, so much so that they are considered as excellent candidates for applications under extreme conditions, such as the divertor of fusion reactors. The excellent thermomechanical stability of molybdenum at high temperatures (400-1000$^{\rm o}$C) has also been detected through nanoindentation, pointing towards connections to emergent local dislocation mechanisms related to defect nucleation. In this work, we carry out a computational study of the effects of high temperature on the mechanical deformation properties of single crystalline Mo under nanoindentation. Molecular dynamics (MD) simulations of spherical nanoindentation are performed at two indenter tip diameters and crystalline sample orientations [100], [110], and [111], for the temperature range of 10-1000K. We investigate how the increase of temperature influences the nanoindentation process, modifying dislocation densities, mechanisms, atomic displacements and also, hardness, in agreement with reported experimental measurements. Our results suggest that the characteristic formation and high-temperature stability of [001] dislocation junctions in Mo during nanoindentation, in contrast to other BCC metals, may be the cause of the persistent thermomechanical stability of Mo. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Propagating bands of plastic deformation in a metal alloy as critical avalanches

Authors: Tero Mäkinen, Pasi Karppinen, Markus Ovaska, Lasse Laurson, Mikko J. Alava

Abstract: The plastic deformation of metal alloys localizes in the Portevin--Le Chatelier effect in bands of different types, including propagating, or type 'A' bands, usually characterised by their width and a typical propagation velocity. This plastic instability arises from collective dynamics of dislocations interacting with mobile solute atoms but the resulting sensitivity to the strain-rate lacks fund… ▽ More The plastic deformation of metal alloys localizes in the Portevin--Le Chatelier effect in bands of different types, including propagating, or type 'A' bands, usually characterised by their width and a typical propagation velocity. This plastic instability arises from collective dynamics of dislocations interacting with mobile solute atoms but the resulting sensitivity to the strain-rate lacks fundamental understanding. Here we show, by employing high-resolution imaging in tensile deformation experiments of an aluminum alloy, that the band velocities exhibit large fluctuations. Each band produces a velocity signal reminiscent of crackling noise bursts observed in numerous driven avalanching systems from propagating cracks in fracture to the Barkhausen effect in ferromagnets. The statistical features of these velocity bursts including their average shapes and size distributions obey predictions of a simple mean-field model of critical avalanche dynamics. Our results thus reveal a novel paradigm of criticality in the localization of deformation. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Journal ref: Science Advances 6(41), eabc7350 (2020)

Probing the transition from dislocation jamming to pinning by machine learning

Authors: Henri Salmenjoki, Lasse Laurson, Mikko J. Alava

Abstract: Collective motion of dislocations is governed by the obstacles they encounter. In pure crystals, dislocations form complex structures as they become jammed by their anisotropic shear stress fields. On the other hand, introducing disorder to the crystal causes dislocations to pin to these impeding elements and, thus, leads to a competition between dislocation-dislocation and dislocation-disorder in… ▽ More Collective motion of dislocations is governed by the obstacles they encounter. In pure crystals, dislocations form complex structures as they become jammed by their anisotropic shear stress fields. On the other hand, introducing disorder to the crystal causes dislocations to pin to these impeding elements and, thus, leads to a competition between dislocation-dislocation and dislocation-disorder interactions. Previous studies have shown that, depending on the dominating interaction, the mechanical response and the way the crystal yields change. Here we employ three-dimensional discrete dislocation dynamics simulations with varying density of fully coherent precipitates to study this phase transition $-$ from jamming to pinning $-$ using unsupervised machine learning. By constructing descriptors characterizing the evolving dislocation configurations during constant loading, a confusion algorithm is shown to be able to distinguish the systems into two separate phases. These phases agree well with the observed changes in the relaxation rate during the loading. Our results also give insights on the structure of the dislocation networks in the two phases. △ Less

Submitted 28 September, 2020; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: 16 pages, 9 figures. Accepted for publication in Materials Theory

Plastic yielding and deformation bursts in the presence of disorder from coherent precipitates

Authors: Henri Salmenjoki, Arttu Lehtinen, Lasse Laurson, Mikko J. Alava

Abstract: Alloying metals with other elements is often done to improve the material strength or hardness. A key microscopic mechanism is precipitation hardening, where precipitates impede dislocation motion, but the role of such obstacles in determining the nature of collective dislocation dynamics remains to be understood. Here, three-dimensional discrete dislocation dynamics simulations of FCC single crys… ▽ More Alloying metals with other elements is often done to improve the material strength or hardness. A key microscopic mechanism is precipitation hardening, where precipitates impede dislocation motion, but the role of such obstacles in determining the nature of collective dislocation dynamics remains to be understood. Here, three-dimensional discrete dislocation dynamics simulations of FCC single crystals are performed with fully coherent spherical precipitates from zero precipitate density upto $ρ_p = 10^{21}\,\text{m}^{-3}$ and at various dislocation-precipitate interaction strengths. When the dislocation-precipitate interactions do not play a major role the yielding is qualitatively as for pure crystals, i.e., dominated by "dislocation jamming", resulting in glassy dislocation dynamics exhibiting critical features at any stress value. We demonstrate that increasing the precipitate density and/or the dislocation-precipitate interaction strength creates a true yield or dislocation assembly depinning transition, with a critical yield stress. This is clearly visible in the statistics of dislocation avalanches observed when quasistatically ramping up the external stress, and is also manifested in the response of the system to constant applied stresses. The scaling of the yielding with precipitates is discussed in terms of the Bacon-Kocks-Scattergood relation. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Comments: 7 pages, 5 figures. Submitted to Physical Review Materials

Journal ref: Phys. Rev. Materials 4, 083602 (2020)

Machine learning and predicting the time dependent dynamics of local yielding in dry foams

Authors: Leevi Viitanen, Jonatan R. Mac Intyre, Juha Koivisto, Antti Puisto, Mikko Alava

Abstract: The yielding of dry foams is enabled by small elementary yield events on the bubble scale, "T1"s. We study the large scale detection of these in an expanding 2D flow geometry using artificial intelligence (AI) and nearest neighbour analysis. A good level of accuracy is reached by the AI approach using only a single frame, with the maximum score for vertex centered images highlighting the important… ▽ More The yielding of dry foams is enabled by small elementary yield events on the bubble scale, "T1"s. We study the large scale detection of these in an expanding 2D flow geometry using artificial intelligence (AI) and nearest neighbour analysis. A good level of accuracy is reached by the AI approach using only a single frame, with the maximum score for vertex centered images highlighting the important role the vertices play in the local yielding of foams. We study the predictability of T1s ahead of time and show that this is possible on a timescale related to the waiting time statistics of T1s in local neighborhoods. The local T1 event predictability development is asymmetric in time, and measures the variation of the local property to yielding and similarly the existence of a relaxation timescale post local yielding. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Modelling aerosol transport and virus exposure with numerical simulations in relation to SARS-CoV-2 transmission by inhalation indoors

Authors: Ville Vuorinen, Mia Aarnio, Mikko Alava, Ville Alopaeus, Nina Atanasova, Mikko Auvinen, Nallannan Balasubramanian, Hadi Bordbar, Panu Erästö, Rafael Grande, Nick Hayward, Antti Hellsten, Simo Hostikka, Jyrki Hokkanen, Ossi Kaario, Aku Karvinen, Ilkka Kivistö, Marko Korhonen, Risto Kosonen, Janne Kuusela, Sami Lestinen, Erkki Laurila, Heikki Nieminen, Petteri Peltonen, Juho Pokki , et al. (5 additional authors not shown)

Abstract: We provide research findings on the physics of aerosol dispersion relevant to the hypothesized aerosol transmission of SARS-CoV-2. We utilize physics-based modeling at different levels of complexity, and literature on coronaviruses, to investigate the possibility of airborne transmission. The previous literature, our 0D-3D simulations by various physics-based models, and theoretical calculations,… ▽ More We provide research findings on the physics of aerosol dispersion relevant to the hypothesized aerosol transmission of SARS-CoV-2. We utilize physics-based modeling at different levels of complexity, and literature on coronaviruses, to investigate the possibility of airborne transmission. The previous literature, our 0D-3D simulations by various physics-based models, and theoretical calculations, indicate that the typical size range of speech and cough originated droplets (d < 20microns) allows lingering in the air for O(1h) so that they could be inhaled. Consistent with the previous literature, numerical evidence on the rapid drying process of even large droplets, up to sizes O(100microns), into droplet nuclei/aerosols is provided. Based on the literature and the public media sources, we provide evidence that the infected individuals could have been exposed to aerosols/droplet nuclei by inhaling them in significant numbers e.g. O(100). By 3D computational fluid dynamics (CFD) simulations, we give examples on the transport and dilution of aerosols (d<20microns) over distances O(10m) in generic environments. We study susceptible and infected individuals in generic public places by Monte-Carlo modeling. The model accounts for the locally varying aerosol concentration levels which the susceptible accumulate via inhalation. The introduced concept, 'exposure time' to virus containing aerosols is proposed to complement the traditional 'safety distance' thinking. We show that the exposure time to inhale O(100) aerosols could range from O(1s) to O(1min) or even to O(1h) depending on the situation. The Monte Carlo analysis provides clear quantitative insight to the exposure time in different public indoor environments. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: Preprint, 58 pages, 15 figures

Machine learning plastic deformation of crystals

Authors: Henri Salmenjoki, Mikko J. Alava, Lasse Laurson

Abstract: Plastic deformation of micron-scale crystalline solids exhibits stress-strain curves with significant sample-to-sample variations. It is a pertinent question if this variability is purely random or to some extent predictable. Here we show, by employing machine learning techniques such as regression neural networks and support vector machines that deformation predictability evolves with strain and… ▽ More Plastic deformation of micron-scale crystalline solids exhibits stress-strain curves with significant sample-to-sample variations. It is a pertinent question if this variability is purely random or to some extent predictable. Here we show, by employing machine learning techniques such as regression neural networks and support vector machines that deformation predictability evolves with strain and crystal size. Using data from discrete dislocations dynamics simulations, the machine learning models are trained to infer the mapping from features of the pre-existing dislocation configuration to the stress-strain curves. The predictability vs strain relation is non-monotonic and exhibits a system size effect: larger systems are more predictable. Stochastic deformation avalanches give rise to fundamental limits of deformation predictability for intermediate strains. However, the large-strain deformation dynamics of the samples can be predicted surprisingly well. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: Published version from https://www.nature.com/articles/s41467-018-07737-2. Supplementary material can also be found from the link

Journal ref: Nat Commun 9, 5307 (2018)

The effect of pressure pulsing on the mechanical dewatering of nanofiber suspensions

Authors: Marko Korhonen, Antti Puisto, Mikko Alava, Thaddeus Maloney

Abstract: Dewatering processes are invariably encountered in the chemical manufacturing and processing of various bioproducts. In this study, Computational Fluid Mechanics (CFD) simulations and theory are utilized to model and optimize the dewatering of commercial nanofiber suspensions. The CFD simulations are based on the volume-averaged Navier-Stokes equations while the analytical model is deduced from th… ▽ More Dewatering processes are invariably encountered in the chemical manufacturing and processing of various bioproducts. In this study, Computational Fluid Mechanics (CFD) simulations and theory are utilized to model and optimize the dewatering of commercial nanofiber suspensions. The CFD simulations are based on the volume-averaged Navier-Stokes equations while the analytical model is deduced from the empirical Darcy's law for dewatering flows. The results are then successfully compared to experimental data on commercial cellulose suspensions obtained with a Dynamic Drainage Analyzer (DDA). Both the CFD simulations and the analytical model capture the dewatering flow profiles of the commercial suspensions in an experiment utilizing a constant pressure profile. However, a temporally varying pressure profile offers a superior dewatering performance, as indicated by both the simulations and the analytical model. Finally, the analytical model also predicts an optimized number of pressure pulses, minimizing the time required to completely dewater the suspension. △ Less

Submitted 30 October, 2019; v1 submitted 29 October, 2019; originally announced October 2019.

Entropy-driven impurity-induced nematic-isotropic transition of liquid crystals

Authors: Pritam Kumar Jana, Julien Lam, Nagma Parveen, Mikko J Alava, Lasse Laurson

Abstract: Phase behavior of liquid crystals is of long-standing interest due to numerous applications, with one of the key issues being how the presence of impurities affects the liquid crystalline order. Here we study the orientational order of 4-cyano-4$^{'}$-pentylbiphenyl (5CB) and 4-cyano-4$^{'}$-hexylbiphenyl (6CB) nematic liquid crystals in the presence of varying concentrations of water and n-hexane… ▽ More Phase behavior of liquid crystals is of long-standing interest due to numerous applications, with one of the key issues being how the presence of impurities affects the liquid crystalline order. Here we study the orientational order of 4-cyano-4$^{'}$-pentylbiphenyl (5CB) and 4-cyano-4$^{'}$-hexylbiphenyl (6CB) nematic liquid crystals in the presence of varying concentrations of water and n-hexane molecules serving as impurities, by carrying out both fully atomistic simulations and experiments. Our results reveal that mixing of the impurities (in case of hexane) with the host liquid crystals causes a nematic-to-isotropic phase transition with hexane concentration as the control parameter while demixing (in case of water) results in only weak impurity-induced perturbations to the nematic liquid-crystalline order. We develop a coarse-grained model illustrating the general nature and entropic origin of the mixing-induced phase transition. △ Less

Submitted 12 October, 2019; originally announced October 2019.

Nanoscale Liquid Crystal Lubrication Controlled by Surface Structure and Film Composition

Authors: Pritam Kumar Jana, Wei Chen, Mikko J. Alava, Lasse Laurson

Abstract: Liquid crystals have emerged as potential candidates for next-generation lubricants due to their tendency to exhibit long-range ordering. Here, we construct a full atomistic model of 4-cyano-4-hexylbiphenyl (6CB) nematic liquid crystal lubricants mixed with hexane and confined by mica surfaces. We explore the effect of the surface structure of mica, as well as lubricant composition and thickness,… ▽ More Liquid crystals have emerged as potential candidates for next-generation lubricants due to their tendency to exhibit long-range ordering. Here, we construct a full atomistic model of 4-cyano-4-hexylbiphenyl (6CB) nematic liquid crystal lubricants mixed with hexane and confined by mica surfaces. We explore the effect of the surface structure of mica, as well as lubricant composition and thickness, on the nanoscale friction in the system. Our results demonstrate the key role of the structure of the mica surfaces, specifically the positions of potassium ($\mathrm{K}^+$) ions, in determining the nature of sliding friction with monolayer lubricants, including the presence or absence of stick-slip dynamics. With the commensurate setup of confining surfaces, when the grooves created between the periodic $\mathrm{K}^+$ ions are parallel to the sliding direction we observe a lower friction force as compared to the perpendicular situation. Random positions of ions exhibit even smaller friction forces with respect to the previous two cases. For thicker lubrication layers the surface structure becomes less important and we observe a good agreement with the experimental data on bulk viscosity of 6CB and the additive hexane. In case of thicker lubrication layers, friction may still be controlled by tuning the relative concentrations of 6CB and hexane in the mixture. △ Less

Submitted 23 December, 2018; originally announced December 2018.

Journal ref: Phys. Chem. Chem. Phys., 2018,20, 18737-18743

Mimicking complex dislocation dynamics by interaction networks

Authors: Henri Salmenjoki, Mikko J. Alava, Lasse Laurson

Abstract: Two-dimensional discrete dislocation models exhibit complex dynamics in relaxation and under external loading. This is manifested both in the time-dependent velocities of individual dislocations and in the ensemble response, the strain rate. Here we study how well this complexity may be reproduced using so-called Interaction Networks, an Artificial Intelligence method for learning the dynamics of… ▽ More Two-dimensional discrete dislocation models exhibit complex dynamics in relaxation and under external loading. This is manifested both in the time-dependent velocities of individual dislocations and in the ensemble response, the strain rate. Here we study how well this complexity may be reproduced using so-called Interaction Networks, an Artificial Intelligence method for learning the dynamics of complex interacting systems. We test how to learn such networks using creep data, and show results on reproducing individual and collective dislocation velocities. The quality of reproducing the interaction kernel is discussed. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Journal ref: Eur. Phys. J. B, 91 11 (2018) 275

Excitation spectra in crystal plasticity

Authors: Markus Ovaska, Arttu Lehtinen, Mikko J. Alava, Lasse Laurson, Stefano Zapperi

Abstract: Plastically deforming crystals exhibit scale-free fluctuations that are similar to those observed in driven disordered elastic systems close to depinning, but the nature of the yielding critical point is still debated. Here, we study the marginal stability of ensembles of dislocations and compute their excitation spectrum in two and three dimensions. Our results show the presence of a singularity… ▽ More Plastically deforming crystals exhibit scale-free fluctuations that are similar to those observed in driven disordered elastic systems close to depinning, but the nature of the yielding critical point is still debated. Here, we study the marginal stability of ensembles of dislocations and compute their excitation spectrum in two and three dimensions. Our results show the presence of a singularity in the distribution of {\it excitation stresses}, i.e., the stress needed to make a localized region unstable, that is remarkably similar to the one measured in amorphous plasticity and spin glasses. These results allow us to understand recent observations of extended criticality in bursty crystal plasticity and explain how they originate from the presence of a pseudogap in the excitation spectrum. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: 5 pages, 4 figures, to appear in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 119, 265501 (2017)

Mechanics of disordered auxetic metamaterials

Authors: Maryam Hanifpour, Charlotte F. Petersen, Mikko J. Alava, Stefano Zapperi

Abstract: Auxetic materials are of great engineering interest not only because of their fascinating negative Poisson's ratio, but also due to their increased toughness and indentation resistance. These materials are typically synthesized polyester foams with a very heterogeneous structure, but the role of disorder in auxetic behavior is not fully understood. Here, we provide a systematic theoretical and exp… ▽ More Auxetic materials are of great engineering interest not only because of their fascinating negative Poisson's ratio, but also due to their increased toughness and indentation resistance. These materials are typically synthesized polyester foams with a very heterogeneous structure, but the role of disorder in auxetic behavior is not fully understood. Here, we provide a systematic theoretical and experimental investigation in to the effect of disorder on the mechanical properties of a paradigmatic auxetic lattice with a re-entrant hexagonal geometry. We show that disorder has a marginal effect on the Poisson's ratio unless the lattice topology is altered, and in all cases examined the disorder preserves the auxetic characteristics. Depending on the direction of loading applied to these disordered auxetic lattices, either brittle or ductile failure is observed. It is found that brittle failure is associated with a disorder-dependent tensile strength, whereas in ductile failure disorder does not affect strength. Our work thus provides general guidelines to optimize elasticity and strength of disordered auxetic metamaterials. △ Less

Submitted 11 October, 2017; v1 submitted 4 April, 2017; originally announced April 2017.

Comments: 5 pages, 5 figures

Friction controls submerged granular flows

Authors: Juha Koivisto, Marko Korhonen, Mikko J. Alava, Carlos P. Ortiz, Douglas J. Durian, Antti Puisto

Abstract: We investigate the coupling between interstitial medium and granular particles by studying the hopper flow of dry and submerged system experimentally and numerically. In accordance with earlier studies, we find, that the dry hopper empties at a constant rate. However, in the submerged system we observe the surging of the flow rate. We model both systems using the discrete element method, which we… ▽ More We investigate the coupling between interstitial medium and granular particles by studying the hopper flow of dry and submerged system experimentally and numerically. In accordance with earlier studies, we find, that the dry hopper empties at a constant rate. However, in the submerged system we observe the surging of the flow rate. We model both systems using the discrete element method, which we couple with computational fluid dynamics in the case of a submerged hopper. We are able to match the simulations and the experiments with good accuracy. To do that, we fit the particle-particle contact friction for each system separately, finding that submerging the hopper changes the particle-particle contact friction from $μ_{vacuum}=0.15$ to $μ_{sub}=0.13$, while all the other simulation parameters remain the same. Furthermore, our experiments find a particle size dependence to the flow rate, which is comprehended based on arguments on the terminal velocity and drag. These results jointly allow us to conclude that at the large particle limit, the interstitial medium does not matter, in contrast to small particles. The particle size limit, where this occurs depends on the viscosity of the interstitial fluid. △ Less

Submitted 7 April, 2017; v1 submitted 28 March, 2017; originally announced March 2017.

Comments: 8 pages, 8 figures, 2 Supplementary videos at https://www.youtube.com/playlist?list=PLlB0dcWeUNfvpPvM6nAU1KsGv_4wbMKPw

Start-up inertia as an origin for heterogeneous flow

Authors: Marko Korhonen, Mikael Mohtaschemi, Antti Puisto, Xavier Illa, Mikko J. Alava

Abstract: For quite some time non-monotonic flow curve was thought to be a requirement for shear banded flows in complex fluids. Thus, in simple yield stress fluids shear banding was considered to be absent. Recent spatially resolved rheological experiments have found simple yield stress fluids to exhibit shear banded flow profiles. One proposed mechanism for the initiation of such transient shear banding p… ▽ More For quite some time non-monotonic flow curve was thought to be a requirement for shear banded flows in complex fluids. Thus, in simple yield stress fluids shear banding was considered to be absent. Recent spatially resolved rheological experiments have found simple yield stress fluids to exhibit shear banded flow profiles. One proposed mechanism for the initiation of such transient shear banding process has been a small stress heterogeneity rising from the experimental device geometry. Here, using Computational Fluid Dynamics methods, we show that transient shear banding can be initialized even under homogeneous stress conditions by the fluid start-up inertia, and that such mechanism indeed is present in realistic experimental conditions. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Journal ref: Phys. Rev. E 95, 022608 (2017)

Inter-event correlations from avalanches hiding below the detection threshold

Authors: Sanja Janićević, Lasse Laurson, Knut Jørgen Måløy, Stéphane Santucci, Mikko J. Alava

Abstract: Numerous systems ranging from deformation of materials to earthquakes exhibit bursty dynamics, which consist of a sequence of events with a broad event size distribution. Very often these events are observed to be temporally correlated or clustered, evidenced by power-law distributed waiting times separating two consecutive activity bursts. We show how such inter-event correlations arise simply be… ▽ More Numerous systems ranging from deformation of materials to earthquakes exhibit bursty dynamics, which consist of a sequence of events with a broad event size distribution. Very often these events are observed to be temporally correlated or clustered, evidenced by power-law distributed waiting times separating two consecutive activity bursts. We show how such inter-event correlations arise simply because of a finite detection threshold, created by the limited sensitivity of the measurement apparatus, or used to subtract background activity or noise from the activity signal. Data from crack propagation experiments and numerical simulations of a non-equilibrium crack line model demonstrate how thresholding leads to correlated bursts of activity by separating the avalanche events into sub-avalanches. The resulting temporal sub-avalanche correlations are well-described by our general scaling description of thresholding-induced correlations in crackling noise. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: 6 pages, 4 figures, Supplemental Material as an ancillary file, accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 117, 230601 (2016)

Bursts of activity in collective cell migration

Authors: Oleksandr Chepizhko, Costanza Giampietro, Eleonora Mastrapasqua, Mehdi Nourazar, Miriam Ascagni, Michela Sugni, Umberto Fascio, Livio Leggio, Chiara Malinverno, Giorgio Scita, Stephane Santucci, Mikko J. Alava, Stefano Zapperi, Caterina A. M. La Porta

Abstract: Dense monolayers of living cells display intriguing relaxation dynamics, reminiscent of soft and glassy materials close to the jamming transition, and migrate collectively when space is available, as in wound healing or in cancer invasion. Here we show that collective cell migration occurs in bursts that are similar to those recorded in the propagation of cracks, fluid fronts in porous media and f… ▽ More Dense monolayers of living cells display intriguing relaxation dynamics, reminiscent of soft and glassy materials close to the jamming transition, and migrate collectively when space is available, as in wound healing or in cancer invasion. Here we show that collective cell migration occurs in bursts that are similar to those recorded in the propagation of cracks, fluid fronts in porous media and ferromagnetic domain walls. In analogy with these systems, the distribution of activity bursts displays scaling laws that are universal in different cell types and for cells moving on different substrates. The main features of the invasion dynamics are quantitatively captured by a model of interacting active particles moving in a disordered landscape. Our results illustrate that collective motion of living cells is analogous to the corresponding dynamics in driven, but inanimate, systems. △ Less

Submitted 25 October, 2016; originally announced October 2016.

Journal ref: PNAS 113 (41), 11408-11413 (2016)

Glassy features of crystal plasticity

Authors: Arttu Lehtinen, Giulio Costantini, Mikko J. Alava, Stefano Zapperi, Lasse Laurson

Abstract: Crystal plasticity occurs by deformation bursts due to the avalanche-like motion of dislocations. Here we perform extensive numerical simulations of a three-dimensional dislocation dynamics model under quasistatic stress-controlled loading. Our results show that avalanches are power-law distributed, and display peculiar stress and sample size dependence: The average avalanche size grows exponentia… ▽ More Crystal plasticity occurs by deformation bursts due to the avalanche-like motion of dislocations. Here we perform extensive numerical simulations of a three-dimensional dislocation dynamics model under quasistatic stress-controlled loading. Our results show that avalanches are power-law distributed, and display peculiar stress and sample size dependence: The average avalanche size grows exponentially with the applied stress, and the amount of slip increases with the system size. These results suggest that intermittent deformation processes in crystalline materials exhibit an extended critical-like phase in analogy to glassy systems, instead of originating from a non-equilibrium phase transition critical point. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 6 pages, 4 figures, Supplemental Material as an ancillary file, accepted for publication in Phys. Rev. B

Multi-scale modeling of dislocation-precipitate interactions in Fe: from molecular dynamics to discrete dislocations

Authors: Arttu Lehtinen, Fredric Granberg, Lasse Laurson, Kai Nordlund, Mikko J. Alava

Abstract: The stress-driven motion of dislocations in crystalline solids, and thus the ensuing plastic deformation process, is greatly influenced by the presence or absence of various point-like defects such as precipitates or solute atoms. These defects act as obstacles for dislocation motion and hence affect the mechanical properties of the material. Here we combine molecular dynamics studies with three-d… ▽ More The stress-driven motion of dislocations in crystalline solids, and thus the ensuing plastic deformation process, is greatly influenced by the presence or absence of various point-like defects such as precipitates or solute atoms. These defects act as obstacles for dislocation motion and hence affect the mechanical properties of the material. Here we combine molecular dynamics studies with three-dimensional discrete dislocation dynamics simulations in order to model the interaction between different kinds of precipitates and a $\frac{1}{2}\langle 1 1 1\rangle$ $\{1 1 0\}$ edge dislocation in BCC iron. We have implemented immobile spherical precipitates into the ParaDis discrete dislocation dynamics code, with the dislocations interacting with the precipitates via a Gaussian potential, generating a normal force acting on the dislocation segments. The parameters used in the discrete dislocation dynamics simulations for the precipitate potential, the dislocation mobility, shear modulus and dislocation core energy are obtained from molecular dynamics simulations. We compare the critical stresses needed to unpin the dislocation from the precipitate in molecular dynamics and discrete dislocation dynamics simulations in order to fit the two methods together, and discuss the variety of the relevant pinning/depinning mechanisms. △ Less

Submitted 9 February, 2016; originally announced February 2016.

Journal ref: Phys.Rev.E 93, 013309, 2016

Avalanches in wood compression

Authors: Tero Mäkinen, Amandine Miksic, Markus Ovaska, Mikko J. Alava

Abstract: Wood is a multi-scale material exhibiting a complex viscoplastic response. We study avalanches in small wood samples in compression. "Woodquakes" measured by acoustic emission are surprisingly similar to earthquakes and crackling noise in rocks and laboratory tests on brittle materials. Both the distributions of event energies and of waiting (silent) times follow power-laws. The stress- strain res… ▽ More Wood is a multi-scale material exhibiting a complex viscoplastic response. We study avalanches in small wood samples in compression. "Woodquakes" measured by acoustic emission are surprisingly similar to earthquakes and crackling noise in rocks and laboratory tests on brittle materials. Both the distributions of event energies and of waiting (silent) times follow power-laws. The stress- strain response exhibits clear signatures of localization of deformation to "weak spots" or softwood layers, as identified using Digital Image Correlation. Even though material structure-dependent localization takes place, the avalanche behavior remains scale-free. △ Less

Submitted 15 June, 2015; originally announced June 2015.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 115, 055501 (2015)

Stick-Slip Control in Nanoscale Boundary Lubrication by Surface Wettability

Authors: Wei Chen, Adam S. Foster, Mikko J. Alava, Lasse Laurson

Abstract: We study the effect of atomic scale surface-lubricant interactions on nanoscale boundary-lubricated friction, by considering two example surfaces - hydrophilic mica and hydrophobic graphene - confining thin layers of water in molecular dynamics simulations. We observe stick-slip dynamics for thin water films confined by mica sheets, involving periodic breaking-reforming transitions of atomic scale… ▽ More We study the effect of atomic scale surface-lubricant interactions on nanoscale boundary-lubricated friction, by considering two example surfaces - hydrophilic mica and hydrophobic graphene - confining thin layers of water in molecular dynamics simulations. We observe stick-slip dynamics for thin water films confined by mica sheets, involving periodic breaking-reforming transitions of atomic scale capillary water bridges formed around the potassium ions of mica. However, only smooth sliding without stick-slip events is observed for water confined by graphene, as well as for thicker water layers confined by mica. Thus, our results illustrate how atomic scale details affect the wettability of the confining surfaces, and consequently control the presence or absence of stick-slip dynamics in nanoscale friction. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: 5 pages, 5 figures, accepted for publication in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 114, 095502 (2015)

Quenched pinning and collective dislocation dynamics

Authors: Markus Ovaska, Lasse Laurson, Mikko J. Alava

Abstract: Several experiments show that crystalline solids deform in a bursty and intermittent fashion. Power-law distributed strain bursts in compression experiments of micron-sized samples, and acoustic emission energies from larger-scale specimens, are the key signatures of the underlying critical-like collective dislocation dynamics - a phenomenon that has also been seen in discrete dislocation dynamics… ▽ More Several experiments show that crystalline solids deform in a bursty and intermittent fashion. Power-law distributed strain bursts in compression experiments of micron-sized samples, and acoustic emission energies from larger-scale specimens, are the key signatures of the underlying critical-like collective dislocation dynamics - a phenomenon that has also been seen in discrete dislocation dynamics (DDD) simulations. Here we show, by performing large-scale two-dimensional DDD simulations, that the character of the dislocation avalanche dynamics changes upon addition of sufficiently strong randomly distributed quenched pinning centres, present e.g. in many alloys as immobile solute atoms. For intermediate pinning strength, our results adhere to the scaling picture of depinning transitions, in contrast to pure systems where dislocation jamming dominates the avalanche dynamics. Still stronger disorder quenches the critical behaviour entirely. △ Less

Submitted 1 June, 2015; v1 submitted 24 January, 2014; originally announced January 2014.

Comments: 16 pages, 5 figures

Journal ref: Scientific Reports 5, 10580 (2015)

Financial interaction networks inferred from traded volumes

Authors: Hongli Zeng, Rémi Lemoy, Mikko Alava

Abstract: In order to use the advanced inference techniques available for Ising models, we transform complex data (real vectors) into binary strings, by local averaging and thresholding. This transformation introduces parameters, which must be varied to characterize the behaviour of the system. The approach is illustrated on financial data, using three inference methods -- equilibrium, synchronous and async… ▽ More In order to use the advanced inference techniques available for Ising models, we transform complex data (real vectors) into binary strings, by local averaging and thresholding. This transformation introduces parameters, which must be varied to characterize the behaviour of the system. The approach is illustrated on financial data, using three inference methods -- equilibrium, synchronous and asynchronous inference -- to construct functional connections between stocks. We show that the traded volume information is enough to obtain well known results about financial markets, which use however the presumably richer price information: collective behaviour ("market mode") and strong interactions within industry sectors. Synchronous and asynchronous Ising inference methods give results which are coherent with equilibrium ones, and more detailed since the obtained interaction networks are directed. △ Less

Submitted 15 November, 2013; originally announced November 2013.

Comments: 14 pages, 6 figures

A novel local search based on variable-focusing for random K-SAT

Authors: Rémi Lemoy, Mikko Alava, Erik Aurell

Abstract: We introduce a new local search algorithm for satisfiability problems. Usual approaches focus uniformly on unsatisfied clauses. The new method works by picking uniformly random variables in unsatisfied clauses. A Variable-based Focused Metropolis Search (V-FMS) is then applied to random 3-SAT. We show that it is quite comparable in performance to the clause-based FMS. Consequences for algorithmic… ▽ More We introduce a new local search algorithm for satisfiability problems. Usual approaches focus uniformly on unsatisfied clauses. The new method works by picking uniformly random variables in unsatisfied clauses. A Variable-based Focused Metropolis Search (V-FMS) is then applied to random 3-SAT. We show that it is quite comparable in performance to the clause-based FMS. Consequences for algorithmic design are discussed. △ Less

Submitted 12 December, 2013; v1 submitted 9 October, 2013; originally announced November 2013.

Comments: 7 pages, 3 figures

Damage accumulation in quasi-brittle fracture

Authors: Claudio Manzato, Mikko J. Alava, Stefano Zapperi

Abstract: The strength of quasi-brittle materials depends on the ensemble of defects inside the sample and on the way damage accumulates before failure. Using large scale numerical simulations of the random fuse model, we investigate the evolution of the microcrack distribution that is directly related to the strength distribution and its size effects. We show that the broadening of the distribution tail or… ▽ More The strength of quasi-brittle materials depends on the ensemble of defects inside the sample and on the way damage accumulates before failure. Using large scale numerical simulations of the random fuse model, we investigate the evolution of the microcrack distribution that is directly related to the strength distribution and its size effects. We show that the broadening of the distribution tail originates from the dominating microcracks in each sample and is related to a tendency of crack coalescence that increases with system size. We study how the observed behavior depends on the disorder present in the sample. △ Less

Submitted 12 August, 2013; originally announced August 2013.

Comments: 5 pages, 7 figures

Avalanches in 2D Dislocation Systems: Plastic Yielding is not Depinning

Authors: Péter Dusán Ispánovity, Lasse Laurson, Michael Zaiser, István Groma, Stefano Zapperi, Mikko J. Alava

Abstract: We study the properties of strain bursts (dislocation avalanches) occurring in two-dimensional discrete dislocation dynamics models under quasistatic stress-controlled loading. Contrary to previous suggestions, the avalanche statistics differs fundamentally from predictions obtained for the depinning of elastic manifolds in quenched random media. Instead, we find an exponent τ=1 of the power-law d… ▽ More We study the properties of strain bursts (dislocation avalanches) occurring in two-dimensional discrete dislocation dynamics models under quasistatic stress-controlled loading. Contrary to previous suggestions, the avalanche statistics differs fundamentally from predictions obtained for the depinning of elastic manifolds in quenched random media. Instead, we find an exponent τ=1 of the power-law distribution of slip or released energy, with a cut-off that increases exponentially with the applied stress and diverges with system size at all stresses. These observations demonstrate that the avalanche dynamics of 2D dislocation systems is scale-free at every applied stress and, therefore, can not be envisaged in terms of critical behavior associated with a depinning transition. △ Less

Submitted 12 July, 2013; originally announced July 2013.

Comments: 6 pages, 4 figures

Role of the sample thickness in planar crack propagation

Authors: Pallab Barai, Phani K. V. V. Nukala, Mikko J. Alava, Stefano Zapperi

Abstract: We study the effect of the sample thickness in planar crack front propagation in a disordered elastic medium using the random fuse model. We employ different loading conditions and we test their stability with respect to crack growth. We show that the thickness induces characteristic lengths in the stress enhancement factor in front of the crack and in the stress transfer function parallel to the… ▽ More We study the effect of the sample thickness in planar crack front propagation in a disordered elastic medium using the random fuse model. We employ different loading conditions and we test their stability with respect to crack growth. We show that the thickness induces characteristic lengths in the stress enhancement factor in front of the crack and in the stress transfer function parallel to the crack. This is reflected by a thickness-dependent crossover scale in the crack front morphology that goes from from multi-scaling to self-affine with exponents in agreement with line depinning models and experiments. Finally, we compute the distribution of crack avalanches which is shown to depend on the thickness and the loading mode. △ Less

Submitted 4 July, 2013; originally announced July 2013.

Comments: 8 pages, 11 figures

Witness of unsatisfiability for a random 3-satisfiability formula

Authors: Lu-Lu Wu, Hai-Jun Zhou, Mikko Alava, Erik Aurell, Pekka Orponen

Abstract: The random 3-satisfiability (3-SAT) problem is in the unsatisfiable (UNSAT) phase when the clause density $α$ exceeds a critical value $α_s \approx 4.267$. However, rigorously proving the unsatisfiability of a given large 3-SAT instance is extremely difficult. In this paper we apply the mean-field theory of statistical physics to the unsatisfiability problem, and show that a specific type of UNSAT… ▽ More The random 3-satisfiability (3-SAT) problem is in the unsatisfiable (UNSAT) phase when the clause density $α$ exceeds a critical value $α_s \approx 4.267$. However, rigorously proving the unsatisfiability of a given large 3-SAT instance is extremely difficult. In this paper we apply the mean-field theory of statistical physics to the unsatisfiability problem, and show that a specific type of UNSAT witnesses (Feige-Kim-Ofek witnesses) can in principle be constructed when the clause density $α> 19$. We then construct Feige-Kim-Ofek witnesses for single 3-SAT instances through a simple random sampling algorithm and a focused local search algorithm. The random sampling algorithm works only when $α$ scales at least linearly with the variable number $N$, but the focused local search algorithm works for clause densty $α> c N^{b}$ with $b \approx 0.59$ and prefactor $c \approx 8$. The exponent $b$ can be further decreased by enlarging the single parameter $S$ of the focused local search algorithm. △ Less

Submitted 10 March, 2013; originally announced March 2013.

Comments: 9 pages, 7 figures included. Submitted to Physical Review E

Journal ref: Physical Review E 87, 052807 (2013)

Transient Shear Banding in Time-dependent Fluids

Authors: Xavier Illa, Antti Puisto, Arttu Lehtinen, Mikael Mohtaschemi, Mikko J. Alava

Abstract: We study the dynamics of shear-band formation and evolution using a simple rheological model. The description couples the local structure and viscosity to the applied shear stress. We consider in detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation to obtain the time-evolution of the local velocity and viscosity fields. It is found that the underlying… ▽ More We study the dynamics of shear-band formation and evolution using a simple rheological model. The description couples the local structure and viscosity to the applied shear stress. We consider in detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation to obtain the time-evolution of the local velocity and viscosity fields. It is found that the underlying reason for dynamic effects is the non-homogeneous shear distribution, which is amplified due to a positive feedback between the flow field and the viscosity response of the shear thinning fluid. This offers a simple explanation for the recent observations of Transient Shear Banding in time-dependent fluids. Extensions to more complicated rheological systems are considered. △ Less

Submitted 14 February, 2013; originally announced February 2013.

Comments: 6 pages, 5 figures

Maximum likelihood reconstruction for Ising models with asynchronous updates

Authors: Hong-Li Zeng, Mikko Alava, Erik Aurell, John Hertz, Yasser Roudi

Abstract: We describe how the couplings in an asynchronous kinetic Ising model can be inferred. We consider two cases, one in which we know both the spin history and the update times and one in which we only know the spin history. For the first case, we show that one can average over all possible choices of update times to obtain a learning rule that depends only on spin correlations and can also be derived… ▽ More We describe how the couplings in an asynchronous kinetic Ising model can be inferred. We consider two cases, one in which we know both the spin history and the update times and one in which we only know the spin history. For the first case, we show that one can average over all possible choices of update times to obtain a learning rule that depends only on spin correlations and can also be derived from the equations of motion for the correlations. For the second case, the same rule can be derived within a further decoupling approximation. We study all methods numerically for fully asymmetric Sherrington-Kirkpatrick models, varying the data length, system size, temperature, and external field. Good convergence is observed in accordance with the theoretical expectations. △ Less

Submitted 23 April, 2013; v1 submitted 11 September, 2012; originally announced September 2012.

Dynamic hysteresis in cyclic deformation of crystalline solids

Authors: Lasse Laurson, Mikko J. Alava

Abstract: The hysteresis or internal friction in the deformation of crystalline solids stressed cyclically is studied from the viewpoint of collective dislocation dynamics. Stress-controlled simulations of a dislocation dynamics model at various loading frequencies and amplitudes are performed to study the stress - strain rate hysteresis. The hysteresis loop areas exhibit a maximum at a characteristic frequ… ▽ More The hysteresis or internal friction in the deformation of crystalline solids stressed cyclically is studied from the viewpoint of collective dislocation dynamics. Stress-controlled simulations of a dislocation dynamics model at various loading frequencies and amplitudes are performed to study the stress - strain rate hysteresis. The hysteresis loop areas exhibit a maximum at a characteristic frequency and a power law frequency dependence in the low frequency limit, with the power law exponent exhibiting two regimes, corresponding to the jammed and the yielding/moving phases of the system, respectively. The first of these phases exhibits non-trivial critical-like viscoelastic dynamics, crossing over to intermittent viscoplastic deformation for higher stress amplitudes. △ Less

Submitted 10 September, 2012; originally announced September 2012.

Comments: 5 pages, 4 figures, to appear in Physical Review Letters

Journal ref: Phys. Rev. Lett 109, 155504 (2012)

Droplet spreading and pinning on heterogeneous substrates

Authors: Mikko J. Alava, Martin Dube

Abstract: The contact angle of a fluid droplet on an heterogeneous surface is analysed using the statistical dynamics of the spreading contact line. The statistical properties of the final droplet radius and contact angle are obtained through applications of depinning transitions of contact lines with non-local elasticity and features of pinning-depinning dynamics. Such properties not only depend on disorde… ▽ More The contact angle of a fluid droplet on an heterogeneous surface is analysed using the statistical dynamics of the spreading contact line. The statistical properties of the final droplet radius and contact angle are obtained through applications of depinning transitions of contact lines with non-local elasticity and features of pinning-depinning dynamics. Such properties not only depend on disorder strength and surface details, but also on the droplet volume and disorder correlation length. Deviations from Wenzel or Cassie/Baxter behaviour are particularly apparent in the case of small droplet volumes and small contact angles. △ Less

Submitted 10 July, 2012; originally announced July 2012.

Comments: accepted to Phys. Rev. E, 5 figures