Auxetic B4N Monolayer: A promising 2D material ... - ACS Publications

fundamental mechanical responses of solids against external loads.1 Most materials ..... minimum value of 80 N/m at α=45o from a value of 153 N/m at ...
0 downloads 0 Views 1MB Size
Subscriber access provided by Nottingham Trent University

Functional Nanostructured Materials (including low-D carbon)

Auxetic B4N Monolayer: A promising 2D material with In-Plane Negative Poisson’s Ratio and Large Anisotropic Mechanics Bing Wang, Qisheng Wu, Yehui Zhang, Liang Ma, and Jinlan Wang ACS Appl. Mater. Interfaces, Just Accepted Manuscript • Publication Date (Web): 22 Aug 2019 Downloaded from pubs.acs.org on August 22, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

Auxetic B4N Monolayer: A Promising 2D Material with In-Plane Negative Poisson’s Ratio and Large Anisotropic Mechanics Bing Wang,† Qisheng Wu, † Yehui Zhang, Liang Ma,* and Jinlan Wang* School of Physics, Southeast University, Nanjing 211189, P. R. China

ABSTRACT: Auxetic materials, known with negative Poisson’s ratio, are highly desirable for many advanced applications but the candidates are rather scarce, especially at low dimension. Motivated by the re-entrant structure that often exposes negative Poisson’s ratio, we predict a two-dimensional (2D) planar B4N monolayer as a promising auxetic material with unusual in-plane negative Poisson’s ratio within the framework of density functional theory calculations. B4N monolayer also exhibits a highly emerged mechanical anisotropy, characterized by Young’s modulus and Poisson’s ratio. In addition, this monolayer shows superior mechanical flexibility in ideal tensile strength and critical strain values. The phonon dispersion calculations and ab initio molecular dynamics simulations further demonstrate that this monolayer also owns excellent dynamical and thermal stabilities. The fantastic mechanical properties coupled with robust structural stability render the auxetic B4N monolayer promising for future nanomechanical devices. KEYWORDS: Auxetic materials, re-entrant structure, mechanical property, B4N monolayer, first-principles method.

1

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 18

INTRODUCTION Poisson’s ratio (ν) of a material is defined as the ratio of the transverse contraction strain to the longitudinal extension strain, which measures the fundamental mechanical responses of solids against external loads.1 Most materials expose positive Poisson’s ratio, which means that when experiencing a longitudinal tension (or compression) they will become thinner (or thicker) in the lateral directions.2-3 In contrast, a small group of materials with negative Poisson’s ratio (NPR), namely auextics or auxetic materials, exhibit a counterintuitive behavior in contrast to that of conventional materials, which expand when stretched and contract when compressed.4 The first auxetic material was discovered in 1982,5 and since then, auxetic materials have attracted world-wide attention because of their unique mechanical properties5-8, and highly demanding in many advanced applications, such as tissue engineering9, medicine10, tougher composites11, fasteners12, bulletproof vests, aircraft13, national security, and so on 14-16. Within the last 30 years, a variety of NPR materials and structures have been discovered, synthesized or fabricated, ranging from macroscopic down to molecular levels.17 Although most reported auxetic materials are in bulk form, the negative Poisson’s ratio is possible in two-dimensional (2D) structures according to the theory of Gibson.5 As the blooming of 2D family, several systems with NPR have been predicted theoretically18-29, and only one of them, black phosphorus19-20, has been observed experimentally. To date, the intrinsic 2D auxetic materials are still rather limited, far from meeting the wide range of needs in practice. Seeking for new 2D auextic materials is still a pressing task for modern nanoscale devices. Re-entrant structure is one of the typical auxetic structures, which has been widely used to model the deformation of auxetic foam by many researchers.2,

5, 15, 30-36

Lakes et al.

33

reported the first man-made auxetic material with a re-entrant foam structure. Wang et al.30 designed a novel re-entrant auxetic honeycomb with the in-plane impact responses computationally. Jin et al.

34

and Chang et al.

35

used the re-entrant

structure as sandwich structure core to improve the blast resistance and dynamic 2

ACS Paragon Plus Environment

Page 3 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

response. Although this structure is allowable in traditional engineered materials, to our knowledge, 2D auxetic material with re-entrant structure has yet to come. Very recently, various 2D boron monolayers (BMs) were predicted theoretically 37-41,

and the χ3 phase (χ3-BM) was successfully synthesized on Ag(111) substrate by

