The diffuse plasma inside clusters of galaxies has X-ray emitting temperatures of a few keV. The physical mechanisms that heat this intracluster medium (ICM) to such temperatures include the accretion shock at the periphery of a galaxy cluster, the shocks driven by merger events, as well as a somewhat overlooked mechanism -- the dissipation of intracluster turbulent motions. We study the relative role of these heating mechanisms using galaxy clusters in Lagrangian tracer particle re-simulations of the Omega500 cosmological simulation. We adopt a novel analysis method of decomposing the temperature increase at each time step into the contribution from dissipative heating and that from adiabatic heating. In the high-resolution spatial-temporal map of these heating rates, merger tracks are clearly visible, demonstrating the dominant role of merger events in heating the ICM. The dissipative heating contributed by each merger event is extended in time and also occurs in the rarefaction regions, suggesting the importance of heating by the dissipation of merger-induced turbulence. Quantitative analysis shows that turbulence heating, rather than direct heating at merger shocks, dominates the temperature increase of the ICM, especially at inner radii $r < r_{\rm 500c}$. In addition, we find that many merger shocks can propagate with almost constant velocity to very large radii $r \gg r_{\rm 500c}$, some even reach and join with the accretion shock and becoming the outer boundary of the ICM. Altogether, these results suggest that the ICM is heated more in an `inside-out' fashion rather than `outside-in' as depicted in the classical smooth accretion picture.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

The co-evolution between black holes and galaxies suggests that feedback of active galactic nuclei influence host galaxies through ejecting radiative and kinetic energies to surroundings. Larger scale outflow in local universe are frequently observed by spatially resolved spectroscopy, while smaller scale outflow cannot be directly resolved by current observations. At the scale of the dusty torus, radiative and kinetic energies ejected from the central active nucleus interact with the materials. However, observations of such outflow are rarely reported due to the lack detection of unambiguously gas emission. Here we report the detection of clear and rich emission lines origin from the scale of dusty tours in an partially obscured quasar. The lines share a common intermediate width with full width at half maximum about 1900 \kmps\ and are shown in two systems: a major system is unshifted and a minor system has a blue-shifts of about 2600 \kmps. The line intensity ratios, combining photo-ionization simulations, indicates an ultradense line-emitting region with the density as high as $\sim$ $10^{13}~\rm cm^{-3}$. We interpret this as the lines being excited by a shock induced by the high-density and high-temperature gases at the scale of dusty torus, rather than photo-ionized by the central accretion disk. We speculate that the outflow, launched from the accretion disk, collides onto the inner wall of the dusty torus and shock-heat the gases to cause the major emission lines. The outflowing gases may also collide onto surrounding isolated clouds, and give rise to blue-shifted minor emission lines.

ASAS J174406+2446.8 was originally found as a $\delta$ Scuti-type pulsating star with the period P=0.189068 $days$ by ASAS survey. However, the LAMOST stellar parameters reveal that it is far beyond the red edge of pulsational instability strip on the $\log g-T$ diagram of $\delta$ Scuti pulsating stars. To understand the physical properties of the variable star, we observed it by the 1.0-m Cassegrain reflecting telescope at Yunnan Observatories. Multi-color light curves in B, V, R$_{c}$ and I$_{c}$ bands were obtained and are analyzed by using the W-D program. It is found that this variable star is a shallow-contact binary with an EB-type light curve and an orbital period of 0.3781\,days rather than a $\delta$ Scuti star. It is a W-subtype contact binary with a mass ratio of $1.135(\pm0.019)$ and a fill-out factor of $10.4(\pm5.6)\,\%$. The situation of ASAS J174406+2446.8 resembles those of other EB-type marginal-contact binaries such as UU Lyn, II Per and GW Tau. All of them are at a key evolutionary phase from a semi-detached configuration to a contact system predicted by the thermal relaxation oscillation theory. The linear ephemeris was corrected by using 303 new determined times of light minimum. It is detected that the O - C curve shows a sinusoidal variationthat could be explained by the light-travel-time effect via the presence of a cool red dwarf. The present investigation reveals that some of the $\delta$ Scuti-type stars beyond the red edge of pulsating instability strip on the $\log g-T$ diagram are misclassified eclipsing binaries. To understand their structures and evolutionary states, more studies are required in the future.

The goal of Project GAUSS is to return samples from the dwarf planet Ceres. Ceres is the most accessible ocean world candidate and the largest reservoir of water in the inner solar system. It shows active cryovolcanism and hydrothermal activities in recent history that resulted in minerals not found in any other planets to date except for Earth's upper crust. The possible occurrence of recent subsurface ocean on Ceres and the complex geochemistry suggest possible past habitability and even the potential for ongoing habitability. Aiming to answer a broad spectrum of questions about the origin and evolution of Ceres and its potential habitability, GAUSS will return samples from this possible ocean world for the first time. The project will address the following top-level scientific questions: 1) What is the origin of Ceres and the origin and transfer of water and other volatiles in the inner solar system? 2) What are the physical properties and internal structure of Ceres? What do they tell us about the evolutionary and aqueous alteration history of icy dwarf planets? 3) What are the astrobiological implications of Ceres? Was it habitable in the past and is it still today? 4) What are the mineralogical connections between Ceres and our current collections of primitive meteorites? GAUSS will first perform a high-resolution global remote sensing investigation, characterizing the geophysical and geochemical properties of Ceres. Candidate sampling sites will then be identified, and observation campaigns will be run for an in-depth assessment of the candidate sites. Once the sampling site is selected, a lander will be deployed on the surface to collect samples and return them to Earth in cryogenic conditions that preserves the volatile and organic composition as well as the original physical status as much as possible.

IRAS F11119+3257 is a quasar-dominated Ultra-Luminous InfraRed Galaxy, with a partially obscured narrow-line seyfert 1 nucleus. In this paper, we present the NIR spectroscopy of F11119+3257, in which we find unusual Paschen emission lines, and metastable He I* $\lambda$10830 absorption associated with the previously reported atomic sodium and molecular OH mini-BAL (Broad Absorption Line) outflow. Photo-ionization diagnosis confirms previous findings that the outflows are at kilo-parsec scales. Such large-scale outflows should produce emission lines. We indeed find that high-ionization emission lines ([O III], [Ne III], and [Ne V]) are dominated by blueshifted components at similar speeds to the mini-BALs. The blueshifted components are also detected in some low-ionization emission lines, such as [O II] $\lambda$3727 and some Balmer lines (H$\alpha$, H$\beta$, and H$\gamma$), even though their cores are dominated by narrow ($FWHM_{\rm NEL} = 570\pm40$km s$^{-1}$) or broad components at the systemic redshift of $z=0.18966\pm0.00006$. The mass flow rate (230-730$~M_\odot \rm yr^{-1}$) and the kinetic luminosity ($\dot{E}_k \sim 10^{43.6-44.8} $erg s$^{-1}$) are then inferred jointly from the blueshifted emission and absorption lines. In the NIR spectrum of F11119+3257, we also find that the Paschen emission lines are unique, in which a very narrow ($FWHM=260\pm20~$km s$^{-1}$) component is shown in only Pa$\alpha$. This narrow component most probably comes from heavily obscured star formation. Based on the Pa$\alpha$ and Pa$\beta$ emissions, we obtain an extinction at the $H$ band, $A_H~>~2.1$ (or a reddenning of $E_{B-V}~>~$3.7), and a star formation rate of $SFR~>~130\rm M_\odot yr^{-1}$ that resembles the estimates inferred from the FIR emissions ($SFR_{\rm FIR} = 190\pm90$ $M_\odot$ yr$^{-1}$).

We extract the multiple-epoch Balmer-line profiles of the heavily obscured quasar 3C 445 from the spectral curves in the literature, and analyze the emission-line profiles of the H$\alpha$ and H$\beta$ lines and the profile variability in the H$\alpha$ line in the large time interval of more than three decades. The profile comparison between the H$\alpha$ and H$\beta$ lines shows that both Balmer lines share the profile with the same form, while the blue system of the H$\beta$ line is seriously weaker than that of the H$\alpha$ line. Moreover, the blue system of the H$\alpha$ line suddenly disappeared completely and then did not appear again, however the other two components did not exhibit significant variation in the velocity or the amplitude. These findings suggest that the blue system of 3C 445, as with SDSS J153636.22+044127.0 and its analogs, is probable the result of the shock-heated outflowing gases. The observation angle of almost edge-on which the previous studies suggested can easily produce the high-speed and high-temperature shock in the collision between the massive outflow and the inner surface of the dusty torus.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

Proton-proton interaction and photo-hadronic interaction in cosmic accelerators are the two main channels for the production of cosmic ultra-high energy photons and neutrinos (TeV-PeV). In this Letter, we use FWW approach to obtain the production of cosmic ultra-high energy photons and neutrinos from the collision between UHE proton with magnetic field which could be considered as the virtualphoton in the rest frame of UHE proton. We name this as $pB$ process. The threshold for the occurrence of the $pB$ process is that the combination of the Lorentz factor of proton and the strength of the magnetic field is about $\gamma_p B \simeq 5\times 10^{18}$Gauss. Beyond this threshold, the rate of energy loss of proton due to the $pB$ process is about three orders higher than that due to the synchrotron radiation of proton in the same magnetic field. The $pB$ process might potentially happen in the atmosphere of white dwarfs, neutron stars or even that of stellar massive black holes.

Turbulence evolution in a density-stratified medium differs from that of homogeneous isotropic turbulence described by the Kolmogorov picture. We evaluate the degree of this effect in the intracluster medium (ICM) with hydrodynamical simulations. We find that the buoyancy effect induced by ICM density stratification introduces qualitative changes to the turbulence energy evolution, morphology, and the density fluctuation - turbulence Mach number relation, and likely explains the radial dependence of the ICM turbulence amplitude as found previously in cosmological simulations. A new channel of energy flow between the kinetic and the potential energy is opened up by buoyancy. When the gravitational potential is kept constant with time, this energy flow leaves oscillations to the energy evolution, and leads to a balanced state of the two energies where both asymptote to power-law time evolution with slopes shallower than that for the turbulence kinetic energy of homogeneous isotropic turbulence. We discuss that the energy evolution can differ more significantly from that of homogeneous isotropic turbulence when there is a time variation of the gravitational potential. Morphologically, ICM turbulence can show a layered vertical structure and large horizontal vortical eddies in the central regions with the greatest density stratification. In addition, we find that the coefficient in the linear density fluctuation - turbulence Mach number relation caused by density stratification is in general a variable with position and time.

