We study the structure of accretion disks around supermassive black holes in the radial range $30\sim 100$ gravitational radii, using a three dimensional radiation magneto-hydrodynamic simulation. For typical conditions in this region of Active Galactic Nuclei (AGN), the Rosseland mean opacity is expected to be larger than the electron scattering value. We show that the iron opacity bump causes the disk to be convective unstable. Turbulence generated by convection puffs up the disk due to additional turbulent pressure support and enhances the local angular momentum transport. This also results in strong fluctuations in surface density and heating of the disk. The opacity drops with increasing temperature and convection is suppressed. The disk cools down and the whole process repeats again. This causes strong oscillations of the disk scale height and luminosity variations by more than a factor of $\approx 3-6$ over a few years' timescale. Since the iron opacity bump will move to different locations of the disk for black holes with different masses and accretion rates, we suggest that this is a physical mechanism that can explain the variability of AGN with a wide range of amplitudes over a time scale of years to decades.

Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events.

A luminous body embedded in an accretion disk can generate asymmetric density perturbations that lead to a net torque and thus orbital migration of the body. Linear theory has shown that this heating torque gives rise to a migration term linear in the body's mass that can oppose or even reverse that arising from the sum of gravitational Lindblad and co-orbital torques. We use high-resolution local simulations in an unstratified disk to assess the accuracy and domain of applicability of the linear theory. We find agreement between analytic and simulation results to better than 10\% in the appropriate regime (low luminosity, low thermal conductivity), but measure deviations in the non-linear (high luminosity) regime and in the high thermal conductivity regime. In the non-linear regime, linear theory overpredicts the acceleration due to the heating torque, which we find to be due to the neglect of non-linear terms in the heat flux. In the high thermal conductivity regime linear theory underpredicts the acceleration, although here both non-linear and computational constraints play a role. We discuss the impact of the heating torque for the evolution of low-mass planets in protoplanetary disks, and for massive stars or accreting compact objects embedded in AGN disks. For the latter case, we show that the thermal torque is likely to be the dominant physical effect at disk radii where the optical depth drops below a critical value.

An observation-based Galactic Cosmic Ray (GCR) spectral model for heavy nuclei is developed. Zhao and Qin (J. Geophys. Res. Space Phys.118, 1837(2013)) proposed an empirical elemental GCR spectra model for nuclear charge 5-28 over the energy range from 30 to 500 MeV/nuc, which is proved to be successful in predicting yearly averaged GCR heavy nuclei spectra.Based on the latest highly statistically precise measurements from ACE/CRIS,a further elemental GCR model with monthly averaged spectra is presented. The model can reproduce the past and predict the futureGCR intensity monthly by correlating model parameters with thecontinuous sunspot number (SSN) record. The effects of solar activity on GCR modulation are considered separately for odd and even solar cycles. Compared with other comprehensive GCR models, our modeling results are satisfyingly consistent with the GCR spectral measurements from ACE/SIS and IMP-8, and have comparable prediction accuracy as the Badhwar & O'Neill 2014 model.A detailed error analysis is also provided.Finally, the GCR carbon and iron nuclei fluxes for the subsequent two solar cycles (SC 25 and 26) are predicted and they show a potential trend in reduced flux amplitude, which is suspected to be relevant to possible weak solar cycles.

Site-testing is crucial for achieving the goal of scientific research and analysis of meteorological and optical observing conditions is one of the basic tasks of it. As one of three potential sites to host 12-meter Large Optical/infrared Telescope (LOT), Muztagh-ata site which is located on the Pamirs Plateau in west China's Xinjiang began its site-testing task in the spring of 2017. In this paper, we firstly start with an introduction to the site and then present a statistical analysis of the ground-level meteorological properties such as air temperature, barometric pressure, relative humidity, wind speed and direction, recorded by automatic weather station with standard meteorological sensors for two-year long. We also show the monitoring results of sky brightness during this period.

In this article, we present a detailed analysis of the statistical properties of seeing for the Muztagh-ata site which is the candidate site for hosting future Chinese Large Optical/infrared Telescope (LOT) project. The measurement was obtained with Differential Image Motion Monitor (DIMM) from April 2017 to November 2018 at different heights during different periods. The median seeing at 11 meters and 6 meters are very close but different significantly from that on the ground. We mainly analyzed the seeing at 11 meters monthly and hourly, having found that the best season for observing was from late autumn to early winter and seeing tended to improve during the night only in autumn. The analysis of the dependence on temperature inversion, wind speed, direction also was made and the best meteorological conditions for seeing is given.

Solar system materials are variably depleted in moderately volatile elements (MVEs) relative to the proto-solar composition. To address the origin of this MVE depletion, we conducted a systematic study of high-precision K isotopic composition on 16 carbonaceous chondrites (CCs) of types CM1-2, CO3, CV3, CR2, CK4-5 and CH3 and 28 ordinary chondrites (OCs) covering petrological types 3 to 6 and chemical groups H, L, and LL. We observed significant overall K isotope (delta41K) variations (-1.54 to 0.70 permil). The K isotope compositions of CCs are largely higher than the Bulk Silicate Earth (BSE) value, whereas OCs show typically lower values than BSE. Neither CCs nor OCs show resolvable correlations between K isotopes and chemical groups, petrological types, shock levels, exposure ages, fall or find occurrence, or terrestrial weathering. The lack of a clear trend between K isotopes and K content indicates that the K isotope fractionations were decoupled from the relative elemental K depletions. The range of K isotope variations in the CCs is consistent with a four-component (chondrule, refractory inclusion, matrix and water) mixing model that is able to explain the bulk elemental and isotopic compositions of the main CC groups, but requires a fractionation in K isotopic compositions in chondrules. We propose that the major control of the isotopic compositions of group averages is condensation or vaporization in nebular environments that is preserved in the compositional variation of chondrules. Parent-body processes (aqueous alteration, thermal metamorphism, and metasomatism) can mobilize K and affect the K isotopes in individual samples. In the case of the OCs, the full range of K isotopic variations can only be explained by the combined effects of the size and relative abundances of chondrules, parent-body aqueous and thermal alteration.

FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the "surface accretion" mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3.

We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equations on an unstructured Voronoi mesh. We then solve the source part of the radiative transfer equations following Jiang et al (2014) using an implicit solver, and couple this to the hydrodynamic equations. The use of an implicit solver ensure reliable convergence and preserves the conservation properties of these equations even in situations where the source terms are stiff due to the small coupling timescales between radiation and matter. We present the results of a limited number of test cases (energy conservation, momentum conservation, dynamic diffusion, linear waves, crossing beams, and multiple shadows) to show convergence with analytic results and numerical stability. We also show that it produces qualitatively the correct results in the presence of multiple sources in the optically thin case.

The Large Optical/infrared Telescope (LOT) is a ground-based 12m diameter optical/infrared telescope which is proposed to be built in the western part of China in the next decade. Based on satellite remote sensing data, along with geographical, logistical and political considerations, three candidate sites were chosen for ground-based astronomical performance monitoring. These sites include: Ali in Tibet, Daocheng in Sichuan, and Muztagh Ata in Xinjiang. Up until now, all three sites have continuously collected data for two years. In this paper, we will introduce this site testing campaign, and present its monitoring results obtained during the period between March 2017 and March 2019.

A comprehensive study of multifrequency correlations can shed light on the nature of variation for blazars. In this work, we collect the long-term radio, optical and $\gamma$-ray light curves of PKS 1502+106. After performing the localized cross-correlation function analysis, we find that correlations between radio and $\gamma$-ray or $V$ band are beyond the $3\sigma$ significance level. The lag of the $\gamma$-ray relative to 15 GHz is $-60^{+5}_{-10}$ days, translating to a distance $3.18^{+0.50}_{-0.27}$ parsec (pc) between them. Within uncertainties, the locations of the $\gamma$-ray and optical emitting regions are roughly the same, and are away from the jet base within $1.2$ pc. The derived magnetic field in optical and $\gamma$-ray emitting regions is about $0.36$ G. The logarithm of $\gamma$-ray flux is significantly linearly correlated with that of $V$ band fluxes, which can be explained by the synchrotron self-Compton (SSC) process, the external Compton (EC) processes, or the combination of them. We find a significant linear correlation in the plot of $\log\prod$ (polarization degree) versus $\log \nu F_{\nu}$ at $V$ band, and use the empirical relation $\Pi \sim \sin^n \theta'$ ($\theta'$ is the observing angle in the comoving frame blob) to explain it. The behaviors of color index (generally redder when brighter at the active state) and $\gamma$-ray spectral index (softer when brighter) could be well explained by the twisted jet model. These findings suggest that the curvature effect (mainly due to the change of the viewing angle) is dominant in the variation phenomena of fluxes, spectral indices, and polarization degrees for PKS 1502+106.

We use multi-frequency radiation hydrodynamics (rad-HD) to simulate radiative acceleration of a spherically symmetric stellar wind. We demonstrate the rad-HD capabilities of Athena++ for a series of test problems with multi-group radiation transfer. We then model the radiative transfer of a single spectral line through a spherically symmetric, isothermal, "CAK"-type line driven wind. We find that correctly accounting for the Doppler shift of the absorbed radiation, the force is well described by the analytic Sobolev line transfer in the supersonic parts of the solution where the flow is stationary and the effects of Abbott waves is negligible. Unlike in the analytic, steady-state solution re-radiation is important and leads to non-trivial radiation energy density and fluxes in the outer parts of the wind. We discuss a variety of applications to these multi-group methods that are currently computationally tractable.

The typical blazar S5 0716$+$714 is very interesting due to its rapid and large amplitude variability and high duty cycle of micro-variability in optical band. We analyze the observations in I, R and V bands obtained with the $1.0m$ telescope at Weihai observatory of Shandong University from 2011 to 2018. The model of synchrotron radiation from turbulent cells in a jet has been proposed as a mechanism for explaining micro-variability seen in blazar light curves. Parameters such as the sizes of turbulent cells, the enhanced particle densities, and the location of the turbulent cells in the jet can be studied using this model. The model predicts a time lag between variations as observed in different frequency bands. Automatic model fitting method for micro-variability is developed, and the fitting results of our multi-frequency micro-variability observations support the model. The results show that both the amplitude and duration of flares decomposed from the micro-variability light curves confirm to the log-normal distribution. The turbulent cell size is within the range of about 5 to 55 AU, and the time lags of the micro-variability flares between the I-R and R-V bands should be several minutes. The time lags obtained from the turbulence model are consistent with the fitting statistical results, and the time lags of flares are correlated with the time lags of the whole light curve.

