Sputtering of Graphite by Hydrogen Isotopes in the Fusion Environment: A Molecular Dynamics Simulation Study

2018 ◽  
Vol 4 (4) ◽  
Author(s):  
Qiang Zhao ◽  
Yang Li ◽  
Zheng Zhang ◽  
Xiaoping Ouyang

The sputtering of graphite due to the bombardment of hydrogen isotopes is crucial to successfully using graphite in the fusion environment. In this work, we use molecular dynamics to simulate the sputtering using the large-scale atomic/molecular massively parallel simulator (lammps). The calculation results show that the peak values of the sputtering yield are between 25 eV and 50 eV. When the incident energy is greater than the energy corresponding to the peak value, a lower carbon sputtering yield is obtained. The temperature that is most likely to sputter is approximately 800 K for hydrogen, deuterium, and tritium. Below the 800 K, the sputtering yields increase with temperature. By contrast, above the 800 K, the yields decrease with increasing temperature. Under the same temperature and incident energy, the sputtering rate of tritium is greater than that of deuterium, which in turn is greater than that of hydrogen. When the incident energy is 25 eV, the sputtering yield at 300 K increases below an incident angle at 30 deg and remains steady after that.

Author(s):  
Qiang Zhao ◽  
Yang Li ◽  
Zheng Zhang ◽  
Xiaoping Ouyang

The sputtering of graphite due to the bombardment of hydrogen isotopes is one of the critical issues in successfully using graphite in the fusion environment. In this work, we use molecular dynamics method to simulate the sputtering by using the LAMMPS. Calculation results show that the peak values of the sputtering yield are located between 25 eV to 50 eV. After the energy of 25 eV, the higher incident energy cause the lower carbon sputtering yield. The temperature which is most likely to sputter is about 800 K for hydrogen, deuterium and tritium. Before the 800 K, the sputtering rates increase when the temperature increase. After the 800 K, they decrease with the temperature increase. Under the same temperature and energy, the sputtering rate of tritium is bigger than that of deuterium, the sputtering rate of deuterium is bigger than that of hydrogen.


2021 ◽  
Vol 11 (20) ◽  
pp. 9702
Author(s):  
Cheng-Hsuan Ho ◽  
Cha’o-Kuang Chen ◽  
Chieh-Li Chen

Interconnects are significant elements in integrated circuits (ICs), as they connect individual components of the circuit into a functioning whole. To form a void-free interconnect, a thin and uniform copper seed layer must be deposited as a basis for electroplating. In this paper, process parameters of sputtering including incident energy, incident angle, substrate temperature, and deposition rate were studied to form a uniform copper seed layer. Different liner/barrier materials and properties including crystal planes were also studied to enhance the quality of the copper seed layer. The study was carried out by molecular dynamics simulation. It revealed that increasing the incident energy and substrate temperature during the sputtering process increases their diffusivity but results in poorer uniformity and larger alloy percentage. By decreasing the deposition rate, the Ostwald ripening effect becomes dominant and increases the uniformity. An adequate incident angle could increase necking and uniformity. Among the sputtering process parameters and material properties discussed in this study, surface diffusion barrier energy of different crystal planes is the most decisive factor, which leads to good uniformity.


2019 ◽  
Author(s):  
Liqun Cao ◽  
Jinzhe Zeng ◽  
Mingyuan Xu ◽  
Chih-Hao Chin ◽  
Tong Zhu ◽  
...  

Combustion is a kind of important reaction that affects people's daily lives and the development of aerospace. Exploring the reaction mechanism contributes to the understanding of combustion and the more efficient use of fuels. Ab initio quantum mechanical (QM) calculation is precise but limited by its computational time for large-scale systems. In order to carry out reactive molecular dynamics (MD) simulation for combustion accurately and quickly, we develop the MFCC-combustion method in this study, which calculates the interaction between atoms using QM method at the level of MN15/6-31G(d). Each molecule in systems is treated as a fragment, and when the distance between any two atoms in different molecules is greater than 3.5 Å, a new fragment involved two molecules is produced in order to consider the two-body interaction. The deviations of MFCC-combustion from full system calculations are within a few kcal/mol, and the result clearly shows that the calculated energies of the different systems using MFCC-combustion are close to converging after the distance thresholds are larger than 3.5 Å for the two-body QM interactions. The methane combustion was studied with the MFCC-combustion method to explore the combustion mechanism of the methane-oxygen system.


2014 ◽  
Vol 1700 ◽  
pp. 61-66
Author(s):  
Guttormur Arnar Ingvason ◽  
Virginie Rollin

