The radial acceleration relation (RAR) in galaxies describes a tight empirical scaling law between the total acceleration $g_\mathrm{tot}(r)=GM_\mathrm{tot}(<r)/r^2$ observed in galaxies and that expected from their baryonic mass $g_\mathrm{bar}(r)=GM_\mathrm{bar}(<r)/r^2$, with a characteristic acceleration scale of $g_\dagger\simeq 1.2\times 10^{-10}$ms$^{-2}$. Here, we examine if such a correlation exists in galaxy clusters using weak-lensing, strong-lensing, and X-ray data sets available for 20 high-mass clusters targeted by the CLASH survey. By combining our CLASH data with stellar mass estimates for the brightest cluster galaxies (BCGs) and accounting for the stellar baryonic component in clusters, we determine, for the first time, an RAR on BCG--cluster scales. The resulting RAR is well described by a tight power-law relation, $g_\mathrm{tot}\propto g_\mathrm{bar}^{0.51^{+0.04}_{-0.05}}$, with lognormal intrinsic scatter of $14.7^{+2.9}_{-2.8}\%$. The slope is consistent with the low acceleration limit of the RAR in galaxies, $g_\mathrm{tot}=\sqrt{g_\dagger\,g_\mathrm{bar}}$, whereas the intercept implies a much higher acceleration scale of $g_\ddagger = (2.02\pm0.11)\times 10^{-9}$ms$^{-2}$, indicating that there is no universal RAR that holds on all scales from galaxies to clusters. We find that the observed RAR in CLASH clusters is consistent with predictions from a semi-analytical model developed in the standard $\Lambda$CDM framework. Our results also predict the presence of a baryonic Faber--Jackson relation ($\sigma_v^4\propto M_\mathrm{bar}$) on cluster scales.

We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbinary disk in a highly anisotropic manner. Stream-like SO emission around the circumbinary disk reveals the complicated shocks caused by the interactions between the disk, accretion flows and outflows. A wall-like structure is discovered south of VLA1623B. The discovery of two outflow cavity walls at the same position traveling at different velocities suggests the two outflows from both VLA1623A and VLA1623B overlays on top of each other in the plane of sky. Our detailed flat and flared disk modeling shows that Cycle 2 C18O J = 2-1 data is inconsistent with the combined binary mass of 0.2 Msun as suggested by early Cycle 0 studies. The combined binary mass for VLA1623A should be modified to 0.3 ~ 0.5 Msun.

We study the properties of strange quark stars by employing a 3-flavor Nambu-Jona-Lasinio model with both scalar-isovector and vector-isovector interactions. Using the constraint on the vector-isoscalar interaction strength obtained from the elliptic flow splitting between particles and their antiparticles in relativistic heavy-ion collisions, we investigate the dependence of the properties of strange quark stars on the vector-isovector and the scalar-isovector interactions, and compare the results with the state-of-art astrophysical constraints on the compact star radius and mass as well as its tidal deformability from the GW170817 event. Results from our study reinforce the prospect of using both heavy-ion collisions and astrophysical observations to provide constraints on the isovector coupling strength in quark matter and thus the quark matter equation of state as well as the QCD phase structure at finite isospin chemical potentials.

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

Recently, from the new Spitzer Photometry and Accurate Rotation Curves (SPARC) data, McGaugh et al. (2016) reported a tight Radial Acceleration Relation (RAR) between the observed total acceleration and the acceleration produced by baryons in spiral galaxies. The relation can be fitted by different functions. However, these functions can be discerned if we express the data in the form of halo acceleration relation (HAR). The data reveals a maximum in the halo acceleration. We examined NFW (cusp) and Burkert (core) profiles in the context of dark matter and different parameter families of the interpolating function in the framework of Modified Newtonian Dynamics (MOND).

Processes of dissipation of energy released during capture and disruption stars by Galactic central Black hole. This process should periodically happen in the Galactic center and affect Galactic center environment as well as Galactic halo. Routine star disruptions by the central supermassive black hole BH provide enough cumulative energy to form the Fermi Bubbles. A single star of super-Eddington eruption provide an subsonic outflow gas whose total energy is about $10^{52}$ erg. The average rate of disruption events is expected to be $10^{-4}-10^{-5}$yr$^{-1}$ that provides the average power of energy release in the GC $\dot{W}\sim 3\times 10^{41}$ erg s$^{-1}$, just needed to support Fermi Bubbles in the halo. In the exponential atmosphere the energy from the GC propagates perpendicular to the Galactic plane and provides symmetric features in the halo.

Molecular gas in Arches cloud located near the Arches cluster is one of the emitters of K-$\alpha$ line of neutral iron and X-ray continuum in the Galactic center (GC). Similarly to the cloud Sgr B2, another well-known emitter of the iron line in the GC, the Arches cloud demonstrates temporal decline of the X-ray emission. The most natural origin of this emission is irradiation of primary photons of an X-ray flare from a distant source, most likely Sgr A*. However, recent observations of the Arches cloud discovered variations of equivalent width of the 6.4 keV iron line, which indicated that the X-ray emission from the cloud is a combination of two components with different origin and different equivalent width, one of which is time-variable, while the other is stationary during the period of observations. We considered two different scenarios: a) this emission is formed by reflection from two clouds, which are at some distance from each other, when they are irradiated by two different flares; and b) the other scenario assumes a combination of X-ray fluxes produced in the same cloud by reflection of primary photons and by subrelativistic cosmic rays. We present restrictions for both model and conditions at which these scenarios can be realized. Although none of the models can be completely ruled out, we find that the X-ray reflection model requires less assumption and therefore is the most viable.

A flux of cosmic rays (CRs) propagating through a diffuse ionized gas can excite MHD waves, thus generating magnetic disturbances. We propose a generic model of CR penetration into molecular clouds through their diffuse envelopes, and identify the leading physical processes controlling their transport on the way from a highly ionized interstellar medium to a dense interior of the cloud. The model allows us to describe a transition between a free streaming of CRs and their diffusive propagation, determined by the scattering on the self-generated disturbances. A self consistent set of equations, governing the diffusive transport regime in an envelope and the MHD turbulence generated by the modulated CR flux, is essentially characterized by two dimensionless numbers. We demonstrate a remarkable mutual complementarity of different mechanisms leading to the onset of the diffusive regime, which results in a universal energy spectrum of the modulated CRs. In conclusion, we briefly discuss implications of our results for several fundamental astrophysical problems, such as the spatial distribution of CRs in the Galaxy as well as the ionization, heating, and chemistry in dense molecular clouds.

We study the Mass Discrepancy-Acceleration Relation (MDAR) of 57 elliptical galaxies by their Einstein rings from the Sloan Lens ACS Survey (SLACS). The mass discrepancy between the lensing mass and the baryonic mass derived from population synthesis is larger when the acceleration of the elliptical galaxy lenses is smaller. The MDAR is also related to surface mass density discrepancy. At the Einstein ring, these lenses belong to high-surface-mass density galaxies. Similarly, we find that the discrepancy between the lensing and stellar surface mass density is small. It is consistent with the recent discovery of dynamical surface mass density discrepancy in disk galaxies where the discrepancy is smaller when surface density is larger. We also find relativistic modified Newtonian dynamics (MOND) can naturally explain the MDAR and surface mass density discrepancy in 57 Einstein rings. Moreover, the lensing mass, the dynamical mass and the stellar mass of these galaxies are consistent with each other in relativistic MOND.

