Prognostic value of biochemical liver parameters in neonatal sepsis-associated cholestasis

Aims The aim of the study was to evaluate the significance of total bilirubin, aspartate transaminase (AST), alanine transaminase and gamma-glutamyltransferase (GGT) for predicting outcome in sepsis-associated cholestasis. Methods: A retrospective cohort review of the hospital records was performed in 181 neonates admitted to the Neonatal Care Unit. A comparison was performed between subjects with low and high liver values based on cut-off values from ROC analysis. We defined poor prognosis to be when a subject had prolonged cholestasis of more than 3.5 months, developed severe sepsis, septic shock or had a fatal outcome. Results: The majority of the subjects were male (56%), preterm (56%) and had early onset sepsis (73%). The poor prognosis group had lower initial values of GGT compared with the good prognosis group (P = 0.003). Serum GGT (cut-off value of 85.5 U/L) and AST (cut-off value of 51 U/L) showed significant correlation with the outcome following multivariate analysis. The odds ratio (OR) of low GGT and high AST were OR 4.3 (95% CI:1.6 to11.8) and OR 2.9 (95% CI:1.1 to 8), respectively, for poor prognosis. In subjects with normal AST values, those with low GGT value had relative risk of 2.52 (95% CI:1.4 to 3.5) for poorer prognosis compared with those with normal or high GGT. Conclusion: Serum GGT and AST values can be used to predict the prognosis of patients with sepsis-associated cholestasis

Approximate solution of u'=-A2u using a relation between exp(-tA2) and exp(itA)

Let A be a positive definite operator in a Hilbert space and consider the initial value problem for u(t) = -A(2)u. Using a representation of the semigroup exp(-A(2)t) in terms of the group exp(iAt) we express u in terms of the solution of the standard heat equation w(t) = W-yy, with initial values v solving the initial value problem for v(y) = iAv. This representation is used to construct a method for approximating u in terms of approximations of v. In the case that A is a 2(nd) order elliptic operator the method is combined with finite elements in the spatial variable and then reduces the solution of the 4(th) order equation for u to that of the 2(nd) order equation for v, followed by the solution of the heat equation in one space variable.

Beyond fractional derivatives: local approximation of other convolution integrals

Dynamic systems involving convolution integrals with decaying kernels, of which fractionally damped systems form a special case, are non-local in time and hence infinite dimensional. Straightforward numerical solution of such systems up to time t needs O(t(2)) computations owing to the repeated evaluation of integrals over intervals that grow like t. Finite-dimensional and local approximations are thus desirable. We present here an approximation method which first rewrites the evolution equation as a coupled in finite-dimensional system with no convolution, and then uses Galerkin approximation with finite elements to obtain linear, finite-dimensional, constant coefficient approximations for the convolution. This paper is a broad generalization, based on a new insight, of our prior work with fractional order derivatives (Singh & Chatterjee 2006 Nonlinear Dyn. 45, 183-206). In particular, the decaying kernels we can address are now generalized to the Laplace transforms of known functions; of these, the power law kernel of fractional order differentiation is a special case. The approximation can be refined easily. The local nature of the approximation allows numerical solution up to time t with O(t) computations. Examples with several different kernels show excellent performance. A key feature of our approach is that the dynamic system in which the convolution integral appears is itself approximated using another system, as distinct from numerically approximating just the solution for the given initial values; this allows non-standard uses of the approximation, e. g. in stability analyses.

Primordial Perturbations from a Self-interacting Curvaton

