Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamical evolution and stability”

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 163 records · Page 9

Structure, Bonding, and Vibrational Dynamics of a Triamine High Energy Density Material under Pressure

High energy density materials have complex intermolecular interactions which influence their stability and performance. We used a combination of synchrotron X-ray diffraction, synchrotron infrared spectroscopy, and Raman vibrational spectroscopy, supplemented by density functional theory calculations, to probe pressure-induced changes in structure and intermolecular interactions of 1H,4'H-[3,3'-bis(1,2,4-triazole)]-4',5,5'-triamine as a model high energy density material up to 40 GPa. We find that compression of the triamine is accompanied by increased intermolecular interactions that give rise to an interesting evolution of the structure and bonding with pressure. Analysis of the equation of state determined from the X-ray diffraction indicates a change in compression mechanism near 19 GPa consistent with changes in vibrational spectra that provide evidence for a structural rearrangement associated with changes in hydrogen bonding near that pressure. As a result, the overall compressional behavior calculated theoretically agrees with that observed experimentally though differences are found that indicate the need for improved treatment of the intermolecular interactions including hydrogen bonding under pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic response of a freely rotating butterfly valve in the advanced test reactor – dynamic fluid-body interaction modeling

To regulate primary coolant flow in the Advanced Test Reactor (ATR), a butterfly valve was installed between the primary coolant pumps and the reactor core. If the mechanical connection between the valve's disk and its shaft ever fails, the disk will rotate freely. Rapid disk rotation may induce pressure surges that could damage upstream pipes. In the present work, the rotational trajectory and pressure evolution during a disk free-rotation scenario were analyzed via the dynamic fluid-body interaction (DFBI) approach in STAR-CCM+, with the movement of a solid being driven by the forces and moment/torque imposed by its surrounding fluid. Assuming a large initial opening angle, the disk accelerates slowly, but swiftly passes the closed position. As a result of the sudden valve closure, a pressure surge occurs in the upstream pipes, exceeding their maximum allowable pressure. Furthermore, the disk does not stabilize at the closed position but continually oscillates due to the unsteady nature of the coolant flow. Because of the significant and continuous water hammer effect, a fix to the butterfly valve is being implemented to prevent rapid valve closure due to failure at the valve's disk/shaft mechanical connection.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Vortex splitting in two-dimensional fluids and non-neutral electron plasmas with smooth vorticity profiles

Initially, elliptical, quasi-two-dimensional (2D) fluid vortices can split into multiple pieces if the aspect ratio is sufficiently large due to the growth and saturation of perturbations known as Love modes on the vortex edge. Presented here are experiments and numerical simulations, showing that the aspect ratio threshold for vortex splitting is significantly higher for vortices with realistic, smooth edges than that predicted by a simple “vortex patch” model, where the vorticity is treated as piecewise constant inside a deformable boundary. The experiments are conducted by exploiting the E × B drift dynamics of collisionless, pure electron plasmas in a Penning–Malmberg trap, which closely model 2D vortex dynamics due to an isomorphism between the Drift–Poisson equations describing the plasmas and the Euler equations describing ideal fluids. The simulations use a particle-in-cell method to model the evolution of a set of point vortices. The aspect ratio splitting threshold ranges up to about twice as large as the vortex patch prediction and depends on the edge vorticity gradient. This is thought to be due to spatial Landau damping, which decreases the vortex aspect ratio over time and, thus, stabilizes the Love modes. Near the threshold, asymmetric splitting events are observed in which one of the split products contains much less circulation than the other. These results are relevant to a wide range of quasi-2D fluid systems, including geophysical fluids, astrophysical disks, and drift-wave eddies in tokamak plasmas.

Physics↗

Tri‐Metallic Catalyst for Oxygen Evolution Reaction Enables Continuous Operation of Anion Exchange Membrane Electrolyzer at 1A cm −2 for Hundreds of Hours

Abstract Although numerous efforts are made to synthesize active electrocatalysts for green hydrogen production; catalyst stability, and facile synthesis to scale up the production are still challenging. Herein, the production of novel non‐PGM catalysts for the oxygen reduction reaction (OER) in an alkaline aqueous medium is reported, which is based on the synthesis of a trimetallic metal–organic framework (MOF) precursors. Fine‐tuning of the composition of the metal centers (Ni, Co, and Fe) shows a great effect on OER activity after the MOF undergoes dynamic chemical and structural transformations under OER conditions. In situ characterization reveals the origin of OER activity enhancement as metals’ oxidation state increases, inducing compressive mechanical strain on metal centers, enhancing the electronic conductivity through the formation of oxygen vacancies, and stronger metal–oxygen covalency. Catalysts are used in membrane electrode assembly (MEA) setup within an industrial full‐cell anion exchange membrane electrolyzer (AEMEC), showing a stable performance for 550 h without noticeable decay at 750 and 1000 mA cm −2 industrial level current densities.

