Isotropic Negative Thermal Expansion Metamaterials - ACS Publications

Jun 22, 2016 - thermal expansion are rare and usually unsatisfied in performance. Here, we propose a novel method to achieve two- and three-dimensiona...
1 downloads 0 Views 9MB Size
Download PDF
Research Article www.acsami.org

Isotropic Negative Thermal Expansion Metamaterials Lingling Wu,† Bo Li,‡ and Ji Zhou*,† †

State Key Laboratory of New Ceramics and Fine Processing, School of Materials Science and Engineering, Tsinghua University, Beijing 100084, China ‡ Advanced Materials Institute, Shenzhen Graduate School, Tsinghua University, Shenzhen, 518055, China S Supporting Information *

ABSTRACT: Negative thermal expansion materials are important and desirable in science and engineering applications. However, natural materials with isotropic negative thermal expansion are rare and usually unsatisfied in performance. Here, we propose a novel method to achieve two- and three-dimensional negative thermal expansion metamaterials via antichiral structures. The two-dimensional metamaterial is constructed with unit cells that combine bimaterial strips and antichiral structures, while the threedimensional metamaterial is fabricated by a multimaterial 3D printing process. Both experimental and simulation results display isotropic negative thermal expansion property of the samples. The effective coefficient of negative thermal expansion of the proposed models is demonstrated to be dependent on the difference between the thermal expansion coefficient of the component materials, as well as on the circular node radius and the ligament length in the antichiral structures. The measured value of the linear negative thermal expansion coefficient of the threedimensional sample is among the largest achieved in experiments to date. Our findings provide an easy and practical approach to obtaining materials with tunable negative thermal expansion on any scale. KEYWORDS: negative thermal expansion, antichiral, bimaterial, 3D printing, metamaterial

■

INTRODUCTION Negative thermal expansion materials are highly desired in infrastructure engineering and precision instrument manufacturing.1−6 An important related research is to design composite materials with a coefficient of thermal expansion (CTE) of zero by combining materials with negative- and positive-CTE in an appropriate ratio.7,8 Extensive use of zero-CTE materials could be found in various fields,9 such as optics10,11 and electronics.12 However, only a few materials naturally exhibit negative thermal expansion characteristics because the phenomenon is derived from abnormal mechanisms such as phase transitions,13 shortening of bond lengths,14 rigid unit modes,15 electronic effects,16 and magnetostriction.17 Although several natural negative thermal expansion materials such as perovskite, NaZn13-type La(Fe,Si,Co)13 compounds,18 Ag3[Co(CN)6],19 (Hf, Mg)(WO4)3,20 and the well-known ZrW2O8 family of materials4,21 exist, few are widely used because of their narrow temperature range of negative thermal expansion, low thermal expansion coefficient value, anisotropy of thermal response,22−24 and low design freedom. In the past decade, metamaterials have offered an entirely new route to construct artificial materials with abnormal properties not present in the component materials. Several approaches have been proposed to achieve negative thermal expansion structures.10,25−28 However, because most of the negative thermal expansion cells in these models are © XXXX American Chemical Society

asymmetric, i.e., their thermal expansion property is tailored in only one particular direction, few exhibit isotropic negative thermal expansion. In 2015, a two-dimensional lattice with a large negative CTE was fabricated from bimetallic strips by Ha based on the chiral structure.29 However, the relationship between the effective coefficient of thermal expansion of the model and its size parameters remains uninvestigated. Furthermore, three-dimensional negative thermal expansion metamaterials have not been constructed so far. In this paper, a novel model with an isotropic negative thermal expansion is proposed and demonstrated through the use of two- and three-dimensional antichiral structures. The metamaterials are also demonstrated to possess significant negative thermal expansion coefficients. Both simulation and experimental results show that tailored negative thermal expansion could be achieved by changing the component materials and size parameters of the model. Because the auxetic properties are independent of scale, we believe that the proposed structures could be applied at the macro-, micro-, or even nanoscale.30 Modeling Procedure and Simulation. The unit cell of the model is designed on the basis of antichiral structurea Received: May 12, 2016 Accepted: June 22, 2016

A

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

