The uncertainty of cosmological data generated from complex processes, such as observational Hubble parameter data (OHD) and the Type Ia supernova (SN Ia) data, cannot be accurately modeled by simple analytical probability distributions, e.g. Gaussian distribution. This necessitates the use of likelihood-free inference, which bypasses the direct calculation of likelihood. In this paper, we propose a new procedure to perform likelihood-free cosmological inference using two artificial neural networks (ANN), the Masked Autoregressive Flow (MAF) and denoising autoencoder (DAE). Our procedure is the first to use DAE to extract features from data in order to simplify the structure of MAF needed to estimate the posterior. Tested on simulated Hubble parameter data with a simple Gaussian likelihood, the procedure shows the capability of extracting feature from data and estimating posteriors without the need for tractable likelihood. We demonstrate that it can accurately approximate the real posterior, and achieves a performance comparable to the traditional MCMC method. We also discuss the application of the proposed procedure on OHD and Pantheon SN Ia data, and use them to constrain cosmological parameters from the non-flat $\Lambda$CDM model. With little prior knowledge added, we found constraints on $H_0, \Omega_m, \Omega_\Lambda$ similar to relevant work. In addition, this work is also the first to use Gaussian process in the procedure of OHD simulation.

Galaxy spins can be predicted from the initial conditions in the early Universe through the tidal tensor twist. In simulations, their directions are well preserved through cosmic time, consistent with expectations of angular momentum conservation. We report a $\sim 3 \sigma$ detection of correlation between observed oriented directions of galaxy angular momenta and their predictions based on the initial density field reconstructed from the positions of SDSS galaxies. The detection is driven by a group of spiral galaxies classified by the Galaxy Zoo as (anti-)clockwise, with a modest improvement from adding galaxies from MaNGA and SAMI surveys. This is the first such detection of the oriented galaxy spin direction, which opens a way to use measurements of galaxy spins to probe fundamental physics in the early Universe.

Polycyclic aromatic hydrocarbon (PAH) emission has long been proposed to be a potential star formation rate indicator, as it arises from the photodissociation region bordering the Str\"omgren sphere of young, massive stars. We apply a recently developed technique of mid-infrared spectral decomposition to obtain a uniform set of PAH measurements from Spitzer low-resolution spectra of a large sample of star-forming galaxies spanning a wide range in stellar mass ($M_\star \approx 10^6-10^{11.4}\,M_{\odot}$) and star formation rate ($\sim 0.1- 2000\,M_{\odot}\rm\,yr^{-1}$). High-resolution spectra are also analyzed to measure [Ne II] 12.8 $\mu$m and [Ne III] 15.6 $\mu$m, which effectively trace the Lyman continuum. We present a new relation between PAH luminosity and star formation rate based on the [Ne II] and [Ne III] lines. Calibrations are given for the integrated $5-15$ $\mu$m PAH emission, the individual features at 6.2, 7.7, 8.6, and 11.3 $\mu$m, as well as several mid-infrared bandpasses sensitive to PAH. We confirm that PAH emission is suppressed in low-mass, dwarf galaxies, and we discuss the possible physical origin of this effect.

We present a new method to investigate the effective magnetic field decay of isolated neutron stars, from the analysis of the long-term timing data of a large sample of radio pulsars \citep{2010MNRAS.402.1027H}. There are some differences between the distributions of frequency's second derivatives of the pulsar spins with different effective field decay timescales. Kolmogorov-Smirnov tests are performed to reexamine the consistency of distributions of the simulated and reported data for a series values of decay timescales. \textbf{We show that the timescale of the effective field decay exceeds $\sim 5~\rm {Myr}$ for pulsars with spin-down age $\tau_{\rm C}< 10^{7}~{\rm yr}$ or $\sim 100~\rm {Myr}$ for pulsars with $10^{7}<\tau_{\rm C}< 10^{9}~{\rm yr}$ in the sample. The result does not depend on any specific theories of the field evolution, the inclination decay or the variation in the moment of inertia. It is also found that the extent of the closed line region of the magnetic field is close to the light cylinder $r_{\rm lc}$, i.e., the corotating radius $r_{\rm c}\approx r_{\rm lc}$ is a good approximation for the observed pulsar population.

Flux density measurements provide fundamental observational parameters that describe a pulsar. In the current pulsar catalogue, 27% of radio pulsars have no flux density measurement in the 20 cm observing band. Here, we present the first measurements of the flux densities in this band for 32 pulsars observed using the Parkes radio telescope and provide updated pulse profiles for these pulsars. We have used both archival and new observations to make these measurements. Various schemes exist for measuring flux densities. We show how the flux densities measured vary between these methods and how the presence of radio-frequency-interference will bias flux density measurements

