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.

Dynamic high-resolution computer simulation of isotachophoretic enantiomer separation and zone stability

The development of electrophoretic computer models and their use for simulation of electrophoretic processes has increased significantly during the last few years. Recently, GENTRANS and SIMUL5 were extended with algorithms that describe chemical equilibria between solutes and a buffer additive in a fast 1:1 interaction process, an approach that enables simulation of the electrophoretic separation of enantiomers. For acidic cationic systems with sodium and H3 0(+) as leading and terminating components, respectively, acetic acid as counter component, charged weak bases as samples, and a neutral CD as chiral selector, the new codes were used to investigate the dynamics of isotachophoretic adjustment of enantiomers, enantiomer separation, boundaries between enantiomers and between an enantiomer and a buffer constituent of like charge, and zone stability. The impact of leader pH, selector concentration, free mobility of the weak base, mobilities of the formed complexes and complexation constants could thereby be elucidated. For selected examples with methadone enantiomers as analytes and (2-hydroxypropyl)-β-CD as selector, simulated zone patterns were found to compare well with those monitored experimentally in capillary setups with two conductivity detectors or an absorbance and a conductivity detector. Simulation represents an elegant way to provide insight into the formation of isotachophoretic boundaries and zone stability in presence of complexation equilibria in a hitherto inaccessible way.

Upconverter materials and upconversion solar-cell devices: simulation and characterization with broad solar spectrum illumination

Upconverter materials and upconverter solar devices were recently investigated with broad-band excitation revealing the great potential of upconversion to enhance the efficiency of solar cell at comparatively low solar concentration factors. In this work first attempts are made to simulate the behavior of the upconverter β-NaYF4 doped with Er3+ under broad-band excitation. An existing model was adapted to account for the lower absorption of broader excitation spectra. While the same trends as observed for the experiments were found in the simulation, the absolute values are fairly different. This makes an upconversion model that specifically considers the line shape function of the ground state absorption indispensable to achieve accurate simulations of upconverter materials and upconverter solar cell devices with broadband excitations, such as the solar radiation.

Comparison of implant stability by means of resonance frequency analysis for flapless and conventionally inserted implants

Background: Resonance frequency analysis (RFA) is a noninvasive technique for the quantitative assessment of implant stability. Information on the implant stability quotient (ISQ) of transmucosally inserted implants is limited. Purpose: The aim of this investigation was to compare the ISQ of conventionally inserted implants by raising a muco-periostal flap with implants inserted using a flapless procedure. Materials and Methods: Forty elderly patients with complete edentulous maxilla were consecutively admitted for treatment with implant-supported prostheses. A computer tomography was obtained for the computer-assisted implant planning. One hundred ten implants were placed conventionally in 23 patients (flap-group) and 85 implants in 17 patients by means of the flapless method (flapless-group) using a stereolithographic template. RFA measurements were performed after implant placement (baseline) and after a healing time of 12 weeks (reentry). Results: All implants exhibited clinically and radiographically successful osseointegration. Bone level did not change significantly neither for genders nor type of surgical protocol. Mean ISQ values of the flapless-group were significantly higher at baseline (p < .001) and at reentry (p < .001) compared with the flap-group. The ISQ values were significantly lower at reentry compared with baseline for the flap-group (p = .028) but not for the flapless-group. This group showed a moderate, but insignificant increase. RFA measurements of males resulted in ISQ values that were thoroughly higher as compared with females at both time-points in both groups. Correlation between RFA and bone level was not found. Conclusions: The flapless procedure showed favorable conditions with regard to implant stability and crestal bone level. Some changes of the ISQ values that represent primary (mechanical) and secondary (bone remodeling) implant stability were observed in slight favor of the flapless method and male patients. In properly planned and well-selected cases, the minimal invasive transmucosal technique using a drill-guide is a safe procedure.

Surveillance and simulation of bovine spongiform encephalopathy and scrapie in small ruminants in Switzerland

BACKGROUND: After bovine spongiform encephalopathy (BSE) emerged in European cattle livestock in 1986 a fundamental question was whether the agent established also in the small ruminants' population. In Switzerland transmissible spongiform encephalopathies (TSEs) in small ruminants have been monitored since 1990. While in the most recent TSE cases a BSE infection could be excluded, for historical cases techniques to discriminate scrapie from BSE had not been available at the time of diagnosis and thus their status remained unclear. We herein applied state-of-the-art techniques to retrospectively classify these animals and to re-analyze the affected flocks for secondary cases. These results were the basis for models, simulating the course of TSEs over a period of 70 years. The aim was to come to a statistically based overall assessment of the TSE situation in the domestic small ruminant population in Switzerland. RESULTS: In sum 16 TSE cases were identified in small ruminants in Switzerland since 1981, of which eight were atypical and six were classical scrapie. In two animals retrospective analysis did not allow any further classification due to the lack of appropriate tissue samples. We found no evidence for an infection with the BSE agent in the cases under investigation. In none of the affected flocks, secondary cases were identified. A Bayesian prevalence calculation resulted in most likely estimates of one case of BSE, five cases of classical scrapie and 21 cases of atypical scrapie per 100'000 small ruminants. According to our models none of the TSEs is considered to cause a broader epidemic in Switzerland. In a closed population, they are rather expected to fade out in the next decades or, in case of a sporadic origin, may remain at a very low level. CONCLUSIONS: In summary, these data indicate that despite a significant epidemic of BSE in cattle, there is no evidence that BSE established in the small ruminant population in Switzerland. Classical and atypical scrapie both occur at a very low level and are not expected to escalate into an epidemic. In this situation the extent of TSE surveillance in small ruminants requires reevaluation based on cost-benefit analysis.

