
We extend to threedimensional space the approximate M_2 model for the slab
geometry studied in our previous paper. The B_2 model therein, as a special
case of the second order extended quadrature method of moments (EQMOM), is
proved to be globally hyperbolic. The model we proposed here extends EQMOM to
multiple dimensions following the idea to approximate the maximum entropy
closure for the slab geometry case. Like the M_2 closure, the ansatz of the new
model has the capacity to capture both isotropic and beamlike solutions, while
the new model has fluxes in closedform, thus is applicable to practical
numerical simulations. The rotational invariance, realizability, and
hyperbolicity of the model are studied.

Using a Bayesian technology we derived distances and extinctions for over
100,000 red giant stars observed by the Apache Point Observatory Galactic
Evolution Experiment (APOGEE) survey by taking into account spectroscopic
constraints from the APOGEE stellar parameters and photometric constraints from
2MASS, as well as a prior knowledge on the Milky Way. Derived distances are
compared with those from four other independent methods, the Hipparcos
parallaxes, star clusters, APOGEE red clump stars, and asteroseismic distances
from APOKASC (Rodrigues et al. 2014) and SAGA Catalogues (Casagrande et al.
2014). These comparisons covers four orders of magnitude in the distance scale
from 0.02 kpc to 20 kpc. The results show that our distances agree very well
with those from other methods: the mean relative difference between our
Bayesian distances and those derived from other methods ranges from 4.2% to
+3.6%, and the dispersion ranges from 15% to 25%. The extinctions toward all
stars are also derived and compared with those from several other independent
methods: the RayleighJeans Color Excess (RJCE) method, Gonzalez's
twodimensional extinction map, as well as threedimensional extinction maps
and models. The comparisons reveal that, overall, estimated extinctions agree
very well, but RJCE tends to overestimate extinctions for cool stars and
objects with low logg.

A maximum likelihood method is used to deal with the combined estimation of
multimeasurements of a branching ratio, where each result can be presented as
an upper limit. The joint likelihood function is constructed using observed
spectra of all measurements and the combined estimate of the branching ratio is
obtained by maximizing the joint likelihood function. The Bayesian credible
interval, or upper limit of the combined branching ratio, is given in cases
both with and without inclusion of systematic error.

We investigate the photocurrent properties of the topological insulator
(Bi$_{0.5}$Sb$_{0.5}$)$_2$Te$_3$ on SrTiO$_3$substrates. We find reproducible,
submicron photocurrent patterns generated by longrange chemical potential
fluctuations, occurring predominantly at the topological insulator/substrate
interface. We fabricate nanoplowed constrictions which comprise single
potential fluctuations. Hereby, we can quantify the magnitude of the disorder
potential to be in the meV range. The results further suggest a dominating
photothermoelectric current generated in the surface states in such nanoscale
constrictions.

Stellar clusters are regularly used to study the evolution of their host
galaxy. Except for a few nearby galaxies, these studies rely on the
interpretation of integrated cluster properties, especially integrated
photometry observed using multiple filters (i.e. the Spectral Energy
Distribution SED). To allow interpretation of such observations, we present a
large set of GALEV cluster models using the realistic approach of adopting
stochasticallysampled stellar IMFs. We provide models for a wide range of
cluster masses (10^3  2 * 10^5$ Msun), metallicities (2.3 <= [Fe/H] <= +0.18
dex), foreground extinction, and 184 regularly used filters. We analyze various
sets of stochastic cluster SEDs by fitting them with nonstochastic models,
which is the procedure commonly used in this field. We identify caveats and
quantify the fitting uncertainties associated with this standard procedure. We
show that this can yield highly unreliable fitting results, especially for
lowmass clusters.

