Source author record

Alessandra Feo

Alessandra Feo 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

22works
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

22 published item(s)

preprint2022arXiv

Investigating the migration of immiscible contaminant fluid flow in homogeneous and heterogeneous aquifers with high-precision numerical simulations

Numerical modeling of the migration of three-phase immiscible fluid flow in variably saturated zones is challenging due to the different behavior of the system between unsaturated and saturated zones. This behavior results in the use of different numerical methods for the numerical simulation of the fluid flow depending on whether it is in the unsaturated or saturated zones. This paper shows that using a high-resolution shock-capturing conservative method to resolve the nonlinear governing coupled partial differential equations of a three-phase immiscible fluid flow allows the numerical simulation of the system through both zones providing a unitary vision (and resolution) of the migration of an immiscible contaminant problem within a porous medium. In particular, using different initial scenarios (including impermeable 'lenses' in heterogeneous aquifers), three-dimensional numerical simulation results are presented on the temporal evolution of the contaminant migration following the saturation profiles of the three-phases fluids flow in variably saturated zones. It is considered either light nonaqueous phase liquid with a density less than the water, or dense nonaqueous phase liquid, which has densities greater than the water initially released in unsaturated dry soil. Our study shows that the fate of the migration of immiscible contaminants in variably saturated zones can be accurately described, using a unique mathematical conservative model, with different evolution depending on the value of the system's physical parameters, including the contaminant density, and accurately tracking the evolution of the sharp (shock) contaminant front.

preprint2020arXiv

Numerical-relativity simulations of long-lived remnants of binary neutron star mergers

We analyze the properties of the gravitational wave signal emitted after the merger of a binary neutron star system when the remnant survives for more than a 80 ms (and up to 140ms). We employ four different piecewise polytropic equations of state supplemented by an ideal fluid thermal component. We find that the post-merger phase can be subdivided into three phases: an early post-merger phase (where the quadrupole mode and a few subdominant features are active), the intermediate post-merger phase (where only the quadrupole mode is active) and the late post-merger phase (where convective instabilities trigger inertial modes). The inertial modes have frequencies somewhat smaller than the quadrupole modes. In one model, we find an interesting association of a corotation of the quadrupole mode in parts of the star with a revival of its amplitude. The gravitational wave emission of inertial modes in the late post-merger phase is concentrated in a narrow frequency region and is potentially detectable by the planned third-generation detectors. This allows for the possibility of probing not only the cold part of the equation of state, but also its dependence on finite temperature. In view of these results, it will be important to investigate the impact of various type of viscosities on the potential excitation of inertial modes in binary neutron star merger remnants.

preprint2016arXiv

Binary neutron star merger simulations with different initial orbital frequency and equation of state

We present results from three-dimensional general relativistic simulations of binary neutron star coalescences and mergers using public codes. We considered equal mass models where the baryon mass of the two Neutron Stars (NS) is $1.4M_{\odot}$, described by four different equations of state (EOS) for the cold nuclear matter (APR4, SLy, H4, and MS1; all parametrized as piecewise polytropes). We started the simulations from four different initial interbinary distances ($40, 44.3, 50$, and $60$ km), including up to the last 16 orbits before merger. That allows to show the effects on the gravitational wave phase evolution, radiated energy and angular momentum due to: the use of different EOSs, the orbital eccentricity present in the initial data and the initial separation (in the simulation) between the two stars. Our results show that eccentricity has a major role in the discrepancy between numerical and analytical waveforms until the very last few orbits, where "tidal" effects and missing high-order post-Newtonian coefficients also play a significant role. We test different methods for extrapolating the gravitational wave signal extracted at finite radii to null infinity. We show that an effective procedure for integrating the Newman-Penrose $ψ_4$ signal to obtain the gravitational wave strain $h$ is to apply a simple high-pass digital filter to $h$ after a time domain integration, where only the two physical motivated integration constants are introduced. That should be preferred to the more common procedures of introducing additional integration constants, integrating in the frequency domain or filtering $ψ_4$ before integration.

preprint2016arXiv

Modeling Equal and Unequal Mass Binary Neutron Star Mergers Using Public Codes

We present three-dimensional simulations of the dynamics of binary neutron star (BNS) mergers 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). We investigate five equal-mass models of total gravitational mass $2.207$, $2.373$, $2.537$, $2.697$ and $2.854 M_\odot$, respectively, and four unequal mass models with $M_{\mathrm{ADM}}\simeq 2.53\ M_\odot$ and $q\simeq 0.94$, $0.88$, $0.82$, and $0.77$ (where $q = M^{(1)}/M^{(2)}$ is the mass ratio). We use a semi-realistic equation of state (EOS) namely, the seven-segment piece-wise polytropic SLyPP with a thermal component given by $Γ_{th} = 1.8$. We have also compared the resulting dynamics (for one model) using both, the BSSN-NOK and CCZ4 methods for the evolution of the gravitational sector, and also different reconstruction methods for the matter sector, namely PPM, WENO and MP5. Our results show agreement and high resolution, but superiority of BSSN-NOK supplemented by WENO reconstruction at lower resolutions. One of the important characteristics of the present investigation is that, for the first time, this has been done using only publicly available open source software, in particular, the Einstein Toolkit code deployed for the dynamical evolution and the LORENE code for the generation of the initial models. All of the source code and parameters used for the simulations have been made publicly available. This not only makes it possible to re-run and re-analyze our data; it also enables others to directly build upon this work for future research.