When photons from distant galaxies and stars pass through our neighboring environment, the wavelengths of the photons would be shifted by our local gravitational potential. This local gravitational redshift effect can potentially have an impact on the measurement of cosmological distance-redshift relation. Using available supernovae data, Wojtak et al [1] found seemingly large biases of cosmological parameters for some extended models (non-flat $\Lambda$CDM, $w$CDM, etc.). Huang [2] pointed out that, however, the biases can be reduced to a negligible level if cosmic microwave background (CMB) data are added to break the strong degeneracy between parameters in the extended models. In this article we forecast the cosmological bias due to local gravitational redshifts for a future WFIRST-like supernovae survey. We find that the local gravitational redshift effect remains negligible, provided that CMB data or some future redshift survey data are added to break the degeneracy between parameters.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Recent detection of superflares on solar-type stars by Kepler mission raised a possibility that they can be associated with energetic coronal mass ejections (CMEs) and energetic particle events (SEPs). These space weather events can impact habitability of exoplanets around these stars. Here we use the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, to model the time intensity profile and spectrum of SEPs accelerated at CME-driven shocks from stars of different ages traced by their rotation rates. We consider a solar-like (G-type) star with 6 different rotation rates varying from 0.5 to 3.0 times rotation speed of the sun. In all 6 cases, a fast CME is launched with the same speed of 1500 km/sec and the resulting time intensity profiles at 3 locations and and energy spectra at 5 locations at 1 AU are obtained. The maximum particle energy at the shock front as a function of r is also shown. Our results suggest that within 0.8 AU the maximum particle energy at the shock front increases with the rotation rate of the star. However, event integrated spectra for the five selected locations along the CME path show complicated patterns. This is because the Parker magnetic field for rapidly rotating stars is more tightly winded. Our results can be used in estimating the radiation environments of terrestrial-type exoplanets around solar-type stars.

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.

SS433 is the only Galactic binary system known to accrete at highly super-critical rates, analogous to tidal disruption events, and needed to explain the mass of some high redshift quasars. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but has not yet been possible due to obscuration. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its super-orbital precession period. Spectral-timing tools have allowed us to confirm that the hard X-ray emission from the inner regions is scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of >= 3x10^37 erg/s and that, if viewed face on, the apparent luminosity would be > 1x10^39 erg/s, confirming its long-suspected nature as an ultraluminous X-ray source (ULX). A lag due to absorption by Fe XXV/XXVI in outflowing material travelling at least 0.14-0.29c matches absorption lines seen in ULXs and - in future - will allow us to map a super-critical outflow for the first time.

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

We collect long term $\gamma$-ray, optical and radio $15$ GHz light curves of quasar object PMN J2345-1555. The correlation analyses between them are performed via the local cross-correlation function (LCCF). We found that all the optical $V$, $R$ band and the infrared $J$ band are correlated with the radio 15 GHz at beyond $3\sigma$ significance level, and the lag times are $-221.81^{+6.26}_{-6.72}$, $-201.38^{+6.42}_{-6.02}$ and $-192.27^{+8.26}_{-7.37}$ days, respectively. The $\gamma$-ray is strongly correlated with optical, but weakly correlated with the radio. We present that time lags between different frequencies can be used as an alternative parameter to derive the core-shift measurement. For this target, the magnetic field and particle density at 1 parsec in jet are derived to be $0.61$ Gauss and $1533/\gamma_{\rm min}$ cm$^{-3}$, respectively. The black hole mass and the 15 GHz core position in jet are estimated to be $10^{8.44} {\rm M}_{\odot}$ and $30$ parsec, respectively. The lag times enable us to derive that the optical and the $\gamma$-ray emitting regions coincide, which are located at $4.26^{+0.83}_{-0.79}$ pc away from 15 GHz core position in jet and beyond the broad line region (BLR). We found that a $3\sigma$ correlation between the color index and the radio light curve, which indicates that opacity may play an important role in the variation. The $\delta V-\delta R$ behaviors are complex, while the $R-J$ shows a bluer when brighter trend. As hinted from radio images, we proposed a positional dependent spectral index model to explain the color index behaviors, which is complementary for the shock in jet model. The curvature effects and contribution from accretion disk may also affect variables of blazars in many aspects.

(Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} \approx 2 - 4 \times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} \approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{\odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.

We investigated 64 pairs of interacting-CME events identified by the simultaneous observations of SOHO and STEREO spacecraft from 2010 January to 2014 August, to examine the relationship between the large SEP events in the energy of 25-60MeV and the properties of the interacting CMEs. We found that during CME interactions the large SEP events in this study were all generated by CMEs with the presence of enhanced type II radio bursts, which also have wider longitudinal distributions comparing to events with the absence of type II radio burst or its enhancement (almost associated with small SEP events). It seems that the signature of type II radio bursts enhancement is a good discriminator between large SEP and small or none SEP event producers during CME interactions. The type II radio burst enhancement is more likely to be generated by CME interactions, with the main CME having larger speed (v), angular width (WD), mass (m) and kinetic energy (E_k), that taking over the preceding CMEs which also have higher v, WD, m and E_k, than those preceding CMEs in CME pairs missing the type II radio bursts or enhancements. Our analysis also revealed that the intensities of associated SEP events correlate negatively with the intersection height of the two CMEs. Most of type-II-enhanced events and SEP events are coincidentally and almost always made by the fast and wide main CMEs that sweeping fully over a relatively slower and narrower preceding CMEs. We suggest that a fast CME with enough energy completely overtaking a relatively narrower preceding CME, especially in low height, can drive more energetic shock signified by the enhanced type II radio bursts. The shock may accelerate ambient particles and lead to large SEP event more easily.

We investigate the equilibrium points and orbits around asteroid 1333 Cevenola by considering the full gravitational potential caused by the 3D irregular shape. The gravitational potential and effective potential of asteroid 1333 Cevenola are calculated. The zero velocity curves for a massless particle orbiting in the gravitational environment have been discussed. The linearized dynamic equation, the characteristic equation, and the conserved quantity of the equilibria for the large size ratio binary asteroid system have been derived. It is found that there are totally five equilibrium points close to 1333 Cevenola. The topological cases of the outside equilibrium points have a staggered distribution. The simulation of orbits in the full gravitational potential caused by the 3D irregular shape of 1333 Cevenola shows that the moonlets orbit is more likely to be stable if the orbit inclination is small.

This study presents a study of equilibrium points, periodic orbits, stabilities, and manifolds in a rotating plane symmetric potential field. It has been found that the dynamical behaviour near equilibrium points is completely determined by the structure of the submanifolds and subspaces. The non-degenerate equilibrium points are classified into twelve cases. The necessary and sufficient conditions for linearly stable, non resonant unstable and resonant equilibrium points are established. Furthermore, the results show that a resonant equilibrium point is a Hopf bifurcation point. In addition, if the rotating speed changes, two non degenerate equilibria may collide and annihilate each other. The theory developed here is lastly applied to two particular cases, motions around a rotating, homogeneous cube and the asteroid 1620 Geographos. We found that the mutual annihilation of equilibrium points occurs as the rotating speed increases, and then the first surface shedding begins near the intersection point of the x axis and the surface. The results can be applied to planetary science, including the birth and evolution of the minor bodies in the Solar system, the rotational breakup and surface mass shedding of asteroids, etc.

The rotational speeds of irregular shaped minor celestial bodies can be changed by the YORP effect. This variation in speed can make the numbers, positions, stabilities, and topological cases of the minor body s relative equilibrium points vary. The numbers of relative equilibrium points can be reduced through the collision and annihilation of relative equilibrium points, or increase through the creation and separation of relative equilibrium points. Here we develop a classification system of multiple annihilation behaviors of the equilibrium points for irregular shaped minor celestial bodies. Most minor bodies have five equilibrium points; there are twice the number of annihilations per equilibrium point when the number of equilibrium points is between one and five. We present the detailed annihilation classifications for equilibria of objects which have five equilibrium points. Additionally, the annihilation classification for the seven equilibria of Kleopatra shaped objects and the nine equilibria of Bennu shaped objects are also discussed. Equilibria of different objects fall into different annihilation classifications. By letting the rotational speeds vary, we studied the annihilations and creations of relative equilibria in the gravitational field of ten minor bodies.

We used binary octahedrons to investigate the dynamical behaviors of binary asteroid systems. The mutual potential of the binary polyhedron method is derived from the fourth order to the sixth order. The irregular shapes, relative orbits, attitude angles, as well as the angular velocities of the binary asteroid system are included in the model. We investigated the relative trajectory of the secondary relative to the primary, the total angular momentum and total energy of the system, the three-axis attitude angular velocity of the binary system, as well as the angular momentum of the two components. The relative errors of total angular momentum and total energy indicate the calculation has a high precision. It is found that the influence of the orbital and attitude motion of the primary from the gravitational force of the secondary is obvious. This study is useful for understanding the complicated dynamical behaviors of the binary asteroid systems discovered in our Solar system.

We are interested in stable periodic orbits for spacecrafts in the gravitational field of minor celestial bodies. The stable periodic orbits around minor celestial bodies are useful not only for the mission design of the deep space exploration, but also for studying the long-time stability of small satellites in the large-size-ratio binary asteroids. The irregular shapes and gravitational fields of the minor celestial bodies are modeled by the polyhedral model. Using the topological classifications of periodic orbits and the grid search method, the stable periodic orbits can be calculated and the topological cases can be determined. Furthermore, we find five different types of stable periodic orbits around minor celestial bodies: A) Stable periodic orbits generated from the stable equilibrium points outside the minor celestial body, B) Stable periodic orbits continued from the unstable periodic orbits around the unstable equilibrium points, C) Retrograde and nearly circular periodic orbits with zero-inclination around minor celestial bodies, D) Resonance periodic orbits, E) Near-surface inclined periodic orbits. We take asteroid 243 Ida, 433 Eros, 6489 Golevka, 101955 Bennu, and the comet 1P/Halley for examples.

We studied the dynamical environment in the vicinity of the primary of the binary asteroid. The gravitational field of the primary is calculated by the polyhedron model with observational data of the irregular shape. The equilibrium points, zero velocity surfaces, as well as Jacobi integral have been investigated. The results show that the deviations of equilibrium points are large from the principal axes of moment of inertia. We take binary asteroid 41 Daphne and S 2008 41 1 for example. The distribution of topological cases of equilibrium points around 41 Daphne is different from other asteroids. The topological cases of the outer equilibrium points E1 E4 are Case 2, Case 5, Case 2, and Case 1. The topological case of the inner equilibrium point E5 is Case 1. Among the four outer equilibrium points E1 E4, E4 is linearly stable and other outer equilibrium points are unstable. Considering the shape variety of the body from Daphne to a sphere, we calculated the zero velocity surfaces and the locations as well as eigenvalues of equilibrium points. It is found that the topological case of the outer equilibrium point E2 change from Case 5 to Case 1, and its stability change from unstable to linearly stable. Using the gravitational force acceleration calculated by the polyhedron model with the irregular shape, we simulated the orbit for the moonlet in the potential of 41 Daphne.

We investigate the orbital stability close to the unique L4-point Jupiter binary Trojan asteroid 624 Hektor. The gravitational potential of 624 Hektor is calculated using the polyhedron model with observational data of 2038 faces and 1021 vertexes. Previous studies have presented three different density values for 624 Hektor. The equilibrium points in the gravitational potential of 624 Hektor with different density values have been studied in detail. There are five equilibrium points in the gravitational potential of 624 Hektor no matter the density value. The positions, Jacobian, eigenvalues, topological cases, stability, as well as the Hessian matrix of the equilibrium points are investigated. For the three different density values the number, topological cases, and the stability of the equilibrium points with different density values are the same. However, the positions of the equilibrium points vary with the density value of the asteroid 624 Hektor. The outer equilibrium points move away from the asteroid s mass center when the density increases, and the inner equilibrium point moves close to the asteroid s mass center when the density increases. There exist unstable periodic orbits near the surface of 624 Hektor. We calculated an orbit near the primary s equatorial plane of this binary Trojan asteroid, the results indicate that the orbit remains stable after 28.8375 d.

We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $\gamma$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_{\gamma}=344.0\pm16.4$ days are obtained for radio 15 GHz and $\gamma$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $\gamma$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $\gamma$-ray is produced by both the SSC and EC processes.

