Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular functionalization”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Wave function analysis with a maximum flow algorithm

An efficient algorithm for computing the maximum-flow path in a network is applied to the identification of the dominant configuration state functions (CSFs) in a graphically contracted function (GCF), configuration interaction, wave function. The flow network is a space of spin-adapted CSFs represented by a Shavitt graph, wherein the nodes correspond to orbital occupations and spin quantum numbers. The graph nodes are connected by arcs, and an arc density is defined as sums of the associated squared CSF coefficients. A max-min approach determines an upper bound to the maximum possible incoming flow for each graph node. A backtracking step generates a candidate walk and is followed by a limited search of alternative branching paths for the dominant CSF. The arc density contributions are removed from the graph, and the algorithm is reapplied to the updated graph. This list of generated walks can be partitioned in order to guarantee that the dominant CSFs have been identified. All of the steps in this algorithm are computationally efficient and do not depend on the potentially large dimension of the underlying linear CSF expansion space. An analysis of low-lying valence states of C-2 illustrates the method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Ab Initio Dynamics of Graphene and Graphyne Electrodes in Vacuum and in the Presence of Electrolytes

Graphene and graphyne electrodes are both relevant for electrochemical energy storage applications due to their unique physical, chemical, and electronic properties. Graphyne has been considered even more attractive than graphene due to its larger specific surface area, greater electronic mobility, and intrinsic band gap. The analyses presented in this work reveal relevant features of graphyne in an interaction with an electrolyte that can be useful for practical applications. Although there is already a considerable volume of work on graphyne, the analysis of its properties taking into account the dynamics of its network and the interaction with the condensed medium has not yet been presented. Here, ab initio molecular dynamics simulations were performed to study the energetic, structural, spectroscopic, and electronic properties of graphene and graphyne electrodes in vacuum and in the presence of an electrolyte (EMIM-BF 4 ). A detailed analysis of the electrode–electrolyte interaction is reported for both systems. Structural features are presented in terms of both radial and spatial distribution functions and angular orientation distributions. Power spectral analyses and molecular orbital maps provided important insights into the intermolecular interactions in graphene and graphyne systems, describing the main differences as well as revealing similarities between them.

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

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

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

Semi-adsorption-controlled growth window for half-Heusler FeVSb epitaxial films

The electronic, magnetic, thermoelectric, and topological properties of Heusler compounds (composition XYZ or X 2 YZ) are highly sensitive to stoichiometry and defects. Here we establish the existence and experimentally map the bounds of a semi-adsorption-controlled growth window for semiconducting half-Heusler FeVSb films, grown by molecular beam epitaxy (MBE). We show that due to the high volatility of Sb, the Sb stoichiometry is self-limiting for a finite range of growth temperatures and Sb fluxes, similar to the growth of III-V semiconductors such as GaSb and GaAs. Films grown within this window are nearly structurally indistinguishable by x-ray diffraction (XRD) and reflection high energy electron diffraction (RHEED). Here, the highest electron mobility and lowest background carrier density are obtained towards the Sb-rich bound of the window, suggesting that Sb vacancies may be a common defect. Similar semi-adsorption-controlled bounds are expected for other ternary intermetallics that contain a volatile species Z = {Sb, As, Bi}, e.g., CoTiSb, LuPtSb, GdPtBi, and NiMnSb. However, outstanding challenges remain in controlling the remaining Fe/V (X/Y) transition metal stoichiometry.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Literary review and modeling outlook for radiation damage in fused-silica optical fibers

The purpose of this document is to distill, from literature, the effects of radiation damage in fused-silica optical fibers. This information will be used to design a theoretical framework capable of predictive modeling that will be published elsewhere. By inspection of the relevant literature, the devised method is calculation of defect amorphous structures by molecular dynamic techniques, and computation of the macroscopic dielectric function using density functional theory. The dielectric function is then used to directly compute the transmittance of the core material before and after irradiation via solution of Maxwell’s equations in a homogeneous medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2020

Viscous flow properties and hydrodynamic diameter of phenothiazine-based redox-active molecules in different supporting salt environments

