BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alain Karma

Alain Karma appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.softphysics.comp-phnlin.PSTissues and OrgansBiological Physicscond-mat.mes-hallmath.DSnlin.CDphysics.app-ph

Catalog footprint

What is connected

18works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

18 published item(s)

preprint2023arXiv

Microstructural Pattern Formation during Far-from-Equilibrium Alloy Solidification

We introduce a new phase-field formulation of rapid alloy solidification that quantitatively incorporates nonequilibrium effects at the solid-liquid interface over a very wide range of interface velocities. Simulations identify a new dynamical instability of dendrite tip growth driven by solute trapping at velocities approaching the absolute stability limit. They also reproduce the formation of the widely observed banded microstructures, revealing how this instability triggers transitions between dendritic and microsegregation-free solidification. Predicted band spacings agree quantitatively with observations in rapidly solidified Al-Cu thin films.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Topological Control of Liquid-Metal-Dealloyed Structures

The past few years have witnessed the rapid development of liquid metal dealloying to fabricate nano-/meso-scale porous and composite structures with ultra-high interfacial area for diverse materials applications. However, this method currently has two important limitations. First, it produces bicontinuous structures with high-genus topologies for a limited range of alloy compositions. Second, structures have a large ligament size due to substantial coarsening during dealloying at high temperature. Here we demonstrate computationally and experimentally that those limitations can be overcome by adding to the metallic melt an element that promotes high-genus topologies by limiting the leakage of the immiscible element during dealloying. We further interpret this finding by showing that bulk diffusive transport of the immiscible element in the liquid melt strongly influences the evolution of the solid fraction and topology of the structure during dealloying. The results shed light on fundamental differences in liquid metal and electrochemical dealloying and establish a new approach to produce liquid-metal-dealloyed structures with desired size and topologies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Isotropic finite-difference approximations for phase-field simulations of polycrystalline alloy solidification

Phase-field models of microstructural pattern formation during alloy solidification are commonly solved numerically using the finite-difference method, which is ideally suited to carry out computationally efficient simulations on massively parallel computer architectures such as Graphic Processing Units. However, one known drawback of this method is that the discretization of differential terms involving spatial derivatives introduces a spurious lattice anisotropy that can influence the solid-liquid interface dynamics. We find that this influence is significant for the case of polycrystalline dendritic solidification, where the crystal axes of different grains do not generally coincide with the reference axes of the finite-difference lattice. In particular, we find that with the commonly used finite-difference implementation of the quantitative phase-field model of binary alloy solidification, both the operating state of the dendrite tip and the dendrite growth orientation are strongly affected by the lattice anisotropy. To circumvent this problem, we use known methods in both real and Fourier space to derive finite-difference approximations of leading differential terms in 2D and 3D that are isotropic at order $h^2$ of the lattice spacing $h$. Importantly, those terms include the divergence of the anti-trapping current that is found to have a critical influence on pattern selection. The 2D and 3D discretizations use an approximated form of the anti-trapping current that facilitates the Fourier-space derivation of the associated isotropic differential operator at $O(h^2)$, but we also derive a 2D discretization of the standard form of this current. Finally, we present 2D and 3D phase-field simulations of alloy solidification, showing that the isotropic finite-difference implementations can dramatically reduce spurious lattice anisotropy effects.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Configurational stability of a crack propagating in a material with mode-dependent fracture energy -- Part II: Drift of fracture facets in mixed-mode I+II+III

