By comparing cumulative linear mass profiles of 12 Gould Belt molecular clouds within 500 pc, we study how the linear mass distributions of molecular clouds vary with the angles between the molecular cloud long axes and the directions of the local magnetic fields (cloud-field direction offsets). We find that molecular clouds with the long axes perpendicular to the magnetic field directions show more even distributions of the linear mass. The result supports that magnetic field orientations can affect the fragmentation of molecular clouds (Li et al. 2017).

The evolution of a large-scale poloidal magnetic field in an accretion disc is an important problem because it determines the launching of winds and the feasibility of the magnetorotational instability to generate turbulence or channel flows. Recent studies, both semi-analytical calculations and numerical simulations, have highlighted the crucial role non-ideal MHD effects (Ohmic resistivity, Hall drift and ambipolar diffusion), relevant in the protoplanetary disc context, might play in magnetic flux evolution in the disc. We investigate the flux transport in discs through the use of two one-dimensional semi-analytic models in the vertical direction, exploring regimes where different physical source terms and effects dominate. The governing equations for both models are derived by performing an asymptotic expansion in the limit of a thin disc, with the different regimes isolated through setting the relative order of the leading terms between variables. Flux transport rates and vertical structure profiles are calculated for a range of diffusivities and disc magnetisations. We found that Ohmic and ambipolar diffusivities drive radially outward flux transport with an outwardly inclined field. A wind outflow drives inward flux transport, which is significantly enhanced in the presence of Hall drift in the positive polarity case, $\eta_H (\boldsymbol{B}_z \cdot \boldsymbol{\Omega}) > 0$, an effect which has only been briefly noted before. Coupled only with outward inclination, the Hall effect reduces the flux transport given by a background Ohmic and/or ambipolar diffusivity, but drives no flux transport when it is the only non-ideal effect present.

Observations of synchrotron polarization at multiple frequencies in the presence of Faraday rotation can provide a way to reconstruct the 3D magnetic field distribution. In this paper we compare the well known Faraday Tomography (FT, Burn 1966) technique to the new approach named Synchrotron Polarization Gradients (SPG, Lazarian & Yuen 2018b). We compare the strengths and limitations of the two techniques, and describe their synergy. In particular, we show that in the situations when FT technique fails, e.g. due to insufficient frequency coverage, the SPG can still trace the 3D structure of magnetic field.

(Abridged) It is well known that the tidal deformability of a compact star carries important information about the interior equation-of-state (EOS) of the star. The first gravitational-wave event GW170817 from a binary compact star merger observed by the LIGO/VIRGO detectors have already put limits on the tidal deformability and provided constraints on the ultra-high nuclear density EOS. In view of this ground breaking discovery, we revisit and extend our previous work [Phys. Rev. D 95, 101302(R) (2017)] which found that taking the effect of elasticity into account in the calculation of the tidal deformability of compact star models composed of crystalline color-superconducting (CCS) quark matter can break the universal I-Love relation discovered for fluid compact stars. In this paper, we present our formulation in detail and provide more analysis to complement our previous findings. We focus and extend the study of the screening effect on the tidal deformability, which we found previously for hybrid star models, to various theoretical two-layer compact star models. We show that the screening effect of these two-layer models in general depends on the thickness of the envelope and the ratio between the density gap and the core density at the core-envelope interface. However, for models with a fluid envelope and a vanishing small density gap, the screening effect remains strong even as the thickness of the envelops tends to zero if the quark matter core has a fairly uniform density. The relevance of our study to GW170817 is also discussed. Our study advocates that the tidal deformability not only provides us information on the EOS, but may also give insights into the multi-layer structure and elastic properties of compact star models composed of CCS quark matter.

The physical significance of tidal deformation in astronomical systems has long been known. The recently discovered universal I-Love-Q relations, which connect moment of inertia, quadrupole tidal Love number, and spin-induced quadrupole moment of compact stars, also underscore the special role of tidal deformation in gravitational wave astronomy. Motivated by the observation that such relations also prevail in Newtonian stars and crucially depend on the stiffness of a star, we consider the tidal Love numbers of Newtonian polytropic stars whose stiffness is characterised by a polytropic index $n$. We first perturbatively solve the Lane-Emden equation governing the profile of polytropic stars through the application of the scaled delta expansion method and then formulate perturbation series for the multipolar tidal Love number about the two exactly solvable cases with $n=0$ and $n=1$, respectively. Making use of these two series to form a two-point Pad\'e approximant, we find an approximate expression of the quadrupole tidal Love number, whose error is less than $2.5 \times 10^{-5}$ per cent (0.39 per cent) for $n\in[0,1]$ ($n\in[0,3]$). Similarly, we also determine the mass moments for polytropic stars accurately. Based on these findings, we are able to show that the I-Love-Q relations are in general stationary about the incompressible limit irrespective of the equation of state (EOS) of a star. Moreover, for the I-Love-Q relations, there is a secondary stationary point near $n \approx 0.4444$, thus showing the insensitivity to $n$ for $n\in[0,1]$. Our investigation clearly tracks the universality of the I-Love-Q relations from their validity for stiff stars such as neutron stars to their breakdown for soft stars.

