We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

Quenched massive spiral galaxies have attracted great attention recently, as more data is available to constrain their environment and cold gas content. However, the quenching mechanism is still uncertain, as it depends on the mass range and baryon budget of the galaxy. In this letter, we report the identification of a rare population of very massive, quenched spiral galaxies with stellar mass $\gtrsim10^{11}{\rm~M_\odot}$ and halo mass $\gtrsim10^{13}{\rm~M_\odot}$ from the Sloan Digital Sky Survey at redshift $z\sim0.1$. Our CO observations using the IRAM-30m telescope show that these galaxies contain only a small amount of molecular gas. Similar galaxies are also seen in the state-of-the-art semi-analytical models and hydro-dynamical simulations. It is found from these theoretical models that these quenched spiral galaxies harbor massive black holes, suggesting that feedback from the central black holes has quenched these spiral galaxies. This quenching mechanism seems to challenge the popular scenario of the co-evolution between massive black holes and massive bulges.

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 Search for Extraterrestrial Intelligence (SETI) attempts to address the possibility of the presence of technological civilizations beyond the Earth. Benefiting from high sensitivity, large sky coverage, an innovative feed cabin for China's Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed the SETI first observations with FAST's newly commisioned 19-beam receiver; we report preliminary results in this paper. Using the data stream produced by the SERENDIP VI realtime multibeam SETI spectrometer installed at FAST, as well as its off-line data processing pipelines, we identify and remove four kinds of radio frequency interference(RFI): zone, broadband, multi-beam, and drifting, utilizing the Nebula SETI software pipeline combined with machine learning algorithms. After RFI mitigation, the Nebula pipeline identifies and ranks interesting narrow band candidate ET signals, scoring candidates by the number of times candidate signals have been seen at roughly the same sky position and same frequency, signal strength, proximity to a nearby star or object of interest, along with several other scoring criteria. We show four example candidates groups that demonstrate these RFI mitigation and candidate selection. This preliminary testing on FAST data helps to validate our SETI instrumentation techniques as well as our data processing pipeline.

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored.

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Filament eruptions occurring at different places within a relatively short time internal, but with a certain physical causal connection are usually known as sympathetic eruption. Studies on sympathetic eruptions are not uncommon. However, in the existed reports, the causal links between sympathetic eruptions remain rather speculative. In this work, we present detailed observations of a sympathetic filament eruption event, where an identifiable causal link between two eruptive filaments is observed. On 2015 November 15, two filaments (F1 in the north and F2 in the south) were located at the southwestern quadrant of solar disk. The main axes of them were almost parallel to each other. Around 22:20 UT, F1 began to erupt, forming two flare ribbons. The southwestern ribbon apparently moved to southwest and intruded southeast part of F2. This continuous intrusion caused F2's eventual eruption. Accompanying the eruption of F2, flare ribbons and post-flare loops appeared in northwest region of F2. Meanwhile, neither flare ribbons nor post-flare loops could be observed in southeastern area of F2. In addition, the nonlinear force-free field (NLFFF) extrapolations show that the magnetic fields above F2 in the southeast region are much weaker than that in the northwest region. These results imply that the overlying magnetic fields of F2 were not uniform. So we propose that the southwest ribbon formed by eruptive F1 invaded F2 from its southeast region with relatively weaker overlying magnetic fields in comparison with its northwest region, disturbing F2 and leading F2 to erupt eventually.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

Discovering pulsars is a significant and meaningful research topic in the field of radio astronomy. With the advent of astronomical instruments such as he Five-hundred-meter Aperture Spherical Telescope (FAST) in China, data volumes and data rates are exponentially growing. This fact necessitates a focus on artificial intelligence (AI) technologies that can perform the automatic pulsar candidate identification to mine large astronomical data sets. Automatic pulsar candidate identification can be considered as a task of determining potential candidates for further investigation and eliminating noises of radio frequency interferences or other non-pulsar signals. It is very hard to raise the performance of DCNN-based pulsar identification because the limited training samples restrict network structure to be designed deep enough for learning good features as well as the crucial class imbalance problem due to very limited number of real pulsar samples. To address these problems, we proposed a framework which combines deep convolution generative adversarial network (DCGAN) with support vector machine (SVM) to deal with imbalance class problem and to improve pulsar identification accuracy. DCGAN is used as sample generation and feature learning model, and SVM is adopted as the classifier for predicting candidate's labels in the inference stage. The proposed framework is a novel technique which not only can solve imbalance class problem but also can learn discriminative feature representations of pulsar candidates instead of computing hand-crafted features in preprocessing steps too, which makes it more accurate for automatic pulsar candidate selection. Experiments on two pulsar datasets verify the effectiveness and efficiency of our proposed method.

The rate of fast radio bursts (FRBs) in the direction of nearby galaxy clusters is expected to be higher than the mean cosmological rate if intrinsically faint FRBs are numerous. In this paper, we describe a targeted search for faint FRBs near the core of the Virgo cluster using the Australian Square Kilometer Array Pathfinder telescope. During 300 hr of observations, we discovered one burst, FRB 180417, with dispersion measure DM $=474.8$ cm$^{-3}$pc. The FRB was promptly followed up by several radio telescopes for 27 h, but no repeat bursts were detected. An optical follow-up of FRB 180417 using the PROMPT5 telescope revealed no new sources down to an $R$-band magnitude of 20.1. We argue that FRB 180417 is likely behind the Virgo cluster as the Galactic and intracluster DM contribution are small compared to the DM of the FRB, and there are no galaxies in the line of sight. The non-detection of FRBs from Virgo constrains the faint-end slope, $\alpha<1.52$ (at 68\% confidence limit), and the minimum luminosity, $L_{\rm min}\gtrsim 2\times 10^{40}$ erg s$^{-1}$ (at 68\% confidence limit), of the FRB luminosity function assuming cosmic FRB rate of $10^4$ FRBs sky$^{-1}$ day$^{-1}$ with flux above 1 Jy located out to redshift of 1. Further FRB surveys of galaxy clusters with high-sensitivity instruments will tighten the constraints on the faint end of the luminosity function and, thus, are strongly encouraged.

The Large Synoptic Survey Telescope will complete its survey in 2022 and produce terabytes of imaging data each night. To work with this massive onset of data, automated algorithms to classify astronomical light curves are crucial. Here, we present a method for automated classification of photometric light curves for a range of astronomical objects. Our approach is based on the gradient boosting of decision trees, feature extraction and selection, and augmentation. The solution was developed in the context of The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) and achieved one of the top results in the challenge.

With the accumulation of big data of CME observations by coronagraphs, automatic detection and tracking of CMEs has proven to be crucial. The excellent performance of convolutional neural network in image classification, object detection and other computer vision tasks motivates us to apply it to CME detection and tracking as well. We have developed a new tool for CME Automatic detection and tracking with MachinE Learning (CAMEL) techniques. The system is a three-module pipeline. It is first a supervised image classification problem. We solve it by training a neural network LeNet with training labels obtained from an existing CME catalog. Those images containing CME structures are flagged as CME images. Next, to identify the CME region in each CME-flagged image, we use deep descriptor transforming to localize the common object in an image set. A following step is to apply the graph cut technique to finely tune the detected CME region. To track the CME in an image sequence, the binary images with detected CME pixels are converted from cartesian to polar coordinate. A CME event is labeled if it can move in at least two frames and reach the edge of coronagraph field of view. For each event, a few fundamental parameters are derived. The results of four representative CMEs with various characteristics are presented and compared with those from four existing automatic and manual catalogs. We find that CAMEL can detect more complete and weaker structures, and has better performance to catch a CME as early as possible.

If dark matter was produced in the early Universe by the decoupling of its annihilations into known particles, there is a sharp experimental target for the size of its coupling. We show that if dark matter was produced by inelastic scattering against a lighter particle from the thermal bath, then its coupling can be exponentially smaller than the coupling required for its production from annihilations. As an application, we demonstrate that dark matter produced by inelastic scattering against electrons provides new thermal relic targets for direct detection and fixed target experiments.

