Biomechanical model of human cornea based on stromal microstructure

The optical characteristics of the human cornea depends on the mechanical balance between the intra-ocular pressure and intrinsic tissue stiffness. A wide range of ophthalmic surgical procedures alter corneal biomechanics to induce local or global curvature changes for the correction of visual acuity. Due to the large number of surgical interventions performed every day, a deeper understanding of corneal biomechanics is needed to improve the safety of these procedures and medical devices. The aim of this study is to propose a biomechanical model of the human cornea, based on stromal microstructure. The constitutive mechanical law includes collagen fiber distribution based on X-ray scattering analysis, collagen cross-linking, and fiber uncrimping. Our results showed that the proposed model reproduced inflation and extensiometry experimental data [Elsheikh et al., Curr. Eye Res., 2007; Elsheikh et al., Exp. Eye Res., 2008] successfully. The mechanical properties obtained for different age groups demonstrated an increase in collagen cross-linking for older specimens. In future work such a model could be used to simulate non-symmetric interventions, and provide better surgical planning.

Numerical model of the human cornea based on stromal microstructure: identification of mechanical properties for two age groups

The optical quality of the human eye mainly depends on the refractive performance of the cornea. The shape of the cornea is a mechanical balance between intraocular pressure and tissue intrinsic stiffness. Several surgical procedures in ophthalmology alter the biomechanics of the cornea to provoke local or global curvature changes for vision correction. Legitimated by the large number of surgical interventions performed every day, the demand for a deeper understanding of corneal biomechanics is rising to improve the safety of procedures and medical devices. The aim of our work is to propose a numerical model of corneal biomechanics, based on the stromal microstructure. Our novel anisotropic constitutive material law features a probabilistic weighting approach to model collagen fiber distribution as observed on human cornea by Xray scattering analysis (Aghamohammadzadeh et. al., Structure, February 2004). Furthermore, collagen cross-linking was explicitly included in the strain energy function. Results showed that the proposed model is able to successfully reproduce both inflation and extensiometry experimental data (Elsheikh et. al., Curr Eye Res, 2007; Elsheikh et. al., Exp Eye Res, May 2008). In addition, the mechanical properties calculated for patients of different age groups (Group A: 65-79 years; Group B: 80-95 years) demonstrate an increased collagen cross-linking, and a decrease in collagen fiber elasticity from younger to older specimen. These findings correspond to what is known about maturing fibrous biological tissue. Since the presented model can handle different loading situations and includes the anisotropic distribution of collagen fibers, it has the potential to simulate clinical procedures involving nonsymmetrical tissue interventions. In the future, such mechanical model can be used to improve surgical planning and the design of next generation ophthalmic devices.

Projected 21st century decrease in marine productivity: a multi-model analysis

Changes in marine net primary productivity (PP) and export of particulate organic carbon (EP) are projected over the 21st century with four global coupled carbon cycle-climate models. These include representations of marine ecosystems and the carbon cycle of different structure and complexity. All four models show a decrease in global mean PP and EP between 2 and 20% by 2100 relative to preindustrial conditions, for the SRES A2 emission scenario. Two different regimes for productivity changes are consistently identified in all models. The first chain of mechanisms is dominant in the low- and mid-latitude ocean and in the North Atlantic: reduced input of macro-nutrients into the euphotic zone related to enhanced stratification, reduced mixed layer depth, and slowed circulation causes a decrease in macro-nutrient concentrations and in PP and EP. The second regime is projected for parts of the Southern Ocean: an alleviation of light and/or temperature limitation leads to an increase in PP and EP as productivity is fueled by a sustained nutrient input. A region of disagreement among the models is the Arctic, where three models project an increase in PP while one model projects a decrease. Projected changes in seasonal and interannual variability are modest in most regions. Regional model skill metrics are proposed to generate multi-model mean fields that show an improved skill in representing observation-based estimates compared to a simple multi-model average. Model results are compared to recent productivity projections with three different algorithms, usually applied to infer net primary production from satellite observations.

Global Search for Minimum Energy (H2O)n Clusters, n = 3−5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative ΔG° values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3−5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20 736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

Economic Policies to Address the Environmental Consequences of Global Reuse

This paper summarizes a two-country model that solves for optimal tax rates to achieve efficiency in an economy with international trade in used consumer electronics. If only the developed nation can tax the disposal of e-waste, then the global Pareto Optimum can be obtained by either imposing an import tariff on used consumer electronics or subsidizing the return of e-waste for disposal in the developed country. The global Pareto Optimum can also be obtained by reducing the disposal tax in the developed country to a level below the external marginal cost of disposal should no other policy option be available.