Inflation is a period of accelerated expansion in the very early universe, which has the appealing aspect that it can create primordial perturbations via quantum fluctuations. These primordial perturbations have been observed in the cosmic microwave background, and these perturbations also function as the seeds of all large-scale structure in the universe. Curvaton models are simple modifications of the standard inflationary paradigm, where inflation is driven by the energy density of the inflaton, but another field, the curvaton, is responsible for producing the primordial perturbations. The curvaton decays after inflation as ended, where the isocurvature perturbations of the curvaton are converted into adiabatic perturbations. Since the curvaton must decay, it must have some interactions. Additionally realistic curvaton models typically have some self-interactions. In this work we consider self-interacting curvaton models, where the self-interaction is a monomial in the potential, suppressed by the Planck scale, and thus the self-interaction is very weak. Nevertheless, since the self-interaction makes the equations of motion non-linear, it can modify the behaviour of the model very drastically. The most intriguing aspect of this behaviour is that the final properties of the perturbations become highly dependent on the initial values. Departures of Gaussian distribution are important observables of the primordial perturbations. Due to the non-linearity of the self-interacting curvaton model and its sensitivity to initial conditions, it can produce significant non-Gaussianity of the primordial perturbations. In this work we investigate the non-Gaussianity produced by the self-interacting curvaton, and demonstrate that the non-Gaussianity parameters do not obey the analytically derived approximate relations often cited in the literature. Furthermore we also consider a self-interacting curvaton with a mass in the TeV-scale. Motivated by realistic particle physics models such as the Minimally Supersymmetric Standard Model, we demonstrate that a curvaton model within the mass range can be responsible for the observed perturbations if it can decay late enough.

A parallel multistep predictor-corrector algorithm for solving ordinary differential equations

In this paper we give a generalized predictor-corrector algorithm for solving ordinary differential equations with specified initial values. The method uses multiple correction steps which can be carried out in parallel with a prediction step. The proposed method gives a larger stability interval compared to the existing parallel predictor-corrector methods. A method has been suggested to implement the algorithm in multiple processor systems with efficient utilization of all the processors.

Surface Chemistry of Thiobacillus ferrooxidans Relevant to Adhesion on Mineral Surfaces

Thiobacillus ferrooxidans cells grown on sulfur, pyrite, and chalcopyrite exhibit greater hydrophobicity than ferrous ion-grown cells. The isoelectric points of sulfur-, pyrite-, and chalcopyrite-grown cells were observed to be at a pH higher than that for ferrous ion-grown cells. Microbe-mineral interactions result in change in the surface chemistry of the organism as well as that of the minerals with which it has interacted. Sulfur, pyrite, and chalcopyrite after interaction with T. ferrooxidans exhibited a significant shift in their isoelectric points from the initial values exhibited by uninteracted minerals. With antibodies raised against sulfur-grown T. ferrooxidans, pyrite- and chalcopyrite-grown cells showed immunoreactivity, whereas ferrous ion-grown cells failed to do so. Fourier transform infrared spectroscopy of sulfur-grown cells suggested that a proteinaceous new cell surface appendage synthesized in mineral-grown cells brings about adhesion to the solid mineral substrates. Such an appendage was found to be absent in ferrous ion-grown cells as it is not required during growth in liquid substrates.

Accelerated gradient based diffuse optical tomographic image reconstruction

Purpose: Fast reconstruction of interior optical parameter distribution using a new approach called Broyden-based model iterative image reconstruction (BMOBIIR) and adjoint Broyden-based MOBIIR (ABMOBIIR) of a tissue and a tissue mimicking phantom from boundary measurement data in diffuse optical tomography (DOT). Methods: DOT is a nonlinear and ill-posed inverse problem. Newton-based MOBIIR algorithm, which is generally used, requires repeated evaluation of the Jacobian which consumes bulk of the computation time for reconstruction. In this study, we propose a Broyden approach-based accelerated scheme for Jacobian computation and it is combined with conjugate gradient scheme (CGS) for fast reconstruction. The method makes explicit use of secant and adjoint information that can be obtained from forward solution of the diffusion equation. This approach reduces the computational time many fold by approximating the system Jacobian successively through low-rank updates. Results: Simulation studies have been carried out with single as well as multiple inhomogeneities. Algorithms are validated using an experimental study carried out on a pork tissue with fat acting as an inhomogeneity. The results obtained through the proposed BMOBIIR and ABMOBIIR approaches are compared with those of Newton-based MOBIIR algorithm. The mean squared error and execution time are used as metrics for comparing the results of reconstruction. Conclusions: We have shown through experimental and simulation studies that Broyden-based MOBIIR and adjoint Broyden-based methods are capable of reconstructing single as well as multiple inhomogeneities in tissue and a tissue-mimicking phantom. Broyden MOBIIR and adjoint Broyden MOBIIR methods are computationally simple and they result in much faster implementations because they avoid direct evaluation of Jacobian. The image reconstructions have been carried out with different initial values using Newton, Broyden, and adjoint Broyden approaches. These algorithms work well when the initial guess is close to the true solution. However, when initial guess is far away from true solution, Newton-based MOBIIR gives better reconstructed images. The proposed methods are found to be stable with noisy measurement data. (C) 2011 American Association of Physicists in Medicine. DOI: 10.1118/1.3531572]

