pith. sign in

arxiv: 2605.11012 · v2 · pith:T27QYCIHnew · submitted 2026-05-10 · ❄️ cond-mat.soft

Inverse Design of Metainterfaces for Static Friction Control: Beyond the Hertzian Limit

Pith reviewed 2026-05-20 22:58 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords inverse designstatic frictionmetainterfacesaxisymmetric asperitiesnonlinear contactdifferentiable physicstribological surfacesboundary element method
0
0 comments X

The pith

General axisymmetric asperities enable inverse design of metainterfaces that achieve arbitrary nonlinear static friction laws.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes a computational method to program static friction by designing the shapes of a few asperities in each periodic cell of an interface. Standard Hertzian contacts produce only linear area-load scaling, which limits functional range, but general axisymmetric profiles can produce nonlinear responses. The approach places a fully differentiable contact solver inside an optimization loop that uses physical gradients to find topographies matching any chosen target friction curve. Designs are then checked against high-fidelity boundary-element simulations to confirm they remain valid.

Core claim

A fully differentiable contact mechanics engine is embedded in a neural-network-plus-optimizer pipeline that automatically discovers non-Hertzian axisymmetric asperity topographies whose collective real contact area reproduces prescribed nonlinear friction-force versus load relations, using only a small number of asperities per unit cell and with results validated by independent boundary-element simulations.

What carries the argument

Regularized physical gradients through a differentiable contact solver that computes real contact area for arbitrary axisymmetric shapes and feeds the result directly into a quadratic optimizer.

If this is right

  • Surfaces can be engineered to exhibit any desired friction-load curve rather than being restricted to linear scaling.
  • Complex target behaviors become achievable with only a few asperities per periodic cell.
  • The same differentiable-physics loop can be reused for other contact quantities such as stiffness or thermal resistance.
  • Scale-invariant designs become possible because the optimization operates on normalized geometry and load.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The framework could be extended to include adhesion or viscoelasticity once the underlying solver incorporates those effects.
  • Fabrication tolerances on asperity height and curvature would need to be quantified to determine how closely real devices match the predicted curves.
  • Similar inverse-design loops might apply to controlling other interface properties such as electrical contact resistance or fluid leakage.

Load-bearing premise

The differentiable contact engine must correctly compute real contact area and its dependence on load for any axisymmetric shape so that the optimized profiles remain physically realizable.

What would settle it

Fabricate the optimized asperity arrays and measure the static friction force as a function of normal load to test whether the experimental curve matches the target nonlinear law within measurement uncertainty.

Figures

Figures reproduced from arXiv: 2605.11012 by Arnav Singhal, Ga\"etan Cortes, Jacopo Bilotto, Jean-Fran\c{c}ois Molinari, Joaquin Garcia-Suarez, Lucas Fourel.

