Temperature- and strain-rate-dependent microfabric evolution in monomineralic mylonite: evidence from in situ deformation of norcamphor

Norcamphor (C7H10O) was subjected to plane strain simple shear in a see-through deformation rig at four different strain rate and temperature conditions. Two transient stages in the microfabric evolution to steady state are distinguished. The grain scale mechanisms associated with the microstructural and textural evolution vary with the applied temperature, strain rate and strain. In high-temperature-low-strain-rate experiments, computer integrated polarization microscopy reveals that the texture evolution is closely related to the crystallographic rotation paths and rotation rates of individual grains. High c-axis rotation rates at low to intermediate shear strains are related to the development of a symmetrical c-axis cross girdle by the end of the first transient stage (γ = 1.5 to 2). During the second transient stage (γ = 1.5 to 6), the cross girdle yields to an oblique c-axis single girdle as c-axis rotation rates decrease and the relative activity of grain boundary migration recrystallization increases. Steady state (γ > 8) is characterized by a stable end orientation of the sample texture and the cyclic growth, rotation and consumption of individual grains within the aggregate.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Dynamic changes of cardiac conduction during rapid pacing

Slow conduction and unidirectional conduction block (UCB) are key mechanisms of reentry. Following abrupt changes in heart rate, dynamic changes of conduction velocity (CV) and structurally determined UCB may critically influence arrhythmogenesis. Using patterned cultures of neonatal rat ventricular myocytes grown on microelectrode arrays, we investigated the dynamics of CV in linear strands and the behavior of UCB in tissue expansions following an abrupt decrease in pacing cycle length (CL). Ionic mechanisms underlying rate-dependent conduction changes were investigated using the Pandit-Clark-Giles-Demir model. In linear strands, CV gradually decreased upon a reduction of CL from 500 ms to 230-300 ms. In contrast, at very short CLs (110-220 ms), CV first decreased before increasing again. The simulations suggested that the initial conduction slowing resulted from gradually increasing action potential duration (APD), decreasing diastolic intervals, and increasing postrepolarization refractoriness, which impaired Na(+) current (I(Na)) recovery. Only at very short CLs did APD subsequently shorten again due to increasing Na(+)/K(+) pump current secondary to intracellular Na(+) accumulation, which caused recovery of CV. Across tissue expansions, the degree of UCB gradually increased at CLs of 250-390 ms, whereas at CLs of 180-240 ms, it first increased and subsequently decreased. In the simulations, reduction of inward currents caused by increasing intracellular Na(+) and Ca(2+) concentrations contributed to UCB progression, which was reversed by increasing Na(+)/K(+) pump activity. In conclusion, CV and UCB follow intricate dynamics upon an abrupt decrease in CL that are determined by the interplay among I(Na) recovery, postrepolarization refractoriness, APD changes, ion accumulation, and Na(+)/K(+) pump function.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Volume dependence of high-frequency respiratory mechanics in healthy adults

Respiratory system input impedance (Zrs) at low to medium frequencies below 100 Hz, and study of its volume dependence, have been used extensively to quantify airway and tissue mechanics. Zrs at high oscillation frequencies including the first antiresonant frequency (far,1) may contain important information about airway mechanics. Changes in high-frequency Zrs with lung volume have not been studied. The volume-dependent behavior of high-frequency Zrs, specifically far,1 and respiratory system resistance at first antiresonance (Rrs(far,1)), was characterized in 16 healthy adults. Zrs was measured with a forced oscillation signal (5-302.5 Hz) through a wavetube setup. To track Zrs, subjects performed slow deep inspiratory and expiratory maneuvers over 30-s measurements, during which average impedance was calculated over 0.4-s intervals, with successive overlapping estimates every 0.156 s. Flow was measured using a pneumotachometer and integrated to obtain volume. Transpulmonary pressure dependence (Ptp) of Zrs was separately determined in five subjects. Both far,1 and Rrs(far,1) decreased with increasing lung volume and Ptp, consistent with an increase in airway caliber and decreased airway wall compliance as volume increased. These characterizations provide insight into airway mechanics, and are furthermore a necessary first step toward determining whether volume dependence of the first antiresonance is altered in disease.

Ultrastructural quantification of cell death after injurious compression of bovine calf articular cartilage

