First-order viscous relativistic hydrodynamics has long been thought to be unsta- ble and acausal. This is not true; it is only with certain definitions of the hydrody- namic variables that the equations of motion display these properties. It is possible to define the hydrodynamic variables such that a fluid is both stable and causal at first order. This thesis does so for both uncharged and charged fluids, mostly for fluids at rest. Work has also been done in limited cases on fluids in motion. A class of stable and causal theories is identified via constraints on transport coefficients derived from linearized perturbations of the equilibrium state. Causality conditions are also derived for the full non-linear hydrodynamic equations. / Graduate

This Thesis puts its focus on the physics of neutron stars and its description with methods of numerical relativity. In the first step, a new numerical framework the Whisky2D code will be developed, which solves the relativistic equations of hydrodynamics in axisymmetry. Therefore we consider an improved formulation of the conserved form of these equations. The second part will use the new code to investigate the critical behaviour of two colliding neutron stars. Considering the analogy to phase transitions in statistical physics, we will investigate the evolution of the entropy of the neutron stars during the whole process. A better understanding of the evolution of thermodynamical quantities, like the entropy in critical process, should provide deeper understanding of thermodynamics in relativity. More specifically, we have written the Whisky2D code, which solves the general-relativistic hydrodynamics equations in a flux-conservative form and in cylindrical coordinates. This of course brings in 1/r singular terms, where r is the radial cylindrical coordinate, which must be dealt with appropriately. In the above-referenced works, the flux operator is expanded and the 1/r terms, not containing derivatives, are moved to the right-hand-side of the equation (the source term), so that the left hand side assumes a form identical to the one of the three-dimensional (3D) Cartesian formulation. We call this the standard formulation. Another possibility is not to split the flux operator and to redefine the conserved variables, via a multiplication by r. We call this the new formulation. The new equations are solved with the same methods as in the Cartesian case. From a mathematical point of view, one would not expect differences between the two ways of writing the differential operator, but, of course, a difference is present at the numerical level. Our tests show that the new formulation yields results with a global truncation error which is one or more orders of magnitude smaller than those of alternative and commonly used formulations. The second part of the Thesis uses the new code for investigations of critical phenomena in general relativity. In particular, we consider the head-on-collision of two neutron stars in a region of the parameter space where two final states a new stable neutron star or a black hole, lay close to each other. In 1993, Choptuik considered one-parameter families of solutions, S[P], of the Einstein-Klein-Gordon equations for a massless scalar field in spherical symmetry, such that for every P > P⋆, S[P] contains a black hole and for every P < P⋆, S[P] is a solution not containing singularities. He studied numerically the behavior of S[P] as P → P⋆ and found that the critical solution, S[P⋆], is universal, in the sense that it is approached by all nearly-critical solutions regardless of the particular family of initial data considered. All these phenomena have the common property that, as P approaches P⋆, S[P] approaches a universal solution S[P⋆] and that all the physical quantities of S[P] depend only on |P − P⋆|. The first study of critical phenomena concerning the head-on collision of NSs was carried out by Jin and Suen in 2007. In particular, they considered a series of families of equal-mass NSs, modeled with an ideal-gas EOS, boosted towards each other and varied the mass of the stars, their separation, velocity and the polytropic index in the EOS. In this way they could observe a critical phenomenon of type I near the threshold of black-hole formation, with the putative solution being a nonlinearly oscillating star. In a successive work, they performed similar simulations but considering the head-on collision of Gaussian distributions of matter. Also in this case they found the appearance of type-I critical behaviour, but also performed a perturbative analysis of the initial distributions of matter and of the merged object. Because of the considerable difference found in the eigenfrequencies in the two cases, they concluded that the critical solution does not represent a system near equilibrium and in particular not a perturbed Tolmann-Oppenheimer-Volkoff (TOV) solution. In this Thesis we study the dynamics of the head-on collision of two equal-mass NSs using a setup which is as similar as possible to the one considered above. While we confirm that the merged object exhibits a type-I critical behaviour, we also argue against the conclusion that the critical solution cannot be described in terms of equilibrium solution. Indeed, we show that, in analogy with what is found in, the critical solution is effectively a perturbed unstable solution of the TOV equations. Our analysis also considers fine-structure of the scaling relation of type-I critical phenomena and we show that it exhibits oscillations in a similar way to the one studied in the context of scalar-field critical collapse. / Diese Arbeit legt seinen Schwerpunkt auf die Physik von Neutronensternen und deren Beschreibung mit Methoden der numerischen Relativitätstheorie. Im ersten Schritt wird eine neue numerische Umgebung, der Whisky2D Code entwickelt, dieser löst die relativistischen Gleichungen der Hydrodynamik in Axialymmetrie. Hierzu betrachten wir eine verbesserte Formulierung der sog. "flux conserved formulation" der Gleichungen. Im zweiten Teil wird der neue Code verwendet / um das kritische Verhalten zweier kollidierenden Neutronensternen zu untersuchen. In Anbetracht der Analogie, um Übergänge in der statistischen Physik Phase werden wir die Entwicklung der Entropie der Neutronensterne während des gesamten Prozesses betrachten. Ein besseres Verständnis der Evolution von thermodynamischen Größen, wie der Entropie in kritischer Prozess, sollte zu einem tieferen Verständnis der relativistischen Thermodynamik führen. Der Whisky2D Code, zur Lösung Gleichungen relativistischer Hydrodynamik wurde in einer „flux conserved form“ und in zylindrischen Koordinaten geschrieben. Hierdurch entstehen 1 / r singuläre Terme, wobei r der ist, die entsprechend behandelt werden müssen. In früheren Arbeiten, wird der Operator expandiert und die 1 / r spezifisch Therme auf die rechte Seite geschrieben, so dass die linke Seite eine Form annimmt, die identisch ist mit der kartesischen Formulierung. Wir nennen dies die Standard-Formulierung. Eine andere Möglichkeit ist, die Terme nicht zu expandieren, den und den 1/r Term in die Gleichung hinein zu ziehen. Wir nennen dies die Neue-Formulierung. Die neuen Gleichungen werden mit den gleichen Verfahren wie im kartesischen Fall gelöst. Aus mathematischer Sicht ist keine Unterschiede zwischen den beiden Formulierungen zu erwarten, erst die numerische Sicht zeigt die Unterschiede auf. Versuche zeigen, dass die Neue-Formulierung numerische Fehler um mehrere Größenordnungen reduziert. Der zweite Teil der Dissertation verwendet den neuen Code für die Untersuchung kritischer Phänomene in der allgemeinen Relativitätstheorie. Insbesondere betrachten wir die Kopf-auf-Kollision zweier Neutronensterne in einem Bereich des Parameter Raums, deren zwei mögliche Endzustände entweder einen neuen stabilen Neutronenstern oder ein Schwarzes Loch darstellen. Im Jahr 1993, betrachtete Choptuik Ein-Parameter-Familien von Lösungen, S [P], der Einstein-Klein-Gordon-Gleichung für ein masseloses Skalarfeld in sphärischer Symmetrie, so dass für jedes P> P ⋆, S[P] ein Schwarzes Loch enthalten ist und jedes P <P ⋆, S [P] eine Lösung ohne eine Singularität. Er studierte das numerisch Verhalten von S [P] als P → P ⋆ und stellte fest, dass die kritische Lösung, S[P ⋆], universell in dem Sinne ist, dass die Lösung nahe des kritischen Wertes unabhängig von den Anfangsdaten ist. Alle diese Phänomene haben die gemeinsame Eigenschaft, dass sobald sich P P ⋆ annähert, auch S[P] S [P ⋆] annähert und dass alle physikalischen Größen von S [P] nur noch von | P - P ⋆ | abhängen. Die erste Studie der kritischen Phänomene über den Frontalzusammenstoß von zwei Neutronensternen wurde von Jin und Suen im Jahr 2007 durchgeführt. Insbesondere untersuchten sie eine Reihe von Anfangsdaten gleicher Neutronensternmasse, mit einer idealen EOS, die aufeinander zu beschleunigt werden. Variiert wurden die Massen der Sterne, ihr Abstand, die Geschwindigkeit und die polytropen Index der EOS. Auf diese Weise konnten sie kritische Phänomen des Typ I beobachten. In weiteren Versuchen, führten sie ähnliche Simulationen frontal kollidierender Materie in Gauß Verteilungen durch. Auch in diesem Fall fanden sie Typ-I-kritisches Verhalten. Zudem führten sie eine störungstheoretische Analyse der Anfangsobjekte als auch der stabilen Endobjekte durch. Wegen der beträchtlichen Unterschiede in den Eigenfrequenzen in beiden Fallen, schlossen sie daraus, dass die kritische Lösung keine linear-gestörten Tolmann-Oppenheimer-Volkoff (TOV) Sterne im Gleichgewicht darstellen. In dieser Arbeit untersuchen wir die Dynamik der Frontalzusammenstoß zweier Neutronensterne gleicher Masse mit ähnlichem Setup wie oben besprochen. Während wir bestätigen, dass die erzeugten Objekte ebenfalls ein Typ-I-kritische Verhalten aufweisen, wiedersprechen wir der der Aussage, dass sich die kritische Lösung nicht als Gleichgewichtslösung dargestellt werden kann. In der Tat zeigen wir, dass die kritische Lösung als linear-gestörte instabile Lösung eines TOV-Sterns dargestellt werden kann. Unsere Analyse berücksichtigt auch die Feinstruktur der Skalenverhältnisse Typ-I-kritischer Phänomene und wir zeigen ebenfalle, dass hier Oszillationen auftreten, die bereits in ähnlichen Studien zum Kollaps kritischer Skalar-Feld gefunden wurden.

