Jingxiang Xu

Multi-nanoparticle model simulations of the porosity effect on sintering processes in Ni/YSZ and Ni/ScSZ by the molecular dynamics method

Understanding the sintering mechanism in porous anodes is necessary for developing durable anodes suitable for use in solid oxide fuel cells. A multi-nanoparticle sintering simulation method based on molecular dynamics (MD) calculation was developed for this purpose [J. Xu et al., J. Phys. Chem. C, 2013, 117, 9663–9672]. The method can be used to calculate the effect of the porous structure properties, such as the porosity and framework structure, on the sintering, unlike previous sintering simulations with conventional nanoparticle models. We revealed that, in a Ni/YSZ porous anode, the YSZ nanoparticle framework suppresses sintering of Ni nanoparticles by disrupting the growth of the neck between two Ni nanoparticles. In this paper, we used our method to reveal the effect of ceramic type on the sintering processes. We investigated the difference between the sintering and degradation processes in Ni/YSZ and Ni/ScSZ anodes. In the simulation, the degree of sintering of the Ni nanoparticles in Ni/ScSZ was smaller than that in Ni/YSZ. The stronger adhesion of Ni to ScSZ nanoparticles than to YSZ nanoparticles prevented the Ni nanoparticles from approaching each other in the Ni/ScSZ anode, inhibiting sintering. Our multi-nanoparticle sintering MD simulations revealed the different sintering processes for Ni nanoparticles in Ni/YSZ and Ni/ScSZ anodes. We also investigated the effect of sintering on degradation. The hydrogen adsorption sites and electrochemical reaction sites of the hydrogen oxidation decreased as the degree of sintering increased. A low degradation of the Ni/ScSZ anode relative to that of the Ni/YSZ anode was observed. Furthermore, we showed the effect of porosity on degradation induced by sintering in Ni/YSZ and Ni/ScSZ, and found an optimal porosity. These findings cannot be obtained by conventional two- or three-nanoparticle sintering MD simulations. Our multi-nanoparticle sintering simulation method is useful for revealing the types of ceramics suitable for inhibiting sintering and degradation in anodes, and can be used to design durable anodes.

Communication: Different behavior of Young's modulus and fracture strength of CeO2: Density functional theory calculations

In this Communication, we use density functional theory (DFT) to examine the fracture properties of ceria (CeO2), which is a promising electrolyte material for lowering the working temperature of solid oxide fuel cells. We estimate the stress-strain curve by fitting the energy density calculated by DFT. The calculated Youngs modulus of 221.8 GPa is of the same order as the experimental value, whereas the fracture strength of 22.7 GPa is two orders of magnitude larger than the experimental value. Next, we combine DFT and Griffith theory to estimate the fracture strength as a function of a crack length. This method produces an estimated fracture strength of 0.467 GPa, which is of the same order as the experimental value. Therefore, the fracture strength is very sensitive to the crack length, whereas the Youngs modulus is not.

Tight-Binding Quantum Chemical Molecular Dynamics Study on the Friction and Wear Processes of Diamond-Like Carbon Coatings: Effect of Tensile Stress

Diamond-like carbon (DLC) coatings have attracted much attention as an excellent solid lubricant due to their low-friction properties. However, wear is still a problem for the durability of DLC coatings. Tensile stress on the surface of DLC coatings has an important effect on the wear behavior during friction. To improve the tribological properties of DLC coatings, we investigate the friction process and wear mechanism under various tensile stresses by using our tight-binding quantum chemical molecular dynamics method. We observe the formation of C-C bonds between two DLC substrates under high tensile stress during friction, leading to a high friction coefficient. Furthermore, under high tensile stress, C-C bond dissociation in the DLC substrates is observed during friction, indicating the atomic-level wear. These dissociations of C-C bonds are caused by the transfer of surface hydrogen atoms during friction. This work provides atomic-scale insights into the friction process and the wear mechanism of DLC coatings during friction under tensile stress.

Parallel Large-Scale Molecular Dynamics Simulation Opens New Perspective to Clarify the Effect of a Porous Structure on the Sintering Process of Ni/YSZ Multiparticles

