SEARCH · Search NASA

Results for “Equilibrium constant”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Kinetic Pathways to Diblock Copolymer Micelles via Thin Film Dissolution and Cosolvent Variation

We investigate how processing pathways control micelle formation during dissolution of lamellar diblock copolymers using dissipative particle dynamics (DPD) simulations and DPD-selfconsistent field theory (DPD-SCFT). Comparing thin film (TF) and cosolvent-assisted (CS) dissolution reveals fundamentally different micellization mechanisms and final states. In the TF pathway, dissolution proceeds through a metastable cylindrical intermediate that dictates final micelle dimensions. The lamellar melt first transforms into cylindrical micelles, which then fragment into spherical micelles. The final spherical radius is intrinsically linked to the cylinder radius by a constant ratio, Rsphere ≈ 1.4Rcylinder. DPD-SCFT confirms that both cylinders and spheres are near their equilibrium dimensions. In contrast, CS dissolution produces kinetically trapped micelles with sizes significantly below equilibrium. Micellization is hindered by high free energy barriers to single-chain exchange and fusion. Growth occurs primarily via fusion, but increasing corona repulsion suppresses further coarsening, preventing equilibration. These results demonstrate that, when starting from a bulk copolymer phase, the metastable cylindrical intermediate plays a central role in determining final micelle size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Drag exerted by drift-acoustic radiation on a magnetic island

The Island Equilibrium and Transport (IslET) model is extended to allow for the effects of ion compressibility albeit for a constant electron temperature. The model describes fully relaxed, saturated, steady-state magnetic islands arising from the evolution of a tearing mode caused either by an instability or an externally imposed resonant perturbation. It assumes that the collision frequency lies in the semi-collisional regime. A salient new feature of the extended model is an equation describing the transport of parallel momentum. The model is applied to evaluate the drag exerted on a thin island by the emission of drift-acoustic waves.

Physics↗

DOE OSTI2024

CFD unified approach under Eulerian–Lagrangian framework for methanol and gasoline direct injection sprays in evaporative and flash boiling conditions

Innovative synthetic fuels for advanced propulsion systems, such as methanol and ammonia, and synthetic blended fuels (E00, E10, and E30), known for their high volatility, are often injected directly into combustion chambers. It follows that Eulerian–Lagrangian spray models need to accurately capture the spray collapse as a consequence of flash boiling onset and be capable of proficiently handling the preferential evaporation of multi-component fuels in evaporative scenarios. So, we performed the assessment of an Eulerian–Lagrangian CFD code for simulating methanol and E00 gasoline blend sprays in both early and late injection conditions involving flash boiling conditions and preferential evaporation. The adoption of an effervescent breakup model and of a non-equilibrium phase transition model for the discrete phase allows the adoption of a setup that is almost completely free from specific constant tuning, especially for what concerns the breakup model. We validated the simulations using experimental PLV maps of methanol and E00 sprays issued from the ECN Spray M injector. The results highlight a significantly different morphology of the methanol spray compared to the E00 one under late injection conditions. Under stratified combustion, low-volatile fuels are likely to be ignited first, and the flame propagates toward the high-volatile fuels. In conclusion, the spray collapse was also correctly reproduced, inducing the presence of a low-pressure zone and modifying the spray morphology.

E00↗

DOE OSTI2024

Open Circuit Potential of a Au Catalyst during Selective Oxidation of Glycerol

It has been of great interest to understand relationships between the potential of a metal catalyst and its thermochemical catalytic activity, especially in an aqueous phase environment. In the literature, there are correlations of open circuit potential with reaction rates or surface concentrations of reaction intermediates. In this study, we measured the open circuit potential (OCP) of a Au gauze during selective oxidation of glycerol to glyceric acid and hydrogen peroxide in a basic solution, as a function of hydroxide and alcoholate concentrations. It is found that applying a potential to the Au catalyst has no influence on the reaction rate. Although a rough correlation appears to exist between OCP and reaction rate, the data are better fit to an equation which assumes that the potential of the metal (i.e. the OCP) is in equilibrium with the electrochemical potential of the solution, which is defined by the thermodynamic activities of the oxidizing and reducing species. The equation: OCP = constant + Σ ox [(RT/FZ ox )ln[Ox]] - Σ red [(RT/FZ red )ln[Red]] implicitly assumes that the Au metal functions as a probing electrode. It is further found that this equation also applies to the Au-catalyzed H 2 O 2 oxidation/decomposition reaction in a basic medium, and possibly to formic acid dehydrogenation. Here, we postulate that the apparent correlation between OCP and reaction rate is due to the fact that the reaction rate is proportional to the concentrations of the reaction products and/or reactants that define the electrochemical potential of the solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Degrees of Rate Control in Interconnected Reaction Networks