The interstellar medium is a key ingredient that governs star formation in galaxies. We present a detailed study of the infrared (~ 1-500 micron) spectral energy distributions of a large sample of 193 nearby (z ~ 0.088) luminous infrared galaxies (LIRGs) covering a wide range of evolutionary stages along the merger sequence. The entire sample has been observed uniformly by 2MASS, WISE, Spitzer, and Herschel. We perform multi-component decomposition of the spectra to derive physical parameters of the interstellar medium, including the intensity of the interstellar radiation field and the mass and luminosity of the dust. We also constrain the presence and strength of nuclear dust heated by active galactic nuclei. The radiation field of LIRGs tends to have much higher intensity than in quiescent galaxies, and it increases toward advanced merger stages as a result of central concentration of the interstellar medium and star formation. The total gas mass is derived from the dust mass and the galaxy stellar mass. We find that the gas fraction of LIRGs is on average ~ 0.3 dex higher than that of main-sequence star-forming galaxies, rising moderately toward advanced merger stages. All LIRGs have star formation rates that place them above the galaxy star formation main sequence. Consistent with recent observations and numerical simulations, the global star formation efficiency of the sample spans a wide range, filling the gap between normal star-forming galaxies and extreme starburst systems.

Interstellar dust spans a wide range in size distribution, ranging from ultrasmall grains of a few Angstroms to micrometer-size grains. While the presence of nanometer-size dust grains in the Galactic interstellar medium was speculated six decades ago and was previously suggested based on early infrared observations, systematic and direct analysis of their properties over a wide range of environments has been lacking. Here we report the detection of nanometer-size dust grains that appear to be universally present in a wide variety of astronomical environments, from Galactic high-latitude clouds to nearby star-forming galaxies and galaxies with low levels of nuclear activity. The prevalence of such a grain population is revealed conclusively as prominent mid-infrared continuum emission at wavelength short than 10\mu m seen in the Spitzer/IRS data, characterized by temperatures of ~300-400K that are significantly higher than the equilibrium temperatures of common, submicron-size grains in typical galactic environments. We propose that the optimal carriers of this pervasive, featureless hot dust component are very small carbonaceous (e.g., graphite) grains of nanometer size that are transiently heated by single-photon absorption. This grain population accounts for ~1.4% of the total infrared emission at ~5-3000\mu m and ~0.4% of the total interstellar dust mass.

Weak equivalence principle (WEP) is, for the first time, tested by astrometry on quasars in the sky measured in two wavelengths. Compared to previous WEP tests based on the Shapiro time delay of massless particles, this one has profound superiority that nearly 1,700 quasars with best measured positions commonly in the optical and radio bands are available. It ensures that, among the tests with photons, this one can give the most significantly robust bound on possible violation of WEP.

We present a new method to decompose the emission features of polycyclic aromatic hydrocarbons (PAHs) from mid-infrared spectra using theoretical PAH templates in conjunction with modified blackbody components for the dust continuum and an extinction term. The primary goal is to obtain robust measurements of the PAH features, which are sensitive to the star formation rate, in a variety of extragalactic environments. We demonstrate the effectiveness of our technique, starting with the simplest Galactic high-latitude clouds to extragalactic systems of ever-increasing complexity, from normal star-forming galaxies to low-luminosity active galaxies, quasars, and heavily obscured infrared-luminous galaxies. In addition to providing accurate measurements of the PAH emission, including the upper limits thereof, our fits can reproduce reasonably well the overall continuum shape and constrain the line-of-sight extinction. Our new PAH line flux measurements differ systematically and significantly from those of previous methods by ~15% to as much as a factor of ~6. The decomposed PAH spectra show remarkable similarity among different systems, suggesting a uniform set of conditions responsible for their excitation.

The interstellar medium is crucial to understanding the physics of active galaxies and the coevolution between supermassive black holes and their host galaxies. However, direct gas measurements are limited by sensitivity and other uncertainties. Dust provides an efficient indirect probe of the total gas. We apply this technique to a large sample of quasars, whose total gas content would be prohibitively expensive to measure. We present a comprehensive study of the full (1 to 500 micron) infrared spectral energy distributions of 87 redshift <0.5 quasars selected from the Palomar-Green sample, using photometric measurements from 2MASS, WISE, and Herschel, combined with Spitzer mid-infrared (5 to 40 micron) spectra. With a newly developed Bayesian Markov Chain Monte Carlo fitting method, we decompose various overlapping contributions to the integrated spectral energy distribution, including starlight, warm dust from the torus, and cooler dust on galaxy scales. This procedure yields a robust dust mass, which we use to infer the gas mass, using a gas-to-dust ratio constrained by the host galaxy stellar mass. Most (90%) quasar hosts have gas fractions similar to those of massive, star-forming galaxies, although a minority (10%) seem genuinely gas-deficient, resembling present-day massive early-type galaxies. This result indicates that "quasar mode" feedback does not occur or is ineffective in the host galaxies of low-redshift quasars. We also find that quasars can boost the interstellar radiation field and heat dust on galactic scales. This cautions against the common practice of using the far-infrared luminosity to estimate the host galaxy star formation rate.

