Anisotropic buckling of few-layer black phosphorus

Andres Castellanos-Gomez; Luis Vaquero-Garzon; Riccardo Frisenda

arxiv: 1907.07371 · v1 · pith:AIP4UMVHnew · submitted 2019-07-17 · ⚛️ physics.app-ph · cond-mat.mes-hall· cond-mat.mtrl-sci

Anisotropic buckling of few-layer black phosphorus

Luis Vaquero-Garzon , Riccardo Frisenda , Andres Castellanos-Gomez This is my paper

Pith reviewed 2026-05-24 20:06 UTC · model grok-4.3

classification ⚛️ physics.app-ph cond-mat.mes-hallcond-mat.mtrl-sci

keywords black phosphorusbucklinganisotropyYoung's modulusripplesfew-layercompression2D materials

0 comments

The pith

Compression of few-layer black phosphorus on a substrate produces ripples whose period is 40 percent longer along the zig-zag crystal direction than along the armchair direction.

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

The paper establishes that the puckered atomic structure of black phosphorus causes its buckling under compression to be direction-dependent. When flakes adhered to a compliant substrate are compressed, the resulting periodic ripples have different wavelengths depending on whether the force is applied along the armchair or zig-zag axis. Analysis of those wavelengths with the thin-film buckling formula gives Young's modulus values of 35.1 GPa armchair and 93.3 GPa zig-zag. A reader would care because the result supplies a direct optical route to the material's anisotropic stiffness without separate mechanical tests. The finding also shows how the atomic puckering produces measurable macroscopic differences in elastic response.

Core claim

When black phosphorus flakes adhered onto a compliant substrate are subjected to compression they buckle into periodic ripples. The ripple period is 40 percent longer when compression is applied along the zig-zag crystal direction than along the armchair direction. This anisotropy originates in the puckered honeycomb crystal structure. Quantitative analysis of the observed ripple periods with the standard continuum buckling model determines the Young's modulus of few-layer black phosphorus as 35.1 ± 6.3 GPa along the armchair direction and 93.3 ± 21.8 GPa along the zig-zag direction.

What carries the argument

The periodic rippling instability under directional compression, whose wavelength depends on the directional Young's modulus through the thin-film buckling relation on a compliant substrate.

If this is right

The buckling wavelength supplies a non-contact method to extract the two principal Young's moduli of few-layer black phosphorus.
The zig-zag direction is stiffer than the armchair direction by a factor of approximately 2.7.
The atomic puckering of the lattice is the direct cause of the observed difference in ripple periods.
Any mechanical model of black phosphorus devices under strain must incorporate this directional stiffness difference.

Where Pith is reading between the lines

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

The same ripple-period measurement could be applied to other layered puckered materials to obtain their anisotropic moduli without separate test equipment.
Device designers working with black phosphorus under compression can predict the direction of easiest buckling from the crystal orientation alone.
Strain sensors or actuators built from black phosphorus could exploit the 40 percent difference in ripple wavelength for directional response.

Load-bearing premise

The standard continuum thin-film buckling model applies directly to few-layer black phosphorus with ripple wavelength fixed only by directional modulus, thickness, and substrate properties.

What would settle it

Independent measurement of Young's modulus along each crystal direction in the same few-layer flakes by nanoindentation or atomic-force bending, followed by direct numerical comparison to the reported values 35.1 GPa and 93.3 GPa.

Figures

Figures reproduced from arXiv: 1907.07371 by Andres Castellanos-Gomez, Luis Vaquero-Garzon, Riccardo Frisenda.

**Figure 1.** Figure 1: a) Schematic diagram of the process employed to fabricate the samples, where the flakes are transferred onto a stretched elastomeric substrate. When the stress is released, the flakes are subjected to compressive stress that produces ripples with a certain thickness-dependent period. b) Transmission mode optical microscopy images of a black phosphorus multilayer flake after releasing the stress on the elas… view at source ↗