We perform two-dimensional and three-dimensional radiation hydrodynamic simulations to study cold clouds accelerated by radiation pressure on dust in the environment of rapidly star-forming galaxies dominated by infrared flux. We utilize the reduced speed of light approximation to solve the frequency-averaged, time-dependent radiative transfer equation. We find that radiation pressure is capable of accelerating the clouds to hundreds of kilometers per second while remaining dense and cold, consistent with observations. We compare these results to simulations where acceleration is provided by entrainment in a hot wind, where the momentum injection of the hot flow is comparable to the momentum in the radiation field. We find that the survival time of the cloud accelerated by the radiation field is significantly longer than that of a cloud entrained in a hot outflow. We show that the dynamics of the irradiated cloud depends on the initial optical depth, temperature of the cloud, and the intensity of the flux. Additionally, gas pressure from the background may limit cloud acceleration if the density ratio between the cloud and background is $\lesssim 10^{2}$. In general, a 10 pc-scale optically thin cloud forms a pancake structure elongated perpendicular to the direction of motion, while optically thick clouds form a filamentary structure elongated parallel to the direction of motion. The details of accelerated cloud morphology and geometry can also be affected by other factors, such as the cloud lengthscale, the reduced speed of light approximation, spatial resolution, initial cloud structure, and the dimensionality of the run, but these have relatively little affect on the cloud velocity or survival time.

Numerical solutions of the cosmic-ray (CR) magneto-hydrodynamic equations are dogged by a powerful numerical instability, which arises from the constraint that CRs can only stream down their gradient. The standard cure is to regularize by adding artificial diffusion. Besides introducing ad-hoc smoothing, this has a significant negative impact on either computational cost or complexity and parallel scalings. We describe a new numerical algorithm for CR transport, with close parallels to two moment methods for radiative transfer under the reduced speed of light approximation. It stably and robustly handles CR streaming without any artificial diffusion. It allows for both isotropic and field-aligned CR streaming and diffusion, with arbitrary streaming and diffusion coefficients. CR transport is handled explicitly, while source terms are handled implicitly. The overall time-step scales linearly with resolution (even when computing CR diffusion), and has a perfect parallel scaling. It is given by the standard Courant condition with respect to a constant maximum velocity over the entire simulation domain. The computational cost is comparable to that of solving the ideal MHD equation. We demonstrate the accuracy and stability of this new scheme with a wide variety of tests, including anisotropic streaming and diffusion tests, CR modified shocks, CR driven blast waves, and CR transport in multi-phase media. The new algorithm opens doors to much more ambitious and hitherto intractable calculations of CR physics in galaxies and galaxy clusters. It can also be applied to other physical processes with similar mathematical structure, such as saturated, anisotropic heat conduction.

The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation.

We present the analysis and computational results for the inclination relative effect of moonlets of triple asteroidal systems. Perturbations on moonlets due to the primarys non-sphericity gravity, the solar gravity, and moonlets relative gravity are discussed. The inclination vector for each moonlet follows a periodic elliptical motion; the motion period depends on the moonlets semi-major axis and the primarys J2 perturbations. Perturbation on moonlets from the Solar gravity and moonlets relative gravity makes the motion of the x component of the inclination vector of moonlet 1 and the y component of the inclination vector of moonlet 2 to be periodic.The mean motion of x component and the y component of the inclination vector of each moonlet forms an ellipse. However, the instantaneous motion of x component and the y component of the inclination vector may be an elliptical disc due to the coupling effect of perturbation forces. Furthermore, the x component of the inclination vector of moonlet 1 and the y component of the inclination vector of moonlet 2 form a quasi-periodic motion. Numerical calculation of dynamical configurations of two triple asteroidal systems (216) Kleopatra and (153591) 2001 SN263 validates the conclusion.

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equibria or the stable periodic orbits will not escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric froce, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we clacluated the strengths of different forces.

We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5\times 10^8M_{\odot}$ black hole with accretion rates varying from $\sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $\sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $\sim 0.1-0.4c$. When the accretion rate is smaller than $\sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $\sim 5\%-7\%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1\%$. We always find the kinetic energy luminosity associated with the outflow is only $\sim 15\%-30\%$ of the radiative luminosity. The mass flux lost in the outflow is $\sim 15\%-50\%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.

The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observations of Sgr A* with Markov Chain Monte Carlo sampling, self-consistently modelling the 2-D inflow-outflow solution for the first time. We find the temperature and density at flow onset are consistent with the origin of the gas in the stellar winds of massive stars in the vicinity of Sgr A*. We place the first observational constraints on the angular momentum of the gas and estimate the centrifugal radius, r$_c$ $\approx$ 0.056 r$_b$ $\approx8\times10^{-3}$ pc, where r$_b$ is the Bondi radius. Less than 1\% of the inflowing gas accretes onto the SMBH, the remainder being ejected in a polar outflow. For the first time, we decouple the quiescent point-like emission from the spatially extended flow. We find this point-like emission, accounting for $\sim4$\% of the quiescent flux, is spectrally too steep to be explained by unresolved flares, nor Bremsstrahlung, but is likely a combination of a relatively steep synchrotron powerlaw and the high energy tail of inverse-Compton emission. With this self-consistent model of the accretion flow structure, we make predictions for the flow dynamics and discuss how future X-ray spectroscopic observations can further our understanding of the Sgr A* accretion flow.

We study the time lags between the continuum emission of quasars at different wavelengths, based on more than four years of multi-band ($g$, $r$, $i$, $z$) light-curves in the Pan-STARRS Medium Deep Fields. As photons from different bands emerge from different radial ranges in the accretion disk, the lags constrain the sizes of the accretion disks. We select 240 quasars with redshifts $z \approx 1$ or $z \approx 0.3$ that are relatively emission line free. The light curves are sampled from day to month timescales, which makes it possible to detect lags on the scale of the light crossing time of the accretion disks. With the code JAVELIN, we detect typical lags of several days in the rest frame between the $g$ band and the $riz$ bands. The detected lags are $\sim 2-3$ times larger than the light crossing time estimated from the standard thin disk model, consistent with the recently measured lag in NGC5548 and micro-lensing measurements of quasars. The lags in our sample are found to increase with increasing luminosity. Furthermore, the increase in lags going from $g-r$ to $g-i$ and then to $g-z$ is slower than predicted in the thin disk model, particularly for high luminosity quasars. The radial temperature profile in the disk must be different from what is assumed. We also find evidence that the lags decrease with increasing line ratios between ultraviolet FeII lines and MgII, which may point to changes in the accretion disk structure at higher metallicity.

LX Ser is a deeply eclipsing cataclysmic variable with an orbital period of $0.^d 1584325$. Sixty two new eclipse times were determined by our observations and the AAVSO International Data base. Combining all available eclipse times, we analyzed the O-C behavior of LX Ser. We found that the O-C diagram of LX Ser shows a sinusoidal oscillation with a period of 22.8 yr and an amplitude of 0.00035 days. Two mechanisms (i.e., the Applegate mechanism and the light travel time effect) are applied to explain the cyclic modulation. We found that the Applegate mechanism is difficult to explain the cyclic oscillation in the orbital period. Therefore, the cyclic period change is most likely to be caused by the light travel time effect due to the presence of a third body. The mass of the tertiary component was determined to be $M_3\sim7.5 M_{Jup}$. We supposed that the tertiary companion is plausible a giant planet. The stability of the giant planet was checked, and we found that the multiple system is stable.

We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfv\'en velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $\approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

The discovery of binary and triple asteroids in addition to the execution of space missions to minor celestial bodies in the past several years have focused increasing attention on periodic orbits around irregular-shaped celestial bodies. In the present work, we adopt a polyhedron shape model for providing an accurate representation of irregular-shaped bodies, and employ the model to calculate their corresponding gravitational and effective potentials. We also investigate the characteristics of periodic orbit families and the continuation of periodic orbits. We prove a fact, which provides a conserved quantity that permits restricting the number of periodic orbits in a fixed energy curved surface about an irregular-shaped body. The collisions of Floquet multipliers are maintained during the continuation of periodic orbits around the comet 1P/Halley. Multiple bifurcations in the periodic orbit families about irregular-shaped bodies are also discussed. Three bifurcations in the periodic orbit family have been found around the asteroid 216 Kleopatra, which include two real saddle bifurcations and one period-doubling bifurcation.

We show that in a multi-Higgs model in which one Higgs fits the LHC 125 GeV state, one or more of the other Higgs bosons can mediate DM-nucleon interactions with maximal DM isospin violation being possible for appropriate Higgs-quark couplings, independent of the nature of DM. We then consider the explicit example of a Type II two-Higgs-doublet model, identifying the h or H as the 125 GeV state while the H or h, respectively, mediates DM-nucleon interactions. Finally, we show that if a stable scalar, S, is added then it can be a viable light DM candidate with correct relic density while obeying all direct and indirect detection limits.

This manuscript considers the main features of the nonlinear dynamics of multiple irregular celestial body systems. The gravitational potential, static electric potential, and magnetic potential are considered. Based on the three established potentials, we show that three conservative values exist for this system, including a Jacobi integral. The equilibrium conditions for the system are derived and their stability analyzed. The equilibrium conditions of a celestial system comprised of n irregular bodies are reduced to 12n minus 9 equations. The dynamical results are applied to simulate the motion of multiple-asteroid systems. The simulation is useful for the study of the stability of multiple irregular celestial body systems and for the design of spacecraft orbits to triple asteroid systems discovered in the solar system. The dynamical configurations of the five triple-asteroid systems 45 Eugenia, 87 Sylvia, 93 Minerva, 216 Kleopatra, and 136617 1994CC, and the six-body system 134340 Pluto are calculated and analyzed.

We study the motion and equilibria of the grains on the surface of the irregular celestial body (hereafter called irregular bodies). Motions for the grains on the smooth and unsmooth surfaces are discussed, respectively. The linearized equations of motion relative to a surface equilibrium point and its characteristic equations are presented. Considering the stick-slip effect, the damping forces and the spring forces for the grain are calculated, then the linearized equations of motion and the characteristic equations relative to the surface equilibrium points are derived. The number of non-degenerate surface equilibria is an even number. We compute the motion of a grain released above three different regions relative to the irregular asteroid 6489 Golevka, including the flat surface, the concave region, and the convex region. Following the grain release and initial bounce, three kinds of motions exist, the orbital motion, the impact motion and the surface motion. We find that the maximum height of the next hop may be bigger than the maximum height of the current hopping. We also used Monte Carlo simulations to calculate 100 grains hopping motions, the results shows that the stable surface equilibria are on the concave region and flat surface of the asteroid.

We investigate the multiple bifurcations in periodic orbit families in the potential field of a highly irregular-shaped celestial body. Topological cases of periodic orbits and four kinds of basic bifurcations in periodic orbit families are studied. Multiple bifurcations in periodic orbit families consist of four kinds of basic bifurcations. We found both binary period-doubling bifurcations and binary tangent bifurcations in periodic orbit families around asteroid 433 Eros. The periodic orbit family with binary period-doubling bifurcations is nearly circular, with almost zero inclination, and is reversed relative to the body of the asteroid 433 Eros. This implies that there are two stable regions separated by one unstable region for the motion around this asteroid. In addition, we found triple bifurcations which consist of two real saddle bifurcations and one period-doubling bifurcation. A periodic orbit family generated from an equilibrium point of asteroid 433 Eros has five bifurcations, which are one real saddle bifurcation, two tangent bifurcations, and two period-doubling bifurcations.

