SEARCH · Engineering Papers

Results for “hydrogen binding energy”

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 OSTI2020

In silico discovery of active, stable, CO-tolerant and cost-effective electrocatalysts for hydrogen evolution and oxidation

Various databases of density functional theory (DFT) calculations for materials and adsorption properties are currently available. Using the Materials Project and GASpy databases of material stability and binding energies (H* and CO*), respectively, we evaluate multiple aspects of catalysts to discover active, stable, CO-tolerant, and cost-effective hydrogen evolution and oxidation catalysts. Finally, we suggest a few candidate materials for future experimental validations. We highlight that the stability analysis is easily obtainable but provides invaluable information to assess thermodynamic and electrochemical stability, bridging the gap between simulations and experiments. Furthermore, it reduces the number of expensive DFT calculations required to predict catalytic activities of surfaces by filtering out unstable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Fluorinated Covalent Organic Frameworks: A Novel Pathway to Enhance Hydrogen Sorption and Control Isosteric Heats of Adsorption; HyMARC Seed Project Final Report

Organic framework materials, as a unique sub-class of carbon-based sorbents, have gained increasing attention for promising attributes toward gas storage. Calculations of the hydrogen capacity of optimized frameworks show potential to achieve greater than 60 g/L storage of hydrogen, placing metal-organic frameworks (MOFs) and COFs near the top of the class of porous materials. However, at the ensemble level, poor stacking creates a quasi-amorphous material with low structural integrity and low effective surface area. The chemical versatility of COFs allows for additional methods for producing long-range order that specifically target the inter-layer interactions in COFs but that leave key pore-accessible linker sites open. These modifications (either during synthesis or post-synthetic) enable a host of strategies to both improve crystalline order for more stable and higher surface, as well as add metals with open coordination for enhanced H 2 binding enthalpy.

08 HYDROGEN↗

NASA NTRS1985

Adsorption and recombination of hydrogen atoms on a model graphite surface

The adsorption and recombination of atomic hydrogen on a model graphite grain have been examined in a series of calculations in which a modified, iterative, extended Hueckel program was used. The hydrogen atom is found to be chemisorbed at a site with a zero-point binding energy of 0.7 eV and at an equilibrium distance of 2.25 A above the site. Despite a barrier of about 0.4 eV between adjacent sites, calculations suggest that at temperatures as low as 10 K, an H atom will tunnel through to adjacent sites in less than one nanosecond. However, a potential barrier to the recombination of two hydrogen atoms has been found which displays high sensitivity to the mutual arrangement of the two hydrogen atoms with respect to the graphite surface. Results show that at very low temperatures, recombinations can occur only by tunneling. Consistent with experiment, the region in which H2 begins to form exhibits a repulsive potential with respect to possible chemisorption of the incipient H2 entity.

Aronowitz, S.↗

DOE OSTI2026

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong↗

DOE OSTI2020

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Influence of Heat Treatments in H 2 and Ar on the E 1 center in β -Gs 2 O 3

In this report the influence of heat treating n-type bulk β-Ga 2 O 3 in hydrogen (H2) and argon (Ar) gases on the presence of the defect level commonly labeled as E 1 was studied. Fourier transform-infrared spectroscopy confirms that hydrogen (H) is incorporated into β-Ga 2 O 3 during H 2 annealing at 900 °C. Deep-level transient spectroscopy measurements reveal that the concentration of the E 1 level is promoted by the introduction of H, in contrast to what is observed in samples heat-treated in an Ar flow. We further find the E 1 level to be stable against heat treatments at 650 K, both with and without an applied reverse-bias voltage. Potential candidates for the defect origin of E 1 are investigated using hybrid-functional calculations, and three types of defect complexes involving H are found to exhibit charge-state transition levels compatible with E 1 , including substitutional H at one of the threefold coordinated O sites, Ga-substitutional shallow donor impurities passivated by H, and certain configurations of singly hydrogenated Ga–O divacancies. Among these types, only the latter exhibit H binding energies that are consistent with the observed thermal stability of E 1 .

36 MATERIALS SCIENCE↗

DOE OSTI2023

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2–25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Electrochemical CO 2 Reduction Reaction over Cu Nanoparticles with Tunable Activity and Selectivity Mediated by Functional Groups in Polymeric Binder

