Surface dosimetry for very small X-ray beams

Introduction The dose to skin surface is an important factor for many radiotherapy treatment techniques. It is known that TPS predicted surface doses can be significantly different from actual ICRP skin doses as defined at 70 lm. A number of methods have been implemented for the accurate determination of surface dose including use of specific dosimeters such as TLDs and radiochromic film as well as Monte Carlo calculations. Stereotactic radiosurgery involves delivering very high doses per treatment fraction using small X-ray fields. To date, there has been limited data on surface doses for these very small field sizes. The purpose of this work is to evaluate surface doses by both measurements and Monte Carlo calculations for very small field sizes. Methods All measurements were performed on a Novalis Tx linear accelerator which has a 6 MV SRS X-ray beam mode which uses a specially thin flattening filter. Beam collimation was achieved by circular cones with apertures that gave field sizes ranging from 4 to 30 mm at the isocentre. The relative surface doses were measured using Gafchromic EBT3 film which has the active layer at a depth similar to the ICRP skin dose depth. Monte Carlo calculations were performed using the BEAMnrc/EGSnrc Monte Carlo codes (V4 r225). The specifications of the linear accelerator, including the collimator, were provided by the manufacturer. Optimisation of the incident X-ray beam was achieved by an iterative adjustment of the energy, spatial distribution and radial spread of the incident electron beam striking the target. The energy cutoff parameters were PCUT = 0.01 MeV and ECUT = 0.700 - MeV. Directional bremsstrahlung splitting was switched on for all BEAMnrc calculations. Relative surface doses were determined in a layer defined in a water phantom of the same thickness and depth as compared to the active later in the film. Results Measured surface doses using the EBT3 film varied between 13 and 16 % for the different cones with an uncertainty of 3 %. Monte Carlo calculated surface doses were in agreement to better than 2 % to the measured doses for all the treatment cones. Discussion and conclusions This work has shown the consistency of surface dose measurements using EBT3 film with Monte Carlo predicted values within the uncertainty of the measurements. As such, EBT3 film is recommended for in vivo surface dose measurements.

Importance of orientational rearrangement during vitrification of hydrocarbons: dependence on molecular shape

On the basis of Monte Carlo calculations of 2,2-dimethylpropane (neopentane), n-pentane, and 2,2-dimethylbutane (neohexane) at several temperatures, thermodynamic properties and radial distribution functions as well as dimerization and bonding energy distribution functions are reported for both liquid and glassy states. Changes in the radial distribution functions on cooling depend on whether the groups are accessible (peripheral) or inaccessible. Peaks in the radial distribution functions corresponding to peripheral groups do not shift to lower distances on cooling and at times display a large increase in the intensity of the first peak. The peaks due to inaccessible groups, on the other hand, shift to lower distances on cooling. The magnitude of the reorientational contribution in determining the resulting structure of the glass is estimated for the different hydrocarbon molecules investigated. The reorientational contribution is highest for neopentane (26%) followed by isopentane (24%), neohexane (22%), and n-pentane (0%). It appears that molecular geometry has an important role in determining the magnitude of the reorientational contribution to the structure of the glass.

Sensory-organ-like response determines the magnetism of zigzag-edged honeycomb nanoribbons

We present an analytical effective theory for the magnetic phase diagram for zigzag-edge terminated honeycomb nanoribbons described by a Hubbard model with an interaction parameter U. We show that the edge magnetic moment varies as ln U and uncover its dependence on the width W of the ribbon. The physics of this owes its origin to the sensory-organ-like response of the nanoribbons, demonstrating that considerations beyond the usual Stoner-Landau theory are necessary to understand the magnetism of these systems. A first-order magnetic transition from an antiparallel orientation of the moments on opposite edges to a parallel orientation occurs upon doping with holes or electrons. The critical doping for this transition is shown to depend inversely on the width of the ribbon. Using variational Monte Carlo calculations, we show that magnetism is robust to fluctuations. Additionally, we show that the magnetic phase diagram is generic to zigzag-edge terminated nanostructures such as nanodots. Furthermore, we perform first-principles modeling to show how such magnetic transitions can be realized in substituted graphene nanoribbons. DOI: 10.1103/PhysRevB.87.085412

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

Consistent nucleon-nucleon potentials and three-body forces

We construct microscopic three-nucleon forces consistent with the Bonn and Nijmegen two-nucleon potentials, and including , Roper, and nucleon-antinucleon excitations. Recent results for the choice of the meson parameters are discussed. The forces are used in Brueckner calculations and the saturation properties of nuclear matter are determined.

