The development of a Monte Carlo system to verify radiotherapy treatment dose calculations

Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the “gold standard” for predicting dose deposition in the patient. In this study, software has been developed that enables the transfer of treatment plan information from the treatment planning system to a Monte Carlo dose calculation engine. A database of commissioned linear accelerator models (Elekta Precise and Varian 2100CD at various energies) has been developed using the EGSnrc/BEAMnrc Monte Carlo suite. Planned beam descriptions and CT images can be exported from the treatment planning system using the DICOM framework. The information in these files is combined with an appropriate linear accelerator model to allow the accurate calculation of the radiation field incident on a modelled patient geometry. The Monte Carlo dose calculation results are combined according to the monitor units specified in the exported plan. The result is a 3D dose distribution that could be used to verify treatment planning system calculations. The software, MCDTK (Monte Carlo Dicom ToolKit), has been developed in the Java programming language and produces BEAMnrc and DOSXYZnrc input files, ready for submission on a high-performance computing cluster. The code has been tested with the Eclipse (Varian Medical Systems), Oncentra MasterPlan (Nucletron B.V.) and Pinnacle3 (Philips Medical Systems) planning systems. In this study the software was validated against measurements in homogenous and heterogeneous phantoms. Monte Carlo models are commissioned through comparison with quality assurance measurements made using a large square field incident on a homogenous volume of water. This study aims to provide a valuable confirmation that Monte Carlo calculations match experimental measurements for complex fields and heterogeneous media.

First principle study of hydrogenation of MgB2 : an important step toward reversible hydrogen storage in the coupled LiBH4/MgH2 system

Recent experiments [F. E. Pinkerton, M. S. Meyer, G. P. Meisner, M. P. Balogh, and J. J. Vajo, J. Phys. Chem. C 111, 12881 (2007) and J. J. Vajo and G. L. Olson, Scripta Mater. 56, 829 (2007)] demonstrated that the recycling of hydrogen in the coupled LiBH4/MgH2 system is fully reversible. The rehydrogenation of MgB2 is an important step toward the reversibility. By using ab initio density functional theory calculations, we found that the activation barrier for the dissociation of H2 are 0.49 and 0.58 eV for the B and Mg-terminated MgB2(0001) surface, respectively. This implies that the dissociation kinetics of H2 on a MgB2 (0001) surface should be greatly improved compared to that in pure Mg materials. Additionally, the diffusion of dissociated H atom on the Mg-terminated MgB2(0001) surface is almost barrier-less. Our results shed light on the experimentally-observed reversibility and improved kinetics for the coupled LiBH4/MgH2 system.

Exact solutions of coupled multispecies linear reaction–diffusion equations on a uniformly growing domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.

Direct estimation of an element of order matrix from H-1 NMR spectra of strongly dipolar coupled spins

NMR spectra of molecules oriented in thermotropic liquid crystalline media provide information on the molecular structure and order. The spins are generally strongly dipolar coupled and the spectral analyse require the tedious and time consuming numerical iterative calculations. The present study demonstrates the application of multiple quantum spin state selective detection of single quantum transitions for mimicking the homonuclear decoupling and the direct estimation of an element of ordering matrix. This information is utilized to estimate the nearly accurate starting dipolar couplings for iterative calculations. The studies on the spectra of strongly dipolar coupled five and six interacting spin systems are reported.

Feasibility of a stepped wedge cluster RCT and concurrent observational sub-study to evaluate the effects of modified ward night lighting on inpatient fall rates and sleep quality: A protocol for a pilot trial

Background: Falls among hospitalised patients impose a considerable burden on health systems globally and prevention is a priority. Some patient-level interventions have been effective in reducing falls, but others have not. An alternative and promising approach to reducing inpatient falls is through the modification of the hospital physical environment and the night lighting of hospital wards is a leading candidate for investigation. In this pilot trial, we will determine the feasibility of conducting a main trial to evaluate the effects of modified night lighting on inpatient ward level fall rates. We will test also the feasibility of collecting novel forms of patient level data through a concurrent observational sub-study. Methods/design: A stepped wedge, cluster randomised controlled trial will be conducted in six inpatient wards over 14 months in a metropolitan teaching hospital in Brisbane (Australia). The intervention will consist of supplementary night lighting installed across all patient rooms within study wards. The planned placement of luminaires, configurations and spectral characteristics are based on prior published research and pre-trial testing and modification. We will collect data on rates of falls on study wards (falls per 1000 patient days), the proportion of patients who fall once or more, and average length of stay. We will recruit two patients per ward per month to a concurrent observational sub-study aimed at understanding potential impacts on a range of patient sleep and mobility behaviour. The effect on the environment will be monitored with sensors to detect variation in light levels and night-time room activity. We will also collect data on possible patient-level confounders including demographics, pre-admission sleep quality, reported vision, hearing impairment and functional status. Discussion: This pragmatic pilot trial will assess the feasibility of conducting a main trial to investigate the effects of modified night lighting on inpatient fall rates using several new methods previously untested in the context of environmental modifications and patient safety. Pilot data collected through both parts of the trial will be utilised to inform sample size calculations, trial design and final data collection methods for a subsequent main trial.