Figure 1. (A) Existing two- and three-dimensional systems of antichiral models. (B) Deformation schematic of the anti-tetrachiral model under thermal stress. (C) Deformation schematic of the anti-trichiral model under thermal stress.

Figure 2. (A) Diagram of the two-dimensional bimaterial anti-tetrachiral element and the anti-trichiral element. (B) Relationship between the circular node radius (r) and the effective CTE of antichiral elements with different lengths of ligament (l). (C) Relationship between the ligament length (l) and the effective CTE of antichiral elements with different node radii (r). (D) Simulated deformation of a two-dimensional anti-tetrachiral element with r = 40 mm, t = 1 mm, and l = 100 mm. (E) Simulated deformation of a two-dimensional anti-trichiral element with r = 40 mm, l = 100 mm, and t = 1 mm. (F) Relationship between the effective CTE of antichiral elements and the difference in CTE of their component materials. In the simulation of (F), l = 100 mm, r = 25 mm, t = 1 mm, α1 = 1.85 × 10−5 1/K, and α2 varies from 1.85 × 10−5 to 9.25 × 10−5 1/K in increments of 1.85 × 10−6 1/K. In the simulation of (B) and (C), r changes from 20 to 40 mm in increments of 5 mm and l varies from 100 to 140 mm in increments of 10 mm.

fact, the antichiral structures contain a series of similar systems when the number of attached ligaments (henceforth denoted as N) changes.36 The unit cell is rotationally symmetrical of order N, where N has been demonstrated to be restricted to only 3 (anti-trichiral) or 4 (anti-tetrachiral) in two-dimensional

well-established mechanical structure with an isotropic negative Poisson’s ratio.31−33 As proposed by Wojciechowski34 and initially realized by Lakes and Sigmound,31,35 the antichiral structure is constructed as a central point (henceforth referred to as a “node”) with several tangentially attached ligaments. In B

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

Figure 3. (A) Simulated deformation of a three-dimensional anti-tetrachiral element with r = 25 mm, l = 100 mm and t = 2 mm. (B) Relationship curves between the effective CTE along three orthogonal directions and the node radius for the three-dimensional anti-tetrachiral element.

systems.37 In a three-dimensional mechanism, only an antitetrachiral system can be constructed because a threedimensional anti-trichiral structure does not geometrically exist.38 The existing patterns of antichiral systems are shown in Figure 1A. To obtain a negative thermal expansion effect, we introduced bimaterial strips into the above antichiral structures. Previous research has demonstrated by finite elements that bimaterial elements in anti-tetrachiral honeycombs presented negative thermal properties.33 In this paper, the linear coefficient of thermal expansion αL was used to define the thermal expansion of the model, which is calculated by the following equation αL =

1 ∂L · L ∂T

varied from 20 to 40 mm in increments of 5 mm. The twodimensional anti-tetrachiral and anti-trichiral models are shown in Figure 2A. To study the relationship between the effective CTE of the elements and the CTE of the component materials, we set the CTE of material 1 (α1) to 1.85 × 10−5 1/K and increased the CTE of material 2 (α2) from 1.85 × 10−5 to 9.25 × 10−5 1/K in increments of 1.85 × 10−6 1/K. The simulated temperature range was 303.15−773.15 K. The simulated relationship curves between the effective CTEs of the antichiral models and their size parameters are plotted in Figure 2B,C, which show that, for both the antitrichiral and anti-tetrachiral structures, the effective CTE is proportional to the node radius r and ligament length l. The deformation of the two-dimensional antichiral structures with the temperature increasing are shown in Figure 2D,E. As evident in these figures, when the temperature was increased from 303.15 to 773.15 K, both structures shrank isotopically. The curve in Figure 2F indicates that the effective CTE of the bimaterial antichiral structures is proportional to the difference between the CTEs of the component materials. Furthermore, the relationship curves for the two-dimensional anti-trichiral and anti-tetrachiral structures are exactly the same when the circular node radius (r) and ligament length (l) are constant. The simulated results indicate that tunable negative thermal expansion can be achieved by carefully selecting the component materials and adjusting the geometrical parameters of the bimaterial antichiral models. The simulated deformation of a three-dimensional anti-tetrachiral element with r = 25 mm, l = 100 mm, and t = 2 mm is shown in Figure 3A, and the relationship curves between the effective CTE along three orthogonal directions and the node radius for the threedimensional anti-tetrachiral element are plotted in Figure 3B, which shows an isotropic negative thermal expansion of the model. The complete simulation results for the two-and threedimensional antichiral models are presented in the Supporting Information (Figure S1).

