Search

We have performed time-dependent density functional theory simulations for the neutralization process of a hydrogen ion injected at 100eV onto the (110) surface of tungsten material. We have also proposed a method for evaluating the detection probability of electrons in a small region. This probability is interpreted as that of detecting hydrogen in each state: positive ion, neutral atom, and negative ion. As a result, the probabilities of detecting hydrogen after a collision as a positive ion, neutral atom, and negative ion were approximately 30 percent, 50 percent, and 20 percent, respectively., Comment: 10 pages, 6 figures, 1 table, 26-th International Conference on Plasma Surface Interaction

The neutralization process of an ion particle on a surface is key issue of the plasma-wall interaction (PWI). We investigated helium (He) ion injection onto a tungsten surface using time-dependent density functional theory (TDDFT) simulation. We developed the TDDFT code QUMASUN for this study, simulating the process of electron transfer from the surface to the He nucleus by solving the time evolution of electron wave function and the classical motion of nuclei simultaneously. Our results show that the probabilities of injected $\text{He}^{2+}$ changing into $\text{He}^{1+}$ and $\text{He}^{0}$ on the surface are approximately 40 percent and 13 percent, respectively. The captured electrons by $\text{He}^{1+}$ and $\text{He}^{0}$ predominantly occupy the 2s and 2p orbitals, corresponding to excited states. In addition, we stated some challenges for applying TDDFT to plasma-wall interactions., Comment: 4 figures, 1 table, resented in 26th International conference on plasma-surface interaction

The total mass estimate of molecular clouds suffers from the uncertainty in the H$_2$-CO conversion factor, the so-called $X_{\rm CO}$ factor, which is used to convert the $^{12}$CO (1--0) integrated intensity to the H$_2$ column density. We demonstrate the machine learning's ability to predict the H$_2$ column density from the $^{12}$CO, $^{13}$CO, and C$^{18}$O (1--0) data set of four star-forming molecular clouds; Orion A, Orion B, Aquila, and M17. When the training is performed on a subset of each cloud, the overall distribution of the predicted column density is consistent with that of the Herschel column density. The total column density predicted and observed is consistent within 10\%, suggesting that the machine learning prediction provides a reasonable total mass estimate of each cloud. However, the distribution of the column density for values $> \sim 2 \times 10^{22}$ cm$^{-2}$, which corresponds to the dense gas, could not be predicted well. This indicates that molecular line observations tracing the dense gas are required for the training. We also found a significant difference between the predicted and observed column density when we created the model after training the data on different clouds. This highlights the presence of different $X_{\rm CO}$ factors between the clouds, and further training in various clouds is required to correct for these variations. We also demonstrated that this method could predict the column density toward the area not observed by Herschel if the molecular line and column density maps are available for the small portion, and the molecular line data are available for the larger areas., Comment: Accepted for publication in MNRAS

By using molecular dynamics simulation, formation mechanisms of amorphous carbon in particular sp${}^3$ rich structure was researched. The problem that reactive empirical bond order potential cannot represent amorphous carbon properly was cleared in the transition process from graphite to diamond by high pressure and the deposition process of amorphous carbon thin films. Moreover, the new potential model which is based on electron distribution simplified as a point charge was developed by using downfolding method. As a result, the molecular dynamics simulation with the new potential could demonstrate the transition from graphite to diamond at the pressure of 15 GPa corresponding to experiment and the deposition of sp${}^3$ rich amorphous carbon., Comment: International Symposium ISPlasma 2012, submitted to Jpn. J. Appl. Phys

By molecular dynamics simulation, the chemical vapor deposition of amorphous carbon onto graphite and diamond surfaces was studied. In particular, we investigated the effect of source H/C ratio, which is the ratio of the number of hydrogen atoms to the number of carbon atoms in a source gas, on the deposition process. In the present simulation, the following two source gas conditions were tested: one was that the source gas was injected as isolated carbon and hydrogen atoms, and the other was that the source gas was injected as hydrocarbon molecules. Under the former condition, we found that as the source H/C ratio increases, the deposition rate of carbon atoms decreases exponentially. This exponential decrease in the deposition rate with increasing source H/C ratio agrees with experimental data. However, under the latter molecular source condition, the deposition rate did not decrease exponentially because of a chemical reaction peculiar to the type of hydrocarbon in the source gas., Comment: accepted by Jpn. J. Appl. Phys. (2008)

Molecular dynamics (MD) simulation with modified Brenner's reactive empirical bond order (REBO) potential is a powerful tool to investigate plasma wall interaction on divertor plates in a nuclear fusion device. However, MD simulation box's size is less than several nm for the performance of a computer. To extend the size of the MD simulation, we develop a hybrid simulation code between MD code using REBO potential and binary collision approximation (BCA) code. Using the BCA code instead of computing all particles with a high kinetic energy for every step in the MD simulation, considerable computation time is saved. By demonstrating a hydrogen atom injection on a graphite by the hybrid simulation code, it is found that the hybrid simulation code works efficiently in a large simulation box., Comment: 5 pages, 5 figures

We have investigated plasma-surface interactions with molecular dynamics (MD) simulations. It, however, is high cost computation and is limited to simulations for materials of nanometer order. In order to overcome the limitation, a complementary model based on binary collision approximation (BCA) can be established. We employed a BCA-based simulation code ACAT and extended to handle any structure involving crystalline and amorphous. The extended code, named "ACaT", stores all positions of projectile and target atoms and velocities of recoil atoms, so it can be combined with the MD code. It also holds the potential to reproduce channeling phenomena. Thus it is expected to be useful for evaluation of channeling effects.