Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sequential analysis”

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 595 records · Page 33

Development of the VOLT-A Shuttle experiment

The NASA Lewis Research Center (LeRC) is investigating potential problems associated with the operation of high voltage solar cell arrays in the space plasma environment. At high voltages, interactions between the solar array and the space plasma could result in unacceptable levels of electrical discharge (arcing) and/or parasitic losses (current drains from the array to the plasma). The objective of the Voltage Operating Limit Tests (VOLT-A) Shuttle bay experiment is to characterize space plasma/solar cell panel interactions in low earth orbit. VOLT-A consists of an experiment plate subassembly which contains four solar panels, an electronics subassembly and a Langmuir probe subassembly mounted on an MPESS carrier. During a given 8.25 hour data taking period (5-1/2 continuous orbits), the solar panels, which represent state-of-the-art solar cell technologies, will be sequentially subjected to bias voltages in steps ranging from minus 626 V to plus 313 V. Appropriate measurements will be made at each voltage to characterize arcing and parasitic losses. Corresponding measurements of the plasma environment (plasma density, electron temperature and neutral density) will also be made. Data will be recorded on an on-board tape recorder for subsequent data reduction and analysis.

Bifano, William J.↗

Enhanced orbit determination filter: Inclusion of ground system errors as filter parameters

The theoretical aspects of an orbit determination filter that incorporates ground-system error sources as model parameters for use in interplanetary navigation are presented in this article. This filter, which is derived from sequential filtering theory, allows a systematic treatment of errors in calibrations of transmission media, station locations, and earth orientation models associated with ground-based radio metric data, in addition to the modeling of the spacecraft dynamics. The discussion includes a mathematical description of the filter and an analytical comparison of its characteristics with more traditional filtering techniques used in this application. The analysis in this article shows that this filter has the potential to generate navigation products of substantially greater accuracy than more traditional filtering procedures.

Masters, W. C.↗

TDRSS Onboard Navigation System (TONS) flight qualification experiment

The National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) is currently developing an operational Tracking and Data Relay Satellite (TDRS) System (TDRSS) Onboard Navigation System (TONS) to provide realtime, autonomous, high-accuracy navigation products to users of TDRSS. A TONS experiment was implemented on the Explorer Platform/Extreme Ultraviolet Explorer (EP/EUVE) spacecraft, launched June 7, 1992, to flight qualify the TONS operational system using TDRSS forward-link communications services. This paper provides a detailed evaluation of the flight hardware, an ultrastable oscillator (USO) and Doppler extractor (DE) card in one of the TDRSS user transponders and the ground-based prototype flight software performance, based on the 1 year of TONS experiment operation. The TONS experiment results are used to project the expected performance of the TONS 1 operational system. TONS 1 processes Doppler data derived from scheduled forward-link S-band services using a sequential estimation algorithm enhanced by a sophisticated process noise model to provide onboard orbit and frequency determination and time maintenance. TONS 1 will be the prime navigation system on the Earth Observing System (EOS)-AM1 spacecraft, currently scheduled for launch in 1998. Inflight evaluation of the USO and DE short-term and long-term stability indicates that the performance is excellent. Analysis of the TONS prototype flight software performance indicates that realtime onboard position accuracies of better than 25 meters root-mean-square are achievable with one tracking contact every one to two orbits for the EP/EUVE 525-kilometer altitude, 28.5 degree inclination orbit. The success of the TONS experiment demonstrates the flight readiness of TONS to support the EOS-AM1 mission.

Gramling, C. J.↗

Novel software for analysis of root gravitropism: comparative response patterns of Arabidopsis wild-type and axr1 seedlings

