Planning system for the optimization of electric field delivery using implanted electrodes for brain tumor control

BACKGROUND: The use of non-ionizing electric fields from low-intensity voltage sources (\u3c 10 V) to control malignant tumor growth is showing increasing potential as a cancer treatment modality. A method of applying these low-intensity electric fields using multiple implanted electrodes within or adjacent to tumor volumes has been termed as intratumoral modulation therapy (IMT). PURPOSE: This study explores advancements in the previously established IMT optimization algorithm, and the development of a custom treatment planning system for patient-specific IMT. The practicality of the treatment planning system is demonstrated by implementing the full optimization pipeline on a brain phantom with robotic electrode implantation, postoperative imaging, and treatment stimulation. METHODS: The integrated planning pipeline in 3D Slicer begins with importing and segmenting patient magnetic resonance images (MRI) or computed tomography (CT) images. The segmentation process is manual, followed by a semi-automatic smoothing step that allows the segmented brain and tumor mesh volumes to be smoothed and simplified by applying selected filters. Electrode trajectories are planned manually on the patient MRI or CT by selecting insertion and tip coordinates for a chosen number of electrodes. The electrode tip positions and stimulation parameters (phase shift and voltage) can then be optimized with the custom semi-automatic IMT optimization algorithm where users can select the prescription electric field, voltage amplitude limit, tissue electrical properties, nearby organs at risk, optimization parameters (electrode tip location, individual contact phase shift and voltage), desired field coverage percent, and field conformity optimization. Tables of optimization results are displayed, and the resulting electric field is visualized as a field-map superimposed on the MR or CT image, with 3D renderings of the brain, tumor, and electrodes. Optimized electrode coordinates are transferred to robotic electrode implantation software to enable planning and subsequent implantation of the electrodes at the desired trajectories. RESULTS: An IMT treatment planning system was developed that incorporates patient-specific MRI or CT, segmentation, volume smoothing, electrode trajectory planning, electrode tip location and stimulation parameter optimization, and results visualization. All previous manual pipeline steps operating on diverse software platforms were coalesced into a single semi-automated 3D Slicer-based user interface. Brain phantom validation of the full system implementation was successful in preoperative planning, robotic electrode implantation, and postoperative treatment planning to adjust stimulation parameters based on actual implant locations. Voltage measurements were obtained in the brain phantom to determine the electrical parameters of the phantom and validate the simulated electric field distribution. CONCLUSIONS: A custom treatment planning and implantation system for IMT has been developed in this study and validated on a phantom brain model, providing an essential step in advancing IMT technology toward future clinical safety and efficacy investigations

Coherent instabilities in a semiconductor laser with fast gain recovery

We report the observation of a coherent multimode instability in quantum cascade lasers (QCLs), which is driven by the same fundamental mechanism of Rabi oscillations as the elusive Risken-Nummedal-Graham-Haken (RNGH) instability predicted 40 years ago for ring lasers. The threshold of the observed instability is significantly lower than in the original RNGH instability, which we attribute to saturable-absorption nonlinearity in the laser. Coherent effects, which cannot be reproduced by standard laser rate equations, can play therefore a key role in the multimode dynamics of QCLs, and in lasers with fast gain recovery in general.Comment: 5 pages, 4 figure

An Exact Algorithm for Side-Chain Placement in Protein Design

Computational protein design aims at constructing novel or improved functions on the structure of a given protein backbone and has important applications in the pharmaceutical and biotechnical industry. The underlying combinatorial side-chain placement problem consists of choosing a side-chain placement for each residue position such that the resulting overall energy is minimum. The choice of the side-chain then also determines the amino acid for this position. Many algorithms for this NP-hard problem have been proposed in the context of homology modeling, which, however, reach their limits when faced with large protein design instances. In this paper, we propose a new exact method for the side-chain placement problem that works well even for large instance sizes as they appear in protein design. Our main contribution is a dedicated branch-and-bound algorithm that combines tight upper and lower bounds resulting from a novel Lagrangian relaxation approach for side-chain placement. Our experimental results show that our method outperforms alternative state-of-the art exact approaches and makes it possible to optimally solve large protein design instances routinely

Multipolarity of quasicontinuum γ-rays from collective high-spin states in 152Dy

Measured internal conversion coefficients for quasicontinuum transitions in 152Dy in the spin range of 30–50 establish their predominantly stretched E2 character. Those transitions are attributed to triaxial bands near the yrast line as calculated in terms of the cranking approximation using the Woods-Saxon potential

Typical equilibrium state of an embedded quantum system