Overall reactions in interconnected networks exhibit net, forward, and reverse rates that are governed by both constitutive elementary steps in the pathway of interest and branching elementary steps that lead to alternative products. Accordingly, steps in branching pathways exhibit negative net, forward, and reverse degrees of rate control, as they reduce reaction flux to the desired product. We here contextualize the forward and reverse degrees of rate control in terms of kinetic resistances (inverse of rates) and leverage the additive nature of kinetic resistance to decouple kinetic driving forces contributed by constitutive elementary steps and branching points (nodal species) in interconnected networks. Regardless of the network connectivity, forward and reverse degrees of rate control are shown to converge at equilibrium. Away from equilibrium, we identify two critical features of interconnected networks: stoichiometric regularity─condition where all stoichiometric numbers are unity─and pathway symmetry around nodal species─condition where branching pathways share the same rate constants, stoichiometry, and species concentrations/activities─that result in (i) equal forward, reverse, and consequently net degrees of rate control and (ii) forward and reverse degrees of rate control that exhibit constant offsets, respectively, across all extents of reaction. Furthermore, our discourse further provides a mathematical description for the influence of stoichiometric irregularity and pathway asymmetry on forward and reverse degrees of rate control. Altogether, the presented work details the effects of network (inter)connectivity and stoichiometry on reaction kinetics and, in doing so, establishes general protocols for capturing these effects as additive terms in the formulation of forward and reverse degrees of rate control.

10 SYNTHETIC FUELS↗

DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Toward Rapid Actinium-225 Purification via Membrane Adsorbers with Covalently Tethered Diglycolamide Ligands

Extractive diglycolamide (DGA) resins are used in several state-of-the-art techniques for purifying 225 Ac, a promising radiometal for targeted alpha therapy. Unfortunately, separation processes that rely on resins are often limited to slow flow rates, high elution volumes, and long processing times. Membrane adsorbers functionalized with DGA ligands are an alternative separation material that may overcome these challenges. This work presents (1) the synthesis of an aminated tetrahexyldiglycolamide ligand, (2) the covalent tethering of the ligand to electrospun poly(vinylbenzyl chloride) fiber mats, and (3) the adsorption and desorption of La(III) and 225 Ac. Chemical and physical characterization supports the covalent tethering of the ligand to the fiber mat, as well as the preservation of the fiber surface area and porosity after functionalization. Equilibrium adsorption experiments were performed with stable La(III) and radioactive 225 Ac. Trends in affinity are consistent between commercial resins and the synthesized membrane adsorbers; however, the Langmuir constants and the maximum binding capacity of the membrane adsorbers were generally lower than the resins. Despite these differences, the modeled selectivity for an equimolar solution of La(III)/ 225 Ac in 10 M nitric acid is 57. Furthermore, 225 Ac is rapidly desorbed from the fibers in 10 M nitric acid (<20 min). The La(III)/ 225 Ac selectivity and rapid 225 Ac desorption indicate this class of materials is promising for rapid radioanalytical separations.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗

DOE OSTI2025

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Visualizing the Active Site Oxyanion Loop Transition Upon Ensitrelvir Binding and Transient Dimerization of SARS-CoV-2 Main Protease