The large majority of the accreting supermassive black holes in the Universe are obscured by large columns of gas and dust. The location and evolution of this obscuring material have been the subject of intense research in the past decades, and are still highly debated. A decrease in the covering factor of the circumnuclear material with increasing accretion rates has been found by studies carried out across the electromagnetic spectrum. The origin of this trend has been suggested to be driven either by the increase in the inner radius of the obscuring material with incident luminosity due to the sublimation of dust; by the gravitational potential of the black hole; by radiative feedback; or by the interplay between outflows and inflows. However, the lack of a large, unbiased and complete sample of accreting black holes, with reliable information on gas column density, luminosity and mass, has left the main physical mechanism regulating obscuration unclear. Using a systematic multi-wavelength survey of hard X-ray-selected black holes, here we show that radiation pressure on dusty gas is indeed the main physical mechanism regulating the distribution of the circumnuclear material. Our results imply that the bulk of the obscuring dust and gas in these objects is located within the sphere of influence of the black hole (i.e., a few to tens of parsecs), and that it can be swept away even at low radiative output rates. The main physical driver of the differences between obscured and unobscured accreting black holes is therefore their mass-normalized accretion rate.

The tori around active galactic nuclei (AGN) are potential formation sites for large amounts of dust, and they may help resolve the so-called dust budget crisis at high redshift. We investigate the dust composition in 53 of the 87 Palomar Green (PG) quasars showing the 9.7 micron silicate feature in emission. By simultaneously fitting the mid-infrared spectroscopic features and the underlying continuum, we estimate the mass fraction in various amorphous and crystalline dust species. We find that the dust consists predominantly of alumina and amorphous silicates, with a small fraction in crystalline form. The mean crystallinity is 8 +/- 6%, with more than half of the crystallinities greater than 5%, well above the upper limit determined for the Galaxy. Higher values of crystallinity are found for higher oxide fractions and for more luminous sources.

The unification theory of active galactic nuclei (AGNs) hypothesizes that all AGNs are surrounded by an anisotropic dust torus and are essentially the same objects but viewed from different angles. However, little is known about the dust which plays a central role in the unification theory. There are suggestions that the AGN dust extinction law appreciably differs from that of the Galaxy. Also, the silicate emission features observed in type 1 AGNs appear anomalous (i.e., their peak wavelengths and widths differ considerably from that of the Galaxy). In this work, we explore the dust properties of 147 AGNs of various types at redshifts z<0.5, with special attention paid to 93 AGNs which exhibit the 9.7 and 18 $\mu$m silicate emission features. We model their silicate emission spectra obtained with the Infrared Spectrograph aboard the Spitzer Space Telescope. We find that 60/93 of the observed spectra can be well explained with "astronomical silicate", while the remaining sources favor amorphous olivine or pyroxene. Most notably, all sources require the dust to be $\mu$m-sized (with a typical size of ~1.5$\pm$0.1 $\mu$m), much larger than sub-$\mu$m-sized Galactic interstellar grains, implying a flat or "gray" extinction law for AGNs. We also find that, while the 9.7 $\mu$m emission feature arises predominantly from warm silicate dust of temperature T~270 K, the ~5--8 $\mu$m continuum emission is mostly from carbon dust of T~640 K. Finally, the correlations between the dust properties (e.g., mass, temperature) and the AGN properties (e.g., luminosity, black hole mass) have also been investigated.

We investigate the dust properties in three spectroscopically anomalous galaxies (IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948). Their Spitzer/IRS spectra are characterized by a steep ~5-8 micron emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent 10 micron silicate emission. The steep ~5-8 micron continuum and strong PAH emission features suggest the presence of starbursts, while the silicate emission is indicative of significant heating from AGNs. The simultaneous detection of these two observational properties has rarely been reported on galactic scale. We employ the PAHFIT software to estimate their starlight contributions, and the CLUMPY model for the components contributed by the AGN tori. We find that the CLUMPY model is generally successful in explaining the overall dust infrared emission, although it appears to emit too flat at the ~5-8 micron continuum to be consistent with that observed in IRAS F10398+1455 and IRAS F21013-0739. The flat ~5-8 micron continuum calculated from the CLUMPY model could arise from the adopted specific silicate opacity of Ossenkopf et al. (1992) which exceeds that of the Draine & Lee (1984) "astronomical silicate" by a factor up to 2 in the ~5-8 micron wavelength range. Future models with a variety of dust species incorporated in the CLUMPY radiation transfer regime are needed for a thorough understanding of the dust properties of these spectroscopically anomalous galaxies.

