A simulation-based algorithm for solving the resource-assignment problem in satellite telecommunication networks

This paper proposes an heuristic for the scheduling of capacity requests and the periodic assignment of radio resources in geostationary (GEO) satellite networks with star topology, using the Demand Assigned Multiple Access (DAMA) protocol in the link layer, and Multi-Frequency Time Division Multiple Access (MF-TDMA) and Adaptive Coding and Modulation (ACM) in the physical layer.

Analytical expressions for the flame front speed in the downward combustion of thin solid fuels and comparison to experiments

We derive analytical expressions for the propagation speed of downward combustion fronts of thin solid fuels with a background flow initially at rest. The classical combustion model for thin solid fuels that consists of five coupled reaction-convection-diffusion equations is here reduced into a single equation with the gas temperature as the single variable. For doing so we apply a two-zone combustion model that divides the system into a preheating region and a pyrolyzing region. The speed of the combustion front is obtained after matching the temperature and its derivative at the location that separates both regions.We also derive a simplified version of this analytical expression expected to be valid for a wide range of cases. Flame front velocities predicted by our analyticalexpressions agree well with experimental data found in the literature for a large variety of cases and substantially improve the results obtained from a previous well-known analytical expression

The Effect Of Induced Intraocular Straylight (cataract Simulation) On Sensitivity Measures Of Standard Automated Perimetry, Pulsar Perimetry And Moorfields Motion Displacement Test

Purpose: To investigate the effect of incremental increases in intraocular straylight on threshold measurements made by three modern forms of perimetry: Standard Automated Perimetry (SAP) using Octopus (Dynamic, G-Pattern), Pulsar Perimetry (PP) (TOP, 66 points) and the Moorfields Motion Displacement Test (MDT) (WEBS, 32 points).Methods: Four healthy young observers were recruited (mean age 26yrs [25yrs, 28yrs]), refractive correction [+2 D, -4.25D]). Five white opacity filters (WOF), each scattering light by different amounts were used to create incremental increases in intraocular straylight (IS). Resultant IS values were measured with each WOF and at baseline (no WOF) for each subject using a C-Quant Straylight Meter (Oculus, Wetzlar, Germany). A 25 yr old has an IS value of ~0.85 log(s). An increase of 40% in IS to 1.2log(s) corresponds to the physiological value of a 70yr old. Each WOFs created an increase in IS between 10-150% from baseline, ranging from effects similar to normal aging to those found with considerable cataract. Each subject underwent 6 test sessions over a 2-week period; each session consisted of the 3 perimetric tests using one of the five WOFs and baseline (both instrument and filter were randomised).Results: The reduction in sensitivity from baseline was calculated. A two-way ANOVA on mean change in threshold (where subjects were treated as rows in the block and each increment in fog filters was treated as column) was used to examine the effect of incremental increases in straylight. Both SAP (p<0.001) and Pulsar (p<0.001) were significantly affected by increases in straylight. The MDT (p=0.35) remained comparatively robust to increases in straylight.Conclusions: The Moorfields MDT measurement of threshold is robust to effects of additional straylight as compared to SAP and PP.

What Was I Thinking? Eye-Tracking Experiments Underscore the Bias that Architecture Exerts on Nuclear Grading in Prostate Cancer.

We previously reported that nuclear grade assignment of prostate carcinomas is subject to a cognitive bias induced by the tumor architecture. Here, we asked whether this bias is mediated by the non-conscious selection of nuclei that "match the expectation" induced by the inadvertent glance at the tumor architecture. 20 pathologists were asked to grade nuclei in high power fields of 20 prostate carcinomas displayed on a computer screen. Unknown to the pathologists, each carcinoma was shown twice, once before a background of a low grade, tubule-rich carcinoma and once before the background of a high grade, solid carcinoma. Eye tracking allowed to identify which nuclei the pathologists fixated during the 8 second projection period. For all 20 pathologists, nuclear grade assignment was significantly biased by tumor architecture. Pathologists tended to fixate on bigger, darker, and more irregular nuclei when those were projected before kigh grade, solid carcinomas than before low grade, tubule-rich carcinomas (and vice versa). However, the morphometric differences of the selected nuclei accounted for only 11% of the architecture-induced bias, suggesting that it can only to a small part be explained by the unconscious fixation on nuclei that "match the expectation". In conclusion, selection of « matching nuclei » represents an unconscious effort to vindicate the gravitation of nuclear grades towards the tumor architecture.