Chemistry↗

The Role of Hydrophobic Nodes in the Dynamics of Class A β-Lactamases

Class A β-lactamases are known for being able to rapidly gain broad spectrum catalytic efficiency against most β-lactamase inhibitor combinations as a result of elusively minor point mutations. The evolution in class A β-lactamases occurs through optimisation of their dynamic phenotypes at different timescales. At long-timescales, certain conformations are more catalytically permissive than others while at the short timescales, fine-grained optimisation of free energy barriers can improve efficiency in ligand processing by the active site. Free energy barriers, which define all coordinated movements, depend on the flexibility of the secondary structural elements. The most highly conserved residues in class A β-lactamases are hydrophobic nodes that stabilize the core. To assess how the stable hydrophobic core is linked to the structural dynamics of the active site, we carried out adaptively sampled molecular dynamics (MD) simulations in four representative class A β-lactamases (KPC-2, SME-1, TEM-1, and SHV-1). Using Markov State Models (MSM) and unsupervised deep learning, we show that the dynamics of the hydrophobic nodes is used as a metastable relay of kinetic information within the core and is coupled with the catalytically permissive conformation of the active site environment. Finally, our results collectively demonstrate that the class A enzymes described here, share several important dynamic similarities and the hydrophobic nodes comprise of an informative set of dynamic variables in representative class A β-lactamases.

59 BASIC BIOLOGICAL SCIENCES↗

Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation

Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Field-resolved observation of exciton coherence in a van der Waals magnet

The emergence of coherence among electronic quasiparticles underlies collective quantum phenomena from superconductivity to superradiance. In semiconductors, exciton coherence is generally thought to decay rapidly due to scattering and dephasing, limiting its persistence on ultrafast timescales. Here, in this work, we demonstrate a light-field-driven mechanism that creates and stabilizes exciton coherence in the layered antiferromagnet CrSBr. We directly record the coherent optical field emitted by excitons and track in real time how a deterministic phase, imprinted by the excitation laser, drives incoherent excitons to synchronize into a collective state. This ensemble remains phase coherent for more than 2 ps, whereas its resonance energy undergoes an ultrafast modulation mediated by spin and lattice interactions. The time-resolved field evolution indicates that the multiple peaks seen in conventional spectra originate from a single excitonic resonance subject to dynamic energy modulation. Our findings establish optical phase imprinting as a mechanism to control and sustain collective order in semiconducting magnets, bridging light-driven dynamics with excitonic and magnetic correlations in layered quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Investigation of the Impact of Flow on MHD Perturbations with the NIMROD Code (Final Report)

This work will enable improved understanding of tokamak edge flows and impurity-species dynamics and their impact on MHD stability as it relates to ELMs, RMPs and QH-mode by leveraging computational developments of the NIMROD code. Experimentally, it is established that the tokamak flows and impurities species can have a substantial effect on MHD stability. This is particularly true for the edge plasma where, for example, large flow shear is correlated with the occurrence of the Quiescent H-mode (QH- mode) state as opposed to operation with edge localized modes (ELMs). Understanding plasma flows is challenging because of the many physical effects that come into play: particle orbits loss, neutrals dynamics, and interaction between multiple ion species with decoupled motion. We will develop a model that incorporates multiple ion species into the MHD framework that includes self-consistent magnetic-field evolution. With this model we will understand how the presence of impurities impacts the tokamak edge-pedestal flows. Finally, we will quantify the impact of the new multiple species modeling on 3D QH-mode simulations. This work will have broader impacts on all MHD studies when multiple species are present (as is the case in modern tokamaks), as well as contributing to the edge and transport communities who also are interested in the studies of multiple-ion-species flows in the edge pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural adaptation of oxygen tolerance in 4-hydroxybutyrl-CoA dehydratase, a key enzyme of archaeal carbon fixation

