Source author record

Flyura Djurabekova

Flyura Djurabekova appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

Anisotropic Kinetics of Ion-Irradiation-Induced Phase Transition in Gallium Oxide

Radiation-tolerant semiconductors have traditionally been engineered by the principle of suppressing defect accumulation and amorphization, based on the assumption that radiation damage is inherently stochastic. Here we show that, in monoclinic $β$-\ce{Ga2O3}, a promising ultrawide-bandgap semiconductor, surface crystallographic orientation deterministically governs radiation tolerance through highly anisotropic kinetics of the $β$-to-$γ$ phase transition. Using machine-learning molecular dynamics coupled with a local configurational-entropy descriptor, we quantitatively map anisotropic $β$-to-$γ$ transition kinetics, showing that the critical dose, transition-layer depth, and kinetic stability of the $γ$-phase are fundamentally governed by surface orientation. Under ion irradiation, non-channeling surfaces such as (100), (001), and (-201) undergo severe surface amorphization, whereas the strongly channeling (010) surface resists damage accumulation and promotes subsurface $γ$-phase nucleation. During thermal annealing recovery process, these initial states follow two distinct recovery pathways: the channeling (010) surface reverts directly from $γ$-to-$β$, whereas non-channeling surfaces follow a sequential amorphous-to-$γ$-to-$β$ transition pathway. This work establishes surface orientation as a fundamental design principle for achieving radiation tolerance through controlled polymorphic transitions, providing a universal framework for engineering functional materials capable of withstanding extreme irradiation environments.

preprint2026arXiv

Insights Into Radiation Damage in YBa$_2$Cu$_3$O$_{7-δ}$ From Machine-Learned Interatomic Potentials

Accurate prediction of radiation damage in YBa$_2$Cu$3$O${7-δ}$ (YBCO) is essential for assessing the performance of high-temperature superconducting (HTS) tapes in compact fusion reactors. Existing empirical interatomic potentials have been used to model radiation damage in stoichiometric YBCO, but fail to describe oxygen-deficient compositions, which are ubiquitous in industrial Rare-Earth Barium Copper Oxide conductors and strongly influence superconducting properties. In this work, we demonstrate that modern machine-learned interatomic potentials enable predictive modelling of radiation damage in YBCO across a wide range of oxygen stoichiometries, with higher fidelity than previous empirical models. We employ two recently developed approaches: an Atomic Cluster Expansion (ACE) potential and a tabulated Gaussian Approximation Potential (tabGAP). Both models accurately reproduce Density Functional Theory (DFT) energies, forces, and threshold displacement energy distributions, providing a reliable description of atomic-scale collision processes. Molecular dynamics simulations of 5 keV cascades predict enhanced peak defect production and recombination relative to a widely used empirical potential, indicating different cascade evolution. By explicitly varying oxygen deficiency, we show that total defect production depends only weakly on stoichiometry, offering insight into the robustness of radiation damage processes in oxygen-deficient YBCO. Finally, fusion-relevant 300 keV cascade simulations reveal amorphous regions with dimensions comparable to the superconducting coherence length, consistent with electron microscopy observations of neutron-irradiated HTS tapes. These results establish machine-learned interatomic potentials as efficient and predictive tools for investigating radiation damage in YBCO across relevant compositions and irradiation conditions.

preprint2022arXiv

Insights into nanoparticle shape transformation by energetic ions using atomistic simulations

The shape of metal nanoparticles embedded in dielectric matrices influences the optical properties of the composite material. Swift heavy ion irradiation can induce a controllable shape transformation in gold and other metals embedded in amorphous silicon dioxide, where the particles elongate along the direction of the ion beam. The details of this transformation are not fully understood, but it is presumably related to nanometer-scale phase transitions induced by individual ion impacts. The phenomenon has been reproduced using atomistic simulations, although the time scale limitations and the lack of accurate interatomic models within the metal-silica interface lead unavoidably to severe simplifications. We improve the realism in the simulations with an accurate model for surface adhesion between gold and silica and by simulating the processes in the matrix between impacts. The simulations with correct adhesion show that the nanoparticles can grow in aspect ratio in the molten state even after silicon dioxide solidifies. Moreover, we demonstrate the active role of the matrix: without explicitly modeling processes in the matrix between impacts, the elongation is limited and does not reach significant aspect ratios seen in experiments. These results significantly improve the theoretical understanding of processes developing in embedded nanoparticles under swift heavy ion irradiation. The knowledge brings forward the ion beam technology as a precise tool for shaping of embedded nanostructures for various optical applications.

preprint2022arXiv