It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the $\Lambda$CDM model. One type of anisotropy recently detected in the SDSS (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called "lopsided". Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set up when comparing pairs of galaxies in $\Lambda$CDM simulations indicating that such a set up is not inconsistent with $\Lambda$CDM, the origin has yet to be explained. Here we examine the origin of such lopsided setups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at $z=0$ display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards $z=0$. Second, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the $z=0$ signal is driven primarily by satellites that are on first approach, who have yet to experience a "flyby". This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.

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.

It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called `shape alignment' has been observed in isolated Milky Way-like galaxies but not in `Local Group'-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the {\it facing} region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the {\it away} region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at $\sim5\sigma$ levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

The Commensal Radio Astronomy Five-hundred-meter Aperture Spherical radio Telescope (FAST) Survey (CRAFTS) utilizes the novel drift-scan commensal survey mode of FAST and can generate billions of pulsar candidate signals. The human experts are not likely to thoroughly examine these signals, and various machine sorting methods are used to aid the classification of the FAST candidates. In this study, we propose a new ensemble classification system for pulsar candidates. This system denotes the further development of the pulsar image-based classification system (PICS), which was used in the Arecibo Telescope pulsar survey, and has been retrained and customized for the FAST drift-scan survey. In this study, we designed a residual network model comprising 15 layers to replace the convolutional neural networks (CNNs) in PICS. The results of this study demonstrate that the new model can sort >96% of real pulsars to belong the top 1% of all candidates and classify >1.6 million candidates per day using a dual--GPU and 24--core computer. This increased speed and efficiency can help to facilitate real-time or quasi-real-time processing of the pulsar-search data stream obtained from CRAFTS. In addition, we have published the labeled FAST data used in this study online, which can aid in the development of new deep learning techniques for performing pulsar searches.

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$.

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

The correlation between the spins of dark matter halos and the large-scale structure (LSS) has been studied in great detail over a large redshift range, while investigations of galaxies are still incomplete. Motivated by this point, we use the state-of-the-art hydrodynamic simulation, Illustris-1, to investigate mainly the spin--LSS correlation of galaxies at redshift of $z=0$. We mainly find that the spins of low-mass, blue, oblate galaxies are preferentially aligned with the slowest collapsing direction ($e_3$) of the large-scale tidal field, while massive, red, prolate galaxy spins tend to be perpendicular to $e_3$. The transition from a parallel to a perpendicular trend occurs at $\sim10^{9.4} M_{\odot}/h$ in the stellar mass, $\sim0.62$ in the g-r color, and $\sim0.4$ in triaxiality. The transition stellar mass decreases with increasing redshifts. The alignment was found to be primarily correlated with the galaxy stellar mass. Our results are consistent with previous studies both in N-body simulations and observations. Our study also fills the vacancy in the study of the galaxy spin--LSS correlation at $z=0$ using hydrodynamical simulations and also provides important insight to understand the formation and evolution of galaxy angular momentum.

Context. Pulsars scintillate. Dynamic spectra show brightness variation of pulsars in the time and frequency domain. Secondary spectra demonstrate the distribution of fluctuation power in the dynamic spectra. Aims. Dynamic spectra strongly depend on observational frequencies, but were often observed at frequencies lower than 1.5 GHz. Scintillation observations at higher frequencies help to constrain the turbulence feature of the interstellar medium over a wide frequency range and can detect the scintillations of more distant pulsars. Methods. Ten pulsars were observed at 2250 MHz (S-band) with the Jiamusi 66 m telescope to study their scintillations. Their dynamic spectra were first obtained, from which the decorrelation bandwidths and time scales of diffractive scintillation were then derived by autocorrelation. Secondary spectra were calculated by forming the Fourier power spectra of the dynamic spectra. Results. Most of the newly obtained dynamic spectra are at the highest frequency or have the longest time span of any published data for these pulsars. For PSRs B0540+23, B2324+60 and B2351+61, these were the first dynamic spectra ever reported. The frequencydependence of scintillation parameters indicates that the intervening medium can rarely be ideally turbulent with a Kolmogorov spectrum. The thin screen model worked well at S-band for the scintillation of PSR B1933+16. Parabolic arcs were detected in the secondary spectra of three pulsars, PSRs B0355+54, B0540+23 and B2154+40, all of which were asymmetrically distributed. The inverted arclets of PSR B0355+54 were seen to evolve along the main parabola within a continuous observing session of 12 hours, from which the angular velocity of the pulsar was estimated that was consistent with the measurement by very long baseline interferometry (VLBI).

Pulsar dispersion measure (DM) accounts for the total electron content between a pulsar and us. High-precision observations for projects of pulsar timing arrays show temporal DM variations of millisecond pulsars. The aim of this paper is to decompose the DM variations of 30 millisecond pulsars by using Hilbert-Huang Transform (HHT) method, so that we can determine the general DM trends from interstellar clouds and the annual DM variation curves from solar wind, interplanetary medium and/or ionosphere. We find that the decomposed annual variation curves of 22 pulsars exhibit quasi-sinusoidal, one component and double components features of different origins. The amplitudes and phases of the curve peaks are related to ecliptic latitude and longitude of pulsars, respectively.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system.

Having achieved 'first-light' right before the opening ceremony on September 25, 2016, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) is being busily commissioned. Its innovative design requires ~1000 points to be measured and driven instead of just the two axes of motion, e.g. Azimuth and Elevation for most of the conventional antennae, to realize pointing and tracking. We have devised a survey plan to utilized the full sensitivity of FAST, while minimizing the complexities in operation the system. The 19-beam L band focal plan array will be rotated to specific angles and taking continuous data streams while the surface shape and the focal cabin stay fixed. Such a survey will cover the northern sky in about 220 full days. Our aim is to obtain data for pulsar search, HI (neutral hydrogen) galaxies, HI imaging, and radio transients, simultaneously, through multiple backends. These data sets could be a significant contribution to all related fields in radio astronomy and remain relevant for decades.

Both simulations and observations have confirmed that the spin of haloes/galaxies is correlated with the large scale structure (LSS) with a mass dependence such that the spin of low-mass haloes/galaxies tend to be parallel with the LSS, while that of massive haloes/galaxies tend to be perpendicular with the LSS. It is still unclear how this mass dependence is built up over time. We use N-body simulations to trace the evolution of the halo spin-LSS correlation and find that at early times the spin of all halo progenitors is parallel with the LSS. As time goes on, mass collapsing around massive halo is more isotropic, especially the recent mass accretion along the slowest collapsing direction is significant and it brings the halo spin to be perpendicular with the LSS. Adopting the $fractional$ $anisotropy$ (FA) parameter to describe the degree of anisotropy of the large-scale environment, we find that the spin-LSS correlation is a strong function of the environment such that a higher FA (more anisotropic environment) leads to an aligned signal, and a lower anisotropy leads to a misaligned signal. In general, our results show that the spin-LSS correlation is a combined consequence of mass flow and halo growth within the cosmic web. Our predicted environmental dependence between spin and large-scale structure can be further tested using galaxy surveys.

In this paper we use high-resolution cosmological simulations to study halo intrinsic alignment and its dependence on mass, formation time and large-scale environment. In agreement with previous studies using N-body simulations, it is found that massive halos have stronger alignment. For given mass, older halos have stronger alignment than younger ones. By identifying the cosmic environment of halo using Hessian matrix, we find that for given mass, halos in cluster regions also have stronger alignment than those in filament. The existing theory has not addressed these dependencies explicitly. In this work we extend the linear alignment model with inclusion of halo bias and find that the halo alignment with its mass and formation time dependence can be explained by halo bias. However, the model can not account for the environment dependence, as it is found that halo bias is lower in cluster and higher in filament. Our results suggest that halo bias and environment are independent factors in determining halo alignment. We also study the halo alignment correlation function and find that halos are strongly clustered along their major axes and less clustered along the minor axes. The correlated halo alignment can extend to scale as large as $100h^{-1}$Mpc where its feature is mainly driven by the baryon acoustic oscillation effect.

The statistical properties of the soft gamma repeater SGR J1550--5418 are investigated carefully. We find that the cumulative distributions of fluence, peak flux and duration can be well fitted by a bent power law, while the cumulative distribution of waiting time follows a simple power law. In particular, the probability density functions of fluctuations of fluence, peak flux, and duration have a sharp peak and fat tails, which can be well fitted by a $q$-Gaussian function. The $q$ values keep approximately steady for different scale intervals, indicating a scale-invariant structure of soft gamma repeaters. Those results support that the origin of soft gamma repeaters is crustquakes of neutron stars with extremely strong magnetic fields.

Both simulations and observations have found that the spin of halo/galaxy is correlated with the large scale environment, and particularly the spin of halo flips in filament. A consistent picture of halo spin evolution in different environments is still lacked. Using N-body simulation we find that halo spin with its environment evolves continuously from sheet to cluster, and the flip of halo spin happens both in filament and nodes. For the flip in filament can be explained by halo formation time and the migrating time when its environment changes from sheet to filament. For low-mass haloes, they form first in sheets and migrate into filaments later, so their mass and spin growth inside filament are lower, and the original spin is still parallel to filament. For massive haloes, they migrate into filaments first, and most of their mass and spin growth are obtained in filaments, so the resulted spin is perpendicular to filament. Our results well explain the overall evolution of cosmic web in the cold dark matter model and can be tested using high-redshift data. The scenario can also be tested against alternative models of dark matter, such as warm/hot dark matter, where the structure formation will proceed in a different way.

Highly polarized components of pulse profiles are investigated by analyzing observational data and simulating the emission processes. The highly polarized components appear at the leading or trailing part of a pulse profile, which preferably have a flat spectrum and a flat polarization angle curve compared with the low polarized components. By considering the emission processes and propagation effects, we simulate the distributions of wave modes and fractional linear polarization within the entire pulsar emission beam. We show that the highly polarized components can appear at the leading, central, and/or trailing parts of pulse profiles in the models, depending on pulsar geometry. The depolarization is caused by orthogonal modes or scattering. When a sight line cuts across pulsar emission beam with a small impact angle, the detected highly polarized component will be of the O mode, and have a flat polarization angle curve and/or a flat spectrum as observed. Otherwise, the highly polarized component will be of the X mode and have a steep polarization angle curve.