two different groups.42 It is interesting to note that a re-entrant-like structure can be extracted from the χ3-BM. However, the expected NPR is actually absent without negative inclined and vertical strut angle. In this work, we predict a new planar auxetic 2D B4N monolayer, which can be viewed as nitrogen embedded χ3-BM. Our density functional theory (DFT) computations show that the B4N monolayer exposes large anisotropic NPR and high in-plane Young’s modulus. It also owns ultrahigh tensile strains (19% and 18% along x and y directions, respectively), which are comparable to those of graphene. Moreover, the planar B4N monolayer shows excellent thermal, dynamic, and mechanical stabilities as well. RESULTS AND DISCUSSION The typical hexagonal and re-entrant honeycomb structures, as shown in Figure 1a and 1b, respectively, can be numerically distinguished by the length of the horizontal strut h, the length of the inclined strut l, and the angle θ (ϵ[-π/2, π/2]) between the inclined and vertical strut. In such a representation, the Poisson's ratio along x direction (νx) for both hexagonal and re-entrant honeycomb structure with mechanical loading along y-direction, as an example, can be expressed by the Eq. (1):43-45 𝜈𝑥 = ―

𝜀𝑦 𝜀𝑥

=

(ℎ + 𝑙sin 𝜃)sin 𝜃 𝑙(cos 𝜃)2

(1)

Markedly, the sign of νx is determined by the sinθ term in the numerator, and negative angle θ (ϵ [-π/2, 0]) will produce negative νx, that is, the NPR effect. As can be clearly seen from Figure 1, the angle θ is always positive for the typical hexagonal structure but negative for the re-entrant structure, which means that the NPR effect would definitely appear in the re-entrant structure of B4N monolayer but not the regular hexagonal structure. 3

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 18

Figure 1. Sketch for the hexagonal (a) and re-entrant (b) honeycomb structure. The lengths (in angstroms) of the representative chemical bonds of χ3-BM (c), and the atomic configuration for B4N monolayer (d). The unit cells are denoted with the black lines. The blue and green balls represent N and B atoms, respectively. The synthesized χ3-BM consists of waved and double-chained boron atoms that connect each other with regular hexagonal honeycomb hole in-between (Figure 1c). Geometrically, the angle θ for χ3-BM is always positive due to its hexagonal honeycomb structure. Thus, only positive Poisson’s ratio is expected for the χ3-BM monolayer. It is known that B and N atoms can also form strong B-N bonds but their lengths (~1.45 Å in h-BN) would be much shorter than the B-B bonds. It would be very interesting to figure out whether the honeycomb structure can be altered by introducing shorter B-N bonds into the 2D χ3-BM monolayer and even give rise to the negative angle θ with charming NPR. To this end, we design a planar 2D re-entrant structure by injecting N atoms into χ3-BM, namely B4N monolayer, as shown in Figure 1d. After structure relaxation, it is shown that the length of B-N bonds (1.33 Å) in B4N monolayer is much shorter than the B-B bonds in χ3-BM (1.64 ~ 1.72 Å, Figure 1c) and even shorter than that in h-BN ~1.45 Å, which is able to convert the honeycomb structure into a stable re-entrant structure with NPR. More detailed 4

ACS Paragon Plus Environment

Page 5 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

discussion will be given in the following sections. The structure exposes Cmmm symmetry (space group No. 65) with orthogonal lattice. The optimized lattice parameters are a=2.9693 Å and b=10.7035 Å with B1 at 4j (0, 0.2212, 1/2), B2 at 4j (1/2, 0.1246, 1/2), and N at 2c (1/2, 0, 1/2), respectively. N atoms act as ‘bridges’ by bonding with two B atoms in B4N monolayer. Young’s modulus (E), reflecting the flexibility or rigidity of materials, and Poisson’s ratio (ν), reflecting the mechanical responses of solids against external loads, are two main mechanical parameters of a material. For 2D linear elastic solid materials, E and ν along x and y directions can be derived from the elastic constants as below: Ex = (C11C22-C12C21)/C22,

(2)

νx = C21/C22,

(4)

Ey = (C11C22-C12C21)/C11, νy = C12/C11,

(3) (5)

Table 1. The elastic constants Cij (N/m), Young’s modulus E (N/m), and Poisson’s ratio ν of B4N, χ3-BM, graphene, h-BN, and MoS2. The present results are also compared with previous literature.

structures

reference

C11

C22

C12

C66

Ex

Ey

νx

νy

B4N

Our work

153

267

-4.9

28

153

268

-0.018

-0.032

χ3-BM

Ref. 38

196

208

0.11

0.12

334

334

0.206

0.206

340

340

332

332

Graphene

h-BN MoS2

Our work

349

349

72

138

Ref. 46 Ref. 26

343

343

62

Our work

290

290

64

112

276

276

0.220

0.220

Ref. 47

289.8

