Beamlet based direct aperture optimization for MERT using a photon MLC

PURPOSE A beamlet based direct aperture optimization (DAO) for modulated electron radiotherapy (MERT) using photon multileaf collimator (pMLC) shaped electron fields is developed and investigated. METHODS The Swiss Monte Carlo Plan (SMCP) allows the calculation of dose distributions for pMLC shaped electron beams. SMCP is interfaced with the Eclipse TPS (Varian Medical Systems, Palo Alto, CA) which can thus be included into the inverse treatment planning process for MERT. This process starts with the import of a CT-scan into Eclipse, the contouring of the target and the organs at risk (OARs), and the choice of the initial electron beam directions. For each electron beam, the number of apertures, their energy, and initial shape are defined. Furthermore, the DAO requires dose-volume constraints for the structures contoured. In order to carry out the DAO efficiently, the initial electron beams are divided into a grid of beamlets. For each of those, the dose distribution is precalculated using a modified electron beam model, resulting in a dose list for each beamlet and energy. Then the DAO is carried out, leading to a set of optimal apertures and corresponding weights. These optimal apertures are now converted into pMLC shaped segments and the dose calculation for each segment is performed. For these dose distributions, a weight optimization process is launched in order to minimize the differences between the dose distribution using the optimal apertures and the pMLC segments. Finally, a deliverable dose distribution for the MERT plan is obtained and loaded back into Eclipse for evaluation. For an idealized water phantom geometry, a MERT treatment plan is created and compared to the plan obtained using a previously developed forward planning strategy. Further, MERT treatment plans for three clinical situations (breast, chest wall, and parotid metastasis of a squamous cell skin carcinoma) are created using the developed inverse planning strategy. The MERT plans are compared to clinical standard treatment plans using photon beams and the differences between the optimal and the deliverable dose distributions are determined. RESULTS For the idealized water phantom geometry, the inversely optimized MERT plan is able to obtain the same PTV coverage, but with an improved OAR sparing compared to the forwardly optimized plan. Regarding the right-sided breast case, the MERT plan is able to reduce the lung volume receiving more than 30% of the prescribed dose and the mean lung dose compared to the standard plan. However, the standard plan leads to a better homogeneity within the CTV. The results for the left-sided thorax wall are similar but also the dose to the heart is reduced comparing MERT to the standard treatment plan. For the parotid case, MERT leads to lower doses for almost all OARs but to a less homogeneous dose distribution for the PTV when compared to a standard plan. For all cases, the weight optimization successfully minimized the differences between the optimal and the deliverable dose distribution. CONCLUSIONS A beamlet based DAO using multiple beam angles is implemented and successfully tested for an idealized water phantom geometry and clinical situations.

A critical evaluation of secondary cancer risk models applied to Monte Carlo dose distributions of 2-dimensional, 3-dimensional conformal and hybrid intensity-modulated radiation therapy for breast cancer

The comparison of radiotherapy techniques regarding secondary cancer risk has yielded contradictory results possibly stemming from the many different approaches used to estimate risk. The purpose of this study was to make a comprehensive evaluation of different available risk models applied to detailed whole-body dose distributions computed by Monte Carlo for various breast radiotherapy techniques including conventional open tangents, 3D conformal wedged tangents and hybrid intensity modulated radiation therapy (IMRT). First, organ-specific linear risk models developed by the International Commission on Radiological Protection (ICRP) and the Biological Effects of Ionizing Radiation (BEIR) VII committee were applied to mean doses for remote organs only and all solid organs. Then, different general non-linear risk models were applied to the whole body dose distribution. Finally, organ-specific non-linear risk models for the lung and breast were used to assess the secondary cancer risk for these two specific organs. A total of 32 different calculated absolute risks resulted in a broad range of values (between 0.1% and 48.5%) underlying the large uncertainties in absolute risk calculation. The ratio of risk between two techniques has often been proposed as a more robust assessment of risk than the absolute risk. We found that the ratio of risk between two techniques could also vary substantially considering the different approaches to risk estimation. Sometimes the ratio of risk between two techniques would range between values smaller and larger than one, which then translates into inconsistent results on the potential higher risk of one technique compared to another. We found however that the hybrid IMRT technique resulted in a systematic reduction of risk compared to the other techniques investigated even though the magnitude of this reduction varied substantially with the different approaches investigated. Based on the epidemiological data available, a reasonable approach to risk estimation would be to use organ-specific non-linear risk models applied to the dose distributions of organs within or near the treatment fields (lungs and contralateral breast in the case of breast radiotherapy) as the majority of radiation-induced secondary cancers are found in the beam-bordering regions.