Stability of the hard and soft tissue profile after mandibular advancement in sagittal split osteotomies: a longitudinal and long-term follow-up study

The aim of the study was to conduct a long-term follow-up investigation of the stability of hard and soft tissues after bilateral sagittal split osteotomy (BSSO) with rigid internal (RIF) fixation to advance the mandible. Sixteen consecutive patients (12 females and 4 males, mean age 21.4 years) were available for re-examination 12.7 years (T5) after surgery. The preceding follow-ups were before (T1), and 5 days (T2), 7.3 months (T3), and 13.9 months (T4) after surgery. Lateral cephalograms were traced by hand, digitized, and evaluated with the Dentofacial Planner program. The x-axis for the system of co-ordinates ran through sella (point zero) and the line NSL -7 degrees. Thus, the program determined the x- and y-values of each variable and the usual angles and distances. Statistical analysis was carried out using Wilcoxon's matched-pair signed-ranks test with Bonferroni adjustments. The relationships between the examined variables were analysed by Spearman rank correlation coefficients. The backward relapse at point B (T5) was 2.42 mm, or 50 per cent, and at pogonion 3.21 mm, or 60 per cent of the initial advancement. The mean net effect at T5 on the labial fold (soft tissue point B) was 94 per cent of the advancement at point B. For the soft tissue chin (soft tissue pogonion), it was 119 per cent of the advancement at pogonion. The net effect on the lower lip (labrale inferior) was 55 per cent of the advancement at incision inferior. The amount of the surgical advancement of the mandible was correlated with the long-term relapse in point B. Among possible reasons for this relapse are the initial soft tissue profile, the initial growth direction, and the remodelling processes of the hard tissue.

Altered perception of apparent motion in schizophrenia spectrum disorder

Apparent motion (AM), the Gestalt perception of motion in the absence of physical motion, was used to study perceptual organization and neurocognitive binding in schizophrenia. Associations between AM perception and psychopathology as well as meaningful subgroups were sought. Circular and stroboscopic AM stimuli were presented to 68 schizophrenia spectrum patients and healthy participants. Psychopathology was measured using the Positive and Negative Syndrome Scale (PANSS). Psychopathology was related to AM perception differentially: Positive and disorganization symptoms were linked to reduced gestalt stability; negative symptoms, excitement and depression had opposite regression weights. Dimensions of psychopathology thus have opposing effects on gestalt perception. It was generally found that AM perception was closely associated with psychopathology. No difference existed between patients and controls, but two latent classes were found. Class A members who had low levels of AM stability made up the majority of inpatients and control subjects; such participants were generally young and male, with short reaction times. Class B typically contained outpatients and some control subjects; participants in class B were older and showed longer reaction times. Hence AM perceptual dysfunctions are not specific for schizophrenia, yet AM may be a promising stage marker.

Investigation of young water inflow in karst aquifers using SF6–CFC–3H/He–85Kr–39Ar and stable isotope components

Karst aquifers are known for their wide distribution of water transfer velocities. From this observation, a multiple geochemical tracer approach seems to be particularly well suited to provide a significant assessment of groundwater flows, but the choice of adapted tracers is essential. In this study, several common tracers in karst aquifers such as physicochemical parameters, major ions, stable isotopes, and d13C to more specific tracers such as dating tracers – 14C, 3H, 3H–3He, CFC-12, SF6 and 85Kr, and 39Ar – were used, in a fractured karstic carbonated aquifer located in Burgundy (France). The information carried by each tracer and the best sampling strategy are compared on the basis of geochemical monitoring done during several recharge events and over longer time periods (months to years). This study’s results demonstrate that at the seasonal and recharge event time scale, the variability of concentrations is low for most tracers due to the broad spectrum of groundwater mixings. The tracers used traditionally for the study of karst aquifers, i.e., physicochemical parameters and major ions, efficiently describe hydrological processes such as the direct and differed recharge, but require being monitored at short time steps during recharge events to be maximized. From stable isotopes, tritium, and Cl� contents, the proportion of the fast direct recharge by the largest porosity was estimated using a binary mixing model. The use of tracers such as CFC-12, SF6, and 85Kr in karst aquifers provides additional information, notably an estimation of apparent age, but they require good preliminary knowledge of the karst system to interpret the results suitably. The CFC-12 and SF6 methods efficiently determine the apparent age of baseflow, but it is preferable to sample the groundwater during the recharge event. Furthermore, these methods are based on different assumptions such as regional enrichment in atmospheric SF6, excess air, and flow models among others. 85Kr and 39Ar concentrations can potentially provide a more direct estimation of groundwater residence time. Conversely, the 3H–3He method is inefficient in the karst aquifer for dating due to 3He degassing.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic to total work ratio were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

