SPH calculations of asteroid disruptions: The role of pressure dependent failure models

We present recent improvements of the modeling of the disruption of strength dominated bodies using the Smooth Particle Hydrodynamics (SPH) technique. The improvements include an updated strength model and a friction model, which are successfully tested by a comparison with laboratory experiments. In the modeling of catastrophic disruptions of asteroids, a comparison between old and new strength models shows no significant deviation in the case of targets which are initially non-porous, fully intact and have a homogeneous structure (such as the targets used in the study by Benz and Asphaug, 1999). However, for many cases (e.g. initially partly or fully damaged targets and rubble-pile structures) we find that it is crucial that friction is taken into account and the material has a pressure dependent shear strength. Our investigations of the catastrophic disruption threshold (27, as a function of target properties and target sizes up to a few 100 km show that a fully damaged target modeled without friction has a Q(D)*:, which is significantly (5-10 times) smaller than in the case where friction is included. When the effect of the energy dissipation due to compaction (pore crushing) is taken into account as well, the targets become even stronger (Q(D)*; is increased by a factor of 2-3). On the other hand, cohesion is found to have an negligible effect at large scales and is only important at scales less than or similar to 1 km. Our results show the relative effects of strength, friction and porosity on the outcome of collisions among small (less than or similar to 1000 km) bodies. These results will be used in a future study to improve existing scaling laws for the outcome of collisions (e.g. Leinhardt and Stewart, 2012). (C) 2014 Elsevier Ltd. All rights reserved.

Exact moment calculations for genetic models with migration, mutation, and drift

Using properties of moment stationarity we develop exact expressions for the mean and covariance of allele frequencies at a single locus for a set of populations subject to drift, mutation, and migration. Some general results can be obtained even for arbitrary mutation and migration matrices, for example: (1) Under quite general conditions, the mean vector depends only on mutation rates, not on migration rates or the number of populations. (2) Allele frequencies covary among all pairs of populations connected by migration. As a result, the drift, mutation, migration process is not ergodic when any finite number of populations is exchanging genes. in addition, we provide closed form expressions for the mean and covariance of allele frequencies in Wright's finite-island model of migration under several simple models of mutation, and we show that the correlation in allele frequencies among populations can be very large for realistic rates of mutation unless an enormous number of populations are exchanging genes. As a result, the traditional diffusion approximation provides a poor approximation of the stationary distribution of allele frequencies among populations. Finally, we discuss some implications of our results for measures of population structure based on Wright's F-statistics.

Verification of intensity modulated stereotactic radiotherapy using Monte Carlo calculations and EPID dosimetry

The purpose of this work was to develop a comprehensive IMSRT QA procedure that examined, using EPID dosimetry and Monte Carlo (MC) calculations, each step in the treatment planning and delivery process. These steps included verification of the field shaping, treatment planning system (RTPS) dose calculations, and patient dose delivery. Verification of each step in the treatment process is assumed to result in correct dose delivery to the patient. ^ The accelerator MC model was verified against commissioning data for field sizes from 0.8 × 0.8 cm 2 to 10 × 10 cm 2. Depth doses were within 2% local percent difference (LPD) in low gradient regions and 1 mm distance to agreement (DTA) in high gradient regions. Lateral profiles were within 2% LPD in low gradient regions and 1 mm DTA in high gradient regions. Calculated output factors were within 1% of measurement for field sizes ≥1 × 1 cm2. ^ The measured and calculated pretreatment EPID dose patterns were compared using criteria of 5% LPD, 1 mm DTA, or 2% of central axis pixel value with ≥95% of compared points required to pass for successful verification. Pretreatment field verification resulted in 97% percent of the points passing. ^ The RTPS and Monte Carlo phantom dose calculations were compared using 5% LPD, 2 mm DTA, or 2% of the maximum dose with ≥95% of compared points required passing for successful verification. RTPS calculation verification resulted in 97% percent of the points passing. ^ The measured and calculated EPID exit dose patterns were compared using criteria of 5% LPD, 1 mm DTA, or 2% of central axis pixel value with ≥95% of compared points required to pass for successful verification. Exit dose verification resulted in 97% percent of the points passing. ^ Each of the processes above verified an individual step in the treatment planning and delivery process. The combination of these verification steps ensures accurate treatment delivery to the patient. This work shows that Monte Carlo calculations and EPID dosimetry can be used to quantitatively verify IMSRT treatments resulting in improved patient care and, potentially, improved clinical outcome. ^

Monte Carlo and analytical dose calculations for ocular proton therapy

Uveal melanoma is a rare but life-threatening form of ocular cancer. Contemporary treatment techniques include proton therapy, which enables conservation of the eye and its useful vision. Dose to the proximal structures is widely believed to play a role in treatment side effects, therefore, reliable dose estimates are required for properly evaluating the therapeutic value and complication risk of treatment plans. Unfortunately, current simplistic dose calculation algorithms can result in errors of up to 30% in the proximal region. In addition, they lack predictive methods for absolute dose per monitor unit (D/MU) values. ^ To facilitate more accurate dose predictions, a Monte Carlo model of an ocular proton nozzle was created and benchmarked against measured dose profiles to within ±3% or ±0.5 mm and D/MU values to within ±3%. The benchmarked Monte Carlo model was used to develop and validate a new broad beam dose algorithm that included the influence of edgescattered protons on the cross-field intensity profile, the effect of energy straggling in the distal portion of poly-energetic beams, and the proton fluence loss as a function of residual range. Generally, the analytical algorithm predicted relative dose distributions that were within ±3% or ±0.5 mm and absolute D/MU values that were within ±3% of Monte Carlo calculations. Slightly larger dose differences were observed at depths less than 7 mm, an effect attributed to the dose contributions of edge-scattered protons. Additional comparisons of Monte Carlo and broad beam dose predictions were made in a detailed eye model developed in this work, with generally similar findings. ^ Monte Carlo was shown to be an excellent predictor of the measured dose profiles and D/MU values and a valuable tool for developing and validating a broad beam dose algorithm for ocular proton therapy. The more detailed physics modeling by the Monte Carlo and broad beam dose algorithms represent an improvement in the accuracy of relative dose predictions over current techniques, and they provide absolute dose predictions. It is anticipated these improvements can be used to develop treatment strategies that reduce the incidence or severity of treatment complications by sparing normal tissue. ^

Pencil -beam redefinition algorithm dose calculations for electron therapy treatment planning

The electron pencil-beam redefinition algorithm (PBRA) of Shiu and Hogstrom has been developed for use in radiotherapy treatment planning (RTP). Earlier studies of Boyd and Hogstrom showed that the PBRA lacked an adequate incident beam model, that PBRA might require improved electron physics, and that no data existed which allowed adequate assessment of the PBRA-calculated dose accuracy in a heterogeneous medium such as one presented by patient anatomy. The hypothesis of this research was that by addressing the above issues the PBRA-calculated dose would be accurate to within 4% or 2 mm in regions of high dose gradients. A secondary electron source was added to the PBRA to account for collimation-scattered electrons in the incident beam. Parameters of the dual-source model were determined from a minimal data set to allow ease of beam commissioning. Comparisons with measured data showed 3% or better dose accuracy in water within the field for cases where 4% accuracy was not previously achievable. A measured data set was developed that allowed an evaluation of PBRA in regions distal to localized heterogeneities. Geometries in the data set included irregular surfaces and high- and low-density internal heterogeneities. The data was estimated to have 1% precision and 2% agreement with accurate, benchmarked Monte Carlo (MC) code. PBRA electron transport was enhanced by modeling local pencil beam divergence. This required fundamental changes to the mathematics of electron transport (divPBRA). Evaluation of divPBRA with the measured data set showed marginal improvement in dose accuracy when compared to PBRA; however, 4% or 2mm accuracy was not achieved by either PBRA version for all data points. Finally, PBRA was evaluated clinically by comparing PBRA- and MC-calculated dose distributions using site-specific patient RTP data. Results show PBRA did not agree with MC to within 4% or 2mm in a small fraction (<3%) of the irradiated volume. Although the hypothesis of the research was shown to be false, the minor dose inaccuracies should have little or no impact on RTP decisions or patient outcome. Therefore, given ease of beam commissioning, documentation of accuracy, and calculational speed, the PBRA should be considered a practical tool for clinical use. ^

Immobilisation dose rates and photogrammetry based mass calculations for southern elephant seals during expedition JUB2010

Immobilization and anaesthesia of adult male southern elephant seals (Mirounga leonina) is potentially risky for animals and scientists. A tiletamine/zolazepam injection is considered the most appropriate drug combination for field application in this species. Since appropriate dosages are difficult to assess due to uncertainties in weight estimation, we used photogrammetry-derived weight estimates to ensure precise post hoc calculations of dosages. We report on 15 intramuscular tiletamine/zolazepam immobilizations of post-moult males of the upper weight class at King George Island/Isla 25 de Mayo, in April 2010. Initial injections were made using blowpipe syringes. Mean tiletamine/zolazepam combined dosages of 0.71 mg/kg (SD ± 0.16) ranged between 0.46 and 1.01 mg/kg. In four cases, ketamine was added in dosages between 0.96 and 2.61 mg/kg. Mean induction period was 23 min (± 15), and the mean duration of the procedures from first injection to release of the animals required 96 min (± 51). Four seals exhibited periods of apnoea, and one case of an extended, repetitive, and potentially critical apnoea (> 25 and 8 min) required intervention in order to successfully re-initiate spontaneous respiration. All procedures resulted in proper immobilizations allowing for the deployment of the satellite tags on the seals' heads. The fact that even substantial deviations between the initial weight estimates and the photogrammetry-derived weight estimates had no apparent effect on the course of the immobilization underlines the drugs' wide safety margin in this species.

Sedimentation rate calculations and isotopic analysis on sediment cores along the IODP Expedition 311 transect

The oceanographic and tectonic conditions of accretionary margins are well-suited for several potential processes governing methane generation, storage and release. To identify the relevant methane evolution pathways in the northern Cascadia accretionary margin, a four-site transect was drilled during Integrated Ocean Drilling Program Expedition 311. The d13C values of methane range from a minimum value of -82.2 per mil on an uplifted ridge of accreted sediment near the deformation front (Site U1326, 1829 mbsl, meters below sea level) to a maximum value of -39.5 per mil at the most landward location within an area of steep canyons near the shelf edge (Site U1329, 946 mbsl). An interpretation based solely on methane isotope values might conclude the 13C-enrichment of methane indicates a transition from microbially- to thermogenically-sourced methane. However, the co-existing CO2 exhibits a similar trend of 13C-enrichment along the transect with values ranging from -22.5 per mil to +25.7 per mil. The magnitude of the carbon isotope separation between methane and CO2 (Ec = 63.8 ± 5.8) is consistent with isotope fractionation during microbially mediated carbonate reduction. These results, in conjunction with a transect-wide gaseous hydrocarbon content composed of > 99.8% (by volume) methane and uniform dDCH4 values (-172 per mil ± 8) that are distinct from thermogenic methane at a seep located 60 km from the Expedition 311 transect, suggest microbial CO2 reduction is the predominant methane source at all investigated sites. The magnitude of the intra-site downhole 13C-enrichment of CO2 within the accreted ridge (Site U1326) and a slope basin nearest the deformation front (Site U1325, 2195 mbsl) is ~ 5 per mil. At the mid-slope site (Site U1327, 1304 mbsl) the downhole 13C-enrichment of the CO2 is ~ 25 per mil and increases to ~ 40 per mil at the near-shelf edge Site U1329. This isotope fractionation pattern is indicative of more extensive diagenetic alteration at sites with greater 13C-enrichment. The magnitude of the 13C-enrichment of CO2 correlates with decreasing sedimentation rates and a diminishing occurrence of stratigraphic gas hydrate. We suggest the decreasing sedimentation rates increase the exposure time of sedimentary organic matter to aerobic and anaerobic degradation, during burial, thereby reducing the availability of metabolizable organic matter available for methane production. This process is reflected in the occurrence and distribution of gas hydrate within the northern Cascadia margin accretionary prism. Our observations are relevant for evaluating methane production and the occurrence of stratigraphic gas hydrate within other convergent margins.

(Table 1) CaCO3, terrigenous matter and Corg calculations from DSDP Hole 75-532B

Samples taken at 10 cm intervals from DSDP Core 532B-17 contain variations in carbonate, opal, organic carbon, and terrigenous components that correlate with light-dark cycles in sediment color. The core site, at 1300 m water depth, is well above the CCD, yet the color variations appear to result largely from cyclical fluctuations in carbonate dissolution, which was greater during glacial periods. Higher concentrations of organic carbon and of terrigenous sediment components correlate with enhanced carbonate dissolution, but opal concentrations inversely correlate and suggest that biological productivity at this site diminished during glacial periods. A complicated glacial-interglacial picture emerges from the data. In interglacial times, upwelling associated with the Benguela Current produced abundant opaline material, organic matter was fairly well preserved, and carbonate was only moderately dissolved. In glacial times, the upwelling core shifted as sea level fell and winds intensified. Productivity in the waters over Site 532 decreased, but lateral supply of oxidizable organic matter enhanced carbonate dissolution, giving rise to light-dark cycles in these sediments.

Paleo-sea surface temperature calculations in the equatorial east Atlantic

We present two ~270 kyr paleo-sea surface temperature (SST) records from the Equatorial Divergence and the South Equatorial Current derived from Mg/Ca ratios in the planktic foraminifer Globigerinoides sacculifer. The present study suggests that the magnesium signature of G. sacculifer provides a seasonal SST estimate from the upper ~50 m of the water column generated during upwelling in austral low-latitude fall/winter. Common to both down-core records is a glacial-interglacial amplitude of ~3°-3.5°C for the last climatic changes and lower Holocene and glacial oxygen isotope stage 2 temperatures compared with interglacial stage 5.5 and glacial stage 6 temperatures, respectively. The comparison to published SST estimates from alkenones, oxygen isotopes, and foraminiferal transfer function from the same core material pinpoints discrepancies and conformities between methods.