Source author record

M. Krause

M. Krause appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

Mirror production for the Cherenkov telescopes of the ASTRI Mini-Array and of the MST project for the Cherenkov Telescope Array

The Cherenkov Telescope Array (CTA) is the next ground-based $γ$-ray observatory in the TeV $γ$-ray spectral region operating with the Imaging Atmospheric Cherenkov Technique. It is based on almost 70 telescopes of different class diameters - LST, MST and SST of 23, 12, and 4 m, respectively - to be installed in two sites in the two hemispheres (at La Palma, Canary Islands, and near Paranal, Chile). Several thousands of reflecting mirror tiles larger than 1 m$^2$ will be produced for realizing the segmented primary mirrors of a so large number of telescopes. Almost in parallel, the ASTRI Mini-Array (MA) is being implemented in Tenerife (Canary Islands), composed of nine 4 m diameter dual-mirror Cherenkov telescopes (very similar to the SSTs). We completed the mirror production for all nine telescopes of the ASTRI MA and two MST telescopes (400 segments in total) using the cold glass slumping replication technology. The results related to the quality achieved with a so large-scale production are presented, also discussing the adopted testing methods and approaches. They will be very useful for the adoption and optimization of the quality assurance process for the huge production (almost 3000 m$^2$ of reflecting surface) of the MST and SST CTA telescopes.

preprint2020arXiv

A decade of multi-wavelength observations of the TeV blazar 1ES 1215+303: Extreme shift of the synchrotron peak frequency and long-term optical-gamma-ray flux increase

Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multi-wavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z=0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by VLBA at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.

preprint2020arXiv

Evidence for proton acceleration up to TeV energies based on VERITAS and Fermi-LAT observations of the Cas A SNR

We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $γ$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $γ$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{μG}$ is deduced.

preprint2020arXiv

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

preprint2020arXiv

VERITAS Detection of LS 5039 and HESS J1825-137

With 8 hours of observations, VERITAS confirms the detection of two very high energy gamma-ray sources. The gamma-ray binary LS 5039 is detected with a statistical significance of $8.8σ$. The measured flux above 1 TeV is $(2.5 \pm 0.4) \times 10^{-12} \rm \, cm^{-2} \, s^{-1}$ near inferior conjunction and $(7.8 \pm 2.8) \times 10^{-13} \rm \, cm^{-2} \, s^{-1}$ near superior conjunction. The pulsar wind nebula HESS J1825-137 is detected with a statistical significance of $6.7σ$ and a measured flux above 1 TeV of $(3.9 \pm 0.8) \times 10^{-12} \rm \, cm^{-2} \, s^{-1}$.

preprint2020arXiv

VVV-WIT-01: highly obscured classical nova or protostellar collision?

A search of the first Data Release of the VISTA Variables in the Via Lactea (VVV) Survey discovered the exceptionally red transient VVV-WIT-01 ($H-K_s=5.2$). It peaked before March 2010, then faded by $\sim$9.5 mag over the following two years. The 1.6--22 $μ$m spectral energy distribution in March 2010 was well fit by a highly obscured black body with $T \sim 1000$ K and $A_{K_s} \sim 6.6$ mag. The source is projected against the Infrared Dark Cloud (IRDC) SDC G331.062$-$0.294. The chance projection probability is small for any single event ($p \approx 0.01$ to 0.02) which suggests a physical association, e.g. a collision between low mass protostars. However, black body emission at $T \sim 1000$ K is common in classical novae (especially CO novae) at the infrared peak in the light curve, due to condensation of dust $\sim$30--60 days after the explosion. Radio follow up with the Australia Telescope Compact Array (ATCA) detected a fading continuum source with properties consistent with a classical nova but probably inconsistent with colliding protostars. Considering all VVV transients that could have been projected against a catalogued IRDC raises the probability of a chance association to $p=0.13$ to 0.24. After weighing several options, it appears likely that VVV-WIT-01 was a classical nova event located behind an IRDC.

preprint2019arXiv

ewN2HDECAY - A program for the Calculation of Electroweak One-Loop Corrections to Higgs Decays in the Next-to-Minimal Two-Higgs-Doublet Model Including State-of-the-Art QCD Corrections