Ni sintering in the Ni/YSZ porous anode of a solid oxide fuel cell changes the porous structure, leading to degradation. Preventing sintering and degradation during operation is a great challenge. Usually, a sintering molecular dynamics (MD) simulation model consisting of two particles on a substrate is used; however, the model cannot reflect the porous structure effect on sintering. In our previous study, a multi-nanoparticle sintering modeling method with tens of thousands of atoms revealed the effect of the particle framework and porosity on sintering. However, the method cannot reveal the effect of the particle size on sintering and the effect of sintering on the change in the porous structure. In the present study, we report a strategy to reveal them in the porous structure by using our multi-nanoparticle modeling method and a parallel large-scale multimillion-atom MD simulator. We used this method to investigate the effect of YSZ particle size and tortuosity on sintering and degradation in the Ni/YSZ anodes. Our parallel large-scale MD simulation showed that the sintering degree decreased as the YSZ particle size decreased. The gas fuel diffusion path, which reflects the overpotential, was blocked by pore coalescence during sintering. The degradation of gas diffusion performance increased as the YSZ particle size increased. Furthermore, the gas diffusion performance was quantified by a tortuosity parameter and an optimal YSZ particle size, which is equal to that of Ni, was found for good diffusion after sintering. These findings cannot be obtained by previous MD sintering studies with tens of thousands of atoms. The present parallel large-scale multimillion-atom MD simulation makes it possible to clarify the effects of the particle size and tortuosity on sintering and degradation.

Computational study on low friction mechanism of diamond-like carbon induced by oxidation reaction

Water lubrication has been attracting attention for environment-friendly society due to low CO2 emission. Furthermore, carbon-based materials such as diamond-like carbon (DLC) show the low friction properties in water lubrication due to the oxidation reaction on the surface in pre-sliding. However, the influence of oxidation reactions on low friction mechanism is still unclear. In this study, we clarify the structure change of DLC with the oxidation reaction in the pre-sliding using first-principles calculation, which suggests the low friction mechanism of DLC in water lubrication. The results show the structure change from sp3 carbon (Csp3) to sp2 carbon (Csp2) by the oxidation reaction on the surface. Furthermore, the Csp2 rich surface in water lubrication indicates the smooth sliding. We suggest that the structure change from Csp3 to Csp2 would affect low friction properties of DLC in water lubrication.

A theoretical study on sintering of Ni nanoparticles in the anode of solid oxide fuel cell under water vapor environment

Sintering of Ni particles in the Ni-based anode is a major obstacle to the widespread use of solid oxide fuel cell because the sintering induces the degradation in the anode. The large amount of water vapor in the fuel is known to accelerate the degradation during the operation. However, the detailed accelerated sintering mechanism is unclear. In this study, to clear the accelerated sintering mechanism during the initial stage of the sintering by water vapor, we investigated the adsorption and dissociation of water molecules on the Ni surface as well as the effect of the terminations of O, H, and OH on the interaction between the Ni clusters by density functional theory because the sintering of Ni particles is started by the Ni particles approaching each other due to the attractive interaction between Ni particles at the initial stage. In our adsorption and dissociation calculations, increasing the amount of water vapor facilitates the adsorption of H2O molecule on the Ni surface due to the H-bond interaction. Meanwhile, the activation energy for the dissociation of H2O molecules on the Ni surface is also lowered with increasing the amount of H2O molecules. Then, we calculated the interaction between Ni cluster, and O terminated, OH terminated, and H terminated Ni clusters. We observed that the attractive interaction between Ni cluster and O terminated Ni cluster is larger than that between two Ni clusters in vacuum. Enhancing the attractive interaction is not observed in the H terminated and OH terminated Ni clusters. It indicates that two Ni clusters approaching each other is faster under the water vapor environment than that in vacuum due to the strong attractive interaction caused by the termination of O. Thus, we suggest that the generation of O termination plays an important role in the sintering at the initial stage under the water vapor environment.