Multiscale machine-learning interatomic potentials for ferromagnetic and liquid iron

We develop and compare four interatomic potentials for iron: a simple machine-learned embedded atom method (EAM) potential, a potential with machine-learned two- and three-body-dependent terms, a potential with machine-learned EAM and three-body terms, and a Gaussian approximation potential with the SOAP descriptor. All potentials are trained to the same diverse database of body-centered cubic and liquid structures computed with density functional theory. The four presented potentials represent different levels of complexity and span three orders of magnitude in computational cost. The first three potentials are tabulated and evaluated efficiently using cubic spline interpolations, while the fourth one is implemented without additional optimization. We compare and discuss the advantages of each implementation, transferability and applicability in terms of the balance between required accuracy versus computational cost.

preprint2022arXiv

Polarization characteristics of adatoms self-diffusing on metal surfaces under high electric fields

Although atomic diffusion on metal surfaces under high electric fields has been studied theoretically and experimentally since the 1970s, its accurate and quantitative theoretical description remains a significant challenge. In our previous work, we developed a theoretical framework that describes the atomic dynamics on metal surfaces in the presence of an electric field in terms of the local polarization characteristics of the surface at the vicinity of a moving atom. Here, we give a deeper analysis of the physics underlying this framework, introducing and rigorously defining the concept of the effective polarization characteristics (permanent dipole moment $μ$ and polarizability $α$) of a moving atom on a metal surface, which are shown to be the relevant atomic quantities determining the dynamics of a moving atom via a compact equation. We use density functional theory (DFT) to calculate $μ$ and $α$ of a W adatom moving on a W {110} surface, where additional adatoms are present in its vicinity. We analyze the dependence of $μ$ and $α$ and hence the migration barriers under electric fields on the local atomic environments (LAE) of an adatom. We find that the LAE significantly affects $μ$ and $α$ of a moving atom in the limited cases we studied, which implies that further systematic DFT calculations are needed to fully parameterize surface diffusion in terms of energy barriers for long-term large scale simulations, such as our recently developed Kinetic Monte Carlo model for surface diffusion under electric field.

preprint2022arXiv

Simple machine-learned interatomic potentials for complex alloys

Developing data-driven machine-learning interatomic potentials for materials containing many elements becomes increasingly challenging due to the vast configuration space that must be sampled by the training data. We study the learning rates and achievable accuracy of machine-learning interatomic potentials for many-element alloys with different combinations of descriptors for the local atomic environments. We show that for a five-element alloy system, potentials using simple low-dimensional descriptors can reach meV/atom-accuracy with modestly sized training datasets, significantly outperforming the high-dimensional SOAP descriptor in data efficiency, accuracy, and speed. In particular, we develop a computationally fast machine-learned and tabulated Gaussian approximation potential (tabGAP) for Mo-Nb-Ta-V-W alloys with a combination of two-body, three-body, and a new simple scalar many-body density descriptor based on the embedded atom method.

preprint2020arXiv

Classification of vacuum arc breakdowns in a pulsed DC system

Understanding the microscopic phenomena behind vacuum arc ignition and generation is crucial for being able to control the breakdown rate, thus improving the effectiveness of many high-voltage applications where frequent breakdowns limit the operation. In this work, statistical properties of various aspects of breakdown, such as the number of pulses between breakdowns, breakdown locations and crater sizes are studied independently with almost identical Pulsed DC Systems at the University of Helsinki and in CERN. In high-gradient experiments, copper electrodes with parallel plate capacitor geometry, undergo thousands of breakdowns. The results support the classification of the events into primary and secondary breakdowns, based on the distance and number of pulses between two breakdowns. Primary events follow a power law on the log--log scale with the slope $α\approx 1.33$, while the secondaries are highly dependent on the pulsing parameters.

preprint2020arXiv

Crystal Defects: A Portal To Dark Matter Detection

We propose to use the defect creation energy loss in commonly used high energy physics solid state detectors as a tool to statistically identify dark matter signal from background. We simulate the energy loss in the process of defect creation using density functional theory and molecular dynamics methods and calculate the corresponding expected dark matter spectra. We show that in phonon-mediated solid state detectors, the energy loss due to defect creation convolved with the expected dark matter interaction signal results in a significant change in the expected spectra for common detector materials. With recent progress towards $\sim$10 eV threshold low-mass dark matter searches, this variation in expected dark matter spectrum can be used as a direct signature of dark matter interactions with atomic nuclei.

preprint2020arXiv

Dynamic coupling between particle-in-cell and atomistic simulations

