[Paper Review] Anisotropic Young's Modulus for Single-Layer Black Phosphorus: The Third Principle Direction Besides Armchair and Zigzag

Jin-Wu Jiang|arXiv (Cornell University)|Mar 1, 2015

2D Materials and Applications3 citations

TL;DR

This paper derives an analytic formula for the directional dependence of Young’s modulus in single-layer black phosphorus using the valence force field model, revealing a previously unreported third principal direction at φₜₚ = 0.268π where the Young’s modulus reaches a maximum of 111.4 Nm⁻¹—exceeding values in both armchair (52.2 Nm⁻¹) and zigzag (85.4 Nm⁻¹) directions. The result highlights a new anisotropic mechanical axis critical for understanding mechanical and electronic properties in this 2D material.

ABSTRACT

We derive an analytic formula for the Young's modulus in single-layer black phosphorus using the valence force field model. By analyzing the directional dependence for the Young's modulus, we explore the third principle direction with direction angle phi_tp = 0.268pi besides armchair and zigzag directions. The maximum Young's modulus value is in the third principle direction. More specifically, the Young's modulus is 52.2 N/m, 85.4 N/m, and 111.4 N/m in the armchair direction, zigzag direction, and the third principle direction, respectively. This new principle direction is of significance for future discussions of other anisotropic properties in the single-layer black phosphorus.

Motivation & Objective

• To identify and characterize the full directional dependence of Young’s modulus in single-layer black phosphorus beyond the conventional armchair and zigzag directions.
• To resolve the limitation in prior studies that only compared two principal directions, missing a third maximum in stiffness.
• To derive an analytic expression for Young’s modulus using the valence force field model to enable explicit understanding of anisotropic mechanical behavior.
• To establish the significance of the third principal direction for future studies on mechanical, electronic, and optical anisotropy in black phosphorus.

Proposed method

• Employed the valence force field model (VFFM) with nine potential terms to describe bond stretching, angle bending, and dihedral interactions in single-layer black phosphorus.
• Derived an analytic formula for directional Young’s modulus by computing the second derivative of the potential energy with respect to strain in a given direction.
• Calculated geometric coefficients (e.g., α₁, α₃, xᵢⱼₚ) from atomic coordinates and bond parameters to express the strain energy in terms of direction angle φ.
• Solved for extreme points of Young’s modulus by setting the derivative of the modulus expression to zero, identifying three critical directions.
• Used experimental lattice parameters: a₁ = 4.376 Å, a₂ = 3.314 Å, d₁ = 2.2449 Å, d₂ = 2.2340 Å, θ₁ = 0.535π, and θ₂ = 0.567π.
• Validated results by comparing contributions from individual VFFM terms and confirming the dominant role of the bond-stretching term (Vᵣ) in determining the directional anisotropy.

Experimental results

Research questions

• RQ1 What is the full directional dependence of Young’s modulus in single-layer black phosphorus beyond the armchair and zigzag directions?
• RQ2 Does a third principal direction exist where the Young’s modulus reaches a maximum, and if so, what is its orientation?
• RQ3 How does the third principal direction relate to the bond angles and atomic geometry, particularly in relation to θ₁/2?
• RQ4 Why do previous studies report only two principal directions, and what is the physical origin of the missing maximum stiffness direction?
• RQ5 Can an analytic formula derived from the valence force field model accurately predict the mechanical anisotropy in black phosphorus?

Key findings

• The Young’s modulus reaches a maximum value of 111.4 Nm⁻¹ in the third principal direction at φₜₚ = 0.268π, exceeding both armchair (52.2 Nm⁻¹) and zigzag (85.4 Nm⁻¹) directions.
• The third principal direction is not aligned with the bond direction r₂₃, despite being very close, and only coincides with θ₁/2 when θ₁ = 0.5π, which is not the case for real black phosphorus (θ₁ = 0.535π).
• The armchair direction (φ = 0) exhibits the minimum Young’s modulus, contradicting the common assumption that the minimum occurs at the direction of maximum strain.
• The maximum stiffness direction is primarily governed by the bond-stretching potential term (Vᵣ), which dominates the directional anisotropy in Young’s modulus.
• The derived analytic expression for Young’s modulus enables explicit, closed-form understanding of mechanical anisotropy without relying on numerical simulations.
• The third principal direction is a critical axis for future studies on mechanical, electronic, and optical properties in single-layer black phosphorus, as it represents the stiffest mechanical response.