preprint2014arXiv

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

We present results about the effect of the use of a stiffer equation of state, namely the ideal-fluid $Γ=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 $β$ for the emergence of the instability when the adiabatic index $Γ$ 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 $β_c$. We also extend the analysis to lower values of $β$ and show that low-beta shear instabilities are present also in the case of matter described by a simple polytropic equation of state.

preprint2014arXiv

Stiffness effects on the dynamics of the bar-mode instability of Neutron Stars in full General Relativity

We present results on the effect of the stiffness of the equation of state on the dynamical bar-mode instability in rapidly rotating polytropic models of neutron stars in full General Relativity. We determine the change in the threshold for the emergence of the instability for a range of the adiabatic $Γ$ index from 2.0 to 3.0, including two values chosen to mimic more realistic equations of state at high densities.

preprint2012arXiv

Thermodynamics of SU(N) Yang-Mills theories in 2+1 dimensions II - The deconfined phase

We present a non-perturbative study of the equation of state in the deconfined phase of Yang-Mills theories in D=2+1 dimensions. We introduce a holographic model, based on the improved holographic QCD model, from which we derive a non-trivial relation between the order of the deconfinement phase transition and the behavior of the trace of the energy-momentum tensor as a function of the temperature T. We compare the theoretical predictions of this holographic model with a new set of high-precision numerical results from lattice simulations of SU(N) theories with N=2, 3, 4, 5 and 6 colors. The latter reveal that, similarly to the D=3+1 case, the bulk equilibrium thermodynamic quantities (pressure, trace of the energy-momentum tensor, energy density and entropy density) exhibit nearly perfect proportionality to the number of gluons, and can be successfully compared with the holographic predictions in a broad range of temperatures. Finally, we also show that, again similarly to the D=3+1 case, the trace of the energy-momentum tensor appears to be proportional to T^2 in a wide temperature range, starting from approximately 1.2 T_c, where T_c denotes the critical deconfinement temperature.

preprint2011arXiv

Thermodynamics of SU(N) Yang-Mills theories in 2+1 dimensions I - The confining phase

We compute the equation of state in the confining phase of SU(N) Yang-Mills theories with N=2, 3, 4, 5 and 6 colors in 2+1 dimensions, via lattice simulations. At low enough temperatures, the results are accurately described by a gas of glueballs, including all known states below the two-particle threshold. Close to the deconfinement temperature, however, this prediction underestimates the numerical results, and the contribution from heavier glueballs has to be included. We show that the spectral density of the latter can be accurately described using a simple bosonic string model.

preprint2011arXiv

Universal signatures of the effective string in finite temperature lattice gauge theories

We study the behaviour of the interquark potential in lattice gauge theories at high temperature, but still in the confining phase, and propose a new observable which could play in this regime the same role played by the Luscher term in the low temperature limit. This quantity is related to the exponent of the power prefactor in the effective string partition function or, equivalently, to the coefficient of the logarithmic correction in the interquark potential and, as for the usual Luscher term, its value does not depend on the particular gauge group under consideration or on the form of the effective string action used to model the flux tube. In this respect it can be considered as a universal signature of the effective string behaviour of the flux tube. As a test of our proposal we studied this quantity with a set of high-precision numerical simulations in the (2+1) dimensional SU(2), SU(3) and SU(4) Yang-Mills theories and in the Z_2 gauge model, always finding a perfect agreement with the predicted values.

preprint2010arXiv

Species Doublers as Super Multiplet Partners in Lattice Supersymmetry

We propose a new lattice superfield formalism in momentum representation which accommodates species doublers of the lattice fermions and their bosonic counterparts as super multiplets. We explicitly show that one dimensional $N=2$ model with interactions has exact supersymmetry on the lattice for all super charges with lattice momentum. In coordinate representation the finite difference operator is made to satisfy Leibnitz rule by introducing a non local product, the "star" product, and the exact lattice supersymmetry is realized. Supersymmetric Ward identities are shown to be satisfied at one loop level.

preprint2010arXiv

Species Doublers as Super Multiplets in Lattice Supersymmetry: Exact Supersymmetry with Interactions for D=1 N=2

