Showing 1–50 of 56 results for author: Lenz, M

Achieving specific yet transient bonds between anisotropic colloids

Authors: Muraleedharapai Mayarani, Martin Lenz, Olivia du Roure, Julien Heuvingh

Abstract: Self-assembly of colloidal particles is a promising avenue to control the shape and dynamics of larger aggregates. However, achieving the necessary fine control over the dynamics and specificity of the bonds between such particles remains a challenge. Here we demonstrate such control in bonds mediated by depletion interactions between anisotropic colloids that we 3D-print in the shape of half disk… ▽ More Self-assembly of colloidal particles is a promising avenue to control the shape and dynamics of larger aggregates. However, achieving the necessary fine control over the dynamics and specificity of the bonds between such particles remains a challenge. Here we demonstrate such control in bonds mediated by depletion interactions between anisotropic colloids that we 3D-print in the shape of half disks with sub-micron resolution. When brought together by diffusion, the particles interact in different configurations but the interaction through the flat faces is by far the longest-lasting. All bonds are flexible and transient, and we demonstrate control over their life time through the depletant concentration in quantitative agreement with a simple physical model. This basic design could be extended to manufacture particles with multiple binding sites to engineer directional assembly with multiple particles. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Comments: Total submission: 23 pages (Article: 15 pages followed by Supporting Information: 8 pages )

How do particles with complex interactions self-assemble?

Authors: Lara Koehler, Pierre Ronceray, Martin Lenz

Abstract: In living cells, proteins self-assemble into large functional structures based on specific interactions between molecularly complex patches. Due to this complexity, protein self-assembly results from a competition between a large number of distinct interaction energies, of the order of one per pair of patches. Current self-assembly models however typically ignore this aspect, and the principles by… ▽ More In living cells, proteins self-assemble into large functional structures based on specific interactions between molecularly complex patches. Due to this complexity, protein self-assembly results from a competition between a large number of distinct interaction energies, of the order of one per pair of patches. Current self-assembly models however typically ignore this aspect, and the principles by which it determines the large-scale structure of protein assemblies are largely unknown. Here, we use Monte-Carlo simulations and machine learning to start to unravel these principles. We observe that despite widespread geometrical frustration, aggregates of particles with complex interactions fall within only a few categories that often display high degrees of spatial order, including crystals, fibers, and oligomers. We then successfully identify the most relevant aspect of the interaction complexity in predicting these outcomes, namely the particles' ability to form periodic structures. Our results provide a first extensive characterization of the rich design space associated with identical particles with complex interactions, and could inspire engineered self-assembling nanoobjects as well as help understand the emergence of robust functional protein structures. △ Less

Submitted 29 July, 2024; v1 submitted 11 December, 2023; originally announced December 2023.

Comments: 12 pages, 10 figures. Supplementary Information: 23 pages, 19 figures

Collective deformation modes promote fibrous self-assembly in protein-like particles

Authors: Hugo Le Roy, M. Mert Terzi, Martin lenz

Abstract: The self-assembly of particles into organized structures is a key feature of living organisms and a major engineering challenge. While it may proceed through the binding of perfectly matched, puzzle-pieces-like particles, many other instances involve ill-fitting particles that must deform to fit together. These include some pathological proteins, which have a known propensity to form fibrous aggre… ▽ More The self-assembly of particles into organized structures is a key feature of living organisms and a major engineering challenge. While it may proceed through the binding of perfectly matched, puzzle-pieces-like particles, many other instances involve ill-fitting particles that must deform to fit together. These include some pathological proteins, which have a known propensity to form fibrous aggregates. Despite this observation, the general relationship between the individual characteristics of the particles and the overall structure of the aggregate is not understood. To elucidate it, we analytically and numerically study the self-assembly of two-dimensional, deformable ill-fitting particles. We find that moderately sticky particles tend to form equilibrium self-limited aggregates whose size is set by an elastic boundary layer associated with collective deformations that may extend over many particles. Particles with a soft internal deformation mode thus give rise to large aggregates. Besides, when the particles are incompressible, their aggregates tend to be anisotropic and fiber-like. Our results are preserved in a more complex particle model with randomly chosen elastic properties. This indicates that generic protein characteristics such as allostery and incompressibility could favor the formation of fibers in protein aggregation, and suggests design principles for artificial self-assembling structures. △ Less

Submitted 5 April, 2024; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: 21 pages, 12 figures

Gender mobility in the labor market with skills-based matching models

Authors: Ajaya Adhikari, Steven Vethman, Daan Vos, Marc Lenz, Ioana Cocu, Ioannis Tolios, Cor J. Veenman

Abstract: Skills-based matching promises mobility of workers between different sectors and occupations in the labor market. In this case, job seekers can look for jobs they do not yet have experience in, but for which they do have relevant skills. Currently, there are multiple occupations with a skewed gender distribution. For skills-based matching, it is unclear if and how a shift in the gender distributio… ▽ More Skills-based matching promises mobility of workers between different sectors and occupations in the labor market. In this case, job seekers can look for jobs they do not yet have experience in, but for which they do have relevant skills. Currently, there are multiple occupations with a skewed gender distribution. For skills-based matching, it is unclear if and how a shift in the gender distribution, which we call gender mobility, between occupations will be effected. It is expected that the skills-based matching approach will likely be data-driven, including computational language models and supervised learning methods. This work, first, shows the presence of gender segregation in language model-based skills representation of occupations. Second, we assess the use of these representations in a potential application based on simulated data, and show that the gender segregation is propagated by various data-driven skills-based matching models.These models are based on different language representations (bag of words, word2vec, and BERT), and distance metrics (static and machine learning-based). Accordingly, we show how skills-based matching approaches can be evaluated and compared on matching performance as well as on the risk of gender segregation. Making the gender segregation bias of models more explicit can help in generating healthy trust in the use of these models in practice. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: This paper was presented during the AAAI Spring Symposium 2023 (AI Trustworthiness Assessment (AITA) track)

Transverse fluctuations control the assembly of semiflexible filaments

Authors: Valerio Sorichetti, Martin Lenz

Abstract: The kinetics of the assembly of semiflexible filaments through end-to-end annealing is key to the structure of the cytoskeleton, but is not understood. We analyze this problem through scaling theory and simulations, and uncover a regime where filaments ends find each other through bending fluctuations without the need for the whole filament to diffuse. This results in a very substantial speed-up o… ▽ More The kinetics of the assembly of semiflexible filaments through end-to-end annealing is key to the structure of the cytoskeleton, but is not understood. We analyze this problem through scaling theory and simulations, and uncover a regime where filaments ends find each other through bending fluctuations without the need for the whole filament to diffuse. This results in a very substantial speed-up of assembly in physiological regimes, and could help understand the dynamics of actin and intermediate filaments in biological processes such as wound healing and cell division. △ Less

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

Comments: Main: 5 pages, 4 figures. Supplementary: 10 pages, 9 figures. To be published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 131, 228401 (2023)

Runaway Transition in Irreversible Polymer Condensation with Cyclisation

Authors: Maria Panoukidou, Simon Weir, Valerio Sorichetti, Yair Gutierrez Fosado, Martin Lenz, Davide Michieletto

Abstract: The process of polymer condensation, i.e. the formation of bonds between reactive end-groups, is ubiquitous in both industry and biology. Here we study generic systems undergoing polymer condensation in competition with cyclisation. Using a generalised Smoluchowski theory, molecular dynamics simulations and experiments using DNA and T4 ligase, we find that this system displays a transition, from a… ▽ More The process of polymer condensation, i.e. the formation of bonds between reactive end-groups, is ubiquitous in both industry and biology. Here we study generic systems undergoing polymer condensation in competition with cyclisation. Using a generalised Smoluchowski theory, molecular dynamics simulations and experiments using DNA and T4 ligase, we find that this system displays a transition, from a regime with finite-length chains at infinite time and dominated by rings to one dominated by linear polymers that grow in time. Finally, we show that fluids prepared close to the transition may have profoundly different compositions and rheology at large condensation times. △ Less

Submitted 3 June, 2024; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: 6 pages, 7 figures (main)

Electro-optic sampling based characterization of broad-band high efficiency THz-FEL

Authors: Maximilian Lenz, Andrew Fisher, Youna Park, Alexander Ody, Pietro Musumeci

Abstract: Extremely high beam-to-radiation energy conversion efficiencies can be obtained in a THz FEL using a strongly tapered helical undulator at the zero-slippage resonant condition, where a circular waveguide is used to match the radiation group velocity to the electron beam longitudinal velocity. In this paper we report on the first electro-optic sampling (EOS) based measurements of the broadband THz… ▽ More Extremely high beam-to-radiation energy conversion efficiencies can be obtained in a THz FEL using a strongly tapered helical undulator at the zero-slippage resonant condition, where a circular waveguide is used to match the radiation group velocity to the electron beam longitudinal velocity. In this paper we report on the first electro-optic sampling (EOS) based measurements of the broadband THz FEL radiation pulses emitted in this regime. The THz field waveforms are reconstructed in the spatial and temporal domains using multi-shot and single-shot EOS schemes respectively. The measurements are performed varying the input electron beam energy in the undulator providing insights on the complex dynamics in a waveguide FEL. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Generic Stress Rectification in Nonlinear Elastic Media

