Appendix (Supplementary Material)

Similar documents
Lecture Note October 1, 2009 Nanostructure characterization techniques

Mapping the mechanical stiffness of live cells with the scanning ion conductance microscope

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).

Measuring Young s modulus of 20LP10L20-LLA40 Microspheres and Gelatin-Methacrylamide (GelMA) Hydrogel using nanoindentation

Characterisation Programme Polymer Multi-scale Properties Industrial Advisory Group 22 nd April 2008

ADVANCED DYNAMIC MECHANICAL ANALYSIS OF A TIRE SAMPLE BY NANOINDENTATION

Improving the accuracy of Atomic Force Microscope based nanomechanical measurements. Bede Pittenger Bruker Nano Surfaces, Santa Barbara, CA, USA

MECHANICAL PROPERTIES OF HYDROGEL USING NANOINDENTATION

SIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING

Elasticity Models for the Spherical Indentation of Gels and Soft Biological Tissues

Revealing bending and force in a soft body through a plant root inspired. approach. Lucia Beccai 1* Piaggio 34, Pontedera (Italy)

Basic Laboratory. Materials Science and Engineering. Atomic Force Microscopy (AFM)

Final Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13

On Temporal Instability of Electrically Forced Axisymmetric Jets with Variable Applied Field and Nonzero Basic State Velocity

MSE 383, Unit 3-3. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.

Intensity (a.u.) Intensity (a.u.) Raman Shift (cm -1 ) Oxygen plasma. 6 cm. 9 cm. 1mm. Single-layer graphene sheet. 10mm. 14 cm

ME 323 Examination #2

Micro-rheology of cells and soft matter with the NanoTracker

Supplementary Information. Text S1:

Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses

Exact Shape Functions for Timoshenko Beam Element

Slanted Functional Gradient Micropillars for Optimal Bioinspired Dry

Analysis of contact deformation between a coated flat plate and a sphere and its practical application

Time-Dependent Nanomechanics of Cartilage

ERASMUS UNIVERSITY ROTTERDAM Information concerning the Entrance examination Mathematics level 2 for International Business Administration (IBA)

Lecture 6 mechanical system modeling equivalent mass gears

Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis

SUPPLEMENTARY INFORMATION 1

Math 216 Second Midterm 28 March, 2013

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

Determining the Elastic Modulus and Hardness of an Ultrathin Film on a Substrate Using Nanoindentation

Mechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection

64. Measurement of vibrations of rotating elements in a folding machine

Linear viscoelastic behavior

Introductory guide to measuring the mechanical properties of nanoobjects/particles

Dynamic Mechanical Analysis (DMA) of Polymers by Oscillatory Indentation

Supporting Information

Application Note #149 Improving the Accuracy of Nanomechanical Measurements with Force-Curve-Based AFM Techniques

Atomic Force Microscopy imaging and beyond

Exploring Nonlinear Oscillator Models for the Auditory Periphery

Lecture 7 Constitutive Behavior of Asphalt Concrete

Torsional Buckling of Double-Walled Carbon Nanotubes

SUPPLEMENTARY NOTES Supplementary Note 1: Fabrication of Scanning Thermal Microscopy Probes

Synthesis, Characterization and Tribological Evaluation of SDS stabilized. Antiwear Lubricant Additives

Unconfined Compression of Articular Cartilage: Nonlinear Behavior and Comparison With a Fibril-Reinforced Biphasic Model

Supplementary Figures

Frontiers of Fracture Mechanics. Adhesion and Interfacial Fracture Contact Damage

Strain Measurements. Isaac Choutapalli

P = ρ{ g a } + µ 2 V II. FLUID STATICS

height trace of a 2L BN mechanically exfoliated on SiO 2 /Si with pre-fabricated micro-wells. Scale bar 2 µm.

202 Index. failure, 26 field equation, 122 force, 1

Determining thermal noise limiting properties of thin films

Exact Free Vibration of Webs Moving Axially at High Speed

Appendix A Glossary of Mathematical Symbols

Lateral Nanomechanics of Cartilage Aggrecan Macromolecules

Measurement Techniques for Engineers. Motion and Vibration Measurement

9 MECHANICAL PROPERTIES OF SOLIDS

COMPRESSION INSTRUMENT FOR TISSUE EXPERIMENTS (CITE) AT THE MESO-SCALE: DEVICE VALIDATION

Micro-Rheology Measurements with the NanoTracker

Supplementary Information. Synthesis of soft colloids with well controlled softness

PERIYAR CENTENARY POLYTECHNIC COLLEGE PERIYAR NAGAR - VALLAM THANJAVUR. DEPARTMENT OF MECHANICAL ENGINEERING QUESTION BANK

I. Comparison between rheological data obtained by using the shear rate sweep and frequency

Regupol. Vibration

EE C247B / ME C218 INTRODUCTION TO MEMS DESIGN SPRING 2014 C. Nguyen PROBLEM SET #4

