Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Preliminary Phase Field Computational Model Development

This interim report presents progress towards the development of meso-scale models of magnetic behavior that incorporate microstructural information. Modeling magnetic signatures in irradiated materials with complex microstructures (such as structural steels) is a significant challenge. The complexity is addressed incrementally, using the monocrystalline Fe (i.e., ferrite) film as model systems to develop and validate initial models, followed by polycrystalline Fe films, and by more complicated and representative alloys. In addition, the modeling incrementally addresses inclusion of other major phases (e.g., martensite, austenite), minor magnetic phases (e.g., carbides, FeCr precipitates), and minor nonmagnetic phases (e.g., Cu precipitates, voids). The focus of the magnetic modeling is on phase-field models. The models are based on the numerical solution to the Landau-Lifshitz-Gilbert equation. From the computational standpoint, phase-field modeling allows the simulation of large enough systems that relevant defect structures and their effects on functional properties like magnetism can be simulated. To date, two phase-field models have been generated in support of this work. First, a bulk iron model with periodic boundary conditions was generated as a proof-of-concept to investigate major loop effects of single versus polycrystalline bulk iron and effects of single non-magnetic defects. More recently, to support the experimental program herein using iron thin films, a new model was generated that uses finite boundary conditions representing surfaces and edges. This model has provided key insights into the domain structures observed in magnetic force microscopy (MFM) measurements. Simulation results for single crystal thin-film iron indicate the feasibility of the model for determining magnetic domain wall thickness and mobility in an externally applied field. Because the phase-field model dimensions are limited relative to the size of most specimens used in experiments, special experimental methods were devised to create similar boundary conditions in the iron films. Preliminary MFM studies conducted on single and polycrystalline iron films with small sub-areas created with focused ion beam have correlated quite well qualitatively with phase-field simulations. However, phase-field model dimensions are still small relative to experiments thus far. We are in the process of increasing the size of the models and decreasing specimen size so both have identical dimensions. Ongoing research is focused on validation of the phase-field model. Validation is being accomplished through comparison with experimentally obtained MFM images (in progress), and planned measurements of major hysteresis loops and first order reversal curves. Extrapolation of simulation sizes to represent a more stochastic bulk-like system will require sampling of various simulations (i.e., with single non-magnetic defect, single magnetic defect, single grain boundary, single dislocation, etc.) with distributions of input parameters. These outputs can then be compared to laboratory magnetic measurements and ultimately to simulate magnetic Barkhausen noise signals.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Exploring the Potential of Fulvalene Dimetals as Platforms for Molecular Solar Thermal Energy Storage: Computations, Syntheses, Structures, Kinetics, and Catalysis

A study of the scope and limitations of varying the ligand framework around the dinuclear core of FvRu(2) in its function as a molecular solar thermal energy storage framework is presented. It includes DFT calculations probing the effect of substituents, other metals, and CO exchange for other ligands on Delta H-storage. Experimentally, the system is shown to be robust in as much as it tolerates a number of variations, except for the identity of the metal and certain substitution patterns. Failures include 1,1',3,3'-tetra-tert-butyl (4), 1,2,2',3'-tetraphenyl (9), diiron (28), diosmium (24), mixed iron-ruthenium (27), dimolybdenum (29), and di-tungsten (30) derivatives. An extensive screen of potential catalysts for the thermal reversal identified AgNO3-SiO2 as a good candidate, although catalyst decomposition remains a challenge.

Borjesson, K↗
DOE OSTI2012

Simulation of current-filament dynamics and relaxation in the Pegasus Spherical Tokamak