We propose a method to directly couple molecular dynamics, finite element method and particle-in-cell techniques to simulate metal surface response to high electric fields. We use this method to simulate the evolution of a field emitting tip under thermal runaway by fully including the 3D space-charge effects. We also present a comparison of the runaway process between the two tip geometries of different widths. The results show with high statistical significance that in case of sufficiently narrow field emitters, the thermal runaway occurs in cycles where intensive neutral evaporation alternates with the cooling periods. The comparison with previous works shows, that the evaporation rate in the regime of intensive evaporation is sufficient to ignite a plasma arc above the simulated field emitters.

preprint2020arXiv

Fundamental Phenomena and Applications of Swift Heavy Ion Irradiations

This review concentrates on the specific properties and characteristics of damage structures generated with high-energy ions in the electronic energy loss regime. Irradiation experiments with so-called swift heavy ions (SHI) find applications in many different fields, with examples presented in ion-track nanotechnology, radiation hardness analysis of functional materials, and laboratory tests of cosmic radiation. The basics of the SHI-solid interaction are described with special attention to processes in the electronic subsystem. The broad spectrum of damage phenomena is exemplified for various materials and material classes, along with a description of typical characterization techniques. The review also presents state-of-the-art modeling efforts that try to account for the complexity of the coupled processes of the electronic and atomic subsystems. Finally, the relevance of SHI phenomena for effects induced by fission fragments in nuclear fuels and how this knowledge can be applied to better estimate damage risks in nuclear materials is discussed.

preprint2020arXiv

Growth mechanism for nanotips in high electric fields

In this work we show using atomistic simulations that the biased diffusion in high electric field gradients creates a mechanism whereby nanotips may start growing from small surface asperities. It has long been known that atoms on a metallic surface have biased diffusion if electric fields are applied and that microscopic tips may be sharpened using fields, but the exact mechanisms have not been well understood. Our Kinetic Monte Carlo simulation model uses a recently developed theory for how the migration barriers are affected by the presence of an electric field. All parameters of the model are physically motivated and no fitting parameters are used. The model has been validated by reproducing characteristic faceting patterns of tungsten surfaces that have in previous experiments been observed to only appear in the presence of strong electric fields. The growth effect is found to be enhanced by increasing fields and temperatures.

preprint2020arXiv

Tungsten migration energy barriers for surface diffusion: a parameterization for KMC simulations

We have calculated the migration barriers for surface diffusion on Tungsten. Our results form a self-sufficient parameterization for Kinetic Monte Carlo simulations of arbitrarily rough atomic tungsten surfaces, as well as nanostructures such as nanotips and nanoclusters. The parameterization includes first- and second-nearest neighbour atom jump processes, as well as a third-nearest neighbour exchange process. The migration energy barriers of all processes are calculated with the Nudged Elastic Band method. The same attempt frequency for all processes is found sufficient and the value is fitted to Molecular Dynamics simulations. The model is validated by correctly simulating with Kinetic Monte Carlo the energetically favourable W nanocluster shapes, in good agreement with Molecular Dynamics simulations.

preprint2019arXiv

Vacuum electrical breakdown conditioning study in a parallel plate electrode pulsed DC system

Conditioning of a metal structure in a high-voltage system is the progressive development of resistance to vacuum arcing over the operational life of the system. This is, for instance, seen during the initial operation of radio frequency (rf) cavities in particle accelerators. It is a relevant topic for any technology where breakdown limits performance, and where conditioning continues for a significant duration of system runtime. Projected future linear accelerators require structures with accelerating gradients of up to 100 MV/m. Currently, this performance level is only achievable after a multi-month conditioning period. In this work, a pulsed DC system applying voltage pulses over parallel disk electrodes was used to study the conditioning process, with the objective of obtaining insight into its underlying mechanics, and ultimately, to find ways to shorten the conditioning process. Two kinds of copper electrodes were tested: As-prepared machine-turned electrodes ("hard" copper), and electrodes that additionally had been subjected to high temperature treatments ("soft" copper). The conditioning behaviour of the soft electrodes was found to be similar to that of comparably treated accelerating structures, indicating a similar conditioning process. The hard electrodes reached the same ultimate performance as the soft electrodes much faster, with a difference of more than an order of magnitude in the number of applied voltage pulses. Two distinctly different distributions of breakdown locations were observed on the two types of electrodes. Considered together, our results support the crystal structure dislocation theory of breakdown, and suggest that the conditioning of copper in high field systems such as rf accelerating structures is dominated by material hardening.

preprint2016arXiv

A general computational method for electron emission and thermal effects in field emitting nanotips

