Source author record

Alain Goriely

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

22 published item(s)

preprint2026arXiv

A whole-brain model of amyloid beta accumulation and cerebral hypoperfusion in Alzheimer's disease

Accumulation of amyloid beta proteins is a defining feature of Alzheimer's disease, and is usually accompanied by cerebrovascular pathology. Evidence suggests that amyloid beta and cerebrovascular pathology are mutually reinforcing; in particular, amyloid beta suppresses perfusion by constricting capillaries, and hypoperfusion promotes the production of amyloid beta. Here, we propose a whole-brain model coupling amyloid beta and blood vessel through a hybrid model consisting of a reaction-diffusion system for the protein dynamics and porous-medium model of blood flow within and between vascular networks: arterial, capillary and venous. We discretize the resulting parabolic--elliptic system of PDEs by means of a high-order discontinuous Galerkin method in space and an implicit Euler scheme in time. Simulations in realistic brain geometries demonstrate the emergence of multistability, implying that a sufficiently large pathogenic protein seeds is necessary to trigger disease outbreak. Motivated by the "two-hit vascular hypothesis" of Alzheimer's disease that hypoperfusive vascular damage triggers amyloid beta pathology, we also demonstrate that localized hypoperfusion, in response to injury, can destabilize the healthy steady state and trigger brain-wide disease outbreak.

preprint2022arXiv

Barcodes distinguish morphology of neuronal tauopathy

The geometry of neurons is known to be important for their functions. Hence, neurons are often classified by their morphology. Two recent methods, persistent homology and the topological morphology descriptor, assign a morphology descriptor called a barcode to a neuron equipped with a given function, such as the Euclidean distance from the root of the neuron. These barcodes can be converted into matrices called persistence images, which can then be averaged across groups. We show that when the defining function is the path length from the root, both the topological morphology descriptor and persistent homology are equivalent. We further show that persistence images arising from the path length procedure provide an interpretable summary of neuronal morphology. We introduce {topological morphology functions}, a class of functions similar to Sholl functions, that can be recovered from the associated topological morphology descriptor. To demonstrate this topological approach, we compare healthy cortical and hippocampal mouse neurons to those affected by progressive tauopathy. We find a significant difference in the morphology of healthy neurons and those with a tauopathy at a postsymptomatic age. We use persistence images to conclude that the diseased group tends to have neurons with shorter branches as well as fewer branches far from the soma.

preprint2020arXiv

A plate theory for nematic liquid crystalline solids

We derive a Föppl-von Kármán-type constitutive model for solid liquid crystalline plates where the nematic director may or may not rotate freely relative to the elastic network. To obtain the reduced two-dimensional model, we rely on the deformation decomposition of a nematic solid into an elastic deformation and a natural shape change. The full solution to the resulting equilibrium equations consists of both the deformation displacement and stress fields. The model equations are applicable to a wide range of thin nematic bodies subject to optothermal stimuli and mechanical loads. For illustration, we consider certain reversible natural shape changes in simple systems which are stress free, and their counterparts, where the natural deformations are blocked and internal stresses appear. More general problems can be addressed within the same framework.

preprint2020arXiv

Dynamic buckling of an inextensible elastic ring: Linear and nonlinear analyses

Slender elastic objects such as a column tend to buckle under loads. While static buckling is well understood as a bifurcation problem, the evolution of shapes during dynamic buckling is much harder to study. Elastic rings under normal pressure have emerged as a theoretical and experimental paradigm for the study of dynamic buckling with controlled loads. Experimentally, an elastic ring is placed within a soap film. When the film outside the ring is removed, surface tension pulls the ring inward, mimicking an external pressurization. Here we present a theoretical analysis of this process by performing a post-bifurcation analysis of an elastic ring under pressure. This analysis allows us to understand how inertia, material properties, and loading affect the observed shape. In particular, we combine direct numerical solutions with a post-bifurcation asymptotic analysis to show that inertia drives the system towards higher modes that cannot be selected in static buckling. Our theoretical results explain experimental observations that cannot be captured by a standard linear stability analysis.