In this paper we investigate the impact of different star formation histories (SFHs) on the relation between stellar mass M$_{*}$ and star formation rate (SFR) using a sample of galaxies with reliable spectroscopic redshift zspec>2 drawn from the VIMOS Ultra-Deep Survey (VUDS). We produce an extensive database of dusty model galaxies, calculated starting from the new library of single stellar population (SSPs) models presented in Cassara' et al. 2013 and weighted by a set of 28 different SFHs based on the Schmidt function, and characterized by different ratios of the gas infall time scale $\tau_{infall}$ to the star formation efficiency $\nu$. The treatment of dust extinction and re-emission has been carried out by means of the radiative transfer calculation. The spectral energy distribution (SED) fitting technique is performed by using GOSSIP+, a tool able to combine both photometric and spectroscopic information to extract the best value of the physical quantities of interest, and to consider the Intergalactic Medium (IGM) attenuation as a free parameter. We find that the main contribution to the scatter observed in the $SFR-M_{*}$ plane is the possibility of choosing between different families of SFHs in the SED fitting procedure, while the redshift range plays a minor role. The majority of the galaxies, at all cosmic times, are best-fit by models with SFHs characterized by a high $\tau_{\rm infall}/\nu$ ratio. We discuss the reliability of the presence of a small percentage of dusty and highly star forming galaxies, in the light of their detection in the FIR.

We report the discovery of two new millisecond pulsars (PSRs J0024$-$7204aa and J0024$-$7204ab) in the globular cluster 47\,Tucanae (NGC 104). Our results bring the total number of pulsars in 47\,Tucanae to 25. These pulsars were discovered by reprocessing archival observations from the Parkes radio telescope. We reprocessed the data using a standard search procedure based on the PRESTO software package as well as using a new method in which we incoherently added the power spectra corresponding to $\sim$1100\,hr of observations. The newly discovered PSR~J0024$-$7204aa, has a pulse frequency of $\rm \sim$541\,Hz (corresponding to a $\rm \sim$1.84 ms period), which is higher than any other pulsars currently known in the cluster and ranks 12$^{\rm{th}}$ amongst all the currently known pulsars. The dispersion measure of this pulsar, 24.941(7)\,cm$^{-3}$ pc, is the highest in the cluster. The second discovered pulsar, PSR~J0024$-$7204ab, is an isolated pulsar with a pulse frequency of $\rm \sim$270\,Hz (corresponding to a period of $\rm \sim$3.70 ms).

Based on the basic definition of Fermi energy of degenerate and relativistic electrons, we obtain a special solution to electron Fermi energy, $E_{\rm F}(e)$, and express $E_{\rm F}(e)$ as a function of electron fraction, $Y_{e}$, and matter density, $\rho$. Several useful analytical formulae for $Y_{e}$ and $\rho$ within classical models and the work of Dutra et al. 2014 (Type-2) in relativistic mean field theory are obtained using numerically fitting. When describing the mean-field Lagrangian, density, we adopt the TMA parameter set, which is remarkably consistent with with the updated astrophysical observations of neutron stars. Due to the importance of the density dependence of the symmetry energy, $S$, in nuclear astrophysics, a brief discussion on the symmetry parameters $S_v$ and $L$ (the slope of $S$) is presented. Combining these fit formulae with boundary conditions for different density regions, we can evaluate the value of $E_{\rm F}(e)$ in any given matter density, and obtain a schematic diagram of $E_{\rm F}(e)$ as a continuous function of $\rho$. Compared with previous study on the electron Fermi energy in other models, our methods of calculating $E_{\rm F}(e)$ are more simple and convenient, and can be universally suitable for the relativistic electron regions in the circumstances of common neutron stars. We have deduced a general expression of $E_{\rm F}(e)$ and $n_{e}$, which could be used to indirectly test whether one EoS of a NS is correct in our future studies on neutron star matter properties. Since URCA reactions are expected in the center of a massive star due to high-value electron Fermi energy and electron fraction, this study could be useful in the future studies on the NS thermal evolution.

This paper describes the first data release (DR1) of the VIMOS Ultra Deep Survey (VUDS). The DR1 includes all low-resolution spectroscopic data obtained in 276.9 arcmin2 of the CANDELS-COSMOS and CANDELS-ECFDS survey areas, including accurate spectroscopic redshifts z_spec and individual spectra obtained with VIMOS on the ESO-VLT. A total of 698 objects have a measured redshift, with 677 galaxies, two type-I AGN and a small number of 19 contaminating stars. The targets of the spectroscopic survey are selected primarily on the basis of their photometric redshifts to ensure a broad population coverage. About 500 galaxies have z_spec>2, 48 with z_spec>4, and the highest reliable redshifts reach beyond z_spec=6. This dataset approximately doubles the number of galaxies with spectroscopic redshifts at z>3 in these fields. We discuss the general properties of the sample in terms of the spectroscopic redshift distribution, the distribution of Lyman-alpha equivalent widths, and physical properties including stellar masses M_star and star formation rates (SFR) derived from spectral energy distribution fitting with the knowledge of z_spec. We highlight the properties of the most massive star-forming galaxies, noting the large range in spectral properties, with Lyman-alpha in emission or in absorption, and in imaging properties with compact, multi-component or pair morphologies. We present the catalogue database and data products. All data are publicly available and can be retrieved from a dedicated query-based database available at http://cesam.lam.fr/vuds.

The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($\beta$) as well as Ly$\alpha$ equivalent widths (EW). On average, the spectroscopically measured $\beta$ (-1.36$\pm$0.02), is comparable to the photometrically measured $\beta$ (-1.32$\pm$0.02), and has smaller measurement uncertainties. The positive correlation of $\beta$ with the Spectral Energy Distribution (SED)-based measurement of dust extinction, E$_{\rm s}$(B-V), emphasizes the importance of $\beta$ as an alternative dust indicator at high redshifts. To make a proper comparison, we divide these SFGs into three subgroups based on their rest-frame Ly$\alpha$ EW: SFGs with no Ly$\alpha$ emission (SFG$_{\rm N}$; EW$\le$0\AA), SFGs with Ly$\alpha$ emission (SFG$_{\rm L}$; EW$>$0\AA), and Ly$\alpha$ emitters (LAEs; EW$\ge$20\AA). The fraction of LAEs at these redshifts is $\sim$10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SFG$_{\rm N}$, SFG$_{\rm L}$ and LAE samples. For the luminosities probed here ($\sim$L$^*$), we find that galaxies with and without Ly$\alpha$ in emission have small but significant differences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant compared to the SFR and E$_{\rm s}$(B-V). [abridged]

In quantum gravity, a foamy structure of space-time leads to Lorentz invariance violation (LIV). As the most energetic astrophysical processes in the Universe, gamma-ray bursts (GRBs) provide an effective way to probe quantum gravity effects. We use continuous spectra of 20 short GRBs detected by the Swift satellite to give a conservative lower limit of quantum gravity energy scale $M_\textrm{QG} $. Due to the LIV effect, photons with different energy have different velocities. This will lead to the delayed arrival of high energy photons relative to the low energy ones. Based on the fact that the LIV-induced time delay can't be longer than the duration of a GRB, we present the most conservative estimation of the quantum gravity energy scales from 20 short GRBs. The most strict constraint, $M_\textrm{QG}>5.05\times10^{14}$ GeV, is from GRB 140622A.

The distribution of galaxies displays anisotropy on different scales and it is often referred as galaxy alignment. To understand the origin of galaxy alignments on small scales, one must investigate how galaxies were accreted in the early universe and quantify their primordial anisotropic at the time of accretion. In this paper we use N-body simulations to investigate the accretion of dark matter subhaloes, focusing on their alignment with the host halo shape and the orientation of mass distribution on large scale, defined using the hessian matrix of the density field. The large/small (e1/e3) eigenvalues of the hessian matrix define the fast/slow collapse direction of dark matter on large scale. We find that: 1) the halo major axis is well aligned with the e3 (slow collapse) direction, and it is stronger for massive haloes; 2) subhaloes are predominately accreted along the major axis of the host halo, and the alignment increases with the host halo mass. Most importantly, this alignment is universal; 3) accretion of subhaloes with respect to the e3 direction is not universal. In massive haloes, subhaloes are accreted along the e3 (even stronger than the alignment with the halo major axis), but in low-mass haloes subhaloes are accreted perpendicular to the e3. The transit mass is lower at high redshift. The last result well explains the puzzled correlation (both in recent observations and simulations) that massive galaxies/haloes have their spin perpendicular to the filament, and the spin of low-mass galaxies/haloes is slightly aligned with the filament, under the assumption that the orbital angular momentum of subhaloes is converted to halo spin.