Modelling Nd-isotopes with a coarse resolution ocean circulation model: Sensitivities to model parameters and source/sink distributions

The neodymium (Nd) isotopic composition (Nd) of seawater is a quasi-conservative tracer of water mass mixing and is assumed to hold great potential for paleoceanographic studies. Here we present a comprehensive approach for the simulation of the two neodymium isotopes 143Nd, and 144Nd using the Bern3D model, a low resolution ocean model. The high computational efficiency of the Bern3D model in conjunction with our comprehensive approach allows us to systematically and extensively explore the sensitivity of Nd concentrations and Nd to the parametrisation of sources and sinks. Previous studies have been restricted in doing so either by the chosen approach or by computational costs. Our study thus presents the most comprehensive survey of the marine Nd cycle to date. Our model simulates both Nd concentrations as well as Nd in good agreement with observations. Nd covaries with salinity, thus underlining its potential as a water mass proxy. Results confirm that the continental margins are required as a Nd source to simulate Nd concentrations and Nd consistent with observations. We estimate this source to be slightly smaller than reported in previous studies and find that above a certain magnitude its magnitude affects Nd only to a small extent. On the other hand, the parametrisation of the reversible scavenging considerably affects the ability of the model to simulate both, Nd concentrations and Nd. Furthermore, despite their small contribution, we find dust and rivers to be important components of the Nd cycle. In additional experiments, we systematically varied the diapycnal diffusivity as well as the Atlantic-to-Pacific freshwater flux to explore the sensitivity of Nd concentrations and its isotopic signature to the strength and geometry of the overturning circulation. These experiments reveal that Nd concentrations and Nd are comparatively little affected by variations in diapycnal diffusivity and the Atlantic-to-Pacific freshwater flux. In contrast, an adequate representation of Nd sources and sinks is crucial to simulate Nd concentrations and Nd consistent with observations. The good agreement of our results with observations paves the way for the evaluation of the paleoceanographic potential of Nd in further model studies.

Measurement of distributions sensitive to the underlying event in inclusive Z-boson production in pp collisions at √ s=7 TeV with the ATLAS detector

A measurement of charged-particle distributions sensitive to the properties of the underlying event is presented for an inclusive sample of events containing a Z-boson, decaying to an electron or muon pair. The measurement is based on data collected using the ATLAS detector at the LHC in proton–proton collisions at a centre-of-mass energy of 7 TeV with an integrated luminosity of 4.6fb−1. Distributions of the charged particle multiplicity and of the charged particle transverse momentum are measured in regions of azimuthal angle defined with respect to the Z-boson direction. The measured distributions are compared to similar distributions measured in jet events, and to the predictions of various Monte Carlo generators implementing different underlying event models.

Measurement of the electroweak production of dijets in association with a Z-boson and distributions sensitive to vector boson fusion in proton-proton collisions at √s = 8 TeV using the ATLAS detector

