Viral dynamics in vivo: limitations on estimates of intracellular delay and virus decay.

Anti-viral drug treatment of human immunodeficiency virus type I (HIV-1) and hepatitis B virus (HBV) infections causes rapid reduction in plasma virus load. Viral decline occurs in several phases and provides information on important kinetic constants of virus replication in vivo and pharmacodynamical properties. We develop a mathematical model that takes into account the intracellular phase of the viral life-cycle, defined as the time between infection of a cell and production of new virus particles. We derive analytic solutions for the dynamics following treatment with reverse transcriptase inhibitors, protease inhibitors, or a combination of both. For HIV-1, our results show that the phase of rapid decay in plasma virus (days 2-7) allows precise estimates for the turnover rate of productively infected cells. The initial quasi-stationary phase (days 0-1) and the transition phase (days 1-2) are explained by the combined effects of pharmacological and intracellular delays, the clearance of free virus particles, and the decay of infected cells. Reliable estimates of the first three quantities are not possible from data on virus load only; such estimates require additional measurements. In contrast with HIV-1, for HBV our model predicts that frequent early sampling of plasma virus will lead to reliable estimates of the free virus half-life and the pharmacological properties of the administered drug. On the other hand, for HBV the half-life of infected cells cannot be estimated from plasma virus decay.

Existence and uniqueness of Q-processes with a given finite μ-invariant measure

Let Q be a stable and conservative Q-matrix over a countable state space S consisting of an irreducible class C and a single absorbing state 0 that is accessible from C. Suppose that Q admits a finite mu-subinvariant measure in on C. We derive necessary and sufficient conditions for there to exist a Q-process for which m is mu-invariant on C, as well as a necessary condition for the uniqueness of such a process.

Stress corrosion cracking and hydrogen embrittlement of an Al-Zn-Mg-Cu alloy

The age hardening, stress corrosion cracking (SCC) and hydrogen embrittlement (HE) of an Al-Zn-Mg-Cu 7175 alloy were investigated experimentally. There were two peak-aged states during ageing. For ageing at 413 K, the strength of the second peak-aged state was slightly higher than that of the first one, whereas the SCC susceptibility was lower, indicating that it is possible to heat treat 7175 to high strength and simultaneously to have high SCC resistance. The SCC susceptibility increased with increasing Mg segregation at the grain boundaries. Hydrogen embrittlement (HE) increased with increased hydrogen charging and decreased with increasing ageing time for the same hydrogen charging conditions. Computer simulations were carried out of (a) the Mg grain boundary segregation using the embedded atom method and (b) the effect of Mg and H segregation on the grain boundary strength using a quasi-chemical approach. The simulations showed that (a) Mg grain boundary segregation in Al-Zn-Mg-Cu alloys is spontaneous, (b) Mg segregation decreases the grain boundary strength, and (c) H embrittles the grain boundary more seriously than does Mg. Therefore, the SCC mechanism of Al-Zn-Mg Cu alloys is attributed to the combination of HE and Mg segregation induced grain boundary embrittlement. (C) 2004 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Extinction times for a general birth, death and catastrophe process

The birth, death and catastrophe process is an extension of the birth-death process that incorporates the possibility of reductions in population of arbitrary size. We will consider a general form of this model in which the transition rates are allowed to depend on the current population size in an arbitrary manner. The linear case, where the transition rates are proportional to current population size, has been studied extensively. In particular, extinction probabilities, the expected time to extinction, and the distribution of the population size conditional on nonextinction (the quasi-stationary distribution) have all been evaluated explicitly. However, whilst these characteristics are of interest in the modelling and management of populations, processes with linear rate coefficients represent only a very limited class of models. We address this limitation by allowing for a wider range of catastrophic events. Despite this generalisation, explicit expressions can still be found for the expected extinction times.

A stochastic metapopulation model accounting for habitat dynamics

