Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPLEXITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The solution structures and relative stability constants of lanthanide–EDTA complexes predicted from computation

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Evaluating the photophysical and photochemical characteristics of green-emitting cerium(III) mono-cyclooctatetraenide complexes

The photophysical and photochemical reactivity of the organometallic mono(cyclooctatetraenide)-Ce(III) complexes: [(C 8 H 8 )Ce(μ-X)(THF) 2 ] 2 X = O 3 SCF 3 - (1), Cl - (2), were studied. Here, in the course of these studies, a new polymorph of [(C 8 H 8 )Ce(μ-O 3 SCF 3 )(THF) 2 ] 2 (1) was described. Photoluminescence (PL) studies of 1 and 2 showed characteristic excitation and emission bands in the visible region, with bright green light emission under light irradiation from states corresponding to the 5d → 4f transition. Complexes 1 and 2 exhibit relatively long lifetime of 205.4 ± 0.2 and 145.8 ± 0.4 ns respectively. And the quantum yields (Φ) for 1 and 2 are 0.180 and 0.068 respectively. Electrochemical studies were performed on complexes 1 and 2 with a reversible Ce(III)/Ce(IV) redox couple recorded at $E_{1/2}$ = -1.53 V versus the Fc/Fc + for 2. Complex 1 shows an irreversible Ce(III/IV) oxidation wave. Complexes 1 and 2 revealed strongly-reducing, estimated, excited state reduction potentials ($E^*_{1/2}$) of -3.31 and -4.02 V versus the Fc/Fc + respectively, with small Stokes shifts of 0.12 eV. With their associated relatively long lifetimes, small Stokes shifts and large, negative $E^*_{1/2}$ values, both complexes were evaluated as potential photosensitisers for halogen atom transfer (XAT) using a test reaction of the dehalogenation of benzyl chloride.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me 2 -1.1.1, and 3 6 Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3) 2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6–2.0 eV) and hence these complexes belong to the superalkali category. Further, the ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3) 2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Hydrogenation of ethylene over molybdenum–sulfur complexes supported on UiO-66

Development of supported single-site catalysts using small metal sulfide complexes could significantly help in the development of cost-effective catalytic materials to drive selective hydrogenation and hydrogenolysis. The goal of this study is to contribute to the development of metal sulfide catalysts by calculating the thermodynamics of a catalytic cyclic involving a metal organic framework functionalized by insertion of metal sulfide. Anchored metal sulfide complexes can potentially be designed with ligands with distinctly different electronic and catalytic properties for specific catalytic applications. Here we examine the hydrogenation of ethylene as a model. We use density functional theory to investigate molybdenum–sulfur complexes as active catalysts anchored on the metal–organic framework UiO-66 as a stable support. Our calculations show that the anchored complexes with more than two sulfur ligands are unfavorable for ethylene adsorption, so we study complexes with one or two sulfur ligands. Hydrogenation of the unsaturated carbon double bond requires the transfer of two hydrogen atoms, which can occur via heterolytic activation of hydrogen to form a Mo-hydride and a protonated sulfur – either by hydride transfer followed by proton transfer or via proton transfer followed by hydride transfer, and we find that both mechanisms proceed via two-state reactivity involving two spin states along the reaction path. Of the two catalysts studied in gas the phase, the MoS single-sulfur–ligand complex with lower oxidation states produces thermodynamically more favorable intermediates along the pathway for the first hydrogen transfer for both the hydride-first mechanism and the proton-first mechanism. As a result, the quantum mechanical calculations provide experimentally inaccessible partial atomic charges and geometries of the various intermediates encountered along the steps of the reaction mechanisms.