Abstract Autotrophic microorganisms that convert inorganic carbon into organic matter were key players in the evolution of life on Earth. As the early atmosphere became oxygenated, microorganisms needed to develop mechanisms for oxygen protection, especially those relying on enzymes containing oxygen-sensitive metal clusters (e.g., Fe-S). Here we investigated how 4-hydroxybutyryl-CoA dehydratase (4HBD) - the key enzyme of the 3-hydroxypropionate/4-hydroxybutyrate (HP/HB) cycle for CO 2 -fixation - adapted as conditions shifted from anoxic to oxic. 4HBD is found in both anaerobic bacteria and aerobic ammonia-oxidizing archaea (AOA). The oxygen-sensitive bacterial 4HBD and oxygen-tolerant archaeal 4HBD share 59 % amino acid identity. To examine the structural basis of oxygen tolerance in archaeal 4HBD, we determined the atomic resolution structure of the enzyme. Two tunnels providing access to the canonical [4Fe-4S] cluster in oxygen-sensitive bacterial 4HBD were closed with four conserved mutations found in all aerobic AOA and other archaea. Further biochemical experiments and molecular dynamics simulations support our findings that restricting access to the active site is the key to oxygen tolerance, explaining how active site evolution drove a major evolutionary transition. Significance statement Autotrophy (primary production) was the first life strategy on Earth. Before photosynthesis (using solar energy to fix carbon dioxide), life relied on chemical reactions for energy. These chemosynthetic reactions are present in all domains of life, including archaea possessing the most energy-efficient carbon fixation pathway - the 3-hydroxypropionate/4-hydroxybutyrate cycle. This efficiency results from enzyme modifications, including enhanced enzyme stability and catalysis of multiple reactions. We reveal the first structure of aerobic 4-hydroxybutyryl-CoA dehydratase (4HBD) from ammonia-oxidizing archaea. These archaea are among the most abundant organisms on the planet, and their 4HBD active site evolved oxygen tolerance to support aerobic metabolism. This modification can provide further insight into enzyme evolution on early earth, as photosynthesis developed and began oxygenating the atmosphere.

DeMirci, Hasan (ORCID:0000000291355397)↗

Tools Assessing Performance (FY2022 Quarter 4 Report)

We were finally successful in acquiring TTU’s 200-meter tower data for the duration of the Isolated Building Wake Experiment (IBWE). This data will provide an independent measurement of the undisturbed wind profile including heights far above those that are typical of meteorological sensor networks. These measurements will characterize the mean and turbulence throughout the atmospheric surface layer, which will be particularly useful for validation of highfidelity computational fluid dynamics models. This dataset can be shared freely among the TAP team but not outside of it as TTU sells this data to the wind power industry. The journal article overview has been submitted to LANL’s publication review and will be submitted to the Journal of Wind Engineering and Industrial Aerodynamics once the review is complete. We are nearing completion on another two articles. The first discusses analyses of the evolution of the measurements of mean wind speed along the centerline of the building as a function of atmospheric stability. Rather than using individual time steps, all of the wind speed measurements are normalized using the upwind measurement at 2 building heights (H) and filtered by wind direction and atmospheric stability. The second discusses the development of the diffusive wake model and uses the averaged profiles from the first article to validate the behavior of the model. The portion of the section below on the diffusive wake model discussing comparisons with the averaged profiles, will be included in the validation section of this paper.

17 WIND ENERGY↗

Cyclodextrin-Derived Porous Liquids Enabled by In Situ Solvation Shell Formation

Porous liquids (PLs) represent a unique platform for molecular separations by combining permanent porosity with liquid-phase mobility. However, it remains a formidable challenge to construct and stabilize PLs with sub-5 Å pores using readily available porous host and liquid media. Here, we report the construction of cyclodextrin (CD)-derived PLs enabled by in situ solvation shell formation. The acid–base neutralization reaction between CD and an organic base was leveraged to generate a thin ionic solvation shell around the CD host, effectively liquefying CD and preventing its segregation in the liquid base medium while preserving accessible molecular-scale cavities. Spectroscopic analysis, neutron scattering, density functional theory calculations, and molecular dynamics simulations collectively confirm the structural evolution and existence of abundant internal porosity in PLs. The unique architectures of CD-derived PLs enable highly selective encapsulation of fluorinated alkanes and significantly enhanced uptake of inert gases. This facile and generalizable strategy enables construction of high-quality PLs with engineered ultramicroporosity to facilitate molecular separations.

Cavities↗

Flash Formation of I-Rich Clusters during Multistage Halide Segregation Studied in MAPbI 1.5 Br 1.5