Figure 1
Figure 1. Figure 1: (a) Illustration of axisymmetric asperities in power-law form [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The physics-informed inverse modeling pipeline based on a deep neural network. The input [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Topographical test-set sub-domains and corresponding neural reconstructions of contact [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Training dynamics and hyperparameter schedules for the neural surrogate. (a) Convergence [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The left panels, (a) and (c), illustrate the macroscopic constitutive response, contact area [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Inverse design performance on out-of-distribution targets: a saturating limit, a bilinear [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (Left) The macroscopic contact response, [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Additional examples of BEM validation plots for representative surfaces from the testing [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
read the original abstract

Programming the static friction of mechanical interfaces is critical for soft robotics, haptics, and precision gripping. Static friction is governed by the real contact area, and standard rough surfaces exhibit a linear area-load scaling inherent to classical Archard and Greenwood-Williamson models, severely restricting their functional range. Here, we propose a framework for the inverse design of tribological metainterfaces engineered for programmable contact behaviors. By utilizing general axisymmetric asperities, we unlock nonlinear macroscopic responses unattainable by standard Hertzian contacts. To solve the inverse problem, we embed a fully differentiable contact mechanics engine within a neural network and a quadratic optimizer. We leverage regularized physical gradients to automatically discover non-standard topographies that reproduce complex target friction laws, with only a few asperities in unit cells. The predicted designs are strictly validated against high-fidelity Boundary Element Method (BEM) simulations. This framework bridges data-driven optimization and rigorous physics, offering a scale-invariant pathway for discovering functional tribological surfaces.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper presents a framework for inverse design of tribological metainterfaces to program static friction beyond the linear area-load scaling of Hertzian contacts. It uses general axisymmetric asperities in unit cells, embeds a fully differentiable contact mechanics engine inside a neural network and quadratic optimizer, and employs regularized physical gradients to discover non-standard topographies that match complex target friction laws. Final designs are validated against high-fidelity BEM simulations.

Significance. If the central results hold, the work offers a scale-invariant, physics-constrained route to custom nonlinear friction responses with minimal asperities per cell. This could impact soft robotics, haptics, and precision gripping by moving beyond Archard/Greenwood-Williamson restrictions. The combination of differentiable physics with optimization is a methodological strength that enables automated discovery while retaining physical interpretability.

major comments (1)
  1. [Method (differentiable contact mechanics engine) and Results (validation)] The central claim requires that the embedded differentiable contact mechanics engine produces real-contact-area predictions whose load dependence is sufficiently faithful to BEM for arbitrary (non-spherical) axisymmetric profiles. The abstract states that final designs are “strictly validated” against BEM, yet no quantitative benchmark (error metrics, load-area curves, or direct comparison) is provided for the non-standard shapes that appear in the optimized solutions. If the surrogate deviates systematically from BEM on those profiles, the optimizer can converge to topographies that satisfy the surrogate but fail the true physics. Please add a dedicated validation subsection with direct engine-vs-BEM comparisons on representative optimized asperity shapes.
minor comments (1)
  1. [Abstract] The abstract refers to “regularized physical gradients” without specifying the regularization form or strength; a brief parenthetical or reference would improve clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript. The feedback on validation of the differentiable contact mechanics engine is well taken, and we have prepared a point-by-point response below. We agree that strengthening the quantitative comparison for optimized non-standard profiles will improve the paper.

read point-by-point responses
  1. Referee: [Method (differentiable contact mechanics engine) and Results (validation)] The central claim requires that the embedded differentiable contact mechanics engine produces real-contact-area predictions whose load dependence is sufficiently faithful to BEM for arbitrary (non-spherical) axisymmetric profiles. The abstract states that final designs are “strictly validated” against BEM, yet no quantitative benchmark (error metrics, load-area curves, or direct comparison) is provided for the non-standard shapes that appear in the optimized solutions. If the surrogate deviates systematically from BEM on those profiles, the optimizer can converge to topographies that satisfy the surrogate but fail the true physics. Please add a dedicated validation subsection with direct engine-vs-BEM comparisons on representative optimized asperity shapes.

    Authors: We agree that explicit quantitative validation for the non-standard axisymmetric profiles discovered by the optimizer is necessary to fully support the central claims. Although the manuscript states that final designs are validated against high-fidelity BEM, we acknowledge that a dedicated subsection presenting error metrics, load-area curves, and direct comparisons specifically for the optimized (non-Hertzian) shapes was not included in the main text. We will add a new subsection in the Results section that reports these benchmarks, including relative errors in predicted contact area, overlaid load-area curves, and representative profile comparisons between the differentiable engine and BEM. This revision will demonstrate that the surrogate remains faithful for the relevant topographies and that the optimizer does not exploit systematic deviations. revision: yes

Circularity Check

0 steps flagged

No significant circularity; optimization and validation remain independent

full rationale

The paper embeds a differentiable contact mechanics engine inside an optimizer to discover axisymmetric asperity topographies that match prescribed target friction laws, then validates the resulting designs against separate high-fidelity BEM simulations. No equation or step in the abstract or described workflow shows the target law or optimized topography reducing to a fitted parameter by construction, nor does any load-bearing premise rest on a self-citation chain. The use of regularized physical gradients and external BEM checks keeps the central claim self-contained against independent benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The framework rests on standard tribological assumptions about real contact area controlling friction and on the accuracy of the differentiable engine; no new free parameters or invented entities are explicitly introduced in the abstract.

axioms (1)
  • domain assumption Static friction is governed by the real contact area.
    Stated in the opening sentence of the abstract as the governing principle.

pith-pipeline@v0.9.0 · 5725 in / 1171 out tokens · 39501 ms · 2026-05-20T22:58:46.221973+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages · 1 internal anchor

  1. [1]

    Introduction, in: Allgower, E.L., Georg, K

    Allgower, E.L., Georg, K., 1990. Introduction, in: Allgower, E.L., Georg, K. (Eds.), Numerical Continuation Methods: An Introduction. Springer, Berlin, Heidelberg, pp. 1–6. URL:https: //doi.org/10.1007/978-3-642-61257-2_1, doi:10.1007/978-3-642-61257-2_1

  2. [2]

    Elastic deformation and the laws of friction

    Archard, J.F., Allibone, T.E., 1997. Elastic deformation and the laws of friction. Pro- ceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 243, 190–205. URL:https://royalsocietypublishing.org/doi/10.1098/rspa.1957.0214, doi:10. 1098/rspa.1957.0214

  3. [3]

    End-to-End Differentiable Physics for Learning and Control, in: Advances in Neural Information Processing Systems, Curran Associates, Inc

    de Avila Belbute-Peres, F., Smith, K., Allen, K., Tenenbaum, J., Kolter, J.Z., 2018. End-to-End Differentiable Physics for Learning and Control, in: Advances in Neural Information Processing Systems, Curran Associates, Inc. URL:https://papers.nips.cc/paper_files/paper/2018/ hash/842424a1d0595b76ec4fa03c46e8d755-Abstract.html

  4. [4]

    Designing metainterfaces with specified friction laws

    Aymard, A., Delplanque, E., Dalmas, D., Scheibert, J., 2024. Designing metainterfaces with specified friction laws. Science 383, 200–204. URL:https://www.science.org/doi/10.1126/ science.adk4234, doi:10.1126/science.adk4234

  5. [5]

    Curriculum learning, in: Proceedings of the 26th Annual International Conference on Machine Learning, Association for Computing Machinery, New York, NY, USA

    Bengio, Y., Louradour, J., Collobert, R., Weston, J., 2009. Curriculum learning, in: Proceedings of the 26th Annual International Conference on Machine Learning, Association for Computing Machinery, New York, NY, USA. pp. 41–48. URL:https://dl.acm.org/doi/10.1145/1553374. 1553380, doi:10.1145/1553374.1553380

  6. [6]

    Bessa, M.A., Bostanabad, R., Liu, Z., Hu, A., Apley, D.W., Brinson, C., Chen, W., Liu, W.K.,

  7. [8]

    Fluid-mediated impact of soft solids

    Bilotto, J., Kolinski, J.M., Lecampion, B., Molinari, J.F., Subhash, G., Garcia-Suarez, J., 2024. Fluid-mediated impact of soft solids. Journal of Fluid Mechanics 997, A35. URL:https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/ fluidmediated-impact-of-soft-solids/E208CFA7E3181F48FB66DCE0DF20A182, doi:10.1017/ jfm.2024.820

  8. [9]

    Masset, R

    Bonari, J., Paggi, M., Dini, D., 2022. A new finite element paradigm to solve contact prob- lems with roughness. International Journal of Solids and Structures 253, 111643. URL: https://www.sciencedirect.com/science/article/pii/S0020768322001640, doi:10.1016/j. ijsolstr.2022.111643

  9. [10]

    Application des potentiels à l’étude de l’équilibre et du mouve- ment des solides élastiques

    Boussinesq, J., 1885. Application des potentiels à l’étude de l’équilibre et du mouve- ment des solides élastiques. Gauthier-Villars. URL:https://gallica.bnf.fr/ark:/12148/ bpt6k9651115r

  10. [11]

    Oxford Univer- sity Press (2018).https://doi.org/10.1093/oso/9780198814788.001.0001

    Bowden, F.P., Tabor, D., 2001. Area of Contact Between Solids, in: Bowden, F.P., Tabor, D. (Eds.), The Friction and Lubrication of Solids. Oxford University Press, p. 0. URL:https:// doi.org/10.1093/oso/9780198507772.003.0002, doi:10.1093/oso/9780198507772.003.0002. 13

  11. [12]

    The role of adhesion in contact mechanics

    Ciavarella, M., Joe, J., Papangelo, A., Barber, J.R., 2019. The role of adhesion in contact mechanics. Journal of The Royal Society Interface 16, 20180738. URL:https://doi.org/10. 1098/rsif.2018.0738, doi:10.1098/rsif.2018.0738

  12. [13]

    Inverse design and flexible parameterization of meta-optics using algorithmic differentiation

    Colburn, S., Majumdar, A., 2021. Inverse design and flexible parameterization of meta-optics using algorithmic differentiation. Communications Physics 4, 65. URL:https://www.nature. com/articles/s42005-021-00568-6, doi:10.1038/s42005-021-00568-6

  13. [14]

    Squeaking at soft–rigid frictional interfaces

    Djellouli, A., Albertini, G., Wilt, J., Tournat, V., Weitz, D., Rubinstein, S., Bertoldi, K., 2026. Squeaking at soft–rigid frictional interfaces. Nature 650, 891–897. URL:https://www.nature. com/articles/s41586-026-10132-3, doi:10.1038/s41586-026-10132-3

  14. [15]

    JOSS5(51), 2121 (2020) https://doi

    Frérot, L., Anciaux, G., Rey, V., Pham-Ba, S., Molinari, J.F., 2020. Tamaas: a library for elastic-plastic contact of periodic rough surfaces. Journal of Open Source Software 5, 2121. URL: https://joss.theoj.org/papers/10.21105/joss.02121, doi:10.21105/joss.02121

  15. [16]

    Elastische Beanspruchung des Erdbodens unter Fundamenten

    Föppl, L., 1941. Elastische Beanspruchung des Erdbodens unter Fundamenten. Forschung auf dem Gebiet des Ingenieurwesens A 12, 31–39. URL:https://doi.org/10.1007/BF02593958, doi:10.1007/BF02593958

  16. [17]

    A Matter of Shape: Contact Area Optimization in Soft Lubricated Impact

    Garcia-Suarez, J., 2026. A Matter of Shape: Contact Area Optimization in Soft Lubricated Impact. Tribology Letters 74, 15. URL:https://doi.org/10.1007/s11249-026-02108-1, doi:10.1007/s11249-026-02108-1

  17. [18]

    A., & Williamson, J

    Greenwood, J.A., Williamson, J.B.P., Bowden, F.P., 1997. Contact of nominally flat surfaces. Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 295, 300–319. URL:https://royalsocietypublishing.org/doi/10.1098/rspa.1966.0242, doi:10. 1098/rspa.1966.0242

  18. [19]

    Elastic-Plastic Contact Analysis of a Sphere and a Rigid Flat

    Kogut, L., Etsion, I., 2002. Elastic-Plastic Contact Analysis of a Sphere and a Rigid Flat. Journal of Applied Mechanics 69, 657–662. URL:https://doi.org/10.1115/1.1490373, doi:10.1115/ 1.1490373

  19. [20]

    Inverse-designed spinodoid meta- materials

    Kumar, S., Tan, S., Zheng, L., Kochmann, D.M., 2020. Inverse-designed spinodoid meta- materials. npj Computational Materials 6, 1–10. URL:https://www.nature.com/articles/ s41524-020-0341-6, doi:10.1038/s41524-020-0341-6

  20. [21]

    Training Deep Neural Networks for the Inverse Design of Nanophotonic Structures

    Liu, D., Tan, Y., Khoram, E., Yu, Z., 2018. Training Deep Neural Networks for the Inverse Design of Nanophotonic Structures. ACS Photonics 5, 1365–1369. URL:https://doi.org/10.1021/ acsphotonics.7b01377, doi:10.1021/acsphotonics.7b01377

  21. [22]

    Decoupled Weight Decay Regularization

    Loshchilov, I., Hutter, F., 2019. Decoupled Weight Decay Regularization. URL:http://arxiv. org/abs/1711.05101, doi:10.48550/arXiv.1711.05101. arXiv:1711.05101 [cs]

  22. [23]

    A Review of Advances in Fabrication Methods and Assistive Technologies of Micro-Structured Surfaces

    Ma, Y., Zhang, G., Cao, S., Huo, Z., Han, J., Ma, S., Huang, Z., 2023. A Review of Advances in Fabrication Methods and Assistive Technologies of Micro-Structured Surfaces. Processes 11,

  23. [24]

    URL:https://www.mdpi.com/2227-9717/11/5/1337, doi:10.3390/pr11051337

  24. [25]

    A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code

    McKay, M.D., Beckman, R.J., Conover, W.J., 1979. A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code. Technometrics 21, 239–245. URL:https://www.jstor.org/stable/1268522, doi:10.2307/1268522

  25. [26]

    Friction on demand: A generative framework for the inverse design of metainterfaces

    Mouton, V., Mélot, A., 2026. Friction on demand: A generative framework for the inverse design of metainterfaces. Tribology International 218, 111724. URL:https://www.sciencedirect.com/ science/article/pii/S0301679X26000666, doi:10.1016/j.triboint.2026.111724

  26. [27]

    Tailoring Frictional Properties of Surfaces Using Diffusion Models

    Nordhagen, E., Sveinsson, H.A., Malthe-Sørenssen, A., 2025. Tailoring Frictional Properties of Surfaces Using Diffusion Models. The Journal of Physical Chemistry C 129, 14559–14564. URL: https://pubs.acs.org/doi/10.1021/acs.jpcc.5c02768, doi:10.1021/acs.jpcc.5c02768. 14

  27. [28]

    Handbook of Contact Mechanics: Exact Solutions of Axisymmetric Contact Problems

    Popov, V.L., Heß, M., Willert, E., 2019. Handbook of Contact Mechanics: Exact Solutions of Axisymmetric Contact Problems. Springer Berlin Heidelberg, Berlin, Heidelberg. URL:http: //link.springer.com/10.1007/978-3-662-58709-6, doi:10.1007/978-3-662-58709-6

  28. [29]

    Frictional Contact of Soft Polymeric Shells

    Sahli, R., Mikkelsen, J., Boye, M.S., Dias, M.A., Aghababaei, R., 2024. Frictional Contact of Soft Polymeric Shells. Physical Review Letters 133, 106202. URL:https://link.aps.org/doi/10. 1103/PhysRevLett.133.106202, doi:10.1103/PhysRevLett.133.106202

  29. [30]

    Why soft contacts are stickier when breaking than when making them

    Sanner, A., Kumar, N., Dhinojwala, A., Jacobs, T.D.B., Pastewka, L., 2024. Why soft contacts are stickier when breaking than when making them. Science Advances 10, eadl1277. URL:https: //www.science.org/doi/full/10.1126/sciadv.adl1277, doi:10.1126/sciadv.adl1277

  30. [31]

    Zur Frage der Druckverteilung unter elastisch gelagerten Tragwerken

    Schubert, G., 1942. Zur Frage der Druckverteilung unter elastisch gelagerten Tragwerken. Ingenieur-Archiv 13, 132–147. URL:https://doi.org/10.1007/BF02095912, doi:10.1007/ BF02095912

  31. [32]

    Real Area of Contact in a Soft Transparent Interface by Particle Exclusion Microscopy

    Schulze, K.D., Bennett, A.I., Marshall, S., Rowe, K.G., Dunn, A.C., 2016. Real Area of Contact in a Soft Transparent Interface by Particle Exclusion Microscopy. Journal of Tribology 138. URL: https://doi.org/10.1115/1.4032822, doi:10.1115/1.4032822

  32. [33]

    The relation between load and penetration in the axisymmetric boussinesq problem for a punch of arbitrary profile

    Sneddon, I.N., 1965. The relation between load and penetration in the axisymmetric boussinesq problem for a punch of arbitrary profile. International Journal of Engineering Science 3, 47–

  33. [34]

    1016/0020-7225(65)90019-4

    URL:https://www.sciencedirect.com/science/article/pii/0020722565900194, doi:10. 1016/0020-7225(65)90019-4

  34. [35]

    The nonlinear nature of friction

    Urbakh, M., Klafter, J., Gourdon, D., Israelachvili, J., 2004. The nonlinear nature of friction. Nature 430, 525–528. URL:https://www.nature.com/articles/nature02750, doi:10.1038/ nature02750

  35. [36]

    Persson’s theory of purely normal elastic rough surface contact: A tutorial based on stochastic process theory

    Xu, Y., Li, X., Chen, Q., Zhou, Y., 2024. Persson’s theory of purely normal elastic rough surface contact: A tutorial based on stochastic process theory. International Journal of Solids and Structures 290, 112684. URL:https://www.sciencedirect.com/science/article/pii/ S0020768324000416, doi:10.1016/j.ijsolstr.2024.112684

  36. [37]

    JAX-FEM: A differen- tiable GPU-accelerated 3D finite element solver for automatic inverse design and mechanistic data science

    Xue, T., Liao, S., Gan, Z., Park, C., Xie, X., Liu, W.K., Cao, J., 2023. JAX-FEM: A differen- tiable GPU-accelerated 3D finite element solver for automatic inverse design and mechanistic data science. Computer Physics Communications 291, 108802. URL:https://www.sciencedirect. com/science/article/pii/S0010465523001479, doi:10.1016/j.cpc.2023.108802

  37. [38]

    The Contact of Elastic Regular Wavy Surfaces Revisited

    Yastrebov, V.A., Anciaux, G., Molinari, J.F., 2014. The Contact of Elastic Regular Wavy Surfaces Revisited. TribologyLetters56, 171–183. URL:https://doi.org/10.1007/s11249-014-0395-z, doi:10.1007/s11249-014-0395-z

  38. [39]

    Friction Anisotropy with Respect to Topographic Orientation

    Yu, C., Wang, Q.J., 2012. Friction Anisotropy with Respect to Topographic Orientation. Scientific Reports 2, 988. URL:https://www.nature.com/articles/srep00988, doi:10.1038/srep00988

  39. [40]

    Normal contact of metainterfaces: The roles of finite size and microcontact interactions

    Zeka, D., Blal, N., Fekak, F.E., Duval, A., Gravouil, A., Scheibert, J., 2026. Normal contact of metainterfaces: The roles of finite size and microcontact interactions. Journal of the Mechanics and Physics of Solids 214, 106646. URL:https://www.sciencedirect.com/science/article/ pii/S0022509626001468, doi:10.1016/j.jmps.2026.106646. A Details on Neural Ne...