Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Grain Boundary‐Derived Cu + /Cu 0 Interfaces in CuO Nanosheets for Low Overpotential Carbon Dioxide Electroreduction to Ethylene

Abstract Electrochemical CO 2 reduction reaction can be used to produce value‐added hydrocarbon fuels and chemicals by coupling with clean electrical energy. However, highly active, selective, and energy‐efficient CO 2 conversion to multicarbon hydrocarbons, such as C 2 H 4 , remains challenging because of the lack of efficient catalysts. Herein, an ultrasonication‐assisted electrodeposition strategy to synthesize CuO nanosheets for low‐overpotential CO 2 electroreduction to C 2 H 4 is reported. A high C 2 H 4 Faradaic efficiency of 62.5% is achieved over the CuO nanosheets at a small potential of −0.52 V versus a reversible hydrogen electrode, corresponding to a record high half‐cell cathodic energy efficiency of 41%. The selectivity toward C 2 H 4 is maintained for over 60 h of continuous operation. The CuO nanosheets are prone to in situ restructuring during CO 2 reduction, forming abundant grain boundaries (GBs). Stable Cu + /Cu 0 interfaces are derived from the low‐coordinated Cu atoms in the reconstructed GB regions and act as highly active sites for CO 2 reduction at low overpotentials. In situ Raman spectroscopic analysis and density functional theory computation reveal that the Cu + /Cu 0 interfaces offer high *CO surface coverage and lower the activation energy barrier for *CO dimerization, which, in synergy, facilitates CO 2 reduction to C 2 H 4 at low overpotentials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electronic structure, magnetic properties, and exchange splitting of gadolinium intermetallics

Here, we report here the electronic structure calculations of the magnetic moment, magneto-crystalline anisotropy energy (MAE), and density of states (DOS) of Gd-metal, GdAl 2 , GdAl 3 , GdCo 2 , GdCo 3 , and GdCo 5 . From the DOS, the exchange splitting energy is derived to map ferrimagnetic to paramagnetic phase transition (Curie temperature, T C ). The employed theoretical method here is based on advanced density functional theory, specifically the linearized augmented plane wave method (LAPW) within the local spin density approximation (LSDA), including spin–orbit coupling (SOC) and Hubbard model parameter (U). DOS analysis of both Gd-Al and Gd-Co systems reveals a strong correlation between exchange splitting energy and T C . The Gd-Co system exhibits an enormous exchange splitting energy among Co-3d electrons, vastly overwhelming the vanishing Gd 5d exchange splitting. Finally, these findings suggest that exchange splitting of Co-3d electrons plays a key role in the ground state magnetism and magnetic transition temperatures in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic structure of LaNiO 2 and CaCuO 2 from a self-consistent vertex-corrected GW approach

In this work, the electronic structure of one of the nickelates (LaNiO 2 ) and one of the cuprates (CaCuO 2 ) is studied with three self-consistent GW-based methods: scGW, sc(GW+vertex), and quasiparticle self-consistent GW. Low-energy features obtained in our study are in many respects similar to the features reported in previous density functional theory plus dynamical mean-field theory (DFT+DMFT) studies. Consistent with the DFT+DMFT conclusion, we find LaNiO 2 to be more correlated than CaCuO 2 . However, correlation effects included in our study change the DFT Fermi surface near the Γ point differently from that reported in DMFT studies. Features that are a few electronvolts away from the Fermi level are broader in our calculations than in the DFT+DMFT, which reflects the differences between the DFT and the GW methods. Our results are in qualitative agreement with previous G0W0 results, but the self-consistency results in quantitative differences. Generally, correlation effects are found to be sufficiently weak in both materials, which allows one to use totally ab initio diagrammatic approaches such as sc(GW+Vertex) and to avoid the methods with adjustable parameters (DFT+U or DFT+DMFT). However, the possibility of some strong correlations at low energy that cannot be captured by perturbative methods cannot be completely excluded. For instance, differences in the Fermi surface should be resolved, thus experimental studies are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Resonantly Enhanced Electromigration Forces for Adsorbates on Graphene