We present in this paper our new program package ewN2HDECAY for the calculation of the partial decay widths and branching ratios of the Higgs bosons of the Next-to-Minimal 2-Higgs Doublet Model (N2HDM). The N2HDM is based on a general CP-conserving 2HDM which is extended by a real scalar singlet field. The program computes the complete electroweak one-loop corrections to all non-loop-induced two-body on-shell Higgs boson decays in the N2HDM and combines them with the state-of-the-art QCD corrections that are already implemented in the existing program N2HDECAY. Most of the independent input parameters of the electroweak sector of the N2HDM are renormalized in an on-shell scheme. The soft-$\mathbb{Z}_2$-breaking squared mass scale $m_{12}^2$ and the vacuum expectation value $v_S$ of the $SU(2)_L$ singlet field, however, are renormalized with $\overline{\text{MS}}$ conditions, while for the four scalar mixing angles $α_i$ ($i=1,2,3$) and $β$ of the N2HDM, several different renormalization schemes are applied. By giving out the leading-order and the loop-corrected partial decay widths separately from the branching ratios, the program ewN2HDECAY not only allows for phenomenological analyses of the N2HDM at highest precision, it can also be used for a study of the impact of the electroweak corrections and the remaining theoretical uncertainty due to missing higher-order corrections based on a change of the renormalization scheme. The input parameters are then consistently calculated with a parameter conversion routine when switching from one renormalization scheme to the other. The latest version of the program ewN2HDECAY can be downloaded from the URL \href{https://github.com/marcel-krause/ewN2HDECAY}{https://github.com/marcel-krause/ewN2HDECAY}.

preprint2019arXiv

Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $γ$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

preprint2019arXiv

Gauge Dependences of Higher-Order Corrections to NMSSM Higgs Boson Masses and the Charged Higgs Decay $H^\pm \to W^\pm h_i$

In this paper we compute the electroweak corrections to the charged Higgs boson decay into a $W$ boson and a neutral Higgs boson in the CP-conserving NMSSM. We calculate the process in a general $R_ξ$ gauge and investigate the dependence of the loop-corrected decay width on the gauge parameter $ξ$. The gauge dependence arises from the mixing of different loop orders. Phenomenology requires the inclusion of mass and mixing corrections to the external Higgs bosons in order to match the experimentally measured mass values. As a result, we move away from a strict one-loop calculation and consequently mix orders in perturbation theory. Moreover, determination of the loop-corrected masses in an iterative procedure also results in the mixing of different loop orders. Gauge dependence then arises from the mismatch with tree-level Goldstone boson couplings that are applied in the loop calculation, and from the gauge dependence of the loop-corrected masses themselves. We find that the gauge dependence is significant.

preprint2019arXiv

Impact of Electroweak Corrections on Neutral Higgs Boson Decays in Extended Higgs Sectors

Precision predictions play an important role in the search for indirect New Physics effects in the Higgs sector itself. For the electroweak (EW) corrections of the Higgs bosons in extended Higgs sectors several renormalization schemes have been worked out that provide gauge-parameter-independent relations between the input parameters and the computed observables. Our recently published program codes 2HDECAY and ewN2HDECAY allow for the computation of the EW corrections to the Higgs decay widths and branching ratios of the Two-Higgs-Doublet Model (2HDM) and the Next-to-Minimal-2HDM (N2HDM) for different renormalization schemes of the scalar mixing angles. In this paper, we present a comprehensive and complete overview over the relative size of the EW corrections to the branching ratios of the 2HDM and N2HDM neutral Higgs bosons for different applied renormalization schemes. We quantify the size of the EW corrections of Standard Model(SM)- and non-SM-like Higgs bosons and moreover also identify renormalization schemes that are well-behaved and do not induce unnaturally large corrections. We furthermore pin down decays and parameter regions that feature large EW corrections and need further treatment in order to improve the predictions. Our study sets the scene for future work in the computation of higher-order corrections to the decays of non-minimal Higgs sectors.

preprint2019arXiv

Measurement of the extragalactic background light spectral energy distribution with VERITAS

The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the Universe's history. Spectral measurements of blazars at very high energies ($>$100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the gamma-ray photons with the EBL. The spectra of 14 VERITAS-detected blazars are included in a new measurement of the EBL SED that is independent of EBL SED models. The resulting SED covers an EBL wavelength range of 0.56--56 $μ$m, and is in good agreement with lower limits obtained by assuming that the EBL is entirely due to radiation from cataloged galaxies.

preprint2019arXiv

Probing the Properties of the Pulsar Wind in the Gamma-Ray Binary HESS J0632+057 with NuSTAR and VERITAS Observations

HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about $315$ days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in November and December 2017. These observations correspond to the orbital phases $ϕ\approx0.22$ and $0.3$, where the fluxes are rising towards the first light-curve peak. A significant variation of the spectral index from 1.77$\pm$0.05 to 1.56$\pm$0.05 is observed in the X-ray data. The multi-wavelength spectral energy distributions (SED) derived from the observations are interpreted in terms of a leptonic model, in which the compact object is assumed to be a pulsar and non-thermal radiation is emitted by high-energy electrons accelerated at the shock formed by the collision between the stellar and pulsar wind. The results of the SED fitting show that our data can be consistently described within this scenario, and allow us to estimate the magnetization of the pulsar wind at the location of the shock formation. The constraints on the pulsar-wind magnetization provided by our results are shown to be consistent with those obtained from other systems.

preprint2018arXiv

The extreme HBL behaviour of Markarian 501 during 2012

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

preprint2016arXiv

A Search for Brief Optical Flashes Associated with the SETI Target KIC 8462852

The F-type star KIC 8462852 has recently been identified as an exceptional target for SETI (search for extraterrestrial intelligence) observations. We describe an analysis methodology for optical SETI, which we have used to analyse nine hours of serendipitous archival observations of KIC 8462852 made with the VERITAS gamma-ray observatory between 2009 and 2015. No evidence of pulsed optical beacons, above a pulse intensity at the Earth of approximately 1 photon per m^2, is found. We also discuss the potential use of imaging atmospheric Cherenkov telescope arrays in searching for extremely short duration optical transients in general.

preprint2016arXiv

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.

preprint2016arXiv

A Search for Very High-Energy Gamma Rays from the Missing Link Binary Pulsar J1023+0038 with VERITAS

The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259--63/LS 2883, making it an ideal candidate for the study of high-energy non-thermal emission. It has been the subject of multi-wavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high-energy gamma-ray observations carried out by VERITAS before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that very high-energy gamma rays are produced via an inverse-Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than $\sim$2 G before the disappearance of the radio pulsar and greater than $\sim$10 G afterwards.

preprint2016arXiv

Commissioning of the vacuum system of the KATRIN Main Spectrometer

The KATRIN experiment will probe the neutrino mass by measuring the beta-electron energy spectrum near the endpoint of tritium beta-decay. An integral energy analysis will be performed by an electro-static spectrometer (Main Spectrometer), an ultra-high vacuum vessel with a length of 23.2 m, a volume of 1240 m^3, and a complex inner electrode system with about 120000 individual parts. The strong magnetic field that guides the beta-electrons is provided by super-conducting solenoids at both ends of the spectrometer. Its influence on turbo-molecular pumps and vacuum gauges had to be considered. A system consisting of 6 turbo-molecular pumps and 3 km of non-evaporable getter strips has been deployed and was tested during the commissioning of the spectrometer. In this paper the configuration, the commissioning with bake-out at 300°C, and the performance of this system are presented in detail. The vacuum system has to maintain a pressure in the 10^{-11} mbar range. It is demonstrated that the performance of the system is already close to these stringent functional requirements for the KATRIN experiment, which will start at the end of 2016.

preprint2016arXiv

Exceptionally bright TeV flares from the binary LS I +61$^\circ$ 303

The TeV binary system LS I +61$^\circ$ 303 is known for its regular, non-thermal emission pattern which traces the orbital period of the compact object in its 26.5 day orbit around its B0 Ve star companion. The system typically presents elevated TeV emission around apastron passage with flux levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV). In this article, VERITAS observations of LS I +61$^\circ$ 303 taken in late 2014 are presented, during which bright TeV flares around apastron at flux levels peaking above 30% of the Crab Nebula flux were detected. This is the brightest such activity from this source ever seen in the TeV regime. The strong outbursts have rise and fall times of less than a day. The short timescale of the flares, in conjunction with the observation of 10 TeV photons from LS I +61$^\circ$ 303 during the flares, provides constraints on the properties of the accelerator in the source.