**Figure 2.** Figure 2: a) Optical microscopy images of ripples with different periods for 4 black phosphorus flakes of different thickness. b) Wavelength vs thickness graph for several black phosphorus flakes with different thickness, their error bars with 95% confidence curves. The solid blue line is the linear fit and the shaded area around indicates the uncertainty of the fit (95% confidence) [PITH_FULL_IMAGE:figures/full_fi… view at source ↗

**Figure 3.** Figure 3: a) Process of ripples formation in the two directions (ZZ and AC) in a controlled way two manipulators to apply the compression along both directions. b) Optical microscopy images of black phosphorus compressed in both directions. The period of the ripple pattern sizeably depends on the direction where the compression is applied. The coloured arrows indicate the ZZ (blue) and AC (red) compression direction… view at source ↗

read the original abstract

When a two-dimensional material, adhered onto a compliant substrate, is subjected to compression it can undertake a buckling instability yielding to a periodic rippling. Interestingly, when black phosphorus flakes are compressed along the zig-zag crystal direction the flake buckles forming ripples with a 40% longer period than that obtained when the compression is applied along the armchair direction. This anisotropic buckling stems from the puckered honeycomb crystal structure of black phosphorus and a quantitative analysis of the ripple period allows us to determine the Youngs's modulus of few-layer black phosphorus along the armchair direction (EbP_AC = 35.1 +- 6.3 GPa) and the zig-zag direction (EbP_ZZ = 93.3 +- 21.8 GPa).

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.

Desk Editor's Note private letter to a colleague

The paper measures 40% anisotropic ripple periods in few-layer black phosphorus and extracts directional Young's moduli of 35 GPa AC and 93 GPa ZZ, but the continuum buckling model needs explicit checks for few-layer validity.

read the letter

The main thing to know is that they compressed few-layer BP flakes on a substrate and saw ripples with a 40% longer period along zig-zag than armchair, then used the wavelength difference to report Young's moduli of 35.1 ± 6.3 GPa armchair and 93.3 ± 21.8 GPa zig-zag. Those specific numbers with uncertainties are the new experimental result here. The anisotropy itself follows from the puckered lattice, as expected from earlier theory, but the quantitative values for few-layer material are what stand out. The experiment is straightforward and gives concrete numbers instead of just direction-dependent behavior, which is useful for device modeling. The analysis applies the standard thin-film buckling relation between ripple period, film modulus, thickness, and substrate properties. That step is where the soft spot sits. For only a few layers the homogeneous-plate assumption could miss discrete atomic effects or direction-dependent adhesion, and nothing in the abstract shows they tested those or ran sensitivity checks. The reported uncertainties are large, around 18-23%, which at least signals the precision limits. No circularity issue since the moduli are outputs from the fit rather than inputs. This is the sort of paper that device engineers working with black phosphorus or 2D material mechanics would want to read for the numbers. It deserves peer review so referees can examine the raw ripple data and the model application details.

Referee Report

2 major / 2 minor

Summary. The paper claims that few-layer black phosphorus flakes on a compliant substrate exhibit anisotropic buckling under compression, with ripple periods 40% longer along the zig-zag direction than the armchair direction due to the puckered honeycomb crystal structure. A quantitative analysis of the observed ripple periods is used to extract the directional Young's moduli as EbP_AC = 35.1 ± 6.3 GPa and EbP_ZZ = 93.3 ± 21.8 GPa.

Significance. If the underlying model holds, the work demonstrates an experimental route to extract anisotropic in-plane moduli of 2D materials from buckling wavelengths, which could complement other techniques such as AFM indentation. The reported values are consistent with the expected strong anisotropy arising from the puckered structure and provide concrete numbers that can be compared against first-principles calculations or other measurements.

major comments (2)

[Abstract] Abstract: the moduli are reported with uncertainties but the abstract (and by extension the quantitative analysis) provides no raw ripple-period data, no explicit form of the buckling relation used, and no validation that the continuum model reproduces the measured wavelengths; this is load-bearing for the central claim that the periods directly yield the moduli.
[Quantitative analysis (results section)] The extraction implicitly applies the standard thin-film buckling wavelength formula (or its anisotropic generalization) with effective thickness t and substrate modulus E_s as the only other parameters; the manuscript does not discuss or test corrections for discrete-layer effects, direction-dependent adhesion, or the validity of replacing the atomic puckered lattice by a homogeneous anisotropic plate, which directly affects whether the reported numbers can be taken as the intrinsic directional moduli.

