The High Density Region of QCD from an Effective Model

We study the high density region of QCD within an effective model obtained in the frame of the hopping parameter expansion and choosing Polyakov-type loops as the main dynamical variables representing the fermionic matter. This model still shows the so-called sign problem, a difficulty peculiar to non-zero chemical potential, but it permits the development of algorithms which ensure a good overlap of the simulated Monte Carlo ensemble with the true one. We review the main features of the model and present results concerning the dependence of various observables on the chemical potential and on the temperature, in particular of the charge density and the Polykov loop susceptibility, which may be used to characterize the various phases expected at high baryonic density. In this way, we obtain information about the phase structure of the model and the corresponding phase transitions and cross over regions, which can be considered as hints about the behaviour of non-zero density QCD.Comment: 7 pages, 5 figures, talk presented at the XXVth International Symposium on Lattice Field Theory, July 30 - August 4 2007, Regensburg, German

Adaptive Optimization of Wave Functions for Lattice Field Models

The accuracy of Green Function Monte Carlo (GFMC) simulations can be greatly improved by a clever choice of the approximate ground state wave function that controls configuration sampling. This trial wave function typically depends on many free parameters whose fixing is a non trivial task. Here, we discuss a general purpose adaptive algorithm for their non-linear optimization. As a non trivial application we test the method on the two dimensional Wess-Zumino model, a relativistically invariant supersymmetric field theory with interacting bosonic and fermionic degrees of freedom.Comment: 12 pages, 5 EPS figures, Contribution to the Proceedings of the "Quantum Monte Carlo" meeting (Trento, Italy, July 3-6, 2001

Exact lattice Ward-Takahashi identity for the N=1 Wess-Zumino model

We consider a lattice formulation of the four dimensional N=1 Wess-Zumino model that uses the Ginsparg-Wilson relation. This formulation has an exact supersymmetry on the lattice. We show that the corresponding Ward-Takahashi identity is satisfied, both at fixed lattice spacing and in the continuum limit. The calculation is performed in lattice perturbation theory up to order $g^2$ in the coupling constant. We also show that this Ward-Takahashi identity determines the finite part of the scalar and fermion renormalization wave functions which automatically leads to restoration of supersymmetry in the continuum limit. In particular, these wave functions coincide in this limit.Comment: 19 pages, 6 figure

Neutron Star instabilities in full General Relativity using a Γ=2.75\Gamma=2.75 ideal fluid

We present results about the effect of the use of a stiffer equation of state, namely the ideal-fluid $\Gamma=2.75$ ones, on the dynamical bar-mode instability in rapidly rotating polytropic models of neutron stars in full General Relativity. We determine the change on the critical value of the instability parameter $\beta$ for the emergence of the instability when the adiabatic index $\Gamma$ is changed from 2 to 2.75 in order to mimic the behavior of a realistic equation of state. In particular, we show that the threshold for the onset of the bar-mode instability is reduced by this change in the stiffness and give a precise quantification of the change in value of the critical parameter $\beta_c$. We also extend the analysis to lower values of $\beta$ and show that low-beta shear instabilities are present also in the case of matter described by a simple polytropic equation of state.Comment: 16 pages, 16 figure

Modeling Mergers of Known Galactic Systems of Binary Neutron Stars

We present a study of the merger of six different known galactic systems of binary neutron stars (BNS) of unequal mass with a mass ratio between $0.75$ and $0.99$. Specifically, these systems are J1756-2251, J0737-3039A, J1906+0746, B1534+12, J0453+1559 and B1913+16. We follow the dynamics of the merger from the late stage of the inspiral process up to $\sim$ 20 ms after the system has merged, either to form a hyper-massive neutron star (NS) or a rotating black hole (BH), using a semi-realistic equation of state (EOS), namely the seven-segment piece-wise polytropic SLy with a thermal component. For the most extreme of these systems ($q=0.75$, J0453+1559), we also investigate the effects of different EOSs: APR4, H4, and MS1. Our numerical simulations are performed using only publicly available open source code such as, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. We show results on the gravitational wave signals, spectrogram and frequencies of the BNS after the merger and the BH properties in the two cases in which the system collapse within the simulated time.Comment: 13 pages, 10 figure

Spectral analysis of gravitational waves from binary neutron star merger remnants

In this work we analyze the gravitational wave signal from hypermassive neutron stars formed after the merger of binary neutron star systems, focusing on its spectral features. The gravitational wave signals are extracted from numerical relativity simulations of models already considered by De Pietri et al. [Phys. Rev. D 93, 064047 (2016)], Maione et al. [Classical Quantum Gravity 33, 175009 (2016)], and Feo et al. [Classical Quantum Gravity 34, 034001 (2017)], and allow us to study the effect of the total baryonic mass of such systems (from $2.4 M_{\odot}$ to $3 M_{\odot}$), the mass ratio (up to $q = 0.77$), and the neutron star equation of state, both in equal and highly unequal mass binaries. We use the peaks we find in the gravitational spectrum as an independent test of already published hypotheses of their physical origin and empirical relations linking them with the characteristics of the merging neutron stars. In particular, we highlight the effects of the mass ratio, which in the past was often neglected. We also analyze the temporal evolution of the emission frequencies. Finally, we introduce a modern variant of Prony's method to analyze the gravitational wave postmerger emission as a sum of complex exponentials, trying to overcome some drawbacks of both Fourier spectra and least-squares fitting. Overall, the spectral properties of the postmerger signal observed in our simulation are in agreement with those proposed by other groups. More specifically, we find that the analysis of Bauswein and Stergioulas [Phys. Rev. D 91, 124056 (2015)] is particularly effective for binaries with very low masses or with a small mass ratio and that the mechanical toy model of Takami et al. [Phys. Rev. D 91, 064001 (2015)] provides a comprehensive and accurate description of the early stages of the postmerger.Comment: 19 pages, 6 figure

Study of dynamical supersymmetry breaking for the two dimensional lattice Wess-Zumino model

A new approach to the study of the transition point in a class of two dimensional Wess-Zumino models is presented. The method is based on the calculation of rigorous lower bounds on the ground state energy density in the infinite lattice limit. Such bounds are useful in the discussion of supersymmetry phase transition. The transition point is then determined and compared with recent results based on large-scale Green Function Monte Carlo simulations with good agreement.Comment: 4 pages, 1 figure, Talk presented at Quark Confinement and the Hadron Spectrum VI, Villasimius, Sardinia, Italy, September 21-25, 200

Dynamical supersymmetry breaking and phase diagram of the lattice Wess-Zumino model

We study dynamical supersymmetry breaking and the transition point by non-perturbative lattice techniques in a class of two-dimensional N=1 Wess-Zumino model. The method is based on the calculation of rigorous lower bounds on the ground state energy density in the infinite-lattice limit. Such bounds are useful in the discussion of supersymmetry phase transition. The transition point is determined with this method and then compared with recent results based on large-scale Green Function Monte Carlo simulations with good agreement.Comment: 6 pages, 4 figs. Lattice2005 (theory), uses Po

The high density region of QCD in a large mass and chemical potential model

We study the high density region of QCD within an effective model obtained in the frame of the hopping parameter expansion. The model still acknowledges the sign problem peculiar to non-zero chemical potential, but it permits the development of refined algorithms which ensure a good overlap of the Monte Carlo ensemble with the true one. We review the main features of the model, including the most explicit form of the resumed expansion, and present calculations concerning the dependence of various observables on the chemical potential and on the temperature, in particular of the charge density and the diquark susceptibility, which may be used to characterize the various phases expected at high baryonic density.Comment: 6 pages, 5 figures, Lattice'05 proceedings, PoS styl