ASIC1a oligomerization structure-function of its extracellular vestibule and functional characterization of αS243Pβɣ ENaC

Quatre cristaux du canal ASIC1a ont été publiés et soutiennent une stoechiométrie trimérique. Cependant, ces données contredisant de précédentes analyses fonctionnelles effectuées sur des canaux de la même famille, notre intérêt fut porté sur l'oligomérisation d'ASIC1a. Dans ce sens, un nouvel essai couplant la méthode d'analyse par substitution de cystéines (SCAM) avec l'utilisation de réactifs sulfhydryls bifonctionnels (crosslinkers) a été mis en place. Le but étant de stabiliser, puis sélectionner les canaux fonctionnels, pour ensuite les séparer selon leur taille par SDS-PAGE. Grâce à cette technique, nous avons démontré que le complexe stabilisé a une taille coïncidant avec une organisation tétramérique. En plus de son oligomérisation, le chemin emprunté par les ions pour traverser le canal n'est pas clairement défini dans ces structures. De ce fait, utilisant une approche électrophysiologique, nous avons étudié le lien entre la structure et la fonction du vestibule extracellulaire d'ASIC1a. Dans ce but, nous nous sommes intéressés l'accessibilité de cystéines spécifiques localisées dans ce vestibule pour des réactifs méthanethiosulfonates (MTS). Ainsi, nous avons pu corréler les cinétiques de modification de ces cystéines par les MTS avec les effets sur le courant sodique, et donc avoir des informations supplémentaires sur la voie empruntée par les ions. De plus, la simulation informatique de liaison de ces réactifs illustre le remplissage total de ce vestibule. Fonctionnellement, cette interaction ne perturbe pas le passage de ions, c'est pourquoi il nous apparaît probable que le vestibule présente une taille plus large que celle illustrée par les cristaux. Dans un deuxième temps, notre intérêt fut porté sur ENaC. Ce canal est composé des trois sous-unités (a, ß et y) et est exprimé dans divers épithéliums, dont les tubules des reins. Il participe à l'homéostasie sodique et est essentiellement régulé par voie hormonale via l'aldostérone et la Vasopressine, mais également par des sérines protéases ou le Na+. Nous avons étudié la répercussion fonctionnelle de la mutation aS243P, découverte chez un nouveau-né prématuré atteint de pseudohypoaldostéronisme de type 1. Cette maladie autosomale récessive se caractérise, généralement, par une hyponatrémie liée à d'importantes pertes de sel dans les urines, une hyperkaliémie, ainsi qu'un niveau élevé d'aldostérone. Tout d'abord aucune des expériences biochimiques et électrophysiologiques n'a pu démontrer un défaut d'expression ou une forte diminution de l'activité soutenant les données cliniques. Cependant, en challengeant aS243PßyENaC avec une forte concentration de Na+ externe, une hypersensibilité de canal fut observée. En effet, ni les phénomènes régulateurs de « feedback inhibition » ou de « Na+ self-inhibition » n'étaient semblables au canal sauvage. De ce fait, ils apparaissaient exacerbés en présence de la mutation, amenant ainsi à une diminution de la réabsorption de Na+. Ceci corrobore entièrement l'hyponatrémie diagnostiquée. Le rein d'un prématuré étant immature, la quantité de Na+ atteignant la partie distale du néphron est plus élevée, du fait que les autres mécanismes de réabsorption en amont ne sont probablement pas encore en place. Cette hypothèse est renforcée par l'existence d'un frère présentant la même mutation, mais qui, né à terme, ne présentait aucun signe d'hyponatrémie. - The main topic of my thesis is the structure-function relationship of the ENaC/Deg family of ion channels, namely the Acid-Sensing Ion Channel ASIC1a and the Epithelial Na Channel ENaC. The primary part of this research is dedicated to the structure of ASIC1a. Four channel crystals have been published, which support a trimeric stoichiometry, although these data contradict previous functional experiments on other ENaC/Deg members. We are therefore interested in ASIC1a oligomerization and have set up a new assay combining the Substituted- Cysteine Accessibility Method (SCAM) with Afunctional sulfhydryl reagents (crosslinkers) allowing its study. The aim was to first stabilize the channels, then select those that are functional and then resolve them according to their size on SDS-PAGE. We demonstrated that the stabilized complex has a molecular weight corresponding to a tetrameric stoichiometry. In addition to our interest in the oligomerization of the ENaC/Deg family of ion channels, we also wanted to investigate the thus far undefined way of permeation for these channels. Therefore, taking the advantage of a more electrophysiological approach, we studied the accessibility of specific cysteines for methanethiosulfonate reagents (MTS) and were able to correlate the MTS association kinetics on cysteine residues with Na+ currents. These results have given us an insight into ion permeation and our functional evidence indicates that the extracellular is larger than that depicted by the crystal structures. As a side project, we focused on ENaC, which is made up of three subunits (a, ß and y) and is expressed in various epithelia, especially in the distal nephron of the kidneys. It plays a role in Na+ homeostasis and is essentially regulated by hormones via aldosterone and vasopressin, but also by serine proteases or Na+. We have studied the functional impact of the aS243P mutation, discovered in a premature baby suffering from pseudohypoaldosteronism of type 1. This autosomal recessive disease is characterized by hyponatremia, hyperkalemia and high aldosterone levels. Firstly, neither biochemical nor electrophysiological experiments indicated an expression defect or a strong decrease in activity. However, challenging aS243PßyENaC with increased external Na+ concentration showed channel hypersensitivity. Indeed, both the "feedback inhibition" and the "Na+ self-inhibition" regulatory mechanisms are impaired, leading to a decrease in Na+ reabsorption, entirely supports the diagnosis. The kidneys in preterm infants are immature and Na+ levels reaching the distal nephron are higher than normally observed. We hypothesize that the upstream reabsorption machinery is unlikely to be sufficiently matured and this assumption is supported by an asymptomatic sibling carrying the same mutation, but born at term. - La cellule, unité fonctionnelle du corps humain, est délimitée par une membrane plasmique servant de barrière biologique entre les milieux intra et extracellulaires. Une communication entre cellules est indispensable pour un fonctionnement adéquat. Sa survie dépend, entre autres, du maintien de la teneur en ions dans chacun des milieux qui doivent pouvoir être réabsorbés, ou sécrétés, selon les besoins. Les protéines insérées dans la membrane forment un canal et sont un moyen de communication permettant spécifiquement à des ions tel que le sodium (Na+) de traverser. Le Na+ se trouve dans la plupart des aliments et le sel, et est spécifiquement réabsorbé au niveau des reins grâce au canal sodique épithélial ENaC. Cette réabsorption se fait de l'urine primaire vers l'intérieur de la cellule, puis est transporté vers le sang. Pour maintenir un équilibre, une régulation de ce canal est nécessaire. En effet, des dysfonctionnements impliquant la régulation ou l'activité d'ENaC lui-même sont à l'origine de maladies telles que la mucoviscidose, l'hypertension ou encore, le pseudohypoaldostéronisme (PHA). Cette maladie est caractérisée, notamment, par d'importantes pertes de sel dans les urines. Des pédiatres ont diagnostiqué un PHA chez un nouveau-né, ce dernier présentant une modification du canal ENaC, nous avons recréé cette protéine afin d'étudier l'impact de ce changement sur son activité. Nous avons démontré que la régulation d'ENaC était effectivement perturbée, conduisant ainsi à une forte réduction de la réabsorption sodique. Afin de développer des molécules capables de moduler l'activité de protéines. Il est nécessaire d'en connaître la structure. Celle du canal sodique sensible à l'acidification ASIC1, un canal cousin d'ENaC, est connue. Ces données structurales contredisant cependant les analyses fonctionnelles, nous nous sommes penchés une nouvelle fois sur ASIC1. Une protéine est une macromolécule biologique composée d'une chaîne d'acides aminés (aa). De l'enchaînement d'aa à la protéine fonctionnelle, quatre niveaux de structuration existent. Chaque aa donne une indication quant au repliement et plus particulièrement la cystéine. Arborant un groupe sulfhydryle (SH) capable de former une liaison spécifique et stable avec un autre SH, celle-ci est souvent impliquée dans la structure tridimensionnelle de la protéine. Ce type de liaison intervient également dans la stabilisation de la structure quaternaire, qui est l'association de plusieurs protéines identiques (homomère), ou pas (hétéromère). Dans cette partie, nous avons remplacé des aa par des cystéines à des endroits spécifiques. Le but était de stabiliser plusieurs homomères d'ASICl ensemble avec des réactifs créant des ponts entre deux SH. Ainsi, nous avons pu déterminer le nombre de protéines ASIC1 participant à la formation d'un canal fonctionnel. Nos résultats corroborent les données fonctionnelles soutenant un canal tétramérique. Nous avons également étudié l'accessibilité de ces nouvelles cystéines afin d'obtenir des informations supplémentaires sur la structure du chemin emprunté par le Na+ à travers ASIC1 et plus particulièrement du vestibule extracellulaire.