289.8

63.7

113.1

275.8

275.8

0.220

0.220

Our work

128.8

127.9

28.8

50

122.3

121.5

0.224

0.225

123

123

0.250

Ref. 48 Ref

46:

Experiment; Ref

26, 38, 47, 48:

Theory.

The calculated mechanical properties of B4N monolayer are shown in Table 1, together with those of other well-studied 2D materials for comparison. The calculated results for known 2D materials are in good agreement with previous experimental/theoretical reports 26, 38, 46-48, demonstrating the reliability of our present 5

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 18

theoretical method. Based on eqs. (2) and (3), the in-plane Young’s moduli along x (Ex) and y (Ey) directions of B4N monolayer are calculated to be 153 N/m and 268 N/m, respectively, exposing strong anisotropy. To gain a comprehensive description of the anisotropic mechanical properties of B4N monolayer, the in-plane Young’s moduli as functions of the arbitrary direction α (α is the angle relative to the positive x direction in this monolayer) were calculated and presented in 2D polar representation curve (Figure 2a). Obviously, the Young’s moduli of B4N monolayer are highly anisotropic in the whole plane. The Young’s modulus E(α) first decreases to its minimum value of 80 N/m at α=45o from a value of 153 N/m at α=0°(x direction) and then increases to a maximum value of 268 N/m at α=90o (y direction). The maximum value (268 N/m) is comparable to that of h-BN monolayer (275.8 N/m)47 and graphene (340 N/m)46 and is about two times of that of MoS2 (123 N/m) 48.

Figure 2. Calculated orientation-dependent (a) Young’s modulus E(α) in N/m and (b) Poisson’s ratio ν(α). (c) Total energy with respect to lattice response of the other direction when the B4N monolayer lattice is under 5% tensile strain along x and y directions, respectively. The arrows indicate the equilibrium magnitude of εx and εy, respectively. (d) The stress-strain relationship of B4N/graphene along x (0o)/zigzag and y (90o)/armchair directions. The arrows denote the strain at the maximum stress. 6

ACS Paragon Plus Environment

Page 7 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

The diagram of calculated Poisson’s ratio as functions of the arbitrary direction shows that the Poisson’s ratio of B4N monolayer is also spatially varying and highly anisotropic (Figure S1). More remarkably, an intriguing unconventional NPR phenomenon is observed in this Figure for B4N monolayer along both x and y directions. Due to the small values of negative Poisson’s ratio, a zoom of Figure S1 near the center point is plotted in Figure 2b. According to eqs. (4) and (5), the NPRs of B4N monolayer along x and y direction are -0.018 and -0.032, respectively, due to the negative value of C12. The Poisson’s ratio calculated according to the uniaxial strain method is also in good agreement with the results obtained through elastic constants (Figure S2). To confirm this interesting point, we also applied a uniaxial strain of 5% in x and y direction of B4N monolayer, respectively. Just as expected, we find that the equilibrium lattice constants of y and x directions are expanded by ~0.28% and 0.52% (Figure 2c), respectively. The calculated values of NPR are comparable to that of Pmmn-borophene21 (νx= -0.04, νy= -0.02) and black phosphorus19 (−0.027 along z direction). Along the directions other than x and y, B4N monolayers expose the normal positive Poisson’s ratio. The NPR value indicates that B4N monolayer will actually expand in the transverse direction when stretched along x/y direction. In addition to Young’s modulus and Poisson’s ratio, the ideal strength is another important mechanical property of material. We thus investigate the ideal strength of B4N monolayer by calculating the stress-strain curves of x and y directions as depicted in Figure 2d. With applied small strains, the sheet exhibits linear stress–strain relationship with notable elastic anisotropy. As the strain increases, the stress–strain behavior becomes more and more nonlinear and shows enhanced anisotropy. The maximum stress along y direction is higher than that along x direction. The peak stress reaches 25 N/m along the ridges at εy = 18% and 11 N/m at εx = 19% across the ridges, which is secondary only to graphene (34 N/m)

38

but higher than black phosphorus

(16 N/m)49. For comparison, our calculated maximum tensile strains for graphene is 19% and 25% along armchair and zigzag directions, respectively, in good agreement 7

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 18

with the previous report50 (19.4% and 26.6% along armchair and zigzag directions), further confirming the reliability of our computational method.