No description available.

Neste trabalho usamos um modelo hidrodinâmico com condições iniciais flutuantes evento a evento para estudar alguns observáveis em colisões de íons pesados relativísticos. Mostramos que a introdução de flutuações, caracterizadas por tubos longitudinais de alta densidade de energia, torna o cálculo dos observáveis mais próximo dos dados, em comparação à abordagem tradicional, na qual as condições iniciais são parametrizadas através de funções suaves e simétricas. No caso da distribuição de partículas em função do momento transversal, observamos um aumento do número de partículas na região onde p_t>2.0GeV, como consequência da expansão rápida dos tubos de alta densidade de energia, principalmente os tubos próximos à superfície da matéria, que chamamos de tubos periféricos. No caso de v_2 em função do momento transversal, observamos uma redução do fluxo elíptico na mesma região, devido à expansão isotrópica dos tubos periféricos. No caso de v_2 em função da pseudo-rapidez, observamos, com o uso das condições iniciais flutuantes, uma curva com um perfil triangular muito próximo dos dados. No caso da função de correlação de duas partículas, observamos que os tubos periféricos dão origem às chamadas estruturas de Ridge. Estas estruturas representam correlações de longo alcance, na direção longitudinal, entre as partículas trigger e as partículas associadas. Obtivemos uma função de correlação com três Ridges, respectivamente, nas posições: \\Delta\\phi=0, \\Delta\\phi=2 e \\Delta\\phi=-2. Mostramos, também, que o modelo reproduz o chamado efeito in-plane/out-of-plane, que está relacionado com a dependência da função de correlação de duas partículas em relação ao ângulo azimutal do momento da partícula trigger. De modo geral, observamos que, com o ajuste correto dos parâmetros do modelo, como por exemplo, a distribuição de densidade de energia inicial média e a temperatura de desacoplamento, todos os observáveis discutidos nesta tese apresentam um razoável acordo com os dados. / In this work we apply a hydrodynamic model with fluctuating initial conditions in a event by event basis to study some observables in relativistic heavy ion collisions. We show that the introduction of fluctuations, characterized by longitudinal high energy density tubes, turn the computation of the observables closer to data, in comparison to the traditional approach, in which the initial conditions are parametrized by using smooth and symmetrical functions. In the case of the transverse momentum distribution, we observe an enhancement of the number of particles in the region where p_t>2.0GeV, as a consequence of the fast expansion of the high energy density tubes, mainly the tubes that are close to the surface of the matter, called peripheral tubes. In the case of v_2 as a function of the transverse momentum, we observe a reduction of the elliptic flow at the same region, due to the isotropic expansion of the peripheral tubes. In the case of v_2 as a function of the pseudorapidity, we observe, by using the fluctuating initial conditions, a curve with a triangular shape very close to data. In the case of the two particle correlation function, we observed that the peripheral tubes give rise to the Ridge structures. These structures represent a long range correlation, in the longitudinal direction, between the trigger particles and the associated particles. We got a two particle correlation function with three Ridges, respectively, at the positions: \\Delta\\phi=0, \\Delta\\phi=2 and \\Delta\\phi=-2. In addition, we show that the model reproduces the in-plane/out-of-plane effect, which is related with the two particle correlation function dependence on the azimuthal angle of the momentum of the trigger particle. Generally, we observed that, with the proper adjustment of the parameters of the model, as, for example, the average initial energy density and the freeze-out temperature, all observables discussed in this thesis show a reasonable agreement with data.