Authors: Félix Benoist, Guglielmo Saggiorato, Martin Lenz

Abstract: Stress propagation in nonlinear media is crucial in cell biology, where molecular motors exert anisotropic force dipoles on the fibrous cytoskeleton. While the force dipoles can be either contractile or expansile, a medium made of fibers which buckle under compression rectifies these stresses towards a biologically crucial contraction. A general understanding of this rectification phenomenon as a… ▽ More Stress propagation in nonlinear media is crucial in cell biology, where molecular motors exert anisotropic force dipoles on the fibrous cytoskeleton. While the force dipoles can be either contractile or expansile, a medium made of fibers which buckle under compression rectifies these stresses towards a biologically crucial contraction. A general understanding of this rectification phenomenon as a function of the medium's elasticity is however lacking. Here we use theoretical continuum elasticity to show that rectification is actually a very general effect in nonlinear materials subjected to anisotropic internal stresses. We analytically show that both bucklable and constitutively linear materials subjected to geometrical nonlinearities rectify small forces towards contraction, while granular-like materials rectify towards expansion. Using simulations, we moreover show that these results extend to larger forces. Beyond fiber networks, these results could shed light on the propagation of stresses in brittle or granular materials following a local plastic rearrangement. △ Less

Submitted 11 June, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Main: 5 pages, 3 figures; SI: 8 pages, 5 figures

Valence can control the nonexponential viscoelastic relaxation of multivalent reversible gels

Authors: Hugo Le Roy, Jake Song, David Lundberg, leksandr V. Zhukhovitskiy, Jeremiah A. Johnson, Gareth H. McKinley, Niels Holten-Andersen, Martin Lenz

Abstract: Gels made of telechelic polymers connected by reversible crosslinkers are a versatile design platform for biocompatible viscoelastic materials. Their linear response to a step strain displays a fast, near-exponential relaxation when using low valence crosslinkers, while larger supramolecular crosslinkers bring about much slower dynamics involving a wide distribution of time scales whose physical o… ▽ More Gels made of telechelic polymers connected by reversible crosslinkers are a versatile design platform for biocompatible viscoelastic materials. Their linear response to a step strain displays a fast, near-exponential relaxation when using low valence crosslinkers, while larger supramolecular crosslinkers bring about much slower dynamics involving a wide distribution of time scales whose physical origin is still debated. Here, we propose a model where the relaxation of polymer gels in the dilute regime originates from elementary events in which the bonds connecting two neighboring crosslinkers all disconnect. Larger crosslinkers allow for a greater average number of bonds connecting them, but also generate more heterogeneity. We characterize the resulting distribution of relaxation time scales analytically, and accurately reproduce stress relaxation measurements on metal-coordinated hydrogels with a variety of crosslinker sizes including ions, metal-organic cages, and nanoparticles. Our approach is simple enough to be extended to any crosslinker size and could thus be harnessed for the rational design of complex viscoelastic materials. △ Less

Submitted 10 October, 2023; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 6 pages 5 figures 1 table for the main text and 9 pages 7 figures for the supplement

Creep properties and deformation mechanisms of single-crystalline $γ^\prime$-strengthened superalloys in dependence of the Co/Ni ratio

Authors: Nicklas Volz, Christopher H. Zenk, Nicolas Karpstein, Malte Lenz, Erdmann Spiecker, Mathias Göken, Steffen Neumeier

Abstract: Co-base superalloys are considered as promising high temperature materials besides the well-established Ni-base superalloys. However, Ni appears to be an indispensable alloying element also in Co-base superalloys. To address the influence of the base elements on the deformation behavior, high-temperature compressive creep experiments were performed on a single crystal alloy series that was designe… ▽ More Co-base superalloys are considered as promising high temperature materials besides the well-established Ni-base superalloys. However, Ni appears to be an indispensable alloying element also in Co-base superalloys. To address the influence of the base elements on the deformation behavior, high-temperature compressive creep experiments were performed on a single crystal alloy series that was designed to exhibit a varying Co/Ni ratio and a constant Al, W and Cr content. Creep tests were performed at 900 °C and 250 MPa and the resulting microstructures and defect configurations were characterized via electron microscopy. The minimum creep rates differ by more than one order of magnitude with changing Co/Ni ratio. An intermediate CoNi-base alloy exhibits the overall highest creep strength. Several strengthening contributions like solid solution strengthening of the $γ$ phase, effective diffusion coefficients or stacking fault energies were quantified. Precipitate shearing mechanisms differ significantly when the base element content is varied. While the Ni-rich superalloys exhibit SISF and SESF shearing, the Co-rich alloys develop extended APBs when the $γ^\prime$ phase is cut. This is mainly attributed to a difference in planar fault energies, caused by a changing segregation behavior. As result, it is assumed that the shearing resistivity and the occurring deformation mechanisms in the $γ^\prime$ phase are crucial for the creep properties of the investigated alloy series. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Image reconstruction in light-sheet microscopy: spatially varying deconvolution and mixed noise

Authors: Bogdan Toader, Jerome Boulanger, Yury Korolev, Martin O. Lenz, James Manton, Carola-Bibiane Schonlieb, Leila Muresan

Abstract: We study the problem of deconvolution for light-sheet microscopy, where the data is corrupted by spatially varying blur and a combination of Poisson and Gaussian noise. The spatial variation of the point spread function (PSF) of a light-sheet microscope is determined by the interaction between the excitation sheet and the detection objective PSF. First, we introduce a model of the image formation… ▽ More We study the problem of deconvolution for light-sheet microscopy, where the data is corrupted by spatially varying blur and a combination of Poisson and Gaussian noise. The spatial variation of the point spread function (PSF) of a light-sheet microscope is determined by the interaction between the excitation sheet and the detection objective PSF. First, we introduce a model of the image formation process that incorporates this interaction, therefore capturing the main characteristics of this imaging modality. Then, we formulate a variational model that accounts for the combination of Poisson and Gaussian noise through a data fidelity term consisting of the infimal convolution of the single noise fidelities, first introduced in L. Calatroni et al. "Infimal convolution of data discrepancies for mixed noise removal", SIAM Journal on Imaging Sciences 10.3 (2017), 1196-1233. We establish convergence rates in a Bregman distance under a source condition for the infimal convolution fidelity and a discrepancy principle for choosing the value of the regularisation parameter. The inverse problem is solved by applying the primal-dual hybrid gradient (PDHG) algorithm in a novel way. Finally, numerical experiments performed on both simulated and real data show superior reconstruction results in comparison with other methods. △ Less

Submitted 8 August, 2021; originally announced August 2021.

Comments: 34 pages, 13 figures

Twist-induced local curvature of filaments in DNA toroids

Authors: Luca Barberi, Martin Lenz

Abstract: DNA toroidal bundles form upon condensation of one or multiple DNA filaments. DNA filaments in toroidal bundles are hexagonally packed, and collectively twist around the center line of the toroid. In a previous study, we and our coworkers argue that the filaments' curvature locally correlates with their density in the bundle, with the filaments less closely packed where their curvature appears to… ▽ More DNA toroidal bundles form upon condensation of one or multiple DNA filaments. DNA filaments in toroidal bundles are hexagonally packed, and collectively twist around the center line of the toroid. In a previous study, we and our coworkers argue that the filaments' curvature locally correlates with their density in the bundle, with the filaments less closely packed where their curvature appears to be higher. We base our claim on the assumption that twist has a negligible effect on the local curvature of filaments in DNA toroids. However, this remains to be proven. We fill this gap here, by calculating the distribution of filaments' curvature in a geometric model of twisted toroidal bundle, which we use to describe DNA toroids by an appropriate choice of parameters. This allows us to substantiate our previous study and suggest directions for future experiments. △ Less

Submitted 18 January, 2021; originally announced January 2021.

Towards an Argument Mining Pipeline Transforming Texts to Argument Graphs

Authors: Mirko Lenz, Premtim Sahitaj, Sean Kallenberg, Christopher Coors, Lorik Dumani, Ralf Schenkel, Ralph Bergmann

Abstract: This paper targets the automated extraction of components of argumentative information and their relations from natural language text. Moreover, we address a current lack of systems to provide complete argumentative structure from arbitrary natural language text for general usage. We present an argument mining pipeline as a universally applicable approach for transforming German and English langua… ▽ More This paper targets the automated extraction of components of argumentative information and their relations from natural language text. Moreover, we address a current lack of systems to provide complete argumentative structure from arbitrary natural language text for general usage. We present an argument mining pipeline as a universally applicable approach for transforming German and English language texts to graph-based argument representations. We also introduce new methods for evaluating the results based on existing benchmark argument structures. Our results show that the generated argument graphs can be beneficial to detect new connections between different statements of an argumentative text. Our pipeline implementation is publicly available on GitHub. △ Less

Submitted 28 September, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Chiral active hexatics: Giant number fluctuations, waves and destruction of order

Authors: Ananyo Maitra, Martin Lenz, Raphael Voituriez

