Institutes
Refine
Year of publication
Document Type
- Article (401)
- Doctoral Thesis (222)
- Preprint (43)
- Conference Proceeding (17)
- Master's Thesis (12)
- Bachelor Thesis (9)
- Contribution to a Periodical (9)
- Book (5)
- Habilitation (3)
- Other (1)
Has Fulltext
- yes (724)
Is part of the Bibliography
- no (724) (remove)
Keywords
- Relativistic heavy-ion collisions (6)
- ALICE (5)
- Black holes (5)
- Heavy-ion collisions (5)
- LHC (5)
- Toroidales Magnetfeld (5)
- Equation of state (4)
- Fluctuations (4)
- HADES (4)
- QCD (4)
Institute
In this thesis, we use lattice QCD to study a part of the QCD phase diagram, specifically the QCD phase transition at mu=0, where the QCD matter changes from hadron gas to quark-gluon plasma (QGP) with increasing temperature.
This phase transition takes place as a crossover, but when theoretically changing the masses of the quarks, the order of the phase transition changes as well.
We focus on the region of heavy quark masses with Nf=2 flavours, where we investigate the critical quark mass at the second order phase transition in the form of a Z2 point between the first-order and the crossover region.
The first-order region is positioned at infinitely heavy quarks. As the quark masses decrease, the associated Z3 centre symmetry breaks explicitly, causing the first-order phase transition to weaken until it turns into the Z2 point and finally into a crossover.
We study this Z2 point using simulations at Nf=2 and lattices of the sizes Nt = {6, 8, 10, 12}, partially building on previous work, in which the simulations for Nt = {6, 8, 10} were started.
The simulations for Nt=12 are not finished yet though, but we were able to draw some preliminary conclusions. These simulations are run on GPUs and CPUs, using the codes Cl2QCD and open-QCD-FASTSUM, respectively. Afterwards, the data goes through a first analysis step in the form of the Python program PLASMA, preparing it for the two techniques we use to analyse the nature of the phase transition.
As a first, reliable analysis method, we perform a finite size scaling analysis of the data to find the location of the Z2 point. Since we are using lattice QCD, performing a continuum extrapolation is necessary to reach the continuum result.
In regard to this, the finite size scaling analysis method is hampered by the excessive amount of simulated data that is needed regarding statistics and the total number of simulations, which is why this thesis is only an intermediate step towards the continuum limit.
This also leads to the second analysis technique we explore in this thesis.
We start to design a Landau theory which describes the phase boundary for heavy masses at Nf=2 based on the simulated data.
We develop a Landau functional for every Nt we have simulation data for.
Albeit the results are not at the same precision as the ones from the finite size scaling analysis, we are able to reproduce the position of the Z2 point for every Nt.
Even though we are not able to take a continuum extrapolation right now, after more development takes place in future works, this approach might, in the long run, lead to a continuum result that won't need as many simulations as the finite size scaling analysis.
Im Rahmen dieser Arbeit wurde ein verbessertes Buncher-System für Hochfrequenzbeschleuniger mit niedrigem und mittlerem Ionenstrom entwickelt. Die entwickelte Methodik hat ermöglicht, ein effektives, vereinfachtes Buncher-System für die Injektion in HF-Beschleuniger wie RFQs, Zyklotrons, DTLs usw. zu entwerfen, welches kleine Ausgangsemittanzen und beträchtliche Strahltransmissionen erzielt. Um einen mono-energetischen und kontinuierlichen Strahl aus einer Ionenquelle für den Einschuss in eine Hochfrequenz-Beschleunigerstruktur anzupassen, wird eine Energiemodulation benötigt, die im weiteren Verlauf (Driftstrecke) zur Längsfokussierung des Strahls führt. Durch eine Sägezahnwellenform wird die ideale Energiemodulation aufgrund der linearen Abhängigkeit zwischen der Energie der Teilchen und ihren relativen Phasen erreicht. Dies ist jedoch technologisch nicht möglich, da Teilchenbeschleuniger Spannungsniveaus im Bereich kV bis 100 kV benötigen. Dagegen ist für eine solche Zielsetzung eine räumliche Trennung der sinusförmigen Anregung mit der Grundfrequenz und höheren Harmonischen möglich.
Daher wurde in dieser Arbeit ein verbesserter harmonischer Buncher, der sogenannte „Double Drift Harmonic Buncher - DDHB“ entwickelt, welcher zahlreiche Vorteile hat. Eine geringe longitudinale Emittanz sowie finanzielle Aspekte sprechen für diesen Lösungsansatz. Die Hauptelemente eines DDHB Systems sind zwei Kavitäten, die durch eine Driftlänge L1 getrennt sind, wobei der erste Resonator mit der Grundfrequenz bei -90° synchroner Phase und angelegter Spannung V1 und der zweite Resonator bei der zweiten harmonischen Frequenz mit +90 synchroner Phase und angelegter Spannung V2 betrieben werden. Schließlich ist eine zweite Drift L2 am Ende des Arrays für eine longitudinale Strahlfokussierung am Hauptbeschleunigereingang erforderlich. Somit erfüllt ein solcher Aufbau das angestrebte Ziel einer hohen Einfangseffizienz und einer kleinen longitudinalen Emittanz durch Anpassen der vier Designparameter V1, L1, V2 und L2.
Das Verständnis der Fokussierung, ausgehend von einem Gleichstromstrahl, einschließlich der Raumladungskräfte, ist einer der wesentlichen Bestandteile der Strahlphysik. Viele kommerzielle Codes bieten Simulationsmöglichkeiten in diesem Anwendungsbereich. Ihre Ansätze bleiben jedoch dem Anwender meist verborgen, oder es fehlen wichtige Details zur genauen Abbildung des vorliegenden Konzepts. Daher bestand eine Hauptaufgabe dieser Arbeit darin, einen speziellen Multi-Particle-Tracking-Beam-Dynamics-Code (BCDC) zu entwickeln, bei dem der Raumladungseffekt während des Bunch-Vorgangs, ausgehend von einem DC-Strahl berechnet wird. Der BCDC - Code enthält elementare Routinen wie Drift und Beschleunigungsspalt oder magnetische Linse für die transversale Strahlfokussierung und Raumladungsberechnungen unter Berücksichtigung der Auswirkungen der nächsten Nachbar-Bunche (NNB). Der Raumladungsalgorithmus in BCDC basiert auf einer direkten Coulomb- Gitter-Gitter-Wechselwirkung und Berechnungen des elektrischen Feldes durch Lokalisierung der Ladungsdichte auf einem kartesischen Gitter. Um Genauigkeit zu erreichen, werden die Feldberechnungen in Längsrichtung symmetrisch um das zentrale Bucket (βλ-Größe) erweitert, so dass das Simulationsfeld dreimal so groß ist. Die zentrale Teilchenverteilung wird dann nach jedem Schritt in die benachbarten Buckets kopiert. Anschließend werden die resultierenden Felder im Hauptgitterfeld neu berechnet, indem die elektrischen Felder im Hauptgitterfeld mit denen aus den benachbarten Regionen überlagert werden. Ohne diese Methode würde z. B. ein kontinuierlicher Strahl, welcher jedoch in der Simulation nur innerhalb einer Zelle der Länge βλ definiert ist, zu einer resultierenden Raumladungsfeldkomponente Ez an beiden Rändern der Zelle führen. Ein solches unphysikalisches Ergebnis konnte durch die Anwendung der NNB-Technik bereits weitgehend eliminiert werden. Zusätzlich zum NNB-Feature verfügt das BCDC über eine weitere Besonderheit nämlich die sogenannte Raumladungskompensation (SCC). Aufgrund der Ionisierung des Restgases kommt es entlang des Niederenergiestrahltransports zu einer teilweisen Raumladungskompensation, und zwar am und hinter dem Bunchersystem mit unterschiedlichen Prozentsätzen. Eines der Hauptziele des DDHB-Konzepts besteht darin, es für Hochstromstrahlanwendungen zu entwickeln. Dabei ermöglicht die teilweise Raumladungskompensation, dass das Design in der Praxis höhere Stromniveaus erreicht. Dadurch ist das BCDC-Programm ein leistungsstarkes Werkzeug für Simulationen in künftigen, stromstarken Projekten. Proof-of-Principle-Designs wurden in dieser Arbeit entwickelt.
Supermassive black hole binaries (SMBHBs) are among the most powerful known sources of gravitational waves (GWs). Accordingly, these systems could dominate the stochastic gravitational wave background (GWB) in the micro- and millihertz frequency range. The time until the merger of two SMBHs in the nucleus of a galaxy can be shortened through dynamical friction due to the presence of dark matter (DM) spikes around the SMBHs. To calculate the orbital evolution of individual SMBHBs within the Newtonian approximation, the SMBHBpy code is developed. This work confirms that the GW signals from SMBHBs with DM spikes can be clearly distinguished from those without any matter. Making use of the upper limit on the characteristic strain of the GWB derived from the data of the Cassini spacecraft mission in 2001/2002, a lower limit on the matter density around SMBHBs is derived in this study. The result is subsequently compared with the theoretical density profiles for cold dark matter and self-interacting dark matter spikes.
Binary neutron star mergers represent unique observational phenomena because all four fundamental interactions play an important role at various stages of their evolution by leaving imprints in astronomical observables. This makes their accurate numerical modeling a challenging multiphysics problem that promises to increase our understanding of the high-energy astrophysics at play, thereby providing constraints for the underlying fundamental theories such as the gravitational interaction or the strong interaction of dense matter. For example, the first and so far only multi-messenger observation of the binary neutron star merger GW170817 resulted in numerous bounds on the parameters of isolated non-rotating neutron stars, e.g., their maximum mass or their distribution in radii, which can be directly used to constrain the equation of state of cold nuclear matter. While many of these results stem from the observation of the inspiral gravitational-wave signal, the postmerger phase of binary neutron star mergers encodes even more details about the extreme physics of hot and dense neutron star matter. In this Thesis we focus on the exploration of dissipative and shearing effects in binary neutron star mergers in order to identify novel approaches to constrain hot and dense neutron star matter.
The first effect is the well-motivated dissipation of energy due to the bulk viscosity which arises from violations of weak chemical equilibrium. We start by exploring the impact of bulk viscosity on black-hole accretion. This simplified problem gives us the opportunity to develop a test case for future codes taking into account the effects of dissipation in a fully general-relativistic setup and build intuition in the physics of relativistic dissipation. Next, we move on to isolated neutron stars and binary neutron star mergers by developing a robust implementation of bulk-viscous dissipation for numerical relativity simulations. We test our implementation by calculating the damping of eigenmodes of isolated neutron stars and the violent migration scenario. Finally, we present the first results on the impact of bulk viscosity on binary neutron star mergers. We identify a number of ways how bulk viscosity impacts the postmerger phase, out of which the suppression of gravitational-wave emission and dynamical mass ejection are the most notable ones.
In the last part of this Thesis we investigate how the shearing dynamics at the beginning of the merger affects the amplification of different initial magnetic-field topologies. We explore the hypothesis that magnetic fields which are located only in a small region near the stellar surface prior to merger lead to a weaker magnetic-field amplification. We show first evidence which confirms this hypothesis and discuss possible implications for constraining the physics of superconduction in cold neutron stars.
We continue previous investigations of the (inhomogeneous) phase structure of the Gross-Neveu model in a noninteger number of spatial dimensions (1≤d<3) in the limit of an infinite number of fermion species (N→∞) at (non)zero chemical potential μ. In this work, we extend the analysis from zero to nonzero temperature T.
The phase diagram of the Gross-Neveu model in 1≤d<3 spatial dimensions is well known under the assumption of spatially homogeneous condensation with both a symmetry broken and a symmetric phase present for all spatial dimensions. In d=1 one additionally finds an inhomogeneous phase, where the order parameter, the condensate, is varying in space. Similarly, phases of spatially varying condensates are also found in the Gross-Neveu model in d=2 and d=3, as long as the theory is not fully renormalized, i.e., in the presence of a regulator. For d=2, one observes that the inhomogeneous phase vanishes, when the regulator is properly removed (which is not possible for d=3 without introducing additional parameters).
In the present work, we use the stability analysis of the symmetric phase to study the presence (for 1≤d<2) and absence (for 2≤d<3) of these inhomogeneous phases and the related moat regimes in the fully renormalized Gross-Neveu model in the μ,T-plane. We also discuss the relation between "the number of spatial dimensions" and "studying the model with a finite regulator" as well as the possible consequences for the limit d→3.
Inhomogeneous condensation in the Gross-Neveu model in noninteger spatial dimensions 1 ≤ d < 3
(2023)
The Gross-Neveu model in the N→∞ limit in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to noninteger spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to noninteger spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability toward an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various (2+1)-dimensional four-fermion and Yukawa models. All models are studied at nonzero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially nonuniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave-function renormalization is negative, in these models is ruled out.
A synchrotron is a particular type of cyclic particle accelerator and the first accelerator concept to enable the construction of large-scale facilities [10], such as the largest particle accelerator in the world, the 27-kilometre-circumference Large Hadron Collider (LHC) by CERN near Geneva, Switzerland, the European Synchrotron Radiation Facility (ESRF) in Grenoble, France for the synchrotron radiation, the superconducting, heavy ion synchrotron SIS100 under construction for the FAIR facility at GSI, Darmstadt, Germany and so on. Unlike a cyclotron, which can accelerate particles starting at low kinetic energy, a synchrotron needs a pre-acceleration facility to accelerate particles to an appropriate initial value before synchrotron injection. A pre-acceleration can be realized by a chain of other accelerator structures like a linac, a microtron in case of electrons, for example, Proton and ion injectors Linac 4 and Linac 3 for the LHC, UNLAC as the injector for the SIS18 in GSI and in future the SIS18 as injector for the SIS100. The linac is a commonly used injector for the ion synchrotron and consists of some key components. The three main parts of a linac are: An ion source creating the particles, a buncher system or an RFQ followed by the main drift tube accelerator DTL. In order to meet the energy and the beam current requirement of a synchrotron injector linac, its cost is a remarkable percentage of the total facility costs.
However, the normal conducting linac operation at cryogenic temperatures can be a promising solution in improving the efficiency and reducing the costs of a linac. Synchrotron injectors operate at very low duty factor with beam pulse lengths in 1 micros to 100 micros range, as most of the time is needed to perform the synchrotron cycle. Superconducting linacs are not convenient, as they cannot efficiently operate at low duty factor and high beam currents.
The cryogenic operation of ion linacs is discussed and investigated at IAP in Frankfurt since around 2012 [1, 37]. The motivation was to develop very compact synchrotron injectors at reduced overall linac costs per MV of acceleration voltage. As the needed beam currents for new facilities are increasing as well, the new technology will also allow an efficient realization of higher injector linac energies, which is needed in that case. Operating normal conducting structures at cryogenic temperature exploits the significantly higher conductivity of copper at temperatures of liquid nitrogen and below. On the other hand, the anomalous skin effect reduces the gain in shunt impedance quite a bit[25, 31, 9]. Some intense studies and experiments were performed recently, which are encouraging with respect to increased field levels at linac operation temperatures between 30 K and 70 K [17, 24, 4, 23, 5, 8]. While these studies are motivated by applications in electron acceleration at GHz-frequencies, the aim of this paper is to find applications in the 100 to 700 MHz range, typical for proton and ion acceleration. At these frequencies, a higher impact in saving RF power is expected due to the larger skin depth, which is proportional to the frequency to the power of negative half with respect to the normal skin effect. On the other hand, it is assumed that the improvement in maximum surface field levels will be similar to what was demonstrated already for electron accelerator cavities. This should allow to find a good compromise between reduced RF power needs for achieving a given accelerator voltage and a reduced total linac length to save building costs.
A very important point is the temperature stability of the cavity surface during the RF pulse. This is of increasing importance the lower the operating temperature is chosen: the temperature dependence of the electric conductivity in copper gets rather strong below 80 K, as long as the RRR - value of the copper is adequate. It is very clear, that this technology is suited for low duty cycle operated cavities only - with RF pulse lengths below one millisecond. At longer pulses the cavity surface will be heated within the pulse to temperatures, where the conductivity advantage is reduced substantially. These conditions fit very well to synchrotron injectors or to pulsed beam power applications.
H – Mode structures of the IH – and of the CH – type are well-known to have rather small cavity diameters at a given operating frequency. Moreover, they can achieve effective acceleration voltage gains above 10 MV/m even at low beam energies, and already at room temperature operation[29]. With the new techniques of 3d – printing of stainless steel and copper components one can reduce cavity sizes even further – making the realization of complex cooling channels much easier.
Another topic are copper components in superconducting cavities – like power couplers. It is of great importance to know exactly the thermal losses at these surfaces, which can’t be cooled efficiently in an easy way.
In this work we study the 3+1-dimensional Nambu-Jona-Lasinio (NJL) model in the mean field-approximation. We carry out calculations using five different regularization schemes (two continuum and three lattice regularization schemes) with particular focus on inhomogeneous phases and condensates. The regularization schemes lead to drastically different inhomogeneous regions. We provide evidence that inhomogeneous condensates appear for all regularization schemes almost exclusively at values of the chemical potential and with wave numbers, which are of the order of or even larger than the corresponding regulators. This can be interpreted as indication that inhomogeneous phases in the 3+1-dimensional NJL model are rather artifacts of the regularization and not a consequence of the NJL Lagrangian and its symmetries.
In the last twenty years, a variety of unexpected resonances had been observed within the charmonium mass region. Although the existence of unconventional states has been predicted by the quantum chromodynamics (QCD), a quantum field theory describing the strong force, a clear evidence was missing. The Y(4260) is such an unexpected and supernummerary state, first observed at BaBar in 2005, and aroused great interest, because it couples much stronger to hidden charm decays (charm-anticharm states like J/Psi or h_c) instead of open charm decays (D meson pairs). This is unusual for states with masses above the D anti-D threshold. Furthermore, it decays into a charged exotic state Y(4260)->Z_c(3900)^+- pi^-+. The charge of the Z_c(3900)^+- is an indication that it comprises of two more quarks than the charm-anticharm pair, and could therefore be assumed to be a four-quark state. Due to these still not understood properties of these QCD-allowed states, they are referred to as exotic XYZ states to emphasize their particularity.
In 2017, the collaboration of the Beijing Spectrometer III (BESIII) investigated the production reaction of the Y(4260) resonance based on a high-luminosity data set. This significantly improved precision of the measurement of the cross-section sigma(e+e- -> J/Psi pi^+ pi^-) permitted a resolution into two resonances, the Y(4230) and the Y(4360). The Z_c(3900)^+- had been discovered by the BESIII collaboration in 2013, thus this experiment at the Beijing Electron-Positron Collider II (BEPCII) is a top-performing facililty to study exotic charmonium-like states.
In this work, an inclusive reconstruction of the strange hyperon Lambda in the charmonium mass region is performed to study possible decays of Y states in order to provide further insight into their nature. Finding more states or new decay channels may provide crucial hints to understand the strong interaction beyond nonperturbative approaches.
Three resonances are observed in the energy dependent cross-section: the first with a mass of (4222.01 +- 5.68) MeV and a width of (154.26 +- 28.16) MeV, the second with a mass of (4358.88 +- 4.97) MeV and a width of (49.58 +- 13.54) MeV and the third with a mass of (4416.41 +- 2.37) MeV and a width of (23.88 +- 7.18) MeV. These resonances, with a statistical significance Z > 5sigma, can be interpreted as the states Y(4230), Y(4360) and psi(4415).
Additionally, a proton momentum-dependent analysis strategy has been used in terms of the inclusiveness of the reconstruction and to address the momentum discrepancies between generic MC and measured data.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various 2+1-dimensional four-fermion and Yukawa models. All models are studied at non-zero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially non-uniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave function renormalization is negative, in these models is ruled out.
Inhomogeneous condensation in the Gross-Neveu model in noninteger spatial dimensions 1 ≤ d < 3
(2023)
The Gross-Neveu model in the N→∞ approximation in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to non-integer spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to non-integer spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability towards an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
Artificial intelligence in heavy-ion collisions : bridging the gap between theory and experiments
(2023)
Artificial Intelligence (AI) methods are employed to study heavy-ion collisions at intermediate collision energies, where high baryon density and moderate temperature QCD matter is produced. The experimental measurements of various conventional observables such as collective flow, particle number fluctuations, etc. are usually compared with expensive model calculations to infer the physics governing the evolution of the matter produced in the collisions. Various experimental effects and processing algorithms can greatly affect the sensitivity of these observables. AI methods are used to bridge this gap between theory and experiments of heavy-ion collisions. The problems with conventional methods of analyzing experimental data are illustrated in a comparative study of the Glauber MC model and the UrQMD transport model. It is found that the centrality determination and the estimated fluctuations of the number of participant nucleons suffer from strong model dependencies for Au-Au collisions at 1.23 AGeV. This can bias the results of the experimental analysis if the number of participant nucleons used is not consistent throughout the analysis and in the final model-to-data comparison. The measurable consequences of this model dependence of the number of participant nucleons are also discussed. In this context, PointNet-based AI models are developed to accurately reconstruct the impact parameter or the number of participant nucleons in a collision event from the hits and/or reconstructed track of particles in 10 AGeV Au-Au collisions at the CBM experiment. In the last part of the thesis, different AI methods to study the equation of state (EoS) at high baryon densities are discussed. First, a Bayesian inference is performed to constrain the density dependence of the EoS from the available experimental measurements of elliptical flow and mean transverse kinetic energy of mid rapidity protons in intermediate energy collisions. The UrQMD model was augmented to include arbitrary potentials (or equivalently the EoSs) in the QMD part to provide a consistent treatment of the EoS throughout the evolution of the system. The experimental data constrain the posterior constructed for the EoS for densities up to four times saturation density. However, beyond three times saturation density, the shape of the posterior depends on the choice of observables used. There is a tension in the measurements at a collision energy of about 4 GeV. This could indicate large uncertainties in the measurements, or alternatively the inability of the underlying model to describe the observables with a given input EoS. Tighter constraints and fully conclusive statements on the EoS require accurate, high statistics data in the whole beam energy range of 2-10 GeV, which will hopefully be provided by the beam energy scan programme of STAR-FXT at RHIC, the upcoming CBM experiment at FAIR, and future experiments at HIAF and NICA. Finally, it is shown that the PointNet-based models can also be used to identify the equation of state in the CBM experiment. Despite the uncertainties due to limited detector acceptance and biases in the reconstruction algorithms, the PointNet-based models are able to learn the features that can accurately identify the underlying physics of the collision. The PointNet-based models are an ideal AI tool to study heavy-ion collisions, not only to identify the geometric event features, such as the impact parameter or the number of participant nucleons, but also to extract abstract physical features, such as the EoS, directly from the detector outputs.
Human feline leukaemia virus subgroup C receptor-related proteins 1 and 2 (FLVCR1 and 2) are major facilitator superfamily transporters from the solute carrier family 49. Dysregulation of these ubiquitous transporters has been linked to various haematological and neurological disorders. While both FLVCRs were initially proposed to hold a physiological function in heme transport, subsequent studies questioned this notion. Here, we used structural, computational and biochemical methods and conclude that these two FLVCRs function as human choline transporters. We present cryo-electron microscopy structures of FLVCRs in different inward- and outward-facing conformations, captured in the apo state or in complex with choline in their translocation pathways. Our findings provide insights into the molecular framework of choline coordination and transport, largely mediated by conserved cation-π interactions, and further illuminate the conformational dynamics of the transport cycle. Moreover, we identified a heme binding site on the protein surface of the FLVCR2 N-domain, and observed that heme actively drives the conformational transitions of the protein. This auxiliary binding site might indicate a potential regulatory role of heme in the FLVCR2 transport mechanisms. Our work resolves the contested substrate specificity of the FLVCRs, and sheds light on the process of maintaining cellular choline homeostasis at the molecular level.
Vertebrate life depends on renal function to filter excess fluid and remove low-molecular-weight waste products. An essential component of the kidney filtration barrier is the slit diaphragm (SD), a specialized cell-cell junction between podocytes. Although the constituents of the SD are largely known, its molecular organization remains elusive. Here, we use super-resolution correlative light and electron microscopy to quantify a linear rate of reduction in albumin concentration across the filtration barrier under no-flow conditions. Next, we use cryo-electron tomography of vitreous lamellae from high-pressure frozen native glomeruli to analyze the molecular architecture of the SD. The resulting densities resemble a fishnet pattern. Fitting of Nephrin and Neph1, the main constituents of the SD, results in a complex interaction pattern with multiple contact sites between the molecules. Using molecular dynamics simulations, we construct a blueprint of the SD that explains its molecular architecture. Our architectural understanding of the SD reconciles previous findings and provides a mechanistic framework for the development of novel therapies to treat kidney dysfunction.
Vertebrate life depends on renal function to filter excess fluid and remove low-molecular-weight waste products. An essential component of the kidney filtration barrier is the slit diaphragm (SD), a specialized cell-cell junction between podocytes. Although the constituents of the SD are largely known, its molecular organization remains elusive. Here, we use super-resolution correlative light and electron microscopy to quantify a linear rate of reduction in albumin concentration across the filtration barrier. Next, we use cryo-electron tomography of vitreous lamellae from high-pressure frozen native glomeruli to analyze the molecular architecture of the SD. The resulting densities resemble a fishnet pattern. Fitting of Nephrin and Neph1, the main constituents of the SD, results in a complex interaction pattern with multiple contact sites between the molecules. Using molecular dynamics flexible fitting, we construct a blueprint of the SD, where we describe all interactions. Our architectural understanding of the SD reconciles previous findings and provides a mechanistic framework for the development of novel therapies to treat kidney dysfunction.
In this work, inhomogeneous chiral phases are studied in a variety of Four-Fermion and Yukawa models in 2+1 dimensions at zero and non-zero temperature and chemical potentials. Employing the mean-field approximation, we do not find indications for an inhomogeneous phase in any of the studied models. We show that the homogeneous phases are stable against inhomogeneous perturbations. At zero temperature, full analytic results are presented.
Inhomogeneous condensation in the Gross-Neveu model in non-integer spatial dimensions 1 ≤ d < 3
(2023)
he Gross-Neveu model in the N→∞ approximation in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to non-integer spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to non-integer spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability towards an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
We investigate the space-time dependence of electromagnetic fields produced by charged participants in an expanding fluid. To address this problem, we need to solve the Maxwell's equations coupled to the hydrodynamics conservation equation, specifically the relativistic magnetohydrodynamics (RMHD) equations, since the charged participants move with the flow. To gain analytical insight, we approximate the problem by solving the equations in a fixed background Bjorken flow, onto which we solve Maxwell's equations. The dynamical electromagnetic fields interact with the fluid's kinematic quantities such as the shear tensor and the expansion scalar, leading to additional non-trivial coupling. We use mode decomposition of Green's function to solve the resulting non-linear coupled wave equations. We then use this function to calculate the electromagnetic field for two test cases: a point source and a transverse charge distribution. The results show that the resulting magnetic field vanishes at very early times, grows, and eventually falls at later times.
Baryonic models of ultra-low-mass compact stars for the central compact object in HESS J1731-347
(2023)
The recent attempt on mass and radius inference of the central compact object within the supernova remnant HESS J1731-347 suggests for this object an unusually low mass of M=0.77−0.17+0.20M⊙ and a small radius of R=10.4−0.78+0.86km. We explore the ways such a result can be accommodated within models of dense matter with heavy baryonic degrees of freedom which are constrained by the multi-messenger observations. We find that to do so using only purely nucleonic models, one needs to assume a rather small value of the slope of symmetry energy Lsym. Once heavy baryons are included higher values of the slope Lsym become acceptable at the cost of a slightly reduced maximum mass of static configuration. These two scenarios are distinguished by the particle composition and will undergo different cooling scenarios. In addition, we show that the universalities of the I-Love-Q relations for static configurations can be extended to very low masses without loss in their accuracy.
We study the production of the light JPC=1−+ hybrid isoscalars η1′(1855) and the yet undiscovered η1(1660) as intermediate states in the radiative decays of the charmonium (J/ψ) to two conventional mesons using a flavor symmetric Lagrangian. For this purpose, we use the J/ψ→γη1′(1855)→γηη′ process as the reference. We find that some of the decay channels have branching ratios similar to or larger than that of the γηη′ channel and are sensitive to the mixing between the hybrid isoscalars. We propose that relatively stable γηf1(1285) channel be explored for the presence of the light hybrid isoscalar η1hyb(1660). We also exploit the strong decay channels containing at least one vector meson to study the radiative decays of the whole hybrid nonet {π1(1600),K1hyb(1750),η1hyb(1660),η1′(1855)}. We find that the hybrids cannot radiatively decay into the I=0 pseudoscalars. Furthermore, the vector decay channels ((ρ/ω/ϕ)γ) of the hybrid isoscalars are sensitive to the strangeness content of the hybrids. We also provide estimates for the branching fractions for the radiative production and partial widths for the radiative decays of the hybrids.
We investigate the space-time dependence of electromagnetic fields produced by charged participants in an expanding fluid. To address this problem, we need to solve the Maxwell's equations coupled to the hydrodynamics conservation equation, specifically the relativistic magnetohydrodynamics (RMHD) equations, since the charged participants move with the flow. To gain analytical insight, we approximate the problem by solving the equations in a fixed background Bjorken flow, onto which we solve Maxwell's equations. The dynamical electromagnetic fields interact with the fluid's kinematic quantities such as the shear tensor and the expansion scalar, leading to additional non-trivial coupling. We use mode decomposition of Green's function to solve the resulting non-linear coupled wave equations. We then use this function to calculate the electromagnetic field for two test cases: a point source and a transverse charge distribution. The results show that the resulting magnetic field vanishes at very early times, grows, and eventually falls at later times.
This work focuses on the investigation of K+, K- and ϕ-meson production in Ag(1.58 A GeV)+Ag collisions. The energetically cheapest channel for direct K+ production in binary NN-collisions NN→NΛK+ lies at exactly this energy. For the remaining K- and ϕ-mesons, an excess energy of 0.31 GeV and 0.34 GeV in the centre of mass system has to be provided by the system. This makes these particles an excellent probe for effects inside the medium.
K+ and K- mesons can be reconstructed directly as they possess a cτ of approximately 3.7 m. Using the approximately 3 billion recorded Ag(1.58 A GeV)+Ag 0-30% most central collision events, all reconstructed K+ and K- within the detector acceptance are investigated for their kinematic properties and their particle production rates compared to a selection of existing models.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various 2+1-dimensional four-fermion and Yukawa models. All models are studied at non-zero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially non-uniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave function renormalization is negative, in these models is ruled out.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √sNN = 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconstructed via its charged two-body mesonic decay channel (3ΛH → 3He + π− and the charge-conjugate process). The measured values τ = [253 ± 11 (stat) ± 6 (syst)] ps and BΛ = [102 ± 63 (stat) ± 67 (syst)] keV are compatible with predictions from effective field theories and confirm that the 3ΛH structure is consistent with a weakly bound system.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √= 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconsNN structed via its charged two-body mesonic decay channel (3ΛH→ 3He + π− and the charge-conjugate process). The measured values τ=[253±11 (stat.)±6 (syst.)] ps and BΛ=[102±63 (stat.)±67 (syst.)] keV are compatible with predictions from effective field theories and confirm that the 3ΛH structure is consistent with a weakly-bound system.
The Heidelberg Ion-Beam Therapy Centre (HIT) provides proton, helium, and carbon-ion beams with different energies and intensities for cancer treatment and oxygen-ion beams for experiments. For several experiments and possible future applications, such as helium ion beam radiography, a low-intensity ion beam monitor integrated into the dose delivery feedback system for the accelerator control is a necessary pre-requisite. The updated 2D prototype for this purpose consists of scintillating fibres with enhanced radiation hardness, silicon photomultipliers (SiPMs) to amplify the emitted light, and a dedicated front-end readout system (FERS) to process and record the generated signals. This setup was tested successfully on monitoring ion-beam position and profile horizontally and vertically, as well as the beam intensity, for all four ion types with energies from 50 to 430 MeV/u and intensities from 1E2 to 1E7 ions/s. Additionally, time-of-arrival (ToA) measurements on single ions have been successfully performed for a limited intensity range, allowing for ion tracking in a further update. This will reduce noise, and will also improve the accuracy and usability of ion radiography.
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.
The equation of state (EoS) of matter at extremely high temperatures and densities is currently not fully understood, and remains a major challenge in the field of nuclear physics. Neutron stars harbor such extreme conditions and therefore serve as celestial laboratories for constraining the dense matter EoS. In this thesis, we present a novel algorithm that utilizes the idea of Bayesian analysis and the computational efficiency of neural networks to reconstruct the dense matter equation of state from mass-radius observations of neutron stars. We show that the results are compatible with those from earlier works based on conventional methods, and are in agreement with the limits on tidal deformabilities obtained from the gravitational wave event, GW170817. We also observe that the resulting squared speed of sound from the reconstructed EoS features a peak, indicating a likely convergence to the conformal limit at asymptotic densities, as expected from quantum chromodynamics. The novel algorithm can also be applied across various fields faced with computational challenges in solving inverse problems. We further examine the efficiency of deep learning methods for analyzing gravitational waves from compact binary coalescences in this thesis. In particular, we develop a deep learning classifier to segregate simulated gravitational wave data into three classes: signals from binary black hole mergers, signals from binary neutron star mergers, or white noise without any signals. A second deep learning algorithm allows for the regression of chirp mass and combined tidal deformability from simulated binary neutron star mergers. An accurate estimation of these parameters is crucial to constrain the underlying EoS. Lastly, we explore the effects of finite temperatures on the binary neutron star merger remnant from GW170817. Isentropic EoSs are used to infer the frequencies of the rigidly rotating remnant and are noted to be significantly lower compared to previous estimates from zero temperature EoSs. Overall, this thesis presents novel deep learning methods to constrain the neutron star EoS, which will prove useful in future, as more observational data is expected in the upcoming years.
Efficient modeling and mitigation of quadrupole errors in synchrotrons and their beam transfer lines
(2023)
This thesis investigates the problem of estimating quadrupole errors on synchrotrons as well as how to minimize the influence of quadrupole errors for beam transfer lines (beamlines). It emphasizes the importance to treat possible error sources in all parts of an accelerator in order to provide constantly high beam quality to the experimental stations. While the presented methods have been investigated by using the example of the SIS18 synchrotron and the HEST beamlines at GSI Helmholtz Centre for Heavy Ion Research, they are equally relevant for the future synchrotrons and beamlines of the Facility for Antiproton and Ion Research in Europe (FAIR).
Part 1 discusses the problem of estimating quadrupole errors via orbit response measurements at synchrotrons. An emphasis is put on investigating the influence of the availability of steerer magnets and beam position monitors (BPMs) on the solvability of the inverse problem as well as on the propagation of measurement uncertainty for the estimation of quadrupole errors. The problem is approached via analytical considerations as well as via dedicated simulation studies. By developing an analytical expression for the Jacobian matrix, the theoretical boundaries for the solvability of the inverse problem are derived. Moreover, it is shown that the analytical expressions for the Jacobian matrix can be used during the fitting procedure to achieve a significant improvement in the computational efficiency by a factor $N_{steerers} \times N_{quadrupoles}$, where $N$ denotes the number of lattice elements of the respective type. The presented results are tested via dedicated measurements at the SIS18 synchrotron.
Part 2 discusses – complementary to part 1 – the influence of quadrupole errors in beam transfer lines with respect to the beam quality requirements given by the experimental stations. A preventive approach is presented which allows to minimize the influence of possible quadrupole errors on the degradation of beam quality. By identifying and selecting robust quadrupole configurations, a stable operation of the beamline can be enabled and the time needed by operators to readjust the beamline parameters can be reduced. The concept of beamline robustness is developed and is studied with the help of dedicated simulations. The simulation results are used to identify certain properties that distinguish robust from nonrobust quadrupole configurations. Also, various methods for improving the computational process of identifying robust quadrupole configurations are presented. The methods and results are tested via dedicated measurements at two different beamlines at GSI Helmholtz Centre for Heavy Ion Research and at Forschungszentrum Jülich.
The theoretical and experimental investigation of exotic hadrons like tetraquarks is an important branch of modern elementary particle physics. In this thesis I investigate different four-quark systems using lattice QCD and search for evidence of stable tetraquark states or resonances.
Lattice QCD as a non-perturbative approach to QCD allows an accurate and reliable determination of the masses of strongly bound hadrons.
However, most tetraquarks appear as weakly bound states or resonances, which makes a theoretical investigation using lattice QCD difficult due to the finite spatial volume. A rigorous treatment of such systems is feasible using the so-called Lüscher method. This allows to calculate the scattering amplitude based on the finite-volume energy spectrum determined in a lattice QCD calculation. Similarly to the analysis of experimental data, this scattering amplitude can be used to determine the binding energies of bound states or the masses and decay widths of resonances in the infinite volume.
In my work I calculate the low-energy energy spectra of different four-quark systems and use - if necessary - the Lüscher method to determine the masses of potential tetraquark states.
I focus on systems consisting of two heavy antiquarks and two light quarks, where at least one of the heavy antiquarks is a bottom quark.
Even though such tetraquarks have not yet been experimentally detected, they are considered promising candidates for particles that are stable with respect to the strong interaction.
A decisive step for successfully calculating low-lying energy levels for such four-quark systems is a carefully chosen set of creation operators, which represent the physical states most accurately. In addition to operators that generate a local structure where all four quarks are located at the same space-time point, I also use so-called scattering operators that resemble two spatially separated mesons. These scattering operators turned out to be relevant for successfully determining the lowest energy levels and are therefore essential, especially if a Lüscher analysis is carried out.
In my work, I considered two different lattice setups to study the four-quark systems $\bar{b}\bar{b}ud$ with $I(J^P)=0(1^+) $, $\bar{b}\bar{b}us$ with $J^P=1^+ $ and $\bar{b}\bar{c}ud$ with $I(J^P)=0(0^+) $ and $I(J^P)=0(1^+) $ and to predict potential tetraquark states. In both setups, I considered scattering operators. While in the first setup I used them only as annihilation operators, in the second setup they were included both as creation and annihilation operators. Additionally, in the second lattice setup, I performed a simplified investigation of the $\bar{b}\bar{b}ud$ system with $I(J^P)=0(1^-) $, which is a potential candidate for a tetraquark resonance. The results of the investigation of the mentioned four-quark systems can be summarized as follows:
For the $ \bar{b}\bar{b}ud $ four-quark system with $ I(J^P)=0(1^+) $ I found a deeply bound ground state slightly more than $ 100\,\textrm{MeV} $ below the lowest meson-meson threshold. The existence of a corresponding $\bar{b}\bar{b}ud$ tetraquark in the infinite volume was confirmed using a Lüscher analysis and possible systematic errors due to the use of lattice QCD were taken into account.
Similar results were obtained for the $ \bar{b}\bar{b}us $ four-quark system with $ J^P=1^+ $. Again, I found a ground state well below the lowest meson-meson threshold, but slightly weaker bound than for the $ \bar{b}\bar{b}ud $ system. Effects due to the finite volume turned out to be negligible for this system, as already predicted for the $ \bar{b}\bar{b}ud $ system. \item For the $ \bar{b}\bar{c}ud $ four-quark systems with $ (J^P)=0(0^+) $ and $ (J^P)=0(1^+) $ I was able to rule out the existence of a deeply bound tetraquark states based on the energy spectrum in the finite volume. However, by means of a scattering analysis using the Lüscher method, I found evidence a broad resonance for both channels.
In the case of the $ \bar{b}\bar{b}ud $ four-quark system with $ I(J^P)=0(1^-) $, I could neither confirm the existence of a resonance, nor rule out its existence with certainty.
In particular, my investigations showed that the results of the two different lattice simulations are consistent. The theoretical prediction of the bound tetraquark states $\bar{b}\bar{b}ud$ and $\bar{b}\bar{b}us$ as well as the tetraquark resonances in the $\bar{b}\bar{c}ud$ system in this work represent an important contribution to the future experimental search for exotic hadrons and can support the discovery of previously unobserved particles.
DNA binding redistributes activation domain ensemble and accessibility in pioneer factor Sox2
(2023)
More than 1600 human transcription factors orchestrate the transcriptional machinery to control gene expression and cell fate. Their function is conveyed through intrinsically disordered regions (IDRs) containing activation or repression domains but lacking quantitative structural ensemble models prevents their mechanistic decoding. Here we integrate single-molecule FRET and NMR spectroscopy with molecular simulations showing that DNA binding can lead to complex changes in the IDR ensemble and accessibility. The C-terminal IDR of pioneer factor Sox2 is highly disordered but its conformational dynamics are guided by weak and dynamic charge interactions with the folded DNA binding domain. Both DNA and nucleosome binding induce major rearrangements in the IDR ensemble without affecting DNA binding affinity. Remarkably, interdomain interactions are redistributed in complex with DNA leading to variable exposure of two activation domains critical for transcription. Charged intramolecular interactions allowing for dynamic redistributions may be common in transcription factors and necessary for sensitive tuning of structural ensembles.
Neutron stars are unique laboratories for the investigation of the high density properties of bulk matter. In this work, the astrophysical constraints for a phase transition from hadronic matter to deconfined quark matter are examined thoroughly. A scheme for relating known astrophysical observables such as mass, radius and tidal deformability to the parameter space of such a transition is devised and applied to the set of data currently available.
In order to span a wide parameter space, a highly parameterizable relativistic mean field equation in compliance with chiral effective field theory results is used, where the stiffness of the equation of state can be varied via the effective mass at saturation density. The phase transitions are modelled using a Maxwell construction and assumed to be of first order, with a constant speed of sound quark matter model. The resulting equations of state are analyzed and divided into four categories, which can be used to constrain the parameter space that allows phase transition. It is highlighted, that a subset of this parameter space would even be detectable without the need of higher precision measurements. A phase transition at high densities is shown to be particularly promising in this regard. Finally, the groundwork is laid to apply the equation of state used in this work for supernova or merger simulations, by extending it to non-zero temperatures.
Asymptotic giant branch (AGB) stars are responsible for the production of the main component of the solar s-process distribution. Despite enormous progress in the theoretical modeling of these objects over the last few decades, many uncertainties remain. The still-unknown mechanism leading to the production of 13C neutron source is one example. The nucleosynthetic signature of AGB stars can be examined in a number of stellar sources, from spectroscopic observations of intrinsic and extrinsic stars to the heavy-element isotopic composition of presolar grains found in meteorites. The wealth of available observational data allows for constraining the processes occurring in AGB interiors. In this view, we discuss recent results from new AGB models including the effects of mixing triggered by magnetic fields, and show comparisons of the related s-process nucleosynthesis with available observations.
We investigate the development of the directed, v1, and elliptic flow, v2, in heavy ion collisions in mid-central Au+Au reactions at Elab=1.23A GeV. We demonstrate that the elliptic flow of hot and dense matter is initially positive (v2>0) due to the early pressure gradient. This positive v2 transfers its momentum to the spectators, which leads to the creation of the directed flow v1. In turn, the spectator shadowing of the in-plane expansion leads to a preferred decoupling of hadrons in the out-of-plane direction and results in a negative v2 for the observable final state hadrons. We propose a measurement of v1−v2 flow correlations and of the elliptic flow of dileptons as methods to pin down this evolution pattern. The elliptic flow of the dileptons allows then to determine the early-state EoS more precisely, because it avoids the strong modifications of the momentum distribution due to shadowing seen in the protons. This opens the unique opportunity for the HADES and CBM collaborations to measure the Equation-of-State directly at 2-3 times nuclear saturation density.
Transient receptor potential (TRP) ion channels are among the most well-studied classes of temperature-sensing molecules. Yet, the molecular mechanism and thermodynamic basis for the temperature sensitivity of TRP channels remains to this day poorly understood. One hypothesis is that the temperature-sensing mechanism can simply be described by a difference in heat capacity between the closed and open channel states. While such a two-state model may be simplistic it nonetheless has descriptive value, in the sense that it can be used to to compare overall temperature sensitivity between different channels and mutants. Here, we introduce a mathematical framework based on the two-state model to reliably extract temperature-dependent thermodynamic potentials and heat capacities from measurements of equilibrium constants at different temperatures. Our framework is implemented in an open-source data analysis package that provides a straightforward way to fit both linear and nonlinear van ‘t Hoff plots, thus avoiding some of the previous, potentially erroneous, assumptions when extracting thermodynamic variables from TRP channel electrophysiology data.
The amplification of magnetic fields plays an important role in explaining numerous astrophysical phenomena associated with binary neutron star mergers, such as mass ejection and the powering of short gamma-ray bursts. Magnetic fields in isolated neutron stars are often assumed to be confined to a small region near the stellar surface, while they are normally taken to fill the whole star in numerical modeling of mergers. By performing high-resolution, global, and high-order general-relativistic magnetohydrodynamic simulations, we investigate the impact of a purely crustal magnetic field and contrast it with the standard configuration consisting of a dipolar magnetic field with the same magnetic energy but filling the whole star. While the crust configurations are very effective in generating strong magnetic fields during the Kelvin–Helmholtz-instability stage, they fail to achieve the same level of magnetic-field amplification of the full-star configurations. This is due to the lack of magnetized material in the neutron-star interiors to be used for further turbulent amplification and to the surface losses of highly magnetized matter in the crust configurations. Hence, the final magnetic energies in the two configurations differ by more than 1 order of magnitude. We briefly discuss the impact of these results on astrophysical observables and how they can be employed to deduce the magnetic topology in merging binaries.
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.
The thermodynamic properties of the interacting particle–antiparticle boson system at high temperatures and densities were investigated within the framework of scalar and thermodynamic mean-field models. We assume isospin (charge) density conservation in the system. The equations of state and thermodynamic functions are determined after solving the self-consistent equations. We study the relationship between attractive and repulsive forces in the system and the influence of these interactions on the thermodynamic properties of the bosonic system, especially on the development of the Bose–Einstein condensate. It is shown that under “weak” attraction, the boson system has a phase transition of the second order, which occurs every time the dependence of the particle density crosses the critical curve or even touches it. It was found that with a “strong” attractive interaction, the system forms a Bose condensate during a phase transition of the first order, and, despite the finite value of the isospin density, these condensate states are characterized by a zero chemical potential. That is, such condensate states cannot be described by the grand canonical ensemble since the chemical potential is involved in the conditions of condensate formation, so it cannot be a free variable when the system is in the condensate phase.
This thesis is concerned with the investigation of static and dynamic properties of quantum Heisenberg paramagnets in the absence of a magnetic field and therefore for vanishing magnetization. For this purpose a new formulation of the spin functional renormalization group (SFRG) is employed. The first manifestations of the SFRG were developed by Krieg and Kopietz, motivated by the FRG approach to ordinary field theories and the older works of Vaks, Larkin and Pikin on diagrammatic methods for spin operators.
The main idea is to study quantum spin systems by considering the evolution of correlation functions under a continuous deformation of the interaction between magnetic moments, starting from a solvable limit. This leads to nonperturbative results for quantities like the spin-spin correlation function. After a basic introduction to the phenomena and concomitant problems discussed in this thesis, a detailed description of the SFRG method in its initial formulation is given in the second chapter. We start with the generating functional of connected imaginary-time spin-correlation functions GΛ [h], for which an exact flow equation is derived. A particular issue, already pointed out by Krieg and Kopietz, arises here, namely the singular non-interacting limit of its subtracted Legendre transform ΓΛ [m]. As a consequence the initial condition of that functional does not have a proper series expansion in powers of m. This prevents us from working directly within a pure one-particle irreducible (1-PI) parametrization of the correlation functions, as is often done in the context of field theories. Thus motivated, we develop a workaround explicitly tailored to paramagnets, which provides us with a functional that has a well-behaved Legendre transform. The new approach is based on a different treatment of fluctuations at zero and finite frequencies, analogous to a previous hybrid formulation for the symmetry-broken phase. Certain properties, considered to be highly relevant for isotropic paramagnets, as well as previous observations, already made in the study of simpler spin systems like the Ising model, serve as additional justifications for choosing this construction.
In the third chapter our new method is assessed by calculating the dynamic susceptibility G(k, iω) and thus the dynamic structure factor S(k, ω) in the symmetric phase. For this purpose an approximate integral equation for the dynamic polarization function Π̃(k, iω) was derived. This equation results from a truncation of the hierarchy of flow equations and contains static quantities, that are assumed to be known from another source. Our first application is the high-temperature limit T → ∞ in d ≤ 3 dimensions. Salient features, believed to be part of the spin dynamics in isotropic Heisenberg magnets are also exhibited by our solution, like (anomalous) diffusion in a suitable hydrodynamic limit. Moreover we obtain the same order of magnitude for the diffusion coefficient D as in experiments and other theoretical calculations. Other aspects do not entirely agree with previous approaches.
Afterwards we continue by investigating systems close to the critical point Tc. Dynamic scaling forms for Π̃(k, iω) and S(k, ω), which, like spin diffusion, are postulated on the basis of quite general physical arguments, are reproduced. Agreement of the line-shapes 2with neutron scattering experiments at T = Tc is found to be satisfying, with deviations for ω → 0, that may be attributed to the simplicity of the approximation, like at infinite temperature.
Finally, we focus our attention on the thermodynamic properties of isotropic Heisenberg paramagnets by calculating the static susceptibility G(k). For this purpose we employ simple truncation schemes of the flow equations for the static self-energy ΣΛ (k) and four-spin vertex ΓΛ , together with a basic ansatz for the dynamic polarization Π̃(k, iω) in quantum systems. As a result we obtain transition temperatures Tc of three-dimensional nonfrustrated magnets within an accuracy of 5 percent compared to established benchmark values from Quantum Monte Carlo and high temperature expansion series. We conclude this chapter by giving an outlook on the application of our method to frustrated systems, which may require a combined non-trivial calculation of static and dynamic properties.
The main focus of this thesis is the application of the nonperturbative Functional Renormalization Group (FRG) to the study of low-energies effective models for Quantum Chromodynamics (QCD). The study of effective field theories and models is crucial for our understanding of physics, especially when we deal with fundamental interaction theories like QCD. In particular, the ultimate goal is the understanding of the critical properties of these models in such a way that we can have an insight on the actual critical phenomena of QCD, with a special focus on its chiral phase transition. The choice of the FRG method derives from the fact that it belongs to the class of functional non-perturbative methods and has also the advantage of linking physics at different energy scales. These features make FRG perfectly compatible with the task of studying non-perturbative phenomena and in particular phase transitions, like the ones expected for strongly interacting matter. However, the functional nature of the FRG approach and of the Wetterich equation has a consequence that its exact resolution is hardly possible, and an ansatz for the effective action is generally needed. In this work we choose to adopt the local-potential approximation (LPA), which prescribes to stop at zeroth order in the expansion in derivative operators of the quantum effective action, including only the quantum effective potential. In this work we exploited the key observation that the FRG flow equation can be cast, for specific models and truncation schemes, in the form of an advection-diffusion, possibly with a source term. This type of equation belongs to the class of problems faced in the context of viscous hydrodynamics. Therefore, an innovative approach to the solution of the FRG flow equation consists in the choice of a method developed specifically for the resolution of this class of hydrodynamic equations. In particular, the Kurganov-Tadmor finite-volume scheme is adopted. Throughout this work we apply this scheme to the study of different physical systems, showing the reliability and the flexibility of this approach.
In the first part of the thesis, we discuss the well-known O(N) model, using the hydrodynamic formulation to solve the FRG flow equation in the LPA truncation. We focus on the study of the critical behaviour of the system and calculate the corresponding critical exponents. Particular attention is given to the error estimation in the extraction of critical exponents, which is a needed and not widely explored aspect. The results are well compatible with others in the literature, obtained with different perturbative and nonperturbative methods, which validates the procedure. In the second part of the thesis, we introduce the quark-meson model as a low-energy effective model for QCD, with a specific focus on its chiral symmetry-breaking pattern and the subsequent dynamical quark-mass generation. The LPA flow equation is of the advection-diffusion type, with an extra source contribution which is due to the inclusion of fermionic degrees of freedom. We thus adopt the developed numerical techniques to derive the phase diagram of the model, which is in agreement with the one obtained with other techniques in the literature.
We also follow another possible way for the study of the critical properties of the quark-meson model: the so-called thermodynamic geometry. This approach is based on the interpretation of the parameter space of the system as a differential manifold. One can then obtain relevant information about the phase transitions from the Ricci scalar. We studied the chiral crossover investigating the behavior of the Ricci scalar up to the critical point, featuring a peaking behavior in the presence of the crossover. We then repeated this analysis in the chiral limit, where the phase transition is expected to be of second order. Via this geometric technique it is possible to have a different view on the chiral phase transition of QCD. This is the case since this approach is based on the calculation of quantities which are influenced by higher-order momenta of the thermodynamic potential, thus allowing for a more comprehensive analysis of the phase transition.
Finally, we exploit the numerical advancement to face the issue of the regulator choice in the FRG calculations. This is one of the most delicate issues which arise when using approximations to solve the FRG flow equation and deserves extensive investigation. In particular, we performed a vacuum parameter study and used the RG consistency requirement to determine the impact of the choice of the regulator on the physical observables and on the phase diagram of the model. Via this study we develop a systematic method to comparison the results obtained via different regulators. We show the importance of the choice of an appropriate UV cutoff in the determination of UV-independent IR observables and, consequently, the impact on the latter that the truncation of the effective average action and the choice of the regulator have.
Die Arbeit behandelt die Messung von Photonen mit Teilchendetektoren, die auf digitalen Silizium-Pixelsensoren basieren. Diskutiert werden zwei wesentliche Schritte in den Upgrade-Programmen des ALICE-Experiments am CERN-LHC:
1. FOCAL-Detektor-Upgrade (2027): Untersuchung der Detektorantwort des elektromagnetischen Pixel-Kalorimeters EPICAL-2 und der Form elektromagnetischer Schauer durch Teststrahl-Messungen und Monte Carlo Simulationen.
2. ALICE 3-Upgrade (2035): Simulationsstudien zum Untergrund in der Messung von Photonen mit sehr kleinem Transversalimpuls.
Teil 1: Performance des elektromagnetischen Pixel-Kalorimeters EPICAL-2
Detektordesign und Testmessungen: EPICAL-2, ein SiW-Sandwich-Design-Kalorimeter mit ALPIDE Sensoren, besitzt eine Tiefe von ca. 20 Strahlungslängen und etwa 25 Millionen Pixel. Testmessungen wurden an der Universität Utrecht (kosmische Myonen) sowie am DESY und CERN-SPS (Elektronen) durchgeführt.
Simulation und Validierung: Das EPICAL-2 wird im Simulationspaket Allpix2 implementiert, um die Testmessungen zu validieren und das Detektorverhalten zu untersuchen. Systematische Variationen bestätigen die Stabilität und Reproduzierbarkeit der Simulation.
Datenaufbereitung und Schauerprofile: Im Rahmen der Datenanalyse werden fehlerhafte Pixel ausgeschlossen, Pixel-Treffer zu Clustern gruppiert, Chips kalibriert und der Strahlwinkel korrigiert. Das longitudinale Profil elektromagnetischer Schauer zeigt, dass das Schauermaximum in der Simulation etwas tiefer liegt als in den Testdaten, was auf zusätzliches Material oder eine unvollständige Beschreibung der Schauerentwicklung in der Simulation zurückzuführen sein könnte. Das laterale Profil zeigt, dass eine Schauertrennung im Millimeter-Bereich möglich ist.
Energieantwort und -auflösung: Die nicht-lineare Energieantwort wird sowohl in Testdaten als auch in Simulationen beobachtet. Die Energieauflösung des EPICAL-2 für Cluster ist besser als für Pixeltreffer und vergleichbar mit dem analogen CALICE-Prototypen. Simulationen ohne Strahlenergie-Fluktuationen zeigen eine bessere Energieauflösung als in den Testdaten.
Teil 2: Untergrund in der Messung von Photonen in ALICE 3
Simulationssetup: Die ALICE 3-Detektorgeometrie wird in GEANT4 implementiert, um den Untergrund in der Messung weicher Photonen zu untersuchen. Simulationen mit PYTHIA und GEANT4 zeigen, dass der Untergrund hauptsächlich aus Zerfallsphotonen und Photonen aus externer Bremsstrahlung besteht.
Ergebnisse der Untergrundstudien: Der Untergrund durch Photonen aus externer Bremsstrahlung dominiert und liegt im Akzeptanzbereich des FCT um einen Faktor von 5 bis 10 über dem theoretischen Signal weicher Photonen. In der Simulation wird das Material zu 8%—14% X0 in ALICE 3 bestimmt, wobei bereits bei 5% X0 der Untergrund genauso stark ist wie das erwartete Signal.
Möglichkeiten zur Untergrundreduzierung: Untersuchungen zeigen, dass ein Elektron-Veto das Signal-zu-Untergrund-Verhältnis um den Faktor 30 verbessern und eine Materialreduktion durch ein optimiertes Strahlrohr um den Faktor 7.
Die Ergebnisse des ersten Teils dieser Arbeit demonstrieren insgesamt die gute Performance des EPICAL-2 in Bezug auf die Energiemessung und die Bestimmung der Schauerform. Darüber hinaus unterstützen sie den Einsatz digitaler Kalorimeter im FOCAL-Upgrade des ALICE-Experiments und zeigen das Potenzial der digitalen Kalorimetertechnologie für zukünftige Hochenergiephysik-Experimente.
Die Ergebnisse des zweiten Teils dieser Arbeit liefern einen wesentliche Beitrag zum geplanten ALICE 3-Upgrade. Weiterhin veranschaulichen sie, wie ein Elektron-Veto und die Reduzierung des Materials zusammen eine vielversprechende Messstrategie bilden können.
IHMCIF: an extension of the PDBx/mmCIF data standard for integrative structure determination methods
(2024)
IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.
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 assess 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 modeling of the X-ray emission from the surface of rotating stars.
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.
The spike protein of SARS-CoV-2 is a highly flexible membrane receptor that triggers the translocation of the virus into cells by attaching to the human receptors. Like other type I membrane receptors, this protein has several extracellular domains connected by flexible hinges. The presence of these hinges results in high flexibility, which consequently results in challenges in defining the conformation of the protein. Here, We developed a new method to define the conformational space based on a few variables inspired by the robotic field’s methods to determine a robotic arm’s forward kinematics. Using newly performed atomistic molecular dynamics (MD) simulations and publicly available data, we found that the Denavit-Hartenberg (DH) parameters can reliably show the changes in the local conformation. Furthermore, the rotational and translational components of the homogenous transformation matrix constructed based on the DH parameters can identify the changes in the global conformation of the spike and also differentiate between the conformation with a similar position of the spike head, which other types of parameters, such as spherical coordinates, fail to distinguish between such conformations. Finally, the new method will be beneficial for looking at the conformational heterogeneity in all other type I membrane receptors.
Im Rahmen dieser Doktorarbeit werden drei Schwerpunkte behandelt: 1) Die hocheffektive Beschleunigung von Elektronen und Protonen durch die Wechselwirkung von relativistischen Laserpulsen mit Schäumen. 2) Die Erzeugung und Messung hochintensiver Betatronstrahlung von direkt laserbeschleunigten (DLA-) Elektronen. 3) Die Anwendung von DLA-Elektronen für den biologischen FLASH-Effekt mit einer rekordbrechenden Dosisrate.
Die direkte Laserbeschleunigung von Elektronen wurde durch die Wechselwirkung eines sub-ps-Laserpulses mit einer Intensität von ~ 10^19 W/cm^2 mit einem Plasma nahe kritischer Elektronendichte (NCD) untersucht. Ein sub-mm langes NCD-Plasma wurde durch Erhitzen eines Schaums mit einer niedrigen Dichte mit einem ns-Puls von 10^13-10^14 W/cm^2 erzeugt. Die Experimente wurden an der PHELIX-Anlage (Petawatt Hoch- Energie Laser für Schwerionenexperimente) in den Jahren 2019 – 2023 durchgeführt. Während der Suche nach optimalen Bedingungen für die Beschleunigung von Elektronen und Protonen wurden die Parameter des ns-Pulses variiert und verschiedene Targets verwendet. Es wurde gezeigt, dass das Plasma im Schaum gute Voraussetzungen für die Erzeugung gerichteter, ultrarelativistischer DLA-Elektronen mit Energien von bis zu 100 MeV bietet. Die Elektronen weisen eine Boltzmann-ähnliche Energieverteilung mit einer Temperatur von 10-20 MeV auf.
Optimale Bedingungen für eine effektive Beschleunigung von DLA-Elektronen wurden bei der Kombination eines CHO-Schaums mit einer Dichte von 2 mg/cm3 und einer Dicke von 300-500 µm mit einer Metallfolie erreicht. Die Gesamtladung der detektierten Elektronen mit Energien über 1,5 MeV erreichte 0,5-1 µC mit der Umwandlungseffizienz der Laserenergie von ~ 20-30%.
Außerdem wird die Beschleunigung von Protonen durch DLA-Elektronen anders verursacht als bei typischer Target Normal Sheath Acceleration (TNSA). Für die Untersuchung der lokalen Protonenenergieverteilung wurden Magnetspektrometer unter verschiedenen Winkeln zur Laserachse verwendet. Dafür wurde eine Filtermethode entwickelt, welche es ermöglicht, Spektren von Protonen mit Energien von bis zu 100 MeV zu rekonstruieren. Es wurde gezeigt, dass am PHELIX durch die Kombination von einem ~ 300-400 µm dicken CHO-Schaum mit einer Dichte von 2 mg/cm^3 und einer 10 µm dicken Au-Folie bei einer Intensität des sub-ps-Pulses von ~ 10^19 W/cm^2 und unter Verwendung eines optimierten ns-Vorpulses eine optimale Protonenbeschleunigung erreicht wurde. Es wurde ein TNSA-ähnliches Regime mit einer maximalen Cut-off-Energie von 34±0,5 MeV beobachtet. Im Vergleich dazu wurde bei der typischen TNSA unter Verwendung einer 10 µm dicken Au-Folie als Target und derselben Laserintensität eine maximale Cut-off-Energie von 24±0,5 MeV gemessen. Darüber hinaus beobachteten wir einen sehr schwachen Abfall der Protonenanzahl in Abhängigkeit von der Protonenenergie (anders als bei der typischen TNSA) und eine sehr regelmäßige Protonenstrahlverteilung in einem breiten Winkelbereich bis zu hohen Energien. Dies könnte zur Verbesserung der Qualität der Protonenradiographie von Plasmafeldern genutzt werden.
Beim DLA-Prozess (im NCD-Plasma) entsteht Betatronstrahlung durch die Oszillationen von Elektronen in quasi-statischen elektrischen und magnetischen Feldern des Plasmakanals. Um diese Strahlung zu untersuchen, wurde ein neues modifiziertes Magnetspektrometer (X-MS) konstruiert. Das X-MS ermöglicht die 1D-Auflösung mehrerer Quellen. Dank dieser Spezifikation war es möglich, Betatronstrahlung von Bremsstrahlung der ponderomotorischen Elektronen im Metallhalter zu trennen und zu messen.
Im Experiment mit einem CHO-Schaum mit einer Dichte von 2 mg/cm^3 und einer Dicke von ~ 800 µm als Target wurde die von den optimierten DLA-Elektronen erzeugte Betatronstrahlung gemessen. Bei einer Peak-Intensität des dreieckigen ns-Pulses von ~ 3·10^13 W/cm^2 und des sub-ps-Pulses von ~ 10^19 W/cm^2, welcher 4±0,5 ns gegenüber dem ns-Puls verzögert war, betrug der Halbwinkel im FWHM-Bereich des Elektronenstrahls 17±2°. Unter diesen Bedingungen war die Betatronstrahlung mit einem Halbwinkel im FWHM-Bereich von 11±2° für die Photonen mit Energien über 10 keV ebenfalls gerichtet. Die Photonenanzahl mit Energien über 10 keV wurde auf etwa 3·10^10 / 3·10^11 (gerichtete Photonen / Photonen im Halbraum entlang der Laserstrahlrichtung) abgeschätzt. Die maximale Photonenanzahl pro Raumwinkel betrug ~2·10^11 photons/sr. Die Brillanz der registrierten Betatronstrahlung erreichte ~ 2·10^20 photons/s/mm^2/mrad^2/(0.1% BW) bei 10 keV.
Die Verwendung eines Hochstromstrahls aus DLA-Elektronen für die FLASH-Strahlentherapie ermöglicht das Erreichen einer Dosis von bis zu 50-70 Gy während eines sub-ps-Laserpulses. Im Jahr 2021, während der P213-Strahlzeit am PHELIX wurde der Sauerstoffkonzentrationsabfall bei der Bestrahlung von Medien (Wasser und andere biologische Medien) mit DLA-Elektronen in Abhängigkeit von der Dosis untersucht. Die Strahlendosis wurde hierbei indirekt gemessen. Hierfür wurde eine Rekonstruktionsmethode entwickelt, die es ermöglicht, die Dosis innerhalb des „Wasser-Containers“ auf Basis von Messungen außerhalb des Containers mit einem untersuchten Medium zu ermitteln. Es wurde eine gute Übereinstimmung zwischen dem Experiment und einer Monte-Carlo-Simulation für Wasser gezeigt. Die registrierte Dosisrate erreichte einen Rekordwert von ~ 70 TGy/s.
Das Heidelberger Ionenstrahl-Therapiezentrum (HIT) stellt Protonen-, Helium- und Kohlenstoff-Ionenstrahlen unterschiedlicher Energie und Intensität für die Krebsbehandlung und Sauerstoff-Ionenstrahlen für Experimente zur Verfügung. Der hierfür verwendete Beschleuniger ist darüber hinaus in der Lage auch Ionenstrahlintensitäten unterhalb der für Therapien verwendeten bereitzustellen. Allerdings ist das derzeit installierte Strahldiagnosesystems nicht in der Lage, das Strahlprofil bei solchen geringen Intensitäten (< 10^5 Ionen/s) zu messen. Dabei existieren mögliche medizinische Anwendung für diese niederintensiven Ionen-strahlen, wie beispielsweise eine neuartige und potentiell klinisch vorteilhafte Bildgebung: die Ionenradiographie. Eine essentielle Voraussetzung für diese und andere Anwendungen ist ein System zur Überwachung von Ionenstrahlen niedriger Intensität. Ein solches System wurde im Rahmen dieser Arbeit konzipiert, realisiert, getestet und optimiert.
Das Funktionsprinzip basiert auf szintillierenden Fasern, insbesondere solchen mit erhöhter Strahlungshärte für die Möglichkeit einer dauerhaften Platzierung im Therapiestrahl. Ein diese Fasern durchlaufendes Ion regt den darin enthaltenen Szintillator durch Stoßprozesse kurzzeitig an. Die dabei deponierte Energie wird anschließend in Form von Photonen wieder emittiert. Silizium-Photomultiplier sind an den Enden der Fasern montiert und wandeln die Photonensignale in verstärkte elektrische Impulse um. Diese Impulse werden von einer neuartigen und dedizierten Ausleseelektronik aufgezeichnet und verarbeitet. Ein Prototypaufbau, bestehend aus den genannten Teilen, wurde im Strahl getestet und kann das transversale Strahlprofil erfolgreich im Intensitätsbereich von 10^7 Ionen/s bis hinunter zu 10^2 Ionen/s aufzeichnen. Darüber hinaus konnte, durch die erfolgreiche Ankunftszeitmessung von einzelnen Ionen bis zu Intensitäten von 5*10^4 Ionen/s, ein Machbarkeitsnachweis für die Messung der Spur von einzelnen Teilchen erbracht werden.