minor comments (2)

[Abstract] Abstract: 'Youngs's modulus' contains a typographical error and should read 'Young's modulus'.
[Abstract] The abstract states that the anisotropic buckling 'stems from the puckered honeycomb crystal structure' but does not cite or show the corresponding atomic model or simulation that would make this link quantitative.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. We address each major comment below and will revise the manuscript to improve the clarity and completeness of the quantitative analysis.

read point-by-point responses

Referee: [Abstract] Abstract: the moduli are reported with uncertainties but the abstract (and by extension the quantitative analysis) provides no raw ripple-period data, no explicit form of the buckling relation used, and no validation that the continuum model reproduces the measured wavelengths; this is load-bearing for the central claim that the periods directly yield the moduli.

Authors: The abstract serves as a concise summary and space constraints preclude inclusion of raw data or equations. Raw ripple periods (with means and standard deviations for both crystal directions) are reported in the results section, Figure 2, and Table S1. The explicit buckling relation is the standard thin-film wrinkling formula generalized to anisotropic in-plane moduli, λ = 2π t (E / 3 E_s)^{1/3}, and is stated in the Methods. Validation consists of substituting the extracted directional moduli back into the formula and recovering the observed wavelengths within experimental error. We will add a brief statement of the formula and a note on this consistency check to the abstract and ensure the relation is highlighted in the results section. revision: yes
Referee: [Quantitative analysis (results section)] The extraction implicitly applies the standard thin-film buckling wavelength formula (or its anisotropic generalization) with effective thickness t and substrate modulus E_s as the only other parameters; the manuscript does not discuss or test corrections for discrete-layer effects, direction-dependent adhesion, or the validity of replacing the atomic puckered lattice by a homogeneous anisotropic plate, which directly affects whether the reported numbers can be taken as the intrinsic directional moduli.

Authors: We agree that an explicit discussion of model assumptions would strengthen the paper. The continuum approximation is standard for buckling wavelengths much larger than the atomic scale, and we use an effective thickness t equal to the number of layers times the interlayer spacing. No direction-dependent adhesion was observed in our data. We will add a short subsection in the discussion addressing the applicability of the homogeneous anisotropic plate model to few-layer black phosphorus, citing supporting literature on continuum treatments of puckered lattices, and noting that discrete-layer corrections are expected to be negligible for the observed wavelengths. revision: yes

Circularity Check

0 steps flagged

No circularity: directional moduli extracted from measured ripple periods via external continuum buckling formula

full rationale

The paper reports experimental ripple periods under directional compression and applies the standard thin-film buckling wavelength relation (known from continuum mechanics on compliant substrates) to solve for the two Young's moduli. No equation in the provided text or abstract reduces the reported EbP values back to fitted inputs by construction, nor does any load-bearing step rely on self-citation chains or ansatzes imported from the authors' prior work. The derivation is self-contained against external benchmarks and the moduli are genuine outputs rather than renamed inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of the standard thin-film buckling wavelength formula to few-layer BP and on the assumption that measured periods directly encode only the directional Young's moduli.

free parameters (2)

Young's modulus (armchair)
Extracted from ripple period analysis; reported as 35.1 GPa with uncertainty.
Young's modulus (zig-zag)
Extracted from ripple period analysis; reported as 93.3 GPa with uncertainty.

axioms (1)

domain assumption Buckling ripple wavelength is governed by the standard continuum mechanics formula relating film modulus, thickness, and substrate compliance.
Invoked to convert observed periods into modulus values.

pith-pipeline@v0.9.0 · 5665 in / 1335 out tokens · 23210 ms · 2026-05-24T20:06:12.528033+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

anisotropic buckling stems from the puckered honeycomb crystal structure

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

[1]

1 L. Li, Y. Yu, G. J. Ye, Q. Ge, X. Ou, H. Wu, D. Feng, X. H. Chen and Y. Zhang, Nat. Nanotechnol., 2014, 9, 372–7. 2 F. Xia, H. Wang and Y. Jia, Nat. Commun., 2014, 5,

work page 2014
[2]

