Engineering Papers⌕ Search

Engineering topics

Weber, Adam Z

Publications and source records attributed to Weber, Adam Z.

At least 19 records

How Does Water Dissociation Work in Bipolar Membranes?

Bipolar membranes (BPMs) create counteracting spatial gradients of pH and electrostatic potential in electrochemical systems, enabling applications in pH regulation, electrocatalysis, and separations. At the polarized junction of a BPM the water dissociation (WD, 2H2O ⇌ H3O+ + OH-) reaction can be driven, but it remains poorly understood. In this Perspective, we integrate molecular insights from bulk-water autoionization and the associated field effects with continuum descriptions of BPM electrostatics and experimental WD kinetic analyses to describe possible mechanisms of voltage-driven WD. Pristine BPM junctions highlight both the limits of primarily electric-field-driven WD and the practical challenges of junction stability at extreme reverse bias. Introducing heterogeneous catalyst layers, commonly metal oxides and graphene oxides, accelerates WD by orders of magnitude through hypothesized coupled effects in which surface acid-base functionality and high-density hydroxyl sites mediate proton-transfer steps, and catalyst mobile electronic/ionic charges redistribute the junction electric potential drop to shape the local electric fields and reactive microenvironments. Kinetic analyses suggest two regimes of heterogeneous WD mechanism, including field-driven ordering of interfacial water and a Second-Wien-Effect dissociation-barrier lowering. We conclude by defining the key unknown variables (local pH, electrostatic potential, catalyst charge state and relationships among mechanisms) and outlining experimental and multiscale modeling strategies needed for predictive WD catalysis and for controlling related ion-transfer reactions.

Wu, Yifan↗

PFSA-Ionomer Adsorption to C and Pt/C Particles in Fuel-Cell Inks

Catalyst inks used to make fuel-cell electrodes consist of Pt/C catalyst particles and a perfluorosulfonic acid (PFSA) ionomer dispersed in water/alcohol solvent mixtures. PFSA ionomer in the ink adsorbs to the surface of the catalyst particles, dictating the dispersion colloid properties. Following adsorption, the subsequent distribution of excess nonadsorbed ionomer in the ink then governs the final structure of the electrode. Here, we characterize the adsorption of the PFSA ionomer onto Pt/C catalyst particles. PFSA adsorption is largely irreversible. Adsorbed sulfonic-acid moieties impart a negative charge on the catalyst surface, causing electrostatic repulsion between the free ionomer in solution and the ionomer-covered Pt/C particle surface. The amount of adsorption is limited by the resulting electrostatic charge that grows as more ionomer adsorbs, and the catalyst surface becomes more negatively charged. Attenuating electrostatic repulsion by increasing the ink ionic strength promotes ionomer adsorption. Electrostatically limited adsorption is observed, irrespective of the solvent water/n-propanol ratio or the catalyst particle porosity and Pt loading. Experimentally measured ionomer adsorption isotherms are well predicted by a Smoluchowski-based kinetic adsorption model, in which the electrostatic energy barrier for adsorption is predicted from DLVO theory. These findings help to unravel the complex phenomena within these colloidal dispersions, allowing for subsequent tailoring of inks to optimize fuel-cell electrode structure and performance.

Adsorption↗

Achieving the hydrogen shot: Interrogating ionomer interfaces

The aim of this study is to enable the hydrogen economy and decarbonize various sectors in our environment that requires less expensive and more durable water electrolyzers, which can meet the Hydrogen-Shot target. The key is to improve the ionomer interfaces in low-temperature water electrolyzers as rapidly as possible, but to do so, it requires a systematic and holistic campaign combining both experiments and theory. In this perspective, we discuss the issues of electrolyzers and needs for translational science. We then present the approach that the Energy EarthShot Research Center: Center for Ionomer-based Water Electrolysis is taking in hopes of inspiring the community with this approach that can be leveraged to multiple problems and technologies.Graphical abstractHighlightsOne way to achieve the Hydrogen-Shot goal of low-cost, clean hydrogen, is advancing research and development on the interfaces of water electrolyzers for both performance and lifetime. The Center for Ionomer-based Water Electrolysis is exploring new techniques and strategies to not only interrogate interfacial phenomena in water electrolyzers to increase efficiency and durability, but also a new paradigm related to synergistic, cojoined experimental and theoretical research.DiscussionCatalyst\ionomer interfaces are complex and not fully understood, but through investigating different interfaces and utilizing digital and physical twins, we can elucidate key mechanisms and understanding.Understanding the dynamic double layer in electrochemical systems that use solid electrolytes is crucial to identifying and mitigating the controlling phenomena to enable increased performance and durability at the technology level.Studying the time and length scales of interfacial changes can be a powerful tool to understand reaction mechanisms and changes in the electrolyzer performance and durability.

