
We introduce an affine generalization of counter automata, and analyze their
ability as well as affine finite automata. Our contributions are as follows. We
show that there is a language that can be recognized by exact realtime affine
counter automata but by neither 1way deterministic pushdown automata nor
realtime deterministic kcounter automata. We also show that a certain promise
problem, which is conjectured not to be solved by twoway quantum finite
automata in polynomial time, can be solved by Las Vegas affine finite automata.
Lastly, we show that how a counter helps for affine finite automata by showing
that the language MANYTWINS, which is conjectured not to be recognized by
affine, quantum or classical finite state models in polynomial time, can be
recognized by affine counter automata with onesided boundederror in realtime.

We present an extensive study of the BL Lac object Mrk 501 based on a data
set collected during the multiinstrument campaign spanning from 2009 March 15
to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS,
Whipple 10m, FermiLAT, RXTE, Swift, GASPWEBT and VLBA. We find an increase
in the fractional variability with energy, while no significant interband
correlations of flux changes are found in the acquired data set. The higher
variability in the very high energy (>100 GeV, VHE) gammaray emission and the
lack of correlation with the Xray emission indicate that the highestenergy
electrons that are responsible for the VHE gammarays do not make a dominant
contribution to the ~1 keV emission. Alternatively, there could be a very
variable component contributing to the VHE gammaray emission in addition to
that coming from the synchrotron selfCompton (SSC) scenarios. The space of SSC
model parameters is probed following a dedicated gridscan strategy, allowing
for a wide range of models to be tested and offering a study of the degeneracy
of modeltodata agreement in the individual model parameters. We find that
there is some degeneracy in both the onezone and the twozone SSC scenarios
that were probed, with several combinations of model parameters yielding a
similar modeltodata agreement, and some parameters better constrained than
others. The SSC model gridscan shows that the flaring activity around 2009 May
22 cannot be modeled adequately with a onezone SSC scenario, while it can be
suitably described within a twoindependentzone SSC scenario. The observation
of an electric vector polarization angle rotation coincident with the gammaray
flare from 2009 May 1 resembles those reported previously for low frequency
peaked blazars, hence suggesting that there are many similarities in the
flaring mechanisms of blazars with different jet properties.

We study the multiband variability and correlations of the TeV blazar Mrk
421 on year time scales, which can bring additional insight on the processes
responsible for its broadband emission. We observed Mrk 421 in the very high
energy (VHE) gammaray range with the Cherenkov telescope MAGICI from March
2007 to June 2009 for a total of 96 hours of effective time after quality cuts.
The VHE flux variability is quantified with several methods, including the
Bayesian Block algorithm, which is applied to data from Cherenkov telescopes
for the first time. The 2.3 year long MAGIC light curve is complemented with
data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASPWEBT, OVRO,
and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an
excellent characterisation of the multiband variability and correlations over
year time scales. Mrk 421 was found in different gammaray emission states
during the 2.3 year long observation period. Flares and different levels of
variability in the gammaray light curve could be identified with the Bayesian
Block algorithm. The same behaviour of a quiet and active emission was found in
the Xray light curves measured by Swift/BAT and the RXTE/ASM, with a direct
correlation in time. The behaviour of the optical light curve of GASPWEBT and
the radio light curves by OVRO and Mets\"ahovi are different as they show no
coincident features with the higher energetic light curves and a less variable
emission. The fractional variability is overall increasing with energy. The
comparable variability in the Xray and VHE bands and their direct correlation
during both high and lowactivity periods spanning many months show that the
electron populations radiating the Xray and gammaray photons are either the
same, as expected in the SynchrotronSelfCompton mechanism, or at least
strongly correlated, as expected in electromagnetic cascades.

Completely positive quantum operations are frequently discussed in the
contexts of statistical mechanics and quantum information. They are customarily
given by maps forming positive operatorvalues measures. To intuitively
understand physical meanings of such abstract operations, the method of
phasespace representations is examined. This method enables one to grasp the
operations in terms of the classical statistical notions. As an example of
physical importance, here, the phasespace representation of the completely
positive quantum operation arising from the singlemode subdynamics of the
twomode squeezed vacuum state, which maps from the vacuum state at vanishing
temperature to mixed states with perfect decoherence including the thermal
state, is studied. It is found in the $P$ representation that remarkably this
operation is invertible, implying that coherence lost by the quantum operation
can be recovered.