3 H. Liu, A. T. Neal, Z. Zhu, Z. Luo, X. Xu, D. Tománek and P. D. Ye, ACS Nano, 2014, 8, 4033–41. 4 A. Castellanos-Gomez, L. Vicarelli, E. Prada, J. O. Island, K. L. Narasimha-Acharya, S. I. Blanter, D. J. Groenendijk, M. Buscema, G. A. Steele, J. V. Alvarez, H. W. Zandbergen, J. J. Palacios and H. S. J. van der Zant, 2D Mater., 2014, 1, 025001. 5 M. Busc...

work page 2014
[3]

Gillgren, D

8 N. Gillgren, D. Wickramaratne, Y. Shi, T. Espiritu, J. Yang, J. Hu, J. Wei, X. Liu, Z. Mao, K. Watanabe, T. Taniguchi, M. Bockrath, Y. Barlas, R. K. Lake and C. Ning Lau, 2D Mater., 2014, 2, 011001. 9 L. Li, F. Yang, G. J. Ye, Z. Zhang, Z. Zhu, W.-K. Lou, X. Zhou, L. Li, K. Watanabe, T. Taniguchi, K. Chang, Y. Wang, X. H. Chen and Y. Zhang, Nat. Nanotec...

work page 2014
[4]

11 L. Li, J. Kim, C. Jin, G. J. Ye, D. Y. Qiu, H. Felipe, Z. Shi, L. Chen, Z. Zhang and F. Yang, Nat. Nanotechnol., 2017, 12,

work page 2017
[6]

22 Y. Deng, Z. Luo, N. J. Conrad, H. Liu, Y. Gong, S. Najmaei, P. M. Ajayan, J. Lou, X. Xu and P. D. Ye, ACS Nano, 2014, 8, 8292–8299. 23 H. Yuan, X. Liu, F. Afshinmanesh, W. Li, G. Xu, J. Sun, B. Lian, G. Ye, Y. Hikita, Z. Shen, S.-C. Zhang, X. Chen, M. Brongersma, H. Y. Hwang and Y. Cui, . 24 H. Chen, P. Huang, D. Guo and G. Xie, J. Phys. Chem. C, 2016,...

work page 2014
[7]

Favron, E

32 A. Favron, E. Gaufrès, F. Fossard, A.-L. Phaneuf-L’Heureux, N. Y.-W. Tang, P. L. Lévesque, A. Loiseau, R. Leonelli, S. Francoeur and R. Martel, Nat. Mater., 2015, 14, 826–32. 33 A. A. Kistanov, Y. Cai, K. Zhou, S. V Dmitriev and Y.-W. Zhang, 2D Mater., 2016, 4, 15010. 34 A. A. Kistanov, Y. Cai, K. Zhou, S. V Dmitriev and Y.-W. Zhang, J. Phys. Chem. C, ...

work page 2015
[9]

Iguiñiz, R

45 N. Iguiñiz, R. Frisenda, R. Bratschitsch and A. Castellanos‐Gomez, Adv. Mater., 2019, 1807150. 46 P. Feicht, R. Siegel, H. Thurn, J. W. Neubauer, M. Seuss, T. Szabó, A. V Talyzin, C. E. Halbig, S. Eigler and D. A. Kunz, Carbon N. Y., 2017, 114, 700–705. 47 D. A. Kunz, J. Erath, D. Kluge, H. Thurn, B. Putz, A. Fery and J. Breu, ACS Appl. Mater. Interfac...

work page doi:10.1021/acs.nanolett.5b04670 2019
[12]

Method and conditions (environment, thickness) E [GPa] Ref

Summary of the reported valued (both theoretical and experimental) in the literature, indicating the method and conditions employed to obtain them. Method and conditions (environment, thickness) E [GPa] Ref. AC ZZ Theory Density functional theory 41.3 106.4 56 Density functional theory 37-44 159-166 55 Density functional theory 52.3 191.9 54 Molecular dyn...

work page doi:10.1039/c9nr03009c 2019

[1] [1]

1 L. Li, Y. Yu, G. J. Ye, Q. Ge, X. Ou, H. Wu, D. Feng, X. H. Chen and Y. Zhang, Nat. Nanotechnol., 2014, 9, 372–7. 2 F. Xia, H. Wang and Y. Jia, Nat. Commun., 2014, 5,