On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

The binary emission-line system, particularly the binary broad-line emission system, is considered the most effective indicator of the black hole binary. A plausible example of such a system, SDSS J153636.22+044127.0, was reported as the first known object with two hydrogen Balmer broad-line systems, which are interpreted to be the result of broad-line regions around a pair of black holes (Boroson \& Lauer 2009). Here, we show the follow-up optical and near-infrared spectral observations of SDSS J153636.22+044127.0 and its analogues. In these objects, the broad hydrogen Balmer and Paschen, He I and Mg II lines share the same peculiar emission-line profile (including a blue system, a red system and a double-peaked disk-line component); however, the invariance in the large time interval, the absence of the blue system in He I $\lambda$10830 profile and the abnormally strong emission of the hydrogen Pa$\beta$ blue system oppose the binary proposal. We suggest that these unique broad lines arise from the AGN emission-line region and the shock-heated outflowing gases rather than a binary system of two active black holes.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma production rates and their relation to the nucleus characteristics, and inherent spatial resolution of the data. This study relies on a detailed 3D nucleus shape model, illumination conditions and the pointing information of the viewing geometry. With these parameters, we can evaluate the relative contribution of water density originating from facets directly inside the MIRO beam as well outside of the beam as a function of distance along the MIRO line-of-sight. We also calculate the ratio of in-beam versus out-of-beam number density. We demonstrate that despite the rather small MIRO field-of-view there is only a small fraction of molecules that originate from facets within the MIRO beam. This is true for nadir, but also translated into the limb observing geometry. The MIRO instrument cannot discriminate active from non-active regions directly from observations. This study also suggests that the beam averaged solar incidence angle, local time and mean normal vectors are not necessary related to molecules within the MIRO beam. These results also illustrate why the 1D spherical Haser model can be applied with relative success for analyzing the MIRO data (and generally any Rosetta measurements). The future possibilities of constraining gas activity distribution on the surface should use 3D codes extracting information from the MIRO spectral line shapes which contain additional information. The presented results are in fact applicable to all relevant instruments onboard Rosetta.

In April 2016, the Rosetta spacecraft performed a low-altitude low-phase-angle flyby over the Imhotep-Khepry transition of 67P/Churyumov-Gerasimenko's nucleus. The OSIRIS/Narrow-Angle-Camera (NAC) acquired 112 images with mainly 3 broadband filters in the visible at a resolution of up to 0.53 m/px and for phase angles between 0.095{\deg} and 62{\deg}. Using those images, we have investigated the morphological and spectrophotometrical properties of this area. We assembled the images into coregistered color cubes. Using a 3D shape model, we produced the illumination conditions and georeference for each image. We projected the observations on a map to investigate its geomorphology. Observations were photometrically corrected using the Lommel-Seeliger disk law. Spectrophotometric analyses were performed on the coregistered color cubes. These data were used to estimate the local phase reddening. This region of the nucleus hosts numerous and varied types of terrains and features. We observe an association between a feature's nature, its reflectance, and its spectral slope. Fine material deposits exhibit an average reflectance and spectral slope, while terrains with diamictons, consolidated material, degraded outcrops, or features such as somber boulders, present a lower-than-average reflectance and higher-than-average spectral slope. Bright surfaces present here a spectral behavior consistent with terrains enriched in water-ice. We find a phase-reddening slope of 0.064{\pm}0.001{\%}/100nm/{\deg} at 2.7 au outbound, similarly to the one obtained at 2.3 au inbound during the February 2015 flyby. Identified as the source region of multiple jets and a host of water-ice material, the Imhotep-Khepry transition appeared in April 2016, close to the frost line, to further harbor several potential locations with exposed water-ice material among its numerous different morphological terrain units.

The third-generation gravitational wave detectors are under development by operating the detector in cryogenic temperature to reduce the thermal noise. Silicon and sapphire are promising candidate materials for the test masses and suspension elements due to their remarkable mechanical and thermal properties at cryogenic temperature. Here we present the performances of the cryogenic thermal cycling and strength testing on hydroxide catalysis bonding between sapphire and silicon. Our results suggest that although these two materials have very different coefficients of thermal expansion, but if the flatness and the thermally grown $\mathrm{SiO_2}$ oxidation layer on the silicon surface are controlled well, the bonded samples can still survive thermal cycling from room temperature to 5.5 K. A breaking strength of 3.6$\pm 0.6$ MPa is measured for the bonds between sapphire and silicon with a 190 nm silicon oxidation thickness after cooling cycle. We construct a hybrid sapphire-silicon Fabry-Perot cavity with the developing bonding technique in our lab. The measurement results reveal that the cavity can survive repeated thermal cycling while maintaining a good finesse.

Matching astronomical catalogs in crowded regions of the sky is challenging both statistically and computationally due to the many possible alternative associations. Budav\'ari and Basu (2016) modeled the two-catalog situation as an Assignment Problem and used the famous Hungarian algorithm to solve it. Here we treat cross-identification of multiple catalogs by introducing a different approach based on integer linear programming. We first test this new method on problems with two catalogs and compare with the previous results. We then test the efficacy of the new approach on problems with three catalogs. The performance and scalability of the new approach is discussed in the context of large surveys.

The phase function of the dust coma of comet 67P has been determined from Rosetta/OSIRIS images \citep{Bertini17}. This function show a deep minimum at phase angles near 100$^\circ$, and a strong backscattering enhancement. These two properties cannot be reproduced by regular models of cometary dust, most of them based on wavelength-sized and randomly-oriented aggregate particles. We show, however, that an ensamble of oriented elongated particles of a wide variety of aspect ratios, with radii $r \gtrsim$10 $\mu$m, and whose long axes are perpendicular to the direction of the solar radiation, are capable of reproducing the observed phase function. These particles must be absorbing, with an imaginary part of the refractive index of about 0.1 to match the expected geometric albedo, and with porosity in the 60-70\% range.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

The diffuse hot medium inside clusters of galaxies typically exhibits turbulent motions whose amplitude increases with radius, as revealed by cosmological hydrodynamical simulations. However, its physical origin remains unclear. It could either be due to an excess injection of turbulence at large radii, or faster turbulence dissipation at small radii. We investigate this by studying the time evolution of turbulence in the intracluster medium (ICM) after major mergers, using the Omega500 non-radiative hydrodynamical cosmological simulations. By applying a novel wavelet analysis to study the radial dependence of the ICM turbulence spectrum, we discover that faster turbulence dissipation in the inner high density regions leads to the increasing turbulence amplitude with radius. We also find that the ICM turbulence at all radii decays in two phases after a major merger: an early fast decay phase followed by a slow secular decay phase. The buoyancy effects resulting from the ICM density stratification becomes increasingly important during turbulence decay, as revealed by a decreasing turbulence Froude number $Fr \sim \mathcal{O}(1)$. Our results indicate that the stronger density stratification and smaller eddy turn-over time are the likely causes of the faster turbulence dissipation rate in the inner regions of the cluster.

We identified a peculiar proximate sub-damped Lya absorption system (sub-DLA) at z=3.234 with a neutral hydrogen column density of $N_{H} =10^{19.70\pm0.15}~\rm cm^{-2}$ toward the quasar SDSS J125903.26+621211.5 in two epoch optical spectra of the Sloan Digital Sky Survey. We detected Lya residue in the proximate sub-DLA trough at $> 8\sigma$ level. To our surprise, extreme ultraviolet (EUV) continuum emission is significantly ($> 4\sigma$) detected in the corresponding Lymann limit absorption edge at both of the FUV and NUV bands by the Galaxy Evolution Explorer. The high neutral hydrogen column density should allow a negligible transmission of both the Lya line photons and EUV continuum photons due to the high optical depth of the gas. The possible scenarios of foreground galaxy contamination, partial coverage, emission from the quasar host galaxy, and extended Lya emission are excluded in turn, and we speculate that the residual Lya and EUV emission is due to photons scattering (broad Lya and the continuum emission) of electrons residing at a spatial scale larger than that of the proximate sub-DLA. Electron scattering is wavelength-independent, the scattered light is therefore a copy of the incident spectrum that might originate in the accretion disk. With the assistances from the neutral hydrogen absorbers as the "natural coronagraph" and the scatterers as the "natural mirror", we yielded a very hard EUV spectral index of $\alpha = 0.65\pm 0.25$ ($F_{\nu}\propto \nu^{\alpha}$), consistent with the standard picture of the locally heated accretion disk in the inner EUV-emitting radii, as well as in the outer near-infrared-emitting radii suggested by Kishimoto et al. (2008).

We report the discovery of highly reddening and hydrogen Balmer and metastable helium broad absorption lines in the quasar SDSS J145057.28+530007.6, based on the optical and near-infrared spectra taken from the SDSS-III/BOSS and the TripleSpec observations. The nuclear continuum, Balmer decrement and absorption-line depth analyses suggest that (1) the accretion disk is completely obscured and the covering factor of the broad-line region (BLR) is only $0.39\pm0.03$, (2) the power-law continuum is reddened by the SMC extinction law of $E(B-V )=0.72\pm 0.01$ mag and the dusty materials are mainly associated with \ion{Ca}{2} H and K rather than the Balmer and \ion{He}{1}* absorption-line system, (3) the unsaturated Balmer (H$\beta$, H$\gamma$, and H$\delta$) and \ion{He}{1}* $\lambda3889$ absorption lines have same two-Gaussian profiles with the shifts of $-931\pm33$ and $-499\pm39$ km s$^{-1}$ and the widths of $121\pm28$ and $196\pm37$ km s$^{-1}$, respectively. Constrained mutually by the Balmer, \ion{He}{1}* absorption lines and undetected \ion{Fe}{2}* $\lambda5169$ in the photoionization simulations, the physical properties of the outflow gases are derived as follows: ionization parameter $10^{-1.4} \lesssim U \lesssim 10^{-0.8}$, density $10^{8.2\pm0.4} \lesssim n_{\rm H} \lesssim 10^{9.0\pm0.4}$ cm$^{-3}$, and column density $10^{22.0\pm0.2} \lesssim N_{\rm H} \lesssim 10^{22.2-22.3}$ cm$^{-2}$. We propose that the ultra-dense outflow gases appear in the vicinity of the surface of the BLR or are located at most 3.12 pc away from the engine. That probably implies that the outflow originates from the BLR, and this kind of ultra-dense BLR scale outflow gases can effectively test the physical properties of the outer gases of the BLR.

