Search NASA⌕ Search

SEARCH · Search NASA

Results for “Facets”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Graphically Contracted Function Construction with the Recursive Pairwise Merge Algorithm

A new implementation of a recursive pairwise merge algorithm to construct a GCF from a list of CSF expansion coefficients is presented. The essential new feature is the preallocation of some work arrays used within the intermediate steps of the merge procedure. This results in roughly an order of magnitude improvement in overall efficiency and also approximately eliminates a factor of n, the molecular orbital dimension, from the original implementation. Initial application of this merge procedure to a series of H m molecules shows that the GCF wave functions can be represented well both with delocalized canonical Hartree-Fock orbitals and with localized molecular orbitals. Finally, for a given wave function complexity, as measured by the average facet count, $\bar{\text {f}}$, the delocalized Hartree-Fock orbitals show smaller errors for small $\bar{\text {f}}$ values, while the localized orbitals show smaller errors for larger $\bar{\text {f}}$ values.

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

Self-Assembly of a Linear Alkylamine Bilayer around a Cu Nanocrystal: Molecular Dynamics

Copper nanocrystals are often grown with the help of alkylamine capping agents, which direct the nanocrystal shape. However, the role of these molecules is still unclear. We characterized the assembly of aqueous tetradecylamine (TDA) around a Cu nanocrystal and found that TDA exhibits a temperature-dependent bilayer structure. The bilayer involves an inner layer, in which TDA binds to Cu via the amine group and tends to orient the alkyl tail perpendicular to the surface, and an outer layer whose structure depends on temperature. At low temperatures, alkylamines in the inner layer form bundles with no apparent relation to the crystal facets. Alkylamines in the outer layer tend to orient their long axes perpendicular to the Cu surfaces, with interdigitation into the inner layer. At high temperatures, alkylamines in the inner layer lose their bundle structure, and outer-layer alkylamines tend to orient themselves tangential to the Cu surfaces, forming a “web” above inner-layer TDA. TDA exhibits a rapid interlayer exchange at typical synthesis temperatures, consistent with experiment. Here, the variety in the assemblies seen here and in other studies of alkanethiols around gold nanocrystals indicates a richness in the assemblies that can be achieved by modulating the interaction between the strongly binding end group and the surface.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Surface Hydration and Hydroxyl Configurations of Gibbsite and Boehmite Nanoplates

Discontinuation of the bulk structure at the interface between metal oxide particles and water leads to altered bonding characteristics and unique facet-dependent molecular environments. Surface hydration and consequent hydroxylation adds further complexity to the interface, details that for metal (oxy)hydroxides are especially difficult to isolate from the background signal of bulk structural hydroxyls. Here we probe for the first time the hydroxylation and hydration structures on basal surfaces of gibbsite (?-Al(OH)3) and boehmite (?-AlOOH) nanoplates at ambient conditions, using the interface-sensitive technique vibrational sum frequency generation spectroscopy (VSFG). Hydroxyl stretching modes at the interfaces with adsorbed water layers compared directly to Raman and infrared bulk modes show that while gibbsite surface frequencies were sharp and nearly identical to those in the bulk, boehmite surface hydroxyls displayed a very different broad spectrum of states. Ab initio molecular dynamics simulations of both basal surfaces with and without hydration waters reveal that gibbsite surface hydroxyls interact only weakly with overlying hydration waters remaining essentially unperturbed, whereasthose on boehmite can interact more strongly facilitated by higher configurational degrees of freedom at the interface. The findings clearly unveil substantial differences in the hydrated interfacial dynamics of these two otherwise similar materials, with implications for their interfacial chemistry and wettability.

Sassi, Michel JPC↗
DOE OSTI2020

Density Functional Theory Investigation of Oxidation Intermediates on Gold and Gold–Silver Surfaces

Gold and gold-silver alloys can be active and selective oxidation catalysts. Previous work has suggested that O 2 dissociation occurs at bimetallic step sites on gold-silver alloys, but the site responsible for the rest of the reaction steps has not been studied. As a first step in gaining insight into this issue, in this work we investigated the adsorption of oxygen and other oxidation intermediates on the (111) and (211) facets of gold-silver alloys using density functional theory. Oxygen and silver coverage effects were analyzed, and different model structures were compared. We also examined the energy barriers for the diffusion of atomic oxygen to gain insight into O migration and spillover. On (111) surfaces, O adsorption is much stronger at low O coverage (less than 0.22 ML), while on (211) surfaces O is strongly bound at both high and low O coverage. O diffusion across the step is faster than diffusion along the step. Ag stabilizes O, both when directly bound to it and when in an adjacent site. Ag also reduces repulsive O-O interactions at low O coverage. Our calculated reaction barriers for O-assisted CH 3 O dehydrogenation suggest that reaction is faster on steps than on terraces. Overall, our findings suggest that spillover of O from Ag-rich steps to Au-rich terraces does not occur, and that oxidation reactions on gold-silver alloys occurs on step sites. More specifically, oxidation likely occurs either on Ag-rich step sites or on Au-rich step sites that are adjacent to Ag-rich step sites.

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

