Theoretical considerations to the verification of dynamic multileaf collimated fields with a SLIC-type portal image detector

The verification possibilities of dynamically collimated treatment beams with a scanning liquid ionization chamber electronic portal image device (SLIC-EPID) are investigated. The ion concentration in the liquid of a SLIC-EPID and therefore the read-out signal is determined by two parameters of a differential equation describing the creation and recombination of the ions. Due to the form of this equation, the portal image detector describes a nonlinear dynamic system with memory. In this work, the parameters of the differential equation were experimentally determined for the particular chamber in use and for an incident open 6 MV photon beam. The mathematical description of the ion concentration was then used to predict portal images of intensity-modulated photon beams produced by a dynamic delivery technique, the sliding window approach. Due to the nature of the differential equation, a mathematical condition for 'reliable leaf motion verification' in the sliding window technique can be formulated. It is shown that the time constants for both formation and decay of the equilibrium concentration in the chamber is in the order of seconds. In order to guarantee reliable leaf motion verification, these time constants impose a constraint on the rapidity of the image-read out for a given maximum leaf speed. For a leaf speed of 2 cm s(-1), a minimum image acquisition frequency of about 2 Hz is required. Current SLIC-EPID systems are usually too slow since they need about a second to acquire a portal image. However, if the condition is fulfilled, the memory property of the system can be used to reconstruct the leaf motion. It is shown that a simple edge detecting algorithm can be employed to determine the leaf positions. The method is also very robust against image noise.

Next-to-Next-to-Next-to-Leading Order QCD Prediction for the Top Antitop S-Wave Pair Production Cross Section Near Threshold in e + e − Annihilation

We present the third-order QCD prediction for the production of top antitop quark pairs in electron-positron collisions close to the threshold in the dominant S-wave state. We observe a significant reduction of the theoretical uncertainty and discuss the sensitivity to the top quark mass and width.

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

Dual-energy CT-cholangiography in potential donors for living-related liver transplantation: Improved biliary visualization by intravenous morphine co-medication

PURPOSE: To prospectively evaluate whether intravenous morphine co-medication improves bile duct visualization of dual-energy CT-cholangiography. MATERIALS AND METHODS: Forty potential donors for living-related liver transplantation underwent CT-cholangiography with infusion of a hepatobiliary contrast agent over 40min. Twenty minutes after the beginning of the contrast agent infusion, either normal saline (n=20 patients; control group [CG]) or morphine sulfate (n=20 patients; morphine group [MG]) was injected. Forty-five minutes after initiation of the contrast agent, a dual-energy CT acquisition of the liver was performed. Applying dual-energy post-processing, pure iodine images were generated. Primary study goals were determination of bile duct diameters and visualization scores (on a scale of 0 to 3: 0-not visualized; 3-excellent visualization). RESULTS: Bile duct visualization scores for second-order and third-order branch ducts were significantly higher in the MG compared to the CG (2.9±0.1 versus 2.6±0.2 [P<0.001] and 2.7±0.3 versus 2.1±0.6 [P<0.01], respectively). Bile duct diameters for the common duct and main ducts were significantly higher in the MG compared to the CG (5.9±1.3mm versus 4.9±1.3mm [P<0.05] and 3.7±1.3mm versus 2.6±0.5mm [P<0.01], respectively). CONCLUSION: Intravenous morphine co-medication significantly improved biliary visualization on dual-energy CT-cholangiography in potential donors for living-related liver transplantation.

Intraoperative three-dimensional fluoroscopic cholangiography

Precise intraoperative assessment of the architecture of the biliary tree could reduce lesions to intra- or extrahepatic bile ducts. The aim of this study was to test feasibility of intraoperative three-dimensional imaging during liver resections. Isocentric C-arm fluoroscopy acquires three-dimensional images during a 190 degrees orbital rotation. The bile ducts were displayed three-dimensionally by realtime rotational projections or multiplanar reconstructions. The technique was established ex vivo in a preserved cadaveric human liver. Intraoperative three-dimensional cholangiography was performed in five patients with centrally located liver malignancies. Complete data acquisition in 3 patients depicted precise anatomical details of the architecture of the biliary tree up to third order divisions. Biliary imaging can be improved by the application of real-time intraoperative three-dimensional cholangiography. For the development of computer-aided navigation in hepatobiliary procedures, this technique could be an important prerequisite for defining landmarks of the liver in a three-dimensional space.

High-resolution sequence stratigraphic correlation in the Upper Jurassic (Kimmeridgian)–Upper Cretaceous (Cenomanian) peritidal carbonate deposits (Western Taurides, Turkey)