Figure 3. (a) The lengths (in angstroms) of the representative chemical bonds of B4N monolayer. (b) Charge density of B4N monolayer with the iso-surface of 0.99 e/Å3. (c) Explanation of negative Poisson’s ratio via hinging of ribs in response to applied stretching in x direction. (d) Band structure and 2D Brillouin zone of B4N monolayer. To probe the origin of the NPR, we further analyze the bond lengths of B4N monolayer (Figure 3a). Take the left part of B4N monolayer as an example, the bond lengths of B4-B5 (1.65 Å), B2-B5 (1.61 Å) and B2-B7 (1.65 Å) are shorter than those of B4-B2 (1.81 Å) and B5-B7 (1.81 Å) by ~9% - 11%. The shorter lengths indicate the B4-B5, B2-B5, and B2-B7 B-B bonds are stronger than those B4-B2 and B5-B7 B-B bonds. These bonding strength differences can give rise to the typical re-entrant honeycomb-like structure with negative angle  and the NPR phenomenon (ν < 0). The charge density of B4N monolayer is also plotted to elucidate its bond properties shown in Figure 3b. The higher charge density between B4-B5, B5-B2, and B2-B7 also indicate stronger bond strength than that between B4-B2 and B5-B7 bonds, which is in good agreement with the analysis of bond lengths. Figure 3c displays the detailed structure deformation along the hinging of the diagonal ribs for re-entrant structure in response to an applied uniaxial load in x direction. Under a global tensile loading in x direction (red arrows), the atoms B5 and B6 of this re-entrant structure 8

ACS Paragon Plus Environment

Page 9 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

move apart, while atoms B2 and B1 maintain their separation. Meanwhile, the tensile strain (red arrows) in x direction of B4N monolayer causes the ribs (B2-B5-B2 and B1-B6-B1 chains) aligned along y direction to separate more (black arrows). As a result, the overall structure expands in y direction, and its area is also apparently expanded while the structure is stretched, leading to negative Poisson’s ratio. Our analyses show that the auxeticity of our novel B4N monolayer is mainly originated from its particular atomic arrangement. We also calculate the projected band structures and partial density of states to probe the electronic properties of B4N monolayer. As shown in Figure 3d, the band lines of B4N monolayer travel across its Fermi level, indicating that it is a nonmagnetic metal. The PDOS analysis shows that the B-p states and N-p states dominate around the Fermi level, which is the principal cause of metallicity (Figure S3). The revealed apparent hybridization between both B-p and N-p states are consistent with the strong B-N bonds. The orbital electronic band structure shows the main contribution to the states around the Fermi level arises from B_pz and N_px orbitals. The metallic property is also maintained and the band structures change slightly under the different external strains (Figure S4).

9

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 18

Figure 4. (a) Phonon dispersion of B4N monolayer. (b) Fluctuation of total energy during AIMD simulation of B4N monolayer at 500 K. The inset is the structure of B4N monolayer at the end of the AIMD simulation after 5 ps. (c) Schematic illustration of preparing B4N. To assess the dynamical stability of B4N monolayer, we examined its phonon dispersion as shown in Figure 4a. The absence of imaginary modes confirms dynamically stability of the B4N monolayer. The highest frequency of B4N monolayer reaches up to 50.84 THz (1695 cm-1), which is comparable to that of graphene (1586 cm-1), 51 indicating the robust B-B and B-N interactions in B4N monolayer. Ab initio molecular dynamics (AIMD) simulations were also conducted to probe the thermal stability of B4N monolayer. Figure 4b shows that the average value of the total potential energy during the AIMD simulations is oscillating within a narrow range and the integral configuration of B4N monolayer is well maintained after 5 ps AIMD simulation, confirming its thermal stability. The mechanical stability was also investigated by using the strain-stress method. For a mechanically stable 10

ACS Paragon Plus Environment

Page 11 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

orthorhombic 2D sheet, the elastic constants need to satisfy that C11C22-C12C21 > 0 and C66 > 0. 47 Considering the symmetry of B4N, we have C12=C21. Table 1 shows that the calculated elastic constants of B4N monolayer satisfy stability criteria, confirming the mechanical stability of this sheet. Therefore, we can conclude that the B4N monolayer has very good dynamical, thermal, and mechanical stability. Although the demonstrated dynamical and thermal stabilities have already suggested great feasibility to synthesize B4N monolayers, it would be much great to propose any practical synthesis approach. The χ3-BM was successfully synthesized on Ag(111) surface by molecular beam epitaxy.27 In the experimental process, one-dimensional boron chains or nanoribbons were also realized in the process of BM synthesis.17,

39, 41, 21, 52

If the N atoms are injected among the boron chains or