preprint2016arXiv

TeV Gamma-ray Observations of The Galactic Center Ridge By VERITAS

The Galactic Center Ridge has been observed extensively in the past by both GeV and TeV gamma-ray instruments revealing a wealth of structure, including a diffuse component as well as the point sources G0.9+0.1 (a composite supernova remnant) and Sgr A* (believed to be associated with the supermassive black hole located at the center of our Galaxy). Previous very high energy (VHE) gamma-ray observations with the H.E.S.S. experiment have also detected an extended TeV gamma-ray component along the Galactic plane in the >300 GeV gamma-ray regime. Here we report on observations of the Galactic Center Ridge from 2010-2014 by the VERITAS telescope array in the >2 TeV energy range. From these observations we 1.) provide improved measurements of the differential energy spectrum for Sgr A* in the >2 TeV gamma-ray regime, 2.) provide a detection in the >2 TeV gamma-ray emission from the composite SNR G0.9+0.1 and an improved determination of its multi-TeV gamma-ray energy spectrum, 3.) report on the detection of VER J1746-289, a localized enhancement of >2 TeV gamma-ray emission along the Galactic plane.

preprint2016arXiv

Timing calibration and spectral cleaning of LOFAR time series data

We describe a method for spectral cleaning and timing calibration of short voltage time series data from individual radio interferometer receivers. It makes use of the phase differences in Fast Fourier Transform (FFT) spectra across antenna pairs. For strong, localized terrestrial sources these are stable over time, while being approximately uniform-random for a sum over many sources or for noise. Using only milliseconds-long datasets, the method finds the strongest interfering transmitters, a first-order solution for relative timing calibrations, and faulty data channels. No knowledge of gain response or quiescent noise levels of the receivers is required. With relatively small data volumes, this approach is suitable for use in an online system monitoring setup for interferometric arrays. We have applied the method to our cosmic-ray data collection, a collection of measurements of short pulses from extensive air showers, recorded by the LOFAR radio telescope. Per air shower, we have collected 2 ms of raw time series data for each receiver. The spectral cleaning has a calculated optimal sensitivity corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well sufficient for our application. Timing calibration across individual antenna pairs has been performed at 0.4 ns precision; for calibration of signal clocks across stations of 48 antennas the precision is 0.1 ns. Monitoring differences in timing calibration per antenna pair over the course of the period 2011 to 2015 shows a precision of 0.08 ns, which is useful for monitoring and correcting drifts in signal path synchronizations. A cross-check method for timing calibration is presented, using a pulse transmitter carried by a drone flying over the array. Timing precision is similar, 0.3 ns.

