Researcher profile

Ermal Rrapaj

Ermal Rrapaj contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
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

12 published item(s)

preprint2024arXiv

Correlations and Distinguishability Challenges in Supernova Models: Insights from Future Neutrino Detectors

This paper explores core-collapse supernovae as crucial targets for neutrino telescopes, addressing uncertainties in their simulation results. We comprehensively analyze eighteen modern simulations and discriminate among supernova models using realistic detectors and interactions. A significant correlation between the total neutrino energy and cumulative counts, driven by massive lepton neutrinos and oscillations, is identified, particularly noticeable with the DUNE detector. Bayesian techniques indicate strong potential for model differentiation during a Galactic supernova event, with HK excelling in distinguishing models based on equation of state, progenitor mass, and mixing scheme.

preprint2022arXiv

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales $t\proptoμ^{-1}$, with $μ$ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number $N$ of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales $τ\approxμ^{-1}\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales $τ\approxμ^{-1}\log(N)$. We present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

preprint2022arXiv

Improving Schrödinger Equation Implementations with Gray Code for Adiabatic Quantum Computers

We reformulate the continuous space Schrödinger equation in terms of spin Hamiltonians. For the kinetic energy operator, the critical concept facilitating the reduction in model complexity is the idea of position encoding. Binary encoding of position produces a Heisenberg-like model and yields exponential improvement in space complexity when compared to classical computing. Encoding with a binary reflected Gray code, and a Hamming distance 2 Gray code yields the additional effect of reducing the spin model down to the XZ and transverse Ising model respectively. We also identify the bijective mapping between diagonal unitaries and the Walsh series, producing the mapping of any real potential to a series of $k$-local Ising models through the fast Walsh transform. Finally, in a finite volume, we provide some numerical evidence to support the claim that the total time needed for adiabatic evolution is protected by the infrared cutoff of the system. As a result, initial state preparation from a free-field wavefunction to an interacting system is expected to exhibit polynomial time complexity with volume and constant scaling with respect to lattice discretization for all encodings. For the Hamming distance 2 Gray code, the evolution starts with the transverse Hamiltonian before introducing penalties such that the low lying spectrum reproduces the energy levels of the Laplacian. The adiabatic evolution of the penalty Hamiltonian is therefore sensitive to the ultraviolet scale. It is expected to exhibit polynomial time complexity with lattice discretization, or exponential time complexity with respect to the number of qubits given a fixed volume.

preprint2022arXiv

Sphaleron Transition Rates and the Chiral Magnetic Effect

The chiral magnetic effect is a novel quantum phenomenon proposed for high-energy nuclear collisions but which has yet to be observed. We quantify the axial charge relaxation time, due to sphalerons, which enters in simulations of this effect. An extrapolation of weak coupling calculations of the sphaleron rate yields rather different relaxation times than strong coupling AdS/CFT calculations. The AdS/CFT relaxation time is the larger one of the two by an order of magnitude, but the weak coupling relaxation time may not be reliable because it is only marginally bigger than the microscopic thermalization time. The role of quark masses has yet to be accurately assessed.

preprint2021arXiv

Exact representations of many body interactions with RBM neural networks

Restricted Boltzmann Machines (RBM) are simple statistical models defined on a bipartite graph which have been successfully used in studying more complicated many-body systems, both classical and quantum. In this work, we exploit the representation power of RBMs to provide an exact decomposition of many-body contact interactions into one-body operators coupled to discrete auxiliary fields. This construction generalizes the well known Hirsch's transform used for the Hubbard model to more complicated theories such as Pionless EFT in nuclear physics, which we analyze in detail. We also discuss possible applications of our mapping for quantum annealing applications and conclude with some implications for RBM parameter optimization through machine learning.

preprint2021arXiv

Lossy compression of statistical data using quantum annealer

We present a new lossy compression algorithm for statistical floating-point data through a representation learning with binary variables. The algorithm finds a set of basis vectors and their binary coefficients that precisely reconstruct the original data. The optimization for the basis vectors is performed classically, while binary coefficients are retrieved through both simulated and quantum annealing for comparison. A bias correction procedure is also presented to estimate and eliminate the error and bias introduced from the inexact reconstruction of the lossy compression for statistical data analyses. The compression algorithm is demonstrated on two different datasets of lattice quantum chromodynamics simulations. The results obtained using simulated annealing show 3.5 times better compression performance than the algorithms based on a neural-network autoencoder and principal component analysis. Calculations using quantum annealing also show promising results, but performance is limited by the integrated control error of the quantum processing unit, which yields large uncertainties in the biases and coupling parameters. Hardware comparison is further studied between the previous generation D-Wave 2000Q and the current D-Wave Advantage system. Our study shows that the Advantage system is more likely to obtain low-energy solutions for the problems than the 2000Q.