preprint2020arXiv

Dynamic buckling of elastic rings in a soap film

Dynamic buckling may occur when a load is rapidly applied to, or removed from, an elastic object at rest. In contrast to its static counterpart, dynamic buckling offers a wide range of accessible patterns depending on the parameters of the system and the dynamics of the load. To study these effects, we consider experimentally the dynamics of an elastic ring in a soap film when part of the film is suddenly removed. The resulting change in tension applied to the ring creates a range of interesting patterns that cannot be easily accessed in static experiments. Depending on the aspect ratio of the ring's cross section, high-mode buckling patterns are found in the plane of the remaining soap film or out of the plane. Paradoxically, while inertia is required to observe these non-trivial modes, the selected pattern does not depend on inertia itself. The evolution of this pattern beyond the initial instability is studied experimentally and explained through theoretical arguments linking dynamics to pattern selection and mode growth. We also explore the influence of dynamic loading and show numerically that by imposing a rate of loading that competes with the growth rate of instability, the observed pattern can be selected and controlled.

preprint2020arXiv

Likely cavitation and radial motion of stochastic elastic spheres

The cavitation of solid elastic spheres is a classical problem of continuum mechanics. Here, we study this problem within the context of "stochastic elasticity" where the constitutive parameters are characterised by probability density functions. We consider homogeneous spheres of stochastic neo-Hookean material, composites with two concentric stochastic neo-Hookean phases, and inhomogeneous spheres of locally neo-Hookean material with a radially varying parameter. In all cases, we show that the material at the centre determines the critical load at which a spherical cavity forms there. However, while under dead-load traction, a supercritical bifurcation, with stable cavitation, is obtained in a static sphere of stochastic neo-Hookean material, for the composite and radially inhomogeneous spheres, a subcritical bifurcation, with snap cavitation, is also possible. For the dynamic spheres, oscillatory motions are produced under suitable dead-load traction, such that a cavity forms and expands to a maximum radius, then collapses again to zero periodically, but not under impulse traction. Under a surface impulse, a subcritical bifurcation is found in a static sphere of stochastic neo-Hookean material and also in an inhomogeneous sphere, whereas in composite spheres, supercritical bifurcations can occur as well. Given the non-deterministic material parameters, the results can be characterised in terms of probability distributions.

preprint2020arXiv

Likely cavitation and radial motion of stochastic elastic spheres 2: Impulse driven

Cavitation in solids can be caused by tensile dead-load traction or impulse traction. The two different types of boundary conditions lead to different static and dynamic solutions. In addition, if the material is stochastic, i.e., the model parameters are represented by probability distributions, the expected behaviour is more complicated to describe. Here, following the first instalment of this work, we examine the static and dynamic cavitation of a stochastic material under a uniform tensile impulse traction in different spherical geometries. We find that the critical load at which a cavity forms at the centre of the sphere is the same as for the homogeneous sphere composed entirely of the material found at its centre, while the post-cavitation radial motion is non-oscillatory. However, there are some important differences in the nonlinear elastic responses. Specifically, subcritical bifurcation, with snap cavitation, is obtained in a static sphere of stochastic neo-Hookean material and also in a radially inhomogeneous sphere, whereas for composite spheres, a supercritical bifurcation, with stable cavitation, is possible as well. Given the non-deterministic material parameters, the results are characterised by probability distributions.

preprint2020arXiv

Likely striping in stochastic nematic elastomers