In an earlier study (Evans, Ishikawa & Estelle 1994, Planta 194, 215-222) we used a video digitizer system to compare the kinetics of auxin action on root elongation in wild-type seedlings and seedlings of auxin response mutants of Arabidopsis thaliana (L.) Heynh. We have since modified the system software to allow determination of elongation on opposite sides of vertical or gravistimulated roots and to allow continuous measurement of the angle of orientation of sequential subsections of the root during the response. We used this technology to compare the patterns of differential growth that generate curvature in roots of the Columbia ecotype and in the mutants axr1-3, axr1-12 and axr2, which show reduced gravitropic responsiveness and reduced sensitivity to inhibition by auxin. The pattern of differential growth during gravitropism differed in roots of wild-type and axr1 seedlings. In wild-type roots, initial curvature resulted from differential inhibition of elongation in the distal elongation zone (DEZ). This was followed by an acceleration of elongation along the top side of the DEZ. In roots of axr1-3, curvature resulted from differential stimulation of elongation whereas in roots of axr1-12 the response was variable. Roots of axr2 did not exhibit gravitropic curvature. The observation that the pattern of differential growth causing curvature is dramatically altered by a change in sensitivity to auxin is consistent with the classical Cholodny-Went theory of gravitropism which maintains that differential growth patterns induced by gravistimulation are mediated primarily by gravi-induced shifts in auxin distribution. The new technology introduced with this report allows automated determination of stimulus response patterns in the small but experimentally popular roots of Arabidopsis.

Non-NASA Center↗

Data from a throughfall exclusion experiment: Fine root dynamics, morphology, chemistry, and AMF colonization across four lowland Panamanian forests

Fine roots regulate forest nutrient, carbon, and water cycling, yet their variation within and among tropical forests remains under-characterized. We quantified root productivity, disappearance, and stocks to 1 m using minirhizotron imaging, and we measured morphology, elemental composition [root carbon (C), root nitrogen (N), root phosphorus (P)], and arbuscular mycorrhizal fungi (AMF) colonization to 20 cm using ingrowth cores and sequential coring. Sampling took place in four distinct lowland Panamanian forests (32 plots; 8 per forest) from 2018 through 2022 under control and throughfall-exclusion (drought) treatments in the Panama Rainforest Changes with Experimental Drying (PARCHED) experiment.The dataset is presented as an Excel workbook with six tabs. The first tab is the data dictionary. Tab S1 contains ingrowth-core production and mortality, morphology and soil moisture. Tab S2 contains sequential-coring standing stocks with associated morphology and soil moisture. Tab S3 contains minirhizotron row data records to 1 m depth, including per-frame root length and diameter, normalized length metrics, and session timing. Tab S4 contains AMF colonization. Tab S5 contains fine-root chemistry at 0–10 cm, reporting %P, %C, %N, and C:N for samples collected via ingrowth cores and sequential-coring standing stocks. CSV mirrors for each tab are provided, and a KML file supplies coordinates for all 32 plots.Key variables span live and dead fine-root biomass (and coarse fractions where applicable), specific root length (SRL) and area (SRA), diameter, root tissue density (RTD), soil moisture, AMF colonization, root %N, %C, %P, and C:N, along with minirhizotron root length and diameter. Depth, season, treatment, and plot/site identifiers are included to support cross-tab integration and analysis from 0–100 cm (minirhizotron) and 0–20 cm (cores).Units are reported in-column and missing values are coded as NA. No special software is required to open or use the files (Excel, CSV, and KML compatible).

54 ENVIRONMENTAL SCIENCES↗

Understanding Solvent-Induced Glass Transition in Polymer Thin Films Using Absorption–Desorption Isotherms

The fundamental thermodynamic and mechanical underpinnings of polymer thin films exposed to solvent vapor are critical for the development of advanced nanolithography and high-performance coatings. This work investigates the solvent− polymer interactions of glassy thin films by using the solvent absorption−desorption isotherms. An analogous relationship to the Flory−Fox equation was observed between solvent−induced glass transition, swelling, Flory−Huggins interaction parameter, and molecular weight. Isothermal swelling measurements revealed that the glass transition trends are more robust in the absorption curve compared to desorption, contrary to previous reports. Excess osmotic pressure analysis of the isotherm provides a measure of the degree of physical aging in thin films annealed below the glass transition. This is further validated in the ordering of block copolymer (BCP) films annealed at low solvent activity. In agreement with the thermal analysis, free-surface plasticization effects become the most prominent below 100 nm. However, solvent annealing is largely dependent on solvent mass transport, as made evident by the strong dependence on solvent viscosity. From these observations, four general types of isotherms are identified that graphically capture distinct solvent−polymer interaction regimes. More broadly, these results inform solvent vapor annealing-induced self-assembly, sequential infiltration synthesis, membrane-based separations, adsorptive processes, and swelling-based responsive materials design.