Nonlinear numerical computation is used to investigate the relaxation of non-axisymmetric current-channels from washer-gun plasma sources into 'tokamak-like' plasmas in the Pegasus toroidal experiment [Eidietis et al. J. Fusion Energy 26, 43 (2007)]. Resistive MHD simulations with the NIMROD code [Sovinec et al. Phys. Plasmas 10(5), 1727-1732 (2003)] utilize ohmic heating, temperature-dependent resistivity, and anisotropic, temperature-dependent thermal conduction corrected for regions of low magnetization to reproduce critical transport effects. Adjacent passes of the simulated current-channel attract and generate strong reversed current sheets that suggest magnetic reconnection. With sufficient injected current, adjacent passes merge periodically, releasing axisymmetric current rings from the driven channel. The current rings have not been previously observed in helicity injection for spherical tokamaks, and as such, provide a new phenomenological understanding for filament relaxation in Pegasus. After large-scale poloidal-field reversal, a hollow current profile and significant poloidal flux amplification accumulate over many reconnection cycles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Skyrmion pinball and directed motion on obstacle arrays

Abstract We examine skyrmions interacting with a square array of obstacles under ac drives applied in one or two directions. For a single direction of ac driving, we find that the Magnus force in conjunction with the obstacle interactions can create elliptical skyrmion orbits of increasing size, leading to localized phases, chaotic phases, and translating or ratcheting orbits. Under two ac drives that are out of phase by 90° and applied in two directions, the skyrmions form localized commensurate orbits that encircle an integer number of obstacles, similar to the electron pinball effect observed for electrons in antidot lattices. As a function of ac amplitude, Magnus force strength, and obstacle size, we find that chaotic scattering regimes and directed motion can emerge even in the absence of asymmetry in the substrate. The directed motion follows different symmetry axes of the periodic substrate, and we observe a variety of reversed ratchet effects. The Magnus force in the skyrmion system produces a significantly larger number of directed motion regimes than are exhibited by overdamped systems. We discuss how these results could be used to move skyrmions in a controlled way for possible applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Rapid estimation of earthquake locations using waveform traveltimes

SUMMARY We introduce a new approach for locating earthquakes using arrival times derived from waveforms. The most costly computational step of the algorithm scales as the number of stations in the active seismographic network. In this approach, a variation on existing grid search methods, a series of full waveform simulations are conducted for all receiver locations, with sources positioned successively at each station. The traveltime field over the region of interest is calculated by applying a phase picking algorithm to the numerical wavefields produced from each simulation. An event is located by subtracting the stored traveltime field from the arrival time at each station. This provides a shifted and time-reversed traveltime field for each station. The shifted and time-reversed fields all approach the origin time of the event at the source location. The mean or median value at the source location thus approximates the event origin time. Measures of dispersion about this mean or median time at each grid point, such as the sample standard error and the average deviation, are minimized at the correct source position. Uncertainty in the event position is provided by the contours of standard error defined over the grid. An application of this technique to a synthetic data set indicates that the approach provides stable locations even when the traveltimes are contaminated by additive random noise containing a significant number of outliers and velocity model errors. It is found that the waveform-based method out-performs one based upon the eikonal equation for a velocity model with rapid spatial variations in properties due to layering. A comparison with conventional location algorithms in both a laboratory and field setting demonstrates that the technique performs at least as well as existing techniques.

58 GEOSCIENCES↗
DOE OSTI2017

Structure of a pentameric virion-associated fiber with a potential role in Orsay virus entry to host cells

Despite the wide use of Caenorhabditis elegans as a model organism, the first virus naturally infecting this organism was not discovered until six years ago. The Orsay virus and its related nematode viruses have a positive-sense RNA genome, encoding three proteins: CP, RdRP, and a novel δ protein that shares no homology with any other proteins. δ can be expressed either as a free δ or a CP-δ fusion protein by ribosomal frameshift, but the structure and function of both δ and CP-δ remain unknown. Using a combination of electron microscopy, X-ray crystallography, computational and biophysical analyses, here we show that the Orsay δ protein forms a ~420-Å long, pentameric fiber with an N-terminal α-helical bundle, a β-stranded filament in the middle, and a C-terminal head domain. The pentameric nature of the δ fiber has been independently confirmed by both mass spectrometry and analytical ultracentrifugation. Recombinant Orsay capsid containing CP-δ shows protruding long fibers with globular heads at the distal end. Mutant viruses with disrupted CP-δ fibers were generated by organism-based reverse genetics. These viruses were found to be either non-viable or with poor infectivity according to phenotypic and qRT-PCR analyses. Furthermore, addition of purified δ proteins to worm culture greatly reduced Orsay infectivity in a sequence-specific manner. Based on the structure resemblance between the Orsay CP-δ fiber and the fibers from reovirus and adenovirus, we propose that CP-δ functions as a cell attachment protein to mediate Orsay entry into worm intestine cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Understanding Reaction Networks through Controlled Approach to Equilibrium Experiments Using Transient Methods

