ASU Scholarship Showcase

Dehalococcoides mccartyi strains are of particular importance for bioremediation due to their unique capability of transforming perchloroethene (PCE) and trichloroethene (TCE) to non-toxic ethene, through the intermediates cis-dichloroethene (cis-DCE) and vinyl chloride (VC). Despite the widespread environmental distribution of Dehalococcoides, biostimulation sometimes fails to promote dechlorination beyond cis-DCE. In our study, microcosms established with garden soil and mangrove sediment also stalled at cis-DCE, albeit Dehalococcoides mccartyi containing the reductive dehalogenase genes tceA, vcrA and bvcA were detected in the soil/sediment inocula. Reductive dechlorination was not promoted beyond cis-DCE, even after multiple biostimulation events with fermentable substrates and a lengthy incubation.

However, transfers from microcosms stalled at cis-DCE yielded dechlorination to ethene with subsequent enrichment cultures containing up to 10⁹ Dehalococcoides mccartyi cells mL^-1. Proteobacterial classes which dominated the soil/sediment communities became undetectable in the enrichments, and methanogenic activity drastically decreased after the transfers. We hypothesized that biostimulation of Dehalococcoides in the cis-DCE-stalled microcosms was impeded by other microbes present at higher abundances than Dehalococcoides and utilizing terminal electron acceptors from the soil/sediment, hence, outcompeting Dehalococcoides for H₂. In support of this hypothesis, we show that garden soil and mangrove sediment microcosms bioaugmented with their respective cultures containing Dehalococcoides in high abundance were able to compete for H² for reductive dechlorination from one biostimulation event and produced ethene with no obvious stall. Overall, our results provide an alternate explanation to consolidate conflicting observations on the ubiquity of Dehalococcoides mccartyi and occasional stalling of dechlorination at cis-DCE; thus, bringing a new perspective to better assess biological potential of different environments and to understand microbial interactions governing bioremediation.

We investigate high-dimensional nonlinear dynamical systems exhibiting multiple resonances under adiabatic parameter variations. Our motivations come from experimental considerations where time-dependent sweeping of parameters is a practical approach to probing and characterizing the bifurcations of the system. The question is whether bifurcations so detected are faithful representations of the bifurcations intrinsic to the original stationary system. Utilizing a harmonically forced, closed fluid flow system that possesses multiple resonances and solving the Navier-Stokes equation under proper boundary conditions, we uncover the phenomenon of the early effect. Specifically, as a control parameter, e.g., the driving frequency, is adiabatically increased from an initial value, resonances emerge at frequency values that are lower than those in the corresponding stationary system. The phenomenon is established by numerical characterization of physical quantities through the resonances, which include the kinetic energy and the vorticity field, and a heuristic analysis based on the concept of instantaneous frequency. A simple formula is obtained which relates the resonance points in the time-dependent and time-independent systems. Our findings suggest that, in general, any true bifurcation of a nonlinear dynamical system can be unequivocally uncovered through adiabatic parameter sweeping, in spite of a shift in the bifurcation point, which is of value to experimental studies of nonlinear dynamical systems.

Most previous works on complete synchronization of chaotic oscillators focused on the one-channel interaction scheme where the oscillators are coupled through only one variable or a symmetric set of variables. Using the standard framework of master-stability function (MSF), we investigate the emergence of complex synchronization behaviors under all possible configurations of two-channel coupling, which include, for example, all possible cross coupling schemes among the dynamical variables. Utilizing the classic Rössler and Lorenz oscillators, we find a rich variety of synchronization phenomena not present in any previously extensively studied, single-channel coupling configurations. For example, in many cases two coupling channels can enhance or even generate synchronization where there is only weak or no synchronization under only one coupling channel, which has been verified in a coupled neuron system. There are also cases where the oscillators are originally synchronized under one coupling channel, but an additional synchronizable coupling channel can, however, destroy synchronization. Direct numerical simulations of actual synchronization dynamics verify the MSF-based predictions. Our extensive computation and heuristic analysis provide an atlas for synchronization of chaotic oscillators coupled through two channels, which can be used as a systematic reference to facilitate further research in this area.