(1)

where L is the length of the unit cell along the studied direction, and T is the temperature. A bimaterial structure composed of two materials with different CTEs has been demonstrated to bend with a curvature of 1/ρ under a thermal stress proportional to the difference of the CTEs of the component materials39

1 ∝ (α2 − α1) ρ

(2)

where α1 and α2 are the CTEs of the component materials. Combining the bimaterial strips with the antichiral structures results in new negative thermal expansion metamaterials. Equation 2 is known to be valid for bimaterial strips made from thin ligaments where the end effects are neglected,40 and both thermal stress and applied loading stress can cause material deformation. Therefore, the bending of ligaments in the models due to thermal stress will result in a shortening of the node−node distance and a tendency to contract the structure in both in-plane principal directions, as shown in Figure 1B,C. As a result, the entire system will exhibit isotropic negative thermal expansion. To verify this hypothesis, theoretical simulation work was performed using the Solid Mechanics module in COMSOL Multiphysics commercial software. Models with two-dimensional anti-tetrachiral and anti-trichiral structures and a threedimensional anti-tetrachiral structure constructed with two materials with different CTEs were considered. To simplify the simulation work, the thickness of both component materials (t) was set to 1 mm, the ligament length (l) was varied from 100 to 140 mm in increments of 10 mm, and the node radius (r) was

■

METHODS Preparation of the Two-Dimensional Antichiral Samples. To experimentally verify the negative thermal expansion property of the models, we first fabricated twodimensional samples of bimetal antichiral metamaterials with a ligament length l of 100 mm and a variable node radius r. The model requires two component materials with different thermal expansion coefficients. We selected a bimetal plate with C

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

Figure 4. Two-dimensional measurement setup. (A) View of the measurement setup. (B) Fabricated anti-tetrachiral sample placed on the heating platform. (C) Fabricated parts for the two-dimensional antichiral samples.

Figure 5. Three-dimensional experiment. (A) Measured relationship between the CTE of the component materials of the three-dimensional antitetrachiral samples and the temperature. (B) Measured relationship between node radius and the effective CTE of the three-dimensional antitetrachiral samples. (C) Sample with r = 20 mm placed in a vacuum drying oven and the measured lengths of the sample. (D) Deformation of the tested sample with r = 20 mm after the increasing of temperature.