Fornaciari, Julie C↗

Modeling Framework for the Assessment of a Sustainable Hydrogen Production and Supply Chain Network in California

The cost-effective and sustainable deployment of hydrogen supply and demand networks, especially in large economic regions like California, can be challenging considering the spatial-temporal availability and variability of the different actors across the network such as production processes, distribution modes, and end-users. In this presentation, we will provide an overview and demonstration of a modeling framework used to assess the environmental, economic, and human health impacts of plausible hydrogen production and supply chain networks in California. Scenarios focus on green hydrogen production pathways using water electrolysis and biomass gasification. End-use applications included in the model are transit, medium and heavy-duty trucking, port authorities, and power and aviation companies that currently consume natural gas, diesel, and aviation fuel for their day-to-day operation. Representative locations for hydrogen production and end-use are based on recent projections of the hydrogen economy in California. All mass and energy flows, as well as estimated emissions, are based on H2A process model designs and projections of technology performance, literature review, and LBNL process, economic and life cycle modeling, and not on company data for the sake of this presentation. Human health impacts are included following methodologies developed for the University of California Irvine HyDeal project. Life cycle phases associated with hydrogen production include feedstock preparation (water and biomass), energy production and consumption (renewable, grid, and combination of renewable and grid electricity), maintenance (chemical utilization in electrolysis and natural gas combustion in gasification), carbon sequestration, hydrogen storage (compression and liquefaction), and distribution (truck and pipeline). We apply the framework utilizing California specific emission factors, financial data, and human health damages and explore the impact of network characteristics on results. Example variations include: the inclusion of policy incentives or not, different representations of the electricity grid and source, electrolysis versus gasification versus combinations of both for production, liquefaction versus compression based on producer capacity cutoffs, transportation truck versus pipeline based on existing infrastructure, and ultimate end use. Comparison of these different scenarios can help inform future projects by demonstrating the trade-offs among environmental, economic, and human health impacts. This model, automated in R, is a starting platform upon which new analysis, modeling capabilities, locations, and emission factors can be rapidly tested and integrated.

Zaki, Mohammed Tamim↗

Asymmetric Bipolar Membrane for High Current Density Electrodialysis Operation with Exceptional Stability

Bipolar membranes (BPMs) enable isolated acidic/alkaline regions in electrochemical devices, facilitating optimized environments for electrochemical separations and catalysis. For economic viability, BPMs must attain stable, high current density operation with low overpotentials in a freestanding configuration. We report an asymmetric, graphene oxide (GrOx)-catalyzed BPM capable of freestanding electrodialysis operation at 1 A cm–2 with overpotentials <250 mV. Use of a thin anion-exchange layer improves water transport while maintaining near unity Faradaic efficiency for acid and base generation. Voltage stability exceeding 1100 h with an average drift of 70 μV/h at 80 mA cm–2 and 100 h with an average drift of −300 μV/h at 500 mA cm–2 and implementation in an electrodialysis stack demonstrate real-world applicability. Continuum modeling reveals that water dissociation in GrOx BPMs is both catalyzed and electric-field enhanced, where low pK a moieties on GrOx enhance local electric fields and high pK a moieties serve as active sites for surface-catalyzed water dissociation. These results establish commercially viable BPM electrodialysis and provide fundamental insight to advance design of next-generation devices.

Lucas, Éowyn↗

Unraveling the core of fuel cell performance: engineering the ionomer/catalyst interface