We explain some phenomena existing generally in the timing residuals: amplitude and sign of the second derivative of a pulsar's spin-frequency ($\ddot\nu$), some sophisticated residual patterns, which also change with the time span of data segments. The sample is taken from Hobbs et al.\,(2010), in which the pulsar's spin-frequency and its first derivative have been subtracted from the timing solution fitting. We first classify the timing residual patterns into different types based on the sign of $\ddot\nu$. Then we use the magnetic field oscillation model developed in our group \citep{zhang12a} to fit successfully the different kinds of timing residuals with the Markov Chain Monte Carlo method. Finally, we simulate the spin evolution over 20 years for a pulsar with typical parameters and analyze the data with the conventional timing solution fitting. By choosing different segments of the simulated data, we find that most of the observed residual patterns can be reproduced successfully. This is the first time that the observed residual patterns are fitted by a model and reproduced by simulations with very few parameters. From the distribution of the different residual patterns in the $P-\dot P$ diagram, we argue that (1) a single magnetic field oscillation mode exists commonly in all pulsars throughout their lifetimes; (2) there may be a transition period over the lifetimes of pulsars, in which multiple magnetic field oscillation modes exist.

The \textit{Spitzer}/Infrared Spectrograph spectra of three spectroscopically anomalous galaxies (IRAS~F10398+1455, IRAS~F21013-0739 and SDSS~J0808+3948) are modeled in terms of a mixture of warm and cold silicate dust, and warm and cold carbon dust. Their unique infrared (IR) emission spectra are characterized by a steep $\simali$5--8$\mum$ emission continuum, strong emission bands from polycyclic aromatic hydrocarbon (PAH) molecules, and prominent silicate emission. The steep $\simali$5--8$\mum$ emission continuum and strong PAH emission features suggest the dominance of starbursts, while the silicate emission is indicative of significant heating from active galactic nuclei (AGNs). With warm and cold silicate dust of various compositions ("astronomical silicate," amorphous olivine, or amorphous pyroxene) combined with warm and cold carbon dust (amorphous carbon, or graphite), we are able to closely reproduce the observed IR emission of these %spectroscopically anomalous galaxies. We find that the dust temperature is the primary cause in regulating the steep $\sim$5--8$\mum$ continuum and silicate emission, insensitive to the exact silicate or carbon dust mineralogy and grain size $a$ as long as $a\simlt1\mum$. More specifically, the temperature of the $\simali$5--8$\mum$ continuum emitter (which is essentially carbon dust) of these galaxies is $\sim$250--400$\K$, much lower than that of typical quasars which is $\sim$640$\K$. Moreover, it appears that larger dust grains are preferred in quasars. The lower dust temperature and smaller grain sizes inferred for these three galaxies compared with that of quasars could be due to the fact that they may harbor a young/weak AGN which is not maturely developed yet.

We model the evolution of the spin frequency's second derivative $\ddot\nu$ and the braking index $n$ of radio pulsars with simulations within the phenomenological model of their surface magnetic field evolution, which contains a long-term power-law decay modulated by short-term oscillations. For the pulsar PSR B0329+54, a model with three oscillation components can reproduce its $\ddot\nu$ variation. We show that the "averaged" $n$ is different from the instantaneous $n$, and its oscillation magnitude decreases abruptly as the time span increases, due to the "averaging" effect. The simulated timing residuals agree with the main features of the reported data. Our model predicts that the averaged $\ddot\nu$ of PSR B0329+54 will start to decrease rapidly with newer data beyond those used in Hobbs et al.. We further perform Monte Carlo simulations for the distribution of the reported data in $|\ddot\nu|$ and $|n|$ versus characteristic age $\tau_{\rm c}$ diagrams. It is found that the magnetic field oscillation model with decay index $\alpha=0$ can reproduce the distributions quite well. Compared with magnetic field decay due to the ambipolar diffusion ($\alpha=0.5$) and the Hall cascade ($\alpha=1.0$), the model with no long term decay ($\alpha=0$) is clearly preferred for old pulsars by the p-values of the two-dimensional Kolmogorov-Smirnov test.

