Geometric Approach to Pontryagin's Maximum Principle

Since the second half of the 20th century, Pontryagin's Maximum Principle has been widely discussed and used as a method to solve optimal control problems in medicine, robotics, finance, engineering, astronomy. Here, we focus on the proof and on the understanding of this Principle, using as much geometric ideas and geometric tools as possible. This approach provides a better and clearer understanding of the Principle and, in particular, of the role of the abnormal extremals. These extremals are interesting because they do not depend on the cost function, but only on the control system. Moreover, they were discarded as solutions until the nineties, when examples of strict abnormal optimal curves were found. In order to give a detailed exposition of the proof, the paper is mostly self\textendash{}contained, which forces us to consider different areas in mathematics such as algebra, analysis, geometry.Comment: Final version. Minors changes have been made. 56 page

IBD Serology and Disease Outcomes in African Americans with Crohn's Disease

Backgrounds: Recent studies have identified the role of serologic markers in characterizing disease phenotype, location, complications, and severity among Northern Europeans (NE) with Crohn's disease (CD). However, very little is known about the role of serology in CD among African Americans (AA). Our study explored the relationship between serology and disease phenotype in AA with CD, while controlling for genetic ancestry. Methods: AAs with CD were enrolled as participants through multicenter collaborative efforts. Serological levels of IgA anti-Saccharomyces cervisiae antibody (ASCA), IgG ASCA, E. coli outermembrane porin C, anti-CBir1, and ANCA were measured using enzyme-linked immunosorbent assays. Genotyping was performed using Illumina immunochip technology; an admixture rate was calculated for each subject. Multiple imputation by chained equations was performed to account for data missing at random. Logistic regression was used to calculate adjusted odds ratio (OR) for associations between serological markers and both complicated disease and disease requiring surgery. Results: A total of 358 patients were included in the analysis. The majority of our patients had inflammatory, noncomplicated disease (58.4%), perianal disease (55.7%), and documented colonic inflammation (86.8%). On multivariable analysis, both IgG ASCA and OmpC were associated with complicated disease (OR, 2.67; 95% CI, 1.67-4.28; OR, 2.23; 95% CI, 1.41-3.53, respectively) and disease requiring surgery (OR, 2.51; 95% CI, 1.49-4.22; OR, 3.57; 95% CI, 2.12-6.00). NE admixture to the African genome did not have any associations or interactions in relation to clinical outcome. Conclusions: Our study comprises the largest cohort of AAs with CD. The utility of serological markers for the prognosis of CD in NE applies equally to AA populations

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)

Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

Relationship of edge localized mode burst times with divertor flux loop signal phase in JET

A phase relationship is identified between sequential edge localized modes (ELMs) occurrence times in a set of H-mode tokamak plasmas to the voltage measured in full flux azimuthal loops in the divertor region. We focus on plasmas in the Joint European Torus where a steady H-mode is sustained over several seconds, during which ELMs are observed in the Be II emission at the divertor. The ELMs analysed arise from intrinsic ELMing, in that there is no deliberate intent to control the ELMing process by external means. We use ELM timings derived from the Be II signal to perform direct time domain analysis of the full flux loop VLD2 and VLD3 signals, which provide a high cadence global measurement proportional to the voltage induced by changes in poloidal magnetic flux. Specifically, we examine how the time interval between pairs of successive ELMs is linked to the time-evolving phase of the full flux loop signals. Each ELM produces a clear early pulse in the full flux loop signals, whose peak time is used to condition our analysis. The arrival time of the following ELM, relative to this pulse, is found to fall into one of two categories: (i) prompt ELMs, which are directly paced by the initial response seen in the flux loop signals; and (ii) all other ELMs, which occur after the initial response of the full flux loop signals has decayed in amplitude. The times at which ELMs in category (ii) occur, relative to the first ELM of the pair, are clustered at times when the instantaneous phase of the full flux loop signal is close to its value at the time of the first ELM

