Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface 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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Probing electronic dead layers in homoepitaxial n-SrTiO 3 (001) films

We combine state-of-the-art oxide epitaxial growth by hybrid molecular beam epitaxy with transport, x-ray photoemission, and surface diffraction, along with classical and first-principles quantum mechanical modeling to investigate the nuances of insulating layer formation in otherwise high-mobility homoepitaxial n-SrTiO 3 (001) films. Our analysis points to charge immobilization at the buried n-SrTiO 3 /undoped SrTiO 3 (001) interface as well as within the surface contamination layer resulting from air exposure as the drivers of electronic dead-layer formation. As Fermi level equilibration occurs at the surface and the buried interface, charge trapping reduces the sheet carrier density (n 2D ) and renders the n-STO film insulating if n 2D falls below the critical value for the metal-to-insulator transition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scanning tunneling microscope tip-induced formation of Bi bilayers on Bi 2 Te 3

We report the formation of Bi(111) bilayer islands and crater structures on Bi 2 Te 3 (111) surfaces induced by voltage pulses from a scanning tunneling microscope tip. Pulses above a threshold voltage (+3 V) produce craters ∼0.5μm in diameter, similar to the size of the tip. Redeposited material self-assembles into a network of atomically ordered islands with a lattice constant identical to the underlying Bi2Te3 surface. The island size monotonically decreases over several micrometers from the pulse site, until the pristine Bi 2 Te 3 surface is recovered. We assign these islands to Bi bilayer based on atomic resolution images, analysis of step heights, and tunneling spectroscopy. Here, the dependence of bilayer formation on bias polarity and the evidence for defect diffusion together suggest a mechanism driven by the interplay of field evaporation and tunneling-current-induced Joule heating.

Bi bilayer↗
DOE OSTI2024

On‐Surface Synthesis and Determination of the Open‐Shell Singlet Ground State of Tridecacene**

Abstract The character of the electronic structure of acenes has been the subject of longstanding discussion. However, convincing experimental evidence of their open‐shell character has so far been missing. Here, we present the on‐surface synthesis of tridecacene molecules by thermal annealing of octahydrotridecacene on a Au(111) surface. We characterized the electronic structure of the tridecacene by scanning probe microscopy, which reveals the presence of an inelastic signal at 126 meV. We attribute the inelastic signal to spin excitation from the singlet diradical ground state to the triplet excited state. To rationalize the experimental findings, we carried out many‐body ab initio calculations as well as model Hamiltonians to take into account the effect of the metallic substrate. Moreover, we provide a detailed analysis of how the dynamic electron correlation and virtual charge fluctuation between the molecule and metallic surface reduces the singlet‐triplet band gap. Thus, this work provides the first experimental confirmation of the magnetic character of tridecacene.

Chemistry↗
DOE OSTI2022

Domain architecture and catalysis of the Staphylococcus aureus fatty acid kinase

Fatty acid kinase (Fak) is a two-component enzyme that generates acyl-phosphate for phospholipid synthesis. Fak consists of a kinase domain protein (FakA) that phosphorylates a fatty acid enveloped by a fatty acid binding protein (FakB). The structural basis for FakB function has been established, but little is known about FakA. Here, we used limited proteolysis to define three separate FakA domains: the amino terminal FakA_N, the central FakA_L, and the carboxy terminal FakA_C. The isolated domains lack kinase activity, but activity is restored when FakA_N and FakA_L are present individually or connected as FakA_NL. The X-ray structure of the monomeric FakA_N captures the product complex with ADP and two Mg 2+ ions bound at the nucleotide site. The FakA_L domain encodes the dimerization interface along with conserved catalytic residues Cys240, His282, and His284. AlphaFold analysis of FakA_L predicts the catalytic residues are spatially clustered and pointing away from the dimerization surface. Furthermore, the X-ray structure of FakA_C shows that it consists of two subdomains that are structurally related to FakB. Analytical ultracentrifugation demonstrates that FakA_C binds FakB, and site-directed mutagenesis confirms that a positively charged wedge on FakB meshes with a negatively charged groove on FakA_C. Finally, small angle X-ray scattering analysis is consistent with freely rotating FakA_N and FakA_C domains tethered by flexible linkers to FakA_L. These data reveal specific roles for the three independently folded FakA protein domains in substrate binding and catalysis.

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

The interplay between surface facet and reconstruction on isopropanol conversion over SrTiO 3 nanocrystals