preprint2016arXiv

Upper limits from five years of blazar observations with the VERITAS Cherenkov telescopes

Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E>100 GeV) γ-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of ~570 hours. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog which are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have a spectroscopic distance estimate. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data-set, which shows a 4 σ excess.

preprint2016arXiv

VERITAS and Multiwavelength Observations of the BL Lacertae Object 1ES 1741+196

We present results from multiwavelength observations of the BL Lacertae object 1ES 1741+196, including results in the very-high-energy $γ$-ray regime using the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The VERITAS time-averaged spectrum, measured above 180 GeV, is well-modelled by a power law with a spectral index of $2.7\pm0.7_{\mathrm{stat}}\pm0.2_{\mathrm{syst}}$. The integral flux above 180 GeV is $(3.9\pm0.8_{\mathrm{stat}}\pm1.0_{\mathrm{syst}})\times 10^{-8}$ m$^{-2}$ s$^{-1}$, corresponding to 1.6% of the Crab Nebula flux on average. The multiwavelength spectral energy distribution of the source suggests that 1ES 1741+196 is an extreme-high-frequency-peaked BL Lacertae object. The observations analysed in this paper extend over a period of six years, during which time no strong flares were observed in any band. This analysis is therefore one of the few characterizations of a blazar in a non-flaring state.

preprint2016arXiv

Very High Energy outburst of Markarian 501 in May 2009

The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $γ$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.

preprint2016arXiv

Very-high-energy observations of the binaries V 404 Cyg and 4U 0115+634 during giant X-ray outbursts

Transient X-ray binaries produce major outbursts in which the X-ray flux can increase over the quiescent level by factors as large as $10^7$. The low-mass X-ray binary V 404 Cyg and the high-mass system 4U 0115+634 underwent such major outbursts in June and October 2015, respectively. We present here observations at energies above hundreds of GeV with the VERITAS observatory taken during some of the brightest X-ray activity ever observed from these systems. No gamma-ray emission has been detected by VERITAS in 2.5 hours of observations of the microquasar V 404 Cyg from 2015, June 20-21. The upper flux limits derived from these observations on the gamma-ray flux above 200 GeV of F $< 4.4\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ correspond to a tiny fraction (about $10^{-6}$) of the Eddington luminosity of the system, in stark contrast to that seen in the X-ray band. No gamma rays have been detected during observations of 4U 0115+634 in the period of major X-ray activity in October 2015. The flux upper limit derived from our observations is F $< 2.1\times 10^{-12}$ cm$^{-2}$ s$^{-1}$ for gamma rays above 300 GeV, setting an upper limit on the ratio of gamma-ray to X-ray luminosity of less than 4%.

preprint2015arXiv

Calibrating the absolute amplitude scale for air showers measured at LOFAR

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

preprint2015arXiv

First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsähovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

preprint2015arXiv

Science Highlights from VERITAS

The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is a ground-based array located at the Fred Lawrence Whipple Observatory in southern Arizona and is one of the world's most sensitive gamma-ray instruments at energies of 85 GeV to $>$30 TeV. VERITAS has a wide scientific reach that includes the study of extragalactic and Galactic objects as well as the search for astrophysical signatures of dark matter and the measurement of cosmic rays. In this paper, we will summarize the current status of the VERITAS observatory and present some of the scientific highlights from the last two years, focusing in particular on those results shown at the 2015 ICRC in The Hague, Netherlands.

preprint2015arXiv

VERITAS detection of $γ$-ray flaring activity from the BL Lac object 1ES 1727+502 during bright moonlight observations