Physics 210: Worksheet 26 Name:

Lecture #2: Split Hopkinson Bar Systems

Regufoam. Vibration 270 Plus.

Size dependence of the mechanical properties of ZnO nanobelts

DYNAMIC MECHANICAL ANALYZER DMA 2980

1 Introduction IPICSE-2016

Nonlinear Dynamic Systems Homework 1

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads

An Experimental Investigation of the Elastic and Viscoelastic Properties of Small

STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS

Keysight Technologies Instrumented Indentation Testing with the Keysight Nano Indenter G200. Application Note

DETERMINATION OF THE ADHESION PROPERTIES OF MICA VIA ATOMIC FORCE SPECTROSCOPY

Characterization of Fixture-Workpiece Static Friction

Characterisation of the mechanical behaviour of a polyurethane elastomer based on indentation and tensile creep experiments

FRICTION FORCES ON THE LUBRICATED SURFACES IN MICRO AND NANO SCALE

AERO 214. Lab II. Measurement of elastic moduli using bending of beams and torsion of bars

AFM Imaging In Liquids. W. Travis Johnson PhD Agilent Technologies Nanomeasurements Division

589. Investigation of motion control of piezoelectric unimorph for laser shutter systems

EE 5344 Introduction to MEMS CHAPTER 6 Mechanical Sensors. 1. Position Displacement x, θ 2. Velocity, speed Kinematic

Name : Applied Physics II Exam One Winter Multiple Choice ( 7 Points ):

External Pressure... Thermal Expansion in un-restrained pipeline... The critical (buckling) pressure is calculated as follows:

BioAFM spectroscopy for mapping of Young s modulus of living cells

3.052 Nanomechanics of Materials and Biomaterials Thursday 02/15/07 Prof. C. Ortiz, MIT-DMSE I LECTURE 4: FORCE-DISTANCE CURVES

CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS

MSC Elastomers Seminar Some Things About Elastomers

Keysight Technologies Carbon Nanotube Tips for MAC Mode AFM Measurements in Liquids. Application Note

Long-base Hydrostatic Incline Meter. gravitational antenna

: APPLIED MECHANICS & STRENGTH OF MATERIALS COURSE CODE : 4021 COURSE CATEGORY : A PERIODS/ WEEK : 5 PERIODS/ SEMESTER : 75 CREDIT : 5 TIME SCHEDULE

EFFECTIVE SIMULATION APPROACH FOR STUDY OF CARBON NANOTUBE MECHANICAL PROPERTIES

Finite Element Analysis of Permeation Tests on Articular Cartilage under Different Testing Conditions Using COMSOL Multiphysics.

Presented By: EAS 6939 Aerospace Structural Composites

Features of static and dynamic friction profiles in one and two dimensions on polymer and atomically flat surfaces using atomic force microscopy

ERASMUS UNIVERSITY ROTTERDAM

8. More about calculus in physics

Figure 1.1: Flaccid (a) and swollen (b) red blood cells being drawn into a micropipette. The scale bars represent 5 µm. Figure adapted from [2].

Numerical analysis of three-lobe journal bearing with CFD and FSI

Transcription:

Appendix (Supplementary Material) Appendix A. Real and Imaginary Parts of the complex modulus E*(f), and E*(f) for the Pyramidal Proe Tip The magnitude of the complex modulus, E*(f), and the phase angle, (f), etween the measured force and the applied sinusoidal displacement (Fig. 1d) can e expressed in terms of (the real) storage modulus, E'(f), and the (imaginary) complex modulus, E''(f), as follows: E * (f ) E (f ) je (f ) (A1) * E (f ) E (f ) E (f ) (A) E(f ) tan (f ) (A3) E (f ) For the pyramidal proe tip, E*(f), was calculated using a Taylor series expansion of the Hertzian contact model (Radmacher 1997; Mahaffy 00): * 1 F(f ) E pyramidal (f ) (A) d0 tan d(f ) where : d 0 is indentation depth at zero frequency, Poisson's ratio is assumed to e 0. (Freeman 199), the pyramidal tip opening angle α is 35 (Ng 007), F(f) is the magnitude of the force and d(f) is the magnitude of the applied displacement at frequency f during dynamic oscillatory compression.