nanoribbons, the strong interaction between B and N atoms may form the B4N monolayer with re-entrant structure. A possible synthesis route for 2D B4N monolayer is therefore proposed in Figure 4c. To evaluate the relative stability between the χ3 B sheet and B4N monolayer, we calculated the formation energy, which is defined as Ef = (EB4N - 4μχ3B - 1/2μN2)/5. EB4N is the total energy per formula unit, and μχ3B and μN2 are chemical potentials of 2D χ3-BM and nitrogen, respectively. The calculated formation energy for B4N monolayer is -0.16 eV per atom. The negative formation energy suggests that the synthesis of the B4N monolayer is exothermic. We also performed a global search for the lowest-energy stable structure for B4N monolayer by using the particle-swarm optimization method as implemented in CALYPSO code.53 The results show that our predicted 2D monolayer is the metastable state which has very close energy to the most stable one (only 18 meV/atom high). More details are given in Supporting Information (Figure S5). Considering the excellent dynamical and thermal stability, the experimental fabrication of B4N monolayer might be feasible although it would be very challenging. CONCLUSIONS To summarize, we have successfully predicted a 2D auxetic B4N monolayer material with attractive negative Poisson’s ratio by using first-principles calculations. 11

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 18

It is revealed that the auxeticity of B4N monolayer mainly originates from its particular re-entrant structure. The intriguing robust in-plane NPR and high Young’s modulus of auxetic B4N monolayer are both highly anisotropic. Besides, B4N monolayer can sustain tensile strains of up to 19% and 18% with maximum stresses of 11 and 25 N/m in the x and y directions, respectively, which is comparable to that of graphene. The thermal, dynamical, and mechanical stability evaluations demonstrate the feasibility of the synthesis of this 2D B4N monolayer and a viable experimental route for the synthesis of B4N monolayer is proposed as well. In view of the novel properties, we wish this new monolayer could be synthesized in near future and applied to the nanomechanical devices. COMPUTATIONAL METHODS. The structure relaxations and self-consistent energy calculations were carried out by using Perdew-Burke-Ernzerhof parameterization of the generalized gradient approximation (PBE-GGA)54 as implemented in the Vienna ab initio simulation package (VASP)55. The projector-augmented plane wave (PAW) approach56 was used to describe the ion-electron interaction. The energy cutoff of the plane wave was set to 500 eV with the energy precision of 10−6 eV and the force precision of 10−2 eV/Å. The Brillouin zone was sampled with a 21 × 7 × 1 Γ-centered k point grid. Phonon dispersions were calculated based on density functional perturbation theory (DFPT) as embedded in phonopy program.57 Van der Waals interactions were considered by using the DFT-D2 method.58 The vacuum length was set to 20 Å. The ab initio molecular dynamics (AIMD) simulations were carried out to evaluate the thermal stabilities of B4N monolayer. At each temperature, AIMD simulation lasted for 5 ps in NVT ensemble with a time step of 1 fs. The temperature was controlled by using the Nosé-Hoover method59. The orientation-dependent Young’s modulus E(α) and Poisson’s ratio ν(α) were calculated as,28, 60, 61 𝐶11𝐶22 ― 𝐶212

𝐸(α) = 4

4

𝐶11𝑠 + 𝐶22𝑐 +

(

𝐶11𝐶22 ― 𝐶212 𝐶44

)

― 2𝐶12 𝑐2𝑠2 12

ACS Paragon Plus Environment

Page 13 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

ν(α) =

𝐶11𝐶12 ― 𝐶212 2 2 (𝐶11 + 𝐶22 ― )𝑐 𝑠 ― 𝐶12(𝑐4 + 𝑠4) 𝐶44 4

4

𝐶11𝑠 + 𝐶22𝑐 +

(

𝐶11𝐶22 ― 𝐶212 𝐶44

)

― 2𝐶12 𝑐2𝑠2

where c = cosα, and s = sinα. ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI:XXX The

details

about

the

derivation

process

of

Equation

1;

alculated

orientation-dependent the Poisson’s ratio ν(θ) in xy plane (Figure S1); Correlations between applied strain along the x (y) direction and the resultant transverse strains along the y (x) direction (Figure S2); The projected of density of state for B4N monolayer (Figure S3); Band structure of B4N monolayer under different strains (Figure S4); The structure ponon dispersion and B4N monolayer (Figure S5). (PDF) AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]; [email protected], Author Contributions †Bing Wang and Qisheng Wu contributed equally to this work. Notes The authors declare no competing financial interest. ACKNOWLEDGMENTS This work is supported by the National Key Research and Development Program of China (2017YFA0204800), Natural Science Funds of China (21773027,21525311), Jiangsu 333 project (BRA2016353), the Fundamental Research Funds for the Central Universities (2242019R10021), the Scientific Research Foundation of Graduate School of Southeast University (YBJJ1773). The authors thank the computational resources at the SEU and National Supercomputing Center in Tianjin.