OBJECTIVE: It has been suggested that chondrocyte death by apoptosis may play a role in the pathogenesis of cartilage destruction in osteoarthritis, but the results of in-vivo and in-vitro investigations have been conflicting. To investigate further the cell death in our in-vitro model for traumatic joint injury, we performed a quantitative analysis by electron microscopy (EM) of cell morphology after injurious compression. For comparison, the TUNEL assay was also performed. DESIGN: Articular cartilage explant disks were harvested from newborn calf femoropatellar groove. The disks were subjected to injurious compression (50% strain at a strain rate of 100%/s), incubated for 3 days, and then fixed for quantitative morphological analysis. RESULTS: By TUNEL, the cell apoptosis rate increased from 7 +/- 2% in unloaded controls to 33 +/- 6% after injury (P=0.01; N=8 animals). By EM, the apoptosis rate increased from 5 +/- 1% in unloaded controls to 62 +/- 10% in injured cartilage (P=0.02, N=5 animals). Analysis by EM also identified that of the dead cells in injured disks, 97% were apoptotic by morphology. CONCLUSIONS: These results confirm a significant increase in cell death after injurious compression and suggest that most cell death observed here was by an apoptotic process.

Context-Oriented Programming

Context-dependent behavior is becoming increasingly important for a wide range of application domains, from pervasive computing to common business applications. Unfortunately, mainstream programming languages do not provide mechanisms that enable software entities to adapt their behavior dynamically to the current execution context. This leads developers to adopt convoluted designs to achieve the necessary runtime flexibility. We propose a new programming technique called Context-oriented Programming (COP) which addresses this problem. COP treats context explicitly, and provides mechanisms to dynamically adapt behavior in reaction to changes in context, even after system deployment at runtime. In this paper we lay the foundations of COP, show how dynamic layer activation enables multi-dimensional dispatch, illustrate the application of COP by examples in several language extensions, and demonstrate that COP is largely independent of other commitments to programming style.

Seismic anisotropy in the Morcles nappe shear zone: Implications for seismic imaging of crustal scale shear zones

Microstructures and textures of calcite mylonites from the Morcles nappe large-scale shear zone in southwestern Switzerland develop principally as a function of 1) extrinsic physical parameters including temperature, stress, strain, strain rate and 2) intrinsic parameters, such as mineral composition. We collected rock samples at a single location from this shear zone, on which laboratory ultrasonic velocities, texture and microstructures were investigated and quantified. The samples had different concentration of secondary mineral phases (< 5 up to 40 vol.%). Measured seismic P wave anisotropy ranges from 6.5% for polyphase mylonites (~ 40 vol.%) to 18.4% in mylonites with < 5 vol.% secondary phases. Texture strength of calcite is the main factor governing the seismic P wave anisotropy. Measured S wave splitting is generally highest in the foliation plane, but its origin is more difficult to explain solely by calcite texture. Additional texture measurements were made on calcite mylonites with low concentration of secondary phases (≤ 10 vol.%) along the metamorphic gradient of the shear zone (15 km distance). A systematic increase in texture strength is observed moving from the frontal part of the shear zone (anchimetamorphism; 280 °C) to the higher temperature, basal part (greenschist facies; 350–400 °C). Calculated P wave velocities become increasingly anisotropic towards the high-strain part of the nappe, from an average of 5.8% in the frontal part to 13.2% in the root of the basal part. Secondary phases raise an additional complexity, and may act either to increase or decrease seismic anisotropy of shear zone mylonites. In light of our findings we reinterpret the origin of some seismically reflective layers in the Grône–Zweisimmen line in southwestern Switzerland (PNR20 Swiss National Research Program). We hypothesize that reflections originate in part from the lateral variation in textural and microstructural arrangement of calcite mylonites in shear zones.

Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches

Since no single experimental or modeling technique provides data that allow a description of transport processes in clays and clay minerals at all relevant scales, several complementary approaches have to be combined to understand and explain the interplay between transport relevant phenomena. In this paper molecular dynamics simulations (MD) were used to investigate the mobility of water in the interlayer of montmorillonite (Mt), and to estimate the influence of mineral surfaces and interlayer ions on the water diffusion. Random Walk (RW) simulations based on a simplified representation of pore space in Mt were used to estimate and understand the effect of the arrangement of Mt particles on the meso- to macroscopic diffusivity of water. These theoretical calculations were complemented with quasielastic neutron scattering (QENS) measurements of aqueous diffusion in Mt with two pseudo-layers of water performed at four significantly different energy resolutions (i.e. observation times). The size of the interlayer and the size of Mt particles are two characteristic dimensions which determine the time dependent behavior of water diffusion in Mt. MD simulations show that at very short time scales water dynamics has the characteristic features of an oscillatory motion in the cage formed by neighbors in the first coordination shell. At longer time scales, the interaction of water with the surface determines the water dynamics, and the effect of confinement on the overall water mobility within the interlayer becomes evident. At time scales corresponding to an average water displacement equivalent to the average size of Mt particles, the effects of tortuosity are observed in the meso- to macroscopic pore scale simulations. Consistent with the picture obtained in the simulations, the QENS data can be described using a (local) 3D diffusion at short observation times, whereas at sufficiently long observation times a 2D diffusive motion is clearly observed. The effects of tortuosity measured in macroscopic tracer diffusion experiments are in qualitative agreement with RW simulations. By using experimental data to calibrate molecular and mesoscopic theoretical models, a consistent description of water mobility in clay minerals from the molecular to the macroscopic scale can be achieved. In turn, simulations help in choosing optimal conditions for the experimental measurements and the data interpretation. (C) 2014 Elsevier B.V. All rights reserved.

Evolution of mylonitic microfabrics (EMM), a computer application for educational purposes

