Source author record

Ermal Rrapaj

Ermal Rrapaj 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

16works
9topics
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

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

preprint2016arXiv

Microscopically constrained mean field models from chiral nuclear thermodynamics

We explore the use of mean field models to approximate microscopic nuclear equations of state derived from chiral effective field theory across the densities and temperatures relevant for simu- lating astrophysical phenomena such as core-collapse supernovae and binary neutron star mergers. We consider both relativistic mean field theory with scalar and vector meson exchange as well as energy density functionals based on Skyrme phenomenology and compare to thermodynamic equa- tions of state derived from chiral two- and three-nucleon forces in many-body perturbation theory. Quantum Monte Carlo simulations of symmetric nuclear matter and pure neutron matter are used to determine the density regimes in which perturbation theory with chiral nuclear forces is valid. Within the theoretical uncertainties associated with the many-body methods, we find that select mean field models describe well microscopic nuclear thermodynamics. As an additional consistency requirement, we study as well the single-particle properties of nucleons in a hot/dense environment, which affect e.g., charged-current weak reactions in neutron-rich matter. The identified mean field models can be used across a larger range of densities and temperatures in astrophysical simulations than more computationally expensive microscopic models.

preprint2016arXiv

Nucleon-nucleon bremsstrahlung of dark gauge bosons and revised supernova constraints

We calculate the rate of production of hypothetical light vector bosons (LVBs) from nucleon-nucleon bremsstrahlung reactions in the soft radiation limit directly in terms of the measured nucleon-nucleon elastic cross sections. We use these results and the observation of neutrinos from supernova SN1987a to deduce constraints on the couplings of vector bosons with masses $\lesssim 200$ MeV to either electric charge (dark photons) or to baryon number. We establish for the first time strong constraints on LVB that couple only to baryon number, and revise earlier constraints on the dark photon. For the latter, we find that the excluded region of parameter space is diminished by about a factor of 10.

preprint2015arXiv

Charged-current reactions in the supernova neutrino-sphere

We calculate neutrino absorption rates due to charged-current reactions $ν_e+n \rightarrow e^- + p $ and $\barν_e+p \rightarrow e^+ + n $ in the outer regions of a newly born neutron star called the neutrino-sphere. To improve on recent work which has shown that nuclear mean fields enhance the $ν_e$ cross-section and suppress the $\barν_e$ cross-section, we employ realistic nucleon-nucleon interactions that fit measured scattering phase shifts. Using these interactions we calculate the momentum-, density-, and temperature-dependent nucleon self-energies in the Hartree-Fock approximation. A potential derived from chiral effective field theory and a pseudo-potential constructed to reproduce nucleon-nucleon phase shifts at the mean-field level are used to study the equilibrium proton fraction and the charged-current rates are studied in detail. We compare our results to earlier calculations obtained using phenomenological mean-field models and to those obtained in the virial expansion valid at low density. We find that for typical ambient conditions in the neutrino-sphere, $T=5-10$ MeV and $ρ=10^{11}-10^{13}$ g/cm$^3$, the difference between the $ν_{e}$ and $\barν_e$ absorption rates is much larger than predicted earlier. Our results have important implications for heavy-element nucleosynthesis in supernovae and for supernova neutrino detection.

preprint2014arXiv

Electron-neutron scattering and transport properties of neutron stars

We show that electrons can couple to the neutron excitations in neutron stars and find that this can limit their contribution to the transport properties of dense matter, especially the shear viscosity. The coupling between electrons and neutrons is induced by protons in the core, and by ions in the crust. We calculate the effective electron-neutron interaction for the kinematics of relevance to the scattering of degenerate electrons at high density. We use this interaction to calculate the electron thermal conductivity, electrical conductivity, and shear viscosity in the neutron star inner crust, and in the core where we consider both normal and superfluid phases of neutron-rich matter. In some cases, particularly when protons are superconducting and neutrons are in their normal phase, we find that electron-neutron scattering can be more important than the other scattering mechanisms considered previously.