Mrk 501 is one of the brightest blazars at TeV energies and has been
extensively studied since its first VHE detection in 1996. Our goal is to
characterize in detail the source gammaray emission, together with the
radiotoXray emission, during the nonflaring (low) activity, which is less
often studied than the occasional flaring (high) activity. We organized a
multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This
multiinstrument effort included the most sensitive VHE gammaray instruments
in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes
MAGIC and VERITAS, as well as Swift, RXTE, the FGAMMA, GASPWEBT, and other
collaborations and instruments. Mrk 501 was found to be in a low state of
activity during the campaign, with a VHE flux in the range of 10%20% of the
Crab nebula flux. Nevertheless, significant flux variations were detected with
various instruments, with a trend of increasing variability with energy. The
broadband spectral energy distribution during the two different emission states
of the campaign can be adequately described within the homogeneous onezone
synchrotron selfCompton model, with the (slightly) higher state described by
an increase in the electron number density. This agrees with previous studies
of the broadband emission of this source during flaring and nonflaring states.
We report for the first time a tentative XraytoVHE correlation during a low
VHE activity. Although marginally significant, this positive correlation
between Xray and VHE, which has been reported many times during flaring
activity, suggests that the mechanisms that dominate the Xray/VHE emission
during nonflaringactivity are not substantially different from those that are
responsible for the emission during flaring activity.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV)
gammarays, only three belong to the subclass of Flat Spectrum Radio Quasars
(FSRQs). MAGIC observed FSRQ PKS 1510089 in FebruaryApril 2012 during a high
activity state in the high energy (HE, E>100 MeV) gammaray band observed by
AGILE and Fermi. MAGIC observations result in the detection of a source with
significance of 6.0 sigma. In agreement with the previous VHE observations of
the source, we find no statistically significant variability during the MAGIC
observations in daily, weekly or monthly time scales. The other two known VHE
FSRQs have shown daily scale to subhour variability.
We study the multifrequency behaviour of the source at the epoch of MAGIC
observation, collecting quasisimultaneous data at radio and optical (GASPWEBT
and FGamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA
telescopes), Xray (Swift satellite) and HE gammaray frequencies. The
gammaray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no
hint of a break. The multifrequency light curves suggest a common origin for
the millimeter radio and HE gammaray emission and the HE gammaray flaring
starts when the new component is ejected from the 43GHz VLBA core. The
quasisimultaneous multifrequency SED is modelled with a onezone inverse
Compton model. We study two different origins of the seed photons for the
inverse Compton scattering, namely the infrared torus and a slow sheath
surrounding the jet around the VLBA core. Both models fit the data well.
However, the fast HE gammaray variability requires that within the modelled
large emitting region, there must exist more compact regions. We suggest that
these observed signatures would be most naturally explained by a turbulent
plasma flowing at a relativistic speed down the jet and crossing a standing
conical shock.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV)
gammaray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17,
triggered by the high activity detected by the Fermi Large Area Telescope (LAT)
in high energy (HE; E > 100 MeV) gammarays, poses intriguing questions on the
location of the gammaray emitting region in this flat spectrum radio quasar.
We present multifrequency data of 4C +21.35 collected from centimeter to VHE
during 2010 to investigate the properties of this source and discuss a possible
emission model. The first hint of detection at VHE was observed by MAGIC on
2010 May 3, soon after a gammaray flare detected by FermiLAT that peaked on
April 29. The same emission mechanism may therefore be responsible for both the
HE and VHE emission during the 2010 flaring episodes. Two optical peaks were
detected on 2010 April 20 and June 30, close in time but not simultaneous with
the two gammaray peaks, while no clear connection was observed between the
Xray an gammaray emission. An increasing flux density was observed in radio
and mm bands from the beginning of 2009, in accordance with the increasing
gammaray activity observed by FermiLAT, and peaking on 2011 January 27 in the
mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C
+21.35 for the two periods of the VHE detection and a quiescent state, using a
onezone model with the emission coming from a very compact region outside the
broad line region. The three SEDs can be fit with a combination of synchrotron
selfCompton and external Compton emission of seed photons from a dust torus,
changing only the electron distribution parameters between the epochs. The fit
of the optical/UV part of the spectrum for 2010 April 29 seems to favor an
inner disk radius of <6 gravitational radii, as one would expect from a
prograderotating Kerr black hole.