Hendeniya, Nayanathara [Iowa State Univ., Ames, IA↗

Genetic algorithms as global random search methods

Genetic algorithm behavior is described in terms of the construction and evolution of the sampling distributions over the space of candidate solutions. This novel perspective is motivated by analysis indicating that that schema theory is inadequate for completely and properly explaining genetic algorithm behavior. Based on the proposed theory, it is argued that the similarities of candidate solutions should be exploited directly, rather than encoding candidate solution and then exploiting their similarities. Proportional selection is characterized as a global search operator, and recombination is characterized as the search process that exploits similarities. Sequential algorithms and many deletion methods are also analyzed. It is shown that by properly constraining the search breadth of recombination operators, convergence of genetic algorithms to a global optimum can be ensured.

Peck, Charles C.↗

Genetic algorithms as global random search methods

Genetic algorithm behavior is described in terms of the construction and evolution of the sampling distributions over the space of candidate solutions. This novel perspective is motivated by analysis indicating that the schema theory is inadequate for completely and properly explaining genetic algorithm behavior. Based on the proposed theory, it is argued that the similarities of candidate solutions should be exploited directly, rather than encoding candidate solutions and then exploiting their similarities. Proportional selection is characterized as a global search operator, and recombination is characterized as the search process that exploits similarities. Sequential algorithms and many deletion methods are also analyzed. It is shown that by properly constraining the search breadth of recombination operators, convergence of genetic algorithms to a global optimum can be ensured.

Peck, Charles C.↗

Elevation Change of the Southern Greenland Ice Sheet from Satellite Radar Altimeter Data

Long-term changes in the thickness of the polar ice sheets are important indicators of climate change. Understanding the contributions to the global water mass balance from the accumulation or ablation of grounded ice in Greenland and Antarctica is considered crucial for determining the source of the about 2 mm/yr sea-level rise in the last century. Though the Antarctic ice sheet is much larger than its northern counterpart, the Greenland ice sheet is more likely to undergo dramatic changes in response to a warming trend. This can be attributed to the warmer Greenland climate, as well as a potential for amplification of a global warming trend in the polar regions of the Northern Hemisphere. In collaboration with Drs. Curt Davis and Craig Kluever of the University of Missouri, we are using data from satellite radar altimeters to measure changes in the elevation of the Southern Greenland ice sheet from 1978 to the present. Difficulties with systematic altimeter measurement errors, particularly in intersatellite comparisons, beset earlier studies of the Greenland ice sheet thickness. We use altimeter data collected contemporaneously over the global ocean to establish a reference for correcting ice-sheet data. In addition, the waveform data from the ice-sheet radar returns are reprocessed to better determine the range from the satellite to the ice surface. At JPL, we are focusing our efforts principally on the reduction of orbit errors and range biases in the measurement systems on the various altimeter missions. Our approach emphasizes global characterization and reduction of the long-period orbit errors and range biases using altimeter data from NASA's Ocean Pathfinder program. Along-track sea-height residuals are sequentially filtered and backwards smoothed, and the radial orbit errors are modeled as sinusoids with a wavelength equal to one revolution of the satellite. The amplitudes of the sinusoids are treated as exponentially-correlated noise processes with a time-constant of six days. Measurement errors (e.g., altimeter range bias) are simultaneously recovered as constant parameters. The corrections derived from the global ocean analysis are then applied over the Greenland ice sheet. The orbit error and measurement bias corrections for different missions are developed in a single framework to enable robust linkage of ice-sheet measurements from 1978 to the present. In 1998, we completed our re-evaluation of the 1978 Seasat and 1985-1989 Geosat Exact Repeat Mission data. The estimates of ice thickness over Southern Greenland (south of 72N and above 2000 m) from 1978 to 1988 show large regional variations (+/-18 cm/yr), but yield an overall rate of +1.5 +/- 0.5 cm/yr (one standard error). Accounting for systematic errors, the estimate may not be significantly different from the null growth rate. The average elevation change from 1978 to 1988 is too small to assess whether the Greenland ice sheet is undergoing a long-term change.

Haines, Bruce J.↗

Organic amendments change soil organic C structure and microbial community but not total organic matter on sub-decadal scales

Organic C has many benefits for soil, but it is depleted by tillage and crop harvest, and especially so for biofuel crops. Accordingly, strategies such as partially retaining stover or planting a cover crop can help ameliorate the negative effect of C removal. We used a long-term field experiment to study the impacts of stover retention and planting cover crop on soil organic matter (SOM), its extractable components and the soil microbial community. SOM chemical composition characterization was determined by electrospray ionization (ESI) coupled with Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) in sequential water, methanol (MeOH), and chloroform (CHCl3) extracts. The characteristics of the soil bacterial community were measured by phospholipid fatty acid (PLFA), real-time quantitative PCR, and 16S rRNA gene sequence. The variations in total SOM content, total microbial biomass, and bacterial population were slight among treatments, but SOM chemical compounds, arbuscular mycorrhizal fungi (AMF) biomass, and bacterial structure changed significantly, and especially so in the coupled application of stover retention and cover crop. Specifically, stover retention enriched more lignin-like compounds in soil, whereas cover crop enriched more condensed hydrocarbons, and had more compounds with an aromaticity index (AI) >0.5. The bacterial community was not altered by the cover crop, but the corn stover retention increased the relative abundances of Myxococcales (Deltaproteobacteria) and decreased that of Actinobacteria. Redundancy analysis (RDA) further unveiled that the bacterial community in the stover treatments had significantly positive association with CHCl3-extracted chemical classes, i.e. unsaturated hydrocarbon and lipid with the coupled application (stover and cover crop), and lignin and protein with the single corn stover treatment. Taken together, our study shows how different C addition practices influence the molecular composition of SOM and the structure of soil microbial communities.