Characterization and simulation of a CdTe detector for use in PET

The Voxel Imaging PET (VIP) Path nder project got the 4 year European Research Council FP7 grant in 2010 to prove the feasibility of using CdTe detectors in a novel conceptual design of PET scanner. The work presented in this thesis is a part of the VIP project and consists of, on the one hand, the characterization of a CdTe detector in terms of energy resolution and coincidence time resolution and, on the other hand, the simulation of the setup with the single detector in order to extend the results to the full PET scanner. An energy resolution of 0.98% at 511 keV with a bias voltage of 1000 V/mm has been measured at low temperature T=-8 ºC. The coincidence time distribution of two twin detectors has been found to be as low as 6 ns FWHM for events with energies above 500 keV under the same temperature and bias conditions. The measured energy and time resolution values are compatible with similar ndings available in the literature and prove the excellent potential of CdTe for PET applications. This results have been presented in form of a poster contribution at the IEEE NSS/MIC & RTSD 2011 conference in October 2011 in Valencia and at the iWoRID 2012 conference in July 2012 in Coimbra, Portugal. They have been also submitted for publication to "Journal of Instrumentation (JINST)" in September 2012.

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested

Improved earthquake response via simulation and integrated space- and ground-based technologies: the TREMOR proposal

Earthquakes occurring around the world each year cause thousands ofdeaths, millions of dollars in damage to infrastructure, and incalculablehuman suffering. In recent years, satellite technology has been asignificant boon to response efforts following an earthquake and itsafter-effects by providing mobile communications between response teamsand remote sensing of damaged areas to disaster management organizations.In 2007, an international team of students and professionals assembledduring theInternational Space University’s Summer Session Program in Beijing, Chinato examine how satellite and ground-based technology could be betterintegrated to provide an optimised response in the event of an earthquake.The resulting Technology Resources for Earthquake MOnitoring and Response(TREMOR) proposal describes an integrative prototype response system thatwill implement mobile satellite communication hubs providing telephone anddata links between response teams, onsite telemedicine consultation foremergency first-responders, and satellite navigation systems that willlocate and track emergency vehicles and guide search-and-rescue crews. Aprototype earthquake simulation system is also proposed, integratinghistorical data, earthquake precursor data, and local geomatics andinfrastructure information to predict the damage that could occur in theevent of an earthquake. The backbone of these proposals is a comprehensiveeducation and training program to help individuals, communities andgovernments prepare in advance. The TREMOR team recommends thecoordination of these efforts through a centralised, non-governmentalorganization.

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

