|Home | About | Journals | Submit | Contact Us | Français|
Free-standing nanomembranes are two-dimensional materials with nanometer thickness but can have macroscopic lateral dimensions. We develop a fracture model to evaluate a pre-stretched free standing circular ultrathin nanomembrane and establish a relation between the energy release rate of a circumferential interface crack and the pre-strain in the membrane. Our results demonstrate that detachment cannot occur when the radius of the membrane is smaller than a critical size. This critical radius is inversely proportional to the Young’s modulus and square of the pre-strain of the membrane.
Free-standing ultrathin nanomembranes are a new class of two-dimensional materials that possess nanoscale thickness across macroscopic dimensions. Such nanomembranes are not only ultra-lightweight but also robust and flexible. It has been reported that elastic moduli of ultrathin nanomembranes can be 1–10 GPa with ultimate strengths of up to 100 MPa [1-6]. These striking properties of free-standing nanomembranes have resulted in a broad spectrum of applications in separation, sensing, biomedicine and energy harvesting [7-10].
These free-standing nanomembranes are usually subjected to residual tensions induced during formation or by thermal mismatch. As a result, they may detach from the suspending micro-holes and lose their functionality [11,12]. Experiments have demonstrated that smaller free standing membranes are structurally and thermally more stable than larger ones [9,10]. This size dependence is also observed by Cheng et al.  who fabricated free standing membranes composed of highly ordered arrays of gold nanoparticles linked by single strand DNA molecules. Briefly, these membranes are formed using the following process: first, droplets of a solution consisting of 13 nm diameter gold nanoparticles capped with 5′-thiolated single strand DNA (ssDNA) is dried on a thin silicon nitrate substrate with micron size holes in it. Satellite micro-droplets were formed and trapped in these microholes due to pinning of contact lines by the hole edges. As solvent evaporated, these pinned satellite micro-droplets thinned and formed monolayer nanoparticle membranes. Pre-tensions were built up on the membranes during this drying stage. It was found that these membranes tended to detach from the edges of larger microholes, but not smaller ones (Fig. (Fig.11).
We propose a model to understand the size effect on the detachment of membranes due to pre-stretch. Micro-holes used in the experiments are usually rectangular or circular. For simplicity, we assume that the micro-holes are circular with radius R. The thin substrate with micro-holes is assumed to be rigid. Let h R denote the thickness of the membrane. For a circular membrane without defect, the pre-stretch is equi-biaxial and is spatially uniform. Since the out of plane stresses are identically zero, the membrane deforms under plane stress conditions. With respect to a polar coordinate system (r, θ) with origin at the center of the circular membrane, the non-zero pre-strains in a membrane without defect are where are the axial and hoop strains, respectively. This equi-biaxial stretch state can be achieved mechanically by imposing a radial displacement of on the circumference of the membrane. The strain energy density of the pre-stretched membrane is where E and ν are the Young’s modulus and the Poisson’s ratio of the membrane, respectively. The total elastic strain energy Γ of the membrane without defect is
A simple scaling argument shows why membranes with larger radius are more likely to fail. Specifically, we envision two failure modes: adhesive or cohesive. In adhesive failure, a crack grows along the edge of the micro-hole. In cohesive failure, crack growth occurs inside a region of the membrane that is adjacent to the interface (e.g. see Fig. Fig.1a).1a). To quantify these failures, we define W as the energy required to detach or break a unit area of membrane in the interface region. This region includes the edge of the micro-hole as well as a thin layer of membrane adjacent to it. If adhesive failure occurs, W is identified as adhesion energy per unit area between the membrane and the substrate, while in cohesive failure, W is the facture toughness of the membrane. The energy needed to detach the entire membrane from the edge is 2πRhW. The elastic energy stored in the membrane Γ will not be large enough to detach the entire membrane if
Using (1), the inequality in (2) can be rewritten as
Equation (3) shows why sheets over smaller holes (smaller R) is less likely to fail. This conclusion is also valid for non-circular shaped micro-holes; in this case R should be replaced by the characteristic length of the micro-hole.
The above argument is too simplistic since it assumes that detachment occurs on the entire interface, whereas in reality, only partial detachment is observed (Fig. (Fig.1a).1a). We ask a more general question: suppose there is a defect on or near the interface between the membrane and the micro-hole edge, will this defect grow? If the defect grows, elastic energy is released locally to debond or break the membrane and part of the membrane is relaxed. To model this process, we consider an interface crack of length a < 2πR as shown in Fig. Fig.2.2. The surface of this crack is traction-free whereas the rest of the membrane circumference is subjected to a radial displacement of According to fracture mechanics , a necessary condition for the growth of this interface crack is
where G is the energy release rate of the crack. Energy release rate is the elastic energy per unit film thickness that would be released if the crack were to extend by a unit distance along the interface. In this analysis, deviation of the crack trajectory from the interface (e.g. in cohesive failure) is assumed to be small so that the energy release rate G can be computed by assuming the crack is right on the interface.
Linearity and dimensional argument show that
where f is an unknown dimensionless function of θo (0 ≤ θo ≤ π), the half angle sustained by the interface crack (see Fig. Fig.2).2). We will call the normalized energy release rate. Its behavior for small cracks (θo 1) can be found using a simple argument. In this limit, the crack can be viewed as a straight crack with length a lying on the interface between the lower half plane (circular membrane) and the upper half plane (rigid substrate). The crack is loaded by a hydrostatic tension at infinity. The energy release rate for this geometry is:
where the dimensionless numerical factor in (6) depends only on the membrane’s Poisson’s ratio, and is an order 1 quantity . Equation (6) shows that f(θoν) is a linear function of θo for small cracks. Applying the crack growth criterion (4) using (6), we found small defects (a/R 1) with
will not grow.
For longer cracks, growth is possible. To study this possibility, we note that (5) and (6) imply that the normalized energy release rate f(θo → 0, ν) → 0. In addition, since the strain energy of the membrane goes to zero as a → 2πR (or θo → π), the energy release rate is expected to vanish in this “long” crack limit. This explains why complete detachment does not occur. The fact that f vanishes at both end points (0 and π) means that f must have an interior maximum f* > 0 at some θo* (0, π). The existence of this maximum and the crack growth criterion (4) imply that no defect can grow if
where f*is a numerical constant which depends only on the Poisson’s ratio of the membrane.
The normalized energy release rate f and its maximum f* are determined numerically using finite element method. Details of the finite element calculation are given in Appendix A. Figure Figure33 plots the normalized energy release rate f(θo, ν) against θo for ν = 0, 0.1, 1/3, 0.4 and 0.49. The energy release rate increases linearly θo for small cracks as predicted by (7), reaches a maximum at θo* and then decreases rapidly to zero. The fact that the energy release rate decreases to zero after implies that the interface crack will always be arrested before it can de-cohere the entire membrane, consistent with experimental observation that the interface crack eventually arrests before the complete detachment of membrane can take place (Fig. (Fig.1).1). According to (8), irrespective of the size of a defect, it cannot grow if . Values of f* for different Poisson’s ratio ν can be approximated by f* = 4.26ν2 + 0.49ν + 0.92 with ν [0, 0.5]. This result suggests that membrane detachment can be prevented when the membrane radius R is smaller than a critical radius Rc:
Equation (9) also implies that stiff membranes are easier to detach than soft membranes for the same W and This result is consistent with our preliminary experimental results. In these experiments, the stiffness of the membrane can be tuned by varying the length of the DNA molecules attached to the gold particles . However, more data is needed to confirm this result.
Our experimental observations suggest that Rc is in the range of 1–4 μm. The Young’s modulus E of a typical membrane is reported to be about 6.5 GPa . Assuming ν is 1/3, the pre-strain is estimated to be 0.12%. This estimate is based on the initial slope of the force–displacement curve of indentation tests , where we have assumed that the force–displacement curve is controlled by the pre-tension for small deflections and approximated the indenter as a point load. There is no direct measurement of W. Using the values of RcE and v listed above, W is estimated to be 15–50 mJ/m2, which is consistent with the strength of van der Waals interaction .
In summary, a fracture mechanics model is used to explain why small free-standing membranes are more resistant to detachment. We show that detachment can be prevented by making the membrane smaller for a given pre-strain and W, which is consistent with our experimental observations. A useful expression for critical radius of the membrane is obtained and may guide future design of free-standing membrane systems.
Numerical calculations were carried out using ABAQUS, a finite element software . Radial displacement of was imposed on the edge of the circular membrane outside the crack while the crack face was left traction free. A typical mesh consisted of approximately 4,000 plane stress four-node quadrilateral elements. Finer mesh was used near the crack tip. The total strain energy of the system Î“ can be found for any normalized crack length Î¸o, and the energy release rate can be computed according to
where Î”Î¸o is a small increment in normalized crack length. We started the simulation with a very small crack and gradually increased its length. The calculation is tedious since each crack length increment requires remeshing and change of boundary conditions. This process was automated using a PYTHON script. The energy release rate was also computed using the J-integral option in ABAQUS. We found practically no difference between the two methods.
C.Y. Hui and R. Long are supported by a grant from the Department of Energy (DE-FG02-07ER46463). W. Cheng, M. Campolongo and D. Luo are partially supported by NYSTAR and the NSF CAREER award (grant number: 0547330).
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.