We report a combined experimental/theoretical approach to studying heterogeneous gas/solid catalytic processes using low-pressure pulse response experiments achieving a controlled approach to equilibrium that combined with quantum mechanics (QM)-based computational analysis provides information needed to reconstruct the role of the different surface reaction steps. We demonstrate this approach using model catalysts for ammonia synthesis/decomposition. Polycrystalline iron and cobalt are studied via low-pressure TAP (temporal analysis of products) pulse response, with the results interpreted through reaction free energies calculated using QM on Fe-BCC(110), Fe-BCC(111), and Co-FCC(111) facets. In TAP experiments, simultaneous pulsing of ammonia and deuterium creates a condition where the participation of reactants and products can be distinguished in both forward and reverse reaction steps. This establishes a balance between competitive reactions for D* surface species that is used to observe the influence of steps leading to nitrogen formation as the nitrogen product remains far from equilibrium. Here, the approach to equilibrium is further controlled by introducing delay timing between NH 3 and D 2 which allows time for surface reactions to evolve before being driven in the reverse direction from the gas phase. The resulting isotopic product distributions for NH 2 D, NHD 2 , and HD at different temperatures and delay times and NH 3 /D 2 pulsing order reveal the role of the N 2 formation barrier in controlling the surface concentration of NH x * species, as well as providing information on the surface lifetimes of key reaction intermediates. Conclusions derived for monometallic materials are used to interpret experimental results on a more complex and active CoFe bimetallic catalyst.

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

The physical foundation of the reconnection electric field

Magnetic reconnection is a key charged particle transport and energy conversion process in environments ranging from astrophysical systems to laboratory plasmas [Yamada et al., Rev. Mod. Phys. 82, 603–664 (2010)]. Magnetic reconnection facilitates plasma transport by establishing new connections of magnetic flux tubes, and it converts, often explosively, energy stored in the magnetic field to kinetic energy of charged particles [J. L. Burch and J. F. Drake, Am. Sci. 97, 392–299 (2009)]. The intensity of the magnetic reconnection process is measured by the reconnection electric field, which regulates the rate of flux tube connectivity changes. The change of magnetic connectivity occurs in the current layer of the diffusion zone, where the plasma transport is decoupled from the transport of magnetic flux. Here we report on computer simulations and analytic theory to provide a self-consistent understanding of the role of the reconnection electric field, which extends substantially beyond the simple change of magnetic connections. Rather, we find that the reconnection electric field is essential to maintain the current density in the diffusion region, which would otherwise be dissipated by a set of processes. Natural candidates for current dissipation are the average convection of current carriers away from the reconnection region by the outflow of accelerated particles, or the average rotation of the current density by the magnetic field reversal in the vicinity. Instead, we show here that the current dissipation is the result of thermal effects, underlying the statistical interaction of current-carrying particles with the adjacent magnetic field. We find that this interaction serves to redirect the directed acceleration of the reconnection electric field to thermal motion. This thermalization manifests itself in form of quasi-viscous terms in the thermal energy balance of the current layer. This collisionless viscosity, found in the pressure evolution equation, dominates near the x-line. These quasi-viscous terms act to increase the average thermal energy. Our predictions regarding current and thermal energy balance are readily amenable to exploration in the laboratory or by satellite missions, in particular, by NASA's Magnetospheric Multiscale mission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Reliability analysis of the F-8 digital fly-by-wire system