Electrochemical carbon dioxide reduction reaction (CO 2 RR) using copper (Cu)-based catalysts has received significant attention mainly because Cu is an element capable of producing hydrocarbons and oxygenates. One possible way to control the CO 2 RR performance at the electrode interface is by modifying catalysts with specific functional groups of different polymeric binders, which are necessary components in the process of electrode fabrication. However, the modification effect of the key functional groups on the CO 2 RR activity and selectivity is poorly understood over Cu-based catalysts. In this work, the role of functional groups (e.g., -COOH and -CF 2 groups) in hydrophilic and hydrophobic polymeric binders on the CO 2 RR of Cu-based catalysts is investigated using a combination of electrochemical measurements, in-situ characterization and density functional theory (DFT) calculations. DFT results reveal that functional groups influence the binding energies of key intermediates involved in both CO 2 RR and the competing hydrogen evolution reaction, consistent with experimental observation of binder-dependent product distributions among formic acid, CO, CH 4 , and H 2 . This study provides a fundamental understanding that the selection of desired polymeric binders is a useful strategy for tuning the CO 2 RR activity and selectivity.

25 ENERGY STORAGE↗

DOE OSTI2021

Hydrogenation and C-S bond activation pathways in thiophene and tetrahydrothiophene reactions on sulfur-passivated surfaces of Ru, Pt, and Re nanoparticles

