Predicting dynamic stability from static features in power grid models using machine learning.

A reliable supply with electric power is vital for our society. Transmission line failures are among the biggest threats for power grid stability as they may lead to a splitting of the grid into mutual asynchronous fragments. New conceptual methods are needed to assess system stability that complement existing simulation models. In this article, we propose a combination of network science metrics and machine learning models to predict the risk of desynchronization events. Network science provides metrics for essential properties of transmission lines such as their redundancy or centrality. Machine learning models perform inherent feature selection and, thus, reveal key factors that determine network robustness and vulnerability. As a case study, we train and test such models on simulated data from several synthetic test grids. We find that the integrated models are capable of predicting desynchronization events after line failures with an average precision greater than 0.996 when averaging over all datasets. Learning transfer between different datasets is generally possible, at a slight loss of prediction performance. Our results suggest that power grid desynchronization is essentially governed by only a few network metrics that quantify the networks' ability to reroute the flow without creating exceedingly high static line loadings.

Optimizing the geometry of transportation networks in the presence of congestion.

Urban transport systems are gaining in importance, as an increasing share of the global population lives in cities and mobility-based carbon emissions must be reduced to mitigate climate change and improve air quality and citizens' health. As a result, public transport systems are prone to congestion, raising the question of how to optimize them to cope with this challenge. In this paper, we analyze the optimal design of urban transport networks to minimize the average travel time in monocentric as well as in polycentric cities. We suggest an elementary model for congestion and introduce a numerical method to determine the optimal shape among a set of predefined geometries considering different models for the behavior of individual travelers. We map out the optimal shape of fundamental network geometries with a focus on the impact of congestion.

Dual communities in spatial networks.

Both human-made and natural supply systems, such as power grids and leaf venation networks, are built to operate reliably under changing external conditions. Many of these spatial networks exhibit community structures. Here, we show that a relatively strong connectivity between the parts of a network can be used to define a different class of communities: dual communities. We demonstrate that traditional and dual communities emerge naturally as two different phases of optimized network structures that are shaped by fluctuations and that they both suppress failure spreading, which underlines their importance in understanding the shape of real-world supply networks.

Network isolators inhibit failure spreading in complex networks.

In our daily lives, we rely on the proper functioning of supply networks, from power grids to water transmission systems. A single failure in these critical infrastructures can lead to a complete collapse through a cascading failure mechanism. Counteracting strategies are thus heavily sought after. In this article, we introduce a general framework to analyse the spreading of failures in complex networks and demostrate that not only decreasing but also increasing the connectivity of the network can be an effective method to contain damages. We rigorously prove the existence of certain subgraphs, called network isolators, that can completely inhibit any failure spreading, and we show how to create such isolators in synthetic and real-world networks. The addition of selected links can thus prevent large scale outages as demonstrated for power transmission grids.

Discontinuous transition to loop formation in optimal supply networks.

The structure and design of optimal supply networks is an important topic in complex networks research. A fundamental trait of natural and man-made networks is the emergence of loops and the trade-off governing their formation: adding redundant edges to supply networks is costly, yet beneficial for resilience. Loops typically form when costs for new edges are small or inputs uncertain. Here, we shed further light on the transition to loop formation. We demonstrate that loops emerge discontinuously when decreasing the costs for new edges for both an edge-damage model and a fluctuating sink model. Mathematically, new loops are shown to form through a saddle-node bifurcation. Our analysis allows to heuristically predict the location and cost where the first loop emerges. Finally, we unveil an intimate relationship among betweenness measures and optimal tree networks. Our results can be used to understand the evolution of loop formation in real-world biological networks.

Order parameter allows classification of planar graphs based on balanced fixed points in the Kuramoto model.

Phase balanced states are a highly underexplored class of solutions of the Kuramoto model and other coupled oscillator models on networks. So far, coupled oscillator research focused on phase synchronized solutions. Yet, global constraints on oscillators may forbid synchronized state, rendering phase balanced states as the relevant stable state. If, for example, oscillators are driving the contractions of a fluid filled volume, conservation of fluid volume constrains oscillators to balanced states as characterized by a vanishing Kuramoto order parameter. It has previously been shown that stable, balanced patterns in the Kuramoto model exist on circulant graphs. However, which noncirculant graphs first of all allow for balanced states and what characterizes the balanced states is unknown. Here, we derive rules of how to build noncirculant, planar graphs allowing for balanced states from the simple cycle graph by adding loops or edges to it. We thereby identify different classes of small planar networks allowing for balanced states. Investigating the balanced states' characteristics, we find that the variance in basin stability scales linearly with the size of the graph for these networks. We introduce the balancing ratio as an order parameter based on the basin stability approach to classify balanced states on networks and evaluate it analytically for a subset of the network classes. Our results offer an analytical description of noncirculant graphs supporting stable, balanced states and may thereby help to understand the topological requirements on oscillator networks under global constraints.