Strontium titanate (SrTiO 3 ) is an extensively investigated perovskite for various applications due to its optical, electrical and chemical properties. To gain an in-depth understanding of the active sites involved in heterogeneous catalysis over the broadly used SrTiO 3 (STO), we studied a model reaction, isopropanol conversion, on three differently shape-controlled nanocrystals: cube, truncated cube and dodecahedra. SEM, XRD and XPS confirmed the morphology, phase and composition of STO shapes. Low energy ion scattering (LEIS) revealed the occurrence of surface reconstruction over STO shapes during O 2 pretreatment at different temperatures. Based on the catalytic activities, scanning transmission electron microscopy images and density functional theory calculations, the step sites on STO derived from surface reconstruction were proposed to be the active sites for isopropanol conversion. This was further confirmed by steady state isotopic kinetic analysis (SSITKA) which demonstrated similar intrinsic turnover frequencies (TOFs) for the differently reconstructed STO shapes. It is concluded that the crystal facets impose an indirect effect on the catalysis of STO via controlling the degrees of surface reconstruction: the less stable STO (1 1 0) facet (dodecahedra) leads to more step sites after reconstruction and hence higher overall reaction rate than the more stable (1 0 0) facet (cube). This work highlights the important interplay between the crystal facet and surface reconstruction in controlling the nature and density of active sties and thus catalysis over complex oxides.

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

NETL Electrical Conductivity Relaxation(ECR) Analysis Tool

This application was developed within the solid oxide fuel cell (SOFC) work plan. It calculates the oxygen surface exchange coefficient (k) and oxygen diffusion coefficient (D) for an oxide sample using electrical conductivity relaxation (ECR) data. The tool also produces a heat map showing the associated uncertainty in fitting the ECR data for given pairs of k and D values, so the user can determine the most likely range of k and D values for the tested values. Multiple datasets for the same material can be plotted to decrease uncertainty. The tool can also simulate ECR curves for inputted k and D values and overlay the simulated curves with experimental data for comparison.

electrical conductivity relaxation↗
DOE OSTI2024

Multi-Fidelity Scheme for Accelerating 4D Finite Element Analysis for Mircoreactors

Next-generation microreactors are currently being designed to be operated terrestrial and extraterrestrial for remote surface power production. These systems will provide an alternative source of carbon-free energy that is versatile and can be utilized for various applications. These applications of microreactors have prompted the usage of high-fidelity unstructured finite element (FE) based approaches to provide time-dependent solutions for multiple physics fields. Computing these high-fidelity (4-D) solutions for multiple design iterations, physics fields, and transient events requires an immense computational cost. These high-fidelity solutions are computational expensive and the cost can be reduced through the implementation of less accurate low-fidelity solutions. Unlike the current fleet of commercial nuclear reactors, these next-gen systems present challenges due to the material and physical limitations required. Due to these constraints, the brute force technique of parameterizing important system characteristics determines whether a design meets the project objective. A system such as a nuclear reactor could have thousands of design parameters that affect system performance. In order to analyze the entire parameter space of such a complex system would require millions of CPU hours and countless design iterations. This costly approach is not practical due to regulatory and budget limitations. In this proposal, an approach that utilizes hybrid high-fidelity and low-fidelity FE models to reduce the computational cost of evaluating these applications in 4-D will be presented. The accuracy of the high-fidelity model and the computational efficiency of the low-fidelity model are taken advantage of to produce a solution that closely resembles the full order high-fidelity solution. By utilizing an unconverged coarse FE mesh, operation limits such as temperature, structural loading, etc., can be evaluated in an accelerated fashion and then can be spatially interpolated onto a finer mesh. The resulting error arising from the coarse mesh can be mitigated with a discrepancy function that actively quantifies and corrects the error in the coarse solution. This discrepancy function can be periodically updated with high-fidelity calculations across the temporal domain thus requiring less iterations on the fine mesh. As a result, the computational cost can be reduced for the evaluation and design iteration of microreactors. The proposal for this research contains three sections: Section 2 provides a literature review, Section 3 outlines the methodology for the proposed multi-fidelity scheme, and Section 4 displays preliminary results of the proposed research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Novel Lightweight, Low-Cost Heliostat for Concentrating Solar Power (SBIR Phase I Final Report)