Nonhyperbolicity, as characterized by the coexistence of Kolmogorov-Arnold-Moser (KAM) tori and chaos in the phase space, is generic in classical Hamiltonian systems. An open but fundamental question in physics concerns the relativistic quantum manifestations of nonhyperbolic dynamics. We choose the mushroom billiard that has been mathematically proven to be nonhyperbolic, and study the resonant tunneling dynamics of a massless Dirac fermion. We find that the tunneling rate as a function of the energy exhibits a striking "clustering" phenomenon, where the majority of the values of the rate concentrate on a narrow region, as a result of the chaos component in the classical phase space. Relatively few values of the tunneling rate, however, spread outside the clustering region due to the integrable component. Resonant tunneling of electrons in nonhyperbolic chaotic graphene systems exhibits a similar behavior. To understand these numerical results, we develop a theoretical framework by combining analytic solutions of the Dirac equation in certain integrable domains and physical intuitions gained from current understanding of the quantum manifestations of chaos. In particular, we employ a theoretical formalism based on the concept of self-energies to calculate the tunneling rate and analytically solve the Dirac equation in one dimension as well as in two dimensions for a circular-ring-type of tunneling systems exhibiting integrable dynamics in the classical limit. Because relatively few and distinct classical periodic orbits are present in the integrable component, the corresponding relativistic quantum states can have drastically different behaviors, leading to a wide spread in the values of the tunneling rate in the energy-rate plane. In contrast, the chaotic component has embedded within itself an infinite number of unstable periodic orbits, which provide far more quantum states for tunneling. Due to the nature of chaos, these states are characteristically similar, leading to clustering of the values of the tunneling rate in a narrow band. The appealing characteristic of our work is a demonstration and physical understanding of the "mixed" role played by chaos and regular dynamics in shaping relativistic quantum tunneling dynamics.

Background: Buffering to achieve pH control is crucial for successful trichloroethene (TCE) anaerobic bioremediation. Bicarbonate (HCO3−) is the natural buffer in groundwater and the buffer of choice in the laboratory and at contaminated sites undergoing biological treatment with organohalide respiring microorganisms. However, HCO3− also serves as the electron acceptor for hydrogenotrophic methanogens and hydrogenotrophic homoacetogens, two microbial groups competing with organohalide respirers for hydrogen (H2). We studied the effect of HCO3− as a buffering agent and the effect of HCO3−-consuming reactions in a range of concentrations (2.5-30 mM) with an initial pH of 7.5 in H2-fed TCE reductively dechlorinating communities containing Dehalococcoides, hydrogenotrophic methanogens, and hydrogenotrophic homoacetogens.

Results: Rate differences in TCE dechlorination were observed as a result of added varying HCO3− concentrations due to H2-fed electrons channeled towards methanogenesis and homoacetogenesis and pH increases (up to 8.7) from biological HCO3− consumption. Significantly faster dechlorination rates were noted at all HCO3− concentrations tested when the pH buffering was improved by providing 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) as an additional buffer. Electron balances and quantitative PCR revealed that methanogenesis was the main electron sink when the initial HCO3− concentrations were 2.5 and 5 mM, while homoacetogenesis was the dominant process and sink when 10 and 30 mM HCO3− were provided initially.

Conclusions: Our study reveals that HCO3− is an important variable for bioremediation of chloroethenes as it has a prominent role as an electron acceptor for methanogenesis and homoacetogenesis. It also illustrates the changes in rates and extent of reductive dechlorination resulting from the combined effect of electron donor competition stimulated by HCO3− and the changes in pH exerted by methanogens and homoacetogens.