Abstract: Active materials, composed of internally driven particles, have properties that are qualitatively distinct from matter at thermal equilibrium. However, the most spectacular departures from equilibrium phase behaviour are thought to be confined to systems with polar or nematic asymmetry. In this paper, we show that such departures are also displayed in more symmetric phases such as hexatics if, in… ▽ More Active materials, composed of internally driven particles, have properties that are qualitatively distinct from matter at thermal equilibrium. However, the most spectacular departures from equilibrium phase behaviour are thought to be confined to systems with polar or nematic asymmetry. In this paper, we show that such departures are also displayed in more symmetric phases such as hexatics if, in addition, the constituent particles have chiral asymmetry. We show that chiral active hexatics whose rotation rate does not depend on density have giant number fluctuations. If the rotation rate depends on density, the giant number fluctuations are suppressed due to a novel orientation-density sound mode with a linear dispersion which propagates even in the overdamped limit. However, we demonstrate that beyond a finite but large lengthscale, a chirality and activity-induced relevant nonlinearity invalidates the predictions of the linear theory and destroys the hexatic order. In addition, we show that activity modifies the interactions between defects in the active chiral hexatic phase, making them non-mutual. Finally, to demonstrate the generality of a chiral active hexatic phase we show that it results from the melting of chiral active crystals in finite systems. △ Less

Submitted 28 December, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Design and Beam Test Results for the 2D Projective sPHENIX Electromagnetic Calorimeter Prototype

Authors: C. A. Aidala, S. Altaf, R. Belmont, S. Boose, D. Cacace, M. Connors, E. Desmond, J. Frantz, E. A. Gamez, N. Grau, J. S. Haggerty, A. Hodges, J. Huang, Y. Kim, M. D. Lenz, W. Lenz, N. A. Lewis, E. J. Mannel, J. D. Osborn, D. V. Perepelitsa, M. Phipps, R. Pisani, S. Polizzo, A. Pun, M. L. Purschke , et al. (13 additional authors not shown)

Abstract: sPHENIX is a new experiment under construction for the Relativistic Heavy Ion Collider at Brookhaven National Laboratory which will study the quark-gluon plasma to further the understanding of QCD matter and interactions. A prototype of the sPHENIX electromagnetic calorimeter (EMCal) was tested at the Fermilab Test Beam Facility in Spring 2018 as experiment T-1044. The EMCal prototype corresponds… ▽ More sPHENIX is a new experiment under construction for the Relativistic Heavy Ion Collider at Brookhaven National Laboratory which will study the quark-gluon plasma to further the understanding of QCD matter and interactions. A prototype of the sPHENIX electromagnetic calorimeter (EMCal) was tested at the Fermilab Test Beam Facility in Spring 2018 as experiment T-1044. The EMCal prototype corresponds to a solid angle of $ Δη\times Δφ= 0.2 \times 0.2$ centered at pseudo-rapidity $η= 1$. The prototype consists of scintillating fibers embedded in a mix of tungsten powder and epoxy. The fibers project back approximately to the center of the sPHENIX detector, giving 2D projectivity. The energy response of the EMCal prototype was studied as a function of position and input energy. The energy resolution of the EMCal prototype was obtained after applying a position dependent energy correction and a beam profile correction. Two separate position dependent corrections were considered. The EMCal energy resolution was found to be $σ(E)/\langle E\rangle = 3.5(0.1) \oplus 13.3(0.2)/\sqrt{E}$ based on the hodoscope position dependent correction, and $σ(E)/\langle E\rangle = 3.0(0.1) \oplus 15.4(0.3)/\sqrt{E}$ based on the cluster position dependent correction. These energy resolution results meet the requirements of the sPHENIX physics program. △ Less

Submitted 15 March, 2021; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 9 pages, 9 figures, 3 tables. Published in IEEE Transactions on Nuclear Science (vol. 68, no. 2, pp. 173-181, Feb. 2021). Revisions over previous version in response to referee and editor comments

Journal ref: IEEE Transactions on Nuclear Science, vol. 68, no. 2, pp. 173-181, Feb. 2021

Geometry of martensite needles in shape memory alloys

Authors: Sergio Conti, Martin Lenz, Nora Lüthen, Martin Rumpf, Barbara Zwicknagl

Abstract: We study the geometry of needle-shaped domains in shape-memory alloys. Needle-shaped domains are ubiquitously found in martensites around macroscopic interfaces between regions which are laminated in different directions, or close to macroscopic austenite/twinned-martensite interfaces. Their geometry results from the interplay of the local nonconvexity of the effective energy density with long-ran… ▽ More We study the geometry of needle-shaped domains in shape-memory alloys. Needle-shaped domains are ubiquitously found in martensites around macroscopic interfaces between regions which are laminated in different directions, or close to macroscopic austenite/twinned-martensite interfaces. Their geometry results from the interplay of the local nonconvexity of the effective energy density with long-range (linear) interactions mediated by the elastic strain field, and is up to now poorly understood. We present a two-dimensional shape optimization model based on finite elasticity and discuss its numerical solution. Our results indicate that the tapering profile of the needles can be understood within finite elasticity, but not with linearized elasticity. The resulting tapering and bending reproduce the main features of experimental observations on NiAl. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: 9 pages, 6 figures

MSC Class: 74N15 (Primary) 74P10; 74S05 (Secondary)

Parametrically Constrained Geometry Relaxations for High-Throughput Materials Science

Authors: Maja-Olivia Lenz, Thomas A. R. Purcell, David Hicks, Stefano Curtarolo, Matthias Scheffler, Christian Carbogno

Abstract: Reducing parameter spaces via exploiting symmetries has greatly accelerated and increased the quality of electronic-structure calculations. Unfortunately, many of the traditional methods fail when the global crystal symmetry is broken, even when the distortion is only a slight perturbation (e.g. Jahn-Teller like distortions). Here we introduce a flexible and generalizable parametric relaxation sch… ▽ More Reducing parameter spaces via exploiting symmetries has greatly accelerated and increased the quality of electronic-structure calculations. Unfortunately, many of the traditional methods fail when the global crystal symmetry is broken, even when the distortion is only a slight perturbation (e.g. Jahn-Teller like distortions). Here we introduce a flexible and generalizable parametric relaxation scheme, and implement it in the all-electron code FHI-aims. This approach utilizes parametric constraints to maintain symmetry at any level. After demonstrating the method's ability to relax metastable structures, we highlight its adaptability and performance over a test set of 359 materials, across thirteen lattice prototypes. Finally we show how these constraints can reduce the number of steps needed to relax local lattice distortions by an order of magnitude. The flexibility of these constraints enables a significant acceleration of the high-throughput searches for novel materials for numerous applications. △ Less

Submitted 30 October, 2019; v1 submitted 5 August, 2019; originally announced August 2019.

Swimmer suspensions on substrates: anomalous stability and long-range order

Authors: Ananyo Maitra, Pragya Srivastava, M. Cristina Marchetti, Sriram Ramaswamy, Martin Lenz

Abstract: We present a comprehensive theory of the dynamics and fluctuations of a two-dimensional suspension of polar active particles in an incompressible fluid confined to a substrate. We show that, depending on the sign of a single parameter, a state with polar orientational order is anomalously stable (or anomalously unstable), with a nonzero relaxation (or growth) rate for angular fluctuations at zero… ▽ More We present a comprehensive theory of the dynamics and fluctuations of a two-dimensional suspension of polar active particles in an incompressible fluid confined to a substrate. We show that, depending on the sign of a single parameter, a state with polar orientational order is anomalously stable (or anomalously unstable), with a nonzero relaxation (or growth) rate for angular fluctuations at zero wavenumber. This screening of the broken-symmetry mode in the stable state does lead to conventional rather than giant number fluctuations as argued by Bricard et al., Nature ${\bf 503}$, 95 (2013), but their bend instability in a splay-stable flock does not exist and the polar phase has long-range order in two dimensions. Our theory also describes confined three-dimensional thin-film suspensions of active polar particles as well as dense compressible active polar rods, and predicts a flocking transition without a banding instability △ Less

Submitted 4 January, 2019; originally announced January 2019.

Journal ref: Phys. Rev. Lett. 124, 028002 (2020)

Fiber plucking by molecular motors yields large emergent contractility in stiff biopolymer networks

Authors: Pierre Ronceray, Chase P. Broedersz, Martin Lenz

Abstract: The mechanical properties of the cell depend crucially on the tension of its cytoskeleton, a biopolymer network that is put under stress by active motor proteins. While the fibrous nature of the network is known to strongly affect the transmission of these forces to the cellular scale, our understanding of this process remains incomplete. Here we investigate the transmission of forces through the… ▽ More The mechanical properties of the cell depend crucially on the tension of its cytoskeleton, a biopolymer network that is put under stress by active motor proteins. While the fibrous nature of the network is known to strongly affect the transmission of these forces to the cellular scale, our understanding of this process remains incomplete. Here we investigate the transmission of forces through the network at the individual filament level, and show that active forces can be geometrically amplified as a transverse motor-generated force force "plucks" the fiber and induces a nonlinear tension. In stiff and densely connnected networks, this tension results in large network-wide tensile stresses that far exceed the expectation drawn from a linear elastic theory. This amplification mechanism competes with a recently characterized network-level amplification due to fiber buckling, suggesting that that fiber networks provide several distinct pathways for living systems to amplify their molecular forces. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Comments: 6 pages, 4 figures