After thirty years of the discovery of the fundamental plane, explanations to the tilt of the fundamental plane with respect to the virial plane still suffer from the need of fine-tuning. In this paper, we try to explore the origin of this tilt from the perspective of modified Newtonian dynamics (MOND) by applying the 16 Coma galaxies available in Thomas et al.[1]. Based on the mass models that can reproduce de Vaucouleurs' law closely, we find that the tilt of the traditional fundamental plane is naturally explained by the simple form of the MONDian interpolating function, if we assume a well motivated choice of anisotropic velocity distribution, and adopt the Kroupa or Salpeter stellar mass-to-light ratio. Our analysis does not necessarily rule out a varying stellar mass-to-light ratio.

We give a short review of processes of stochastic acceleration in the Galaxy. We discuss: how to estimate correctly the number of accelerated particles, and at which condition the stochastic mechanism is able to generate power-law nonthermal spectra. We present an analysis of stochastic acceleration in the Galactic halo and discuss whether this mechanism can be responsible for production of high energy electrons there, which emit gamma-ray and microwave emission from the giant Fermi bubbles. Lastly, we discuss whether the effects of stochastic acceleration can explain the CR distribution in the Galactic disk (CR gradient).

We study a possible connection between different non-thermal emissions from the inner few parsecs of the Galaxy. We analyze the origin of the gamma-ray source 2FGL J1745.6-2858 (or 3FGL J1745.6-2859c) in the Galactic Center and the diffuse hard X-ray component recently found by NuSTAR, as well as the radio emission and processes of hydrogen ionization from this area. We assume that a source in the GC injected energetic particles with power-law spectrum into the surrounding medium in the past or continues to inject until now. The energetic particles may be protons, electrons or a combination of both. These particles diffuse to the surrounding medium and interact with gas, magnetic field and background photons to produce non-thermal emissions. We study the spectral and spatial features of the hard X-ray emission and gamma-ray emission by the particles from the central source. Our goal is to examine whether the hard X-ray and gamma-ray emissions have a common origin. Our estimations show that in the case of pure hadronic models the expected flux of hard X-ray emission is too low. Despite protons can produce a non-zero contribution in gamma-ray emission, it is unlikely that they and their secondary electrons can make a significant contribution in hard X-ray flux. In the case of pure leptonic models it is possible to reproduce both X-ray and gamma-ray emissions for both transient and continuous supply models. However, in the case of continuous supply model the ionization rate of molecular hydrogen may significantly exceed the observed value.