A strategy for scheduling tightly coupled parallel applications on clusters

Although various strategies have been developed for scheduling parallel applications with independent tasks, very little work exists for scheduling tightly coupled parallel applications on cluster environments. In this paper, we compare four different strategies based on performance models of tightly coupled parallel applications for scheduling the applications on clusters. In addition to algorithms based on existing popular optimization techniques, we also propose a new algorithm called Box Elimination that searches the space of performance model parameters to determine the best schedule of machines. By means of real and simulation experiments, we evaluated the algorithms on single cluster and multi-cluster setups. We show that our Box Elimination algorithm generates up to 80% more efficient schedule than other algorithms. We also show that the execution times of the schedules produced by our algorithm are more robust against the performance modeling errors.

Stability of multi-cluster synchronization

A system of many coupled oscillators on a network can show multicluster synchronization. We obtain existence conditions and stability bounds for such a multicluster synchronization. When the oscillators are identical, we obtain the interesting result that network structure alone can cause multicluster synchronization to emerge even when all the other parameters are the same. We also study occurrence of multicluster synchronization when two different types of oscillators are coupled.

Synthesis and Characterization of a Class of Donor-Acceptor Conjugated Molecules: Experiments and Theoretical Calculations

A class of conjugated molecules containing donor (thiophene) and acceptor (malononitrile) is synthesized by Knoevenagel condensation reaction between 2-(2,6-dimethy1-4H-pyran-4-ylidene) malononitrile and thiophene carbaldehyde containing two and three thiophene units. The resulting molecules are characterized by H-1 and C-13 NMR. We have performed UV-vis absorption, fluorescence, and cyclic voltammetry measurements on these materials. The spectroscopic and electrochemical measurements proved beyond doubt that these materials possess lowexcitation gap and are suitable for being an active material in various electronic devices. We have also performed electronic structure calculations using density functional theory (DFT) and INDO/SCI methods to characterize the ground and excited states of this class of molecules. These donor-acceptor molecules show a strong charge transfercharacter that increases with the increase in the number of thiophene rings coupled to the malononitrile acceptor moiety. We have also calculated the pi-coherence length, Stoke's shift, and effect of solvents on excited states for this class of molecules, Our theoretical values agree well with experimental results.

Electronic structure and bonding of beta-rhombohedral boron using cluster fragment approach

Theoretical studies using density functional theory are carried out to understand the electronic structure and bonding and electronic properties of elemental beta-rhombohedral boron. The calculated band structure of ideal beta-rhombohedral boron (B-105) shows valence electron deficiency and depicts metallic behavior. This is in contrast to the experimental result that it is a semiconductor. To understand this ambiguity we discuss the electronic structure and bonding of this allotrope with cluster fragment approach using our recently proposed mno rule. This helps us to comprehend in greater detail the structure of B-105 and materials which are closely related to beta-rhombohedral boron. The molecular structures B12H12-2, B28H21+1, BeB27H21, LiB27H21-1, CB27H21+2, B57H36+3, Be3B54H36, and Li2CB54H36, and corresponding solids Li8Be3B102 and Li10CB102 are arrived at using these ideas and studied using first principles density functional theory calculations.

Synthesis, structure and magnetic properties of the polyoxovanadate cluster Zn-2(NH2(CH2)(2)NH2)(5)]{Zn(NH2(CH2)(2)NH2)(2)}(2){V18O42(H2O)\ ]center dot xH(2)O (x similar to 12), possessing a layered structure