During moonlit nights, observations with ground-based Cherenkov telescopes at very high energies (VHE, $E>100$ GeV) are constrained since the photomultiplier tubes (PMTs) in the telescope camera are extremely sensitive to the background moonlight. Observations with the VERITAS telescopes in the standard configuration are performed only with a moon illumination less than 35$\%$ of full moon. Since 2012, the VERITAS collaboration has implemented a new observing mode under bright moonlight, by either reducing the voltage applied to the PMTs (reduced-high-voltage configuration, RHV), or by utilizing UV-transparent filters. While these operating modes result in lower sensitivity and increased energy thresholds, the extension of the available observing time is useful for monitoring variable sources such as blazars and sources requiring spectral measurements at the highest energies. In this paper we report the detection of $γ$-ray flaring activity from the BL Lac object 1ES 1727+502 during RHV observations. This detection represents the first evidence of VHE variability from this blazar. The integral flux is $(1.1\pm0.2)\times10^{-11}\mathrm{cm^{-2}s^{-1}}$ above 250 GeV, which is about five times higher than the low-flux state. The detection triggered additional \veritas\ observations during standard dark-time. Multiwavelength observations with the FLWO 48" telescope, and the Swift and Fermi satellites are presented and used to produce the first spectral energy distribution (SED) of this object during $γ$-ray flaring activity. The SED is then fitted with a standard synchrotron-self-Compton model, placing constraints on the properties of the emitting region and of the acceleration mechanism at the origin of the relativistic particle population in the jet.

preprint2014arXiv

LORA: A scintillator array for LOFAR to measure extensive air showers

The measurement of the radio emission from extensive air showers, induced by high-energy cosmic rays is one of the key science projects of the LOFAR radio telescope. The LOfar Radboud air shower Array (LORA) has been installed in the core of LOFAR in the Netherlands. The main purpose of LORA is to measure the properties of air showers and to trigger the read-out of the LOFAR radio antennas to register extensive air showers. The experimental set-up of the array of scintillation detectors and its performance are described.

preprint2014arXiv

Polarized radio emission from extensive air showers measured with LOFAR

We present LOFAR measurements of radio emission from extensive air showers. We find that this emission is strongly polarized, with a median degree of polarization of nearly $99\%$, and that the angle between the polarization direction of the electric field and the Lorentz force acting on the particles, depends on the observer location in the shower plane. This can be understood as a superposition of the radially polarized charge-excess emission mechanism, first proposed by Askaryan and the geomagnetic emission mechanism proposed by Kahn and Lerche. We calculate the relative strengths of both contributions, as quantified by the charge-excess fraction, for $163$ individual air showers. We find that the measured charge-excess fraction is higher for air showers arriving from closer to the zenith. Furthermore, the measured charge-excess fraction also increases with increasing observer distance from the air shower symmetry axis. The measured values range from $(3.3\pm 1.0)\%$ for very inclined air showers at $25\, \mathrm{m}$ to $(20.3\pm 1.3)\%$ for almost vertical showers at $225\, \mathrm{m}$. Both dependencies are in qualitative agreement with theoretical predictions.

preprint2014arXiv

The shape of the radio wavefront of extensive air showers as measured with LOFAR

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

preprint2013arXiv

A Detailed Study of the Radio--FIR Correlation in NGC6946 with Herschel-PACS/SPIRE from KINGFISH

We derive the distribution of the synchrotron spectral index across NGC6946 and investigate the correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the KINGFISH Herschel PACS and SPIRE data. The radio--FIR correlation is studied as a function of star formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. The synchrotron emission follows both star-forming regions and the so-called magnetic arms present in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms ($α_n \sim 1$), while it is flat in places of giant H{\sc ii} regions and in the center of the galaxy ($α_n \sim 0.6-0.7$). The map of $α_n$ provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in the disk of the galaxy. Variations in the synchrotron--FIR correlation across the galaxy are shown to be a function of both star formation and magnetic fields. We find that the synchrotron emission correlates better with cold rather than with warm dust emission, when the interstellar radiation field is the main heating source of dust. The synchrotron--FIR correlation suggests a coupling between the magnetic field and the gas density. NGC6946 shows a power-law behavior between the total (turbulent) magnetic field strength B and the star formation rate surface density $Σ_{\rm SFR}$ with an index of 0.14\,(0.16)$\pm$0.01. This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence expected in actively star forming regions. The scale-by-scale analysis of the synchrotron--FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic ray electrons, resulting in a diffusion coefficient of $D_0=4.6\times 10^{28}$\,cm$^2$\,s$^{-1}$ for 2.2\,GeV CREs.

preprint2013arXiv

Detecting cosmic rays with the LOFAR radio telescope

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

preprint2013arXiv

Pathway to the Square Kilometre Array - The German White Paper -

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