Recently it has been suggested that the "twin-CME" scenario Li.etal2012 may be a very effective mechanism in causing extreme Solar Energetic Particle (SEP) events and in particular Ground Level Enhancement (GLE) events. Ding.etal2013 performed a statistical examination of the twin-CME scenario with a total of $126$ fast and wide western Coronal Mass Ejections (CMEs). They found that CMEs having a preceding CME with a speed $>$ 300 $km/s$ within $9$ hours from the same active region have larger probability of leading to large SEP events than CMEs that do not have preceding CMEs. The choice of $9$ hours being the time lag $\tau$ between the preceding CME and the main CME was based on some crude estimates of the decay time of the turbulence downstream of the shock driven by the preceding CME. In this work, we examine this choice. For the $126$ fast wide CMEs examined in Ding.etal2013, we vary the time lag $\tau$ from $1$ hour to $24$ hours with an increment of $1$ hour. By considering three quantities whose values depend on the choice of this time lag $\tau$, we show that the choice of $13$ hours for $\tau$ is more appropriate. Our study confirms our earlier result that twin CMEs are more likely to lead to large SEP events than single fast CMEs. The results shown here are of great relevance to space weather studies.

Some asteroids contain water ice, and a space mission landing on an asteroid may take liquid to the surface of the asteroid. Gas pressure is very weak on the surface of asteroids. Here we consider the capillary action in a crack on the surface of irregular asteroids. The crack is modelled as a capillary which has a fixed radius. An asteroid s irregular gravitational potential influences the height of the liquid in the capillary. The height of the liquid in the capillary on the surface of such asteroids is derived from the asteroid s irregular gravitational potential. Capillary mechanisms are expected to produce an inhomogeneaous distribution of emergent liquid on the surface. This result is applied to asteroid 433 Eros, which has an irregular, elongated, and concave shape. Two cases are considered 1) we calculate the height of the liquid in the capillary when the direction of the capillary is perpendicular to the local surface of the asteroid; 2) we calculate the height of the liquid in the capillary when the direction of the capillary is parallel to the vector from the center of mass to the surface position. The projected height in the capillary on the local surface of the asteroid seems to depend on the assumed direction of the capillary.

Stream-stream collisions play an important role for the circularization of highly eccentric streams resulting from tidal disruption events (TDEs). We perform three dimensional radiation hydrodynamic simulations to show that stream collisions can contribute significant optical and ultraviolet light to the flares produced by TDEs, and can sometimes explain the majority of the observed emission. Our simulations focus on the region near the radiation pressure dominated shock produced by a collision and track how the kinetic energy of the stream is dissipated by the associated shock. When the mass flow rate of the stream $\dot{M}$ is a significant fraction of the Eddington accretion rate, $\gtrsim2\%$ of the initial kinetic energy is converted to radiation directly as a result of the collision. In this regime, the collision redistributes the specific kinetic energy into the downstream gas and more than $16\%$ of the mass can become unbound. The fraction of unbound gas decreases rapidly as $\dot{M}$ drops significantly below the Eddington limit, with no unbound gas being produced when $\dot{M}$ drops to $1\%$ of Eddington; we find however that the radiative efficiency increases slightly to $\lesssim 8\%$ in these low $\dot{M}$ cases. The effective radiation temperature and size of the photosphere is determined by the stream velocity and $\dot{M}$, which we find to be a few times $10^4$~K and $10^{14}$~cm in our calculations, comparable to the inferred values of some TDE candidates. The photosphere size is directly proportional to $\dot{M}$, which can explain the rapidly changing photosphere sizes seen in TDE candidates such as PS1-10jh.

We have assembled the historical variability data of Mrk 421 at radio 15 GHz, X-ray and gamma-ray bands, spanning about 6.3, 10.3 and 7.5 yr, respective- ly. We analyzed the variability by using three methods. The results indicated that there is a period of 287.6 +- 4.4 days for 15 GHz, 309.5 +- 5.8 days for X- ray and 283.4+- 4.7 days for gamma-ray, respectively. This period can be reasonably explained by the nonballistic helical motion of the emitting material. The cor- relation analysis suggested that the variabilities of radio 15 GHz, X-ray and 0-ray are remarkable correlated, and the emission of radio 15 GHz lags behind that of X-ray, and the X-ray flux lags behind the gamma-ray. This suggests that the gamma-ray derives from inverse Compton (IC) scattering of the synchrotron photons, supporting the synchrotron self-Compton (SSC) model. Moreover, the time delay between different wavebands could be explained by the shock-in-jet models, in which a moving emission region produces the radio to gamma-ray activity, implying that the emission region of gamma-ray is closer to the center than ones of X-ray and radio emission.

Accretion disks around supermassive black holes have regions where the Rosseland mean opacity can be much larger than the electron scattering opacity primarily due to the large number of bound-bound transitions in iron. We study the effects of this iron opacity "bump" on the thermal stability and vertical structure of radiation pressure dominated accretion disks, utilizing three dimensional radiation magneto-hydrodynamic simulations in the local shearing box approximation. The simulations self-consistently calculate the heating due to MHD turbulence caused by magneto-rotational instability and radiative cooling by using the radiative transfer module based on a variable Eddington tensor in Athena. For a $5\times 10^8$ solar mass black hole with $\sim 3\%$ of the Eddington luminosity, a model including the iron opacity bump maintains its structure for more than $10$ thermal times without showing significant signs of thermal runaway. In contrast, if only electron scattering and free-free opacity are included as in the standard thin disk model, the disk collapses on the thermal time scale. The difference is caused by a combination of 1) an anti-correlation between the total optical depth and the temperature, and 2) enhanced vertical advective energy transport. These results suggest that the iron opacity bump may have a strong impact on the stability and structure of AGN accretion disks, and may contribute to a dependence of AGN properties on metallicity. Since this opacity is relevant primarily in UV emitting regions of the flow, it may help to explain discrepancies between observation and theory that are unique to AGNs.

In this paper, new determined BVR_cI_c light curves and radial velocities of V781 Tau are presented. By analyzing the light curves and radial velocities simultaneously, we found that V781 Tau is a W-subtype medium contact binary star with a mass ratio of q=2.207+-0.005 and a contact degree of f=21.6(+-1.0)%. The difference between the two light maxima was explained by a dark spot on the less massive primary component. The orbital period change of V781 Tau was also investigated. A secular decrease at a rate of $-6.01(+-2.28)*10^{-8} d/yr and a cyclic modulation with a period of 44.8+-5.7 yr and an amplitude of 0.0064+-0.0011 day were discovered. The continuous period decrease may be caused by angular momentum loss due to magnetic stellar wind. Applegate mechanism failed to explain the cyclic modulation. It is highly possible that the cyclic oscillation is the result of the light travel time effect by a third companion.

We are interested in the periodic motion and bifurcations near the surface of an asteroid. The gravity field of an irregular asteroid and the equation of motion of a particle near the surface of an asteroid are studied. The periodic motions around the major body of triple asteroid 216 Kleopatra and the OSIRIS REx mission target asteroid 101955 Bennu are discussed. We find that motion near the surface of an irregular asteroid is quite different from the motion near the surface of a homoplastically spheroidal celestial body. The periodic motions around the asteroid 101955 Bennu and 216 Kleopatra indicate that the geometrical shapes of the orbits are probably very sophisticated. There exist both stable periodic motions and unstable periodic motions near the surface of the same irregular asteroid. This periodic motion which is unstable can be resonant or non resonant. The period doubling bifurcation and pseudo period doubling bifurcation of periodic orbits coexist in the same gravity field of the primary of the triple asteroid 216 Kleopatra. It is found that both of the period doubling bifurcations of periodic orbits and pseudo period-doubling bifurcation of periodic orbits have four different paths. The pseudo period doubling bifurcation found in the potential field of primary of triple asteroid 216 Kleopatra shows that there exist stable periodic orbits near the primary s equatorial plane, which gives an explanation for the motion stability of the triple asteroid 216 Kleopatra s two moonlets, Alexhelios and Cleoselene.

On 13 December 2012, Chang'e-2 completed a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 meters from the asteroid's surface. The observations show that Toutatis has an irregular surface and its shape resembles a ginger-root of a smaller lobe (head) and a larger lobe (body). Such bilobate shape is indicative of a contact binary origin for Toutatis. In addition, the high-resolution images better than 3 meters provide a number of new discoveries about this asteroid, such as an 800-meter depression at the end of the large lobe, a sharply perpendicular silhouette near the neck region, boulders, indicating that Toutatis is probably a rubble-pile asteroid. Chang'e-2 observations have significantly revealed new insights into the geological features and the formation and evolution of this asteroid. In final, we brief the future Chinese asteroid mission concept.

Boulders are ubiquitously found on the surfaces of small rocky bodies in the inner solar system and their spatial and size distributions give insight into the geological evolution and collisional history of the parent bodies. Using images acquired by the Chang'e-2 spacecraft, more than 200 boulders have been identified over the imaged area of the near-Earth asteroid Toutatis. The cumulative boulder size frequency distribution (SFD) shows a steep slope of -4.4 $\pm$ 0.1, which is indicative of a high degree of fragmentation. Similar to Itokawa, Toutatis probably has a rubble-pile structure, as most boulders on its surface cannot solely be explained by impact cratering. The significantly steeper slope for Toutatis' boulder SFD compared to Itokawa may imply a different preservation state or diverse formation scenarios. In addition, the cumulative crater SFD has been used to estimate a surface crater retention age of approximately 1.6 $\pm$ 0.3 Gyr.

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $\tau_0$, to $\tau_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $\tau_0 \gg \tau_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $\tau_0 \ll \tau_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

The order and chaos of the motion near equilibrium points in the potential of a rotating highly irregular-shaped celestial body are investigated from point of view of the dynamical system theory. The positions of the non-degenerate equilibrium points vary continuously when the parameter changes. The topological structures in the vicinity of equilibrium points are classified into several different cases. Bifurcations at equilibrium points and the topological transfers between different cases for equilibrium points are also discussed. The conclusions can be applied to all kinds of rotating celestial bodies, simple-shaped or highly irregular-shaped, including asteroids, comets, planets and satellites of planets to help one to understand the dynamical behaviors around them. Applications to asteroids 216 Kleopatra, 2063 Bacchus, and 25143 Itokawa are significant and interesting: eigenvalues affiliated to the equilibrium points for the asteroid 216 Kleopatra move and always belong to the same topological cases; while eigenvalues affiliated to two different equilibrium points for the asteroid 2063 Bacchus and 25143 Itokawa move through the resonant cases of equilibrium points, and the collision of eigenvalues in the complex plane occurs. Poincar\'e sections in the potential of the asteroid 216 Kleopatra show the chaos behaviors of the orbits in large scale.