On a galactic scale the 9.7um silicate emission is usually only seen in type 1 active galactic nuclei (AGNs). They usually also display a flat emission continuum at ~5--8um and the absence of polycyclic aromatic hydrocarbon (PAH) emission bands. In contrast, starburst galaxies, luminous infrared (IR) galaxies (LIRGs), and ultraluminous IR galaxies (ULIRGs) exhibit a red 5--8um emission continuum, strong 10um and 18um silicate absorption features, and strong PAH emission bands. Here we report the detection of anomalous dust properties by Spitzer/Infrared Spectrograph in three galaxies (IRASF10398+1455, IRASF21013-0739 and SDSSJ0808+3948) which are characterized by the simultaneous detection of a red 5--8um emission continuum, the 9.7 and 18um silicate emission features as well as strong PAH emission bands. These apparently contradictory dust IR emission properties are discussed in terms of iron-poor silicate composition, carbon dust deficit, small grain size and low dust temperature in the young AGN phase of these three galaxies.

The observed long-term spin-down evolution of isolated radio pulsars cannot be explained by the standard magnetic dipole radiation with a constant braking torque. However, how and why the torque varies still remains controversial, which is a major issue in understanding neutron stars. Many pulsars have been observed with significant long-term changes of their spin-down rates modulated by quasi-periodic oscillations. Applying the phenomenological model of pulsar timing noise we developed recently to the observed precise pulsar timing data of isolated neutron stars, here we show that the observed long-term evolutions of their spin-down rates and quasi-periodic modulations can be explained by Hall effects in their crusts. Therefore the evolution of their crustal magnetic fields, rather than that in their cores, dominates the observed long term spin-down evolution of these young pulsars. Understanding of the nature of pulsar timing noise not only reveals the interior physics of neutron stars, but also allows physical modeling of pulsar spin-down and thus improves the sensitivity of gravitational wave detections with pulsars.

Ground-based and space-borne observatories studying exoplanetary transits now and in the future will considerably increase the number of known exoplanets and the precision of the measured times of transit minima. Variations in the transit times can not only be used to infer the presence of additional planets, but might also provide opportunities for testing new physics in the places beyond the Solar system. In this work, we take deviation from the inverse-square law of gravity as an example, focus on the fifth-force-like Yukawa-type correction to the Newtonian gravitational force which parameterizes this deviation, investigate its effects on the secular transit timing variations and analyze their observability in exoplanetary systems. It is found that the most optimistic values of Yukawa-type secular transit timing variations are at the level of $\sim 0.1$ seconds per year. Those values unfortunately appear only in rarely unique cases and, most importantly, they are still at least two orders of magnitude below the current capabilities of observations. Such a deviation from the inverse-square law of gravity is likely too small to detect for the foreseeable future. Meanwhile, systematic uncertainties, such as the presence of additional and unknown planets, will likely be exceptionally difficult to remove from a signal that should be seen.

The observed long-term spin-down evolution of isolated radio pulsars cannot be explained by the standard magnetic dipole radiation with a constant braking torque. However how and why the torque varies still remains controversial, which is an outstanding problem in our understanding of neutron stars. We have constructed a phenomenological model of the evolution of surface magnetic fields of pulsars, which contains a long-term decay modulated by short-term oscillations; a pulsar's spin is thus modified by its magnetic field evolution. The predictions of this model agree with the precisely measured spin evolutions of several individual pulsars; the derived parameters suggest that the Hall drift and Hall waves in the NS crusts are probably responsible for the long-term change and short-term quasi-periodical oscillations, respectively. Many statistical properties of the timing noise of pulsars can be well re-produced with this model, including correlations and the distributions of the observed braking indices of the pulsars, which span over a range of more than 100 millions. We have also presented a phenomenological model for the recovery processes of classical and slow glitches, which can successfully model the observed slow and classical glitch events without biases.