Direct Cartesian-Space Solutions of Generalized Bloch Equations in the Rotating Frame

Analytical solutions of the generalized Bloch equations for an arbitrary set of initial values of the x, y, and z magnetization components are given in the rotating frame. The solutions involve the decoupling of the three coupled differential equations such that a third-order differential equation in each magnetization variable is obtained. In contrast to the previously reported solutions given by Torrey, the present attempt paves the way for more direct physical insight into the behavior of each magnetization component. Special cases have been discussed that highlight the utility of the general solutions. Representative trajectories of magnetization components are given, illustrating their behavior with respect to the values of off-resonance and initial conditions. (C) 1995 Academic Press, Inc.

Wavelet based spectral finite element modelling and detection of de-lamination in composite beams

In this paper, a model for composite beam with embedded de-lamination is developed using the wavelet based spectral finite element (WSFE) method particularly for damage detection using wave propagation analysis. The simulated responses are used as surrogate experimental results for the inverse problem of detection of damage using wavelet filtering. The WSFE technique is very similar to the fast fourier transform (FFT) based spectral finite element (FSFE) except that it uses compactly supported Daubechies scaling function approximation in time. Unlike FSFE formulation with periodicity assumption, the wavelet-based method allows imposition of initial values and thus is free from wrap around problems. This helps in analysis of finite length undamped structures, where the FSFE method fails to simulate accurate response. First, numerical experiments are performed to study the effect of de-lamination on the wave propagation characteristics. The responses are simulated for different de-lamination configurations for both broad-band and narrow-band excitations. Next, simulated responses are used for damage detection using wavelet analysis.

Estimation of Parameters for Battery Storage Models

This paper presents a comprehensive and robust strategy for the estimation of battery model parameters from noise corrupted data. The deficiencies of the existing methods for parameter estimation are studied and the proposed parameter estimation strategy improves on earlier methods by working optimally for low as well as high discharge currents, providing accurate estimates even under high levels of noise, and with a wide range of initial values. Testing on different data sets confirms the performance of the proposed parameter estimation strategy.

Propiedades básicas de materiales soporte y parámetros de control en biofiltración.