Stable Surface Terminations of a Perovskite Oxyhydride from First-Principles

Successful synthesis of some perovskite oxyhydrides and their unique catalytic properties have recently attracted researchers’ attention. However, their surface structure remains unclear. Here we identify stable surface terminations of a prototypical perovskite oxyhydride, BaTiO 2.5 H 0.5 , under catalytically relevant temperatures and pressures by using first-principles thermodynamics based on density functional theory. Five low-index facets, including (100), (010), (210), (011), and (211), and their various terminations for a total of 47 different surfaces have been examined for relative stability at different temperatures (700, 500, 300 K) and gas environments (10 –15 ≤ P O2 ≤ 1 atm, 10–15 ≤ PH2 ≤ 100 atm). The most stable ones are found to be (010)-Ba 2 O 2 , (210)-Ti 2 O 2 , and (211)-Ba 2 O 4 H surface terminations. These polar surfaces are stabilized by charge compensation. This work provides important insights into the stable surfaces of perovskite oxyhydrides for future studies of their catalytic properties.

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

Conversion of Formic Acid on Single- and Nano-Crystalline Anatase TiO 2 (101)

Understanding thermochemical transformations of formic acid (FA) on metal oxide surfaces is important for many catalytical reactions. Here we study thermally induced reactions of FA on a single-crystalline and nanocrystalline anatase TiO 2 (101). We employ a combination of scanning tunneling microscopy (STM), temperature-programmed desorption (TPD), infrared reflection absorption spectroscopy (IRAS), diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and density functional theory (DFT) to follow the FA surface intermediates and reaction products above room temperature. We find that the primary reaction products desorbing at about 300, 480, and 515 K are molecular water, carbon monoxide, and formaldehyde, respectively. Bidentate (BD) formate and bridging hydroxyl (HO b ) are identified as central intermediates in the FA transformations. Bridging oxygen vacancies (V O ) are also likely participants despite their low stability at the surface. In conclusion, the parallel studies on single crystals and faceted TiO 2 (101) nanoparticles reveal the spectroscopic commonalities of surface species and of the thermal conversion of molecular and deprotonated forms of FA.

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

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mechanistic Study of the Li–Air Battery with a Co 3 O 4 Cathode and Dimethyl Sulfoxide Electrolyte

The lithium–air battery, a powerful competitor to replace the traditional lithium-ion battery, has attracted increasing attention due to its extremely high theoretical energy density. However, its development is limited by the cathode and electrolyte properties, which should include high stability, conductivity, and electrocatalytic properties in oxygen-rich environments. Here, we employ a systematic first-principles study of Li–O 2 discharge and charge reactions on the Co 3 O 4 -based cathode with the assistance of dimethyl sulfoxide (DMSO) electrolyte. The structure, stability, and electronic properties of different surface reconstructions of the Co 3 O 4 (100) facet are investigated. In addition, the mechanisms and thermodynamic overpotentials of multi-step reactions between Li + /e – and O 2 are provided, where lithium suboxide products (Li 2 O 2 or Li 3 O 2 ) are formed on the different Co 3 O 4 (100) terminations. The solvation shell of Li + components in explicit DMSO solvent is investigated through ab initio molecular dynamics simulations. In general, we find that the Co 3 O 4 (100)-O (oxidized) surface is the most stable one under standard conditions, and the stable Li + solvation structure is found in a tetrahedral Li(DMSO) 4 + shell in the DMSO-based electrolyte. Moreover, in the system of the Co 3 O 4 (100)-O cathode and DMSO electrolyte, the solution model pathway is energetically favorable for the Li–O 2 discharge reaction. It provides a low constant overpotential of 0.17 V during a long-term discharging process, thus causing the final toroid Li 2 O 2 formation on the cathode. During the charging process, an overpotential of 0.36 V is required to rapidly decompose Li 2 O 2 .

25 ENERGY STORAGE↗
DOE OSTI2022

Trends in Formic Acid Electro-Oxidation on Transition Metals Alloyed with Platinum and Palladium