`Evolution of mylonitic microfabrics' (EMM) is an interactive Filemaker Pro 3.0 application that documents a series of see-through deformation experiments on polycrystalline norcamphor. The application comprises computer animations, graphics and text explanations designed to give students and researchers insight into the interaction and dynamic nature of small-scale, mylonitic processes like intracrystalline glide, dynamic recrystallization and strain localization (microshearing). EMM shows how mylonitic steady state is achieved at different strain rates and temperatures. First, rotational mechanisms like glide-induced vorticity, subgrain rotation recrystallization and rigid-body rotation bring grains' crystal lattices into orientations that are favorable for intracrystalline glide. In a second stage, selective elimination of grains whose lattices are poorly oriented for glide involves grain boundary migration. This strengthens the texture. Temperature and strain rate affect both the relative activity of different strain accommodation mechanisms and the rate of microfabric change. Steady-state microfabrics are characterized by stable texture, grain size and shape-preferred orientations of grains and domains. This involves the cyclical generation and elimination of dynamically recrystallized grains and microshear zones.

pH-dependent distribution of chlorin e6 derivatives across phospholipid bilayers probed by NMR spectroscopy

The pH-dependent membrane adsorption and distribution of three chlorin derivatives, chlorin e6 (CE), rhodin G7 (RG), and monoaspartyl-chlorin e6 (MACE), in the physiological pH range (pH 6-8) were probed by NMR spectroscopy. Unilamellar vesicles consisting of dioleoyl-phosphatidyl-choline (DOPC) were used as membrane models. The chlorin derivatives were characterized with respect to their aggregation behavior, the pK(a) values of individual carboxylate groups, the extent of membrane adsorption, and their flip-flop rates across the bilayer membrane for pH 6-8. External membrane adsorption was found to be lower for RG than for CE and MACE. Both electrostatic interactions and the extent of aggregation seemed to be the main determinants of membrane adsorption. Rate constants for chlorin transfer across the membrane were found to correlate strongly with the pH of the surrounding medium, in particular, for CE and RG. In acidic solution, CE and RG transfer across the membrane was strongly accelerated, and in basic solution, all compounds were retained, mostly in the outer monolayer. In contrast, MACE flip-flop across the membrane remained very low even at pH 6. The protonation of ionizable groups is suggested to be a major determinant of chlorin transfer rates across the bilayer. pK(a) values of CE and RG were found to be between 6 and 8, and two of the carboxylate groups in MACE had pK(a) values below 6. For CE and RG, the kinetic profiles at acidic pH indicated that the initial fast membrane distribution was followed by secondary steps that are discussed in this article.

Pericellular fibronectin is required for RhoA-dependent responses to cyclic strain in fibroblasts

To test the hypothesis that the pericellular fibronectin matrix is involved in mechanotransduction, we compared the response of normal and fibronectin-deficient mouse fibroblasts to cyclic substrate strain. Normal fibroblasts seeded on vitronectin in fibronectin-depleted medium deposited their own fibronectin matrix. In cultures exposed to cyclic strain, RhoA was activated, actin-stress fibers became more prominent, MAL/MKL1 shuttled to the nucleus, and mRNA encoding tenascin-C was induced. By contrast, these RhoA-dependent responses to cyclic strain were suppressed in fibronectin knockdown or knockout fibroblasts grown under identical conditions. On vitronectin substrate, fibronectin-deficient cells lacked fibrillar adhesions containing alpha5 integrin. However, when fibronectin-deficient fibroblasts were plated on exogenous fibronectin, their defects in adhesions and mechanotransduction were restored. Studies with fragments indicated that both the RGD-synergy site and the adjacent heparin-binding region of fibronectin were required for full activity in mechanotransduction, but not its ability to self-assemble. In contrast to RhoA-mediated responses, activation of Erk1/2 and PKB/Akt by cyclic strain was not affected in fibronectin-deficient cells. Our results indicate that pericellular fibronectin secreted by normal fibroblasts is a necessary component of the strain-sensing machinery. Supporting this hypothesis, induction of cellular tenascin-C by cyclic strain was suppressed by addition of exogenous tenascin-C, which interferes with fibronectin-mediated cell spreading.

Role of RhoA/ROCK-dependent actin contractility in the induction of tenascin-C by cyclic tensile strain

In chick embryo fibroblasts, the mRNA for extracellular matrix protein tenascin-C is induced 2-fold by cyclic strain (10%, 0.3 Hz, 6 h). This response is attenuated by inhibiting Rho-dependent kinase (ROCK). The RhoA/ROCK signaling pathway is primarily involved in actin dynamics. Here, we demonstrate its crucial importance in regulating tenascin-C expression. Cyclic strain stimulated RhoA activation and induced fibroblast contraction. Chemical activators of RhoA synergistically enhanced the effects of cyclic strain on cell contractility. Interestingly, tenascin-C mRNA levels perfectly matched the extent of RhoA/ROCK-mediated actin contraction. First, RhoA activation by thrombin, lysophosphatidic acid, or colchicine induced tenascin-C mRNA to a similar extent as strain. Second, RhoA activating drugs in combination with cyclic strain caused a super-induction (4- to 5-fold) of tenascin-C mRNA, which was again suppressed by ROCK inhibition. Third, disruption of the actin cytoskeleton with latrunculin A abolished induction of tenascin-C mRNA by chemical RhoA activators in combination with cyclic strain. Lastly, we found that myosin II activity is required for tenascin-C induction by cyclic strain. We conclude that RhoA/ROCK-controlled actin contractility has a mechanosensory function in fibroblasts that correlates directly with tenascin-C gene expression. Previous RhoA/ROCK activation, either by chemical or mechanical signals, might render fibroblasts more sensitive to external tensile stress, e.g., during wound healing.

Cell shape-dependent early responses of fibroblasts to cyclic strain

Randomly spread fibroblasts on fibronectin-coated elastomeric membranes respond to cyclic strain by a varying degree of focal adhesion assembly and actin reorganization. We speculated that the individual shape of the cells, which is linked to cytoskeletal structure and pre-stress, might tune these integrin-dependent mechanotransduction events. To this aim, fibronectin circles, squares and rectangles of identical surface area (2000μm(2)) were micro-contact printed onto elastomeric substrates. Fibroblasts plated on these patterns occupied the corresponding shapes. Cyclic 10% equibiaxial strain was applied to patterned cells for 30min, and changes in cytoskeleton and cell-matrix adhesions were quantified after fluorescence staining. After strain, megakaryocytic leukemia-1 protein translocated to the nucleus in most cells, indicating efficient RhoA activation independently of cell shape. However, circular and square cells (with radial symmetry) showed a significantly greater increase in the number of actin stress fibers and vinculin-positive focal adhesions after cyclic strain than rectangular (bipolar) cells of identical size. Conversely, cyclic strain induced larger changes in pY397-FAK positive focal complexes and zyxin relocation from focal adhesions to stress fibers in bipolar compared to symmetric cells. Thus, radially symmetric cells responded to cyclic strain with a larger increase in assembly, whereas bipolar cells reacted with more pronounced reorganization of actin stress fibers and matrix contacts. We conclude that integrin-mediated responses to external mechanical strain are differentially modulated in cells that have the same spreading area but different geometries, and do not only depend on mere cell size.