We study the tidal deformability of bare quark stars and hybrid compact stars composed of a quark matter core in general relativity, assuming that the deconfined quark matter exists in a crystalline color superconducting phase. We find that taking the elastic property of crystalline quark matter into account in the calculation of the tidal deformability can break the universal I-Love relation discovered for fluid compact stars, which connects the moment of inertia and tidal deformability. Our result suggests that measurements of the moment of inertia and tidal deformability can in principle be used to test the existence of solid quark stars, despite our ignorance of the high density equation of state (EOS). Assuming that the moment of inertia can be measured to 10% level, one can then distinguish a 1.4 (1) $M_\odot$ solid quark star described by our quark matter EOS model with a gap parameter $\Delta=25$ MeV from a fluid compact star if the tidal deformability can be measured to about 10% (45%) level. On the other hand, we find that the nuclear matter fluid envelope of a hybrid star can screen out the effect of the solid core significantly so that the resulting I-Love relation for hybrid stars still agrees with the universal relation for fluid stars to about 1% level.

Under suitable scaling, the structure of self-gravitating polytropes is described by the standard Lane-Emden equation (LEE), which is characterised by the polytropic index $n$. Here we use the known exact solutions of the LEE at $n=0$ and $1$ to solve the equation perturbatively. We first introduce a scaled LEE (SLEE) where polytropes with different polytropic indices all share a common scaled radius. The SLEE is then solved perturbatively as an eigenvalue problem. Analytical approximants of the polytrope function, the radius and the mass of polytropes as a function of $n$ are derived. The approximant of the polytrope function is well-defined and uniformly accurate from the origin down to the surface of a polytrope. The percentage errors of the radius and the mass are bounded by $8.1 \times 10^{-7}$ per cent and $8.5 \times 10^{-5}$ per cent, respectively, for $n\in[0,1]$. Even for $n\in[1,5)$, both percentage errors are still less than $2$ per cent.

The emergence of the I-Love-Q relations, revealing that the moment of inertia, the tidal Love number (deformability) and the spin-induced quadrupole moment of compact stars are, to high accuracy, interconnected in a universal way disregarding the wide variety of equations of state (EOSs) of dense matter, has attracted much interest recently. However, the physical origin of these relations is still a debatable issue. In the present paper, we focus on the I-Love relation for self-bound stars (SBSs) such as incompressible stars and quark stars. We formulate perturbative expansions for the moment of inertia, the tidal Love number (deformability) and the I-Love relation of SBSs. By comparing the respective I-Love relations of incompressible stars and a specific kind of SBSs, we show analytically that the I-Love relation is, to relevant leading orders in stellar compactness, stationary with respect to changes in the EOS about the incompressible limit. Hence, the universality of the I-Love relation is indeed attributable to the proximity of compact stars to incompressible stars, and the stationarity of the relation as unveiled here. We also discover that the moment of inertia and the tidal deformability of a SBS with finite compressibility are, to leading order in compactness, equal to their counterparts of an incompressible star with an adjusted compactness, thus leading to a novel explanation for the I-Love universal relation.

Most molecular clouds are filamentary or elongated. Among those forming low-mass stars, their long axes tend to be either parallel or perpendicular to the large-scale (10-100 pc) magnetic field (B-field) in the surrounding inter cloud medium. This arises because, along the dynamically dominant B-fields, the competition between self-gravity and turbulent pressure will shape the cloud to be elongated either perpendicular or parallel to the fields. Recent study also suggested that, on the scales of 0.1-0.01 pc, fields are dynamically important within cloud cores forming massive stars. But whether the core field morphologies are inherited from the inter cloud medium or governed by cloud turbulence is under vigorous debate, so is the role played by B-fields in cloud fragmentation at 10 - 0.1 pc scales. Here we report B-field maps covering 100-0.01 pc scales inferred from polarimetric observations of a massive-star forming region, NGC 6334. First, the main filament also lies perpendicular to the ambient field. NGC 6334 hosts young star-forming sites where fields are not severely affected by stellar feedback, and their directions do not change significantly over the entire scale range. This means that the fields are dynamically important. At various scales, we find that the hourglass-shaped field lines are pinched where the gas column density peaks and the field strength is proportional to the 0.4-power of the density. We conclude that B-fields play a crucial role in the fragmentation of NGC 6334.