The Lyman continuum (LyC) flux escaping from high-z galaxies into the IGM is a fundamental quantity to understand the physical processes involved in the reionization epoch. We have investigated a sample of star-forming galaxies at z~3.3 in order to search for possible detections of LyC photons escaping from galaxy halos. UV deep imaging in the COSMOS field obtained with the prime focus camera LBC at the LBT telescope was used together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 45 galaxies at z~3.3 with L>0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 3.27<z<3.40 both in the R and deep U bands. A sub-sample of 10 galaxies apparently shows escape fractions>28% but a detailed analysis of their properties reveals that, with the exception of two marginal detections (S/N~2) in the U band, all the other 8 galaxies are most likely contaminated by the UV flux of low-z interlopers located close to the high-z targets. The average escape fraction derived from the stacking of the cleaned sample was constrained to fesc_rel<2%. The implied HI photo-ionization rate is a factor two lower than that needed to keep the IGM ionized at z~3, as observed in the Lyman forest of high-z QSO spectra or by the proximity effect. These results support a scenario where high redshift, relatively bright (L>0.5L*) star-forming galaxies alone are unable to sustain the level of ionization observed in the cosmic IGM at z~3. Star-forming galaxies at higher redshift and at fainter luminosities (L<<L*) can be the major contributors to the reionization of the Universe only if their physical properties are subject to rapid changes from z~3 to z~6-10. Alternatively, ionizing sources could be discovered looking for fainter sources among the AGN population at high-z.

The cosmic web contains a large fraction of the total gas mass in the universe but is difficult to detect at most wavelengths. Synchrotron emission from shock-accelerated electrons may offer the chance of imaging the cosmic web at radio wavelengths. In this work we use 3D cosmological ENZO-MHD simulations (combined with a post-processing renormalisation of the magnetic field to bracket for missing physical ingredients and resolution effects) to produce models of the radio emission from the cosmic web. In post-processing we study the capabilities of 13 large radio surveys to detect this emission. We find that surveys by LOFAR, SKA1-LOW and MWA have a chance of detecting the cosmic web, provided that the magnetisation level of the tenuous medium in filaments is of the order of 1% of the thermal gas energy.

We report the discovery of 31 low-luminosity (-14.5 > M_{AB}(B) > -18.8), extreme emission line galaxies (EELGs) at 0.2 < z < 0.9 identified by their unusually high rest-frame equivalent widths (100 < EW[OIII] < 1700 A) as part of the VIMOS Ultra Deep Survey (VUDS). VIMOS optical spectra of unprecedented sensitivity ($I_{AB}$ ~ 25 mag) along with multiwavelength photometry and HST imaging are used to investigate spectrophotometric properties of this unique sample and explore, for the first time, the very low stellar mass end (M* < 10^8 M$_{\odot}$) of the luminosity-metallicity (LZR) and mass-metallicity (MZR) relations at z < 1. Characterized by their extreme compactness (R50 < 1 kpc), low stellar mass and enhanced specific star formation rates (SFR/M* ~ 10^{-9} - 10^{-7} yr^{-1}), the VUDS EELGs are blue dwarf galaxies likely experiencing the first stages of a vigorous galaxy-wide starburst. Using T_e-sensitive direct and strong-line methods, we find that VUDS EELGs are low-metallicity (7.5 < 12+log(O/H) < 8.3) galaxies with high ionization conditions, including at least three EELGs showing HeII 4686A emission and four EELGs of extremely metal-poor (<10% solar) galaxies. The LZR and MZR followed by EELGs show relatively large scatter, being broadly consistent with the extrapolation toward low luminosity and mass from previous studies at similar redshift. However, we find evidences that galaxies with younger and more vigorous star formation -- as characterized by their larger EWs, ionization and sSFR -- tend to be more metal-poor at a given stellar mass.

Frequency dependence of pulsar linear polarization is investigated by simulations of emission and propagation processes. Linearly polarized waves are generated through curvature radiation by relativistic particles streaming along curved magnetic field lines, which have ordinary mode (O-mode) and extra-ordinary mode (X-mode) components. As emitted waves propagate outwards, two mode components are separated due to re- fraction of the O mode, and their polarization states are also modified. According to the radius to frequency mapping, low frequency emission is generated from higher magnetosphere, where significant rotation effect leads the X and O modes to be sepa- rated. Hence, the low frequency radiation has a large fraction of linear polarization. As the frequency increases, emission is generated from lower heights, where the rotation effect becomes weaker and the distribution regions of two modes are more overlapped. Hence, more significant depolarization appears for emission at higher frequencies. In addition, refraction effect of the O mode becomes serious in very deep magnetosphere, which bends the O mode emission towards outer parts of a pulsar beam and also causes the separation of mode distribution regions and hence the fractional linear polariza- tion increasing with frequency. If emission of different frequencies is generated from a region of the same height, serious O mode refraction can result in the decrease of both profile width and fractional linear polarization. The observed frequency dependence of linear polarization for some pulsars can be naturally explained within the scope of our scenario.

The relation between the galaxy stellar mass M_star and the dark matter halo mass M_h gives important information on the efficiency in forming stars and assembling stellar mass in galaxies. We present the stellar mass to halo mass ratio (SMHR) measurements at redshifts 2<z<5, obtained from the VIMOS Ultra Deep Survey. We use halo occupation distribution (HOD) modelling of clustering measurements on ~3000 galaxies with spectroscopic redshifts to derive the dark matter halo mass M_h, and SED fitting over a large set of multi-wavelength data to derive the stellar mass M_star and compute the SMHR=M_star/M_h. We find that the SMHR ranges from 1% to 2.5% for galaxies with M_star=1.3x10^9 M_sun to M_star=7.4x10^9 M_sun in DM halos with M_h=1.3x10^{11} M_sun} to M_h=3x10^{11} M_sun. We derive the integrated star formation efficiency (ISFE) of these galaxies and find that the star formation efficiency is a moderate 6-9% for lower mass galaxies while it is relatively high at 16% for galaxies with the median stellar mass of the sample ~7x10^9 M_sun. The lower ISFE at lower masses may indicate that some efficient means of suppressing star formation is at work (like SNe feedback), while the high ISFE for the average galaxy at z~3 is indicating that these galaxies are efficiently building-up their stellar mass at a key epoch in the mass assembly process. We further infer that the average mass galaxy at z~3 will start experiencing star formation quenching within a few hundred millions years.

The influence of cosmological constant type dark energy in the early universe is investigated. This is accommodated by a new dispersion relation in de Sitter spacetime. We perform a global fitting to explore the cosmological parameters space by using the CosmoMC package with the recently released Planck TT and WMAP Polarization datasets. Using the results from global fitting, we compute a new CMB temperature-temperature spectrum. The obtained TT spectrum has lower power compared with the one based on $\Lambda$CDM model at large scales.

A number of scientific numerical codes can currently exploit GPUs with remarkable performance. In astrophysics, Enzo and Ramses are prime examples of such applications. The two codes have been ported to GPUs adopting different strategies and programming models, Enzo adopting CUDA and Ramses using OpenACC. We describe here the different solutions used for the GPU implementation of both cases. Performance benchmarks will be presented for Ramses. The results of the usage of the more mature GPU version of Enzo, adopted for a scientific project within the CHRONOS programme, will be summarised.

Using new spectroscopic observations obtained as part of the VIMOS Ultra-Deep Survey (VUDS), we perform a systematic search for overdense environments in the early universe ($z>2$) and report here on the discovery of Cl J0227-0421, a massive protocluster at $z=3.29$. This protocluster is characterized by both the large overdensity of spectroscopically confirmed members, $\delta_{gal}=10.5\pm2.8$, and a significant overdensity in photometric redshift members. The halo mass of this protocluster is estimated, by a variety of methods, to be roughly $3\times10^{14}$ $\mathcal{M}_{\odot}$ at $z\sim3.3$, which, evolved to $z=0$ results in a halo mass rivaling or exceeding that of the Coma cluster. The properties of 19 spectroscopically confirmed member galaxies are compared with a large sample of VUDS/VVDS galaxies in lower density field environments at similar redshifts. We find tentative evidence for an excess of redder, brighter, and more massive galaxies within the confines of the protocluster relative to the field population, which suggests that we may be observing the beginning of environmentally-induced quenching. The properties of these galaxies are investigated, including a discussion of the brightest protocluster galaxy which appears to be undergoing vigorous coeval nuclear and starburst activity. The remaining member galaxies appear to have characteristics which are largely similar to the field population. Though we find weaker evidence of the suppression of the median star formation rates amongst and differences in stacked spectra of member galaxies with respect to the field, we defer any conclusions of these trends to future work with the ensemble of protostructures that are found in the full VUDS sample.