In earlier papers (Leblond et.al., 2011, 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a "double'' propagation criterion combining Griffith (1920)'s energetic condition and Goldstein and Salganik (1974)'s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the "threshold'' value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Crack Path Selection in Orientationally Ordered Composites

While cracks in isotropic homogeneous materials propagate straight, perpendicularly to the tensile axis, cracks in natural and synthetic composites deflect from a straight path, often increasing the toughness of the material. Here we combine experiments and simulations to identify materials properties that predict whether cracks propagate straight or kink on a macroscale larger than the composite microstructure. Those properties include the anisotropy of the fracture energy, which we vary several folds by increasing the volume fraction of orientationally ordered alumina platelets inside a polymer matrix, and a microstructure-dependent process zone size that is found to modulate the additional stabilizing or destabilizing effect of the non-singular stress acting parallel to the crack. Those properties predict the existence of an anisotropy threshold for crack kinking and explain the surprisingly strong dependence of this threshold on sample geometry and load distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Vulnerable Window of Yield Strength for Swelling-Driven Fracture of Phase-Transforming Battery Materials

Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith's criterion. Here we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our results reveal the existence of a vulnerable window of yield strength inside which pillars fracture during lithiation. They identify two different modes of fracture inside that window with and without surface localization of plastic deformation prior to fracture for lower and higher yield strength, respectively, and highlight the importance of taking into account this localization to accurately predict the onset of fracture within Griffith theory. The results further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing mechanically stable phase-transforming battery materials undergoing large volume changes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Elastically-mediated interactions between grain boundaries and precipitates in two-phase coherent solids

We investigate analytically and numerically the interaction between grain boundaries and second phase precipitates in two-phase coherent solids in the presence of misfit strain. Our numerical study uses amplitude equations that describe the interaction of composition and stress [R. Spatschek and A. Karma, Phys. Rev. B {\bf 81}, 214201 (2010)] and free-energies corresponding to two-dimensional hexagonal and three-dimensional BCC crystal structures that exhibit isotropic and anisotropic elastic properties, respectively. We consider two experimentally motivated geometries where (i) a lamellar precipitate nucleates along a planar grain boundary that is centered inside the precipitate, and (ii) a circular precipitate nucleates inside a grain at a finite distance to an initially planar grain boundary. For the first geometry, we find that the grain boundary becomes morphologically unstable due to the combination of long-range elastic interaction between the grain boundary and compositional domain boundaries, and shear-coupled grain boundary motion. We characterize this instability analytically by extending the linear stability analysis carried out recently [P.-A. Geslin, Y.-C. Xu, and A. Karma, Phys. Rev. Lett. {\bf 114}, 105501 (2015)] to the more general case of elastic anisotropy. The analysis predicts that elastic anisotropy hinders but does not suppress the instability. Simulations also reveal that, in a well-developed non-linear regime, this instability can lead to the break-up of low-angle grain boundaries when the misfit strain exceeds a threshold that depends on the grain boundary misorientation. For the second geometry, simulations show that the elastic interaction between an initially planar grain boundary and an adjacent circular precipitate causes the precipitate to migrate to and anchor at the grain boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Crack Front Segmentation and Facet Coarsening in Mixed-Mode Fracture

A planar crack generically segments into an array of "daughter cracks" shaped as tilted facets when loaded with both a tensile stress normal to the crack plane (mode I) and a shear stress parallel to the crack front (mode III). We investigate facet propagation and coarsening using in-situ microscopy observations of fracture surfaces at different stages of quasi-static mixed-mode crack propagation and phase-field simulations. The results demonstrate that the bifurcation from propagating planar to segmented crack front is strongly subcritical, reconciling previous theoretical predictions of linear stability analysis with experimental observations. They further show that facet coarsening is a self- similar process driven by a spatial period-doubling instability of facet arrays with a growth rate dependent on mode mixity. Those results have important implications for understanding the failure of a wide range of materials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ginzburg-Landau theory of the bcc-liquid interface kinetic coefficient

We extend the Ginzburg-Landau (GL) theory of atomically rough bcc-liquid interfaces [Wu {\it et al.}, Phys. Rev. B \textbf{73}, 094101 (2006)] outside of equilibrium. We use this extension to derive an analytical expression for the kinetic coefficient, which is the proportionality constant $μ(\hat n)$ between the interface velocity along a direction $\hat n$ normal to the interface and the interface undercooling. The kinetic coefficient is expressed as a spatial integral along the normal direction of a sum of gradient square terms corresponding to different nonlinear density wave profiles. Anisotropy arises naturally from the dependence of those profiles on the angles between the principal reciprocal lattice vectors $\vec K_i$ and $\hat n$. Values of the kinetic coefficient for the$(100)$, $(110)$ and $(111)$ interfaces are compared quantitatively to the prediction of linear Mikheev-Chernov (MC) theory [J. Cryst. Growth \textbf{112}, 591 (1991)] and previous molecular dynamics (MD) simulation studies of crystallization kinetics for a classical model of Fe. Additional MD simulations are carried out here to compute the relaxation time of density waves in the liquid in order to make this comparison free of fit parameter. The GL theory predicts a similar expression for $μ$ as the MC theory but yields a better agreement with MD simulations for both its magnitude and anisotropy due to a fully nonlinear description of density wave profiles across the solid-liquid interface. GL theory is also used to derive an inverse relation between $μ$ and the solid-liquid interfacial free-energy. The general methodology used here to derive an expression for $μ(\hat n)$ also applies to amplitude equations derived from the phase-field-crystal model, which only differ from GL theory by the choice of cubic and higher order nonlinearities in the free-energy density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatiotemporal Dynamics of Calcium-Driven Cardiac Alternans

We investigate the dynamics of spatially discordant alternans (SDA) driven by an instability of intracellular calcium cycling using both amplitude equations [P. S. Skardal, A. Karma, and J. G. Restrepo, Phys. Rev. Lett. 108, 108103 (2012)] and ionic model simulations. We find that, close to the alternans bifurcation, SDA is manifested as a smooth wavy modulation of the amplitudes of both repolarization and calcium transient (CaT) alternans, similarly to the well-studied case of voltage-driven alternans. In contrast, further away from the bifurcation, the amplitude of CaT alternans jumps discontinuously at the nodes separating out-of-phase regions, while the amplitude of repolarization alternans remains smooth. We show that node motion of discontinuous SDA patterns is strongly hysteretic even in homogeneous tissue due to the novel phenomenon of "unidirectional pinning": node movement can only be induced towards, but not away from, the pacing site in response to a change of pacing rate or physiological parameter. In addition, we show that the wavelength of discontinuous SDA patterns scales linearly with the conduction velocity restitution length scale, in contrast to the wavelength of smooth patterns that scales sub-linearly with this length scale. Those results are also shown to be robust against cell-to-cell fluctuations owing to the property that unidirectional node motion collapses multiple jumps accumulating in nodal regions into a single jump. Amplitude equation predictions are in good overall agreement with ionic model simulations. Finally, we briefly discuss physiological implications of our findings. In particular, we suggest that due to the tendency of conduction blocks to form near nodes, the presence of unidirectional pinning makes calcium-driven alternans potentially more arrhythmogenic than voltage-driven alternans.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unified Theoretical Framework for Polycrystalline Pattern Evolution

The rate of curvature-driven grain growth in polycrystalline materials is well-known to be limited by interface dissipation. We show analytically and by simulations that, for systems forming modulated phases or non-equilibrium patterns with crystal ordering, growth is limited by bulk dissipation associated with lattice translation, which dramatically slows down grain coarsening. We also show that bulk dissipation is reduced by thermal noise so that those systems exhibit faster coarsening behavior dominated by interface dissipation for high Peierls barrier and high noise. Those results provide a unified theoretical framework for understanding and modeling polycrystalline pattern evolution in diverse systems over a broad range of length and time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase-field-crystal study of grain boundary premelting and shearing in bcc iron

We use the phase-field-crystal (PFC) method to investigate the equilibrium premelting and nonequilibrium shearing behaviors of $[001]$ symmetric tilt grain boundaries (GBs) at high homologous temperature over the complete range of misorientation $0<θ<90^\circ$ in classical models of bcc Fe. We characterize the dependence of the premelted layer width $W$ as a function of temperature and misorientation and compute the thermodynamic disjoining potential whose derivative with respect to $W$ represents the structural force between crystal-melt interfaces due to the spatial overlap of density waves. The disjoining potential is also computed by molecular dynamics (MD) simulations, for quantitative comparison with PFC simulations, and coarse-grained amplitude equations (AE) derived from PFC that provide additional analytical insights. We find that, for GBs over an intermediate range of misorientation ($θ_{\rm min}<θ<θ_{\rm max}$), $W$ diverges as the melting temperature is approached from below, corresponding to a purely repulsive disjoining potential, while for GBs outside this range ($θ<θ_{\rm min}$ or $θ_{\rm max}<θ<90^\circ$), $W$ remains finite at the melting point, with its value corresponding to a shallow attractive minimum of the disjoining potential. In response to a shear stress parallel to the GB plane, GBs in PFC simulations exhibit coupled motion normal to this plane, with a discontinuous change of the coupling factor as a function of $θ$ that reflects a transition between two coupling modes, and/or sliding (shearing of the two grains). The coupling factor for the two coupling modes is in excellent quantitative agreement with previous theoretical predictions [J. W. Cahn, Y. Mishin, and A. Suzuki, Acta Mater. 54, 4953 (2006)].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structural short-range forces between solid-melt interfaces

We predict the structural interaction of crystalline solid-melt interfaces using amplitude equations which are derived from classical density functional theory or phase-field-crystal modeling. The solid ordering decays exponentially on the scale of the interface thickness at solid-melt interfaces; the overlap of two such profiles leads to a short range interaction, which is mainly carried by the longest-range density waves, which can facilitate grain boundary premelting. We calculate the tail of these interactions, depending on the relative translation of the two crystals fully analytically and predict the interaction potential, and compare it to numerical simulations. For grain boundaries the interaction is predicted to decay twice faster as for two crystals without misorientation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unidirectional Pinning and Hysteresis of Spatially Discordant Alternans in Cardiac Tissue

Spatially discordant alternans is a widely observed pattern of voltage and calcium signals in cardiac tissue that can precipitate lethal cardiac arrhythmia. Using spatially coupled iterative maps of the beat-to-beat dynamics, we explore this pattern's dynamics in the regime of a calcium-dominated period-doubling instability at the single cell level. We find a novel nonlinear bifurcation associated with the formation of a discontinuous jump in the amplitude of calcium alternans at nodal lines separating discordant regions. We show that this jump unidirectionally pins nodal lines by preventing their motion away from the pacing site following a pacing rate decrease, but permitting motion towards this site following a rate increase. This unidirectional pinning leads to strongly history-dependent nodal line motion that is strongly arrhythmogenic.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Amplitude equations for polycrystalline materials with interaction between composition and stress

We investigate the ability of frame-invariant amplitude equations [G. H. Gunaratne, Q. Ouyang, and H. Swinney, Phys. Rev. E {\bf 50}, 2802 (1994)] to describe quantitatively the evolution of polycrystalline microstructures and we extend this approach to include the interaction between composition and stress. Validations for elemental materials include studies of the Asaro-Tiller-Grinfeld morphological instability of a stressed crystal surface, polycrystalline growth from the melt, grain boundary energies over a wide range of misorientation, and grain boundary motion coupled to shear deformation. Amplitude equations with accelerated strain relaxation in the solid are shown to model accurately the Asaro-Tiller-Grinfeld instability. Polycrystalline growth is also well described. However, the survey of grain boundary energies shows that the approach is only valid for a restricted range of misorientations as a direct consequence of an amplitude expansion. This range covers approximately half the complete range allowed by crystal symmetry for some fixed reference set of density waves used in the expansion. Over this range, coupled motion to shear is well described by known geometrical rules and a transition from coupling to sliding motion is also reproduced. Amplitude equations for alloys are derived phenomenologically in a Ginzburg-Landau spirit. Vegard's law is shown to be naturally described by seeking a gauge invariant form of those equations under a transformation that corresponds to a lattice expansion and deviations from Vegard's law can be easily incorporated. Those equations realistically describe the dilute alloy limit and have the same flexibility as conventional phase-field models for incorporating arbitrary free-energy/composition curves...

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase-field-crystal model for fcc ordering

We develop and analyze a two-mode phase-field-crystal model to describe fcc ordering. The model is formulated by coupling two different sets of crystal density waves corresponding to <111> and <200> reciprocal lattice vectors, which are chosen to form triads so as to produce a simple free- energy landscape with coexistence of crystal and liquid phases. The feasibility of the approach is demonstrated with numerical examples of polycrystalline and (111) twin growth. We use a two-mode amplitude expansion to characterize analytically the free-energy landscape of the model, identifying parameter ranges where fcc is stable or metastable with respect to bcc. In addition, we derive analytical expressions for the elastic constants for both fcc and bcc. Those expressions show that a non-vanishing amplitude of [200] density waves is essential to obtain mechanically stable fcc crystals with a non-vanishing tetragonal shear modulus (C11 - C12)/2. We determine the model parameters for specific materials by fitting the peak liquid structure factor properties and solid density wave amplitudes following the approach developed for bcc [K.-A. Wu and A. Karma, Phys. Rev. B 76, 184107 (2007)]. This procedure yields reasonable predictions of elastic constants for both bcc Fe and fcc Ni using input parameters from molecular dynamics simulations. The application of the model to two-dimensional square lattices is also briefly examined.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structural disjoining potential for grain boundary premelting and grain coalescence from molecular-dynamics simulations

We describe a molecular dynamics framework for the direct calculation of the short-ranged structural forces underlying grain-boundary premelting and grain-coalescence in solidification. The method is applied in a comparative study of (i) a Sigma 9 <115> 120 degress twist and (ii) a Sigma 9 <110> {411} symmetric tilt boundary in a classical embedded-atom model of elemental Ni. Although both boundaries feature highly disordered structures near the melting point, the nature of the temperature dependence of the width of the disordered regions in these boundaries is qualitatively different. The former boundary displays behavior consistent with a logarithmically diverging premelted layer thickness as the melting temperature is approached from below, while the latter displays behavior featuring a finite grain-boundary width at the melting point. It is demonstrated that both types of behavior can be quantitatively described within a sharp-interface thermodynamic formalism involving a width-dependent interfacial free energy, referred to as the disjoining potential. The disjoining potential for boundary (i) is calculated to display a monotonic exponential dependence on width, while that of boundary (ii) features a weak attractive minimum. The results of this work are discussed in relation to recent simulation and theoretical studies of the thermodynamic forces underlying grain-boundary premelting.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Phase-Field Formulation for Quantitative Modeling of Alloy Solidification

A phase-field formulation is introduced to simulate quantitatively microstructural pattern formation in alloys. The thin-interface limit of this formulation yields a much less stringent restriction on the choice of interface thickness than previous formulations and permits to eliminate non-equilibrium effects at the interface. Dendrite growth simulations with vanishing solid diffusivity show that both the interface evolution and the solute profile in the solid are well resolved.

0 citations0 reviewsNo code linkedOpen access