Electron emission from nanometric size emitters becomes of increasing interest due to its involvement to sharp electron sources, vacuum breakdown phenomena and various other vacuum nanoelectronics applications. The most commonly used theoretical tools for the calculation of electron emission are still nowadays the Fowler-Nordheim and the Richardson-Laue-Dushman equations although it has been shown since the 1990's that they are inadequate for nanometrically sharp emitters or in the intermediate thermal-field regime. In this paper we develop a computational method for the calculation of emission currents and Nottingham heat, which automatically distinguishes among different emission regimes, and implements the appropriate calculation method for each. Our method covers all electron emission regimes (thermal, field and intermediate), aiming to maximize the calculation accuracy while minimizing the computational time. As an example, we implemented it in atomistic simulations of the thermal evolution of Cu nanotips under strong electric fields and found that the predicted behaviour of such nanotips by the developed technique differs significantly from estimations obtained based on the Fowler-Nordheim equation. Finally, we show that our tool can be also successfully applied in the analysis of experimental I-V data.

preprint2016arXiv

Electrodynamics - molecular dynamics simulations of the stability of Cu nanotips under high electric field

The shape memory effect and pseudoelasticity in Cu nanowires is one possible pair of mechanisms that prevents high aspect ratio nanosized field electron emitters to be stable at room temperature and permits their growth under high electric field. By utilizing hybrid electrodynamics molecular dynamics simulations we show that a global electric field of 1 GV/m or more significantly increases the stability and critical temperature of spontaneous reorientation of nanosized <100> Cu field emitters. We also show that in the studied tips the stabilizing effect of an external applied electric field is an order of magnitude greater than the destabilization caused by the field emission current. We detect the critical temperature of spontaneous reorientation using the tool that spots the changes in crystal structure. The method is compatible with techniques that consider the change in potential energy, has a wider range of applicability and allows pinpointing different stages in the reorientation processes.

preprint2016arXiv

Laser-induced asymmetric faceting and growth of nano-protrusion on a tungsten tip

Irradiation of a sharp tungsten tip by a femtosecond laser and exposed to a strong DC electric field led to gradual and reproducible surface modifications. By a combination of field emission microscopy and scanning electron microscopy, we observed asymmetric surface faceting with sub-ten nanometer high steps. The presence of well pronounced faceted features mainly on the laser-exposed side implies that the surface modification was driven by a laser-induced transient temperature rise -- on a scale of a couple of picoseconds -- in the tungsten tip apex. Moreover, we identified the formation of a nano-tip a few nanometers high located at one of the corners of a faceted plateau. The results of simulations emulating the experimental conditions, are consistent with the experimental observations. The presented conditions can be used as a new method to fabricate nano-tips of few nm height, which can be used in coherent electron pulses generation. Besides the direct practical application, the results also provide insight into the microscopic mechanisms of light-matter interaction. The apparent growth mechanism of the features may also help to explain the origin of enhanced electron field emission, which leads to vacuum arcs, in high electric-field devices such as radio-frequency particle accelerators.

preprint2016arXiv

Long-term stability of Cu surface nanotips

Sharp nanoscale tips on metal surfaces of electrodes enhance locally applied electric fields. Strongly enhanced electric fields trigger electron field emission and atom evaporation from the apexes of the nanotips. Combined together, these processes may explain electric discharges in form of small local arcs observed near metal surfaces in the presence of electric fields even in ultra high vacuum conditions. In the present work we investigate the stability of nanoscale tips by means of computer simulations of surface diffusion processes on copper (Cu), the main material of high voltage electronics. We study the stability and life-time of thin Cu surface nanotips at different temperatures in terms of diffusion processes. For this purpose, we have developed a surface Kinetic Monte Carlo (KMC) model where the jump processes are described by tabulated precalculated energy barriers. We show that tall surface features with high aspect ratios can be fairly stable at room temperature. However, the stability was found to depend strongly on the temperature: 13 nm nanotips with the major axes in the <110> crystallographic directions were found to flatten down to half of the original height in less than 100 ns at temperatures close to the melting point, whereas no significant change in the height of these nanotips was observed after 10 $μ$s at room temperature. Moreover, the nanotips built up along the <110> crystallographic directions were found significantly more stable than those oriented in the <100> or <111> crystallographic directions. The proposed KMC model has been found well suited for simulating atomic surfaces processes and was validated against Molecular Dynamics (MD) simulation results via the comparison of the flattening times obtained by both methods. We also note that the KMC simulations were two orders of magnitude computationally faster than the corresponding MD calculations.