Multistability in lossy power grids and oscillator networks.

Networks of phase oscillators are studied in various contexts, in particular, in the modeling of the electric power grid. A functional grid corresponds to a stable steady state such that any bifurcation can have catastrophic consequences up to a blackout. Also, the existence of multiple steady states is undesirable as it can lead to transitions or circulatory flows. Despite the high practical importance there is still no general theory of the existence and uniqueness of steady states in such systems. Analytic results are mostly limited to grids without Ohmic losses. In this article, we introduce a method to systematically construct the solutions of the real power load-flow equations in the presence of Ohmic losses and explicitly compute them for tree and ring networks. We investigate different mechanisms leading to multistability and discuss the impact of Ohmic losses on the existence of solutions.

Mechanical Model of Nuclei Ordering in Drosophila Embryos Reveals Dilution of Stochastic Forces.

During the initial development of syncytial embryos, nuclei go through cycles of nuclear division and spatial rearrangement. The arising spatial pattern of nuclei is important for subsequent cellularization and morphing of the embryo. Although nuclei are contained within a common cytoplasm, cytoskeletal proteins are nonuniformly packaged into regions around every nucleus. In fact, cytoskeletal elements like microtubules and their associated motor proteins exert stochastic forces between nuclei, actively driving their rearrangement. Yet, it is unknown how the stochastic forces are balanced to maintain nuclear order in light of increased nuclear density upon every round of divisions. Here, we investigate the nuclear arrangements in Drosophila melanogaster over the course of several nuclear divisions starting from interphase 11. We develop a theoretical model in which we distinguish long-ranged passive forces due to the nuclei as inclusions in the elastic matrix, namely the cytoplasm, and active, stochastic forces arising from the cytoskeletal dynamics mediated by motor proteins. We perform computer simulations and quantify the observed degree of orientational and spatial order of nuclei. Solely doubling the nuclear density upon nuclear division, the model predicts a decrease in nuclear order. Comparing results to experimental recordings of tracked nuclei, we make contradictory observations, finding an increase in nuclear order upon nuclear divisions. Our analysis of model parameters resulting from this comparison suggests that overall motor protein density as well as relative active-force amplitude has to decrease by a factor of about two upon nuclear division to match experimental observations. We therefore expect a dilution of cytoskeletal motors during the rapid nuclear division to account for the increase in nuclear order during syncytial embryo development. Experimental measurements of kinesin-5 cluster lifetimes support this theoretical finding.

Dosimetry intercomparison of four proton therapy institutions in Germany employing spot scanning.

AIM: To verify the consistency of dose and range measurement in an interinstitution comparison among proton therapy institutions in Germany which use the pencil-beam scanning technique. METHODS: Following a peer-to-peer approach absorbed dose and range have been intercompared in several missions at two hosting centers with two or three visiting physics teams of participating institutions using their own dosimetry equipment. A meta-analysis has been performed integrating the results of the individual missions. Dose has been determined with ionization chambers according to the dosimetry protocol IAEA TRS-398. For determination of the depth of the distal 80% dose the teams used either a scanning water phantom, a variable water column or a multi-layer ionization chamber. RESULTS: The systematic deviation between measured doses of the participating institutions is less than 1%. Ranges differ systematically less than 0.4mm. CONCLUSIONS: The match of measured dose and range is better than expected from the respective uncertainties. As all physics teams agree on the assessment of absorbed dose and range, an important prerequisite for a start of joint clinical studies is fulfilled.

Fusion of freehand SPECT and ultrasound: First experience in preoperative localization of sentinel lymph nodes.