Space debris attitude simulation - ιOTA (In-Orbit Tumbling Analysis)

Today, there is little knowledge on the attitude state of decommissioned intact objects in Earth orbit. Observational means have advanced in the past years, but are still limited with respect to an accurate estimate of motion vector orientations and magnitude. Especially for the preparation of Active Debris Removal (ADR) missions as planned by ESA’s Clean Space initiative or contingency scenarios for ESA spacecraft like ENVISAT, such knowledge is needed. ESA's “Debris Attitude Motion Measurements and Modelling” project (ESA Contract No. 40000112447), led by the Astronomical Institute of the University of Bern (AIUB), addresses this problem. The goal of the project is to achieve a good understanding of the attitude evolution and the considerable internal and external effects which occur. To characterize the attitude state of selected targets in LEO and GTO, multiple observation methods are combined. Optical observations are carried out by AIUB, Satellite Laser Ranging (SLR) is performed by the Space Research Institute of the Austrian Academy of Sciences (IWF) and radar measurements and signal level determination are provided by the Fraunhofer Institute for High Frequency Physics and Radar Techniques (FHR). The In-Orbit Tumbling Analysis tool (ιOTA) is a prototype software, currently in development by Hyperschall Technologie Göttingen GmbH (HTG) within the framework of the project. ιOTA will be a highly modular software tool to perform short-(days), medium-(months) and long-term (years) propagation of the orbit and attitude motion (six degrees-of-freedom) of spacecraft in Earth orbit. The simulation takes into account all relevant acting forces and torques, including aerodynamic drag, solar radiation pressure, gravitational influences of Earth, Sun and Moon, eddy current damping, impulse and momentum transfer from space debris or micro meteoroid impact, as well as the optional definition of particular spacecraft specific influences like tank sloshing, reaction wheel behaviour, magnetic torquer activity and thruster firing. The purpose of ιOTA is to provide high accuracy short-term simulations to support observers and potential ADR missions, as well as medium-and long-term simulations to study the significance of the particular internal and external influences on the attitude, especially damping factors and momentum transfer. The simulation will also enable the investigation of the altitude dependency of the particular external influences. ιOTA's post-processing modules will generate synthetic measurements for observers and for software validation. The validation of the software will be done by cross-calibration with observations and measurements acquired by the project partners.

Inclination-dependent changes of the critical shoulder angle significantly influence superior glenohumeral joint stability.

BACKGROUND The critical shoulder angle combines the acromion index and glenoid inclination and has potential to discriminate between shoulders at risk for rotator cuff tear or osteoarthritis and those that are asymptomatic. However, its biomechanics, and particularly the role of the glenoid inclination, are not yet fully understood. METHODS A shoulder simulator was used to analyze the independent influence of glenoid inclination during abduction from 0 to 60°. Spindle motors transferred tension forces by a cable-pulley on human cadaveric humeri. A six-degree-of-freedom force transducer was mounted directly behind the polyethylene glenoid to measure shear and compressive joint reaction force and calculate the instability ratio (ratio of shear and compressive joint reaction force) with the different force ratios of the deltoid and supraspinatus muscles (2:1 and 1:1). A stepwise change in the inclination by 5° increments allowed simulation of a critical shoulder angle range of 20° to 45°. FINDINGS Tilting the glenoid to cranial (increasing the critical shoulder angle) increases the shear joint reaction force and therefore the instability ratio. A balanced force ratio (1:1) between the deltoid and the supraspinatus allowed larger critical shoulder angles before cranial subluxation occurred than did the deltoid-dominant ratio (2:1). INTERPRETATION Glenoid inclination-dependent changes of the critical shoulder angle have a significant impact on superior glenohumeral joint stability. The increased compensatory activity of the rotator cuff to keep the humeral head centered may lead to mechanical overload and could explain the clinically observed association between large angles and degenerative rotator cuff tears.

3D Direct Simulation Monte Carlo Modelling of the Inner Gas Coma of Comet 67P/Churyumov-Gerasimenko: A Parameter Study

Direct Simulation Monte Carlo (DSMC) is a powerful numerical method to study rarefied gas flows such as cometary comae and has been used by several authors over the past decade to study cometary outflow. However, the investigation of the parameter space in simulations can be time consuming since 3D DSMC is computationally highly intensive. For the target of ESA's Rosetta mission, comet 67P/Churyumov-Gerasimenko, we have identified to what extent modification of several parameters influence the 3D flow and gas temperature fields and have attempted to establish the reliability of inferences about the initial conditions from in situ and remote sensing measurements. A large number of DSMC runs have been completed with varying input parameters. In this work, we present the simulation results and conclude on the sensitivity of solutions to certain inputs. It is found that among cases of water outgassing, the surface production rate distribution is the most influential variable to the flow field.