Measurements of fiducial cross sections for the electroweak production of two jets in association with a Z-boson are presented. The measurements are performed using 20.3 fb−1 of proton-proton collision data collected at a centre-of-mass energy of p s = 8TeV by the ATLAS experiment at the Large Hadron Collider. The electroweak component is extracted by a fit to the dijet invariant mass distribution in a fiducial region chosen to enhance the electroweak contribution over the dominant background in which the jets are produced via the strong interaction. The electroweak cross sections measured in two fiducial regions are in good agreement with the Standard Model expectations and the background-only hypothesis is rejected with significance above the 5ơ level. The electroweak process includes the vector boson fusion production of a Z-boson and the data are used to place limits on anomalous triple gauge boson couplings. In addition, measurements of cross sections and differential distributions for inclusive Z-boson-plus-dijet production are performed in five fiducial regions, each with different sensitivity to the electroweak contribution. The results are corrected for detector effects and compared to predictions from the Sherpa and Powheg event generators.

A class of exact Navier-Stokes solutions for homogeneous flat-plate boundary layers and their linear stability

We introduce a new boundary layer formalism on the basis of which a class of exact solutions to the Navier–Stokes equations is derived. These solutions describe laminar boundary layer flows past a flat plate under the assumption of one homogeneous direction, such as the classical swept Hiemenz boundary layer (SHBL), the asymptotic suction boundary layer (ASBL) and the oblique impingement boundary layer. The linear stability of these new solutions is investigated, uncovering new results for the SHBL and the ASBL. Previously, each of these flows had been described with its own formalism and coordinate system, such that the solutions could not be transformed into each other. Using a new compound formalism, we are able to show that the ASBL is the physical limit of the SHBL with wall suction when the chordwise velocity component vanishes while the homogeneous sweep velocity is maintained. A corresponding non-dimensionalization is proposed, which allows conversion of the new Reynolds number definition to the classical ones. Linear stability analysis for the new class of solutions reveals a compound neutral surface which contains the classical neutral curves of the SHBL and the ASBL. It is shown that the linearly most unstable Görtler–Hämmerlin modes of the SHBL smoothly transform into Tollmien–Schlichting modes as the chordwise velocity vanishes. These results are useful for transition prediction of the attachment-line instability, especially concerning the use of suction to stabilize boundary layers of swept-wing aircraft.

Influence of polymers on microstructure and adhesive strength of cementitious tile adhesive mortars

The impact of polymer modification on the physical properties of cementitious mortars is investigated using a multimethod approach. Special emphasis is put on the identification and quantification of different polymer components within the cementitious matrix. With respect to thin-bed applications, particularly tile adhesives, the spatial distributions of latex, cellulose ether (CE), polyvinyl alcohol (PVA), and cement hydration products can be quantified. It is shown that capillary forces and evaporation induce water fluxes in the interconnected part of the pore system, which transport CE, PVA, and cement ions to the mortar interfaces. In contrast, the distribution of latex remains homogeneous. In combination with results from qualitative experiments, the quantitative findings allow reconstruction of the evolution from fresh to hardened mortar, including polymer film formation, cement hydration, and water migration. The resulting microstructure and the failure modes can be correlated with the final adhesive strength of the tile adhesive. The results demonstrate that skinning prior to tile inlaying can strongly reduce wetting properties of the fresh mortar and lower final adhesive strength.

Constraining age distributions of groundwater from public supply wells in diverse hydrogeological settings in Scania, Sweden

Twenty-five public supply wells throughout the hydrogeologically diverse region of Scania, southern Sweden are subjected to environmental tracer analysis (³H–³He,⁴He, CFCs, SF₆ and for one well only also ⁸⁵Kr and ³⁹Ar) to study well and aquifer vulnerability and evaluate possibilities of groundwater age distribution assessment. We find CFC and SF₆ concentrations well above solubility equilibrium with modern atmosphere, indicating local contamination, as well as indications of CFC degradation. The tracer-specific complications considerably constrain possibilities for sound quantitative regional ground- water age distribution assessment and demonstrate the importance of initial qualitative assessment of tracer-specific reliability, as well a need for additional, complementary tracers (e.g. ⁸⁵Kr,³⁹Ar and potentially also ¹⁴C). Lumped parameter modelling yields credible age distribution assessments for representative wells in four type aquifers. Pollution vulnerability of the aquifer types was based on the selected LPM models and qualitative age characterisation. Most vulnerable are unconfined dual porosity and fractured bedrock aquifers, due to a large component of very young groundwater. Unconfined sedimentary aquifers are vulnerable due to young groundwater and a small pre-modern component. Less vulnerable are semi-confined sedimentary or dual-porosity aquifers, due to older age of the modern component and a larger pre-modern component. Confined aquifers appear least vulnerable, due an entirely pre-modern groundwater age distribution (recharged before 1963). Tracer complications aside, environmental tracer analyses and lumped parameter modelling aid in vulnerability assessment and protection of regional groundwater resources.