Though individual stellar parameters of compact stars usually demonstrate obvious dependence on the equation of state (EOS), EOS-insensitive universal formulas relating these parameters remarkably exist. In the present paper, we explore the interrelationship between two such formulas, namely the $f$-$I$ relation connecting the $f$-mode quadrupole oscillation frequency $\omega_2$ and the moment of inertia $I$, and the $I$-Love-$Q$ relations relating $I$, the quadrupole tidal deformability $\lambda_2$, and the quadrupole moment $Q$, which have been proposed by Lau, Leung, and Lin [Astrophys. J. {\bf 714}, 1234 (2010)] and Yagi and Yunes [Science {\bf 341}, 365 (2013)], respectively. A relativistic universal relation between $\omega_l$ and $\lambda_l$ with the same angular momentum $l=2,3,\ldots$, the so-called "diagonal $f$-Love relation" that holds for realistic compact stars and stiff polytropic stars, is unveiled here. An in-depth investigation in the Newtonian limit is further carried out to pinpoint its underlying physical mechanism and hence leads to a unified $f$-$I$-Love relation. We reach the conclusion that these EOS-insensitive formulas stem from a common physical origin --- compact stars can be considered as quasiincompressible when they react to slow time variations introduced by $f$-mode oscillations, tidal forces and rotations.

The recent discovery of the universal I-Love-Q relations connecting the moment of inertia, tidal deformability, and the spin-induced quadrupole moment of compact stars is intriguing and totally unexpected. In this paper, we provide numerical evidence showing that the universality can be attributed to the incompressible limit of the I-Love-Q relations. The fact that modern equations of state are stiff, with an effective adiabatic index larger than about two, above the nuclear density range is the key to establishing the universality for neutron stars and quark stars with typical compactness from about 0.1 to 0.3. On the other hand, the I-Love-Q relations of low-mass neutron stars near the minimum mass limit depend more sensitively on the underlying equation of state because these stars are composed mainly of softer matter at low densities. However, the I-Love-Q relations for low-mass quark stars can still be represented accurately by the incompressible limit. We also study the I-Love relation connecting the moment of inertia and tidal deformability analytically in Newtonian gravity and show why the I-Love-Q relation is weakly dependent on the underlying equation of state and can be attributed to its incompressible limit.

In spite of the diversity in the equations of state of nuclear matter, the recently discovered I-Love-Q relations [Yagi and Yunes, Science {\bf 341}, 365 (2013)], which relate the moment of inertia, tidal Love number (deformability) and the spin-induced quadrupole moment of compact stars, hold for various kinds of realistic neutron stars and quark stars. While the physical origin of such universality is still a current issue, the observation that the I-Love-Q relations of incompressible stars can well approximate those of realistic compact stars hints at a new direction to approach the problem. In this paper, by establishing recursive post-Minkowskian expansion for the moment of inertia and the tidal deformability of incompressible stars, we analytically derive the I-Love relation for incompressible stars and show that the so obtained formula can be used to accurately predict the behavior of realistic compact stars from the Newtonian limit to the maximum mass limit.

Due to our ignorance of the equation of state (EOS) beyond nuclear density, there is still no unique theoretical model for neutron stars (NSs). It is therefore surprising that universal EOS-independent relations connecting different physical quantities of neutron stars can exist. Lau et al. [ApJ, 714, 1234 (2010)] found that the frequency of the $f$-mode oscillation, the mass, and the moment of inertia are connected by universal relations. More recently, Yagi and Yunes [Science, 341, 365 (2013)] discovered the I-Love-Q universal relations among the mass, the moment of inertia, the Love number, and the quadrupole moment. In this paper, we study these universal relations in the Eddington-inspired Born-Infeld (EiBI) gravity. This theory differs from general relativity (GR) significantly only at high densities due to the nonlinear coupling between matter and gravity. It thus provides us an ideal case to test how robust the universal relations of NSs are with respect to the change of the gravity theory. Thanks to the apparent EOS formulation of EiBI gravity developed recently by Delsate and Steinhoff [Phys. Rev. Lett., 109, 021101 (2012)], we are able to study the universal relations in EiBI gravity using the same techniques as those in GR. We find that the universal relations in EiBI gravity are essentially the same as those in GR. Our work shows that, within the currently viable coupling constant, there exists at least one modified gravity theory that is indistinguishable from GR in view of the unexpected universal relations.

