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.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

Present state of global wetland extent and wetland methane modelling: conclusions from a model inter-comparison project (WETCHIMP)

Global wetlands are believed to be climate sensitive, and are the largest natural emitters of methane (CH4). Increased wetland CH4 emissions could act as a positive feedback to future warming. The Wetland and Wetland CH4 Inter-comparison of Models Project (WETCHIMP) investigated our present ability to simulate large-scale wetland characteristics and corresponding CH4 emissions. To ensure inter-comparability, we used a common experimental protocol driving all models with the same climate and carbon dioxide (CO2) forcing datasets. The WETCHIMP experiments were conducted for model equilibrium states as well as transient simulations covering the last century. Sensitivity experiments investigated model response to changes in selected forcing inputs (precipitation, temperature, and atmospheric CO2 concentration). Ten models participated, covering the spectrum from simple to relatively complex, including models tailored either for regional or global simulations. The models also varied in methods to calculate wetland size and location, with some models simulating wetland area prognostically, while other models relied on remotely sensed inundation datasets, or an approach intermediate between the two. Four major conclusions emerged from the project. First, the suite of models demonstrate extensive disagreement in their simulations of wetland areal extent and CH4 emissions, in both space and time. Simple metrics of wetland area, such as the latitudinal gradient, show large variability, principally between models that use inundation dataset information and those that independently determine wetland area. Agreement between the models improves for zonally summed CH4 emissions, but large variation between the models remains. For annual global CH4 emissions, the models vary by ±40% of the all-model mean (190 Tg CH4 yr−1). Second, all models show a strong positive response to increased atmospheric CO2 concentrations (857 ppm) in both CH4 emissions and wetland area. In response to increasing global temperatures (+3.4 °C globally spatially uniform), on average, the models decreased wetland area and CH4 fluxes, primarily in the tropics, but the magnitude and sign of the response varied greatly. Models were least sensitive to increased global precipitation (+3.9 % globally spatially uniform) with a consistent small positive response in CH4 fluxes and wetland area. Results from the 20th century transient simulation show that interactions between climate forcings could have strong non-linear effects. Third, we presently do not have sufficient wetland methane observation datasets adequate to evaluate model fluxes at a spatial scale comparable to model grid cells (commonly 0.5°). This limitation severely restricts our ability to model global wetland CH4 emissions with confidence. Our simulated wetland extents are also difficult to evaluate due to extensive disagreements between wetland mapping and remotely sensed inundation datasets. Fourth, the large range in predicted CH4 emission rates leads to the conclusion that there is both substantial parameter and structural uncertainty in large-scale CH4 emission models, even after uncertainties in wetland areas are accounted for.

Genetically determined heterogeneity of lung disease in a mouse model of airway mucus obstruction

Mucus clearance is an important airway innate defense mechanism. Airway-targeted overexpression of the epithelial Na(+) channel β-subunit [encoded by sodium channel nonvoltage gated 1, beta subunit (Scnn1b)] in mice [Scnn1b-transgenic (Tg) mice] increases transepithelial Na(+) absorption and dehydrates the airway surface, which produces key features of human obstructive lung diseases, including mucus obstruction, inflammation, and air-space enlargement. Because the first Scnn1b-Tg mice were generated on a mixed background, the impact of genetic background on disease phenotype in Scnn1b-Tg mice is unknown. To explore this issue, congenic Scnn1b-Tg mice strains were generated on C57BL/6N, C3H/HeN, BALB/cJ, and FVB/NJ backgrounds. All strains exhibited a two- to threefold increase in tracheal epithelial Na(+) absorption, and all developed airway mucus obstruction, inflammation, and air-space enlargement. However, there were striking differences in neonatal survival, ranging from 5 to 80% (FVB/NJspace enlargement and postnatal death were only present when Na(+) hyperabsorption occurred early, and 2) inflammation and mucus obstruction developed whenever Na(+) hyperabsorption was expressed. In summary, the genetic context and timing of airway innate immune dysfunction critically determines lung disease phenotype. These mouse strains may be useful to identify key modifier genes and pathways.