[ES]En la problemática medioambiental generada por la contaminación atmosférica hay tres aspectos que marcan las directrices de actuación institucional: la presión social, la legislación vigente y la tecnología disponible. En cuanto a este último aspecto, la biofiltración es una tecnología eficaz, asequible y sostenible basada en la actividad biodegradadora de microorganismos específicos adheridos a la superficie de un material soporte que constituye el lecho del biofiltro. La elección de un material soporte adecuado es de especial importancia para asegurar el correcto funcionamiento de los biofiltros. Esta decisión está basada en las propiedades intrínsecas del material que deben ser analizadas previamente a su uso. En este proyecto se ha seleccionado cuales son estas propiedades básicas a partir de una revisión bibliográfica, destacando la capacidad de retención de humedad, superficie específica, porosidad, y estabilidad física y química del material. En este trabajo, también se han fijado los parámetros de control que deben ser medidos de forma rutinaria en un biofiltro para asegurar la eficacia del tratamiento de descontaminación. En base a la información bibliográfica recopilada, se ha concluido que los parámetros básicos son pH, temperatura, contenido de humedad del lecho y pérdidas de carga. A nivel experimental, se han medido las pérdidas de carga generadas en biofiltros empacados con tres materiales soportes que son de especial interés para una investigación posterior a desarrollar por el grupo Biofiltración de la Universidad del País Vasco (UPV/EHU). Durante el período de arranque inicial de los tres biofiltros, las pérdidas de carga medidas fueron muy bajas en los tres casos, aunque algo superiores cuando la alimentación era en sentido ascendente frente al flujo descendente. Sin embargo, cuando se midieron las pérdidas de carga bajo condiciones de encharcamiento del lecho, que es una de las situaciones más problemáticas en un biofiltro, todos los soportes presentaron un aumento puntual de la pérdida de carga aunque la tendencia en los días posteriores fue claramente descendente, recuperando valores de operación habituales. La comparativa del comportamiento frente a las pérdidas de carga permitió seleccionar el soporte más idóneo de los tres analizados, aunque los otros dos podrían ser alternativas viables en caso de sustitución.

A model of the pelagic ecossystem in the eastern tropical Pacific Ocean

English: Recent calls for a more holistic approach to fisheries management have motivated development of trophic mass-balance models of ecosystems that underlie fisheries production. We developed a model hypothesis of the pelagic ecosystem in the eastern tropical Pacific Ocean (ETP) to gain insight into the relationships among the various species in the system and to explore the ecological implications of alternative methods of harvesting tunas. We represented the biomasses of and fluxes between the principal elements in the ecosystem with Ecopath, and examined the ecosystem's dynamic, time-series behavior with Ecosim. We parameterized the model for 38 species or groups of species, and described the sources, justifications, assumptions, and revisions of our estimates of the various parameters, diet relations, fisheries landings, and fisheries discards in the model. We conducted sensitivity analyses with an intermediate version of the model, for both the Ecopath mass-balance and the dynamic trajectories predicted by Ecosim. The analysis showed that changes in the basic parameters for two components at middle trophic levels, Cephalopods and Auxis spp., exert the greatest influence on the system. When the Cephalopod Q/B and Auxis spp. P/B were altered from their initial values and the model was rebalanced, the trends of the biomass trajectories predicted by Ecosim were not sensitive, but the scaling was sensitive for several components. We described the review process the model was subjected to, which included reviews by the IATTC Purse-seine Bycatch Working Group and by a working group supported by the National Center for Ecological Analysis and Synthesis. We fitted the model to historical time series of catches per unit of effort and mortality rates for yellowfin and bigeye tunas in simulations that incorporated historical fishing effort and a climate driver to represent the effect of El Niño-Southern Oscillation-scale variation on the system. The model was designed to evaluate the possible ecological implications of fishing for tunas in various ways. We recognize that a model cannot possibly represent all the complexity of a pelagic ocean ecosystem, but we believe that the ETP model provides insight into the structure and function of the pelagic ETP. Spanish: Llamamientos recientes hacia un enfoque más holístico al ordenamiento de la pesca han motivado el desarrollo de modelos tróficos de balance de masas de los ecosistemas que sostienen la producción pesquera. Desarrollamos una hipótesis modelo del ecosistema pelágico en el Océano Pacífico oriental tropical (POT) con miras a mejorar los conocimientos de las relaciones entre las distintas especies en el sistema y explorar las implicaciones ecológicas de métodos alternativos de capturar atunes. Con Ecopath representamos las biomasas de los elementos principales en el ecosistema, y los flujos entre los mismos, y con Ecosim examinamos el comportamiento dinámico del ecosistema con el tiempo. Parametrizamos el modelo para 38 especies o grupos de especies (denominados “componentes” del modelo), y describimos las fuentes, justificaciones, supuestos, y revisiones de nuestras estimaciones de los distintos parámetros, relaciones basadas en dieta, capturas retenidas de las pesquerías, y descartes de las mismas en el modelo. Realizamos análisis de sensibilidad con una versión intermedia del modelo, para el balance de masas de Ecopath y las trayectorias dinámicas predichas por Ecosim también. El análisis demostró que cambios en los parámetros básicos para dos componentes en niveles tróficos medianos, Cefalópodos y Auxis spp., ejercieron la mayor influencia sobre el sistema. Cuando se alteraron el Q/B de los Cefalópodos y el P/B de los Auxis spp. de sus valores iniciales y se balanceó el modelo de nuevo, las tendencias de las trayectorias de la biomasa predichas por Ecosim no fueron sensibles, pero la escala fue sensible para varios componentes. Describimos el proceso de revisión al que fue sujeto el modelo, inclusive revisiones por el Grupo de Trabajo sobre Captura Incidental de la CIAT y un grupo de trabajo apoyado por el Centro Nacional para Síntesis y Análisis Ecológicos. Ajustamos el modelo a series de tiempo históricas de capturas por unidad de esfuerzo y tasas de mortalidad de atunes aleta amarilla y patudo en simulaciones que incorporaron esfuerzo de pesca histórico e impulsos climáticos para representar el efecto de variaciones a escala de El Niño-Oscilación del Sur sobre el sistema. El modelo fue diseñado para evaluar las posibles implicaciones ecológicas de la pesca atunera de varias formas. Reconocemos la imposibilidad de que el modelo represente toda la complejidad de un ecosistema oceánico pelágico, pero creemos que el modelo del POT mejora los conocimientos de la estructura y función del POT pelágico.