In this work, we investigate the bifurcations of relative equilibria in the gravitational potential of asteroids. A theorem concerning a conserved quantity, which is about the eigenvalues and number of relative equilibria, is presented and proved. The conserved quantity can restrict the number of non-degenerate equilibria in the gravitational potential of an asteroid. It is concluded that the number of non-degenerate equilibria in the gravitational field of an asteroid varies in pairs and is an odd number. In addition, the conserved quantity can also restrict the kinds of bifurcations of relative equilibria in the gravitational potential of an asteroid when the parameter varies. Furthermore, studies have shown that there exist transcritical bifurcations, quasi-transcritical bifurcations, saddle-node bifurcations, saddle-saddle bifurcations, binary saddle-node bifurcations, supercritical pitchfork bifurcations, and subcritical pitchfork bifurcations for the relative equilibria in the gravitational potential of asteroids. It is found that for the asteroid 216 Kleopatra, when the rotation period varies as a parameter, the number of relative equilibria changes from 7 to 5 to 3 to 1, and the bifurcations for the relative equilibria are saddle-node bifurcations and saddle-saddle bifurcations.

The eruption of the filament with the kink fashion is often regarded as a signature of the kink instability. However, the kink instability threshold for the filament magnetic structure has been not widely understood. Using the H-alpha observation from the New Vacuum Solar Telescope (NVST), we present a partial eruptive filament. In the eruption, a filament thread appeared to split from the middle portion of the filament and to break out in a kinklike fashion. During this period, the left filament material remained below, which erupted without the kinking motion later on. The coronal magnetic field lines associated with the filament are obtained from the nonlinear force-free field (NLFFF) extrapolations using the 12 minutes cadence vector magnetograms of the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). We studied the extrapolated field lines passing through the magnetic dips that are in good agreement with the observed filament. The field lines are non-uniformly twisted and appear to be made up by two twisted flux ropes winding about each other. One of them has higher twist than the other, and the highly twisted one has its dips aligned with the kinking eruptive thread at the beginning of its eruption. Before the eruption, moreover, the highly twisted flux rope was found to expand with the approximately constant field twist. In addition, the helicity flux maps deduced from the HMI magnetograms show that some helicity is injected into the overlying magnetic arcade, but no significant helicity is injected into the flux ropes. Accordingly, we suggest that the highly twisted flux rope became kink unstable when the instability threshold declined with the expansion of the flux rope.

We compiled the radio, optical, and X-ray data of blazars from the Sloan Digital Sky Survey (SDSS) database, and presented the distribution of luminosities and broad band spectral indices. The distribution of luminosities shows that the averaged luminosity of flat-spectral radio quasars (FSRQs) is larger than that of BL Lacs objects. On the other hand, the broad band spectral energy distribution reveals that FSRQs and low energy peaked BL Lac objects (LBLs) objects have similar spectral properties, but high energy peaked BL Lac objects (HBLs) have a distinct spectral property. This may be due to that different subclasses of blazars have different intrinsic environments and are at different cooling levels. Even so, a unified scheme also is revealed from the color-color diagram, which hints that there are similar physical processes operating in all objects under a range of intrinsic physical conditions or beaming parameter.

The equilibrium points of the gravitational potential field of minor celestial bodies, including asteroids, comets, and irregular satellites of planets, are studied. In order to understand better the orbital dynamics of massless particles moving near celestial minor bodies and their internal structure, both internal and external equilibrium points of the potential field of the body are analyzed. In this paper, the location and stability of the equilibrium points of 23 minor celestial bodies are presented. In addition, the contour plots of the gravitational effective potential of these minor bodies are used to point out the differences between them. Furthermore, stability and topological classifications of equilibrium points are discussed, which clearly illustrate the topological structure near the equilibrium points and help to have an insight into the orbital dynamics around the irregular-shaped minor celestial bodies. The results obtained here show that there is at least one equilibrium point in the potential field of a minor celestial body, and the number of equilibrium points could be one, five, seven, and nine, which are all odd integers. It is found that for some irregular-shaped celestial bodies, there are more than four equilibrium points outside the bodies while for some others there are no external equilibrium points. If a celestial body has one equilibrium point inside the body, this one is more likely linearly stable.

In this study, equilibrium points and periodic orbits in the potential field of asteroids are investigated. We present the linearized equations of motion relative to the equilibrium points and characteristic equations. We find that the distribution of characteristic multipliers of periodic orbits around the equilibrium point and the distribution of eigenvalues of the equilibrium point correspond to each other. The distribution of eigenvalues of the equilibrium point confirms the topology and the stability of periodic orbits around the equilibrium point.

This paper studies the distribution of characteristic multipliers, the structure of submanifolds, the phase diagram, bifurcations and chaotic motions in the potential field of rotating highly irregular-shaped celestial bodies (hereafter called irregular bodies). The topological structure of the submanifolds for the orbits in the potential field of an irregular body is shown to be classified into 34 different cases, including 6 ordinary cases, 3 collisional cases, 3 degenerate real saddle cases, 7 periodic cases, 7 period-doubling cases, 1 periodic and collisional case, 1 periodic and degenerate real saddle case, 1 period-doubling and collisional case, 1 period-doubling and degenerate real saddle case, and 4 periodic and period-doubling cases. The different distribution of the characteristic multipliers has been shown to fix the structure of the submanifolds, the type of orbits, the dynamical behaviour and the phase diagram of the motion. Classifications and properties for each case are presented. Moreover, tangent bifurcations, period-doubling bifurcations, Neimark-Sacker bifurcations and the real saddle bifurcations of periodic orbits in the potential field of an irregular body are discovered. Submanifolds appear to be Mobius strips and Klein bottles when the period-doubling bifurcation occurs.

Under the steady state condition, the spectrum of electrons is investigated by solving the continuity equation under the complex radiation of both the synchrotron and Compton processes. The resulted GRB spectrum is a broken power law in both the fast and slow cooling phases. On the basis of this electron spectrum, the spectral indices of the Band function in four different phases are presented. In the complex radiation frame, the detail investigation on physical parameters reveals that both the reverse shock photosphere model and the forward shock with strong coupling model can answer the $\alpha \sim -1$ problem. A possible marginal to fast cooling phase transition in GRB 080916C is discussed. The time resolved spectra in different pulses of GRB 100724B, GRB 100826A and GRB 130606B are investigated. We found that the flux is proportional to the peak energy in almost all pulses. The phases for different pulses are determined according to the spectral index revolution. We found the strong correlations between spectral indices and the peak energy in GRB 100826A, which can be explained by the Compton effect in the fast cooling phase. However, the complex scenario predicts a steeper index for the injected electrons, which challenges the acceleration mechanism in GRBs.

The polarization measurement is an important tool to probe the prompt emission mechanism in gamma-ray bursts (GRBs). The synchrotron photons can be scattered by cold electrons in the outflow via Compton scattering processes. The observed polarization depends on both the photon energy and the viewing angle. With the typical bulk Lorentz factor $\Gamma \sim 200$, photons with energy $E>10$ MeV tend to have smaller polarization than photons with energy $E<1$ MeV. At the right viewing angle, i.e. $\theta \sim \Gamma^{-1}$, the polarization achieves its maximal value, and the polarization angle changes $90^{\circ}$ relative to the initial polarization direction. Thus, the synchrotron radiation plus Compton scattering model can naturally explain the $90^{\circ}$ change of the polarization angle in GRB 100826A.

We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ~220L_edd/c^2 and forms a radiation driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ~20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ~10L_edd. This yields a radiative efficiency ~4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

When observing the blazars, we identified a new eclipsing binary GSC 03517-00663. The light curves of GSC 03517-00663 are typical EW-type light curves. Based on the observation using the 1m telescope at Weihai Observatory of Shandong University, complete VRI light curves were determined. Then, we analyzed the multiple light curves using the W-D program. It is found that GSC 03517-00663 has a mass ratio of q=0.164 and a contact degree of f=69.2%. GSC 03517-00663 is a deep, low mass ratio overcontact binary. The light curves of GSC 03517-00663 show strong O'Connell effect, it was explained by employing a dark spot on the secondary component.

The first photometric analysis of IR Cas was carried out based on the new observed BVRI light curves. The symmetric light curves and nearly flat secondary minimum indicate that very precise photometric results can be determined. We found that IR Cas is a near contact binary with the primary component filling its Roche lobe. An analysis of the O-C diagram based on all available times of light minimum reveals evidence for a periodic change with a semiamplitude of 0.0153 days and a period of 39.7 years superimposed on a secular decrease at a rate of dp/dt=-1.28(\pm0.09)\times10^{-7} d yr^{-1}. The most reasonable explanation for the periodic change is the light time-travel effect due to a third body. The period decrease may be caused by mass transfer from the primary component to the secondary. With the decreasing period, IR Cas would eventually evolve into a contact system.