Dominance of chemokine ligand 2 and matrix metalloproteinase-2 and -9 and suppression of pro-inflammatory cytokines in the epidural compartment after intervertebral disc extrusion in a canine model

BACKGROUND CONTEXT In canine intervertebral disc (IVD) disease, a useful animal model, only little is known about the inflammatory response in the epidural space. PURPOSE To determine messenger RNA (mRNA) expressions of selected cytokines, chemokines, and matrix metalloproteinases (MMPs) qualitatively and semiquantitatively over the course of the disease and to correlate results to neurologic status and outcome. STUDY DESIGN/SETTING Prospective study using extruded IVD material of dogs with thoracolumbar IVD extrusion. PATIENT SAMPLE Seventy affected and 13 control (24 samples) dogs. OUTCOME MEASURES Duration of neurologic signs, pretreatment, neurologic grade, severity of pain, and outcome were recorded. After diagnostic imaging, decompressive surgery was performed. METHODS Messenger RNA expressions of interleukin (IL)-1β, IL-2, IL-4, IL-6, IL-8, IL-10, tumor necrosis factor (TNF), interferon (IFN)γ, MMP-2, MMP-9, chemokine ligand (CCL)2, CCL3, and three housekeeping genes was determined in the collected epidural material by Panomics 2.0 QuantiGene Plex technology. Relative mRNA expression and fold changes were calculated. Relative mRNA expression was correlated statistically to clinical parameters. RESULTS Fold changes of TNF, IL-1β, IL-2, IL-4, IL-6, IL-10, IFNγ, and CCL3 were clearly downregulated in all stages of the disease. MMP-9 was downregulated in the acute stage and upregulated in the subacute and chronic phase. Interleukin-8 was upregulated in acute cases. MMP-2 showed mild and CCL2 strong upregulation over the whole course of the disease. In dogs with severe pain, CCL3 and IFNγ were significantly higher compared with dogs without pain (p=.017/.020). Dogs pretreated with nonsteroidal anti-inflammatory drugs revealed significantly lower mRNA expression of IL-8 (p=.017). CONCLUSIONS The high CCL2 levels and upregulated MMPs combined with downregulated T-cell cytokines and suppressed pro-inflammatory genes in extruded canine disc material indicate that the epidural reaction is dominated by infiltrating monocytes differentiating into macrophages with tissue remodeling functions. These results will help to understand the pathogenic processes representing the basis for novel therapeutic approaches. The canine IVD disease model will be rewarding in this process.

The rabbit blood shunt subarachnoid haemorrhage model.

The recently introduced rabbit blood shunt subarachnoid haemorrhage model is based on the two standard procedures of subclavian artery cannulation and transcutaneous cisterna magna puncture. An extracorporeal shunt placed in between the arterial system and the subarachnoid space allows examiner-independent SAH in a closed cranium. Despite its straightforwardness, it is worth examining some specific features and characteristics of the model. We outline technical considerations to successfully perform the model with minimal mortality and morbidity. In addition, we discuss outcome measures, advantages and limitations, and the applicability of the model for the study of early brain injury and delayed cerebral vasospasm after SAH.

Articulated Statistical Shape Model-Based 2D-3D Reconstruction of a Hip Joint

In this paper, reconstruction of three-dimensional (3D) patient-specific models of a hip joint from two-dimensional (2D) calibrated X-ray images is addressed. Existing 2D-3D reconstruction techniques usually reconstruct a patient-specific model of a single anatomical structure without considering the relationship to its neighboring structures. Thus, when those techniques would be applied to reconstruction of patient-specific models of a hip joint, the reconstructed models may penetrate each other due to narrowness of the hip joint space and hence do not represent a true hip joint of the patient. To address this problem we propose a novel 2D-3D reconstruction framework using an articulated statistical shape model (aSSM). Different from previous work on constructing an aSSM, where the joint posture is modeled as articulation in a training set via statistical analysis, here it is modeled as a parametrized rotation of the femur around the joint center. The exact rotation of the hip joint as well as the patient-specific models of the joint structures, i.e., the proximal femur and the pelvis, are then estimated by optimally fitting the aSSM to a limited number of calibrated X-ray images. Taking models segmented from CT data as the ground truth, we conducted validation experiments on both plastic and cadaveric bones. Qualitatively, the experimental results demonstrated that the proposed 2D-3D reconstruction framework preserved the hip joint structure and no model penetration was found. Quantitatively, average reconstruction errors of 1.9 mm and 1.1 mm were found for the pelvis and the proximal femur, respectively.