The biggest obstacle to the widespread implementation of polymer electrolyte membrane fuel cells (PEMFCs) is their cost, primarily due to the use of platinum catalysts. The high intrinsic catalyst activity exhibited on a rotating disk electrode (RDE) is rarely realized in a membrane electrode assembly (MEA), which is a long-standing challenge for PEMFCs and a cause of low catalyst utilization. To translate the high RDE performance of a catalyst into a MEA, the design of an ideal ionomer/catalyst interface is proposed: a thin, conformal ionomer film covers the maximum surface of a Pt nanoparticle and thus simultaneously maximizes catalyst utilization, (i.e., high mass activity and electrochemically active surface area) and O 2 diffusion rate (i.e., high current density performance) without compromising proton conduction. Building such an interface is a long-standing challenge due to the lack of interaction between the ionomer and catalyst particles, resulting in large ionomer agglomerates and inhomogeneous ionomer coverage over the catalyst nanoparticle, with consequent poor fuel cell performance. In this work, this ionomer/catalyst interface has been engineered, utilizing the electrostatic attraction between positively charged catalyst and negatively charged ionomer particles in a catalyst ink and preserved in a solid catalyst layer. As a result, this interface leads to previously unachieved proton exchange membrane fuel cell performance in terms of both catalyst utilization (75% vs. 45%) and peak/rated power density (i.e., 1.430/0.930 W cm -2 , H 2 /air, cathode Pt loading: 0.1 mgPt cm -2 ) for pure Pt catalysts, even better than those of Pt alloy catalysts. This work demonstrates the formation of an interface in the liquid phase (using ultra-small-angle X-ray scattering in combination with cryo-TEM, isothermal–titration–calorimetry) and the preservation of the interface in the solid catalyst layer (using TEM) and estimates the effective coverage and thickness of the ionomer film (using limiting current density, RDE and fuel cell performance).

25 ENERGY STORAGE↗

Using Membranes with Internal Microchannels to Prevent Drying-out during CO 2 Electrolysis

Scaling up CO 2 electrolysis is a vital aspect in the transition to manufacturing sustainable fuels and chemical compounds, satisfying the demand for chemicals and demand for storing renewable electricity. Depending on the employed catalyst, different products can be produced, such as carbon monoxide and ethylene, by applying a voltage on a CO 2 and H 2 O fed electrolyzer. CO is a desirable product according to techno-economic analysis, because it can be produced selectively using a silver catalyst and is a precursor for hydrocarbons in the Fischer–Tropsch process. In a state-of-the-art CO 2 electrolyzer, two electrodes are directly pressed against an ion-exchange membrane – this is called a zero-gap configuration. Therefore, the membrane is a crucial component for the system, since it has the role of providing a conductive medium between the electrodes. One of the challenges in CO 2 electrolysis is that water is consumed in the reaction. At high current density, this may cause the membrane's surface near the cathode to dry out, lowering efficiency or perhaps stopping the process entirely, since there is no longer a conductive medium. As a result, water management is critical for this process. In this work, we approach the drying-out challenge by studying the novel concept of a membrane with internal microchannels. These channels allow the circulation of water or an electrolyte inside the membrane, which reduces the water diffusion path and affects the membrane’s conductivity. The effects of channel geometry, location, and concentrations of electrolyte inside on water content, conductivity and overall performance are studied in a 2D COMSOL model. In addition, the effect of internal concentration of electrolyte on the membrane’s resistance, on the process performance and the K + cross-over to cathode side were investigated experimentally. Our modeling results prove that the presence of the channels can keep the membrane hydrated. The highest current densities are observed when the channel is closest to the cathode, and with smaller pores. Smaller pores are advantageous due to the trade-off between enhanced membrane conductivity and the lower conductivity of the liquid itself. If water is circulated in a large channel, it increases the membrane’s ionomer conductivity due to hydration but the overall conductivity is decreased since water is not highly conductive. Nonetheless, the results also show that a higher concentration of electrolyte inside the microchannels can significantly increase the total conductivity of the membrane, and therefore the energy efficiency of the process. These effects are most significant at higher current densities. In the experimental results, we’ve observed similar effects in terms of membrane conductivity and current density of the process – the higher the electrolyte concentration the higher the current density. Furthermore, a low concentration decreases the amount of potassium which crosses over to the cathode side, inhibiting salt deposition. We’ve concluded that a small channel, up to 90 µm wide, close to the pore with an electrolyte with a concentration of up to 10 mM could be very beneficial for the water management and energy efficiency of the process. This helps to keep the membrane hydrated at higher current densities, improves the conductivity of the membrane, and it doesn’t have significant impacts on the salt deposition.

Petrov, Kostadin Veselinov↗