Given a complex geospatial network with nodes distributed in a two-dimensional region of physical space, can the locations of the nodes be determined and their connection patterns be uncovered based solely on data? We consider the realistic situation where time series/signals can be collected from a single location. A key challenge is that the signals collected are necessarily time delayed, due to the varying physical distances from the nodes to the data collection centre. To meet this challenge, we develop a compressive-sensing-based approach enabling reconstruction of the full topology of the underlying geospatial network and more importantly, accurate estimate of the time delays. A standard triangularization algorithm can then be employed to find the physical locations of the nodes in the network. We further demonstrate successful detection of a hidden node (or a hidden source or threat), from which no signal can be obtained, through accurate detection of all its neighbouring nodes. As a geospatial network has the feature that a node tends to connect with geophysically nearby nodes, the localized region that contains the hidden node can be identified.

We develop a framework to uncover and analyse dynamical anomalies from massive, nonlinear and non-stationary time series data. The framework consists of three steps: preprocessing of massive datasets to eliminate erroneous data segments, application of the empirical mode decomposition and Hilbert transform paradigm to obtain the fundamental components embedded in the time series at distinct time scales, and statistical/scaling analysis of the components. As a case study, we apply our framework to detecting and characterizing high-frequency oscillations (HFOs) from a big database of rat electroencephalogram recordings. We find a striking phenomenon: HFOs exhibit on–off intermittency that can be quantified by algebraic scaling laws. Our framework can be generalized to big data-related problems in other fields such as large-scale sensor data and seismic data analysis.

Recent works revealed that the energy required to control a complex network depends on the number of driving signals and the energy distribution follows an algebraic scaling law. If one implements control using a small number of drivers, e.g. as determined by the structural controllability theory, there is a high probability that the energy will diverge. We develop a physical theory to explain the scaling behaviour through identification of the fundamental structural elements, the longest control chains (LCCs), that dominate the control energy. Based on the LCCs, we articulate a strategy to drastically reduce the control energy (e.g. in a large number of real-world networks). Owing to their structural nature, the LCCs may shed light on energy issues associated with control of nonlinear dynamical networks.

A challenging problem in network science is to control complex networks. In existing frameworks of structural or exact controllability, the ability to steer a complex network toward any desired state is measured by the minimum number of required driver nodes. However, if we implement actual control by imposing input signals on the minimum set of driver nodes, an unexpected phenomenon arises: due to computational or experimental error there is a great probability that convergence to the final state cannot be achieved. In fact, the associated control cost can become unbearably large, effectively preventing actual control from being realized physically. The difficulty is particularly severe when the network is deemed controllable with a small number of drivers. Here we develop a physical controllability framework based on the probability of achieving actual control. Using a recently identified fundamental chain structure underlying the control energy, we offer strategies to turn physically uncontrollable networks into physically controllable ones by imposing slightly augmented set of input signals on properly chosen nodes. Our findings indicate that, although full control can be theoretically guaranteed by the prevailing structural controllability theory, it is necessary to balance the number of driver nodes and control cost to achieve physical control.

Network reconstruction is a fundamental problem for understanding many complex systems with unknown interaction structures. In many complex systems, there are indirect interactions between two individuals without immediate connection but with common neighbors. Despite recent advances in network reconstruction, we continue to lack an approach for reconstructing complex networks with indirect interactions. Here we introduce a two-step strategy to resolve the reconstruction problem, where in the first step, we recover both direct and indirect interactions by employing the Lasso to solve a sparse signal reconstruction problem, and in the second step, we use matrix transformation and optimization to distinguish between direct and indirect interactions. The network structure corresponding to direct interactions can be fully uncovered. We exploit the public goods game occurring on complex networks as a paradigm for characterizing indirect interactions and test our reconstruction approach. We find that high reconstruction accuracy can be achieved for both homogeneous and heterogeneous networks, and a number of empirical networks in spite of insufficient data measurement contaminated by noise. Although a general framework for reconstructing complex networks with arbitrary types of indirect interactions is yet lacking, our approach opens new routes to separate direct and indirect interactions in a representative complex system.