preprint2013arXiv

The relation between magnetic and material arms in models for spiral galaxies

Context. Observations of polarized radio emission show that large-scale (regular) magnetic fields in spiral galaxies are not axisymmetric, but generally stronger in interarm regions. In some nearby galaxies such as NGC 6946 they are organized in narrow magnetic arms situated between the material spiral arms. Aims. The phenomenon of magnetic arms and their relation to the optical spiral arms (the material arms) call for an explanation in the framework of galactic dynamo theory. Several possibilities have been suggested but are not completely satisfactory; here we attempt a consistent investigation. Methods. We use a 2D mean-field dynamo model in the no-z approximation and add injections of small-scale magnetic field, taken to result from supernova explosions, to represent the effects of dynamo action on smaller scales. This injection of small scale field is situated along the spiral arms, where star-formation mostly occurs. Results. A straightforward explanation of magnetic arms as a result of modulation of the dynamo mechanism by material arms struggles to produce pronounced magnetic arms, at least with realistic parameters, without introducing new effects such as a time lag between Coriolis force and α-effect. In contrast, by taking into account explicitly the small-scale magnetic field that is injected into the arms by the action of the star forming regions that are concentrated there, we can obtain dynamo models with magnetic structures of various forms that can be compared with magnetic arms. (abbrev). Conclusions. We conclude that magnetic arms can be considered as coherent magnetic structures generated by large-scale dynamo action, and associated with spatially modulated small-scale magnetic fluctuations, caused by enhanced star formation rates within the material arms.

preprint2013arXiv

The role of massive stars in the turbulent infancy of Galactic globular clusters: Feedback on the intracluster medium, and detailed timeline

A major paradigm shift has recently revolutionized our picture of globular clusters (GC) that were long thought to be simple systems of coeval stars born out of homogeneous material. Indeed, detailed abundance studies of GC long-lived low-mass stars performed with 8-10m class telescopes, together with high-precision photometry of Galactic GCs obtained with HST,have brought compelling clues on the presence of multiple stellar populations in individual GCs. These stellar subgroups can be recognized thanks to their different chemical properties (more precisely by abundance differences in light elements from carbon to aluminium; see Bragaglia, this volume) and by the appearance of multimodal sequences in the colour-magnitude diagrams (see Piotto, this volume). This has a severe impact on our understanding of the early evolution of GCs, and in particular of the possible role that massive stars played in shaping the intra-cluster medium (ICM) and in inducing secondary star formation. Here we summarize the detailed timeline we have recently proposed for the first 40 Myrs in the lifetime of a typical GC following the general ideas of our so-called "Fast Rotating Massive stars scenario" (FRMS, Decressin et al. 2007b) and taking into account the dynamics of interstellar bubbles produced by stellar winds and supernovae. More details can be found in Krause et al. (2012, 2013).

preprint2011arXiv

Modeling the total and polarized emission in evolving galaxies: "spotty" magnetic structures

Future radio observations with the SKA and its precursors will be sensitive to trace spiral galaxies and their magnetic field configurations up to redshift $z\approx3$. We suggest an evolutionary model for the magnetic configuration in star-forming disk galaxies and simulate the magnetic field distribution, the total and polarized synchrotron emission, and the Faraday rotation measures for disk galaxies at $z\la 3$. Since details of dynamo action in young galaxies are quite uncertain, we model the dynamo action heuristically relying only on well-established ideas of the form and evolution of magnetic fields produced by the mean-field dynamo in a thin disk. We assume a small-scale seed field which is then amplified by the small-scale turbulent dynamo up to energy equipartition with kinetic energy of turbulence. The large-scale galactic dynamo starts from seed fields of 100 pc and an averaged regular field strength of 0.02\,$μ$G, which then evolves to a "spotty" magnetic field configuration in about 0.8\,Gyr with scales of about one kpc and an averaged regular field strength of 0.6\,$μ$G. The evolution of these magnetic spots is simulated under the influence of star formation, dynamo action, stretching by differential rotation of the disk, and turbulent diffusion. The evolution of the regular magnetic field in a disk of a spiral galaxy, as well as the expected total intensity, linear polarization and Faraday rotation are simulated in the rest frame of a galaxy at 5\,GHz and 150\,MHz and in the rest frame of the observer at 150\,MHz. We present the corresponding maps for several epochs after disk formation. (abridged)

preprint2011arXiv

The sub-mm morphology of the interacting galaxy NGC 3627

We present sub-mm continuum and heterodyne data of the interacting galaxy NGC 3627, obtained with the HHT and APEX. We find significant changes in the molecular line ratios over small scales in the southeastern part of this galaxy. The kinematics of the CO(2-1) line, as well as the morphology in the 870um continuum emission, suggest a continuation of the western spiral arm on the east side of the galaxy's central area. This continued spiral arm crosses the normal eastern spiral arm in an area which shows an unusual magnetic field configuration. This spiral arm crossing, independently of if it's physical or just projected, may help in the understanding of the observed magnetic field configuration.