dimensional antichiral models with different node radii are shown in the Supporting Information (Figure S2). Preparation of the Three-Dimensional Anti-tetrachiral Samples. The three-dimensional anti-tetrachiral model cannot be fabricated by a traditional manufacturing process since it has a very complex bimaterial internal structure. Therefore, we fabricated the three-dimensional tetrachiral samples with l = 15 mm and t = 1.5 mm by a multi-3D printing process (Objet350

aluminum and copper as the raw materials. Aluminum has a relatively high thermal expansion coefficient of 2.3 × 10−5 1/K, whereas copper has a lower coefficient of 1.85 × 10−5 1/K. A cutting machine was used to cut the bimetal plates into strips with a thickness of 1 mm. The circular node was designed as a ring with grooved ears, and both the grooved ears and bimetal strips were drilled so that they could be firmly tightened with a screw, as shown in Figure 4. The fabricated samples of the twoD

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

Figure 6. (A, D) Experimentally measured results of an anti-tetrachiral sample with r = 40 mm, t = 1 mm, and l = 100 mm at 303.15 and 773.15 K, respectively. (C) Comparison of the simulated and experimentally measured relationship between the node radius (r) and the effective CTEs along two orthogonal directions for the two-dimensional anti-tetrachiral samples. (B, E) Experimentally measured results of an anti-trichiral sample with r = 40 mm, t = 1 mm, and l = 100 mm at 303.15 and 773.15 K, respectively. (F) Comparison of simulated and experimentally measured relationship between the node radius (r) and the effective CTE for the two-dimensional anti-trichiral samples. For both the simulation and experiments, α1 = 1.85 × 10−5 1/K and α2 = 2.3 × 10−5 1/K. The measurement unit in (A, B) and (D, E) is millimeters. U.S.) was used to determine the temperature on the surface of the sample in real time for all the samples. Since the thermal deformation temperature for the component materials of the three-dimensional samples is lower than 373.15 K, we carried out the three-dimensional experiment in a vacuum drying oven (DZF-6050, Shanghai Jinghong technology Corp., China.). The temperature in the oven was uniformly increased from 293.15 to 343.15 K (the CTE of the rubber-like material is higher than the rigid opaque material in this temperature range, as shown in Figure 5A) with a heating rate of 4.5 K/min. The computer-controlled CCD camera was fixed on the viewing window of the oven to take digital photographs from the same height and angle for all measurements. To accurately measure the deformation of the samples in the pictures recorded by the CCD camera, we used commercial software (Photoshop) to determine the length of the samples at different temperatures. Rulers placed in the three- and two-dimensional measurement setup were used as a reference for measurement. It must be mentioned that, to avoid the measurement error produced by the expansion of the ruler, all of the measurement for each sample in Photoshop was carried out by referring to the ruler in the picture recorded at room temperature (293.15 K). To improve the measurement accuracy, we made two measurements along each orthogonal direction for the two-dimensional anti-tetrachiral samples and calculated the average values, as shown in Figure 6A,D. The same process was used for the two-dimensional anti-trichiral samples, except that we measured three lengths along three directions with equal included angles because the anti-trichiral model exhibits 3-fold rotational symmetry, as shown in Figure 6B,E. The complete measurement photograph for the two-dimensional anti-tetrachiral and anti-trichiral samples are shown in Figures S4 and S5 of the

connex2, Stratasys Ltd.). The node radius r is varied from 15 to 25 mm with an increment of 2.5 mm. Two materials with different CTEs were chosen as lower-CTE (rigid opaque material: VeroWhitePlus RGD835) and higher-CTE (rubberlike material: TangoPlus FLX930) materials to fabricate the three-dimensional models. The fabricated model with a node radius of 20 mm is shown in Figure 5. The complete photograph of the three-dimensional samples with different node radii is shown in the Supporting Information (Figure S3). The relationship between the CTE of the selected materials and the temperature was measured by a dilatometer (Netzsch Dilatometer DIL402PC), and the results are plotted in Figure 5A. From the figure, we can see that the rubber-like material (TangoPlus FLX930) has a higher CTE than the rigid opaque material (VeroWhitePlus RGD835) in the testing temperature range.

■

EXPERIMENTAL SECTION