A hydrothermal reaction of a mixture of ZnCl2, V2O5, ethylenediamine and water gave rise to a layered poly oxovanadate material. clusters. These clusters, with all the vanadium ions in the +4 state, are connected together through Zn(NH2(CH2)(2)NH2)(2) linkers forming a two-dimensional structure. The layers are also separated by distorted trigonal bipyramidal [Zn-2(NH2(CH2)(2)NH2)(5)] complexes. The Structure, thus, presents a dual role for the Zn-ethylenediamine complex. The magnetic susceptibility studies indicate that the interactions between the V centres in I are predominantly antiferromagnetic in nature and the compound shows highly frustrated behaviour. The magnetic properties are compared to the theoretical calculations based oil the Heisenberg model, in addition to correlating to the structure. Crystal data for the complexes are presented.

MNDO study of coupled strained hydrocarbons: Additivity of exocyclic C---C bond contraction induced by different ring systems

MNDO geometry optimizations were carried out on a series of symmetrically and unsymmetrically coupled strained ring hydrocarbons, R1-R1 and R1–R2 (R1=methyl, cyclopropyl, 1-bicyclo[1.1.0]butyl, 1-bicyclo[1.1.1]pentyl, prismyl, cubyl, 6-tricyclo [3.1.1.03,6]heptyl, and tetrahedryl groups; R2=methyl and cyclopropyl). The remarkable contraction of the C---C bond connecting the strained rings found experimentally in a few cases was reproduced correctly by the calculations. A linear correlation was found between the bond length shortening and the bond angle widening at the corresponding carbon atoms for all the structures considered. The reduction in C---C bond lengths due to various ring systems is additive. The additivity indicates that inter-ring interactions which effect the central bond length are absent and confirms the common electronic origin of bond contraction in these systems, viz. enhanced s-character in the exocyclic bonds of strained rings.

Effect of Sorbate-Zeolite Interaction on Cluster Lifetime and Size of Sorbates: Xe in NaY

Molecular dynamics calculations are reported for Xe in sodium Y zeolite with varying strengths of sorbate-zeolite dispersion interaction. In the absence of any dispersion interaction between the sorbate and the zeolite, the presence of the zeolite has a purely geometrical role. Increase in the strength of the sorbate-zeolite interaction increases the monomer population and decreases the population of dimers and higher sized clusters. The lifetime of the monomers as well as dimers increases with the strength of the dispersion interaction. The observed variations in the lifetime and the population of the different sized clusters is explained in terms of the changes in the potential energy surface caused by the increase in the strength of the dispersion interaction.

Coupled phonons, magnetic excitations, and ferroelectricity in AlFeO3: Raman and first-principles studies

We determine the nature of coupled phonons and magnetic excitations in AlFeO3 using inelastic light scattering from 5 to 315 K covering a spectral range from 100 to 2200 cm(-1) and complementary first-principles density functional theory-based calculations. A strong spin-phonon coupling and magnetic ordering-induced phonon renormalization are evident in (1) anomalous temperature dependence of many modes with frequencies below 850 cm(-1), particularly near the magnetic transition temperature T-c approximate to 250 K, and (2) distinct changes in band positions of high-frequency Raman bands between 1100 and 1800 cm(-1); in particular, a broad mode near 1250 cm(-1) appears only below T-c, attributed to the two-magnon Raman scattering. We also observe weak anomalies in the mode frequencies similar to 100 K due to a magnetically driven ferroelectric phase transition. Understanding of these experimental observations has been possible on the basis of first-principles calculations of the phonons' spectrum and their coupling with spins.

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

Finite differences and a coupled analytic technique with applications to explosions and earthquakes

An analytic technique is developed that couples to finite difference calculations to extend the results to arbitrary distance. Finite differences and the analytic result, a boundary integral called two-dimensional Kirchhoff, are applied to simple models and three seismological problems dealing with data. The simple models include a thorough investigation of the seismologic effects of a deep continental basin. The first problem is explosions at Yucca Flat, in the Nevada test site. By modeling both near-field strong-motion records and teleseismic P-waves simultaneously, it is shown that scattered surface waves are responsible for teleseismic complexity. The second problem deals with explosions at Amchitka Island, Alaska. The near-field seismograms are investigated using a variety of complex structures and sources. The third problem involves regional seismograms of Imperial Valley, California earthquakes recorded at Pasadena, California. The data are shown to contain evidence of deterministic structure, but lack of more direct measurements of the structure and possible three-dimensional effects make two-dimensional modeling of these data difficult.