Enabling Efficient Water Splitting with Advanced Materials Designed for High pH Membrane Interface

This research focuses on durable, high-performance materials and interfaces for advanced water splitting, enabling a clear pathway for achieving <$2/KgH2 (on scale) via anion exchange membrane (AEM)-based electrolysis. We aim to advance this goal via an improved fundamental understanding of both hydrogen and oxygen evolution reactions (HER/OER) leading to platinum group metal (PGM)-free catalyst materials. Here, we use NiFeCo and NiMo as OER and HER catalysts, respectively, in a full alkaline electrolysis cell. A thermally stable, multi-cation AEM is used to operate the cell at elevated temperatures thereby lowering the operating potential by increasing the kinetics of the respective catalytic reactions. The transition metal catalysts, paired with this novel AEM, have been used to achieve an operating potential of 1.79 V at 1 A/cm2. This potential is 300 mV lower than similar electrolysis cells equipped with PGM-based catalysts. We have also demonstrated the efficacy of the less caustic potassium carbonate solution as a replacement for potassium hydroxide as an electrolyte. By using a carbonate-based electrolyte, the hydroxide ion required for the anodic reaction is continuously replenished by re-establishing an equilibrium with water. This consideration, coupled with the reduced alkalinity of the carbonate solution, yields an electrolysis cell capable of sustained operation with minimal increase in operating potential.

Kendrick, Ian↗

(Invited) Understanding Multi-Component Transport and Reaction Phenomena in Bipolar Membranes Used for Electrosynthesis

Bipolar membranes (BPMs), which have long seen usage in electrodialysis reactors for the generation of acid and base, have recently demonstrated potential to become critical components in electrochemical synthesis devices. Because they can operate under large pH gradients, BPMs enable favorable environments for electrocatalysis at the individual electrodes. Critical to the implementation of BPMs in these devices is understanding the kinetics of water dissociation that occurs within the BPM junction as well as the co- and counter-ion crossover through the BPM, which both present significant obstacles to developing efficient and stable BPM-devices for electrosynthesis applications. Prior work has modeled ion transport in bipolar membranes in neutral salt solutions for electrodialysis. However, no model exists for the BPM under the harsh applied pH gradients that would be present in electrosynthesis, and there is significant need to explore the effects of the internal hydration on the lifetime and performance of BPMs in such environments. Additionally, a mechanistic understanding of water dissociation catalysis will be required to develop interfacial catalysts that enable the high current density operation required for scalable electrosynthesis of fuels. In this talk, we discuss modeling methodologies and physics inherent in BPMs and present our recent model of ion transport and water dissociation catalysis in BPMs across the pH scale. Specifically, we simulate multi-ion transport for a BPM with various electrolyte combinations on each side of the membrane, demonstrating the significance of co- and counter-ion crossover in BPMs operating under harsh pH gradients. We then investigate effects caused by hydration gradients that occur due to internal ion-exchange and examine potential methods for improving performance and mitigating crossover. Finally, we examine the impact of the interfacial water dissociation catalyst and perform sensitivity analysis on the key properties (catalyst point of zero charge and pKa) that dictate catalyst performance. These results provide information that is critical to developing a comprehensive understanding of multi-component phenomena in BPMs and to informing the design and implementation of BPMs in next-generation devices for the numerous electrosynthesis chemistries that benefit from operation under an applied pH gradient.

Bui, Justin C↗

Sensitivity Analysis and Effective Parametrization of PEM Fuel Cell Models