PURPOSE: Radioguided sentinel lymph node biopsy (SLNB) is the standard of care in breast cancer and melanoma. Additional preoperative Single-photon emission computed tomography (SPECT/CT) for improved anatomical co-registration of the SLNs causes additional radiation exposure and is time-consuming and expensive. The aim of this prospective study was to evaluate a novel approach involving real-time fusion of freehand SPECT (fhSPECT) and ultrasound (US) for anatomical co-registration of SLNs. METHODS: From February 2015 to February 2016, 153 patients were included in this prospective study. All patients underwent lymphoscintigraphy according to practical guidelines and 151 (118 cases of breast cancer, 30 cutaneous malignancies, and three cases of vulvar cancer) of the 153 patients were additionally investigated with fhSPECT-US. FhSPECT connected to a hand-held gamma detector generates three-dimensional images of the radioactivity distribution in the scanned area. For co-registration and real-time fusion of fhSPECT and subsequently performed US, an infrared stereo tracking system was used. RESULTS: In all patients an SLN was found on lymphoscintigraphy, and the fhSPECT detected corresponding hotspots in all but one patient. In 72 % of patients and 73 % of lymph node basins, real-time anatomical co-registration with US was feasible. The rate of success in achieving good co-registration increased from 60 to 75 % after training by a radiologist specialized in breast imaging. A higher co-registration rate (78 %) was observed in patients with only one SLN than in those with two SLNs (68 %) or three or more SLNs (0 %). CONCLUSIONS: Real-time fusion of fhSPECT and US for preoperative anatomical co-registration of SLNs is feasible. However, before this approach can completely replace preoperative lymphatic imaging, further technical developments are needed.

Assessment of the 18F-labeled PET tracer LMI1195 for imaging norepinephrine handling in rat hearts.

UNLABELLED: A novel (18)F-labeled tracer, LMI1195 (N-[3-bromo-4-(3-(18)F-fluoro-propoxy)-benzyl]-guanidine), is being developed for sympathetic nerve imaging; its high specificity for neural uptake-1 mechanism has previously been demonstrated in cell associative studies and in rabbit and nonhuman primate studies assessing heart uptake. The aim of this study was to investigate the mechanisms of (18)F-LMI1195 cardiac uptake in the rat, which is known to contain norepinephrine uptake mechanisms beyond uptake-1. METHODS: Tracer accumulation in the heart was studied over time after intravenous administration of (18)F-LMI1195 in healthy male Wistar rats by quantitative in vivo PET imaging. The uptake mechanism was assessed by pretreatment with the nonselective norepinephrine uptake-1 and norepinephrine uptake-2 inhibitor phenoxybenzamine (50 mg/kg intravenously; n = 4), the selective norepinephrine uptake-1 inhibitor desipramine (2 mg/kg intravenously; n = 4), or saline control (intravenously; n = 4). RESULTS: (18)F-LMI1195 produced high and sustained heart uptake allowing clear delineation of the left ventricular wall over 60 min after tracer administration. Pretreatment with phenoxybenzamine markedly reduced the (18)F-LMI1195 cardiac uptake when compared with controls. In contrast, there was preserved (18)F-LMI1195 uptake after desipramine pretreatment. CONCLUSION: In rats, cardiac uptake of (18)F-LMI1195 was significantly inhibited by phenoxybenzamine but not desipramine, suggesting (18)F-LMI1195 is a substrate for the uptake-2 mechanism and is consistent with the rat heart having a dominant level of the mechanism.

Insulin-mediated oxidative stress and DNA damage in LLC-PK1 pig kidney cell line, female rat primary kidney cells, and male ZDF rat kidneys in vivo.

Hyperinsulinemia, a condition with excessively high insulin blood levels, is related to an increased cancer incidence. Diabetes mellitus is the most common of several diseases accompanied by hyperinsulinemia. Because an elevated kidney cancer risk was reported for diabetic patients, we investigated the induction of genomic damage by insulin in LLC-PK1 pig kidney cells, rat primary kidney cells, and ZDF rat kidneys. Insulin at a concentration of 5nM caused a significant increase in DNA damage in vitro. This was associated with the formation of reactive oxygen species (ROS). In the presence of antioxidants, blockers of the insulin, and IGF-I receptors, and a phosphatidylinositol 3-kinase inhibitor, the insulin-mediated DNA damage was reduced. Phosphorylation of protein kinase B (PKB or AKT) was increased and p53 accumulated. Inhibition of the mitochondrial and nicotinamide adenine dinucleotide phosphatase oxidase-related ROS production reduced the insulin-mediated damage. In primary rat cells, insulin also induced genomic damage. In kidneys from healthy, lean ZDF rats, which were infused with insulin to yield normal or high blood insulin levels, while keeping blood glucose levels constant, the amounts of ROS and the tumor protein (p53) were elevated in the high-insulin group compared with the control level group. ROS and p53 were also elevated in diabetic obese ZDF rats. Overall, insulin-induced oxidative stress resulted in genomic damage. If the same mechanisms are active in patients, hyperinsulinemia might cause genomic damage through the induction of ROS contributing to the increased cancer risk, against which the use of antioxidants and/or ROS production inhibitors might exert protective effects.

