Researcher profile

Flyura Djurabekova

Flyura Djurabekova contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 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.