Wide band gap mixed halide perovskites such as MAPb(I 1–x Br x ) 3 sparked great research interest because of their outstanding optoelectronic properties, ease of fabrication, and band gap tunability. Their application thus far is however limited by light-induced halide segregation in which microscopic clusters with a high iodide content are formed and act as recombination centers. The key mechanism(s) underlying this halide segregation process are still debated. Here, we present a study on the photoluminescence evolution in MAPb(I 1.5 Br 1.5 ) perovskites with varying microstructures under constant illumination at room temperature and at elevated temperature. Our findings reveal a more complicated picture of the segregation mechanism occurring in three stages instead of two as commonly reported. The process starts with a flash formation of I-rich domains. Following is a rapid blue shift before the gradual and typically observed red shift occurs. The evolution of the three stages is fully reversible in the dark and is also present at elevated temperatures (50 °C). We explain the existence of multiple stages during light-induced halide segregation by natural compositional fluctuations of the halides and the formation of halide clusters with a dynamically changing distribution in I–Br content. The variation in the I–Br ratio depends on the grain size and film heterogeneity. Overall, these findings add further details in the quest of unraveling the underlying segregation mechanism(s), which need to be identified to stabilize halides in wide band gap perovskites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Antecedent Hydrologic Conditions Reflected in Stream Lithium Isotope Ratios During Storms

Antecedent hydrological conditions are recorded through the evolution of dissolved lithium isotope signatures (δ 7 Li) by juxtaposing two storm events in an upland watershed subject to a Mediterranean climate. Discharge and δ 7 Li are negatively correlated in both events,but mean δ 7 Li ratios and associated ranges of variation are distinct between them. We apply a previously developed reactive transport model (RTM) for the site to these event-scale flow perturbations, but observed shifts in stream δ 7 Li are not reproduced. To reconcile the stability of the subsurface solute weathering profile with our observations of dynamic stream δ 7 Li signatures, we couple the RTM to a distribution of fluid transit times that evolve based on storm hydrographs. The approach guides appropriate flux-weighting of fluid from the RTM over a range of flow path lengths, or equivalently fluid residence times. This flux-weighted RTM approach accurately reproduces dynamic storm δ 7 Li-discharge patterns distinguished by the antecedent conditions of the watershed.

58 GEOSCIENCES↗

Pt Particles on a Dynamic TiO 2 Support in Near-Ambient Conditions−Disentangling Size, Pressure, and Support Effects

Platinum particles on reducible oxides are known to form complex and highly dynamic catalyst systems at elevated pressures and temperatures, often adopting active structures that differ from those found at room temperature and under ultrahigh vacuum (UHV). Here, we study the oxidation and structural evolution of subnanometer Pt clusters and nanoparticles supported on rutile TiO 2 (110) across an oxygen pressure range from UHV to 0.1 mbar, using near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS), scanning tunneling microscopy (STM) under UHV and NAP conditions, and low-energy ion scattering (LEIS). Our results reveal distinct differences in oxidation behavior and thermal stability between Pt nanoparticles and clusters, which are further modulated by the support stoichiometry and oxygen pressure. Small Pt clusters become oxidized even at room temperature but are susceptible to accelerated sintering in 0.1 mbar O 2 at elevated temperatures. In contrast, well-crystallized Pt nanoparticles on near-stoichiometric TiO 2 show weaker oxidation. On a reduced, defective TiO 2 support, Pt instead quickly becomes deeply buried by new titania layers, which are formed during support reoxidation. This process appears to result primarily from interactions of the support with the gas phase, unlike the classical, selflimited encapsulation that is induced by the strong metal−support interaction (SMSI). Finally, we address the full complexity of real catalysts in a direct side-by-side comparison of the single-crystalline model system with a Pt-loaded TiO 2 powder catalyst (P25). We conclude that the stoichiometry of the model supports must be carefully chosen and controlled to accurately reproduce the expected state of powder supports during redox reactions.

metal nanoparticles↗

Visualizing the SEI formation between lithium metal and solid-state electrolyte

The solid electrolyte interphase (SEI) is regarded as the most important factor affecting the durability of lithium-metal anode in all-solid-state batteries (ASSBs). Despite its significance, the nucleation and growth mechanism of SEI is not yet well understood. Here, we elucidate the thermodynamics and kinetics governing SEI formation at the Li|β-Li 3 PS 4 interface at the atomic scale via thermodynamic phase equilibrium analysis and machine-learning-potential-assisted molecular dynamics (MD) simulations. An accurate moment tensor potential using the machine-learning method is developed for a reactive model of Li|β-Li 3 PS 4 . This potential enabled us to perform large-scale MD simulations with the model size expanded to the experimental dimensions (~40 nm) while maintaining the same level of accuracy as density functional theory calculations. The results reveal a four-stage evolution process at the Li|β-Li 3 PS 4 interface, namely (i) fast ion diffusion, (ii) nucleation, (iii) Li 2 S growth, and (iv) stabilization. Notably, we demonstrate that the SEI can be categorized into crystalline and amorphous regions. The simulated SEI thickness, structure, and composition closely match experimental findings, validating the accuracy of the MD simulations. We further disclose the significant impact of ion diffusion kinetic limitations on the phase formation and crystallization of interfacial products. Furthermore, we shed light on the detailed potential energy (PE) distribution of lithium along the direction perpendicular to the interface. Finally, this information is crucial for better understanding interfacial ion mobility.