A measurement setup was constructed to verify the thermal expansion properties of the proposed two-dimensional antichiral models, as shown in Figure 4A. The experimental setup was composed of three parts: the heating platform and its controller, the CCD camera, and a computer. The sample was placed on the heating platform and covered with a quartz jar to prevent air convection. The computer-controlled CCD camera was positioned above the heating platform to observe the deformation of the sample in real time. To ensure soft lighting during the experiment, two light-compensating lamps were placed above the measurement equipment. An infrared camera (TiS65, Fluke Corp. E

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

ACS Applied Materials & Interfaces

■

Supporting Information, respectively. For the three-dimensional tetrachiral sample, we made four measurements in each photograph recorded from the front view of the sample, as shown in Figure 5C, and then calculated the average values.

Research Article

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.6b05717. Figure S1 shows the complete simulation results of the two- and three-dimensional antichiral structures. Figure S2 represents the photographs of the fabricated twodimensional antichiral samples with different node radii. Figure S3 presents the photograph of the fabricated three-dimensional anti-tetrachiral samples with different node radii. Figures S4 and S5 present the complete measured results of the two-dimensional anti-tetrachiral model and anti-trichiral model with different node radii, respectively (PDF)

■

RESULTS AND DISCUSSION The experimentally measured relationship between the effective CTE and the node radius (r) for the two-dimensional antitetrachiral and anti-trichiral samples is plotted in Figure 6C,F to provide a better comparison with the simulation results. As evident in these figures, both the simulation and experimental results show that the structures exhibit isotropic negative thermal expansion and the values of effective CTEs are proportional to the node radius (r). Reasonable agreement was achieved between the simulation and experimental results, with the experimental data generally showing lower values of the thermal expansion coefficient compared to the simulation results. This discrepancy mainly resulted from the friction between the heating platform and the sample. In Figure 6C, the measured value of the effective CTE for the anti-tetrachiral system reaches −68.1 × 10−6 1/K over a large operating temperature range from 303.15 to 773.15 K. The measurement result for the three-dimensional antitetrachiral sample is shown in Figure 5B, which indicates the relationship between the node radius (r) and the effective CTE of the sample. From the figure, we can see that the samples present a negative thermal expansion property from 293.15 to 343.15 K, and the magnitude of the effective CTE shows a proportional trend with the increase of the node radius of the sample. The discrepancies of the experiment were mainly caused by the gravity of the nodes. The three-dimensional experiment was not carried out as accurately as the twodimensional experiment due to the limitation of the measurement instrument and the precision of the multimaterial 3D printing fabrication process. However, it should be noted that, although the experiment result is imperfect, it provides a new way to realize the artificial three-dimensional NTE material and experimentally demonstrates its feasibility. We believe that, in the short future, more 3D bimaterial anti-tetrachiral models with a finer internal structure could be constructed with the development of precision manufacturing. This method proved to be a cheap and feasible way to promote the practical applications of the artificial negative thermal expansion materials. Efforts will be made on further improving the performance of the three-dimensional NTE material.

■

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Author Contributions

J.Z. conceived the idea. L.W. designed experiments and carried out numerical calculations. B.L. revised the full paper. All authors contributed to scientific discussion and critical revision of the article. Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China under Grant Nos.11274198 and 51532004 and the Science and Technology Plan of Shenzhen City under grant JCYJ 20150827165038323.

■

REFERENCES

(1) Lakes, R. Cellular Solids with Tunable Positive or Negative Thermal Expansion of Unbounded Magnitude. Appl. Phys. Lett. 2007, 90, 221905. (2) Ito, T.; Suganuma, T.; Wakashima, K. Glass fiber/polypropylene Composite Laminates with Negative Coefficients of Thermal Expansion. J. Mater. Sci. Lett. 1999, 18, 1363−1365. (3) Miller, W.; Smith, C. W.; Mackenzie, D. S.; Evans, K. E. Negative Thermal Expansion: A Review. J. Mater. Sci. 2009, 44, 5441−5451. (4) Mary, T. A.; Evans, J. S. O.; Vogt, T.; Sleight, A. W. Negative Thermal Expansion from 0.3 to 1050 K in ZrW2O8. Science 1996, 272, 90−92. (5) Miller, W.; Mackenzie, D. S.; Smith, C. W.; Evans, K. E. A Generalised Scale-independent Mechanism for Tailoring of Thermal Expansivity: Positive and Negative. Mech. Mater. 2008, 40, 351−361. (6) Sleight, A. W. Isotropic Negative Thermal Expansion. Annu. Rev. Mater. Sci. 1998, 28, 29−43. (7) Steeves, C. A.; Lucato, S. L. d. S. e; He, M.; Antinucci, E.; Hutchinson, J. W.; Evans, A. G. Concepts for Structurally Robust Materials that Combine Low Thermal Expansion with High Stiffness. J. Mech. Phys. Solids 2007, 55, 1803−1822. (8) Romao, C. P.; Perras, F. A.; Werner-Zwanziger, U.; Lussier, J. A.; Miller, K. J.; Calahoo, C. M.; Zwanziger, J. W.; Bieringer, M.; Marinkovic, B. A.; Bryce, D. L.; White, M. A. Zero Thermal Expansion in ZrMgMo3O12: NMR Crystallography Reveals Origins of Thermoelastic Properties. Chem. Mater. 2015, 27, 2633−2646. (9) Sleight, A. W. Negative Thermal Expansion Materials. Curr. Opin. Solid State Mater. Sci. 1998, 3, 128−131. (10) Kelly, A.; McCartney, L. N.; Clegg, W. J.; Stearn, R. J. Controlling Thermal Expansion to Obtain Negative Expansivity Using Laminated Composites. Compos. Sci. Technol. 2005, 65, 47−59.

■

CONCLUSIONS In conclusion, we proposed two- and three-dimensional artificial materials with substantial negative thermal expansion property and demonstrated their fabrication by bimaterial strips and a multimaterial 3D printing process. For both the twodimensional and three-dimensional samples, the relationship between the negative thermal expansion of the metamaterials and the node radius of the models was investigated. Large and tailorable negative CTEs were obtained with this approach. The mechanism proposed here is theoretically scale independent, meaning that the concept of manipulation of thermal expansion could be extended to the micro- or nanoscale if appropriate component materials and fabrication processes are available. These models have extensive potential applications in areas such as aerospace engineering, architectural engineering, and other domains which involve great temperature ranges that might cause thermal fatigue and cracking. F

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces (11) Lind, C. Two Decades of Negative Thermal Expansion Research: Where Do We Stand? Materials 2012, 5, 1125−1154. (12) Varghese, J.; Joseph, T.; Sebastian, M. T. ZrSiO4 Ceramics for Microwave Integrated Circuit Applications. Mater. Lett. 2011, 65, 1092−1094. (13) Panda, M. K.; Runcevski, T.; Sahoo, S. C.; Belik, A. A.; Nath, N. K.; Dinnebier, R. E.; Naumov, P. Colossal Positive and Negative Thermal Expansion and Thermosalient Effect in a Pentamorphic Organometallic Martensite. Nat. Commun. 2014, 5, 4811. (14) Welche, P. R. L.; Heine, V.; Dove, M. T. Negative Thermal Expansion in Beta-quartz. Phys. Chem. Miner. 1998, 26 (1), 63−77. (15) Pryde, A. K. A.; Hammonds, K. D.; Dove, M. T.; Heine, V.; Gale, J. D.; Warren, M. C. Origin of the Negative Thermal Expansion in ZrW2O8 and ZrV2O7. J. Phys.: Condens. Matter 1996, 8, 10973− 10982. (16) Arvanitidis, J.; Papagelis, K.; Margadonna, S.; Prassides, K.; Fitch, A. N. Temperature-induced Valence Transition and Associated Lattice Collapse in Samarium Fulleride. Nature 2003, 425, 599−602. (17) Hao, Y.; Gao, Y.; Wang, B.; Qu, J.; Li, Y.; Hu, J.; Deng, J. Negative Thermal Expansion and Magnetic Properties of Y2Al3Fe14−xMnx Compounds. Appl. Phys. Lett. 2001, 78, 3277−3279. (18) Huang, R.; Liu, Y.; Fan, W.; Tan, J.; Xiao, F.; Qian, L.; Li, L. Giant Negative Thermal Expansion in NaZn13-Type La (Fe, Si, Co) 13 Compounds. J. Am. Chem. Soc. 2013, 135, 11469−11472. (19) Goodwin, A. L.; Calleja, M.; Conterio, M. J.; Dove, M. T.; Evans, J. S. O.; Keen, D. A.; Peters, L.; Tucker, M. G. Colossal Positive and Negative Thermal Expansion in the Framework Material Ag3[Co(CN)6]. Science 2008, 319, 794−797. (20) Suzuki, T.; Omote, A. Negative Thermal Expansion in (Hf,Mg) (WO4)3. J. Am. Ceram. Soc. 2004, 87 (7), 1365−1367. (21) Evans, J. S. O.; David, W. I. F.; Sleight, A. W. Structural Investigation of the Negative Thermal Expansion Material ZrW2O8. Acta Crystallogr., Sect. B: Struct. Sci. 1999, 55, 333−340. (22) Das, D.; Jacobs, T.; Barbour, L. J. Exceptionally Large Positive and Negative Anisotropic Thermal Expansion of an Organic Crystalline Material. Nat. Mater. 2010, 9, 36−39. (23) Senn, M. S.; Bombardi, A.; Murray, C. A.; Vecchini, C.; Scherillo, A.; Luo, X.; Cheong, S. W. Negative Thermal Expansion in Hybrid Improper Ferroelectric Ruddlesden-Popper Perovskites by Symmetry Trapping. Phys. Rev. Lett. 2015, 114, 035701. (24) Fortes, A. D.; Suard, E.; Knight, K. S. Negative Linear Compressibility and Massive Anisotropic Thermal Expansion in Methanol Monohydrate. Science 2011, 331, 742. (25) Grima, J. N.; Oliveri, L.; Ellul, B.; Gatt, R.; Attard, D.; Cicala, G.; Recca, G. Adjustable and Negative Thermal Expansion from Multilayered Systems. Phys. Status Solidi RRL 2010, 4, 133−135. (26) Grima, J. N.; Ellul, B.; Gatt, R.; Attard, D. Negative Thermal Expansion from Disc, Cylindrical, and Needle Shaped Inclusions. Phys. Status Solidi B 2013, 250, 2051−2056. (27) Grima, J. N.; Ellul, B.; Attard, D.; Gatt, R.; Attard, M. Composites with Needle-like Inclusions Exhibiting Negative Thermal Expansion: A Preliminary Investigation. Compos. Sci. Technol. 2010, 70, 2248−2252. (28) Landert, M.; Kelly, A.; Stearn, J. R.; Hine, P. J. Negative Thermal Expansion of Laminates. J. Mater. Sci. 2004, 39, 3563−3567. (29) Ha, C. S.; Hestekin, E.; Li, J.; Plesha, M. E.; Lakes, R. S. Controllable Thermal Expansion of Large Magnitude in Chiral Negative Poisson’s Ratio Lattices. Phys. Status Solidi B 2015, 252, 1431−1434. (30) Gatt, R.; Attard, D.; Farrugia, P. S.; Azzopardi, K. M.; Mizzi, L.; Brincat, J. P.; Grima, J. N. A Realistic Generic Model for Antitetrachiral Systems. Phys. Status Solidi B 2013, 250, 2012−2019. (31) Lakes, R. Deformation Mechanism in Negative Poisson’s Ratio Materials: Structural Aspects. J. Mater. Sci. 1991, 26, 2287−2292. (32) Clausen, A.; Wang, F.; Jensen, J. S.; Sigmund, O.; Lewis, J. A. Topology Optimized Architectures with Programmable Poisson’s Ratio over Large Deformations. Adv. Mater. 2015, 27, 5523−5527. (33) Gatt, R.; Grima, J. N. Negative Compressibility. Phys. Status Solidi RRL 2008, 2, 236−238.

(34) Wojciechowski, K. W. Constant Thermodynamic Tension Monte-Carlo Studies of Elastic Properties of a Two-dimensional System of Hard Cyclic Hexamers. Mol. Phys. 1987, 61, 1247−1258. (35) Sigmund, O.; Torquato, S.; Aksay, I. A. On the Design of 1−3 Piezocomposites Using Topology Optimization. J. Mater. Res. 1998, 13, 1038−1048. (36) Alderson, A.; Alderson, K. L.; Attard, D.; Evans, K. E.; Gatt, R.; Grima, J. N.; Miller, W.; Ravirala, N.; Smith, C. W.; Zied, K. Elastic Constants of 3-, 4- and 6-connected Chiral and Anti-chiral Honeycombs Subject to Uniaxial In-plane Loading. Compos. Sci. Technol. 2010, 70, 1042−1048. (37) Grima, J. N. Ph.D. Thesis, University of Exeter, Exeter, U.K., 2000. (38) Coxeter, H. S. M. Regular Skew Polyhedra in Three and Four Dimensions, and Their Topological Analogues. Proceedings of the London Mathematical Society 1937, s2-43, 33−62. (39) Timoshenko, S. Analysis of Bi-metal Thermostats. J. Opt. Soc. Am. 1925, 11, 233−255. (40) Suhir, E. Stresses in Bi-metal Thermostats. J. Appl. Mech. 1986, 53, 657−660.

G

DOI: 10.1021/acsami.6b05717 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX