Search

On the path toward developing the first fusion energy devices, plasma simulations have become indispensable tools for supporting the design and development of fusion machines. Among these critical simulation tools, BIT1 is an advanced Particle-in-Cell code with Monte Carlo collisions, specifically designed for modeling plasma-material interaction and, in particular, analyzing the power load distribution on tokamak divertors. The current implementation of BIT1 relies exclusively on MPI for parallel communication and lacks support for GPUs. In this work, we address these limitations by designing and implementing a hybrid, shared-memory version of BIT1 capable of utilizing GPUs. For shared-memory parallelization, we rely on OpenMP and OpenACC, using a task-based approach to mitigate load-imbalance issues in the particle mover. On an HPE Cray EX computing node, we observe an initial performance improvement of approximately 42%, with scalable performance showing an enhancement of about 38% when using 8 MPI ranks. Still relying on OpenMP and OpenACC, we introduce the first version of BIT1 capable of using GPUs. We investigate two different data movement strategies: unified memory and explicit data movement. Overall, we report BIT1 data transfer findings during each PIC cycle. Among BIT1 GPU implementations, we demonstrate performance improvement through concurrent GPU utilization, especially when MPI ranks are assigned to dedicated GPUs. Finally, we analyze the performance of the first BIT1 GPU porting with the NVIDIA Nsight tools to further our understanding of BIT1 computational efficiency for large-scale plasma simulations, capable of exploiting current supercomputer infrastructures., Comment: Accepted by ICCS 2024 (The 24th International Conference on Computational Science), prepared in English, formatted according to the Springer LNCS templates and consists of 15 pages, which includes the main text, references, and figures

We revisit the topic of probability measures on CAT(0) cube complexes and prove that an amenable group acting on a CAT(0) cube complex, regardless of dimension, necessarily preserves an interval in the Roller compactification. In the finite dimensional case, we prove that there must be an orbit of cardinality $2^N$, where $N$ is bounded by the dimension. This is a slight extension of the author's previous Tits' Alternative., Comment: 11 pages, 11 figures, in english

Atmospheric muons are the most frequently observed form of cosmic radiation. Despite this, the existence of the muon flux component produced in decays of short-lived parent particles, called prompt muon flux still awaits experimental confirmation. This contribution to the muon flux is expected to start dominating at high energies, around PeV, since many of the prompt parent particles are heavy hadrons, containing charm and strange quarks. The aim of this thesis was two-fold: to evaluate the possibility of observing the prompt muon flux and to validate the performance of the KM3NeT detectors. [shortened abstract, the full version is contained in the thesis], Comment: PhD Thesis, 210 pages, Abstract in English and Polish, extensive appendix

Non-unitary $3 \times 3$ mixing matrices are present in many extensions of the Standard Model. One of the most simple extensions is the addition of $n_R$ right-handed neutrinos without imposing Lepton Number Conservation ($n_R$$\nu$SM). In the region of the parameter space where the usual seesaw approximation is valid, deviations from $3 \times 3$ unitarity of the leptonic mixing matrix are negligible, while one has three light neutrinos and two or three heavy neutrinos with masses close to the GUT scale. This gap between the mass scales occurs while having order 1 Yukawa couplings, providing a natural explanation for the smallness of light neutrino masses - hence the name seesaw. Other regions of the parameter space, with lighter heavy neutrinos - with an eV or keV mass, for instance - and bigger deviations from unitarity of the leptonic mixing matrix, are still allowed by experiment. In fact, whenever heavy neutrino masses are many orders of magnitude below the GUT scale there are interesting phenomenological implications at low energies. To analyse such regions with precision, one needs an exact parameterisation, since approximations that were valid in the usual seesaw case now fail. This thesis will consist of a study and classification of regions of the 3$\nu$SM parameter space in the light of a newly-developed exact parameterisation. Furthermore, it will also include a study of models with vector-like quarks, another simple extension of the SM, that contains a $3 \times 3$ non-unitary quark mixing matrix. Studying the most interesting regions of the parameter space of vector-like quarks models also benefits from the use of the aforementioned exact parameterisation. These regions are phenomenologically rich and may explain certain open questions like the CKM unitarity problem/Cabibbo angle anomaly. A chapter discussing the experimental implications of these models is also included., Comment: PhD Thesis, 186 pages, Abstract in English and Portuguese

If $G$ is a finite group, then the spectrum $\omega(G)$ is the set of all element orders of $G$. The prime spectrum $\pi(G)$ is the set of all primes belonging to $\omega(G)$. A simple graph $\Gamma(G)$ whose vertex set is $\pi(G)$ and in which two distinct vertices $r$ and $s$ are adjacent if and only if $rs \in \omega(G)$ is called the Gruenberg-Kegel graph or the prime graph of $G$. In this paper, we prove that if $G$ is a group of even order, then the set of vertices which are non-adjacent to $2$ in $\Gamma(G)$ form a union of cliques. Moreover, we decide when a strongly regular graph is isomorphic to the Gruenberg-Kegel graph of a finite group. Besides this, we prove that a complete bipartite graph with each part of size at least $3$ can not be isomorphic to the Gruenberg-Kegel graph of a finite group., Comment: The authors of this paper are ordered with respect to alphabet ordering in English

Many countries around the world, including Korea, use the school choice lottery system. However, this method has a problem in that many students are assigned to less-preferred schools based on the lottery results. In addition, the task of finding a good assignment with ties often has a time complexity of NP, making it a very difficult problem to improve the quality of the assignment. In this paper, we prove that the problem of finding a stable matching that maximizes the student-oriented preference utility in a two-sided market with one-sided preference can be solved in polynomial time, and we verify through experiments that the quality of assignment is improved. The main contributions of this paper are as follows. We found that stable student-oriented allocation in a two-sided market with one-sided preferences is the same as stable allocation in a two-sided market with symmetric preferences. In addition, we defined a method to quantify the quality of allocation from a preference utilitarian perspective. Based on the above two, it was proven that the problem of finding a stable match that maximizes the preference utility in a two-sided market with homogeneous preferences can be reduced to an allocation problem. In this paper, through an experiment, we quantitatively verified that optimal student assignment assigns more students to schools of higher preference, even in situations where many students are assigned to schools of low preference using the existing assignment method., Comment: 17 pages, 6 figures, in English; 15 pages, 6 figures