We report viscous flow properties of a redox-active organic molecule, N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT), a candidate for non-aqueous redox flow batteries, and two of its radical cation salts. A microfluidic viscometer enabled the use of small sample volumes in determining viscosity as a function of shear rate and concentration in the non-aqueous solvent, acetonitrile, both with and without supporting salts. All solutions tested show Newtonian behavior over shear rates of up to 30 000 s−1, which was rationalized by scaling arguments for the diffusion-based relaxation time of a single MEEPT molecule without aggregation. Neat MEEPT is flowable but with a large viscosity (412 mPa⋅s at room temperature), which is ∼1000 times larger than that of acetonitrile. MEEPT solutions in acetonitrile have low viscosities; at concentrations up to 0.5 M, the viscosity increases by less than a factor of two. From concentration-dependent viscosity measurements, molecular information was inferred from intrinsic viscosity (hydrodynamic diameter) and the Huggins coefficient (interactions). Model fit credibility was assessed using the Bayesian Information Criterion. It is found that the MEEPT and its charged cations are “flowable” and do not flocculate at concentrations up to 0.5 M. MEEPT has a hydrodynamic diameter of around 8.5 Å, which is almost insensitive to supporting salt and state of charge. This size is comparable to molecular dimensions of single molecules obtained from optimized structures using density functional theory calculations. The results suggest that MEEPT is a promising candidate for redox flow batteries in terms of its viscous flow properties.

Wang, Yilin (ORCID:0000000316102995)↗
DOE OSTI2023

Direct Synthesis of Ammonia from Nitrate on Amorphous Graphene with Near 100% Efficiency

Abstract Ammonia is an indispensable commodity in the agricultural and pharmaceutical industries. Direct nitrate‐to‐ammonia electroreduction is a decentralized route yet challenged by competing side reactions. Most catalysts are metal‐based, and metal‐free catalysts with high nitrate‐to‐ammonia conversion activity are rarely reported. Herein, it is shown that amorphous graphene synthesized by laser induction and comprising strained and disordered pentagons, hexagons, and heptagons can electrocatalyze the eight‐electron reduction of NO 3 − to NH 3 with a Faradaic efficiency of ≈100% and an ammonia production rate of 2859 µg cm −2 h −1 at −0.93 V versus reversible hydrogen electrode. X‐ray pair‐distribution function analysis and electron microscopy reveal the unique molecular features of amorphous graphene that facilitate NO 3 − reduction. In situ Fourier transform infrared spectroscopy and theoretical calculations establish the critical role of these features in stabilizing the reaction intermediates via structural relaxation. The enhanced catalytic activity enables the implementation of flow electrolysis for the on‐demand synthesis and release of ammonia with >70% selectivity, resulting in significantly increased yields and survival rates when applied to plant cultivation. The results of this study show significant promise for remediating nitrate‐polluted water and completing the NO x cycle.

Huang, Libei↗
DOE OSTI2020

Selectively Depolymerizable Polyurethanes from Unsaturated Polyols Cleavable by Olefin Metathesis

This communication describes a novel series of linear and crosslinked polyurethanes (PUs) and their selective depolymerization under mild conditions. Two unique polyols are synthesized bearing unsaturated units in a configuration designed to favor ring-closing metathesis (RCM) to five- and six-membered cycloalkenes. These polyols are co-polymerized with toluene diisocyanate to generate linear PUs and trifunctional hexamethylene- and diphenylmethane-based isocyanates to generate crosslinked PUs. The polyol design is such that the RCM reaction cleaves the backbone of the polymer chain. Upon exposure to dilute solutions of Grubbs’ catalyst under ambient conditions, the PUs are rapidly depolymerized to low molecular weight, soluble products bearing vinyl and cycloalkene functionalities. These functionalities enable further re-polymerization by traditional strategies for polymerization of double bonds. We anticipate that this general approach can be expanded to develop a range of chemically recyclable condensation polymers that are readily depolymerized by orthogonal metathesis chemistry.

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

Reviewing computational studies of defect formation and behaviors in carbon fiber structural units

