Monte Carlo simulation of a dynamic MLC based on a multiple source model

Detailed knowledge of the characteristics of the radiation field shaped by a multileaf collimator (MLC) is essential in intensity modulated radiotherapy (IMRT). A previously developed multiple source model (MSM) for a 6 MV beam was extended to a 15 MV beam and supplemented with an accurate model of an 80-leaf dynamic MLC. Using the supplemented MSM and the MC code GEANT, lateral dose distributions were calculated in a water phantom and a portal water phantom. A field which is normally used for the validation of the step and shoot technique and a field from a realistic IMRT treatment plan delivered with dynamic MLC are investigated. To assess possible spectral changes caused by the modulation of beam intensity by an MLC, the energy spectra in five portal planes were calculated for moving slits of different widths. The extension of the MSM to 15 MV was validated by analysing energy fluences, depth doses and dose profiles. In addition, the MC-calculated primary energy spectrum was verified with an energy spectrum which was reconstructed from transmission measurements. MC-calculated dose profiles using the MSM for the step and shoot case and for the dynamic MLC case are in very good agreement with the measured data from film dosimetry. The investigation of a 13 cm wide field shows an increase in mean photon energy of up to 16% for the 0.25 cm slit compared to the open beam for 6 MV and of up to 6% for 15 MV, respectively. In conclusion, the MSM supplemented with the dynamic MLC has proven to be a powerful tool for investigational and benchmarking purposes or even for dose calculations in IMRT.

Simulation in training: one-year experience using an efficiency index to assess interventional radiology fellow training status

PURPOSE: To compare objective fellow and expert efficiency indices for an interventional radiology renal artery stenosis skill set with the use of a high-fidelity simulator. MATERIALS AND METHODS: The Mentice VIST simulator was used for three different renal artery stenosis simulations of varying difficulty, which were used to grade performance. Fellows' indices at three intervals throughout 1 year were compared to expert baseline performance. Seventy-four simulated procedures were performed, 63 of which were captured as audiovisual recordings. Three levels of fellow experience were analyzed: 1, 6, and 12 months of dedicated interventional radiology fellowship. The recordings were compiled on a computer workstation and analyzed. Distinct measurable events in the procedures were identified with task analysis, and data regarding efficiency were extracted. Total scores were calculated as the product of procedure time, fluoroscopy time, tools, and contrast agent volume. The lowest scores, which reflected efficient use of tools, radiation, and time, were considered to indicate proficiency. Subjective analysis of participants' procedural errors was not included in this analysis. RESULTS: Fellows' mean scores diminished from 1 month to 12 months (42,960 at 1 month, 18,726 at 6 months, and 9,636 at 12 months). The experts' mean score was 4,660. In addition, the range of variance in score diminished with increasing experience (from a range of 5,940-120,156 at 1 month to 2,436-85,272 at 6 months and 2,160-32,400 at 12 months). Expert scores ranged from 1,450 to 10,800. CONCLUSIONS: Objective efficiency indices for simulated procedures can demonstrate scores directly comparable to the level of clinical experience.

High-resolution computer simulation of electrophoretic mobilization in isoelectric focusing