Neste trabalho usamos um modelo hidrodinâmico com condições iniciais flutuantes evento a evento para estudar alguns observáveis em colisões de íons pesados relativísticos. Mostramos que a introdução de flutuações, caracterizadas por tubos longitudinais de alta densidade de energia, torna o cálculo dos observáveis mais próximo dos dados, em comparação à abordagem tradicional, na qual as condições iniciais são parametrizadas através de funções suaves e simétricas. No caso da distribuição de partículas em função do momento transversal, observamos um aumento do número de partículas na região onde p_t>2.0GeV, como consequência da expansão rápida dos tubos de alta densidade de energia, principalmente os tubos próximos à superfície da matéria, que chamamos de tubos periféricos. No caso de v_2 em função do momento transversal, observamos uma redução do fluxo elíptico na mesma região, devido à expansão isotrópica dos tubos periféricos. No caso de v_2 em função da pseudo-rapidez, observamos, com o uso das condições iniciais flutuantes, uma curva com um perfil triangular muito próximo dos dados. No caso da função de correlação de duas partículas, observamos que os tubos periféricos dão origem às chamadas estruturas de Ridge. Estas estruturas representam correlações de longo alcance, na direção longitudinal, entre as partículas trigger e as partículas associadas. Obtivemos uma função de correlação com três Ridges, respectivamente, nas posições: \\Delta\\phi=0, \\Delta\\phi=2 e \\Delta\\phi=-2. Mostramos, também, que o modelo reproduz o chamado efeito in-plane/out-of-plane, que está relacionado com a dependência da função de correlação de duas partículas em relação ao ângulo azimutal do momento da partícula trigger. De modo geral, observamos que, com o ajuste correto dos parâmetros do modelo, como por exemplo, a distribuição de densidade de energia inicial média e a temperatura de desacoplamento, todos os observáveis discutidos nesta tese apresentam um razoável acordo com os dados. / In this work we apply a hydrodynamic model with fluctuating initial conditions in a event by event basis to study some observables in relativistic heavy ion collisions. We show that the introduction of fluctuations, characterized by longitudinal high energy density tubes, turn the computation of the observables closer to data, in comparison to the traditional approach, in which the initial conditions are parametrized by using smooth and symmetrical functions. In the case of the transverse momentum distribution, we observe an enhancement of the number of particles in the region where p_t>2.0GeV, as a consequence of the fast expansion of the high energy density tubes, mainly the tubes that are close to the surface of the matter, called peripheral tubes. In the case of v_2 as a function of the transverse momentum, we observe a reduction of the elliptic flow at the same region, due to the isotropic expansion of the peripheral tubes. In the case of v_2 as a function of the pseudorapidity, we observe, by using the fluctuating initial conditions, a curve with a triangular shape very close to data. In the case of the two particle correlation function, we observed that the peripheral tubes give rise to the Ridge structures. These structures represent a long range correlation, in the longitudinal direction, between the trigger particles and the associated particles. We got a two particle correlation function with three Ridges, respectively, at the positions: \\Delta\\phi=0, \\Delta\\phi=2 and \\Delta\\phi=-2. In addition, we show that the model reproduces the in-plane/out-of-plane effect, which is related with the two particle correlation function dependence on the azimuthal angle of the momentum of the trigger particle. Generally, we observed that, with the proper adjustment of the parameters of the model, as, for example, the average initial energy density and the freeze-out temperature, all observables discussed in this thesis show a reasonable agreement with data.

