Monte Carlo modeling of polarized light propagation: Stokes vs Jones - Part I

This bipartite comparative study aims at inspecting the similarities and differences between the Jones and Stokes–Mueller formalisms when modeling polarized light propagation with numerical simulations of the Monte Carlo type. In this first part, we review the theoretical concepts that concern light propagation and detection with both pure and partially/totally unpolarized states. The latter case involving fluctuations, or “depolarizing effects,” is of special interest here: Jones and Stokes–Mueller are equally apt to model such effects and are expected to yield identical results. In a second, ensuing paper, empirical evidence is provided by means of numerical experiments, using both formalisms.

Validation of the Swiss Monte Carlo Plan for a static and dynamic 6 MV photon beam

Monte Carlo (MC) based dose calculations can compute dose distributions with an accuracy surpassing that of conventional algorithms used in radiotherapy, especially in regions of tissue inhomogeneities and surface discontinuities. The Swiss Monte Carlo Plan (SMCP) is a GUI-based framework for photon MC treatment planning (MCTP) interfaced to the Eclipse treatment planning system (TPS). As for any dose calculation algorithm, also the MCTP needs to be commissioned and validated before using the algorithm for clinical cases. Aim of this study is the investigation of a 6 MV beam for clinical situations within the framework of the SMCP. In this respect, all parts i.e. open fields and all the clinically available beam modifiers have to be configured so that the calculated dose distributions match the corresponding measurements. Dose distributions for the 6 MV beam were simulated in a water phantom using a phase space source above the beam modifiers. The VMC++ code was used for the radiation transport through the beam modifiers (jaws, wedges, block and multileaf collimator (MLC)) as well as for the calculation of the dose distributions within the phantom. The voxel size of the dose distributions was 2mm in all directions. The statistical uncertainty of the calculated dose distributions was below 0.4%. Simulated depth dose curves and dose profiles in terms of [Gy/MU] for static and dynamic fields were compared with the corresponding measurements using dose difference and γ analysis. For the dose difference criterion of ±1% of D(max) and the distance to agreement criterion of ±1 mm, the γ analysis showed an excellent agreement between measurements and simulations for all static open and MLC fields. The tuning of the density and the thickness for all hard wedges lead to an agreement with the corresponding measurements within 1% or 1mm. Similar results have been achieved for the block. For the validation of the tuned hard wedges, a very good agreement between calculated and measured dose distributions was achieved using a 1%/1mm criteria for the γ analysis. The calculated dose distributions of the enhanced dynamic wedges (10°, 15°, 20°, 25°, 30°, 45° and 60°) met the criteria of 1%/1mm when compared with the measurements for all situations considered. For the IMRT fields all compared measured dose values agreed with the calculated dose values within a 2% dose difference or within 1 mm distance. The SMCP has been successfully validated for a static and dynamic 6 MV photon beam, thus resulting in accurate dose calculations suitable for applications in clinical cases.

Phase I clinical and pharmacokinetic evaluation of the vascular-disrupting agent OXi4503 in patients with advanced solid tumors

Preclinical studies show that OXi4503 (combretastatin A1 diphosphate, CA1P) is more potent than other clinically evaluated vascular-disrupting agents.

Ankle arthrodesis versus total ankle replacement: how do I decide?

Taking into account numerous individual criteria, the correct indication substantially influences the outcome of patients with end-stage ankle arthritis treated by ankle arthrodesis or total ankle replacement. The purpose of this report is to assist the foot and ankle surgeon or orthopedic surgeon involved in choosing ankle arthrodesis or total ankle replacement in decision-making. Balancing the criteria that are discussed in consideration of the recent relevant literature and evidence available, the surgeon is directed to the correct individual decision.

Phase I trial of combretastatin A4 phosphate (CA4P) in combination with bevacizumab in patients with advanced cancer

The vascular disrupting agent (VDA) combretastatin A4 phosphate (CA4P) induces significant tumor necrosis as a single agent. Preclinical models have shown that the addition of an anti-VEGF antibody to a VDA attenuates the revascularization of the surviving tumor rim and thus significantly increases antitumor activity.

More than just Good Luck. How to Explain the Success of Monte Carlo Simulations

Monte Carlo simulation is a powerful method in many natural and social sciences. But what sort of method is it? And where does its power come from? Are Monte Carlo simulations experiments, theories or something else? The aim of this talk is to answer these questions and to explain the power of Monte Carlo simulations. I provide a classification of Monte Carlo techniques and defend the claim that Monte Carlo simulation is a sort of inference.

Monte Carlo modeling of polarized light propagation: Stokes vs Jones - Part II