Size-dependent predation by Dugesia lugubris (Turbellaria) on Physa acuta (Gastropoda): experiments and model

1. We investigated experimentally predation by the flatworm Dugesia lugubris on the snail Physa acuta in relation to predator body length and to prey morphology [shell length (SL) and aperture width (AW)]. 2. SL and AW correlate strongly in the field, but display significant and independent variance among populations. In the laboratory, predation by Dugesia resulted in large and significant selection differentials on both SL and AW. Analysis of partial effects suggests that selection on AW was indirect, and mediated through its strong correlation with SL. 3. The probability P(ij) for a snail of size category i (SL) to be preyed upon by a flatworm of size category j was fitted with a Poisson-probability distribution, the mean of which increased linearly with predator size (i). Despite the low number of parameters, the fit was excellent (r2 = 0.96). We offer brief biological interpretations of this relationship with reference to optimal foraging theory. 4. The largest size class of Dugesia (>2 cm) did not prey on snails larger than 7 mm shell length. This size threshold might offer Physa a refuge against flatworm predation and thereby allow coexistence in the field. 5. Our results are further discussed with respect to previous field and laboratory observations on P acuta life-history patterns, in particular its phenotypic variance in adult body size.

Population PK-guided simulation study in children exposed to drugs in human milk

Risks of significant infant drug exposure through human milk arepoorly defined due to lack of large-scale PK data. We propose to useBayesian approach based on population PK (popPK)-guided modelingand simulation for risk prediction. As a proof-of-principle study, weexploited fluoxetine milk concentration data from 25 women. popPKparameters including milk-to-plasma ratio (MP ratio) were estimatedfrom the best model. The dose of fluoxetine the breastfed infant wouldreceive through mother's milk, and infant plasma concentrations wereestimated from 1000 simulated mother-infant pairs, using randomassignment of feeding times and milk volume. A conservative estimateof CYP2D6 activity of 20% of the allometrically-adjusted adult valuewas assumed. Derived model parameters, including MP ratio were consistentwith those reported in the literature. Visual predictive check andother model diagnostics showed no signs of model misspecifications.The model simulation predicted that infant exposure levels to fluoxetinevia mother's milk were below 10% of weight-adjusted maternal therapeuticdoses in >99% of simulated infants. Predicted median ratio ofinfant-mother serum levels at steady state was 0.093 (range 0.033-0.31),consistent with literature reported values (mean=0.07; range 0-0.59).Predicted incidence of relatively high infant-mother ratio (>0.2) ofsteady-state serum fluoxetine concentrations was <1.3%. Overall, ourpredictions are consistent with clinical observations. Our approach maybe valid for other drugs, allowing in silico prediction of infant drugexposure risks through human milk. We will discuss application of thisapproach to another drug used in lactating women.