ABSTRACTAdding single walled carbon nanotubes (SWCNT) to a polymer matrix can improve the delamination properties of the composite. Due to the complexity of polymer molecules and the curing process, few 3-D Molecular Dynamics (MD) simulations of a polymer-SWCNT composite have been run. Our model runs on the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), with a COMPASS (Condensed phase Optimized Molecular Potential for Atomistic Simulations Studies) potential. This potential includes non-bonded interactions, as well as bonds, angles and dihedrals to create a MD model for a SWCNT and EPON 862/DETDA (Diethyltoluenediamine) polymer matrix. Two simulations were performed in order to test the implementation of the COMPASS parameters. The first one was a tensile test on a SWCNT, leading to a Young’s modulus of 1.4 TPa at 300K. The second one was a pull-out test of a SWCNT from an originally uncured EPON 862/DETDA matrix.


2021 ◽  
Author(s):  
Kai Xu ◽  
Lei Yan ◽  
Bingran You

Force field is a central requirement in molecular dynamics (MD) simulation for accurate description of the potential energy landscape and the time evolution of individual atomic motions. Most energy models are limited by a fundamental tradeoff between accuracy and speed. Although ab initio MD based on density functional theory (DFT) has high accuracy, its high computational cost prevents its use for large-scale and long-timescale simulations. Here, we use Bayesian active learning to construct a Gaussian process model of interatomic forces to describe Pt deposited on Ag(111). An accurate model is obtained within one day of wall time after selecting only 126 atomic environments based on two- and three-body interactions, providing mean absolute errors of 52 and 142 meV/Å for Ag and Pt, respectively. Our work highlights automated and minimalistic training of machine-learning force fields with high fidelity to DFT, which would enable large-scale and long-timescale simulations of alloy surfaces at first-principles accuracy.


Metals ◽  
2018 ◽  
Vol 8 (10) ◽  
pp. 770 ◽  
Author(s):  
Qinghua Chang ◽  
Jingpei Xie ◽  
Aixia Mao ◽  
Wenyan Wang

Large scale Atomic/Molecular dynamic Parallel Simulator (LAMMPS) molecular dynamics simulation software was used to simulate the copper and aluminum atom diffusion and changes of interface during heating and cooling process of copper and aluminum composite panels. The structures of the interface were characterized through scanning electron microscope (SEM), X-ray diffraction (XRD), and transmission electron microscope (TEM), and the mechanical properties were also tested. The simulation results show that the diffusion rate of copper atom is higher than that of aluminum atom, and that the CuAl2 radial distribution function of the interface at 300 K is consistent with that of pure CuAl2 at room temperature. At 930 K, t = 50 ps Cu atoms spread at a distance of approximately four Al lattice constants around the Al layer, and Al atoms spread to about half a lattice constant distance to the Cu layer. The experimental results show that the thickness of the interface in copper–aluminum composite plate is about 1 μm, and only one kind of CuAl2 with tetragonal phase structure is generated in the interface, which corresponds with the result of molecular dynamics simulation.


2021 ◽  
Vol 12 (6) ◽  
pp. 7239-7248

The novel coronavirus, recognized as COVID-19, is the cause of an infection outbreak in December 2019. The effect of temperature and pH changes on the main protease of SARS-CoV-2 were investigated using all-atom molecular dynamics simulation. The obtained results from the root mean square deviation (RMSD) and root mean square fluctuations (RMSF) analyses showed that at a constant temperature of 25℃ and pH=5, the conformational change of the main protease is more significant than that of pH=6 and 7. Also, by increasing temperature from 25℃ to 55℃ at constant pH=7, a remarkable change in protein structure was observed. The radial probability of water molecules around the main protease was decreased by increasing temperature and decreasing pH. The weakening of the binding energy between the main protease and water molecules due to the increasing temperature and decreasing pH has reduced the number of hydrogen bonds between the main protease and water molecules. Finding conditions that alter the conformation of the main protease could be fundamental because this change could affect the virus’s functionality and its ability to impose illness.


2016 ◽  
Vol 18 (26) ◽  
pp. 17461-17469 ◽  
Author(s):  
Z. Y. Hou ◽  
K. J. Dong ◽  
Z. A. Tian ◽  
R. S. Liu ◽  
Z. Wang ◽  
...  

The effect of the cooling rate on the solidification process of liquid aluminium is studied using a large-scale molecular dynamics method.


Sign in / Sign up

Export Citation Format

Share Document