Development and transition of spiral wave in the coupled Hindmarsh-Rose neurons in two-dimensional space

The dynamics and the transition of spiral waves in the coupled Hindmarsh-Rose (H-R) neurons in two-dimensional space are investigated in the paper. It is found that the spiral wave can be induced and developed in the coupled HR neurons in two-dimensional space, with appropriate initial values and a parameter region given. However, the spiral wave could encounter instability when the intensity of the external current reaches a threshold value of 1.945. The transition of spiral wave is found to be affected by coupling intensity D and bifurcation parameter r. The spiral wave becomes sparse as the coupling intensity increases, while the spiral wave is eliminated and the whole neuronal system becomes homogeneous as the bifurcation parameter increases to a certain threshold value. Then the coupling action of the four sub-adjacent neurons, which is described by coupling coefficient D', is also considered, and it is found that the spiral wave begins to breakup due to the introduced coupling action from the sub-adjacent neurons (or sites) and together with the coupling action of the nearest-neighbour neurons, which is described by the coupling intensity D.

The study of banded spherulite patterns by coupled logistic map lattice model

Banded spherulite patterns are simulated in two dimensions by means of a coupled logistic map lattice model. Both target pattern and spiral pattern which have been proved to be existent experimentally in banded spherulite are obtained by choosing suitable parameters in the model. The simulation results also indicate that the band spacing is decreased with the increase of parameter mu in the logistic map and increased with the increase of the coupling parameter epsilon, which is quite similar to the results in some experiments. Moreover, the relationship between the parameters and the corresponding patterns is obtained, and the target patterns and spiral patterns are distinguished for a given group of initial values, which may guide the study of banded spherulite.

Estimation of individual phases from the four-phase structure invariants in the single isomorphous replacement case

The probability distribution of the four-phase invariants in the case of single isomorphous replacement has been developed to estimate some individual phases. An example of its application to obtain the phases having special values of 0, pi or +/-pi /2 is given for a known protein structure in space group P2(1)2(1)2(1). The phasing procedure includes the determination of starting phases and an iterative calculation. The initial values of starting phases, which are required by the formula, can be obtained from the estimate of one-phase seminvariants and by specifying the origin and enantiomorph. In addition, the calculations lead to two sets of possible phases for each type of reflection by assigning arbitrarily an initial phase value. The present method provides a possibility for the multisolution technique to increase greatly the number of known phases while keeping the number of the trials quite small.