Solid-state carbon materials, such as graphite and graphene, are at the forefront of materials research because of their unique electronic, vibrational, and mechanical properties, leading to a broad range of potential and realized applications. One key application is their role as basic structural units of carbon fiber (CF), a lightweight alternative to steel. In CF, a delicate relationship exists between ultimate material strength and the atomic-scale density of defects contained at the inter- and intra-subunit levels. Computational studies provide insight into the stability of various types of defects that can form in these systems and connect with experimental observables such as bandgap and spectroscopic measurements. Therefore, the literature contains many computational studies that focus on changes induced by defects, including vacancies and Dienes transformations (Stone-Wales and Thrower defects). However, wide-ranging methods and cell sizes have been used, and property-specific information is often lacking. Furthermore, this review summarizes key literature findings, paying particular attention to changes in the electronic, vibrational, and mechanical properties induced by defects in graphitic materials relevant to CF.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shaping Nanobodies and Intrabodies against Proteoforms

Proteoforms expand genomic diversity and direct developmental processes. While high-resolution mass spectrometry has accelerated characterization of proteoforms, molecular techniques working to bind, and disrupt the function of specific proteoforms have lagged behind. In this study, we worked to develop intrabodies capable of binding specific proteoforms. Herein we employed a synthetic camelid nanobody library expressed in yeast to identify nanobody binders of different SARS-CoV-2 receptor binding domain (RBD) proteoforms. Importantly, employment of the positive and negative selection mechanism inherent to the synthetic system allowed for amplification of nanobody-expressing yeast that bind to the original (Wuhan strain RBD) but not the E484K (Beta variant) mutation. Nanobodies raised against specific RBD proteoforms were validated by yeast-2-hybrid analysis and sequence comparisons. These results provide a clear framework for development of proteoform-specific nanobodies and intrabodies.

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

Switching of Hybrid Improper Ferroelectricity in Oxide Double Perovskites

In ABO 3 -type perovskite oxides with Pnma symmetry, rotation (Q R+ , a 0 a 0 c + ) and tilt (Q T , a – a – c 0 ) of BO 6 octahedra are the two primary order parameters. These order parameters establish an inherent trilinear coupling with anti-ferroelectric A-site displacement (Q AFE ) to form the low-symmetry phase. The symmetry is further lowered in double perovskite oxides (DPOs) due to A/A' cation ordering. It in turn makes these systems polar via hybrid improper ferroelectric mechanism, primarily driven by Q R+ and Q T . Naturally, it has been believed that functionalities such as polarization can also be switched by tuning these primary order parameters. However, mystery around finding switching mechanism still remains. Our study based on density functional theory calculations combined with finite-temperature molecular dynamics simulations shows that the polarization switching is a two-step process, driven by out-of-phase rotation (Q R– , a 0 a 0 c – when Q T = 0 or, a – a – b – when Q T ≠ 0). A series of polar DPOs such as KLnFeOsO 6 [Ln = Sm, Gd, Dy, Tm (lanthanides) and Y (rare earth)], all belonging to P2 1 symmetry, are considered in this investigation. The polarization switching P ($\overrightarrow{P}$) occurs at a very high temperature of ~1150 K through a phase transition, from a polar (P2 1 ) phase with $\overrightarrow{P}$(+) to $\overrightarrow{P}$(-) via a non-polar P4/n phase. The switching itself is metastable in nature. The switching (both polarization and spin state) is only observed for a very short period of time (~23 ps) that poses limitation on using such a mechanism in memory device realization. We demonstrate a concurrent heating–cooling procedure to overcome such shortcoming. In conclusion, simulations conducted at 600 K further imply that long lasting switching can be achieved, at least for 1.2 ns for 600 K, and ideally for an infinite time, if the material is heated just above the T c followed by rapid cooling to a temperature below T c .

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

Y(III) Sorption at the Orthoclase (001) Surface Measured by X-ray Reflectivity

