530 Physik
Refine
Year of publication
Document Type
- Article (1960)
- Preprint (1307)
- Doctoral Thesis (592)
- Conference Proceeding (242)
- diplomthesis (101)
- Bachelor Thesis (75)
- Master's Thesis (61)
- Contribution to a Periodical (56)
- Part of Periodical (36)
- Diploma Thesis (34)
Keywords
- Kollisionen schwerer Ionen (47)
- heavy ion collisions (44)
- LHC (25)
- Quark-Gluon-Plasma (25)
- Heavy Ion Experiments (21)
- equation of state (19)
- quark-gluon plasma (19)
- BESIII (16)
- QCD (16)
- Relativistic heavy-ion collisions (16)
Institute
- Physik (4269)
- Frankfurt Institute for Advanced Studies (FIAS) (1391)
- Informatik (978)
- Präsidium (70)
- MPI für Biophysik (43)
- ELEMENTS (40)
- Biochemie und Chemie (18)
- Helmholtz International Center for FAIR (12)
- Biochemie, Chemie und Pharmazie (11)
- Geowissenschaften (9)
Post-merger gravitational-wave signal from neutron-star binaries: a new look at an old problem
(2023)
The spectral properties of the post-merger gravitational-wave signal from a binary of neutron stars encodes a variety of information about the features of the system and of the equation of state describing matter around and above nuclear saturation density. Characterizing the properties of such a signal is an “old” problem, which first emerged when a number of frequencies were shown to be related to the properties of the binary through “quasiuniversal” relations. Here we take a new look at this old problem by computing the properties of the signal in terms of the Weyl scalar ψ4. In this way, and using a database of more than 100 simulations, we provide the first evidence for a new instantaneous frequency, y f0 4, associated with the instant of quasi-time-symmetry in the dynamics, and which also follows a quasi-universal relation. We also derive a new quasi-universal relation for the merger frequency f h mer, which provides a description of the data that is 4 times more accurate than previous expressions while requiring fewer fitting coefficients. Finally, consistent with the findings of numerous studies before ours, and using an enlarged ensemble of binary systems, we point out that the ℓ = 2, m = 1 gravitational-wave mode could become comparable with the traditional ℓ = 2, m = 2 mode on sufficiently long timescales, with strain amplitudes in a ratio |h21|/|h22| ∼ 0.1–1 under generic orientations of the binary, which could be measured by present detectors for signals with a large signal-to-noise ratio or by third-generation detectors for generic signals should no collapse occur.
We use holography to study the dynamics of a strongly-coupled gauge theory in four-dimensional de Sitter space with Hubble rate H. The gauge theory is non-conformal with a characteristic mass scale M. We solve Einstein’s equations numerically and determine the time evolution of homogeneous gauge theory states. If their initial energy density is high compared with H4 then the early-time evolution is well described by viscous hydrodynamics with a non-zero bulk viscosity. At late times the dynamics is always far from equilibrium. The asymptotic late-time state preserves the full de Sitter symmetry group and its dual geometry is a domain-wall in AdS5. The approach to this state is characterised by an emergent relation of the form P = w ℰ that is different from the equilibrium equation of state in flat space. The constant w does not depend on the initial conditions but only on H/M and is negative if the ratio H/M is close to unity. The event and the apparent horizons of the late-time solution do not coincide with one another, reflecting its non-equilibrium nature. In between them lies an “entanglement horizon” that cannot be penetrated by extremal surfaces anchored at the boundary, which we use to compute the entanglement entropy of boundary regions. If the entangling region equals the observable universe then the extremal surface coincides with a bulk cosmological horizon that just touches the event horizon, while for larger regions the extremal surface probes behind the event horizon.
One-photon and multi-photon absorption, spontaneous and stimulated photon emission, resonance Raman scattering and electron transfer are important molecular processes that commonly involve combined vibrational-electronic (vibronic) transitions. The corresponding vibronic transition profiles in the energy domain are usually determined by Franck-Condon factors (FCFs), the squared norm of overlap integrals between vibrational wavefunctions of different electronic states. FC profiles are typically highly congested for large molecular systems and the spectra usually become not well-resolvable at elevated temperatures. The (theoretical) analyses of such spectra are even more difficult when vibrational mode mixing (Duschinsky) effects are significant, because contributions from different modes are in general not separable, even within the harmonic approximation. A few decades ago Doktorov, Malkin and Man'ko [1979 J. Mol. Spectrosc. 77, 178] developed a coherent state-based generating function approach and exploited the dynamical symmetry of vibrational Hamiltonians for the Duschinsky relation to describe FC transitions at zero Kelvin. Recently, the present authors extended the method to incorporate thermal, single vibronic level, non-Condon and multi-photon effects in energy, time and probability density domains for the efficient calculation and interpretation of vibronic spectra. Herein, recent developments and corresponding generating functions are presented for single vibronic levels related to fluorescence, resonance Raman scattering and anharmonic transition.
As cryo-EM approaches the physical resolution limits imposed by electron optics and radiation damage, it becomes increasingly urgent to address the issues that impede high-resolution structure determination of biological specimens. One of the persistent problems has been beam-induced movement, which occurs when the specimen is irradiated with high-energy electrons. Beam-induced movement results in image blurring and loss of high-resolution information. It is particularly severe for biological samples in unsupported thin films of vitreous water. By controlled devitrification of conventionally plunge-frozen samples, the suspended film of vitrified water was converted into cubic ice, a polycrystalline, mechanically stable solid. It is shown that compared with vitrified samples, devitrification reduces beam-induced movement in the first 5 e Å−2 of an exposure by a factor of ∼4, substantially enhancing the contribution of the initial, minimally damaged frames to a structure. A 3D apoferritin map reconstructed from the first frames of 20 000 particle images of devitrified samples resolved undamaged side chains. Devitrification of frozen-hydrated specimens helps to overcome beam-induced specimen motion in single-particle cryo-EM, as a further step towards realizing the full potential of cryo-EM for high-resolution structure determination.
CryoEM at IUCRJ: a new era
(2016)
We carry out an in-depth analysis of the prompt-collapse behaviour of binary neutron star (BNS) mergers. To this end, we perform more than 80 general relativistic BNS merger simulations using a family of realistic Equations of State (EOS) with different stiffness, which feature a first order deconfinement phase transition between hadronic and quark matter. From these simulations we infer the critical binary mass Mcrit that separates the prompt from the non-prompt collapse regime. We show that the critical mass increases with the stiffness of the EOS and obeys a tight quasi-universal relation, Mcrit/MTOV ≈ 1.41 ± 0.06, which links it to the maximum mass MTOV of static neutron stars, and therefore provides a straightforward estimate for the total binary mass beyond which prompt collapse becomes inevitable. In addition, we introduce a novel gauge independent definition for a one-parameter family of threshold masses in terms of curvature invariants of the Riemann tensor which characterizes the development toward a more rapid collapse with increasing binary mass. Using these diagnostics, we find that the amount of matter remaining outside the black hole sharply drops in supercritical mass mergers compared to subcritical ones and is further reduced in mergers where the black hole collapse is induced by the formation of a quark matter core. This implies that Mcrit, particularly for merger remnants featuring quark matter cores, imposes a strict upper limit on the emission of any detectable electromagnetic counterpart in BNS mergers.
Post-merger gravitational-wave signal from neutron-star binaries: a new look at an old problem
(2023)
The spectral properties of the post-merger gravitational-wave signal from a binary of neutron stars encodes a variety of information about the features of the system and of the equation of state describing matter around and above nuclear saturation density. Characterising the properties of such a signal is an “old” problem, which first emerged when a number of frequencies were shown to be related to the properties of the binary through “quasi-universal” relations. Here we take a new look at this old problem by computing the properties of the signal in terms of the Weyl scalar ψ4. In this way, and using a database of more than 100 simulations, we provide the first evidence for a new instantaneous frequency, f ψ4 0, associated with the instant of quasi timesymmetry in the postmerger dynamics, and which also follows a quasi-universal relation. We also derive a new quasi-universal relation for the merger frequency f h mer, which provides a description of the data that is four times more accurate than previous expressions while requiring fewer fitting coefficients. Finally, consistently with the findings of numerous studies before ours, and using an enlarged ensamble of binary systems we point out that the ℓ = 2, m = 1 gravitational-wave mode could become comparable with the traditional ℓ = 2, m = 2 mode on sufficiently long timescales, with strain amplitudes in a ratio |h 21|/|h 22| ∼ 0.1 − 1 under generic orientations of the binary, which could be measured by present detectors for signals with large signal-to-noise ratio or by third-generation detectors for generic signals should no collapse occur.
A considerable effort has been dedicated recently to the construction of generic equations of state (EOSs) for matter in neutron stars. The advantage of these approaches is that they can provide model-independent information on the interior structure and global properties of neutron stars. Making use of more than 106 generic EOSs, we asses the validity of quasi-universal relations of neutron star properties for a broad range of rotation rates, from slow-rotation up to the mass-shedding limit. In this way, we are able to determine with unprecedented accuracy the quasi-universal maximum-mass ratio between rotating and nonrotating stars and reveal the existence of a new relation for the surface oblateness, i.e., the ratio between the polar and equatorial proper radii. We discuss the impact that our findings have on the imminent detection of new binary neutron-star mergers and how they can be used to set new and more stringent limits on the maximum mass of nonrotating neutron stars, as well as to improve the modelling of the X-ray emission from the surface of rotating stars.
We have investigated the systematic differences introduced when performing a Bayesian-inference analysis of the equation of state of neutron stars employing either variable- or constant-likelihood functions. The former have the advantage that it retains the full information on the distributions of the measurements, making an exhaustive usage of the data. The latter, on the other hand, have the advantage of a much simpler implementation and reduced computational costs. In both approaches, the EOSs have identical priors and have been built using the sound-speed parameterization method so as to satisfy the constraints from X-ray and gravitationalwaves observations, as well as those from Chiral Effective Theory and perturbative QCD. In all cases, the two approaches lead to very similar results and the 90%-confidence levels are essentially overlapping. Some differences do appear, but in regions where the probability density is extremely small and are mostly due to the sharp cutoff set on the binary tidal deformability Λ˜ ≤ 720 employed in the constant-likelihood analysis. Our analysis has also produced two additional results. First, a clear inverse correlation between the normalized central number density of a maximally massive star, nc,TOV/ns, and the radius of a maximally massive star, RTOV. Second, and most importantly, it has confirmed the relation between the chirp mass Mchirp and the binary tidal deformability Λ˜. The importance of this result is that it relates a quantity that is measured very accurately, Mchirp, with a quantity that contains important information on the micro-physics, Λ˜. Hence, once Mchirp is measured in future detections, our relation has the potential of setting tight constraints on Λ˜.
We use the quantum null energy condition in strongly coupled two-dimensional field theories (QNEC2) as diagnostic tool to study a variety of phase structures, including crossover, second and first order phase transitions. We find a universal QNEC2 constraint for first order phase transitions with kinked entanglement entropy and discuss in general the relation between the QNEC2-inequality and monotonicity of the Casini-Huerta c-function. We then focus on a specific example, the holographic dual of which is modelled by three-dimensional Einstein gravity plus a massive scalar field with one free parameter in the self-interaction potential. We study translation invariant stationary states dual to domain walls and black branes. Depending on the value of the free parameter we find crossover, second and first order phase transitions between such states, and the c-function either flows to zero or to a finite value in the infrared. We present evidence that evaluating QNEC2 for ground state solutions allows to predict the existence of phase transitions at finite temperature.
Holography has provided valuable insights into the time evolution of strongly coupled gauge theories in a fixed spacetime. However, this framework is insufficient if this spacetime is dynamical. We present a scheme to evolve a four-dimensional, strongly interacting gauge theory coupled to four-dimensional dynamical gravity in the semiclassical regime. As in previous work, we use holography to evolve the quantum gauge theory stress tensor, whereas the four-dimensional metric evolves according to Einstein’s equations coupled to the expectation value of the stress tensor. The novelty of our approach is that both the boundary and the bulk spacetimes are constructed dynamically, one time step at a time. We focus on Friedmann-Lemaître-Robertson-Walker geometries and evolve far-from-equilibrium initial states that lead to asymptotically expanding, flat or collapsing Universes.
According to the inflationary theory of cosmology, most elementary particles in the current Universe were created during a period of reheating after inflation. In this Letter, we self-consistently couple the Einstein-inflaton equations to a strongly coupled quantum field theory as described by holography. We show that this leads to an inflating universe, a reheating phase, and finally a universe dominated by the quantum field theory in thermal equilibrium.
Using full 3+1 dimensional general-relativistic hydrodynamic simulations of equal- and unequal-mass neutron-star binaries with properties that are consistent with those inferred from the inspiral of GW170817, we perform a detailed study of the quark-formation processes that could take place after merger. We use three equations of state consistent with current pulsar observations derived from a novel finite-temperature framework based on V-QCD, a non-perturbative gauge/gravity model for Quantum Chromodynamics. In this way, we identify three different post-merger stages at which mixed baryonic and quark matter, as well as pure quark matter, are generated. A phase transition triggered collapse already ≲10ms after the merger reveals that the softest version of our equations of state is actually inconsistent with the expected second-long post-merger lifetime of GW170817. Our results underline the impact that multi-messenger observations of binary neutron-star mergers can have in constraining the equation of state of nuclear matter, especially in its most extreme regimes.
Using full 3+1 dimensional general-relativistic hydrodynamic simulations of equal- and unequal-mass neutron-star binaries with properties that are consistent with those inferred from the inspiral of GW170817, we perform a detailed study of the quark-formation processes that could take place after merger. We use three equations of state consistent with current pulsar observations derived from a novel finite-temperature framework based on V-QCD, a non-perturbative gauge/gravity model for Quantum Chromodynamics. In this way, we identify three different post-merger stages at which mixed baryonic and quark matter, as well as pure quark matter, are generated. A phase transition triggered collapse already ≲ 10 ms after the merger reveals that the softest version of our equations of state is actually inconsistent with the expected second-long post-merger lifetime of GW170817. Our results underline the impact that gravitational wave observations of binary neutron-star mergers can have in constraining the equation of state of nuclear matter, especially in its most extreme regimes.
We present a novel framework for the equation of state of dense and hot quantum chromodynamics (QCD), which focuses on the region of the phase diagram relevant for neutron star mergers and core-collapse supernovae. The model combines predictions from the gauge/gravity duality with input from lattice field theory, QCD perturbation theory, chiral effective theory, and statistical modeling. It is therefore, by construction, in good agreement with theoretical constraints both at low and high densities and temperatures. The main ingredients of our setup are the nonperturbative V-QCD model based on the gauge/gravity duality, a van der Waals model for nucleon liquid, and the DD2 version of the Hempel-Schaffner-Bielich statistical model of nuclear matter. By consistently combining these models, we also obtain a description for the nuclear to quark matter phase transition and its critical end point. The parameter dependence of the model is represented by three (soft, intermediate, and stiff) variants of the equation of state, all of which agree with observational constraints from neutron stars and their mergers. We discuss resulting constraints for the equation of state, predictions for neutron stars, and the location of the critical point.
Holography has provided valuable insights into the time evolution of strongly coupled gauge theories in a fixed spacetime. However, this framework is insufficient if this spacetime is dynamical. We present a novel scheme to evolve a four-dimensional, strongly interacting gauge theory coupled to four-dimensional dynamical gravity in the semiclassical regime. We use holography to evolve the quantum gauge theory stress tensor. The four-dimensional metric evolves according to the four-dimensional Einstein equations coupled to the expectation value of the stress tensor. We focus on Friedmann-Lemaître-Robertson-Walker geometries and evolve far-from-equilibrium initial states that lead to asymptotically expanding, flat or collapsing Universes.
We use the quantum null energy condition in strongly coupled two-dimensional field theories (QNEC2) as diagnostic tool to study a variety of phase structures, including crossover, second and first order phase transitions. We find a universal QNEC2 constraint for first order phase transitions with kinked entanglement entropy and discuss in general the relation between the QNEC2-inequality and monotonicity of the Casini-Huerta c-function. We then focus on a specific example, the holographic dual of which is modelled by three-dimensional Einstein gravity plus a massive scalar field with one free parameter in the self-interaction potential. We study translation invariant stationary states dual to domain walls and black branes. Depending on the value of the free parameter we find crossover, second and first order phase transitions between such states, and the c-function either flows to zero or to a finite value in the infrared. Strikingly, evaluating QNEC2 for ground state solutions allows to predict the existence of phase transitions at finite temperature.
We use the quantum null energy condition in strongly coupled two-dimensional field theories (QNEC2) as diagnostic tool to study a variety of phase structures, including crossover, second and first order phase transitions. We find a universal QNEC2 constraint for first order phase transitions with kinked entanglement entropy and discuss in general the relation between the QNEC2-inequality and monotonicity of the Casini-Huerta c-function. We then focus on a specific example, the holographic dual of which is modelled by three-dimensional Einstein gravity plus a massive scalar field with one free parameter in the self-interaction potential. We study translation invariant stationary states dual to domain walls and black branes. Depending on the value of the free parameter we find crossover, second and first order phase transitions between such states, and the c-function either flows to zero or to a finite value in the infrared. Strikingly, evaluating QNEC2 for ground state solutions allows to predict the existence of phase transitions at finite temperature.
We use the quantum null energy condition in strongly coupled two-dimensional field theories (QNEC2) as diagnostic tool to study a variety of phase structures, including crossover, second and first order phase transitions. We find a universal QNEC2 constraint for first order phase transitions with kinked entanglement entropy and discuss in general the relation between the QNEC2-inequality and monotonicity of the Casini-Huerta c-function. We then focus on a specific example, the holographic dual of which is modelled by three-dimensional Einstein gravity plus a massive scalar field with one free parameter in the self-interaction potential. We study translation invariant stationary states dual to domain walls and black branes. Depending on the value of the free parameter we find crossover, second and first order phase transitions between such states, and the c-function either flows to zero or to a finite value in the infrared. Strikingly, evaluating QNEC2 for ground state solutions allows to predict the existence of phase transitions at finite temperature.
Holography has provided valuable insights into the time evolution of strongly coupled gauge theories in a fixed spacetime. However, this framework is insufficient if this spacetime is dynamical. We present a scheme to evolve a four-dimensional, strongly interacting gauge theory coupled to four-dimensional dynamical gravity in the semiclassical regime. As in previous work, we use holography to evolve the quantum gauge theory stress tensor, whereas the four-dimensional metric evolves according to Einstein's equations coupled to the expectation value of the stress tensor. The novelty of our approach is that both the boundary and the bulk spacetimes are constructed dynamically, one time step at a time. We focus on Friedmann-Lemaître-Robertson-Walker geometries and evolve far-from-equilibrium initial states that lead to asymptotically expanding, flat or collapsing Universes.
We present a novel framework for the equation of state of dense and hot Quantum Chromodynamics (QCD), which focuses on the region of the phase diagram relevant for neutron star mergers and core-collapse supernovae. The model combines predictions from the gauge/gravity duality with input from lattice field theory, QCD perturbation theory, chiral effective theory and statistical modeling. It is therefore, by construction, in good agreement with theoretical constraints both at low and high densities and temperatures. The main ingredients of our setup are the non-perturbative V-QCD model based on the gauge/gravity duality, a van der Waals model for nucleon liquid, and the DD2 version of the Hempel-Schaffner-Bielich statistical model of nuclear matter. By consistently combining these models, we also obtain a description for the nuclear to quark matter phase transition and its critical endpoint. The parameter dependence of the model is represented by three (soft, intermediate and stiff) variants of the equation of state, all of which agree with observational constraints from neutron stars and their mergers. We discuss resulting constraints for the equation of state, predictions for neutron stars and the location of the critical point.
According to the inflationary theory of cosmology, most elementary particles in the current universe were created during a period of reheating after inflation. In this work we self-consistently couple the Einstein-inflaton equations to a strongly coupled quantum field theory (QFT) as described by holography. We show that this leads to an inflating universe, a reheating phase and finally a universe dominated by the QFT in thermal equilibrium.
We use holography to study the dynamics of a strongly-coupled gauge theory in four-dimensional de Sitter space with Hubble rate H. The gauge theory is non-conformal with a characteristic mass scale M. We solve Einstein’s equations numerically and determine the time evolution of homogeneous gauge theory states. If their initial energy density is high compared with H4 then the early-time evolution is well described by viscous hydrodynamics with a non-zero bulk viscosity. At late times the dynamics is always far from equilibrium. The asymptotic late-time state preserves the full de Sitter symmetry group and its dual geometry is a domain-wall in AdS5. The approach to this state is characterised by an emergent relation of the form P = w E that is different from the equilibrium equation of state in flat space. The constant w does not depend on the initial conditions but only on H/M and is negative if the ratio H/M is close to unity. The event and the apparent horizons of the late-time solution do not coincide with one another, reflecting its non-equilibrium nature. In between them lies an “entanglement horizon” that cannot be penetrated by extremal surfaces anchored at the boundary, which we use to compute the entanglement entropy of boundary regions. If the entangling region equals the observable universe then the extremal surface coincides with a bulk cosmological horizon that just touches the event horizon, while for larger regions the extremal surface probes behind the event horizon.
We present the first holographic simulations of non-equilibrium steady state formation in strongly coupled N=4 SYM theory in 3+1 dimensions. We initially join together two thermal baths at different temperatures and chemical potentials and compare the subsequent evolution of the combined system to analytic solutions of the corresponding Riemann problem and to numeric solutions of ideal and viscous hydrodynamics. The time evolution of the energy density that we obtain holographically is consistent with the combination of a shock and a rarefaction wave: A shock wave moves towards the cold bath, and a smooth broadening wave towards the hot bath. Between the two waves emerges a steady state with constant temperature and flow velocity, both of which are accurately described by a shock+rarefaction wave solution of the Riemann problem. In the steady state region, a smooth crossover develops between two regions of different charge density. This is reminiscent of a contact discontinuity in the Riemann problem. We also obtain results for the entanglement entropy of regions crossed by shock and rarefaction waves and find both of them to closely follow the evolution of the energy density.
We present the first holographic simulations of non-equilibrium steady state formation in strongly coupled N=4 SYM theory in 3+1 dimensions. We initially join together two thermal baths at different temperatures and chemical potentials and compare the subsequent evolution of the combined system to analytic solutions of the corresponding Riemann problem and to numeric solutions of ideal and viscous hydrodynamics. The time evolution of the energy density that we obtain holographically is consistent with the combination of a shock and a rarefaction wave: A shock wave moves towards the cold bath, and a smooth broader wave towards the hot bath. Between the two waves emerges a steady state with constant temperature and flow velocity, both of which are accurately described by a shock+rarefaction wave solution of the Riemann problem. In the steady state region develops a smooth crossover between two regions of different charge densities that diffuses on a timescale proportional to t√ and is reminiscent of a contact discontinuity in the Riemann problem. We also obtain results for the entanglement entropy of regions crossed by shock and rarefaction waves and find both of them to closely follow the evolution of the energy density.
Biological as well as advanced artificial intelligences (AIs) need to decide which goals to pursue. We review nature's solution to the time allocation problem, which is based on a continuously readjusted categorical weighting mechanism we experience introspectively as emotions. One observes phylogenetically that the available number of emotional states increases hand in hand with the cognitive capabilities of animals and that raising levels of intelligence entail ever larger sets of behavioral options. Our ability to experience a multitude of potentially conflicting feelings is in this view not a leftover of a more primitive heritage, but a generic mechanism for attributing values to behavioral options that can not be specified at birth. In this view, emotions are essential for understanding the mind. For concreteness, we propose and discuss a framework which mimics emotions on a functional level. Based on time allocation via emotional stationarity (TAES), emotions are implemented as abstract criteria, such as satisfaction, challenge and boredom, which serve to evaluate activities that have been carried out. The resulting timeline of experienced emotions is compared with the “character” of the agent, which is defined in terms of a preferred distribution of emotional states. The long-term goal of the agent, to align experience with character, is achieved by optimizing the frequency for selecting individual tasks. Upon optimization, the statistics of emotion experience becomes stationary.
The core of this work is represented by the investigation of the chiral phase transition, using Monte Carlo simulations and unimproved staggered fermions, both in the weak and strong coupling regimes of Quantum Chromodynamics. Based on recent results from Monte Carlo simulations, both using unimproved staggered fermions and Wilson fermions, the chiral phase transition in the continuum and chiral limit shows compatibility with a second-order phase transition for Nf (number of flavours) in range [2:7], at zero baryon chemical potential. This achievement relies on the analytic continuation of Nf to non-integer values on the lattice, which allows to make use of extrapolation techniques to the chiral limit, where simulations are not possible. Furthermore, these results provide a resolution to the ambiguous scenario for Nf = 2 in the chiral limt. The first part of this thesis is devoted to the investigation of the chiral phase transition when a non-zero imaginary baryon chemical potential is involved, whose value corresponds to the 81% of the Roberge-Weiss one. Using the same extrapolation techniques aforementioned, the order of the chiral phase transition in the continuum and chiral limit shows compatibility with a second-order phase transition for Nf in range [2:6], highlighting a lack of dependence of the order of the chiral phase transition on the imaginary baryon chemical potential value. The second part of this thesis is about the study of the extension of the first-order chiral region in the strong coupling regime, at zero baryon chemical potential. Using Monte Carlo techniques, this can be done by investigating the Z2 boundary on a coarse lattice, whose temporal extent reads Nt = 2, and simulations are realised for Nf = 4, 8. The results in the weak coupling regime show, for $Nt = 8, 6, 4 and fixed Nf value, an inflating first-order chiral region. As in the strong coupling limit a second-order chiral phase transition is expected, the first-order chiral region has to shrink as the strong coupling regime is approached, resulting in a non-monotonic behaviour of the Z2 boundary. For Nf = 8, a critical mass on the Z2 boundary has been obtained, confirming the expected non-monotonic behaviour. For Nf = 4 the results do not provide a unique conclusion: Either a Z2 boundary at extremely low bare quark mass or a second-order chiral phase transition in the O(2) universality class in the chiral limit can take place. In addition to the two main topics, the performances of the second-order minimum norm integrator (2MN) and the fourth-order minimum norm integrator (4MN) have been compared, after implementing the 4MN one in the CL2QCD code used to realise our simulations. The 2MN integrator had already been implemented in the code since the first version was released. The two integrators belong to the class of symplectic integrators and represent an essential component of the RHMC algorithm, involved in our investigation. This step is extremely important, in order to guarantee the best quality when collecting data from simulations, and the results of the comparison suggested to favor the 2MN integrator, for both the topics.
A powerful technique to distinguish the enantiomers of a chiral molecule is the Coulomb Explosion Imaging (CEI). This technique allows us to determine the handedness of a single molecule. In CEI, the molecule becomes charged by losing many electrons in a very short period of time by interacting with the light. The repulsion forces between the positive charged particles of the molecule leads the molecule to break into parts-fragments. By measuring the three vector momentum of (at least) four fragments, the handedness observable can be determined. In this thesis, CEI is induced by absorption of a single high energy photon, which creates an inner-shell hole (K shell) of the molecule. The subsequent cascade of Auger decays lead to fragmentation. We decided to work with the formic acid molecule in this thesis. Two different experiments were conducted. The first experiment focused on exciting electrons to different energy states, while the second experiment focused on extracting directly a photoelectron to the continuum and measure the angular distribution of the photoelectron in the molecular frame. The primary goal was to search for chiral signal in a pure achiral planar molecule under the previous electron processes. The results of these findings were further implemented to two more molecules.
In the framework of the LHC Injectors Upgrade Project (LIU), the CERN Proton Synchrotron Booster (PSB) went through major upgrades resulting in new effects to study, challenges to overcome and new parameter regimes to explore. To assess the achievable beam brightness limit of the machine, a series of experimental and computational studies in the transverse planes were performed. In particular, the new injection scheme induces optics perturbations that are strongly enhanced near the half-integer resonance. In this thesis, methods for dynamically measuring and correcting these perturbations and their impact on the beam performance will be presented. Additionally, the quality of the transverse beam distributions and strategies for improvement will be addressed. Finally, the space charge effects when dynamically crossing the half-integer resonance will be characterized. The results of these studies and their broader significance beyond the PSB will be discussed.
Measurements of the pT-dependent flow vector fluctuations in Pb-Pb collisions at sNN−−−√=5.02 TeV using azimuthal correlations with the ALICE experiment at the LHC are presented. A four-particle correlation approach [1] is used to quantify the effects of flow angle and magnitude fluctuations separately. This paper extends previous studies to additional centrality intervals and provides measurements of the pT-dependent flow vector fluctuations at sNN−−−√=5.02 TeV with two-particle correlations. Significant pT-dependent fluctuations of the V⃗ 2 flow vector in Pb-Pb collisions are found across different centrality ranges, with the largest fluctuations of up to ∼15% being present in the 5% most central collisions. In parallel, no evidence of significant pT-dependent fluctuations of V⃗ 3 or V⃗ 4 is found. Additionally, evidence of flow angle and magnitude fluctuations is observed with more than 5σ significance in central collisions. These observations in Pb-Pb collisions indicate where the classical picture of hydrodynamic modeling with a common symmetry plane breaks down. This has implications for hard probes at high pT, which might be biased by pT-dependent flow angle fluctuations of at least 23% in central collisions. Given the presented results, existing theoretical models should be re-examined to improve our understanding of initial conditions, quark--gluon plasma (QGP) properties, and the dynamic evolution of the created system.
The intense photon fluxes from relativistic nuclei provide an opportunity to study photonuclear interactions in ultraperipheral collisions. The measurement of coherently photoproduced π+π−π+π− final states in ultraperipheral Pb-Pb collisions at sNN−−−√=5.02 TeV is presented for the first time. The cross section, dσ/dy, times the branching ratio (ρ→π+π+π−π−) is found to be 47.8±2.3 (stat.)±7.7 (syst.) mb in the rapidity interval |y|<0.5. The invariant mass distribution is not well described with a single Breit-Wigner resonance. The production of two interfering resonances, ρ(1450) and ρ(1700), provides a good description of the data. The values of the masses (m) and widths (Γ) of the resonances extracted from the fit are m1=1385±14 (stat.)±3 (syst.) MeV/c2, Γ1=431±36 (stat.)±82 (syst.) MeV/c2, m2=1663±13 (stat.)±22 (syst.) MeV/c2 and Γ2=357±31 (stat.)±49 (syst.) MeV/c2, respectively. The measured cross sections times the branching ratios are compared to recent theoretical predictions.
Measurement of beauty-quark production in pp collisions at √s = 13 TeV via non-prompt D mesons
(2024)
The pT-differential production cross sections of non-prompt D0, D+, and D+s mesons originating from beauty-hadron decays are measured in proton−proton collisions at a centre-of-mass energy s√ of 13 TeV. The measurements are performed at midrapidity, |y|<0.5, with the data sample collected by ALICE from 2016 to 2018. The results are in agreement with predictions from several perturbative QCD calculations. The fragmentation fraction of beauty quarks to strange mesons divided by the one to non-strange mesons, fs/(fu+fd), is found to be 0.114±0.016 (stat.)±0.006 (syst.)±0.003 (BR)±0.003 (extrap.). This value is compatible with previous measurements at lower centre-of-mass energies and in different collision systems in agreement with the assumption of universality of fragmentation functions. In addition, the dependence of the non-prompt D meson production on the centre-of-mass energy is investigated by comparing the results obtained at s√=5.02 and 13 TeV, showing a hardening of the non-prompt D-meson pT-differential production cross section at higher s√. Finally, the bb¯¯¯ production cross section per unit of rapidity at midrapidity is calculated from the non-prompt D0, D+, D+s, and Λ+c hadron measurements, obtaining dσ/dy=75.2±3.2 (stat.)±5.2 (syst.)+12.3−3.2 (extrap.) μb.
The two-particle momentum correlation functions between charm mesons (D∗± and D±) and charged light-flavor mesons (π± and K±) in all charge-combinations are measured for the first time by the ALICE Collaboration in high-multiplicity proton–proton collisions at a center-of-mass energy of √s = 13 TeV. For DK and D∗K pairs, the experimental results are in agreement with theoretical predictions of the residual strong interaction based on quantum chromodynamics calculations on the lattice and chiral effective field theory. In the case of Dπ and D∗π pairs, tension between the calculations including strong interactions and the measurement is observed. For all particle pairs, the data can be adequately described by Coulomb interaction only, indicating a shallow interaction between charm and light-flavor mesons. Finally, the scattering lengths governing the residual strong interaction of the Dπ and D∗π systems are determined by fitting the experimental correlation functions with a model that employs a Gaussian potential. The extracted values are small and compatible with zero.
The influence of an ac current of arbitrary amplitude and frequency on the mixed-state dc-voltage-ac-drive tiltingratchet response of a superconducting film with uniaxial cosine pinning potential at finite temperature is theoretically investigated. The results are obtained in the single-vortex approximation, within the frame of an exact solution of the Langevin equation for non-interacting vortices. Both experimentally achievable, the dc ratchet response and absorbed ac power are predicted to demonstrate a pronounced filter-like behavior at microwave frequencies. Based on our findings, we propose a cut-off filter and discuss its operating curves as functions of the driving parameters, i.e, ac amplitude, frequency, and dc bias. The predicted results can be examined, e.g, on superconducting films with a washboard pinning potential landscape.
LICE is one of the four major LHC experiments at CERN. When the accelerator enters the Run 3 data-taking period, starting in 2021, ALICE expects almost 100 times more Pb-Pb central collisions than now, resulting in a large increase of data throughput. In order to cope with this new challenge, the collaboration had to extensively rethink the whole data processing chain, with a tighter integration between Online and Offline computing worlds. Such a system, code-named ALICE O2, is being developed in collaboration with the FAIR experiments at GSI. It is based on the ALFA framework which provides a generalized implementation of the ALICE High Level Trigger approach, designed around distributed software entities coordinating and communicating via message passing.
We will highlight our efforts to integrate ALFA within the ALICE O2 environment. We analyze the challenges arising from the different running environments for production and development, and conclude on requirements for a flexible and modular software framework. In particular we will present the ALICE O2 Data Processing Layer which deals with ALICE specific requirements in terms of Data Model. The main goal is to reduce the complexity of development of algorithms and managing a distributed system, and by that leading to a significant simplification for the large majority of the ALICE users.
This thesis provides a detailed derivation of dissipative spin hydrodynamics from quantum field theory for systems composed of spin-0, spin-1/2, or spin-1 particles.
The Wigner function formalism is introduced for quantum fields in the respective representations of the Poincaré group, and the conserved currents, i.e., the energy-momentum tensor and the total angular momentum tensor, in various so-called pseudogauges are derived. An expansion around the semiclassical limit in powers of the Planck constant is performed.
Subsequently, kinetic equations are obtained for binary elastic scattering, using both the de Groot-van Leeuwen-van Weert and Kadanoff-Baym method, with the latter retaining the effect of quantum statistics. The resulting collision term features both local and nonlocal contributions, with the latter providing a relaxation mechanism for the spin degrees of freedom of the quasiparticles. The local-equilibrium distribution function is derived from the requirement that the local part of the collision term vanishes.
From quantum kinetic theory, dissipative spin hydrodynamics is then constructed via the method of moments, extended to particles with spin. The system of moment equations is closed via the Inverse-Reynolds Dominance (IReD) approach, resulting in a set of equations of motion describing the evolution of both ideal and dissipative degrees of freedom. The application to polarization phenomena relevant to heavy-ion collisions is discussed.
Hadron lists based on experimental studies summarized by the Particle Data Group (PDG) are a crucial input for the equation of state and thermal models used in the study of strongly-interacting matter produced in heavy-ion collisions. Modeling of these strongly-interacting systems is carried out via hydrodynamical simulations, which are followed by hadronic transport codes that also require a hadronic list as input. To remain consistent throughout the different stages of modeling of a heavy-ion collision, the same hadron list with its corresponding decays must be used at each step. It has been shown that even the most uncertain states listed in the PDG from 2016 are required to reproduce partial pressures and susceptibilities from Lattice Quantum Chromodynamics with the hadronic list known as the PDG2016+. Here, we update the hadronic list for use in heavy-ion collision modeling by including the latest experimental information for all states listed in the Particle Data Booklet in 2021. We then compare our new list, called PDG2021+, to Lattice Quantum Chromodynamics results and find that it achieves even better agreement with the first principles calculations than the PDG2016+ list. Furthermore, we develop a novel scheme based on intermediate decay channels that allows for only binary decays, such that PDG2021+ will be compatible with the hadronic transport framework SMASH. Finally, we use these results to make comparisons to experimental data and discuss the impact on particle yields and spectra.
This paper reports on Monte Carlo simulation results for future measurements of the moduli of time-like proton electromagnetic form factors, |GE | and |GM|, using the ¯pp → μ+μ− reaction at PANDA (FAIR). The electromagnetic form factors are fundamental quantities parameterizing the electric and magnetic structure of hadrons. This work estimates the statistical and total accuracy with which the form factors can be measured at PANDA, using an analysis of simulated data within the PandaRoot software framework. The most crucial background channel is ¯pp → π+π−,due to the very similar behavior of muons and pions in the detector. The suppression factors are evaluated for this and all other relevant background channels at different values of antiproton beam momentum. The signal/background separation is based on a multivariate analysis, using the Boosted Decision Trees method. An expected background subtraction is included in this study, based on realistic angular distribuations of the background contribution. Systematic uncertainties are considered and the relative total uncertainties of the form factor measurements are presented.
Fossil dental remains are an archive of unique information for paleobiological studies. Computed microtomography based on X-ray microfocus sources (X-μCT) and Synchrotron Radiation (SR-μCT) allow subtle quantification at the micron and sub-micron scale of the meso- and microstructural signature imprinted in the mineralized tissues, such as enamel and dentine, through high-resolution “virtual histology”. Nonetheless, depending on the degree of alterations undergone during fossilization, X-ray analyses of tooth tissues do not always provide distinct imaging contrasts, thus preventing the extraction of essential morphological and anatomical details. We illustrate here by three examples the successful application of neutron microtomography (n-μCT) in cases where X-rays have previously failed to deliver contrasts between dental tissues of fossilized specimen.
How long does it take to emit an electron from an atom? This question has intrigued scientists for decades. As such emission times are in the attosecond regime, the advent of attosecond metrology using ultrashort and intense lasers has re-triggered strong interest on the topic from an experimental standpoint. Here, we present an approach to measure such emission delays, which does not require attosecond light pulses, and works without the presence of superimposed infrared laser fields. We instead extract the emission delay from the interference pattern generated as the emitted photoelectron is diffracted by the parent ion’s potential. Targeting core electrons in CO, we measured a 2d map of photoelectron emission delays in the molecular frame over a wide range of electron energies. The emission times depend drastically on the photoelectrons’ emission directions in the molecular frame and exhibit characteristic changes along the shape resonance of the molecule.
Nuclear magnetic resonance (NMR) spectroscopy is a powerful and popular technique for probing the molecular structures, dynamics and chemical properties. However the conventional NMR spectroscopy is bottlenecked by its low sensitivity. Dynamic nuclear polarization (DNP) boosts NMR sensitivity by orders of magnitude and resolves this limitation. In liquid-state this revolutionizing technique has been restricted to a few specific non-biological model molecules in organic solvents. Here we show that the carbon polarization in small biological molecules, including carbohydrates and amino acids, can be enhanced sizably by in situ Overhauser DNP (ODNP) in water at room temperature and at high magnetic field. An observed connection between ODNP 13C enhancement factor and paramagnetic 13C NMR shift has led to the exploration of biologically relevant heterocyclic compound indole. The QM/MM MD simulation underscores the dynamics of intermolecular hydrogen bonds as the driving force for the scalar ODNP in a long-living radical-substrate complex. Our work reconciles results obtained by DNP spectroscopy, paramagnetic NMR and computational chemistry and provides new mechanistic insights into the high-field scalar ODNP.
Upon antibiotic stress Gram-negative pathogens deploy resistance-nodulation-cell division-type tripartite efflux pumps. These include a H+/drug antiporter module that recognizes structurally diverse substances, including antibiotics. Here, we show the 3.5 Å structure of subunit AdeB from the Acinetobacter baumannii AdeABC efflux pump solved by single-particle cryo-electron microscopy. The AdeB trimer adopts mainly a resting state with all protomers in a conformation devoid of transport channels or antibiotic binding sites. However, 10% of the protomers adopt a state where three transport channels lead to the closed substrate (deep) binding pocket. A comparison between drug binding of AdeB and Escherichia coli AcrB is made via activity analysis of 20 AdeB variants, selected on basis of side chain interactions with antibiotics observed in the AcrB periplasmic domain X-ray co-structures with fusidic acid (2.3 Å), doxycycline (2.1 Å) and levofloxacin (2.7 Å). AdeABC, compared to AcrAB-TolC, confers higher resistance to E. coli towards polyaromatic compounds and lower resistance towards antibiotic compounds.
Single-particle tracking enables the analysis of the dynamics of biomolecules in living cells with nanometer spatial and millisecond temporal resolution. This technique reports on the mobility of membrane proteins and is sensitive to the molecular state of a biomolecule and to interactions with other biomolecules. Trajectories describe the mobility of single particles over time and provide information such as the diffusion coefficient and diffusion state. Changes in particle dynamics within single trajectories lead to segmentation, which allows to extract information on transitions of functional states of a biomolecule. Here, mean-squared displacement analysis is developed to classify trajectory segments into immobile, confined diffusing, and freely diffusing states, and to extract the occurrence of transitions between these modes. We applied this analysis to single-particle tracking data of the membrane receptor MET in live cells and analyzed state transitions in single trajectories of the un-activated receptor and the receptor bound to the ligand internalin B. We found that internalin B-bound MET shows an enhancement of transitions from freely and confined diffusing states into the immobile state as compared to un-activated MET. Confined diffusion acts as an intermediate state between immobile and free, as this state is most likely to change the diffusion state in the following segment. This analysis can be readily applied to single-particle tracking data of other membrane receptors and intracellular proteins under various conditions and contribute to the understanding of molecular states and signaling pathways.
Gradient-consistent enrichment of finite element spaces for the DNS of fluid-particle interaction
(2019)
Highlights
• Monolithic scheme for particulate flows preventing an oscillating pressure along the interface.
• The choice of enriching shape functions is driven by the properties of its gradient instead of its value.
• The choice of enriching shape functions inherits a natural stabilization on small cut elements.
Abstract
We present gradient-consistent enriched finite element spaces for the simulation of free particles in a fluid. This involves forces being exchanged between the particles and the fluid at the interface. In an earlier work [23] we derived a monolithic scheme which includes the interaction forces into the Navier-Stokes equations by means of a fictitious domain like strategy. Due to an inexact approximation of the interface oscillations of the pressure along the interface were observed. In multiphase flows oscillations and spurious velocities are a common issue. The surface force term yields a jump in the pressure and therefore the oscillations are usually resolved by extending the spaces on cut elements in order to resolve the discontinuity. For the construction of the enriched spaces proposed in this paper we exploit the Petrov-Galerkin formulation of the vertex-centered finite volume method (PG-FVM), as already investigated in [23]. From the perspective of the finite volume scheme we argue that wrong discrete normal directions at the interface are the origin of the oscillations. The new perspective of normal vectors suggests to look at gradients rather than values of the enriching shape functions. The crucial parameter of the enrichment functions therefore is the gradient of the shape functions and especially the one of the test space. The distinguishing feature of our construction therefore is an enrichment that is based on the choice of shape functions with consistent gradients. These derivations finally yield a fitted scheme for the immersed interface. We further propose a strategy ensuring a well-conditioned system independent of the location of the interface. The enriched spaces can be used within any existing finite element discretization for the Navier-Stokes equation. Our numerical tests were conducted using the PG-FVM. We demonstrate that the enriched spaces are able to eliminate the oscillations.
Rotational test spaces for a fully-implicit FVM and FEM for the DNS of fluid-particle interaction
(2019)
The paper presents a fully-implicit and stable finite element and finite volume scheme for the simulation of freely moving particles in a fluid. The developed method is based on the Petrov-Galerkin formulation of a vertex-centered finite volume method (PG-FVM) on unstructured grids. Appropriate extension of the ansatz and test spaces lead to a formulation comparable to a fictitious domain formulation. The purpose of this work is to introduce a new concept of numerical modeling reducing the mathematical overhead which many other methods require. It exploits the identification of the PG-FVM with a corresponding finite element bilinear form. The surface integrals of the finite volume scheme enable a natural incorporation of the interface forces purely based on the original bilinear operator for the fluid. As a result, there is no need to expand the system of equations to a saddle-point problem. Like for fictitious domain methods the extended scheme treats the particles as rigid parts of the fluid. The distinguishing feature compared to most existing fictitious domain methods is that there is no need for an additional Lagrange multiplier or other artificial external forces for the fluid-solid coupling. Consequently, only one single solve for the derived linear system for the fluid together with the particles is necessary and the proposed method does not require any fractional time stepping scheme to balance the interaction forces between fluid and particles. For the linear Stokes problem we will prove the stability of both schemes. Moreover, for the stationary case the conservation of mass and momentum is not violated by the extended scheme, i.e. conservativity is accomplished within the range of the underlying, unconstrained discretization scheme. The scheme is applicable for problems in two and three dimensions.
We investigate the applicability of the well-known multilevel Monte Carlo (MLMC) method to the class of density-driven flow problems, in particular the problem of salinisation of coastal aquifers. As a test case, we solve the uncertain Henry saltwater intrusion problem. Unknown porosity, permeability and recharge parameters are modelled by using random fields. The classical deterministic Henry problem is non-linear and time-dependent, and can easily take several hours of computing time. Uncertain settings require the solution of multiple realisations of the deterministic problem, and the total computational cost increases drastically. Instead of computing of hundreds random realisations, typically the mean value and the variance are computed. The standard methods such as the Monte Carlo or surrogate-based methods are a good choice, but they compute all stochastic realisations on the same, often, very fine mesh. They also do not balance the stochastic and discretisation errors. These facts motivated us to apply the MLMC method. We demonstrate that by solving the Henry problem on multi-level spatial and temporal meshes, the MLMC method reduces the overall computational and storage costs. To reduce the computing cost further, parallelization is performed in both physical and stochastic spaces. To solve each deterministic scenario, we run the parallel multigrid solver ug4 in a black-box fashion.
Collective flow phenomena are a sensitive probe for the properties of extreme QCD matter. However, their interpretation relies on the understanding of the initial conditions e.g. the eccentricity of the nuclear overlap region. HADES [1] provides a large acceptance combined with a high mass-resolution and therefore allows to study di-electron and hadron production in heavy-ion collisions with unprecedented precision. In this contribution, the capability of HADES to study flow harmonics by utilizing multi-particle azimuthal correlation techniques is discussed. Due to the high statistics of seven billion Au+Au collisions at 1.23 AGeV collected in 2012, a systematic study of higher-order flow harmonics, the differentiation between collective and non-flow effects, and as well the multi-differential (pt, rapidity, centrality) analysis is possible.
In this thesis, the flow coefficients vn of the orders n = 1 − 6 are studied for protons and light nuclei in Au+Au collisions at Ebeam = 1.23 AGeV, equivalent to a center-of-mass energy in the nucleon-nucleon system of √sNN = 2.4 GeV. The detailed multi-differential measurement is performed with the HADES experiment at SIS18/GSI. HADES, with its large acceptance, covering almost full azimuth angle, combined with its high mass-resolution and good particle-identification capability, is well equipped to study the azimuthal flow pattern not only for protons, deuterons, and tritons but also for charged pions, kaons, the φ-mesons, electrons/positrons, as well as light nuclei like helions and alphas. The high statistics of more than seven billion Au-Au collisions recorded in April/May 2012 with HADES enables for the first time the measurement of higher order flow coefficients up to the 6th harmonic. Since the Fourier coefficient of 7th and 8th order are beyond the statistical significance only an upper bound is given. The Au+Au collision system is the largest reaction system with the highest particle multiplicities, which was measured so far with HADES. A dedicated correction method for the flow measurement had to be developed to cope with the reconstruction in-efficiencies due to occupancies of the detector system. The systematical bias of the flow measurement is studied and several sources of uncertainties identified, which mainly arise from the quality selection criteria applied to the analyzed tracks, the correction procedure for reconstruction inefficiencies, the procedures for particle identification (PID) and the effects of an azimuthally non-uniform detector acceptance. The systematic point-to-point uncertainties are determined separately for each particle type (proton, deuteron and triton), the order of the flow harmonics vn, and the centrality class. Further, the validity of the results is inspected in the range of their evaluated systematic uncertainties with several consistency checks. In order to enable meaningful comparisons between experimental observations and predictions of theoretical models, the classification of events should be well defined and in sufficiently narrow intervals of impact parameter. Part of this work included the implementation of the procedure to determine the centrality and orientation of the reaction.
In the conclusion the experimental results are discussed, including various scaling properties of the flow harmonics. It is found that the ratio v4/v2 for protons and light nuclei (deuterons and tritons) at midrapidity for all centrality classes approaches values close to 0.5 at high transverse momenta, which was suggested to be indicative for an ideal hydrodynamic behaviour. A remarkable scaling is observed in the pt dependence of v2 (v4) at mid-rapidity of the three hydrogen isotopes, when dividing by their nuclear mass number A (A^2) and pt by A. This is consistent with naive expectations from nucleon coalescence, butraises the question whether this mass ordering can also be explained by a hydrodynamical-inspired approach, like the blast-wave model. The relation of v2 and v4 to the shape of the initial eccentricity of the collision system is studied. It is found that v2 is independent of centrality for all three particle species after dividing it by the averaged second order participant eccentricity v2/⟨ε2⟩. A similar scaling is shown for v4 after division by ⟨ε2⟩^2.
This thesis contains three theoretical works about certain aspects of the interplay of electronic correlations and topology in the Hubbard model.
In the first part of this thesis, the applicability of elementary band representations (EBRs) to diagnose interacting topological phases, that are protected by spatial symmetries and time-reversal-symmetry, in terms of their single-particle Matsubara Green’s functions is investigated. EBRs for the Matsubara Green’s function in the zero-temperature limit can be defined via the topological Hamiltonian. It is found that the Green’s function EBR classification can only change by (i) a gap closing in the spectral function at zero frequency, (ii) the Green’s function becoming singular i.e. having a zero eigenvalue at zero frequency or (iii) the Green’s function breaking a protecting symmetry. As an example, the use of the EBRs for Matsubara Green’s functions is demonstrated on the Su-Schriefer-Heeger model with exact diagonalization.
In the second part the Two-Particle Self-Consistent approach (TPSC) is extended to include spin-orbit coupling (SOC). Time-reversal symmetry, that is preserved in the presence of SOC, is used to derive new TPSC self-consistency equations including SOC. SOC breaks spin rotation symmetry which leads to a coupling of spin and charge channel. The local and constant TPSC vertex then consists of three spin vertices and one charge vertex. As a test case to study the interplay of Hubbard interaction and SOC, the Kane-Mele-Hubbard model is studied. The antiferromagnetic spin fluctuations are the leading instability which confirms that the Kane-Mele-Hubbard model is an XY antiferromagnet at zero temperature. Mixed spin-charge fluctuations are found to be small. Moreover, it is found that the transversal spin vertices are more strongly renormalized than the longitudinal spin vertex, SOC leads to a decrease of antiferromagnetic spin fluctuations and the self-energy shows dispersion and sharp features in momentum space close to the phase transition.
In the third part TPSC with SOC is used to calculate the spin Hall conductivity in the Kane-Mele-Hubbard model at finite temperature. The spin Hall conductivity is calculated once using the conductivity bubble and once including vertex corrections. Vertex corrections for the spin Hall conductivity within TPSC corresponds to the analogues of the Maki-Thompson contributions which physically correspond to the excitation and reabsorption of a spin, a charge or a mixed spin-charge excitation by an electron. At all temperatures, the vertex corrections show a large contribution in the vicinity of the phase transition to the XY antiferromagnet where antiferromagnetic spin fluctuations are large. It is found that vertex corrections are crucial to recover the quantized value of −2e^2/h in the zero-temperature limit. Further, at non-zero temperature, increasing the Hubbard interaction leads to a decrease of the spin Hall conductivity. The results indicate that scattering of electrons off antiferromagnetic spin fluctuations renormalize the band gap. Decreasing the gap can be interpreted as an effective increase of temperature leading to a decrease of the spin Hall conductivity.
Many Polyakov loop models can be written in a dual formulation which is free of sign problem even when a non-vanishing baryon chemical potential is introduced in the action. Here, results of numerical simulations of a dual representation of one such effective Polyakov loop model at finite baryon density are presented. We compute various local observables such as energy density, baryon density, quark condensate and describe in details the phase diagram of the model. The regions of the first order phase transition and the crossover, as well as the line of the second order phase transition, are established. We also compute several correlation functions of the Polyakov loops.