Institut für Theoretische Physik und Astrophysik
Refine
Has Fulltext
- yes (138)
Is part of the Bibliography
- yes (138)
Year of publication
Document Type
- Doctoral Thesis (138) (remove)
Keywords
- Monte-Carlo-Simulation (12)
- Supersymmetrie (12)
- Topologischer Isolator (12)
- Blazar (9)
- Aktiver galaktischer Kern (8)
- Elementarteilchenphysik (8)
- Astrophysik (7)
- Gammastrahlung (7)
- Hubbard-Modell (7)
- Starke Kopplung (7)
Institute
Sonstige beteiligte Institutionen
Strong correlations caused by interaction in systems of electrons can bring about unusual physical phenomena due to many-body quantum effects that cannot properly be captured by standard electronic structure methods like density functional theory. In this thesis, we apply the state-of-the-art continuous-time quantum Monte Carlo algorithm in hybridization expansion (CT-HYB) for the strongly correlated multi-orbital Anderson impurity model (AIM) to the solution of models of magnetic impurities on metallic surfaces and, via dynamical mean-field theory (DMFT), to the solution of a lattice model, the multi-orbital Hubbard model with Hund's coupling.
A concise introduction to the theoretical background focuses on information directly relevant to the understanding of applied models, methods, and the interpretation of results. It starts with a discussion of the AIM with its parameters and its solution in the path integral formalism, the basis of the CT-HYB algorithm. We consider its derivation and implementation in some detail before reviewing the DMFT approach to correlated lattice models and the interpretation of the single-particle Green's function.
We review two algorithmic developments for the CT-HYB algorithm that help to increase the performance of calculations especially in case of a complex structure of the interaction matrix and allow the precise calculation of self-energies and vertex functions also at intermediate and higher frequencies.
Our comparative analysis of Kondo screening in the cobalt on copper impurity system points out the importance of an accurate interaction matrix for qualitatively correct Kondo temperatures and the relevance of all d-orbitals in that case. Theoretical modeling of cobalt impurities in copper "atomic wires" fails to reproduce variations and partial absence of Kondo resonances depending on the wire size. We analyze the dependence of results on parameters and consider possible reasons for the discrepancy. Different Kondo temperatures of iron adatoms adsorbed on clean or oxygen-reconstructed niobium in the normal state are qualitatively reproduced, with the adsorption distance identified as major factor and implications for the superconducting state pointed out.
Moving on to lattice problems, we demonstrate the connection between Hund's coupling, shown to cause first-order character of the interaction-driven Mott transition at half-filling in the two-orbital Hubbard model, and a phase separation zone ending in a quantum critical point at finite doping. We touch on similarities in realistic models of iron-pnictide superconductors. We analyze the manifestation of the compressibility divergence at the finite-temperature critical points away from half-filling in the eigenbasis of the two-particle generalized susceptibility. A threshold for impurity susceptibility eigenvalues that indicates divergence of the DMFT lattice compressibility and distinguishes thermodynamic stability and instability of DMFT solutions is determined.
Next to the emergence of nearly isolated quantum systems such as ultracold atoms with unprecedented experimental tunability, the conceptualization of the eigenstate thermalization hypothesis (ETH) by Deutsch and Srednicki in the late 20th century has sparked exceptional interest in the mechanism of quantum thermalization. The ETH conjectures that the expectation value of a local observable within the quantum state of an isolated, interacting quantum system converges to the thermal equilibrium value at large times caused by a loss of phase coherence, referred to as dephasing. The thermal behavior within the quantum expectation value is traced back to the level of individual eigenstates, who locally act as a thermal bath to subsystems of the full quantum system and are hence locally indistinguishable to thermal states. The ETH has important implications for the understanding of the foundations of statistical mechanics, the quantum-to-classical transition, and the nature of quantum entanglement. Irrespective of its theoretical success, a rigorous proof has remained elusive so far. $$ \ $$
An alternative approach to explain thermalization of quantum states is given by the concept of typicality. Typicality deals with typical states \(\Psi\) chosen from a subspace of Hilbert space with energy \(E\) and small fluctuations \(\delta\) around it. It assumes that the possible microstates of this subspace of Hilbert space are uniformly distributed random vectors. This is inspired by the microcanonical ensemble in classical statistical mechanics, which assumes equal weights for all accessible microstates with energy \(E\) within an energy allowance \(\delta\). It follows from the ergodic hypothesis, which states that the time spent in each part of phase space is proportional to its volume leading to large time averages being equated to ensemble averages. In typicality, the Hilbert space of quantum mechanics is hence treated as an analogue of classical phase space where statistical and thermodynamic properties can be defined. Since typicality merely shifts assumptions of statistical mechanics to the quantum realm, it does not provide a complete understanding of the emergence of thermalization on a fundamental microscopic level. $$ \ $$
To gain insights on quantum thermalization and derive it from a microscopic approach, we exclusively consider the fundamental laws of quantum mechanics. In the joint work with T. Hofmann, R. Thomale and M. Greiter, on which this thesis reports, we explore the ETH in generic local Hamiltonians in a two-dimensional spin-\(1/2\) lattice with random nearest neighbor spin-spin interactions and random on-site magnetic fields. This isolated quantum system is divided into a small subsystem weakly coupled to the remaining part, which is assumed to be large and which we refer to as bath. Eigenstates of the full quantum system as well as the action of local operators on those can then be decomposed in terms of a product basis of eigenstates of the small subsystem and the bath. Central to our analysis is the fact that the coupling between the subsystem and the bath, represented in terms of the uncoupled product eigenbasis, is given by an energy dependent random band matrix, which is obtained from both analytical and numerical considerations. $$ \ $$
Utilizing the methods of Dyson-Brownian random matrix theory for random band matrices, we analytically show that the overlaps of eigenstates of the full quantum system with the uncoupled product eigenbasis are described by Cauchy-Lorentz distributions close to their respective peaks. The result is supported by an extensive numerical study using exact diagonalization, where the numerical parameters for the overlap curve agree with the theoretical calculation. The information on the decomposition of the eigenstates of the full quantum system enables us to derive the reduced density matrix within the small subsystem given the pure density matrix of a single eigenstate. We show that in the large bath limit the reduced density matrix converges to a thermal density matrix with canonical Boltzmann probabilities determined by renormalized energies of the small subsystem which are shifted from their bare values due the influence of the coupling to the bath. The behavior of the reduced density matrix is confirmed through a finite size scaling analysis of the numerical data. Within our calculation, we make use of the pivotal result, that the density of states of a local random Hamiltonian is given by a Gaussian distribution under very general circumstances. As a consequence of our analysis, the quantum expectation value of any local observable in the subsystem agrees with its thermal expectation value, which proves the validity of the ETH in the equilibrium phase for the considered class of random local Hamiltonians and elevates it from hypothesis to theory. $$ \ $$
Our analysis of quantum thermalization solely relies on the application of quantum mechanics to large systems, locality and the absence of integrability. With the self-averaging property of large random matrices, random matrix theory does not entail a statistical assumption, but is rather applied as a mathematical tool to extract information about the behavior of large quantum systems. The canonical distribution of statistical mechanics is derived without resorting to statistical assumptions such as the concepts of ergodicity or maximal entropy, nor assuming any characteristics of quantum states such as in typicality. In future research, with this microscopic approach it may become possible to exactly pinpoint the origin of failure of quantum thermalization, e.g. in systems that exhibit many body localization or many body quantum scars. The theory further enables the systematic investigation of equilibration, i.e. to study the time scales on which thermalization takes place.
In this thesis, I establish new relations between quantum information measures in a two-dimensional CFT and geometric objects in a three-dimensional AdS space employing the AdS/CFT correspondence. I focus on two quantum information measures: the computational cost of quantum circuits in a CFT and Berry phases in two entangled CFTs. In particular, I show that these quantities are associated with geometric objects in the dual AdS space.
The hunt for topological materials is one of the main topics of recent research in condensed matter physics. We analyze the 4-band Luttinger model, which considers the total angular momentum \(j = 3/2\) hole states of many semiconductors. Our analysis shows that this model hosts a wide array of topological phases and allows analytical calculations of the related topological surface states. The existence of these surface states is highly desired due to their strong protection against perturbations.
In the first part of the thesis, we predict the existence of either one or two two-dimensional (2D) surface states of topological origin in the three-dimensional (3D) quadratic-node semimetal phase of the Luttinger model, called the Luttinger semimetal phase. We associate the origin of these states with the inverted order of s and p-orbital states in the band structure and approximate chiral symmetry around the node. Hence, our findings are essential for many materials, including HgTe, α-Sn, and iridate compounds. Such materials are often modified with strain engineering by growing the crystal on a substrate with a different lattice constant, which adds a deformation potential to the electrons. While tensile strain is often used to drive such materials into a gapped topological insulator regime, we apply compressive strain to induce a topological semimetal regime. Here, we differentiate between Dirac and Weyl semimetals based on inversion and time-reversal symmetry being simultaneously present or not. One major part of this thesis is the theoretical study of the evolution of the Luttinger semimetal surface states in these topological semimetal phases.
The relative strength of the compressive strain and typical bulk inversion asymmetry (BIA) terms allow the definition of a symmetry hierarchy in the system. The cubic symmetric \(O_h\) Luttinger model is the highest symmetry low-energy parent model. Since the BIA terms in the Weyl semimetal phase are small in most materials, we find a narrow energy and momentum range around the Weyl points where the surface states form Fermi arcs between two Weyl nodes with opposite chirality. Consequently, we see 2D momentum planes between the Weyl points, which can be considered as effective 2D Chern insulators with chiral edge states connecting the valence and conduction band in the bulk gap. Exceeding the range of the BIA terms, the compressive strain becomes dominating, and the system behaves like a Dirac semimetal with two doubly degenerate linear Dirac nodes in the band structure. For energies larger than the compressive strain strength, the quadratic terms in the Luttinger model dominate and surface band structure is indistinguishable from an unperturbed Luttinger semimetal. To conclude this symmetry hierarchy, we analyze the limit of the Luttinger model when the remote \(j = 1/2\)
electron states show a considerable hybridization with the \(j = 3/2\) hole states around the Fermi level. Here, the Luttinger model is not valid anymore and one needs to consider more complicated models, like the 6-band Kane Hamiltonian.
In the second part of this thesis, we analyze theoretically two different setups for s-wave superconductivity proximitized \(j = 3/2\) particles in Luttinger materials under a magnetic field. First, we explore a one-dimensional wire setup, where the intrinsic BIA of inversion asymmetric crystals opens a topological gap in the bulk states. In contrast to wires, modeled by a quadratic dispersion with Rashba or Dresselhaus spin-orbit coupling, we find two topological phase transitions due to the different effects of magnetic fields to \(|j_z| = 3/2\) heavy-hole (HH) and \(|j_z| = 1/2\) light-hole (LH) states. Second, we discuss a two-dimensional Josephson junction setup, where we find Andreev-bound states inside the superconducting gap. Here, the intrinsic spin-orbit coupling of the Luttinger model is sufficient to open a topological gap even in the presence of inversion symmetry. This originates from the hybridization of the light and heavy-hole bands in combination with the superconducting pairing.
Consequently, both setups can form Majorana-bound states at the boundaries of the system.
The existence of these states are highly relevant in the scientific community due to their nonabelian braiding statistics and stability against decoherence, making them a prime candidate for the realization of topological quantum computation. Majorana-bound states form at zero energy and are protected by the topological gap. We predict that our findings of the topological superconductor phase of the Luttinger model are valid for both semimetal and metal phases. Hence, our study is additionally relevant for metallic systems, like p-doped GaAs. This opens a new avenue for the search for topological superconductivity.
Explaining the baryon asymmetry of the Universe has been a long-standing problem of particle physics, with the consensus being that new physics is required as the Standard Model (SM) cannot resolve this issue. Beyond the Standard Model (BSM) scenarios would need to incorporate new sources of \(CP\) violation and either introduce new departures from thermal equilibrium or modify the existing electroweak phase transition. In this thesis, we explore two approaches to baryogenesis, i.e. the generation of this asymmetry.
In the first approach, we study the two-particle irreducible (2PI) formalism as a means to investigate non-equilibrium phenomena. After arriving at the renormalised equations of motions (EOMs) to describe the dynamics of a phase transition, we discuss the techniques required to obtain the various counterterms in an on-shell scheme. To this end, we consider three truncations up to two-loop order of the 2PI effective action: the Hartree approximation, the scalar sunset approximation and the fermionic sunset approximation. We then reconsider the renormalisation procedure in an \(\overline{\text{MS}}\) scheme to evaluate the 2PI effective potential for the aforementioned truncations. In the Hartree and the scalar sunset approximations, we obtain analytic expressions for the various counterterms and subsequently calculate the effective potential by piecing together the finite contributions. For the fermionic sunset approximation, we obtain similar equations for the counterterms in terms of divergent parts of loop integrals. However, these integrals cannot be expressed in an analytic form, making it impossible to evaluate the 2PI effective potential with the fermionic contribution. Our main results are thus related to the renormalisation programme in the 2PI formalism: \( (i) \)the procedure to obtain the renormalised EOMs, now including fermions, which serve as the starting point for the transport equations for electroweak baryogenesis and \( (ii) \) the method to obtain the 2PI effective potential in a transparent manner.
In the second approach, we study baryogenesis via leptogenesis. Here, an asymmetry in the lepton sector is generated, which is then converted into the baryon asymmetry via the sphaleron process in the SM. We proceed to consider an extension of the SM along the lines of a scotogenic framework. The newly introduced particles are charged odd under a \(\mathbb{Z}_2\) symmetry, and masses for the SM neutrinos are generated radiatively. The \(\mathbb{Z}_2\) symmetry results in the lightest BSM particle being stable, allowing for a suitable dark matter (DM) candidate. Furthermore, the newly introduced heavy Majorana fermionic singlets provide the necessary sources of \(CP\) violation through their Yukawa interactions and their out-of-equilibrium decays produce a lepton asymmetry. This model is constrained from a wide range of observables, such as consistency with neutrino oscillation data, limits on branching ratios of charged lepton flavour violating decays, electroweak observables and obtaining the observed DM relic density. We study leptogenesis in this model in light of the results of a Markov chain Monte Carlo scan, implemented in consideration of the aforementioned constraints. Successful leptogenesis in this model, to account for the baryon asymmetry, then severely constrains the available parameter space.
Relativistic effects crucially influence the fundamental properties of many quantum materials. In the accelerated reference frame of an electron, the electric field of the nuclei is transformed into a magnetic field that couples to the electron spin. The resulting interaction between an electron spin and its orbital angular momentum, known as spin-orbit coupling (SOC), is hence fundamental to the physics of many condensed matter phenomena. It is particularly important quantitatively in low-dimensional quantum systems, where its coexistence with inversion symmetry breaking can lead to a splitting of spin degeneracy and spin momentum locking. Using the paradigm of Landau Fermi liquid theory, the physics of SOC can be adequately incorporated in an effective single particle picture. In a weak coupling approach, electronic correlation effects beyond single particle propagator renormalization can trigger Fermi surface instabilities such as itinerant magnetism, electron nematic phases, superconductivity, or other symmetry broken states of matter.
In this thesis, we use a weak coupling-based approach to study the effect of SOC on Fermi surface instabilities and, in particular, superconductivity. This encompasses a weak coupling renormalization group formulation of unconventional superconductivity as well as the random phase approximation. We propose a unified formulation for both of these two-particle Green’s function approaches based on the notion of a generalized susceptibility.
In the half-Heusler semimetal and superconductor LuPtBi, both SOC and electronic correlation
effects are prominent, and thus indispensable for any concise theoretical description. The metallic and weakly dispersive surface states of this material feature spin momentum locked Fermi surfaces, which we propose as a possible domain for the onset of unconventional surface superconductivity. Using our framework for the analysis of Fermi surface instability and combining it with ab-initio density functional theory calculations, we analyse the surface band structure of LuPtBi, and particularly its propensity towards Cooper pair formation. We study how the presence of strong SOC modifies the classification of two-electron wave functions as well as the screening of electron-electron interactions. Assuming an electronic mechanism, we identify a chiral superconducting condensate featuring Majorana edge modes to be energetically favoured over a wide range of model parameters.
The last years have witnessed an exciting scientific quest for intriguing topological phenomena in time-dependent quantum systems. A key to many manifestations of topology in dynamical systems relies on the effective dimensional extension by time-periodic drives. An archetypal example is provided by the Thouless pump in one spatial dimension, where a robust and quantized charge transport can be described in terms of an integer quantum Hall effect upon interpreting time as an extra dimension. Generalizing this fundamental concept to multifrequency driving, a variety of higher-dimensional topological models can be engineered in dynamical synthetic dimensions, where the underlying topological classification leads to quantized pumping effects in the associated lower-dimensional time-dependent systems.
In this Thesis, we explore how correlations profoundly impact the topological features of dynamical synthetic quantum materials. More precisely, we demonstrate that the interplay of interaction and dynamical synthetic dimension gives rise to striking topological phenomena that go beyond noninteracting implementations. As a starting point, we exploit the Floquet counterpart of an integer quantum Hall scenario, namely a two-level system driven by two incommensurate frequencies. In this model, the topologically quantized response translates into a process in which photons of different frequencies are exchanged between the external modes, referred to as topological frequency conversion. We extend this prototypical setup to an interacting version, focusing on the minimal case of two correlated spins equally exposed to the external drives. We show that the topological invariant determining the frequency conversion can be changed by odd integers, something explicitly forbidden in the noninteracting limit of two identical spins. This correlated topological feature may, in turn, result in an enhancement of the quantized response.
Robust response signals, such as those predicted for the topological frequency converter, are of fundamental interest for potential technological applications of topological quantum matter. Based on an open quantum system implementation of the frequency converter, we propose a novel mechanism of topological quantization coined ''topological burning glass effect''. Remarkably, this mechanism amplifies the local response of the driven two-level system by an integer that is proportional to the number of environmental degrees of freedom to which the system is strongly coupled. Specifically, our findings are illustrated by the extension of the frequency converter to a central spin model. There, the local energy transfer mediated exclusively by the central spin is significantly enhanced by the collective motion of the surrounding spins. In this sense, the central spin adopts the topological nature of the total system in its non-unitary dynamics, taking into account the correlations with the environment.
Magnetism is a phenomenon ubiquitously found in everyday life. Yet, together with superconductivity and superfluidity, it is among the few macroscopically realized quantum states. Although well-understood on a quasi-classical level, its microscopic description is still far from being solved. The interplay of strong interactions present in magnetic condensed-matter systems and the non-trivial commutator structure governing the underlying spin algebra prevents most conventional approaches in solid-state theory to be applied.
On the other hand, the quantum limit of magnetic systems is fertile land for the development of exotic phases of matter called spin-liquids. In these states, quantum fluctuations inhibit the formation of magnetic long-range order down to the lowest temperatures. From a theoretical point of view, spin-liquids open up the possibility to study their exotic properties, such as fractionalized excitations and emergent gauge fields. However, despite huge theoretical and experimental efforts, no material realizing spin-liquid properties has been unambiguously identified with a three-dimensional crystal structure. The search for such a realization is hindered by the inherent difficulty even for model calculations. As most numerical techniques are not applicable due to the interaction structure and dimensionality of these systems, a methodological gap has to be filled.
In this thesis, to fill this void, we employ the pseudo-fermion functional renormalization group (PFFRG), which provides a scheme to investigate ground state properties of quantum magnetic systems even in three spatial dimensions.
We report the status quo of this established method and extend it by alleviating some of its inherent approximations. To this end, we develop a multi-loop formulation of PFFRG, including hitherto neglected terms in the underlying flow equations consistently, rendering the outcome equivalent to a parquet approximation. As a necessary prerequisite, we also significantly improve the numerical accuracy of our implementation of the method by switching to a formulation respecting the asymptotic behavior of the vertex functions as well as employing state-of-the-art numerical algorithms tailored towards PFFRG. The resulting codebase was made publicly accessible in the open-source code PFFRGSolver.jl.
We subsequently apply the technique to both model systems and real materials. Augmented by a classical analysis of the respective models, we scan the phase diagram of the three-dimensional body-centered cubic lattice up to third-nearest neighbor coupling and the Pyrochlore lattice up to second-nearest neighbor. In both systems, we uncover in addition to the classically ordered phases, an extended parameter regime, where a quantum paramagnetic phase appears, giving rise to the possibility of a quantum spin liquid.
Additionally, we also use the nearest-neighbor antiferromagnet on the Pyrochlore lattice as well as the simple cubic lattice with first- and third-nearest neighbor couplings as a testbed for multi-loop PFFRG, demonstrating, that the inclusion of higher loop orders has quantitative effects in paramagnetic regimes and that the onset of order can be signaled by a lack of loop convergence.
Turning towards material realizations, we investigate the diamond lattice compound MnSc\(_2\)S\(_4\), explaining on grounds of ab initio couplings the emergence of a spiral spin liquid at low temperatures, but above the ordering transition.
In the Pyrochlore compound Lu\(_2\)Mo\(_2\)O\(_5\)N\(_2\), which is known to not magnetically order down to lowest temperatures, we predict a spin liquid state displaying a characteristic gearwheel pattern in the spin structure factor.
Context. In active galaxies, matter is accreted onto super massive black holes (SMBH). This accretion process causes a region roughly the size of our solar system to outshine the entire host galaxy, forming an active galactic nucleus (AGN). In some of these active galaxies, highly relativistic particle jets are formed parallel to the rotation axis of the super massive black hole. A fraction of these sources is observed under a small inclination angle between the pointing direction of the jet and the observing line of sight. These sources are called blazars. Due to the small inclination angle and the highly relativistic speeds of the particles in the jet, beaming effects occur in the radiation of these particles. Blazars can be subdivided into the high luminosity flat spectrum radio quasars (FSRQs) and the low luminosity BL Lacertae objects (BL Lacs). As all AGN, blazars are broadband emitters and therefore observable from the longest wavelengths in the radio regime to the shortest wavelengths in the gamma-ray regime. In this thesis I will analyze blazars at these two extremes with respect to their parsec-scale properties in the radio and their time evolution properties in gamma-ray flux.
Method. In the radio regime the technique of very long baseline interferometry (VLBI) can be used in order to spatially resolve the synchrotron radiation coming from those objects down to sub-parsec scales. This information can be used to observe the time evolution of the structure of such sources. This is done in large monitoring programs such as the MOJAVE (15 GHz) and the Boston University blazar monitoring program (43 GHz). In this thesis I utilize data of 28 sources from these monitoring programs spanning 10 years of observation from 2003 to 2013, resulting in over 1800 observed epochs, to study the brightness temperature and diameter gradients of these jets. I conduct a search for systematic geometry transitions in the radio jets. The synchrotron cooling time in the radio core of the jets is used to determine the magnetic field strength in the radio core. Considering the jet geometry, these magnetic field strengths are scaled to the ergosphere of the SMBH in order to obtain the distance of the radio core to the SMBH.
In the gamma-regime these blazars cannot be spatially resolved. Due to this, it is hard to put strong constrains onto where the gamma-ray emitting region is. Blazars have shown to be variable at high energies on time scales down to minutes. The nature of this variability can be studied in order to put constrains on the particle acceleration mechanism and possibly the region and size of the gamma-ray emitting region. The variability of blazars in the energy range between 0.1 GeV and 1 GeV can for example be observed with the pair-conversion telescope on board the Fermi satellite. I use 10 years of data from the Fermi-LAT (LAT: Large Area Telescope) satellite in order to study the variability of a large sample of blazars (300-800, depending on the used significance filters for data points). I quantify this variability with the Ornstein-Uhlenbeck (OU) parameters and the power spectral density (PSD) slopes. The same procedure is applied to 20 light curves available for the radio sample.
Results. The diameter evolution along the jet axis of the radio sources suggests, that FSRQs feature flatter gradients than BL Lacs. Fitting these gradients, it is revealed that BL Lacs are systematically better described by a simple single power law than FSRQs. I found 9 sources with a strongly constrained geometry transition. The sources are 0219+421, 0336-019, 0415+379, 0528+134, 0836+710, 1101+384, 1156+295, 1253-055 and 2200+420. In all of these sources, the geometry transition regions are further out in the jet than the Bondi sphere. The magnetic field strengths of BL Lacs is systematically larger than that of FSRQs. However the scaling of these fields suggest that the radio cores of BL Lac objects are closer to the SMBHs than the radio cores of FSRQs. Analyzing the variability of Fermi-LAT light curves yields consistent results for all samples. FSRQs show systematically steeper PSD slopes and feature OU parameters more favorable to strong variability than BL Lacs. The Fermi-LAT light curves of the sub-sample of radio jets, suggest an anticorrelation between the jet complexity from the radio observations and the OU-parameters as well as the PSD slopes from the gamma-ray observations.
Conclusion.
The flatter diameter gradients of FSRQs suggest that these sources are more collimated further down the jet than BL Lacs. The systematically better description of the diameter and brightness temperature gradient by a single power law of BL Lacs, suggest that FSRQs are more complex with respect to the diameter evolution along the jet and the surface brightness distribution than BL Lac objects. FSRQs often feature regions where recollimation can occur in distinct knots within the jets. For the sources where a geometry transition could be constrained, the Bondi radius, being systematically smaller than the position of the transition region along the jet axis, suggest that changing pressure gradients are not the sole cause for these systematic geometry transitions. Nevertheless they may be responsible for recollimation regions, found typically downstream the jet, beyond the Bondi radius and the transition zone. The difference in the distance of the radio cores between FSRQs and BL Lacs is most likely due to the combination of differences in SMBH masses and systematically smaller jet powers in BL Lacs. The variability in energy ranges above 100 MeV and above 1 GeV-regime suggest that many light curves of BL Lac objects are more likely to be white noise while the PSD slopes and the OU parameters of FSRQ gamma-ray light curves favor stronger variability on larger time scales with respect to the time binning of the analyzed light curve. Although the anticorrelation of the jet complexity acquired from the radio observations and the PSD slopes and OU parameters from the gamma-observations suggest that more complex sources favor OU parameters and PSD slopes resulting in more variability (not white noise) it is beyond the scope of this thesis to pinpoint whether this correlation results from causation. The question whether a complex jet causes more gamma-ray variability or more gamma-ray variability causes more complex jets cannot be answered at this point. Nevertheless the computed correlation measures suggest that this dependence is most likely not linear and therefore an indication that these effects might even interact.
The extragalactic gamma-ray sky is dominated by blazars, active galactic nuclei (AGN) with a relativistic jet that is closely aligned with the line of sight. Galaxies develop an active nucleus if the central supermassive black hole (BH) accretes large amounts of ambient matter and magnetic flux. The inflowing mass accumulates around the plane perpendicular to the accretion flow's angular momentum. The flow is heated through viscous friction and part of the released energy is radiated as blackbody or non-thermal radiation, with luminosities that can dominate the accumulated stellar luminosity of the host galaxy. A fraction of the accretion flow luminosity is reprocessed in a surrounding field of ionised gas clouds. These clouds, revolving around the central BH, emit Doppler-broadened atomic emission lines. The region where these broad-line-emitting clouds are located is called broad-line region (BLR).
About one in ten AGN forms an outflow of radiation and relativistic particles, called a relativistic jet. According to the Blandford-Znajek mechanism, this is facilitated through electromagnetic processes in the magnetosphere of a spinning BH. The latter induces a magnetospheric poloidal current circuit, generating a decelerating torque on the BH and inducing a toroidal magnetic field. Consequently, rotational energy of the BH is converted to Poynting flux streaming away mainly along the rotational axis and starting the jet. One possibility for particle acceleration near the jet base is realised by magnetospheric vacuum gaps, regions temporarily devoid of plasma, such that an intermittent electric field arises parallel to the magnetic field lines, enabling particle acceleration and contributing to the mass loading of the jets.
Magnetised structures, containing bunches of relativistic electrons, propagate away from the galactic nucleus along the jets. Assuming that these electrons emit synchrotron radiation and that they inverse-Compton (IC) up-scatter abundant target photons, which can either be the synchrotron photons themselves or photons from external emitters, the emitted spectrum can be theoretically determined. Additionally taking into account that these emission regions move relativistically themselves and that the emission is Doppler-boosted and beamed in forward direction, the typical two-hump spectral energy distribution (SED) of blazars is recovered.
There are however findings that challenge this well-established model. Short-time variability, reaching down to minute scales at very high energy gamma rays, is today known to be a widespread phenomenon of blazars, calling for very compact emission regions. In most models of such optically thick emission regions, the gamma-ray flux is usually pair-absorbed exponentially, without considering the cascade evolving from the pair-produced electrons. From the observed flux, it is often concluded that emission emanates from larger distances where the region is optically thin, especially from outside of the BLR. Only in few blazars gamma-ray attenuation associated with pair absorption in the BLR was clearly reported.
With the advent of sophisticated high-energy or very high energy gamma-ray detectors, like the Fermi Large Area Telescope or the Major Atmospheric Gamma-ray Imaging Cherenkov telescopes, besides the extraordinarily fast variability spectral features have been found that cannot be explained by conventional models reproducing the two-hump SED. Two such narrow spectral features are discussed in this work. For the nearby blazar Markarian 501, hints to a sharp peak around 3 TeV have been reported from a multi-wavelength campaign carried out in July 2014, while for 3C 279 a spectral dip was found in 2018 data, that can hardly be described with conventional fitting functions. In this work it is examined whether these spectral peculiarities of blazar jet emission can be explained, if the full radiation reprocessing through an IC pair cascade is accounted for.
Such a cascade is the multiple concatenation of IC scattering events and pair production events. In the cascades generally considered in this work, relativistic electrons and high-energy photons are injected into a fixed soft target photon field. A mathematical description for linear IC pair cascades with escape terms is delivered on the basis of preliminary works. The steady-state kinetic equations for the electrons and for the photons are determined, whereby it is paid attention to an explicit formulation and to motivating the correct integration borders of all integrals from kinematic constraints. In determining the potentially observable gamma-ray flux, both the attenuated injected flux and the flux evolving as an effect of IC up-scattering, pair absorption and escape are incorporated, giving the emerging spectra very distinct imprints.
Much effort is dedicated to the numerical solution of the electrons' kinetic equation via iterative schemes. It is explained why pointwise iteration from higher to lower Lorentz factors is more efficient than iterating the whole set of sampling points. The algorithm is parallelised at two positions. First, several workers can perform pointwise iterations simultaneously. Second, the most demanding integral is cut into a number of part integrals which can be determined by multiple workers. Through these measures, the Python code can be readily applied to simulate steady-state IC pair cascades with escape.
In the case of Markarian 501 the developed framework is as follows. The AGN hosts an advection-dominated accretion flow with a normalised accretion rate of several \(10^{-4}\) and an electron temperature near \(10^{10}\) K. On the one hand, the accretion flow illuminates the few ambient gas clouds with approximate radius \(10^{11}\) m, which reprocess a fraction 0.01 of the luminosity into hydrogen and helium emission lines. On the other hand, the gamma rays from the accretion flow create electrons and positrons in a sporadically active vacuum gap in the BH magnetosphere. In the active gap, a power of roughly 0.001 of the Blandford-Znajek power is extracted from the rotating BH through a gap potential drop of several \(10^{18}\) V, generating ultra-relativistic electrons, which subsequently are multiplied by a factor of about \(10^6\) through interaction with the accretion flow photons. This electron beam propagates away from the central engine and encounters the photon field of one passing ionised cloud. The resulting IC pair cascade is simulated and the evolving gamma-ray spectrum is determined. Just above the absorption troughs due to the hydrogen lines, the spectrum exhibits a narrow bump around 3 TeV. When the cascaded emission is added to the emission generated at larger distances, the observed multi-wavelength SED including the sharp peak at 3 TeV is reproduced, underlining that radiation processes beyond conventional models are motivated by distinct spectral features.
The dip in the spectrum of 3C 279 is addressed by a similar cascade model. Three types of injection are considered, varying in the ratio of the photon density to the electron density and varying in the spectral shape. The IC pair cascade is assumed to happen either in the dense BLR photon field with a luminosity of several \(10^{37}\) W and a radial size of few \(10^{14}\) m or in the diluted photon field outside of the BLR. The latter scenario is however rejected as the spectral slope around several 100 MeV and the dip at few 10 GeV cannot be reconciled within this model. The radiation cascaded in the BLR can explain the observational data, irrespective of the assumed injected rate. It is therefore concluded that for this period of gamma-ray emission, the radiation production happens at the edge of the BLR of 3C 279.
Both investigations show that IC pair cascades can account for fine structure seen in blazar SEDs. It is insufficient to restrict the radiation transport to pure exponential absorption of an injection term. Pair production and IC up-scattering by all generations of photons and electrons in the optically thick regime critically shape the emerging spectra. As the advent of future improved detectors will provide more high-precision spectra, further observations of narrow spectral features can be expected. It seems therefore recommendable to incorporate cascading into conventional radiation production models or to extend the model developed in this work by synchrotron radiation.