For monodomain nematic elastomers, we construct generalised elastic-nematic constitutive models combining purely elastic and neoclassical-type strain-energy densities. Inspired by recent developments in stochastic elasticity, we extend these models to stochastic-elastic-nematic forms where the model parameters are defined by spatially-independent probability density functions at a continuum level. To investigate the behaviour of these systems and demonstrate the effects of the probabilistic parameters, we focus on the classical problem of shear striping in a stretched nematic elastomer for which the solution is given explicitly. We find that, unlike in the neoclassical case where the inhomogeneous deformation occurs within a universal interval that is independent of the elastic modulus, for the elastic-nematic models, the critical interval depends on the material parameters. For the stochastic extension, the bounds of this interval are probabilistic, and the homogeneous and inhomogeneous states compete in the sense that both have a a given probability to occur. We refer to the inhomogeneous pattern within this interval as "likely striping".

preprint2020arXiv

Neuronal Oscillations on Evolving Networks: Dynamics, Damage, Degradation, Decline, Dementia, and Death

Neurodegenerative diseases, such as Alzheimer's or Parkinson's disease, show characteristic degradation of structural brain networks. This degradation eventually leads to changes in the network dynamics and degradation of cognitive functions. Here, we model the progression in terms of coupled physical processes: The accumulation of toxic proteins, given by a nonlinear reaction-diffusion transport process, yields an evolving brain connectome characterized by weighted edges on which a neuronal-mass model evolves. The progression of the brain functions can be tested by simulating the resting-state activity on the evolving brain network. We show that while the evolution of edge weights plays a minor role in the overall progression of the disease, dynamic biomarkers predict a transition over a period of 10 years associated with strong cognitive decline.

preprint2018arXiv

The mechanics of human brain organoids

Organoids are prototypes of human organs derived from cultured human stem cells. They provide a reliable and accurate experimental model to study the physical mechanisms underlying the early developmental stages of human organs morphogenesis and, in particular, the early morphogenesis of the cortex. Here, we propose a mathematical model to elucidate the role played by two mechanisms which have been experimentally proven to be crucial in shaping human brain organoids: the contraction of the inner core of the organoid and the microstructural remodeling of the outer cortex. Our results show that both mechanisms are crucial for the final shape of the organoid and can explain the origin of brain pathologies such as lissencephaly (smooth brain).

preprint2016arXiv

Capacitance-Power-Hysteresis Trilemma in Nanoporous Supercapacitors

Nanoporous supercapacitors are an important player in the field of energy storage that fill the gap between dielectric capacitors and batteries. The key challenge in the development of supercapacitors is the perceived trade-off between capacitance and power delivery. Current efforts to boost the capacitance of nanoporous supercapacitors focus on reducing the pore size so that they can only accommodate a single layer of ions. However, this tight packing compromises the charging dynamics and hence power density. We show via an analytical theory and Monte Carlo simulations that charging is sensitively dependent on the affinity of ions to the pores, and that high capacitances can be obtained for ionophobic pores of widths significantly larger than the ion diameter. Our theory also predicts that charging can be hysteretic with a significant energy loss per cycle for intermediate ionophilicities. We use these observations to explore the parameter regimes in which a capacitance-power-hysteresis trilemma may be avoided.

preprint2016arXiv

Design and stability of a family of deployable structures

A large family of deployable filamentary structures can be built by connecting two elastic rods along their length. The resulting structure has interesting shapes that can be stabilized by tuning the material properties of each rod. To model this structure and study its stability, we show that the equilibrium equations describing unloaded states can be derived from a variational principle. We then use a novel geometric method to study the stability of the resulting equilibria. As an example we apply the theory to establish the stability of all possible equilibria of the Bristol ladder.

preprint2016arXiv

Growth, collapse, and stalling in a mechanical model for neurite motility

Neurites, the long cellular protrusions that form the routes of the neuronal network are capable to actively extend during early morphogenesis or to regenerate after trauma. To perform this task, they rely on their cytoskeleton for mechanical support. In this paper, we present a three-component active gel model that describes neurites in the three robust mechanical states observed experimentally: collapsed, static, and motile. These states arise from an interplay between the physical forces driven by growth of the microtubule-rich inner core of the neurite and the acto-myosin contractility of its surrounding cortical membrane. In particular, static states appear as a mechanical traction/compression balance of these two parallel structures. The model predicts how the response of a neurite to a towing force depends on the force magnitude and recovers the response of neurites to several drug treatments that modulate the cytoskeleton active and passive properties.