The F-8 Digital Fly-by-Wire (DFBW) flight test program intended to provide the technology for advanced control systems, giving aircraft enhanced performance and operational capability is addressed. A detailed analysis of the experimental system was performed to estimated the probabilities of two significant safety critical events: (1) loss of primary flight control function, causing reversion to the analog bypass system; and (2) loss of the aircraft due to failure of the electronic flight control system. The analysis covers appraisal of risks due to random equipment failure, generic faults in design of the system or its software, and induced failure due to external events. A unique diagrammatic technique was developed which details the combinatorial reliability equations for the entire system, promotes understanding of system failure characteristics, and identifies the most likely failure modes. The technique provides a systematic method of applying basic probability equations and is augmented by a computer program written in a modular fashion that duplicates the structure of these equations.

Brock, L. D.↗
DOE OSTI2016

Y-box-binding protein-1 (YB-1) promotes cell proliferation, adhesion and drug resistance in diffuse large B-cell lymphoma

YB-1 is a multifunctional protein, which has been shown to correlate with resistance to treatment of various tumor types. This study investigated the expression and biologic function of YB-1 in diffuse large B-cell lymphoma (DLBCL). Immunohistochemical analysis showed that the expression statuses of YB-1 and pYB-1{sup S102} were reversely correlated with the clinical outcomes of DLBCL patients. In addition, we found that YB-1 could promote the proliferation of DLBCL cells by accelerating the G1/S transition. Ectopic expression of YB-1 could markedly increase the expression of cell cycle regulators cyclin D1 and cyclin E. Furthermore, we found that adhesion of DLBCL cells to fibronectin (FN) could increase YB-1 phosphorylation at Ser102 and pYB-1{sup S102} nuclear translocation. In addition, overexpression of YB-1 could increase the adhesion of DLBCL cells to FN. Intriguingly, we found that YB-1 overexpression could confer drug resistance through cell-adhesion dependent and independent mechanisms in DLBCL. Silencing of YB-1 could sensitize DLBCL cells to mitoxantrone and overcome cell adhesion-mediated drug resistance (CAM-DR) phenotype in an AKT-dependent manner. - Highlights: • The expression statuses of YB-1 and pYB-1{sup S102} are reversely correlated with outcomes of DLBCL patients. • YB-1 promotes cell proliferation by accelerating G1/S transition in DLBCL. • YB-1 confers drug resistance to mitoxantrone in DLBCL.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

High pressure suppression of plasticity due to an overabundance of shear embryo formation

Abstract High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is known to occur in a variety of these materials, the governing dynamics of the mechanisms are not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands, “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high-strain rate application materials such as those for high explosives and hypersonic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2013

New intrinsic mechanism on gum-like superelasticity of multifunctional alloys

Ti-Nb-based Gum Metals exhibit extraordinary superelasticity with ultralow elastic modulus, superior strength and ductility, and a peculiar dislocation-free deformation behavior, most of which challenge existing theories of crystal strength. Additionally, this kind of alloys actually displays even more anomalous mechanical properties, such as the non-linear superelastic behavior, accompanied by a pronounced tension-to-compression asymmetry, and large ductility with a low Poisson’s ratio. Two main contradictory arguments exist concerning the deformation mechanisms of those alloys, i.e., formation of reversible nanodisturbance and reversible martensitic transformation. Herein we used the in-situ synchrotron high-energy X-ray scattering technique to reveal the novel intrinsic physical origin of all anomalous mechanical properties of the Ti-24Nb-4Zr-8Sn-0.10O alloy, a typical gum-like metal. Our experiments provide direct evidence on two different kinds of interesting, stress-induced, reversible nanoscale martensitic transitions, i.e., the austenitic regions with B2 structure transform to a0 martensite and those with BCC structure transform to d martensite.

COMPUTATIONAL METHODS↗
DOE OSTI2017

Examination of the α-ω two-phase shock-induced microstructure in zirconium and titanium

