Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct Sequence”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Unraveling the mechanism of recognition of the 3’ splice site of the adenovirus major late promoter intron by the alternative splicing factor PUF60

Pre-mRNA splicing is critical for achieving required amounts of a transcript at a given time and for regulating production of encoded protein. A given pre-mRNA may be spliced in many ways, or not at all, giving rise to multiple gene products. Numerous splicing factors are recruited to pre-mRNA splice sites to ensure proper splicing. One such factor, the 60 kDa poly(U)-binding splicing factor (PUF60), is recruited to sites that are not always spliced, but rather function as alternative splice sites. In this study, we characterized the interaction of PUF60 with a splice site from the adenovirus major late promoter (the AdML 3' splice site, AdML3’). We found that the PUF60–AdML3’ dissociation constants are in the micromolar range, with the binding affinity predominantly provided by PUF60’s two central RNA recognition motifs (RRMs). A 1.95 Å crystal structure of the two PUF60 RRMs in complex with AdML3’ revealed a dimeric organization placing two stretches of nucleic acid tracts in opposing directionalities, which can cause looping of nucleic acid and explain how PUF60 affects pre-mRNA geometry to effect splicing. Solution characterization of this complex by light-scattering and UV/Vis spectroscopy suggested a potential 2:1 (PUF60 2 :AdML3’) stoichiometry, consistent with the crystal structure. This work defines the sequence specificity of the alternative splicing factor PUF60 at the pre-mRNA 3’ splice site. Our observations suggest that control of pre-mRNA directionality is important in the early stage of spliceosome assembly, and advance our understanding of the molecular mechanism by which alternative and constitutive splicing factors differentiate among 3’ splice sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Automated Generation of Weather and Climate Analysis Products

Wind roses are an important part to site operations as they depict the wind speed and direction percentage over a time period. In this project, I automated the production of wind roses for SRS meteorological towers. I developed a script to sequence through the dataset by period and create a wind rose for each of those period sets. Graphics were generated for every 4-hour period of the day, for each of the 4 heights of the instruments, for every half of each month. A 10-year climatological period consisted of wind speed and directions measurements at 15-minute intervals. We used the years 2014-2024 as the wind instruments on the tower were upgraded to sonic anemometers in early 2014. I then compared these wind roses to those from a previous study done in 2003 to analyze the differences and similarities in the wind patterns. The wind roses are also used to understand environmental transport conditions at SRS. When comparing the different heights of the anemometers to the 2003 report, there is a similarity in the fact that the direction is relatively the same across the levels, with the wind speeds increasing as you get higher in the air. When comparing the season of our graphs against the 2003 report, we noticed that the winter and summer months appear to have very similar wind directions. However, in the spring we noticed that there were more southerly winds compared to the 2003 report which had more westerly winds. Finally, in the fall months, we noticed that there was a higher percentage of northeasterly winds, while the 2003 report showed more southeasterly winds. We can use these summaries to estimate the directional dependence of dose that the surrounding areas receive throughout the year according to time of day.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Symmetry-Breaking Dendrimer Synthons in Colloidal Crystal Engineering with DNA

Breaking symmetry in colloidal crystals is challenging due to the inherent chemical and structural isotropy of many nanoscale building blocks. If a non-particle component could be used to anisotropically encode such building blocks with orthogonal recognition properties, one could expand the scope of structural and compositional possibilities of colloidal crystals beyond what is thus far possible with purely particle-based systems. In this work, we report the synthesis and characterization of novel DNA dendrimers that function as symmetry-breaking synthons, capable of programming anisotropic and orthogonal interactions within colloidal crystals. When the DNA dendrimers have identical sticky ends, they hybridize with DNA-functionalized nanoparticles to yield three distinct colloidal crystals, dictated by dendrimer size, including a structure not previously reported in the field of colloidal crystal engineering, Si 2 Sr. When used as symmetry-breaking synthons (when the sticky ends deliberately consist of orthogonal sequences), the synthesis of binary and ternary colloidal alloys with structures that can only be realized through directional interactions is possible. Furthermore, by modulating the extent of shape anisotropy within the DNA dendrimers, the local distribution of the nanoparticles within the crystals can be directed.

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

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Borgs are giant genetic elements with potential to expand metabolic capacity

Anaerobic methane oxidation exerts a key control on greenhouse gas emissions1, yet factors that modulate the activity of microorganisms performing this function remain poorly understood. Here we discovered extraordinarily large, diverse DNA sequences that primarily encode hypothetical proteins through studying groundwater, sediments and wetland soil where methane production and oxidation occur. Four curated, complete genomes are linear, up to approximately 1 Mb in length and share genome organization, including replichore structure, long inverted terminal repeats and genome-wide unique perfect tandem direct repeats that are intergenic or generate amino acid repeats. We infer that these are highly divergent archaeal extrachromosomal elements with a distinct evolutionary origin. Gene sequence similarity, phylogeny and local divergence of sequence composition indicate that many of their genes were assimilated from methane-oxidizing Methanoperedens archaea. We refer to these elements as ‘Borgs’. We identified at least 19 different Borg types coexisting with Methanoperedens spp. in four distinct ecosystems. Borgs provide methane-oxidizing Methanoperedens archaea access to genes encoding proteins involved in redox reactions and energy conservation (for example, clusters of multihaem cytochromes and methyl coenzyme M reductase). These data suggest that Borgs might have previously unrecognized roles in the metabolism of this group of archaea, which are known to modulate greenhouse gas emissions, but further studies are now needed to establish their functional relevance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Reconstruction of the Room D, B, G, and Q Closure Histories at the Waste Isolation Pilot Plant

Experimental measurements of room closure in salt repositories are valuable for understanding the evolution of the underground and for validating geomechanical models. Room closure was measured during a number of experiments at the Waste Isolation Pilot Plant (WIPP) during the 1980's and 1990's. Most rooms were excavated using a multi-pass mining sequence, where each pass necessarily destroyed some of the mining sequence closure measurement points. These destroyed points were promptly reinstalled to capture the closure after the mining pass. After the room was complete, the mining sequence closure measurement stations were supplemented with remotely read closure measurement stations. Although many aspects of these experiments were thoroughly documented, the digital copies of the closure data were inadvertently destroyed, the non-trivial process of zeroing and shifting the raw closure measurements after each mining pass was not precisely described, the various closure measurements within a given room were not directly compared on the same plot, and the measurements were collected for several years longer than previously reported. Consequently, the hand-written mining sequence closure measurements for Rooms D, B, G, and Q were located in the WIPP archives, digitized, and reanalyzed for this report. The process of reconstructing the mining sequence closure histories was documented in detail and the raw data can be found in the appendices. Within the mid-section of a given room, the reconstructed closure histories were largely consistent with other mining sequence and remotely read closure histories, which builds confidence in the experiments and suggests that plane strain is an appropriate modeling assumption. The reconstructed closure histories were also reasonably consistent with previously published results, except in one notable case: the reconstructed Room Q closure histories 30 days after excavation were about 45 % less than the corresponding closures reported in Munson's 1997 capstone paper.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Toward Robust and Routine Determination of Mw for Small Earthquakes: Application to the 2020 Mw 5.7 Magna, Utah, Seismic Sequence

To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.

58 GEOSCIENCES↗
DOE OSTI2023

DeCOIL: Optimization of Degenerate Codon Libraries for Machine Learning-Assisted Protein Engineering

With advances in machine learning (ML)-assisted protein engineering, models based on data, biophysics, and natural evolution are being used to propose informed libraries of protein variants to explore. Synthesizing these libraries for experimental screens is a major bottleneck, as the cost of obtaining large numbers of exact gene sequences is often prohibitive. Degenerate codon (DC) libraries are a cost-effective alternative for generating combinatorial mutagenesis libraries where mutations are targeted to a handful of amino acid sites. However, existing computational methods to optimize DC libraries to include desired protein variants are not well suited to design libraries for ML-assisted protein engineering. To address these drawbacks, we present DEgenerate Codon Optimization for Informed Libraries (DeCOIL), a generalized method which directly optimizes DC libraries to be useful for protein engineering: to sample protein variants that are likely to have both high fitness and high diversity in the sequence search space. Using computational simulations and wet-lab experiments, we demonstrate that DeCOIL is effective across two specific case studies, with potential to be applied to many other use cases. DeCOIL offers several advantages over existing methods, as it is direct, easy-to-use, generalizable, and scalable. With accompanying software, DeCOIL can be readily implemented to generate desired informed libraries.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Phase evolution and amorphous stability upon solid-state reaction in superlattice like Ge–Sb–Te combinatorial thin-film

In this paper, the superlattice-like (SLL) Ge–Sb–Te combinatorial thin films were prepared by using a high-throughput ion beam sputtering system. The phase evolution and amorphous stability of such films undergoing heat treatment as a function of the coating sequence and modulation period were systematically studied. The composition structure diagram was constructed via an automated process of data obtained by high-throughput synchrotron micro-X-ray diffraction and lab-based micro-X-ray fluorescence. Furthermore, the element distribution and microstructure in the depth direction of the SLL thin films were characterized with time-of-flight secondary ion mass spectrometry and transmission electron microscopy, respectively. These studies demonstrated that the coating sequence has a significant effect on the element distribution in the as-deposited SLL thin films and the structure of the final product upon solid-state reaction. Reducing the modulation period of the SLL thin film improves the stability of the amorphous Ge–Sb–Te phase. This work lays a solid foundation for the rational design of SLL Ge–Sb–Te thin films to improve their performance.

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

scMicrobe PTA: near complete genomes from single bacterial cells

Microbial genomes produced by standard single-cell amplification methods are largely incomplete. Here, we show that primary template-directed amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard multiple displacement amplification-based approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Impact of control signal phase noise on qubit fidelity

As qubit decoherence times are increased and readout technologies are improved, nonidealities in the drive signals, such as phase noise, are going to represent a growing limitation to the fidelity achievable at the end of complex control pulse sequencies. Here we study the impact on fidelity of phase noise affecting reference oscillators with the help of numerical simulations, which allow us to directly take into account the interaction between the phase fluctuations in the control signals and the evolution of the qubit state. Our method is based on the generation of phase noise realizations consistent with a given power spectral density, that are then applied to the pulse carrier in simulations, with Qiskit-Dynamics, of the qubit temporal evolution. By comparing the final state obtained at the end of a noisy pulse sequence with that in the ideal case and averaging over multiple noise realizations, we estimate the resulting degradation in fidelity, and exploiting an approximate analytical representation of a carrier affected by phase fluctuations, we discuss the contributions of the different spectral components of phase noise.

Barsotti, Agata [Pisa U.]↗
DOE OSTI2026

Methyl formate oxidation kinetics up to 100 atm

Methyl formate (MF, CH3OCHO), the simplest ester, is a representative oxygenated fuel with high oxygen content, and low sooting tendency. However, its oxidation behavior under high-pressure and intermediate-temperature conditions remains insufficiently understood, especially where low-temperature peroxy radical chemistry, methanol chemistry, and pressure-dependent reaction pathways play a critical role. In this study, MF oxidation experiments were conducted in the Princeton supercritical-pressure jet-stirred reactor (SP-JSR) at 20 and 100 atm over the temperature range of 400–950 K under both fuel-lean and fuel-rich conditions. Based on the experimental results, an updated HP-Mech was developed by incorporating previous MF sub-mechanisms, expanded low-temperature peroxy pathways, and evaluated pressure-dependent decomposition kinetics. The newly updated HP-Mech shows greatly improved performance in predicting the onset temperature, the key intermediate species fractions, methanol formation, and the progression of MF oxidation across all the experimental conditions. Path flux analysis indicates that MF consumption at the onset stage is dominated by H-abstraction at the methyl site, forming CH2OCHO radicals that lead to the formation and isomerization of O2CH2OCHO, driving low-temperature chain propagation. Moreover, H-abstraction at the formate site forms CH3OCO radicals that preferentially decompose to CH3, initiating the methanol formation pathway linked to CH3O2 and HO2 chemistry. At the same time, HO2 formation is strongly coupled to MF oxidation through multiple MF-derived radical pathways. HCO originates from MF oxidation and acts as a key coupling species linking fuel consumption to HO2 buildup, especially under high-pressure and intermediate-temperature conditions. In addition to this dominant channel, supplementary HO2 formation pathways involving CH3, CH3O, CH2OH, and CH3O2 reacting with O2 further connect methanol chemistry and oxygenated radical chemistry to the HO2 pool, indicating the central role of HO2 in governing MF oxidation. Sensitivity analysis identifies MF with OH/HO2/CH3O2 reactions and the HO2/H2O2/OH sequence as the key factors controlling reactivity in the high-pressure and intermediate-temperature regime. MF directly reacts with OH/HO2/CH3O2 to consume the fuel and produce reactive radicals like CH2OCHO and CH3OCO that undergo subsequent oxidation pathways. Moreover, HO2 recombination suppresses oxidation at lower temperatures, while thermal decomposition of H2O2 accelerates OH production and promotes fuel consumption as temperature increases. The direct formation of active OH from HO2 radicals further completes the mechanism, improving its prediction especially during the oxidation onset stage.

Low-temperature Chemistry↗
DOE OSTI2023

Gas-phase microactuation using kinetically controlled surface states of ultrathin catalytic sheets

Biological systems convert chemical energy into mechanical work by using protein catalysts that assume kinetically controlled conformational states. Synthetic chemomechanical systems using chemical catalysis have been reported, but they are slow, require high temperatures to operate, or indirectly perform work by harnessing reaction products in liquids (e.g., heat or protons). Here, we introduce a bioinspired chemical strategy for gas-phase chemomechanical transduction that sequences the elementary steps of catalytic reactions on ultrathin (<10 nm) platinum sheets to generate surface stresses that directly drive microactuation (bending radii of 700 nm) at ambient conditions (T = 20 °C; P total = 1 atm). When fueled by hydrogen gas and either oxygen or ozone gas, we show how kinetically controlled surface states of the catalyst can be exploited to achieve fast actuation (600 ms/cycle) at 20 °C. Finally, we also show that the approach can integrate photochemically controlled reactions and can be used to drive the reconfiguration of microhinges and complex origami- and kirigami-based microstructures.

42 ENGINEERING↗
DOE OSTI2023

Plug & play directed evolution of proteins with gradient-based discrete MCMC

Abstract A long-standing goal of machine-learning-based protein engineering is to accelerate the discovery of novel mutations that improve the function of a known protein. We introduce a sampling framework for evolving proteins in silico that supports mixing and matching a variety of unsupervised models, such as protein language models, and supervised models that predict protein function from sequence. By composing these models, we aim to improve our ability to evaluate unseen mutations and constrain search to regions of sequence space likely to contain functional proteins. Our framework achieves this without any model fine-tuning or re-training by constructing a product of experts distribution directly in discrete protein space. Instead of resorting to brute force search or random sampling, which is typical of classic directed evolution, we introduce a fast Markov chain Monte Carlo sampler that uses gradients to propose promising mutations. We conduct in silico directed evolution experiments on wide fitness landscapes and across a range of different pre-trained unsupervised models, including a 650 M parameter protein language model. Our results demonstrate an ability to efficiently discover variants with high evolutionary likelihood as well as estimated activity multiple mutations away from a wild type protein, suggesting our sampler provides a practical and effective new paradigm for machine-learning-based protein engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Pentose Phosphate Pathway Reactions in Photosynthesizing Cells

The pentose phosphate pathway (PPP) is divided into an oxidative branch that makes pentose phosphates and a non-oxidative branch that consumes pentose phosphates, though the non-oxidative branch is considered reversible. A modified version of the non-oxidative branch is a critical component of the Calvin–Benson cycle that converts CO2 into sugar. The reaction sequence in the Calvin–Benson cycle is from triose phosphates to pentose phosphates, the opposite of the typical direction of the non-oxidative PPP. The photosynthetic direction is favored by replacing the transaldolase step of the normal non-oxidative PPP with a second aldolase reaction plus sedoheptulose-1,7-bisphosphatase. This can be considered an anabolic version of the non-oxidative PPP and is found in a few situations other than photosynthesis. In addition to the strong association of the non-oxidative PPP with photosynthesis metabolism, there is recent evidence that the oxidative PPP reactions are also important in photosynthesizing cells. These reactions can form a shunt around the non-oxidative PPP section of the Calvin–Benson cycle, consuming three ATP per glucose 6-phosphate consumed. A constitutive operation of this shunt occurs in the cytosol and gives rise to an unusual labeling pattern of photosynthetic metabolites while an inducible shunt in the stroma may occur in response to stress.

13C↗
DOE OSTI2021

Kinetic Measurements in Heterogeneous Catalysis

This contribution is about the experimental determination of the rate of a heterogeneous catalytic reaction and the analysis of kinetic data. In this case, the reaction rate can be defined as the frequency at which the closed sequence of elementary steps transforming reactants into products is occurring. However, the rate of chemical reaction is not directly observed; rather, one records the rate of substance change. The rate of chemical reaction is calculated based on the rate of substance change and assumed stoichiometry of the reaction. In order to make this an intensive quantity (i.e., independent of the volume in which the reaction is performed, this frequency is divided by the latter or by a quantity proportional to the latter). As the focus of what follows is on heterogeneously catalyzed reactions this can be the catalyst mass, the catalyst surface area or the total number of active sites present in the reaction volume confined by the walls of a chemical reactor. For example, the rate may be reported as one of the following: Net rate of consumption/production of component i mol m-3 s-1 Net specific rate of consumption/production of component i mol kg-1 s-1 Rate of substance change (per unit volume of catalyst) mol m-3 s-1 Specific reaction rate (per unit mass of catalyst) mol kg-1 s-1 Net rate of consumption/production of component i mol m-2 s-1 The turnover frequency (TOF) is a characteristic originally introduced by Boudart [1] is expressed as TOF=R/G_tot , where R is the steady-state rate of reactant consumption or product generation and G_totis the total areal density of active sites [mol m-2]. The latter is typically taken from experimental chemisorption data obtained at low temperature. The advancement of a chemical reaction leads to changes in the amounts of reactants and products which, for a stoichiometric single reaction, are connected by the stoichiometric coefficients. In many cases, for the sake of simplicity, we will assume that the investigated reaction proceeds according to a single reaction pathway, so that the net production rate of a reaction component is directly proportional to the reaction rate and that the latter equals the net production rate of any of the involved components when divided by the appropriate stoichiometric coefficient. For complex stoichiometric reactions the net production rates of the involved components are linear combinations of the reaction rates. Different goals for kinetic measurements can be formulated: Catalyst testing, i.e., obtaining the kinetic dependences for the development and discrimination of efficient catalytic materials. Precise kinetic characterization of active materials via high-throughput screening is a part of such activity. Detailed kinetics, i.e., revealing the detailed mechanism of complex catalytic reaction via systematic kinetic studies, both steady-state and non-steady-state Industrial kinetics, i.e., obtaining data and relationships for describing and predicting the behaviour of catalytic reactors and processes at industrial scale Finally, mathematical modeling and analysis where kinetic measurements provide data for understanding complex kinetic phenomena, e.g., oscillations, non-linear self-organization, etc.

36 MATERIALS SCIENCE↗
DOE OSTI2022

AF2Complex predicts direct physical interactions in multimeric proteins with deep learning

Abstract Accurate descriptions of protein-protein interactions are essential for understanding biological systems. Remarkably accurate atomic structures have been recently computed for individual proteins by AlphaFold2 (AF2). Here, we demonstrate that the same neural network models from AF2 developed for single protein sequences can be adapted to predict the structures of multimeric protein complexes without retraining. In contrast to common approaches, our method, AF2Complex, does not require paired multiple sequence alignments. It achieves higher accuracy than some complex protein-protein docking strategies and provides a significant improvement over AF-Multimer, a development of AlphaFold for multimeric proteins. Moreover, we introduce metrics for predicting direct protein-protein interactions between arbitrary protein pairs and validate AF2Complex on some challenging benchmark sets and the E. coli proteome. Lastly, using the cytochrome c biogenesis system I as an example, we present high-confidence models of three sought-after assemblies formed by eight members of this system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An automated procedure built on MTEX for reconstructing deformation twin hierarchies from electron backscattered diffraction datasets of heavily twinned microstructures

Here we present a set of algorithms built on the MTEX and MATLAB graph toolboxes for automatic reconstruction of deformation twin hierarchies from Electron Backscatter Diffraction (EBSD) datasets with a focus on developing methods for heavily twinned microstructures (twin fractions >0.5). The algorithms address key issues arising at large strains, mainly: missing twin relationships, grouping of heavily deformed grain fragments into families of similar orientation originating from a single initial grain, identification of parent fragments for large twin volume fractions, and classification of families having twin relationships with multiple families. To facilitate the development of these algorithms, large-grained ultra-high purity α-Ti deformed in compression along two directions is investigated. Graphs are utilized to handle non-local geometric merging and to represent relationships throughout the reconstruction process. When determining if a grain fragment is from the undeformed microstructure, the combined metrics of the fragment's orientation volume fraction in the initial texture and the directed graph centrality measure of out-closeness (the number of nodes reached in a graph from a given node) are essential. To address automation in reconstructing the sequence of twinning and relating fragments originating from a single grain in the initial microstructure, the twin family tree is formulated as a minimum spanning tree emanating from the initial grain family. A scheme constructing the distances associated with twin relationship comprising the spanning tree is developed, and a novel quasi-directional Prim spanning tree algorithm is used to determine the twin family tree. The procedure is demonstrated to significantly improve the level of automation in reconstructing twin hierarchies in heavily twinned microstructure compared to other methodologies in literature. The procedure can readily be applied to analyses of twinning in metals, as well as provide an approach for routinely extracting twin statistics at larger deformation levels than previously possible. Significantly, the procedure is demonstrated to be capable of identifying third generation twinning in α-Ti microstructures.

36 MATERIALS SCIENCE↗