25 ENERGY STORAGE↗

Investigation of the Use of Dynamic Probabilistic Risk Assessment Methodologies for Identifying Digital I&C System Common Cause Failures

Digital Instrumentation and Control (I&C) systems have a key role in nuclear power plants in the upgrade of aging analog systems. Digital systems improve plant safety and reliability through features such as increased hardware reliability and stability and improved failure detection capability. There is no consensus on which of the current probabilistic risk assessment methods are most suitable for use in the reliability analysis of digital I&C systems. While the traditional event-tree/fault-tree (ET/FT) approach is still used for their reliability modeling, there are concerns regarding this approach in properly accounting for dynamic interactions among system components since potentially significant dependencies among failure events may not be identified and/or their likelihood may not be properly quantified. Dynamic methodologies are expected to provide a much more accurate representation of probabilistic evolution of the I&C systems in time due to their capability to more properly account for complex interactions than the static approach. The applicability of dynamic PRA methodologies for digital I&C system is investigated using the criteria presented in the NUREG/CR-6901, and the comparisons made in NUREG/CR-6901 are updated in light of the latest studies. The Dynamic Event Tree (DET) approach has been identified as one of the top dynamic methods when evaluated against the requirements for the reliability modeling of digital I&C systems. The DET method is a strong candidate for integration into existing PRA studies, as it bears many similarities to the traditional ET approach. In this study, the DET approach has been applied to the Plant Protection System of the APR1400 design, and the results are compared to results from its available traditional ET/FT analysis. Possible approaches to evaluate and quantify the effects of common cause failures on system safety using dynamic methods are also examined.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗

Visualizing Size-Dependent Dynamics of CeO 2-δ {100}-Supported CoO x Nanoparticles Under CO 2 Hydrogenation Conditions

Carbon dioxide is a major greenhouse gas. In order to optimize processes focused on its chemical valorization, one needs detailed information about the effects of CO 2 and/or CO 2 /H 2 mixtures on the structure and morphology of metal/oxide catalysts. Here, in this study, the evolution of a catalyst with cobalt supported on CeO 2 -cube nanostructures under CO 2 hydrogenation conditions was investigated by using a set of in situ characterization techniques (X-ray absorption fine structure, X-ray diffraction, diffuse reflectance infrared Fourier transform spectroscopy, and environmental transmission electron microscopy (TEM)). The {100} facets of the ceria support displayed an unexpectedly high stability due to strong interactions with the aggregates of cobalt oxide. A significant influence of interfacial bonding between CoO x and CeO 2-δ {100} is evident through a clear preference in the orientation of CoO x nanoparticles (NPs) with respect to the substrate. For initially reduced Co/CeO 2 -cube nanostructures, a kinetically controlled oxidation of cobalt upon the introduction of CO 2 was observed during the early stages of CO 2 hydrogenation. Environmental TEM revealed the size-dependent morphological behavior of cobalt oxide NPs due to strong interactions with the CeO 2 {100} surface. When the environment was switched from H 2 to a mixture of H 2 and CO 2 at 250 °C, small CoO x NPs (in the largest dimension < 2.5 nm) rapidly transformed from a pyramidal three-dimensional (3D) form to a planar, monatomic layer attached to the concurrently oxidized CeO 2-δ {100} surface. This maximizes the number of sites available for the binding of CO 2 or reaction intermediates. The shape transformation reflected the oxophilic character of cobalt and strong metal–support interactions. The removal of CO 2 from the gas phase led to a reduction of the cobalt oxide NPs by hydrogen and a reversible two-dimensional → 3D transformation. In contrast, no significant morphological changes, apart from further oxidation, were observed for big CoO x NPs (in the largest dimension > 3 nm). These trends are not seen for nanoparticles of noble metals. The observed morphological and structural changes in the small CoO x NPs affected the stability of reaction intermediates and modified the selectivity of the CoO x /CeO 2 catalyst system for methane production.

36 MATERIALS SCIENCE↗