Xu, Jiangbin↗

Data assimilation in operator algebras

We develop an algebraic framework for sequential data assimilation of partially observed dynamical systems. In this framework, Bayesian data assimilation is embedded in a nonabelian operator algebra, which provides a representation of observables by multiplication operators and probability densities by density operators (quantum states). In the algebraic approach, the forecast step of data assimilation is represented by a quantum operation induced by the Koopman operator of the dynamical system. Moreover, the analysis step is described by a quantum effect, which generalizes the Bayesian observational update rule. Projecting this formulation to finite-dimensional matrix algebras leads to computational schemes that are i) automatically positivity-preserving and ii) amenable to consistent data-driven approximation using kernel methods for machine learning. Moreover, these methods are natural candidates for implementation on quantum computers. Applications to the Lorenz 96 multiscale system and the El Niño Southern Oscillation in a climate model show promising results in terms of forecast skill and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗

Shoulder fillet effects in strength distributions of microelectromechanical system components

The failure forces and fracture strengths of polysilicon microelectromechanical system (MEMS) components in the form of stepped tensile bars with shoulder fillets were measured using a sequential failure chain methodology. Approximately 150 specimens for each of four fillet geometries with different stress concentration factors were tested. The resulting failure force and strength distributions of the four geometries were related by a common sidewall flaw population existing within different effective stressed lengths. The failure forces, strengths, and flaw population were well described by a weakest-link based analytical framework. Finite element analysis was used to verify body-force based expressions for the stress concentration factors and to provide insight into the variation of specimen effective length with fillet geometry. Monte Carlo simulations of flaw size and location, based on the strength measurements, were also used to provide insight into fillet shape and size effects. The successful description of the shoulder fillet specimen strengths provides further empirical support for application of the strength and flaw framework in MEMS fabrication and design optimization.

36 MATERIALS SCIENCE↗

A Concept for the Deployment of a Large Lunar Crater Radio Telescope Using Teams of Tethered Robots