Theoretical spectra of Xray bursting neutron star (NS) model atmospheres are
widely used to determine the basic NS parameters such as their masses and
radii. We construct accurate NS atmosphere models using for the first time an
exact treatment of Compton scattering via the integral relativistic kinetic
equation. We also compare the results with the previous calculations based on
the Kompaneets operator. We solve the radiation transfer equation together with
the hydrostatic equilibrium equation accounting exactly for the radiation
pressure by electron scattering. We thus construct a new set of planeparallel
atmosphere models in LTE for hot NSs. The models were computed for six chemical
compositions (pure H, pure He, solar H/He mix with various heavy elements
abundances Z = 1, 0.3, 0.1, and 0.01 Z_sun, and three log g = 14.0, 14.3, and
14.6. For each chemical composition and log g, we compute more than 26 model
atmospheres with various luminosities relative to the Eddington luminosity
L_Edd computed for the Thomson crosssection. The maximum relative luminosities
L/L_Edd reach values of up to 1.1 for high gravity models. The emergent spectra
of all models are redshifted and fitted by diluted blackbody spectra in the
320 keV energy range appropriate for the RXTE/PCA. We also compute the color
correction factors f_c. The radiative acceleration g_rad in our luminous,
hotatmosphere models is significantly smaller than in corresponding models
based on the Kompaneets operator, because of the KleinNishina reduction of the
electron scattering crosssection, and therefore formally "superEddington"
model atmospheres do exist. The differences between the new and old model
atmospheres are small for L / L_Edd < 0.8. For the same g_rad / g, the new f_c
are slightly larger (by approximately 1%) than the old values.

All the neutron star (NS) atmosphere models published so far have been
calculated in the "cold plasma approximation", which neglects the relativistic
effects in the radiative processes, such as cyclotron emission/absorption at
harmonics of cyclotron frequency. Here we present new NS atmosphere models
which include such effects. We calculate a set of models for effective
temperatures T_eff =13 MK and magnetic fields B \sim 10^{10}10^{11} G,
typical for the socalled central compact objects (CCOs) in supernova remnants,
for which the electron cyclotron energy E_{c,e} and its first harmonics are in
the observable soft Xray range. Although the relativistic parameters, such as
kT_eff /(m_e c^2) and E_{c,e} /(m_e c^2), are very small for CCOs, the
relativistic effects substantially change the emergent spectra at the cyclotron
resonances, E \approx sE_{c,e} (s=1, 2,...). Although the cyclotron absorption
features can form in a cold plasma due to the quantum oscillations of the
freefree opacity, the shape and depth of these features change substantially
if the relativistic effects are included. In particular, the features acquire
deep Doppler cores, in which the angular distribution of the emergent intensity
is quite different from that in the cold plasma approximation. The relative
contributions of the Doppler cores to the equivalent widths of the features
grow with increasing the quantization parameter b_eff = E_{c,e}/kT_eff and
harmonic number s. The total equivalent widths of the features can reach \sim
150250 eV; they increase with growing b_eff and are smaller for higher
harmonics.

We suggest a new method of calculation of the equilibrium correlation
functions of an arbitrary order for the interacting Fermigas model in the
frame of the static fluctuation approximation (SFA) method. This method based
only on the single and controllable approximation allows to obtain the
socalled fardistance equations (FDEs). These equations connecting the quantum
states of a Fermi particle with variables of the local field operator contains
all necessary information related to calculation of the desired correlation
functions and basic thermodynamic parameters of the manybody system
considered. The basic expressions for the mean energy and heat capacity for
electron gas at low temperatures in the limit of high density were obtained.
All expressions are given in the units of r_s,where r_s determines the ratio of
a mean distance between electrons to the Bohr radius a_0. In these expressions
we calculated the terms of the order r_s and r_s^2, correspondingly. It was
shown also that the SFA allows to find the terms related with high orders of
the decomposition with respect to the parameter r_s.

Xray bursting neutron stars in low mass Xray binaries constitute an
appropriate source class to constrain masses and radii of neutron stars, but a
sufficiently extended set of corresponding model atmospheres is necessary for
these investigations. We computed such a set of model atmospheres and emergent
spectra in a planeparallel, hydrostatic, and LTE approximation with Compton
scattering taken into account. The models were calculated for six different
chemical compositions: pure hydrogen and pure helium atmospheres, and
atmospheres with solar mix of hydrogen and helium, and various heavy element
abundances Z = 1, 0.3, 0.1, and 0.01 Z_sun. For each chemical composition the
models are computed for three values of surface gravity, log g =14.0, 14.3, and
14.6, and for 20 values of the luminosity in units of the Eddington luminosity,
L/L_Edd, in the range 0.0010.98. The emergent spectra of all models are
redshifted and fitted by a diluted blackbody in the RXTE/PCA 320 keV energy
band, and corresponding values of the color correction (hardness factors) f_c
are presented. Theoretical dependences f_c  L/L_Edd can fitted to the observed
dependence K^{1/4}  F of the blackbody normalization K on flux during cooling
stages of Xray bursts to determine the Eddington flux and the ratio of the
apparent neutron star radius to the source distance. If the distance is known,
these parameters can be transformed to the constraints on neutron star mass and
radius. The theoretical atmosphere spectra can also be used for direct
comparison with the observed Xray burst spectra.