When observing the transiting extrasolar planets, we discovered a new W UMa type binary GSC 03553-00845. Following observation was carried out using the 1m telescope at Weihai Observatory of Shandong University. Complete BVR light curves were determined. Using the W-D program, we analyzed the light curves. Photometric solution reveals that GSC 03553-00845 is a W-subtype W UMa type binary with a mass ratio of q=2.904, it is an overcontact binary system by a contact degree of f=29.5% with a small temperature difference between the components (\Delta T=206 K indicating a good thermal contact between the components. More observation of GSC 03553-00845 is needed in order to analyze the light curve variation and orbital period change.

(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the Athena code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers. The method is ideally suited for problems in which characteristic velocities are non-relativistic, but still within a few percent or more of the speed of light. The method is an intermediate step towards algorithms for fully relativistic flows.

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

This study investigates the different novel forms of the dynamical equations of a particle orbiting a rotating asteroid and the effective potential, the Jacobi integral, etc. on different manifolds. Nine new forms of the dynamical equations of a particle orbiting a rotating asteroid are presented, and the classical form of the dynamical equations has also been found. The dynamical equations with the potential and the effective potential in scalar form in the arbitrary body-fixed frame and the special body-fixed frame are presented and discussed. Moreover, the simplified forms of the effective potential and the Jacobi integral have been derived. The dynamical equation in coefficient-matrix form has been derived. Other forms of the dynamical equations near the asteroid are presented and discussed, including the Lagrange form, the Hamilton form, the symplectic form, the Poisson form, the Poisson-bracket form, the cohomology form, and the dynamical equations on the K\"ahler manifold and another complex manifold. Novel forms of the effective potential and the Jacobi integral are also presented. The dynamical equations in scalar form and coefficient-matrix form can aid in the study of the dynamical system, the bifurcation, and the chaotic motion of the orbital dynamics of a particle near a rotating asteroid. The dynamical equations of a particle near a rotating asteroid are presented on several manifolds, including the symplectic manifold, the Poisson manifold, and complex manifolds, which may lead to novel methods of studying the motion of a particle in the potential field of a rotating asteroid.

We study the orbits and manifolds near the equilibrium points of a rotating asteroid. The linearised equations of motion relative to the equilibrium points in the gravitational field of a rotating asteroid, the characteristic equation and the stable conditions of the equilibrium points are derived and discussed. First, a new metric is presented to link the orbit and the geodesic of the smooth manifold. Then, using the eigenvalues of the characteristic equation, the equilibrium points are classified into 8 cases. A theorem is presented and proved to describe the structure of the submanifold as well as the stable and unstable behaviours of a massless particle near the equilibrium points. The linearly stable, the non-resonant unstable, and the resonant equilibrium points are discussed. There are three families of periodic orbits and four families of quasi-periodic orbits near the linearly stable equilibrium point. For the non-resonant unstable equilibrium points, there are four cases; for the periodic orbit and the quasi-periodic orbit, the structures of the submanifold and the subspace near the equilibrium points are studied for each case. For the resonant equilibrium points, the dimension of the resonant manifold is greater than four, and we find at least 1 family of periodic orbits near the resonant equilibrium points. Besides, this theory is applied to asteroids 216 Kleopatra, 1620 Geographos, 4769 Castalia, and 6489 Golevka.

Synchrotron radiation and Compton scattering are widely accepted as the most likely emission mechanisms of some astrophysical phenomena, such as gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). The measurement on polarization of photons provides a useful tool to distinguish different emission mechanisms and structures of the emission region. Based on the differential cross section of a polarized photon scattered by an unpolarized electron of any initial momentum, we derive analytical formula of polarization for beamed photons scattered by isotropic electrons with a power law distribution. Numerical calculations are carried out in four special cases: electrons at rest, Thomson limit, head-on collision and monochromatic electrons. It is found that the maximum polarization can be as high as $100\%$ for low energy photons, if the electrons are at rest. Although polarization is highly suppressed due to the isotropic electrons, a maximum value of $\sim 10\% \-- 20\%$ can still be achieved. Compton scattering process can be used to explain the polarization of GRB 041219A and GRB 100826A.

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magneto-hydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the alpha-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from AGNs, and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

Based on the quantum electrodynamics, we present a generic formalism of the polarization for beamed monochromatic photons scattered by electrons in any spectral distribution. The formulae reduce to the components of the Fano matrix when electrons are at rest. We mainly investigate the polarization in three scenarios, i.e., electrons at rest, isotropic electrons with a power law spectrum and thermal electrons. If the incident beam is polarized, the polarization is reduced significantly by isotropic electrons at large viewing angles, and the degree of polarization due to thermal electrons is about one times less than that of electrons in a power law. If the incident bean is unpolarized, soft $\gamma$-rays can lead to about 15% polarization at viewing angles around $\pi/4$. For isotropic electrons, one remarkable feature is that the polarization as a function of the incident photon energy always peaks roughly at 1 MeV, this is valid for both the thermal and power law cases. This feature can be used to distinguish the model of the inverse Compton scattering from that of the synchrotron radiation.

We report on the first phase of our study of cloud irradiation. We study irradiation by means of numerical, two-dimensional time-dependent radiation-hydrodynamic simulations of a cloud irradiated by a strong radiation. We adopt a very simple treatment of the opacity, neglect photoionization and gravity, and instead focus on assessing the role of the type and magnitude of the opacity on the cloud evolution. Our main result is that even relatively dense clouds that are radiatively heated (i.e., with significant absorption opacity) do not move as a whole instead they undergo a very rapid and major evolution in its shape, size and physical properties. In particular, the cloud and its remnants become optical thin within less than one sound crossing time and before they can travel over a significant distance (a distance of a few radii of the initial cloud). We also found that a cloud can be accelerated as a whole under quite extreme conditions, e.g., the opacity must be dominated by scattering. However, the acceleration due to the radiation force is relatively small and unless the cloud is optically thin the cloud quickly changes its size and shape. We discuss implications for the modelling and interpetation broad line regions of active galactic nuclei.

Clusters of galaxies have a huge mass which can act as gravitational lenses. Galaxies behind clusters can be distorted to form arcs in images by the lenses. Herein a search was done for giant lensed arcs by galaxy clusters using the SDSS data. By visually inspecting SDSS images of newly identified clusters in the SDSS DR8 and Stripe 82 data, we discover 8 strong lensing clusters together with additional 3 probable and 6 possible cases. The lensed arcs show bluer colors than the member galaxies of clusters. The masses and optical luminosities of galaxy clusters interior to the arcs are calculated, and the mass-to-light ratios are found to be in the range of a few tens of M_Solar/L_Solar, consistent with the distribution of previously known lensing clusters.

We study the long-term thermal stability of radiation dominated disks in which the vertical structure is determined self-consistently by the balance of heating due to dissipation of MHD turbulence driven by the magneto-rotational instability (MRI), and cooling due to radiation emitted at the photosphere. The calculations adopt the local shearing box approximation, and utilize the recently developed radiation transfer module in the Athena MHD code based on a variable Eddington tensor rather than an assumed local closure. After saturation of the MRI, in many cases the disk maintains a steady vertical structure for many thermal times. However, in every case in which the box size in the horizontal directions is at least one pressure scale height, fluctuations associated with MRI turbulence and dynamo action in the disk eventually trigger a thermal runaway which causes the disk to either expand or contract until the calculation must be terminated. During runaway, the dependence of the heating and cooling rates on total pressure satisfy the simplest criterion for classical thermal instability. We identify several physical reasons why the thermal runaway observed in our simulations differ from the standard alpha disk model, for example the advection of radiation contributes a non-negligible fraction to the vertical energy flux at the largest radiation pressure, most of the dissipation does not happen in the disk mid-plane, and the change of dissipation scale height with mid-plane pressure is slower than the change of density scale height. We discuss how and why our results differ from those published previously. Such thermal runaway behavior might have important implications for interpreting temporal variability in observed systems, but fully global simulations are required to study the saturated state before detailed predictions can be made.

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches 1% of the rest mass energy density of the gas), we find Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of 8 compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find radiative shear viscosity is negligible compared to the Maxwell and Reynolds stress in the flow. This may have important implications for the structure of radiation dominated accretion disks.

Possible Lorentz invariance violation (LIV) has been investigated for a long time based on observations of GRBs . These arguments relied on the assumption that photons with different energy are emitted at the same place and time. In this work, we try to take account of the intrinsic time delay $\Delta t_{\rm int}$ between emissions of low and high energy photons by using the magnetic jet model. The possible LIV effects are discussed in a unified scenario both for long and short {\it Fermi}-detected GRBs. This leads to a unique quantum gravity energy scale $M_1c^2 \sim 1.0 \times 10^{20}$ GeV respecting the linear dispersion relation.

Recent observations by the Fermi-LAT showed that there are delayed arrivals of GeV photons relative to the onset of MeV photons in some GRBs. In order to avoid a large optical depth, the minimal value of the Lorentz factor has been estimated to be higher than 1000 in some brightest bursts. In this paper, we present a detailed calculation of the time delay between the MeV and GeV photons in the framework of the magnetic-dominated jet model. We find that the time delay strongly depends on the saturated bulk Lorentz factor of the jet. Inspired by this fact, we use this model to calculate the Lorentz factors of the four brightest Fermi bursts. The results indicate that the Lorentz factors are much smaller than that obtained from the "single-zone" scenario. The short GRB 090510 has a minimal Lorentz factor 385, while the three long bursts GRB 080916c, GRB090902b and GRB 090926 have almost the same Lorentz factors, with an average value near 260. Another interesting result is that, for long bursts, GeV photons are emitted after the bulk Lorentz factor saturates. For the short GRB, however, MeV and GeV photons are emitted at the same phase, i.e., either at the expansion phase or at the coast phase.

Light axion-like particles (ALPs) that couple to two photons are allowed in a number of proposed extensions to the Standard Model of elementary particles. Of particular interest from a theoretical and observational standpoint is the energy regime near 0.1 meV. We present results from a pilot experiment to search for a signal from a 0.14 meV scalar ALP by way of its coupling to two photons. Using a copper resonant cavity cooled to four degrees Kelvin while immersed in a seven Tesla magnetic field, and coupled to a low noise cryogenic amplifier and room temperature receiver, we exclude an ALP-driven excess of 34 GHz photons with g$>10^{-8}$/GeV with 5$\sigma$ confidence. We discuss the ramifications of this initial measurement as well as planned modifications to the experiment for increased sensitivity.

The non-linear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small scale structures are also suppressed in this case. In the non-linear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a variable Eddington tensor (VET) versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the non-linear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth rate of RTI is smaller than acceleration time scale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.

We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow.

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)

We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.

It is shown that the extremal Reissner-Nordstr\"{o}m black hole, the non-extremal one with multiple scattering particles, and the Schwarzschild black hole with radial head-on particles are stable under the collision of the particles near the horizon, if the back-reaction effect and the effect generated by gravity of particles are involved. Moreover, the collision near Reissner-Nordstr\"{o}m black holes with astrophysically typical mass can not generate the Planck-scale center-of-mass energy. However, the head-on collision near the typical primordial black hole could just occur at the Planck-energy scale.

We report on the rare eruption of a miniature H_alpha filament that took a surge form. The filament first underwent a full development within 46 minutes and then began to erupt 9 minutes later, followed by a compact, impulsive X-ray class M2.2 flare with a two ribbon nature only at the early eruption phase. During the eruption, its top rose, whereas the two legs remained rooted in the chromosphere and swelled little perpendicular to the rising direction. This led to a surge-like eruption with a narrow angular extent. Similar to the recent observations for standard and blowout X-ray jets by Moore et al., we thus define it as a "blowout H_alpha surge". Furthermore, our observations showed that the eruption was associated with (1) a coronal mass ejection guided by a preexisting streamer, (2) abrupt, significant, and persistent changes in the photospheric magnetic field around the filament, and (3) sudden disappearance of a small pore. These observations thus provide evidence that blowout surge is a small-scale version of large-scale filament eruption in many aspects. Our observations further suggest that at least part of H_alpha surges belong to blowout-type ones, and exact distinction between standard and blowout H_alpha surges is important in understanding their different origins and associated eruptive phenomena.

We identify new strong lensing clusters of galaxies from the Sloan Digital Sky Survey III (SDSS DR8) by visually inspecting color images of a large sample of clusters of galaxies. We find 68 new clusters showing giant arcs in addition to 30 known lensing systems. Among 68 cases, 13 clusters are "almost certain" lensing systems with tangential giant arcs, 22 clusters are "probable" and 31 clusters are "possible" lensing systems. We also find two exotic systems with blue rings. The giant arcs have angular separations of 2.0"-25.7" from the bright central galaxies. We note that the rich clusters are more likely to be lensing systems, and the separations between arcs and the central galaxies increase with cluster richness.

We study the possibility of getting infinite energy in the center of mass frame of colliding charged particles in a general stationary charged black hole. For black holes with two-fold degenerate horizon, it is found that arbitrary high center-of-mass energy can be attained, provided that one of the particle has critical angular momentum or critical charge, and the remained parameters of particles and black holes satisfy certain restriction. For black holes with multiple-fold degenerate event horizons, the restriction is released. For non-degenerate black holes, the ultra-high center-of-mass is possible to be reached by invoking the multiple scattering mechanism. We obtain a condition for the existence of innermost stable circular orbit with critical angular momentum or charge on any-fold degenerate horizons, which is essential to get ultra-high center-of-mass energy without fine-tuning problem. We also discuss the proper time spending by the particle to reach the horizon and the duality between frame dragging effect and electromagnetic interaction. Some of these general results are applied to braneworld small black holes.

We study the scaling between bulge magnitude and central black hole (BH) mass in galaxies with virial BH masses < 10^6 solar mass. Based on careful image decomposition of a snapshot Hubble Space Telescope I-band survey, we found that these BHs are found predominantly in galaxies with pseudobulges. Here we show that the \mbulge\ relation for the pseudobulges at low mass is significantly different from classical bulges with BH masses >10^7 solar mass. Specfically, bulges span a much wider range of bulge luminosity, and on average the luminosity is larger, at fixed black hole mass. The trend holds both for the active galaxies from Bentz et al. and the inactive sample of Gultekin et al. and cannot be explained by differences in stellar populations, as it persists when we use dynamical bulge masses. Put another way, the ratio between bulge and BH mass is much larger than $\sim 1000$ for our sample. This is consistent with recent suggestions that black hole mass does not scale with the pseudobulge luminosity. The low-mass scaling relations appear to flatten, consistent with predictions from Volonteri & Natarajan for massive seed BHs.