Detection of Internal Defects in Concrete Members Using Global Vibration Characteristics

Rock-pocket and honeycomb defects impair overall stiffness, accelerate aging, reduce service life, and cause structural problems in hardened concrete members. Traditional methods for detecting such deficient volumes involve visual observations or localized nondestructive methods, which are labor-intensive, time-consuming, highly sensitive to test conditions, and require knowledge of and accessibility to defect locations. The authors propose a vibration response-based nondestructive technique that combines experimental and numerical methodologies for use in identifying the location and severity of internal defects of concrete members. The experimental component entails collecting mode shape curvatures from laboratory beam specimens with size-controlled rock pocket and honeycomb defects, and the numerical component entails simulating beam vibration response through a finite element (FE) model parameterized with three defect-identifying variables indicating location (x, coordinate along the beam length) and severity of damage (alpha, stiffness reduction and beta, mass reduction). Defects are detected by comparing the FE model predictions to experimental measurements and inferring the low number of defect-identifying variables. This method is particularly well-suited for rapid and cost-effective quality assurance for precast concrete members and for inspecting concrete members with simple geometric forms.

Global Reuse and Optimal Waste Policy

Electronic waste generated from the consumption of durable goods in developed countries is often exported to underdeveloped countries for reuse, recycling and disposal with unfortunate environmental consequences. The lack of efficient disposal policies within developing nations coupled with global free trade agreements make it difficult for consumers to internalize these costs. This paper develops a two-country model, one economically developed and the other underdeveloped, to solve for optimal tax policies necessary to achieve the efficient allocation of economic resources in an economy with a durable good available for global reuse without policy measures in the underdeveloped country. A tax in the developed country on purchases of the new durable good combined with a waste tax set below the domestic external cost of disposal is sufficient for global efficiency. The implication of allowing free global trade in electronic waste is also examined, where optimal policy resembles a global deposit-refund system.

Constraining global methane emissions and uptake by ecosystems

Natural methane (CH4) emissions from wet ecosystems are an important part of today's global CH4 budget. Climate affects the exchange of CH4 between ecosystems and the atmosphere by influencing CH4 production, oxidation, and transport in the soil. The net CH4 exchange depends on ecosystem hydrology, soil and vegetation characteristics. Here, the LPJ-WHyMe global dynamical vegetation model is used to simulate global net CH4 emissions for different ecosystems: northern peatlands (45°–90° N), naturally inundated wetlands (60° S–45° N), rice agriculture and wet mineral soils. Mineral soils are a potential CH4 sink, but can also be a source with the direction of the net exchange depending on soil moisture content. The geographical and seasonal distributions are evaluated against multi-dimensional atmospheric inversions for 2003–2005, using two independent four-dimensional variational assimilation systems. The atmospheric inversions are constrained by the atmospheric CH4 observations of the SCIAMACHY satellite instrument and global surface networks. Compared to LPJ-WHyMe the inversions result in a~significant reduction in the emissions from northern peatlands and suggest that LPJ-WHyMe maximum annual emissions peak about one month late. The inversions do not put strong constraints on the division of sources between inundated wetlands and wet mineral soils in the tropics. Based on the inversion results we diagnose model parameters in LPJ-WHyMe and simulate the surface exchange of CH4 over the period 1990–2008. Over the whole period we infer an increase of global ecosystem CH4 emissions of +1.11 Tg CH4 yr−1, not considering potential additional changes in wetland extent. The increase in simulated CH4 emissions is attributed to enhanced soil respiration resulting from the observed rise in land temperature and in atmospheric carbon dioxide that were used as input. The long-term decline of the atmospheric CH4 growth rate from 1990 to 2006 cannot be fully explained with the simulated ecosystem emissions. However, these emissions show an increasing trend of +3.62 Tg CH4 yr−1 over 2005–2008 which can partly explain the renewed increase in atmospheric CH4 concentration during recent years.

A Computational Model for Processing Lexical Information"