N-terminal autoprocessing from its polyprotein precursor enables creating the mature-like stable dimer interface of SARS-CoV-2 main protease (MPro), concomitant with the active site oxyanion loop equilibrium transitioning to the active conformation (E*) and onset of catalytic activity. Here, through mutagenesis of critical interface residues and evaluating noncovalent inhibitor (ensitrelvir, ESV) facilitated dimerization through its binding to MPro, we demonstrate that residues extending from Ser1 through Glu14 are critical for dimerization. Combined mutations G11A, E290A and R298A (MPro™) restrict dimerization even upon binding of ESV to monomeric MPro™ with an inhibitor dissociation constant of 7.4 ± 1.6 µM. Contrasting the covalent inhibitor NMV or GC373 binding to monomeric MPro, ESV binding enabled capturing the transition of the oxyanion loop conformations in the absence of a reactive warhead and independent of dimerization. Characterization of complexes by room-temperature X-ray crystallography reveals ESV bound to the E* state of monomeric MPro as well as an intermediate approaching the inactive state (E). It appears that the E* to E equilibrium shift occurs initially from G138-F140 residues, leading to the unwinding of the loop and formation of the 3 10 -helix. Finally, we describe a transient dimer structure of the MPro precursor held together through interactions of residues A5-G11 with distinct states of the active sites, E and E*, likely representing an intermediate in the autoprocessing pathway.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Selective electron-phonon coupling strength from nonequilibrium optical spectroscopy: The case of MgB 2

The coupling between quasiparticles and bosonic excitations rules the energy transfer pathways in condensed matter systems. The possibility of inferring the strength of specific coupling channels from their characteristic timescales measured in nonequilibrium experiments is still an open question. Here, in this study, we investigate MgB 2 , in which conventional superconductivity at temperatures as high as 39 K is mediated by the strong coupling between the conduction electrons and the 𝐸 2⁢𝑔 phonon mode. By means of broadband time-resolved optical spectroscopy, we show that this selective electron-phonon coupling dictates the nonequilibrium optical response of MgB 2 at early times (<100 fs) after photoexcitation. Furthermore, based on an effective temperature model analysis, we estimate its contribution to the total electron-boson coupling function extracted from complementary equilibrium spectroscopy approaches, namely, optical reflectivity and angle-resolved photoemission spectroscopy. The coupling strength with the 𝐸 2⁢𝑔 phonon modes is thus estimated to be 𝜆 ≃ 0.56, which is approximately half of the total coupling constant, in agreement with 𝑎𝑏 initio calculations from the literature. As a benchmark, broadband time-resolved optical spectroscopy is performed also on the isostructural and nonsuperconducting compound AlB 2 , showing that the nonequilibrium optical response relaxes on a slower timescale due to the lack of strongly coupled phonon modes. Our findings demonstrate the possibility to resolve and quantify selective electron-phonon coupling from nonequilibrium optical spectroscopy.