A hybrid multi-compartment model of granuloma formation and T cell priming in tuberculosis

Tuberculosis is a worldwide health problem with 2 billion people infected with Mycobacterium tuberculosis (Mtb, the bacteria causing TB). The hallmark of infection is the emergence of organized structures of immune cells forming primarily in the lung in response to infection. Granulomas physically contain and immunologically restrain bacteria that cannot be cleared. We have developed several models that spatially characterize the dynamics of the host-mycobacterial interaction, and identified mechanisms that control granuloma formation and development. In particular, we published several agent-based models (ABMs) of granuloma formation in TB that include many subtypes of T cell populations, macrophages as well as key cytokine and chemokine effector molecules. These ABM studies emphasize the important role of T-cell related mechanisms in infection progression, such as magnitude and timing of T cell recruitment, and macrophage activation. In these models, the priming and recruitment of T cells from the lung draining lymph node (LN) was captured phenomenologically. In addition to these ABM studies, we have also developed several multi-organ models using ODEs to examine trafficking of cells between, for example, the lung and LN. While we can predict temporal dynamic behaviors, those models are not coupled to the spatial aspects of granuloma. To this end, we have developed a multi-organ model that is hybrid: an ABM for the lung compartment and a non-linear system of ODE representing the lymph node compartment. This hybrid multi-organ approach to study TB granuloma formation in the lung and immune priming in the LN allows us to dissect protective mechanisms that cannot be achieved using the single compartment or multi-compartment ODE system. The main finding of this work is that trafficking of important cells known as antigen presenting cells from the lung to the lymph node is a key control mechanism for protective immunity: the entire spectrum of infection outcomes can be regulated by key immune cell migration rates. Our hybrid multi-organ implementation suggests that effector CD4+ T cells can rescue the system from a persistent infection and lead to clearance once a granuloma is fully formed. This could be effective as an immunotherapy strategy for latently infected individuals. © 2011 Elsevier Ltd

Dynamics of HIV infection of CD4 T cells

We examine a model for the interaction of HIV with CD4+ T cells that considers four populations: uninfected T cells, latently infected T cells, actively infected T cells, and free virus. Using this model we show that many of the puzzling quantitative features of HIV infection can be explained simply. We also consider effects of AZT on viral growth and T-cell population dynamics. The model exhibits two steady states, an uninfected state in which no virus is present and an endemically infected state, in which virus and infected T cells are present. We show that if N, the number of infectious virions produced per actively infected T cell, is less a critical value, Ncrit, then the uninfected state is the only steady state in the nonnegative orthant, and this state is stable. For N > Ncrit, the uninfected state is unstable, and the endemically infected state can be either stable, or unstable and surrounded by a stable limit cycle. Using numerical bifurcation techniques we map out the parameter regimes of these various behaviors. oscillatory behavior seems to lie outside the region of biologically realistic parameter values. When the endemically infected state is stable, it is characterized by a reduced number of T cells compared with the uninfected state. Thus T-cell depletion occurs through the establishment of a new steady state. The dynamics of the establishment of this new steady state are examined both numerically and via the quasi-steady-state approximation. We develop approximations for the dynamics at early times in which the free virus rapidly binds to T cells, during an intermediate time scale in which the virus grows exponentially, and a third time scale on which viral growth slows and the endemically infected steady state is approached. Using the quasi-steady-state approximation the model can be simplified to two ordinary differential equations the summarize much of the dynamical behavior. We compute the level of T cells in the endemically infected state and show how that level varies with the parameters in the model. The model predicts that different viral strains, characterized by generating differing numbers of infective virions within infected T cells, can cause different amounts of T-cell depletion and generate depletion at different rates. Two versions of the model are studied. In one the source of T cells from precursors is constant, whereas in the other the source of T cells decreases with viral load, mimicking the infection and killing of T-cell precursors