The group analysed some syntactic and phonological phenomena that presuppose the existence of interrelated components within the lexicon, which motivate the assumption that there are some sublexicons within the global lexicon of a speaker. This result is confirmed by experimental findings in neurolinguistics. Hungarian speaking agrammatic aphasics were tested in several ways, the results showing that the sublexicon of closed-class lexical items provides a highly automated complex device for processing surface sentence structure. Analysing Hungarian ellipsis data from a semantic-syntactic aspect, the group established that the lexicon is best conceived of being as split into at least two main sublexicons: the store of semantic-syntactic feature bundles and a separate store of sound forms. On this basis they proposed a format for representing open-class lexical items whose meanings are connected via certain semantic relations. They also proposed a new classification of verbs to account for the contribution of the aspectual reading of the sentence depending on the referential type of the argument, and a new account of the syntactic and semantic behaviour of aspectual prefixes. The partitioned sets of lexical items are sublexicons on phonological grounds. These sublexicons differ in terms of phonotactic grammaticality. The degrees of phonotactic grammaticality are tied up with the problem of psychological reality, of how many degrees of this native speakers are sensitive to. The group developed a hierarchical construction network as an extension of the original General Inheritance Network formalism and this framework was then used as a platform for the implementation of the grammar fragments.

Effects of thoracic epidural anesthesia on hemodynamics and global oxygen transport in ovine endotoxemia

Besides providing effective analgesia, thoracic epidural anesthesia (TEA) has been shown to decrease perioperative morbidity and mortality. Because of its vasodilatory properties in association with the sympathetic blockade, however, TEA may potentially aggravate cardiovascular dysfunctions resulting from sepsis and systemic inflammatory response syndrome. The objective of the present study was to assess the effects of TEA on hemodynamics, global oxygen transport, and renal function in ovine endotoxemia. After a baseline measurement in healthy sheep (n = 18), Salmonella typhosa endotoxin was centrally infused at incremental doses to induce and maintain a hypotensive-hypodynamic circulation using an established protocol. The animals were then randomly assigned to one of two groups. In the treatment group, continuous TEA was initiated with 0.1 mL.kg of 0.125% bupivacaine at the onset of endotoxemia and maintained with 0.1 mL.kg.h. In the control group, the same amount of isotonic sodium chloride solution was injected through the epidural catheter. In the animals surviving the entire experiment (n = 7 per group), cardiac index and mean arterial pressure decreased in a dose-dependent manner during endotoxin infusion. In the TEA group, neither systemic hemodynamics nor global oxygen transport were impaired beyond the changes caused by endotoxemia itself. Urinary output was increased in the TEA group as compared with the control group (P < 0.05). In this model of endotoxic shock, TEA improved renal perfusion without affecting cardiopulmonary hemodynamics and global oxygen transport.

A Functional-Based Distribution Diagnostic for a Linear Model with Correlated Outcomes: Technical Report

Despite the widespread popularity of linear models for correlated outcomes (e.g. linear mixed modesl and time series models), distribution diagnostic methodology remains relatively underdeveloped in this context. In this paper we present an easy-to-implement approach that lends itself to graphical displays of model fit. Our approach involves multiplying the estimated marginal residual vector by the Cholesky decomposition of the inverse of the estimated marginal variance matrix. Linear functions or the resulting "rotated" residuals are used to construct an empirical cumulative distribution function (ECDF), whose stochastic limit is characterized. We describe a resampling technique that serves as a computationally efficient parametric bootstrap for generating representatives of the stochastic limit of the ECDF. Through functionals, such representatives are used to construct global tests for the hypothesis of normal margional errors. In addition, we demonstrate that the ECDF of the predicted random effects, as described by Lange and Ryan (1989), can be formulated as a special case of our approach. Thus, our method supports both omnibus and directed tests. Our method works well in a variety of circumstances, including models having independent units of sampling (clustered data) and models for which all observations are correlated (e.g., a single time series).

Drug-induced respiratory depression: an integrated model of drug effects on the hypercapnic and hypoxic drive

Drug-induced respiratory depression is a common side effect of the agents used in anesthesia practice to provide analgesia and sedation. Depression of the ventilatory drive in the spontaneously breathing patient can lead to severe cardiorespiratory events and it is considered a primary cause of morbidity. Reliable predictions of respiratory inhibition in the clinical setting would therefore provide a valuable means to improve the safety of drug delivery. Although multiple studies investigated the regulation of breathing in man both in the presence and absence of ventilatory depressant drugs, a unified description of respiratory pharmacodynamics is not available. This study proposes a mathematical model of human metabolism and cardiorespiratory regulation integrating several isolated physiological and pharmacological aspects of acute drug-induced ventilatory depression into a single theoretical framework. The description of respiratory regulation has a parsimonious yet comprehensive structure with substantial predictive capability. Simulations relative to the synergistic interaction of the hypercarbic and hypoxic respiratory drive and the global effect of drugs on the control of breathing are in good agreement with published experimental data. Besides providing clinically relevant predictions of respiratory depression, the model can also serve as a test bed to investigate issues of drug tolerability and dose finding/control under non-steady-state conditions.