We propose a new lattice superfield formalism in momentum representation which accommodates species doublers of the lattice fermions and their bosonic counterparts as super multiplets. We explicitly show that one dimensional N=2 model with interactions has exact Lie algebraic supersymmetry on the lattice for all super charges. In coordinate representation the finite difference operator is made to satisfy Leibnitz rule by introducing a non local product, the ``star'' product, and the exact lattice supersymmetry is realized. The standard momentum conservation is replaced on the lattice by the conservation of the sine of the momentum, which plays a crucial role in the formulation. Half lattice spacing structure is essential for the one dimensional model and the lattice supersymmetry transformation can be identified as a half lattice spacing translation combined with alternating sign structure. Invariance under finite translations and locality in the continuum limit are explicitly investigated and shown to be recovered. Supersymmetric Ward identities are shown to be satisfied at one loop level. Lie algebraic lattice supersymmetry algebra of this model suggests a close connection with Hopf algebraic exactness of the link approach formulation of lattice supersymmetry.

preprint2010arXiv

Thermodynamics of SU(N) gauge theories in 2+1 dimensions in the $T <\ T_c$ regime

We present Monte Carlo results for the thermodynamics of pure SU(N) gauge theories with $N=2,...,6$ in 2+1 dimensions. We focus on the confined phase region $T<T_c$ and study thermodynamics variables such as the trace of the energy-momentum tensor, pressure, energy and entropy density using the integral method. We also investigate scaling properties with $N$ of the different observables. We compare our results with a gas of free glueballs and the bosonic string predictions for the Hagedorn spectrum.

preprint2009arXiv

Lattice Supersymmetry: Some Ideas from Low Dimensional Models

In the framework of the so called link approach we study exact lattice supersymmetry for the simplest supersymmetric model: N=1 supersymmetry in D=1. The model is described by a lattice with spacing a/2, thus containing twice as many sites as the conventional one. The boson and fermion are related through a 2pi/a momentum shift, which can provide an interpretation of them being species doublers to each other. An exactly supersymmetric lattice action can be written within this scheme in momentum representation, which however turns out to be non local in coordinate space.

preprint2005arXiv

Density Matrix Renormalization Group Study of the Disorder Line in the Quantum ANNNI Model

We apply Density Matrix Renormalization Group methods to study the phase diagram of the quantum ANNNI model in the region of low frustration where the ferromagnetic coupling is larger than the next-nearest-neighbor antiferromagnetic one. By Finite Size Scaling on lattices with up to 80 sites we locate precisely the transition line from the ferromagnetic phase to a paramagnetic phase without spatial modulation. We then measure and analyze the spin-spin correlation function in order to determine the disorder transition line where a modulation appears. We give strong numerical support to the conjecture that the Peschel-Emery one-dimensional line actually coincides with the disorder line. We also show that the critical exponent governing the vanishing of the modulation parameter at the disorder transition is $β_q = 1/2$.

preprint2005arXiv

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.

preprint2003arXiv

Gauge Theories on a 2+2 Anisotropic Lattice

The implementation of gauge theories on a four-dimensional anisotropic lattice with two distinct lattice spacings is discussed, with special attention to the case where two axes are finely and two axes are coarsely discretized. Feynman rules for the Wilson gauge action are derived and the renormalizability of the theory and the recovery of the continuum limit are analyzed. The calculation of the gluon propagator and the restoration of Lorentz invariance in on-shell states is presented to one-loop order in lattice perturbation theory for $SU(N_c)$ on both 2+2 and 3+1 lattices.

preprint2001arXiv

Hamiltonian Monte Carlo simulation of the two-dimensional Wess-Zumino model

We study a Hamiltonian lattice version of the two-dimensional Wess-Zumino model. Preliminary results obtained by Quantum Monte Carlo with a many-parameter guiding wave function are presented. We analyze the pattern of supersymmetry breaking by measuring the ground state energy and a set of supersymmetric Ward identities. The algorithm is quite effective and allows very precise measurements.

preprint2001arXiv

SUSY Ward identities in 1-loop perturbation theory

We present preliminary results of a study of the supersymmetric (SUSY) Ward identities (WIs) for the N=1 SU(2) SUSY Yang-Mills theory in the context of one-loop lattice perturbation theory. The supersymmetry on the lattice is explicitly broken by the gluino mass and the lattice artifacts. However, the renormalization of the supercurrent can be carried out in a scheme that restores the nominal continuum WIs. The perturbative calculation of the renormalization constants and mixing coefficients for the local supercurrent is presented.

preprint2000arXiv

On the 1-loop lattice perturbation theory of the supersymmetric Ward identities

The one loop corrections to the supersymmetric Ward identities (WIs) in the discretized N=1 SU(2) supersymmetric Yang-Mills theory can be investigated by means of lattice perturbation theory. The supersymmetry (SUSY) is explicitly broken by the lattice discretization as well as by the introduction of Wilson fermions. However, the renormalization of the supercurrent can be carried out in a scheme that restores the nominal continuum WIs. We present our work in progress which is concerned with the 1-loop renormalization of the local supercurrent, i.e. with the perturbative computation of the corresponding renormalization constants and mixing coefficients.