We study how generic phase transitions taking place in compact stars constructed in the framework of the Eddington-inspired Born-Infeld (EiBI) gravity can lead to anomalous behavior of these stars. For the case with first-order phase transitions, compact stars in EiBI gravity with a positive coupling parameter $\kappa $ exhibit a finite region with constant pressure, which is absent in general relativity. However, for the case with a negative $\kappa $, an equilibrium stellar configuration cannot be constructed. Hence, EiBI gravity seems to impose stricter constraints on the microphysics of stellar matter. Besides, in the presence of spatial discontinuities in the sound speed $c_s$ due to phase transitions, the Ricci scalar is spatially discontinuous and contains $\delta$-function singularities proportional to the jump in $c_s^2$ acquired in the associated phase transition.

We study the hydrostatic equilibrium structure of compact stars in the Eddington-inspired Born-Infeld gravity recently proposed by Banados and Ferreira [Phys. Rev. Lett. 105, 011101 (2010)]. We also develop a framework to study the radial perturbations and stability of compact stars in this theory. We find that the standard results of stellar stability still hold in this theory. The frequency square of the fundamental oscillation mode vanishes for the maximum-mass stellar configuration. The dependence of the oscillation mode frequencies on the coupling parameter \kappa of the theory is also investigated. We find that the fundamental mode is insensitive to the value of \kappa, while higher order modes depend more strongly on \kappa.

We propose here a robust scheme to infer the physical parameters of compact stars from their f-mode gravitational wave signals. We first show that the frequency and the damping rate of f-mode oscillation of compact stars can be expressed in terms of universal functions of stellar mass and moment of inertia. By employing the universality in the f-mode one can then infer accurate values of the mass, the moment of inertia and the radius of a compact star. In addition, we demonstrate that our new scheme works well for both realistic neutron stars and quark stars, and hence provides a unifying way to infer the physical parameters of compact stars.

We review results from general relativistic axisymmetric magnetohydrodynamic simulations of accretion in Sgr A*. We use general relativistic radiative transfer methods and to produce a broad band (from millimeter to gamma-rays) spectrum. Using a ray tracing scheme we also model images of Sgr A* and compare the size of image to the VLBI observations at 230 GHz. We perform a parameter survey and study radiative properties of the flow models for various black hole spins, ion to electron temperature ratios, and inclinations. We scale our models to reconstruct the flux and the spectral slope around 230 GHz. The combination of Monte Carlo spectral energy distribution calculations and 230 GHz image modeling constrains the parameter space of the numerical models. Our models suggest rather high black hole spin ($a_*\approx 0.9$), electron temperatures close to the ion temperature ($T_i/T_e \sim 3$) and high inclination angles ($i \approx 90 \deg$).

Using flow models based on axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations, we construct radiative models for sgr A*. Spectral energy distributions that include the effects of thermal synchrotron emission and absorption, and Compton scattering, are calculated using a Monte Carlo technique. Images are calculated using a ray-tracing scheme. All models are scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters are the dimensionless black hole spin a*, the inclination i, and the ion-to-electron temperature ratio Ti/Te. We find that: (1) models with Ti/Te=1 are inconsistent with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly increasing function of a*; (3) the X-ray flux is a strongly increasing function of i; (4) 230 GHz image size is a complicated function of i, a*, and Ti/Te, but the Ti/Te = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with Ti/Te=10 and i=85 deg the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) in all models whose SEDs are consistent with observations the event horizon is uncloaked at 230 GHz; (7) the models that are most consistent with the observations have a* \sim 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

It has been well established that the $f$-mode of relativistic ordinary-fluid neutron stars displays a universal scaling behavior. Here we study whether the "ordinary" $f_{\rm o}$- and "superfluid" $f_{\rm s}$-modes of superfluid neutron stars also show similar universal behavior. We first consider a simple case where the neutron superfluid and normal fluid are decoupled, and with each fluid modeled by a polytropic equation of state. We find that the $f_{\rm o}$-mode obeys the same scaling laws as established for the $f$-mode of orindary-fluid stars. However, the oscillation frequency of the $f_{\rm s}$-mode obeys a different scaling law, which can be derived analytically from a homogenous two-fluid stellar model in Newtonian gravity. Next the coupling effect between the two fluids is studied via a parameterized model of entrainment. We find that the coupling in general breaks the universal behavior seen in the case of decoupled fluids. Based on a relativistic variational principle, an approximated expression is derived for the first-order shift of the $f_{\rm s}$-mode squared frequency due to the entrainment.