Search for neutral Higgs bosons of the minimal supersymmetric standard model in pp collisions at sqrts = 8 TeV with the ATLAS detector

A search for the neutral Higgs bosons predicted by the Minimal Supersymmetric Standard Model (MSSM) is reported. The analysis is performed on data from proton-proton collisions at a centre-of-mass energy of 8 TeV collected with the ATLAS detector at the Large Hadron Collider. The samples used for this search were collected in 2012 and correspond to integrated luminosities in the range 19.5-20.3 fb−1. The MSSM Higgs bosons are searched for in the τ τ final state. No significant excess over the expected background is observed, and exclusion limits are derived for the production cross section times branching fraction of a scalar particle as a function of its mass. The results are also interpreted in the MSSM parameter space for various benchmark scenarios.

Measurement of the electroweak production of dijets in association with a Z-boson and distributions sensitive to vector boson fusion in proton-proton collisions at √s = 8 TeV using the ATLAS detector

Measurements of fiducial cross sections for the electroweak production of two jets in association with a Z-boson are presented. The measurements are performed using 20.3 fb−1 of proton-proton collision data collected at a centre-of-mass energy of p s = 8TeV by the ATLAS experiment at the Large Hadron Collider. The electroweak component is extracted by a fit to the dijet invariant mass distribution in a fiducial region chosen to enhance the electroweak contribution over the dominant background in which the jets are produced via the strong interaction. The electroweak cross sections measured in two fiducial regions are in good agreement with the Standard Model expectations and the background-only hypothesis is rejected with significance above the 5ơ level. The electroweak process includes the vector boson fusion production of a Z-boson and the data are used to place limits on anomalous triple gauge boson couplings. In addition, measurements of cross sections and differential distributions for inclusive Z-boson-plus-dijet production are performed in five fiducial regions, each with different sensitivity to the electroweak contribution. The results are corrected for detector effects and compared to predictions from the Sherpa and Powheg event generators.

CODE’s new solar radiation pressure model for GNSS orbit determination

The Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), which was developed in the early 1990s, is widely used in the International GNSS Service (IGS) community. For a rather long time, spurious spectral lines are known to exist in geophysical parameters, in particular in the Earth Rotation Parameters (ERPs) and in the estimated geocenter coordinates, which could recently be attributed to the ECOM. These effects grew creepingly with the increasing influence of the GLONASS system in recent years in the CODE analysis, which is based on a rigorous combination of GPS and GLONASS since May 2003. In a first step we show that the problems associated with the ECOM are to the largest extent caused by the GLONASS, which was reaching full deployment by the end of 2011. GPS-only, GLONASS-only, and combined GPS/GLONASS solutions using the observations in the years 2009–2011 of a global network of 92 combined GPS/GLONASS receivers were analyzed for this purpose. In a second step we review direct solar radiation pressure (SRP) models for GNSS satellites. We demonstrate that only even-order short-period harmonic perturbations acting along the direction Sun-satellite occur for GPS and GLONASS satellites, and only odd-order perturbations acting along the direction perpendicular to both, the vector Sun-satellite and the spacecraft’s solar panel axis. Based on this insight we assess in the third step the performance of four candidate orbit models for the future ECOM. The geocenter coordinates, the ERP differences w. r. t. the IERS 08 C04 series of ERPs, the misclosures for the midnight epochs of the daily orbital arcs, and scale parameters of Helmert transformations for station coordinates serve as quality criteria. The old and updated ECOM are validated in addition with satellite laser ranging (SLR) observations and by comparing the orbits to those of the IGS and other analysis centers. Based on all tests, we present a new extended ECOM which substantially reduces the spurious signals in the geocenter coordinate z (by about a factor of 2–6), reduces the orbit misclosures at the day boundaries by about 10 %, slightly improves the consistency of the estimated ERPs with those of the IERS 08 C04 Earth rotation series, and substantially reduces the systematics in the SLR validation of the GNSS orbits.