The overall goal of this Department of Energy (DOE) SBIR Phase I project was to demonstrate the feasibility of the L.Garde, Inc. heliostat mirror facet for use in concentrated solar power applications. Prior to this Phase I project, L.Garde, Inc. developed a design for a lightweight and low cost mirror facet and successfully fabricated small scale (0.667 m 2 ) with a mirror surface flatness of 0.7 mrad. Therefore, the focus of the Phase I effort was to (1) demonstrate scalability by fabricating larger size mirrors while retaining the same mirror surface flatness, (2) demonstrate the design can withstand operating environment conditions through field testing, and (3) perform cost analysis to show the heliostat mirror facet design can meet DoE goals of < 12 kg/m2 and <$30/m 2 . During the Phase I project, L.Garde successfully fabricated large area (1.395 m 2 ) heliostat mirror facets, with a mass per unit area of 8 kg/m 2 and a cost of $29.80/m 2 . Surface flatness measurements (using contact probe technique) were performed on the mirror surface and results showed a flatness of 0.68 mrad, which were in agreement with previously fabricated small scale mirror facets. Slope error measurements was also performed on the mirror surface and results showed a slope error of 0.77 mrad which also was in agreement with previously fabricated small scale mirror facets. One large area mirror facet was subjected to a three (3) month field test. Surface flatness and slope error measurements were performed on the mirror surface after the field test and results were 0.663 mrad and 0.58 mrad respectively. The large scale heliostat mirror facets were also subjected to accelerated aging. However, during this exposure we observed failure due to the breaking of the mirror facets. Further investigation into the source of the failure showed that the mirror component alone showed signs of degradation during the same accelerated aging test. We believe that moisture ingress is occurring when subjected to high humidity, leading to buildup of water with the layers of protective paints and resulting in mirror breakage when water expands in freezing conditions. In order to mitigate or eliminate this issue, a new weatherproofing strategy may be implemented or a new mirror component may be used that is able to survive accelerated aging conditions.

14 SOLAR ENERGY↗
DOE OSTI2026

In Situ Analysis of Manganese Antimonate Oxygen Evolution Electrocatalysts via Ambient Pressure X-ray Photoelectron Spectroscopy

Here, manganese antimonates are earth-abundant potential alternatives to precious metal catalysts for the oxygen-evolution reaction (OER). Herein, X-ray photoelectron spectroscopy was used to determine the surface chemistry of a manganese antimonate catalyst for the OER under ultra-high vacuum, ambient pressure, and in situ reaction conditions. Ex situ and in situ analysis revealed the oxidation states of surficial species as a function of material stoichiometry, water content, and applied potential. In situ XPS measurements in 1.0 M KOH(aq) indicated that relative to the rest state, the surface Mn(III) partially oxidized while effecting the OER, with approximately 15% of the Mn signal attributable to Mn(IV) and the remainder attributable to Mn(III).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Protein Surface Printer for Exploring Protein Domains

The surface of proteins is vital in determining protein functions. Herein, a program, Protein Surface Printer (PSP), is built that performs multiple functions in quantifying protein surface domains. Two proteins, PETase and cytochrome P450, are used to validate that the program supports atomistic simulations with different combinations of programs and force fields. A case study is conducted on the structural analysis of the spike proteins of SARS-CoV-2 and SARS-CoV and the human cell receptor ACE2. Although the surface domains of both spike proteins are highly similar, their receptor-binding domains (RBDs) and the O-linked glycan domains are structurally different. The O-linked glycan domain of SARS-CoV-2 is highly positively charged, which may promote binding to negatively charged human cells.

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

Plasma-assisted catalysis for ammonia synthesis in a dielectric barrier discharge reactor: key surface reaction steps and potential causes of low energy yield

Ammonia synthesis experiments were carried out in a coaxial dielectric barrier discharge (DBD) reactor packed with several different supports and metal catalysts. There was a marked increase in the reaction rate, over that obtained in an empty DBD plasma reactor, upon introduction of a packed bed of γ-Al 2 O 3 , Ru/γ-Al 2 O 3 , or SiO 2 particles. The difference in the reaction rates over γ-Al 2 O 3 and Ru/γ-Al 2 O 3 was minimal. Complementary zero-dimensional plasma kinetic model analysis was also performed using inputs from experimental data. This kinetic analysis allowed for gas phase reactions, Eley–Rideal (E–R) reactions, and direct adsorption of radical species on the γ-Al 2 O 3 surface. On the metal surface, dissociative adsorption of N 2 and H 2 , and Langmuir–Hinshelwood reactions were also included. This analysis revealed that, under the conditions of our experiments, ammonia synthesis proceeds principally by the formation of reactive radicals in the gas phase, which then adsorb and participate in E–R reactions on both the metal and support material surfaces. Furthermore, this finding illustrates a challenge for substantially increasing the energy yield for plasma-assisted ammonia synthesis in typical DBD reactors containing packed catalyst beads.

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

Thermal-Hydraulic Analysis of the AFIP-7 Irradiation Test in the Advanced Test Reactor - Oxide Growth Prediction and Correlation