hp-DGFEM for second-order mixed elliptic problems polyhedra

We prove exponential rates of convergence of hp-version discontinuous Galerkin (dG) interior penalty finite element methods for second-order elliptic problems with mixed Dirichlet-Neumann boundary conditions in axiparallel polyhedra. The dG discretizations are based on axiparallel, σ-geometric anisotropic meshes of mapped hexahedra and anisotropic polynomial degree distributions of μ-bounded variation. We consider piecewise analytic solutions which belong to a larger analytic class than those for the pure Dirichlet problem considered in [11, 12]. For such solutions, we establish the exponential convergence of a nonconforming dG interpolant given by local L 2 -projections on elements away from corners and edges, and by suitable local low-order quasi-interpolants on elements at corners and edges. Due to the appearance of non-homogeneous, weighted norms in the analytic regularity class, new arguments are introduced to bound the dG consistency errors in elements abutting on Neumann edges. The non-homogeneous norms also entail some crucial modifications of the stability and quasi-optimality proofs, as well as of the analysis for the anisotropic interpolation operators. The exponential convergence bounds for the dG interpolant constructed in this paper generalize the results of [11, 12] for the pure Dirichlet case.

Materials properties from electron density distributions: Molecular magnetism and linear optical properties

The general goal of this thesis is correlating observable properties of organic and metal-organic materials with their ground-state electron density distribution. In a long-term view, we expect to develop empirical or semi-empirical approaches to predict materials properties from the electron density of their building blocks, thus allowing to rationally engineering molecular materials from their constituent subunits, such as their functional groups. In particular, we have focused on linear optical properties of naturally occurring amino acids and their organic and metal-organic derivatives, and on magnetic properties of metal-organic frameworks. For analysing the optical properties and the magnetic behaviour of the molecular or sub-molecular building blocks in materials, we mostly used the more traditional QTAIM partitioning scheme of the molecular or crystalline electron densities, however, we have also investigated a new approach, namely, X-ray Constrained Extremely Localized Molecular Orbitals (XC-ELMO), that can be used in future to extracted the electron densities of crystal subunits. With the purpose of rationally engineering linear optical materials, we have calculated atomic and functional group polarizabilities of amino acid molecules, their hydrogen-bonded aggregates and their metal-organic frameworks. This has enabled the identification of the most efficient functional groups, able to build-up larger electric susceptibilities in crystals, as well as the quantification of the role played by intermolecular interactions and coordinative bonds on modifying the polarizability of the isolated building blocks. Furthermore, we analysed the dependence of the polarizabilities on the one-electron basis set and the many-electron Hamiltonian. This is useful for selecting the most efficient level of theory to estimate susceptibilities of molecular-based materials. With the purpose of rationally design molecular magnetic materials, we have investigated the electron density distributions and the magnetism of two copper(II) pyrazine nitrate metal-organic polymers. High-resolution X-ray diffraction and DFT calculations were used to characterize the magnetic exchange pathways and to establish relationships between the electron densities and the exchange-coupling constants. Moreover, molecular orbital and spin-density analyses were employed to understand the role of different magnetic exchange mechanisms in determining the bulk magnetic behaviour of these materials. As anticipated, we have finally investigated a modified version of the X-ray constrained wavefunction technique, XC-ELMOs, that is not only a useful tool for determination and analysis of experimental electron densities, but also enables one to derive transferable molecular orbitals strictly localized on atoms, bonds or functional groups. In future, we expect to use XC-ELMOs to predict materials properties of large systems, currently challenging to calculate from first-principles, such as macromolecules or polymers. Here, we point out advantages, needs and pitfalls of the technique. This work fulfils, at least partially, the prerequisites to understand materials properties of organic and metal-organic materials from the perspective of the electron density distribution of their building blocks. Empirical or semi-empirical evaluation of optical or magnetic properties from a preconceived assembling of building blocks could be extremely important for rationally design new materials, a field where accurate but expensive first-principles calculations are generally not used. This research could impact the community in the fields of crystal engineering, supramolecular chemistry and, of course, electron density analysis.