Efficiency of frost-cracking processes through space and time: An example from the eastern Italian Alps

It is widely accepted that climate has a strong impact and exerts important feedbacks on erosional processes and sediment transport mechanisms. However, the extent at which climate influences erosion is still a matter of debate. In this paper we test whether frost-cracking processes and related temperature variations can influence the sediment production and surface erosion in a small catchment situated in the eastern Italian Alps. To this extent, we first present a geomorphic map of the region that we complement with published 10Be-based denudation rates. We then apply a preexisting heat-flow model in order to analyze the variations of the frost-cracking intensity (FCI) in the study area, which could have controlled the sediment production in the basin. Finally, we compare the model results with the pattern of denudation rates and Quaternary deposits in the geomorphic map. The model results, combined with field observations, mapping, and quantitative geomorphic analyses, reveal that frost-cracking processes have had a primary role in the production of sediment where the intensity of sediment supply has been dictated and limited by the combined effect of temperature variations and conditions of bedrock preservation. These results highlight the importance of a yet poorly understood process for the production of sediment in mountain areas.

Temperature Prediction Model for Bone Drilling Based on Density Distribution and In Vivo Experiments for Minimally Invasive Robotic Cochlear Implantation

Surgical robots have been proposed ex vivo to drill precise holes in the temporal bone for minimally invasive cochlear implantation. The main risk of the procedure is damage of the facial nerve due to mechanical interaction or due to temperature elevation during the drilling process. To evaluate the thermal risk of the drilling process, a simplified model is proposed which aims to enable an assessment of risk posed to the facial nerve for a given set of constant process parameters for different mastoid bone densities. The model uses the bone density distribution along the drilling trajectory in the mastoid bone to calculate a time dependent heat production function at the tip of the drill bit. Using a time dependent moving point source Green's function, the heat equation can be solved at a certain point in space so that the resulting temperatures can be calculated over time. The model was calibrated and initially verified with in vivo temperature data. The data was collected in minimally invasive robotic drilling of 12 holes in four different sheep. The sheep were anesthetized and the temperature elevations were measured with a thermocouple which was inserted in a previously drilled hole next to the planned drilling trajectory. Bone density distributions were extracted from pre-operative CT data by averaging Hounsfield values over the drill bit diameter. Post-operative [Formula: see text]CT data was used to verify the drilling accuracy of the trajectories. The comparison of measured and calculated temperatures shows a very good match for both heating and cooling phases. The average prediction error of the maximum temperature was less than 0.7 °C and the average root mean square error was approximately 0.5 °C. To analyze potential thermal damage, the model was used to calculate temperature profiles and cumulative equivalent minutes at 43 °C at a minimal distance to the facial nerve. For the selected drilling parameters, temperature elevation profiles and cumulative equivalent minutes suggest that thermal elevation of this minimally invasive cochlear implantation surgery may pose a risk to the facial nerve, especially in sclerotic or high density mastoid bones. Optimized drilling parameters need to be evaluated and the model could be used for future risk evaluation.

Space debris attitude simulation - ιOTA (In-Orbit Tumbling Analysis)