Frequency of transferrin receptor positive reticulocytes (TF-Ret) in blood as an indicator of total-body radiation exposure: a pilot study in nuclear medicine patients.

Approximately 3-20% of all reticulocytes in blood of healthy persons are immature and transferrin receptor positive (Tf-Ret). Tf-Ret were measured by flow cytometry in 27 patients treated with three different radiopharmaceuticals labeled with (131)I and in 25 healthy controls. Patients were chronically exposed within 6 days to blood doses from 0.18-1.89 Gy (D6). Typically, two-thirds of D6 was administered within the first day (D1). The study had to be confined to intra-subject investigations due to high biological variability of Tf-Ret counts. A significant radiation-induced decline was found in patients D1 doses that were ≥0.5 Gy. Tf-Ret frequency declined during the first 4 to 5 days of nuclear therapy to about 30-60% of its initial value, and increased in the following 3 days without reaching the initial value. At the time of nadir, the relative frequency of Tf-Ret was more depressed than that of reticulocytes and lymphocytes. The relative Tf-Ret frequency at nadir could be fitted to the equation: %-Tf-Ret=exp-(D1/D(o)). D(o) was found to be 1.0 ± 0.4 Gy (Mean ± SEM). The study shows that Tf-Ret frequency in blood might be a good parameter for estimation of the radiation dose to red marrow.

A comparison of different experimental methods for general recombination correction for liquid ionization chambers.

Radiation dosimetry of highly modulated dose distributions requires a detector with a high spatial resolution. Liquid filled ionization chambers (LICs) have the potential to become a valuable tool for the characterization of such radiation fields. However, the effect of an increased recombination of the charge carriers, as compared to using air as the sensitive medium has to be corrected for. Due to the presence of initial recombination in LICs, the correction for general recombination losses is more complicated than for air-filled ionization chambers. In the present work, recently published experimental methods for general recombination correction for LICs are compared and investigated for both pulsed and continuous beams. The experimental methods are all based on one of two approaches: either measurements at two different dose rates (two-dose-rate methods), or measurements at three different LIC polarizing voltages (three-voltage methods). In a comparison with the two-dose-rate methods, the three-voltage methods fail to achieve accurate corrections in several instances, predominantly at low polarizing voltages and dose rates. However, for continuous beams in the range of polarizing voltages recommended by the manufacturer of the LICs used, the agreement between the different methods is generally within the experimental uncertainties. For pulsed beams, the agreement between the methods is poor. The inaccuracies found in the results from the three-voltage methods are associated with numerical difficulties in solving the resulting equation systems, which also make these methods sensitive to small variations in the experimental data. These issues are more pronounced for the case of pulsed beams. Furthermore, the results suggest that the theoretical modelling of initial recombination used in the three-voltage methods may be a contributing factor to the deviating results observed.

The estrogen receptor-α is required and sufficient to maintain physiological glucose uptake in the mouse heart.

Estrogens attenuate cardiac hypertrophy and increase cardiac contractility via their cognate estrogen receptors (ERs) ERα and ERß. Because female sex hormones enhance global glucose use and because myocardial function and mass are tightly linked to cardiac glucose metabolism, we tested the hypothesis that expression and activation of the ERα might be required and sufficient to maintain physiological cardiac glucose uptake in the murine heart. Cardiac glucose uptake quantified in vivo by 18F-fluorodeoxyglucose positron emission tomography was strongly impaired in ovariectomized compared with gonadal intact female C57BL/6JO mice. The selective ERα agonist 16α-LE2 and the nonselective ERα and ERß agonist 17ß-estradiol completely restored cardiac glucose uptake in ovariectomized mice. Cardiac 18F-fluorodeoxyglucose uptake was strongly decreased in female ERα knockout mice compared with wild-type littermates. Analysis of cardiac mRNA accumulation by quantitative RT-PCR revealed an upregulation of genes involved in glycolisis and tricarboxylic acid cycle by ERα treatment. In conclusion, systemic activation of ERα is sufficient, and its expression is required to maintain physiological glucose uptake in the murine heart, which is likely to contribute to known cardioprotective estrogen effects.

Initial recombination in the track of heavy charged particles: numerical solution for air filled ionization chambers.