A stochastic metapopulation model accounting for habitat dynamics is presented. This is the stochastic SIS logistic model with the novel aspect that it incorporates varying carrying capacity. We present results of Kurtz and Barbour, that provide deterministic and diffusion approximations for a wide class of stochastic models, in a form that most easily allows their direct application to population models. These results are used to show that a suitably scaled version of the metapopulation model converges, uniformly in probability over finite time intervals, to a deterministic model previously studied in the ecological literature. Additionally, they allow us to establish a bivariate normal approximation to the quasi-stationary distribution of the process. This allows us to consider the effects of habitat dynamics on metapopulation modelling through a comparison with the stochastic SIS logistic model and provides an effective means for modelling metapopulations inhabiting dynamic landscapes.

Stochastic models for mainland-island metapopulations in static and dynamic landscapes

This paper has three primary aims: to establish an effective means for modelling mainland-island metapopulations inhabiting a dynamic landscape: to investigate the effect of immigration and dynamic changes in habitat on metapopulation patch occupancy dynamics; and to illustrate the implications of our results for decision-making and population management. We first extend the mainland-island metapopulation model of Alonso and McKane [Bull. Math. Biol. 64:913-958,2002] to incorporate a dynamic landscape. It is shown, for both the static and the dynamic landscape models, that a suitably scaled version of the process converges to a unique deterministic model as the size of the system becomes large. We also establish that. under quite general conditions, the density of occupied patches, and the densities of suitable and occupied patches, for the respective models, have approximate normal distributions. Our results not only provide us with estimates for the means and variances that are valid at all stages in the evolution of the population, but also provide a tool for fitting the models to real metapopulations. We discuss the effect of immigration and habitat dynamics on metapopulations, showing that mainland-like patches heavily influence metapopulation persistence, and we argue for adopting measures to increase connectivity between this large patch and the other island-like patches. We illustrate our results with specific reference to examples of populations of butterfly and the grasshopper Bryodema tuberculata.

A novel signal processing method for intraoperative neurophysiological monitoring in spinal surgeries

Intraoperative neurophysiologic monitoring is an integral part of spinal surgeries and involves the recording of somatosensory evoked potentials (SSEP). However, clinical application of IONM still requires anywhere between 200 to 2000 trials to obtain an SSEP signal, which is excessive and introduces a significant delay during surgery to detect a possible neurological damage. The aim of this study is to develop a means to obtain the SSEP using a much less, twelve number of recordings. The preliminary step involved was to distinguish the SSEP with the ongoing brain activity. We first establish that the brain activity is indeed quasi-stationary whereas an SSEP is expected to be identical every time a trial is recorded. An algorithm was developed using Chebychev time windowing for preconditioning of SSEP trials to retain the morphological characteristics of somatosensory evoked potentials (SSEP). This preconditioning was followed by the application of a principal component analysis (PCA)-based algorithm utilizing quasi-stationarity of EEG on 12 preconditioned trials. A unique Walsh transform operation was then used to identify the position of the SSEP event. An alarm is raised when there is a 10% time in latency deviation and/or 50% peak-to-peak amplitude deviation, as per the clinical requirements. The algorithm shows consistency in the results in monitoring SSEP in up to 6-hour surgical procedures even under this significantly reduced number of trials. In this study, the analysis was performed on the data recorded in 29 patients undergoing surgery during which the posterior tibial nerve was stimulated and SSEP response was recorded from scalp. This method is shown empirically to be more clinically viable than present day approaches. In all 29 cases, the algorithm takes 4sec to extract an SSEP signal, as compared to conventional methods, which take several minutes. The monitoring process using the algorithm was successful and proved conclusive under the clinical constraints throughout the different surgical procedures with an accuracy of 91.5%. Higher accuracy and faster execution time, observed in the present study, in determining the SSEP signals provide a much improved and effective neurophysiological monitoring process.

Energy levels and radiative rates for Cr-like Cu VI and Zn VII

Energy levels and radiative rates (. A-values) for transitions in Cr-like Cu VI and Zn VII are reported. These data are determined in the quasi-relativistic approach (QR), by employing a very large configuration interaction (CI) expansion which is highly important for these ions. No radiative rates are available in the literature to compare with our results, but our calculated energies are in close agreement with those compiled by NIST and other available theoretical data, for a majority of the levels. The A-values (and resultant lifetimes) are listed for all significantly contributing E1, E2 and M1 radiative transitions among the energetically lowest 322 levels of each ion.