In this second part of our comparative study inspecting the (dis)similarities between “Stokes” and “Jones,” we present simulation results yielded by two independent Monte Carlo programs: (i) one developed in Bern with the Jones formalism and (ii) the other implemented in Ulm with the Stokes notation. The simulated polarimetric experiments involve suspensions of polystyrene spheres with varying size. Reflection and refraction at the sample/air interfaces are also considered. Both programs yield identical results when propagating pure polarization states, yet, with unpolarized illumination, second order statistical differences appear, thereby highlighting the pre-averaged nature of the Stokes parameters. This study serves as a validation for both programs and clarifies the misleading belief according to which “Jones cannot treat depolarizing effects.”

Minimum carbon payment along an aridity gradient for dryland afforestation in non annex I countries with stochastic weather and prices

Recent findings demonstrate that trees in deserts are efficient carbon sinks. It remains however unknown whether the Clean Development Mechanism will accelerate the planting of trees in Non Annex I dryland countries. We estimated the price of carbon at which a farmer would be indifferent between his customary activity and the planting of trees to trade carbon credits, along an aridity gradient. Carbon yields were simulated by means of the CO2FIX v3.1 model for Pinus halepensis with its respective yield classes along the gradient (Arid – 100mm to Dry Sub Humid conditions – 900mm). Wheat and pasture yields were predicted on somewhat similar nitrogen-based quadratic models, using 30 years of weather data to simulate moisture stress. Stochastic production, input and output prices were afterwards simulated on a Monte Carlo matrix. Results show that, despite the high levels of carbon uptake, carbon trading by afforesting is unprofitable anywhere along the gradient. Indeed, the price of carbon would have to raise unrealistically high, and the certification costs would have to drop significantly, to make the Clean Development Mechanism worthwhile for non annex I dryland countries farmers. From a government agency's point of view the Clean Development Mechanism is attractive. However, such agencies will find it difficult to demonstrate “additionality”, even if the rule may be somewhat flexible. Based on these findings, we will further discuss why the Clean Development Mechanism, a supposedly pro-poor instrument, fails to assist farmers in Non Annex I dryland countries living at minimum subsistence level.

Stent-Retriever Thrombectomy after Intravenous t-PA vs. t-PA Alone in Stroke

Background Among patients with acute ischemic stroke due to occlusions in the proximal anterior intracranial circulation, less than 40% regain functional independence when treated with intravenous tissue plasminogen activator (t-PA) alone. Thrombectomy with the use of a stent retriever, in addition to intravenous t-PA, increases reperfusion rates and may improve long-term functional outcome. Methods We randomly assigned eligible patients with stroke who were receiving or had received intravenous t-PA to continue with t-PA alone (control group) or to undergo endovascular thrombectomy with the use of a stent retriever within 6 hours after symptom onset (intervention group). Patients had confirmed occlusions in the proximal anterior intracranial circulation and an absence of large ischemic-core lesions. The primary outcome was the severity of global disability at 90 days, as assessed by means of the modified Rankin scale (with scores ranging from 0 [no symptoms] to 6 [death]). Results The study was stopped early because of efficacy. At 39 centers, 196 patients underwent randomization (98 patients in each group). In the intervention group, the median time from qualifying imaging to groin puncture was 57 minutes, and the rate of substantial reperfusion at the end of the procedure was 88%. Thrombectomy with the stent retriever plus intravenous t-PA reduced disability at 90 days over the entire range of scores on the modified Rankin scale (P<0.001). The rate of functional independence (modified Rankin scale score, 0 to 2) was higher in the intervention group than in the control group (60% vs. 35%, P<0.001). There were no significant between-group differences in 90-day mortality (9% vs. 12%, P=0.50) or symptomatic intracranial hemorrhage (0% vs. 3%, P=0.12). Conclusions In patients receiving intravenous t-PA for acute ischemic stroke due to occlusions in the proximal anterior intracranial circulation, thrombectomy with a stent retriever within 6 hours after onset improved functional outcomes at 90 days.

3D Direct Simulation Monte Carlo Modelling of the Inner Gas Coma of Comet 67P/Churyumov-Gerasimenko: A Parameter Study

Direct Simulation Monte Carlo (DSMC) is a powerful numerical method to study rarefied gas flows such as cometary comae and has been used by several authors over the past decade to study cometary outflow. However, the investigation of the parameter space in simulations can be time consuming since 3D DSMC is computationally highly intensive. For the target of ESA's Rosetta mission, comet 67P/Churyumov-Gerasimenko, we have identified to what extent modification of several parameters influence the 3D flow and gas temperature fields and have attempted to establish the reliability of inferences about the initial conditions from in situ and remote sensing measurements. A large number of DSMC runs have been completed with varying input parameters. In this work, we present the simulation results and conclude on the sensitivity of solutions to certain inputs. It is found that among cases of water outgassing, the surface production rate distribution is the most influential variable to the flow field.