Using HST observations of 147 host galaxies of low-mass black holes (BHs), we systematically study the structures and scaling relations of these active galaxies. Our sample is selected to have central BHs with virial masses ~10^5-10^6 solar mass. The host galaxies have total I-band magnitudes of -23.2<M_I<-18.8 mag and bulge magnitudes of -22.9<M_I<-16.1 mag. Detailed bulge-disk-bar decompositions with GALFIT show that 93% of the galaxies have extended disks, 39% have bars and 5% have no bulges at all at the limits of our observations. Based on the Sersic index and bulge-to-total ratio, we conclude that the majority of the galaxies with disks are likely to contain pseudobulges and very few of these low-mass BHs live in classical bulges. The fundamental plane of our sample is offset from classical bulges and ellipticals in a way that is consistent with the scaling relations of pseudobulges. The sample has smaller velocity dispersion at fixed luminosity in the Faber-Jackson plane, compared with classical bulges and elliptical galaxies. The galaxies without disks are structurally more similar to spheroidals than to classical bulges according to their positions in the fundamental plane, especially the Faber-Jackson projection. Overall, we suggest that BHs with mass < 10^6 solar mass live in galaxies that have evolved secularly over the majority of their history. A classical bulge is not a prerequisite to host a black hole.

We present new measurements of stellar velocity dispersions, using spectra obtained with the Keck Echellette Spectrograph and Imager (ESI) and the Magellan Echellette (MagE), for 76 Seyfert 1 galaxies from the recent catalogue of Greene & Ho. These objects were selected from the Sloan Digital Sky Survey (SDSS) to have estimated black hole (BH) masses below 2\times10^6 M\odot. Combining our results with previous ESI observations of similar objects, we obtain an expanded sample of 93 galaxies and examine the relation between BH mass and velocity dispersion (the M-Sigma relation) for active galaxies with low BH masses. The low-mass active galaxies tend to follow the extrapolation of the M-Sigma relation of inactive galaxies. Including results for active galaxies of higher BH mass from the literature, we find a zero point{\alpha}= 7.68\pm0.08 and slope of {\beta}= 3.32\pm0.22 for the M-Sigma relation [in the form log Mbh={\alpha}+{\beta}log({\sigma}*/200 km s-1)], with intrinsic scatter of 0.46\pm0.03 dex. This result is consistent, within the uncertainties, with the slope of the M-Sigma relation for reverberation-mapped active galaxies with BH masses from 10^6 to 10^9 M\odot. For the subset of our sample having morphological information from Hubble Space Telescope images, we examine the slope of the M-Sigma relation separately for subsamples of barred and unbarred host galaxies, and find no significant evidence for a difference in slope. We do find a mild offset between low-inclination and high-inclination disk galaxies, such that more highly inclined galaxies tend to have larger {\sigma}* at a given value of BH mass, presumably due to the contribution of disk rotation within the spectroscopic aperture. We also find that the velocity dispersion of the ionized gas trace the stellar velocity dispersion well for this large sample of low-mass Seyfert 1 galaxies.

We present simultaneous observations of three recurring jets in EUV and soft X-ray (SXR), which occurred in an active region on 2007 June 5. By comparing their morphological and kinematic characteristics in these two different wavelengths, we found that EUV and SXR jets had similar locations, directions, sizes and velocities. We also analyzed their spectral properties by using six spectral lines from the EUV Imaging Spectrometer (EIS) on board Hinode, and found that these jets had temperatures from 0.05 to 2.0 MK and maximum electron densities from 6.6$\times10^{9}$ to 3.4$\times10^{10}$ cm$^{-3}$. For each jet, an elongated blue-shifted component and a red-shifted component at the jet base were simultaneously observed in Fe{\sc xii} $\lambda$195 and He{\sc ii} $\lambda$256 lines. The three jets had maximum Doppler velocities from 25 to 121 km s$^{-1}$ in Fe{\sc xii} $\lambda$195 line and from 115 to 232 km s $^{-1}$ in He{\sc ii} $\lambda$256 line. They had maximum non-thermal velocities from 98 to 181 km s$^{-1}$ in Fe{\sc xii} $\lambda$195 line and from 196 to 399 km s$^{-1}$ in He{\sc ii} $\lambda$256 line. We also examined the relationship between averaged Doppler velocities and maximum ionization temperatures of these three jets, and found that averaged Doppler velocities decreased with the increase of maximum ionization temperatures. In the photosphere, magnetic flux emergences and cancellations continuously took place at the jet base. These observational results were consistent with the magnetic reconnection jet model that magnetic reconnection between emerging magnetic flux and ambient magnetic field occurred in the lower atmosphere.

Using a version of the ZEUS code, we carry out two-dimensional simulations of self-gravitating shearing sheets, with application to QSO accretion disks at a few thousand Schwarzschild radii, corresponding to a few hundredths of a parsec for a 10^8 solar-mass black hole. Radiation pressure and optically thick radiative cooling are implemented via vertical averages. We determine dimensionless versions of the maximum surface density, accretion rate, and effective viscosity that can be sustained by density-wave turbulence without fragmentation. Where fragments do form, we study the final masses that result. The maximum Shakura-Sunyaev viscosity parameter is approximately 0.4. Fragmentation occurs when the cooling time is less than about twice the shearing time, as found by Gammie and others, but can also occur at very long cooling times in sheets that are strongly radiation-pressure dominated. For accretion at the Eddington rate onto a 10^8 solar-mass black hole, fragmentation occurs beyond four thousand Schwarzschild radii, r_s. Near this radius, initial fragment masses are several hundred suns, consistent with estimates from linear stability; final masses after merging increase with the size of the sheet, reaching several thousand suns in our largest simulations. With increasing black-hole mass at a fixed Eddington ratio, self-gravity prevails to smaller multiples of r_s, where radiation pressure is more important and the cooling time is longer compared to the dynamical time; nevertheless, fragmentation can occur and produces larger initial fragment masses. We also find energy conservation is likely to be a challenge for all eulerian codes in self-gravitating regimes where radiation pressure dominates.

We study possible correlations between ultrahigh energy cosmic rays (UHECRs), observed by Auger, AGASA and Yakutsk, and nearby active galactic nuclei (AGNs) and $Fermi$ sources. We consider the deflection effects by a Galactic magnetic field (GMF) model constrained by the most updated measurements. We found that the average deflection angles of UHECRs by the Galactic magnetic fields are less than $4^\circ$. A correlation between the Auger cosmic-ray events and nearby AGNs with a significance level of $\sim 4\sigma$ was found for the Auger UHECR data sets with or without deflection correction. No correlation was found between the AGASA/Yakutsk events with nearby AGNs. Marginal correlations between the Auger events and the $Fermi$ sources, and between AGASA events and $Fermi$ AGNs were found when the deflections calculated by the GMF model were considered. However, no correlation was found between the Yakutsk data and $Fermi$ sources. Some $Fermi$ sources are close to the arrival directions of UHECR events detected by Auger, AGASA, and Yakutsk, most of which are probably chance coincidence rather than objects producing UHECRs in the nearby Universe. Four $Fermi$ sources, NGC 4945, ESO 323-G77, NGC 6951, and Cen A, within 100~Mpc have UHECR events within $3.1^{\circ}$ from their positions, which could potentially be cosmic ray accelerators. However, the association can only be confirmed if more UHECRs are preferably detected in these directions.

Both radiative and mechanical feedback from Active Galactic Nuclei have been found to be important for the evolution of elliptical galaxies. We compute how a shock may be driven from a central black hole into the gaseous envelope of an elliptical galaxy by mechanical as well as radiative feedback (in the form of nuclear winds) using high resolution 1-D hydrodynamical simulations. We calculate the synchrotron emission from the electron cosmic rays accelerated by the shocks (not the jets) in the central part of elliptical galaxies, and we also study the synchrotron spectrum's evolution using the standard diffusive shock acceleration mechanism, which is routinely applied to the scaled volume case of supernova remnants. We find good agreement quantitatively between the synchrotron radio emission produced via this mechanism with extant observations of elliptical galaxies which are undergoing outbursts. Additionally, we also find that synchrotron optical and X-ray emission can co-exist inside elliptical galaxies during a certain phase of evolution subsequent to central outbursts. In fact, our calculations predict a synchrotron luminosity of $\sim 1.3\times 10^6 L_{\odot}$ at the frequency 5 GHz (radio band), of $\sim 1.1\times 10^6 L_{\odot}$ at $4.3\times10^{14}$ Hz (R band, corresponding to the absolute magnitude -10.4), and of $\sim 1.5\times 10^{7} L_{\odot}$ at $2.4\times10^{17}$ Hz (soft X-ray, 0.5 -- 2.0 keV band).

We study the orbital evolution of wide binary stars in the solar neighborhood due to gravitational perturbations from passing stars. We include the effects of the Galactic tidal field and continue to follow the stars after they become unbound. For a wide variety of initial semi-major axes and formation times, we find that the number density (stars per unit logarithmic interval in projected separation) exhibits a minimum at a few times the Jacobi radius r_J, which equals 1.7 pc for a binary of solar-mass stars. The density peak interior to this minimum arises from the primordial distribution of bound binaries, and the exterior density, which peaks at \sim 100--300 pc separation, arises from formerly bound binaries that are slowly drifting apart. The exterior peak gives rise to a significant long-range correlation in the positions and velocities of disk stars that should be detectable in large astrometric surveys such as GAIA that can measure accurate three-dimensional distances and velocities.

It is well-recognized that the presence of magnetic fields will lead to anisotropic energy cascade and dissipation of astrophysical turbulence. With the diffusion approximation and linear dissipation rates, we study the cascade and damping of Alfv\'en-cyclotron fluctuations in solar plasmas numerically. For an isotropic case the steady-state turbulence spectra are nearly isotropic in the inertial range and can be fitted by a single power-law function with a spectral index of -3/2, similar to the Iroshnikov-Kraichnan phenomenology. Beyond the MHD regime the kinetic effects make the spectrum softer at higher wavenumbers. In the dissipation range the turbulence spectrum cuts off at the wavenumber, where the damping rate becomes comparable to the cascade rate, and the cutoff wavenumber changes with the wave propagation direction. The angle averaged turbulence spectrum of the isotropic model resembles a broken power-law. Taking into account the Doppler effects, the model naturally reproduces the broken power-law turbulence spectra observed in the solar wind and predicts that a higher break frequency always comes along with a softer dissipation range spectrum that may be caused by the increase of the turbulence intensity, the reciprocal of the plasma \beta, and/or the angle between the solar wind velocity and the mean magnetic field. These predictions can be tested by detailed comparisons with more accurate observations.