O Plasma de Quarks e Glúons é um novo estado da matéria criado brevemente em colisões nucleares de alta energia no RHIC (acelerador localizado no Brookhaven National Laboratory) e no LHC (acelerador do CERN, na Suíça). Para pesquisar suas propriedades e estado inicial, várias análises e métodos são usados. Neste mestrado, foi pesquisada a recém proposta Análise em Componentes Principais para multiplicidade e fluxo anisotrópico. É discutida sua sensibilidade ao tamanho das inomogeneidades no estado inicial bem como (de maneira indireta) ao tipo de condições iniciais e energia da colisão. Foi feita também uma comparação com dados experimentais recém obtidos pela colaboração CMS do LHC. O resultado talvez mais promissor é a diferença qualitativa entre dados de flutuações de multiplicidade e resultados de simulações hidrodinâmicas com certas condições iniciais. / The Quark Gluon Plasma is a new state of matter created briefly in high energy nuclear collisions at RHIC (accelerator located at Brookhaven National Laboratory) and at LHC (accelerator at CERN, Switzerland). To research its properties and initial state, many analysis and methods are used. In this master, the recently proposed Principal Component Analysis was studied for multiplicity and anisotropic flow. Its sensitivity to the size of inhomogeneities in the initial state as well as (indirectly) to the kind of initial conditions and collision energy is discussed. Also, a comparison was done with experimental data recently obtained by the CMS collaboration from LHC. Maybe, the most promising result is the qualitative difference between data from multiplicity fluctuations and results from hydrodynamic simulations with different initial conditions.

