Refine
Year of publication
Document Type
- Article (15)
- Doctoral Thesis (3)
- Postprint (3)
Language
- English (21)
Is part of the Bibliography
- yes (21)
Keywords
- magnetic fields (21) (remove)
Institute
We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGNs), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma-rays from an AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background photons via inverse-Compton scattering to produce gamma-rays. Due to the deflection of the electron- positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218 vertical bar 304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around 10(-14) G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES. 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.
Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated and, on account of weak energy losses, can accumulate, suggesting that clusters may also be sources of very high energy (VHE; E > 100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cerenkov telescopes, with complementing Fermi Large Area Telescope observations at GeV energies. No significant gamma-ray emission from the Coma Cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5) x 10(-8) photonsm(-2) s(-1) (VERITAS, >220 GeV) and similar to 2 x 10(-6) photonsm(-2) s(-1) (Fermi, 1-3GeV), respectively. We use the gamma-ray upper limits to constrain cosmic rays (CRs) and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and <1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of similar to(2-5.5) mu G, depending on the radial magnetic field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this renders the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <sigma nu >.
During an unusually massive filament eruption on 7 June 2011, SDO/AIA imaged for the first time significant EUV emission around a magnetic reconnection region in the solar corona. The reconnection occurred between magnetic fields of the laterally expanding CME and a neighbouring active region. A pre-existing quasi-separatrix layer was activated in the process. This scenario is supported by data-constrained numerical simulations of the eruption. Observations show that dense cool filament plasma was re-directed and heated in situ, producing coronal-temperature emission around the reconnection region. These results provide the first direct observational evidence, supported by MHD simulations and magnetic modelling, that a large-scale re-configuration of the coronal magnetic field takes place during solar eruptions via the process of magnetic reconnection.
Aims. We study an evolving bipolar active region that exhibits flux cancellation at the internal polarity inversion line, the formation of a soft X-ray sigmoid along the inversion line and a coronal mass ejection. The aim is to investigate the quantity of flux cancellation that is involved in flux rope formation in the time period leading up to the eruption.
Methods. The active region is studied using its extreme ultraviolet and soft X-ray emissions as it evolves from a sheared arcade to flux rope configuration. The evolution of the photospheric magnetic field is described and used to estimate how much flux is reconnected into the flux rope.
Results. About one third of the active region flux cancels at the internal polarity inversion line in the 2.5 days leading up to the eruption. In this period, the coronal structure evolves from a weakly to a highly sheared arcade and then to a sigmoid that crosses the inversion line in the inverse direction. These properties suggest that a flux rope has formed prior to the eruption. The amount of cancellation implies that up to 60% of the active region flux could be in the body of the flux rope. We point out that only part of the cancellation contributes to the flux in the rope if the arcade is only weakly sheared, as in the first part of the evolution. This reduces the estimated flux in the rope to similar to 30% or less of the active region flux. We suggest that the remaining discrepancy between our estimate and the limiting value of similar to 10% of the active region flux, obtained previously by the flux rope insertion method, results from the incomplete coherence of the flux rope, due to nonuniform cancellation along the polarity inversion line. A hot linear feature is observed in the active region which rises as part of the eruption and then likely traces out the field lines close to the axis of the flux rope. The flux cancellation and changing magnetic connections at one end of this feature suggest that the flux rope reaches coherence by reconnection immediately before and early in the impulsive phase of the associated flare. The sigmoid is destroyed in the eruption but reforms quickly, with the amount of cancellation involved being much smaller than in the course of its original formation.
Supernovae are known to be the dominant energy source for driving turbulence in the interstellar medium. Yet, their effect on magnetic field amplification in spiral galaxies is still poorly understood. Analytical models based on the uncorrelated-ensemble approach predicted that any created field will be expelled from the disk before a significant amplification can occur. By means of direct simulations of supernova-driven turbulence, we demonstrate that this is not the case. Accounting for vertical stratification and galactic differential rotation, we find an exponential amplification of the mean field on timescales of 100Myr. The self-consistent numerical verification of such a “fast dynamo” is highly beneficial in explaining the observed strong magnetic fields in young galaxies. We, furthermore, highlight the importance of rotation in the generation of helicity by showing that a similar mechanism based on Cartesian shear does not lead to a sustained amplification of the mean magnetic field. This finding impressively confirms the classical picture of a dynamo based on cyclonic turbulence.
It has been suggested that coronal mass ejections (CMEs) remove the magnetic he-licity of their coronal source region from the Sun. Such removal is often regarded to be necessary due to the hemispheric sign preference of the helicity, which inhibits a simple annihilation by reconnection between volumes of opposite chirality. Here we monitor the relative magnetic he-licity contained in the coronal volume of a simulated flux rope CME, as well as the upward flux of relative helicity through horizontal planes in the simulation box. The unstable and erupting flux rope carries away only a minor part of the initial relative helicity; the major part remains in the volume. This is a consequence of the requirement that the current through an expanding loop must decrease if the magnetic energy of the configuration is to decrease as the loop rises, to provide the kinetic energy of the CME.
We quantitatively address the conjecture that magnetic helicity must be shed from the Sun by eruptions launching coronal mass ejections in order to limit its accumulation in each hemisphere. By varying the ratio of guide and strapping field and the flux rope twist in a parametric simulation study of flux rope ejection from approximately marginally stable force-free equilibria, different ratios of self- and mutual helicity are set and the onset of the torus or helical kink instability is obtained. The helicity shed is found to vary over a broad range from a minor to a major part of the initial helicity, with self helicity being largely or completely shed and mutual helicity, which makes up the larger part of the initial helicity, being shed only partly. Torus-unstable configurations with subcritical twist and without a guide field shed up to about two-thirds of the initial helicity, while a highly twisted, kink-unstable configuration sheds only about one-quarter. The parametric study also yields stable force-free flux rope equilibria up to a total flux-normalized helicity of 0.25, with a ratio of self- to total helicity of 0.32 and a ratio of flux rope to external poloidal flux of 0.94. These results numerically demonstrate the conjecture of helicity shedding by coronal mass ejections and provide a first account of its parametric dependence. Both self- and mutual helicity are shed significantly; this reduces the total initial helicity by a fraction of ∼0.4--0.65 for typical source region parameters.
Stellar magnetic fields, as a crucial component of star formation and evolution, evade direct observation at least with current and near future instruments. However investigating whether magnetic fields are generated by a dynamo process or represent relics from the formation process, or whether they show a behavior similar to the sun or something very different, it is essential to investigate their structure and temporal evolution. Fortunately nature provides us with the possibility to indirectly observe surface topologies on distant stars by means of Doppler shift and polarization of light, though not without its challenges. Based on the mentioned effects, the so called Zeeman-Doppler Imaging technique is a powerful method to retrieve magnetic fields from rapid rotating stars based on measurements of spectropolarimetric observations in terms of Stokes profiles. In recent years, a large number of stellar magnetic field distributions could be reconstructed by Zeeman-Doppler Imaging (ZDI). However, the implementation of this method often relies on many approximations because, as an inversion method, it entails enormous computational requirements. The aim of this thesis is to develop methods for a ZDI, designed to invert time-resolved spectropolarimetric data of active late type stars, and to account for the expected complex and small scale magnetic fields on these stars. In order to reliably reconstruct the detailed field orientation and strength, the inversion method is employed to be able to use of all four Stokes components. Furthermore it is based on fully polarized radiative transfer calculations to account for the intricate interplay between temperature and magnetic field. Finally, the application of a newly developed ZDI code to Stokes I and V observations of II Pegasi (short: II Peg) was supposed to deliver the first magnetic surface maps for this highly active star. To accomplish the high computational burden of a radiative transfer based ZDI, we developed a novel approximation method to speed up the inversion process. It is based on Principal Component Analysis and Artificial Neural Networks. The latter approximate the functional mapping between atmospheric parameters and the corresponding local Stokes profiles. Inverse problems, as we are dealing with, are potentially ill-posed and require a regularization method. We propose a new regularization scheme, which implements a local entropy function that accounts for the peculiarities of the reconstruction of localized magnetic fields. To deal with the relatively large noise that is always present in polarimetric data, we developed a multi-line denoising technique based on Principal Component Analysis. In contrast to other multi-line techniques that extract from a large number of spectral lines a sort of mean profile, this method allows to extract individual spectral lines and thus allows for an inversion on the basis of specific lines. All these methods are incorporated in our newly developed ZDI code iMap, which is based on a conjugated gradient method. An in depth validation of our new synthesis method demonstrates the reliability and accuracy of this approach as well as a gain in computation time by almost three orders of magnitude relative to the conventional radiative transfer calculations. We investigated the influence of the different Stokes components (IV / IVQU) on the ability to reconstruct a known synthetic field configuration. In doing so we validate the capability of our inversion code, and we also assess limitations of magnetic field inversions in general. In a first application to II Peg, a K2 IV subgiant, we derived temperature and magnetic field surface distributions from spectropolarimetric data obtained in 2004 and 2007. It gives for the first time simultaneously the temporal evolution of the surface temperature and magnetic field distribution on II Peg.
The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.
Magnetic agnetic field generation in core-sheath jets via the kinetic Kelvin-Helmholtz instability
(2014)
We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas, we find generation of strong large-scale DC currents and magnetic fields that extend over the entire shear surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas, we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates a transverse structure similar to that produced by the Weibel instability.
The architecture of exoplanetary systems is often different from the solar system, with some exoplanets being in close orbits around their host stars and having orbital periods of only a few days. In analogy to interactions between stars in close binary systems, one may expect interactions between the star and the exoplanet as well. From theoretical considerations, effects on the host star through tidal and magnetic interaction with the exoplanet are possible; for the exoplanet, some interesting implications are the evaporation of the planetary atmosphere and potential effects on the planetary magnetism. In this review, several possible interaction pathways and their observational prospects and existing evidence are discussed. A particular emphasis is put on observational opportunities for these kinds of effects in the high-energy regime.
In order to explain the variance of the solar rotation law during the activity minima and maxima, the angular momentum transport by rotating magnetoconvection is simulated in a convective box penetrated by an inclined azimuthal magnetic field. Turbulence-induced kinetic and magnetic stresses and the Lorentz force of the large-scale magnetic background field are the basic transporters of angular momentum. Without rotation, the sign of the magnetic stresses naturally depends on the signs of the field components as positive (negative) B theta B phi transport the angular momentum poleward (equatorward). For fast enough rotation, however, the turbulence-originated Reynolds stresses start to dominate the transport of the angular momentum flux. The simulations show that positive ratios of the two meridional magnetic field components to the azimuthal field reduce the inward radial as well as the equatorward latitudinal transport, which result from hydrodynamic calculations. Only for B theta B phi>0 (generated by solar-type rotation laws with an accelerated equator) does the magnetic-influenced rotation at the solar surface prove to be flatter than the nonmagnetic profile together with the observed slight spin-down of the equator. The latter phenomenon does not appear for antisolar rotation with polar vortex as well as for rotation laws with prevailing radial shear.
A conducting Taylor-Couette flow with quasi-Keplerian rotation law containing a toroidal magnetic field serves as a mean-field dynamo model of the Tayler-Spruit type. The flows are unstable against non-axisymmetric perturbations which form electromotive forces defining a effect and eddy diffusivity. If both degenerated modes with m = +/- 1 are excited with the same power then the global a effect vanishes and a dynamo cannot work. It is shown, however, that the Tayler instability produces finite alpha effects if only an isolated mode is considered but this intrinsic helicity of the single-mode is too low for an alpha(2) dynamo. Moreover, an alpha Omega dynamo model with quasi-Keplerian rotation requires a minimum magnetic Reynolds number of rotation of Rm similar or equal to 2000 to work. Whether it really works depends on assumptions about the turbulence energy. For a steeper-than-quadratic dependence of the turbulence intensity on the magnetic field, however, dynamos are only excited if the resulting magnetic eddy diffusivity approximates its microscopic value, eta(T) similar or equal to eta. By basically lower or larger eddy diffusivities the dynamo instability is suppressed.
We revisit the effect of nonlinear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in k-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that nonlinear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wave numbers where wave driving is weak or absent. The relaxation time also strongly depends on the intergalactic medium temperature, T-IGM, and for T-IGM << 10 eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small k, that nonlinearly drains the wave energy at the resonant k of the pair-beam instability. Any other dissipation process operating at small k would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other nonlinear processes, such as the modulation instability, could provide additional damping of the nonresonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.
This work revisits the electrostatic instability for blazar-induced pair beams propagating through the intergalactic medium (IGM) using linear analysis and PIC simulations. We study the impact of the realistic distribution function of pairs resulting from the interaction of high-energy gamma-rays with the extragalactic background light. We present analytical and numerical calculations of the linear growth rate of the instability for the arbitrary orientation of wave vectors. Our results explicitly demonstrate that the finite angular spread of the beam dramatically affects the growth rate of the waves, leading to the fastest growth for wave vectors quasi-parallel to the beam direction and a growth rate at oblique directions that is only a factor of 2-4 smaller compared to the maximum. To study the nonlinear beam relaxation, we performed PIC simulations that take into account a realistic wide-energy distribution of beam particles. The parameters of the simulated beam-plasma system provide an adequate physical picture that can be extrapolated to realistic blazar-induced pairs. In our simulations, the beam looses only 1% of its energy, and we analytically estimate that the beam would lose its total energy over about 100 simulation times. An analytical scaling is then used to extrapolate the parameters of realistic blazar-induced pair beams. We find that they can dissipate their energy slightly faster by the electrostatic instability than through inverse-Compton scattering. The uncertainties arising from, e.g., details of the primary gamma-ray spectrum are too large to make firm statements for individual blazars, and an analysis based on their specific properties is required.
The general electromagnetic fluctuation theory is a powerful tool to analyze the magnetic fluctuation spectrum of a plasma. Recent works utilizing this theory for a magnetized non-relativistic isotropic Maxwellian electron-proton plasma have demonstrated that the equilibrium ratio of vertical bar delta B vertical bar/B-0 can be as high as 10(-12). This value results from the balance between spontaneous emission of fluctuations and their damping, and it is considerably smaller than the observed value vertical bar delta B vertical bar/B-0 in the solar wind at 1 au, where 10(-3) less than or similar to vertical bar delta B vertical bar/B-0 less than or similar to 10(-1). In the present manuscript, we consider an anisotropic bi-Maxwellian distribution function to investigate the effect of plasma instabilities on the magnetic field fluctuations. We demonstrate that these instabilities strongly amplify the magnetic field fluctuations and provide a sufficient mechanism to explain the observed value of vertical bar delta B vertical bar/B-0 in the solar wind at 1 au.
Context. For magnetically driven events, the magnetic energy of the system is the prime energy reservoir that fuels the dynamical evolution. In the solar context, the free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. A trustworthy estimation of the magnetic energy is therefore needed in three-dimensional (3D) models of the solar atmosphere, e. g., in coronal fields reconstructions or numerical simulations.
Aims. The expression of the energy of a system as the sum of its potential energy and its free energy (Thomson's theorem) is strictly valid when the magnetic field is exactly solenoidal. For numerical realizations on a discrete grid, this property may be only approximately fulfilled. We show that the imperfect solenoidality induces terms in the energy that can lead to misinterpreting the amount of free energy present in a magnetic configuration.
Methods. We consider a decomposition of the energy in solenoidal and nonsolenoidal parts which allows the unambiguous estimation of the nonsolenoidal contribution to the energy. We apply this decomposition to six typical cases broadly used in solar physics. We quantify to what extent the Thomson theorem is not satisfied when approximately solenoidal fields are used.
Results. The quantified errors on energy vary from negligible to significant errors, depending on the extent of the nonsolenoidal component of the field. We identify the main source of errors and analyze the implications of adding a variable amount of divergence to various solenoidal fields. Finally, we present pathological unphysical situations where the estimated free energy would appear to be negative, as found in some previous works, and we identify the source of this error to be the presence of a finite divergence.
Conclusions. We provide a method of quantifying the effect of a finite divergence in numerical fields, together with detailed diagnostics of its sources. We also compare the efficiency of two divergence-cleaning techniques. These results are applicable to a broad range of numerical realizations of magnetic fields.
Microfabricated solid-state surfaces, also called atom chip', have become a well-established technique to trap and manipulate atoms. This has simplified applications in atom interferometry, quantum information processing, and studies of many-body systems. Magnetic trapping potentials with arbitrary geommetries are generated with atom chip by miniaturized current-carrying conductors integrated on a solid substrate. Atoms can be trapped and cooled to microKelvin and even nanoKelvin temperatures in such microchip trap. However, cold atoms can be significantly perturbed by the chip surface, typically held at room temperature. The magnetic field fluctuations generated by thermal currents in the chip elements may induce spin flips of atoms and result in loss, heating and decoherence. In this thesis, we extend previous work on spin flip rates induced by magnetic noise and consider the more complex geometries that are typically encountered in atom chips: layered structures and metallic wires of finite cross-section. We also discuss a few aspects of atom chips traps built with superconducting structures that have been suggested as a means to suppress magnetic field fluctuations. The thesis describes calculations of spin flip rates based on magnetic Green functions that are computed analytically and numerically. For a chip with a top metallic layer, the magnetic noise depends essentially on the thickness of that layer, as long as the layers below have a much smaller conductivity. Based on this result, scaling laws for loss rates above a thin metallic layer are derived. A good agreement with experiments is obtained in the regime where the atom-surface distance is comparable to the skin depth of metal. Since in the experiments, metallic layers are always etched to separate wires carrying different currents, the impact of the finite lateral wire size on the magnetic noise has been taken into account. The local spectrum of the magnetic field near a metallic microstructure has been investigated numerically with the help of boundary integral equations. The magnetic noise significantly depends on polarizations above flat wires with finite lateral width, in stark contrast to an infinitely wide wire. Correlations between multiple wires are also taken into account. In the last part, superconducting atom chips are considered. Magnetic traps generated by superconducting wires in the Meissner state and the mixed state are studied analytically by a conformal mapping method and also numerically. The properties of the traps created by superconducting wires are investigated and compared to normal conducting wires: they behave qualitatively quite similar and open a route to further trap miniaturization, due to the advantage of low magnetic noise. We discuss critical currents and fields for several geometries.
Spectroscopic observations play essential roles in astrophysics. They are crucial for determining physical parameters in our Universe, providing information about the chemistry of various astronomical environments. The proper execution of the spectroscopic analysis requires accounting for all the physical effects that are compatible to the signal-to-noise ratio. We find in this paper the influence on spectroscopy from the atomic/ground state alignment owing to anisotropic radiation and modulated by interstellar magnetic field, has significant impact on the study of interstellar gas. In different observational scenarios, we comprehensively demonstrate how atomic alignment influences the spectral analysis and provide the expressions for correcting the effect. The variations are even more pronounced for multiplets and line ratios. We show the variation of the deduced physical parameters caused by the atomic alignment effect, including alpha-to-iron ratio ([X/Fe]) and ionization fraction. Synthetic observations are performed to illustrate the visibility of such effect with current facilities. A study of Photodissociation regions in rho Ophiuchi cloud is presented to demonstrate how to account for atomic alignment in practice. Our work has shown that due to its potential impact, atomic alignment has to be included in an accurate spectroscopic analysis of the interstellar gas with current observational capability.