Kilometer-scale craters on the far side of the Moon have unique potential as future locations for large radio telescopes, which can observe the universe at wavelengths and frequencies (> 10 m, < 30 MHz) not possible with conventional Earth or orbital-based approaches. Distinct advantages of building a Lunar Crater Radio Telescope (LCRT) on the far side include i) isolation from radio noise due to the Earth’s ionosphere, orbiting satellites, and the Sun, ii) days of uninterrupted dark/cold sky viewing during lunar night, and iii) terrain geometry naturally suited for constructing the largest mesh antenna structure in the Solar System. A key challenge to constructing LCRT on the Moon is related to the complexity of deploying a 1-km diameter antenna and hanging receiver within a lunar crater whose diameter, depth, and slope are 3-5 km, 1 km, and ~30 degrees respectively. In this paper, we first evaluate the trade space for deploying a large, complex structure within a crater, and then provide a more detailed concept evaluation of our favored approach, which employs coordinated teams of tethered rovers to extract and suspend a folded antenna from a lander at the base of a crater. NASA’s Jet Propulsion Laboratory in collaboration with California Institute of Technology (Caltech) have developed a novel robotic system for accessing extremely steep terrains; the Axel rover is a two-wheeled rugged terrain vehicle that is supported by an electro-mechanical tether that provides power, data, and tensile support from a top-side anchor location. Recently, a pair of Axel robots have been used in a DuAxel configuration that allows for four-wheel driving and repeated passive anchoring at different locations. The DuAxel system has unique advantages for deploying an LCRT antenna, including the ability to deploy from a lander near a crater, drive a distance to the crater rim to deploy an Axel, and later, retract the deployed Axel in order to sequentially lift up sections of the antenna. Our proposed concept involves delivering a packaged antenna and receiver to the bottom-center of a crater floor on a lander, then later sending a team of multiple DuAxel rovers to retrieve guide wires from the lander, which are pulled to the top of the crater. We explore this concept in detail and provide some initial quantitative analysis to demonstrate the feasibility of our system with respect to the spatial and mass properties of the antenna as juxtaposed to DuAxel capabilities. Finally, we outline next steps towards validating our concept on the way to a future lunar deployment opportunity.

Hallinan, Greg↗

A Reversed Impedance-Based Stability Criterion for IBR Grids: Preprint

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗

Comprehensive process and environmental impact analysis of integrated DBD plasma steam methane reforming

Utilization of electricity generated from renewable sources to obtain hydrogen, H 2 , is of critical importance to decrease the overall carbon footprint. Here in this work, integration of a dielectric discharge barrier (DBD) plasma reactor to convert low calorific value gas, such as landfill gas or coal mine gas into hydrogen, into the existing steam methane reforming (SMR) technology was evaluated using process design considerations. In particular, a DBD-enhanced catalytic SMR reactor was modeled to operate at near atmospheric pressure and 500 °C sequentially with the conventional reformer to obtain ~ 65 kmol/hr H2 for distributed production. This allowed decreasing the size of the conventional reformer albeit at the increased overall electricity consumption. Calculated process economics showed that only at an electricity cost of less than $0.004/kWh does the hybrid DBD plasma process derived H 2 price become competitive with that of the conventional SMR. A Life Cycle Assessment framework was used to compare environmental impacts from the conventional SMR, hybrid DBD SMR and hybrid DBD SMR utilizing only onshore wind-derived electricity. Larger environmental impacts in the plasma reformer were obtained due to the use of electricity for the plasma reforming operation, which was modeled as coming from the typical U.S. grid mix. Utilizing only 100% wind-derived electricity provided certain environmental benefits, except for the ecotoxicity impact where the wind power scenario modeled here only reduced ecotoxicity impacts associated with electricity by 30%.

08 HYDROGEN↗

The effects of sequential hydrothermal-mechanical refining pretreatment on cellulose structure changes and sugar recoveries

The recalcitrance of lignocellulosic biomass necessitates an efficient pretreatment protocol for operating a successful cellulosic biorefinery. It is critical to improve cellulose accessibility for hydrolysis and fermentation by altering the plant cell wall’s physical structure and chemical composition. Sequential hydrothermal-mechanical refining pretreatment (HMR) allows efficient recovery of cellulosic sugars without utilizing any hazardous chemicals. HMR has been successfully applied to Liberty switchgrass, a bioenergy cultivar released by the USDA, and now it is being applied to oilcane, a recently developed transgenic sugarcane variety engineered to accumulate lipids in its vegetative tissues. Sugar yields of oilcane bagasse (OCB) and switchgrass (SG) treated with HMR are 96.4% and 75.4%, respectively. This study sought to correlate cellulosic sugar yields with structural changes within the cell wall caused by HMR on two distinct bioenergy crops. Simon’s staining technique for the specific surface area analysis showed that HMR increased the specific surface area of pretreated biomass residues by 80–112%. In addition, ATR-FTIR was performed to determine the effects of HMR on physical structures based on the total crystallinity index (TCI) and hydrogen bonding intensity (HBI). Irrespective of biomass type, HMR decreased the initial crystalline cellulose contents of untreated biomass residues by 3.5% and reduced TCI and HBI by 7–13%. Here, the study found that sugar yields were negatively correlated to reducing values of hydrogen bonding intensity, crystalline cellulose content, and total crystallinity index.