No description available.

Cette thèse porte sur l'étude de différents aspects de la physique nucléaire depuis les basses énergies jusqu'aux énergies intermédiaires. Pour les basses énergies, où la matière nucléaire est essentiellement constituée de nucléons en interaction, la partie I traite de la fusion-fission des noyaux super-lourds, et la partie II des règles de somme associées aux interactions de type Skyrme. Pour les énergies intermédiaires, la matière nucléaire étant alors considérée comme une phase hadronique principalement constituée de pions, la partie III se focalise sur l'hydrodynamique relativiste de la matière nucléaire avec brisure spontanée de symétrie chirale. Dans la partie I, on s'intéresse à la formation puis à la désexcitation des noyaux super-lourds. On étudie donc la formation du noyau composé avec effets de mémoire. Pour la désexcitation d'un noyau super-lourd, l'existence d'un puits isomérique dans la barrière de potentiel change la dynamique de désexcitation et augmente les temps de fission. Cette dernière étude pourrait être utile à l'étude de la dynamique de la fission des actinides. Dans la partie II, les règles de somme M1 et M3 associées aux potentiels phénoménologiques de type Skyrme sont calculées à partir de leurs définitions intrinsèques. On détermine alors M1 jusqu'au niveau tensoriel et M3 avec potentiel central. Dans la partie III, pour le traitement hydrodynamique de la matière hadronique appliqué aux collisions d'ions lourds on peut, en première approximation, écarter les modifications induites par la brisure spontanée de symétrie chirale mais pas celles dues à l'aspect dissipatif. / This study focuses on different aspects of nuclear physics from low energies to intermediate ones. For the low energies, the nuclear matter is essentially constituted from interacting nucleons. Part I is on the fusion-fission of super-heavy elements, while Part II is on the Skyrme interactions-associated sum rules. In the case of the intermediate energies, where the nuclear matter is considered as being an hadronic phase mainly constituted from pions, Part III is focused on nuclear matter relativistic hydrodynamics with spontaneous chiral symmetry breaking. In Part I, the formation and the desexcitation of super-heavy nuclei are being studied. We analyzed the formation of compound nuclei including the memory effects. For super-heavy nuclei desexcitation, the existence of isomeric state within the potential barrier modifies the desexcitation dynamics and increases the fission time. This latter study could be useful for the study of the actinides fission. In Part II, the phenomenological Skyrme effective interactions- associated M1 and M3 sum rules are being calculated based on their intrinsic definitions. We identify then M1 up to the tensorial level and M3 with central potential. In Part III, as for the hadronic matter hydrodynamics being applied to heavy ions collisions, and as a first approach only, we can neglect spontaneous chiral symmetry but certainly not the dissipative impact.