13

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 18

References 1. Greaves, G. N.; Greer, A. L.; Lakes, R. S.; Rouxel, T., Poisson's Ratio and Modern materials. Nat. Mater. 2011, 10, 823-837. 2. Alderson, K. E. E. a. A., Auxetic Materials: Functional Materials and Structures from Lateral Thinking. Adv. Mater. 2000, 12, 617-628. 3. Lakes, R., Advances in Negative Poisson’s Ratio Materials Material. Adv. Muter. 1993, 4, 293-296. 4. Lakes, R., Foam Structures with a Negative Poisson's Ratio. Science 1987, 235, 1038-1040. 5. Gibson, L. J.; Ashby, M.; Schajer, G.; Robertson, C., the Mechanics of Two-Dimensional Cellular Materials, Proc. R. Soc. Lond. A. 1982, 382, 25-42. 6. Choi, J. B.; Lakes, R. S., Non-linear Properties of Metallic Cellular Materials with a Negative Poisson's Ratio. J. Mater. Sci. 1992, 27, 5375-5381. 7. Lakes, R. S.; Elms, K., Indentability of Conventional and Negative Poisson’s Ratio Foams. J. Comp. Mater. 1993, 27, 1193-1202. 8. Choi, J. B.; Lakes, R. S., Fracture Toughness of Re-entrant Foam Materials with A Negative Poisson's Ratio: Experiment and Analysis. Int. J. Fracture 1996, 80, 73-83. 9. Park, Y. J.; Kim, J. K., The Effect of Negative Poisson’s Ratio Polyurethane Scaffolds for Articular Cartilage Tissue Engineering Applications. Adv. Mater. Sci. Eng. 2013, 2013, 1-5. 10. Scarpa, F., Auxetic Materials for Bioprostheses. IEEE Signal Proc. Mag. 2008, 25, 128-126. 11. Sun, Y.; Pugno, N., Hierarchical Fibers with A Negative Poisson’s Ratio for Tougher Composites. Materials 2013, 6, 699-712. 12. Xu, X.; Zhang, Q.; Hao, M.; Hu, Y.; Lin, Z.; Peng, L.; Wang, T.; Ren, X.; Wang, C.; Zhao, Z., Double-Negative-Index Ceramic Aerogels for Thermal Superinsulation. Science 2019, 363, 723-727. 13. Baughman, R. H.; Shacklette, J. M.; Zakhidov, A. A.; Stafström, S., Negative Poisson’s Ratios As a Common Feature of Cubic Metals. Nature 1998, 392, 362-365. 14. Yang, W.; Li, Z. M.; Shi, W.; Xie, B. H.; Yang, M. B., Review on Auxetic Materials. J. Mater. Sci. 2004, 39, 3269-3279. 15. Prawoto, Y., Seeing Auxetic Materials from the Mechanics Point of View: A Structural Review on the Negative Poisson’s Ratio. Comp. Mater. Sci. 2012, 58, 140-153. 16. Yin, H.; Liu, C.; Zheng, G. P.; Wang, Y.; Ren, F., Ab Initio Simulation Studies on the Room-temperature Ferroelectricity in Two-Dimensional β-Phase GeS. Appl. Phys. Lett. 2019, 114, 192903. 17. Evans, K. E.; Alderson, A., Auxetic Materials: Functional Materials and Structures from Lateral Thinking! Adv. Mater. 2000, 12, 617-628. 18. Zhang, C.; He, T.; Matta, S. K.; Liao, T.; Kou, L.; Chen, Z.; Du, A., Predicting Novel 2D MB2 (M= Ti, Hf, V, Nb, Ta) Monolayers with Ultrafast Dirac Transport Channel and Electron-Orbital Controlled Negative Poisson’s Ratio. J. Phys. Chem. Lett. 2019, 10, 2567-2573. 19. Jiang, J. W.; Park, H. S., Negative Poisson's Ratio in Single-Layer Black Phosphorus. Nat. Commun 2014, 5, 4727. 20. Du, Y.; Maassen, J.; Wu, W.; Luo, Z.; Xu, X.; Ye, P. D., Auxetic Black Phosphorus: A 14

ACS Paragon Plus Environment

Page 15 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