preprint2010arXiv

Exploring the magnetic fields in local and distant galaxies

(abridged) We investigate the possibility to recognize the magnetic field structures in nearby galaxies and to test the cosmological evolution of their large- and small-scale magnetic fields with the SKA and its precursors. We estimate the required density of the background polarized sources detected with the SKA for reliable reconstruction and reconstruction of magnetic field structures in nearby spiral galaxies. The dynamo theory is applied to distant galaxies to explore the evolution of magnetic fields in distant galaxies in the context of a hierarchical dark matter cosmology. Under favorite conditions, a \emph{recognition} of large-scale magnetic structures in local star-forming disk galaxies (at a distance $\la 100$ Mpc) is possible from $\ga 10$ RMs towards background polarized sources. Galaxies with strong turbulence or small inclination need more polarized sources for a statistically reliable recognition. A reliable \emph{reconstruction} of the field structure without precognition needs at least 20 RM values on a cut along the projected minor axis which translates to $\approx1200$ sources towards the galaxy. We demonstrate that early regular fields are already in place at $z \sim 4$ (approximately 1.5 Gyr after the disk formation) in massive gas-rich galaxies ($>10^9$ M$_{\sun}$) which then evolve to Milky-Way type galaxies. Major and minor mergers influence the star formation rate and geometry of the disk which has an effect of shifting the generation of regular fields in disks to later epochs. Predictions of the evolutionary model of regular fields, simulations of the evolution of turbulent and large-scale regular fields, total and polarized radio emission of disk galaxies, as well as future observational tests with the SKA are discussed.

preprint2010arXiv

Gas dynamics of the central few parsec region of NGC 1068 fuelled by the evolving nuclear star cluster

High resolution observations with the NIR adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of three-dimensional high resolution hydrodynamical simulations with the Pluto code, we follow the evolution of such clusters, focusing on stellar mass loss. This leads to clumpy or filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. In order to capture the relevant physics in the inner region, we treat this disc separately by viscously evolving the radial surface density distribution. This enables us to link the tens of parsec scale region (accessible via SINFONI observations) to the (sub-)parsec scale region (observable with the MIDI instrument and via water maser emission). In this work, we concentrate on the effects of a parametrised turbulent viscosity to generate angular momentum and mass transfer in the disc and additionally take star formation into account. Input parameters are constrained by observations of the nearby Seyfert 2 galaxy NGC 1068. At the current age of its nuclear starburst of 250 Myr, our simulations yield disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 solar masses and mass transfer rates of 0.025 solar masses per year through the inner rim of the disc. This shows that our large scale torus model is able to approximately account for the disc size as inferred from interferometric observations in the mid-infrared and compares well to the extent and mass of a rotating disc structure as inferred from water maser observations. Several other observational constraints are discussed as well.

preprint2009arXiv