In this paper, we present surface brightness profiles for 79 globular
clusters in M31, using images observed with {\it Hubble Space Telescope}, some
of which are from new observations. The structural and dynamical parameters are
derived from fitting the profiles to several different models for the first
time. The results show that in the majority of cases, King models fit the M31
clusters as well as Wilson models, and better than S\'{e}rsic models. However,
there are 11 clusters best fitted by S\'{e}rsic models with the S\'{e}rsic
index $n>2$, meaning that they have cuspy central density profiles. These
clusters may be the wellknown corecollapsed candidates. There is a bimodality
in the size distribution of M31 clusters at large radii, which is different
from their Galactic counterparts. In general, the properties of clusters in M31
and the Milky Way fall in the same regions of parameter spaces. The tight
correlations of cluster properties indicate a "fundamental plane" for clusters,
which reflects some universal physical conditions and processes operating at
the epoch of cluster formation.

The PANDA experiment will not use any hardware trigger, i.e. all raw data are
streaming in the data acquisition with a bandwidth of ~280 GB/s. The PANDA
Online System is designed to perform data reduction by a factor of ~800 by
reconstruction algorithms programmed in VHDL (Very High Speed Integrated
Circuit Hardware Description Language) on FPGAs (Field Programmable Gate
Arrays).

The threesite model has been offered as a benchmark for studying the
collider phenomenology of Higgsless models. In this paper we analyze how well
the threesite model performs as a general exemplar of Higgsless models in
describing W_L W_L scattering, and which modifications can make it more
representative. We employ general sum rules relating the masses and couplings
of the KaluzaKlein (KK) modes of the gauge fields in continuum and
deconstructed Higgsless models as a way to compare the different theories. We
show that the size of the fourpoint vertex for the (unphysical)
NambuGoldstone modes and the degree to which the sum rules are saturated by
contributions from the lowestlying KK resonances both provide good measures of
the extent to which a highlydeconstructed theory can accurately describe the
lowenergy physics of a continuum 5d Higgsless model. After comparing the
threesite model to flat and warped continuum models, we analyze extensions of
the threesite model to a longer open linear moose with an additional U(1)
group and to a ring ("BESS' or "hidden local symmetry") model with three sites
and three links. Both cases may be readily analyzed in the framework of the
general sum rules. We demonstrate that W_L W_L scattering in the ring model can
very closely approximate scattering in the continuum models, provided that the
hidden local symmetry parameter `a' is chosen to mimic rhomeson dominance of
pionpion scattering in QCD. The hadron and lepton collider phenomenology of
both extended models is briefly discussed, with a focus on the complementary
information to be gained from precision measurements of the Z' line shape and
ZWW coupling at a highenergy lepton collider.

A deformed relativistic HartreeBogoliubov (DRHB) model is developed aiming
at a proper description of exotic nuclei, particularly deformed ones with large
spatial extension. In order to give an adequate description of both the
contribution of the continuum and the large spatial distribution in exotic
nuclei, the DRHB equations are solved in a WoodsSaxon basis in which the
radial wave functions have proper asymptotic behaviors at large distance from
the nuclear center which is crucial for the formation of halo. The formalism
and the numerical procedure of the DRHB model in a WoodsSaxon basis are
briefly presented.

Episodic activity of quasars is driving growth of supermassive black holes
(SMBHs) via accretion of baryon gas. In this Letter, we develop a simple method
to analyse the duty cycle of quasars up to redshift $z\sim 6$ universe from
luminosity functions (LFs). We find that the duty cycle below redshift $z\sim
2$ follows the cosmic history of star formation rate (SFR) density. Beyond
$z\sim 2$, the evolutionary trends of the duty cycle are just opposite to that
of the cosmic SFR density history, implying the role of feedback from black
hole activity. With the duty cycle, we get the net lifetime of quasars ($z\le
5$) about $\sim 10^9$yrs. Based on the local SMBHs, the mean mass of SMBHs is
obtained at any redshifts and their seeds are of $10^5\sunm$ at the
reionization epoch ($z_{\rm re}$) of the universe through the conservation of
the black hole number density in comoving frame. We find that primordial black
holes ($\sim 10^3\sunm$) are able to grow up to the seeds via a moderate
superEddington accretion of $\sim 30$ times of the critical rate from $z=24$
to $z_{\rm re}$. Highly superEddington accretion onto the primordials is not
necessary.

It is argued that there is a linear correlation between star formation rate
(SFR) and accretion rate for normal bright active galactic nuclei (AGNs).
However, it is still unclear whether this correlation holds for LINERs, of
which the accretion rates are relatively lower than those of normal bright
AGNs. The radiatively inefficient accretion flows (RIAFs) are believed to be
present in these LINERs. In this work, we derive accretion rates for a sample
of LINERs from their hard Xray luminosities based on spectral calculations for
RIAFs. We find that LINERs follow the same correlation between star formation
rate and accretion rate defined by normal bright AGNs, when reasonable
parameters are adopted for RIAFs. It means that the gases feed the black hole
and star formation in these lowluminosity LINERs may follow the same way as
that in normal bright AGNs, which is roughly consistent with recent numerical
simulations on quasar evolution.

Seed black holes formed in the collapse of population III stars have been
invoked to explain the presence of supermassive black holes at high redshift.
It has been suggested that a seed black hole can grow up to $10^{5\sim 6}\sunm$
through highly superEddington accretion for a period of $\sim 10^{6\sim 7}$ yr
between redshift $z=20\sim 24$. We studied the feedback of radiation pressure,
Compton heating and outflow during the seed black hole growth. It is found that
its surrounding medium fueled to the seed hole is greatly heated by Compton
heating. For a supercritical accretion onto a $10^3\sunm$ seed hole, a Compton
sphere (with a temperature $\sim 10^6$K) forms in a timescale of $1.6\times
10^3$yr so that the hole is only supplied by a rate of $10^{3}$ Eddington
limit from the Compton sphere. Beyond the Compton sphere, the kinetic feedback
of the strong outflow heats the medium at large distance, this leads to a
dramatical decrease of the outer Bondi accretion onto the black hole and avoid
the accumulation of the matter. The highly supercritical accretion will be
rapidly halted by the strong feedback. The seed black holes hardly grow up at
the very early universe unless the strong feedback can be avoided.

We model the constraints set on the evolution of dusty starburst galaxies by
the current deep extragalactic surveys performed in the farinfrared with $\it
Spitzer$, and at radio wavelengths with the VLA. Our models fit the number
counts in all the available spectral bands well, and also provide a reasonably
close match to the redshift distribution of the $\it Spitzer$ detections. We
find: 1.) dusty starburst galaxies with infrared burst phases triggered by
galactic interactions at redshift $z \sim 12$ are good candidates to fit the
$\it Spitzer$ results at $24\mu m$, $70 \mu m$ and $160 \mu m$, assuming
plausible strengths for the PAH features for the infrared luminous sources. An
Arp220like spectral energy distribution (SED) for Ultraluminous Infrared
Galaxies (ULIGs) of $L_{\rm ir}>10^{12} L_{\odot}$ and one like that of M82 for
Luminous Infrared Galaxies (LIGs) of $L_{\rm ir}\sim 10^{1112} L_{\odot}$ give
a successful fit to the $\it Spitzer$ $24\mu m$ and $\rm ISOCAM 15\mu m$ number
counts at flux levels of $S_{\nu}< 1 mJy$; 2.) the strong evolution of the
number density of the ULIGs from redshift $z\sim 0$ to $\sim 1$ predicted by
our models is consistent with the current deep $\rm 1.4 GHz$ radio surveys and
accounts for the upturn in the $\rm 1.4 GHz$ differential counts at the submJy
flux level; and 3.) comparing with number counts at near infrared bands, as
well as the background measurements using DIRBE and 2MASS, shows that only a
fraction of the stellar mass in the Universe is included in our models of dusty
starburst mergers at $z \sim 1  2$.

We investigate the radial velocity difference between the narrow
emissionline components of [O III] $\lambda$ 5007 and H$\beta$ in a sample of
150 SDSS NLS1s. Seven "blue outliers" with [O III] blueshifted by more than 250
\kms are found. A strong correlation between the [O III] blueshift and the
Eddington ratio is found for these seven "blue outliers". For the entire
sample, we found a modest correlation between the blueshift and the linewidth
of the narrow component of the [O III] line. The reflected profile of [O III]
indicates two kinematically and physically distinct regions. The [O III]
linewidth depends not only on the bulge stellar gravitational potential, but
also on the central black hole potential.

The missed particleantiparticle degrees of freedom are retrieved and the
corresponding particleantiparticle intrinsic space are introduced to study the
dynamical symmetry of the Dirac particle. As a result, the
particleantiparticle quantum number appears naturally and the Dirac particle
has five quantum numbers instead of four. An antisymmetry (different from the
conventional symmetry) of the Dirac Hamiltonian and a dual symmetry of its
eigenfunctions are explored. The $\hat{\kappa}$ operator of the Dirac equation
in central potentials is found to be the analog of the helicity operator of the
free particlethe alignment of the spin along the angular momentum.

We study the complete K\"{a}hlerEinstein metric of a Hartogs domain
$\widetilde {\Omega}$, which is obtained by inflation of an irreducible bounded
symmetric domain $\Omega $, using a power $N^{\mu}$ of the generic norm of
$\Omega$. The generating function of the K\"{a}hlerEinstein metric satisfies a
complex MongeAmp\`{e}re equation with boundary condition. The domain
$\widetilde {\Omega}$ is in general not homogeneous, but it has a subgroup of
automorphisms, the orbits of which are parameterized by $X\in\lbrack0,1[$. This
allows to reduce the MongeAmp\`{e}re equation to an ordinary differential
equation with limit condition. This equation can be explicitly solved for a
special value $\mu_{0}$ of $\mu$, called the critical exponent. We work out the
details for the two exceptional symmetric domains. The critical exponent seems
also to be relevant for the properties of other invariant metrics like the
Bergman metric; a conjecture is stated, which is proved for the exceptional
domains.

We study the complete K\"{a}hlerEinstein metric of a Hartogs domain
$\widetilde {\Omega}$, which is obtained by inflation of an irreducible bounded
symmetric domain $\Omega $, using a power $N^{\mu}$ of the generic norm of
$\Omega$. The generating function of the K\"{a}hlerEinstein metric satisfies a
complex MongeAmp\`{e}re equation with boundary condition. The domain
$\widetilde {\Omega}$ is in general not homogeneous, but it has a subgroup of
automorphisms, the orbits of which are parameterized by $X\in\lbrack0,1[$. This
allows to reduce the MongeAmp\`{e}re equation to an ordinary differential
equation with limit condition. This equation can be explicitly solved for a
special value $\mu_{0}$ of $\mu$, called the critical exponent. We work out the
details for the two exceptional symmetric domains. The critical exponent seems
also to be relevant for the properties of other invariant metrics like the
Bergman metric; a conjecture is stated, which is proved for the exceptional
domains.

We review the results from the various hydrodynamical and transport models on
the collective flow observables from AGS to RHIC energies. A critical
discussion of the present status of the CERN experiments on hadron collective
flow is given. We emphasize the importance of the flow excitation function from
1 to 50 A$\cdot$GeV: here the hydrodynamic model has predicted the collapse of
the $v_1$flow and of the $v_2$flow at $\sim 10$ A$\cdot$GeV; at 40
A$\cdot$GeV it has been recently observed by the NA49 collaboration. Since
hadronic rescattering models predict much larger flow than observed at this
energy we interpret this observation as evidence for a first order phase
transition at high baryon density $\rho_B$. Moreover, the connection of the
elliptic flow $v_2$ to jet suppression is examined. It is proven experimentally
that the collective flow is not faked by minijet fragmentation. Additionally,
detailed transport studies show that the awayside jet suppression can only
partially ($<$ 50%) be due to hadronic rescattering. Furthermore, the change in
sign of $v_1, v_2$ closer to beam rapidity is related to the occurence of a
high density first order phase transition in the RHIC data at 62.5, 130 and 200
A$\cdot$GeV.

The dynamical evolution of super star clusters (SSCs) moving in the
background of a dark matter halo has been investigated as a possible event
responsible for the formation of bulges in latetype spirals. The underlying
physical processes include sinking of SSCs due to the dynamical friction and
stripping of SSCs on their way to the center. Our model calculations show that
only sinking of circumnuclear SSCs contributes to the formation of galactic
bulges at the early stage. Based on the assumption of a universal density
profile for the dark matter halo, and an isothermal model for the SSCs, our
simulations have yielded bulges that are similar in many aspects to the
observational ones. In particular, the derived surface density profiles can be
well fitted by anexponential structure with nuclear cusps, which is consistent
with HST observations.

Yields of hadrons, their average masses and energies per hadron at the stage
of chemical freezeout in (ultra)relativistic heavyion collisions are analyzed
within the statistical model. The violation of the scaling <E>/<N> = 1 GeV
observed in Au+Au collisions at $\sqrt{s}$ = 130 AGeV is linked to the
formation of resonancerich matter with a considerable fraction of baryons and
antibaryons. The rise of the energyperhadron ratio in baryondominated matter
is discussed. A violation of the scaling condition is predicted for a very
central zone of heavyion collisions at energies around 40 AGeV.

Discrete symmetry breaking and possible restoration at finite temperature $T$
are analysed in 2D GrossNeveu model by the realtime thermal field theory in
the fermion bubble approximation. The dynamical fermion mass $m$ is proven to
be scaleindependent and this fact indicates the equivalence between the
fermion bubble diagram approximation and the mean field approximation used in
the auxialiary scalar field approach. Reproducing of the nonzero critical
temperature $T_c=0.567 m(0)$, ($m(0)$ is the dynamical fermion mass at T=0),
shows the equivalence between the realtime and the imaginarytime thermal
field theory in this problem. However, in the realtime formalism, more results
including absence of scalar bound state, the equation of criticality curve of
chemical potentialtemperature and the $\ln(T_c/T)$ behavior of $m^2$ at
$T\stackrel{<}{\sim} T_c$ can be easily obtained. The last one indicates the
secondorder phase transition feature of the symmetry restoration.

Based on a general analysis of Green functions in the realtime thermal field
theory, we have proven that the fourpoint amputated functions in a NJL model
in the fermion bubble diagram approximation behave like usual twopoint
functions. We expound the thermal transformations of the matrix propagator for
a scalar bound state in the $F\bar{F}$ basis and in the $RA$ basis. The
resulting physical causal, advanced and retarded propagator are respectively
identical to corresponding ones derived in the imaginarytime formalism and
this shows once again complete equivalence of the two formalisms of thermal
field theory on the discussed problem in the NJL model.

The temporal and spectral properties of terrestrial gammaray flashes (TGFs)
are studied. The delay of low energy photons relative to high energy ones in
the gammaray variations of the TGFs with high signal to noise ratio has been
revealed by an approach of correlation analysis in the time domain on different
time scales. The temporal structures of the TGFs in high energy band are
usually narrower than that in low energy band. The spectral hardness has a
general trend of decreasing with time in a flash. The observed temporal and
spectral characteristics give constraints and valuable hints on the flash
production mechanism.

The SchwingerDyson equation of fermion selfenergy in the linearization
approximation is solved exactly in a theory with gauge and effective
fourfermion interactions. Different expressions for the indepedent solutions
which respectively submit to irregular and regular ultraviolet boundary
condition are derived and expounded.

By rigorous reanalysis of the results, we have proven that the propagators at
finite temperature for scalar bound states in onegeneration fermion condensate
scheme of electroweak symmetry breaking are in fact identical in the
imaginarytime and the realtime formalism. This dismisses the doubt about
possible discrepancy between the two formalisms in this problem. Identity of
the derived thermal transformation matrices of the realtime matrix propagators
for scalar bound states without and with chemical potential and the ones for
corresponding elementary scalar particles shows similarity of thermodynamic
property between the two types of particles. Only one former inference is
modified, i.e. when the two flavors of fermions have unequal nonzero masses,
the amplitude of the composite Higgs particle will decay instead grow in time.