Omega (ω) phase is formed in alpha (α) zirconium during dynamic loading and can be retained in the recovered material. The pathway for the α to ω change (or the reverse transformation) is not well understood. Zirconium was shock-loaded and the resulting two-phase microstructure was examined. Electron backscatter diffraction (EBSD) was used to characterize the orientation relationships and habit planes between phases to understand the pathways between α and ω phases and compare to Molecular Dynamics (MD) simulations. Based on key microstructural features, a significant amount of α phase appears to have originated from the reverse transformation from ω-Zr on unloading. Results of microstructural analysis will be discussed, along with implications toward phase transformation pathways.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A Scalable Algorithm for Radiative Heat Transfer Using Reverse Monte Carlo Ray Tracing

Radiative heat transfer is an important mechanism in a class of challenging engineering and research problems. A direct all-to-all treatment of these problems is prohibitively expensive on large core counts due to pervasive all- to-all MPI communication. The massive heat transfer problem arising from the next generation of clean coal boilers being modeled by the Uintah framework has radiation as a dominant heat transfer mode. Reverse Monte Carlo ray tracing (RMCRT) can be used to solve for the radiative-flux divergence while accounting for the effects of participating media. The ray tracing approach used here replicates the geometry of the boiler on a multi-core node and then uses an all-to-all communication phase to distribute the results globally. The cost of this all-to-all is reduced by using an adaptive mesh approach in which a fine mesh is only used locally, and a coarse mesh is used elsewhere. Here, a model for communication and computation complexity is used to predict performance of this new method. We show this model is consistent with observed results and demonstrate excellent strong scaling to 262K cores on the DOE Titan system on problem sizes that were previously computationally intractable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Covalent inhibition of hAChE by organophosphates causes homodimer dissociation through long-range allosteric effects

Acetylcholinesterase (EC 3.1.1.7), a key acetylcholine-hydrolyzing enzyme in cholinergic neurotransmission, is present in a variety of states in situ, including monomers, C-terminally disulfide-linked homodimers, homotetramers, and up to three tetramers covalently attached to structural subunits. Could oligomerization that ensures high local concentrations of catalytic sites necessary for efficient neurotransmission be affected by environmental factors? Using small-angle X-ray scattering (SAXS) and cryo-EM, we demonstrate that homodimerization of recombinant monomeric human acetylcholinesterase (hAChE) in solution occurs through a C-terminal four-helix bundle at micromolar concentrations. We show that diethylphosphorylation of the active serine in the catalytic gorge or isopropylmethylphosphonylation by the R P enantiomer of sarin promotes a 10-fold increase in homodimer dissociation. We also demonstrate the dissociation of organophosphate (OP)-conjugated dimers is reversed by structurally diverse oximes 2PAM, HI6, or RS194B, as demonstrated by SAXS of diethylphosphoryl-hAChE. However, binding of oximes to the native ligand-free hAChE, binding of high-affinity reversible ligands, or formation of an S P -sarin-hAChE conjugate had no effect on homodimerization. Dissociation monitored by time-resolved SAXS occurs in milliseconds, consistent with rates of hAChE covalent inhibition. OP-induced dissociation was not observed in the SAXS profiles of the double-mutant Y337A/F338A, where the active center gorge volume is larger than in wildtype hAChE. These observations suggest a key role of the tightly packed acyl pocket in allosterically triggered OP-induced dimer dissociation, with the potential for local reduction of acetylcholine-hydrolytic power in situ. Computational models predict allosteric correlated motions extending from the acyl pocket toward the four-helix bundle dimerization interface 25 Å away.

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

Bi{sub 1−x}Nb{sub x}O{sub 1.5+x} (x=0.0625, 0.12) fast ion conductors: Structures, stability and oxide ion migration pathways