Mor, Selene [Università Cattolica del Sacro Cuore,↗

DOE OSTI2026

Improved Accuracy in Semi-Experimental Structure Determination by Resolving Problems Associated with Rotation of Principal Inertial Axes of Isotopologues: Structures of 1,3-Oxazole ( c -C 3 H 3 NO)

The rotational spectrum of the normal isotopologue of 1,3-oxazole (c-C 3 H 3 NO) was observed from 43 to 750 GHz. Over 3900 transitions for the ground vibrational state are measured, assigned, and least-squares fit to sextic centrifugally distorted-rotor Hamiltonians. The measured frequencies and resulting spectroscopic constants from this extended spectral range, combined with previous measurements of the nuclear quadrupole coupling constants, will facilitate astronomical searches for oxazole across the majority of the range of modern radiotelescopes. Spectra for a set of 30 oxazole isotopologues, which include multiple isotopic substitutions of each atom, are used to determine the first semi-experimental equilibrium ($r$$^{SE}_{e}$) structure and semi-experimental substitution structure ($r$$^{SE}_{e}$), each using CCSD(T) computed values for the vibration–rotation interaction and electron-mass corrections. The large number of isotopologues, including 21 isotopologues observed for the first time, and the redundant substitutions of each atom provide sufficient spectroscopic information to determine the $r$$^{SE}_{e}$ structure with the expected high level of accuracy and precision (0.0001 or 0.0002 Å in bond distances and 0.013 to 0.025° in bond angles). In the course of this study, we analyzed a known issue for some $r$$^{SE}_{e}$ structure determinations of near-oblate asymmetric tops in which inclusion of individual isotopologues degrades the structure determination. We demonstrate that this problem primarily arises from the difference in the values of the computed vibration–rotation interaction corrections as evaluated at the computed re geometry vs the $r$$^{SE}_{e}$ geometry of the “real” molecule. Our solution to this problem substantially improves the $r$$^{SE}_{e}$ structure of oxazole and likely can be generalized to many other molecules.

Chemical structure↗

DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Temperature-Dependent Sn Incorporation and Defect Formation in Pseudomorphic SiSn Layers on Si (001) via Molecular Beam Epitaxy

SiSn alloys have attracted growing interest for group-IV bandgap engineering, although their epitaxial growth remains challenging due to the extremely low equilibrium solubility of Sn in Si. In this work, fully strained (pseudomorphic) SiSn epitaxial layers were grown on Si (001) substrates by means of molecular beam epitaxy. A systematic investigation reveals a strong inverse correlation between growth temperature and Sn incorporation efficiency. Despite a constant Sn flux, the incorporated Sn composition decreases from 5.5% to 3.2% as the growth temperature increases, indicating a pronounced temperature dependence of Sn incorporation. Reflection high-energy electron diffraction indicates a gradual transition of the growth from two-dimensional to three-dimensional with increasing film thickness. Structural characterization by means of X-ray diffraction, atomic force microscopy, and transmission electron microscopy confirms the pseudomorphic growth and smooth surface morphology and reveals twins and stacking faults near the surface region. These results establish a quantitative reference for SiSn growth kinetics and provide guidance for future studies of SiSn and SiGeSn alloys in silicon-compatible electronic and optoelectronic applications.

SiSn alloys↗

DOE OSTI2025

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗

DOE OSTI2025

Non-equilibrium dynamics in geometrically frustrated spin glass Bi 2 Fe 3 GaO 9 with a Cairo lattice

We have explored the magnetic relaxation of the spin glass phase in Bi 2 Fe 3 GaO 9 , a geometrically frustrated magnet with a unique Cairo lattice made of pentagonal building blocks. Using dc magnetization measurements, we find the relaxation behavior in single crystals follows the Arrhenius law as a function of temperature but remains relatively constant under different magnetic fields. Through a carefully designed protocol, we observe significant rejuvenation and memory effects, reminiscent of classical dilute spin glasses rather than spin jamming systems. Our finding suggests that Bi 2 Fe 3 GaO 9 is probably located in the crossover between conventional spin glass and perfect spin jam and thus provides an excellent platform for investigating the evolution of underlying magnetic structure, paving the way for further neutron scattering research.

aging↗

DOE OSTI2025

Revolutionizing Methane Transformation with the Dual Production of Aromatics and Electricity in a Protonic Ceramic Electrocatalytic Membrane Reactor

Reducing the energy and carbon intensity of the conventional chemical processing industry can be achieved by electrochemically transforming natural gases into higher-value chemicals with higher efficiency and near-zero emissions. In this work, the direct conversion of methane to aromatics and electricity has been achieved in a protonic ceramic electrocatalytic membrane reactor through the integration of a proton-conducting membrane assembly and a trimetallic Pt–Cu/Mo/ZSM-5 catalyst for the nonoxidative methane dehydro-aromatization reaction. In this integrated system, a remarkable 15.6% single-pass methane conversion with an 11.4% benzene yield has been demonstrated, while a peak power density of 276 mW cm –2 is obtained at 700 °C. The enhanced 15.7% increase in conversion and 16.0% improvement in the yield are observed when compared with the thermochemical process, which is attributed to the shift of reaction equilibrium by the removal of hydrogen through the protonic membrane. Concurrently, the faster H2 removal at a higher electrical current gave rise to a higher methane conversion and benzene yield. Furthermore, the catalyst can be efficiently regenerated by eliminating carbon deposition. A stable cell potential is maintained for 45 h under a constant current load of 0.13 A cm –2 . Lastly, the dual production of aromatics and electricity in the electrocatalytic membrane reactor has been demonstrated to be an attractive approach for decarbonizing chemical processing.

aromatic compounds↗