We investigate the electromigration forces for weakly bonded adsorbates on graphene by using density-functional based calculations. We find that the nature of electromigration forces on an adsorbate critically depends on the energy level alignment between the adsorbate state and the Fermi level of the graphene. For a resonant adsorbate, whose frontier orbitals lie close to the Fermi level, the electromigration force is dominated by the electron wind force that is strongly enhanced along the electron flow direction, irrespective of the sign of the adsorbate charge. For a nonresonant adsorbate, the electromigration force is essentially the direct force that depends on the adsorbate charge. We also show that the magnitude of electromigration forces can be continuously tunable through electrostatic gating for resonant adsorbates. Our results provide new insight for understanding and controlling how nanoscale objects behave in or on host materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

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

From Design to Device: Challenges and Opportunities in Computational Discovery of p -Type Transparent Conductors

A high-performance p -type transparent conductor (TC) does not yet exist but could lead to advances in a wide range of optoelectronic applications and enable new architectures for, e.g., next-generation photovoltaic (PV) devices. High-throughput computational material screenings have been a promising approach to filter databases and identify new p -type TC candidates and some of these predictions have been experimentally validated. However, most of these predicted candidates do not have experimentally achieved properties on par with n -type TCs used in solar cells and therefore have not yet been used in commercial devices. Thus, there is still a significant divide between transforming predictions into results that are actually achievable in the laboratory and an even greater lag in scaling predicted materials into functional devices. In this perspective, we outline some of the major disconnects in this materials discovery process—from scaling computational predictions into synthesizable crystals and thin films in the laboratory to scaling laboratory-grown films into real-world solar devices—and share insights to inform future strategies for TC discovery and design. Published by the American Physical Society 2024

14 SOLAR ENERGY↗
DOE OSTI2023

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward using collective x-ray Thomson scattering to study C–H demixing and hydrogen metallization in warm dense matter conditions

The insulator–metal transition in liquid hydrogen is an important phenomenon to understand the interiors of gas giants, such as Jupiter and Saturn, as well as the physical and chemical behavior of materials at high pressures and temperatures. Here, the path toward an experimental approach is detailed based on spectrally resolved x-ray scattering, tailored to observe and characterize hydrogen metallization in dynamically compressed hydrocarbons in the regime of carbon–hydrogen phase separation. With the help of time-dependent density functional theory calculations and scattering spectra from undriven carbon samples collected at the European x-ray Free-Electron Laser Facility (EuXFEL), we demonstrate sufficient data quality for observing C–H demixing and investigating the presence of liquid metallic hydrogen in future experiments using the reprated drive laser systems at EuXFEL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probing Mg Intercalation in the Tetragonal Tungsten Bronze Framework V 4 Nb 18 O 55

While commercial Li-ion batteries offer the highest energy densities of current rechargeable battery technologies, their energy storage limit has almost been achieved. Therefore, there is considerable interest in Mg batteries, which could offer increased energy densities in comparison to Li-ion batteries if a high-voltage electrode material, such as a transition-metal oxide, can be developed. However, there are currently very few oxide materials which have demonstrated reversible and efficient Mg 2+ insertion and extraction at high voltages; this is thought to be due to poor Mg 2+ diffusion kinetics within the oxide structural framework. Here, in this paper, the authors provide conclusive evidence of electrochemical insertion of Mg 2+ into the tetragonal tungsten bronze V 4 Nb 18 O 55 , with a maximum reversible electrochemical capacity of 75 mA h g –1 , which corresponds to a magnesiated composition of Mg 4 V 4 Nb 18 O 55 . Experimental electrochemical magnesiation/demagnesiation revealed a large voltage hysteresis with charge/discharge (1.12 V vs Mg/Mg 2+ ); when magnesiation is limited to a composition of Mg 2 V 4 Nb 18 O 55 , this hysteresis can be reduced to only 0.5 V. Hybrid-exchange density functional theory (DFT) calculations suggest that a limited number of Mg sites are accessible via low-energy diffusion pathways, but that larger kinetic barriers need to be overcome to access the entire structure. The reversible Mg 2+ intercalation involved concurrent V and Nb redox activity and changes in crystal structure, as confirmed by an array of complementary methods, including powder X-ray diffraction, X-ray absorption spectroscopy, and energy-dispersive X-ray spectroscopy. Consequently, it can be concluded that the tetragonal tungsten bronzes show promise as intercalation electrode materials for Mg batteries.

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