Cationic and anionic electrophoretic mobilization for focusing of hemoglobins (Hb's) in the presence of 100 carrier ampholytes covering a pI range of 6.00-7.98 was studied by computer simulation at a constant current density of 300 A/m(2). Electropherograms that would be produced by whole column imaging and by single detectors placed at different locations along the focusing column are presented. Upon mobilization, peak heights of the Hb zones decrease, but the zones retain a relatively sharp constant profile and are migrating at a constant velocity. A further peak decrease occurs during readjustment at the locations of the original buffer/column interfaces, indicating that detection sensitivity is the lowest at these locations. An anionic carrier ampholyte mobility smaller than that of its cationic species produces a cathodic drift which is smaller than the transport rate used for electrophoretic mobilization. Compared to the case with equal mobilities of carrier ampholyte species, a small increase (decrease) is predicted for the cationic (anionic) mobilization rate within the focusing column. Simulation data suggest that electrophoretic mobilization after focusing and focusing with concurrent electrophoretic mobilization are comparable isotachophoretic processes that occur when there is an uninterrupted flux of an ion through the focusing column. Cathodic drift caused by unequal mobilities of the species of carrier ampholytes, electrophoretic mobilization, and decomposition occurring at the pH gradient edges are related electrophoretic processes.

Range of motion after computed tomography-based simulation of intertrochanteric corrective osteotomy in cases of slipped capital femoral epiphysis: comparison of uniplanar flexion osteotomy and multiplanar flexion, valgisation, and rotational osteotomies

BACKGROUND: Various osteotomy techniques have been developed to correct the deformity caused by slipped capital femoral epiphysis (SCFE) and compared by their clinical outcomes. The aim of the presented study was to compare an intertrochanteric uniplanar flexion osteotomy with a multiplanar osteotomy by their ability to improve postoperative range of motion as measured by simulation of computed tomographic data in patients with SCFE. METHODS: We examined 19 patients with moderate or severe SCFE as classified based on slippage angle. A computer program for the simulation of movement and osteotomy developed in our laboratory was used for study execution. According to a 3-dimensional reconstruction of the computed tomographic data, the physiological range was determined by flexion, abduction, and internal rotation. The multiplanar osteotomy was compared with the uniplanar flexion osteotomy. Both intertrochanteric osteotomy techniques were simulated, and the improvements of the movement range were assessed and compared. RESULTS: The mean slipping and thus correction angles measured were 25 degrees (range, 8-46 degrees) inferior and 54 degrees (range, 32-78 degrees) posterior. After the simulation of multiplanar osteotomy, the virtually measured ranges of motion as determined by bone-to-bone contact were 61 degrees for flexion, 57 degrees for abduction, and 66 degrees for internal rotation. The simulation of the uniplanar flexion osteotomy achieved a flexion of 63 degrees, an abduction of 36 degrees, and an internal rotation of 54 degrees. CONCLUSIONS: Apart from abduction, the improvement in the range of motion by a uniplanar flexion osteotomy is comparable with that of the multiplanar osteotomy. However, the improvement in flexion for the simulation of both techniques is not satisfactory with regard to the requirements of normal everyday life, in contrast to abduction and internal rotation. LEVEL OF EVIDENCE: Level III, Retrospective comparative study.

A 700-year paleoecological record of boreal ecosystem responses to climatic variation from Alaska

Recent observations and model simulations have highlighted the sensitivity of the forest - tundra ecotone to climatic forcing. In contrast, paleoecological studies have not provided evidence of tree-line fluctuations in response to Holocene climatic changes in Alaska, suggesting that the forest - tundra boundary in certain areas may be relatively stable at multicentennial to millennial time scales. We conducted a multiproxy study of sediment cores from an Alaskan lake near the altitudinal limits of key boreal-forest species. Paleoecological data were compared with independent climatic reconstructions to assess ecosystem responses of the forest - tundra boundary to Little Ice Age (LIA) climatic. uctuations. Pollen, diatom, charcoal, macrofossil, and magnetic analyses provide the first continuous record of vegetation -. re - climate interactions at decadal to centennial time scales during the past 700 years from southern Alaska. Boreal-forest diebacks characterized by declines of Picea mariana, P. glauca, and tree Betula occurred during the LIA ( AD 1500 - 1800), whereas shrubs ( Alnus viridis, Betula glandulosa/nana) and herbaceous taxa (Epilobium, Aconitum) expanded. Marked increases in charcoal abundance and changes in magnetic properties suggest increases in. re importance and soil erosion during the same period. In addition, the conspicuous reduction or disappearance of certain aquatic ( e. g., Isoetes, Nuphar, Pediastrum) and wetland ( Sphagnum) plants and major shifts in diatom assemblages suggest pronounced lake-level. uctuations and rapid ecosystem reorganization in response to LIA climatic deterioration. Our results imply that temperature shifts of 1 - 2 degrees C, when accompanied by major changes in moisture balance, can greatly alter high-altitudinal terrestrial, wetland, and aquatic ecosystems, including conversion between boreal-forest tree line and tundra. The climatic and ecosystem variations in our study area appear to be coherent with changes in solar irradiance, suggesting that changes in solar activity contributed to the environmental instability of the past 700 years.

Evaluation of a Bayesian MCMC Random Effects Inference Methodology for Capture-Mark-Recapture Data

Monte Carlo simulation was used to evaluate properties of a simple Bayesian MCMC analysis of the random effects model for single group Cormack-Jolly-Seber capture-recapture data. The MCMC method is applied to the model via a logit link, so parameters p, S are on a logit scale, where logit(S) is assumed to have, and is generated from, a normal distribution with mean μ and variance σ2 . Marginal prior distributions on logit(p) and μ were independent normal with mean zero and standard deviation 1.75 for logit(p) and 100 for μ ; hence minimally informative. Marginal prior distribution on σ2 was placed on τ2=1/σ2 as a gamma distribution with α=β=0.001 . The study design has 432 points spread over 5 factors: occasions (t) , new releases per occasion (u), p, μ , and σ . At each design point 100 independent trials were completed (hence 43,200 trials in total), each with sample size n=10,000 from the parameter posterior distribution. At 128 of these design points comparisons are made to previously reported results from a method of moments procedure. We looked at properties of point and interval inference on μ , and σ based on the posterior mean, median, and mode and equal-tailed 95% credibility interval. Bayesian inference did very well for the parameter μ , but under the conditions used here, MCMC inference performance for σ was mixed: poor for sparse data (i.e., only 7 occasions) or σ=0 , but good when there were sufficient data and not small σ .

A real time simulation model of glucose-insulin metabolism for type 1 diabetes patients

In this paper, a simulation model of glucose-insulin metabolism for Type 1 diabetes patients is presented. The proposed system is based on the combination of Compartmental Models (CMs) and artificial Neural Networks (NNs). This model aims at the development of an accurate system, in order to assist Type 1 diabetes patients to handle their blood glucose profile and recognize dangerous metabolic states. Data from a Type 1 diabetes patient, stored in a database, have been used as input to the hybrid system. The data contain information about measured blood glucose levels, insulin intake, and description of food intake, along with the corresponding time. The data are passed to three separate CMs, which produce estimations about (i) the effect of Short Acting (SA) insulin intake on blood insulin concentration, (ii) the effect of Intermediate Acting (IA) insulin intake on blood insulin concentration, and (iii) the effect of carbohydrate intake on blood glucose absorption from the gut. The outputs of the three CMs are passed to a Recurrent NN (RNN) in order to predict subsequent blood glucose levels. The RNN is trained with the Real Time Recurrent Learning (RTRL) algorithm. The resulted blood glucose predictions are promising for the use of the proposed model for blood glucose level estimation for Type 1 diabetes patients.

Palaeohydrology, fires and vegetation succession in the southern Baltic during the last 7500 years reconstructed from a raised bog based on multi-proxy data

We present the first 7500 yr long multi-proxy record from a raised bog located at the southern Baltic coast, Poland. Testate amoebae, plant macrofossils, pollen and microscopic charcoal were used to reconstruct environmental changes in Pomerania (northern Poland, Kaszuby Lakeland) from a 7-m thick peat archive of Stążki bog dated 5500 BC–AD 1250. We obtained a record of proxies representing different spatial scales: regional vegetation changed simultaneously with local vegetation, and testate amoebae showed a pattern of change similar to that of pollen and plant macrofossils. On the basis of the combined proxies, we distinguished three hydroclimatic stages: moist conditions 5500–3450 BC, drier conditions with regionally increased fires up to 600 BC, and again moist conditions from 600 BC onward. During the drier interval, a first climatic shift to wetter conditions at 1700 BC is indicated by regional pollen as the replacement of Corylus by Carpinus, and locally by, e.g., the increase of Hyalosphenia elegans and mire plants such as Sphagnum sec. Cuspidata. Furthermore, we observed a correlation since 600 BC among the re-expansion of Carpinus (after a sudden decline ca. 950 BC), increased peat accumulation, increase of Hyalosphenia species, and fewer fires, suggesting lower evapotranspiration and a stable high water table in the bog. Fagus started to expand after AD 810 gradually replacing Carpinus, which was possibly due to a gradually more oceanic climate, though we cannot exclude human impact on the forests. Peat accumulation, determined by radiocarbon dating, varied with bog surface wetness. The hydroclimatic phases found in Stążki peatland are similar to moisture changes recorded in other sites from Poland and Europe. This is the first detailed record of hydroclimatic change during the Holocene in the southern Baltic region, so it forms a reference site for further studies on other southern Baltic bogs that are in progress.

Climatic and human impacts on mountain vegetation at Lauenensee (Bernese Alps, Switzerland) during the last 14,000 years

Lake sediments from Lauenensee (1381 m a.s.l.), a small lake in the Bernese Alps, were analysed to reconstruct the vegetation and fire history. The chronology is based on 11 calibrated radiocarbon dates on terrestrial plant macrofossils suggesting a basal age of 14,200 cal. BP. Pollen and macrofossil data imply that treeline never reached the lake catchment during the Bølling–Allerød interstadial. Treeline north of the Alps was depressed by c. 300 altitudinal meters, if compared with southern locations. We attribute this difference to colder temperatures and to unbuffered cold air excursions from the ice masses in northern Europe. Afforestation started after the Younger Dryas at 11,600 cal. BP. Early-Holocene tree-Betula and Pinus sylvestris forests were replaced by Abies alba forests around 7500 cal. BP. Continuous high-resolution pollen and macrofossil series allow quantitative assessments of vegetation dynamics at 5900–5200 cal. BP (first expansion of Picea abies, decline of Abies alba) and 4100–2900 cal. BP (first collapse of Abies alba). The first signs of human activity became noticeable during the late Neolithic c. 5700–5200 cal. BP. Cross-correlation analysis shows that the expansion of Alnus viridis and the replacement of Abies alba by Picea abies after c. 5500 cal. BP was most likely a consequence of human disturbance. Abies alba responded very sensitively to a combination of fire and grazing disturbance. Our results imply that the current dominance of Picea abies in the upper montane and subalpine belts is a consequence of anthropogenic activities through the millennia.

Investigating the consistency between proxy-based reconstructions and climate models using data assimilation: a mid-Holocene case study

The mid-Holocene (6 kyr BP; thousand years before present) is a key period to study the consistency between model results and proxy-based reconstruction data as it corresponds to a standard test for models and a reasonable number of proxy-based records is available. Taking advantage of this relatively large amount of information, we have compared a compilation of 50 air and sea surface temperature reconstructions with the results of three simulations performed with general circulation models and one carried out with LOVECLIM, a model of intermediate complexity. The conclusions derived from this analysis confirm that models and data agree on the large-scale spatial pattern but the models underestimate the magnitude of some observed changes and that large discrepancies are observed at the local scale. To further investigate the origin of those inconsistencies, we have constrained LOVECLIM to follow the signal recorded by the proxies selected in the compilation using a data-assimilation method based on a particle filter. In one simulation, all the 50 proxy-based records are used while in the other two only the continental or oceanic proxy-based records constrain the model results. As expected, data assimilation leads to improving the consistency between model results and the reconstructions. In particular, this is achieved in a robust way in all the experiments through a strengthening of the westerlies at midlatitude that warms up northern Europe. Furthermore, the comparison of the LOVECLIM simulations with and without data assimilation has also objectively identified 16 proxy-based paleoclimate records whose reconstructed signal is either incompatible with the signal recorded by some other proxy-based records or with model physics.

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P.a.axillaris, P.a.parodii and P.a.subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P.axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100000years, which is compatible with a divergence time between 35000 and 107000years ago between P.a.axillaris and P.a.parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P.a.axillaris and P.a.parodii in their contact zone.

Relaxation of species-specific neighborhood effects in Bornean rain forest under climatic perturbation

Evidence of negative conspecific density dependence (NDD) operating on seedling survival and sapling recruitment has accumulated recently. In contrast, evidence of NDD operating on growth of trees has been circumstantial at best. Whether or not local NDD at the level of individual trees leads to NDD at the level of the community is still an open question. Moreover, whether and how perturbations interfere with these processes have rarely been investigated. We applied neighborhood models to permanent plot data from a Bornean dipterocarp forest censused over two 10-11 year periods. Although the first period was only lightly perturbed, a moderately strong El Nino event causing severe drought occurred in the first half of the second period. Such events are an important component of the environmental stochasticity affecting the region. We show that local NDD on growth of small-to-medium-sized trees may indeed translate to NDD at the level of the community. This interpretation is based on increasingly negative effects of bigger conspecific neighbors on absolute growth rates of individual trees with increasing basal area across the 18 most abundant overstory species in the first period. However, this relationship was much weaker in the second period. We interpreted this relaxation of local and community-level NDD as a consequence of increased light levels at the forest floor due to temporary leaf and twig loss of large trees in response to the drought event. Mitigation of NDD under climatic perturbation acts to decrease species richness, especially in forest overstory and therefore has an important role in determining species relative abundances at the site.

Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

Flow disturbances in stent-related coronary evaginations: a computational fluid-dynamic simulation study

Aims: Angiographic ectasias and aneurysms in stented segments have been associated with late stent thrombosis. Using optical coherence tomography (OCT), some stented segments show coronary evaginations reminiscent of ectasias. The purpose of this study was to explore, using computational fluid-dynamic (CFD) simulations, whether OCT-detected coronary evaginations can induce local changes in blood flow. Methods and results: OCT-detected evaginations are defined as outward bulges in the luminal vessel contour between struts, with the depth of the bulge exceeding the actual strut thickness. Evaginations can be characterised cross ectionally by depth and along the stented segment by total length. Assuming an ellipsoid shape, we modelled 3-D evaginations with different sizes by varying the depth from 0.2-1.0 mm, and the length from 1-9 mm. For the flow simulation we used average flow velocity data from non-diseased coronary arteries. The change in flow with varying evagination sizes was assessed using a particle tracing test where the particle transit time within the segment with evagination was compared with that of a control vessel. The presence of the evagination caused a delayed particle transit time which increased with the evagination size. The change in flow consisted locally of recirculation within the evagination, as well as flow deceleration due to a larger lumen - seen as a deflection of flow towards the evagination. Conclusions: CFD simulation of 3-D evaginations and blood flow suggests that evaginations affect flow locally, with a flow disturbance that increases with increasing evagination size.