(arXiv abridged abstract) The observed UV rest-frame spectra of distant galaxies are the result of their intrinsic emission combined with absorption along the line of sight produced by the inter-galactic medium (IGM). Here we analyse the evolution of the mean IGM transmission Tr(Ly_alpha) and its dispersion along the line of sight for 2127 galaxies with 2.5<z<5.5 in the VIMOS Ultra Deep Survey (VUDS). We fit model spectra combined with a range of IGM transmission to the galaxy spectra using the spectral fitting algorithm GOSSIP+. We use these fits to derive the mean IGM transmission towards each galaxy for several redshift slices from z=2.5 to z=5.5. We find that the mean IGM transmission defined as Tr(Ly_alpha)=e^{-tau} (with tau the HI optical depth) is 79%, 69%, 59%, 55% and 46% at redshifts 2.75, 3,22, 3.70, 4.23, 4.77, respectively. We compare these results to measurements obtained from quasars lines of sight and find that the IGM transmission towards galaxies is in excellent agreement with quasar values up to redshift z~4. We find tentative evidence for a higher IGM transmission at z>= 4 compared to results from QSOs, but a degeneracy between dust extinction and IGM prevents to draw firm conclusions if the internal dust extinction for star-forming galaxies at z>4 takes a mean value significantly in excess of E(B-V)>0.15. Most importantly, we find a large dispersion of IGM transmission along the lines of sight towards distant galaxies with 68% of the distribution within 10 to 17% of the median value in delta z=0.5 bins, similar to what is found on the LOS towards QSOs. We demonstrate the importance of taking into account this large range of IGM transmission when selecting high redshift galaxies based on their colour properties (e.g. LBG or photometric redshift selection) or otherwise face a significant incompleteness in selecting high redshift galaxy populations.

We investigate the evolution of galaxy clustering for galaxies in the redshift range 2.0<$z$<5.0 using the VIMOS Ultra Deep Survey (VUDS). We present the projected (real-space) two-point correlation function $w_p(r_p)$ measured by using 3022 galaxies with robust spectroscopic redshifts in two independent fields (COSMOS and VVDS-02h) covering in total 0.8 deg$^2$. We quantify how the scale dependent clustering amplitude $r_0$ changes with redshift making use of mock samples to evaluate and correct the survey selection function. Using a power-law model $\xi(r) = (r/r_0)^{-\gamma}$ we find that the correlation function for the general population is best fit by a model with a clustering length $r_0$=3.95$^{+0.48}_{-0.54}$ h$^{-1}$Mpc and slope $\gamma$=1.8$^{+0.02}_{-0.06}$ at $z$~2.5, $r_0$=4.35$\pm$0.60 h$^{-1}$Mpc and $\gamma$=1.6$^{+0.12}_{-0.13}$ at $z$~3.5. We use these clustering parameters to derive the large-scale linear galaxy bias $b_L^{PL}$, between galaxies and dark matter. We find $b_L^{PL}$ = 2.68$\pm$0.22 at redshift $z$~3 (assuming $\sigma_8$ = 0.8), significantly higher than found at intermediate and low redshifts. We fit an HOD model to the data and we obtain that the average halo mass at redshift $z$~3 is $M_h$=10$^{11.75\pm0.23}$ h$^{-1}$M$_{\odot}$. From this fit we confirm that the large-scale linear galaxy bias is relatively high at $b_L^{HOD}$ = 2.82$\pm$0.27. Comparing these measurements with similar measurements at lower redshifts we infer that the star-forming population of galaxies at $z$~3 should evolve into the massive and bright ($M_r$<-21.5) galaxy population which typically occupy haloes of mass $\langle M_h\rangle$ = 10$^{13.9}$ h$^{-1}$ $M_{\odot}$ at redshift $z$=0.

We study the evolution of the star formation rate (SFR) - stellar mass (M_star) relation and specific star formation rate (sSFR) of star forming galaxies (SFGs) since a redshift z~5.5 using 2435 (4531) galaxies with highly reliable (reliable) spectroscopic redshifts in the VIMOS Ultra-Deep Survey (VUDS). It is the first time that these relations can be followed over such a large redshift range from a single homogeneously selected sample of galaxies with spectroscopic redshifts. The log(SFR) - log(M_star) relation for SFGs remains roughly linear all the way up to z=5 but the SFR steadily increases at fixed mass with increasing redshift. We find that for stellar masses M_star>3.2 x 10^9 M_sun the SFR increases by a factor ~13 between z=0.4 and z=2.3. We extend this relation up to z=5, finding an additional increase in SFR by a factor 1.7 from z=2.3 to z=4.8 for masses M_star > 10^10 M_sun. We observe a turn-off in the SFR-M_star relation at the highest mass end up to a redshift z~3.5. We interpret this turn-off as the signature of a strong on-going quenching mechanism and rapid mass growth. The sSFR increases strongly up to z~2 but it grows much less rapidly in 2<z<5. We find that the shape of the sSFR evolution is not well reproduced by cold gas accretion-driven models or the latest hydrodynamical models. Below z~2 these models have a flatter evolution (1+z)^{Phi} with Phi=2-2.25 compared to the data which evolves more rapidly with Phi=2.8+-0.2. Above z~2, the reverse is happening with the data evolving more slowly with Phi=1.2+-0.1. The observed sSFR evolution over a large redshift range 0<z<5 and our finding of a non linear main sequence at high mass both indicate that the evolution of SFR and M_star is not solely driven by gas accretion. The results presented in this paper emphasize the need to invoke a more complex mix of physical processes {abridge}

The amplification of primordial magnetic fields via a small-scale turbulent dynamo during structure formation might be able to explain the observed magnetic fields in galaxy clusters. The magnetisation of more tenuous large-scale structures such as cosmic filaments is more uncertain, as it is challenging for numerical simulations to achieve the required dynamical range. In this work, we present magneto-hydrodynamical cosmological simulations on large uniform grids to study the amplification of primordial seed fields in the intracluster medium (ICM) and in the warm-hot-intergalactic medium (WHIM). In the ICM, we confirm that turbulence caused by structure formation can produce a significant dynamo amplification, even if the amplification is smaller than what is reported in other papers. In the WHIM inside filaments, we do not observe significant dynamo amplification, even though we achieve Reynolds numbers of $R_{\rm e} \sim 200-300$. The maximal amplification for large filaments is of the order of $\sim 100$ for the magnetic energy, corresponding to a typical field of a few $\sim \rm nG$ starting from a primordial weak field of $10^{-10}$ G (comoving). In order to start a small-scale dynamo, we found that a minimum of $\sim 10^2$ resolution elements across the virial radius of galaxy clusters was necessary. In filaments we could not find a minimum resolution to set off a dynamo. This stems from the inefficiency of supersonic motions in the WHIM in triggering solenoidal modes and small-scale twisting of magnetic field structures. Magnetic fields this small will make it hard to detect filaments in radio observations.

Aims. The aim of this work is to constrain the evolution of the fraction of Lya emitters among UV selected star forming galaxies at 2<z<6, and to measure the stellar escape fraction of Lya photons over the same redshift range. Methods. We exploit the ultradeep spectroscopic observations collected by the VIMOS Ultra Deep Survey (VUDS) to build an unique, complete and unbiased sample of 4000 spectroscopically confirmed star forming galaxies at 2<z<6. Our galaxy sample UV luminosities brighter than M* at 2<z<6, and luminosities down to one magnitude fainter than M* at 2<z<3.5. Results. We find that 80% of the star forming galaxies in our sample have EW0(Lya)<10A, and correspondingly fesc(Lya)<1%. By comparing these results with literature, we conclude that the bulk of the Lya luminosity at 2<z<6 comes from galaxies that are fainter in the UV than those we sample in this work. The strong Lya emitters constitute, at each redshift, the tail of the distribution of the galaxies with extreme EW0(Lya) and fesc(Lya) . This tail of large EW0 and fesc(Lya) becomes more important as the redshift increases, and causes the fraction of Lya with EW0> 25A to increase from 5% at z=2 to 30% at z=6, with the increase being relatively stronger beyond z=4. We observe no difference, for the narrow range of UV luminosities explored in this work, between the fraction of strong Lya emitters among galaxies fainter or brighter than M*, although the fraction for the FUV faint galaxies evolves faster, at 2<z<3.5, than for the bright ones. We do observe an anticorrelation between E(B-V) and fesc(Lya): generally galaxies with high fesc(Lya) have also small amounts of dust (and viceversa). However, when the dust content is low (E(B-V)<0.05) we observe a very broad range of fesc(Lya), ranging from 10^-3 to 1. This implies that the dust alone is not the only regulator of the amount of escaping Lya photons.