Knowing the thickness of the oxide layer on the surface of aluminum fuel cladding is vitally important for predicting fuel temperature due to the low thermal conductivity of the oxide layer. Several correlation models for predicting oxide growth can be found in the literature. In previous research, the correlations were combined with heat transfer simulations in Abaqus, a finite element analysis code, to forecast the oxide growth. However, this approach requires heat transfer coefficients for modeling heat exchanges with the external flow field, and such coefficients were obtained through empirical equations. Since different empirical equations yield varying heat transfer coefficients, the cladding temperature and predicted oxide thickness both carry a high degree of uncertainty. Here, this research develops a new approach that integrates the fluid flow, fluid and solid heat transfer, and oxide growth correlation(s) into a single computational fluid dynamics model. We demonstrate this approach’s ability to predict oxide development on the AFIP-7 plates during two Advanced Test Reactor (ATR) irradiation cycles. The projected oxide thickness falls within the experimental measurements taken during post irradiation examination.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Super-Resolved Single-Molecule Tracking Studies of Rhodamine B Accumulation on Fresh and Aged Polyethylene Terephthalate

It is well-known that toxic organic micropollutants (OMs) accumulate on the surfaces of microplastics. However, much remains to be learned about the exact molecular level mechanisms of OM accumulation and how these evolve as the plastics age. In this work, super-resolved single-molecule tracking (SMT) is used for the first time to investigate the accumulation of Rhodamine B (RhB) dye on fresh and artificially aged polyethylene terephthalate (PET) surfaces. PET thin films serve as models for microplastics, while RhB serves as a proxy for the OMs they accumulate. Artificial aging of the films is accomplished by exposing them in a UV-ozone chamber. Water contact angle, spectroscopic ellipsometry, and carbonyl index measurements reveal a gradual decrease in film hydrophobicity, thickness, and carbonyl content with age. Atomic force microscopy (AFM) data reveal an increase in surface roughness and confirm that the films remain largely intact and continuous across the aging times explored. In SMT experiments, wide-field fluorescence videos acquired from the water/PET interface under 7.5 pM RhB reveal both mobile and immobile dye molecules. Measurements of the frame-to-frame displacements of the dye show that diffusion occurs by a desorption-mediated mechanism and that the diffusion rate varies with PET film age. The surface density of mobile dye molecules decreases with increasing PET age, while the population of immobile molecules becomes relatively larger, suggesting an age-dependent transformation of the mechanism(s) by which the dye is accumulated. SMT data reveal that both mobile and immobile molecules repeatedly adsorb over the same surface sites, consistent with the emergence of nanoscale PET surface heterogeneity also revealed by AFM. Estimates of the adsorption coefficients are obtained using a nearest-neighbor analysis, giving values from 9.9 × 10 5 to 2.1 × 10 6 M –1 for immobile molecules and from 1.8 × 10 5 to 2.5 × 10 5 M –1 for mobile molecules on fresh and 5 min aged PET, respectively. Here, anomalous age-dependent variations in the velocity of molecular motion on the PET surface and in the population of immobile molecules are shown to correlate with changes in the strength of RhB adsorption.

Adsorption↗
DOE OSTI2022

Random Walks and Sticky Surfaces: Single-Molecule Measurements of Solute Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Nanoporous anodic aluminum oxide (AAO) membranes are now being explored for use in advanced chemical separations, including in the dehydration of biofuels such as ethanol. Optimization of membrane performance requires an in-depth understanding of how solvent mixtures and solutes behave under nanoconfinement. In this work, the diffusion of rhodamine B (RhB) dye through 10 nm and 20 nm AAO nanopores filled with a series of ethanol/water mixtures (0 - 33% water) is explored by fluorescence correlation spectroscopy (FCS). RhB was found to diffuse through the pores by two distinct mechanisms with mean diffusion coefficients, $D_f$ and $D_s$, reflecting fast and slow diffusive motions, respectively, with values that differ by nearly two orders of magnitude. Further, both $D_f$ and $D_s$ increased with pore size and were significantly smaller than $D_b$, the RhB diffusion coefficient in bulk liquid. Mean $D_f$ values follow a composition-dependent trend that closely mimics the viscosity dependence of $D_b$. Additional slowing of fast RhB diffusion is attributed to both hydrodynamic drag and electrostatic interactions with the nanopore surface. The mean $D_s$ values exhibit a different trend with increasing water content, revealing an increase in $D_s$ and a decrease in the contributions of slow diffusion to the observed dynamics. The fluorescence time transient data used in the analysis show that the slow diffusion process is strongly hindered and likely involves frequent adsorption of RhB to the pore surfaces. These results provide new insights into the detailed molecular-level mechanisms of mass transport in nanoporous AAO membranes.

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