Three-body force effect on P-3 F-2 neutron superfluidity in neutron matter, neutron star matter, and neutron stars

We investigate the effect of microscopic three-body forces on the P-3 F-2 neutron superfluidity in neutron matter, beta-stable neutron star matter, and neutron stars by using the BCS theory and the Brueckner-Hartree-Fock approach. We adopt the Argonne V18 potential supplemented with a microscopic three-body force as the realistic nucleon-nucleon interaction. We have concentrated on studying the three-body force effect on the P-3 F-2 neutron pairing gap. It is found that the three-body force effect considerably enhances the P-3 F-2 neutron superfluidity in neutron star matter and neutron stars.

Microscopic three-body forces and kaon condensation in cold neutrino-trapped matter

We investigate the composition and the equation of state of the kaon condensed phase in neutrino-free and neutrino-trapped star matter within the framework of the Brueckner-Hartree-Fock approach with three-body forces. We find that neutrino trapping shifts the onset density of kaon condensation to a larger baryon density, and reduces considerably the kaon abundance. As a consequence, when kaons are allowed, the equation of state of neutrino-trapped star matter becomes stiffer than the one of neutrino free matter. The effects of different three-body forces are compared and discussed. Neutrino trapping turns out to weaken the role played by the symmetry energy in determining the composition of stellar matter, and thus reduces the difference between the results obtained by using different three-body forces.

Measuring distributions of jump rates in disordered metal–hydrogen systems by nuclear magnetic relaxation

Monte Carlo calculations of the nuclear magnetic relaxation rate in a disordered metal–hydrogen system having a distribution of jump rates are reported. The calculations deal specifically with the spin-locked rotating-frame relaxation time T1ρ. The results demonstrate that the temperature variation of the rate is only weakly dependent on the distribution and it is therefore unlikely that the jump rate distribution can be extracted from relaxation measurements in which temperature is the main variable. It is shown that the alternative of measuring the relaxation rate over a wide range of spin-locking field strengths at a constant temperature can lead to an evaluation of the distribution.

Shifts in electron capture to the continuum at low collision energies: Enhanced role of target postcollision interactions

Measurements of electron velocity distributions emitted at 0degrees for collisions of 10- and 20-keV H+ incident ions on H-2 and He show that the electron capture to the continuum cusp formation, which is still possible at these low impact energies, is shifted to lower momenta than its standard position (centered on the projectile velocity), as recently predicted. Classical trajectory Monte Carlo calculations reproduce the observations remarkably well, and indicate that a long-range residual interaction of the electron with the target ion after ionization is responsible for the shifts, which is a general effect that is enhanced at low nuclear velocities.

On the equation of state and atomic mean square displacement of crystals

We have presented a Green's function method for the calculation of the atomic mean square displacement (MSD) for an anharmonic Hamil toni an . This method effectively sums a whole class of anharmonic contributions to MSD in the perturbation expansion in the high temperature limit. Using this formalism we have calculated the MSD for a nearest neighbour fcc Lennard Jones solid. The results show an improvement over the lowest order perturbation theory results, the difference with Monte Carlo calculations at temperatures close to melting is reduced from 11% to 3%. We also calculated the MSD for the Alkali metals Nat K/ Cs where a sixth neighbour interaction potential derived from the pseudopotential theory was employed in the calculations. The MSD by this method increases by 2.5% to 3.5% over the respective perturbation theory results. The MSD was calculated for Aluminum where different pseudopotential functions and a phenomenological Morse potential were used. The results show that the pseudopotentials provide better agreement with experimental data than the Morse potential. An excellent agreement with experiment over the whole temperature range is achieved with the Harrison modified point-ion pseudopotential with Hubbard-Sham screening function. We have calculated the thermodynamic properties of solid Kr by minimizing the total energy consisting of static and vibrational components, employing different schemes: The quasiharmonic theory (QH), ).2 and).4 perturbation theory, all terms up to 0 ().4) of the improved self consistent phonon theory (ISC), the ring diagrams up to o ().4) (RING), the iteration scheme (ITER) derived from the Greens's function method and a scheme consisting of ITER plus the remaining contributions of 0 ().4) which are not included in ITER which we call E(FULL). We have calculated the lattice constant, the volume expansion, the isothermal and adiabatic bulk modulus, the specific heat at constant volume and at constant pressure, and the Gruneisen parameter from two different potential functions: Lennard-Jones and Aziz. The Aziz potential gives generally a better agreement with experimental data than the LJ potential for the QH, ).2, ).4 and E(FULL) schemes. When only a partial sum of the).4 diagrams is used in the calculations (e.g. RING and ISC) the LJ results are in better agreement with experiment. The iteration scheme brings a definitive improvement over the).2 PT for both potentials.