Stress-dependent amplification of active forces in nonlinear elastic media

Authors: Pierre Ronceray, Chase Broedersz, Martin Lenz

Abstract: The production of mechanical stresses in living organisms largely relies on localized, force-generating active units embedded in filamentous matrices. Numerical simulations of discrete fiber networks with fixed boundaries have shown that buckling in the matrix dramatically amplifies the resulting active stresses. Here we extend this result to a bucklable continuum elastic medium subjected to an ar… ▽ More The production of mechanical stresses in living organisms largely relies on localized, force-generating active units embedded in filamentous matrices. Numerical simulations of discrete fiber networks with fixed boundaries have shown that buckling in the matrix dramatically amplifies the resulting active stresses. Here we extend this result to a bucklable continuum elastic medium subjected to an arbitrary external stress, and derive analytical expressions for the active, nonlinear constitutive relations characterizing the full active medium. Inserting these relations into popular "active gel" descriptions of living tissues and the cytoskeleton will enable investigations into nonlinear regimes previously inaccessible due to the phenomenological nature of these theories. △ Less

Submitted 7 May, 2018; originally announced May 2018.

Comments: 7 pages, 5 figures

Spontaneous rotation can stabilise ordered chiral active fluids

Authors: Ananyo Maitra, Martin Lenz

Abstract: Active hydrodynamic theories are a powerful tool to study the emergent ordered phases of internally driven particles such as bird flocks, bacterial suspension and their artificial analogues. While theories of orientationally ordered phases are by now well established, the effect of chirality on these phases is much less studied. In this paper, we present the first complete dynamical theory of orie… ▽ More Active hydrodynamic theories are a powerful tool to study the emergent ordered phases of internally driven particles such as bird flocks, bacterial suspension and their artificial analogues. While theories of orientationally ordered phases are by now well established, the effect of chirality on these phases is much less studied. In this paper, we present the first complete dynamical theory of orientationally ordered chiral particles in two-dimensional incompressible systems. We show that phase-coherent states of rotating chiral particles are remarkably stable in both momentum-conserved and non-conserved systems in contrast to their non-rotating counterparts. Furthermore, defect separation -- which drives chaotic flows in non-rotating active fluids -- is suppressed by intrinsic rotation of chiral active particles. We thus establish chirality as a source of dramatic stabilization in active systems, which could be key in interpreting the collective behaviours of some biological tissues, cytoskeletal systems and collections of bacteria. △ Less

Submitted 4 January, 2019; v1 submitted 26 April, 2018; originally announced April 2018.

Homogenization in magnetic-shape-memory polymer composites

Authors: Sergio Conti, Martin Lenz, Matthäus Pawelczyk, Martin Rumpf

Abstract: Magnetic-shape-memory materials (e.g. specific NiMnGa alloys) react with a large change of shape to the presence of an external magnetic field. As an alternative for the difficult to manifacture single crystal of these alloys we study composite materials in which small magnetic-shape-memory particles are embedded in a polymer matrix. The macroscopic properties of the composite depend strongly on t… ▽ More Magnetic-shape-memory materials (e.g. specific NiMnGa alloys) react with a large change of shape to the presence of an external magnetic field. As an alternative for the difficult to manifacture single crystal of these alloys we study composite materials in which small magnetic-shape-memory particles are embedded in a polymer matrix. The macroscopic properties of the composite depend strongly on the geometry of the microstructure and on the characteristics of the particles and the polymer. We present a variational model based on micromagnetism and elasticity, and derive via homogenization an effective macroscopic model under the assumption that the microstructure is periodic. We then study numerically the resulting cell problem, and discuss the effect of the microstructure on the macroscopic material behavior. Our results may be used to optimize the shape of the particles and the microstructure. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 17 pages, 4 figures

MSC Class: 74Q05; 74N15; 74F15; 74S15

Stability from activity

Authors: Ananyo Maitra, Pragya Srivastava, M. Cristina Marchetti, Juho Lintuvuori, Sriram Ramaswamy, Martin Lenz

Abstract: Suspensions of actively driven anisotropic objects exhibit distinctively nonequilibrium behaviors, and current theories predict that they are incapable of sustaining orientational order at high activity. By contrast, here we show that nematic suspensions on a substrate can display order at arbitrarily high activity due to a previously unreported, potentially stabilizing active force. The resulting… ▽ More Suspensions of actively driven anisotropic objects exhibit distinctively nonequilibrium behaviors, and current theories predict that they are incapable of sustaining orientational order at high activity. By contrast, here we show that nematic suspensions on a substrate can display order at arbitrarily high activity due to a previously unreported, potentially stabilizing active force. The resulting nonequilibrium ordered phase displays robust giant number fluctuations that cannot be suppressed even by an incompressible solvent. Our results apply to virtually all experimental assays used to investigate the active nematic ordering of self-propelled colloids, bacterial suspensions and the cytoskeleton, and have testable implications in interpreting their nonequilibrium behaviors △ Less

Submitted 7 November, 2017; originally announced November 2017.

Extensile actomyosin?

Authors: Martin Lenz

Abstract: Living cells move thanks to assemblies of actin filaments and myosin motors that range from very organized striated muscle tissue to disordered intracellular bundles. The mechanisms powering these disordered structures are debated, and all models studied so far predict that they are contractile. We reexamine this prediction through a theoretical treatment of the interplay of three well-characteriz… ▽ More Living cells move thanks to assemblies of actin filaments and myosin motors that range from very organized striated muscle tissue to disordered intracellular bundles. The mechanisms powering these disordered structures are debated, and all models studied so far predict that they are contractile. We reexamine this prediction through a theoretical treatment of the interplay of three well-characterized internal dynamical processes in actomyosin bundles: actin treadmilling, the attachement-detachment dynamics of myosin and that of crosslinking proteins. We show that these processes enable an extensive control of the bundle's active mechanics, including reversals of the filaments' apparent velocities and the possibility of generating extension instead of contraction. These effects offer a new perspective on well-studied in vivo systems, as well as a robust criterion to experimentally elucidate the underpinnings of actomyosin activity. △ Less

Submitted 24 September, 2017; originally announced September 2017.

Comments: 5 pages, 5 figures

Stanley-Reisner rings for quasi-arithmetic matroids

Authors: Matthias Lenz

Abstract: In this note we define a Stanley-Reisner ring for quasi-arithmetic matroids and more general structures. To this end, we define two types of CW complexes associated with a quasi-arithmetic matroid that generalize independence complexes of matroids. Then we use Stanley's construction of Stanley-Reisner rings for simplicial posets. In this note we define a Stanley-Reisner ring for quasi-arithmetic matroids and more general structures. To this end, we define two types of CW complexes associated with a quasi-arithmetic matroid that generalize independence complexes of matroids. Then we use Stanley's construction of Stanley-Reisner rings for simplicial posets. △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: 14 pages, 4 figures

MSC Class: Primary: 05B35; 13F55; Secondary: 06A07; 06A11; 14N20

Cell contraction induces long-ranged stress stiffening in the extracellular matrix

Authors: Yu Long Han, Pierre Ronceray, Guoqiang Xu, Andrea Malandrino, Roger Kamm, Martin Lenz, Chase P. Broedersz, Ming Guo

Abstract: Animal cells in tissues are supported by biopolymer matrices, which typically exhibit highly nonlinear mechanical properties. While the linear elasticity of the matrix can significantly impact cell mechanics and functionality, it remains largely unknown how cells, in turn, affect the nonlinear mechanics of their surrounding matrix. Here we show that living contractile cells are able to generate a… ▽ More Animal cells in tissues are supported by biopolymer matrices, which typically exhibit highly nonlinear mechanical properties. While the linear elasticity of the matrix can significantly impact cell mechanics and functionality, it remains largely unknown how cells, in turn, affect the nonlinear mechanics of their surrounding matrix. Here we show that living contractile cells are able to generate a massive stiffness gradient in three distinct 3D extracellular matrix model systems: collagen, fibrin, and Matrigel. We decipher this remarkable behavior by introducing Nonlinear Stress Inference Microscopy (NSIM), a novel technique to infer stress fields in a 3D matrix from nonlinear microrheology measurement with optical tweezers. Using NSIM and simulations, we reveal a long-ranged propagation of cell-generated stresses resulting from local filament buckling. This slow decay of stress gives rise to the large spatial extent of the observed cell-induced matrix stiffness gradient, which could form a mechanism for mechanical communication between cells. △ Less

Submitted 3 September, 2017; originally announced September 2017.

Computing the poset of layers of a toric arrangement

Authors: Matthias Lenz