During fall periods in 2002, 2003 and 2004 three major oceanographic expeditions were carried out in Mamala Bay, Hawaii. These were part of the RASP Remote Anthropogenic Sensing Program. Ikonos and Quickbird optical satellite images of sea surface glint revealed ~100 m spectral anomalies in km^2 averaging patches in regions leading from the Honolulu Sand Island Municipal Outfall diffuser to distances up to 20 km. To determine the mechanisms behind this phenomenon, the RASP expeditions monitored the waters adjacent to the outfall with an array of hydrographic, optical and turbulence microstructure sensors in anomaly and ambient background regions. Drogue tracks and mean turbulence parameters for 2x10^4 microstructure patches were analyzed to understand complex turbulence, fossil turbulence and zombie turbulence near-vertical internal wave transport processes. The dominant mechanism appears to be generic to stratified natural fluids including planet and star atmospheres and is termed beamed zombie turbulence maser action (BZTMA). Most of the bottom turbulent kinetic energy is converted to ~100 m fossil turbulence waves. These activate secondary (zombie) turbulence in outfall fossil turbulence patches that transmit heat, mass, chemical species, momentum and information vertically to the sea surface for detection in an efficient maser action. The transport is beamed in intermittent mixing chimneys.

Adopting the Lindblom-Detweiler formalism for polar oscillations of neutron stars, we study the $w$-mode oscillation and find that the Lagrangian change in pressure, measured by the physical quantity $X$, is negligibly small. Based on this observation, we develop the generalised inverse-Cowling approximation (GICA) with the approximation $X=X'=0$, where $X'$ is the derivative of $X$ with respect to the circumferential radius, for $w$-mode oscillations of neutron stars. Under GICA, $w$-mode oscillations are described by a second-order differential system, which can yield accurate frequencies and damping rates of quasi-normal modes.

The radio source Sagittarius A* (Sgr A*) is believed to be a hot, inhomogeneous, magnetized plasma flowing near the event horizon of the 3 million solar mass black hole at the galactic center. At a distance of 8000 parsecs the black hole would be among the largest black holes as judged by angular size. Recent observations are consistent with the idea that the millimeter and sub-millimeter photons are dominated by optically thin, thermal synchrotron emission. Anticipating future Very Long Baseline Interferometry (VLBI) observations of Sgr A* at these wavelengths, we present here the first dynamically self-consistent models of millimeter and sub-millimeter emission from Sgr A* based on general relativistic numerical simulations of the accretion flow. Angle-dependent spectra are calculated assuming a thermal distribution of electrons at the baryonic temperature dictated by the simulation and the accretion rate, which acts as a free parameter in our model. The effects of varying model parameters (black hole spin and inclination of the spin to the line of sight) and source variability on the spectrum are shown. We find that the accretion rate value needed to match our calculated millimeter flux to the observed flux is consistent with constraints on the accretion rate inferred from detections of the rotation measure. We also describe the relativistic jet that is launched along the black hole spin axis by the accretion disk and evolves to scales of 1000 gravitational radii.

We show here how the internal structure of a neutron star can be inferred from its gravitational wave spectrum. Under the premise that the frequencies and damping rates of a few $w$-mode oscillations are found, we apply an inversion scheme to determine its mass, radius and density distribution. In addition, accurate equation of state of nuclear matter can also be determined.