We study the pulsar timing properties and the data analysis methods during glitch recoveries. In some cases one first fits the time-of-arrivals (TOAs) to obtain the "time-averaged" frequency $\nu$ and its first derivative $\dot\nu$, and then fits models to them. However, our simulations show that $\nu$ and $\dot\nu$ obtained this way are systematically biased, unless the time intervals between the nearby data points of TOAs are smaller than about $10^4$ s, which is much shorter than typical observation intervals. Alternatively, glitch parameters can be obtained by fitting the phases directly with relatively smaller biases; but the initial recovery timescale is usually chosen by eyes, which may introduce a strong bias. We also construct a phenomenological model by assuming a pulsar's spin-down law of $\dot\nu\nu^{-3} =-H_0 G(t)$ with $G(t)=1+\kappa e^{-t/\tau}$ for a glitch recovery, where $H_0$ is a constant and $\kappa$ and $\tau$ are the glitch parameters to be found. This model can reproduce the observed data of slow glitches from B1822--09 and a giant classical glitch of B2334+61, with $\kappa<0$ or $\kappa>0$, respectively. We then use this model to simulate TOA data and test several fitting procedures for a glitch recovery. The best procedure is: 1) use a very high order polynomial (e.g. to 50th order) to precisely describe the phase; 2) then obtain $\nu(t)$ and $\dot\nu(t)$ from the polynomial; and 3) the glitch parameters are obtained from $\nu(t)$ or $\dot\nu(t)$. Finally, the uncertainty in the starting time $t_0$ of a classical glitch causes uncertainties to some glitch parameters, but less so to a slow glitch and $t_0$ of which can be determined from data.

We model the evolution of spin frequency's second derivative $\ddot\nu$ and braking index $n$ of radio pulsars with simulations within the phenomenological model of their surface magnetic field evolution, which contains a long-term decay modulated by short-term oscillations. For the pulsar PSR B0329+54, the model can reproduce the main characteristics of its $\ddot\nu$ variation with oscillation periods, predicts another $\sim 50$ yr oscillation component and another recent swing of the sign of $\ddot\nu$. We show that the "averaged" $n$ is different from the instantaneous $n$, and its oscillation magnitude decreases abruptly as the time span increases, due to the "averaging" effect. The simulation predicted timing residuals agree with the main features of the reported data. We further perform Monte Carlo simulations for the distribution of the reported data in $|\ddot\nu|$ versus characteristic age $\tau_{\rm c}$ diagram. The model with a power law index $\alpha=0.5$ can reproduce the slope of the linear fit to pulsars' distributions in the diagrams of $\log|\ddot\nu|-\log\tau_{\rm c}$ and $\log|n|-\log\tau_{\rm c}$, but the oscillations are responsible for the almost equal number of positive and negative values of $\ddot\nu$, in agreement with our previous analytical studies; an oscillation period of about several decades is also preferred. However the range of the oscillation amplitudes is $-11.4\lesssim\log f\lesssim-10.2$, slightly lager than the analytical prediction, $\log f\simeq-11.85$, because the "averaging" effect was not included previously.

Accelerations of both the solar system barycenter (SSB) and stars in the Milky Way cause a systematic observational effect on the stellar proper motions, which was first studied in the early 1990s and developed by J. Kovalevsky (aberration in proper motions, 2003, A&A, 404, 743). This paper intends to extend that work and aims to estimate the magnitude and significance of the aberration in proper motions of stars, especially in the region near the Galactic center. We adopt two models for the Galactic rotation curve to evaluate the aberrational effect on the Galactic plane. Based on the theoretical developments, we show that the effect of aberration in proper motions depends on the galactocentric distance of stars; it is dominated by the acceleration of stars in the central region of the Galaxy. Within 200 pc from the Galactic center, the systematic proper motion can reach an amplitude larger than 1000 uas/yr by applying a flat rotation curve. With a more realistic rotation curve which is linearly rising in the core region of the Galaxy, the aberrational proper motions are limited up to about 150 uas/yr. Then we investigate the applicability of the theoretical expressions concerning the aberrational proper motions, especially for those stars with short period orbits. If the orbital period of stars is only a fraction of the light time from the star to the SSB, the expression proposed by Kovalevsky is not appropriate. With a more suitable formulation, we found that the aberration has no effect on the determination of the stellar orbits on the celestial sphere. The aberrational effect under consideration is small but not negligible with high-accurate astrometry in the future, particularly in constructing the Gaia celestial reference system realized by Galactic stars.

The Yukawa correction to the Newtonian gravitational force is accepted as a parameterization of deviations from the inverse-square law of gravity which might be caused by new physics beyond the standard model of particles and the general theory of relativity. We investigate these effects on the clock onboard a drag-free satellite: dynamics of the satellite and influence on the time transfer link. We find the Yukawa signal in the time transfer is much harder to detect with current state of clocks than those effects on the dynamics, especially the secular change of periastron, by laser ranging in the case of an artificial Earth satellite carrying a frequency standard with an orbit of $a=10^7$\,m and $e=0.01$.

