How to run a campaign: Optimal control of SIS and SIR information epidemics

Information spreading in a population can be modeled as an epidemic. Campaigners (e.g., election campaign managers, companies marketing products or movies) are interested in spreading a message by a given deadline, using limited resources. In this paper, we formulate the above situation as an optimal control problem and the solution (using Pontryagin's Maximum Principle) prescribes an optimal resource allocation over the time of the campaign. We consider two different scenarios-in the first, the campaigner can adjust a direct control (over time) which allows her to recruit individuals from the population (at some cost) to act as spreaders for the Susceptible-Infected-Susceptible (SIS) epidemic model. In the second case, we allow the campaigner to adjust the effective spreading rate by incentivizing the infected in the Susceptible-Infected-Recovered (SIR) model, in addition to the direct recruitment. We consider time varying information spreading rate in our formulation to model the changing interest level of individuals in the campaign, as the deadline is reached. In both the cases, we show the existence of a solution and its uniqueness for sufficiently small campaign deadlines. For the fixed spreading rate, we show the effectiveness of the optimal control strategy against the constant control strategy, a heuristic control strategy and no control. We show the sensitivity of the optimal control to the spreading rate profile when it is time varying. (C) 2014 Elsevier Inc. All rights reserved.

A composition-dependent ``re-entrant'' crystallographic phase transition in the substitutional metal acetylacetonate complex (Cr1-xGax)(acac)(3)

The formation of a complete solid solution between acetylacetonate (acac) complexes of chromium and gallium, (Cr1-x,Ga-x)(acac)(3) for 0.1 composition to have a congruent melting point, making it a substitutional complex. Whereas the pure complexes (i.e. the end members of the solid solution, x = 0 and x = 1) are both centro-symmetric, a composition-dependent crystallographic phase transition to a non-centrosymmetric structure is found to occur for compositions with 0.4 < x < 0.9. Such a ``re-entrant'' crystallographic transition is interpreted to be due to the drive to overcome the disorder present in certain centrosymmetric chromium-rich compositions, by going over to a non-centrosymmetric structure with a doubling of the unit cell. The substitutional complex is shown to lead to a substitutional oxide with the beta-gallate structure. (C) 2014 Elsevier Ltd. All rights reserved.

Contribution of forest fire ash and plant litter decay on stream dissolved composition in a sub-humid tropical watershed (Mule Hole, Southern India)

The current understanding of wildfire effects on water chemistry is limited by the quantification of the elemental dissolution rates from ash and element release rate from the plant litter, as well as quantification of the specific ash contribution to stream water chemistry. The main objective of the study was to provide such knowledge through combination of experimental modelling, field data and end-member mixing analysis (EMMA) of wildfire impact on a watershed scale. The study concerns watershed effects of fire in the Indian subcontinent, a region that is typically not well represented in the fire science literature. In plant litter ash, major elements are either hosted in readily-soluble phases (K, Mg) such as salts, carbonates and oxides or in less-soluble carrier-phases (Si, Ca) such as amorphous silica, quartz and calcite. Accordingly, elemental release rates, inferred from ash leaching experiments in batch reactor, indicated that the element release into solution followed the order K > Mg > Na > Si > Ca. Experiments on plant litter leaching in mixed-flow reactor indicated two dissolution regimes: rapid, over the week and slower over the month. The mean dissolution rates at steady-state (R-ss) indicated that the release of major elements from plant litter followed the order Ca > Si > Cl > Mg > K > Na. R-ss for Si and Ca for tree leaves and herbaceous species are similar to those reported for boreal and European tree species and are higher than that from the dissolution of soil clay minerals. This identifies tropical plant litters as important source of Si and Ca for tropical surface waters. In the wildfire-impacted year 2004, the EMMA indicated that the streamflow composition (Ca, K, Mg, Na, Si, Cl) was controlled by four main sources: rainwater, throughfall, ash leaching and soil solution. The influence of the ash end-member was maximal early in the rainy season (the two first storm events) and decreased later in the rainy season, when the stream was dominated by the throughfall end-member. The contribution of plant litter decay to the streamwater composition for a year not impacted by wildfire is significant with estimated solute fluxes originating from this decay greatly exceed, for most major elements, the annual elemental dissolved fluxes at the Mule Hole watershed outlet. This highlighted the importance of solute retention and vegetation back uptake processes within the soil profile. Overall, the fire increased the mobility and export of major elements from the soils to the stream. It also shifted the vegetation-related contribution to the elemental fluxes at the watershed outlet from long-term (seasonal) to short-term (daily to monthly). (C) 2014 Elsevier B.V. All rights reserved.

Time-Optimal Convergence to a Rectilinear Path in the Presence of Wind

This paper considers the problem of determining the time-optimal path of a fixed-wing Miniature Air Vehicle (MAV), in the presence of wind. The MAV, which is subject to a bounded turn rate, is required to eventually converge to a straight line starting from a known initial position and orientation. Earlier work in the literature uses Pontryagin's Minimum Principle (PMP) to solve this problem only for the no-wind case. In contrast, the present work uses a geometric approach to solve the problem completely in the presence of wind. In addition, it also shows how PMP can be used to partially solve the problem. Using a 6-DOF model of a MAV the generated optimal path is tracked by an autopilot consisting of proportional-integral-derivative (PID) controllers. The simulation results show the path generation and tracking for cases with steady and time-varying wind. Some issues on real-time path planning are also addressed.

Modelling the effects of chorus species composition and caller density on acoustic masking interference in multispecies choruses of crickets and katydids

Natural multispecies acoustic choruses such as the dusk chorus of a tropical rain forest consist of simultaneously signalling individuals of different species whose calls travel through a common shared medium before reaching their `intended' receivers. This causes masking interference between signals and impedes signal detection, recognition and localization. The levels of acoustic overlap depend on a number of factors, including call structure, intensity, habitat-dependent signal attenuation and receiver tuning. In addition, acoustic overlaps should also depend on caller density and the species composition of choruses, including relative and absolute abundance of the different calling species. In this study, we used simulations to examine the effects of chorus species relative abundance and caller density on the levels of effective heterospecific acoustic overlap in multispecies choruses composed of the calls of five species of crickets and katydids that share the understorey of a rain forest in southern India. We found that on average species-even choruses resulted in higher levels of effective heterospecific acoustic overlap than choruses with strong dominance structures. This effect was found consistently across dominance levels ranging from 0.4 to 0.8 for larger choruses of forty individuals. For smaller choruses of twenty individuals, the effect was seen consistently for dominance levels of 0.6 and 0.8 but not 0.4. Effective acoustic overlap (EAO) increased with caller density but the manner and extent of increase depended both on the species' call structure and the acoustic context provided by the composition scenario. The Phaloria sp. experienced very low levels of EAO and was highly buffered to changes in acoustic context whereas other species experienced high FAO across contexts or were poorly buffered. These differences were not simply predictable from call structures. These simulation-based findings may have important implications for acoustic biodiversity monitoring and for the study of acoustic masking interference in natural environments. (C) 2013 Elsevier B.V. All rights reserved.

Effect of post-deposition annealing on composition and electrical properties of dc reactive magnetron sputtered Al2O3 thin films

We have investigated the effect of post- deposition annealing on the composition and electrical properties of alumina (Al2O3) thin films. Al2O3 were deposited on n-type Si < 100 >. substrates by dc reactive magnetron sputtering. The films were subjected to post- deposition annealing at 623, 823 and 1023 K in vacuum. X-ray photoelectron spectroscopy results revealed that the composition improved with post- deposition annealing, and the film annealed at 1023 K became stoichiometric with an O/Al atomic ratio of 1.49. Al/Al2O3/Si metal-oxide-semiconductor (MOS) structures were then fabricated, and a correlation between the dielectric constant epsilon(r) and interface charge density Q(i) with annealing conditions were studied. The dielectric constant of the Al2O3 thin films increased to 9.8 with post- deposition annealing matching the bulk value, whereas the oxide charge density decreased to 3.11 x 10(11) cm(-2.) Studies on current-voltage IV characteristics indicated ohmic and Schottky type of conduction at lower electric fields (<0.16 MV cm(-1)) and space charge limited conduction at higher electric fields.

Synthesis and Mechanism of Composition and Size Dependent Morphology Selection in Nanoparticles of Ag-Cu Alloys Processed by Laser Ablation Under Liquid Medium

We have synthesized Ag-Cu alloy nanoparticles of four different compositions by using the laser ablation technique with the target under aqueous medium. Following this, we report a morphological transition in the nanoparticles from a normal two-phase microstructure to a structure with random segregation and finally a core shell structure at small sizes as a function of Cu concentration. To illustrate the composition dependence of morphology, we report observations carried out on nanoparticles of two different sizes: similar to 5 and similar to 20 nm. The results could be rationalized through the thermodynamic modeling of free energy of phase mixing and wettability of the alloying phases.

Optimal Trajectory Generation for Convergence to a Rectilinear Path

This paper presents a strategy to determine the shortest path of a fixed-wing Miniature Air Vehicle (MAV), constrained by a bounded turning rate, to eventually fly along a given straight line, starting from an arbitrary but known initial position and orientation. Unlike the work available in the literature that solves the problem using the Pontryagin's Minimum Principle (PMP) the trajectory generation algorithm presented here considers a geometrical approach which is intuitive and easy to understand. This also computes the explicit solution for the length of the optimal path as a function of the initial configuration. Further, using a 6-DOF model of a MAV the generated optimal path is tracked by an autopilot consisting of proportional-integral-derivative (PID) controllers. The simulation results show the path generation and tracking for different cases.

Incoherence-based optimal selection of independent measurements in diffuse optical tomography

An optimal measurement selection strategy based on incoherence among rows (corresponding to measurements) of the sensitivity (or weight) matrix for the near infrared diffuse optical tomography is proposed. As incoherence among the measurements can be seen as providing maximum independent information into the estimation of optical properties, this provides high level of optimization required for knowing the independency of a particular measurement on its counterparts. The proposed method was compared with the recently established data-resolution matrix-based approach for optimal choice of independent measurements and shown, using simulated and experimental gelatin phantom data sets, to be superior as it does not require an optimal regularization parameter for providing the same information. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Optimal control of information epidemics modeled as Maki Thompson rumors

We model the spread of information in a homogeneously mixed population using the Maki Thompson rumor model. We formulate an optimal control problem, from the perspective of single campaigner, to maximize the spread of information when the campaign budget is fixed. Control signals, such as advertising in the mass media, attempt to convert ignorants and stiflers into spreaders. We show the existence of a solution to the optimal control problem when the campaigning incurs non-linear costs under the isoperimetric budget constraint. The solution employs Pontryagin's Minimum Principle and a modified version of forward backward sweep technique for numerical computation to accommodate the isoperimetric budget constraint. The techniques developed in this paper are general and can be applied to similar optimal control problems in other areas. We have allowed the spreading rate of the information epidemic to vary over the campaign duration to model practical situations when the interest level of the population in the subject of the campaign changes with time. The shape of the optimal control signal is studied for different model parameters and spreading rate profiles. We have also studied the variation of the optimal campaigning costs with respect to various model parameters. Results indicate that, for some model parameters, significant improvements can be achieved by the optimal strategy compared to the static control strategy. The static strategy respects the same budget constraint as the optimal strategy and has a constant value throughout the campaign horizon. This work finds application in election and social awareness campaigns, product advertising, movie promotion and crowdfunding campaigns. (C) 2014 Elsevier B.V. All rights reserved.

A Stochastic Chemical Dynamic Approach to Correlate Autoimmunity and Optimal Vitamin-D Range

Motivated by several recent experimental observations that vitamin-D could interact with antigen presenting cells (APCs) and T-lymphocyte cells (T-cells) to promote and to regulate different stages of immune response, we developed a coarse grained but general kinetic model in an attempt to capture the role of vitamin-D in immunomodulatory responses. Our kinetic model, developed using the ideas of chemical network theory, leads to a system of nine coupled equations that we solve both by direct and by stochastic (Gillespie) methods. Both the analyses consistently provide detail information on the dependence of immune response to the variation of critical rate parameters. We find that although vitamin-D plays a negligible role in the initial immune response, it exerts a profound influence in the long term, especially in helping the system to achieve a new, stable steady state. The study explores the role of vitamin-D in preserving an observed bistability in the phase diagram (spanned by system parameters) of immune regulation, thus allowing the response to tolerate a wide range of pathogenic stimulation which could help in resisting autoimmune diseases. We also study how vitamin-D affects the time dependent population of dendritic cells that connect between innate and adaptive immune responses. Variations in dose dependent response of anti-inflammatory and pro-inflammatory T-cell populations to vitamin-D correlate well with recent experimental results. Our kinetic model allows for an estimation of the range of optimum level of vitamin-D required for smooth functioning of the immune system and for control of both hyper-regulation and inflammation. Most importantly, the present study reveals that an overdose or toxic level of vitamin-D or any steroid analogue could give rise to too large a tolerant response, leading to an inefficacy in adaptive immune function.

Optimal sequential wireless relay placement on a random lattice path

Our work is motivated by impromptu (or ``as-you-go'') deployment of wireless relay nodes along a path, a need that arises in many situations. In this paper, the path is modeled as starting at the origin (where there is the data sink, e.g., the control center), and evolving randomly over a lattice in the positive quadrant. A person walks along the path deploying relay nodes as he goes. At each step, the path can, randomly, either continue in the same direction or take a turn, or come to an end, at which point a data source (e.g., a sensor) has to be placed, that will send packets to the data sink. A decision has to be made at each step whether or not to place a wireless relay node. Assuming that the packet generation rate by the source is very low, and simple link-by-link scheduling, we consider the problem of sequential relay placement so as to minimize the expectation of an end-to-end cost metric (a linear combination of the sum of convex hop costs and the number of relays placed). This impromptu relay placement problem is formulated as a total cost Markov decision process. First, we derive the optimal policy in terms of an optimal placement set and show that this set is characterized by a boundary (with respect to the position of the last placed relay) beyond which it is optimal to place the next relay. Next, based on a simpler one-step-look-ahead characterization of the optimal policy, we propose an algorithm which is proved to converge to the optimal placement set in a finite number of steps and which is faster than value iteration. We show by simulations that the distance threshold based heuristic, usually assumed in the literature, is close to the optimal, provided that the threshold distance is carefully chosen. (C) 2014 Elsevier B.V. All rights reserved.

Use of Composition-Graded Bi-Electrolyte Cells for Thermodynamic Studies on Lanthanum Aluminates

Innovative bi-electrolyte solid-state cells incorporating single crystal CaF2 and composition-graded solid electrolyte (LaF3) y (CaF2) 1-y (y = 0 to 0.32) were used for measurement of the standard Gibbs energy of formation of hexagonal La0.885Al11.782O19 and cubic LaAlO3 from component binary oxides La2O3 and alpha-Al2O3 in the temperature range from 875 to 1175 K. The cells were designed based on experimentally verified relevant phase relations in the systems La2O3-Al2O3LaF3 and CaF2-LaF3. The results can be summarized as: 5.891 alpha-Al2O3 + 0.4425 La2O3 (A-rare earth)-> La0.885Al11.782O19 (hex), Delta G(f(ox))(degrees)(+/- 2005)/Jmol(-1) = -80982 + 7.313(T/K); 1/2 La2O3 (A-rare earth) + 1/2 a-Al2O3 -> LaAlO3 (cubic), Delta G(f(ox))(degrees)(+/- 2100)/Jmol(-1) = -59810 + 4.51(T/K). Electron probe microanalysis was used to ascertain the non-stoichiometric range of the hexaaluminate phase. The results are critically analyzed in the light of earlier electrochemical measurements. Several imperfections in the electrochemical cells used by former investigators are identified. Data obtained in the study for LaAlO3 are consistent with calorimetric enthalpy of formation and entropy derived from heat capacity data. Estimated are the standard entropy and the standard enthalpy of formation from elements of hexagonal La0.885Al11.782O19 and rhombohedral LaAlO3 at 298.15 K. c 2014 The Electrochemical Society. All rights reserved.

Multimodal therapy for complete regression of malignant melanoma using constrained nonlinear optimal dynamic inversion

Using a realistic nonlinear mathematical model for melanoma dynamics and the technique of optimal dynamic inversion (exact feedback linearization with static optimization), a multimodal automatic drug dosage strategy is proposed in this paper for complete regression of melanoma cancer in humans. The proposed strategy computes different drug dosages and gives a nonlinear state feedback solution for driving the number of cancer cells to zero. However, it is observed that when tumor is regressed to certain value, then there is no need of external drug dosages as immune system and other therapeutic states are able to regress tumor at a sufficiently fast rate which is more than exponential rate. As model has three different drug dosages, after applying dynamic inversion philosophy, drug dosages can be selected in optimized manner without crossing their toxicity limits. The combination of drug dosages is decided by appropriately selecting the control design parameter values based on physical constraints. The process is automated for all possible combinations of the chemotherapy and immunotherapy drug dosages with preferential emphasis of having maximum possible variety of drug inputs at any given point of time. Simulation study with a standard patient model shows that tumor cells are regressed from 2 x 107 to order of 105 cells because of external drug dosages in 36.93 days. After this no external drug dosages are required as immune system and other therapeutic states are able to regress tumor at greater than exponential rate and hence, tumor goes to zero (less than 0.01) in 48.77 days and healthy immune system of the patient is restored. Study with different chemotherapy drug resistance value is also carried out. (C) 2014 Elsevier Ltd. All rights reserved.