Abstract: A toric arrangement is an arrangement of subtori of codimension one in a real or complex torus. The poset of layers is the set of connected components of non-empty intersections of these subtori, partially ordered by reverse inclusion. In this note we present an algorithm that computes this poset in the central case. A toric arrangement is an arrangement of subtori of codimension one in a real or complex torus. The poset of layers is the set of connected components of non-empty intersections of these subtori, partially ordered by reverse inclusion. In this note we present an algorithm that computes this poset in the central case. △ Less

Submitted 24 August, 2017; v1 submitted 22 August, 2017; originally announced August 2017.

Comments: 10 pages, 5 figures, minor corrections (a false statement in the introduction was corrected)

MSC Class: 05B35; 06A07; 06A11; 14N20; 52C35

Geometrical frustration yields fiber formation in self-assembly

Authors: Martin Lenz, Thomas A. Witten

Abstract: Controlling the self-assembly of supramolecular structures is vital for living cells, and a central challenge for engineering at the nano- and microscales. Nevertheless, even particles without optimized shapes can robustly form well-defined morphologies. This is the case in numerous medical conditions where normally soluble proteins aggregate into fibers. Beyond the diversity of molecular mechanis… ▽ More Controlling the self-assembly of supramolecular structures is vital for living cells, and a central challenge for engineering at the nano- and microscales. Nevertheless, even particles without optimized shapes can robustly form well-defined morphologies. This is the case in numerous medical conditions where normally soluble proteins aggregate into fibers. Beyond the diversity of molecular mechanisms involved, we propose that fibers generically arise from the aggregation of irregular particles with short-range interactions. Using a minimal model of ill-fitting, sticky particles, we demonstrate robust fiber formation for a variety of particle shapes and aggregation conditions. Geometrical frustration plays a crucial role in this process, and accounts for the range of parameters in which fibers form as well as for their metastable character. △ Less

Submitted 23 May, 2017; originally announced May 2017.

Comments: 6 pages, 5 figures, 1 ancillary movie; to appear in Nature Physics

Journal ref: Nat. Phys. 13, 1100 (2017)

Representations of weakly multiplicative arithmetic matroids are unique

Authors: Matthias Lenz

Abstract: An arithmetic matroid is weakly multiplicative if the multiplicity of at least one of its bases is equal to the product of the multiplicities of its elements. We show that if such an arithmetic matroid can be represented by an integer matrix, then this matrix is uniquely determined. This implies that the integer cohomology ring of a centred toric arrangement whose arithmetic matroid is weakly mult… ▽ More An arithmetic matroid is weakly multiplicative if the multiplicity of at least one of its bases is equal to the product of the multiplicities of its elements. We show that if such an arithmetic matroid can be represented by an integer matrix, then this matrix is uniquely determined. This implies that the integer cohomology ring of a centred toric arrangement whose arithmetic matroid is weakly multiplicative is determined by its poset of layers. This partially answers a question asked by Callegaro-Delucchi. △ Less

Submitted 4 March, 2018; v1 submitted 27 April, 2017; originally announced April 2017.

Comments: 9 pages, 2 figures, minor corrections

MSC Class: Primary: 05B35; 52C35; Secondary: 14M15

Journal ref: Annals of Combinatorics 23 (2019), no. 2, 335-346

Design and Beam Test Results for the sPHENIX Electromagnetic and Hadronic Calorimeter Prototypes

Authors: C. A. Aidala, V. Bailey, S. Beckman, R. Belmont, C. Biggs, J. Blackburn, S. Boose, M. Chiu, M. Connors, E. Desmond, A. Franz, J. S. Haggerty, X. He, M. M. Higdon, J. Huang, K. Kauder, E. Kistenev, J. LaBounty, J. G. Lajoie, M. Lenz, W. Lenz, S. Li, V. R. Loggins, E. J. Mannel, T. Majoros , et al. (25 additional authors not shown)

Abstract: The super Pioneering High Energy Nuclear Interaction eXperiment (sPHENIX) at the Relativistic Heavy Ion Collider (RHIC) will perform high precision measurements of jets and heavy flavor observables for a wide selection of nuclear collision systems, elucidating the microscopic nature of strongly interacting matter ranging from nucleons to the strongly coupled quark-gluon plasma. A prototype of the… ▽ More The super Pioneering High Energy Nuclear Interaction eXperiment (sPHENIX) at the Relativistic Heavy Ion Collider (RHIC) will perform high precision measurements of jets and heavy flavor observables for a wide selection of nuclear collision systems, elucidating the microscopic nature of strongly interacting matter ranging from nucleons to the strongly coupled quark-gluon plasma. A prototype of the sPHENIX calorimeter system was tested at the Fermilab Test Beam Facility as experiment T-1044 in the spring of 2016. The electromagnetic calorimeter (EMCal) prototype is composed of scintillating fibers embedded in a mixture of tungsten powder and epoxy. The hadronic calorimeter (HCal) prototype is composed of tilted steel plates alternating with plastic scintillator. Results of the test beam reveal the energy resolution for electrons in the EMCal is $2.8\%\oplus~15.5\%/\sqrt{E}$ and the energy resolution for hadrons in the combined EMCal plus HCal system is $13.5\%\oplus 64.9\%/\sqrt{E}$. These results demonstrate that the performance of the proposed calorimeter system satisfies the sPHENIX specifications. △ Less

Submitted 16 December, 2018; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: 19 pages, 31 figures, 3 tables. Published in IEEE Transactions on Nuclear Science (Volume: 65 , Issue: 12 , Dec. 2018 ). Revisions over previous version in response to referee and editor comments, plus updating of references and acknowledgement

Journal ref: IEEE Transactions on Nuclear Science, Volume 65, Issue 12, pp. 2901-2919, December 2018

On powers of Plücker coordinates and representability of arithmetic matroids

Authors: Matthias Lenz

Abstract: The first problem we investigate is the following: given $k\in \mathbb{R}_{\ge 0}$ and a vector $v$ of Plücker coordinates of a point in the real Grassmannian, is the vector obtained by taking the $k$th power of each entry of $v$ again a vector of Plücker coordinates? For $k\neq 1$, this is true if and only if the corresponding matroid is regular. Similar results hold over other fields. We also de… ▽ More The first problem we investigate is the following: given $k\in \mathbb{R}_{\ge 0}$ and a vector $v$ of Plücker coordinates of a point in the real Grassmannian, is the vector obtained by taking the $k$th power of each entry of $v$ again a vector of Plücker coordinates? For $k\neq 1$, this is true if and only if the corresponding matroid is regular. Similar results hold over other fields. We also describe the subvariety of the Grassmannian that consists of all the points that define a regular matroid. The second topic is a related problem for arithmetic matroids. Let $\mathcal{A} = (E, rk, m)$ be an arithmetic matroid and let $k\neq 1 $ be a non-negative integer. We prove that if $\mathcal{A}$ is representable and the underlying matroid is non-regular, then $\mathcal{A}^k := (E, rk, m^k)$ is not representable. This provides a large class of examples of arithmetic matroids that are not representable. On the other hand, if the underlying matroid is regular and an additional condition is satisfied, then $\mathcal{A}^k$ is representable. Bajo-Burdick-Chmutov have recently discovered that arithmetic matroids of type $\mathcal{A}^2$ arise naturally in the study of colourings and flows on CW complexes. In the last section, we prove a family of necessary conditions for representability of arithmetic matroids. △ Less

Submitted 29 September, 2019; v1 submitted 30 March, 2017; originally announced March 2017.

Comments: 36 pages, 1 figure, minor corrections, same content as journal version

MSC Class: Primary: 05B35; 14M15; Secondary: 14T05

The dynamics of filament assembly define cytoskeletal network morphology

Authors: Giulia Foffano, Nicolas Levernier, Martin Lenz

Abstract: The actin cytoskeleton is a key component in the machinery of eukaryotic cells, and it selfassembles out of equilibrium into a wide variety of biologically crucial structures. While the molecular mechanisms involved are well characterized, the physical principles governing the spatial arrangement of actin filaments are not understood. Here we propose that the dynamics of actin network assembly fro… ▽ More The actin cytoskeleton is a key component in the machinery of eukaryotic cells, and it selfassembles out of equilibrium into a wide variety of biologically crucial structures. While the molecular mechanisms involved are well characterized, the physical principles governing the spatial arrangement of actin filaments are not understood. Here we propose that the dynamics of actin network assembly from growing filaments results from a competition between diffusion, bundling, and steric hindrance, and is responsible for the range of observed morphologies. Our model and simulations thus predict an abrupt dynamical transition between homogeneous and strongly bundled networks as a function of the actin polymerization rate. This suggests that cells may effect dramatic changes to their internal architecture through minute modifications of their nonequilibrium dynamics. Our results are consistent with available experimental data. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 9 pages, 7 figures. Submitted to and accepted by Nat. Commun

DETOx: Towards Optimal Software-based Soft-Error Detector Configurations

Authors: Michael Lenz, Horst Schirmeier