In this paper, we present an intercomparison study of two quasars, SDSS J145618.32+340037.2 and SDSS J215331.50-025514.1, which have proximate damped Lya systems (PDLAs) with residual flux upon the Lya absorption trough. Though they both have residual flux as luminous as $10^{43}erg/s$, their PDLAs are quite different in, e.g., HI column density, metal line absorption strength, high-ionization absorption lines as well as residual flux strength. For J1456, the HI column density is $20.6\pm0.2$, with za=2.3138, nearly identical to the quasar redshift (z=2.3142) determined from the [OIII] emission line. The metallicity of this system is typical of DLAs and there is high ionization therein, suggesting that the PDLA system is multiphase, putting it in the quasar environment. For J2153, we measure the HI column density to be $21.5\pm0.1$ at za=3.511, slightly redshifted with respect to the quasar (z=3.490) measured from CIII]. The metallicity of this system is quite low and there is a lack of significant high-ionization absorption lines therein, suggesting that the system is beyond the quasar host galaxy. The residual flux is wide (1000 km/s) in J1456, with a significance of $8\sigma$, while also wide (1500 km/s) but with a smaller significance of $3\sigma$ in J2153. Among many explanations, we find that Lya fuzz or resonant scattering can be used to explain the residual flux in the two sources while partial coverage cannot be excluded for J1456. By comparing these two cases, together with similar cases reported previously, we suggest that the strength of the residual flux is related to properties such as metallicity and high-ionization absorption lines of PDLAs. The residual flux recorded upon the PDLA absorption trough opens a window for us to see the physical conditions and processes of the quasar environment, and their profile and strength further remind us of their spatial scales.

About 3196 EA-type binaries (EAs) were observed by LAMOST by June 16, 2017 and their spectral types were derived. Meanwhile stellar atmospheric parameters of 2020 EAs were determined. In the paper, those EAs are catalogued and their physical properties and evolutionary states are investigated. The period distribution of EAs suggests that the period limit of tidal locking for the close binaries is about 6 days. It is found that the metallicity of EAs is higher than that of EWs indicating that EAs are generally younger than EWs and they are the progenitors of EWs. The metallicities of long-period EWs (0.4 < P < 1 days) are the same as those of EAs with the same periods, while their values of Log (g) are usually smaller than those of EAs. These support the evolutionary process that EAs evolve into long-period EWs through the combination of angular momentum loss (AML) via magnetic braking and case A mass transfer. For short-period EWs, their metallicities are lower than those of EAs, while their gravitational accelerations are higher. These reveal that they may be formed from cool short-period EAs through AML via magnetic braking with little mass transfer. For some EWs with high metallicities, they may be contaminated by material from the evolution of unseen neutron stars and black holes or they have third bodies that may help them to form rapidly through a short timescale of pre-contact evolution. The present investigation suggests that the modern EW populations may be formed through the combination of aforementioned mechanisms.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

The processes that led to the formation of the planetary bodies in the Solar System are still not fully understood. Using the results obtained with the comprehensive suite of instruments on-board ESA's Rosetta mission, we present evidence that comet 67P/Churyumov-Gerasimenko likely formed through the gentle gravitational collapse of a bound clump of mm-sized dust aggregates ("pebbles"), intermixed with microscopic ice particles. This formation scenario leads to a cometary make-up that is simultaneously compatible with the global porosity, homogeneity, tensile strength, thermal inertia, vertical temperature profiles, sizes and porosities of emitted dust, and the steep increase in water-vapour production rate with decreasing heliocentric distance, measured by the instruments on-board the Rosetta spacecraft and the Philae lander. Our findings suggest that the pebbles observed to be abundant in protoplanetary discs around young stars provide the building material for comets and other minor bodies.

LS I +61${^\circ}$303 is a $\gamma$-ray emitting X-ray binary with periodic radio outbursts with time scales of one month. Previous observations have revealed microflares superimposed on these large outbursts with periods ranging from a few minutes to hours. This makes LS I +61${^\circ}$303, along with Cyg X-1, the only TeV emitting X-ray binary exhibiting radio microflares. To further investigate these microflaring activity in LS I +61${^\circ}$303 we observed the source with the 100-m Effelsberg radio telescope at 4.85, 8.35, and 10.45 GHz and performed timing analysis on the obtained data. Radio oscillations of 15 hours time scales are detected at all three frequencies. We also compare the spectral index evolution of radio data to that of the photon index of GeV data observed by Fermi-LAT. We conclude that the observed QPO could result from multiple shocks in a jet.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

We report the detection of heavy reddening and the He I* $\lambda10830$ absorption lines at the AGNs' redshift in three Narrow-line Seyfert 1 galaxies: SDSS J091848.61+211717.0, SDSS J111354.66+124439.0, and SDSS J122749.13+321458.9. They exhibit very red optical to near-infrared colors, narrow Balmer/Paschen broad emission lines and He I* $\lambda10830$ absorption lines. The ultraviolet-optical-infrared nucleus continua are reddened by the SMC extinction law of $E(B-V)\sim 0.74$, 1.17, and 1.24 mag for three objects, which are highly consistent with the values obtained from the broad-line Balmer decrements, but larger than those of narrow emission lines. The reddening analysis suggests the extinction dust simultaneously obscure the accretion disk, the broad emission line region, and the hot dust from the inner edge of the torus. It is possible that the dust obscuring the AGN structures are the dusty torus itself. Furthermore, the Cloudy analysis of the He I* $\lambda10830$ absorption lines propose the distance of the absorption materials to be the extend scale of the torus. That greatly increases probabilities of the obscure and absorption materials being the dusty torus.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

Numerous EWs were discovered by several deep photometric survey and there are about 40785 EW-type binary systems listed in the international variable star index (VSX) by March 13, 2017. 7938 of them were observed by LAMOST by November 30, 2016 and their spectral types were given. Stellar atmospheric parameters of 5363 EW-type binary stars were determined based on good spectroscopic observations. In the paper, those EWs were catalogued and their properties are analyzed. The distributions of the orbital period (P), the effect temperature (T), the gravitational acceleration (Log(g)), the metallicity ([Fe/H]) and the radial velocity (RV) are presented for those observed EW-type systems. It is shown that about 80.6\% sample stars have metallicity below zero indicating that EW-type systems are old stellar population. This is in agreement with the the conclusion that the EW binaries are formed from moderately close binaries through angular momentum loss via magnetic braking that takes a few hundred million to a few billion years. The unusual high metallicities of a few percent of EWs may be caused by contaminating of material from the evolution of unseen neutron stars and black holes in the systems. The correlations between the orbital period and the effect temperature, the gravitational acceleration and the metallicity are presented and their scatters are mainly caused by (i) the presence of the third bodies and (ii) the wrong determined periods sometimes. It is shown that some EW contain evolved component stars and the physical properties of EWs are mainly depending on their orbital periods. It is found that the extremely short-period EWs may be older than their long-period cousins because they have lower metallicities. This reveals that they have a longer timescale of pre-contact evolution and their formation and evolution are mainly driven by angular momentum loss via magnetic braking.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

In this paper we present an analysis of absorption line variability in mini-BAL quasar LBQS 1206+1052. The SDSS spectrum demonstrates that the absorption troughs can be divided into two components of blueshift velocities of $\sim$700 km s$^{-1}$ and $\sim$1400 km s$^{-1}$ relative to the quasar rest-frame. The former component shows rare Balmer absorption, which is an indicator of high density absorbing gas, thus the quasar is worth follow-up spectroscopic observations. Our follow-up optical and near-infrared spectra using MMT, YFOSC, TripleSpec and DBSP reveal that the strengths of the absorption lines vary for both of the two components, while the velocities do not change. We reproduce all of the spectral data by assuming that only the ionization state of the absorbing gas is variable and that all other physical properties are invariable. The variation of ionization is consistent with the variation of optical continuum from the V-band light-curve. Additionally, we can not interpret the data by assuming that the variability is due to a movement of the absorbing gas. Therefore, our analysis strongly indicates that the absorption line variability in LBQS 1206+1052 is photoionization-driven. As shown from photo-ionization simulations, the absorbing gas with blueshift velocity of $\sim$700 km s$^{-1}$ has a density in the range of $10^9$ to $10^{10}$ cm$^{-3}$ and a distance of $\sim$1 pc, and the gas with blueshift velocity of $\sim$1400 km s$^{-1}$ has a density of $10^3$ cm$^{-3}$ and a distance of $\sim$1 kpc.

Broad emission-line outflows of active galactic nuclei (AGNs) have been proposed for many years but are very difficult to quantitatively study because of the coexistence of the gravitationally-bound and outflow emission. We present detailed analysis of a heavily reddened quasar, SDSS J000610.67+121501.2, whose normal ultraviolet (UV) broad emission lines (BELs) are heavily suppressed by the Dusty Torus as a natural "Coronagraph", thus the blueshifted BELs (BBELs) can be reliably measured. The physical properties of the emission-line outflows are derived as follows: ionization parameter $U \sim 10^{-0.5}$, column density $N_{\rm H}\sim 10^{22.0}$ cm$^{-2}$, covering fraction of $\sim 0.1$ and upper limit density of $n_{\rm H}\sim 10^{5.8}$ cm$^{-3}$. The outflow gases are located at least 41 pc away from the central engine, which suggests that they have expanded to the scale of the dust torus or beyond. Besides, Lya shows a narrow symmetric component, to our surprise, which is undetected in any other lines. After inspecting the narrow emission-line region and the starforming region as the origin of the Lya narrow line, we propose the end-result of outflows, diffusing gases in the larger region, acts as the screen of Lya photons. Future high spatial resolution spectrometry and/or spectropolarimetric observation are needed to make a final clarification.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