The cost of proton-exchange-membrane fuel cells (PEMFCs) remains a major hurdle in large-scale commercialization of this technology. To improve their performance and reduce cost, novel materials and electrode designs are continuously envisioned, e.g., non-PGM catalyst layers, ultra-thin Pt/Pt-Ni based catalyst layers, structured ionomer arrays or NSTF catalyst layers.1 Understanding the impact of these improvement strategies can be extremely time and cost intensive due to complex physical phenomena and large design space. We have previously developed a PEMFC modeling framework2 which has been a time and cost effective tool for understanding and optimizing the complex multi-physics phenomena within PEMFCs; however, several of the cell parameters used in the modeling have large spread in measured data.3 Furthermore, several transport parameters such as water adsorption kinetics have not been accurately measured and the approximations are spread over several orders of magnitude. These uncertainties cause problems in ascertaining accuracy of the modeling approach and they reduce the predictive power of the numerical models. The aim of this work is to identify the sensitivity of PEFC numerical model outputs to various input parameters. The previously in-house developed MEA modeling framework2 is used for PEMFC modeling. The sensitivity of the model outputs with respect to inputs parameters is obtained by analyzing the condition numbers for different output-input pairs at varying operating conditions. An example of the sensitivity analysis is shown in Figure 1. The condition numbers are obtained for the entire possible range of input parameters at varying operating conditions to identify the most crucial parameters of the PEMFC model. Based on our preliminary analysis, parameters related to kinetics (exchange current density and ECSA) and heat/water management in electrodes and ionomer (ionomer fraction, thermal conductivity) are most crucial. One of the major goals of this work is to identify the most crucial set of parameters towards which the model shows maximum sensitivity. This will guide future experimentalists to measure these properties with higher accuracy. Furthermore, the sensitivity analysis will also enable us to optimize the PEMFC performance by selectively targeting the most sensitive parameters and thereby making the largest impact. Acknowledgements The work is funded under the Fuel Cell Performance and Durability Consortium (FC-PAD), by the Fuel Cell Technologies Office (FCTO), Office of Energy Efficiency and Renewable Energy (EERE), of the U.S. Department of Energy under contract number DE-AC02-05CH11231. The authors would like to thank Nathan Craig at Robert Bosch LLC for his valuable input in designing the sensitivity analysis. The authors would also like to thank Giovanna Bucci and Matthias Hanauer at Robert Bosch for their valuable inputs and discussion. References P. K. Sinha, W. Gu, A. Kongkanand and E. Thompson, J. Electrochem. Soc., 158, B831 (2011). L. M. Pant, M. R. Gerhardt, N. Macauley, R. Mukundan, R. L. Borup and A. Z. Weber, Electrochim. Acta, 326, 134963 (2019). R. Vetter and J. O. Schumacher, ArXiv181110091 Phys. (2018). Figure 1

Pant, Lalit↗

Influence of Proton Activity Gaps between Electrodes on Open-Circuit Potential of H 2 /H 2 and H 2 /Air Cells

Polymer -electrolyte-fuel-cell open-circuit voltages (OCVs) are exactly defined by equation (1), where cathode and anode proton activities [(aH+)cathode and (aH+)anode, respectively] usually are identical, so the third term in the right-hand side of equation (2) is ignored. OCV=E0+RT/2F*ln(a1/2 O2*(a2 H+)cathode/aH2O)-RT/2F*ln((a2 H+)anode/aH2) (1) =E0+RT/2F*ln(a1/2 O2*aH2/aH2O)+RT/2F*ln((a2 H+)cathode/(a2 H+)anode) (2) Water vapor pressure is a colligative property that fundamentally correlates to electrolyte concentrations in aqueous solutions. Proton activity is a function of acid concentration, such as pH, when electrolytes are acids. In polymer-electrolyte membranes, water vapor pressure and acid concentration are understood as relative humidity (RH) and water uptake (λ), respectively, where λ represents number of water molecules per sulfonic acid molecule. Several investigations have reported the relationship between RH and λ, meaning that proton activities and associated water uptakes are intimately related to RH. In actual fuel cell operation, cathode RH is determined by ambient-atmosphere and/or humidifier RH(s), and anode RH depends on hydrogen-circulator RH. Therefore, RH is not always identical at both electrodes, and the difference between electrode RHs is considerable during dry operation of polymer electrolyte fuel cells. Therefore, the third term in the right-hand side of equation (2) may be significant for dry operation. We measured OCVs when hydrogen was supplied to both electrodes at 80°C. One electrode (A) was fixed at 30% RH, while RH at the other electrode (B) was varied (0, 5, 10, 20, and 30%). Measured OCVs varied from 0 to 75 mV. For fuel cell tests, electrode A was supplied with hydrogen at 30% RH; electrode B, oxygen at 0, 5, 10, 20, and 30% RH. OCVs deviated from that measured when RH at electrode B was 30%, increasing from 0 to 60 mV with decreasing RH at electrode B. Results are also shown in Figure 1. Proton activities of both electrodes were thermodynamically calculated. The Gibbs–Duhem relation was applied to obtain molar Gibbs free energies of water and sulfonic acid, and proton activity coefficient was calculated using the Gibbs free energy of sulfonic acid and the relationship between RH and λ1–4, assuming that protons and sulfonic anions show identical ionic-activity coefficients. OCVs were estimated using the third term in the right-hand side of equation (2). Results are shown in Figure 1. Fuel-cell current–voltage performance was poor when RHs at the anode and cathode were 30 and 20%, respectively. To determine kinetic current, we measured the oxygen-reduction reaction (ORR) using a rotating-disk electrode (RDE) in concentrated-acid aqueous solutions, which modeled catalyst-layer ionomers. Kinetic currents decreased with acid concentrations. References T. A. Zawodzinski, Jr., C. Derouin, S. Radzinski, R. J. Sherman, V. T. Smith, T. E. Springer and S. Gottesfeld , J. Electrochem. Soc., 140,1041 (1993) P. K. Das and A. Z. Weber, Proceedings of the ASME 2013 11th Fuel Cell Science, Engineering and Technology Conference, Fuel Cell 18010 (2013) V. A. Sethuraman, J. W. Weidner, A. T. Haug, S. Motupally,b and L. V. Protsailo, J. Electrochem. Soc., 155, B50 (2008) A. Kusoglu and A. Z. Weber, Chem. Rev., 117, 987 (2017) Figure 1