Prediction of infant drug exposure through breastfeeding: population PK modeling and simulation of fluoxetine

BACKGROUND: Risks of significant infant drug exposurethrough breastmilk are poorly defined for many drugs, and largescalepopulation data are lacking. We used population pharmacokinetics(PK) modeling to predict fluoxetine exposure levels ofinfants via mother's milk in a simulated population of 1000 motherinfantpairs.METHODS: Using our original data on fluoxetine PK of 25breastfeeding women, a population PK model was developed withNONMEM and parameters, including milk concentrations, wereestimated. An exponential distribution model was used to account forindividual variation. Simulation random and distribution-constrainedassignment of doses, dosing time, feeding intervals and milk volumewas conducted to generate 1000 mother-infant pairs with characteristicssuch as the steady-state serum concentrations (Css) and infantdose relative to the maternal weight-adjusted dose (relative infantdose: RID). Full bioavailability and a conservative point estimate of1-month-old infant CYP2D6 activity to be 20% of the adult value(adjusted by weigth) according to a recent study, were assumed forinfant Css calculations.RESULTS: A linear 2-compartment model was selected as thebest model. Derived parameters, including milk-to-plasma ratios(mean: 0.66; SD: 0.34; range, 0 - 1.1) were consistent with the valuesreported in the literature. The estimated RID was below 10% in >95%of infants. The model predicted median infant-mother Css ratio was0.096 (range 0.035 - 0.25); literature reported mean was 0.07 (range0-0.59). Moreover, the predicted incidence of infant-mother Css ratioof >0.2 was less than 1%.CONCLUSION: Our in silico model prediction is consistent withclinical observations, suggesting that substantial systemic fluoxetineexposure in infants through human milk is rare, but further analysisshould include active metabolites. Our approach may be valid forother drugs. [supported by CIHR and Swiss National Science Foundation(SNSF)]