Here we report that the observed braking indices of the 366 pulsars in the sample of Hobbs et al. range from about $-10^8$ to about $+10^8$ and are significantly correlated with their characteristic ages. Using the model of magnetic field evolution we developed previously based on the same data, we derived an analytical expression for the braking index, which agrees with all the observed statistical properties of the braking indices of the pulsars in the sample of Hobbs et al. Our model is, however, incompatible with the previous interpretation that magnetic field growth is responsible for the small values of braking indices ($<3$) observed for "baby" pulsars with characteristic ages of less than $2\times 10^3$ yr. We find that the "instantaneous" braking index of a pulsar may be different from the "averaged" braking index obtained from fitting the data over a certain time span. The close match between our model-predicted "instantaneous" braking indices and the observed "averaged" braking indices suggests that the time spans used previously are usually smaller than or comparable to their magnetic field oscillation periods. Our model can be tested with the existing data, by calculating the braking index as a function of the time span for each pulsar. In doing so, one can obtain for each pulsar all the parameters in our magnetic field evolution model, and may be able to improve the sensitivity of using pulsars to detect gravitational waves.

We propose that a magnetar could be formed during the core collapse of massive stars or coalescence of two normal neutron stars, through collecting and inheriting the magnetic fields magnified by hyperaccreting disk. After the magnetar is born, its dipole magnetic fields in turn have a major influence on the following accretion. The decay of its toroidal field can fuel the persistent X-ray luminosity of either an SGR or AXP; however the decay of only the poloidal field is insufficient to do so.

Under the standard model extension (SME) framework, Lorentz invariance is tested in five binary pulsars: PSR J0737-3039, PSR B1534+12, PSR J1756-2251, PSR B1913+16 and PSR B2127+11C. By analyzing the advance of periastron, we obtain the constraints on a dimensionless combination of SME parameters that is sensitive to timing observations. The results imply no evidence for the break of Lorentz invariance at $10^{-10}$ level, one order of magnitude larger than previous estimation.

We test models for the evolution of neutron star (NS) magnetic fields (B). Our model for the evolution of the NS spin is taken from an analysis of pulsar timing noise presented by Hobbs et al. (2010). We first test the standard model of a pulsar's magnetosphere in which B does not change with time and magnetic dipole radiation is assumed to dominate the pulsar's spin-down. We find this model fails to predict both the magnitudes and signs of the second derivatives of the spin frequencies ($\ddot{\nu}$). We then construct a phenomenological model of the evolution of $B$, which contains a long term decay (LTD) modulated by short term oscillations (STO); a pulsar's spin is thus modified by its B-evolution. We find that an exponential LTD is not favored by the observed statistical properties of $\ddot{\nu}$ for young pulsars and fails to explain the fact that $\ddot{\nu}$ is negative for roughly half of the old pulsars. A simple power-law LTD can explain all the observed statistical properties of $\ddot{\nu}$. Finally we discuss some physical implications of our results to models of the B-decay of NSs and suggest reliable determination of the true ages of many young NSs is needed, in order to constrain further the physical mechanisms of their B-decay. Our model can be further tested with the measured evolutions of $\dot{\nu}$ and $\ddot{\nu}$ for an individual pulsar; the decay index, oscillation amplitude and period can also be determined this way for the pulsar.

We model the evolution of the magnetic fields of neutron stars as consisting of a long term power-law decay modulated by short term small amplitude oscillations. Our model predictions on the timing noise $\ddot\nu$ of neutron stars agree well with the observed statistical properties and correlations of normal radio pulsars. Fitting the model predictions to the observed data, we found that their initial parameter implies their initial surface magnetic dipole magnetic field strength ~ 5E14 G at ~0.4 year old and that the oscillations have amplitude between E-8 to E-5 and period on the order of years. For individual pulsars our model can effectively reduce their timing residuals, thus offering the potential of more sensitive detections of gravitational waves with pulsar timing arrays. Finally our model can also re-produce their observed correlation and oscillations of the second derivative of spin frequency, as well as the "slow glitch" phenomenon.

The observed correlations, between the characteristic ages and dipole surface magnetic field strengths of all pulsars, can be well explained by magnetic field decay with core temperatures of $~2\times10^{8}$ K, $\sim2\times10^{7}$ K, and $\sim10^{5}$ K, for magnetars, normal radio pulsars, and millisecond pulsars, respectively; assuming that their characteristic ages are about two orders of magnitude larger than their true ages, the required core temperatures may be reduced by about a factor of 10. The magnetic decay follows a power-law and is dominated by the solenoidal component of the ambipolar diffusion mode. In this model, all NSs are assumed to have the same initial magnetic field strength, but different core temperature which do not change as the magnetic field decays. This suggests that the key distinguishing property between magnetars and normal pulsars is that magnetars were born much hotter than normal pulsars, and thus have much longer magnetic field decay time scales, resulting in higher surface magnetic field strength even with the same ages of normal pulsars. The above conclusion agrees well with the observed correlations between the surface temperatures of magnetars and other young NSs, which do not agree with the cooling dominated evolution of neutron stars. This suggests a possible scenario that heating, perhaps due to magnetic field decay, balances neutron star cooling for observed pulsars.