Kermani, Maryam Mansoori [University of Minnesota,↗
DOE OSTI2025

Spectroscopic investigation of cation effects in U(VI)-NO3- complexation in aqueous solutions

Understanding and manipulating uranyl speciation in aqueous solutions is critical for advancing chemical separation, sensing, and understanding environmental transport of uranyl. We report on the significant enhancement in the complexation of uranyl with nitrate in aqueous solutions containing quaternary ammonium cations leading to the formation of anionic complexes. We base this on the comparative study of the effect of monovalent cations (Na+, Li+, NH4+, and N(CH3)4+) on the complexation equilibria of uranyl-nitrate in solutions, probed by time-resolved laser induced fluorescence spectroscopy (TRLFS). Lifetime-corrected spectra, obtained by extrapolating the time-resolved spectra to t = 0, were used to study speciation to mitigate the effects of variations in the fluorescence lifetimes that depend on extraneous factors such as dynamic quenching. We demonstrate that the lifetime-corrected spectra can be used to determine uranyl speciation in aqueous solutions where the mono-nitrate complex forms, and in acetonitrile with N(CH3)4NO3 where di-nitrate and tris-nitrato species are formed. Aqueous solutions containing N(CH3)4+ are shown to promote the formation of higher complexes of uranyl compared to other inorganic nitrate salts based on the higher redshift in the spectra, poor fits to the two-component model, and the higher apparent formation constants. Comparing the trends in uranyl speciation in the presence of N(CH3)4+ in aqueous and acetonitrile solutions, it is proposed that the quaternary ammonium cations (quats) promote the formation of anionic complexes of uranyl by the ion-association mechanism. These results provide a basis for designing quat-assisted separation systems that target anionic actinide species in aqueous solutions.

Nayak, Srikanth (ORCID:0000000302135796)↗
DOE OSTI2012

Robust inverse scattering full waveform seismic tomography for imaging complex structure

Seismic tomography becomes important tool recently for imaging complex subsurface. It is well known that imaging complex rich fault zone is difficult. In this paper, The application of time domain inverse scattering wave tomography to image the complex fault zone would be shown on this paper, especially an efficient time domain inverse scattering tomography and their run in cluster parallel computer which has been developed. This algorithm is purely based on scattering theory through solving Lippmann Schwienger integral by using Born's approximation. In this paper, it is shown the robustness of this algorithm especially in avoiding the inversion trapped in local minimum to reach global minimum. A large data are solved by windowing and blocking technique of memory as well as computation. Parameter of windowing computation is based on shot gather's aperture. This windowing technique reduces memory as well as computation significantly. This parallel algorithm is done by means cluster system of 120 processors from 20 nodes of AMD Phenom II. Benchmarking of this algorithm is done by means Marmoussi model which can be representative of complex rich fault area. It is shown that the proposed method can image clearly the rich fault and complex zone in Marmoussi model even though the initial model is quite far from the true model. Therefore, this method can be as one of solution to image the very complex mode.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Investigation of electron-atom/molecule scattering resonances: Two complex multiconfigurational self-consistent field approaches

Resonances are temporarily bound states which lie in the continuum part of the Hamiltonian. If the electronic coordinates of the Hamiltonian are scaled (“dilated”) by a complex parameter, η = αe{sup iθ} (α, θ real), then its complex eigenvalues represent the scattering states (resonant and non-resonant) while the eigenvalues corresponding to the bound states and the ionization and the excitation thresholds remain real and unmodified. These make the study of these transient species amenable to the bound state methods. We developed a quadratically convergent multiconfigurational self-consistent field method (MCSCF), a well-established bound-state technique, combined with a dilated Hamiltonian to investigate resonances. This is made possible by the adoption of a second quantization algebra suitable for a set of “complex conjugate biorthonormal” spin orbitals and a modified step-length constraining algorithm to control the walk on the complex energy hypersurface while searching for the stationary point using a multidimensional Newton-Raphson scheme. We present our computational results for the {sup 2}PBe{sup −} shape resonances using two different computationally efficient methods that utilize complex scaled MCSCF (i.e., CMCSCF). These two methods are to straightforwardly use CMCSCF energy differences and to obtain energy differences using an approximation to the complex multiconfigurational electron propagator. It is found that, differing from previous computational studies by others, there are actually two {sup 2}PBe{sup −} shape resonances very close in energy. In addition, N{sub 2} resonances are examined using one of these methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The roles of SDHAF2 and dicarboxylate in covalent flavinylation of SDHA, the human complex II flavoprotein

Mitochondrial complex II, also known as succinate dehydrogenase (SDH), is an integral-membrane heterotetramer (SDHABCD) that links two essential energy-producing processes, the tricarboxylic acid (TCA) cycle and oxidative phosphorylation. A significant amount of information is available on the structure and function of mature complex II from a range of organisms. However, there is a gap in our understanding of how the enzyme assembles into a functional complex, and disease-associated complex II insufficiency may result from incorrect function of the mature enzyme or from assembly defects. Here, we investigate the assembly of human complex II by combining a biochemical reconstructionist approach with structural studies. We report an X-ray structure of human SDHA and its dedicated assembly factor SDHAF2. Importantly, we also identify a small molecule dicarboxylate that acts as an essential cofactor in this process and works in synergy with SDHAF2 to properly orient the flavin and capping domains of SDHA. This reorganizes the active site, which is located at the interface of these domains, and adjusts the pK a of SDHA R451 so that covalent attachment of the flavin adenine dinucleotide (FAD) cofactor is supported. We analyze the impact of disease-associated SDHA mutations on assembly and identify four distinct conformational forms of the complex II flavoprotein that we assign to roles in assembly and catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Two interaction surfaces between XPA and RPA organize the preincision complex in nucleotide excision repair

The xeroderma pigmentosum protein A (XPA) and replication protein A (RPA) proteins fulfill essential roles in the assembly of the preincision complex in the nucleotide excision repair (NER) pathway. We have previously characterized the two interaction sites, one between the XPA N-terminal (XPA-N) disordered domain and the RPA32 C-terminal domain (RPA32C), and the other with the XPA DNA binding domain (DBD) and the RPA70AB DBDs. Here, we show that XPA mutations that inhibit the physical interaction in either site reduce NER activity in biochemical and cellular systems. Combining mutations in the two sites leads to an additive inhibition of NER, implying that they fulfill distinct roles. Our data suggest a model in which the interaction between XPA-N and RPA32C is important for the initial association of XPA with NER complexes, while the interaction between XPA DBD and RPA70AB is needed for structural organization of the complex to license the dual incision reaction. Integrative structural models of complexes of XPA and RPA bound to single-stranded/double-stranded DNA (ss/dsDNA) junction substrates that mimic the NER bubble reveal key features of the architecture of XPA and RPA in the preincision complex. Most critical among these is that the shape of the NER bubble is far from colinear as depicted in current models, but rather the two strands of unwound DNA must assume a U-shape with the two ss/dsDNA junctions localized in close proximity. Our data suggest that the interaction between XPA and RPA70 is key for the organization of the NER preincision complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mechanism of RanGTP priming H2A-H2B release from Kap114 in an atypical RanGTP•Kap114•H2A-H2B complex

Previously, we showed that the nuclear import receptor Importin-9 wraps around the H2A-H2B core to chaperone and transport it from the cytoplasm to the nucleus. However, unlike most nuclear import systems where RanGTP dissociates cargoes from their importins, RanGTP binds stably to the Importin-9•H2A-H2B complex, and formation of the ternary RanGTP•Importin-9•H2A-H2B complex facilitates H2A-H2B release to the assembling nucleosome. It was unclear how RanGTP and the cargo H2A-H2B can bind simultaneously to an importin, and how interactions of the three components position H2A-H2B for release. Here, we show cryo-EM structures of Importin-9•RanGTP and of its yeast homolog Kap114, including Kap114•RanGTP, Kap114•H2A-H2B, and RanGTP•Kap114•H2A-H2B, to explain how the conserved Kap114 binds H2A-H2B and RanGTP simultaneously and how the GTPase primes histone transfer to the nucleosome. In the ternary complex, RanGTP binds to the N-terminal repeats of Kap114 in the same manner as in the Kap114/Importin-9•RanGTP complex, and H2A-H2B binds via its acidic patch to the Kap114 C-terminal repeats much like in the Kap114/Importin-9•H2A-H2B complex. Ran binds to a different conformation of Kap114 in the ternary RanGTP•Kap114•H2A-H2B complex. Here, Kap114 no longer contacts the H2A-H2B surface proximal to the H2A docking domain that drives nucleosome assembly, positioning it for transfer to the assembling nucleosome or to dedicated H2A-H2B chaperones in the nucleus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Single-domain antibodies pinpoint potential targets within Shigella invasion plasmid antigen D of the needle tip complex for inhibition of type III secretion

Numerous Gram-negative pathogens infect eukaryotes and use the type III secretion system (T3SS) to deliver effector proteins into host cells. One important T3SS feature is an extracellular needle with an associated tip complex responsible for assembly of a pore-forming translocon in the host cell membrane. Shigella spp. cause shigellosis, also called bacillary dysentery, and invade colonic epithelial cells via the T3SS. The tip complex of Shigella flexneri contains invasion plasmid antigen D (IpaD), which initially regulates secretion and provides a physical platform for the translocon pore. The tip complex represents a promising therapeutic target for many important T3SS-containing pathogens. Here, in an effort to further elucidate its function, we created a panel of single-VH domain antibodies (VHHs) that recognize distinct epitopes within IpaD. These VHHs recognized the in situ tip complex and modulated the infectious properties of Shigella. Moreover, structural elucidation of several IpaD–VHH complexes provided critical insights into tip complex formation and function. Of note, one VHH heterodimer could reduce Shigella hemolytic activity by >80%. Our observations along with previous findings support the hypothesis that the hydrophobic translocator (IpaB in Shigella) likely binds to a region within the tip protein that is structurally conserved across all T3SS-possessing pathogens, suggesting potential therapeutic avenues for managing infections by these pathogens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

STAR FORMATION ACROSS THE W3 COMPLEX

We present a multi-wavelength analysis of the history of star formation in the W3 complex. Using deep, near-infrared ground-based images combined with images obtained with Spitzer and Chandra observatories, we identified and classified young embedded sources. We identified the principal clusters in the complex and determined their structure and extension. We constructed extinction-limited samples for five principal clusters and constructed K-band luminosity functions that we compare with those of artificial clusters with varying ages. This analysis provided mean ages and possible age spreads for the clusters. We found that IC 1795, the centermost cluster of the complex, still hosts a large fraction of young sources with circumstellar disks. This indicates that star formation was active in IC 1795 as recently as 2 Myr ago, simultaneous to the star-forming activity in the flanking embedded clusters, W3-Main and W3(OH). A comparison with carbon monoxide emission maps indicates strong velocity gradients in the gas clumps hosting W3-Main and W3(OH) and shows small receding clumps of gas at IC 1795, suggestive of rapid gas removal (faster than the T Tauri timescale) in the cluster-forming regions. We discuss one possible scenario for the progression of cluster formation in the W3 complex. We propose that early processes of gas collapse in the main structure of the complex could have defined the progression of cluster formation across the complex with relatively small age differences from one group to another. However, triggering effects could act as catalysts for enhanced efficiency of formation at a local level, in agreement with previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Energetic lanthanide complexes: coordination chemistry and explosives applications

Here, we report that metals are generally added to organic molecular explosives in a heterogeneous composite to improve overall heat and energy release. In order to avoid creating a mixture that can vary in homogeneity, energetic organic molecules can be directly bonded to high molecular weight metals, forming a single metal complex with Angstrom-scale separation between the metal and the explosive. To probe the relationship between the structural properties of metal complexes and explosive performance, a new series of energetic lanthanide complexes has been prepared using energetic ligands such as NTO (5-nitro-2,4-dihydro-1,2,4-triazole-3-one). These are the first examples of lanthanide NTO complexes where no water is coordinated to the metal, demonstrating novel control of the coordination environment. The complexes have been characterized by X-ray crystallography, NMR and IR spectroscopies, photoluminescence, and sensitivity testing. The structural and energetic properties are discussed in the context of enhanced blast effects and detection. Cheetah calculations have been performed to fine-tune physical properties, creating a systematic method for producing explosives with 'tailor made' characteristics. In conclusion, these new complexes will be benchmarks for further study in the field of metalized high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Modelling the evolution of complex conductivity during calcite precipitation on glass beads

When pH and alkalinity increase, calcite frequently precipitates and hence modifies the petrophysical properties of porous media. The complex conductivity method can be used to directly monitor calcite precipitation in porous media because it is sensitive to the evolution of the mineralogy, pore structure and its connectivity. Here, we have developed a mechanistic grain polarization model considering the electrochemical polarization of the Stern and diffuse layers surrounding calcite particles. Our complex conductivity model depends on the surface charge density of the Stern layer and on the electrical potential at the onset of the diffuse layer, which are computed using a basic Stern model of the calcite/water interface. The complex conductivity measurements of Wu et al. on a column packed with glass beads where calcite precipitation occurs are reproduced by our surface complexation and complex conductivity models. The evolution of the size and shape of calcite particles during the calcite precipitation experiment is estimated by our complex conductivity model. At the early stage of the calcite precipitation experiment, modelled particles sizes increase and calcite particles flatten with time because calcite crystals nucleate at the surface of glass beads and grow into larger calcite grains. Finally, at the later stage of the calcite precipitation experiment, modelled sizes and cementation exponents of calcite particles decrease with time because large calcite grains aggregate over multiple glass beads and only small calcite crystals polarize.

58 GEOSCIENCES↗
DOE OSTI2021

Understanding protein-complex assembly through grand canonical maximum entropy modeling

Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded systems where the abundance of these proteins vary with cell types. While some protein complexes form putative structures that can be visualized with imaging, there are far more protein complexes that are yet to be solved because of their dynamic associations with one another. Nevertheless, it is possible to infer these protein complexes through a physical model. However, it is often not clear to physicists what kind of data from biology is necessary for such a modeling endeavor. Here, we aim to model these clusters of coarse-grained protein assemblies from multiple subunits through the constraints of interactions among the subunits and the chemical potential of each subunit. We obtained the constraints on the interactions among subunits from the known protein structures. We inferred the chemical potential that dictates the particle number distribution of each protein subunit from the knowledge of protein abundance from experimental data. Guided by the maximum entropy principle, we formulated an inverse statistical mechanical method to infer the distribution of particle numbers from the data of protein abundance as chemical potentials for a grand canonical multicomponent mixture. Using grand canonical Monte Carlo simulations, we captured a distribution of high-order clusters in a protein complex of succinate dehydrogenase with four known subunits. The complexity of hierarchical clusters varies with the relative protein abundance of each subunit in distinctive cell types such as lung, heart, and brain. When the crowding content increases, we observed that crowding stabilizes emergent clusters that do not exist in dilute conditions. We, therefore, proposed a testable hypothesis that the hierarchical complexity of protein clusters on a molecular scale is a plausible biomarker of predicting the phenotypes of a cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Charge density studies of 3 d metal (Ni/Cu) complexes with a non-innocent ligand

High-resolution X-ray diffraction experiments and atom-specific X-ray absorption experiments are applied to investigate a series of square planar complexes with the non-innocent ligand of maleonitriledithiolate (mnt), [S 2 C 2 (CN) 2 ] z - , containing M —S bonds. Four complexes of (PyH) z [ M (mnt) 2 ] z - , where M = Ni or Cu, z = 2 or 1 and PyH + = C 5 NH 6 + , were studied in order to clarify whether such one-electron oxidation–reduction, [ M (mnt) 2 ] 2- /[ M (mnt) 2 ] 1- , is taking place at the metal or the ligand site. Combining the techniques of metal K -, L -edge and S K -edge X-ray absorption spectroscopy with high-resolution X-ray charge density studies, it is unambiguously demonstrated that the electron redox reaction is ligand based and metal based for Ni and Cu pairs, respectively. The bonding characters in terms of topological properties associated with the bond critical points are compared between the oxidized form [ ML ] - and the reduced form [ ML ] 2- . In the case of Ni complexes, the formal oxidation state of Ni remains as Ni 2+ and each mnt ligand carries a 2- charge in [Ni(mnt) 2 ] 2- , but only one of the ligands is formally oxidized in [Ni(mnt) 2 ] 1- . In contrast, in the case of Cu complexes, the mnt remains as 2- in both complexes, but the formal oxidation states of the metal are Cu 2+ and Cu 3+ . Bond characterizations and d -orbital populations will be presented. The complementary results of XAS, XRD and DFT calculations will be discussed. The conclusion on the redox reactions in these complexes can be firmly established.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Conforming Morse-Smale Complexes

Morse-Smale (MS) complexes have been gaining popularity as a tool for feature-driven data analysis and visualization. However, the quality of their geometric embedding and the sole dependence on the input scalar field data can limit their applicability when expressing application-dependent features. In this paper we introduce a new combinatorial technique to compute an MS complex that conforms to both an input scalar field and an additional, prior segmentation of the domain. The segmentation constrains the MS complex computation guaranteeing that boundaries in the segmentation are captured as separatrices of the MS complex. We demonstrate the utility and versatility of our approach with two applications. First, we use streamline integration to determine numerically computed basins/mountains and use the resulting segmentation as an input to our algorithm. This strategy enables the incorporation of prior flow path knowledge, effectively resulting in an MS complex that is as geometrically accurate as the employed numerical integration. Our second use case is motivated by the observation that often the data itself does not explicitly contain features known to be present by a domain expert. We introduce edit operations for MS complexes so that a user can directly modify their features while maintaining all the advantages of a robust topology-based representation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Second cancer incidence risk estimates using BEIR VII models for standard and complex external beam radiotherapy for early breast cancer

Purpose: To compare organ specific cancer incidence risks for standard and complex external beam radiotherapy (including cone beam CT verification) following breast conservation surgery for early breast cancer.Method: Doses from breast radiotherapy and kilovoltage cone beam CT (CBCT) exposures were obtained from thermoluminescent dosimeter measurements in an anthropomorphic phantom in which the positions of radiosensitive organs were delineated. Five treatment deliveries were investigated: (i) conventional tangential field whole breast radiotherapy (WBRT), (ii) noncoplanar conformal delivery applicable to accelerated partial beast irradiation (APBI), (iii) two-volume simultaneous integrated boost (SIB) treatment, (iv) forward planned three-volume SIB, and (v) inverse-planned three volume SIB. Conformal and intensity modulated radiotherapy methods were used to plan the complex treatments. Techniques spanned the range from simple methods appropriate for patient cohorts with a low cancer recurrence risk to complex plans relevant to cohorts with high recurrence risk. Delineated organs at risk included brain, salivary glands, thyroid, contralateral breast, left and right lung, esophagus, stomach, liver, colon, and bladder. Biological Effects of Ionizing Radiation (BEIR) VII cancer incidence models were applied to the measured mean organ doses to determine lifetime attributable risk (LAR) for ages at exposure from 35 to 80 yr according to radiotherapy techniques, and included dose from the CBCT imaging. Results: All LAR decreased with age at exposure and were lowest for brain, thyroid, liver, and bladder (<0.1%). There was little dependence of LAR on radiotherapy technique for these organs and for colon and stomach. LAR values for the lungs for the three SIB techniques were two to three times those from WBRT and APBI. Uncertainties in the LAR models outweigh any differences in lung LAR between the SIB methods. Constraints in the planning of the SIB methods ensured that contralateral breast doses and LAR were comparable to WBRT, despite their added complexity. The smaller irradiated volume of the ABPI plan contributed to a halving of LAR for contralateral breast compared with the other plan types. Daily image guided radiotherapy (IGRT) for a left breast protocol using kilovoltage CBCT contributed <10% to LAR for the majority of organs, and did not exceed 22% of total organ dose. Conclusions: Phantom measurements and calculations of LAR from the BEIR VII models predict that complex breast radiotherapy techniques do not increase the theoretical risk of second cancer incidence for organs distant from the treated breast, or the contralateral breast where appropriate plan constraints are applied. Complex SIB treatments are predicted to increase the risk of second cancer incidence in the lungs compared to standard whole breast radiotherapy; this is outweighed by the threefold reduction in 5 yr local recurrence risk for patients of high risk of recurrence, and young age, from the use of radiotherapy. APBI may have a favorable impact on risk of second cancer in the contralateral breast and lung for older patients at low risk of recurrence. Intensive use of IGRTincreased the estimated values of LAR but these are dominated by the effect of the dose from the radiotherapy, and any increase in LAR from IGRT is much lower than the models' uncertainties.

60 APPLIED LIFE SCIENCES↗