Search

The aqueous sodium-ion battery is a promising alternative to the well-known lithium-ion battery owing to the large abundance of sodium ion resources. Although it is safer than the lithium-ion battery, the voltage window of the sodium-ion battery is narrower than that of the lithium-ion battery, thus limiting its practical implementation. Therefore, a highly concentrated electrolyte is required to address this issue. In the present work, the effect of the salt concentration on the transport properties of water molecules is investigated via theoretical analyses at the quantum mechanical level. A molecular dynamics simulation at the quantum mechanical level revealed that as the salt concentration increases, the ion-water interactions became stronger, leading to a lower diffusivity and a lower electronic band gap. These imply that the superconcentrated aqueous-based electrolytes have high potentials for the sodium-ion battery applications.

Ceria (CeO2) is a promising metal-oxide support that is used in three-way catalysis (TWC). The activity of ceria-supported TWC depends on the location and concentration of oxygen vacancies. Oxygen ...

The density-functional tight-binding (DFTB) method is one of the useful quantum chemical methods, which provides a good balance between accuracy and computational efficiency. In this account, we reviewed the basis of the DFTB method, the linear-scaling divide-and-conquer (DC) technique, as well as the parameterization process. We also provide some refinement, modifications, and extension of the existing parameters that can be applicable for lithium-ion battery systems. The diffusion constants of common electrolyte molecules and LiTFSA salt in solution have been estimated using DC-DFTB molecular dynamics simulation with our new parameters. The resulting diffusion constants have good agreement to the experimental diffusion constants.

The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice Ih and Ic) and high-density structures (ice III and melted ice VI). Analysis of the temperature dependence of radial distribution function and bond order parameters served to characterize the distribution and configuration of hundreds of water molecules in a unit cell. We confirmed that ice Ih and Ic possess higher hexagonal symmetries than ice III and melted ice VI. The estimated Grotthuss shuttling diffusion coefficients in ice were larger than that of liquid water, indicating a slower proton diffusion process in high-density structures than in low-density systems. The energy barriers calculated on the basis of the Arrhenius plot of diffusion coefficients were in reasonable agreement with experimental measurement for ice Ih. Furthermore,...

Predicting pKa values for different types of amine species with high accuracy and efficiency is of critical importance for the design of high performance and economical solvents in carbon capture and storage with aqueous amine solutions. In this study, we demonstrate that density-functional tight-binding-based metadynamics simulations are a promising approach to calculate the free energy difference between the protonated and neutral states of amines in aqueous solution with inexpensive computational cost. The calculated pKa values were in satisfactory agreement with the experimental values, the mean absolute deviation being only 0.09 pKa units for 34 amines commonly used in CO2 scrubbing. Such superior reproducibility and correlation compared to estimations by static quantum mechanical calculations highlight the significant effect of dynamical proton transfer processes in explicit solvent molecules for the improvement of the estimation accuracy.

Divide-and-conquer-type density-functional tight-binding molecular dynamics simulations of the CO2 absorption process in monoethanolamine (MEA) solution have been performed for systems containing thousands of atoms. The formation of carbamate anions has been widely investigated for neutral systems via ab initio molecular dynamics simulations, yet the present study is aimed at identifying the role of hydroxide ions in acid-base equilibrium. The structural and electronic analyses reveal that the hydroxide ion approaches, via Grotthuss-type shuttling, the zwitterionic intermediates and abstracts a proton from the nitrogen atom of MEA. We also estimated the fraction of reacted CO2 and carbamate formed at different initial CO2 concentrations that confirm a high absorbed CO2 concentration decreases the fraction of MEA(C) formed due to the abundance of MEA(Z) in the solution.

The diffusion of the hydroxide ion in bulk water was examined by linear-scaling divide-and-conquer density-functional tight-binding molecular dynamics (DC-DFTB-MD) simulations using three different-sized unit cells that contained 522, 1050, and 4999 water molecules as well as one hydroxide ion. The repulsive potential for the oxygen-oxygen pair was improved by iterative Boltzmann inversion, which adjusted the radial distribution function of DFTB-MD simulations to that of the reference density functional theory-MD one. The calculated diffusion coefficients and the Arrhenius diffusion barrier were in good agreement with experimental results. The results of the hydroxide ion coordination number distribution and potential of mean force analyses supported a dynamical hypercoordination diffusion mechanism.

The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice I

Predicting pK

The process of proton diffusion in liquid water was investigated using molecular dynamics (MD) simulations, and the total energy and atomic forces were evaluated by the divide-and-conquer-type density-functional tight-binding (DC-DFTB) method. The effectiveness of this approach was confirmed by comparing the computational time of water clusters with conventional treatments. The unit cell employed herein, which contained 523 water molecules and 1 excess proton, was moderately large in comparison with those used in previous studies. The reasonable accuracy obtained by using this unit cell was confirmed by examining the temperature fluctuation. The diffusion coefficients for the vehicular and Grotthuss processes were accurately reproduced by the DC-DFTB-MD simulations with the unit cell containing 523 water molecules. Furthermore, the energy barriers were evaluated from the temperature dependence of the diffusion coefficient for each process. The calculated barrier for Grotthuss diffusion was in good agreement with the experimental value.