preprint2020arXiv

Exact solution of multi-angle quantum many-body collective neutrino flavor oscillations

I study the flavor evolution of a dense neutrino gas by considering vacuum contributions, matter effects and neutrino self-interactions. Assuming a system of two flavors in a uniform matter background, the time evolution of the many-body system in discretized momentum space is computed. The multi-angle neutrino-neutrino interactions are treated exactly and compared to both the single-angle approximation and mean field calculations. %The time unit chosen is $μ_0^{-1}=(\frac{G_F}{2\sqrt{2}V})^{-1}$. The mono-energetic two neutrino beam scenario is solved analytically. I proceed to solve flavor oscillations for mono-energetic cubic lattices and quadratic lattices of two energy levels. In addition I study various configurations of twelve, sixteen, and twenty neutrinos. I find that when all neutrinos are initially of the same flavor, all methods agree. When both flavors are present, I find collective oscillations and flavor equilibration develop in the many body treatment but not in the mean field method. This difference persists in dense matter with tiny mixing angle and it can be ascribed to non-negligible flavor polarization correlations being present. Entanglement entropy is significant in all such cases. The relevance for supernovae or neutron stars mergers is contingent upon the value of the normalization volume $V$ and the large $N$ dependence of the timescale associated with oscillations. In future work, I intend to study this dependence using larger lattices and also include anti-neutrinos.

preprint2020arXiv

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multi-messenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning (ML) formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two mono-energetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein (MSW) resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results intimate the promise of this "variational annealing" methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

preprint2020arXiv

Spin versus Helicity Equilibration Times and Lagrangian for Strange Quarks in Rotating Quark-Gluon Plasma

Measurements of the net polarization of $Λ$ and $\barΛ$ hyperons at the Relativistic Heavy Ion Collider (RHIC) have stimulated much interest in how strange quarks might align their spin with the vorticity of the matter created in heavy ion collisions. We calculate the Lagrangian in the rest frame of a fluid element undergoing rotation with angular velocity $ω$ including photon and gluon fields. There is an additional coupling between the quarks and the gauge fields proportional to $ω$, but this vertex does not change the spin of the quarks. We also show that the times to equilibrate quark helicity and spin parallel to the vorticity are the same so long as $ω$ is small compared to the temperature.

preprint2020arXiv

Spin-Vorticity Coupling for Massive Vector Mesons

Recent experiments at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) have indicated that hadrons containing strange quarks produced in non-central heavy ion collisions can be polarized. We investigate in detail the coupling of spin and vorticity for electrically neutral, massive vector bosons using the Proca equation, and provide the nonrelativistic reduction of the field equations via a single Foldy--Wouthuysen transformation. We find that the resulting Hamiltonian is not-Hermitian, but ${\cal PT}$ invariant, and involves a spin dependent term $\textstyle{\frac{1}{2}} s_z \hbar ω$ to leading order in vorticity. We also calculate further relativistic and quantum corrections to the Hamiltonian.

preprint2019arXiv

Is Hyperon Polarization in Relativistic Heavy Ion Collisions Connected to Axial U(1) Symmetry Breaking at High Temperature?

Experiments at the Relativistic Heavy Ion Collider (RHIC) have measured the net polarization of $Λ$ and $\barΛ$ hyperons and attributed it to a coupling between their spin and the vorticity of the fluid created in heavy ion collisions, but how the spin comes to equilibrium with vorticity is an open problem. Recently we found that vorticity fluctuations and helicity flip of strange quarks in quark-gluon plasma through perturbative QCD processes resulted in equilibration times far too long to be relevant. Here we consider the Nambu--Jona-Lasinio model with the inclusion of the six-quark Kobayashi--Maskawa--'t Hooft interaction which breaks axial U(1). Using instanton inspired models for the temperature dependence of the axial symmetry breaking, we find that constituent strange quarks can reach spin equilibrium at temperatures below about 170 MeV, just before they hadronize to form hyperons.

preprint2019arXiv

Relaxation Time for Strange Quark Spin in Rotating Quark-Gluon Plasma

Experiments at the Relativistic Heavy Ion Collider (RHIC) have measured the net polarization of $Λ$ and $\barΛ$ hyperons and attributed it to a coupling between their spin and the vorticity of the fluid created in heavy ion collisions. Equipartition of energy is generally assumed, but the dynamical mechanism which polarizes them has yet to be determined. We consider two such mechanisms: vorticity fluctuations and helicity flip in scatterings between strange quarks and light quarks and gluons. With reasonable parameters both mechanisms lead to equilibration times orders of magnitude too large to be relevant to heavy ion collisions. Our conclusion is that strange quark spin or helicity is unchanged from the time they are created to the time they hadronize. A corollary is that vorticity fluctuations do not affect the hyperon spin either.