Isopotential Electron Titration: Hydrogen Adsorbate-Metal Charge Transfer

The extent of charge transfer between an adsorbate and thermocatalytic surface plays a key role in determining catalytic activity, but direct and quantitative measures have remained elusive. Here, we report the method of isopotential electron titration (IET), an approach that directly measures charge transfer between adsorbates and catalytic surfaces. Charge transfer between Pt and adsorbed hydrogen adatoms was investigated using a catalytic condenser, where the Pt surface was separated from a p-type silicon layer by a hafnia dielectric film. By forcing the Pt and Si layers into isopotential conditions, charge transfer between the adsorbate and Pt surface was titrated through an external circuit. Hydrogen atoms donated electrons to Pt upon adsorption, which was quantitatively reversed upon desorption. Across a temperature range of 125–200 °C (surface hydrogen fractional coverages of 80–100%), the charge transferred to Pt by an adsorbed hydrogen atom was measured to be 0.19 ± 0.01% |e|/H. Bader charge analysis of the extent of charge transfer was in agreement with experimental measurements, with a calculated net donation of 0.4% |e|/H. The ability to experimentally quantify surface charge transfer provides an electronic-based approach to characterize catalytic surfaces, the adsorbed moieties residing on them, and the chemical reactions they accelerate.

Adsorption↗
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

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

Nucleation and growth of molybdenum disulfide grown by thermal atomic layer deposition on metal oxides

To enable greater control over thermal atomic layer deposition (ALD) of molybdenum disulfide (MoS 2 ), here we report studies of the reactions of molybdenum hexafluoride (MoF 6 ) and hydrogen sulfide (H 2 S) with metal oxide substrates from nucleation to few-layer films. In situ quartz crystal microbalance experiments performed at 150, 200, and 250 °C revealed temperature-dependent nucleation behavior of the MoF 6 precursor, which is attributed to variations in surface hydroxyl concentration with temperature. In situ Fourier transform infrared spectroscopy coupled with ex situ x-ray photoelectron spectroscopy (XPS) indicated the presence of molybdenum oxide and molybdenum oxyfluoride species during nucleation. Density functional theory calculations additionally support the formation of these species as well as predicted metal oxide to fluoride conversion. Residual gas analysis revealed reaction by-products, and the combined experimental and computational results provided insights into proposed nucleation surface reactions. With additional ALD cycles, Fourier transform infrared spectroscopy indicated steady film growth after ~13 cycles at 200 °C. XPS revealed that higher deposition temperatures resulted in a higher fraction of MoS 2 within the films. Deposition temperature was found to play an important role in film morphology with amorphous films obtained at 200 °C and below, while layered films with vertical platelets were observed at 250 °C. Finally, these results provide an improved understanding of MoS 2 nucleation, which can guide surface preparation for the deposition of few-layer films and advance MoS 2 toward integration into device manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Numerical study of the effect of propped surface area and fracture conductivity on shale gas production: Application for multi-size proppant pumping schedule design

Hydraulic fracturing is a technique extensively used in the oil and gas industry, where water, proppant (sand) and additives are injected into unconventional reservoirs to enhance the recovery of shale hydrocarbon. Although some previous studies have developed pumping schedules that maximize gas production for a single-size proppant, there are very few studies that consider the effect of varying proppant diameters across pumping stages on shale gas production. Motivated by this, we carried out an extensive sensitivity analysis to determine the effect of different proppant diameters on the average fracture conductivity (FC), average propped surface area (PSA) and cumulative shale gas production volume. We found out that the cumulative shale gas production volume depends on both the average PSA and average FC. We also found out that small-diameter proppant resulted in higher average PSA and lower average FC, whereas large-diameter proppant resulted in lower average PSA and higher average FC. Hence, we designed a multi-size proppant pumping schedule considering both of these parameters into account for simultaneously propagating multiple fractures to maximize shale gas production from unconventional reservoirs. Since the size of injected proppant particles determines the average PSA and average FC for the propped hydraulic fractures, we developed a novel framework called Sequentially Interlinked Modeling Structure (SIMS) to predict the average PSA, average FC and cumulative shale gas production volume at the end of 10 years for a given pumping schedule. Then, we used this SIMS framework to obtain a multi-size proppant pumping schedule that maximizes shale gas production. In closing, we demonstrated that the obtained pumping schedule gives a cumulative shale gas production volume greater than the values obtained from the existing pumping schedules.

03 NATURAL GAS↗