Thiophene-H2 reactions proceed via sulfur removal and hydrogenation routes on dispersed metal nanoparticles that become decorated by refractory S-adlayers during catalysis. The identity and kinetic relevance of the required elementary steps are described here based on rates measured at S-chemical potentials set by H2S/H2 ratios similar to those prevalent during practical catalysis on Re, ReSx, Ru, and Pt catalysts. Free energies for S adatom formation (from H2S decomposition and H2 evolution) are strongly exothermic (< -50 kJ mol-1 on Pt(111) and < -150 kJ mol-1 on Re and Ru(0001)), but strong repulsions between S adatoms cause adsorption free energies to increase significantly with coverage on all three surfaces, preventing complete monolayer formation. These adlayers, composed of unreactive S-atoms (S') that cover 1/3–2/3 ML leave residual interstitial spaces (*) that bind S-atoms (S*), intermediates, and transition states reversibly, as required for catalytic turnovers. The number and binding properties of these interstices depend on the identity and chemical state of the nanoparticle bulk phase, which influences S'-binding and coverages and cause large differences in direct desulfurization and hydrogenation turnover rates (per exposed metal atom) on dispersed Re, ReSx, Ru, and Pt. The identity and kinetic relevance of elementary steps for desulfurization (to C4¬ hydrocarbons) and hydrogenation (to tetrahydrothiophene; THT) are similar among these catalysts; they involve the kinetically-relevant formation of a thiophene-derived intermediate (monohydrothiophene on Re and ReSx; dihydrothiophene on Ru and Pt) that either cleaves its C-S bond or “over-hydrogenates” to THT in one surface sojourn. THT then undergoes C-S bond cleavage in secondary reactions that correct such over-hydrogenation to form the more unsaturated species that cleave C-S bonds. THT/C4 product ratios are insensitive to H2S/H2 ratios and thiophene pressure, even though active interstitial spaces are covered by kinetically-detectable coverages of S* and thiophene; therefore, primary and secondary reactions must involve the same active surfaces. The observed increase in THT/C4 ratios with H2 pressure shows that THT formation transition states involve a larger number of H-atoms than for C-S cleavage. The requirement for bound species with intermediate unsaturation (between THT and thiophene) for C-S bond cleavage is reminiscent of the H-shuttling required in C-C and C-O hydrogenolysis, reactions that involve the partial dehydrogenation of alkanes and alkanols, respectively, to weaken such bonds and to increase the formation entropy of the relevant transition states via the evolution of H2(g). These mechanistic details challenge prevalent paradigms about different site requirements for hydrogenation and desulfurization pathways and about how metal-sulfur bond energies act as descriptors of reactivity; in fact, such binding energies merely act to define the refractory S-adlayers that enable the formation of weakly-binding interstices that reversibly bind intermediates and transition states, thus allowing catalytic turnovers.

Yik, Edwin↗

DOE OSTI2019

Structure and regulation of human epithelial cell transforming 2 protein

Epithelial cell transforming 2 (Ect2) protein activates Rho GTPases and controls cytokinesis and many other cellular processes. Dysregulation of Ect2 is associated with various cancers. Here, we report the crystal structure of human Ect2 and complementary mechanistic analyses. The data show the C-terminal PH domain of Ect2 folds back and blocks the canonical RhoA-binding site at the catalytic center of the DH domain, providing a mechanism of Ect2 autoinhibition. Ect2 is activated by binding of GTP-bound RhoA to the PH domain, which suggests an allosteric mechanism of Ect2 activation and a positive-feedback loop reinforcing RhoA signaling. Furthermore this bimodal RhoA binding of Ect2 is unusual and was confirmed with Förster resonance energy transfer (FRET) and hydrogen–deuterium exchange mass spectrometry (HDX-MS) analyses. Several recurrent cancer-associated mutations map to the catalytic and regulatory interfaces, and dysregulate Ect2 in vitro and in vivo. Together, our findings provide mechanistic insights into Ect2 regulation in normal cells and under disease conditions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Catalytic influence of light element incorporation in the lattice of palladium

The alloying of metals with foreign elements is a common approach for manipulating catalytic properties of active sites. For palladium-based catalysts, incorporation of light elements receives far less attention than Pd-metal alloys even though incorporated light element are catalytically significant. Due to their small size, light elements (H, B, C, and O) can occupy octahedral sites in Pd to form an interstitial phase, while palladium adjusts its lattice to accommodate the incorporation of larger light elements (S and P) to form distinct phases. The incorporated light element modifies the electronic structure of palladium, influencing the binding energy of substrates and reactive intermediates. Pd hydride and oxide often form in situ during hydrogenation or oxidation reactions, respectively, representing the (more, less) active phase of the catalyst. Further, the reactive interstitial hydrogen dissolved in Pd tends to hydrogenate reactants completely and unselectively. Pd oxide can be the active phase for oxidation reactions, and the incorporated oxygen or oxygen vacancies may participate in the reaction directly. Incorporated light elements change geometric aspects of Pd ensembles, serving the role of site blocking, ensemble isolation, or serve as an additional active site, endowing Pd with bifunctional character. We review the synthetic approaches for incorporating light elements into the palladium lattice. We summarize the catalytic influence of light element incorporation into the palladium lattice with respect to electronic and geometric modification, in situ incorporated hydrogen and oxygen, and discuss the influence of light element incorporation into Pd on specific reactions (i.e., selective hydrogenation, oxidation, electrocatalysis, and C—C cross coupling).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Development of Density Functional Tight-Binding Parameters Using Relative Energy Fitting and Particle Swarm Optimization

We provide a strategy to optimize density functional tight-binding (DFTB) parameterization for the calculation of the structures and properties of organic molecules consisting of hydrogen, carbon, nitrogen, and oxygen. We utilize an objective function based on similarity measurements and the Particle Swarm Optimization (PSO) method to find an optimal set of parameters. This objective function considers not only the common DFTB descriptors of binding energies and atomic forces but also incorporates relative energies of isomers into the fitting procedure for more chemistry-driven results. The quality in the description of the binding energies and atomic forces is measured based on the Ballester similarity index and relative energies through a similarity index induced by the Levenshtein edit distance to quantify the correct energetic order of isomers. Training and testing datasets were created to include all relevant chemical functional groups. Finally, the accuracy of this strategy is assessed, and its range of applicability is discussed by comparison against our previous parameterization. The improved performance of the new DFTB parameterization is validated with respect to the density functional theory large datasets QM-9 and ANI-1, where excellent agreement is found between the structures and properties available in these datasets, and the ones obtained with DFTB.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1989

Thermonuclear breakup reactions of light nuclei. I - Processes and effects

Temperature and density conditions are considered for the occurrence of breakup reactions of light nuclei in astrophysical plasmas. The proton-induced endothermic process is shown to be the principal mechanism for nuclear breakdown in a plasma. The phenomenon occurs at a temperature of about 1 MeV, which is a fraction of the typical binding energy per nucleon in nuclei. The temperature for breakup of He-4 is about twice as large, because of the higher binding energy. Depending on the temperature attained in the plasma, the initial concentration of elements heavier than hydrogen can be depleted. However, if it attains a temperature of about 1 MeV, breaking up the metals (C, N, O, Ne, Mg) but not He-4, an increase in the He-4 abundance by as much as 10 percent can result, since these elements essentially break down to alpha particles.

Guessoum, Nidhal↗

NASA NTRS1989

A high level Ab initio study of the anionic hydrogen-bonded complexes FH-CN-, FH-NC-, H2O-CN- and H2O-NC-

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

NASA NTRS1989

A high-level ab initio study of the anionic hydrogen-bonded complexes FH-CN(-), FH-NC(-), H2O-CN(-), and H2O-NC(-)

Lee, Timothy J.↗

DOE OSTI2023

Analyzing barium titanate TiO2 surface interactions with tert-butylphosphonic acid using density functional theory

Abstract Barium titanate (BTO) is a widely researched ferroelectric useful for energy storage. While BTO’s surface chemistry is commonly studied using density functional theory, little has been published on the TiO2 surface. Here, we determined that BTO’s surface response can be decoupled from the ferroelectric response by using a pre-optimized ferroelectric slab and allowing only the top three atomic z-layers to respond to ligand binding. Multiple favorable binding modes were identified for hydrogen, hydroxyl, water, and tert-butyl phosphonic acid on BTO’s TiO 2 surface. Of these ligands, tBuPA dominates surface binding with binding energies as low as − 2.61 eV for its nine configurations. Graphical abstract

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Influence of Hydrogen Isotopes on Vacancy Formation and Antisite Defect Diffusion in Palladium and Vanadium Metals

Density functional theory calculations are performed to study the interaction between hydrogen isotopes, vacancy, and antisite defects in Pd and V. Various defect configurations and defect-defect distances are systematically explored. Binding energies and migration energy barriers are calculated and tabulated. The study provides atomistic data for subsequent mesoscale simulations of tritium, vacancy, and antisite diffusion. In Pd, a strong binding is found between a vacancy and tritium (0.16 eV). This tritium binding increases the thermal concentration of vacancies by a factor of ~10 at 500 ?C. The increase in vacancy concentration enhances V antisite diffusion in Pd by a factor of ~6 at 500 ?C. The influence of tritium is even stronger in V, with a tritium and vacancy binding energy of 0.38 eV. Such a strong binding increases the thermal concentration of vacancies by a factor of ~300 at 500 ?C. The increase in vacancy concentration enhances Pd antisite diffusion in V by a factor of ~640 at 500 ?C. Vanadium and Pd exhibit a strong driving force to intermix with a formation energy of -1.57 eV for V antisite in Pd and -1.05 eV for Pd antisite in V. The results suggest vanadium diffusion into Pd is energetically stronger than the reverse. Zero-point-energy corrections are taken into account and calculations for hydrogen and tritium are presented.

Setyawan, Wahyu↗

DOE OSTI2020

Determinants for proton location and electron coupled proton transfer in hydrogen bonded pentafluorophenol-anion clusters

This work aims to reveal the determinant factors for proton locations and electron coupled proton transfer (ECPT) in biologically relevant hydrogen bonded systems. Pentafluorophenol-anion clusters [C6F5O-•H+•A-]- are chosen to model active sites of biological functions, with the anion A- being systematically varied to ensure the proton affinities (PAs) of the anions well cover the referenced PA of C6F5O- from being appreciably smaller, similar, and to significantly larger. Negative ion photoelectron spectroscopy of these clusters provides spectroscopic evidence showing that proton location in the anionic state is largely following the PA prediction, while ECPT is observed only for the clusters with the anion possessing the electron binding energy (EBE) significantly larger than that of the referenced C6F5O-. Theoretical calculations suggest these clusters are stabilized by forming single strong hydrogen bond between donor and acceptor, and the associated charge and MO analyses fully support the experimental observations. The current holistic cluster model study indicates that PA is the right determinant that can be used to predict the proton location and describe hydrogen bonding structures, while both PA and EBE of anionic groups play important roles in facilitating ECPT process. This work was supported by U.S. Department of Energy (DOE), Office of Basic Energy Sciences, Division of Chemical Science, Geosciences, and Biosciences (X.-B.W.), and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. Part of this work was supported by the National Natural Science Foundation of China (No-11604046 for JZ; No. 11727810 for ZRS), as well as “Chenguang Program” by Shanghai Education Development Foundation, Shanghai Municipal Education Commission (No. 16CG38) and Natural Science Foundation of Shanghai (No. 16ZR1448100) (JZ).

Zhang, Jian↗