work page 2014

[2] [2]

3 H. Liu, A. T. Neal, Z. Zhu, Z. Luo, X. Xu, D. Tománek and P. D. Ye, ACS Nano, 2014, 8, 4033–41. 4 A. Castellanos-Gomez, L. Vicarelli, E. Prada, J. O. Island, K. L. Narasimha-Acharya, S. I. Blanter, D. J. Groenendijk, M. Buscema, G. A. Steele, J. V. Alvarez, H. W. Zandbergen, J. J. Palacios and H. S. J. van der Zant, 2D Mater., 2014, 1, 025001. 5 M. Busc...

work page 2014

[3] [3]

Gillgren, D

8 N. Gillgren, D. Wickramaratne, Y. Shi, T. Espiritu, J. Yang, J. Hu, J. Wei, X. Liu, Z. Mao, K. Watanabe, T. Taniguchi, M. Bockrath, Y. Barlas, R. K. Lake and C. Ning Lau, 2D Mater., 2014, 2, 011001. 9 L. Li, F. Yang, G. J. Ye, Z. Zhang, Z. Zhu, W.-K. Lou, X. Zhou, L. Li, K. Watanabe, T. Taniguchi, K. Chang, Y. Wang, X. H. Chen and Y. Zhang, Nat. Nanotec...

work page 2014

[4] [4]

11 L. Li, J. Kim, C. Jin, G. J. Ye, D. Y. Qiu, H. Felipe, Z. Shi, L. Chen, Z. Zhang and F. Yang, Nat. Nanotechnol., 2017, 12,

work page 2017

[5] [6]

22 Y. Deng, Z. Luo, N. J. Conrad, H. Liu, Y. Gong, S. Najmaei, P. M. Ajayan, J. Lou, X. Xu and P. D. Ye, ACS Nano, 2014, 8, 8292–8299. 23 H. Yuan, X. Liu, F. Afshinmanesh, W. Li, G. Xu, J. Sun, B. Lian, G. Ye, Y. Hikita, Z. Shen, S.-C. Zhang, X. Chen, M. Brongersma, H. Y. Hwang and Y. Cui, . 24 H. Chen, P. Huang, D. Guo and G. Xie, J. Phys. Chem. C, 2016,...

work page 2014

[6] [7]

Favron, E

32 A. Favron, E. Gaufrès, F. Fossard, A.-L. Phaneuf-L’Heureux, N. Y.-W. Tang, P. L. Lévesque, A. Loiseau, R. Leonelli, S. Francoeur and R. Martel, Nat. Mater., 2015, 14, 826–32. 33 A. A. Kistanov, Y. Cai, K. Zhou, S. V Dmitriev and Y.-W. Zhang, 2D Mater., 2016, 4, 15010. 34 A. A. Kistanov, Y. Cai, K. Zhou, S. V Dmitriev and Y.-W. Zhang, J. Phys. Chem. C, ...

work page 2015

[7] [9]

Iguiñiz, R

45 N. Iguiñiz, R. Frisenda, R. Bratschitsch and A. Castellanos‐Gomez, Adv. Mater., 2019, 1807150. 46 P. Feicht, R. Siegel, H. Thurn, J. W. Neubauer, M. Seuss, T. Szabó, A. V Talyzin, C. E. Halbig, S. Eigler and D. A. Kunz, Carbon N. Y., 2017, 114, 700–705. 47 D. A. Kunz, J. Erath, D. Kluge, H. Thurn, B. Putz, A. Fery and J. Breu, ACS Appl. Mater. Interfac...

work page doi:10.1021/acs.nanolett.5b04670 2019

[8] [12]

Method and conditions (environment, thickness) E [GPa] Ref

Summary of the reported valued (both theoretical and experimental) in the literature, indicating the method and conditions employed to obtain them. Method and conditions (environment, thickness) E [GPa] Ref. AC ZZ Theory Density functional theory 41.3 106.4 56 Density functional theory 37-44 159-166 55 Density functional theory 52.3 191.9 54 Molecular dyn...

work page doi:10.1039/c9nr03009c 2019