Decoherence models for discrete-time quantum walks and their application to neutral atom experiments

We discuss decoherence in discrete-time quantum walks in terms of a phenomenological model that distinguishes spin and spatial decoherence. We identify the dominating mechanisms that affect quantum-walk experiments realized with neutral atoms walking in an optical lattice. From the measured spatial distributions, we determine with good precision the amount of decoherence per step, which provides a quantitative indication of the quality of our quantum walks. In particular, we find that spin decoherence is the main mechanism responsible for the loss of coherence in our experiment. We also find that the sole observation of ballistic-instead of diffusive-expansion in position space is not a good indicator of the range of coherent delocalization. We provide further physical insight by distinguishing the effects of short- and long-time spin dephasing mechanisms. We introduce the concept of coherence length in the discrete-time quantum walk, which quantifies the range of spatial coherences. Unexpectedly, we find that quasi-stationary dephasing does not modify the local properties of the quantum walk, but instead affects spatial coherences. For a visual representation of decoherence phenomena in phase space, we have developed a formalism based on a discrete analogue of the Wigner function. We show that the effects of spin and spatial decoherence differ dramatically in momentum space.

Influence of the sea state on Mediterranean heavy precipitation: a case-study from HyMeX SOP1

Sea state can influence the turbulent air–sea exchanges, especially the momentum flux, by modifying the sea-surface roughness. The high-resolution non-hydrostatic convection-permitting model MESO-NH is used here to investigate the impact of a more realistic representation of the waves on heavy precipitation during the Intense Observation Period (IOP) 16a of the first HyMeX Special Observation Period (SOP1). Several quasi-stationary mesoscale convective systems developed over the western Mediterranean region, two of them over the sea, and resulted in heavy precipitation on the French and Italian coasts on 26 October 2012. Three different bulk parametrizations are tested in this study: a reference case (NOWAV) without any wave effect, a parametrization taking into account theoretical wave effects (WAV) and a last one with realistic wave characteristics from the MFWAM analyses (WAM). Using a realistic wave representation in WAM significantly increases the roughness length and the friction velocity with respect to NOWAV and WAV. The three MESO-NH sensitivity experiments of the IOP16a show that this surface-roughness increase in WAM generates higher momentum fluxes and directly impacts the low-level dynamics of the atmosphere, with a slowdown of the 10 m wind, when and where the wind speed exceeds 10 m s−1 and the sea state differs from the idealized one. The turbulent heat fluxes are not significantly influenced by the waves, these fluxes being controlled by the moisture content rather than by the wind speed in the simulations. Although the convective activity is globally well reproduced by all the simulations, the difference in the low-level dynamics of the atmosphere influences the localization of the simulated heavy precipitation. Objective evaluation of the daily rainfall amount and of the 10 m wind speed against the observations confirms the positive impact of the realistic wave representation on this simulation of heavy precipitation.

Estudo de modelos simplificados com interações de longo alcance no ensemble microcanônico

Tese (doutorado)Universidade de Brasília, Instituto de Física, Programa de Pós-Graduação em Física, 2015.

A geometric approach to stationary defect solutions in one space dimension

In this manuscript, we consider the impact of a small jump-type spatial heterogeneity on the existence of stationary localized patterns in a system of partial dierential equations in one spatial dimension...

Quasi Monte Carlo approach with uniform separation for solving of the rendering equation, ICNAAM 2007

This paper is directed to the advanced parallel Quasi Monte Carlo (QMC) methods for realistic image synthesis. We propose and consider a new QMC approach for solving the rendering equation with uniform separation. First, we apply the symmetry property for uniform separation of the hemispherical integration domain into 24 equal sub-domains of solid angles, subtended by orthogonal spherical triangles with fixed vertices and computable parameters. Uniform separation allows to apply parallel sampling scheme for numerical integration. Finally, we apply the stratified QMC integration method for solving the rendering equation. The superiority our QMC approach is proved.