Abstract: Application developers often place executable assertions -- equipped with program-specific predicates -- in their system, targeting programming errors. However, these detectors can detect data errors resulting from transient hardware faults in main memory as well. But while an assertion reduces silent data corruptions (SDCs) in the program state they check, they add runtime to the target program t… ▽ More Application developers often place executable assertions -- equipped with program-specific predicates -- in their system, targeting programming errors. However, these detectors can detect data errors resulting from transient hardware faults in main memory as well. But while an assertion reduces silent data corruptions (SDCs) in the program state they check, they add runtime to the target program that increases the attack surface for the remaining state. This article outlines an approach to find an optimal subset of assertions that minimizes the SDC count, without the need to run fault-injection experiments for every possible assertion subset. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: Editor: Gilles Tredan. 12th European Dependable Computing Conference (EDCC 2016), September 5-9, 2016, Gothenburg, Sweden. Fast Abstracts Proceedings -- EDCC 2016

A convolution formula for Tutte polynomials of arithmetic matroids and other combinatorial structures

Authors: Spencer Backman, Matthias Lenz

Abstract: In this note we generalize the convolution formula for the Tutte polynomial of Kook-Reiner-Stanton and Etienne-Las Vergnas to a more general setting that includes both arithmetic matroids and delta-matroids. As corollaries, we obtain new proofs of two positivity results for pseudo-arithmetic matroids and a combinatorial interpretation of the arithmetic Tutte polynomial at infinitely many points in… ▽ More In this note we generalize the convolution formula for the Tutte polynomial of Kook-Reiner-Stanton and Etienne-Las Vergnas to a more general setting that includes both arithmetic matroids and delta-matroids. As corollaries, we obtain new proofs of two positivity results for pseudo-arithmetic matroids and a combinatorial interpretation of the arithmetic Tutte polynomial at infinitely many points in terms of arithmetic flows and colorings. We also exhibit connections with a decomposition of Dahmen-Micchelli spaces and lattice point counting in zonotopes. △ Less

Submitted 21 April, 2017; v1 submitted 8 February, 2016; originally announced February 2016.

Comments: 13 pages, minor corrections, a section on the background was added

MSC Class: Primary: 05B35; Secondary: 05C15; 05C21; 05C31; 05C10; 52C35

A posteriori modeling error estimates in the optimization of two-scale elastic composite materials

Authors: Sergio Conti, Benedict Geihe, Martin Lenz, Martin Rumpf

Abstract: The a posteriori analysis of the discretization error and the modeling error is studied for a compliance cost functional in the context of the optimization of composite elastic materials and a two-scale linearized elasticity model. A mechanically simple, parametrized microscopic supporting structure is chosen and the parameters describing the structure are determined minimizing the compliance obje… ▽ More The a posteriori analysis of the discretization error and the modeling error is studied for a compliance cost functional in the context of the optimization of composite elastic materials and a two-scale linearized elasticity model. A mechanically simple, parametrized microscopic supporting structure is chosen and the parameters describing the structure are determined minimizing the compliance objective. An a posteriori error estimate is derived which includes the modeling error caused by the replacement of a nested laminate microstructure by this considerably simpler microstructure. Indeed, nested laminates are known to realize the minimal compliance and provide a benchmark for the quality of the microstructures. To estimate the local difference in the compliance functional the dual weighted residual approach is used. Different numerical experiments show that the resulting adaptive scheme leads to simple parametrized microscopic supporting structures that can compete with the optimal nested laminate construction. The derived a posteriori error indicators allow to verify that the suggested simplified microstructures achieve the optimal value of the compliance up to a few percent. Furthermore, it is shown how discretization error and modeling error can be balanced by choosing an optimal level of grid refinement. Our two scale results with a single scale microstructure can provide guidance towards the design of a producible macroscopic fine scale pattern. △ Less

Submitted 11 March, 2017; v1 submitted 13 November, 2015; originally announced November 2015.

MSC Class: 35B27; 49M29; 49Q1; 65N30; 65N38; 74P05

Fiber networks amplify active stress

Authors: Pierre Ronceray, Chase Broedersz, Martin Lenz

Abstract: Large-scale force generation is essential for biological functions such as cell motility, embryonic development, and muscle contraction. In these processes, forces generated at the molecular level by motor proteins are transmitted by disordered fiber networks, resulting in large-scale active stresses. While these fiber networks are well characterized macroscopically, this stress generation by micr… ▽ More Large-scale force generation is essential for biological functions such as cell motility, embryonic development, and muscle contraction. In these processes, forces generated at the molecular level by motor proteins are transmitted by disordered fiber networks, resulting in large-scale active stresses. While these fiber networks are well characterized macroscopically, this stress generation by microscopic active units is not well understood. Here we theoretically study force transmission in these networks, and find that local active forces are rectified towards isotropic contraction and strongly amplified as fibers collectively buckle in the vicinity of the active units. This stress amplification is reinforced by the networks' disordered nature, but saturates for high densities of active units. Our predictions are quantitatively consistent with experiments on reconstituted tissues and actomyosin networks, and shed light on the role of the network microstructure in shaping active stresses in cells and tissue. △ Less

Submitted 22 July, 2015; v1 submitted 21 July, 2015; originally announced July 2015.

Comments: 8 pages, 4 figures. Supporting information: 5 pages, 5 figures

Journal ref: Proc. Natl. Acad. Sci. U.S.A. 113, 2827 (2016)

A Reaction-Diffusion Model of the Cadherin-Catenin System: A Possible Mechanism for Contact Inhibition and Implications for Tumorigenesis

Authors: Markus Basan, Timon Idema, Martin Lenz, Jean-François Joanny, Thomas Risler

Abstract: Contact inhibition is the process by which cells switch from a motile growing state to a passive and stabilized state upon touching their neighbors. When two cells touch, an adhesion link is created between them by means of transmembrane E-cadherin proteins. Simultaneously, their actin filaments stop polymerizing in the direction perpendicular to the membrane and reorganize to create an apical bel… ▽ More Contact inhibition is the process by which cells switch from a motile growing state to a passive and stabilized state upon touching their neighbors. When two cells touch, an adhesion link is created between them by means of transmembrane E-cadherin proteins. Simultaneously, their actin filaments stop polymerizing in the direction perpendicular to the membrane and reorganize to create an apical belt that colocalizes with the adhesion links. Here, we propose a detailed quantitative model of the role of the cytoplasmic $β$-catenin and $α$-catenin proteins in this process, treated as a reaction-diffusion system. Upon cell-cell contact, the concentration in $α$-catenin dimers increases, inhibiting actin branching and thereby reducing cellular motility and expansion pressure. This model provides a mechanism for contact inhibition that could explain previously unrelated experimental findings on the role played by E-cadherin, $β$-catenin and $α$-catenin in the cellular phenotype and in tumorigenesis. In particular, we address the effect of a knockout of the adenomatous polyposis coli tumor suppressor gene. Potential direct tests of our model are discussed. △ Less

Submitted 4 May, 2015; originally announced May 2015.

Journal ref: Biophys. J. 98 (12), 2770-2779 (2010)

Hysteresis in Magnetic Shape Memory Composites: Modeling and Simulation

Authors: Sergio Conti, Martin Lenz, Martin Rumpf

Abstract: Magnetic shape memory alloys are characterized by the coupling between a structural phase transition and magnetic one. This permits to control the shape change via an external magnetic field, at least in single crystals. Composite materials with single-crystalline particles embedded in a softer matrix have been proposed as a way to overcome the blocking of the transformation at grain boundaries. W… ▽ More Magnetic shape memory alloys are characterized by the coupling between a structural phase transition and magnetic one. This permits to control the shape change via an external magnetic field, at least in single crystals. Composite materials with single-crystalline particles embedded in a softer matrix have been proposed as a way to overcome the blocking of the transformation at grain boundaries. We investigate hysteresis phenomena for small NiMnGa single crystals embedded in a polymer matrix for slowly varying magnetic fields. The evolution of the microstructure is studied within the rate-independent variational framework proposed by Mielke and Theil (1999). The underlying variational model incorporates linearized elasticity, micromagnetism, stray field and a dissipation term proportional to the volume swept by the phase boundary. The time discretization is based on an incremental minimization of the sum of energy and dissipation. A backtracking approach is employed to approximately ensure the global minimality condition. We illustrate and discuss the influence of the particle geometry (volume fraction, shape, arrangement) and the polymer elastic parameters on the observed hysteresis and compare with recent experimental results. △ Less

Submitted 23 February, 2015; v1 submitted 19 February, 2015; originally announced February 2015.

Comments: 18 pages, 11 figures

MSC Class: 74N30 (Primary); 74F15; 74S15 (Secondary)

Journal ref: J. Mech. Phys. Solids 89 (2016), 272-286

Connecting local active forces to macroscopic stress in elastic media

Authors: Pierre Ronceray, Martin Lenz