We present OSIRIS/NAC observations of decimetre-sized, likely ice-containing aggregates ejected from a confined region on the surface of comet 67P/Churyumov-Gerasimenko. The images were obtained in January 2016 when the comet was at 2 AU from the Sun out-bound from perihelion. We measure the acceleration of individual aggregates through a two-hour image series. Approximately 50% of the aggregates are accelerated away from the nucleus, and 50% towards it, and likewise towards either horizontal direction. The accelerations are up to one order of magnitude stronger than local gravity, and are most simply explained by the combined effect of gas drag accelerating all aggregates upwards, and the recoil force from asymmetric outgassing, either from rotating aggregates with randomly oriented spin axes and sufficient thermal inertia to shift the temperature maximum away from an aggregate's subsolar region, or from aggregates with variable ice content. At least 10% of the aggregates will escape the gravity field of the nucleus and feed the comet's debris trail, while others may fall back to the surface and contribute to the deposits covering parts of the northern hemisphere. The rocket force plays a crucial role in pushing these aggregates back towards the surface. Our observations show the future back fall material in the process of ejection, and provide the first direct measurement of the acceleration of aggregates in the innermost coma (<2km) of a comet, where gas drag is still significant.

The accretion of interstellar medium onto the central super massive black holes is widely accepted as the source of the gigantic energy released by the active galactic nuclei. But few pieces of observational evidence have been confirmed directly demonstrating the existence of the inflows. The absorption line system in the spectra of quasar SDSS J112526.12+002901.3 presents an interesting example, in which the rarely detected hydrogen Balmer and metastable He I absorption lines are found redshifted to the quasar's rest frame along with the low-ionization metal absorption lines Mg II, Fe II, etc. The repeated SDSS spectroscopic observations suggest a transverse velocity smaller than the radial velocity. The motion of the absorbing medium is thus dominated by infall. The He I* lines present a powerful probe to the strength of ionizing flux, while the Balmer lines imply a dense environment. With the help of photoionization simulations, we find the absorbing medium is exposed to the radiation with ionization parameter $U\approx 10^{-1.8}$, and the density is $n(\mathrm{H})\approx 10^9\ \mathrm{cm}^{-3}$. Thus the absorbing medium is located $\sim 4\ \mathrm{pc}$ away from the central engine. According to the similarity in the distance and physical conditions between the absorbing medium and the torus, we strongly propose the absorption line system as a candidate for the accretion inflow which originates from the inner surface of the torus.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