The dynamics of an elliptical galaxy within a couple of effective radii can be probed effectively by stars. However, at larger distances planetary nebulae (PNe) replace stars as the tracer of the dynamics. Making use of the motion of PNe, Romanowsky et al. (2003) measured the dynamics of three luminous elliptical galaxies (NGC821, NGC3379, and NGC4494) at large distances from the galactic center. They found that little dark matter is needed up to 6 effective radii. Milgrom & Sanders (2003) showed that this result can be understood in the framework of MOdified Newtonian Dynamics (MOND). As more data are available in the past decade, we revisit this problem. We combine PNe data (up to 6{8 effective radii) and stellar data from SAURON of 7 elliptical galaxies, including those 3 galaxies in Romanowsky et al. (2003) with updated data and 4 other galaxies which have not been analyzed before. We conclude that the dynamics of these galaxies can be well explained by MOND.

Based on Newtonian dynamics, observations show that the luminous masses of astrophysical objects that are the size of a galaxy or larger are not enough to generate the measured motions which they supposedly determine. This is typically attributed to the existence of dark matter, which possesses mass but does not radiate (or absorb radiation). Alternatively, the mismatch can be explained if the underlying dynamics is not Newtonian. Within this conceptual scheme, Modified Newtonian Dynamics (MOND) is a successful theoretical paradigm. MOND is usually expressed in terms of a nonlinear Poisson equation, which is difficult to analyse for arbitrary matter distributions. We study the MONDian gravitational field generated by slightly non-spherically symmetric mass distributions based on the fact that both Newtonian and MONDian fields are conservative (which we refer to as the compatibility condition). As the non-relativistic version of MOND has two different formulations (AQUAL and QuMOND) and the compatibility condition can be expressed in two ways, there are four approaches to the problem in total. The method involves solving a suitably defined linear deformation potential, which generally depends on the choice of MOND interpolation function. However, for some specific form of the deformation potential, the solution is independent of the interpolation function.

Swift, Chandra and XMM have found a weak but nearly constant X-ray component from Swift J1644+57 that appeared at ~500 days and was visible at least until ~ 1400 days after the stellar capture, which cannot be explained by standard tidal disruption theories. We suggest that this X-ray afterglow component may result from Thomson scattering between the primary X-rays and its surrounding plasma, i.e. the Compton echo effect. Similar phenomena has also been observed from molecular clouds in our Galactic Center, which were caused by the past activity of Srg A*. If this interpretation of Swift J1644+57 afterglow is correct, this is the first Compton Echo effect observed in the cosmological distances.

From the rate of hydrogen ionization and the gamma ray flux, we derived the spectrum of relativistic and subrelativistic cosmic rays (CRs) nearby and inside the molecular cloud Sgr B2 near the Galactic Center (GC). We studied two cases of CR propagation in molecular clouds: free propagation and scattering of particles by magnetic fluctuations excited by the neutral gas turbulence. We showed that in the latter case CR propagation inside the cloud can be described as diffusion with the coefficient $\sim 3\times 10^{27}$ cm$^2$ s$^{-1}$. For the case of hydrogen ionization by subrelativistic protons, we showed that their spectrum outside the cloud is quite hard with the spectral index $\delta>-1$. The energy density of subrelativistic protons ($>50$ eV cm$^{-3}$) is one order of magnitude higher than that of relativistic CRs. These protons generate the 6.4 keV emission from Sgr B2, which was about 30\% of the flux observed by Suzaku in 2013. Future observations for the period after 2013 may discover the background flux generated by subrelativistic CRs in Sgr B2. Alternatively hydrogen ionization of the molecular gas in Sgr B2 may be caused by high energy electrons. We showed that the spectrum of electron bremsstrahlung is harder than the observed continuum from Sgr B2, and in principle this X-ray component provided by electrons could be seen from the INTEGRAL data as a stationary high energy excess above the observed spectrum $E_x^{-2}$.

We present the results obtained from linear stability analysis and 2.5-dimensional magnetohydrodynamic (MHD) simulations of the magnetorotational instability (MRI), including the effects of cosmic rays (CRs). We took into account of the CR diffusion along the magnetic field but neglect the cross-field-line diffusion. Two models are considered in this paper: shearing box model and differentially rotating cylinder model. We studied how MRI is affected by the initial CR pressure (i.e., energy) distribution. In the shearing box model, the initial state is uniform distribution. Linear analysis shows that the growth rate of MRI does not depend on the value of CR diffusion coefficient. In the differentially rotating cylinder model, the initial state is a constant angular momentum polytropic disk threaded by weak uniform vertical magnetic field. Linear analysis shows that the growth rate of MRI becomes larger if the CR diffusion coefficient is larger. Both results are confirmed by MHD simulations. The MHD simulation results show that the outward movement of matter by the growth of MRI is not impeded by the CR pressure gradient, and the centrifugal force which acts to the concentrated matter becomes larger. Consequently, the growth rate of MRI is increased. On the other hand, if the initial CR pressure is uniform, then the growth rate of the MRI barely depends on the value of the CR diffusion coefficient.

We analyse the model of stochastic re-acceleration of electrons, which are emitted by supernova remnants (SNRs) in the Galactic Disk and propagate then into the Galactic halo, in order to explain the origin on nonthermal (radio and gamma-ray) emission from the Fermi Bubbles (FB). We assume that the energy for re-acceleration in the halo is supplied by shocks generated by processes of star accretion onto the central black hole. Numerical simulations show that regions with strong turbulence (places for electron re-acceleration) are located high up in the Galactic Halo about several kpc above the disk. The energy of SNR electrons that reach these regions does not exceed several GeV because of synchrotron and inverse Compton energy losses. At appropriate parameters of re-acceleration these electrons can be re-accelerated up to the energy 10E12 eV which explains in this model the origin of the observed radio and gamma-ray emission from the FB. However although the model gamma-ray spectrum is consistent with the Fermi results, the model radio spectrum is steeper than the observed by WMAP and Planck. If adiabatic losses due to plasma outflow from the Galactic central regions are taken into account, then the re-acceleration model nicely reproduces the Planck datapoints.

Fermi bubbles are giant gamma-ray structures extended north and south of the Galactic center with characteristic sizes of order of 10 kpc recently discovered by Fermi Large Area Telescope. Good correlation between radio and gamma-ray emission in the region covered by Fermi bubbles implies the presence of high-energy electrons in this region. Since it is relatively difficult for relativistic electrons of this energy to travel all the way from the Galactic sources toward Fermi bubbles one can assume that they accelerated in-situ. The corresponding acceleration mechanism should also affect the distribution of the relativistic protons in the Galaxy. Since protons have much larger lifetimes the effect may even be observed near the Earth. In our model we suggest that Fermi bubbles are created by acceleration of electrons on series of shocks born due to periodic star accretions by supermassive black hole Sgr A*. We propose that hadronic CR within the "knee" of the observed CR spectrum are produced by Galactic supernova remnants distributed in the Galactic disk. Reacceleration of these particles in the Fermi Bubble produces CRs beyond the knee. This model provides a natural explanation of the observed CR flux, spectral indexes, and matching of spectra at the knee.

We analyse the origin of the gamma-ray flux from the Fermi Bubbles (FBs) in the framework of the hadronic model in which gamma-rays are produced by collisions of relativistic protons with the protons of background plasma in the Galactic halo. It is assumed in this model that the observed radio emission from the FBs is due to synchrotron radiation of secondary electrons produced by $pp$ collisions. However, if these electrons loose their energy by the synchrotron and inverse-Compton, the spectrum of secondary electrons is too soft, and an additional arbitrary component of primary electrons is necessary in order to reproduce the radio data. Thus, a mixture of the hadronic and leptonic models is required for the observed radio flux. It was shown that if the spectrum of primary electrons is $\propto E_e^{-2}$, the permitted range of the magnetic field strength is within 2 - 7 $\mu$G region. The fraction of gamma-rays produced by $pp$ collisions can reach about 80% of the total gamma-ray flux from the FBs. If magnetic field is <2 $\mu$G or >7 $\mu$G the model is unable to reproduce the data. Alternatively, the electrons in the FBs may lose their energy by adiabatic energy losses if there is a strong plasma outflow in the GC. Then, the pure hadronic model is able to reproduce characteristics of the radio and gamma-ray flux from the FBs. However, in this case the required magnetic field strength in the FBs and the power of CR sources are much higher than those followed from observations.

We analyse processes of electron acceleration in the Fermi Bubbles in order to define parameters and restrictions of the models, which are suggested for the origin of these giant radio and gamma-ray structures. In the case of leptonic origin of the nonthermal radiation from the Bubbles, these electrons should be produced somehow in-situ because of relatively short lifetime of high energy electrons, which lose their energy by synchrotron and inverse Compton processes. It has been suggested that electrons in Bubbles may be accelerated by shocks produced by tidal disruption of star accreting onto the central black hole or a process of re-acceleration of electrons ejected by supernova remnants. These processes will be investigated in subsequent papers. In this paper we focus to study in-situ stochastic (Fermi) acceleration by a hydromagnetic/supersonic turbulence, in which electrons can be directly accelerated from the background plasma. We showed that the acceleration from the background plasma is able to explain the observed fluxes of radio and gamma-ray emission from the Bubbles but the range of permitted parameters of the model is strongly restricted.

The single-particle potentials for nucleons and hyperons in neutron-rich matter generally depends on the density and isospin asymmetry of the medium as well as the momentum and isospin of the particle. It further depends on the temperature of the matter if the latter is in thermal equilibrium. We review here the extension of a Gogny-type isospin- and momentum-dependent interaction in several aspects made in recent years and their applications in studying intermediate-energy heavy ion collisions, thermal properties of asymmetric nuclear matter and properties of neutron stars. The importance of the isospin- and momentum-dependence of the single-particle potential, especially the momentum dependence of the isovector potential, is clearly revealed throughout these studies.

We analyse processes of particle acceleration in the Fermi Bubbles. The goal of our investigations is to obtain restrictions for acceleration mechanisms. Our analysis of the three processes: acceleration from background plasma, re-acceleration of relativistic electrons emitted by supernova remnants, and acceleration by shocks generated by processes of star tidal disruption in the Galactic Center, showed that the model of multi-shock acceleration does not have serious objections at present and therefore seems us more attractive than others.

Time delay in galaxy gravitational lensing systems has been used to determine the value of Hubble constant. As in other dynamical phenomena at the scale of galaxy, dark matter is often invoked in gravitational lensing to account for the "missing mass" (the apparent discrepancy between the dynamical mass and the luminous mass). Alternatively, modified gravity can be used to explain the discrepancy. In this paper we adopt the Tensor-Vector-Scalar gravity (TeVeS), a relativistic version of MOdified Newtonian Dynamics (MOND), to study gravitational lensing phenomena and derive the formulae needed to evaluate the Hubble constant. We test our method on quasar lensing by elliptical galaxies in the literature. We focus on double-image systems with time delay measurement. Three candidates are suitable for our study: HE 2149-2745, FBQ J0951+2635 and SBS 0909+532. The Hubble constant obtained is consistent with the value used in fitting the CMB result in neutrino cosmological model.

The origin of hard X-ray (HXR) excess emission from clusters of galaxies is still an enigma, whose nature is debated. One of the possible mechanism to produce this emission is the bremsstrahlung model. However, previous analytical and numerical calculations showed that in this case the intracluster plasma had to be overheated very fast because suprathermal electrons emitting the HXR excess lose their energy mainly by Coulomb losses, i.e., they heat the background plasma. It was concluded also from these investigations that it is problematic to produce emitting electrons from a background plasma by stochastic (Fermi) acceleration because the energy supplied by external sources in the form of Fermi acceleration is quickly absorbed by the background plasma. In other words the Fermi acceleration is ineffective for particle acceleration. We revisited this problem and found that at some parameter of acceleration the rate of plasma heating is rather low and the acceleration tails of non-thermal particles can be generated and exist for a long time while the plasma temperature is almost constant. We showed also that for some regime of acceleration the plasma cools down instead of being heated up, even though external sources (in the form of external acceleration) supply energy to the system. The reason is that the acceleration withdraws effectively high energy particles from the thermal pool (analogue of Maxwell demon).

The J-Ks color distribution (CD) with a bin size of 0.05 magnitude for the entire Milky Way has been carried out by using the Two Micron All Sky Survey Point Source Catalog (2MASS PSC). The CDs are bimodal, which has a red peak at 0.8 < J-Ks < 0.85 and a blue peak at 0.3 < J-Ks < 0.4. The colors of the red peak are more or less the same for the whole sky, but that of the blue peak depend on Galactic latitude, (J-Ks ~ 0.35 at low Galactic latitudes and 0.35 < J-Ks < 0.4 for other sky areas). The blue peak dominates the bimodal CDs at low Galactic latitudes and becomes comparable with the red peak in other sky regions. In order to explain the bimodal distribution and the global trend shown by the all sky 2MASS CDs, we assemble an empirical HR diagram, which is composed by observational-based near infrared HR diagrams and color magnitude diagrams, and incorporate a Milky Way model. In the empirical HR diagram, the main sequence stars turnoff the thin disk is relatively bluer, (J-Ks)0 = 0.31, when we compare with the thick disk which is (J-Ks)0 = 0.39. The age of the thin/thick disk is roughly estimated to be around 4-5/8-9 Gyr according to the color-age relation of the main sequence turnoff. In general, the 2MASS CDs can be treated as a tool to census the age of stellar population of the Milky Way in a statistical manner and to our knowledge this is a first attempt to measure the age.

The {\it Fermi} Large Area Telescope has recently discovered two giant gamma-ray bubbles which extend north and south of the Galactic center with diameters and heights of the order of $H\sim 10$ kpc. We suggest that the periodic star capture processes by the Galactic supermassive black hole Sgr A$^*$, with a capture rate of $\tau_{\rm cap}^{-1}\sim 3\times 10^{-5}$ yr$^{-1}$ and an energy release of $W\sim 3\times 10^{52}$ erg per capture, can result in hot plasma injecting into the Galactic halo at a wind velocity of $u\sim 10^8$ cm s$^{-1}$. The periodic injection of hot plasma can produce a series of shocks. Energetic protons in the bubble are re-accelerated when they interact with these shocks. We show that for energy larger than $E> 10^{15}$ eV, the acceleration process can be better described by the stochastic second-order Fermi acceleration. We propose that hadronic cosmic rays (CRs) within the ``knee'' of the observed CR spectrum are produced by Galactic supernova remnants distributed in the Galactic disk. Re-acceleration of these particles in the Fermi Bubble produces CRs beyond the knee. With a mean CR diffusion coefficient in this energy range in the bubble $D_B\sim 3\times 10^{30}$ cm$^2$ s$^{-1}$, we can reproduce the spectral index of the spectrum beyond the knee and within. The conversion efficiency from shock energy of the bubble into CR energy is about 10\%. This model provides a natural explanation of the observed CR flux, spectral indices, and matching of spectra at the knee.

Fermi has discovered two giant gamma-ray-emitting bubbles that extend nearly 10 kpc in diameter. We propose that periodic star capture processes by the galactic supermassive black hole, Sgr A*, with a capture rate $<10^{-5}$ yr$^{-1}$ and energy release $\sim 10^{52}$ erg per one capture can produce shocks in the halo, which accelerate electrons to the energy ~ 1 TeV. These electrons generate radio emission via synchrotron radiation, and gamma-rays via inverse Compton scattering with the relic and the galactic soft photons. Estimates of the diffusion coefficient from the observed gamma-ray flux explains consistently the necessary maximum energy of electrons and sharp edges of the bubble.

Fermi LAT has discovered two extended gamma-ray bubbles above and below the galactic plane. We propose that their origin is due to the energy release in the Galactic center (GC) as a result of quasi-periodic star accretion onto the central black hole. Shocks generated by these processes propagate into the Galactic halo and accelerate particles there. We show that electrons accelerated up to ~10 TeV may be responsible for the observed gamma-ray emission of the bubbles as a result of inverse Compton (IC) scattering on the relic photons. We also suggest that the Bubble could generate the flux of CR protons at energies > 10^15 eV because the shocks in the Bubble have much larger length scales and longer lifetimes in comparison with those in SNRs. This may explain the the CR spectrum above the knee.

Parker instability arises from the presence of magnetic fields in a plasma in a gravitational field such as the interstellar medium (ISM), wherein the magnetic buoyant pressure expels the gas and causes the gas to move along the field lines. The process of mixing of this instability in the ISM near the Galactic plane is investigated. The initial ISM is assumed to consist of two fluids: plasma gas and cosmic-ray particles, in hydrostatic equilibrium, coupled with a uniform, azimuthally-aligned magnetic field. The evolution of the instability is explored in two models: an isothermal exponential-declining density model and a two-layered, hyperbolic tangent temperature model. After a small perturbation, the unstable gas aggregates at the bottom of the magnetic loops and forms dense blobs. The growth rate of the instability decreases as the coupling between the cosmic rays and the plasma becomes stronger (meaning a smaller CR diffusion coefficient). The mixing is enhanced by the cosmic-ray diffusion, while the shape of the condensed gas depends sensitively on the initial equilibrium conditions. The hyperbolic tangent temperature model produces a more concentrated and round shape of clumps at the foot points of rising magnetic arches, like the observed giant molecular cloud, whereas the exponential density model gives rise to a filamentary morphology of the clumpy structure. When considering a minimum perpendicular or cross field diffusion of cosmic rays, which is often substantially smaller than the parallel coefficient $\kappa_{\|}$, around $2%-4%$ of $\kappa_{\|}$, the flow speed is significantly increased such that the magnetic loops extend to a greater altitude. We speculate that the galactic wind flow perpendicular to the galactic disk may be facilitated by Parker instabilities through the cross field diffusion of cosmic rays.

The Ks band differential star count of the Two Micron All Sky Survey (2MASS) is used to derive the global structure parameters of the smooth components of the Milky Way. To avoid complication introduced by other fine structures and significant extinction near and at the Galactic plane, we only consider Galactic latitude |b| > 30 degree data. The star count data is fitted with a threecomponent model: double exponential thin disk and thick disk, and a power law decay oblate halo. Using maximum likelihood the best-fit local density of thin disk is n0 = 0.030 +- 0.002 stars/pc^3. The best-fit scale-height and length of the thin disk are Hz1 = 360+-10 pc and Hr1 = 3.7+-1.0 kpc, and those of the thick disk are and Hz2 = 1020+-30 pc and Hr2 = 5.0+-1.0 kpc, the local thick-to-thin disk density ratio is f2 = 7+-1%. The best-fit axis ratio, power law index and local density ratio of the oblate halo are kappa = 0.55+-0.15, p = 2.6+-0.6 and fh = 0.20+-0:10%, respectively. Moreover, we find some degeneracy among the key parameters (e.g. n0,Hz1, f2 and Hz2). Any pair of these parameters are anticorrelated to each other. The 2MASS data can be well-fitted by several possible combinations of parameters. This is probably the reason that there is a wide range of values for the structure parameters in literature similar to this study. Since only medium and high Galactic latitude data are analyzed, the fitting is very insensitive to the scale-lengths of the disks.

We show in the Skyrme-Hartree-Fock approach that unambiguous correlations exist between observables of finite nuclei and nuclear matter properties. Using this correlation analysis to existing data on the neutron skin thickness of Sn isotopes, we find important constraints on the value E_{sym}(rho_0) and density slope L of the nuclear symmetry energy at saturation density. Combining these constraints with those from recent analyses of isospin diffusion and double neutron/proton ratio in heavy ion collisions leads to a value of L=58\pm 18 MeV approximately independent of E_{sym}(\rho_0).

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.

Fermi has discovered two giant gamma-ray-emitting bubbles that extend nearly 10kpc in diameter north and south of the galactic center (GC). The existence of the bubbles was first evidenced in X-rays detected by ROSAT and later WMAP detected an excess of radio signals at the location of the gamma-ray bubbles. We propose that periodic star capture processes by the galactic supermassive black hole, Sgr A$^*$, with a capture rate $3\times 10^{-5}$yr$^{-1}$ and energy release $\sim 3\times 10^{52}$erg per capture can produce very hot plasma $\sim 10$keV with a wind velocity $\sim 10^8$cm/s injected into the halo and heat up the halo gas to $\sim 1$keV, which produces thermal X-rays. The periodic injection of hot plasma can produce shocks in the halo and accelerate electrons to $\sim$TeV, which produce radio emission via synchrotron radiation, and gamma-rays via inverse Compton scattering with the relic and the galactic soft photons.

This study investigates Parker instability in an interstellar medium (ISM) near the Galactic plane using three-dimensional magneto-hydrodynamic simulations. Parker instability arises from the presence of a magnetic field in a plasma, wherein the magnetic buoyant pressure expels the gas and cause the gas to move along the field lines. The process is thought to induce the formation of giant molecular clouds in the Galaxy. In this study, the effects of cosmic-ray (CR) diffusion are examined. The ISM at equilibrium is assumed to comprise a plasma fluid and a CR fluid at various temperatures, with a uniform magnetic field passing through it in the azimuthal direction of the Galactic disk. After a small perturbation, the unstable gas aggregates at the footpoint of the magnetic fields and forms dense blobs. The growth rate of the instability increases with the strength of the CR diffusion. The formation of dense clouds is enhanced by the effect of cosmic rays (CRs), whereas the shape of the clouds depends sensitively on the initial conditions of perturbation.

Describing the hyperonic and quark phases of neutron stars with an isospin- and momentum-dependent effective interaction for the baryon octet and the MIT bag model, respectively, and using the Gibbs conditions to construct the mixed phase, we study the energy release due to the hadron-quark phase transition. Moreover, the frequency and damping time of the first axial $w$-mode of gravitational waves are studied for both hyperonic and hybrid stars. We find that the energy release is much more sensitive to the bag constant than the density dependence of the nuclear symmetry energy. Also, the frequency of the $w$-mode is found to be significantly different with or without the hadron-quark phase transition and depends strongly on the value of the bag constant. Effects of the density dependence of the nuclear symmetry energy become, however, important for large values of the bag constant that lead to higher hadron-quark transition densities.

By assuming that the isospin- and momentum-dependent MDI interaction has a form similar to the Gogny-like effective two-body interaction with a Yukawa finite-range term and the momentum dependence only originates from the finite-range exchange interaction, we determine its parameters by comparing the predicted potential energy density functional in uniform nuclear matter with what has been usually given and used extensively in transport models for studying isospin effects in intermediate-energy heavy-ion collisions as well as in investigating the properties of hot asymmetric nuclear matter and neutron star matter. We then use the density matrix expansion to derive from the resulting finite-range exchange interaction an effective Skyrme-like zero-range interaction with density-dependent parameters. As an application, we study the transition density and pressure at the inner edge of neutron star crusts using the stability conditions derived from the linearized Vlasov equation for the neutron star matter.

Expressing explicitly the parameters of the standard Skyrme interaction in terms of the macroscopic properties of asymmetric nuclear matter, we show in the Skyrme-Hartree-Fock approach that unambiguous correlations exist between observables of finite nuclei and nuclear matter properties. We find that existing data on neutron skin thickness $\Delta r_{np}$ of Sn isotopes give an important constraint on the symmetry energy $E_{sym}({\rho _{0}})$ and its density slope $L$ at saturation density ${\rho _{0}}$. Combining these constraints with those from recent analyses of isospin diffusion and double neutron/proton ratio in heavy-ion collisions at intermediate energies leads to a more stringent limit on $L$ approximately independent of $E_{sym}({\rho _{0}})$. The implication of these new constraints on the $\Delta r_{np}$ of $^{208}$Pb as well as the core-crust transition density and pressure in neutron stars is discussed.

Using strong lensing data Milgrom's MOdified Newtonian Dynamics (MOND) or its covariant TeVeS (Tensor-Vector-Scalar Theory) is being examined here as an alternative to the conventional $\Lambda$CDM paradigm. We examine 10 double-image gravitational lensing systems, in which the lens masses have been estimated by stellar population synthesis models. While mild deviations exist, we do not find out that strong cases for outliers to the TeVeS theory.

Using the momentum-dependent MDI effective interaction for nucleons, we have studied the transition density and pressure at the boundary between the inner crust and liquid core of hot neutron stars. We find that their values are larger in neutrino-trapped neutron stars than in neutrino-free neutron stars. Furthermore, both are found to decrease with increasing temperature of a neutron star as well as increasing slope parameter of the nuclear symmetry energy, except that the transition pressure in neutrino-trapped neutron stars for the case of small symmetry energy slope parameter first increases and then decreases with increasing temperature. We have also studied the effect of the nuclear symmetry energy on the critical temperature above which the inner crust in a hot neutron star disappears and found that with increasing value of the symmetry energy slope parameter, the critical temperature decreases slightly in neutrino-trapped neutron stars but first decreases and then increases in neutrino-free neutron stars.

The isospin- and momentum-dependent MDI interaction, which has been extensively used in intermediate-energy heavy-ion reactions to study the properties of asymmetric nuclear matter, is extended to include the nucleon-hyperon and hyperon-hyperon interactions by assuming same density, momentum and isospin dependence as for the nucleon-nucleon interaction. The parameters in these interactions are determined from the empirical hyperon single-particle potentials in symmetric nuclear matter at saturation density. The extended MDI interaction is then used to study in the mean-field approximation the equation of state of hypernuclear matter and also the properties of hybrid stars by including the phase transition from the hypernuclear matter to the quark matter at high densities. In particular, the effects of attractive and repulsive $\Sigma$N interactions and different values of symmetry energies on the hybrid star properties are investigated.

Using the Hugenholtz-Van Hove theorem, we derive general expressions for the quadratic and quartic symmetry energies in terms of single-nucleon potentials in isospin asymmetric nuclear matter. These analytical relations are useful for gaining deeper insights into the microscopic origins of the uncertainties in our knowledge on nuclear symmetry energies especially at supra-saturation densities. As examples, the formalism is applied to two model single-nucleon potentials that are widely used in transport model simulations of heavy-ion reactions.

We discuss the origin of continuum and line X-ray emission observed in the direction the Galactic Center. We predict a significant flux of de-excitation gamma-ray lines in this direction, which can be produced by subrelativistic protons generated by accretion processes.

The origin and properties of the source of positrons annihilating in the Galactic Center is still a mystery. One of the criterion, which may discriminate between different mechanisms of positron production there, is the positron energy injection. Beacom and Yueksel (2006) suggested a method to estimate this energy from the ratio of the 511 keV line to the MeV in-flight annihilation fluxes. From the COMPTEL data they derived that the maximum injection energy of positron should be about several MeV that cut down significantly a class of models of positron origin in the GC assuming that positrons lose their energy by Coulomb collisions only. However, observations show that the strength of magnetic field in the GC is much higher than in other parts of the Galaxy, and it may range there from 100 $\mu$G to several mG. In these conditions, synchrotron losses of positrons are significant that extends the range of acceptable values of positron injection energy. We show that if positrons injection in the GC is non-stationary and magnetic field is higher than 0.4 mG both radio and gamma-ray restrictions permit their energy to be higher than several GeV.

Using the nuclear symmetry energy that has been recently constrained by the isospin diffusion data in intermediate-energy heavy ion collisions, we have studied the transition density and pressure at the inner edge of neutron star crusts, and they are found to be 0.040 fm$^{-3}$ $\leq \rho_{t}\leq 0.065$ fm$^{-3}$ and 0.01 MeV/fm$^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm$^{3}$, respectively, in both the dynamical and thermodynamical approaches. We have also found that the widely used parabolic approximation to the equation of state of asymmetric nuclear matter gives significantly higher values of core-crust transition density and pressure, especially for stiff symmetry energies. With these newly determined transition density and pressure, we have obtained an improved relation between the mass and radius of neutron stars.

Aims. The accretion of stars onto the central supermassive black hole at the center of the Milky Way is predicted to generate large fluxes of subrelativistic ions in the Galactic center region. We analyze the intensity, shape and spatial distribution of de-excitation gamma-ray lines produced by nuclear interactions of these energetic particles with the ambient medium. Methods. We first estimate the amount and mean kinetic energy of particles released from the central black hole during star disruption. We then calculate from a kinetic equation the energy and spatial distributions of these particles in the Galactic center region. These particle distributions are then used to derive the characteristics of the main nuclear interaction gamma-ray lines. Results. Because the time period of star capture by the supermassive black hole is expected to be shorter than the lifetime of the ejected fast particles against Coulomb losses, the gamma-ray emission is predicted to be stationary. We find that the nuclear de-excitation lines should be emitted from a region of maximum 5$^\circ$ angular radius. The total gamma-ray line flux below 8 MeV is calculated to be $\approx10^{-4}$ photons cm$^{-2}$ s$^{-1}$. The most promising lines for detection are those at 4.44 and $\sim$6.2 MeV, with a predicted flux in each line of $\approx$$10^{-5}$ photons cm$^{-2}$ s$^{-1}$. Unfortunately, it is unlikely that this emission can be detected with the INTEGRAL observatory. But the predicted line intensities appear to be within reach of future gamma-ray space instruments. A future detection of de-excitation gamma-ray lines from the Galactic center region would provide unique information on the high-energy processes induced by the central supermassive black hole and the physical conditions of the emitting region.

Analytical expressions for the saturation density as well as the binding energy and incompressibility at the saturation density of asymmetric nuclear matter are given exactly up to 4th-order in the isospin asymmetry delta =(rho_n - rho_p)/rho using 11 characteristic parameters defined at the normal nuclear density rho_0. Using an isospin- and momentum-dependent modified Gogny (MDI) interaction and the SHF approach with 63 popular Skyrme interactions, we have systematically studied the isospin dependence of the saturation properties of asymmetric nuclear matter, particularly the incompressibility $K_{sat}(\delta )=K_{0}+K_{sat,2}\delta ^{2}+K_{sat,4}\delta ^{4}+O(\delta ^{6})$ at the saturation density. Our results show that the magnitude of the high-order $K_{sat,4}$ parameter is generally small compared to that of the $K_{\sat,2}$ parameter. The latter essentially characterizes the isospin dependence of the incompressibility at the saturation density and can be expressed as $K_{sat,2}=K_{sym}-6L-\frac{J_{0}}{K_{0}}L$, Furthermore, we have constructed a phenomenological modified Skyrme-like (MSL) model which can reasonably describe the general properties of symmetric nuclear matter and the symmetry energy predicted by both the MDI model and the SHF approach. The results indicate that the high-order $J_{0}$ contribution to $K_{sat,2}$ generally cannot be neglected. In addition, it is found that there exists a nicely linear correlation between $K_{sym}$ and $L$ as well as between $J_{0}/K_{0}$ and $K_{0}$. These correlations together with the empirical constraints on $K_{0}$, $L$, $E_{sym}(\rho_{0})$ and the nucleon effective mass lead to an estimate of $K_{sat,2}=-370\pm 120$ MeV.

We analyse new results of Chandra and Suzaku which found a flux of hard X-ray emission from the compact region around Sgr A$^\ast$ (r ~ 100 pc). We suppose that this emission is generated by accretion processes onto the central supermassive blackhole when an unbounded part of captured stars obtains an additional momentum. As a result a flux of subrelativistic protons is generated near the Galactic center which heats the background plasma up to temperatures about 6-10 keV and produces by inverse bremsstrahlung a flux of non-thermal X-ray emission in the energy range above 10 keV.

We showed that if the non-thermal emission from the Galactic center in the range 14-40 keV is due to inverse bremsstrahlung emission of subrelativistic protons, their interactions with hot and cold fractions of the interstellar medium are equally important. Our estimation show that about 30% of the total non-thermal flux from the GC in the range 14-40 keV is generated in regions of cold gas while the rest is produced by proton interaction with hot plasma. From the spatial distribution of 6.7 keV iron line we concluded the spatial distribution of hot plasma is strongly non-uniform that should be taken into account in analysis of protons propagation in the GC. From the Suzaku data we got independent estimates for the diffusion coefficient of subrelativistic protons in the GC, which was in the range $ 10^{26} - 10^{27}$ cm$^2$s$^{-1}$

We analyze the 6.4 keV line and continuum emission from the molecular cloud Sgr B2 and the source HESS J1745-303, which is supposed to be a complex of molecular gas. From the HESS results it follows that Sgr A is a source of high energy protons, which penetrate into molecular clouds producing there a TeV gamma-ray flux. We present arguments that Sgr A may also produce a flux of subrelativistic protons which generate the 6.4 keV line and bremsstrahlung continuum emission from the clouds.

We use N-body simulations to survey the response of embedded star clusters to the dispersal of their parent molecular cloud. The final stages of the clusters can be divided into three classes: the cluster (i) is destroyed, (ii) has a loose structure, and (iii) has a compact core. We are interested in three of the governing parameters of the parent cloud: (i) the mass, (ii) the size, and (iii) the dispersing rate. It is known that the final stage of the cluster is well correlated with the star formation efficiency (SFE) for systems with the same cluster and cloud profile. We deem that the SFE alone is not enough to address systems with clouds of different sizes. Our result shows that the initial cluster-cloud mass ratio at a certain Lagrangian radius, and the initial kinetic energy are better indicators for the survivability of embedded clusters.

Heavy-ion reactions especially those induced by radioactive beams provide useful information about the density dependence of the nuclear symmetry energy, thus the Equation of State of neutron-rich nuclear matter, relevant for many astrophysical studies. The latest developments in constraining the symmetry energy at both sub- and supra-saturation densities from analyses of the isopsin diffusion and the $\pi^-/\pi^+$ ratio in heavy-ion collisions using the IBUU04 transport model are discussed. Astrophysical ramifications of the partially constrained symmetry energy on properties of neutron star crusts, gravitational waves emitted by deformed pulsars and the w-mode oscillations of neutron stars are presented briefly.

Modified Newtonian Dynamics (MOND) and its relativistic version - TeVeS offer us an alternative perspective to understand the universe without the demand of the elusive cold dark matter. This MONDian paradigm is not only competitive with the conventional CDM in a large range of scales, but also even more successful in the galactic scale. Recently, by studying 6 lensing systems, Ferreras et al. (2008) claimed that MOND still needs dark matter even in galactic scales. When we study the same systems, however, we yield an opposite conclusion. In this contribution, we report our result and conclude that MOND does not need dark matter in galactic lensing systems. Furthermore, we extend our study to 22 SLACS (Sloan Lens ACS Survey) lenses, and obtain the same conclusion as well, i.e., no dark matter is needed in elliptical galaxies.

The whole sky differential star counts (DSC) with 1 degree resolution are retrieved from 2MASS online data service. Galaxy with double exponential thin and thick disks and a single power law luminosity function (LF) is used to interpret the 2MASS data. The slope of the DSC appears roughly isotropic over the whole sky, the average value is ~0.32, which corresponds to a power law index ~1.8 of the LF. We find that the scale-length and scale-height the thin disk are ~3.0 kpc and ~245 pc, and those of the thick disk are ~3.0 kpc and ~780 pc. The ratio of the thick disk to the thin disk is ~7%. The location of Sun above the disk is ~15 pc. A comparison of the data and model and their discrepancy are also provided.

Recent analyses of several isospin effects in heavy-ion reactions have allowed us to constrain the density dependence of nuclear symmetry energy at sub-saturation densities within a narrow range. Combined with constraints on the Equation of State (EOS) of symmetric nuclear matter obtained previously from analyzing the elliptic flow in relativistic heavy-ion collisions, the EOS of neutron-rich nuclear matter is thus partially constrained. Here we report effects of the partially constrained EOS of neutron-rich nuclear matter on the mass-radius correlation, moment of inertia, elliptical deformation and gravitational radiation of (rapidly) rotating neutron stars.

The ultimate goal of studying isospin physics via heavy-ion reactions with neutron-rich, stable and/or radioactive nuclei is to explore the isospin dependence of in-medium nuclear effective interactions and the equation of state of neutron-rich nuclear matter, particularly the isospin-dependent term in the equation of state, i.e., the density dependence of the symmetry energy. Because of its great importance for understanding many phenomena in both nuclear physics and astrophysics, the study of the density dependence of the nuclear symmetry energy has been the main focus of the intermediate-energy heavy-ion physics community during the last decade, and significant progress has been achieved both experimentally and theoretically. In particular, a number of phenomena or observables have been identified as sensitive probes to the density dependence of the nuclear symmetry energy. Experimental studies have confirmed some of these interesting isospin-dependent effects and allowed us to constrain relatively stringently the symmetry energy at sub-saturation densities. The impacts of this constrained density dependence of the symmetry energy on the properties of neutron stars have also been studied, and they were found to be very useful for the astrophysical community. With new opportunities provided by the various radioactive beam facilities being constructed around the world, the study of isospin physics is expected to remain one of the forefront research areas in nuclear physics. In this report, we review the major progress achieved during the last decade in isospin physics with heavy ion reactions and discuss future challenges to the most important issues in this field.

Star clusters are formed in molecular clouds which are believed to be the birth places of most stars. From recent observational data, Lada & Lada 2003 estimated that only 4% to 7% of the clusters embedded inside molecular clouds have survived. An important mechanism for the disruption of embedded (bound)-clusters is the dispersion of the parent cloud by UV radiation, stellar winds and/or supernova explosions. In this work we study the effect of this mechanism by N-body simulations. We find that most embedded-clusters survive for more than 30 Myr even when different initial conditions of the cluster may introduce some minor variations, but the general result is rather robust.

Using various relativistic mean-field models, including the nonlinear ones with meson field self-interactions, those with density-dependent meson-nucleon couplings, and the point-coupling models without meson fields, we have studied the isospin-dependent bulk and single-particle properties of asymmetric nuclear matter. In particular, we have determined the density dependence of nuclear symmetry energy from these different relativistic mean-field models and compare the results with the constraints recently extracted from analyses of experimental data on isospin diffusion and isotopic scaling in intermediate-energy heavy ion collisions as well as from measured isotopic dependence of the giant monopole resonances in even-A Sn isotopes. Among the 23 parameter sets in the relativistic mean-filed model that are commonly used for nuclear structure studies, only a few are found to give symmetry energies that are consistent with the empirical constraints. We have also studied the nuclear symmetry potential and the isospin-splitting of the nucleon effective mass in isospin asymmetric nuclear matter. We find that both the momentum dependence of the nuclear symmetry potential at fixed baryon density and the isospin-splitting of the nucleon effective mass in neutron-rich nuclear matter depend not only on the nuclear interactions but also on the definition of the nucleon optical potential.

Heavy-ion reactions provide a unique means to investigate the equation of state (EOS) of neutron-rich nuclear matter, especially the density dependence of the nuclear symmetry energy $E_{sym}(\rho)$. The latter plays an important role in understanding many key issues in both nuclear physics and astrophysics. Recent analyses of heavy-ion reactions have already put a stringent constraint on the $E_{sym}(\rho)$ around the saturation density. This subsequently allowed us to constrain significantly the radii and cooling mechanisms of neutron stars as well as the possible changing rate of the gravitational constant G.

Heavy-ion reactions induced by neutron-rich nuclei provide a unique means to investigate the equation of state of isospin-asymmetric nuclear matter, especially the density dependence of the nuclear symmetry energy. In particular, recent analyses of the isospin diffusion data in heavy-ion reactions have already put a stringent constraint on the nuclear symmetry energy around the nuclear matter saturation density. We review this exciting result and discuss its implications on nuclear effective interactions and the neutron skin thickness of heavy nuclei. In addition, we also review the theoretical progress on probing the high density behaviors of the nuclear symmetry energy in heavy-ion reactions induced by high energy radioactive beams.

Star clusters are formed in molecular clouds which are believed to be the birth places of most stars. From recent observational data, Lada & Lada(2003) estimated that only 4 to 7% of the proto-clusters have survived. Many factors could cause this high infant mortality. Galactic tidal forces, close encounters with molecular clouds and shock heating are among the possible causes but they have a longer timescale than typical lifetime of molecular clouds. Another possible reason is mass loss in very beginning of cluster evolution in the form of UV radiation, stellar winds or supernova explosions. Mass loss is the main factor we study in this work by using N-body simulations. We find that most proto-clusters survive for more than 40 Myr even when the mass loss rate is high.

Using the relativistic impulse approximation with the Love-Franey \textsl{NN} scattering amplitude developed by Murdock and Horowitz, we investigate the low-energy (100 MeV$\leq E_{\mathrm{kin}}\leq 400$ MeV) behavior of the nucleon Dirac optical potential, the Schr\"{o}dinger-equivalent potential, and the nuclear symmetry potential in isospin asymmetric nuclear matter. We find that the nuclear symmetry potential at fixed baryon density decreases with increasing nucleon energy. In particular, the nuclear symmetry potential at saturation density changes from positive to negative values at nucleon kinetic energy of about 200 MeV. Furthermore,the obtained energy and density dependence of the nuclear symmetry potential is consistent with those of the isospin- and momentum-dependent MDI interaction with $x=0$, which has been found to describe reasonably both the isospin diffusion data from heavy-ion collisions and the empirical neutron-skin thickness of $^{208} $Pb.

Recent analysis of the isospin diffusion data from heavy-ion collisions based on an isospin- and momentum-dependent transport model with in-medium nucleon-nucleon cross sections has led to the extraction of a value of $L=88\pm 25$ MeV for the slope of the nuclear symmetry energy at saturation density. This imposes stringent constraints on both the parameters in the Skyrme effective interactions and the neutron skin thickness of heavy nuclei. Among the 21 sets of Skyrme interactions commonly used in nuclear structure studies, the 4 sets SIV, SV, G$_\sigma$, and R$_\sigma$ are found to give $L$ values that are consistent with the extracted one. Further study on the correlations between the thickness of the neutron skin in finite nuclei and the nuclear matter symmetry energy in the Skyrme Hartree-Fock approach leads to predicted thickness of the neutron skin of $0.22\pm 0.04$ fm for $^{208}$Pb, $0.29\pm 0.04$ fm for $^{132}$Sn, and $0.22\pm 0.04$ fm for $^{124}$Sn.

The stochastic acceleration of subrelativistic electrons from a background plasma is studied in order to find a possible explanation of the hard X-ray (HXR) emission detected from the Coma cluster. We calculate the necessary energy supply as a function of the plasma temperature and of the electron energy. We show that, for the same value of the HXR flux, the energy supply changes gradually from its high value (when emitting particle are non-thermal) to lower values (when the electrons are thermal). The kinetic equations we use include terms describing particle thermalization as well as momentum diffusion due to the Fermi II acceleration. We show that the temporal evolution of the particle distribution function has, at its final stationary stage, a rather specific form: it cannot be described by simple exponential or power-law expressions. A broad transfer region is formed by Coulomb collisions at energies between the Maxwellian and power-law parts of the distribution function. In this region the radiative lifetime of a single electron differs greatly from the lifetime of the distribution function as a whole. For a plasma temperature of 8 keV, the particles emitting bremsstrahlung at 20-80 keV lie in this quasi-thermal regime. We show that the energy supply required by quasi-thermal electrons to produce the observed HXR flux from Coma is one or two orders of magnitude smaller than the value derived from the assumption of a nonthermal origin of the emitting particles. This result may solve the problem of rapid cluster overheating by nonthermal electrons. We finally predict the change in Coma's SZ effect caused by the distortions of the Maxwellian electron spectrum, and we show that evidence for acceleration of subrelativistic electrons can be derived from detailed spectral measurements.

Nuclear reactions induced by stable and/or radioactive neutron-rich nuclei provide the opportunity to pin down the equation of state of neutron-rich matter, especially the density ($\rho$) dependence of its isospin-dependent part, i.e., the nuclear symmetry energy $E_{\rm sym}$. A conservative constraint, $32(\rho /\rho_{0})^{0.7} < E_{\rm sym}(\rho ) < 32(\rho /\rho _{0})^{1.1}$, around the nuclear matter saturation density $\rho_0$ has recently been obtained from the isospin diffusion data in intermediate energy heavy-ion collisions. We review this exciting result and discuss its consequences and implications on nuclear effective interactions, radii and cooling mechanisms of neutron stars.

Since Bekenstein's (2004) creation of his Tensor-Vector-Scalar theory (TeVeS), the Modified Newtonian dynamics (MOND) paradigm has been redeemed from the embarrassment of lacking a relativistic version. One primary success of TeVeS is that it provides an enhancement of gravitational lensing, which could not be achieved by other MONDian theories. Following Bekenstein's work, we investigate the phenomena of gravitational lensing including deflection angles, lens equations and time delay. We find that the deflection angle would maintain its value while the distance of closest approach vary in the MOND regime, this coincides with the conclusion of Mortlock and Turner's (2001) intuitional approach. Moreover, the scalar field, which is introduced to enhance the deflection angle in TeVeS, contributes a negative effect on the potential time delay. Unfortunately this phenomenon is unmeasurable in lensing systems where we can only observe the time delay between two images for a given source. However, this measurable time delay offers another constraint on the mass ratio of the DM and MOND scenarios, which in general differs from that given by the deflection angle. In other words, for a lensing system, if two masses, m_gN and m_gM, are mutually alternatives for the deflection angles in their own paradigm, regarding the time delay they are in general in an exclusive relation.

We use linear analysis to examine the effect of cosmic rays (CRs) on the Parker-Jeans instability of magnetized self-gravitating gaseous disks. We adopt a slab equilibrium model in which the gravity (including self-gravity) is perpendicular to the mid-plane, the magnetic field lies along the slab. CR is described as a fluid and only along magnetic field lines diffusion is considered. The linearised equations are solved numerically. The system is susceptible to Parker-Jeans instability. In general the system is less unstable when the CR diffusion coefficient is smaller (i.e., the coupling between the CRs and plasma is stronger). The system is also less unstable if CR pressure is larger. This is a reminiscence of the fact that Jeans instability and Parker instability are less unstable when the gas pressure is larger (or temperature is higher). Moreover, for large CR diffusion coefficient (or small CR pressure), perturbations parallel to the magnetic field are more unstable than those perpendicular to it. The other governing factor on the growth rate of the perturbations in different directions is the thickness of the disk or the strength of the external pressure on the disk. In fact, this is the determining factor in some parameter regimes.

The latest development in determining the density dependence of the nuclear symmetry energy using heavy-ion collisions is reviewed. Within the IBUU04 version of an isospin- and momentum-dependent transport model using a modified Gogny effective interaction, recent experimental data from NSCL/MSU on isospin diffusion are found to be consistent with a nuclear symmetry energy of $E_{sym}(\rho)\approx 31.6(\rho /\rho_{0})^{1.05}$ at subnormal densities. Predictions on several observables sensitive to the density dependence of the symmetry energy at supranormal densities accessible at GSI and the planned Rare Isotope Accelerator (RIA) are also made.

With an isospin- and momentum-dependent transport model, we find that the degree of isospin diffusion in heavy ion collisions at intermediate energies is affected by both the stiffness of the nuclear symmetry energy and the momentum dependence of the nucleon potential. Using a momentum dependence derived from the Gogny effective interaction, recent experimental data from NSCL/MSU on isospin diffusion are shown to be consistent with a nuclear symmetry energy given by $E_{\text{sym}}(\rho)\approx 31.6(\rho /\rho_{0})^{1.05} $ at subnormal densities. This leads to a significantly constrained value of about -550 MeV for the isospin-dependent part of the isobaric incompressibility of isospin asymmetric nuclear matter.

We present the results of linear analysis and two-dimensional local magnetohydrodynamic (MHD) simulations of the Parker instability, including the effects of cosmic rays (CRs), in magnetized gas disks (galactic disks). As an unperturbed state for both the linear analysis and the MHD simulations, we adopted an equilibrium model of a magnetized two-temperature layered disk with constant gravitational acceleration parallel to the normal of the disk. The disk comprises a thermal gas, cosmic rays and a magnetic field perpendicular to the gravitational accelerartion. Cosmic ray diffusion along the magnetic field is considered; cross field-line diffusion is supposed to be small and is ignored. We investigated two cases in our simulations. In the mechanical perturbation case we add a velocity perturbation parallel to the magnetic field lines, while in the explosional perturbation case we add cosmic ray energy into a sphere where the cosmic rays are injected. Linear analysis shows that the growth rate of the Parker instability becomes smaller if the coupling between the CR and the gas is stronger (i.e., the CR diffusion coefficient is smaller). Our MHD simulations of the mechanical perturbation confirm this result. We show that the falling matter is impeded by the CR pressure gradient, this causes a decrease in the growth rate. In the explosional perturbation case, the growth of the magnetic loop is faster when the coupling is stronger in the early stage. However, in the later stage the behavior of the growth rate becomes similar to the mechanical perturbation case.

The coalescence model based on nucleon distribution functions from an isospin-dependent transport model is used to study the production of light clusters such as deuteron, triton, and $^{3}$He from heavy-ion collisions induced by neutron-rich nuclei at intermediate energies. It is found that the emission time of light clusters depends on their masses. For clusters with the same momentum per nucleon, heavier ones are emitted earlier. Both the yield and energy spectrum of light clusters are sensitive to the density dependence of nuclear symmetry energy, with more light clusters produced in the case of a stiff symmetry energy. On the other hand, effects due to the stiffness of the isoscalar part of nuclear equation of state and the medium dependence of nucleon-nucleon cross sections on light cluster production are unimportant. We have also studied the correlation functions of clusters, and they are affected by the density dependence of nuclear symmetry energy as well, with the stiff symmetry energy giving a stronger anti-correlation of light clusters, particularly for those with large kinetic energies. Dependence of light cluster production on the centrality and incident energy of heavy ion collisions as well as the mass of the reaction system is also investigated.

Using an isospin-dependent transport model for heavy-ion collisions induced by neutron-rich nuclei at intermediate energies, we study the production of light clusters such as deuteron, triton, and $^{3}$He via coalescence of nucleons. We find that both the yields and energy spectra of these light clusters are affected significantly by the density dependence of nuclear symmetry energy, with a stiffer symmetry energy giving a larger yield.

Using an isospin-dependent transport model, we study isospin effects on two-nucleon correlation functions in heavy-ion collisions induced by neutron-rich nuclei at intermediate energies. We find that these correlation functions are sensitive to the density dependence of nuclear symmetry energy, but not to the incompressibility of symmetric nuclear matter and the medium dependence of nucleon-nucleon cross sections. This sensitivity is mainly due to effects of nuclear symmetry energy on the emission times of neutrons and protons as well as their relative emission sequence. We also study the variations of the symmetry energy effects on nucleon-nucleon correlations with respect to the impact parameter, incident energy, and mass number of heavy ion collisions.

Using an isospin-dependent transport model, we study the effects of nuclear symmetry energy on two-nucleon correlation functions in heavy ion collisions induced by neutron-rich nuclei. We find that the density dependence of the nuclear symmetry energy affects significantly the nucleon emission times in these collisions, leading to larger values of two-nucleon correlation functions for a symmetry energy that has a stronger density dependence. Two-nucleon correlation functions are thus useful tools for extracting information about the nuclear symmetry energy from heavy ion collisions.