Benchmarking and validation of a Geant4–SHADOW Monte Carlo simulation for dose calculations in microbeam radiation therapy

Microbeam radiation therapy (MRT) is a synchrotron‐based radiotherapy modality that uses high‐intensity beams of spatially fractionated radiation to treat tumours. The rapid evolution of MRT towards clinical trials demands accurate treatment planning systems (TPS), as well as independent tools for the verification of TPS calculated dose distributions in order to ensure patient safety and treatment efficacy. Monte Carlo computer simulation represents the most accurate method of dose calculation in patient geometries and is best suited for the purpose of TPS verification. A Monte Carlo model of the ID17 biomedical beamline at the European Synchrotron Radiation Facility has been developed, including recent modifications, using the Geant4 Monte Carlo toolkit interfaced with the SHADOW X‐ray optics and ray‐tracing libraries. The code was benchmarked by simulating dose profiles in water‐equivalent phantoms subject to irradiation by broad‐beam (without spatial fractionation) and microbeam (with spatial fractionation) fields, and comparing against those calculated with a previous model of the beamline developed using the PENELOPE code. Validation against additional experimental dose profiles in water‐equivalent phantoms subject to broad‐beam irradiation was also performed. Good agreement between codes was observed, with the exception of out‐of‐field doses and toward the field edge for larger field sizes. Microbeam results showed good agreement between both codes and experimental results within uncertainties. Results of the experimental validation showed agreement for different beamline configurations. The asymmetry in the out‐of‐field dose profiles due to polarization effects was also investigated, yielding important information for the treatment planning process in MRT. This work represents an important step in the development of a Monte Carlo‐based independent verification tool for treatment planning in MRT.     

Managing longevity and disability risks in life annuities with long term care

The aim of the paper is twofold. Firstly, it develops a model for risk assessment in a portfolio of life annuities with long term care benefits. These products are usually represented by a Markovian Multi-State model and are affected by both longevity and disability risks. Here, a stochastic projection model is proposed in order to represent the future evolution of mortality and disability transition intensities. Data from the Italian National Institute of Social Security (INPS) and from Human Mortality Database (HMD) are used to estimate the model parameters. Secondly, it investigates the solvency in a portfolio of enhanced pensions. To this aim a risk model based on the portfolio risk reserve is proposed and different rules to calculate solvency capital requirements for life underwriting risk are examined. Such rules are then compared with the standard formula proposed by the Solvency II project.     

Automated detection and localization of bowhead whale sounds in the presence of seismic airgun surveys

An automated procedure has been developed for detecting and localizing frequency-modulated bowhead whale sounds in the presence of seismic airgun surveys. The procedure was applied to four years of data, collected from over 30 directional autonomous recording packages deployed over a 280 km span of continental shelf in the Alaskan Beaufort Sea. The procedure has six sequential stages that begin by extracting 25-element feature vectors from spectrograms of potential call candidates. Two cascaded neural networks then classify some feature vectors as bowhead calls, and the procedure then matches calls between recorders to triangulate locations. To train the networks, manual analysts flagged 219 471 bowhead call examples from 2008 and 2009. Manual analyses were also used to identify 1.17 million transient signals that were not whale calls. The network output thresholds were adjusted to reject 20% of whale calls in the training data. Validation runs using 2007 and 2010 data found that the procedure missed 30%-40% of manually detected calls. Furthermore, 20%-40% of the sounds flagged as calls are not present in the manual analyses; however, these extra detections incorporate legitimate whale calls overlooked by human analysts. Both manual and automated methods produce similar spatial and temporal call distributions.     

Multipulse Phases in k-Mixtures of Bose–Einstein Condensates

For the system

$-\Delta U_i+ U_i=U_i^3-\beta U_i\sum_{j\neq i}U_j^2,\quad i=1,\dots,k,$

(with k ≧ 3), we prove the existence for β large of positive radial solutions on \({\mathbb R^N}\) . We show that as β →  + ∞, the profile of each component U _i separates, in many pulses, from the others. Moreover, we can prescribe the location of such pulses in terms of the oscillations of the changing-sign solutions of the scalar equation  ? ΔW  +  W  =  W³. Within an Hartree–Fock approximation, this provides a theoretical indication of phase separation into many nodal domains for the k-mixtures of Bose–Einstein condensates.

     

Radial Solutions of¶Superlinear Equations on ℝN¶Part II: The Forced Case

We study the effect of a forcing term in the context of the search of multiple nodal solutions u∈h ¹(? ^N ) to a class of elliptic equations of type¶?Δu(x)=f(|x|,u(x))+h(|x|), x∈? ^N ,¶where f(|x|≡0 and f is superlinear but subcritical at infinity.     

An age‐dependent population equation with delayed birth process

Within a semigroup framework, we discuss well posedness and qualitative behaviour of an age‐dependent population equation with delay in the birth process. Using positivity and Perron–Frobenius theory we obtain an explicit stability criterion. Copyright © 2004 John Wiley & Sons, Ltd.     

Hierarchical decomposition of domains with fractures

We consider the efficient and robust numerical solution of elliptic problems with jumping coefficients occurring on a network of thin fractures. We present an iterative solution concept based on a hierarchical separation of the fractures and the surrounding rock matrix. Upper estimates for the convergence rates are independent of the width of the fractures and of the jumps of the coefficients. Inexact solution of the local subproblems is also considered. The theoretical results are illustrated by numerical experiments.

     

New strategies for the synthesis of partially fluorinated acrylic polymers as possible materials for the protection of stone monuments

The conservation of the cultural heritage requires the development of new materials having specific characteristics that encompass particular attention to durability and efficacy. We approached this problem by synthesizing polyacrylic esters containing variable amounts of fluorine in the α-position of the main chain. These products were obtained from the copolymerisation of ammonium 2-fluoroacrylate and acrylic acid. The polyacrylic acids were esterified using different procedures. The polyester characteristics vary in relation to the polymerisation procedures and degree of esterification. The best esterification results, were obtained using a reaction catalyzed by BF₃ or TMSCl. These materials show good properties and are of potential interest for their use as protective agents for stone conservation.