Environmental correlates of species rank − abundance distributions in global drylands

Theoretical models predict lognormal species abundance distributions (SADs) in stable and productive environments, with log-series SADs in less stable, dispersal driven communities. We studied patterns of relative species abundances of perennial vascular plants in global dryland communities to: (i) assess the influence of climatic and soil characteristics on the observed SADs, (ii) infer how environmental variability influences relative abundances, and (iii) evaluate how colonisation dynamics and environmental filters shape abundance distributions. We fitted lognormal and log-series SADs to 91 sites containing at least 15 species of perennial vascular plants. The dependence of species relative abundances on soil and climate variables was assessed using general linear models. Irrespective of habitat type and latitude, the majority of the SADs (70.3%) were best described by a lognormal distribution. Lognormal SADs were associated with low annual precipitation, higher aridity, high soil carbon content, and higher variability of climate variables and soil nitrate. Our results do not corroborate models predicting the prevalence of log-series SADs in dryland communities. As lognormal SADs were particularly associated with sites with drier conditions and a higher environmental variability, we reject models linking lognormality to environmental stability and high productivity conditions. Instead our results point to the prevalence of lognormal SADs in heterogeneous environments, allowing for more evenly distributed plant communities, or in stressful ecosystems, which are generally shaped by strong habitat filters and limited colonisation. This suggests that drylands may be resilient to environmental changes because the many species with intermediate relative abundances could take over ecosystem functioning if the environment becomes suboptimal for dominant species.

Cold-Water Coral Distributions in the Drake Passage Area from Towed Camera Observations - Initial Interpretations

Seamounts are unique deep-sea features that create habitats thought to have high levels of endemic fauna, productive fisheries and benthic communities vulnerable to anthropogenic impacts. Many seamounts are isolated features, occurring in the high seas, where access is limited and thus biological data scarce. There are numerous seamounts within the Drake Passage (Southern Ocean), yet high winds, frequent storms and strong currents make seafloor sampling particularly difficult. As a result, few attempts to collect biological data have been made, leading to a paucity of information on benthic habitats or fauna in this area, particularly those on primarily hard-bottom seamounts and ridges. During a research cruise in 2008 six locations were examined (two on the Antarctic margin, one on the Shackleton Fracture Zone, and three on seamounts within the Drake Passage), using a towed camera with onboard instruments to measure conductivity, temperature, depth and turbidity. Dominant fauna and bottom type were categorized from 200 randomized photos from each location. Cold-water corals were present in high numbers in habitats both on the Antarctic margin and on the current swept seamounts of the Drake Passage, though the diversity of orders varied. Though the Scleractinia (hard corals) were abundant on the sedimented margin, they were poorly represented in the primarily hard-bottom areas of the central Drake Passage. The two seamount sites and the Shackleton Fracture Zone showed high numbers of stylasterid (lace) and alcyonacean (soft) corals, as well as large numbers of sponges. Though data are preliminary, the geological and environmental variability (particularly in temperature) between sample sites may be influencing cold-water coral biogeography in this region. Each area observed also showed little similarity in faunal diversity with other sites examined for this study within all phyla counted. This manuscript highlights how little is understood of these isolated features, particularly in Polar regions.

Seasonal Distributions of Satellite-Measured Phytoplankton Pigment Concentration Along the Chilean Coast