Appendix B. Indentation stiffness of single chondrocytes with GF treatment during time in culture For all of the experiments of Figs -, a force-displacement (F-d) curve was otained during the initial loading phase of each cell (Fig. 1c). This F-d curve was also used to characterize the quasistatic ehavior of each cell prior to force relaxation and susequent dynamic loading. The F-d curves cells in FBS and GF culture as a function of time in culture are shown in Fig. S1 for oth micron-sized and nano-sized proe tips. F-d curves varied with culture duration and tip geometry in oth medium conditions (two-way ANOVA, p<0.05) (Fig. S1). For each proe tip, the stiffness (slope of the F-d curve) depended significantly on the GF treatment and culture duration as well (two-way ANOVA, p<0.05). Chondrocytes cultured in FBS medium (Fig S1a,c) up to day 8 exhiited much lower stiffness than day 0 chondrocytes which had no cell-associated matrix, while chondrocytes treated with GF exhiited stiffness similar to or greater than day 0 cells after day 1 Fig S1,d). These trends are consistent with those reported in a previous study (Ng 007). Appendix C. Characteristic time constants for the poroelastic ehavior of tissue-engineered cell-associated matrix As descried in the main text, we estimated the characteristic poroelastic relaxation time, p ~ [L /(Hk)] to understand the mechanisms underlying the frequency dependence of the measured dynamic stiffness of the cell-associated matrix. The geometric configuration comprised of the cell-associated matrix

sandwiched in-etween the impermeale spherical proe tip and the cell memrane (Fig. 1a) is qualitatively approximated y a uniaxial unconfined compression geometry; the characteristic length scale, L, over which fluid flows is then approximated y the radius of the proe tip,.5 µm. The equilirium modulus H can e approximated y the quasi-static indentation stiffness (~0.1 kpa for day 1 and 0.17 kpa for day 8 in FBS culture (Ng 007). To estimate the hydraulic permeaility, k, it was hypothesized that GAG chains are the most critical determinant of the permeaility of the newly developing matrix, ased on previous measurements and models for the permeaility of native cartilage (Maroudas 1987). A unit cell model was used to calculate k (Happel 1959; Eisenerg 1988); the unit cell is comprised of a central charged GAG rod of radius, a, and its surrounding concentric cylindrical fluid shell having outer radius, (i.e., is the average center-to-center GAG spacing within the matrix (Eisenerg 1988)): k 1 3 1 8 ln 3 a a a 3 ln a 1 a a (D1) where is the fluid viscosity (10-3 Pa s). Outer radius was independently calculated to e ~0-70 nm using the Poisson-Boltzmann cell model (Buschmann 1995) along with the measured total GAG content of the cell-associated matrix (C GAG = 0~130 g/ml for days 7-8, (Ng 007)): 1 / CGAG BN M (D) CS

where B is the distance etween charge groups along the GAG chain (0.6nm), N is the Avogadro's numer (6.0 x 10 3 ), and M cs is the molecular weight of a dissociated chondroitin sulfate disaccharide (58g/mol) (Buschmann 1995). With the estimated value of, equation D1 yields values for the permeaility, k, of the newly developing matrix, which decreased from 6.6 x10-1 m /N s for day 7 to 1.7 x10-1 m /N s for day 8. This range of values for k is much higher than that of the mature chondron (10-17 m /N s (Alexopoulos 005)), suggesting that the newly developed cell-associated matrix is much less resistant to the fluid flow. Based on the aove estimates for L, H, and k, we calculated values for p to e in the range 0.011 0.01 sec for day 1-8 cells in FBS culture. This corresponds to characteristic frequencies f p ~[1/ p ] of 50-90 Hz Appendix D. Analysis and correction for the effects of hydrodynamic drag on amplitude and phase angle of measured force We performed the same dynamic oscillatory compression after the rampand-hold on a silicon sustrate as a control sample, which is elastic and approximated as a hard surface comparing to the cartilage and AFM cantilever (Han 008). Under lower frequencies (< 100 Hz), the phase difference etween the acquired z-piezo movement and the cantilever deflection are negligile; under higher frequencies ( 100 Hz), there is a significant phase difference etween the z-piezo and the cantilever deflection measured on the mica sample, possily due to the data acquisition process of a cantilever and the z-piezo movement and hydrodynamic drags at high frequencies. Hence, the phase angle

difference (e.g. for frequency f, f, control ) was calculated at each tested frequency and we sutracted the corresponding phase angle difference from the value otained on the silicon surface to correct the artificial phase differences from the AFM instrument and hydrodynamic effects: f,real (E1) f,measured f,control

Force (nn) Force (nn) Force (nn) Force (nn) Figure S1 (a).0 10% FBS(n=-5) ().0 IGF-1+OP-1 (n=5) (mean± SEM, n=,5) (n=5) 1.5 Day 0 1.5 Day 8 Day 0 Loading Day 8 Day 1 Day 1 Day 1 Day 1 (c).0 10% FBS (n=3-5) 1.5 Day 0 Day 8 Day 1 Day 1 (d).0 IGF-1+OP-1(n=5) (n=5) 1.5 Day 8 Day 1 Day 0 Day 1 Figure S1 Force-displacement curves of individual chondrocytes and chondrocyte-matrix composites using each proe tip: (a,c) FBS culture and (,d) GF culture (IGF-1+OP-1). n = numer of cells, mean ± SEM, n=3 for day 0 and n = -5 for day 7-8