We report the discovery of 4 strong gravitational lensing systems by visual inspections of the Sloan Digital Sky Survey images of galaxy clusters in Data Release 6 (SDSS DR6). Two of the four systems show Einstein rings while the others show tangential giant arcs. These arcs or rings have large angular separations (>8") from the bright central galaxies and show bluer color compared with the red cluster galaxies. In addition, we found 5 probable and 4 possible lenses by galaxy clusters.

By using H$\alpha$, He I 10830, EUV and soft X-ray (SXR) data, we examined a filament eruption that occurred on a quiet-sun region near the center of the solar disk on 2006 January 12, which disturbed a sigmoid overlying the filament channel observed by the $\emph{GOES-12}$ SXR Imager (SXI), and led to the eruption of the sigmoid. The event was associated with a partial halo coronal mass ejection (CME) observed by the Large Angle and Spectrometric Coronagraphs (LASCO) on board the Solar and Heliospheric Observatory ($\emph{SOHO}$), and resulted in the formation of two flare-like ribbons, post-eruption coronal loops, and two transient coronal holes (TCHs), but there were no significantly recorded $\emph{GOES}$ or H$\alpha$ flares corresponding to the eruption. The two TCHs were dominated by opposite magnetic polarities and were located on the two ends of the eruptive sigmoid. They showed similar locations and shapes in He I 10830, EUV and SXR observations. During the early eruption phase, brightenings first appeared on the locations of the two subsequent TCHs, which could be clearly identified on He I 10830, EUV and SXR images. This eruption event could be explained by the magnetic flux rope model, and the two TCHs were likely to be the feet of the flux rope.

Growing evidence indicate supermassive black holes (SMBHs) in the mass range of $M_{\rm BH}$$\sim 10^6-10^{10}M_{\odot}$ lurking in central bulges of many galaxies. Extensive observations reveal fairly tight power laws of $M_{\rm BH}$ versus the mean stellar velocity dispersion $\sigma$ of the host bulge. The dynamic evolution of a bulge and the formation of a central SMBH should be physically linked by various observational clues. In this contribution, we reproduce the empirical $M_{\rm BH}-\sigma$ power laws based on a self-similar general polytropic quasi-static bulge evolution and a sensible criterion of forming a SMBH surrounding the central density singularity of a general singular polytropic sphere (SPS) \cite{loujiang2008}. Other properties of host bulges and central SMBHs are also examined. Based on our model, we discuss the intrinsic scatter of the $M_{\rm BH}-\sigma$ relation and a scenario for the evolution of SMBHs in different host bulges.

Solar Energetic Particles (SEPs) show a rich variety of spectra and relative abundances of many ionic species and their isotopes. A long standing puzzle has been the extreme enrichments of 3He ions. The most extreme enrichments are observed in low fluence, the so-called impulsive, events which are believed to be produced at the flare site in the solar corona with little scattering and acceleration during transport to the Earth. In two earlier papers (Liu et al. 2004 and 2006) we showed how such extreme enrichments can result in the model developed by Petrosian and Liu (2004), where ions are accelerated stochastically by plasma waves or turbulence. In this paper we address the relative distributions of the fluences of 3He and 4He ions presented by Ho et al. (2005) which show that while the distribution of 4He fluence like many other extensive characteristics of solar flare, is fairly broad, the 3He fluence is limited to a narrow range. Moreover, the ratio of the fluences shows a strong correlation with the 4He fluence. One of the predictions of our model was presence of steep variation of the fluence ratio with the level of turbulence or the rate of acceleration. We show here that this feature of the model can reproduce the observed distribution of the fluences with very few free parameters. The primary reason for the success of the model in both fronts is because fully ionized 3He ion, with its unique charge to mass ratio, can resonantly interact with more plasma modes and accelerate more readily than 4He. Essentially in most flares, all background 3He ions are accelerated to few MeV/nucleon range, while this happens for 4He ions only in very strong events. A much smaller fraction of 4He ions reach such energies in weaker events.

Growing evidence indicate supermassive black holes (SMBHs) in a mass range of $M_{\rm BH}$$\sim 106-10^{10}M_{\odot}$ lurking in central stellar bulges of galaxies.Extensive observations reveal fairly tight power laws of $M_{\rm BH}$ versus the mean stellar velocity dispersion $\sigma$ of the host stellar bulge.Together with evidence for correlations between $M_{\rm BH}$ and other properties of host bulges, the dynamic evolution of a bulge and the formation of a central SMBH should be linked. In this Letter, we reproduce the empirical $M_{\rm BH}-\sigma$ power laws based on our recent theoretical analyses (Lou & Wang; Wang & Lou; Lou, Jiang & Jin) for a self-similar general polytropic quasi-static dynamic evolution of bulges with self-gravity and spherical symmetry and present a sensible criterion of forming a central SMBH. The key result is $M_{\rm BH}={\cal L}\sigma^{1/(1-n)}$ where $2/3<n<1$ and ${\cal L}$ is a proportional coefficient characteristic of different classes of host bulges. By fitting and comparing several empirical $M_{\rm BH}-\sigma$ power laws, we conclude that SMBHs and galactic bulges grow and evolve in a coeval manner and most likely there exist several classes of galactic bulge systems in quasi-static self-similar evolution and that to mix them together can lead to an unrealistic fitting. Based on our bulge-SMBH model, we provide explanations for intrinsic scatter in the relation and a unified scenario for the formation and evolution of SMBHs in different classes of host bulges.

With the diffusion approximation, we study the cascade and damping of Alfv\'{e}n-cyclotron fluctuations in solar plasmas numerically. Motivated by wave-wave couplings and nonlinear effects, we test several forms of the diffusion tensor. For a general locally anisotropic and inhomogeneous diffusion tensor in the wave vector space, the turbulence spectrum in the inertial range can be fitted with power-laws with the power-law index varying with the wave propagation direction. For several locally isotropic but inhomogeneous diffusion coefficients, the steady-state turbulence spectra are nearly isotropic in the absence of damping and can be fitted by a single power-law function. However, the energy flux is strongly polarized due to the inhomogeneity that leads to an anisotropic cascade. Including the anisotropic thermal damping, the turbulence spectrum cuts off at the wave numbers, where the damping rates become comparable to the cascade rates. The combined anisotropic effects of cascade and damping make this cutoff wave number dependent on the wave propagation direction, and the propagation direction integrated turbulence spectrum resembles a broken power-law, which cuts off at the maximum of the cutoff wave numbers or the $^4$He cyclotron frequency. Taking into account the Doppler effects, the model can naturally reproduce the broken power-law wave spectra observed in the solar wind and predicts that a higher break frequency is aways accompanied with a greater spectral index change that may be caused by the increase of the Alfv\'{e}n Mach number, the reciprocal of the plasma beta, and/or the angle between the solar wind velocity and the mean magnetic field. These predictions can be tested by future observations.

A theoretical model framework of spherical symmetry is presented for a composite astrophysical system of two polytropic fluids coupled together by gravity to explore large-scale shocks and flow dynamics in clusters of galaxies or in globular clusters. The existence of such large-scale shocks in clusters of galaxies as inferred by high-resolution X-ray and radio imaging observations implies large-scale systematic flows that are beyond usual static models for clusters of galaxies. Here, we explore self-similar two-fluid flow solutions with shocks for a hot polytropic gas flow in a cluster of galaxies in the presence of a massive dark matter (DM) flow after the initiation of a gravitational core collapse or a central AGN activity or a large-scale merging process. In particular, the possibility of DM shocks or sharp jumps of mass density and of velocity dispersion in dark matter halo is discussed and such DM shocks might be detectable through gravitational lensing effects. To examine various plausible scenarios for clusters of galaxies, we describe three possible classes of shock flows within our model framework for different types of temperature, density and flow speed profiles. Depending upon sensible model parameters and shock locations, the hot ICM and DM halo may have various combinations of asymptotic behaviours of outflow, breeze, inflow, contraction or static envelopes at large radii at a given time. We refer to asymptotic outflows of hot ICM at large radii as the galaxy cluster wind. As a result of such galaxy cluster winds and simultaneous contractions of DM halo during the course of galaxy cluster evolution, there would be less hot ICM within clusters of galaxies as compared to the average baryon fraction in the Universe.

We present data analysis and interpretation of an M1.4-class flare observed with the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) on April 30, 2002. This event, with its footpoints occulted by the solar limb, exhibits a rarely observed, but theoretically expected, double-source structure in the corona. The two coronal sources, observed over the 6-30 keV range, appear at different altitudes and show energy-dependent structures with the higher-energy emission being closer together. Spectral analysis implies that the emission at higher energies in the inner region between the two sources is mainly nonthermal, while the emission at lower energies in the outer region is primarily thermal. The two sources are both visible for about 12 minutes and have similar light curves and power-law spectra above about 20 keV. These observations suggest that the magnetic reconnection site lies between the two sources. Bi-directional outflows of the released energy in the form of turbulence and/or particles from the reconnection site can be the source of the observed radiation. The spatially resolved thermal emission below about 15 keV, on the other hand, indicates that the lower source has a larger emission measure but a lower temperature than the upper source. This is likely the result of the differences in the magnetic field and plasma density of the two sources.

We present analyses of the spatial and spectral evolution of hard X-ray emission observed by {\it RHESSI} during the impulsive phase of an M1.7 flare on 2003 November 13. In general, as expected, the loop top (LT) source dominates at low energies while the footpoint (FP) sources dominate the high energy emission. At intermediate energies, both the LT and FPs may be seen, but during certain intervals emission from the legs of the loop dominates, in contrast to the commonly observed LT and FP emission. The hard X-ray emission tends to rise above the FPs and eventually merge into a single LT source. This evolution starts first at low energies and proceeds to higher energies. The spectrum of the resultant LT source becomes more and more dominated by a thermal component with an increasing emission measure as the flare proceeds. The soft and hard X-rays show a Neupert-type behavior. With a non-thermal bremsstrahlung model the brightness profile along the loop is used to determine the density profile and its evolution, which reveals a gradual increase of the gas density in the loop. These results are evidence for chromospheric evaporation and are consistent with the qualitative features of hydrodynamic simulations of this phenomenon. However, some observed source morphology and its evolution cannot be accounted for by previous simulations. Therefore simulations with more realistic physical conditions are required to explain the results and the particle acceleration and plasma heating processes.

We present a study of the spatial and spectral evolution of the loop-top (LT) sources in a sample of 6 flares near the solar limb observed by {\it RHESSI}. A distinct coronal source, which we identify as the LT source, was seen in each of these flares from the early ``pre-heating'' phase through the late decay phase. Spectral analyses reveal an evident steep power-law component in the pre-heating and impulsive phases, suggesting that the particle acceleration starts upon the onset of the flares. In the late decay phase the LT source has a thermal spectrum and appears to be confined within a small region near the top of the flare loop, and does not spread throughout the loop, as is observed at lower energies. The total energy of this source decreases usually faster than expected from the radiative cooling but much slower than that due to the classical Spitzer conductive cooling along the flare loop. These results indicate the presence of a distinct LT region, where the thermal conductivity is suppressed significantly and/or there is a continuous energy input. We suggest that plasma wave turbulence could play important roles in both heating the plasma and suppressing the conduction during the decay phase of solar flares. With a simple quasi-steady loop model we show that the energy input in the gradual phase can be comparable to that in the impulsive phase and demonstrate how the observed cooling and confinement of the LT source can be used to constrain the wave-particle interaction.

We present data analysis and interpretation of a simple X-class flare observed with RHESSI on November 3, 2003. In contrast to other X-class flares observed previously, this flare shows a very simple morphology with well defined looptop (LT) and footpoint (FP) sources. The almost monotonic upward motion of the LT source and increase in separation of the two FP sources are consistent with magnetic reconnection models proposed for solar flares. In addition, we find that the source motions are relatively slower during the more active phases of hard X-ray emission; the emission centroid of the LT source shifts toward higher altitudes with the increase of energy; the separation between the LT emission centroids at two different photon energies is anti-correlated with the FP flux. Non-uniformity of the reconnecting magnetic fields could be a possible explanation of these features.