Today, there is little knowledge on the attitude state of decommissioned intact objects in Earth orbit. Observational means have advanced in the past years, but are still limited with respect to an accurate estimate of motion vector orientations and magnitude. Especially for the preparation of Active Debris Removal (ADR) missions as planned by ESA’s Clean Space initiative or contingency scenarios for ESA spacecraft like ENVISAT, such knowledge is needed. ESA's “Debris Attitude Motion Measurements and Modelling” project (ESA Contract No. 40000112447), led by the Astronomical Institute of the University of Bern (AIUB), addresses this problem. The goal of the project is to achieve a good understanding of the attitude evolution and the considerable internal and external effects which occur. To characterize the attitude state of selected targets in LEO and GTO, multiple observation methods are combined. Optical observations are carried out by AIUB, Satellite Laser Ranging (SLR) is performed by the Space Research Institute of the Austrian Academy of Sciences (IWF) and radar measurements and signal level determination are provided by the Fraunhofer Institute for High Frequency Physics and Radar Techniques (FHR). The In-Orbit Tumbling Analysis tool (ιOTA) is a prototype software, currently in development by Hyperschall Technologie Göttingen GmbH (HTG) within the framework of the project. ιOTA will be a highly modular software tool to perform short-(days), medium-(months) and long-term (years) propagation of the orbit and attitude motion (six degrees-of-freedom) of spacecraft in Earth orbit. The simulation takes into account all relevant acting forces and torques, including aerodynamic drag, solar radiation pressure, gravitational influences of Earth, Sun and Moon, eddy current damping, impulse and momentum transfer from space debris or micro meteoroid impact, as well as the optional definition of particular spacecraft specific influences like tank sloshing, reaction wheel behaviour, magnetic torquer activity and thruster firing. The purpose of ιOTA is to provide high accuracy short-term simulations to support observers and potential ADR missions, as well as medium-and long-term simulations to study the significance of the particular internal and external influences on the attitude, especially damping factors and momentum transfer. The simulation will also enable the investigation of the altitude dependency of the particular external influences. ιOTA's post-processing modules will generate synthetic measurements for observers and for software validation. The validation of the software will be done by cross-calibration with observations and measurements acquired by the project partners.

Comet 1P/Halley multifluid MHD model for the Giotto fly-by

The interaction of comets with the solar wind has been the focus of many studies including numerical modeling. We compare the results of our multifluid MHD simulation of comet 1P/Halley to data obtained during the flyby of the European Space Agency's Giotto spacecraft in 1986. The model solves the full set of MHD equations for the individual fluids representing the solar wind protons, the cometary light and heavy ions, and the electrons. The mass loading, charge-exchange, dissociative ion-electron recombination, and collisional interactions between the fluids are taken into account. The computational domain spans over several million kilometers, and the close vicinity of the comet is resolved to the details of the magnetic cavity. The model is validated by comparison to the corresponding Giotto observations obtained by the Ion Mass Spectrometer, the Neutral Mass Spectrometer, the Giotto magnetometer experiment, and the Johnstone Plasma Analyzer instrument. The model shows the formation of the bow shock, the ion pile-up, and the diamagnetic cavity and is able to reproduce the observed temperature differences between the pick-up ion populations and the solar wind protons. We give an overview of the global interaction of the comet with the solar wind and then show the effects of the Lorentz force interaction between the different plasma populations.

Granular layer in the periplasmic space of gram-positive bacteria and fine structures of Enterococcus gallinarum and Streptococcus gordonii septa revealed by cryo-electron microscopy of vitreous sections.

High-resolution structural information on optimally preserved bacterial cells can be obtained with cryo-electron microscopy of vitreous sections. With the help of this technique, the existence of a periplasmic space between the plasma membrane and the thick peptidoglycan layer of the gram-positive bacteria Bacillus subtilis and Staphylococcus aureus was recently shown. This raises questions about the mode of polymerization of peptidoglycan. In the present study, we report the structure of the cell envelope of three gram-positive bacteria (B. subtilis, Streptococcus gordonii, and Enterococcus gallinarum). In the three cases, a previously undescribed granular layer adjacent to the plasma membrane is found in the periplasmic space. In order to better understand how nascent peptidoglycan is incorporated into the mature peptidoglycan, we investigated cellular regions known to represent the sites of cell wall production. Each of these sites possesses a specific structure. We propose a hypothetic model of peptidoglycan polymerization that accommodates these differences: peptidoglycan precursors could be exported from the cytoplasm to the periplasmic space, where they could diffuse until they would interact with the interface between the granular layer and the thick peptidoglycan layer. They could then polymerize with mature peptidoglycan. We report cytoplasmic structures at the E. gallinarum septum that could be interpreted as cytoskeletal elements driving cell division (FtsZ ring). Although immunoelectron microscopy and fluorescence microscopy studies have demonstrated the septal and cytoplasmic localization of FtsZ, direct visualization of in situ FtsZ filaments has not been obtained in any electron microscopy study of fixed and dehydrated bacteria.