Direct formic acid fuel cells (DFAFCs) have emerged as a promising power source to meet increased demands for alternative energy sources in the transportation and portable energy storage sectors. Furthermore, these fuel cells utilize formic acid (FA), a nontoxic and carbon-neutral fuel when produced from biomass or via CO 2 reduction. Despite the promise of DFAFCs, the best monometallic catalysts, platinum and palladium, are poisoned by CO through the indirect oxidation pathway and require large overpotentials. By alloying Pt and Pd with other metals, we aim to improve both the activity and selectivity of these catalysts. Here, we present a systematic density functional theory (GGA-PW91) study of FA electro-oxidation (FAO) on the (111) facet of bimetallic Pt (Pt 3 M) or Pd (Pd 3 M) catalysts (M = Au, Ag, Cu, Pt, Pd, Ir, Rh, Ru, or Re) to evaluate the catalytic performance of these surfaces. For each surface, we calculate free energy diagrams and onset potentials of three key reaction mechanisms: direct oxidation of FA via carboxyl (COOH), direct oxidation of FA via formate (HCOO), and the indirect oxidation of FA that first forms CO en route to full oxidation to CO 2 . We then display the trends in the form of phase diagrams that compare the activity of the calculated surfaces against regions of high activity using the free energies of adsorbed CO and OH as descriptors, enabling high-throughput screening and design of improved catalysts, particularly those alloying Pt or Pd with Ir, Ru, or Re.

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

NO Degradation on the Anatase TiO 2 (001) Surface in the Presence of Water

Nitric oxide (NO) is known to degrade to nitric acid (HNO 3 ) on anatase TiO2 facets under visible light irradiation in the presence of water. However, the exact role that water plays in this photoreaction is not fully understood. Here, by employing the density functional theory (DFT) and time-dependent density functional tight binding (TD-DFTB), we show the viability of two suggested degradation pathways involving water. Both reaction pathways are triggered by a charge transfer excitation from the NO molecule to the TiO 2 surface. In one of them, NO interacts with dissociated water molecules adsorbed on the surface to form HONO + , whereas for the second pathway, NO is oxidized to NO 2 + by capturing one oxygen atom from the substrate. Both HONO and NO 2 are known byproducts in the photodegradation of NO, which further react with adsorbed water to finally produce HNO 3 . We also demonstrate by means of nudged elastic band calculations that these reactions are unlikely to happen without illumination.

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

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110): 1. Thermodynamics

Ammonia (NH 3 ) is being considered as a practical means to transport hydrogen (H 2 ) because of its higher volumetric energy density for the same temperature and pressure. Thermodynamics suggest high temperature is needed to decompose NH 3 to nitrogen (N 2 ) and H 2 . Furthermore, overcoming decomposition kinetic barriers requires a catalyst. Via density functional theory, we study this reaction on a model catalyst: the close-packed (110) facet of α-Fe. Specifically, we predict detailed in-operando temperature- and pressure-dependent surface phase diagrams on this benchmark catalyst that offer insights for the design of optimal NH 3 decomposition catalysts. Here, we explore the equilibrium composition(s) of the Fe(110) surface when exposed to NH 3 –H 2 mixtures. Further, we predict that both N and NH partially cover the Fe(110) surface at 300–400 °C (far above the NH 3 decomposition-formation coexistence temperature at standard partial pressures of 1 bar: ~ 180 °C) and 2–4 bar of total reactor pressure. At the equilibrium N/NH coverage, these species inhibit coadsorption of H, indicating that direct H 2 production may occur. However, from thermodynamics alone, removal of N/NH as N 2 (g) is extremely unfavorable even at these elevated temperatures–effectively deactivating the surface toward further NH 3 decomposition. Thus, catalysts that remain active toward N–H bond breaking but with weakened N binding relative to Fe are needed to enable high-turnover catalytic NH 3 decomposition to release H 2 .

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

Role of Oxidizing Conditions in the Dispersion of Supported Platinum Nanoparticles Explored by Ab Initio Modeling

Achieving fine control over the dispersion of supported platinum nanoparticles (Pt) is a promising avenue to enhancing their catalytic activity and selectivity. Experimental observations suggest that exposing ceria-supported Pt nanoparticles to O 2 at 500 °C promotes their dispersion into smaller particles and eventually single atoms. The exact role of oxygen in this process is not yet well understood. Past density functional theory studies of ceria-supported Pt have typically narrowed their scope to single atoms and clusters of a few atoms. Herein, we combine several approaches and types of models in a consistent atomistic framework to evaluate the relative stability of ceria-supported Pt as a function of the degree of oxidation of Pt and of the particle size, ranging from single atoms to nanoparticles 1.5 nm of diameter. Finally, we find that the largest nanoparticles remain the thermodynamically most stable species on the lowest energy facet of ceria even under oxidizing conditions, suggesting that stronger adsorption sites are required to stabilize smaller clusters and single atoms and promote oxidative dispersion.

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

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

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

Evolution and Degradation Patterns of Electrochemical Cells Based on the Analysis of Interfacial Phenomena at Li Metal Anode/Electrolyte Interfaces