We consider an arbitrary quantum system coupled non perturbatively to a large arbitrary and fully quantum environment. In [G. Ithier and F. Benaych-Georges, Phys. Rev. A 96, 012108 (2017)] the typicality of the dynamics of such an embedded quantum system was established for several classes of random interactions. In other words, the time evolution of its quantum state does not depend on the microscopic details of the interaction. Focusing at the long time regime, we use this property to calculate analytically a new partition function characterizing the stationary state and involving the overlaps between eigenvectors of a bare and a dressed Hamiltonian. This partition function provides a new thermodynamical ensemble which includes the microcanonical and canonical ensembles as particular cases. We check our predictions with numerical simulations.Comment: 1 figure, 5 pages. This article supersedes the part on the equilibrium state in arXiv:1510.0435

Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

In this article, we describe an approach for solving partial differential equations with general boundary conditions imposed on arbitrarily shaped boundaries. A continuous function, the domain parameter, is used to modify the original differential equations such that the equations are solved in the region where a domain parameter takes a specified value while boundary conditions are imposed on the region where the value of the domain parameter varies smoothly across a short distance. The mathematical derivations are straightforward and generically applicable to a wide variety of partial differential equations. To demonstrate the general applicability of the approach, we provide four examples herein: (1) the diffusion equation with both Neumann and Dirichlet boundary conditions; (2) the diffusion equation with both surface diffusion and reaction; (3) the mechanical equilibrium equation; and (4) the equation for phase transformation with the presence of additional boundaries. The solutions for several of these cases are validated against corresponding analytical and semi-analytical solutions. The potential of the approach is demonstrated with five applications: surface-reaction-diffusion kinetics with a complex geometry, Kirkendall-effect-induced deformation, thermal stress in a complex geometry, phase transformations affected by substrate surfaces, and a self-propelled droplet.Comment: This document is the revised version of arXiv:0912.1288v

Semi-Hard Scattering Unraveled from Collective Dynamics by Two-Pion Azimuthal Correlations in 158 A GeV/c Pb + Au Collisions

Elliptic flow and two-particle azimuthal correlations of charged hadrons and high-$p_T$ pions ($p_T>$ 1 GeV/$c$) have been measured close to mid-rapidity in 158A GeV/$c$ Pb+Au collisions by the CERES experiment. Elliptic flow ($v_2$) rises linearly with $p_T$ to a value of about 10% at 2 GeV/$c$. Beyond $p_T\approx$ 1.5 GeV/$c$, the slope decreases considerably, possibly indicating a saturation of $v_2$ at high $p_T$. Two-pion azimuthal anisotropies for $p_T>$ 1.2 GeV/$c$ exceed the elliptic flow values by about 60% in mid-central collisions. These non-flow contributions are attributed to near-side and back-to-back jet-like correlations, the latter exhibiting centrality dependent broadening.Comment: Submitted to Phys. Rev. Letters, 4 pages, 5 figure

Recent results from Pb-Au collisions at 158 GeV/c per nucleon obtained with the CERES spectrometer

During the 1996 lead run time, CERES has accumulated 42 million events, corresponding to a factor of 5 more statistics than in 1995 and 2.5 million events of a special photon-run. We report on the results of the low-mass e$^+$e$^-$-pair analysis. Since the most critical item is the poor signal-to-background ratio we also discuss the understanding of this background, in absolute terms, with the help of a detailed Monte Carlo simulation. We show preliminary results of the photon analysis and summarize the results of the hadron analysis preliminarily reported on already at QM'97Comment: 10 pages, 9 figures, Proceedings of the XIV Int. Conf. on Nucleus-Nucleus Collisions,Quark Matter 99, Torino, Italy, May 10 - 15, 199

Low-mass e+e- pair production in 158 A GeV Pb-Au collisions at the CERN SPS, its dependence on multiplicity and transverse momentum

We report a measurement of low-mass electron pairs observed in 158 GeV/nucleon Pb-Au collisions. The pair yield integrated over the range of invariant masses 0.2 < m < 2.0 GeV is enhanced by a factor of 3.5 +/- 0.4 (stat) +/- 0.9 (syst) over the expectation from neutral meson decays. As observed previously in S-Au collisions, the enhancement is most pronounced in the invariant-mass region 300-700 MeV. For Pb-Au we find evidence for a strong increase of the enhancement with centrality. In addition, we show that the enhancement covers a wide range in transverse momentum, but is largest at the lowest observed pt.Comment: 17 pages, 4 figures, submitted to Phys.Lett.