INTRODUCTION: Modern particle therapy facilities enable sub-millimeter precision in dose deposition. Here, also ionization chambers (ICs) are used, which requires knowledge of the recombination effects. Up to now, recombination is corrected using phenomenological approaches for practical reasons. In this study the effect of the underlying dose distribution on columnar recombination, a quantitative model for initial recombination, is investigated. MATERIAL AND METHODS: Jaffé's theory, formulated in 1913 quantifies initial recombination by elemental processes, providing an analytical (closed) solution. Here, we investigate the effect of the underlying charged carrier distribution around a carbon ion track. The fundamental partial differential equation, formulated by Jaffé, is solved numerically taking into account more realistic charge carrier distributions by the use of a computer program (Gascoigne 3D). The investigated charge carrier distributions are based on track structure models, which follow a 1/r(2) behavior at larger radii and show a constant value at small radii. The results of the calculations are compared to the initial formulation and to data obtained in experiments using carbon ion beams. RESULTS: The comparison between the experimental data and the calculations shows that the initial approach made by Jaffé is able to reproduce the effects of initial recombination. The amorphous track structure based charge carrier distribution does not reproduce the experimental data well. A small additional correction in the assessment of the saturation current or charge is suggested by the data. CONCLUSION: The established model of columnar recombination reproduces the experimental data well, whereas the extensions using track structure models do not show such an agreement. Additionally, the effect of initial recombination on the saturation curve (i.e. Jaffé plot) does not follow a linear behavior as suggested by current dosimetry protocols, therefore higher order corrections (such as the investigated ones) might be necessary.

Characterization of the optical properties and stability of Presage™ following irradiation with photons and carbon ions.

BACKGROUND: The on-going development of both intensity-modulated radiotherapy (IMRT), including the more recent intensity-modulated arc therapy, as well as particle beam therapy, has created a clear need for accurate verification of dose distributions in three dimensions (3D). Presage™ is a new 3D dosimetry material that exhibits a radiochromic response when exposed to ionizing radiation. In this study we have 1) developed an improved optical set-up for measurements of changes in OD of Presage™ point dosimeters, 2) investigated the dose response of Presage™ for photons and carbon ions in the therapy range, 3) investigated the dose response of Presage™ for photons in the kGy range and 4) investigated the fading (i.e. bleaching) of Presage™ postirradiation. MATERIALS AND METHODS: Presage™ was examined in 1 × 1 × 4.5 cm(3) optical cuvettes; a cuvette holder assured accurate repositioning, and the optical setup included a reference detector to take into account laser intensity fluctuations. The cuvettes were measured pre- and postirradiation for a two week period. RESULTS: A linear response was observed between dose and optical response between 0 Gy and 100 Gy for γ-radiation from Co-60 and for carbon ions (both plateau and SOBP) from 0 to 20 Gy. The dosimeter was found to have a saturation dose of approximately 100 Gy for photons. A linear energy transfer (LET) effect was not observed in the dose response of different LET radiation. The postirradiation change in optical fading was found to be 0.5% ΔOD/day. CONCLUSIONS: Our study shows that Presage™ remains a dosimeter of interest for radiation therapy with other particles as well as photons in the therapy dose range.

Spatially resolved ballistic optoelectronic transport measured by quantized photocurrent spectroscopy.

GaAs-based quantum point contacts (QPCs) are exploited to spatially resolve and analyze the ballistic, nonequilibrium flow of photogenerated electrons in a nanoscale circuit. Electron-hole pairs are photogenerated in a two-dimensional electron gas (2DEG), and the resulting current through an adjacent QPC is measured as a function of the laser spot position. The transmission of photogenerated electrons through the QPC is governed by the energy dispersion and the quantized momentum values of the electron modes in the QPC.

COTS Silicon diodes as radiation detectors in proton and heavy charged particle radiotherapy 1.

Modern radiotherapy facilities for cancer treatment such as the Heavy Ion Therapy Center (HIT) in Heidelberg, Germany, allow for sub-millimeter precision in dose deposition. For measurement of such dose distributions and characterization of the particle beams, detectors with high spatial resolution are necessary. Here, a detector based on the commercially available COTS photodiode (BPW-34) is presented. When applied in hadronic beams of protons and carbon ions, the detector reproduces dose distribution well, but its response decreases rapidly by radiation damage. However, for MeV photon beams, the detector exhibits a similar behavior as found in diode detectors usually applied in radiotherapy.