Yoshida, Toshihiko↗

Mathematical Modeling of Hydroxide-Exchange-Membrane Water Electrolyzer

Water electrolyzers can transform intermittent renewable energy like solar energy and wind energy into the chemical energy of hydrogen with zero greenhouse-gas emissions. The hydroxide-exchange membrane electrolyzer (HEME) combines the capability to produce pressurized hydrogen with the advantage of being able to use low or non-platinum group metal (PGM) electrocatalysts in the alkaline environment.1 Hydroxide salts, for example, KOH, are added to the HEME water feed on both anode and cathode to improve its performance. However, the specific mechanism of performance improvement still needs to be further understood. In addition, at high current densities, bubble evolution can result in mass-transport limitations, a less well studied phenomena. Mathematical modeling is ideal to explore these issues as it is cost and time efficient and can deconvolute the physics, processes, and observed phenomena and study the applied-voltage breakdown. In this work, we extend our previously developed 1D two-phase continuum model2 to study the varies processes in the HEME and provide insights on performance optimizations. First, the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) kinetics at different hydroxide concentrations have been studied by rotating disk electrodes (RDE) and implemented in the model. Then, the model is calibrated and validated against experimental HEME polarization curves for different KOH concentrations as a liquid electrolyte. The model clearly shows a performance increase with increasing KOH concentrations, which is consistent with the experimental results. The reduced ohmic resistance and increased electrochemical active surface area (ECSA) are the two main reasons for performance increase. The large amount of hydroxide in the liquid electrolyte not only helps to distribute the reactant hydroxide throughout the catalyst layer (CL), which reduces ohmic loss, but also enables reaction at the interface between the liquid electrolyte and electrocatalyst, which increases the ECSA. Applied-voltage breakdown demonstrates that the electrolyzer performance is dominated by anode kinetics and ohmic loss. A comparison with the DI water feed shows a more uniform current distribution in the anode CL when KOH is added, which indicates a higher utilization of the CL. Second, we present modeling on the effects of bubble coverage. As gas evolves, part of the ECSA is minimized due to bubble coverage. To account for this effect, an empirical relationship between the fractional bubble coverage and the current density is implemented in the model.3 The model shows this bubble coverage effect is more pronounced at large current densities with DI water feed. Acknowledgements This work was funded under the HydroGEN Consortium by the Energy Efficiency and Renewable Energy, Hydrogen and Fuel Cell Technologies Office, of the U. S. Department of Energy under contract number DE-AC02-05CH11231. References R. Abbasi, B. P. Setzler, S. Lin, J. Wang, Y. Zhao, H. Xu, B. Pivovar, B. Tian, X. Chen, G. Wu and Y. Yan, 31, 1805876 (2019). L. N. Stanislaw, M. R. Gerhardt and A. Z. Weber, ECS Transactions, 92, 767 (2019). H. Vogt and R. J. Balzer, Electrochimica Acta, 50, 2073 (2005).

Liu, Jiangjin↗