Taming the fixed-node error in diffusion Monte Carlo via range separation

By combining density-functional theory (DFT) and wave function theory via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction known as configuration interaction using a perturbative selection made iteratively (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of mu = 0.5 bohr(-1) provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.

Scemama, Anthony↗
DOE OSTI2025

The electronic structure, crystal fields, and magnetic anisotropy in RECo 5 magnets

The current progress in describing rare-earth-based magnets' electronic structure and magnetic properties is discussed. We use several currently popular electronic structure methods to show the typical values of critical parameters that define the physics of RECo 5 (RE = rare earth atom) materials. The magnetic moments and magnetic anisotropy of 4f atoms are obtained using several approaches, including anisotropic 4f-charge density-constrained DFT and DFT+HI methods. We also suggest the introduction of "penalty" functional for obtaining correct variational total energy in the traditional Hund's rule-constrained DFT-based techniques. The applicability and future extensions are discussed. The proposed combination of methods is potentially suitable for high-throughput computational searches of new rare-earth-containing magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A comparative analysis of different van der Waals treatments for molecular adsorption on the basal plane of 2H-MoS 2

The binding energy of hydrogen sulfide, ammonia, ethane, ethylene, butadiene, benzene, toluene, pyridine, pyrrole, and thiophene on the basal plane of the semi-conducting 2H-molybdenum sulfide (MoS 2 ) was calculated with the following flavors of Density Functional Theory (DFT): GGA-PW91, PBE-D3, vdW-DF, optPBE, optB86b, optB88, vdW-TS, and BEEF-vdW. The GGA-PW91 binding energies are negligible (<0.07 eV in magnitude) in all cases. The predictions with vdW-DF and PBE-D3 are the closest (error <0.05 eV) to the isosteric heats of adsorption calculated from reported temperature programmed desorption data for thiophene and butadiene. For all dispersion flavors examined here, the magnitude of the dispersion contribution to the binding energy increases linearly with the number of heavy atoms in the adsorbate, with each atom contributing 0.05 eV (BEEF-vdW) – 0.09 eV (optB88-vdW). Further, this implies that the calculated adsorption constants of molecules larger than acridine (i.e., comprising > 14 non-heavy atoms) can vary by more than four orders of magnitude at industrial conditions depending on the chosen method of dispersion correction. Further, dispersion effects fall off rapidly (>0.03 eV/ non-hydrogen atom/Å) as the adsorbate-surface distance increases.

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

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fission with Exotic Nuclei (Abbreviated Report)

Nuclear fission is a key mechanism involved in the synthesis of heavy elements in the Cosmos and is the primary explanation for the stability of superheavy elements. Nevertheless, our knowledge of fission remains extremely fragmented. Most experiments have been conducted only on a tiny number of stable actinide nuclei and are often incomplete, leading to gaps in our basic understanding of the process. For many radioactive isotopes, basic fission data such as the charge or mass distribution of the fragments is unknown. These gaps cannot always be filled by simulation alone. Common fission models contain too many free parameters and lack predictive power. In contrast, the fundamental theory of fission under development at LLNL is much more predictive, but its current computational cost is too high to be used extensively for data evaluations. A unique window of opportunity to resolve these limitations has recently opened: the U.S. nuclear science community is ramping up major experimental programs at the Facility for Rare Isotope Beams (FRIB, the DOE flagship facility in low-energy nuclear science), and techniques from machine learning have shown great potential to simplify the use of a fundamental, quantum-mechanical theory of fission. This project has two components. On the experimental side, we acquired and deployed at the HIGS facility a new dual Frisch-Grid ionization chamber to measure correlated fragment-mass, kinetic energy, and angular distributions of fission fragments from induced fission. This new device was used to perform measurements of charge, mass and total kinetic energy of fission fragments in the photofission of 238 U and eight gamma-ray beam energies between 6.2 and 13 MeV, which allowed extracting high-precision independent yields for this reaction. The device was also used to perform measurements of the same quantities in the neutron-induced fission of 234 U with monoenergetic beams of energy between 5 and 8 MeV. In parallel, we collaborated with a team at Commissariat à l’énergie atomique et aux énergies alternatives (CEA) to perform a series of measurements of fission yields in inverse kinematics for the two isotopes of 236 U and 240 Pu. The experiment took place at the Grand Accélérateur National d’Ions Lourds in France in June 2023. The deployment of the VAMOS spectrometer with a new array called PISTA allowed determining the excitation energy of the fissioning system within 1 Mega-electronvolts. The second component of the project involved using deep neural networks to build fast and reliable emulators of our current fission models. In an invited paper published in Frontier in Physics, we showed that autoencoders could successfully compress nuclear wavefunctions in nuclear density functional theory. We achieved a dimensionality reduction of the order of two orders of magnitude while keeping the error in the total energy to less than 0.01%. In a second paper submitted to Physical Review Letters in June 2023 with our collaborators at CEA, we showed that variational autoencoders can learn the collective degrees of freedom driving the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Origin of dislocation structures in an additively manufactured austenitic stainless steel 316L

In this experiment, the origin of dislocation structures in AM stainless steels was systematically investigated by controlling the effect of thermal stress through geometric constraints for the first time. Stainless steel 316L parts were produced in the form of “1D” rods, “2D” walls, and "3D" rectangular prisms to evaluate the effect of constraints to thermal expansion/shrinkage on the development of defect microstructures and to elucidate the origin of additively manufactured (AM) dislocation microstructures. Dislocation density, organization, chemical micro-segregation, precipitate structures, and misorientations were analyzed as a function of increasing constraints around solidifying material in 1D, 2D, and 3D components built using both directed energy deposition (DED) and powder-bed selective laser melting (SLM). In DED parts, the dislocation density was not dependent on local misorientations or micro-segregation patterns, but evolved from approximately ρ ⊥ = 10 12 m -2 in 1D parts to ρ ⊥ = 10 14 m -2 in 3D parts, indicating that it is primarily thermal distortions that produce AM dislocation structures. In DED 3D parts and SLM parts, dislocation densities were highest (ρ ⊥ ≈ 10 14 m -2 ) and corresponded to the formation of dislocation cells approximately 300-450 nm in diameter. Dislocation cells overlapped with dendrite micro-segregation in some but not all cases. The results illustrate that dendritic micro-segregation, precipitates, or local misorientations influence how the dislocations organize during processing, but are not responsible for producing the organized cell structures. Furthermore, this work shows that AM dislocation structures originate due to thermal distortions during printing, which are primarily dictated by constraints surrounding the melt pool and thermal cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bond Selective Photochemistry at Metal Nanoparticle Surfaces: CO Desorption from Pt and Pd

The use of visible photon fluxes to influence catalytic reactions on metal nanoparticle surfaces has attracted attention based on observations of reaction mechanisms and selectivity not observed under equilibrium heating. These observations suggest that photon fluxes can selectively impact the rates of certain elementary steps, creating nonequilibrium energy distributions among various reaction pathways. However, quantitative studies validating these hypotheses on metal nanoparticle surfaces are lacking. Here, we examine the influence of continuous wave visible photon fluxes on the CO desorption rates from 1 to 2 nm diameter Pt and Pd nanoparticle surfaces supported on γ-Al 2 O 3 . Temperature-programmed desorption measurements quantified via diffuse reflectance infrared Fourier transform spectroscopy demonstrate that visible photon fluxes significantly enhanced the rate of CO desorption from Pt nanoparticles in a wavelength-dependent manner. 440 nm photons most efficiently promoted CO desorption from Pt nanoparticle surfaces, aligning with the excitation energy for the interfacial electronic transition within the Pt–CO bond. Conversely, visible photon fluxes had no measurable influence on CO desorption rates from Pd nanoparticle surfaces after accounting for photon-induced heating. Density functional theory calculations demonstrate that the Pt–CO bond exhibits a narrower LUMO resonance, stronger coupling between the photoexcitation and forces induced on the metal–C bond, and vibrational energy dissipation that more effectively couples to desorption as compared to Pd–CO. These results demonstrate the specificity photons provide in facilitating chemical reactions on metal nanoparticle surfaces and substantiate the idea that photon fluxes can steer processes and outcomes of catalytic reactions in ways not achievable by equilibrium heating.

desorption↗
DOE OSTI2024

Adsorbate-induced adatom formation on Au-Cu bimetallic alloys and its possible consequences for CO 2 electroreduction

The adsorbate-induced formation of sub-nanometer clusters on transition-metal single crystals observed in previous high-pressure microscopic studies hinted at the in-situ formation of unique active sites even on large nanoparticle catalysts. We propose that the adatom formation energy can be used as an energetic descriptor for the initial step toward the adsorbate-induced metal-cluster formation process. This descriptor can be efficiently computed using density functional theory (DFT) calculations and applied for screening and identification of metal catalysts where this phenomenon may play an important role in generating active sites in-situ. As a proof of concept, here, we construct an adatom formation energy database for three Au x Cu y alloys (x:y = 3:1, 1:1, or 1:3) and eighteen adsorbates (H, C, N, O, F, S, Cl, Br, I, CH x , NH x (x = 1 – 3), CO, NO, and OH) commonly involved in catalytic reactions. The energetics of adatom formation were examined in all cases where the (111) terrace, (211) step-edge, and (874) kink were the sources of the adatom. We demonstrate that the presence of an adsorbate could alter not only the energetics for adatom formation but also the elemental nature of the preferred adatom being formed. Using our database, we identified promising systems which favor adsorbate-induced adatom formation under near-ambient conditions. Specifically, CO-induced adatom formation on all three Au-Cu alloy surfaces could occur under CO 2 electroreduction (CO 2 RR) conditions. This phenomenon offers a qualitative explanation for the experimentally observed CO 2 RR activity on Au-Cu alloy catalysts. As a result, our methodology offers an easily expandable and efficient approach for large-scale catalyst screening with regards to adatom/cluster formation under reaction conditions and provides insight into the possible nature of active sites on alloy catalysts from a novel perspective.

Active site↗
DOE OSTI2020

Engineering Surface Oxygenated Functionalities on Commercial Carbon toward Ultrafast Sodium Storage in Ether-Based Electrolytes

The pursuit of a high-capacity anode material has been urgently required for commercializing sodium-ion batteries with a high energy density and an improved working safety. In the absence of thermodynamically stable sodium intercalated compounds with graphite, constructing nanostructures with expanded interlayer distances is still the mainstream option for developing high-performance carbonaceous anodes. In this regard, a surface-functionalized and pore-forming strategy through a facile CO 2 thermal etching route was rationally adopted to engineer negligible oxygenated functionalities on commercial carbon for boosting the sodium storage process. Benefitted from the abundant ionic/electronic pathways and more active reaction sites in the microporous structure with noticeable pseudocapacitive behaviors, the functionalized porous carbon could achieve a highly reversible capacity of 505 mA h g -1 at 50 mA g -1 , an excellent rate performance of 181 mA h g -1 at 16,000 mA g -1 , and an exceptional rate cycle stability of 176 mA h g -1 at 3200 mA g -1 over 1000 cycles. These outstanding electrochemical properties should be ascribed to a synergistic mechanism, fully utilizing the graphitic and amorphous structures for synchronous intercalations of sodium ions and solvated sodium ion compounds, respectively. Finally, the controllable generation and evolution of a robust but thin solid electrolyte interphase film with the emergence of obvious capacitive reactions on the defective surface, favoring the rapid migration of sodium ions and solvated species, also contribute to a remarkable electrochemical performance of this porous carbon black.

25 ENERGY STORAGE↗