We study the $\gamma$-ray emission properties of pulsars by using a new self-consistent outer gap model. The outer gap can exist in pulsars with age over million years old if the effect of magnetic inclination angle as well as the average properties of the outer gap are considered. The mature $\gamma$-ray pulsars, whose ages are between 0.3 to 3 million years old, are able to move up to high galactic latitude. Moreover, their $\gamma$-ray luminosity are weaker and their spectra are softer than those of younger $\gamma$-ray pulsars in the galactic plane significantly. We use a Monte Carlo method to simulate the statistical properties of $\gamma$-ray pulsars in the Galaxy as well as in the Gould Belt. We find that $\gamma$-ray pulsars located at $\mid b \mid < 5^{\circ}$ and located at $\mid b \mid > 5^{\circ}$ have very different properties. High galactic latitude $\gamma$-ray pulsars are dominated by mature pulsars with longer periods, weaker fluxes and softer spectra. If the pulsar birth rate in the Galaxy and the Gould Belt are $\sim 10^{-2}yr^{-1}$ and $\sim 2\times 10^{-5}yr^{-1}$ respectively, there are 42 and 35 radio-quiet $\gamma$-ray pulsars for $\mid b \mid < 5^{\circ}$ and $\mid b \mid > 5^{\circ}$ respectively. Radio-quiet $\gamma$-ray pulsars from the Gould Belt are 2 and 13 for $\mid b \mid < 5^{\circ}$ and $\mid b \mid > 5^{\circ}$ respectively. We suggest that a good fraction of unidentified EGRET $\gamma$-ray sources may be these radio-quiet $\gamma$-ray pulsars. Furthermore $\gamma$-ray pulsars located at $\mid b \mid > 5^{\circ}$ satisfies $L_{\gamma} \propto L^{\beta}_{sd}$ whereas $L_{\gamma} \propto L^{\delta}_{sd}$ for $\gamma$-ray pulsars in the galactic plane, where $\beta \sim 0.6$ and $\delta \sim 0.3$ respectively.

We re-examine the outer gap size by taking the geometry of the dipole magnetic field into account. Furthermore, we also consider that instead of taking the gap size at half of the light cylinder radius to represent the entire outer gap it is more appropriate to average the entire outer gap size over the distance. When these two factors are considered, the derived outer gap size f is not only the function of period P and magnetic field B of the neutron star, but also the function of the average radial distance to the neutron star <r>. We use this new outer gap model to study $\gamma$-ray luminosity of pulsars, and to study the death lines of $\gamma$-ray emission of the pulsars. Our model can predict the $\gamma$-ray luminosity of individual pulsar if its $P, B$ and $\alpha$ are known. Since different pulsars have different $\alpha$, this explains why some $\gamma$-ray pulsars have very similar $P$ and $B$ but have very different $\gamma$-ray luminosities. In determining the death line of $\gamma$-ray pulsars, we have used a new criterion based on concrete physical reason. In estimate of the fractional size of the outer gap, two possible X-ray fields are considered: (i) X-rays are produced by the neutron star cooling and polar cap heating, and (ii)X-rays are produced by the bombardment of the relativistic particles from the outer gap on the stellar surface. Since it is very difficult to measure $\alpha$ in general, we use a Monte Carlo method to simulate the properties of $\gamma$-ray pulsars in our galaxy. We find that this new outer gap model predicts many more weak $\gamma$-ray pulsars. For all simulated $\gamma$-ray pulsars with self-sustained outer gaps, $\gamma$-ray luminosity $L_{\gamma}$ satisfies $L_{\gamma}\propto L^{\delta}_{sd}$; where the value of $\delta$ depends on the ensitivity of the $\gamma$-ray detector.

The atomic recombination process leads to a softening of the matter equation of state as reflected by a reduced generalized adiabatic index, with accompanying heat release. We study the effects of this recombination softening and reheating of the cosmic plasma on the ionization history, visibility function, Cold Dark Matter (CDM) transfer function, and the Cosmic Microwave Background (CMB) spectra. The resulting modifications of the CMB spectrm is 1/10 of WMAP's current error and is comparable to PLANCK's error. Therefore, this effect should be considered when data with higher accuracy are analysed.

Turbulence is defined as an eddy-like state of fluid motion where the inertial-vortex forces of the eddies are larger than any of the other forces that tend to damp the eddies out. Energy cascades of irrotational flows from large scales to small are non-turbulent, even if they supply energy to turbulence. Turbulent flows are rotational and cascade from small scales to large, with feedback. Viscous forces limit the smallest turbulent eddy size to the Kolmogorov scale. In stratified fluids, buoyancy forces limit large vertical overturns to the Ozmidov scale and convert the largest turbulent eddies into a unique class of saturated, non-propagating, internal waves, termed fossil-vorticity-turbulence. These waves have the same energy but different properties and spectral forms than the original turbulence patch. The Gibson (1980, 1986) theory of fossil turbulence applies universal similarity theories of turbulence and turbulent mixing to the vertical evolution of an isolated patch of turbulence in a stratified fluid as its growth is constrained and fossilized by buoyancy forces. These theories apply to the dynamics of atmospheric, astrophysical and cosmological turbulence.