09 BIOMASS FUELS↗

Mechanical continuity and reversible chromosome disassembly within intact genomes removed from living cells

Chromatin is thought to be structurally discontinuous because it is packaged into morphologically distinct chromosomes that appear physically isolated from one another in metaphase preparations used for cytogenetic studies. However, analysis of chromosome positioning and movement suggest that different chromosomes often behave as if they were physically connected in interphase as well as mitosis. To address this paradox directly, we used a microsurgical technique to physically remove nucleoplasm or chromosomes from living cells under isotonic conditions. Using this approach, we found that pulling a single nucleolus or chromosome out from interphase or mitotic cells resulted in sequential removal of the remaining nucleoli and chromosomes, interconnected by a continuous elastic thread. Enzymatic treatments of interphase nucleoplasm and chromosome chains held under tension revealed that mechanical continuity within the chromatin was mediated by elements sensitive to DNase or micrococcal nuclease, but not RNases, formamide at high temperature, or proteases. In contrast, mechanical coupling between mitotic chromosomes and the surrounding cytoplasm appeared to be mediated by gelsolin-sensitive microfilaments. Furthermore, when ion concentrations were raised and lowered, both the chromosomes and the interconnecting strands underwent multiple rounds of decondensation and recondensation. As a result of these dynamic structural alterations, the mitotic chains also became sensitive to disruption by restriction enzymes. Ion-induced chromosome decondensation could be blocked by treatment with DNA binding dyes, agents that reduce protein disulfide linkages within nuclear matrix, or an antibody directed against histones. Fully decondensed chromatin strands also could be induced to recondense into chromosomes with pre-existing size, shape, number, and position by adding anti-histone antibodies. Conversely, removal of histones by proteolysis or heparin treatment produced chromosome decondensation which could be reversed by addition of histone H1, but not histones H2b or H3. These data suggest that DNA, its associated protein scaffolds, and surrounding cytoskeletal networks function as a structurally-unified system. Mechanical coupling within the nucleoplasm may coordinate dynamic alterations in chromatin structure, guide chromosome movement, and ensure fidelity of mitosis.

NASA Discipline Cell Biology↗

The hierarchical organization of autocatalytic reaction networks and its relevance to origin of life

Abiogenesis has long been suspected to require that chemical reaction networks contain multiple autocatalytic cores, but little is known about what features of these networks allow for the gradual accretion of complexity. To identify realistic scenarios for the emergence of life-like properties, we develop the concept of a seed-dependent autocatalytic system (SDAS), which is a subnetwork that can autocatalytically self-maintain given a flux of food, but cannot be initiated by food alone. Rather, initiation of SDASs requires the transient introduction of chemical “seeds”. We show that, depending on the topological relationship of SDASs in a chemical reaction network, a food-driven system can accrete complexity in a historically contingent manner, governed by rare seeding events. We develop new algorithms for detecting and analyzing SDASs in chemical reaction databases and describe parallels between multi-SDAS networks and biological ecosystems. Applying our algorithms to both an abiotic reaction network and a biochemical one, each driven by a set of simple food chemicals, we detect SDASs that are organized as trophic tiers, of which the higher tier can be seeded by relatively simple chemicals if the lower tier is already activated. This indicates that sequential activation of trophically organized SDASs by seed chemicals that are not much more complex than what already exist could be a mechanism of gradual complexification from relatively simple abiotic reactions to more complex life-like systems. Interestingly, in both reaction networks, higher-tier SDASs include chemicals that might alter emergent features of chemical systems and could serve as early targets of selection. Our analysis provides computational tools for analyzing very large chemical/biochemical reaction networks and suggests new approaches to studying abiogenesis in the lab.

Zhen Peng↗