2D Material with Negative Poisson's Ratio. Nano Lett. 2016, 16, 6701-6708. 21. Mannix, A. J.; Zhou, X. F.; Kiraly, B.; Wood, J. D.; Alducin, D.; Myers, B. D.; Liu, X.; Fisher, B. L.; Santiago, U.; Guest, J. R.; Yacaman, M. J.; Ponce, A.; Oganov, A. R.; Hersam, M. C.; Guisinger, N. P., Synthesis of Borophenes: Anisotropic, Two-Dimensional Boron Polymorphs. Science 2015, 350, 1513-1516. 22. Özçelik, V. O.; Cahangirov, S.; Ciraci, S., Stable Single-Layer Honeycomblike Structure of Silica. Phys. Rev. Lett. 2014, 112, 246803. 23. Qiao, M.; Wang, Y.; Li, Y.; Chen, Z., Tetra-Silicene: A Semiconducting Allotrope of Silicene with Negative Poisson’s Ratios. J. Phys. Chem. C 2017, 121, 9627−9633. 24. Zhang, S.; Zhou, J.; Wang, Q.; Chen, X.; Kawazoe, Y.; Jena, P., Penta-Graphene: A New Carbon Allotrope. PNAS 2015, 112, 2372-2377. 25. Yagmurcukardes, M.; Sahin, H.; Kang, J.; Torun, E.; Peeters, F. M.; Senger, R. T., Pentagonal Monolayer Crystals of Carbon, Boron Nitride, and Silver Azide. J. Appl. Phys. 2015, 118, 104303. 26. Wang, Y.; Li, F.; Li, Y.; Chen, Z., Semi-Metallic Be5C2 Monolayer Global Minimum with Quasi-Planar Pentacoordinate Carbons and Negative Poisson's Ratio. Nat. Commun. 2016, 7, 11488. 27. Gao, Z.; Dong, X.; Li, N.; Ren, J., Novel Two-Dimensional Silicon Dioxide with In-Plane Negative Poisson's Ratio. Nano Lett. 2017, 17, 772-777. 28. Peng, R.; Ma, Y.; He, Z.; Huang, B.; Kou, L.; Dai, Y., Single-Layer Ag2S: A Two-Dimensional Bidirectional Auxetic Semiconductor. Nano Lett. 2019, 19, 1227-1233. 29. Wang, B.; Yuan, S.; Li, Y.; Shi, L.; Wang, J., A New Dirac Cone Material: A Graphene-Like Be3C2 Monolayer. Nanoscale 2017, 9, 5577-5582. 30. Wang, H.; Lu, Z.; Yang, Z.; Li, X., A Novel Re-entrant Auxetic Honeycomb with Enhanced In-Plane Impact Resistance. Compos. Struct. 2019, 208, 758-770. 31. Li, D.; Yin, J.; Dong, L.; Lakes, R. S., Strong Re-entrant Cellular Structures with Negative Poisson’s Ratio. J. Mater. Sci. 2017, 53, 3493-3499. 32. Hou, X.; Silberschmidt, V. V., Metamaterials with Negative Poisson’s Ratio: A Review of Mechanical Properties and Deformation Mechanisms. 2015, 7, 155-179. 33. Prall, D.; Lakes, R. S., Properties of a Chiral Honeycomb with a Poisson's Ratio of -1. Int. J. Mech. Sci. 1997, 39, 305-314. 34. Jin, X.; Wang, Z.; Ning, J.; Xiao, G.; Liu, E.; Shu, X., Dynamic Response of Sandwich Structures with Graded Auxetic Honeycomb Cores under Blast Loading. Compos. Part B: Eng. 2016, 106, 206-217. 35. Qi, C.; Remennikov, A.; Pei, L. Z.; Yang, S.; Yu, Z. H.; Ngo, T. D., Impact and Close-In Blast Response of Auxetic Honeycomb-Cored Sandwich Panels: Experimental Tests and Numerical Simulations. Compos. Struct. 2017, 180, 161-178. 36. Peng, R.; Ma, Y.; Wu, Q.; Huang, B.; Dai, Y., Two-Dimensional Materials with Intrinsic Auxeticity: Progress and Perspectives. Nanoscale 2019, 11, 11413-11428. 37. Ma, F.; Jiao, Y.; Gao, G.; Gu, Y.; Bilic, A.; Chen, Z.; Du, A., Graphene-like Two-Dimensional Ionic Boron with Double Dirac Cones at Ambient Condition. Nano Lett. 2016, 16, 3022-3028. 38. Zhang, Z.; Yang, Y.; Penev, E. S.; Yakobson, B. I., Elasticity, Flexibility, and Ideal Strength of Borophenes. Adv. Funct. Mater. 2017, 27, 1605059. 15

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 18