In this work, we report the results of a theoretical–computational analysis of the solid electrolyte interphase (SEI) growth and degradation dynamics occurring in lithium metal batteries during cycling. We use ab initio-kinetic Monte Carlo simulations to generate a synthetic data set, which is analyzed by machine learning methods. We aim to determine: (i) how modifications in interfacial interaction energies between solid electrolyte interphase (SEI) blocks and between Li ions and SEI facets impact the Coulombic efficiency (CE) of the battery and (ii) what factors, including reactions, microscopic transport, and other interfacial events, may lead to cell performance “failure” during prolonged charge and discharge cycles, signaled as a sharp decay in the CE over cycling. The demonstration of our approach is done on a cell including a Li metal surface interfacing with a previously introduced state-of-the-art electrolyte, and the idea can be applied to any electrochemical system. Outcomes include the identification of the leading chemical, physical, and structural variables causing cell failure and relating them to the electrolyte formulation, thus paving the way to future more refined analysis and electrolyte design.

batteries↗
DOE OSTI2026

Insights into Nonelectroactive C–C Bond Formation on Cu(100) during Electrochemical CO 2 Reduction from Multiconfigurational Wavefunction Theory

Carbon–carbon (C–C) bond formation is necessary for hydrocarbon (and oxygenate) synthesis beyond methane (and formate/formic acid) during electrochemical CO and CO 2 reduction (ECOR and ECO 2 R). Cu has notable ability to form hydrocarbons compared to other pure metals. In particular, the (100) facet of face-centered cubic Cu forms ethylene competitively with H 2 and methane during both ECOR and ECO 2 R. Past simulations based on density functional theory (DFT) with standard exchange-correlation functional approximations predict fast nonelectroactive C–C bond formation channels involving adsorbed (*) CO together with another *CO, formyl (*CHO), or hydroxymethylidyne (*COH), forming OC*–*CO, OC*–CHO*, and OC*–*COH, respectively. Such simulations support the prevailing hypothesis that emergence of C 2 products is kinetically determined at the early stages of the reduction chemistry. Here we show, via simulations with more accurate many-body, i.e., “correlated”, wavefunction theory (enabled by an embedding scheme), that the coupling of *CO with a *CO or a *COH (previously predicted at the same level of theory to kinetically dominate over *CHO as the one-electron reduction product of *CO) is highly activated (kinetically impeded), with free energy barriers >1 eV, in contradiction to previous DFT-based simulations. Intriguingly, we find that the coupling of two adjacent *COHs incurs only a small barrier (<0.3 eV) and is exoergic (< –1 eV); however, given the predicted low surface mobility of *COH, the emergence of HOC*–*COH is also improbable, at least at low *COH coverages. We therefore conclude that it is highly unlikely for *CO to participate in nonelectroactive C–C bond formation on pristine Cu(100), contrary to conventional wisdom, and that the energetically favorable *COH dimerization may occur only after substantial buildup of *COH on the surface.

adsorption↗
DOE OSTI2026

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2020

Chemisorption-Driven Roughening of Hydrothermally Grown KTa1–xNbxO3 Nanoparticles

The chemical and physical properties of nanoparticle surfaces have significant effects on their growth processes and the resulting morphology. Hydrothermally grown KNbO3, KTaO3, and KTa1–xNbxO3 were studied to examine the complex relationships between surface composition, phase, chemistry, and energetics and how these may be used to model and thereby control nanoparticle growth mechanisms. Two different composition-dependent growth modes were identified, where one type formed smooth surface facets, while the other resulted in roughened nanoparticle morphologies. Electron microscopy characterization, density functional theory calculations, and mathematical growth models were used to illuminate the role of surface properties and chemisorption on the nanoparticle growth morphology. Surface energy reduction by chemisorption can increase the rate of terrace nucleation, driving the roughening of these lower energy nanoparticle surfaces.

Ly, Tiffany↗
DOE OSTI2020

The Prevalence of Anions at Plasmonic Nanojunctions: A Closer Look at P-Nitrothiophenol

We revisit the reductive coupling of p-nitrothiophenol (NTP) to form dimercaptoazobenzene (DMAB), herein monitored through gap-mode tip-enhanced Raman spectroscopy (TERS) and nano-imaging. We employ a plasmonic Au probe (100 nm diameter at its apex) illuminated with a 633 nm laser source (50 ?W/?m2 at the sample position) to image a NTP-coated faceted silver nanoparticle (~70 nm diameter). A detailed analysis of the recorded spectra reveals that anionic NTP species contribute to the recorded spectral images, in addition to the more thoroughly described DMAB product. Notably, the signatures of the anions are more pronounced than those of the DMAB product under our present experimental conditions. Our results thus demonstrate that anions and their spectral signatures must be considered in the analysis of plasmon-enhanced spectra and images.

Wang, Chih-Feng↗