Abstract: In contrast with ordinary materials, living matter drives its own motion by generating active, out-of-equilibrium internal stresses. These stresses typically originate from localized active elements embedded in an elastic medium, such as molecular motors inside the cell or contractile cells in a tissue. While many large-scale phenomenological theories of such active media have been developed, a sy… ▽ More In contrast with ordinary materials, living matter drives its own motion by generating active, out-of-equilibrium internal stresses. These stresses typically originate from localized active elements embedded in an elastic medium, such as molecular motors inside the cell or contractile cells in a tissue. While many large-scale phenomenological theories of such active media have been developed, a systematic understanding of the emergence of stress from the local force-generating elements is lacking. In this paper, we present a rigorous theoretical framework to study this relationship. We show that the medium's macroscopic active stress tensor is equal to the active elements' force dipole tensor per unit volume in both continuum and discrete linear homogeneous media of arbitrary geometries. This relationship is conserved on average in the presence of disorder, but can be violated in nonlinear elastic media. Such effects can lead to either a reinforcement or an attenuation of the active stresses, giving us a glimpse of the ways in which nature might harness microscopic forces to create active materials. △ Less

Submitted 12 November, 2014; originally announced November 2014.

Comments: 9 pages, 4 figures

Splines, lattice points, and arithmetic matroids

Authors: Matthias Lenz

Abstract: Let $X$ be a $(d\times N)$-matrix. We consider the variable polytope $Π_X(u) = \{w \ge 0 : X w = u \}$. It is known that the function $T_X$ that assigns to a parameter $u \in \mathbb{R}^d$ the volume of the polytope $Π_X(u)$ is piecewise polynomial. The Brion-Vergne formula implies that the number of lattice points in $Π_X(u)$ can be obtained by applying a certain differential operator to the func… ▽ More Let $X$ be a $(d\times N)$-matrix. We consider the variable polytope $Π_X(u) = \{w \ge 0 : X w = u \}$. It is known that the function $T_X$ that assigns to a parameter $u \in \mathbb{R}^d$ the volume of the polytope $Π_X(u)$ is piecewise polynomial. The Brion-Vergne formula implies that the number of lattice points in $Π_X(u)$ can be obtained by applying a certain differential operator to the function $T_X$. In this article we slightly improve the Brion-Vergne formula and we study two spaces of differential operators that arise in this context: the space of relevant differential operators (i.e. operators that do not annihilate $T_X$) and the space of nice differential operators (i.e. operators that leave $T_X$ continuous). These two spaces are finite-dimensional homogeneous vector spaces and their Hilbert series are evaluations of the Tutte polynomial of the arithmetic matroid defined by the matrix $X$. They are closely related to the $\mathcal{P}$-spaces studied by Ardila-Postnikov and Holtz-Ron in the context of zonotopal algebra and power ideals. △ Less

Submitted 15 September, 2014; v1 submitted 18 August, 2014; originally announced August 2014.

Comments: 43 pages, 4 figures, minor corrections

MSC Class: Primary: 05B35; 19L10; 52B20; Secondary: 13B25; 14M25; 16S32; 41A15; 47F05; 52B40; 52C35

Journal ref: Journal of Algebraic Combinatorics, March 2016, Volume 43, Issue 2, pp 277-324

Geometrical origins of contractility in disordered actomyosin networks

Authors: Martin Lenz

Abstract: Movement within eukaryotic cells largely originates from localized forces exerted by myosin motors on scaffolds of actin filaments. Although individual motors locally exert both contractile and extensile forces, large actomyosin structures at the cellular scale are overwhelmingly contractile, suggesting that the scaffold serves to favor contraction over extension. While this mechanism is well unde… ▽ More Movement within eukaryotic cells largely originates from localized forces exerted by myosin motors on scaffolds of actin filaments. Although individual motors locally exert both contractile and extensile forces, large actomyosin structures at the cellular scale are overwhelmingly contractile, suggesting that the scaffold serves to favor contraction over extension. While this mechanism is well understood in highly organized striated muscle, its origin in disordered networks such as the cell cortex is unknown. Here we develop a mathematical model of the actin scaffold's local two- or three-dimensional mechanics and identify four competing contraction mechanisms. We predict that one mechanism dominates, whereby local deformations of the actin break the balance between contraction and extension. In this mechanism, contractile forces result mostly from motors plucking the filaments transversely rather than buckling them longitudinally. These findings sheds light on recent $\textit{in vitro}$ experiments, and provides a new geometrical understanding of contractility in the myriad of disordered actomyosin systems found $\textit{in vivo}$. △ Less

Submitted 24 July, 2014; originally announced July 2014.

Comments: 9 pages, 5 figures, 5-page supplement

Journal ref: Phys. Rev. X 4, 041002 (2014)

The PHENIX Forward Silicon Vertex Detector

Authors: C. Aidala, L. Anaya, E. Anderssen, A. Bambaugh, A. Barron, J. G. Boissevain, J. Bok, S. Boose, M. L. Brooks, S. Butsyk, M. Cepeda, P. Chacon, S. Chacon, L. Chavez, T. Cote, C. D'Agostino, A. Datta, K. DeBlasio, L. DelMonte, E. J. Desmond, J. M. Durham, D. Fields, M. Finger, C. Gingu, B. Gonzales , et al. (60 additional authors not shown)

Abstract: A new silicon detector has been developed to provide the PHENIX experiment with precise charged particle tracking at forward and backward rapidity. The Forward Silicon Vertex Tracker (FVTX) was installed in PHENIX prior to the 2012 run period of the Relativistic Heavy Ion Collider (RHIC). The FVTX is composed of two annular endcaps, each with four stations of silicon mini-strip sensors, covering a… ▽ More A new silicon detector has been developed to provide the PHENIX experiment with precise charged particle tracking at forward and backward rapidity. The Forward Silicon Vertex Tracker (FVTX) was installed in PHENIX prior to the 2012 run period of the Relativistic Heavy Ion Collider (RHIC). The FVTX is composed of two annular endcaps, each with four stations of silicon mini-strip sensors, covering a rapidity range of $1.2<|η|<2.2$ that closely matches the two existing PHENIX muon arms. Each station consists of 48 individual silicon sensors, each of which contains two columns of mini-strips with 75 $μ$m pitch in the radial direction and lengths in the $φ$ direction varying from 3.4 mm at the inner radius to 11.5 mm at the outer radius. The FVTX has approximately 0.54 million strips in each endcap. These are read out with FPHX chips, developed in collaboration with Fermilab, which are wire bonded directly to the mini-strips. The maximum strip occupancy reached in central Au-Au collisions is approximately 2.8%. The precision tracking provided by this device makes the identification of muons from secondary vertices away from the primary event vertex possible. The expected distance of closest approach (DCA) resolution of 200 $μ$m or better for particles with a transverse momentum of 5 GeV/$c$ will allow identification of muons from relatively long-lived particles, such as $D$ and $B$ mesons, through their broader DCA distributions. △ Less

Submitted 14 February, 2014; v1 submitted 14 November, 2013; originally announced November 2013.

Comments: version submitted to NIM

Lattice points in polytopes, box splines, and Todd operators

Authors: Matthias Lenz

Abstract: Let $X$ be a list of vectors that is totally unimodular. In a previous article the author proved that every real-valued function on the set of interior lattice points of the zonotope defined by $X$ can be extended to a function on the whole zonotope of the form $p(D)B_X$ in a unique way, where $p(D)$ is a differential operator that is contained in the so-called internal $\Pcal$-space. In this pape… ▽ More Let $X$ be a list of vectors that is totally unimodular. In a previous article the author proved that every real-valued function on the set of interior lattice points of the zonotope defined by $X$ can be extended to a function on the whole zonotope of the form $p(D)B_X$ in a unique way, where $p(D)$ is a differential operator that is contained in the so-called internal $\Pcal$-space. In this paper we construct an explicit solution to this interpolation problem in terms of Todd operators. As a corollary we obtain a slight generalisation of the Khovanskii-Pukhlikov formula that relates the volume and the number of integer points in a smooth lattice polytope. △ Less

Submitted 13 May, 2013; originally announced May 2013.

Comments: 15 pages, 4 figures

MSC Class: Primary: 05B35; 19L10; 52B20; 52B40; Secondary: 13B25; 14M25; 16S32; 41A15; 47F05

Journal ref: International Mathematics Research Notices, Volume 2015, Issue 14, 2015, Pages 5289-5310

Interpolation, box splines, and lattice points in zonotopes

Authors: Matthias Lenz

Abstract: Let $X$ be a totally unimodular list of vectors in some lattice. Let $B_X$ be the box spline defined by $X$. Its support is the zonotope $Z(X)$. We show that any real-valued function defined on the set of lattice points in the interior of $Z(X)$ can be extended to a function on $Z(X)$ of the form $p(D)B_X$ in a unique way, where $p(D)$ is a differential operator that is contained in the so-called… ▽ More Let $X$ be a totally unimodular list of vectors in some lattice. Let $B_X$ be the box spline defined by $X$. Its support is the zonotope $Z(X)$. We show that any real-valued function defined on the set of lattice points in the interior of $Z(X)$ can be extended to a function on $Z(X)$ of the form $p(D)B_X$ in a unique way, where $p(D)$ is a differential operator that is contained in the so-called internal $\Pcal$-space. This was conjectured by Olga Holtz and Amos Ron. We also point out connections between this interpolation problem and matroid theory, including a deletion-contraction decomposition. △ Less

Submitted 6 November, 2012; originally announced November 2012.

Comments: 10 pages, 3 figures