Five years (1979-1983) of Coastal Zone Color Scanner satellite ocean color data are used to examine seasonal patterns of phytoplankton pigment concentration along the Chilean coast from 20 degrees S to 45 degrees S. Four kilometer resolution, 2-4 day composites document the presence of filaments of elevated pigment concentration extending offshore throughout the study area, with maximum offshore extension at higher latitudes. In three years, 1979, 1981, and 1983, sufficient data exist in monthly composites to allow recreation of portions of the seasonal cycle. Data in 1979 are the most complete. Near-shore concentrations and cross-shelf extension of pigment concentrations in 1979 are maximum in austral winter throughout the study area and minimum in summer. Available data from 1981 and 1983 are consistent with this temporal pattern but with concentrations approximately double those of 1979. Seasonal, spatial patterns within 10 km of shore and 50 km offshore indicate a latitudinal discontinuity both in absolute concentration and in the magnitude of the seasonal cycle at approximately 33 degrees S in both 1979 and in the climatological time series. The discontinuity is strongest ill fall-winter and weakest in summer. South of this latitude, concentrations are relatively high (2-3 mg m(-3) in 1979), a strong seasonal cycle is present, and patterns 50 km offshore are correlated with those within 10 km of shore. North of 33 degrees S, concentrations are < 1.5 mg m(-3) (in 1979), and the seasonal cycle within 10 km of shore is present but much weaker and less obviously correlated with that 50 km offshore. The seasonal cycle of pigment concentrations is 180 degrees out of phase with monthly averaged upwelling favorable winds. Noncoincident Pathfinder sea surface temperature data show that over most latitudes, coastal low surface temperatures lag wind forcing by 1-2 months, but these too are out of phase with the pigment seasonal cycle. These data point to control of pigment patterns along the Chilean coast by the interaction of upwelling with circulation patterns unconnected to local wind forcing.

Use of Otolith Strontium:Calcium Ratios for Hindcasting Larval Cod Gadus Morhua Distributions Relative to Water Masses on Georges Bank

The concentration ratios of strontium to calcium in laboratory-reared larval cod otoliths are shown to be related to the water temperature (T) at the time of otolith precipitation. This relationship is curvilinear, and is best described by a simple exponential equation of the form (Sr/Ca x 1000 = a exp(-T/b). We show that when Sr/Ca elemental analyses are related to the daily growth increments in the larval otoliths, relative temperature histories of individual field-caught larvae can be reconstructed from the egg stage to the time of capture. We present preliminary examples of how such reconstructed temperature histories of Atlantic cod Gadus morhua larvae, collected on Georges Bank during April and May 1993, may be interpreted in relation to the broad-scale larval distributions and the hydrography of the Bank.

Distributions and Variability of Particulate Organic Matter in a Coastal Upwelling System

In this study we examined the spatial and temporal variability of particulate organic material (POM) off Oregon during the upwelling season. High-resolution vertical profiling of beam attenuation was conducted along two cross-shelf transects. One transect was located in a region where the shelf is relatively uniform and narrow (off Cascade Head (CH)); the second transect was located in a region where the shelf is shallow and wide (off Cape Perpetua (CP)). In addition, water samples were collected for direct analysis of chlorophyll, particulate organic carbon (POC), and particulate organic nitrogen (PON). Beam attenuation was highly correlated with POC and PON. Striking differences in distribution patterns and characteristics of POM were observed between CH and CP. Off CH, elevated concentrations of chlorophyll and POC were restricted to the inner shelf and were highly variable in time. The magnitude of the observed short-term temporal variability was of the same order as that of the seasonal variability reported in previous studies. Elevated concentrations of nondegraded chlorophyll and POM were observed near the bottom. Downwelling and rapid sinking are two mechanisms by which phytoplankton cells can be delivered to the bottom before being degraded. POM may be then transported across the shelf via the benthic nepheloid layer. Along the CP transect, concentrations of POM were generally higher than they were along the CH transect and extended farther across the shelf. Characteristics of surface POM, namely, C: N ratios and carbon: chlorophyll ratios, differed between the two sites. These differences can be attributed to differences in shelf circulation.