A combined experimental and computational study of Bi{sub 1−x}Nb{sub x}O{sub 1.5+x} (x=0.0625 and 0.12) has been carried out using laboratory X-ray, neutron and electron diffraction, impedance measurements and ab-initio molecular dynamics. We demonstrate that Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625}, previously reported to adopt a cubic fluorite-type superstructure, can form two different polymorphs depending on the synthetic method: a metastable cubic phase is produced by quenching; while slower cooling yields a stable material with a tetragonal √2×√2×1 superstructure, which undergoes a reversible phase transition into the cubic form at ~680 °C on subsequent reheating. Neutron diffraction reveals that the tetragonal superstructure arises mainly from ordering in the oxygen sublattice, with Bi and Nb remaining disordered, although structured diffuse scattering observed in the electron diffraction patterns suggests a degree of short-range ordering. Both materials are oxide ion conductors. On thermal cycling, Bi{sub 0.88}Nb{sub 0.12}O{sub 1.62} exhibits a decrease in conductivity of approximately an order of magnitude due to partial transformation into the tetragonal phase, but still exhibits conductivity comparable to yttria-stabilised zirconia (YSZ). Ab-initio molecular dynamics simulations performed on Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} show that oxide ion diffusion occurs by O{sup 2−} jumps between edge- and corner-sharing OM{sub 4} groups (M=Bi, Nb) via tetrahedral □M{sub 4} and octahedral □M{sub 6} vacancies. - Graphical abstract: Oxide ion migration in tetragonal Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} occurs by O{sup 2−} jumps between edge- and corner-sharing OM{sub 4} groups (M=Bi, Nb) via tetrahedral M{sub 4} and octahedral M{sub 6} vacancies. - Highlights: • Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} adopts a tetragonal √2×√2×1 fluorite superstructure. • Superstructure is due to ordering in the O-sublattice, with Bi/Nb disordered. • Bi{sub 0.9375}Nb{sub 0.0625}O{sub 1.5625} is a good oxide ion conductor. • O{sup 2−} jump between OM{sub 4} groups (M=Bi, Nb) via tetrahedral and octahedral vacancies.

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

Highly active Au/δ-MoC and Au/β-Mo 2 C catalysts for the low-temperature water gas shift reaction: effects of the carbide metal/carbon ratio on the catalyst performance

In this paper, the water gas shift (WGS) reaction catalyzed by orthorhombic β-Mo 2 C and cubic δ-MoC surfaces with and without Au clusters supported thereon has been studied by means of a combination of sophisticated experiments and state-of-the-art computational modeling. Experiments evidence the importance of the metal/carbon ratio on the performance of these systems, where Au/δ-MoC is presented as a suitable catalyst for WGS at low temperatures owing to its high activity, selectivity (only CO 2 and H 2 are detected), and stability (oxycarbides are not observed). Periodic density functional theory-based calculations show that the supported Au clusters and the Au/δ-MoC interface do not take part directly in water dissociation but their presence is crucial to switch the reaction mechanism, drastically decreasing the effect of the reverse WGS reaction and favoring the WGS products desorption, thus leading to an increase in CO 2 and H 2 production. Finally, the present results clearly display the importance of the Mo/C ratio and the synergy with the admetal clusters in tuning the activity and selectivity of the carbide substrate.

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

Network Analysis of Breast Cancer Progression and Reversal Using a Tree-Evolving Network Algorithm

The HMT3522 progression series of human breast cells have been used to discover how tissue architecture, microenvironment and signaling molecules affect breast cell growth and behaviors. However, much remains to be elucidated about malignant and phenotypic reversion behaviors of the HMT3522-T4-2 cells of this series. We employed a ‘‘pan-cell-state’’ strategy, and analyzed jointly microarray profiles obtained from different state-specific cell populations from this progression and reversion model of the breast cells using a tree-lineage multi-network inference algorithm, Treegl. We found that different breast cell states contain distinct gene networks. The network specific to non-malignant HMT3522-S1 cells is dominated by genes involved in normal processes, whereas the T4-2-specific network is enriched with cancer-related genes. The networks specific to various conditions of the reverted T4-2 cells are enriched with pathways suggestive of compensatory effects, consistent with clinical data showing patient resistance to anticancer drugs. We validated the findings using an external dataset, and showed that aberrant expression values of certain hubs in the identified networks are associated with poor clinical outcomes. Thus, analysis of various reversion conditions (including non-reverted) of HMT3522 cells using Treegl can be a good model system to study drug effects on breast cancer.

59 BASIC BIOLOGICAL SCIENCES↗