Interactions of heavy metals with charged mineral surfaces control their mobility in the environment. Here, we investigate the adsorption of Y(III) onto the orthoclase (001) basal plane, the former as a representative of rare earth elements and an analogue of trivalent actinides and the latter as a representative of naturally abundant K-feldspar minerals. In this study, we apply in situ high-resolution X-ray reflectivity to determine the sorption capacity and molecular distribution of adsorbed Y species as a function of the Y 3+ concentration, [Y 3+ ], at pH 7 and 5. With [Y 3+ ] ≥ 1 mM at pH 7, we observe an inner-sphere (IS) sorption complex at a distance of ~1.5 Å from the surface and an outer-sphere (OS) complex at 3–4 Å. Based on the adsorption height of the IS complex, a bidentate, binuclear binding mode, in which Y 3+ binds to two terminal oxygens, is proposed. In contrast, mostly OS sorption is observed at pH 5. The observed maximum Y coverage is ~1.3 Y 3+ /A UC (A UC : area of the unit cell = 111.4 Å 2 ) for all the investigated pH values and Y concentrations, which is in the expected range based on the estimated surface charge of orthoclase (001).

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

Solvation Structure of 237 Np 4+ in a Noncomplexing Environment

Here, the solvation structure of an Np 4+ ion in an aqueous, noncomplexing and nonoxidizing environment of trifluoromethanesulfonic (triflic) acid was investigated with X-ray absorption spectroscopy (XAS) combined with ab initio molecular dynamics (AIMD) and time-dependent density functional theory (TDDFT) calculations. Np L III -edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data were collected for Np 4+ in 1, 3, and 7 M triflic acid using a laboratory-scale spectrometer and separately at a synchrotron facility, producing data sets in excellent agreement. TDDFT calculations revealed a weak pre-edge feature not previously reported for Np L III -edge XANES. AIMD modeling results showed differences in the hydration shell of the Np 4+ ion at different concentrations of triflic acid; these results are supported by the experiment. EXAFS fit models to the experiment resulted in similar coordination of Np 4+ in noncomplexing aqueous media as reported in the literature for 1 M perchloric acid but, together with calculations, revealed more than one distance between Np and O atoms in 7 M triflic acid. These results imply monodentate coordination with sulfonate groups in 7 M triflic acid and suggest the possibility of proto-neptunyl species in relatively low-concentration Np 4+ acid solutions.

Boglaienko, Daria V. [Pacific Northwest National L↗
DOE OSTI2020

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗
DOE OSTI2021

Methane Adsorption on Heteroatom-Modified Maquettes of Porous Carbon Surfaces

Experimental and theoretical studies disagree on the energetics of methane adsorption on carbon materials. However, this information is critical for the rational design and optimization of the structure and composition of adsorbents for natural gas storage. The delicate nature of dispersion interactions, polarization of both the adsorbent and the adsorbate, interplay between H-bonding and tetrel bonding, and induced dipole/Coulomb interactions inherent to methane physisorption require computational treatment at the highest possible level of theory. In this study, we employed the smallest reasonable computational model, a maquette of porous carbon surfaces with a central site for substitution and methane binding. The most accurate predictions of methane adsorption energetics were achieved by electron-correlated molecular orbital theory CCSD(T) and hybrid density functional theory MN15 calculations employing a saturated, all-electron basis set. Here, the characteristic geometry of methane adsorption on a carbon surface (“lander approach”) arises due to bonding interactions of the adsorbent π-system with the proximal H–C bonds of methane, in addition to tetrel bonding between the antibonding orbital of the distal C–H bond and the central atom of the maquette (C, B, or N). The polarization of the electron density, structural deformations, and the comprehensive energetic analysis clearly indicate a ∼3 kJ mol –1 preference for methane binding on the N-substituted maquette. The B-substituted maquette showed a comparable or lower binding energy than the unsubstituted, pure C model, depending on the level of theory employed. The calculated thermodynamic results indicate a strategy for incorporating electron-enriched substitutions (e.g., N) into carbon materials as a way to increase methane storage capacity over electron-deficient (e.g., B) modifications. The thermochemical analysis was revised for establishing a conceptual agreement between the experimental isosteric heat of adsorption and the binding enthalpies from statistical thermodynamics principles.

03 NATURAL GAS↗