Accretion shocks around galaxy clusters mark the position where the infalling diffuse gas is significantly slowed down, heated up, and becomes a part of the intracluster medium (ICM). They play an important role in setting the ICM properties. Hydrodynamical simulations have found an intriguing result that the radial position of this accretion shock tracks closely the position of the `splashback radius' of the dark matter, despite the very different physical processes that gas and dark matter experience. Using the self-similar spherical collapse model for dark matter and gas, we find that an alignment between the two radii happens only for a gas with an adiabatic index of $\gamma \approx 5/3$ and for clusters with moderate mass accretion rates. In addition, we find that some observed ICM properties, such as the entropy slope and the effective polytropic index lying around $\sim 1.1-1.2$, are captured by the self-similar spherical collapse model, and are insensitive to the mass accretion history.

A steepening feature in the outer density profiles of dark matter halos indicating the splashback radius has drawn much attention recently. Possible observational detections have even been made for galaxy clusters. Theoretically, Adhikari et al. have estimated the location of the splashback radius by computing the secondary infall trajectory of a dark matter shell through a growing dark matter halo with an NFW profile. However, since they imposed a shape of the halo profile rather than computing it consistently from the trajectories of the dark matter shells, they could not provide the full shape of the dark matter profile around the splashback radius. We improve on this by extending the self-similar spherical collapse model of Fillmore \& Goldreich to a $\Lambda$CDM universe. This allows us to compute the dark matter halo profile and the trajectories simultaneously from the mass accretion history. Our results on the splashback location agree qualitatively with Adhikari et al. but with small quantitative differences at large mass accretion rates. We present new fitting formulae for the splashback radius $R_{\rm sp}$ in various forms, including the ratios of $R_{\rm sp} / R_{\rm 200c}$ and $R_{\rm sp} / R_{\rm 200m}$. Numerical simulations have made the puzzling discovery that the splashback radius scales well with $R_{\rm 200m}$ but not with $R_{\rm 200c}$. We trace the origin of this to be the correlated increase of $\Omega_{\rm m}$ and the average halo mass accretion rate with an increasing redshift.

SDSS J1634+2049 is a local (z = 0.1293) infrared-luminous quasar with LIR= 10^11.91 Lsun. We present a detailed multiwavelength study of both the host galaxy and the nucleus. The host galaxy demonstrates violent, obscured star formation activities with SFR ~ 140 Msun yr^-1, estimated from either the PAH emission or IR luminosity. The optical to NIR spectra exhibit a blueshifted narrow cuspy component in Hb, HeI5876,10830 and other emission lines consistently with an offset velocity of ~900 km/s, as well as additional blueshifting phenomena in high-ionization lines , while there exist blueshifted broad absorption lines (BALs) in NaID and HeI*3889,10830, indicative of the AGN outflows producing BALs and emission lines. Constrained mutually by the several BALs with CLOUDY, the physical properties of the absorption-line outflow are derived as follows: 10^4 < n_H <= 10^5 cm^-3, 10^-1.3 <= U <= 10^-0.7 and 10^22.5<= N_H <= 10^22.9 cm^-2 , similar to those derived for the emission-line outflows. The similarity suggests a common origin. Taking advantages of both the absorption lines and outflowing emission lines, we find that the outflow gas is located at a distance of 48 - 65 pc from the nucleus, and that the kinetic luminosity of the outflow is 10^44-10^46 erg s^-1. J1634+2049 has a off-centered galactic ring on the scale of ~ 30 kpc that is proved to be formed by a recent head-on collision by a nearby galaxy. Thus this quasar is a valuable object in the transitional phase emerging out of dust enshrouding as depicted by the co-evolution scenario.

SDSS J095253.83+011422.0 (SDSS J0952+0114) was reported by Hall et al. (2004) as an exotic quasar at $z_{\rm em}=3.020$. In contrast to prominent broad metal--line emissions with FWHM~9000 km/s, only a narrow Ly \alpha emission line is present with FWHM~1000 km/s. The absence of broad Ly alpha emission line has been a mystery for more than a decade. In this paper, we demonstrate that this is due to dark Proximate Damped Ly alpha Absorption (PDLA) at $z_{\rm abs}=3.010$ by identifying associated Lyman absorption line series from the damped Ly beta up to Ly9, as well as the Lyman limit absorption edge. The PDLA cloud has a column density of $\log N_{\rm H\,I}({\rm cm}^{-2})=21.8\pm0.2$, a metallicity of [Zn/H]$>-1.0$, and a spatial extent exceeding the Narrow Emission Line Region (NELR) of the quasar. With a luminosity of $L_{{\rm Ly}\alpha}\sim10^{45}$ erg s$^{-1}$, the residual Ly alpha emission superposed on the PDLA trough is of two orders of magnitude stronger than previous reports. This is best explained as re-radiated photons arising from the quasar outflowing gas at a scale larger than the NELR. The PDLA here, acting like a natural coronagraph, provides us with a good insight into the illuminated gas in the vicinity of the quasar, which are usually hard to resolve due to their small size and "seeing fuzz" of bright quasars. Notably, SDSS J0952+0114 analogs might be easily omitted in the spectroscopic surveys of DLAs and PDLAs, as their damped Ly alpha troughs can be fully filled by additional strong Ly alpha emissions. Our preliminary survey shows that such systems are not very rare. They are potentially a unique sample for probing strong quasar feedback phenomena in the early universe.

We report on the discovery of broad Balmer absorption lines variability in the QSO SDSS J125942.80+121312.6, based on the optical and near-infrared spectra taken from the SDSS-I, SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and TripleSpec observations over a timescale of 5.8 years in the QSO's rest-frame. The blueshifted absorption profile of H$\beta$ shows a variation of more than 5$\sigma$ at a high velocity portion ($>3000\ \mathrm{km\ s}^{-1}$) of the trough. We perform a detailed analysis for the physical conditions of the absorber using Balmer lines as well as metastable He I and optical Fe II absorptions ($\lambda 4233$ from b$^4$P$_{5/2}$ level and $\lambda 5169$ from a$^6$S$_{5/2}$) at the same velocity. These Fe II lines are identified in the QSO spectra for the first time. According to the photoionization simulations, we estimate a gas density of $n(\mathrm{H})\approx 10^{9.1}\ \mathrm{cm}^{-3}$ and a column density of $N_{\mathrm{col}}(\mathrm{H})\approx 10^{23}\ \mathrm{cm}^{-2}$ for the BOSS data, but the model fails to predict the variations of ionic column densities between the SDSS and BOSS observations if changes in ionizing flux are assumed. We thus propose transverse motion of the absorbing gas being the cause of the observed broad Balmer absorption line variability. In fact, we find that the changes in covering factors of the absorber can well-reproduce all of the observed variations. The absorber is estimated $\sim 0.94$ pc away from the central engine, which is where the outflow likely experiences deceleration due to the collision with the surrounding medium. This scheme is consistent with the argument that LoBAL QSOs may represent the transition from obscured star-forming galaxies to classic QSOs.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

The existence of intermediate-width emission line regions (IELRs) in active galactic nuclei has been discussed for over two decades. A consensus, however, is yet to be arrived at due to the lack of convincing evidence for their detection. We present a detailed analysis of the broadband spectrophotometry of the partially obscured quasar OI 287. The ultraviolet intermediate-width emission lines (IELs) are very prominent, in high contrast to the corresponding broad emission lines (BELs) which are heavily suppressed by dust reddening. Assuming that the IELR is virialized, we estimated its distance to the central black hole of $\sim 2.9$ pc, similar to the dust sublimation radius of $\sim 1.3$ pc. Photo-ionization calculations suggest that the IELR has a hydrogen density of $\sim 10^{8.8}-10^{9.4} ~ \rm cm^{-3}$, within the range of values quoted for the dusty torus near the sublimation radius. Both its inferred location and physical conditions suggest that the IELR originates from the inner surface of the dusty torus. In the spectrum of this quasar, we identified only one narrow absorption-line system associated with the dusty material. With the aid of photo-ionization model calculations, we found that the obscuring material might originate from an outer region of the dusty torus. We speculate that the dusty torus, which is exposed to the central ionizing source, may produce IELs through photo-ionization processes, while also obscure BELs as a natural "coronagraph". Such a "coronagraph" could be found in a large number of partially obscured quasars and be a useful tool to study IELRs.

Non-thermal pressure in galaxy clusters leads to underestimation of the mass of galaxy clusters based on hydrostatic equilibrium with thermal gas pressure. This occurs even for dynamically relaxed clusters that are used for calibrating the mass-observable scaling relations. We show that the analytical model for non-thermal pressure developed in Shi & Komatsu 2014 can correct for this so-called 'hydrostatic mass bias', if most of the non-thermal pressure comes from bulk and turbulent motions of gas in the intracluster medium. Our correction works for the sample average irrespective of the mass estimation method, or the dynamical state of the clusters. This makes it possible to correct for the bias in the hydrostatic mass estimates from X-ray surface brightness and the Sunyaev-Zel'dovich observations that will be available for clusters in a wide range of redshifts and dynamical states.

We present the discovery of Balmer line absorption from H$\alpha$ to H$\gamma$ in an iron low-ionization broad absorption line (FeLoBAL) quasar SDSS J152350.42+391405.2 (hereafter J1523), by the quasi-simultaneous optical and near-infrared spectroscopy. The Balmer line absorption is at $z_{absor}$ = 0.6039 +/-0.0021 and blueshifted by v=10,353 km/s with respect to the Balmer emission lines. All Balmer BALs have uniform absorption profile with the widths of $\Delta$ v ~12,000 km/s. We also found the absorption trough in He 1* $\lambda$10830 with the same velocity and width in the H-band TripleSpec spectrum of J1523. This object is only the tenth active galactic nucleus known to exhibit non-stellar Balmer absorption, and also the case with the highest velocity and broadest Balmer absorption lines which have ever been found. A CLOUDY analysis shows that the absorbers require an gas density of $log_{10} n_ e (cm^{-3})=9$ and an ionization parameter of $log_{10} U=-1.0$. They locate at a distance of ~0.2 pc from the central ionizing source which is slightly farther than that of BELRs. Furthermore, J1523 is one of the brightest Balmer-BAL quasar ever reported, with unique iron absorption variations, making it as the most promising candidate for follow up high-resolution spectroscopy, multi-band observations, and long-term monitoring.

The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

We present the G-virial method (available at http://gxli.github.io/G-virial/) which aims to quantify (1) the importance of gravity in molecular clouds in the position-position-velocity (PPV) space, and (2) properties of the gas condensations in molecular clouds. Different from previous approaches that calculate the virial parameter for different regions, our new method takes gravitational interactions between all the voxels in 3D PPV data cubes into account, and generates maps of the importance of gravity. This map can be combined with the original data cube to derive relations such as the mass-radius relation. Our method is important for several reasons. First, it offers the the ability to quantify the centrally condensed structures in the 3D PPV data cubes, and enables us to compare them in an uniform framework. Second, it allows us to understand the importance of gravity at different locations in the data cube, and provides a global picture of gravity in clouds. Third, it offers a robust approach to decomposing the data into different regions which are gravitationally coherent. To demonstrate the application of our method we identified regions from the Perseus and Ophiuchus molecular clouds, and analyzed their properties. We found an increase in the importance of gravity towards the centers of the individual molecular condensations. We also quantified the properties of the regions in terms of mass-radius and mass-velocity relations. Through evaluating the virial parameters based on the G-virial, we found that all our regions are almost gravitationally bound. Cluster-forming regions appear are more centrally condensed.

We present the variability study of broad absorption lines (BALs) in a uniformly radio-selected sample of 28 BAL quasars using the archival data from the FIRST Bright Quasar Survey (FBQS) and the Sloan Digital Sky Survey (SDSS), as well as those obtained by ourselves, covering time scales $\sim 1-10$ years in the quasar's rest-frame. To our surprise, 5 quasars showing strong variations are all belong to a special subclass of overlapping iron low ionization BAL (OFeLoBAL) quasars, however, other 4 non-overlapping FeLoBALs (non-OFeLoBALs) are invariable except one case with weak optical depth change. Meanwhile, we also find 6 typical variations of high-ionization and low-ionization BALs in this BAL quasar sample. Photoionization models suggest that OFeLoBALs are formed in a relative dense ($n_e>10^6$ cm$^{-3}$) outflows at a distance from the subparsec to the dozens of parsecs from the continuum source. They differ from those of non-OFeLoBALs, which are likely produced by low-density gas, locating at a distance of hundreds to thousands parsecs. Thus, OFeLoBALs and non-OFeLoBALs, i.e., FeLoBALs with/without strong BAL variations, perhaps represent the bimodality of Fe II absorption, the former is located in the active galactic nucleus environment rather than the host galaxy. We suggest that high density and small distance are the necessary conditions that cause OFeLoBALs. As suggested in the literature, strong BAL variability is possibly due to variability of the covering factor of BAL regions caused by clouds transiting across the line of sight rather than ionization variations.

Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.

Turbulent gas motion inside galaxy clusters provides a non-negligible non-thermal pressure support to the intracluster gas. If not corrected, it leads to a systematic bias in the estimation of cluster masses from X-ray and Sunyaev-Zel'dovich (SZ) observations assuming hydrostatic equilibrium, and affects interpretation of measurements of the SZ power spectrum and observations of cluster outskirts from ongoing and upcoming large cluster surveys. Recently, Shi & Komatsu developed an analytical model for predicting the radius, mass, and redshift dependence of the non-thermal pressure contributed by the kinetic random motions of intracluster gas sourced by the cluster mass growth. In this paper, we compare the predictions of this analytical model to a state-of-the-art cosmological hydrodynamics simulation. As different mass growth histories result in different non-thermal pressure, we perform the comparison on 65 simulated galaxy clusters on a cluster-by-cluster basis. We find an excellent agreement between the modelled and simulated non-thermal pressure profiles. Our results open up the possibility of using the analytical model to correct the systematic bias in the mass estimation of galaxy clusters. We also discuss tests of the physical picture underlying the evolution of intracluster non-thermal gas motions, as well as a way to further improve the analytical modeling, which may help achieve a unified understanding of non-thermal phenomena in galaxy clusters.

Non-thermal pressure in the intracluster gas has been found ubiquitously in numerical simulations, and observed indirectly. In this paper we develop an analytical model for intracluster non-thermal pressure in the virial region of relaxed clusters. We write down and solve a first-order differential equation describing the evolution of non-thermal velocity dispersion. This equation is based on insights gained from observations, numerical simulations, and theory of turbulence. The non-thermal energy is sourced, in a self-similar fashion, by the mass growth of clusters via mergers and accretion, and dissipates with a time-scale determined by the turnover time of the largest turbulence eddies. Our model predicts a radial profile of non-thermal pressure for relaxed clusters. The non-thermal fraction increases with radius, redshift, and cluster mass, in agreement with numerical simulations. The radial dependence is due to a rapid increase of the dissipation time-scale with radii, and the mass and redshift dependence comes from the mass growth history. Combing our model for the non-thermal fraction with the Komatsu-Seljak model for the total pressure, we obtain thermal pressure profiles, and compute the hydrostatic mass bias. We find typically 10% bias for the hydrostatic mass enclosed within $r_{500}$.

With 3rd-order statistics of gravitational shear it will be possible to extract valuable cosmological information from ongoing and future weak lensing surveys which is not contained in standard 2nd-order statistics, due to the non-Gaussianity of the shear field. Aperture mass statistics are an appropriate choice for 3rd-order statistics due to their simple form and their ability to separate E- and B-modes of the shear. However, it has been demonstrated that 2nd-order aperture mass statistics suffer from E-/B-mode mixing because it is impossible to reliably estimate the shapes of close pairs of galaxies. This finding has triggered developments of several new 2nd-order statistical measures for cosmic shear. Whether the same developments are needed for 3rd-order shear statistics is largely determined by how severe this E-/B-mixing is for 3rd-order statistics. We test 3rd-order aperture mass statistics against E-/B-mode mixing, and find that the level of contamination is well-described by a function of $\theta/\theta_{\rm min}$, with $\theta_{\rm min}$ being the cut-off scale. At angular scales of $\theta > 10 \;\theta_{\rm min}$, the decrease in the E-mode signal due to E-/B-mode mixing is smaller than 1 percent, and the leakage into B-modes is even less. For typical small-scale cut-offs this E-/B-mixing is negligible on scales larger than a few arcminutes. Therefore, 3rd-order aperture mass statistics can safely be used to separate E- and B-modes and infer cosmological information, for ground-based surveys as well as forthcoming space-based surveys such as Euclid.

Solar-cycle related variation of differential rotation is investigated through analyzing the rotation rates of magnetic fields, distributed along latitudes and varying with time at the time interval of August 1976 to April 2008. More pronounced differentiation of rotation rates is found to appear at the ascending part of a Schwabe cycle than at the descending part on an average. The coefficient $B$ in the standard form of differential rotation, which represents the latitudinal gradient of rotation, may be divided into three parts within a Schwabe cycle. Part one spans from the start to the $4^{th}$ year of a Schwabe cycle, within which the absolute $B$ is approximately a constant or slightly fluctuates. Part two spans from the $4^{th}$ to the $7^{th}$ year, within which the absolute $B$ decreases. Part three spans from the $7^{th}$ year to the end, within which the absolute $B$ increases. Strong magnetic fields repress differentiation of rotation rates, so that rotation rates show less pronounced differentiation, but weak magnetic fields seem to just reflect differentiation of rotation rates. The solar-cycle related variation of solar differential rotation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation of rotation rates.

The latitudinal distributions of the yearly mean rotation rates measured respectively by Suzuki in 1998 and 2012 and Pulkkinen $\&$ Tuominen in 1998 are utilized to investigate internal-cycle variation of solar differential rotation. The rotation rate at the solar Equator seems to decrease since cycle 10 onwards. The coefficient $B$ of solar differential rotation, which represents the latitudinal gradient of rotation, is found smaller in the several years after the minimum of a solar cycle than in the several years after the maximum time of the cycle, and it peaks several years after the maximum time of the solar cycle. The internal-cycle variation of the solar rotation rates looks similar in profile to that of the coefficient $B$. A new explanation is proposed to address such a solar-cycle related variation of the solar rotation rates. Weak magnetic fields may more effectively reflect differentiation at low latitudes with high rotation rates than at high latitudes with low rotation rates, and strong magnetic fields may more effectively repress differentiation at relatively low latitudes than at high latitudes. The internal-cycle variation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation.

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential rotation law obtained by Makarov, Tlatov, and Callebaut (1997), Javaraiah, Bertello, and Ulrich (2005a, b), and Javaraiah et al. (2009). Our results indicate that the solar surface rotation rate at the Equator (indicated by the A parameter of the standard solar rotation law) shows a secular decrease since cycle 12 onwards, given by about $1\,-\,1.5\times10^{-3}$($deg\ day^{-1} year^{-1}$). The B parameter of the standard differential rotation law seems to also show a secular decrease since cycle 12 onwards, but of weak statistical significance. The rotation rate averaged on latitudes ($0^{o}\,--\,40^{o}$) does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation in the long run.

The daily sunspot numbers of the whole disk as well as the northern and southern hemispheres from January 1, 1945 to December 31, 2010 are used to investigate the temporal variation of the rotational cycle length through the continuous wavelet transformation analysis method. The auto-correlation function analysis of daily hemispheric sunspot numbers shows that the southern hemisphere rotates faster than the northern hemisphere. The results obtained from the wavelet transformation analysis are: there exists no direct relationship between the variation trend of the rotational cycle length and the variation trend of solar activity in the two hemispheres; the rotational cycle length of both hemispheres has no significant period appearing at the 11 years, but has significant period of about 7.6 years. Analysis concerning the solar cycle dependence of the rotational cycle length shows that in the whole disk and the northern hemisphere acceleration seems to appear before the minimum time of solar activity. Furthermore, the cross-correlation study indicates that the rotational cycle length of the two hemispheres has different phases, and the rotational cycle length of the whole disk as well as the northern and southern hemispheres also has phase shifts with the corresponding solar activity. What's more, the temporal variation of North-South (N-S) asymmetry of the rotational cycle length is also studied; it displays the same variation trend as the N-S asymmetry of solar activity in a solar cycle as well as in the considered time interval, and it has two significant periods of 7.7 and 17.5 years. Moreover, the N-S asymmetry of the rotational cycle length and the N-S asymmetry of solar activity are highly correlated. It's inferred that the northern hemisphere should rotate faster at the beginning of solar cycle 24.

With the growing interest in and ability of using weak lensing studies to probe the non-Gaussian properties of the matter density field, there is an increasing need for the study of suitable statistical measures, e.g. shear three-point statistics. In this paper we establish the relations between the three-point configuration space shear and convergence statistics, which are an important missing link between different weak lensing three-point statistics and provide an alternative way of relating observation and theory. The method we use also allows us to derive the relations between other two- and three-point correlation functions. We show the consistency of the relations obtained with already established results and demonstrate how they can be evaluated numerically. As a direct application, we use these relations to formulate the condition for E/B-mode decomposition of lensing three-point statistics, which is the basis for constructing new three-point statistics which allow for exact E/B-mode separation. Our work applies also to other two-dimensional polarization fields such as that of the Cosmic Microwave Background.

Three-point weak lensing statistics provide cosmic information complementary to that of two-point statistics. However, both statistics suffer from intrinsic-shear alignment, which is one of their limiting systematics. The nulling technique is a model-independent method developed to eliminate intrinsic-shear alignment at the two-point level. In this paper we demonstrate that the nulling technique can also be naturally generalized to the three-point level, controlling the corresponding GGI systematics. We show that under the assumption of exact redshift information the intrinsic-shear alignment contamination can be completely eliminated. To show how well the nulling technique performs on data with limited redshift information, we apply the nulling technique to three-point weak lensing statistics from a fictitious survey analogous to a typical future deep imaging survey, in which the three-point intrinsic-shear alignment systematics is generated from a power-law toy model. Using 10 redshift bins, the nulling technique leads to a factor of 10 suppression of the GGI/GGG ratio, and reduces the bias on cosmological parameters to less than the original statistical error. More detailed redshift information allowing for finer redshift bins leads to better reduction of bias. The information loss during the nulling procedure doubles the statistical error on cosmological parameters. A comparison of the nulling technique with an unconditioned compression of the data suggests that part of the information loss can be retained by considering higher-order nulling weights during the nulling procedure. A combined analysis of two- and three-point statistics confirms that the information contained in them is of comparable size and is complementary to each other, both before and after nulling.

To probe cosmological fields beyond the Gaussian level, three-point statistics can be used, all of which are related to the bispectrum. Hence, measurements of CMB anisotropies, galaxy clustering, and weak gravitational lensing alike have to rely upon an accurate theoretical background concerning the bispectrum and its noise properties. If only small portions of the sky are considered, it is often desirable to perform the analysis in the flat-sky limit. We aim at a formal, detailed derivation of the bispectrum covariance in the flat-sky approximation, focusing on a pure two-dimensional Fourier-plane approach. We define an unbiased estimator of the bispectrum, which takes the average over the overlap of annuli in Fourier space, and compute its full covariance. The outcome of our formalism is compared to the flat-sky spherical harmonic approximation in terms of the covariance, the behavior under parity transformations, and the information content. We introduce a geometrical interpretation of the averaging process in the estimator, thus providing an intuitive understanding. Contrary to foregoing work, we find a difference by a factor of two between the covariances of the Fourier-plane and the spherical harmonic approach. We argue that this discrepancy can be explained by the differing behavior with respect to parity. However, in an exemplary analysis it is demonstrated that the Fisher information of both formalisms agrees to high accuracy. Via the geometrical interpretation we are able to link the normalization in the bispectrum estimator to the area enclosed by the triangle configuration at consideration as well as to the Wigner symbol, which leads to convenient approximation formulae for the covariances of both approaches.

The time dependence of the changes in the emission spectra of Comet 9P/Tempel 1 after Deep Impact are derived and discussed. This was a unique event because for the first time it gave astronomers the opportunity to follow the time history of the formation and decay of O(1S), OH, CN, C2, C3, NH, and NH2. Least squares fits of a modified Haser model with constraints using known rate constants were fit to the observed data. In the case of OH a simple two-step Haser model provides a reasonable fit to the observations. Fitting the emissions from O(1S), CN, C2, C3, NH, and NH2 requires the addition of a delayed component to a regular two or three step Haser model. From this information a picture of the Deep Impact encounter emerges where there is an initial formation of gas and dust, which is responsible for the prompt emission that occurs right after impact. A secondary source of gas starts later after impact when the initial dust has dissipated enough so that solar radiation can reach the surface of freshly exposed material. The implications of this and other results are discussed in terms of the implications on the structure and composition of the comet's nucleus.

We assess the effects on Big Bang Nucleosynthesis (BBN) of lepton number generation in the early universe resulting from the two-doublet four-neutrino mass/mixing scheme. It has been argued that this neutrino mass/mixing arrangement gives the most viable fit to the existing data. We study full 4 x 4 mixing matrices and show how possible symmetries in these can affect the BBN He-4 abundance yields. Though there is as yet no consensus on the reliability of BBN calculations with neutrino flavor mixing, we show that, in the case where the sign of the lepton number asymmetry is unpredictable, BBN considerations may pick out specific relationships between mixing angles. In particular, reconciling the observed light element abundances with a \bar\nu_\mu <-> \bar\nu_e oscillation interpretation of LSND would allow unique new constraints on the neutrino mixing angles in this model.

We examine in detail the changes in the production of primordial helium resulting from nonthermal neutrino momentum distributions produced by resonant transformation of electron-type neutrinos to steriles. These transformations, anti-nu_e -> anti-nu_s (nu_e->nu_s), amplify a positive (negative) lepton number asymmetry. We find that the resulting supression relative to a thermal distribution of low energy nu_e reduces n -> p conversion to a greater extent than does the enhancement of n -> p from an identical suppression of anti-nu_e. Thus, equal lepton-number asymmetries of opposite sign have unequal effects on the resulting helium yield in primordial nucleosynthesis.

This is a reply to the preprint "On Neutrino-Mixing-Generated Lepton Asymmetry and the Primordial Helium-4 Abundance" by M. V. Chizhov and D. P. Kirilova (hep-ph/9908525), which criticised our recent publication (X. Shi, G. M. Fuller and K. Abazajian Phys. Rev. D 60, 063002 (1999)). Here we point out factual errors in their description of what our paper says. We also show that their criticisms of our work have no merit.

It has been proposed that an asymmetry in the electron neutrino sector may be generated by resonant active-sterile neutrino transformations during Big Bang Nucleosynthesis (BBN). We calculate the change in the primordial $^4$He yield $Y$ resulting from this asymmetry, taking into account both the time evolution of the $\nu_e$ and $\bar\nu_e$ distribution function and the spectral distortions in these. We calculate this change in two schemes: (1) a lepton asymmetry directly generated by $\nu_e$ mixing with a lighter right-handed sterile neutrino $\nu_s$; and (2) a lepton asymmetry generated by a $\nu_\tau\leftrightarrow\nu_s$ or $\nu_\mu\leftrightarrow\nu_s$ transformation which is subsequently partially converted to an asymmetry in the $\nu_e\bar\nu_e$ sector by a matter-enhanced active-active neutrino transformation. In the first scheme, we find that the percentage change in $Y$ is between -1% and 9% (with the sign depending on the sign of the asymmetry), bounded by the Majorana mass limit $m_{\nu_e}\la 1$ eV. In the second scheme, the maximal percentage reduction in $Y$ is 2%, if the lepton number asymmetry in neutrinos is positive; Otherwise, the percentage increase in $Y$ is $\la 5%$ for $m^2_{\nu_\mu,\nu_\tau}-m^2_{\nu_s}\la 10^4$ eV. We conclude that the change in the primordial $^4$He yield induced by a neutrino-mixing-generated lepton number asymmetry can be substantial in the upward direction, but limited in the downward direction.

We investigate the effects of matter-enhanced (MSW) transformation of neutrinos in the Early Universe on the primordial helium yield (Y_p). We find that Y_p is affected much more by the MSW-induced alterations in the neutrino energy spectra than by the associated change in expansion rate. Specifically, the absence due to transformation of low energy electron neutrinos can significantly affect neutron-proton weak interconversion rates through the lifting of Fermi-blocking of neutron decay at low energies and through halting low-energy neutrino capture on neutrons. We find that the change of Y_p within a causal horizon is -0.005\le \delta Y_p \leq 0.013 for m_{\nu_{\mu,\tau}}^2 - m_{\nu_s}^2 \le 10^4 eV^2 in the case of \nu_{\mu,\tau}-\nu_s-\nu_e mixing, with the lower limit at m_{\nu_{\mu,\tau}}^2 - m_{\nu_s}^2 \approx 100 eV^2, and -0.002\le \delta Y_p \leq 0.020 for m_{\nu_e}^2 - m_{\nu_s}^2 \le 1 eV^2 in the case of \nu_e-\nu_s mixing.

We extend a treatment of the causal structure of space-time to active-sterile neutrino transformation-based schemes for lepton number generation in the early universe. We find that these causality considerations necessarily lead to the creation of spatial domains of lepton number with opposite signs. Lepton number gradients at the domain boundaries can open a new channel for MSW resonant production of sterile neutrinos. The enhanced sterile neutrino production via this new channel allows considerable tightening of Big Bang Nucleosynthesis constraints on active-sterile neutrino mixing, including the proposed $\nu_\mu \to \nu_s$ solution for the Super Kamiokande atmospheric $\nu_\mu$ deficit, and the four-neutrino schemes proposed to simultaneously fit current neutrino experimental results.

We propose a new and unique dark matter candidate: $\sim 100$ eV to $\sim 10$ keV sterile neutrinos produced via lepton number-driven resonant MSW (Mikheyev-Smirnov-Wolfenstein) conversion of active neutrinos. The requisite lepton number asymmetries in any of the active neutrino flavors range from 10$^{-3}$ to 10$^{-1}$ of the photon number - well within primordial nucleosynthesis bounds. The unique feature here is that the adiabaticity condition of the resonance strongly favors the production of lower energy sterile neutrinos. The resulting non-thermal (cold) energy spectrum can cause these sterile neutrinos to revert to non-relativistic kinematics at an early epoch, so that free-streaming lengths at or below the dwarf galaxy scale are possible. Therefore, the main problem associated with light neutrino dark matter candidates can be circumvented in our model.

We suggest that the deceleration of the relativistic shock by a denser part of the interstellar medium off line-of-sight produced the observed radio ``flare'' in the early afterglow of GRB990123. We find that this scenario is consistent with observations if the particle number density of this denser part of the medium is between $\sim 200$ and $\sim 2\times 10^4$ cm$^{-3}$. Because of the premature deceleration of part of the shock, the later stage of the afterglow should decay modestly faster than the powerlaw expected from an isotropic shock propagation.

It has been suggested that the supermassive black holes, at the centers of galaxies and quasars, may initially form in single collapses of relativistic star clusters or supermassive stars built-up during the evolution of dense star clusters. We show that it may be possible for ICECUBE (a planned 1 km^3 neutrino detector in Antarctica) to detect the neutrino bursts associated with those collapses at redshift $z\la 0.2$ with a rate of $\sim$ 0.1 to 1 burst per year. Such detections could give new insights into the formation of structure in the universe, especially when correlated with gravitational wave signatures or even gamma-ray bursts.

We discuss Big Bang Nucleosynthesis constraints on maximal $\nu_\mu\leftrightarrow\nu_s$ mixing. Vacuum $\nu_\mu\leftrightarrow\nu_s$ oscillation has been proposed as one possible explanation of the Super Kamiokande atmospheric neutrino data. Based on the most recent primordial abundance measurements, we find that the effective number of neutrino species for Big Bang Nucleosynthesis (BBN) is $N_\nu\la 3.3$. Assuming that all three active neutrinos are light (with masses $\ll 1$ MeV), we examine BBN constraints on $\nu_\mu\leftrightarrow\nu_s$ mixing in two scenarios: (1) a negligible lepton asymmetry (the standard picture); (2) the presence of a large lepton asymmetry which has resulted from an amplification by $\nu_\tau\leftrightarrow\nu_{s'}$ mixing ($\nu_{s'}$ being $\nu_s$ or another sterile neutrino species). The latter scenario has been proposed recently to reconcile the BBN constraints and large-angle $\nu_\mu\leftrightarrow\nu_s$ mixing. We find that the large-angle $\nu_\mu\leftrightarrow\nu_s$ mixing in the first scenario, which would yield $N_\nu\approx 4$, is ruled out as an explanation of the Super Kamiokande data. It is conceivably possible for the $\nu_\mu\leftrightarrow\nu_s$ solution to evade BBN bounds in the second scenario, but only if 200 eV$^2\la m^2_{\nu_\tau}-m^2_{\nu_{s'}}\la 10^4$ eV$^2$ is satisfied, and if $\nu_\tau$ decays non-radiatively with a lifetime $\la 10^3$ years. This mass-squared difference implies 15 eV$\la m_{\nu_\tau}\la 100$ eV if $\nu_{s'}$ is much lighter than $\nu_\tau$. We conclude that maximal (or near maximal) $\nu_\mu\leftrightarrow\nu_\tau$ mixing is a more likely explanation of the Super Kamiokande data.

In the paper "Big Bang Nucleosynthesis and Active-Sterile Neutrino Mixing: Evidence for Maximal Muon-Neutrino/Sterile-Neutrino Mixing in Super Kamiokande" (astro-ph/9810075), we suggested that to evade the Big Bang Nucleosynthesis exclusion of the muon neutrino to sterile neutrino oscillation explanation of the Super Kamiokande data, the tau neutrino must have a mass over about 15 eV and it must mix with a lighter sterile neutrino. A stable tau neutrino with this mass is inconsistent with cosmological structure formation. In a comment on our paper (astro-ph/9811067), Foot and Volkas argued that our result is incorrect and that the required tau neutrino mass should be much lower. Here we back up our original result with a more detailed calculation. We show that the argument of Foot and Volkas is invalid, most likely due to an insufficient energy resolution in the low energy part of the neutrino spectrum.

The formation of supermassive black holes through the gravitational collapse of supermassive objects ($M \ga 10^4 M_\odot$) has been proposed as a source of cosmological gamma-ray bursts. The major advantage of this model is that such collapses are far more energetic than stellar-remnant mergers. The major drawback of this idea is the severe baryon loading problem in one-dimensional models. We can show that the observed log N - log P (number vs. peak flux) distribution for gamma-ray bursts in the BATSE database is not inconsistent with an identification of supermassive object collapse as the origin of the gamma-ray bursts. This conclusion is valid for a range of plausible cosmological and gamma-ray burst spectral parameters.

We propose that the gravitational collapse of supermassive objects ($ M\ga 10^4 M_\odot$), either as relativistic star clusters or as single supermassive stars (which may result from stellar mergers in dense star clusters), could be a cosmological source of $\gamma$-ray bursts. These events could provide the seeds of the supermassive black holes observed at the center of many galaxies. Collapsing supermassive objects will release a fraction of their huge gravitational binding energy as thermal neutrino pairs. We show that the accompanying neutrino/antineutrino annihilation-induced heating could drive electron/positron ``fireball'' formation, relativistic expansion, and associated $\gamma$-ray emission. The major advantage of this model is its energetics: supermassive object collapses are far more energetic than solar mass-scale compact object mergers; therefore, the conversion of gravitational energy to fireball kinetic energy in the supermassive object scenario need not be highly efficient, nor is it necessary to invoke directional beaming. The major weakness of this model is difficulty in avoiding a baryon loading problem for one dimensional collapse scenarios.

COBE-normalized flat (matter plus cosmological constant) and open Cold Dark Matter (CDM) models are tested by comparing their expected Hubble flow variations and the observed variations in a Type Ia supernova sample and a Tully Fisher cluster sample. The test provides a probe of the CDM power spectrum on scales of $0.02h$ Mpc$^{-1}\la k\la 0.2h$ Mpc$^{-1}$, free of the bias factor $b$. The results favor a low matter content universe, or a flat matter-dominated universe with a very low Hubble constant and/or a very small spectral index $n_{ps}$, with the best fits having $\Omega_0\sim 0.3$ to 0.4. The test is found to be more discriminative to the open CDM models than to the flat CDM models. For example, the test results are found to be compatible with those from the X-ray cluster abundance measurements at smaller length scales, and consistent with the galaxy and cluster correlation analysis of Peacock and Dodds (1994) at similar length scales, if our universe is flat; but the results are marginally incompatible with the X-ray cluster abundance measurements if our universe is open. The open CDM results are consistent with that of Peacock and Dodds only if the matter density of the universe is less than about 60% of the critical density. The shortcoming of the test is discussed, so are ways to minimize it.

We calculate the luminosity and energy spectrum of the neutrino emission from electron-positron pair annihilation during the collapse of a supermassive star (${M} \ga 5\times10^4 {M_\odot}$). We then estimate the cumulative flux and energy spectrum of the resulting neutrino background as a function of the abundance and redshift of supermassive stars and the efficiency of these objects in converting gravitational energy into neutrino energy. We estimate the expected signal in some of the new generation of astrophysical neutrino detectors from both a cumulative background of supermassive stars and single collapse events associated with these objects.

When Cold Dark Matter (CDM) power spectrum normalized by COBE results, its amplitude at smaller scales can be parametrized by $\Gamma\sim\Omega_0h$. The expected variations of Hubble flows in two samples, the sample of 36 clusters in the Mark III catalogue, and the sample of 20 Type Ia supernovae (SNe), are calculated for the power spectrum in critical-density CDM models (including tilted CDM models and vacuum-dominated CDM models). The comparison between the expectations and the real variations in the data offers a bias-free way to constrain the power spectrum. The cluster sample yields $\Gamma\le 0.30-0.88(n_{ps}-1)+1.9 (n_{ps}-1)^2$ at 95% C.L., with best fits $\Gamma = 0.15-0.39(n_{ps}-1)+0.37(n_{ps}-1)^2$, where $n_{ps}$ is the spectral index of the power spectrum. The Type Ia SN sample yields $\Gamma\le 0.25-0.80(n_{ps}-1)+1.6 (n_{ps}-1)^2$ at 95% C.L., strongly favoring lower $\Gamma$'s. The results are inconsistent with a critical-density matter-dominated universe with $n_{ps}\ga 0.8$ and $H_0\ga 50$ km/sec/Mpc.

There are irreducible differences between the Hubble constant measured locally and the global value. They are due to density perturbations and finite sample volume (cosmic variance) and finite number of objects in the sample (sampling variance). We quantify these differences for a suite of COBE-normalized CDM models that are consistent with the observed large-scale structure. For small samples of objects that only extend out to 10,000 km/sec, the variance can approach 4%. For the largest samples of Type Ia supernovae (SNeIa), which include about 40 objects and extend out to almost 40,000 km/sec, the variance is 1-2% and is dominated by sampling variance. Sampling and cosmic variance may be an important consideration in comparing local determinations of the Hubble constant with precision determinations of the global value that will be made from high-resolution maps of CBR anisotropy.

We describe how thermally produced anti-electron neutrinos ($\bar\nu_e$) from the homologously collapsing core of a supermassive star ($M \ga 5\times10^4 M_\odot$) can lead to significant deuterium enrichment in the ejected envelope of such a star. Deuterium-enriched material at high redshift might then serve as a clue to the existence of pregalactic supermassive stars. Conceivably, the ejected deuterium-enriched material could intercept the line of sight to a distant QSO and mimic a Lyman limit absorber. In such a case, the deuterium abundance inferred from absorption lines might not reflect the true primordial abundance of deuterium. We discuss relevant theoretical uncertainties in supermassive star physics as well as potential observational signatures in Lyman $\alpha$ absorber clouds for processing by stars of this kind.

A formalism that simultaneously searches for the monopolar and dipolar peculiar velocities is presented. The formalism is applied to (1) the Mark III catalogue, (2) Lauer and Postman's Abell cluster catalogue, and (3) Riess et al.'s Type Ia supernova sample. The emphasis is drawn to the monopolar peculiar velocities, i.e., peculiar Hubble flows, within these samples. The samples show inconsistent peculiar Hubble flows within a depth of $\sim 60h^{-1}$ Mpc. Beyond a depth of $\sim 80h^{-1}$ Mpc, the Hubble flows of all samples converge to the global Hubble flow to better than 10% at the $2\sigma$ level. The results are compared with theoretical predictions. They at face value disfavor models predicting smaller peculiar velocities such as the tilted Cold Dark Matter model. Limitations of the catalogues are discussed, so are ways to improve the catalogues so that an accurate map of Hubble flows in our local universe can be drawn and be compared with theoretical predictions.

We investigate in detail the parameter space of active-sterile neutrino oscillations that amplifies neutrino chemical potentials at the epoch of Big Bang Nucleosynthesis. We calculate the magnitude of the amplification and show evidences of chaos in the amplification process. We also discuss the implications of the neutrino chemical potential amplification in the Big Bang Nucleosynthesis. It is shown that with a $\sim 1$ eV $\nue$, the amplification of its chemical potential by active-sterile neutrino oscillations can lower the effective number of neutrino species at Big Bang Nucleosynthesis to significantly below 3.

The proximity effect of metal line absorption systems in spectra of high redshift quasars and quasar pairs is investigated through the distribution of \Lyalpha forests. No strong proximity effect is found down to 500 km/sec away from a metal line absorption system for lyman $\alpha$ clouds with $13.3\le\log\NHI <13.8$ at a redshift of $\sim 3$. Its implication on properties of high redshift galaxies and \Lyalpha clouds is discussed. More high resolution spectra of high redshift quasars and pairs of quasars (whose lines of sight are separated by $\sim 0.3h^{-1}$ to 1$h^{-1}$ Mpc) are needed so that the proximity effect can be reliably determined for \Lyalpha clouds with various $\NHI$ and down to a closer distance.

Recent observations of Cepheids in the Virgo cluster have bolstered the evidence that supports a Hubble constant in 70-90 km/s/Mpc range. This evidence, by and large, probes the expansion of the Universe within 100 Mpc. We investigate the possibility that the expansion rate within this region is systematically higher than the true expansion rate due to the presence of a local, large underdense region or void. We begin by calculating the expected deviations between the locally measured Hubble constant and the true Hubble constant for a variety of models. We also discuss the expected correlations between these deviations and mass fluctuation for the sample volume. We find that the fluctuations are small for the standard cold dark matter as well as mixed dark matter models but can be substantial in a number of interesting and viable nonstandard scenarios. However, deviations in the Hubble flow for a region of radius 200 Mpc are small for virtually all reasonable models. Therefore, methods based on supernovae or the Sunyaev-Zel'dovich effect, which can probe 200 Mpc scales, will be essential in determining the true Hubble constant. We discuss, in detail, the fluctuations induced in the cosmic background radiation by voids at the last scattering surface. In addition, we discuss the dipole and quadrupole fluctuations one would expect if the void enclosing us is aspherical or if we lie off-center.

Observations show that there is significant clustering of QSO metal absorption lines within the range of velocity dispersion between 200km/sec and 600 km/sec. With a reasonable supernova rate, it is shown that high velocity gases driven by SNe and/or strong stellar winds could explain the clustered absorptions, provided that QSO metal-line absorbers are galactic halos or dwarf galaxies. Rich clusters of galaxies, on the other hand, cannot yield the observed clustering of QSO metal absorption lines.

The age of globular clusters inferred from observations depends sensitively on assumptions such as the initial helium abundance and the mass loss rate. A high helium abundance (e.g., $Y\approx$0.28), as well as an inclusion of helium diffusion and oxygen-enhancement in stellar models, can lower the current age estimate for metal-poor globular clusters from 14$\pm 1.5$ Gyr to about 11$\pm 1$ Gyr, significantly relaxing the constraints on the Hubble constant, allowing values as high as 60km/sec/Mpc for a universe with the critical density and 90km/sec/Mpc for a baryon-only universe. The uncertainties of a high helium abundance and an instability strip induced mass loss near the turn-off in globular clusters are discussed. Ages lower than 10 Gyr are not possible even with the operation of both of these mechanisms unless the initial helium abundance in globular clusters is $>0.28$, which can hardly be accomodated by indirect inferences of helium abundances in globular clusters.

Without assuming any solar models and neutrino flavor conversions, $^8$B neutrinos seen by the Kamiokande experiment should contribute 2.6$\pm 0.45$ SNU to the chlorine experiment. When this rate is compared with the total event rate of 2.3$\pm 0.2$ SNU observed by the Homestake experiment which should inlude a 0.2 SNU contribution from uncertainty-free pep neutrinos, there may still be a possible evidence that $^7$Be neutrinos are more severely suppressed than the $^8$B neutrinos with respect to the predictions of standard solar models, which cannot be explained by any known astrophysics solution. Given a Kamiokande event rate of larger than 36$\%$ (2$\sigma$) of the prediciton of Bahcall and Pinsonneault's standard solar model, variations of standard solar models yield minimal rates of 3.6 SNU for the Homestake experiment and 114 SNU for GALLEX and SAGE. Therefore, variations of standard solar models as solutions to the solar neutrino problem are so far inconsistent with the Homestake experiment and only marginally allowed by the gallium experiments. If the gallium experiments confirm a flux significantly below 114 SNU, it would imply unconventional physics.

The ages of globular clusters inferred from observations depends sensitively on assumptions like the initial helium abundance and the mass loss rate. A high helium abundance (e.g., $Y\approx$0.28) or a mass loss rate of $\sim$10$^{-11}M_\odot$ yr$^{-1}$ near the main sequence turn-off region lowers the current age estimate from 14 Gyr to about 10--12 Gyr, significantly relaxing the constraints on the Hubble constant, allowing values as high as 60km/sec/Mpc for a universe with the critical density and 90km/sec/Mpc for a baryon-only universe. Possible mechanisms for the helium enhancement in globular clusters are discussed, as are arguments for an instability strip induced mass loss near the turn-off. Ages lower than 10 Gyr are not possible even with the operation of both of these mechanisms unless the initial helium abundance in globular clusters is $>0.30$, which would conflict with indirect measurements of helium abundances in globular clusters.

The role of a resonant $\nu_e$-$\nu_s$ oscillation is discussed in the event of a supernova explosion. It is concluded that a significant $\nu_e$-$\nu_s$ mixing may hinder the ability of the supernova to explode. It may also cool the proto-neutron star too quickly with respect to the observed cooling time of several seconds. The constraints on the $\nu_e$-$\nu_s$ mixing parameters based on the above arguments are calculated.

We discuss in detail the effect of neutrino oscillations in Big Bang nucleosynthesis, between active and sterile neutrinos, as well as between active and active neutrinos. We calculate the constraints on mixings between active and sterile neutrinos from the present observation of the primordial helium abundance and discuss the potential implications on various astrophysical and cosmological problems of such oscillations. In particular, we show that large angle sterile neutrino mixing seems to be excluded as a MSW solution to the solar neutrino situation or a solution to the atmospheric neutrino mixing hinted at in some underground experiments. We show how with this constraint, the next generation of solar neutrino experiments should be able to determine the resolution of the solar neutrino problem. It is also shown how sterile neutrinos remain a viable dark matter candidate.