Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAtomistic-Object Kinetic Monte Carlo simulations of irradiation damage in tungsten
255
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe the development of a new object kinetic Monte Carlo code where the elementary defect objects are off-lattice atomistic configurations. Atomic-level transitions are used to transform and translate objects, to split objects and to merge them together. This gradually constructs a database of atomic configurations -- a set of relevant defect objects and their possible events generated on-the-fly. Elastic interactions are handled within objects with empirical potentials at short distances, and between spatially distinct objects using the dipole tensor formalism. The model is shown to evolve mobile interstitial clusters in tungsten faster than an equivalent molecular dynamics simulation, even at elevated temperatures. We apply the model to the evolution of complex defects generated using molecular dynamics simulations of primary radiation damage in tungsten. We show that we can evolve defect structures formed in cascade simulations to experimentally observable timescales of seconds while retaining atomistic detail. We conclude that the first few nanoseconds of simulation following cascade initiation would be better performed using molecular dynamics, as this will capture some of the near-temperature-independent evolution of small highly-mobile interstitial clusters. We also conclude that, for the 20keV PKA cascades annealing simulations considered here, internal relaxations of sessile objects difficult to capture using conventional object KMC with idealised object geometries establish the conditions for long timescale evolution.
rate research
Read More
High-temperature, high-dose, neutron irradiation of W results in the formation of Re-rich clusters at concentrations one order of magnitude lower than the thermodynamic solubility limit. These clusters may eventually transform into brittle W-Re intermetallic phases, which can lead to high levels of hardening and thermal conductivity losses. Standard theories of radiation enhanced diffusion and precipitation cannot explain the formation of these precipitates and so understanding the mechanism by which nonequilibrium clusters form under irradiation is crucial to predict materials degradation and devise mitigation strategies. Here we carry out a thermodynamic study of W-Re alloys and conduct kinetic Monte Carlo simulations of Re cluster formation in irradiated W-2Re alloys using a generalized Hamiltonian for crystals containing point defects parameterized entirely with electronic structure calculations. Our model incorporates recently-gained mechanistic information of mixed-interstitial solute transport, which is seen to control cluster nucleation and growth by forming quasi-spherical nuclei after an average incubation time of 20 s at 1800 K. These nuclei are seen to grow by attracting more mixed interstitials bringing solute atoms, which in turns attracts vacancies leading to recombination and solute agglomeration. The clusters grow to a maximum size of approximately 4-nm radius, and are not fully dense with Re, containing 50% or less near the center. Our simulations are in reasonable agreement with recent atom probe examinations of ion irradiated W-2Re systems at 773 K.
We present a new efficient method for Monte Carlo simulations of diffusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditional small diffusion hops and propagates the diffusing particles over long distances through a sequence of super-hops, one particle at a time. By partitioning the simulation space into non-overlapping protecting domains each containing only one or two particles, the algorithm factorizes the N-body problem of collisions among multiple Brownian particles into a set of much simpler single-body and two-body problems. Efficient propagation of particles inside their protective domains is enabled through the use of time-dependent Greens functions (propagators) obtained as solutions for the first-passage statistics of random walks. The resulting Monte Carlo algorithm is event-driven and asynchronous; each Brownian particle propagates inside its own protective domain and on its own time clock. The algorithm reproduces the statistics of the underlying Monte-Carlo model exactly. Extensive numerical examples demonstrate that for an important class of diffusion-reaction models the new algorithm is efficient at low particle densities, where other existing algorithms slow down severely.
As impermeable to gas molecules and at the same time transparent to high-energy ions, graphene has been suggested as a window material for separating a high-vacuum ion beam system from targets kept at ambient conditions. However, accumulation of irradiation-induced damage in the graphene membrane may give rise to its mechanical failure. Using atomistic simulations, we demonstrate that irradiated graphene even with a high vacancy concentration does not show signs of such instability, indicating a considerable robustness of graphene windows. We further show that upper and lower estimates for the irradiation damage in graphene can be set using a simple model.
Controlled anisotropic growth of two-dimensional materials provides an approach for the synthesis of large single crystals and nanoribbons, which are promising for applications as low-dimensional semiconductors and in next-generation optoelectronic devices. In particular, the anisotropic growth of transition metal dichalcogenides induced by the substrate is of great interest due to its operability. To date, however, their substrate-induced anisotropic growth is typically driven by the optimization of experimental parameters without uncovering the fundamental mechanism. Here, the anisotropic growth of monolayer tungsten disulfide on an ST-X quartz substrate is achieved by chemical vapor deposition, and the mechanism of substrate-induced anisotropic growth is examined by kinetic Monte Carlo simulations. These results show that, besides the variation of substrate adsorption, the chalcogen to metal (C/M) ratio is a major contributor to the large growth anisotropy and the polarization of undergrowth and overgrowth; either perfect isotropy or high anisotropy can be expected when the C/M ratio equals 2.0 by properly controlling the linear relationship between gas flux and temperature.
The pyramid-to-dome transition in Ge$_{x}$Si$_{1-x}$ on Si(100) initiated by step bunching on pyramidal quantum dots is atomistically simulated using a novel multi-state lattice model incorporating effective surface reconstructions. Results are explained by a simple theory based on a shallow island approximation. Under given deposition conditions in $d$ dimensions, the shape transition is shown to occur at island size $n_c$ following $n_c^{1/d} propto x^{-zeta}$ independent of temperature and deposition rate, where $zetaalt 2$ and $x$ is the actual Ge concentration in the island. The transition has an energy barrier dominated by the facet interface energy. Fast deposition however can out-run and delay the transition to larger island sizes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