It is commonly accepted that a neutron star is produced, when a massive star exhausts its nuclear fuel and ends its life in a core-collapse supernova explosion. This scenario is confirmed by the detection of pulsars, which are believed to be rapidly spinning neutron stars, in the central regions of many supernova remnants. Neutron stars and their associated supernova remnants should therefore have the same ages. As expected, the age of the Crab Pulsar, the first to be connected with a supernova remnant (the Crab Nebula), can be inferred from its current spin period and its derivative and indeed has about the same age of the supernova remnant, that was produced from a historically recorded supernova explosion in 1054. However most neutron stars appear to be much older than the ages of their associated supernova remnants, a puzzle not yet understood. Another puzzle is that so far no convincing evidence has been found in favor of magnetic field decay in neutron stars, that is predicted in most models of neutron stars. Here we show convincing evidence of magnetic field decay in some young neutron stars, and that the magnetic field decay can alter their spinning behaviors significantly such that these neutron stars appear much older than they really are.

We have produced the next generation of quasar spectral energy distributions (SEDs), essentially updating the work of Elvis et al. (1994) by using high-quality data obtained with several space and ground-based telescopes, including NASA's Great Observatories. We present an atlas of SEDs of 85 optically bright, non-blazar quasars over the electromagnetic spectrum from radio to X-rays. The heterogeneous sample includes 27 radio-quiet and 58 radio-loud quasars. Most objects have quasi-simultaneous ultraviolet-optical spectroscopic data, supplemented with some far-ultraviolet spectra, and more than half also have Spitzer mid-infrared IRS spectra. The X-ray spectral parameters are collected from the literature where available. The radio, far-infrared, and near-infrared photometric data are also obtained from either the literature or new observations. We construct composite spectral energy distributions for radio-loud and radio-quiet objects and compare these to those of Elvis et al., finding that ours have similar overall shapes, but our improved spectral resolution reveals more detailed features, especially in the mid and near-infrared.

The luminescent properties of CsI(Na) crystals are studied in this report. By using a TDS3054C oscilloscope with a sampling frequency of 5 GS/s, we find out that nuclear recoil signals are dominated by very fast light pulse with a decay time of ~20 ns, while {\gamma}-ray signals have a decay time of ~600 ns. The wavelength of nuclear recoil and {\gamma}-ray signals are also different. The study of n/{\gamma} separation shows that the rejection factor can reach an order of 10-7 with signal efficiency more than 80% at an equivalent electron recoil energy of 20 keV or more. Such a property makes CsI(Na) an ideal candidate for dark matter searches.

Recent observations and theoretical work on gamma-ray bursts (GRBs) favor the central engine model of a Kerr black hole (BH) surrounded by a magnetized neutrino-dominated accretion flow (NDAF). The magnetic coupling between the BH and disk through a large-scale closed magnetic field exerts a torque on the disk, and transports the rotational energy from the BH to the disk. We investigate the properties of the NDAF with this magnetic torque. For a rapid spinning BH, the magnetic torque transfers enormous rotational energy from BH into the inner disk. There are two consequences: (i) the luminosity of neutrino annihilation is greatly augmented; (ii) the disk becomes thermally and viscously unstable in the inner region, and behaves S-Shape of the surface density versus accretion rate. It turns out that magnetically torqued NDAF can be invoked to interpret the variability of gamma-ray luminosity. In addition, we discuss the possibility of restarting the central engine to produce the X-ray flares with required energy.

Many models of gamma-ray bursts suggest a common central engine: a black hole of several solar masses accreting matter from a disk at an accretion rate from 0.01 to 10 M_\sun / s. The inner region of the disk is cooled by neutrino emission and large amounts of its binding energy were liberated, which could trigger the fireball. We improve the neutrino-dominated accreting flows by considering the effects of the magnetic fields, and find that more than half of the liberating energy can be extracted directly by the large-scale magnetic fields on the disk. And it turns out that the temperature of the disk is a bit lower than the neutrino-dominated accreting flows without magnetic field. Therefore, The outflows are magneticallydominated rather than neutrino dominated. In our model, neutrino mechanism can fuel some GRBs (not the brightest ones), but cannot fuel X-ray flares. However, the magnetic processes (both BZ and electromagnetic luminosity from a disk) are viable mechanisms for most of GRBs and the following X-ray flares.