39. Mannix, A. J.; Zhang, Z.; Guisinger, N. P.; Yakobson, B. I.; Hersam, M. C., Borophene as A Prototype for Synthetic 2D Materials Development. Nat. Nanotechnol. 2018, 13, 444-450. 40. Shi, L.; Ling, C.; Ouyang, Y.; Wang, J., High Intrinsic Catalytic Activity of Two-Dimensional Boron Monolayers for the Hydrogen Evolution Reaction. Nanoscale 2017, 9, 533-537. 41. Zhang, Z.; Penev, E. S.; Yakobson, B. I., Two-Dimensional Boron: Structures, Properties and Applications. Chem. Soc. Rev. 2017, 46, 6746. 42. Feng, B.; Zhang, J.; Zhong, Q.; Li, W.; Li, S.; Li, H.; Cheng, P.; Meng, S.; Chen, L.; Wu, K., Experimental Realization of Two-Dimensional Boron Sheets. Nat. Chem. 2016, 8, 563. 43. Smith, C. W.; Grima, J.; Evans, K., A Novel Mechanism for Generating Auxetic Behaviour in Reticulated Foams: Missing Rib Foam Model. Acta mater. 2000, 48, 4349-4356. 44. Saxena, K. K.; Das, R.; Calius, E. P., Three Decades of Auxetics Research-Materials with Negative Poisson's Ratio: a Review. Adv. Eng. Mater. 2016, 18, 1847-1870. 45. Ingrole, A.; Hao, A.; Liang, R., Design and Modeling of Auxetic and Hybrid Honeycomb Structures for In-Plane Property Enhancement. Mater. Design 2017, 117, 72-83. 46. Lee, C.; Wei, X.; Kysar, J. W.; Hone, J., Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene. Science 2008, 321, 385-388. 47. Andrew, R. C.; Mapasha, R. E.; Ukpong, A. M.; Chetty, N., Mechanical Properties of Graphene and Boronitrene. Phys. Rev. B 2012, 85, 125428. 48. Liu, K.; Yan, Q.; Chen, M.; Fan, W.; Sun, Y.; Suh, J.; Fu, D.; Lee, S.; Zhou, J.; Tongay, S.; Ji, J.; Neaton, J. B.; Wu, J., Elastic Properties of Chemical-Vapor-Deposited Monolayer MoS2, WS2, and Their Bilayer Heterostructures. Nano Lett. 2014, 14, 5097-5103. 49. Wei, Q.; Peng, X., Superior Mechanical Flexibility of Phosphorene and Few-Layer Black Phosphorus. Appl. Phys. Lett. 2014, 104, 251915. 50. Liu, F.; Ming, P.; Li, J., Ab Initio Calculation of Ideal Strength and Phonon Instability of Graphene under Tension. Phys. Rev. B 2007, 76, 064120. 51. Yan, J.-A.; Ruan, W. Y.; Chou, M. Y., Phonon Dispersions and Vibrational Properties of Monolayer, Bilayer, and Trilayer Graphene: Density-Functional Perturbation Theory. Phys. Rev. B 2008, 77, 125041. 52. Zhong, Q.; Kong, L.; Gou, J.; Li, W.; Sheng, S.; Yang, S.; Cheng, P.; Li, H.; Wu, K.; Chen, L., Synthesis of Borophene Nanoribbons on Ag(110) Surface. Phys. Rev. Mater. 2017, 1, 021001(R). 53. Wang, Y.; Lv, J.; Zhu, L.; Ma, Y., Crystal Structure Prediction via Particle-Swarm Optimization. Phys. Rev. B 2010, 82, 094116. 54. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 55. Kresse, G.; Furthmüller, J., Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169-11186. 56. Blöchl, P. E., Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979. 57. Togo, A.; Oba, F.; Tanaka, I., First-Principles Calculations of the Ferroelastic Transition between Rutile-Type and CaCl2-type SiO2 at High Pressures. Phys. Rev. B 2008, 78, 134106. 58. Henkelman, G.; Uberuaga, B. P.; Jónsson, H., A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 16

ACS Paragon Plus Environment

Page 17 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Materials & Interfaces

9901-9904. 59. Martyna, G. J.; Klein, M. L.; Tuckerman, M., Nosé–Hoover chains: The Canonical Ensemble via Continuous Dynamics. J. Chem. Phys. 1992, 97, 2635-2643. 60. Wang, H.; Li, X.; Li, P.; Yang, J., -Phosphorene: A Two Dimensional Material with a Highly Negative Poisson's Ratio. Nanoscale 2017, 9, 850-855. 61. Zhao, Y.; Li, X.; Liu, J.; Zhang, C.; Wang, Q., A New Anisotropic Dirac Cone Material: A B2S Honeycomb Monolayer. J. Phys. Chem. Lett. 2018, 9, 1815-1820.

TOC

17

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

551x340mm (144 x 144 DPI)

ACS Paragon Plus Environment

Page 18 of 18