CosmoMC is a Fortran 95 Markov-Chain Monte-Carlo (MCMC) engine to explore the cosmological parameter space, plus a Python suite for plotting and presenting results (see http://cosmologist.info/cosmomc/). This document describes the installation of the CosmoMC on a Linux system (Ubuntu 14.04.1 LTS 64-bit version). It is written for those who want to use it in their scientific research but without much training on Linux and the program. Besides a step-by-step installation guide, we also give a brief introduction of how to run the program on both a desktop and a cluster. We share our way to generate the plots that are commonly used in the references of cosmology. For more information, one can refer to the CosmoCoffee forum (http://cosmocoffee.info/viewforum.php?f=11) or contact the authors of this document. Questions and comments would be much appreciated.

[Abridged] We characterise a massive proto-cluster at z=2.895 that we found in the COSMOS field using the spectroscopic sample of the VIMOS Ultra-Deep Survey (VUDS). This is one of the rare structures at z~3 not identified around AGNs or radio galaxies, so it is an ideal laboratory to study galaxy formation in dense environments. The structure comprises 12 galaxies with secure spectroscopic redshift in an area of 7'x8', in a z bin of Dz=0.016. The measured galaxy number overdensity is delta_g=12+/-2. This overdensity has total mass of M~8.1x10^(14)M_sun in a volume of 13x15x17 Mpc^3. Simulations indicate that such an overdensity at z~2.9 is a proto-cluster that will collapse in a cluster of total mass M~2.5x10^(15)M_sun at z=0. We compare the properties of the galaxies within the overdensity with a control sample at the same z but outside the overdensity. We did not find any statistically significant difference between the properties (stellar mass, SFR, sSFR, NUV-r, r-K) of the galaxies inside and outside the overdensity. The stacked spectrum of galaxies in the overdensity background shows a significant absorption feature at the wavelength of Lya redshifted at z=2.895 (lambda=4736 A), with a rest frame EW = 4+/- 1.4 A. Stacking only background galaxies without intervening sources at z~2.9 along their line of sight, we find that this absorption feature has a rest frame EW of 10.8+/-3.7 A, with a detection S/N of ~4. These EW values imply a high column density (N(HI)~3-20x10^(19)cm^(-2)), consistent with a scenario where such absorption is due to intervening cold gas streams, falling into the halo potential wells of the proto-cluster galaxies. However, we cannot exclude the hypothesis that this absorption is due to the diffuse gas within the overdensity.

The propagation of polarized emission in pulsar magnetosphere is investigated in this paper. The polarized waves are generated through curvature radiation from the relativistic particles streaming along curved magnetic field lines and co-rotating with the pulsar magnetosphere. Within the 1/{\deg} emission cone, the waves can be divided into two natural wave mode components, the ordinary (O) mode and the extraord nary (X) mode, with comparable intensities. Both components propagate separately in magnetosphere, and are aligned within the cone by adiabatic walking. The refraction of O-mode makes the two components separated and incoherent. The detectable emission at a given height and a given rotation phase consists of incoherent X-mode and O-mode components coming from discrete emission regions. For four particle-density models in the form of uniformity, cone, core and patches, we calculate the intensities for each mode numerically within the entire pulsar beam. If the co-rotation of relativistic particles with magnetosphere is not considered, the intensity distributions for the X-mode and O-mode components are quite similar within the pulsar beam, which causes serious depolarization. However, if the co-rotation of relativistic particles is considered, the intensity distributions of the two modes are very different, and the net polarization of out-coming emission should be significant. Our numerical results are compared with observations, and can naturally explain the orthogonal polarization modes of some pulsars. Strong linear polarizations of some parts of pulsar profile can be reproduced by curvature radiation and subsequent propagation effect.

We present the VIMOS Ultra Deep Survey (VUDS), a spectroscopic redshift survey of ~10.000 very faint galaxies to study the major phase of galaxy assembly 2<z<~6. The survey covers 1 deg^2 in 3 separate fields: COSMOS, ECDFS and VVDS-02h, with targets selection based on an inclusive combination of photometric redshifts and color properties. Spectra covering 3650<lambda<9350 A are obtained with VIMOS on the ESO-VLT with integration times of 14h. Here we present the survey strategy, the target selection, the data processing, as well as the redshift measurement process, emphasizing the specific methods adapted to this high redshift range. The spectra quality and redshift reliability are discussed, and we derive a completeness in redshift measurement of 91%, or 74% for the most reliable measurements, down to i_AB=25, and measurements are performed all the way down to i_AB=27. The redshift distribution of the main sample peaks at z=3-4 and extends over a large redshift range mainly in 2 < z < 6. At 3<z<5, the galaxies cover a large range of luminosities -23< M_U < -20.5, stellar mass 10^9 M_sun< M_star < 10^{11} M_sun, and star formation rates 1 M_sun/yr< SFR < 10^3 M_sun/yr. We discuss the spectral properties of galaxies using individual as well as stacked spectra. The comparison between spectroscopic and photometric redshifts as well as color selection demonstrate the effectiveness of our selection scheme. With ~6000 galaxies with reliable spectroscopic redshifts in 2<z<6 expected when complete, this survey is the largest at these redshifts and offers the opportunity for unprecedented studies of the star-forming galaxy population and its distribution in large scale structures during the major phase of galaxy assembly.

Membership identification is the first step to determine the properties of a star cluster. Low-mass members in particular could be used to trace the dynamical history, such as mass segregation, stellar evaporation, or tidal stripping, of a star cluster in its Galactic environment. We identified member candidates with stellar masses $\sim$0.11--2.4 Msun of the intermediate-age Praesepe cluster (M44), by using Pan-STARRS and 2MASS photometry, and PPMXL proper motions. Within a sky area of 3 deg radius, 1040 candidates are identified, of which 96 are new inclusions. Using the same set of selection criteria on field stars led to an estimate of a false positive rate 16%, suggesting 872 of the candidates being true members. This most complete and reliable membership list allows us to favor the BT-Settl model in comparison with other stellar models. The cluster shows a distinct binary track above the main sequence, with a binary frequency of 20--40%, and a high occurrence rate of similar mass pairs. The mass function is consistent with that of the disk population but shows a deficit of members below 0.3 solar masses. A clear mass segregation is evidenced, with the lowest-mass members in our sample being evaporated from this disintegrating cluster.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Beam radii for cone-dominant pulsars follow a power-law relation with frequency, $\vartheta = (\nu/\nu_0)^k+\vartheta_0$, which has not yet well explained in previous works. We study this frequency dependence of beam radius (FDB) for cone-dominant pulsars by using the curvature radiation mechanism. Considering various density and energy distributions of particles in the pulsar open field line region, we numerically simulate the emission intensity distribution across emission height and rotation phase, and get integrated profiles at different frequencies and obtain the FDB curves. For the density model of a conal-like distribution, the simulated profiles always shrink to one component at high frequencies. In the density model with two separated density patches, the profiles always have two distinct components, and the power-law indices $k$ are found to be in the range from -0.1 to -2.5, consistent with observational results. Energy distributions of streaming particles have significant influence on the frequency-dependence behavior. Radial energy decay of particles are necessary to get proper $\vartheta_0$ in models. We conclude that by using the curvature radiation mechanism, the observed frequency dependence of beam radius for the cone-dominant pulsars can only be explained by the emission model of particles in two density patches with a Gaussian energy distribution and a radial energy loss.

The mass assembly of galaxies can proceed through different physical processes. Here we report on the spectroscopic identification of close physical pairs of galaxies at redshifts 2 <z < 4 and discuss the impact of major mergers in shaping galaxies at these early cosmological times. We aim to identify and characterize close physical pairs of galaxies destined to merge and use their properties to infer the contribution of merging processes to the early mass assembly of galaxies. We search for galaxy pairs with a transverse separation r_p < 25 h^{-1} kpc and a velocity difference Delta v < 500 km s^{-1} using the VIMOS VLT Deep Survey (VVDS) and early data from the VIMOS Ultra Deep Survey (VUDS) that comprise a sample of 1111 galaxies with spectroscopic redshifts measurements at redshifts 1.8 < z < 4 in the COSMOS, ECDFS, and VVDS-02h fields. We analyse their spectra and associated visible and near-infrared photometry to assess the main properties of merging galaxies that have an average stellar mass M_{star}=2.3 x 10^{10} M_{sun} at these redshifts. Using the 12 physical pairs found in our sample we obtain a first estimate of the merger fraction at these rdshifts, fmerg ~(15-20)%. The pair separations indicate that these pairs will merge within 1 Gyr, on average, with each producing a more massive galaxy by the time the cosmic star formation peaks at z~1-2. From the average mass ratio between galaxies in the pairs, the stellar mass of the resulting galaxy after merging will be ~60% larger than the most massive galaxy in the pair before merging. We conclude that major merging of galaxy pairs is on-going at 2<z<4 and significantly contributing to the major assembly phase of galaxies at this early epoch.

Extended quintessence is obtained by coupling a normal scalar field to the Ricci scalar defined in the metric formalism. In this paper, we propose a new extended quintessence dark energy by introducing a non-minimal coupling between the quintessence and gravity, but with the Ricci scalar given from the Palatini formalism rather than the metric one. We find that the equation of state of the new extended quintessence can cross the phantom divide line, and moreover, it oscillates around the -1 line. We also show that the universe driven by the new extended quintessence will enter a dark energy dominated de Sitter phase in the future.

We have investigated the three-body force (TBF) effect on the neutron and proton momentum distributions in asymmetric nuclear matter within the framework of the extended Brueckner-Hartree-Fock approach by adopting the $AV18$ two-body interaction plus a microscopic TBF. In asymmetric nuclear matter, it is shown that the neutron and proton momentum distributions become different from their common distribution in symmetric nuclear matter. The predicted depletion of the proton hole states increases while the neutron one decreases as a function of isospin-asymmetry. The TBF effect on the neutron and proton momentum distributions turns out to be negligibly weak at low densities around and below the normal nuclear density. The TBF effect is found to become sizable only at high densities well above the saturation density, and inclusion of the TBF leads to an overall enhancement of the depletion of the neutron and proton Fermi seas.

We confront two types of phantom dark energy potential with observational data. The models we consider are the power-law potential, $V\propto {\phi}^{\mu}$, and the exponential potential, $V\propto \exp({\lambda}{\phi}/{M_P})$. We fit the models to the latest observations from SN-Ia, CMB and BAO, and obtain tight constraints on parameter spaces. Furthermore, we apply the goodness-of-fit and the information criteria to compare the fitting results from phantom models with that from the cosmological constant and the quintessence models presented in our previous work. The results show that the cosmological constant is statistically most preferred, while the phantom dark energy fits slightly better than the quintessence does.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

We consider the curvature emission properties from relativistic particles streaming along magnetic field lines and co-rotating with pulsar magnetosphere. The co-rotation affects the trajectories of the particles and hence the emission properties, especially the polarization. We consider the modification of the particle velocity and acceleration due to the co-rotation. Curvature radiation from a single particle is calculated using the approximation of a circular path to the particle trajectory. Curvature radiation from particles at a given height actually contains the contributions from particles streaming along all the nearby field lines around the tangential point, forming the emission cone of 1/{\gamma}. The polarization patterns from the emission cone are distorted by the additional rotation, more serious for emission from a larger height. Net circular polarization can be generated by the density gradient in the emission cone. For three typical density models in the form of core, cone and patches, we calculate the polarization profiles for emission generated at a given height. We find that the circular polarization could have a single sign or sign reversal, depending on the density gradient along the rotation phase. The polarization profiles of the total curvature radiation from the whole open field line region, calculated by adding the emission from all possible heights, are similar to that from a dominating emission height. The circular polarization of curvature radiation has sign reversals in the patchy emission, while it has a single sign for the core emission, and is negligible for the cone emission.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

We present a search for Galactic dark matter (DM) satellites using the Large Area Telescope (LAT). N-body simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic DM satellites. These satellites could potentially produce gamma rays through the self-annihilation of DM particles. Some DM satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We searched for LAT sources with these characteristics. We found no candidate DM satellites matching these criteria in one year of LAT data and interpreted this result in the context of N-body simulations.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

We confront tracker field quintessence with observational data. The potentials considered in this paper include $V(\phi)\propto\phi^{-\alpha}$, $\exp(M_{p}/\phi)$, $\exp(M_{p}/\phi)-1$, $\exp(\beta M_{p}/\phi)$ and $\exp(\gamma M_{p}/\phi)-1$; while the data come from the latest SN Ia, CMB and BAO observations. Stringent parameter constraints are obtained. In comparison with the cosmological constant via information criteria, it is found that models with potentials $\exp(M_{p}/\phi)$, $\exp(M_{p}/\phi)-1$ and $\exp(\gamma M_{p}/\phi)-1$ are not supported by the current data.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

Lunar occultations (LO) provide a unique combination of high angular resolution and sensitivity at near-infrared wavelenghts. At the ESO Very Large Telescope, it is possible to achieve about 1 milliarcsecond (mas) resolution and detect sources as faint as K$\approx$12\,mag. We have taken advantage of a passage of the Moon over two crowded and reddened regions in the direction of the inner part of the galactic bulge, in order to obtain a high number of occultation light curves over two half nights. Our goal was to detect and characterize new binary systems, and to investigate highly extincted and relatively unknown infrared sources in search of circumstellar shells and similar peculiarities. Our target list included a significant number of very late-type stars, but in fact the majority of the sources was without spectral classification. A LO event requires the sampling of the light curve at millisecond rates in order to permit a detailed study of the diffraction fringes. For this, we used the so-called burst mode of the ISAAC instrument at the Melipal telescope. Our observing efficiency was ultimately limited by overheads for telescope pointing and data storage, to about one event every three minutes. We could record useful light curves for 184 sources. Of these, 24 were found to be binaries or multiples, all previously unknown. The projected separations are as small as 7.5\,mas, and the magnitude differences as high as $\Delta$K=6.5\,mag. Additionally we could establish, also for the first time, the resolved nature of at least two more stars, with indication of circumstellar emission. We could also put upper limits on the angular size of about 165 unresolved stars, an information that combined with previous and future observations will be very helpful in establishing a list of reliable calibrators for long baseline interferometers.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

Magnetic fields are generally expected to increase the characteristic mass of stars formed in stellar clusters, because they tend to increase the effective Jeans mass. We test this expectation using adaptive mesh refinement (AMR) magnetohydrodynamic simulations of cluster formation in turbulent magnetized clumps of molecular clouds, treating stars as accreting sink particles. We find that, contrary to the common expectation, a magnetic field of strength in the observed range decreases, rather than increases, the characteristic stellar mass. It (1) reduces the number of intermediate-mass stars that are formed through direct turbulent compression, because sub-regions of the clump with masses comparable to those of stars are typically magnetically subcritical and cannot be compressed directly into collapse, and (2) increases the number of low-mass stars that are produced from the fragmentation of dense filaments. The filaments result from mass accumulation along the field lines. In order to become magnetically supercritical and fragment, the filament must accumulate a large enough column density (proportional to the field strength), which yields a high volume density (and thus a small thermal Jeans mass) that is conducive to forming low-mass stars. We find, in addition, that the characteristic stellar mass is reduced further by outflow feedback. The conclusion is that both magnetic fields and outflow feedback are important in shaping the stellar initial mass function (IMF).

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

We describe an implementation of compressible inviscid fluid solvers with block-structured adaptive mesh refinement on Graphics Processing Units using NVIDIA's CUDA. We show that a class of high resolution shock capturing schemes can be mapped naturally on this architecture. Using the method of lines approach with the second order total variation diminishing Runge-Kutta time integration scheme, piecewise linear reconstruction, and a Harten-Lax-van Leer Riemann solver, we achieve an overall speedup of approximately 10 times faster execution on one graphics card as compared to a single core on the host computer. We attain this speedup in uniform grid runs as well as in problems with deep AMR hierarchies. Our framework can readily be applied to more general systems of conservation laws and extended to higher order shock capturing schemes. This is shown directly by an implementation of a magneto-hydrodynamic solver and comparing its performance to the pure hydrodynamic case. Finally, we also combined our CUDA parallel scheme with MPI to make the code run on GPU clusters. Close to ideal speedup is observed on up to four GPUs.

(Abridged) We investigate massive star formation in turbulent, magnetized, parsec-scale clumps of molecular clouds including protostellar outflow feedback using Enzo-based MHD simulations with accreting sink particles and effective resolution $2048^3$. We find that, in the absence of regulation by magnetic fields and outflow feedback, massive stars form readily in a turbulent, moderately condensed clump of $\sim 1,600$ solar masses, along with a cluster of hundreds of lower mass stars. The massive stars are fed at high rates by (1) transient dense filaments produced by large-scale turbulent compression at early times, and (2) by the clump-wide global collapse resulting from turbulence decay at late times. In both cases, the bulk of the massive star's mass is supplied from outside a 0.1 pc-sized "core" that surrounds the star. In our simulation, the massive star is clump-fed rather than core-fed. The need for large-scale feeding makes the massive star formation prone to regulation by outflow feedback, which directly opposes the feeding processes. The outflows reduce the mass accretion rates onto the massive stars by breaking up the dense filaments that feed the massive star formation at early times, and by collectively slowing down the global collapse that fuel the massive star formation at late times. The latter is aided by a moderate magnetic field of strength in the observed range. We conclude that the massive star formation in our simulated turbulent, magnetized, parsec-scale clump is outflow-regulated and clump-fed (ORCF for short). An important implication is that the formation of low-mass stars in a dense clump can affect the formation of massive stars in the same clump, through their outflow feedback on the clump dynamics.

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

Magnetic field strengths inferred for relativistic outflows including gamma-ray bursts (GRB) and active galactic nuclei (AGN) are larger than naively expected by orders of magnitude. We present three-dimensional relativistic magnetohydrodynamics (MHD) simulations demonstrating amplification and saturation of magnetic field by a macroscopic turbulent dynamo triggered by the Kelvin-Helmholtz shear instability. We find rapid growth of electromagnetic energy due to the stretching and folding of field lines in the turbulent velocity field resulting from non-linear development of the instability. Using conditions relevant for GRB internal shocks and late phases of GRB afterglow, we obtain amplification of the electromagnetic energy fraction to $\epsilon_B \sim 5 \times 10^{-3}$. This value decays slowly after the shear is dissipated and appears to be largely independent of the initial field strength. The conditions required for operation of the dynamo are the presence of velocity shear and some seed magnetization both of which are expected to be commonplace. We also find that the turbulent kinetic energy spectrum for the case studied obeys Kolmogorov's 5/3 law and that the electromagnetic energy spectrum is essentially flat with the bulk of the electromagnetic energy at small scales.

We investigate the sensitivity of the Gamma-ray Large Area Space Telescope (GLAST) to indirectly detect weakly interacting massive particles (WIMPs) through the $\gamma$-ray signal that their pair annihilation produces. WIMPs are among the favorite candidates to explain the compelling evidence that about 80% of the mass in the Universe is non-baryonic dark matter (DM). They are serendipitously motivated by various extensions of the standard model of particle physics such as Supersymmetry and Universal Extra Dimensions (UED). With its unprecedented sensitivity and its very large energy range (20 MeV to more than 300 GeV) the main instrument on board the GLAST satellite, the Large Area Telescope (LAT), will open a new window of discovery. As our estimates show, the LAT will be able to detect an indirect DM signature for a large class of WIMP models given a cuspy profile for the DM distribution. Using the current state of the art Monte Carlo and event reconstruction software developed within the LAT collaboration, we present preliminary sensitivity studies for several possible sources inside and outside the Galaxy. We also discuss the potential of the LAT to detect UED via the electron/positron channel. Diffuse background modeling and other background issues that will be important in setting limits or seeing a signal are presented.

Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using both spatially and temporally structured adaptive mesh refinement (AMR). We used the method of lines to discretize the SRHD equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme for time integration. For spatial reconstruction, we have implemented piecewise linear method (PLM), piecewise parabolic method (PPM), third order convex essentially non-oscillatory (CENO) and third and fifth order weighted essentially non-oscillatory (WENO) schemes. Flux is computed using either direct flux reconstruction or approximate Riemann solvers including HLL, modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR part of the code is built on top of the cosmological Eulerian AMR code {\sl enzo}. We discuss the coupling of the AMR framework with the relativistic solvers. Via various test problems, we emphasize the importance of resolution studies in relativistic flow simulations because extremely high resolution is required especially when shear flows are present in the problem. We also present the results of two 3d simulations of astrophysical jets: AGN jets and GRB jets. Resolution study of those two cases further highlights the need of high resolutions to calculate accurately relativistic flow problems.

Using magnetohydrodynamic (MHD) adaptive mesh refinement simulations, we study the formation and early evolution of disk galaxies with a magnetized interstellar medium. For a $10^{10}$ \msun halo with initial NFW dark matter and gas profiles, we impose a uniform $10^{-9}$ G magnetic field and follow its collapse, disk formation and evolution up to 1 Gyr. Comparing to a purely hydrodynamic simulation with the same initial condition, we find that a protogalactic field of this strength does not significantly influence the global disk properties. At the same time, the initial magnetic fields are quickly amplified by the differentially rotating turbulent disk. After the initial rapid amplification lasting $\sim500$ Myr, subsequent field amplification appears self-regulated. As a result, highly magnetized material begin to form above and below the disk. Interestingly, the field strengths in the self-regulated regime agrees well with the observed fields in the Milky Way galaxy both in the warm and the cold HI phase and do not change appreciably with time. Most of the cold phase shows a dispersion of order ten in the magnetic field strength. The global azimuthal magnetic fields reverse at different radii and the amplitude declines as a function of radius of the disk. By comparing the estimated star formation rate (SFR) in hydrodynamic and MHD simulations, we find that after the magnetic field strength saturates, magnetic forces provide further support in the cold gas and lead to a decline of the SFR.

In a hierarchical picture of galaxy formation virialization continually transforms gravitational potential energy into kinetic energies of the baryonic and dark matter. For the gaseous component the kinetic, turbulent energy is transformed eventually into internal thermal energy through shocks and viscous dissipation. Traditionally this virialization and shock heating has been assumed to occur instantaneously allowing an estimate of the gas temperature to be derived from the the virial temperature defined from the embedding dark matter halo velocity dispersion. As the mass grows the virial temperature of a halo grows. Mass accretion hence can be translated into a heating term. We derive this heating rate from the extended Press Schechter formalism and demonstrate its usefulness in semi-analytical models of galaxy formation. Our method explicitly conserves energy unlike the previous impulsive heating assumptions. Our formalism can trivially be applied in all current semi-analytical models as the heating term can be computed directly from the underlying merger trees. Our analytic results for the first cooling halos and the transition from cold to hot accretion are in agreement with numerical simulations.

The Gamma-ray Large Area Space Telescope (GLAST) is an observatory designed to perform gamma-ray astronomy in the energy range 20 MeV to 300 GeV, with supporting measurements for gamma-ray bursts from 10 keV to 25 MeV. GLAST will be launched at the end of 2007, opening a new and important window on a wide variety of high energy astrophysical phenomena . The main instrument of GLAST is the Large Area Telescope (LAT), which provides break-through high-energy measurements using techniques typically used in particle detectors for collider experiments. The LAT consists of 16 identical towers in a four-by-four grid, each one containing a pair conversion tracker and a hodoscopic crystal calorimeter, all covered by a segmented plastic scintillator anti-coincidence shield. The scientific return of the instrument depends very much on how accurately we know its performance, and how well we can monitor it and correct potential problems promptly. We report on a novel technique that we are developing to help in the characterization and monitoring of LAT by using the power of classification trees to pinpoint in a short time potential problems in the recorded data. The same technique could also be used to evaluate the effect on the overall LAT performance produced by potential instrumental problems.

We present a new radiative transfer method Short-N Characteristics that is a hybrid between the standard long characteristic methods and the standard short characteristic methods. We have implemented the numerical method within the SPECT3D imaging and spectral analysis application. Numerical experiments have shown that the short-N method is capable of reproducing the accuracy of the long characteristic methods, while delivering the CPU time efficiency of the short characteristic methods. We apply the short-N method to 2.5D cylindrical coordinates, i.e., we assume a 3D problem under azimuthal symmetry leading to a 2D mathematical problem with r-z as the independent spatial coordinates (2DRZ). We are currently working to extend this new method to other geometries.

We discuss the issue of energy nonconservation in the virialzation process of spherical collapse model with homogeneous dark energy. We propose an approximation scheme to find the virialization radius. By comparing various schemes and estimating the parameter characterizing the ratio of dark energy to dark matter at the turn-around time, we conclude that the problem of energy nonconservation may have sizable effects in fitting models to observations.

We aim at the construction of dark energy models without exotic matter but with a phantom-like equation of state (an effective phantom phase). The first model we consider is decaying vacuum cosmology where the fluctuations of the vacuum are taken into account. In this case, the phantom cosmology (with an effective, observational $\omega$ being less than -1) emerges even for the case of a real dark energy with a physical equation of state parameter $\omega$ larger than -1. The second proposal is a generalized holographic model, which is produced by the presence of an infrared cutoff. It also leads to an effective phantom phase, which is not a transient one as in the first model. However, we show that quantum effects are able to prevent its evolution towards a Big Rip singularity.

We take a phenomenological approach to study the cosmological evolution of decaying vacuum cosmology ($\Lambda(t)$CDM) based on a simple assumption about the form of the modified matter expansion rate. In this framework, almost all the current vacuum decaying models can be unified in a simple manner. We argue that the idea of letting vacuum decay to resolve the fine-tuning problem is inconsistent with cosmological observations. We also discuss some issues in confronting $\Lambda(t)$CDM with observation. Using the effective equation of state formalism, we indicate that $\Lambda(t)$CDM is a possible candidate for phantom cosmology. Moreover, confronted with a possible trouble of effective equation of state formalism, we construct the effect dark energy density. Finally, we discuss the evolution of linear perturbation.

In this paper we derive the Modified Friedmann equation in the Palatini formulation of $R^2$ gravity. Then we use it to discuss the problem of whether in Palatini formulation a $R^2$ term can drive an inflation. We show that the Palatini formulation of $R^2$ gravity cannot lead to gravity-driven inflation. If considering no zero radiation and matter energy densities, we show that only under rather restrictive assumption about the radiation and matter energy densities there will be a mild power-law inflation $a\sim t^2$.

We study the cosmological dynamics of an effective theory for a strongly coupled scalar field in the moduli space of $\mathcal{N}=4$ supersymmetric Yang-Mills theory recently proposed by Silverstein and Tong, called "D-cceleration". We discuss various Energy Conditions in this theory. Then we prove the inflationary attractor property using the Hamilton-Jacobi method and study the phase portrait as well as the cosmological evolution of the scalar field.

We review the Palatini formulation of the higher-derivative gravity of the $L(R)$ form and its applications in cosmology.

The inflationary attractor properties of the canonical scalar field and Born-Infeld field are investigated in the Randall-Sundrum II scenario with a Gauss-Bonnet term in the bulk action. We find that the inflationary attractor property will always hold for both the canonical and Born-Infeld fields for any allowed non-negative Gauss-Bonnet coupling. We also briefly discuss the possibility of explaining the suppressed lower multiples and running scalar spectral index simultaneously in the scenario of Gauss-Bonnet brane inflation.

Recently, corrections of the $L(R)$ type to Einstein-Hilbert action that become important at small curvature are proposed. Those type of models intend to explain the observed cosmic acceleration without dark energy. We derive the full Modified Friedmann equation in the Palatini formulation of those modified gravity model of the $L(R)$ type. Then, we discuss various cosmological predictions of the Modified Friedmann equation.

Recently, corrections to Einstein-Hilbert action that become important at small curvature are proposed. We discuss the first order and second order approximations to the field equations derived by the Palatini variational principle. We work out the first and second order Modified Friedmann equations and present the upper redshift bounds when these approximations are valid. We show that the second order effects can be neglected on the cosmological predictions involving only the Hubble parameter itself, e.g. the various cosmological distances, but the second order effects can not be neglected in the predictions involving the derivatives of the Hubble parameter. Furthermore, the Modified Friedmann equations fit the SN Ia data at an acceptable level.