Upper Jurassic (Kimmeridgian)±Upper Cretaceous (Cenomanian) inner platform carbonates in the Western Taurides are composed of metre-scale upward-shallowing cyclic deposits (parasequences) and important karstic surfaces capping some of the cycles. Peritidal cycles (shallow subtidal facies capped by tidal-¯at laminites or fenestrate limestones) are regressive- and transgressive-prone (upward-deepening followed by upward-shallowing facies trends). Subtidal cycles are of two types and indicate incomplete shallowing. Submerged subtidal cycles are composed of deeper subtidal facies overlain by shallow subtidal facies. Exposed subtidal cycles consist of deeper subtidal facies overlain by shallow subtidal facies that are capped by features indicative of prolonged subaerial exposure. Subtidal facies occur characteristically in the Jurassic, while peritidal cycles are typical for the Lower Cretaceous of the region. Within the foraminiferal and dasyclad algal biostratigraphic framework, four karst breccia levels are recognized as the boundaries of major second-order cycles, introduced for the ®rst time in this study. These levels correspond to the Kimmeridgian±Portlandian boundary, mid-Early Valanginian, mid-Early Aptian and mid-Cenomanian and represent important sea level falls which affected the distribution of foraminiferal fauna and dasyclad ¯ora of the Taurus carbonate platform. Within the Kimmeridgian±Cenomanian interval 26 third-order sequences (types 1 and 2) are recognized. These sequences are the records of eustatic sea level ¯uctuations rather than the records of local tectonic events because the boundaries of the sequences representing 1±4 Ma intervals are correlative with global sea level falls. Third-order sequences and metre-scale cyclic deposits are the major units used for long-distance, high-resolution sequence stratigraphic correlation in the Western Taurides. Metre-scale cyclic deposits (parasequences) in the Cretaceous show genetical stacking patterns within third-order sequences and correspond to fourth-order sequences representing 100±200 ka. These cycles are possibly the E2 signal (126 ka) of the orbital eccentricity cycles of the Milankovitch band. The slight deviation of values, calculated for parasequences, from the mean value of eccentricity cycles can be explained by the currently imprecise geochronology established in the Cretaceous and missed sea level oscillations when the platform lay above fluctuating sea level.

Comparative theoretical analysis of continuous wave laser cutting of metals at 1 and 10 um wavelength

We present a derivation and, based on it, an extension of a model originally proposed by V.G. Niziev to describe continuous wave laser cutting of metals. Starting from a local energy balance and by incorporating heat removal through heat conduction to the bulk material, we find a differential equation for the cutting profile. This equation is solved numerically and yields, besides the cutting profiles, the maximum cutting speed, the absorptivity profiles, and other relevant quantities. Our main goal is to demonstrate the model’s capability to explain some of the experimentally observed differences between laser cutting at around 1 and 10 μm wavelengths. To compare our numerical results to experimental observations, we perform simulations for exactly the same material and laser beam parameters as those used in a recent comparative experimental study. Generally, we find good agreement between theoretical and experimental results and show that the main differences between laser cutting with 1- and 10-μm beams arise from the different absorptivity profiles and absorbed intensities. Especially the latter suggests that the energy transfer, and thus the laser cutting process, is more efficient in the case of laser cutting with 1-μm beams.

Demographic model of the Swiss cattle population for the years 2009-2011 stratified by gender, age and production type

Demographic composition and dynamics of animal and human populations are important determinants for the transmission dynamics of infectious disease and for the effect of infectious disease or environmental disasters on productivity. In many circumstances, demographic data are not available or of poor quality. Since 1999 Switzerland has been recording cattle movements, births, deaths and slaughter in an animal movement database (AMD). The data present in the AMD offers the opportunity for analysing and understanding the dynamic of the Swiss cattle population. A dynamic population model can serve as a building block for future disease transmission models and help policy makers in developing strategies regarding animal health, animal welfare, livestock management and productivity. The Swiss cattle population was therefore modelled using a system of ordinary differential equations. The model was stratified by production type (dairy or beef), age and gender (male and female calves: 0-1 year, heifers and young bulls: 1-2 years, cows and bulls: older than 2 years). The simulation of the Swiss cattle population reflects the observed pattern accurately. Parameters were optimized on the basis of the goodness-of-fit (using the Powell algorithm). The fitted rates were compared with calculated rates from the AMD and differed only marginally. This gives confidence in the fitted rates of parameters that are not directly deductible from the AMD (e.g. the proportion of calves that are moved from the dairy system to fattening plants).

Higgs effects in top anti-top production near threshold in e+e− annihilation

The completion of the third-order QCD corrections to the inclusive top-pair production cross section near threshold demonstrates that the strong dynamics is under control at the few percent level. In this paper we consider the effects of the Higgs boson on the cross section and, for the first time, combine the third-order QCD result with the third-order P-wave, the leading QED and the leading non-resonant contributions. We study the size of the different effects and investigate the sensitivity of the cross section to variations of the top-quark Yukawa coupling due to possible new physics effects.

Optimized CUDA-Based PDE Solver for Reaction Diffusion Systems on Arbitrary Surfaces

Partial differential equation (PDE) solvers are commonly employed to study and characterize the parameter space for reaction-diffusion (RD) systems while investigating biological pattern formation. Increasingly, biologists wish to perform such studies with arbitrary surfaces representing ‘real’ 3D geometries for better insights. In this paper, we present a highly optimized CUDA-based solver for RD equations on triangulated meshes in 3D. We demonstrate our solver using a chemotactic model that can be used to study snakeskin pigmentation, for example. We employ a finite element based approach to perform explicit Euler time integrations. We compare our approach to a naive GPU implementation and provide an in-depth performance analysis, demonstrating the significant speedup afforded by our optimizations. The optimization strategies that we exploit could be generalized to other mesh based processing applications with PDE simulations.