MSC Class: Primary: 05B35; 41A05; 41A15; 52B20; Secondary: 13F20; 41A63; 47F05; 52B40; 52C07

Journal ref: International Mathematics Research Notices, Volume 2014, Issue 20, 2014, Pages 5697-5712

Zonotopal algebra and forward exchange matroids

Authors: Matthias Lenz

Abstract: Zonotopal algebra is the study of a family of pairs of dual vector spaces of multivariate polynomials that can be associated with a list of vectors X. It connects objects from combinatorics, geometry, and approximation theory. The origin of zonotopal algebra is the pair (D(X),P(X)), where D(X) denotes the Dahmen-Micchelli space that is spanned by the local pieces of the box spline and P(X) is a sp… ▽ More Zonotopal algebra is the study of a family of pairs of dual vector spaces of multivariate polynomials that can be associated with a list of vectors X. It connects objects from combinatorics, geometry, and approximation theory. The origin of zonotopal algebra is the pair (D(X),P(X)), where D(X) denotes the Dahmen-Micchelli space that is spanned by the local pieces of the box spline and P(X) is a space spanned by products of linear forms. The first main result of this paper is the construction of a canonical basis for D(X). We show that it is dual to the canonical basis for P(X) that is already known. The second main result of this paper is the construction of a new family of zonotopal spaces that is far more general than the ones that were recently studied by Ardila-Postnikov, Holtz-Ron, Holtz-Ron-Xu, Li-Ron, and others. We call the underlying combinatorial structure of those spaces forward exchange matroid. A forward exchange matroid is an ordered matroid together with a subset of its set of bases that satisfies a weak version of the basis exchange axiom. △ Less

Submitted 31 March, 2016; v1 submitted 17 April, 2012; originally announced April 2012.

Comments: 34 pages, 4 figures, minor corrections, same as journal version (up to layout)

MSC Class: 05A15; 05B35; 13B25; 16S32; 41A15 (Primary); 05C31; 41A63; 47F05 (Secondary)

Journal ref: Advances in Mathematics 294 (May 2016), 819-852

Contractile units in disordered actomyosin bundles arise from F-actin buckling

Authors: Martin Lenz, Todd Thoresen, Margaret L. Gardel, Aaron R. Dinner

Abstract: Bundles of filaments and motors are central to contractility in cells. The classic example is striated muscle, where actomyosin contractility is mediated by highly organized sarcomeres which act as fundamental contractile units. However, many contractile bundles in vivo and in vitro lack sarcomeric organization. Here we propose a model for how contractility can arise in actomyosin bundles without… ▽ More Bundles of filaments and motors are central to contractility in cells. The classic example is striated muscle, where actomyosin contractility is mediated by highly organized sarcomeres which act as fundamental contractile units. However, many contractile bundles in vivo and in vitro lack sarcomeric organization. Here we propose a model for how contractility can arise in actomyosin bundles without sarcomeric organization and validate its predictions with experiments on a reconstituted system. In the model, internal stresses in frustrated arrangements of motors with diverse velocities cause filaments to buckle, leading to overall shortening. We describe the onset of buckling in the presence of stochastic actin-myosin detachment and predict that buckling-induced contraction occurs in an intermediate range of motor densities. We then calculate the size of the "contractile units" associated with this process. Consistent with these results, our reconstituted actomyosin bundles contract at relatively high motor density, and we observe buckling at the predicted length scale. △ Less

Submitted 19 January, 2012; originally announced January 2012.

Comments: 5 pages, 4 figures, Supporting text and movies attached

Journal ref: Phys. Rev. Lett. 108, 238107 (2012)

Coils and transformers - often used but seldomly explained correctly

Authors: Michael Lenz, Torsten Schmidt

Abstract: The devices coil and transformer are subjects of interest in numerous schoolbooks, in introductory scientific textbooks of physics and engineering, and in laboratory courses at universities. Many descriptions, however, draw a somewhat distorted picture of the underlying physical mechanisms and provide half-knowledge or even clear misconceptions that should not be left uncommented and are therefore… ▽ More The devices coil and transformer are subjects of interest in numerous schoolbooks, in introductory scientific textbooks of physics and engineering, and in laboratory courses at universities. Many descriptions, however, draw a somewhat distorted picture of the underlying physical mechanisms and provide half-knowledge or even clear misconceptions that should not be left uncommented and are therefore studied in detail: (1) Primary and secondary voltage at a transformer have a different sign. (2) Electromagnetic induction is the only mechanism of importance for coils and transformers. (3) The terminal voltage at coils and transformers is compensated by the so-called "induced voltage" (emf), which explains why Kirchhoff's voltage law also applies to coils and transformers. (4) The cores of coils and transformers are used for their ability to store energy. Energy is transported from the primary to the secondary coil within the magnetic core. (5) The stray magnetic and electric fields are sencondary effects not having a major effect on energy transport. (6) The higher the load current, the easier a transformer goes into saturation. (7) The higher the number of turns at the primary coil, the larger the magnetic flux in the core. (8) Transformers with cores having an air gap have a lower coupling factor, because the stray inductivity increases. In the paper, the most important characteristics of coil and transformers are derived directly from Maxwell's equation for idealised conditions, and subsequently, the different misconceptions are discussed and corrected. △ Less

Submitted 6 May, 2012; v1 submitted 13 July, 2011; originally announced July 2011.

Comments: Preprint in German language

Matroids and log-concavity

Authors: Matthias Lenz

Abstract: We show that f-vectors of matroid complexes of realisable matroids are log-concave. This was conjectured by Mason in 1972. Our proof uses the recent result by Huh and Katz who showed that the coefficients of the characteristic polynomial of a realisable matroid form a log-concave sequence. We also discuss the relationship between log-concavity of f-vectors and h-vectors of matroids. In the last se… ▽ More We show that f-vectors of matroid complexes of realisable matroids are log-concave. This was conjectured by Mason in 1972. Our proof uses the recent result by Huh and Katz who showed that the coefficients of the characteristic polynomial of a realisable matroid form a log-concave sequence. We also discuss the relationship between log-concavity of f-vectors and h-vectors of matroids. In the last section we explain the connection between zonotopal algebra and f-vectors and characteristic polynomials of matroids. △ Less

Submitted 10 June, 2013; v1 submitted 15 June, 2011; originally announced June 2011.

Comments: 9 pages

MSC Class: Primary: 05A20; 05B35; Secondary: 05E45; 05C31

Requirements for contractility in disordered cytoskeletal bundles

Authors: Martin Lenz, Margaret L. Gardel, Aaron R. Dinner

Abstract: Actomyosin contractility is essential for biological force generation, and is well understood in highly organized structures such as striated muscle. Additionally, actomyosin bundles devoid of this organization are known to contract both in vivo and in vitro, which cannot be described by standard muscle models. To narrow down the search for possible contraction mechanisms in these systems, we inve… ▽ More Actomyosin contractility is essential for biological force generation, and is well understood in highly organized structures such as striated muscle. Additionally, actomyosin bundles devoid of this organization are known to contract both in vivo and in vitro, which cannot be described by standard muscle models. To narrow down the search for possible contraction mechanisms in these systems, we investigate their microscopic symmetries. We show that contractile behavior requires non-identical motors that generate large enough forces to probe the nonlinear elastic behavior of F-actin. This suggests a role for filament buckling in the contraction of these bundles, consistent with recent experimental results on reconstituted actomyosin bundles. △ Less

Submitted 27 January, 2012; v1 submitted 5 January, 2011; originally announced January 2011.

Comments: 10 pages, 6 figures; text shortened

Journal ref: New J. Phys. 14, 033037 (2012)

Hierarchical zonotopal power ideals

Authors: Matthias Lenz

Abstract: Zonotopal algebra deals with ideals and vector spaces of polynomials that are related to several combinatorial and geometric structures defined by a finite sequence of vectors. Given such a sequence X, an integer k>=-1 and an upper set in the lattice of flats of the matroid defined by X, we define and study the associated hierarchical zonotopal power ideal. This ideal is generated by powers of lin… ▽ More Zonotopal algebra deals with ideals and vector spaces of polynomials that are related to several combinatorial and geometric structures defined by a finite sequence of vectors. Given such a sequence X, an integer k>=-1 and an upper set in the lattice of flats of the matroid defined by X, we define and study the associated hierarchical zonotopal power ideal. This ideal is generated by powers of linear forms. Its Hilbert series depends only on the matroid structure of X. Via the Tutte polynomial, it is related to various other matroid invariants, e.g. the shelling polynomial and the characteristic polynomial. This work unifies and generalizes results by Ardila-Postnikov on power ideals and by Holtz-Ron and Holtz-Ron-Xu on (hierarchical) zonotopal algebra. We also generalize a result on zonotopal Cox modules that were introduced by Sturmfels-Xu. △ Less

Submitted 5 January, 2012; v1 submitted 4 November, 2010; originally announced November 2010.

Comments: 22 pages, 1 figure; small changes reflecting referees' comments; to appear in European Journal of Combinatorics

MSC Class: 05A15; 05B35; 05C31; 13B25; 16S32

Journal ref: European Journal of Combinatorics 33 (2012), no. 6, 1120-1141