Nonlinear Schrodinger equation for a superfluid Fermi gas in the BCS-BEC crossover

We introduce a quasianalytic nonlinear Schrodinger equation with beyond mean-field corrections to describe the dynamics of a zero-temperature dilute superfluid Fermi gas in the crossover from the weak-coupling Bardeen-Cooper-Schrieffer (BCS) regime, where k(F)parallel to a parallel to << 1 with a the s-wave scattering length and k(F) the Fermi momentum, through the unitarity limit k(F)a ->+/-infinity to the Bose-Einstein condensation (BEC) regime where k(F)a > 0. The energy of our model is parametrized using the known asymptotic behavior in the BCS, BEC, and the unitarity limits and is in excellent agreement with accurate Green's-function Monte Carlo calculations. The model generates good results for frequencies of collective breathing oscillations of a trapped Fermi superfluid.

Effective pair potentials and structure of water

Classical Monte Carlo calculations have been performed in order to investigate the ability of the TIP4P, SPC, and SPCE water models to reproduce the structural features of liquid water. The simulations were carried out in the NPT ensemble at 4 thermodynamic conditions. The results are compared with recent neutron diffraction data. Essentially, the three models capture equally well the thermodynamic and structural features of water. Although they were parametrized to reproduce the water properties at ambient conditions, the agreement with the experimental pair correlation functions was even better at supercritical conditions. This is because the effective pair potentials have some difficulty to reproduce cooperative interactions, like hydrogen bonds. These interactions are less effective at supercritical conditions, where the liquid behaves roughly like a gas.

S-, P-, and D-wave resonances in positronium-lithium-atom scattering

We investigate ortho-positronium-lithium-atom (Ps-Li) scattering using static-exchange and three-Ps-state coupled-channel calculations. The present three-PS-state scheme, while closely agreeing with the resonance and binding energies in the Ps-H system, predicts S-, P-, and D-wave resonances at 4.25 eV, 4.9 eV, and, 5.25 eV. respectively, in the electronic spin-singlet channel of Ps-Li scattering. The present calculation also yields a Ps-Li binding in this attractive singlet channel with an approximate binding energy of 0.218 eV, which is in adherence with the recent findings of a chemically stable PsLi system using stocastic variational and quantum Monte Carlo calculations. We further report elastic, Ps(2s)-, and Ps(2p)-excitation cross sections at low to medium energies (0.068-30 eV).

Neutron field characteristics of Ciemat's Neutron Standards Laboratory Hector Rene Vega-Carrillo

Monte Carlo calculations were carried out to characterize the neutron field produced by the calibration neutron sources of the Neutron Standards Laboratory at the Research Center for Energy, Environment and Technology (CIEMAT) in Spain. For 241AmBe and 252Cf neutron sources, the neutron spectra, the ambient dose equivalent rates and the total neutron fluence rates were estimated. In the calibration hall, there are several items that modify the neutron field. To evaluate their effects different cases were used, from point-like source in vacuum up to the full model. Additionally, using the full model, the neutron spectra were estimated to different distances along the bench; with these spectra, the total neutron fluence and the ambient dose equivalent rates were calculated. The hall walls induce the largest changes in the neutron spectra and the respective integral quantities. The free-field neutron spectrum is modified due the room return effect.

The effect of protein concentration on ion binding.

The concentration of protein in a solution has been found to have a significant effect on ion binding affinity. It is well known that an increase in ionic strength of the solvent medium by addition of salt modulates the ion-binding affinity of a charged protein due to electrostatic screening. In recent Monte Carlo simulations, a similar screening has been detected to arise from an increase in the concentration of the protein itself. Experimental results are presented here that verify the theoretical predictions; high concentrations of the negatively charged proteins calbindin D9k and calmodulin are found to reduce their affinity for divalent cations. The Ca(2+)-binding constant of the C-terminal site in the Asn-56 --> Ala mutant of calbindin D9k has been measured at seven different protein concentrations ranging from 27 microM to 7.35 mM by using 1H NMR. A 94% reduction in affinity is observed when going from the lowest to the highest protein concentration. For calmodulin, we have measured the average Mg(2+)-binding constant of sites I and II at 0.325, 1.08, and 3.25 mM protein and find a 13-fold difference between the two extremes. Monte Carlo calculations have been performed for the two cases described above to provide a direct comparison of the experimental and simulated effects of protein concentration on metal ion affinities. The overall agreement between theory and experiment is good. The results have important implications for all biological systems involving interactions between charged species.