[The Impact of Nuclear Star Formation on Gas Inflow to AGN

Our adaptive optics observations of nearby AGN at spatial resolutions as small as 0.085arcsec show strong evidence for recent, but no longer active, nuclear star formation. We begin by describing observations that highlight two contrasting methods by which gas can flow into the central tens of parsecs. Gas accumulation in this region will inevitably lead to a starburst, and we discuss the evidence for such events. We then turn to the impact of stellar evolution on the further inflow of gas by combining a phenomenological approach with analytical modelling and hydrodynamic simulations. These complementary perspectives paint a picture in which all the processes are ultimately regulated by the mass accretion rate into the central hundred parsecs, and the ensuing starburst that occurs there. The resulting supernovae delay accretion by generating a starburst wind, which leaves behind a clumpy interstellar medium. This provides an ideal environment for slower stellar outflows to accrete inwards and form a dense turbulent disk on scales of a few parsecs. Such a scenario may resolve the discrepancy between the larger scale structure seen with adaptive optics and the small scale structure seen with VLTI.

preprint2009arXiv

Cosmic rays and the magnetic field in the nearby starburst galaxy NGC 253 I. The distribution and transport of cosmic rays

The diffusive and convective CR transport in NGC 253 from the disk into the halo is investigated using the local CR bulk speed. The connection between the CR transport and the galactic wind is outlined. We observed NGC 253 with the VLA at lambda 6.2 cm in a mosaic with 15 pointings. The missing zero-spacing flux density of the VLA mosaic was filled in using observations with the 100-m Effelsberg telescope. We also obtained a new lambda 3.6 cm map from Effelsberg observations and reproduced VLA maps at lambda 20 cm and lambda 90 cm. We find a thin and a thick radio disk with exponential scaleheights of 0.3 kpc and 1.7 kpc at lambda 6.2 cm. The equipartition total magnetic field strength between 7 micro G and 18 micro G in the disk is remarkably high. We use the spectral aging of the cosmic ray electrons (CREs) seen in the vertical profiles of the spectral index to determine a lower limit for the global CR bulk speed as (170+/-70) km/s. The linear correlation between the scaleheights and the CRE lifetimes, as evident from the dumbbell shaped halo, requires a vertical CR transport with a bulk speed of (300+/-30) km/s in the northeastern halo, similar to the escape velocity of 280 km/s. This shows the presence of a "disk wind" in NGC253. In the southwestern halo, the transport is mainly diffusive with a diffusion coefficient of (2.0+/-0.2) 10E29 cm^2 s^-1. In the northeastern halo, the CR transport is convective and more efficient than in the southwestern halo. The luminous material is transported by the disk wind, which can explain the different amounts of extra-planar HI, Halpha, and soft X-ray emission in the two halo parts. Future low-frequency radio observations will provide the data to analyze the vertical velocity profile of galactic winds.

preprint2009arXiv

Obscuring and feeding supermassive black holes with evolving nuclear star clusters

Recently, high resolution observations with the help of the near-infrared adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of high resolution hydrodynamical simulations with the PLUTO-code, we follow the evolution of such clusters, especially focusing on mass and energy feedback from young stars. This leads to a filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. Here, we concentrate on the long-term evolution of the nuclear disc in NGC 1068 with the help of an effective viscous disc model, using the mass input from the large scale simulations and accounting for star formation in the disc. This two-stage modelling enables us to connect the tens of parsec scale region (observable with SINFONI) with the parsec scale environment (MIDI observations). At the current age of the nuclear star cluster, our simulations predict disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 Msun and mass transfer rates through the inner boundary of 0.025 Msun/yr in good agreement with values derived from observations.

preprint2009arXiv

Very Light Magnetized Jets on Large Scales - I. Evolution and Magnetic Fields

Magnetic fields, which are undoubtedly present in extragalactic jets and responsible for the observed synchrotron radiation, can affect the morphology and dynamics of the jets and their interaction with the ambient cluster medium. We examine the jet propagation, morphology and magnetic field structure for a wide range of density contrasts, using a globally consistent setup for both the jet interaction and the magnetic field. The MHD code NIRVANA is used to evolve the simulation, using the constrained-transport method. The density contrasts are varied between η= 10^{-1} and 10^{-4} with constant sonic Mach number 6. The jets are supermagnetosonic and simulated bipolarly due to the low jet densities and their strong backflows. The helical magnetic field is largely confined to the jet, leaving the ambient medium nonmagnetic. We find magnetic fields with plasma β\sim 10 already stabilize and widen the jet head. Furthermore they are efficiently amplified by a shearing mechanism in the jet head and are strong enough to damp Kelvin-Helmholtz instabilities of the contact discontinuity. The cocoon magnetic fields are found to be stronger than expected from simple flux conservation and capable to produce smoother lobes, as found observationally. The bow shocks and jet lengths evolve self-similarly. The radio cocoon aspect ratios are generally higher for heavier jets and grow only slowly (roughly self-similar) while overpressured, but much faster when they approach pressure balance with the ambient medium. In this regime, self-similar models can no longer be applied. Bow shocks are found to be of low excentricity for very light jets and have low Mach numbers. Cocoon turbulence and a dissolving bow shock create and excite waves and ripples in the ambient gas. Thermalization is found to be very efficient for low jet densities.

preprint2007arXiv

Superposition of Quantum and Classical Rotational Motions in Sc2C2@C84 Fullerite

The superposition of the quantum rotational motion (tunneling) of the encapsulated Sc2C2 complex with the classical rotational motion of the surrounding C84 molecule in a powder crystal of Sc2C2@C84 fullerite is investigated by theory. Since the quantum rotor is dragged along by the C84 molecule, any detection method which couples to the quantum rotor (in casu the C2 bond of the Sc2C2 complex also probes the thermally excited classical motion (uniaxial rotational diffusion and stochastic meroaxial jumps) of the surrounding fullerene. The dynamic rotation-rotation response functions in frequency space are obtained as convolutions of quantum and classical dynamic correlation functions. The corresponding Raman scattering laws are derived, the overall shape of the spectra and the width of the resonance lines are studied as functions of temperature. The results of the theory are confronted with experimental low-frequency Raman spectra on powder crystals of Sc2C2@C84 [M. Krause et al., Phys. Rev. Lett. 93, 137403 (2004)]. The agreement of theory with experiment is very satisfactory in a broad temperature range.