preprint2016arXiv

The surprising dynamics of a chain on a pulley: Lift-off and snapping

The motion of weights attached to a chain or string moving on a frictionless pulley is a classic problem of introductory physics used to understand the relationship between force and acceleration. Here, we consider the dynamics of the chain when one of the weights is removed and, thus, one end is pulled with constant acceleration. This simple change has dramatic consequences for the ensuing motion: at a finite time, the chain `lifts off' from the pulley and the free end subsequently accelerates faster than the end that is pulled. Eventually, the chain undergoes a dramatic reversal of curvature reminiscent of the crack, or snap, of a whip. We combine experiments, numerical simulations, and theoretical arguments to explain key aspects of this dynamical problem.

preprint2015arXiv

Dynamics of Ion Transport in Ionic Liquids

A gap in understanding the link between continuum theories of ion transport in ionic liquids and the underlying microscopic dynamics has hindered the development of frameworks for transport phenomena in these concentrated electrolytes. Here, we construct a continuum theory for ion transport in ionic liquids by coarse graining a simple exclusion process of interacting particles on a lattice. The resulting dynamical equations can be written as a gradient flow with a mobility matrix that vanishes at high densities. This form of the mobility matrix gives rise to a charging behaviour that is different to the one known for electrolytic solutions, but which agrees qualitatively with the phenomenology observed in experiments and simulations.

preprint2014arXiv

Are Room Temperature Ionic Liquids Dilute Electrolytes?

An important question in understanding the structure of ionic liquids is whether ions are truly "free" and mobile which would correspond to a concentrated ionic melt, or are rather "bound" in ion pairs, that is a liquid of ion pairs with a small concentration of free ions. Recent surface force balance experiments from different groups have given conflicting answers to this question. We propose a simple model for the thermodynamics and kinetics of ion pairing in ionic liquids. Our model takes into account screened ion-ion, dipole-dipole and dipole-ion interactions in the mean field limit. The results of this model suggest that almost two thirds of the ions are free at any instant, and ion pairs have a short lifetime comparable to the characteristic timescale for diffusion. These results suggest that there is no particular thermodynamic or kinetic preference for ions residing in pairs. We therefore conclude that ionic liquids are concentrated, rather than dilute, electrolytes.

preprint2014arXiv

Controlled topological transitions in thin film phase separation

In this paper the evolution of a binary mixture in a thin-film geometry with a wall at the top and bottom is considered. By bringing the mixture into its miscibility gap so that no spinodal decomposition occurs in the bulk, a slight energetic bias of the walls towards each one of the constituents ensures the nucleation of thin boundary layers that grow until the constituents have moved into one of the two layers. These layers are separated by an interfacial region where the composition changes rapidly. Conditions that ensure the separation into two layers with a thin interfacial region are investigated based on a phase-field model. Using matched asymptotic expansions a corresponding sharp-interface problem for the location of the interface is established. It is then argued that this newly created two-layer system is not at its energetic minimum but destabilizes into a controlled self-replicating pattern of trapezoidal vertical stripes by minimizing the interfacial energy between the phases while conserving their area. A quantitative analysis of this mechanism is carried out via a thin-film model for the free interfaces, which is derived asymptotically from the sharp-interface model.

preprint2014arXiv

Unravelling Nanoconfined Films of Ionic Liquids

The confinement of an ionic liquid between charged solid surfaces is treated using an exactly solvable 1D Coulomb gas model. The theory highlights the importance of two dimensionless parameters: the fugacity of the ionic liquid, and the electrostatic interaction energy of ions at closest approach relative to thermal energy, in determining how the disjoining pressure exerted on the walls depends on the geometrical confinement. Our theory reveals that thermodynamic fluctuations play a vital role in the "squeezing out" of charged layers as the confinement is increased. The model shows good qualitative agreement with previous experimental data, with all parameters independently estimated without fitting.

preprint2013arXiv

Magnetic chains: From self-buckling to self-assembly

Spherical neodymium-iron-boron magnets are perman-ent magnets that can be assembled into a variety of structures due to their high magnetic strength. A one-dimensional chain of these magnets responds to mechanical loadings in a manner reminiscent of an elastic rod. We investigate the macroscopic mechanical properties of assemblies of ferromagnetic spheres by considering chains, rings, and chiral cylinders of magnets. Based on energy estimates and simple experiments, we introduce an effective magnetic bending stiffness for a chain of magnets and show that, used in conjunction with classic results for elastic rods, it provides excellent estimates for the buckling and vibration dynamics of magnetic chains. We then use this estimate to understand the dynamic self-assembly of a cylinder from an initially straight chain of magnets.

preprint2013arXiv

Nonlinear Correction to the Euler Buckling Formula for Compressed Cylinders with Guided-Guided End Conditions

Euler's celebrated buckling formula gives the critical load $N$ for the buckling of a slender cylindrical column with radius $B$ and length $L$ as \[ N / (π^3 B^2) = (E/4)(B/L)^2, \] where $E$ is Young's modulus. Its derivation relies on the assumptions that linear elasticity applies to this problem, and that the slenderness $(B/L)$ is an infinitesimal quantity. Here we ask the following question: What is the first nonlinear correction in the right hand-side of this equation when terms up to $(B/L)^4$ are kept? To answer this question, we specialize the exact solution of incremental non-linear elasticity for the homogeneous compression of a thick compressible cylinder with lubricated ends to the theory of third-order elasticity. In particular, we highlight the way second- and third-order constants ---including Poisson's ratio--- all appear in the coefficient of $(B/L)^4$.

preprint2013arXiv

Scalar evolution equations for shear waves in incompressible solids: A simple derivation of the Z, ZK, KZK, and KP equations

We study the propagation of two-dimensional finite-amplitude shear waves in a nonlinear pre-strained incompressible solid, and derive several asymptotic amplitude equations in a simple, consistent, and rigorous manner. The scalar Zabolotskaya (Z) equation is shown to be the asymptotic limit of the equations of motion for all elastic generalized neo-Hookean solids (with strain energy depending only on the first principal invariant of Cauchy-Green strain). However, we show that the Z equation cannot be a scalar equation for the propagation of two-dimensional shear waves in general elastic materials (with strain energy depending on the first and second principal invariants of strain). Then we introduce dispersive and dissipative terms to deduce the scalar Kadomtsev-Petviashvili (KP), Zabolotskaya-Khokhlov (ZK) and Khokhlov-Zabolotskaya-Kuznetsov (KZK) equations of incompressible solid mechanics.

preprint2013arXiv

The mechanics of a chain or ring of spherical magnets

Strong magnets, such as neodymium-iron-boron magnets, are increasingly being manufactured as spheres. Because of their dipolar characters, these spheres can easily be arranged into long chains that exhibit mechanical properties reminiscent of elastic strings or rods. While simple formulations exist for the energy of a deformed elastic rod, it is not clear whether or not they are also appropriate for a chain of spherical magnets. In this paper, we use discrete-to-continuum asymptotic analysis to derive a continuum model for the energy of a deformed chain of magnets based on the magnetostatic interactions between individual spheres. We find that the mechanical properties of a chain of magnets differ significantly from those of an elastic rod: while both magnetic chains and elastic rods support bending by change of local curvature, nonlocal interaction terms also appear in the energy formulation for a magnetic chain. This continuum model for the energy of a chain of magnets is used to analyse small deformations of a circular ring of magnets and hence obtain theoretical predictions for the vibrational modes of a circular ring of magnets. Surprisingly, despite the contribution of nonlocal energy terms, we find that the vibrations of a circular ring of magnets are governed by the same equation that governs the vibrations of a circular elastic ring.