Observations of galaxy clustering are made in redshift space, which results in distortions to the underlying isotropic distribution of galaxies. These redshift-space distortions (RSD) not only degrade important features of the matter density field, such as the baryonic acoustic oscillation (BAO) peaks, but also pose challenges for the theoretical modelling of observational probes. Here we introduce an iterative nonlinear reconstruction algorithm to remove RSD effects from galaxy clustering measurements, and assess its performance by using mock galaxy catalogues. The new method is found to be able to recover the real-space galaxy correlation function with an accuracy of $\sim1\%$, and restore the quadrupole accurately to $0$, on scales $s\gtrsim20\Mpch$. It also leads to an improvement in the reconstruction of the initial density field, which could help to accurately locate the BAO peaks. An `internal calibration' scheme is proposed to determine the values of cosmological parameters as a part of the reconstruction process, and possibilities to break parameter degeneracies are discussed. RSD reconstruction can offer a potential way to simultaneously extract the cosmological parameters, initial density field, real-space galaxy positions and large-scale peculiar velocity field (of the real Universe), making it an alternative to standard perturbative approaches in galaxy clustering analysis, bypassing the need for RSD modelling.

Since Sep. 2018, LAMOST started the medium-resolution (R $\sim 7500$) spectral survey (MRS). We proposed the spectral survey of Galactic nebulae, including \hii\ regions, HH objects, supernova remnants, planetary nebulae and the special stars with MRS (LAMOST MRS-N). LAMOST MRS-N covers about 1700 square degrees of the northern Galactic plane within $40^\circ < l < 215^\circ$ and $-5^\circ < b < 5^\circ$. In this 5-year survey, we plan to observe about 500 thousand nebulae spectra. According to the commissioning observations, the nebulae spectra can provide precise radial velocity with uncertainty less than 1 km s$^{-1}$. These high-precision spectral data are of great significance to our understanding of star formation and evolution.

We report on observations of a solar jet propagating along coronal loops taken by the Solar Dynamics Observatory (SDO), the Interface Region Imaging Spectragraph (IRIS) and 1-m New Vacuum Solar Telescope (NVST). The ejecta of the jet consist of multi-thermal components and propagate with a speed greater than 100 km/s. Brightenings are found in the remote footpoints of the coronal loops having compact and round-shape in the Halpha images. The emission peak of the remote brightening in the Atmospheric Imaging Assembly (AIA) 94 \AA passband lags 60 s behind that in the jet base. The brightenings in the remote footpoints are believed to be consequences of heating by nonthermal electrons, MHD waves and/or conduction front generated by the magnetic reconnection processes of the jet. The heating in the remote footpoints leads to extension of the brightening along the loops toward the jet base, which is believed to be the chromospheric evaporation. This apparently acts as a brake on the ejecta, leading to a deceleration in the range from 1.5 to 3 km s$^{-2}$ with an error of $\sim1.0$\,km s$^{-2}$ when the chromospheric evaporation and the ejecta meet at locations near the loop apexes. The dynamics of this jet allows a unique opportunity to diagnose the chromospheric evaporation from the remote footpoints, from which we deduce a velocity in the range of 330--880 km/s.

We investigate the nonlinear growth of large-scale structure in the generalised Proca theory, in which a self-interacting massive vector field plays the role of driving the acceleration of the cosmic expansion. Focusing to the Proca Lagrangian at cubic order -- the cubic vector Galileon model -- we derive the simplified equations for gravity as well as the longitudinal and transverse modes of the vector field under the weak-field and quasi-static approximations, and implement them in a modified version of the ECOSMOG N-body code. Our simulations incorporate the Vainshtein screening effect, which reconciles the fifth force propagated by the longitudinal mode of the cubic vector Galileon model with local tests of gravity. The results confirm that for all scales probed by the simulation, the transverse mode has a negligible impact on structure formation in a realistic cosmological setup. It is well known that in this model the strength of the fifth force is controlled by a free model parameter, which we denote as beta_3. By running a suite of cosmological simulations for different values of beta_3, we show that this parameter also determines the effectiveness of the Vainshtein screening. The model behaves identically to the cubic scalar Galileon for beta_3 going to zero, in which the fifth force is strong in unscreened regions but is efficiently screened in high-density regions. In the opposite limit, beta_3 going to infinity, the model approaches its `quintessence' counterpart, which has a vanishing fifth force but a modified expansion history compared to LambdaCDM. This endows the model with rich phenomenology, which will be investigated in future works.

Using Stergioulas's RNS code for investigating fast pulsars with Equation of States (EOSs) on the causality surface (where the speed of sound equals to that of light) of the high-density EOS parameter space satisfying all known constraints from both nuclear physics and astrophysics, we show that the GW190814's secondary component of mass $(2.50-2.67)$ M$_{\odot}$ can be a super-fast pulsar spinning faster than 971 Hz about 42\% below its Kepler frequency. There is a large and physically allowed EOS parameter space below the causality surface where pulsars heavier than 2.50 M$_{\odot}$ are supported if they can rotate even faster with critical frequencies depending strongly on the high-density behavior of nuclear symmetry energy.

Using an explicitly isospin-dependent parametric Equation of State (EOS) for the core of neutron stars (NSs) within the Bayesian statistical approach, we infer the EOS parameters of super-dense neutron-rich nuclear matter from three sets of imagined mass-radius correlation data representing typical predictions by various nuclear many-body theories, i.e, the radius stays the same, decreases or increases with increasing NS mass within $\pm 15\%$ between 1.4 M$_{\odot}$ and 2.0 M$_{\odot}$. The corresponding average density increases quickly, slowly or slightly decreases as the NS mass increases from 1.4 M$_{\odot}$ to 2.0 M$_{\odot}$. Using the posterior probability distribution functions (PDFs) of EOS parameters inferred from GW170817 and NICER radius data for canonical NSs as references, we investigate how future radius measurements of massive NS will improve our knowledge about the EOS of super-dense neutron-rich nuclear matter, especially its symmetry energy term, compared to what people have already learned from analyzing the GW170817 and NICER data. While the EOS of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding high-density symmetry energies at densities above about $2\rho_0$ are very different, indicating that the radii of massive NSs carry reliable information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS.

We carry out ``full-physics'' hydrodynamical simulations of galaxy formation in the normal-branch Dvali-Gabadadze-Porrati (nDGP) braneworld model using a new modified version of the {\sc Arepo} code and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) which represent, respectively, weak and moderate departures from GR, in boxes of sizes $62\,h^{-1}{\rm Mpc}$ and $25\,h^{-1}{\rm Mpc}$ using $2\times512^3$ dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to $\sim 25$ per cent on large scales. Similar to their impact in the standard cosmological model ($\Lambda$CDM), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to $\sim15$ per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above $\sim0.1$ times the halo radius.

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

A fast radio burst was recently detected to be associated with a hard X-ray burst from the Galactic magnetar SGR 1935+2154. Scenarios involving magnetars for FRBs are hence highly favored. In this work, we suggest that the impact between an asteroid and a magnetar could explain such a detection. According to our calculations, an asteroid of mass $10^{20}$ g will be disrupted at $7 \times 10^9$ cm when approaching the magnetar, the accreted material will flow along the magnetic field lines at Alfv\'en radius $\sim 10^7$ cm. After falling onto the magnetar surface, an instant accretion column will be formed, producing a Comptonized X-ray burst and an FRB in the magnetosphere. We show that all the observational features could be interpreted self-consistently in our scenario. Quasi-periodic oscillations in this specific X-ray burst may help to verify our scenario.

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

Within the minimum model of neutron stars (NS) consisting of neutrons, protons and electrons, a new approach is proposed for inferring the symmetry energy of super-dense neutron-rich nucleonic matter above twice the saturation density $\rho_0$ of nuclear matter directly from integrating iteratively barotropic pressures in both neutron stars and heavy-ion reactions. Simultaneously, the proton fraction of NSs at $\beta$ equilibrium is extracted as a function of baryon density from the same procedure. An application of this approach using the NS pressure from GW170817 and the pressure in cold symmetric nuclear matter (SNM) extracted earlier by analyzing nuclear collective flow data in relativistic heavy-ion collisions provides a useful constraining band for the symmetry energy above $2\rho_0$.

A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects.

The structural property of the magnetic field in flare-bearing solar active regions (ARs) is one of the key aspects for understanding and forecasting solar flares. In this paper, we make a comparative analysis on the chirality and magnetic configurations associated with two X-class two-ribbon flares happening in AR 10930 and AR 11158. The photospheric magnetic fields of the two ARs were observed by space-based instruments, and the corresponding coronal magnetic fields were calculated based on the nonlinear force-free field model. The analysis shows that the electric current in the two ARs was distributed mostly around the main polarity inversion lines (PILs) where the flares happened, and the magnetic chirality (indicated by the signs of force-free factor $\alpha$) along the main PILs is opposite for the two ARs, i.e., left-handed ($\alpha<0$) for AR 10930 and right-handed ($\alpha>0$) for AR 11158. It is found that, for both the flare events, a prominent magnetic connectivity (featured by co-localized strong $\alpha$ and strong current density distributions) was formed along the main PIL before flare and was totally broken after flare eruption. The two branches of the broken magnetic connectivity, combined with the prominent magnetic connectivity before flare, compose the opposite magnetic configurations in the two ARs owing to their opposite chirality, i.e., Z-shaped configuration in AR 10930 with left-handed chirality and inverse Z-shaped configuration in AR 11158 with right-handed chirality. It is speculated that two-ribbon flares can be generally classified to these two magnetic configurations by chirality in the flare source regions of ARs.

We study nonlinear structure formation in the Bound Dark Energy model (BDE), where dark energy (DE) corresponds to a light scalar meson particle $\phi$ dynamically formed at a condensation energy scale $\Lambda_c$. The evolution of this dark-energy meson is determined by the potential $V(\phi)=\Lambda_c^{4+2/3}\phi^{-2/3}$, with a distinguishing phenomenology from other quintessence scenarios. Particularly, the expansion rate of the universe is affected not only at late times, but also when the condensation of $\phi$ occurs, which in linear theory leads to an enhancement (with respect to standard $\Lambda$CDM) of matter perturbations on small scales. We study how much of this signature is still present at late times as well as the properties of dark matter halos in the nonlinear regime through N-body simulations. Our results show that nonlinear corrections wash out this feature from the matter power spectrum even before DE becomes dominant. There is, however, a small but clear suppression of the BDE spectrum of $2\%$ today on the largest scales due to the distinct late-time dynamics of DE. The differences on the clustering power between BDE and $\Lambda$CDM are reflected in the halo mass function, where small halos are more abundant in BDE as opposed to large heavy structures, whose formation is delayed because of the expansion history of the universe. This result is well captured by the semi-analytical Sheth-Tormen formula. However, despite these differences, the halo concentration parameter is essentially the same in both models, which suggest that clustering inside the halos decouple from the general expansion once the halos form.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

Cosmic voids are a key component of the large-scale structure that contain a plethora of cosmological information. Typically, voids are identified from the underlying galaxy distribution, which is a biased tracer of the total matter field. Previous works have shown that 2D voids identified in weak lensing maps -- weak lensing voids -- correspond better to true underdense regions along the line of sight. In this work, we study how the properties of weak lensing voids depend on the choice of void finder, by adapting several popular void finders. We present and discuss the differences between identifying voids directly in the convergence maps, and in the distribution of weak lensing peaks. Particular effort has been made to test how these results are affected by galaxy shape noise, which is a dominant source of noise in weak lensing observations. By studying the signal-to-noise ratios (SNR) for the tangential shear profile of each void finder, we find that voids identified directly in the convergence maps have the highest SNR but are also the ones most affected by galaxy shape noise. Troughs are least affected by noise, but also have the lowest SNR. The tunnel algorithm, which identifies voids in the distribution of weak lensing peaks, represents a good compromise between finding a large tangential shear SNR and mitigating the effect of galaxy shape noise.

We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a minimal distortion gauge, where the linear growth rate is scale-dependent, and the standard method for realising initial particle data is not straightforwardly applicable. A new method is introduced, in which the initial positions of particles are generated from the displacement field realised for a matter power spectrum as usual, but the velocity is calculated by finite-differencing the displacement fields around the initial redshift. In this way, all the information required for setting up the initial conditions is drawn from three consecutive input matter power spectra, and additional assumptions such as scale-independence of the linear growth factor and growth rate are not needed. We implement this method in a modified 2LPTic code, and demonstrate that in a Newtonian setting it can reproduce the velocity field given by the default 2LPTic code with subpercent accuracy. We also show that the matter and velocity power spectra of the initial particle data generated for GRAMSES simulations using this method agree very well with the linear-theory predictions in the particular gauge used by GRAMSES. Finally, we discuss corrections to the finite difference calculation of the velocity when radiation is present, as well as additional corrections implemented in GRAMSES to ensure consistency. This method can be applied in ICs generation for GR simulations in generic gauges, and simulations of cosmological models with scale-dependent linear growth rate.

The experiment to detect the global epoch of reionization signature (EDGES) collaboration reported the detection of a line at 78MHz in the sky-averaged spectrum due to neutral hydrogen (HI) 21-cm hyperfine absorption of cosmic microwave background photons at $z\sim 17$. This requires that the spin temperature of HI be coupled to the kinetic temperature of the gas at this $z$ through the scattering of Ly-$\alpha$ photons emitted by massive stars. To explain the experimental result, star formation needs to be sufficiently efficient at $z\sim 17$ and this can be used to constrain models in which small-scale structure formation is suppressed (DMF models), either due to dark matter free-streaming or non-standard inflationary dynamics. We combine simulations of structure formation with a simple recipe for star formation to investigate whether these models emit enough Ly-$\alpha$ photons to reproduce the experimental signal for reasonable values of the star formation efficiency, $f_\star$. We find that a thermal warm dark matter (WDM) model with mass $m_\mathrm{WDM}\sim 4.3\,\mathrm{keV}$ is consistent with the timing of the signal for $f_\star\lesssim 2\%$. The exponential growth of structure around $z\sim 17$ in such a model naturally generates a sharp onset of the absorption. A warmer model with $m_\mathrm{WDM}\sim 3\,\mathrm{keV}$ requires a higher star formation efficiency, $f_\star\sim 6\%$, which is a factor of few above predictions of current star formation models and observations of satellites in the Milky Way. However, uncertainties in the process of star formation at these $z$ do not allow to derive strong constrains on such models using 21-cm absorption line. The onset of the 21-cm absorption is generally faster in DMF compared to cold dark matter (CDM), unless some process significantly suppresses star formation in halos with masses below $\sim 10^8\,h^{-1}\,\mathrm{M}_\odot$.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

We investigate the consequences of different regularizations and ambiguities in loop cosmological models on the predictions in the scalar and tensor primordial spectrum of the cosmic microwave background using the dressed metric approach. Three models, standard loop quantum cosmology (LQC), and two modified loop quantum cosmologies (mLQC-I and mLQC-II) arising from different regularizations of the Lorentzian term in the classical Hamiltonian constraint are explored for chaotic inflation in spatially-flat Friedmann-Lema\^itre-Robertson-Walker (FLRW) universe. In each model, two different treatments of the conjugate momentum of the scale factor are considered. The first one corresponds to the conventional treatment in dressed metric approach, and the second one is inspired from the hybrid approach which uses the effective Hamiltonian constraint. For these two choices, we find the power spectrum to be scale-invariant in the ultra-violet regime for all three models, but there is at least a $10\%$ relative difference in amplitude in the infra-red and intermediate regimes. All three models result in significant differences in the latter regimes. In mLQC-I, the magnitude of the power spectrum in the infra-red regime is of the order of Planck scale irrespective of the ambiguity in conjugate momentum of the scale factor. The relative difference in the amplitude of the power spectrum between LQC and mLQC-II can be as large as $50\%$ throughout the infra-red and intermediate regimes. Differences in amplitude due to regularizations and ambiguities turn out to be small in the ultra-violet regime.

Solar filaments are dark structures on the solar disk, with an elongated spine and several barbs extending out from the spine. When appearing above the solar limb, a filament is called a prominence, with several feet extending down to the solar surface. It was generally thought that filament barbs are simply the prominence feet veering away from the spine and down to the solar surface. However, it was recently noticed that there might be another dynamic type of barbs, which were proposed to be due to filament thread longitudinal oscillation. If this is the case, the dynamic barbs would not extend down to the solar surface. With the quadrature observations of a filament barb on 2011 June 5 from the {\it Solar Dynamics Observatory} and the {\it STEREO} satellites, we confirm that the filament barb is due to filament thread longitudinal oscillations. Viewed from the side, the filament barb looks like an appendix along the spine of the prominence, and does not extend down to the solar surface as a foot.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

We report pulsar timing observations carried out in L-band with NTSC's 40-meter Haoping Radio Telescope (HRT), which was constructed in 2014. The observations were carried out using the pulsar machine we developed. Timing observations toward millisecond pulsar J0437-4715 obtains a timing residual (r.m.s) of 397ns in the time span of 284 days. And our observations successfully detected Crab pulsar's glitch that happened on July 23rd, 2019.

Similar to the Schwarzschild coordinates for spherical black holes, the Baldwin, Jeffery and Rosen (BJR) coordinates for plane gravitational waves are often singular, and extensions beyond such singularities are necessary, before studying asymptotic properties of such spacetimes at the null infinity of the plane, on which the gravitational waves propagate. The latter is closely related to the studies of memory effects and soft graviton theorems. In this paper, we point out that in the BJR coordinates all the spacetimes are singular physically at the focused point $u = u_s$, except for the two cases: (1) $\alpha =1/2, \; \forall \; \chi_n$; and (2) $\alpha =1, \; \chi_i = 0\; (i = 1, 2, 3)$, where $\chi_n$ are the coefficients in the expansion $\chi \equiv \left[{\mbox{det}}\left(g_{ab}\right) \right]^{1/4} = \left(u - u_s\right)^{\alpha}\sum_{n = 0}^{\infty}\chi_n \left(u - u_s\right)^n$ with $\chi_0 \not= 0$, the constant $\alpha \in (0, 1]$ characterizes the strength of the singularities, and $g_{ab}$ denotes the reduced metric on the two-dimensional plane orthogonal to the propagation direction of the wave. Therefore, the hypersurfaces $u= u_s$ already represent the boundaries of such spacetimes, and the null infinity does not belong to them. As a result, they cannot be used to study properties of plane gravitational waves at null infinities, including memory effects and soft graviton theorems.

Muons in neutron stars (NSs) play especially important roles in addressing several interesting new physics questions associated with detecting as well as understanding interactions and astrophysical effects of muonphilic dark matter particles. The key model inputs for studying the latter are the total muon mass $M_{\mu}$, the muon mass fraction $M_{\mu}/M_{\rm NS}$ over the NS mass $M_{\rm NS}$ and the muon radial density profile $\rho_{\mu}(r)$ in NSs of varying masses. We investigate these quantities within a minimum model for the core of NSs consisting of neutrons, protons, electrons, and muons using an explicitly isospin-dependent parametric Equation of State (EOS) constrained by available nuclear laboratory experiments and the latest astrophysical observations of NS masses, radii and tidal deformabilities. We found that the absolutely maximum muon mass $M_{\mu}$ and its mass fraction $M_{\mu}/M_{\rm NS}$ in the most massive NSs allowed by causality are about 0.025 $M_\odot$ and 1.1\%, respectively. For the most massive NS of mass 2.14 $M_\odot$ observed so far, they reduce to about 0.020 $M_\odot$ and 0.9\%, respectively. We also study respective effects of individual parameters describing the EOS of high-density neutron-rich nucleonic matter on the muon contents in NSs with varying masses. We found that the most important but uncertain nuclear physics ingredient for determining the muon contents in NSs is the high-density nuclear symmetry energy.

Here we propose a new scheme to reconstruct the baryon acoustic oscillations (BAO) signal, with key cosmological information, based on deep convolutional neural networks. After training the network with almost no fine-tuning, in the test set, the network recovers large-scale modes accurately: the correlation coefficient between the ground truth and recovered initial conditions still reach $90\%$ at $k \leq 0.2~ h\mathrm{Mpc}^{-1}$, which significantly improves the BAO signal-to-noise ratio until the scale $k=0.4~ h\mathrm{Mpc}^{-1}$. Furthermore, our scheme is independent of the survey boundary since it reconstructs initial condition based on local density distribution in configuration space, which means that we can gain more information from the whole survey space. Finally, we found our trained network is not sensitive to the cosmological parameters and works very well in those cosmologies close to that of our training set. This new scheme will possibly help us dig out more information from the current, on-going and future galaxy surveys.

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.

Winds from black hole accretion flows are ubiquitous. Previous works mainly focus on the launching of wind in the accretion flow scale. It still remains unclear how far the winds can propagate outward and what is their large-scale dynamics. As the first paper of this series, we study the large-scale dynamics of thermal wind beyond accretion scales via analytical and numerical methods. Boundary conditions, which are crucial to our problem, are analyzed and presented based on the small-scale simulations combined with observations of winds. Both black hole and galaxy potential are taken into account. For winds originated from hot accretion flows, we find that the wind can reach to large scales. The radial profiles of velocity, density, and temperature can be approximated by $v_r\approx v_{r0}, \rho\approx \rho_{0}(r/r_0)^{-2}$, and $T\approx T_0 (r/r_0)^{-2(\gamma-1)}$, where $v_{r0}, \rho_0, T_0$ are the velocity, density, and temperature of winds at the boundary $r_0(\equiv 10^3 r_g)$, $\gamma$ is the polytropic index. During the outward propagation, the enthalpy and the rotational energy compensate the increase of gravitational potential. For thin disks, we find that because the Bernoulli parameter is smaller, winds cannot propagate as far as the hot winds, but stop at a certain radius where the Bernoulli parameter is equal to the potential energy. Before the winds stop, the profiles of dynamical quantities can also be approximated by the above relations. In this case the rotational energy alone compensates the increase of the potential energy.

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments.

We present GRAMSES, a new pipeline for nonlinear cosmological $N$-body simulations in General Relativity (GR). This code adopts the Arnowitt-Deser-Misner (ADM) formalism of GR, with constant mean curvature and minimum distortion gauge fixings, which provides a fully nonlinear and background independent framework for relativistic cosmology. Employing a fully constrained formulation, the Einstein equations are reduced to a set of ten elliptical equations which are solved using multigrid relaxation with adaptive mesh refinements (AMR), and three hyperbolic equations for the evolution of tensor degrees of freedom. The current version of GRAMSES neglects the latter by using the conformal flatness approximation, which allows it to compute the two scalar and two vector degrees of freedom of the metric. In this paper we describe the methodology, implementation, code tests and first results for cosmological simulations in a $\Lambda$CDM universe, while the generation of initial conditions and physical results will be discussed elsewhere. Inheriting the efficient AMR and massive parallelisation infrastructure from the publicly-available $N$-body and hydrodynamic simulation code RAMSES, GRAMSES is ideal for studying the detailed behaviour of spacetime inside virialised cosmic structures and hence accurately quantifying the impact of backreaction effects on the cosmic expansion, as well as for investigating GR effects on cosmological observables using cosmic-volume simulations.

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.

Within the Bayesian framework using available constraining bands on the pressure in symmetric nuclear matter (SNM) derived earlier by others in the density range of 1.3$\rho_0$ to 4.5$\rho_0$ from kaon production and nuclear collective flow data in energetic heavy-ion collisions, we infer the posterior probability distribution functions (PDFs) of SNM incompressibility $K_0$ and skewness $J_0$ using uniform prior PDFs for them in the ranges of $220\leq K_0\leq 260$ MeV and $-800\leq J_0\leq 400$ MeV. The 68\% posterior credible boundaries around the most probable values of $K_0$ and $J_0$ are found to be 222$\pm$2 MeV and -215$\pm$20 MeV, respectively, much narrower than their prior ranges widely used currently in the literature and are consistent with the results of a recent Bayesian analysis of neutron star properties constrained by available X-ray and gravitational wave observations.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

We investigate if, for a fixed number density of targets and redshift, there is an optimal way to select a galaxy sample in order to measure the baryon acoustic oscillation (BAO) scale, which is used as a standard ruler to constrain the cosmic expansion. Using the mock galaxy catalogue built by Smith et al. in the Millennium-XXL N-body simulation with a technique to assign galaxies to dark matter haloes based on halo occupation distribution modelling, we consider the clustering of galaxies selected by luminosity, colour and local density. We assess how well the BAO scale can be extracted by fitting a template to the power spectrum measured for each sample. We find that the BAO peak position is recovered equally well for samples defined by luminosity or colour, while there is a bias in the BAO scale recovered for samples defined by density. The BAO position is contracted to smaller scales for the densest galaxy quartile and expanded to large scales for the two least dense galaxy quartiles. For fixed galaxy number density, density-selected samples have higher uncertainties in the recovered BAO scale than luminosity- or colour-selected samples.

We develop perturbation theory approaches to model the marked correlation function constructed to up-weight low density regions of the Universe, which might help distinguish modified gravity models from the standard cosmology model based on general relativity. Working within Convolution Lagrangian Perturbation Theory, we obtain that weighted correlation functions are expressible as double convolution integrals that we approximate using a combination of Eulerian and Lagrangian schemes. We find that different approaches agree within 1$\%$ on quasi non-linear scales. Compared with {\it N}-body simulations, the perturbation theory is found to provide accurate predictions for the marked correlation function of dark matter fields, dark matter halos as well as Halo Occupation Distribution galaxies down to $30$ Mpc/h. These analytic approaches help to understand the degeneracy between the mark and the galaxy bias and find a way to maximize the differences among various cosmological models.

Loops are one of the fundamental structures that trace the geometry of the magnetic field in the solar atmosphere. Their evolution and dynamics provide a crucial proxy for studying how the magnetized structures are formed and heated in the solar atmosphere. Here, we report on spectroscopic observations of a set of transition region loops taken by the Interface Region Imaging Spectrograph (IRIS) at Si IV 1394 \AA\ with a sit-and-stare mode. The loops are corresponding to the flux emergence at its very late phase when the emerged magentic features in the photosphere have fully developed. We find the transition region loops are still expanding and moving upward with a velocity of a few kilometers per second ($\lesssim$10 km/s) at this stage. The expansion of the loops leads to interactions between themselves and the ambient field, which can drive magnetic reconnection evidenced by multiple intense brightenings, including transition region explosive events and IRIS bombs in the footpoint region associated with the moving polarity. A set of quasi-periodic brightenings with a period of about 130 s is found at the loop apex, from which the Si IV 1394 \AA\ profiles are significantly non-Gaussian with enhancements at both blue and red wings at Doppler velocities of about 50 km/s. We suggest that the transition region loops in the very late phase of flux emergence can be powered by heating events generated by the interactions between the expanding loops and the ambient fields and also by (quasi-)periodic processes, such as oscillation-modulated braiding reconnection.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Modifications to General Relativity (GR) often incorporate screening mechanisms in order to remain compatible with existing tests of gravity. The screening is less efficient in underdense regions, which suggests that cosmic voids can be a useful cosmological probe for constraining modified gravity models. In particular, weak lensing by voids has been proposed as a promising test of such theories. Usually, voids are identified from galaxy distributions, making them biased tracers of the underlying matter field. An alternative approach is to study voids identified in weak lensing maps -- weak lensing voids -- which have been shown to better correspond to true underdense regions. In this paper, we study the ability of weak lensing voids to detect the signatures of modified gravity. Focusing on the void abundance and weak lensing profiles, we find that both statistics are sensitive probes of gravity. These are quantified in terms of the signal-to-noise ratios (SNR) with which an LSST-like survey will be able to distinguish between different gravity models. We find that the tangential shear profiles of weak lensing voids are considerably better than galaxy voids at this, though voids have somewhat lower SNR than weak lensing peaks. The abundances of voids and peaks have respectively $\rm{SNR} = 50$ and $70$ for a popular class of modified gravity in an LSST-like survey.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

Solar energetic particles (SEPs) are an important product of solar activity. They are connected to solar active regions and flares, coronal mass ejections (CMEs), EUV waves, shocks, Type II and III radio emissions, and X-ray bursts. These phenomena are major probes of the partition of energy in solar eruptions, as well as for the organization, dynamics, and relaxation of coronal and interplanetary magnetic fields. Many of these phenomena cause terrestrial space weather, posing multiple hazards for humans and their technology from space to the ground. Since particular flares, shocks, CMEs, and EUV waves produce SEP events but others do not, since propagation effects from the low corona to 1 AU appear important for some events but not others, and since Type II and III radio emissions and X-ray bursts are sometimes produced by energetic particles leaving these acceleration sites, it is necessary to study the whole system with a multi-frequency and multi-instrument perspective that combines both in-situ and remote observations with detailed modelling of phenomena. This article demonstrates this comprehensive approach, and shows its necessity, by analysing a trio of unusual and striking solar eruptions, radio and X-ray bursts, and SEP events that occurred on 4 November 2015. These events show both strong similarities and differences from standard events and each other, despite having very similar interplanetary conditions and only two are sites and CME genesis regions. They are therefore major targets for further in-depth observational studies, and for testing both existing and new theories and models. Based on the very limited modelling available we identify the aspects that are and are not understood, and we discuss ideas that may lead to improved understanding of the SEP, radio, and space-weather events.

We study the statistics of weak lensing convergence peaks, such as their abundance and two-point correlation function (2PCF), for a wide range of cosmological parameters $\Omega_m$ and $\sigma_8$ within the standard $\Lambda$CDM paradigm, focusing on intermediate-height peaks with signal-to-noise ratio (SNR) of $1.5$ to $3.5$. We find that the cosmology dependence of the peak abundance can be described by a one-parameter fitting formula that is accurate to within $\sim3\%$. The peak 2PCFs are shown to feature a self-similar behaviour: if the peak separation is rescaled by the mean inter-peak distance, catalogues with different minimum peak SNR values have identical clustering, which suggests that the peak abundance and clustering are closely interconnected. A simple fitting model for the rescaled 2PCF is given, which together with the peak abundance model above can predict peak 2PCFs with an accuracy better than $\sim5\%$. The abundance and 2PCFs for intermediate peaks have very different dependencies on $\Omega_m$ and $\sigma_8$, implying that their combination can be used to break the degeneracy between these two parameters.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

Recently, GRB 170817A was confirmed to be associated with GW 170817, which was produced by a neutron star - neutron star merger. It indicates that at least some short gamma-ray bursts come from binary neutron star mergers. Theoretically, it is widely accepted that short gamma-ray bursts can be produced by two distinctly different mechanisms, binary neutron star mergers and neutron star - black hole mergers. These two kinds of bursts should be different observationally due to their different trigger mechanisms. Motivated by this idea, we collect a universal data set constituted of 51 short gamma-ray bursts observed by $Swift$/BAT, among which 14 events have extended emission component. We study the observational features of these 51 events statistically. It is found that our samples are consisted of two distinct groups. They clearly show a bimodal distribution when their peak photon fluxes at 15-150 keV band are plotted against the corresponding fluences. Most interestingly, all the 14 short bursts with extended emission lie in a particular region. When the fluences are plotted against the burst durations, short bursts with extended emission again tend to concentrate in the long duration segment. These features strongly indicate that short gamma-ray bursts really may come from two distinct types of progenitors. We argue that those short gamma-ray bursts with extended emission come from the coalescence of neutron stars, while the short gamma-ray bursts without extended emission come from neutron star - black hole mergers.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

We present an analysis of the matter, halo and galaxy clustering in $f(R)$-gravity employing the SHYBONE full-physics hydrodynamical simulation suite. Our analysis focuses on the interplay between baryonic feedback and $f(R)$-gravity in the matter power spectrum, the matter and halo correlation functions, the halo and galaxy-host-halo mass function, the subhalo and satellite-galaxy count and the correlation function of the stars in our simulations. Our studies of the matter power spectrum in full physics simulations in $f(R)$-gravity show, that it will be very difficult to derive accurate fitting formulae for the power spectrum enhancement in $f(R)$-gravity which include baryonic effects. We find that the enhancement of the halo mass function due to $f(R)$-gravity and its suppression due to feedback effects do not show significant back-reaction effects and can thus be estimated from independent GR-hydro and $f(R)$ dark matter only simulations. Our simulations furthermore show, that the number of subhalos and satellite-galaxies per halo is not significantly affected by $f(R)$-gravity. Low mass halos are nevertheless more likely to be populated by galaxies in $f(R)$-gravity. This suppresses the clustering of stars and the galaxy correlation function in the theory compared to standard cosmology.

Future 21 cm intensity mapping surveys such as SKA can provide precise information on the spatial distribution of the neutral hydrogen (HI) in the post-reionization epoch. This information will allow to test the standard LCDM paradigm and with that the nature of gravity. In this work, we employ the SHYBONE simulations, which model galaxy formation in f(R) modified gravity using the IllustrisTNG model, to study the effects of modified gravity on HI abundance and power spectra. We find that the enhanced growth low-mass dark matter halos experience in f(R) gravity at high redshifts alters the HI power spectrum and can be observable through 21 cm intensity mapping. Our results suggest that the HI power spectrum is suppressed by 13\% on scales $k\lesssim2\,h\,\mathrm{Mpc}^{-1}$ at z=2 for F6, a f(R) model which passes most observational constraints. We show that this suppression can be detectable by SKA1-MID with 1000 hours of exposure time, making HI clustering a novel test of gravity at high redshift.

We introduce the LOCal Universe Screening Test Suite (LOCUSTS) project, an effort to create 'screening maps' in the nearby Universe to identify regions in our neighbourhood which are screened, i.e., regions where deviations from General Relativity (GR) are suppressed, in various modified gravity (MG) models. In these models, deviations from the GR force law are often stronger for smaller astrophysical objects, making them ideal test beds of gravity in the local Universe. However, the actual behaviour of the modified gravity force also depends on the environment of the objects, and to make accurate predictions one has to take the latter into account. This can be done approximately using luminous objects in the local Universe as tracers of the underlying dark matter field. Here, we propose a new approach that takes advantage of state-of-the-art Bayesian reconstruction of the mass distribution in the Universe, which allows us to solve the modified gravity equations and predict the screening effect more accurately. This is the first of a series of works, in which we present our methodology and some qualitative results of screening for a specific MG model, $f(R)$ gravity. Applications to test models using observations and extensions to other classes of models will be studied in future works.

It is well-known that ordinary supernovae (SNe) are powered by 56Ni cascade decay. Broad-lined type Ic SNe (SNe Ic-BL) are a subclass of SNe that are not all exclusively powered by 56Ni decay. It was suggested that some SNe Ic-BL are powered by magnetar spin-down. iPTF16asu is a peculiar broad-lined type Ic supernova discovered by the intermediate Palomar Transient Factory. With a rest-frame rise time of only 4 days, iPTF16asu challenges the existing popular models, for example, the radioactive heating (56Ni-only) and the magnetar+56Ni models. Here we show that this rapid rise could be attributed to interaction between the SN ejecta and a pre-existing circumstellar medium ejected by the progenitor during its final stages of evolution, while the late-time light curve can be better explained by energy input from a rapidly spinning magnetar. This model is a natural extension to the previous magnetar model. The mass-loss rate of the progenitor and ejecta mass are consistent with a progenitor that experienced a common envelope evolution in a binary. An alternative model for the early rapid rise of the light curve is the cooling of a shock propagating into an extended envelope of the progenitor. It is difficult at this stage to tell which model (interaction+magnetar+56Ni or cooling+magnetar+56Ni) is better for iPTF16asu. However, it is worth noting that the inferred envelope mass in the cooling+magnetar+56Ni is very high.

The radius $R_{1.4}$ of neutron stars (NSs) with a mass of 1.4 M$_{\odot}$ has been extracted consistently in many recent studies in the literature. Using representative $R_{1.4}$ data, we infer high-density nuclear symmetry energy $E_{\rm{sym}}(\rho)$ and the associated nucleon specific energy $E_0(\rho)$ in symmetric nuclear matter (SNM) within a Bayesian statistical approach using an explicitly isospin-dependent parametric Equation of State (EOS) for nucleonic matter. We found that: (1) The available astrophysical data can already improve significantly our current knowledge about the EOS in the density range of $\rho_0-2.5\rho_0$. In particular, the symmetry energy at twice the saturation density $\rho_0$ of nuclear matter is determined to be $E_{\mathrm{sym}}(2\rho_0)$ =39.2$_{-8.2}^{+12.1}$ MeV at 68\% confidence level. (2) A precise measurement of the $R_{1.4}$ alone with a 4\% 1$\sigma$ statistical error but no systematic error will not improve much the constraints on the EOS of dense neutron-rich nucleonic matter compared to what we extracted from using the available radius data. (3) The $R_{1.4}$ radius data and other general conditions, such as the observed NS maximum mass and causality condition introduce strong correlations for the high-order EOS parameters. Consequently, the high-density behavior of $E_{\rm{sym}}(\rho)$ inferred depends strongly on how the high-density SNM EOS $E_0(\rho)$ is parameterized, and vice versa. (4) The value of the observed maximum NS mass and whether it is used as a sharp cut-off for the minimum maximum mass or through a Gaussian distribution affect significantly the lower boundaries of both the $E_0(\rho)$ and $E_{\rm{sym}}(\rho)$ only at densities higher than about $2.5\rho_0$.

We study the evolution of spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker universe for chaotic and Starobinsky potentials in the framework of modified loop quantum cosmologies. These models result in a non-singular bounce as in loop quantum cosmology, but with far more complex modified Friedmann dynamics with higher order than quadratic terms in energy density. For the kinetic energy dominated bounce, we obtain analytical solutions using different approximations and compare with numerical evolution for various physical variables. The relative error turns out to be less than $0.3\%$ in the bounce regime for both of the potentials. Generic features of dynamics, shared with loop quantum cosmology, are established using analytical and numerical solutions. Detailed properties of three distinct phases in dynamics separating bounce regime, transition stage and inflationary phase are studied. For the potential energy dominated bounce, we qualitatively describe its generic features and confirm by simulations that they all lead to the desired slow-roll phase in the chaotic inflation. However, in the Starobinsky potential, the potential energy dominated bounce cannot give rise to any inflationary phase. Finally, we compute the probability for the desired slow-roll inflation to occur in the chaotic inflation and as in loop quantum cosmology, find a very large probability for the universe to undergo inflation.

We study the extreme-ultraviolet (EUV) emissions modulated by leaky fast sausage modes (FSMs) in solar active region loops and examine their observational signatures via spectrometers like EIS. After computing fluid variables of leaky FSMs with MHD simulations, we forward-model the intensity and spectral properties of the Fe X 185~\AA~and Fe XII 195~\AA~lines by incorporating non-equilibrium ionization (NEI) in the computations of the relevant ionic fractions. The damping times derived from the intensity variations are then compared with the wave values, namely the damping times directly found from our MHD simulations. Our results show that in the equilibrium ionization cases, the density variations and the intensity variations can be either in phase or in anti-phase, depending on the loop temperature. NEI considerably impacts the intensity variations but has only marginal effects on the derived Doppler velocity or Doppler width. We find that the damping time derived from the intensity can largely reflect the wave damping time if the loop temperature is not drastically different from the nominal formation temperature of the corresponding emission line. These results are helpful for understanding the modulations to the EUV emissions by leaky FSMs and hence helpful for identifying FSMs in solar active region loops.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

We investigate cosmological models in which dynamical dark energy consists of a scalar field whose present-day value is controlled by a coupling to the neutrino sector. The behaviour of the scalar field depends on three functions: a kinetic function, the scalar field potential, and the scalar field-neutrino coupling function. We present an analytic treatment of the background evolution during radiation- and matter-domination for exponential and inverse power law potentials, and find a relaxation of constraints compared to previous work on the amount of early dark energy in the exponential case. We then carry out a numerical analysis of the background cosmology for both types of potential and various illustrative choices of the kinetic and coupling functions. By applying bounds from Planck on the amount of early dark energy, we are able to constrain the magnitude of the kinetic function at early times.

It has been found that the Kelvin-Helmholtz instability (KHI) induced by both transverse and torsional oscillations in coronal loops can reinforce the effects of wave heating. In this study, we model a coronal loop as a system of individual strands, and we study wave heating effects by considering a combined transverse and torsional driver at the loop footpoint. We deposit the same energy into the multi-stranded loop and an equivalent monolithic loop, and then observe a faster increase in the internal energy and temperature in the multi-stranded model. Therefore, the multi-stranded model is more efficient in starting the heating process. Moreover, higher temperature is observed near the footpoint in the multi-stranded loop and near the apex in the monolithic loop. The apparent heating location in the multi-stranded loop agrees with the previous predictions and observations. Given the differences in the results from our multi-stranded loop and monolithic loop simulations, and given that coronal loops are suggested to be multi-stranded on both theoretical and observational grounds, our results suggest that the multi-strandedness of coronal loops needs to be incorporated in future wave-based heating mechanisms.

We present a general method to compute the nonlinear matter power spectrum for dark energy and modified gravity scenarios with percent-level accuracy. By adopting the halo model and nonlinear perturbation theory, we predict the reaction of a $\Lambda$CDM matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to $N$-body simulations we demonstrate that with no-free parameters we can recover the nonlinear matter power spectrum for a wide range of different $w_0$-$w_a$ dark energy models to better than 1% accuracy out to $k \approx 1 \, h \, {\rm Mpc}^{-1}$. We obtain a similar performance for both DGP and $f(R)$ gravity, with the nonlinear matter power spectrum predicted to better than 3% accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1%. With a single suite of standard $\Lambda$CDM $N$-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise nonlinear regime.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

We have carried out a set of cosmological hydrodynamical simulations that follow galaxy formation in $f(R)$ modified gravity models. Our simulations employ the Illustris-TNG full physics model and a new modified gravity solver in the AREPO code. For the first time we are able to investigate the degeneracy in the matter power spectrum between the effects of $f(R)$-gravity and feedback from active galactic nuclei (AGN), and the imprint of modified gravity on the properties of galaxies and on the distribution of dark matter, gas and stars in the universe. $f(R)$-gravity has an observable effect on the neutral hydrogen power spectrum at high redshift at a level of 20%. For both the F6 and F5 models, this is significantly larger than the predicted errors for the SKA1-MID survey, making this probe a powerful test of gravity on large scales. A similar effect is present in the power spectrum of the stars at high redshift. We also show that rotationally supported disc galaxies can form in $f(R)$-gravity, even in the partially screened regime. Our simulations indicate that there might be more disc galaxies in F6 compared to GR, and fewer in F5. Finally, we show that the back reaction between AGN feedback and modified gravity in the matter power spectrum is not important in the F6 model but has a sizeable effect in F5.

We investigate the dependence of the galaxy-halo connection and galaxy density field in modified gravity models using the $N-$body simulations for $f(R)$ and nDGP models at $z=0$. Because of the screening mechanisms employed by these models, chameleon and Vainshtein, halos are clustered differently in the non-linear regime of structure formation. We quantify their deviations in the galaxy density field from the standard $\Lambda$CDM model under different environments. We populate galaxies in halos via the (Sub)Halo Abundance Matching. Our main results are: 1) The galaxy-halo connection {\it strongly} depends on the gravity model; a maximum variation of $\sim40\%$ is observed between Halo Occupational Distribution (HOD) parameters; 2) $f(R)$ gravity models predict an excess of galaxies in low density environments of $\sim10\%$ but predict a deficit of $\sim10\%$ at high density environments for $|f_{R0}| = 10^{-4}$ and $10^{-6}$ while $|f_{R0}| = 10^{-5}$ predicts more high density structures; nDGP models are consistent with $\Lambda$CDM; 3) Different gravity models predict different dependences of the galaxy luminosity function (GLF) with the environment, especially in void-like regions we find differences around $\sim10\%$ for the $f(R)$ models while nDPG models remain closer to $\Lambda$CDM for low-luminosity galaxies but there is a deficit of $\sim11\%$ for high-luminosity galaxies in all environments. We conclude that the dependence of the GLF with environment might provide a test to distinguish between gravity models and their screening mechanisms from the $\Lambda$CDM. We provide HOD parameters for the gravity models analyzed in this paper.

Both coronal plumes and network jets are rooted in network lanes. The relationship between the two, however, has yet to be addressed. For this purpose, we perform an observational analysis using images acquired with the Atmospheric Imaging Assembly (AIA) 171{\AA} passband to follow the evolution of coronal plumes, the observations taken by the Interface Region Imaging Spectrograph (IRIS) slit-jaw 1330{\AA} to study the network jets, and the line-of-sight magnetograms taken by the Helioseismic and Magnetic Imager (HMI) to overview the the photospheric magnetic features in the regions. Four regions in the network lanes are identified, and labeled ``R1--R4''. We find that coronal plumes are clearly seen only in ``R1''&''R2'' but not in ``R3''&``R4'', even though network jets abound in all these regions. Furthermore, while magnetic features in all these regions are dominated by positive polarity, they are more compact (suggesting stronger convergence) in ``R1''&``R2'' than that in ``R3''&``R4''. We develop an automated method to identify and track the network jets in the regions. We find that the network jets rooted in ``R1''&``R2'' are higher and faster than that in ``R3''&``R4'',indicating that network regions producing stronger coronal plumes also tend to produce more dynamic network jets. We suggest that the stronger convergence in ``R1''&``R2'' might provide a condition for faster shocks and/or more small-scale magnetic reconnection events that power more dynamic network jets and coronal plumes.

Hemispherical power asymmetry has emerged as a new challenge to cosmology in early universe. While the cosmic microwave background (CMB) measurements indicated the asymmetry amplitude $A \simeq 0.07$ at the CMB scale $k_{\rm CMB}\simeq 0.0045\,{\rm Mpc}^{-1}$, the high-redshift quasar observations found no significant deviation from statistical isotropy. This conflict can be reconciled in some scale-dependent asymmetry models. We put forward a new parameterization of scale-dependent asymmetric power spectrum, inspired by a multi-speed inflation model. The 21-cm power spectrum from the epoch of reionization can be used to constrain the scale-dependent hemispherical asymmetry. We demonstrate that an optimum, multi-frequency observation by the Square Kilometre Array (SKA) Phase 2 can impose a constraint on the amplitude of the power asymmetry anomaly at the level of $\Delta A \simeq 0.2$ at $0.056 \lesssim k_{\rm 21cm} \lesssim 0.15 \,{\rm Mpc}^{-1}$. This limit may be further improved by an order of magnitude as $\Delta A \simeq 0.01$ with a cosmic variance limited experiment such as the Omniscope.

The broadband afterglow of GW170817/GRB 170817A is believed to be from an off-axis structured jet. The central engine of a gamma-ray burst usually launches a pair of outflows that move oppositely, it is reasonable to consider the emission from a double-sided structured jets, with a near-jet moving toward us and a counter-jet moving away from us. Assuming that the two branches of the jet have the same physical parameters, we have calculated their radio emission. It is found that the counter-jet component will emerge in the radio light curves $\sim$ 2500 days post-merger.It typically leads to a plateau in the light curve, thus could be hopefully revealed by high accuracy radio observations. We have also considered the possibilities that both branches have different parameters, and found that if some of the physical parameters of the counter-jet are evaluated more favorably than those of the near-jet, then the counter-jet emission will be enhanced and will clearly show up as an obvious plateau or even a rebrightening. For example, when the circum-burst medium encountered by the counter-jet is assumed to be 100 times denser than that of the near-jet, a remarkable radio plateau will appear at $\sim$ 600 days. However, in X-ray bands, the counter-jet component is generally too faint to be discerned. It is argued that the late radio observations of GW170817/GRB 170817A can help determine the key parameters and diagnose the environment of the event.

Determining the Equation of State (EOS) of dense neutron-rich nuclear matter is a shared goal of both nuclear physics and astrophysics. Except possible phase transitions, the density dependence of nuclear symmetry \esym is the most uncertain part of the EOS of neutron-rich nucleonic matter especially at supra-saturation densities. Much progresses have been made in recent years in predicting the symmetry energy and understanding why it is still very uncertain using various microscopic nuclear many-body theories and phenomenological models. Simultaneously, significant progresses have also been made in probing the symmetry energy in both terrestrial nuclear laboratories and astrophysical observatories. In light of the GW170817 event as well as ongoing or planned nuclear experiments and astrophysical observations probing the EOS of dense neutron-rich matter, we review recent progresses and identify new challenges to the best knowledge we have on several selected topics critical for understanding astrophysical effects of the nuclear symmetry energy.

We study implications of the very recently reported mass $M=2.17^{+0.11}_{-0.10}$M$_\odot$ of PSR~J0740+6620 on the Equation of State (EOS) of super-dense neutron-rich nuclear matter with respect to existing constraints on the EOS based on the mass $M=2.01\pm 0.04$M$_\odot$ of PSR~J0348+0432, the maximum tidal deformability of GW170817 and earlier results of various terrestrial nuclear laboratory experiments. The lower limit of the skewness $J_0$ measuring the stiffness of super-dense isospin-symmetric nuclear matter is raised raised from about -220 MeV to -150 MeV, reducing significantly its current uncertainty range. The lower bound of the high-density symmetry energy also increases appreciably leading to a rise of the minimum proton fraction in neutron stars at $\beta$-equilibrium from about 0 to 5\% around three times the saturation density of nuclear matter. The difficulties for some of the most widely used and previously well tested model EOSs to predict simultaneously both a maximum mass higher than 2.17 M$_\odot$ and a pressure consistent with that extracted from GW170817 present some interesting new challenges for nuclear theories.

We present a novel fitting formula for the halo concentration enhancement in chameleon $f(R)$ gravity relative to General Relativity (GR). The formula is derived by employing a large set of $N$-body simulations of the Hu-Sawicki $f(R)$ model which cover a wide range of model and cosmological parameters, resolutions and simulation box sizes. The complicated dependence of the concentration on halo mass $M$, redshift $z$, and the $f(R)$ and cosmological parameters can be combined into a simpler form that depends only on a rescaled mass $M/10^{p_2}$, with $p_2\equiv1.5\log_{10}\left[|{\bar{f}_R(z)}|/(1+z)\right]+21.64$ and $\bar{f}_R(z)$ the background scalar field at $z$, irrespective of the $f(R)$ model parameter. Our fitting formula can describe the concentration enhancement well for redshifts $z\leq3$, nearly 7 orders of magnitude in $M/10^{p_2}$ and five decades in halo mass. This is part of a series of works which aims to provide a general framework for self-consistent and unbiased tests of gravity using present and upcoming galaxy cluster surveys. The fitting formula, which is the first quantitative model for the concentration enhancement due to chameleon type modified gravity, is an important part in this framework and will allow continuous exploration of the parameter space. It can also be used to model other statistics such as the matter power spectrum.

In this paper, we review observational aspects of three common small-scale energetic events in the solar transition region (TR), namely: TR explosive events, ultraviolet bursts and jets. These events are defined in either (both) spectral or (and) imaging data. The development of multiple instruments capable of observing the TR has allowed researchers to gain numerous insights into these phenomena in recent years. These events have provided a proxy to study how mass and energy are transported between the solar chromosphere and the corona. As the physical mechanisms responsible for these small-scale events might be similar to the mechanisms responsible for large-scale phenomena, such as flares and coronal mass ejections, analysis of these events could also help our understanding of the solar atmosphere from small to large scales. The observations of these small-scale energetic events demonstrate that the TR is extremely dynamic and is a crucial layer in the solar atmosphere between the chromosphere and the corona.

Within a minimum model for neutron stars consisting of nucleons, electrons and muons at $\beta$-equilibrium using about a dozen Equation of States (EOSs) from microscopic nuclear many-body theories and 40,000 EOSs randomly generated using an explicitly isospin-dependent parametric EOS model for high-density neutron-rich nucleonic matter within its currently known uncertainty range, we study correlations among the f-mode frequency, its damping time and the tidal deformability as well as the compactness of neutron stars. Except for quark stars, both the f-mode frequency and damping time of canonical neutron stars are found to scale with the tidal deformability independent of the EOSs used. Applying the constraint on the tidal deformability of canonical neutron stars $\Lambda_{1.4}=190^{+390}_{-120}$ extracted by the LIGO+VIRGO Collaborations from their improved analyses of the GW170817 event, the f-mode frequency and its damping time of canonical neutron stars are limited to 1.67 kHz - 2.18 kHz and 0.155 s - 0.255 s, respectively, providing a useful guidance for the ongoing search for gravitational waves from the f-mode oscillations of isolated neutron stars. Moreover, assuming either or both the f-mode frequency and its damping time will be measured precisely in future observations with advanced gravitational wave detectors, we discuss how information about the mass and/or radius as well as the still rather elusive nuclear symmetry energies at supra-saturation densities may be extracted.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Applying an equation of state (EOS) with its symmetric nuclear matter (SNM) contribution and low-density symmetry energy $E_{sym}(\rho)$ constrained by heavy-ion reaction data, we calculate the dimensionless tidal deformability $\Lambda$ of neutron stars in coalescing binary systems. Corresponding to the partially constrained EOS that previously predicted a radius of 11.5 km $\leq R_{1.4} \leq$ 13.6 km for canonical neutron-star configurations, $\Lambda$ is found to be in the range of 292 $\leq\Lambda_{1.4}\leq$ 680, consistent with the very recent observation of the GW170817 event. We investigate the effect of the high-density behavior of $E_{sym}(\rho)$ on the tidal properties of neutron stars and find that while $\Lambda$ depends strongly on the details of the symmetry energy, different trends of $E_{sym}(\rho)$ lead to very similar values of $\Lambda$. In particular, the transition from stiff/soft to soft/stiff $E_{sym}(\rho)$ could yield the same $\Lambda$. Thus, measuring $\Lambda$ alone may not determine completely the density dependence of the symmetry energy. Coherent analyses of the dense neutron-rich nuclear matter EOS underlying both nuclear laboratory experiments and astrophysical observations are therefore necessary to break this degeneracy and determine precisely the details of the $E_{sym}(\rho)$.

Inspired by recent IRIS observations, we forward model the response of the Fe XXI 1354 \AA\ line to fundamental, standing, linear fast sausage modes (FSMs) in flare loops. Starting with the fluid parameters for an FSM in a straight tube with equilibrium parameters largely compatible with the IRIS measurements, we synthesize the line profiles by incorporating the non-Equilibrium Ionization (NEI) effect in the computation of the contribution function. We find that both the intensity and Doppler shift oscillate at the wave period ($P$). The phase difference between the two differs from the expected value ($90^\circ$) only slightly because NEI plays only a marginal role in determining the ionic fraction of Fe XXI in the examined dense loop. The Doppler width modulations, however, posses an asymmetry in the first and second halves of a wave period, leading to a secondary periodicity at $P/2$ in addition to the primary one at $P$. This behavior results from the competition between the broadening due to bulk flow and that due to temperature variations, with the latter being stronger but not overwhelmingly so. These expected signatures, with the exception of the Doppler width, are largely consistent with the IRIS measurements, thereby corroborating the reported detection of a fundamental FSM. The forward modeled signatures are useful for identifying fundamental FSMs in flare loops from measurements of the Fe \MyRoman{21} 1354 \AA\ line with instruments similar to IRIS, even though a much higher cadence is required for the expected behavior in the Doppler widths to be detected.

By numerically inverting the Tolman-Oppenheimer-Volkov (TOV) equation using an explicitly isospin-dependent parametric Equation of State (EOS) of dense neutron-rich nucleonic matter, a restricted EOS parameter space is established using observational constraints on the radius, maximum mass, tidal polarizability and causality condition of neutron stars (NSs). The constraining band obtained for the pressure as a function of energy (baryon) density is in good agreement with that extracted recently by the LIGO+Virgo Collaborations from their improved analyses of the NS tidal polarizability in GW170817. Rather robust upper and lower boundaries on nuclear symmetry energies are extracted from the observational constraints up to about twice the saturation density $\rho_0$ of nuclear matter. More quantitatively, the symmetry energy at $2\rho_0$ is constrained to $E_{\rm{sym}}(2\rho_0)=46.9\pm10.1$ MeV excluding many existing theoretical predictions scattered between $E_{\rm{sym}}(2\rho_0)=15$ and 100 MeV. Moreover, by studying variations of the causality surface where the speed of sound equals that of light at central densities of the most massive neutron stars within the restricted EOS parameter space, the absolutely maximum mass of neutron stars is found to be 2.40 M$_{\odot}$ approximately independent of the EOSs used. This limiting mass is consistent with findings of several recent analyses and numerical general relativity simulations about the maximum mass of the possible super-massive remanent produced in the immediate aftermath of GW170817.

Currently, FRB 121102 is the only fast radio burst source that was observed to give out bursts repeatedly. It shows a high repeating rate, with more than one hundred bursts being spotted, but with no obvious periodicity in the activities. Thanks to its repetition, the source was well localized with a subarcsecond accuracy, leading to a redshift measurement of about 0.2. FRB 121102 is a unique source that can help us understand the enigmatic nature of fast radio bursts. In this study, we analyze the characteristics of the waiting times between bursts from FRB 121102. It is found that there is a clear bimodal distribution for the waiting times. While most waiting times cluster at several hundred seconds, a small portion of the waiting times are strikingly in the range of 2--40 millisecond. More interestingly, it is found that the waiting time does not correlate with the burst intensity, either for the preceding burst or for the subsequent burst. It strongly indicates that the repeating bursts should be generated by some external mechanisms, but not internal mechanisms. As a result, the models involving collisions between small bodies and neutron stars could be competitive mechanisms for this interesting source.

Marked correlation functions, which are sensitive to the clustering of galaxies in different environments, have been proposed as constraints on modified gravity models. We present measurements of the marked correlation functions of galaxies in redshift space using 361,761 LOWZ ($z_{\rm eff} = 0.32$) galaxies from the Sloan Digital Sky Survey III (SDSS III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) and compare them to $\Lambda$CDM+General Relativity simulations. We apply mass cuts to find the best match between the redshift space autocorrelation function of subhaloes in the simulation and in the observations. We then compare the marked correlation functions, finding no significant evidence for deviations of the marked correlation functions of LOWZ galaxies from $\Lambda$CDM on scales $6 \ h^{-1}$Mpc $ \leq s \leq$ $69 \ h^{-1}$Mpc. The constraining power of marked correlation functions in our analysis is limited by our ability to model the autocorrelation function of galaxies on small scales including the effect of redshift distortions. The statistical errors are well below the differences seen between marked correlation functions of $f$(R) gravity models and $\Lambda$CDM in recent publications (Armijo et al., Hern\'{a}ndez-Aguayo et al.) indicating that improved future theoretical analyses should be able to rule out some models definitively.

We study cosmic voids in the normal-branch Dvali-Gabadadze-Porrati (nDGP) braneworld models, which are representative of a class of modified gravity theories where deviations from General Relativity are usually hidden by the Vainshtein screening in high-density environments. This screening is less efficient away from these environments, which makes voids ideally suited for testing this class of models. We use N-body simulations of $\Lambda$-cold dark matter ($\Lambda$CDM) and nGDP universes, where dark matter haloes are populated with mock galaxies that mimic the clustering and number densities of the BOSS CMASS galaxy sample. We measure the force, density and weak lensing profiles around voids identified with six different algorithms. Compared to $\Lambda$CDM, voids in nDGP are more under-dense due to the action of the fifth force that arises in these models, which leads to a faster evacuation of matter from voids. This leaves an imprint on the weak lensing tangential shear profile around nDGP voids, an effect that is particularly strong for 2D underdensities that are identified in the plane-of-the-sky. We make predictions for the feasibility of distinguishing between nDGP and $\Lambda$CDM using void lensing in upcoming large-scale surveys such as LSST and Euclid. We compare with the analysis of voids in chameleon gravity theories and find that the weak lensing signal for 3D voids is similar to nDGP, whereas for 2D voids the differences with $\Lambda$CDM are much stronger for the chameleon gravity case, a direct consequence of the different screening mechanisms operating in these theories.

Thermal inflation models (which feature two inflationary stages) can display damped primordial curvature power spectra on small scales which generate damped matter fluctuations. For a reasonable choice of parameters, thermal inflation models naturally predict a suppression of the matter power spectrum on galactic and sub-galactic scales, mimicking the effect of warm or interacting dark matter. Matter power spectra in these models are also characterised by an excess of power (w.r.t. the standard $\Lambda$CDM power spectrum) just below the suppression scale. By running a suite of N-body simulations we investigate the non-linear growth of structure in models of thermal inflation. We measure the non-linear matter power spectrum and extract halo statistics, such as the halo mass function, and compare these quantities with those predicted in the standard $\Lambda$CDM model and in other models with damped matter fluctuations. We find that the thermal inflation models considered here produce measurable differences in the matter power spectrum from $\Lambda$CDM at redshifts $z>5$, while the halo mass functions are appreciably different at all redshifts. The halo mass function at $z=0$ for thermal inflation displays an enhancement of around $\sim 20\%$ w.r.t. $\Lambda$CDM and a damping at lower halo masses, with the position of the enhancement depending on the value of the free parameter in the model. The enhancement in the halo mass function (w.r.t. $\Lambda$CDM ) increases with redshift, reaching $\sim 40\%$ at $z=5$. We also study the accuracy of the analytical Press-Schechter approach, using different filters to smooth the density field, to predict halo statistics for thermal inflation. We find that the predictions with the smooth-$k$ filter agree with the simulation results over a wider range of halo masses than is the case with other filters commonly used in the literature.

The reconstruction of the initial conditions of the Universe is an important topic in cosmology, particularly in the context of sharpening the measurement of the baryon acoustic oscillation (BAO) peak. Nonlinear reconstruction algorithms developed in recent years, when applied to late-time matter fields, can recover to a substantial degree the initial density distribution, however, when applied to sparse tracers of the matter field, the performance is poorer. In this paper we apply the Shi et al. non-linear reconstruction method to biased tracers in order to establish what factors affect the reconstruction performance. We find that grid resolution, tracer number density and mass assignment scheme all have a significant impact on the performance of our reconstruction method, with triangular-shaped-cloud (TSC) mass assignment and a grid resolution of ${\sim}1{-}2h^{-1}$ Mpc being the optimal choice. We also show that our method can be easily adapted to include generic tracer biases up to quadratic order in the reconstruction formalism. Applying the reconstruction to halo and galaxy samples with a range of tracer number densities, we find that the linear bias is by far the most important bias term, while including nonlocal and nonlinear biases only leads to marginal improvements on the reconstruction performance. Overall, including bias in the reconstruction substantially improves the recovery of BAO wiggles, down to $k\sim0.25~h\text{Mpc}^{-1}$ for tracer number densities between $2\times10^{-4}$ and $2\times10^{-3}~(h^{-1}\text{Mpc})^{-3}$.

Measurements of redshift space distortions (RSD) provide a means to test models of gravity on large-scales. We use mock galaxy catalogues constructed from large N-body simulations of standard and modified gravity models to measure galaxy clustering in redshift space. We focus our attention on two of the most representative and popular families of modified gravity models: the Hu \& Sawicki $f(R)$ gravity and the normal branch of the DGP model. The galaxy catalogues are built using a halo occupation distribution (HOD) prescription with the HOD parameters in the modified gravity models tuned to match with the number density and the real-space clustering of {\sc boss-cmass} galaxies. We employ two approaches to model RSD: the first is based on linear perturbation theory and the second models non-linear effects on small-scales by assuming standard gravity and including biasing and RSD effects. We measure the monopole to real-space correlation function ratio, the quadrupole to monopole ratio, clustering wedges and multipoles of the correlation function and use these statistics to find the constraints on the distortion parameter, $\beta$. We find that the linear model fails to reproduce the N-body simulation results and the true value of $\beta$ on scales $s < 40\Mpch$, while the non-linear modelling of RSD recovers the value of $\beta$ on the scales of interest for all models. RSD on large scales ($s\gtrsim20$-$40\Mpch$) have been found to show significant deviations from the prediction of standard gravity in the DGP models. However, the potential to use RSD to constrain $f(R)$ models is less promising, due to the different screening mechanism in this model.

We analyse the two-point and marked correlation functions of haloes and galaxies in three variants of the chameleon $f(R)$ gravity model using N-body simulations, and compare to a fiducial $\Lambda$CDM model based on general relativity (GR). Using a halo occupation distribution prescription (HOD) we populate dark matter haloes with galaxies, where the HOD parameters have been tuned such that the galaxy number densities and the real-space galaxy two-point correlation functions in the modified gravity models match those in GR to within $1\sim3\%$. We test the idea that since the behaviour of gravity is dependent on environment, marked correlation functions may display a measurable difference between the models. For this we test marks based on the density field and the Newtonian gravitational potential. We find that the galaxy marked correlation function shows significant differences measured in different models on scales smaller than $r\lesssim 20~h^{-1}$ Mpc. Guided by simulations to identify a suitable mark, this approach could be used as a new probe of the accelerated expansion of the Universe.

The strongly isospin-dependent tensor force leads to short-range correlations (SRC) between neutron-proton (deuteron-like) pairs much stronger than those between proton-proton and neutron-neutron pairs. As a result of the short-range correlations, the single-nucleon momentum distribution develops a high-momentum tail above the Fermi surface. Because of the strongly isospin-dependent short-range correlations, in neutron-rich matter a higher fraction of protons will be depleted from its Fermi sea and populate above the Fermi surface compared to neutrons. This isospin-dependent nucleon momentum distribution may have effects on: (1) nucleon spectroscopic factors of rare isotopes, (2) the equation of state especially the density dependence of nuclear symmetry energy, (3) the coexistence of a proton-skin in momentum space and a neutron-skin in coordinate space (i.e., protons move much faster than neutrons near the surface of heavy nuclei). In this talk, we discuss these features and their possible experimental manifestations. As an example, SRC effects on the nuclear symmetry energy are discussed in detail using a modified Gogny-Hartree-Fock (GHF) energy density functional (EDF) encapsulating the SRC-induced high momentum tail (HMT) in the single-nucleon momentum distribution.

Forward-modeling the emission properties in various passbands is important for confidently identifying magnetohydrodynamic waves in the structured solar corona. We examine how Non-equilibrium Ionization (NEI) affects the Extreme Ultraviolet (EUV) emissions modulated by standing fast sausage modes (FSMs) in coronal loops, taking the Fe IX 171 \AA\ and Fe XII 193 \AA\ emission lines as examples. Starting with the expressions for linear FSMs in straight cylinders, we synthesize the specific intensities and spectral profiles for the two spectral lines by incorporating the self-consistently derived ionic fractions in the relevant contribution functions. We find that relative to the case where Equilibrium Ionization (EI) is assumed, NEI considerably impacts the intensity modulations, but shows essentially no effect on the Doppler velocities or widths. Furthermore, NEI may affect the phase difference between intensity variations and those in Doppler widths for Fe XII 193 \AA\ when the line-of-sight is oblique to the loop axis. While this difference is $180^\circ$ when EI is assumed, it is $\sim 90^\circ$ when NEI is incorporated for the parameters we choose. We conclude that in addition to viewing angles and instrumental resolutions, NEI further complicates the detection of FSMs in spectroscopic measurements of coronal loops in the EUV passband.

Recent numerical studies revealed that transverse motions of coronal loops can induce the Kelvin-Helmholtz Instability (KHI). This process could be important in coronal heating because it leads to dissipation of energy at small spatial-scale plasma interactions. Meanwhile, small amplitude decayless oscillations in coronal loops have been discovered recently in observations of SDO/AIA. We model such oscillations in coronal loops and study wave heating effects, considering a kink and Alfv{\'e}n driver separately and a mixed driver at the bottom of flux tubes. Both the transverse and Alfv{\'e}n oscillations can lead to the KHI. Meanwhile, the Alfv{\'e}n oscillations established in loops will experience phase mixing. Both processes will generate small spatial-scale structures, which can help the dissipation of wave energy. Indeed, we observe the increase of internal energy and temperature in loop regions. The heating is more pronounced for the simulation containing the mixed kink and Alfv{\'e}n driver. This means that the mixed wave modes can lead to a more efficient energy dissipation in the turbulent state of the plasma and that the KHI eddies act as an agent to dissipate energy in other wave modes. Furthermore, we also obtained forward modelling results using the FoMo code. We obtained forward models which are very similar to the observations of decayless oscillations. Due to the limited resolution of instruments, neither Alfv{\'e}n modes nor the fine structures are observable. Therefore, this numerical study shows that Alfv{\'e}n modes probably can co-exist with kink modes, leading to enhanced heating.

Type Ia supernovae (SNIe) are generally accepted to act as standardisable candles, and their use in cosmology led to the first confirmation of the as yet unexplained accelerated cosmic expansion. Many of the theoretical models to explain the cosmic acceleration assume modifications to Einsteinian General Relativity which accelerate the expansion, but the question of whether such modifications also affect the ability of SNIe to be standardisable candles has rarely been addressed. This paper is an attempt to answer this question. For this we adopt a semi-analytical model to calculate SNIe light curves in non-standard gravity. We use this model to show that the average rescaled intrinsic peak luminosity -- a quantity that is assumed to be constant with redshift in standard analyses of Type Ia supernova (SNIa) cosmology data -- depends on the strength of gravity in the supernova's local environment because the latter determines the Chandrasekhar mass -- the mass of the SNIa's white dwarf progenitor right before the explosion. This means that SNIe are no longer standardisable candles in scenarios where the strength of gravity evolves over time, and therefore the cosmology implied by the existing SNIa data will be different when analysed in the context of such models. As an example, we show that the observational SNIa cosmology data can be fitted with both a model where $(\Omega_{\rm M}, \Omega_{\Lambda})=(0.62, 0.38)$ and Newton's constant $G$ varies as $G(z)=G_0(1+z)^{-1/4}$ and the standard model where $(\Omega_{\rm M}, \Omega_{\Lambda})=(0.3, 0.7)$ and $G$ is constant, when the Universe is assumed to be flat.

The intrinsic peak luminosity of Type Ia supernovae (SNIa) depends on the value of Newton's gravitational constant $G$, through the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. If the luminosity distance can be independently determined, the SNIa can be treated as a tracker to constrain the possible time variation of $G$ in different redshift ranges. The gravitational-wave (GW) standard sirens, caused by the coalescence of binary neutron stars, provide a model-independent way to measure the distance of GW events, which can be used to determine the luminosity distances of SNIa by interpolation, provided the GW and SNIa samples have similar redshift ranges. We demonstrate that combining the GW observations of third-generation detectors with SNIa data provides a powerful and model-independent way to measure $G$ in a wide redshift range, which can constrain the ratio $G/G_0$, where $G$ and $G_0$ are respectively the values in the redshift ranges $z>0.1$ and $z<0.1$, at the level of $1.5\%$.

There has been considerable interest in sausage modes in photospheric waveguides like pores and sunspots, and slow surface sausage modes (SSSMs) have been suggested to damp ufficiently rapidly to account for chromospheric heating. Working in the framework of linear resistive magnetohydrodynamics, we examine how efficient electric resistivity and resonant absorption in the cusp continuum can be for damping SSSMs in a photospheric waveguide with equilibrium parameters compatible with recent measurements of a photospheric pore. For SSSMs with the measured wavelength, we find that the damping rate due to the cusp resonance is substantially less strong than theoretically expected with the thin-boundary approximation. The damping-time-to-period ratio ($\tau/P$) we derive for standing modes, equivalent to the damping-length-to-wavelength ratio for propagating modes given the extremely weak dispersion, can reach only $\sim 180$. However, the accepted values for electric resistivity ($\eta$) correspond to a regime where both the cusp resonance and resistivity play a role. The values for $\tau/P$ attained at the largest allowed $\eta$ may reach $\sim 30$. We conclude that electric resistivity can be considerably more efficient than the cusp resonance for damping SSSMs in the pore in question, and it needs to be incorporated into future studies on the damping of SSSMs in photospheric waveguides in general.

The Intermediate Palomar Transient Factory reported the discovery of an unusual type II-P supernova iPTF14hls. Instead of a ~100-day plateau as observed for ordinary type II-P supernovae, the light curve of iPTF14hls has at least five distinct peaks, followed by a steep decline at ~1000 days since discovery. Until 500 days since discovery, the effective temperature of iPTF14hls is roughly constant at 5000-6000K . In this paper we propose that iPTF14hls is likely powered by intermittent fallback accretion. It is found that the light curve of iPTF14hls can be well fit by the usual t^{-5/3} accretion law until ~1000 days post discovery when the light curve transitions to a steep decline. To account for this steep decline, we suggest a power-law density profile for the late accreted material, rather than the constant profile as appropriated for the t^{-5/3} accretion law. Detailed modeling indicates that the total fallback mass is ~0.2M_{sun}, with an ejecta mass M_{ej}~21M_{sun}. We find the third peak of the light curve cannot be well fit by the fallback model, indicating that there could be some extra rapid energy injection. We suggest that this extra energy injection may be a result of a magnetic outburst if the central object is a neutron star. These results indicate that the progenitor of iPTF14hls could be a massive red supergiant.

The recent discovery of gravitational waves marks the culmination of a sequence of successful tests of the general theory of relativity (GR) since its formulation in 1915. Yet these tests remain confined to the scale of stellar systems or the strong gravity regime. A departure from GR on larger, cosmological scales has been advocated by the proponents of modified gravity theories as an alternative to the Cosmological Constant to account for the observed cosmic expansion history. While indistinguishable in these terms by construction, such models on the other hand yield distinct values for the linear growth rate of density perturbations and, as a consequence, for the associated galaxy peculiar velocity field. Measurements of the resulting anisotropy of galaxy clustering, when spectroscopic redshifts are used to derive distances, have thus been proposed as a powerful probe of the validity of GR on cosmological scales. However, despite significant effort in modelling such redshift space distortions, systematic errors remain comparable to current statistical uncertainties. Here, we present the results of a different forward-modelling approach, which fully exploits the sensitivity of the galaxy velocity field to modifications of GR. We use state-of-the-art, high-resolution N-body simulations of a standard GR and a compelling f(R) model, one of GR's simplest variants, to build simulated catalogues of stellar-mass-selected galaxies through a robust match to the Sloan Digital Sky Survey observations. We find that, well within the uncertainty of this technique, f(R) fails to reproduce the observed redshift-space clustering on scales 1-10 Mpc/h. Instead, the standard LCDM GR model agrees impressively well with the data. This result provides a strong confirmation, on cosmological scales, of the robustness of Einstein's general theory of relativity.

What can we learn about the density dependence of nuclear symmetry energy $E_{\rm{sym}}(\rho )$ from precise measurements of the radius ($R_{\rm{1.4}}$) and/or tidal polarizability ($\Lambda_{1.4}$) of canonical neutron stars (NSs) with a mass of 1.4 M$_\odot$? With the $E_{\rm{sym}}(\rho )$ parameterized using three parameters $L$, $K_{\rm{sym}}$, and $J_{\rm{sym}}$ which have the asymptotic meaning of being respectively the slope, curvature, and skewness of symmetry energy at saturation density, we found that, while both the $R_{\rm{1.4}}$ and $\Lambda_{1.4}$ depend strongly on the slope $L$, the $K_{\rm{sym}}$ and $J_{\rm{sym}}$ parameters characterizing the high-density behavior of $E_{\rm{sym}}(\rho )$ also play appreciable roles. Thus, there is not a simple relation between the $\Lambda_{\rm{1.4}}$/$R_{\rm{1.4}}$ and $L$ alone. Precise measurements of just the $\Lambda_{\rm{1.4}}$ and $R_{\rm{1.4}}$ can not completely determine the $E_{\rm{sym}}(\rho )$ but limit combinations of its parameters. In particular, stringent constraints approximately independent of the $J_{\rm{sym}}$ on the $L$-$K_{\rm{sym}}$ correlations can be obtained. However, infinite combinations of the larger (smaller) $L$ and smaller (larger) $K_{\rm{sym}}$ can lead to the same $\Lambda_{\rm{1.4}}$ and $R_{\rm{1.4}}$. Additional observables including those from terrestrial nuclear experiments are thus necessary to break this degeneracy in order to completely determine the density dependence of nuclear symmetry energy $E_{\rm{sym}}(\rho )$.

Qualitative dynamics of three different loop quantizations of spatially flat isotropic and homogeneous models is studied using effective spacetime description of the underlying quantum geometry. These include the standard loop quantum cosmology (LQC), its recently revived modification (referred to as mLQC-I), and another related modification of LQC (mLQC-II) whose dynamics is studied in detail for the first time. Various features of LQC, including quantum bounce and pre-inflationary dynamics, are found to be shared with the mLQC-I and mLQC-II models. We study universal properties of dynamics for chaotic inflation, fractional monodromy inflation, Starobinsky potential, non-minimal Higgs inflation, and an exponential potential. We find various critical points and study their stability, which reveal various qualitative similarities in the post-bounce phase for all these models. The pre-bounce qualitative dynamics of LQC and mLQC-II turns out to be very similar, but is strikingly different from that of mLQC-I. In the dynamical analysis, some of the fixed points turn out to be degenerate for which center manifold theory is used. For all these potentials, non-perturbative quantum gravitational effects always result in a non-singular inflationary scenario with a phase of super-inflation succeeded by the conventional inflation. We show the existence of inflationary attractors, and obtain scaling solutions in the case of the exponential potential. Since all of the models agree with general relativity at late times, our results are also of use in classical theory where qualitative dynamics of some of the potentials has not been studied earlier.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We investigate the nonlinear evolution of structure in variants of the standard cosmological model which display damped density fluctuations relative to cold dark matter (e.g. in which cold dark matter is replaced by warm or interacting DM). Using N-body simulations, we address the question of how much information is retained from different scales in the initial linear power spectrum following the nonlinear growth of structure. We run a suite of N-body simulations with different initial linear matter power spectra to show that, once the system undergoes nonlinear evolution, the shape of the linear power spectrum at high wavenumbers does not affect the non-linear power spectrum, while it still matters for the halo mass function. Indeed, we find that linear power spectra which differ from one another only at wavenumbers larger than their half-mode wavenumber give rise to (almost) identical nonlinear power spectra at late times, regardless of the fact that they originate from different models with damped fluctuations. On the other hand, the halo mass function is more sensitive to the form of the linear power spectrum. Exploiting this result, we propose a two parameter model of the transfer function in generic damped scenarios, and show that this parametrisation works as well as the standard three parameter models for the scales on which the linear spectrum is relevant.

Cosmic voids are an important probe of large-scale structure that can constrain cosmological parameters and test cosmological models. We present a new paradigm for void studies: void detection in weak lensing convergence maps. This approach identifies objects that relate directly to our theoretical understanding of voids as underdensities in the total matter field and presents several advantages compared to the customary method of finding voids in the galaxy distribution. We exemplify this approach by identifying voids using the weak lensing peaks as tracers of the large-scale structure. We find self-similarity in the void abundance across a range of peak signal-to-noise selection thresholds. The voids obtained via this approach give a tangential shear signal up to $\sim40$ times larger than voids identified in the galaxy distribution.

We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions.

The splashback radius is a physical scale in dark matter halos that is set by the gravitational dynamics of recently accreted shells. We use analytical models and N-body simulations to study the dependence of splashback on dark energy and screened modified gravity theories. In modified gravity models, the transition from screened to unscreened regions typically occurs in the cluster outskirts, suggesting potentially observable signatures in the splashback feature. We investigate the location of splashback in both chameleon and Vainshtein screened models and find significant differences compared with $\Lambda$CDM predictions. We also find an interesting interplay between dynamical friction and modified gravity, providing a distinctive signature for modified gravity models in the behavior of the splashback feature as a function of galaxy luminosity.

Fast radio bursts (FRBs) are bright radio pulses from the sky with millisecond durations and Jansky-level flux densities. Their origins are still largely uncertain. Here we suggest a new model for FRBs. We argue that the collision of a white dwarf with a black hole can give birth to a transient accretion disk, from which powerful episodic magnetic blobs will be launched. The collision between two consecutive magnetic blobs can result in a catastrophic magnetic reconnection, which releases a large amount of free magnetic energy and forms a forward shock. The shock propagates through the cold magnetized plasma within the blob in the collision region, radiating through synchrotron maser mechanism, which is responsible for a non-repeating FRB signal. Our calculations show that the theoretical energetics, radiation frequency, duration timescale, and event rate can be well consistent with the observational characteristic of FRBs.

It has been well-established that particular centrophilic orbital families in non-spherical galaxies can, in principle, drive a black hole binary to shrink its orbit through three-body scattering until the black holes are close enough to strongly emit gravitational waves. Most of these studies rely on orbital analysis of a static SMBH-embedded galaxy potential to support this view; it is not clear, however, how these orbits transform as the second SMBH enters the center, so our understanding of which orbits actually interact with a SMBH binary is not ironclad. Here, we analyze two flattened galaxy models, one with a single SMBH and one with a binary, to determine which orbits actually do interact with the SMBH binary and how they compare with the set predicted in single SMBH-embedded models. We find close correspondence between the centrophilic orbits predicted to interact with the binary and those that are actually scattered by the binary, in terms of energy and Lz distribution, where Lz is the z component of a stellar particle's angular momentum. Of minor note: because of the larger mass, the binary SMBH has a radius of influence about 4 times larger than in the single SMBH model, which allows the binary to draw from a larger reservoir of orbits to scatter. Of the prediction particles and scattered particles, nearly half have chaotic orbits, 40% have fx:fy=1:1 orbits, 10% have other resonant orbits.

Two main models have been developed to explain the mechanisms of release, heating and acceleration of the nascent solar wind, the wave-turbulence-driven (WTD) models and reconnection-loop-opening (RLO) models, in which the plasma release processes are fundamentally different. Given that the statistical observational properties of helium ions produced in magnetically diverse solar regions could provide valuable information for the solar wind modelling, we examine the statistical properties of the helium abundance (A_He) and the speed difference between helium ions and protons (v_alpha,p) for coronal holes (CHs), active regions (ARs) and the quiet Sun (QS). We find bimodal distributions in the space of A_He and v_alpha,p/v_A (where v_A is the local Alfven speed)for the solar wind as a whole. The CH wind measurements are concentrated at higher A_He and v_alpha,p/v_A values with a smaller A_He distribution range, while the AR and QS wind is associated with lower A_He and v_alpha,p/v_A, and a larger A_He distribution range. The magnetic diversity of the source regions and the physical processes related to it are possibly responsible for the different properties of A_He and v_alpha,p/v_A. The statistical results suggest that the two solar wind generation mechanisms, WTD and RLO, work in parallel in all solar wind source regions. In CH regions WTD plays a major role, whereas the RLO mechanism is more important in AR and QS.

The brightness of the multi-wavelength afterglow of GRB 170817A is increasing unexpectedly even $\sim$ 160 days after the associated gravitational burst. Here we suggest that the brightening can be caused by a late-time energy injection process. We use an empirical expression to mimic the evolution of the injection luminosity, which is consisted of a power-law rising phase and a power-law decreasing phase. It is found that the power-law indices of the two phases are $0.92$ and $-2.8$, respectively, with the peak time of the injection being $\sim 110$ days. The energy injection could be due to some kind of accretion, with the total accreted mass being $\sim 0.006 M_\odot$. However, normal fall-back accretion, which usually lasts for a much shorter period, cannot provide a natural explanation. Our best-fit decaying index of $-2.8$ is also at odds with the expected value of $-5/3$ for normal fall-back accretion. Noting that the expansion velocities of the kilonova components associated with GW170817 is $0.1-0.3\,c$, we argue that there should also be some ejecta with correspondingly lower velocities during the coalescence of the double neutron star system. They are bound by the gravitational well of the remnant central compact object and might be accreted at a timescale of about 100 days, providing a reasonable explanation for the energy injection. Detailed studies on the long-lasting brightening of GRB 170817A thus may provide useful information on the matter ejection during the merger process of binary neutron stars.

Within the parameter space of equation of state (EOS) of dense neutron-rich matter limited by existing constraints mainly from terrestrial nuclear experiments, we investigate how the neutron star maximum mass $M_{\rm{max}}>2.01\pm0.04$ M$_\odot$, radius $10.62<R_{\rm{1.4}}< 12.83$ km and tidal deformability $\Lambda_{1.4}\leq800$ of canonical neutron stars all together constrain the EOS of dense neutron-rich nucleonic matter. While the 3-D parameter space of $K_{\rm{sym}}$ (curvature of nuclear symmetry energy), $J_{\rm{sym}}$ and $J_0$ (skewness of the symmetry energy and EOS of symmetric nuclear matter, respectively) are narrowed down significantly by the observational constraints, more data are needed to pin down the individual values of $K_{\rm{sym}}$, $J_{\rm{sym}}$ and $J_0$ with quantified uncertainties. The $J_0$ largely controls the maximum mass of neutron stars. While the EOS with $J_0=0$ is sufficiently stiff to support neutron stars as massive as 2.37 M$_{\odot}$, to support the hyperthetical ones as massive as 2.74 M$_{\odot}$ (composite mass of GW170817) requires $J_0$ to be larger than its currently known maximum value of about 400 MeV and beyond the causality limit. The upper limit on the tidal deformability of $\Lambda_{1.4}=800$ from the recent observation of GW170817 is found to provide upper limits on some EOS parameters consistent with but far less restrictive than the existing constraints of other observables studied.

We present a systematic study of the cosmological dynamics resulting from an effective Hamiltonian, recently derived in loop quantum gravity using Thiemann's regularization and earlier obtained in loop quantum cosmology (LQC) by keeping the Lorentzian term explicit in the Hamiltonian constraint. We show that quantum geometric effects result in higher than quadratic corrections in energy density in comparison to LQC causing a non-singular bounce. Dynamics can be described by the Hamilton's or the Friedmann-Raychaudhuri equations, but the map between the two descriptions is not one-to-one. A careful analysis resolves the tension on symmetric versus asymmetric bounce in this model, showing that the bounce must be asymmetric and symmetric bounce is physically inconsistent, in contrast to the standard LQC. In addition, the current observations only allow a scenario where the pre-bounce branch is asymptotically de Sitter, similar to a quantization of the Schwarzschild interior in LQC, and the post-bounce branch yields the classical general relativity. For a quadratic potential, we find that a slow-roll inflation generically happens after the bounce, which is quite similar to what happens in LQC.

We propose a new filter, a smooth-$k$ space filter, to use in the Press-Schechter approach to model the dark matter halo mass function which overcomes shortcomings of other filters. We test this against the mass function measured in N-body simulations. We find that the commonly used sharp-$k$ filter fails to reproduce the behaviour of the halo mass function at low masses measured from simulations of models with a sharp truncation in the linear power spectrum. We show that the predictions with our new filter agree with the simulation results over a wider range of halo masses for both damped and undamped power spectra than is the case with the sharp-$k$ and real-space top-hat filters.

We use high spatial and temporal resolution observations, simultaneously obtained with the New Vacuum Solar Telescope and Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory, to investigate the high-frequency oscillations above a sunspot umbra. A novel time--frequency analysis method, namely the synchrosqueezing transform (SST), is employed to represent their power spectra and to reconstruct the high-frequency signals at different solar atmospheric layers. A validation study with synthetic signals demonstrates that SST is capable to resolving weak signals even when their strength is comparable with the high-frequency noise. The power spectra, obtained from both SST and the Fourier transform, of the entire umbral region indicate that there are significant enhancements between 10 and 14 mHz (labeled as 12 mHz) at different atmospheric layers. Analyzing the spectrum of a photospheric region far away from the umbra demonstrates that this 12~mHz component exists only inside the umbra. The animation based on the reconstructed 12 mHz component in AIA 171 \AA\ illustrates that an intermittently propagating wave first emerges near the footpoints of coronal fan structures, and then propagates outward along the structures. A time--distance diagram, coupled with a subsonic wave speed ($\sim$ 49 km s$^{-1}$), highlights the fact that these coronal perturbations are best described as upwardly propagating magnetoacoustic slow waves. Thus, we first reveal the high-frequency oscillations with a period around one minute in imaging observations at different height above an umbra, and these oscillations seem to be related to the umbral perturbations in the photosphere.

The Interface Region Imaging Spectrograph (IRIS) reveals numerous small-scale (sub-arcsecond) brightenings that appear as bright dots sparkling the solar transition region in active regions. Here, we report a statistical study on these transition region bright dots. We use an automatic approach to identify 2742 dots in a Si IV raster image. We find that the average spatial size of the dots is 0.8 arcsec$^2$ and most of them are located in the faculae area. Their Doppler velocities obtained from the Si IV 1394 {\AA} line range from -20 to 20 km/s. Among these 2742 dots, 1224 are predominantly blue-shifted and 1518 are red-shifted. Their nonthermal velocities range from 4 to 50 km/s with an average of 24 km/s. We speculate that the bright dots studied here are small-scale impulsive energetic events that can heat the active region corona.

Recent follow-up observations of the binary neutron star (NS) merging event GW170817/SGRB 170817A reveal that its X-ray/optical/radio emissions are brightening continuously up to $\sim 100$ days post-merger. This late-time brightening is unexpected from the kilonova model or the off-axis top-hat jet model for gamma-ray burst afterglows. In this paper, by assuming that the merger remnant is a long-lived NS, we propose that the interaction between an electron-positron-pair ($e^+e^-$) wind from the central NS and the jet could produce a long-lived reverse shock, from which a new emission component would rise and can interpret current observations well. The magnetic-field-induced ellipticity of the NS is taken to be $4 \times 10^{-5}$ in our modeling, so that the braking of the NS is mainly through the gravitational wave (GW) radiation rather than the magnetic dipole radiation, and the emission luminosity at early times would not exceed the observational limits. In our scenario, since the peak time of the brightening is roughly equal to the spin-down time scale of the NS, the accurate peak time may help constrain the ellipticity of the remnant NS. We suggest that radio polarization observations of the brightening would help to distinguish our scenario from other scenarios. Future observations on a large sample of short gamma-ray burst afterglows or detections of GW signals from merger remnants would test our scenario.

Impulsively generated sausage wave trains in coronal structures are important for interpreting a substantial number of observations of quasi-periodic signals with quasi-periods of order seconds. We have previously shown that the Morlet spectra of these wave trains in coronal tubes depend crucially on the dispersive properties of trapped sausage waves, the existence of cutoff axial wavenumbers and the monotonicity of the dependence of the axial group speed on the axial wavenumber in particular. This study examines the difference a slab geometry may introduce, for which purpose we conduct a comprehensive eigenmode analysis, both analytically and numerically, on trapped sausage modes in coronal slabs with a considerable number of density profiles. For the profile descriptions examined, coronal slabs can trap sausage waves with longer axial wavelengths, and the group speed approaches the internal Alfv\'en speed more rapidly at large wavenumbers in the cylindrical case. However, common to both geometries, cutoff wavenumbers exist only when the density profile falls sufficiently rapidly at distances far from coronal structures. Likewise, the monotonicity of the group speed curves depends critically on the profile steepness right at the structure axis. Furthermore, the Morlet spectra of the wave trains are shaped by the group speed curves for coronal slabs and tubes alike. Consequently, we conclude that these spectra have the potential for telling the sub-resolution density structuring inside coronal structures, although their detection requires an instrumental cadence of better than $\sim 1$ second.

We present a systematic comparison of several existing and new void finding algorithms, focusing on their potential power to test a particular class of modified gravity models - chameleon $f(R)$ gravity. These models deviate from standard General Relativity (GR) more strongly in low-density regions and thus voids are a promising venue to test them. We use Halo Occupation Distribution (HOD) prescriptions to populate haloes with galaxies, and tune the HOD parameters such that the galaxy two-point correlation functions are the same in both f(R) and GR models. We identify both 3D voids as well as 2D underdensities in the plane-of-the-sky to find the same void abundance and void galaxy number density profiles across all models, which suggests that they do not contain much information beyond galaxy clustering. However, the underlying void dark matter density profiles are significantly different, with f(R) voids being more underdense than GR ones, which leads to f(R) voids having a larger tangential shear signal than their GR analogues. We investigate the potential of each void finder to test f(R) models with near-future lensing surveys such as EUCLID and LSST. The 2D voids have the largest power to probe f(R) gravity, with a LSST analysis of tunnel (which is a new type of 2D underdensity introduced here) lensing distinguishing at 80 and 11$\sigma$ (statistical error) f(R) models with $|f_{R0}|=10^{-5}$ and $10^{-6}$ from GR.

We analytically investigate the Schwinger pair production in the de Sitter background by using {\em the uniform asymptotic approximation method}, and show that the equation of motion in general has two turning points, and the nature of these points could be single, double, real or complex, depending on the choice of the free parameters involved in the theory. Different natures of these points lead to different electric currents. In particular, when $\beta \equiv m^2/H^2-9/4$ is positive, both turning points are complex, and the electric current due to the Schwinger process is highly suppressed, where $m$ and $H$ denote, respectively, the mass of the particle and the Hubble parameter. For the turning points to be real, it is necessary to have $\beta < 0$, and the more negative of $\beta$, the easier to produce particles. In addition, when $\beta < 0$, we also study the particle production when the electric field $E$ is very weak. We find that the electric current in this case is proportional to $E^{1/2 - \sqrt{|\beta|}}$, which is strongly enhanced in the weak electric field limit when $m < \sqrt{2} H$.

We propose a new framework for testing gravity using cluster observations, which aims to provide an unbiased constraint on modified gravity models from Sunyaev Zel'dovich (SZ) and X-ray cluster counts and the cluster gas fraction, among other possible observables. Focusing on a popular $ f(R) $ model of gravity, we propose a novel procedure to recalibrate mass scaling relations from $ \Lambda $CDM to $ f(R) $ gravity for SZ and X-ray cluster observables. We find that the complicated modified gravity effects can be simply modelled as a dependence on a combination of the background scalar field and redshift, $ f_R(z)/(1+z) $, regardless of the $ f(R) $ model parameter. By employing a large suite of N-body simulations, we demonstrate that a theoretically derived tanh fitting formula is in excellent agreement with the dynamical mass enhancement of dark matter haloes for a large range of background field parameters and redshifts. Our framework is sufficiently flexible to allow for tests of other models and inclusion of further observables. The one-parameter description of the dynamical mass enhancement can have important implications on the theoretical modelling of observables and on practical tests of gravity.

In theories of modified gravity with the chameleon screening mechanism, the strength of the fifth force depends on environment. This induces an environment dependence of structure formation, which differs from $\Lambda$CDM. We show that these differences can be captured by the marked correlation function. With the galaxy correlation functions and number densities calibrated to match between $f(R)$ and $\Lambda$CDM models in simulations, we show that the marked correlation functions from using either the local density or halo mass as the marks encode extra information, which can be used to test these theories. We discuss possible applications of these statistics in observations.

We examine the dispersive properties of linear fast standing modes in transversely nonuniform solar coronal slabs with finite gas pressure, or, equivalently, finite plasma beta. We derive a generic dispersion relation governing fast waves in coronal slabs for which the continuous transverse distributions of the physical parameters comprise a uniform core, a uniform external medium, and a transition layer (TL) in between. The profiles in the TL are allowed to be essentially arbitrary. Restricting ourselves to the first several branches of fast modes, which are of most interest from the observational standpoint, we find that a finite plasma beta plays an at most marginal role in influencing the periods ($P$), damping times ($\tau$), and critical longitudinal wavenumbers ($k_{\rm c}$), when both $P$ and $\tau$ are measured in units of the transverse fast time. However, these parameters are in general significantly affected by how the TL profiles are described. We conclude that, for typical coronal structures, the dispersive properties of the first several branches of fast standing modes can be evaluated with the much simpler theory for cold slabs provided that the transverse profiles are properly addressed and the transverse Alfv\'en time in cold MHD is replaced with the transverse fast time.

We present observations of persistent oscillations of some bright features in the upper-chromosphere/transition-region above sunspots taken by IRIS SJ 1400 {\AA} and upward propagating quasi-periodic disturbances along coronal loops rooted in the same region taken by AIA 171 {\AA} passband. The oscillations of the features are cyclic oscillatory motions without any obvious damping. The amplitudes of the spatial displacements of the oscillations are about 1 $^{"}$. The apparent velocities of the oscillations are comparable to the sound speed in the chromosphere, but the upward motions are slightly larger than that of the downward. The intensity variations can take 24-53% of the background, suggesting nonlinearity of the oscillations. The FFT power spectra of the oscillations show dominant peak at a period of about 3 minutes, in consistent with the omnipresent 3 minute oscillations in sunspots. The amplitudes of the intensity variations of the upward propagating coronal disturbances are 10-15% of the background. The coronal disturbances have a period of about 3 minutes, and propagate upward along the coronal loops with apparent velocities in a range of 30-80 km/s. We propose a scenario that the observed transition region oscillations are powered continuously by upward propagating shocks, and the upward propagating coronal disturbances can be the recurrent plasma flows driven by shocks or responses of degenerated shocks that become slow magnetic-acoustic waves after heating the plasma in the coronal loops at their transition-region bases.

We present high-resolution observations of a magnetic reconnection event in the solar atmosphere taken with the New Vacuum Solar Telescope, AIA and HMI. The reconnection event occurred between the threads of a twisted arch filament system (AFS) and coronal loops. Our observations reveal that the relaxation of the twisted AFS drives some of its threads to encounter the coronal loops, providing inflows of the reconnection. The reconnection is evidenced by flared X-shape features in the AIA images, a current-sheet-like feature apparently connecting post-reconnection loops in the \halpha$+$1 \AA\ images, small-scale magnetic cancellation in the HMI magnetograms and flows with speeds of 40--80 km/s along the coronal loops. The post-reconnection coronal loops seen in AIA 94 \AA\ passband appear to remain bright for a relatively long time, suggesting that they have been heated and/or filled up by dense plasmas previously stored in the AFS threads. Our observations suggest that the twisted magnetic system could release its free magnetic energy into the upper solar atmosphere through reconnection processes. While the plasma pressure in the reconnecting flux tubes are significantly different, the reconfiguration of field lines could result in transferring of mass among them and induce heating therein.

We report on high-resolution imaging and spectral observations of eruptions of a spiral structure in the transition region, which were taken with the Interface Region Imaging Spectrometer (IRIS), the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI). The eruption coincided with the appearance of two series of jets, with velocities comparable to the Alfv\'en speeds in their footpoints. Several pieces of evidence of magnetic braiding in the eruption are revealed, including localized bright knots, multiple well-separated jet threads, transition region explosive events and the fact that all these three are falling into the same locations within the eruptive structures. Through analysis of the extrapolated three-dimensional magnetic field in the region, we found that the eruptive spiral structure corresponded well to locations of twisted magnetic flux tubes with varying curl values along their lengths. The eruption occurred where strong parallel currents, high squashing factors, and large twist numbers were obtained. The electron number density of the eruptive structure is found to be $\sim3\times10^{12}$ cm$^{-3}$, indicating that significant amount of mass could be pumped into the corona by the jets. Following the eruption, the extrapolations revealed a set of seemingly relaxed loops, which were visible in the AIA 94 \AA\ channel indicating temperatures of around 6.3 MK. With these observations, we suggest that magnetic braiding could be part of the mechanisms explaining the formation of solar eruption and the mass and energy supplement to the corona.

Various kinds of isovector nucleon effective masses are used in the literature to characterize the momentum/energy dependence of the nucleon symmetry potential or self-energy due to the space/time non-locality of the underlying isovector strong interaction in neutron-rich nucleonic matter. The multifaceted studies on nucleon isovector effective masses are multi-disciplinary in nature. Besides structures, masses and low-lying excited states of nuclei as well as nuclear reactions, studies of the isospin dependence of short-range correlations in nuclei from scatterings of high-energy electrons and protons on heavy nuclei also help understand nucleon effective masses especially the so-called E-mass in neutron-rich matter. A thorough understanding of all kinds of nucleon effective masses has multiple impacts on many interesting issues in both nuclear physics and astrophysics. We review some of the significant progresses made in recent years by the nuclear physics community in resolving some of the hotly debated and longstanding issues regarding nucleon effective masses especially in dense neutron-rich matter.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

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

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

Broad-lined type Ic supernovae (SNe Ic-BL) are a subclass of rare core collapse SNe whose energy source is debated in the literature. Recently a series of investigations on SNe Ic-BL with the magnetar (plus 56Ni) model were carried out. Evidence for magnetar formation was found for the well-observed SNe Ic-BL 1998bw and 2002ap. In this paper we systematically study a large sample of SNe Ic-BL not associated with gamma-ray bursts. We use photospheric velocity data determined in a homogeneous way. We find that the magnetar+56Ni model provides a good description of the light curves and velocity evolution of our sample of SNe Ic-BL, although some SNe (not all) can also be described by the pure-magnetar model or by the two-component pure-56Ni model (3 out of 12 are unlikely explained by two-component model). In the magnetar+56Ni model, the amount of 56Ni required to explain their luminosity is significantly reduced, and the derived initial explosion energy is, in general, in accordance with neutrino heating. Some correlations between different physical parameters are evaluated and their implications regarding magnetic field amplification and the total energy reservoir are discussed.

A theoretically interesting and practically important question in cosmology is the reconstruction of the initial density distribution provided a late-time density field. This is a long-standing question with a revived interest recently, especially in the context of optimally extracting the baryonic acoustic oscillation (BAO) signals from observed galaxy distributions. We present a new efficient method to carry out this reconstruction, which is based on numerical solutions to the nonlinear partial differential equation that governs the mapping between the initial Lagrangian and final Eulerian coordinates of particles in evolved density fields. This is motivated by numerical simulations of the quartic Galileon gravity model, which has similar equations that can be solved effectively by multigrid Gauss-Seidel relaxation. The method is based on mass conservation, and does not assume any specific cosmological model. Our test shows that it has a performance comparable to that of state-of-the-art algorithms which were very recently put forward in the literature, with the reconstructed density field over $\sim80\%$ ($50\%$) correlated with the initial condition at $k\lesssim0.6h/{\rm Mpc}$ ($1.0h/{\rm Mpc}$). With an example, we demonstrate that this method can significantly improve the accuracy of BAO reconstruction.

We investigate the impact of thermal velocities in N-body simulations of structure formation in warm dark matter models. Adopting the commonly used approach of adding thermal velocities, randomly selected from a Fermi-Dirac distribution, to the gravitationally-induced velocities of the simulation particles, we compare the matter and velocity power spectra measured from CDM and WDM simulations, in the latter case with and without thermal velocities. This prescription for adding thermal velocities introduces numerical noise into the initial conditions, which influences structure formation. At early times, the noise affects dramatically the power spectra measured from simulations with thermal velocities, with deviations of the order of $\sim \mathcal{O}(10)$ (in the matter power spectra) and of the order of $\sim \mathcal{O}(10^2)$ (in the velocity power spectra) compared to those extracted from simulations without thermal velocities. At late times, these effects are less pronounced with deviations of less than a few percent. Increasing the resolution of the N-body simulation shifts these discrepancies to higher wavenumbers. We also find that spurious haloes start to appear in simulations which include thermal velocities at a mass that is $\sim$3 times larger than in simulations without thermal velocities.

We present an analysis of galaxy-galaxy weak gravitational lensing (GGL) in chameleon $f(R)$ gravity - a leading candidate of non-standard gravity models. For the analysis we have created mock galaxy catalogues based on dark matter haloes from two sets of numerical simulations, using a halo occupation distribution (HOD) prescription which allows a redshift dependence of galaxy number density. To make a fairer comparison between the $f(R)$ and $\Lambda$CDM models, their HOD parameters are tuned so that the galaxy two-point correlation functions in real space (and therefore the projected two-point correlation functions) match. While the $f(R)$ model predicts an enhancement of the convergence power spectrum by up to $\sim30\%$ compared to the standard $\Lambda$CDM model with the same parameters, the maximum enhancement of GGL is only half as large and less than 5\% on separations above $\sim1$-$2h^{-1}$Mpc, because the latter is a cross correlation of shear (or matter, which is more strongly affected by modified gravity) and galaxy (which is weakly affected given the good match between galaxy auto correlations in the two models) fields. We also study the possibility of reconstructing the matter power spectrum by combination of GGL and galaxy clustering in $f(R)$ gravity. We find that the galaxy-matter cross correlation coefficient remains at unity down to $\sim2$-$3h^{-1}$Mpc at relevant redshifts even in $f(R)$ gravity, indicating joint analysis of GGL and galaxy clustering can be a powerful probe of matter density fluctuations in chameleon gravity. The scale dependence of the model differences in their predictions of GGL can potentially allow to break the degeneracy between $f(R)$ gravity and other cosmological parameters such as $\Omega_m$ and $\sigma_8$.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

In this paper, we study the non-projectable 2d Ho\v{r}ava gravity coupled with a non-relativistic scalar field, where the coupling is in general non-minimal and of the form $f(\phi)R$, where $f(\phi)$ is an arbitrary function of the scalar field $\phi$, and $R$ denotes the 2d Ricci scalar. In particular, we first investigate the Hamiltonian structure, and show that there are two-first and two-second class constraints, similar to the pure gravity case, but now the local degree of freedom is one, due to the presence of the scalar field. Then, we present various exact stationary solutions of this coupled system, and find that some of them represent black holes but now with universal horizons as their boundaries. At these horizons, the Hawking radiations are thermal with temperatures proportional to their surface gravities, which normally depend on the non-linear dispersion relations of the particles radiated, similar to the (3+1)-dimensional case.

We present a detection of 89 candidates of ultra-diffuse galaxies (UDGs) in a 4.9 degree$^2$ field centered on the Hickson Compact Group 95 (HCG 95) using deep $g$- and $r$-band images taken with the Chinese Near Object Survey Telescope. This field contains one rich galaxy cluster (Abell 2588 at $z$=0.199) and two poor clusters (Pegasus I at $z$=0.013 and Pegasus II at $z$=0.040). The 89 candidates are likely associated with the two poor clusters, giving about 50 $-$ 60 true UDGs with a half-light radius $r_{\rm e} > 1.5$ kpc and a central surface brightness $\mu(g,0) > 24.0$ mag arcsec$^{-2}$. Deep $z$'-band images are available for 84 of the 89 galaxies from the Dark Energy Camera Legacy Survey (DECaLS), confirming that these galaxies have an extremely low central surface brightness. Moreover, our UDG candidates are spread over a wide range in $g-r$ color, and $\sim$26% are as blue as normal star-forming galaxies, which is suggestive of young UDGs that are still in formation. Interestingly, we find that one UDG linked with HCG 95 is a gas-rich galaxy with H I mass $1.1 \times 10^{9} M_{\odot}$ detected by the Very Large Array, and has a stellar mass of $M_\star \sim 1.8 \times 10^{8}$ $M_{\odot}$. This indicates that UDGs at least partially overlap with the population of nearly dark galaxies found in deep H I surveys. Our results show that the high abundance of blue UDGs in the HCG 95 field is favored by the environment of poor galaxy clusters residing in H I-rich large-scale structures.

Integrated pulse profiles at 8.6~GHz obtained with the Shanghai Tian Ma Radio Telescope (TMRT) are presented for a sample of 26 pulsars. Mean flux densities and pulse width parameters of these pulsars are estimated. For eleven pulsars these are the first high-frequency observations and for a further four, our observations have a better signal-to-noise ratio than previous observations. For one (PSR J0742-2822) the 8.6~GHz profiles differs from previously observed profiles. A comparison of 19 profiles with those at other frequencies shows that in nine cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping, whereas in the other ten the separation is nearly constant. Different spectral indices of profile components lead to the variation of integrated pulse profile shapes with frequency. In seven pulsars with multi-component profiles, the spectral indices of the central components are steeper than those of the outer components. For the 12 pulsars with multi-component profiles in the high-frequency sample, we estimate the core width using gaussian fitting and discuss the width-period relationship.

How the solar corona is heated to high temperatures remains an unsolved mystery in solar physics. In the present study we analyse observations of 50 whole active-region loops taken with the Extreme-ultraviolet Imaging Spectrometer (EIS) on board the Hinode satellite. Eleven loops were classified as cool (<1 MK) and 39 as warm (1-2 MK) loops. We study their plasma parameters such as densities, temperatures, filling factors, non-thermal velocities and Doppler velocities. We combine spectroscopic analysis with linear force-free magnetic-field extrapolation to derive the three-dimensional structure and positioning of the loops, their lengths and heights as well as the magnetic field strength along the loops. We use density-sensitive line pairs from Fe XII, Fe XIII, Si X and Mg VII ions to obtain electron densities by taking special care of intensity background-subtraction. The emission-measure loci method is used to obtain the loop temperatures. We find that the loops are nearly isothermal along the line-of-sight. Their filling factors are between 8% and 89%. We also compare the observed parameters with the theoretical RTV scaling law. We find that most of the loops are in an overpressure state relative to the RTV predictions. In a followup study, we will report a heating model of a parallel-cascade-based mechanism and will compare the model parameters with the loop plasma and structural parameters derived here.

The morphological properties of large scale structure of the Universe can be fully described by four Minkowski functionals (MFs), which provide important complementary information to other statistical observables such as the widely used 2-point statistics in configuration and Fourier spaces. In this work, for the first time, we present the differences in the morphology of large scale structure caused by modifications to general relativity (to address the cosmic acceleration problem), by measuring the MFs from N-body simulations of modified gravity and general relativity. We find strong statistical power when using the MFs to constrain modified theories of gravity: with a galaxy survey that has survey volume $\sim 0.125 (h^{-1}$Gpc$)^3$ and galaxy number density $\sim 1 / (h^{-1}$Mpc$)^{3}$, the two normal-branch DGP models and the F5 $f(R)$ model that we simulated can be discriminated from $\Lambda$CDM at a significance level >~ 5$\sigma$ with an individual MF measurement. Therefore, the MF of large scale structure is potentially a powerful probe of gravity, and its application to real data deserves active explorations.

The turn-around radii of the galaxy groups show the imprint of a long battle between their self-gravitational forces and the accelerating space. The standard $\Lambda$CDM cosmology based on the general relativity (GR) predicts the existence of an upper bound on the expectation value of the turn-around radius which is rarely violated by individual galaxy groups. We speculate that a deviation of the gravitational law from GR on the cosmological scale could cause an appreciable shift of the mean turn-around radius to higher values and make the occurrence of the bound violation more probable. Analyzing the data from high-resolution N-body simulations for two specific models with modified gravity (MG) and the standard GR+$\Lambda$CDM cosmology, we determine the turn-around radii of the massive Rockstar groups from the peculiar motions of the galactic halos located in the bound zone where the fifth force generated by MG is expected to be at most partially shielded. We detect a $4\sigma$ signal of difference in the odds of the bound violations between a fiducial MG and the GR models, proving that the odds of the bound violations increase with the strength of the fifth force produced by the presence of MG. The advantage of using the odds of the bound violations as a complementary diagnostics to probe the nature of gravity is discussed.

In certain theories of modified gravity, solar system constraints on deviations from general relativity (GR) are satisfied by virtue of a so-called screening mechanism, which enables the theory to revert to GR in regions where the matter density is high or the gravitational potential is deep. In the case of chameleon theories, the screening has two contributions -- self-screening, which is due to the mass of an object itself, and environmental screening, which is caused by the surrounding matter -- which are often entangled, with the second contribution being more crucial for less massive objects. A quantitative understanding of the effect of the environment on the screening can prove critical in observational tests of such theories using systems such as the Local Group and dwarf galaxies, for which the environment may be inferred in various ways. We use the high-resolution {\sc liminality} simulation of Shi et al. (2015) to test the fidelity of different definitions of environment. We find that, although the different ways to define environment in practice do not agree with one another perfectly, they can provide useful guidance, and cross checks about how well a dark matter halo is screened. In addition, the screening of subhaloes in dark matter haloes is primarily determined by the environment, with the subhalo mass playing a minor role, which means that lower-resolution simulations where subhaloes are not well resolved can still be useful for understanding the modification of gravity inside subhaloes.

We consider an alternative cold dark matter candidate, ultralight bosons ($m>10^{-22}$eV) described by a complex scalar field (SFDM) with global U(1) symmetry, with comoving particle number density conserved after particle production during standard reheating. We allow for repulsive self-interaction. In a Lambda-SFDM universe, SFDM starts relativistic, evolving from stiff (w=1) to radiation-like (w=1/3), becoming nonrelativistic (w=0) at late times. Thus, a stiff-SFDM-dominated era precedes the familiar radiation-dominated era. SFDM particle mass $m$ and quartic self-interaction strength \lambda, are therefore constrained by cosmological observables, N_{eff}, the effective number of neutrino species during BBN, and z_{eq}, the matter-radiation equality redshift. Since the stochastic gravitational wave background (SGWB) from inflation is amplified during the stiff-SFDM-dominated era, it can also contribute a radiationlike component large enough to affect these observables. Remarkably, this amplification makes this SGWB detectable by current GW experiments, e.g., aLIGO/Virgo and LISA, for Lambda-SFDM models satisfying cosmological constraints, for a range of reheat temperatures T_{re} and currently allowed values of tensor-to-scalar ratio $r$. For given r and $\lambda/(mc^2)^2$, the marginally-allowed Lambda-SFDM model for each T_{re} has the smallest m that satisfies cosmological constraints. For example, for marginally-allowed models with r=0.01 and $\lambda/(mc^2)^2=10^{-18}$eV$^{-1}$cm$^3$, null detection by the aLIGO O1 run excludes 8.75*10^3<T_{re} (GeV)<1.7*10^5 at 95% confidence, demonstrating that GW experiments already place a new kind of cosmological constraint on SFDM. A wider parameter range should be accessible to aLIGO/Virgo O5, with potential to detect this signature of Lambda-SFDM. For this same illustrative family, 3-sigma detection is predicted for 600<T_{re} (GeV)<10^7.

We refine the mass and environment dependent spherical collapse model of chameleon $f(R)$ gravity by calibrating a phenomenological correction inspired by the parameterized post-Friedmann framework against high-resolution $N$-body simulations. We employ our method to predict the corresponding modified halo mass function, and provide fitting formulas to calculate the fractional enhancement of the $f(R)$ halo abundance with respect to that of General Relativity (GR) within a precision of $\lesssim 5\%$ from the results obtained in the simulations. Similar accuracy can be achieved for the full $f(R)$ mass function on the condition that the modeling of the reference GR abundance of halos is accurate at the percent level. We use our fits to forecast constraints on the additional scalar degree of freedom of the theory, finding that upper bounds competitive with current Solar System tests are within reach of cluster number count analyses from ongoing and upcoming surveys at much larger scales. Importantly, the flexibility of our method allows also for this to be applied to other scalar-tensor theories characterized by a mass and environment dependent spherical collapse.

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG.

We study the imprints that theories of gravity beyond GR can leave on the lensing signal around line of sight directions that are predominantly halo-underdense (called troughs) and halo-overdense. To carry out our investigations, we consider the normal branch of DGP gravity, as well as a phenomenological variant thereof that directly modifies the lensing potential. The predictions of these models are obtained with N-body simulation and ray-tracing methods using the ECOSMOG and Ray-Ramses codes. We analyse the stacked lensing convergence profiles around the underdense and overdense lines of sight, which exhibit, respectively, a suppression and a boost w.r.t. the mean in the field of view. The modifications to gravity in these models strengthen the signal w.r.t. $\Lambda{\rm CDM}$ in a scale-independent way. We find that the size of this effect is the same for both underdense and overdense lines of sight, which implies that the density field along the overdense directions on the sky is not sufficiently evolved to trigger the suppression effects of the screening mechanism. These results are robust to variations in the minimum halo mass and redshift ranges used to identify the lines of sight, as well as to different line of sight aperture sizes and criteria for their underdensity and overdensity thresholds.

Solar filaments/prominences are one of the most common features in the corona, which may lead to energetic coronal mass ejections (CMEs) and flares when they erupt. Filaments are about one hundred times cooler and denser than the coronal material, and physical understanding of their material origin remains controversial. Two types of scenarios have been proposed: one argues that the filament plasma is brought into the corona from photosphere or chromosphere through a siphon or evaporation/injection process, while the other suggests that the material condenses from the surrounding coronal plasma due to thermal instability. The elemental abundance analysis is a reasonable clue to constrain the models, as the siphon or evaporation/injection model would predict that the filament material abundances are close to the photospheric or chromospheric ones, while the condensation model should have coronal abundances. In this letter, we analyze the elemental abundances of a magnetic cloud that contains the ejected filament material. The corresponding filament eruption occurred on 1998 April 29, accompanying an M6.8 class soft X-ray flare located at the heliographic coordinates S18E20 (NOAA 08210) and a fast halo CME with the linear velocity of 1374 km s$^{-1}$ near the Sun. We find that the abundance ratios of elements with low and high First Ionization Potential such as Fe/O, Mg/O, and Si/O are 0.150, 0.050, and 0.070, respectively, approaching their corresponding photospheric values 0.065, 0.081, and 0.066, which does not support the coronal origin of the filament plasma.

Connecting in-situ measured solar-wind plasma properties with typical regions on the Sun can provide an effective constraint and test to various solar wind models. We examine the statistical characteristics of the solar wind with an origin in different types of source regions. We find that the speed distribution of coronal hole (CH) wind is bimodal with the slow wind peaking at ~400 km/s and a fast at ~600 km/s. An anti-correlation between the solar wind speeds and the O7+/O6+ ion ratio remains valid in all three types of solar wind as well during the three studied solar cycle activity phases, i.e. solar maximum, decline and minimum. The N_Fe/N_O range and its average values all decrease with the increasing solar wind speed in different types of solar wind. The N_Fe/N_O range (0.06-0.40, FIP bias range 1-7) for AR wind is wider than for CH wind (0.06-0.20, FIP bias range 1-3) while the minimum value of N_Fe/N_O (0.06) does not change with the variation of speed, and it is similar for all sourceregions. The two-peak distribution of CH wind and the anti-correlation between the speed and O7+/O6+ in all three types of solar wind can be explained qualitatively by both the wave-turbulence-driven (WTD) and reconnection-loop-opening (RLO) models, whereas the distribution features of N_Fe/N_O in different source regions of solar wind can be explained more reasonably by the RLO models.

We investigate the information content of various cosmic shear statistics on the theory of gravity. Focusing on the Hu-Sawicki-type $f(R)$ model, we perform a set of ray-tracing simulations and measure the convergence bispectrum, peak counts and Minkowski functionals. We first show that while the convergence power spectrum does have sensitivity to the current value of extra scalar degree of freedom $|f_{\rm R0}|$, it is largely compensated by a change in the present density amplitude parameter $\sigma_{8}$ and the matter density parameter $\Omega_{\rm m0}$. With accurate covariance matrices obtained from 1000 lensing simulations, we then examine the constraining power of the three additional statistics. We find that these probes are indeed helpful to break the parameter degeneracy, which can not be resolved from the power spectrum alone. We show that especially the peak counts and Minkowski functionals have the potential to rigorously (marginally) detect the signature of modified gravity with the parameter $|f_{\rm R0}|$ as small as $10^{-5}$ ($10^{-6}$) if we can properly model them on small ($\sim 1\, \mathrm{arcmin}$) scale in a future survey with a sky coverage of 1,500 squared degrees. We also show that the signal level is similar among the additional three statistics and all of them provide complementary information to the power spectrum. These findings indicate the importance of combining multiple probes beyond the standard power spectrum analysis to detect possible modifications to General Relativity.

The behavior of the axial group speeds of trapped sausage modes plays an important role in determining impulsively generated wave trains, which have often been invoked to account for quasi-periodic signals with quasi-periods of order seconds in a considerable number of coronal structures. We conduct a comprehensive eigenmode analysis, both analytically and numerically, on the dispersive properties of sausage modes in pressureless tubes with three families of continuous radial density profiles. We find a rich variety of the dependence on the axial wavenumber $k$ of the axial group speed $v_{\rm gr}$. Depending on the density contrast and profile steepness as well as on the detailed profile description, the $v_{\rm gr}-k$ curves either possess or do not possess cutoff wavenumbers, and they can behave in either a monotonical or non-monotonical manner. With time-dependent simulations, we further show that this rich variety of the group speed characteristics heavily influences the temporal evolution and Morlet spectra of impulsively generated wave trains. In particular, the Morlet spectra can look substantially different from "crazy tadpoles" found for the much-studied discontinuous density profiles. We conclude that it is necessary to re-examine available high-cadence data to look for the rich set of temporal and spectral features, which can be employed to discriminate between the unknown forms of the density distributions transverse to coronal structures.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

Neutron-skins in coordinate space and proton-skins in momentum space are predicted to coexist in heavy nuclei and their correlation is governed by Liouville's theorem and Heisenberg's uncertainty principle. An analysis of their correlation within a further extended Thomas-Fermi (ETF$^+$) approximation incorporating effects of nucleon short-range correlations reveals generally that protons move faster than neutrons in neutron-skins of heavy nuclei.

Using the liminality N-body simulations of Shi et. al., we present the first predictions for galaxy clustering in $f(R)$ gravity using subhalo abundance matching. We find that, for a given galaxy density, even for an $f(R)$ model with $f_{R0}=-10^{-6}$, for which the cold dark matter clustering is very similar to $\Lambda$CDM, the predicted clustering of galaxies in the $f(R)$ model is very different from $\Lambda$CDM. The deviation can be as large as $40\%$ for samples with mean densities close to that of $L_*$ galaxies. This large deviation is testable given the accuracy that future large-scale galaxy surveys aim to achieve. Our result demonstrates that galaxy surveys can provide a stringent test of General Relativity on cosmological scales, which is comparable to the tests from local astrophysical observations.

We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest" model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $\leq 10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.

The frequency dependence of the longitudinal group speeds of trapped sausage waves plays an important role in determining impulsively generated wave trains, which have often been invoked to account for quasi-periodic signals in coronal loops. We examine how the group speeds ($v_{\rm gr}$) depend on angular frequency ($\omega$) for sausage modes in pressureless coronal tubes with continuous transverse density distributions by solving the dispersion relation pertinent to the case where the density inhomogeneity of arbitrary form takes place in a transition layer of arbitrary thickness. We find that in addition to the transverse lengthscale $l$ and density contrast $\rho_{\rm i}/\rho_{\rm e}$, the group speed behavior depends also on the detailed form of the density inhomogeneity. For parabolic profiles, $v_{\rm gr}$ always decreases with $\omega$ first before increasing again, as happens for the much studied top-hat profiles. For linear profiles, however, the behavior of the $\omega-v_{\rm gr}$ curves is more complex. When $\rho_{\rm i}/\rho_{\rm e} \lesssim 6$, the curves become monotonical for large values of $l$. On the other hand, for higher density contrasts, a local maximum $v_{\rm gr}^{\rm max}$ exists in addition to a local minimum $v_{\rm gr}^{\rm min}$ when coronal tubes are diffuse. With time-dependent computations, we show that the different behavior of group speed curves, the characteristic speeds $v_{\rm gr}^{\rm min}$ and $v_{\rm gr}^{\rm max}$ in particular, is reflected in the temporal evolution and Morlet spectra of impulsively generated wave trains. We conclude that the observed quasi-periodic wave trains not only can be employed to probe such key parameters as density contrasts and profile steepness, but also have the potential to discriminate between the unknown forms of the transverse density distribution.

While standing fast sausage modes in flare loops are often invoked to interpret quasi-periodic pulsations (QPPs) in solar flares, it is unclear as to how they are influenced by the combined effects of a continuous transverse structuring and a finite internal plasma beta ($\beta_{\rm i}$). We derive a generic dispersion relation (DR) governing linear sausage waves in straight magnetic tubes for which plasma pressure is not negligible and the density and temperature inhomogeneities of essentially arbitrary form take place in a layer of arbitrary width. Focusing on fast modes, we find that $\beta_{\rm i}$ only weakly influences $k_{\rm c}$, the critical longitudinal wavenumber separating the leaky from trapped modes. Likewise, for both trapped and leaky modes, the periods $P$ in units of the transverse fast time depend only weakly on $\beta_{\rm i}$, which is compatible with the fact that the effective wavevectors of fast sausage modes are largely perpendicular to the background magnetic field. However, a weak $\beta_{\rm i}$ dependence of the damping times $\tau$ is seen only when the length-to-radius ratio $L/R$ is $\sim 50\%$ larger than some critical value $\pi/(k_{\rm c} R)$, which itself rather sensitively depends on the density contrast, profile steepness as well as on how the transverse structuring is described. In the context of QPPs, we conclude that the much simpler zero-beta theory can be employed for trapped modes, as long as one sees the deduced internal Alfv\'en speed as actually being the fast speed. In contrast, effects due to a finite beta in flare loops should be considered when leaky modes are exploited.

What is nuclear symmetry energy? Why is it important? What do we know about it? Why is it so uncertain especially at high densities? Can the total symmetry energy or its kinetic part be negative? What are the effects of three-body and/or tensor force on symmetry energy? How can we probe the density dependence of nuclear symmetry energy with terrestrial nuclear experiments? What observables of heavy-ion reactions are sensitive to the high-density behavior of nuclear symmetry energy? How does the symmetry energy affect properties of neutron stars, gravitational waves and our understanding about the nature of strong-field gravity? In this lecture, we try to answer these questions as best as we can based on some of our recent work and/or understanding of research done by others. This note summarizes the main points of the lecture.

Fast Radio Bursts (FRBs) are intense radio flashes from the sky that are characterized by millisecond durations and Jansky-level flux densities. We carried out a statistical analysis on FRBs discovered. Their mean dispersion measure, after subtracting the contribution from the interstellar medium of our Galaxy, is found to be $\sim 660\,\rm pc\,cm^{-3}$, supporting their being from cosmological origin. Their energy released in radio band spans about two orders of magnitude, with a mean value of $\sim 10^{39}$ ergs. More interestingly, although the FRB study is still in a very early phase, the published collection of FRBs enables us to derive a useful intensity distribution function. For the 16 non-repeating FRBs detected by Parkes telescope and the Green Bank Telescope, the intensity distribution can be described as $dN/dF_{\rm obs} = (4.1 \pm 1.3) \times 10^3 \, F_{\rm obs}^{-1.1\pm0.2} \; \rm sky^{-1}\,day^{-1}$, where $F_{\rm obs}$ is the observed radio fluence in units of Jy~ms. Here the power-law index is significantly flatter than the expected value of 2.5 for standard candles distributed homogeneously in a flat Euclidean space. Based on this intensity distribution function, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) will be able to detect about 5 FRBs for every 1000 hours of observation time.

Abell 2146 ($z$ = 0.232) consists of two galaxy clusters undergoing a major merger. The system was discovered in previous work, where two large shock fronts were detected using the $\textit{Chandra X-ray Observatory}$, consistent with a merger close to the plane of the sky, caught soon after first core passage. A weak gravitational lensing analysis of the total gravitating mass in the system, using the distorted shapes of distant galaxies seen with ACS-WFC on $\textit{Hubble Space Telescope}$, is presented. The highest peak in the reconstruction of the projected mass is centred on the Brightest Cluster Galaxy (BCG) in Abell 2146-A. The mass associated with Abell 2146-B is more extended. Bootstrapped noise mass reconstructions show the mass peak in Abell 2146-A to be consistently centred on the BCG. Previous work showed that BCG-A appears to lag behind an X-ray cool core; although the peak of the mass reconstruction is centred on the BCG, it is also consistent with the X-ray peak given the resolution of the weak lensing mass map. The best-fit mass model with two components centred on the BCGs yields $M_{200}$ = 1.1$^{+0.3}_{-0.4}$$\times$10$^{15}$M$_{\odot}$ and 3$^{+1}_{-2}$$\times$10$^{14}$M$_{\odot}$ for Abell 2146-A and Abell 2146-B respectively, assuming a mass concentration parameter of $c=3.5$ for each cluster. From the weak lensing analysis, Abell 2146-A is the primary halo component, and the origin of the apparent discrepancy with the X-ray analysis where Abell 2146-B is the primary halo is being assessed using simulations of the merger.

In this talk, we first briefly review the isospin dependence of the total nucleon effective mass $M^{\ast}_{J}$ inferred from analyzing nucleon-nucleus scattering data within an isospin dependent non-relativistic optical potential model, and the isospin dependence of the nucleon E-mass $M^{\ast,\rm{E}}_{J}$ obtained from applying the Migdal-Luttinger theorem to a phenomenological single-nucleon momentum distribution in nuclei constrained by recent electron-nucleus scattering experiments. Combining information about the isospin dependence of both the nucleon total effective mass and E-mass, we then infer the isospin dependence of nucleon k-mass using the well-known relation $M^{\ast}_{J}=M^{\ast,\rm{E}}_{J}\cdot M^{\ast,\rm{k}}_{J}$. Implications of the results on the nucleon mean free path (MFP) in neutron-rich matter are discussed.

Modified gravity theories are a popular alternative to dark energy as a possible explanation for the observed accelerating cosmic expansion, and their cosmological tests are currently an active research field. Studies in recent years have been increasingly focused on testing these theories in the nonlinear regime, which is computationally demanding. Here we show that, under certain circumstances, a whole class of theories can be ruled out by using background cosmology alone. This is possible because certain classes of models (i) are fundamentally incapable of producing specific background expansion histories, and (ii) said histories are incompatible with local gravity tests. As an example, we demonstrate that a popular class of models, $f(R)$ gravity, would not be viable if observations suggest even a slight deviation of the background expansion history from that of the $\Lambda$CDM paradigm.

We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.

Various small-scale structures abound in the solar atmosphere above active regions, playing an important role in the dynamics and evolution therein. We report on a new class of small-scale transition region structures in active regions, characterized by strong emissions but extremely narrow Si IV line profiles as found in observations taken with the Interface Region Imaging Spectrograph (IRIS). Tentatively named as Narrow-line-width UV bursts (NUBs), these structures are located above sunspots and comprise of one or multiple compact bright cores at sub-arcsecond scales. We found six NUBs in two datasets (a raster and a sit-and-stare dataset). Among these, four events are short-living with a duration of $\sim$10 mins while two last for more than 36 mins. All NUBs have Doppler shifts of 15--18 km/s, while the NUB found in sit-and-stare data possesses an additional component at $\sim$50 km/s found only in the C II and Mg II lines. Given that these events are found to play a role in the local dynamics, it is important to further investigate the physical mechanisms that generate these phenomena and their role in the mass transport in sunspots.

The physical connection among and formation mechanisms of various components of the prominence-horn cavity system remain elusive. Here we present observations of such a system, focusing on a section of the prominence that rises and separates gradually from the main body. This forms a configuration sufficiently simple to yield clues to the above issues. It is characterized by embedding horns, oscillations, and a gradual disappearance of the separated material. The prominence-horn structure exhibits a large amplitude longitudinal oscillation with a period of ~150 minutes and an amplitude of ~30 Mm along the trajectory defined by the concave horn structure. The horns also experience a simultaneous transverse oscillation with a much smaller amplitude (~3 Mm) and shorter period (~10-15 minutes), likely representative of a global mode of the large-scale magnetic structure. The gradual disappearance of the structure indicates that the horn, an observational manifestation of the field-aligned transition region separating the cool and dense prominence from the hot and tenuous corona, is formed due to the heating and diluting process of the central prominence mass, while most previous studies suggest that it is the opposite process, i.e., the cooling and condensation of coronal plasmas, to form the horn. This study also demonstrates how the prominence transports magnetic flux to the upper corona, a process essential for the gradual build-up of pre-eruption magnetic energy.

In this Letter, we report the observational constraints on the Hu-Sawicki $f(R)$ theory derived from weak lensing peak abundances, which are closely related to the mass function of massive halos. In comparison with studies using optical or x-ray clusters of galaxies, weak lensing peak analyses have the advantages of not relying on mass-baryonic observable calibrations. With observations from the Canada-France-Hawaii-Telescope Lensing Survey, our peak analyses give rise to a tight constraint on the model parameter $|f_{R0}|$ for $n=1$. The $95\%$ CL limit is $\log_{10}|f_{R0}| < -4.82$ given WMAP9 priors on $(\Omega_{\rm m}, A_{\rm s})$. With Planck15 priors, the corresponding result is $\log_{10}|f_{R0}| < -5.16$.

The recent direct detection of gravitational waves (GWs) from binary black hole mergers (2016, Phys. Rev. Lett. 116, no. 6, 061102; no. 24, 241103) opens up an entirely new non-electromagnetic window into the Universe making it possible to probe physics that has been hidden or dark to electromagnetic observations. In addition to cataclysmic events involving black holes, GWs can be triggered by physical processes and systems involving neutron stars. Properties of neutron stars are largely determined by the equation of state (EOS) of neutron-rich matter, which is the major ingredient in calculating the stellar structure and properties of related phenomena, such as gravitational wave emission from elliptically deformed pulsars and neutron star binaries. Although the EOS of neutron-rich matter is still rather uncertain mainly due to the poorly known density dependence of nuclear symmetry energy at high densities, significant progress has been made recently in constraining the symmetry energy using data from terrestrial nuclear laboratories. These constraints could provide useful information on the limits of GWs expected from neutron stars. Here after briefly reviewing our previous work on constraining gravitational radiation from elliptically deformed pulsars with terrestrial nuclear laboratory data in light of the recent gravitational wave detection, we estimate the maximum gravitational wave strain amplitude, using an optimistic value for the breaking strain of the neutron star crust, for 15 pulsars at distances 0.16 kpc to 0.91 kpc from Earth, and find it to be in the range of $\sim[0.2-31.1]\times 10^{-24}$, depending on the details of the EOS used to compute the neutron star properties. Implications are discussed.

Fast magnetoacoustic wave is an important tool for inferring solar atmospheric parameters. We numerically simulate the propagation of fast wave pulses in randomly structured plasmas mimicking the highly inhomogeneous solar corona. A network of secondary waves is formed by a series of partial reflections and transmissions. These secondary waves exhibit quasi-periodicities in both time and space. Since the temporal and spatial periods are related simply through the fast wave speed, we quantify the properties of secondary waves by examining the dependence of the average temporal period ($\bar{p}$) on the initial pulse width ($w_0$) as well as the density contrast ($\delta_\rho$) and correlation length ($L_c$) that characterize the randomness of the equilibrium density profiles. For small-amplitude pulses, $\delta_\rho$ does not alter $\bar{p}$ significantly. Large-amplitude pulses, on the other hand, enhance the density contrast when $\delta_\rho$ is small but have a smoothing effect when $\delta_\rho$ is sufficiently large. We found that $\bar{p}$ scales linearly with $L_c$ and that the scaling factor is larger for a narrower pulse. However, in terms of the absolute values of $\bar{p}$, broader pulses generate secondary waves with longer periods, and this effect is stronger in random plasmas with shorter correlation lengths. Secondary waves carry the signatures of both the leading wave pulse and background plasma, our study may find applications in MHD seismology by exploiting the secondary waves detected in the dimming regions after CMEs or EUV waves.

We propose a new method to model cluster scaling relations in modified gravity. Using a suite of non-radiative hydrodynamical simulations, we show that the scaling relations of cumulative gas quantities, such as the Sunyaev Zel'dovich effect (Compton-y parameter) and the x-ray Compton-Y parameter, can be accurately predicted using the known results in the $\Lambda$CDM model with a precision of $\sim3\%$. This method provides a reliable way to analyze the gas physics in modified gravity using the less-demanding and much more efficient pure cold dark matter simulations. Our results therefore have important theoretical and practical implications in constraining gravity using cluster surveys.

The shear measurement is a crucial task in the current and the future weak lensing survey projects. And the reconstruction of the point spread function(PSF) is one of the essential steps. In this work, we present three different methods, including Gaussianlets, Moffatlets and EMPCA to quantify their efficiency on PSF reconstruction using four sets of simulated LSST star images. Gaussianlets and Moffatlets are two different sets of basis functions whose profiles are based on Gaussian and Moffat functions respectively. Expectation Maximization(EM) PCA is a statistical method performing iterative procedure to find principal components of an ensemble of star images. Our tests show that: 1) Moffatlets always perform better than Gaussianlets. 2) EMPCA is more compact and flexible, but the noise existing in the Principal Components (PCs) will contaminate the size and ellipticity of PSF while Moffatlets keeps them very well.

We report Shanghai Tian Ma Radio Telescope detections of several long carbon-chain molecules at C and Ku band, including HC3N, HC5N, HC7N, HC9N, C3S, C6H and C8H toward the starless cloud Serpens South 1a. We detected some transitions (HC9N J=13-12 F=12-11 and F=14-13, H13CCCN J=2-1 F=1-0 and F=1-1, HC13CCN J=2-1 F=2-2, F=1-0 and F=1-1, HCC13CN J=2-1 F=1-0 and F=1-1) and resolved some hyperfine components (HC5N J=6-5 F=5-4, H13CCCN J=2-1 F=2-1) for the first time in the interstellar medium. The column densities of these carbon-chain molecules in a range of 10^{12}-10^{13} cm^{-2} are comparable to two carbon-chain molecule rich sources, TMC-1 and Lupus-1A. The abundance ratios are 1.00:(1.11\pm0.15):(1.47\pm0.18) for [H13CCCN]:[HC13CCN]:[HCC13CN]. This result implies that the 13C isotope is also concentrated in the carbon atom adjacent to the nitrogen atom in HC3N in Serpens south 1a, which is similar to TMC-1. The [HC3N]/[H13CCCN] ratio of 78\pm9, the [HC3N]/[HC13CCN] ratio of 70\pm8, and the [HC3N]/[HCC13CN] ratio of 53\pm4 are also comparable to those in TMC-1. In any case, Serpens South 1a proves a testing ground for understanding carbon-chain chemistry.

The well-known Migdal-Luttinger theorem states that the jump of the single-nucleon momentum distribution at the Fermi surface is equal to the inverse of the nucleon effective E-mass. Recent experiments studying short-range correlations (SRC) in nuclei using electron-nucleus scatterings at the Jefferson National Laboratory (JLAB) together with model calculations constrained significantly the Migdal-Luttinger jump at saturation density of nuclear matter. We show that the corresponding nucleon effective E-mass is consequently constrained to $M_0^{\ast,\rm{E}}/M\approx2.22\pm0.35$ in symmetric nuclear matter (SNM) and the E-mass of neutrons is smaller than that of protons in neutron-rich matter. Moreover, the average depletion of the nucleon Fermi sea increases (decreases) approximately linearly with the isospin asymmetry $\delta$ according to $\kappa_{\rm{p}/\rm{n}}\approx 0.21\pm0.06 \pm (0.19\pm0.08)\delta$ for protons (neutrons). These results will help improve our knowledge about the space-time non-locality of the single-nucleon potential in neutron-rich nucleonic matter useful in both nuclear physics and astrophysics.

We analyze the kinematics of six Virgo cluster dwarf early-type galaxies (dEs) from their globular cluster (GC) systems. We present new Keck/DEIMOS spectroscopy for three of them and reanalyze the data found in the literature for the remaining three. We use two independent methods to estimate the rotation amplitude (Vmax) and velocity dispersion (sigma_GC) of the GC systems and evaluate their statistical significance by simulating non-rotating GC systems with the same number of GC satellites and velocity uncertainties. Our measured kinematics agree with the published values for the three galaxies from the literature and, in all cases, some rotation is measured. However, our simulations show that the null hypothesis of being non-rotating GC systems cannot be ruled out. In the case of VCC1861, the measured Vmax and the simulations indicate that it is not rotating. In the case of VCC1528, the null hypothesis can be marginally ruled out, thus, it might be rotating although further confirmation is needed. In our analysis, we find that, in general, the measured Vmax tends to be overestimated and the measured sigma_GC tends to be underestimated by amounts that depend on the intrinsic Vmax/sigma_GC, the number of observed GCs (N_GC), and the velocity uncertainties. The bias is negligible when N_GC>~20. In those cases where a large N_GC is not available, it is imperative to obtain data with small velocity uncertainties. For instance, errors of <2km/s lead to Vmax<10km/s for a system that is intrinsically not rotating.

Modified theories of gravity provide us with a unique opportunity to generate innovative tests of gravity. In Chameleon f(R) gravity, the gravitational potential differs from the weak-field limit of general relativity (GR) in a mass dependent way. We develop a probe of gravity which compares high mass clusters, where Chameleon effects are weak, to low mass clusters, where the effects can be strong. We utilize the escape velocity edges in the radius/velocity phase space to infer the gravitational potential profiles on scales of 0.3-1 virial radii. We show that the escape edges of low mass clusters are enhanced compared to GR, where the magnitude of the difference depends on the background field value |fR0|. We validate our probe using N-body simulations and simulated light cone galaxy data. For a DESI (Dark Energy Spectroscopic Instrument) Bright Galaxy Sample, including observational systematics, projection effects, and cosmic variance, our test can differentiate between GR and Chameleon f(R) gravity models, |fR0| = 4e-6 (2e-6) at > 5{\sigma} (> 2{\sigma}), more than an order of magnitude better than current cluster-scale constraints.

In the solar wind, magnetic field power spectra usually show several power-laws. In this paper, magnetic field data from the Cluster mission during an undisturbed interval of slow solar wind is analyzed at 0.28Hz, near the spectral break point between the ion inertial and dissipation/dispersion ranges. Assuming Taylor's frozen-in condition, it corresponds to a proton kinetic scale of $kv_A/\Omega_p \sim 0.38$, where $v_A$ and $\Omega_p$ are the Alfv\'en speed and proton angular gyrofrequency, respectively. Data show that the Cluster spacecraft passed through a series of wavepackets. A strong isolated wavepacket is found to be in accordance with the four Cluster satellites crossing an Alfv\'en vortex, a nonlinear solution to the incompressible MHD equations. A strong agreement is seen between the data from four satellites and a model vortex with a radius of the order of $40$ times the local proton gyro-radii. The polarization at different spacecraft is compared and is found to agree with the vortex model, whereas it cannot be explained solely by the linear plane wave approach.

We study intensity disturbances above a solar polar coronal hole seen in the AIA 171 \AA\ and 193 \AA\ passbands, aiming to provide more insights into their physical nature. The damping and power spectra of the intensity disturbances with frequencies from 0.07 mHz to 10.5 mHz are investigated. The damping of the intensity disturbances tends to be stronger at lower frequencies, and their damping behavior below 980" (for comparison, the limb is at 945") is different from what happens above. No significant difference is found between the damping of the intensity disturbances in the AIA 171 \AA\ and that in the AIA 193 \AA. The indices of the power spectra of the intensity disturbances are found to be slightly smaller in the AIA 171 \AA\ than in the AIA 193 \AA, but the difference is within one sigma deviation. An additional enhanced component is present in the power spectra in a period range of 8--40 minutes at lower heights. While the power spectra of spicule is highly correlated with its associated intensity disturbance, it suggests that the power spectra of the intensity disturbances might be a mixture of spicules and wave activities. We suggest that each intensity disturbance in the polar coronal hole is possibly a series of independent slow magnetoacoustic waves triggered by spicular activities.

Significant progress has been made recently in constraining the isospin-dependent parameters characterizing the SRC (short-range correlation)-modified single-nucleon momentum distribution in neutron-rich nucleonic matter using both experimental data and microscopic model calculations. Using the constrained single-nucleon momentum distribution in a nonlinear relativistic mean field (RMF) model, we study the equation of state (EOS) of asymmetric nucleonic matter (ANM), especially the density dependence of nuclear symmetry energy $E_{\rm{sym}}(\rho)$. Firstly, as a test of the model, the average nucleon kinetic energy extracted recently from electron-nucleus scattering experiments using a neutron-proton dominance model is well reproduced by the RMF model incorporating effects of the SRC-induced high momentum nucleons, while it is significantly under predicted by the RMF model using a step function for the single-nucleon momentum distribution as in free Fermi gas (FFG) models. Secondly, the kinetic symmetry energy of quasi-nucleons is found to be $E^{\rm{kin}}_{\rm{sym}}(\rho_0)=-16.94\pm13.66\,\rm{MeV}$ which is dramatically different from the prediction of $E^{\rm{kin}}_{\rm{sym}}(\rho_0)\approx 12.5$ MeV by FFG models at nuclear matter saturation density $\rho_0=0.16\,\rm{fm}^{-3}$. Thirdly, comparing the RMF calculations with and without the high momentum nucleons using two sets of model parameters both reproducing identically all empirically constraints on the EOS of symmetric nuclear matter (SNM) and the symmetry energy of ANM at $\rho_0$, the SRC-modified single-nucleon momentum distribution is found to make the $E_{\rm{sym}}(\rho)$ more concave around $\rho_0$ by softening it significantly at both sub-saturation and supra-saturation densities, leading to an isospin-dependent incompressibility of ANM in better agreement with existing experimental data.

Standing oscillations with multiple periods were found in a number of atmospheric structures on the Sun. The ratio of the period of the fundamental to twice the one of its first overtone, $P_1/2P_2$, is important in applications of solar magneto-seismology. We examine how field-aligned flows impact $P_1/2P_2$ of standing modes in solar magnetic cylinders. For coronal loops, the flow effects are significant for both fast kink and sausage modes. For kink ones, they reduce $P_1/2P_2$ by up to 17\% relative to the static case even when the density contrast between the loop and its surroundings approaches infinity. For sausage modes, the reduction in $P_1/2P_2$ due to flow is typically $\lesssim 5.5\%$ compared with the static case. However, the threshold aspect ratio, only above which can trapped sausage modes be supported, may increase dramatically with the flow magnitude. For photospheric tubes, the flow effect on $P_1/2P_2$ is not as strong. However, when applied to sausage modes, introducing field-aligned flows offers more possibilities in interpreting the multiple periods recently measured. We conclude that field-aligned flows should be taken into account to help better understand what causes the departure of $P_1/2P_2$ from unity.

Standing fast sausage modes in flare loops were suggested to account for a considerable number of quasi-periodic pulsations (QPPs) in the light curves of solar flares. This study continues our investigation into the possibility to invert the measured periods $P$ and damping times $\tau$ of sausage modes to deduce the transverse Alfv\'en time $R/v_{\rm Ai}$, density contrast $\rho_{\rm i}/\rho_{\rm e}$, and the steepness of the density distribution transverse to flare loops. A generic dispersion relation (DR) governing linear sausage modes is derived for pressureless cylinders where density inhomogeneity of arbitrary form takes place within the cylinder. We show that in general the inversion problem is under-determined for QPP events where only a single sausage mode exists, be the measurements spatially resolved or unresolved. While $R/v_{\rm Ai}$ can be inferred to some extent, the range of possible steepness parameters may be too broad to be useful. However, for spatially resolved measurements where an additional mode is present, it is possible to deduce self-consistently $\rho_{\rm i}/\rho_{\rm e}$, the profile steepness, and the internal Alfv\'en speed $v_{\rm Ai}$. We show that at least for a recent QPP event that involves a fundamental kink mode in addition to a sausage one, flare loop parameters are well constrained, even if the specific form of the transverse density distribution remains unknown. We conclude that spatially resolved, multi-mode QPP measurements need to be pursued for inferring flare loop parameters.

Coronal bright points (BPs) are associated with magnetic bipolar features (MBFs) and magnetic cancellation. Here, we investigate how BP-associated MBFs form and how the consequent magnetic cancellation occurs. We analyse longitudinal magnetograms from the Helioseismic and Magnetic Imager to investigate the photospheric magnetic flux evolution of 70 BPs. From images taken in the 193 A passband of the Atmospheric Imaging Assembly (AIA) we dermine that the BPs' lifetimes vary from 2.7 to 58.8 hours. The formation of the BP MBFs is found to involve three processes, namely emergence, convergence and local coalescence of the magnetic fluxes. The formation of a MBF can involve more than one of these processes. Out of the 70 cases, flux emergence is the main process of a MBF buildup of 52 BPs, mainly convergence is seen in 28, and 14 cases are associated with local coalescence. For MBFs formed by bipolar emergence, the time difference between the flux emergence and the BP appearance in the AIA 193 \AA\ passband varies from 0.1 to 3.2 hours with an average of 1.3 hours. While magnetic cancellation is found in all 70 BPs, it can occur in three different ways: (I) between a MBF and small weak magnetic features (in 33 BPs); (II) within a MBF with the two polarities moving towards each other from a large distance (34 BPs); (III) within a MBF whose two main polarities emerge in the same place simultaneously (3 BPs). While a MBF builds up the skeleton of a BP, we find that the magnetic activities responsible for the BP heating may involve small weak fields.

We introduce and demonstrate the power of a method to speed up current iterative techniques for N-body modified gravity simulations. Our method is based on the observation that the accuracy of the final result is not compromised if the calculation of the fifth force becomes less accurate, but substantially faster, in high-density regions where it is weak due to screening. We focus on the nDGP model which employs Vainshtein screening, and test our method by running AMR simulations in which the solutions on the finer levels of the mesh (high density) are not obtained iteratively, but instead interpolated from coarser levels. We show that the impact this has on the matter power spectrum is below $1\%$ for $k < 5h/{\rm Mpc}$ at $z = 0$, and even smaller at higher redshift. The impact on halo properties is also small ($\lesssim 3\%$ for abundance, profiles, mass; and $\lesssim 0.05\%$ for positions and velocities). The method can boost the performance of modified gravity simulations by more than a factor of 10, which allows them to be pushed to resolution levels that were previously hard to achieve.

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field.

The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq \rho$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inflation. Also, we study how the associated enhanced expansion rate due to the stiff phase might facilitate a first-order electroweak symmetry breaking phase transition, in light of the recently measured value of the Higgs boson mass.

Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $\Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150\% while the deviation of the conventional matter power spectrum is less than 50\%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $\Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies. In addition, based on our findings, we also discuss several possible methods of making use of the differences between the conventional and effective dark matter power spectra in $f(R)$ gravity to discriminate the theory from the $\Lambda$CDM model.

We examine the influence of a continuous density structuring transverse to coronal slabs on the dispersive properties of fundamental standing kink and sausage modes supported therein. We derive generic dispersion relations (DRs) governing linear fast waves in pressureless straight slabs with general transverse density distributions, and focus on the cases where the density inhomogeneity takes place in a layer of arbitrary width and in arbitrary form. The physical relevance of the solutions to the DRs is demonstrated by the corresponding time-dependent computations. For all profiles examined, the lowest-order kink modes are trapped regardless of longitudinal wavenumber $k$. A continuous density distribution introduces a difference to their periods of $\lesssim 13\%$ when $k$ is the observed range, relative to the case where the density profile takes a step-function form. Sausage modes and other branches of kink modes are leaky at small $k$, and their periods and damping times are heavily influenced by how the transverse density profile is prescribed, the lengthscale in particular. These modes have sufficiently high quality to be observable only for physical parameters representative of flare loops. We conclude that while the simpler DR pertinent to a step-function profile can be used for the lowest-order kink modes, the detailed information on the transverse density structuring needs to be incorporated into studies of sausage modes and higher-order kink modes.

A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the scalar-tensor subset of alternative gravity models. We focus on effects of the slope $L$ of nuclear symmetry energy at saturation density and the high-density behavior of nuclear symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of nuclear symmetry energy leads to large changes in the binding energy of neutron stars. The difference in predictions using the GR and the scalar-tensor theory appears only for massive neutron stars, and even then is significantly smaller than the difference resulting from variations in the symmetry energy.

In this paper, we report radio observations of the Galactic Center magnetar PSR J1745-2900 at six epochs between June and October, 2014. These observations were carried out using the new Shanghai Tian Ma Radio Telescope at a frequency of 8.6 GHz. Both the flux density and integrated profile of PSR J1745-2900 show dramatic changes from epoch to epoch showing that the pulsar was in its "erratic" phase. On MJD 56836, the flux density of this magnetar was about 8.7 mJy, which was ten times large than that reported at the time of discovery, enabling a single-pulse analysis. The emission is dominated by narrow "spiky" pulses which follow a log-normal distribution in peak flux density. From 1913 pulses, we detected 53 pulses whose peak flux density is ten times greater than that of the integrated profile. They are concentrated in pulse phase at the peaks of the integrated profile. The pulse widths at the 50% level of these bright pulses was between 0.2 to 0.9 deg, much narrower than that of integrated profile (~12 deg). The observed pulse widths may be limited by interstellar scattering. No clear correlation was found between the widths and peak flux density of these pulses and no evidence was found for subpulse drifting. Relatively strong spiky pulses are also detected in the other five epochs of observation, showing the same properties as that detected in MJD 56836. These strong spiky pulses cannot be classified as "giant" pulses but are more closely related to normal pulse emission.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

We explore the Minkowski functionals of weak lensing convergence map to distinguish between $f(R)$ gravity and the general relativity (GR). The mock weak lensing convergence maps are constructed with a set of high-resolution simulations assuming different gravity models. It is shown that the lensing MFs of $f(R)$ gravity can be considerably different from that of GR because of the environmentally dependent enhancement of structure formation. We also investigate the effect of lensing noise on our results, and find that it is likely to distinguish F5, F6 and GR gravity models with a galaxy survey of $\sim3000$ degree$^2$ and with a background source number density of $n_g=30~{\rm arcmin}^{-2}$, comparable to an upcoming survey dark energy survey (DES). We also find that the $f(R)$ signal can be partially degenerate with the effect of changing cosmology, but combined use of other observations, such as the cosmic microwave background (CMB) data, can help break this degeneracy.

We study lensing by voids in Cubic Galileon and Nonlocal gravity cosmologies, which are examples of theories of gravity that modify the lensing potential. We find voids in the dark matter and halo density fields of N-body simulations and compute their lensing signal analytically from the void density profiles, which we show are well fit by a simple analytical formula. In the Cubic Galileon model, the modifications to gravity inside voids are not screened and they approximately double the size of the lensing effects compared to GR. The difference is largely determined by the direct effects of the fifth force on lensing and less so by the modified density profiles. For this model, we also discuss the subtle impact on the force and lensing calculations caused by the screening effects of haloes that exist in and around voids. In the Nonlocal model, the impact of the modified density profiles and the direct modifications to lensing are comparable, but they boost the lensing signal by only $\approx 10\%$, compared with that of GR. Overall, our results suggest that lensing by voids is a promising tool to test models of gravity that modify lensing.

At present, neutron star radii from both observations and model predictions remain very uncertain. Whereas different models can predict a wide range of neutron star radii, it is not possible for most models to predict radii that are smaller than about 10 km, thus if such small radii are established in the future they will be very difficult to reconcile with model estimates. By invoking a new term in the equation of state that enhances the energy density, but leaves the pressure unchanged we simulate the current uncertainty in the neutron star radii. This new term can be possibly due to the exchange of the weakly interacting light U-boson with appropriate in-medium parameters, which does not compromise the success of the conventional nuclear models. The validity of this new scheme will be tested eventually by more precise measurements of neutron star radii.

Standing sausage modes in flare loops are important for interpreting quasi-periodic pulsations (QPPs) in solar flare lightcurves. We propose an inversion scheme that consistently uses their periods $P$ and damping times $\tau$ to diagnose flare loop parameters. We derive a generic dispersion relation governing linear sausage waves in pressure-less straight tubes, for which the transverse density inhomogeneity takes place in a layer of arbitrary width $l$ and is of arbitrary form. We find that $P$ and $\tau$ depend on the combination of $[R/v_{\rm Ai}, L/R, l/R, \rho_{\rm i}/\rho_{\rm e}]$, where $R$ is the loop radius, $L$ is the looplength, $v_{\rm Ai}$ is the internal Alfv\'en speed, and $\rho_{\rm i}/\rho_{\rm e}$ is the density contrast. For all the density profiles examined, $P$ and $\tau$ experience saturation when $L/R \gg 1$, yielding an inversion curve in the $[R/v_{\rm Ai}, l/R, \rho_{\rm i}/\rho_{\rm e}]$ space with a specific density profile when $L/R$ is sufficiently large. When applied to a spatially unresolved QPP event, the scheme yields that $R/v_{\rm Ai}$ is the best constrained, whereas $l/R$ corresponds to the other extreme. For spatially resolved QPPs, while $L/R \gg 1$ cannot be assumed beforehand, an inversion curve remains possible due to additional geometrical constraints. When a spatially resolved QPP event involves another mode, as is the case for a recent event, the full set of $[v_{\rm Ai}, l, \rho_{\rm i}/\rho_{\rm e}]$ can be inferred. We conclude that the proposed scheme provides a useful tool for magneto-seismologically exploiting QPPs.

Effects of the completely unknown symmetry (isovector) potential of the \D on the total and differential \rpi in heavy-ion collisions at beam energies from 100 to 1000 MeV/A are explored within an isospin-dependent transport model IBUU. The effects are found to be negligible at beam energies above the pion production threshold due to the very short lifetimes of less than 2 fm/c for $\Delta$ resonances with masses around $m_{\Delta}=1232$ MeV, leaving the $\pi^-/\pi^+$ ratios of especially the energetic pions still a reliable probe of the high-density behavior of nuclear symmetry energy $E_{sym}(\rho)$. However, as the beam energy becomes deeply sub-threshold for pion production, effects of the $\Delta$ symmetry potential becomes appreciable especially on the \rpi of low-energy pions from the decays of low-mass $\Delta$ resonances which have lived long enough to be affected by their mean-field potentials, providing a useful tool to study the symmetry potential and spectroscopy of $\Delta$ resonances in neutron-rich nuclear matter. Interestingly though, even at the deeply sub-threshold beam energies, the differential \rpi of energetic pions remains sensitive to the \esym at supra saturation densities with little influence from the uncertain symmetry potential of the $\Delta$ resonance.

We propose a new cosmological test of gravity, by using the observed mass fraction of X-ray emitting gas in massive galaxy clusters. The cluster gas fraction, believed to be a fair sample of the average baryon fraction in the Universe, is a well-understood observable, which has previously mainly been used to constrain background cosmology. In some modified gravity models, such as $f(R)$ gravity, gas temperature in a massive cluster is determined by the effective mass of that cluster, which can be larger than its true mass. On the other hand, X-ray luminosity is determined by the true gas density, which in both modified gravity and $\Lambda$CDM models depends mainly on $\Omega_{\rm b}/\Omega_{\rm m}$ and hence the true total cluster mass. As a result, the standard practice of combining gas temperatures and X-ray surface brightnesses of clusters to infer their gas fractions can, in modified gravity models, lead to a larger - in $f(R)$ gravity this can be $1/3$ larger - value of $\Omega_{\rm b}/\Omega_{\rm m}$ than that inferred from other observations such as the CMB. A quick calculation shows that the Hu-Sawicki $n=1$ $f(R)$ model with $|\bar{f}_{R0}|=3\sim5\times10^{-5}$ is in tension with the gas fraction data of the 42 clusters analysed by Allen et al. (2008). We also discuss the implications for other modified gravity models.

We develop and test an algorithm to rescale a simulated dark-matter particle distribution or halo catalogue from a standard gravity model to that of a modified gravity model. This method is based on that of Angulo & White but with some additional ingredients to account for (i) scale-dependent growth of linear density perturbations and (ii) screening mechanisms that are generic features of viable modified gravity models. We attempt to keep the method as general as possible, so that it may plausibly be applied to a wide range of modified theories, although tests against simulations are restricted to a subclass of $f(R)$ models at this stage. We show that rescaling allows the power spectrum of matter to be reproduced at the $\sim 3$ per cent level in both real and redshift space up to $k=0.1h\,\mathrm{Mpc}^{-1}$ if we change the box size and alter the particle displacement field; this limit can be extended to $k=1h\,\mathrm{Mpc}^{-1}$ if we additionally alter halo internal structure. We simultaneously develop an algorithm that can be applied directly to a halo catalogue, in which case the halo mass function and clustering can be reproduced at the $\sim 5$ per cent level. Finally we investigate the clustering of halo particle distributions, generated from rescaled halo catalogues, and find that a similar accuracy can be reached.

Quasi-periodic propagating disturbances (PDs) are ubiquitous in polar coronal holes on the Sun. It remains unclear as to what generates PDs. In this work, we investigate how the PDs are generated in the solar atmosphere by analyzing a fourhour dataset taken by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We find convincing evidence that spicular activities in the solar transition region as seen in the AIA 304 {\AA} passband are responsible for PDs in the corona as revealed in the AIA 171 {\AA} images. We conclude that spicules are an important source that triggers coronal PDs.

We introduce the idea of {\it effective} dark matter halo catalog in $f(R)$ gravity, which is built using the {\it effective} density field. Using a suite of high resolution N-body simulations, we find that the dynamical properties of halos, such as the distribution of density, velocity dispersion, specific angular momentum and spin, in the effective catalog of $f(R)$ gravity closely mimic those in the $\Lambda$CDM model. Thus, when using effective halos, an $f(R)$ model can be viewed as a $\Lambda$CDM model. This effective catalog therefore provides a convenient way for studying the baryonic physics, the galaxy halo occupation distribution and even semi-analytical galaxy formation in $f(R)$ cosmologies.

Sausage modes are important in coronal seismology. Spatially damped propagating sausage waves were recently observed in the solar atmosphere. We examine how wave leakage influences the spatial damping of sausage waves propagating along coronal structures modeled by a cylindrical density enhancement embedded in a uniform magnetic field. Working in the framework of cold magnetohydrodynamics, we solve the dispersion relation (DR) governing sausage waves for complex-valued longitudinal wavenumber $k$ at given real angular frequencies $\omega$. For validation purposes, we also provide analytical approximations to the DR in the low-frequency limit and in the vicinity of $\omega_{\rm c}$, the critical angular frequency separating trapped from leaky waves. In contrast to the standing case, propagating sausage waves are allowed for $\omega$ much lower than $\omega_{\rm c}$. However, while able to direct their energy upwards, these low-frequency waves are subject to substantial spatial attenuation. The spatial damping length shows little dependence on the density contrast between the cylinder and its surroundings, and depends only weakly on frequency. This spatial damping length is of the order of the cylinder radius for $\omega \lesssim 1.5 v_{\rm Ai}/a$, where $a$ and $v_{\rm Ai}$ are the cylinder radius and the Alfv\'en speed in the cylinder, respectively. We conclude that if a coronal cylinder is perturbed by symmetric boundary drivers (e.g., granular motions) with a broadband spectrum, wave leakage efficiently filters out the low-frequency components.

We show that Solar System tests can place very strong constraints on K-mouflage models of gravity, which are coupled scalar field models with nontrivial kinetic terms that screen the fifth force in regions of large gravitational acceleration. In particular, the bounds on the anomalous perihelion of the Moon imposes stringent restrictions on the K-mouflage Lagrangian density, which can be met when the contributions of higher-order operators in the static regime are sufficiently small. The bound on the rate of change of the gravitational strength in the Solar System constrains the coupling strength $\beta$ to be smaller than $0.1$. These two bounds impose tighter constraints than the results from the Cassini satellite and Big Bang Nucleosynthesis. Despite the Solar System restrictions, we show that it is possible to construct viable models with interesting cosmological predictions. In particular, relative to $\Lambda$-CDM, such models predict percent-level deviations for the clustering of matter and the number density of dark matter haloes. This makes these models predictive and testable by forthcoming observational missions.

We report on the first Interface Region Imaging Spectrograph (IRIS) study of cool transition region loops. This class of loops has received little attention in the literature. A cluster of such loops was observed on the solar disk in active region NOAA11934, in the Si IV 1402.8 \AA\ spectral raster and 1400 \AA\ slit-jaw (SJ) images. We divide the loops into three groups and study their dynamics and interaction. The first group comprises relatively stable loops, with 382--626\,km cross-sections. Observed Doppler velocities are suggestive of siphon flows, gradually changing from -10 km/s at one end to 20 km/s at the other end of the loops. Nonthermal velocities from 15 to 25 km/s were determined. These physical properties suggest that these loops are impulsively heated by magnetic reconnection occurring at the blue-shifted footpoints where magnetic cancellation with a rate of $10^{15}$ Mx/s is found. The released magnetic energy is redistributed by the siphon flows. The second group corresponds to two footpoints rooted in mixed-magnetic-polarity regions, where magnetic cancellation occurred at a rate of $10^{15}$ Mx/s and line profiles with enhanced wings of up to 200 km/s were observed. These are suggestive of explosive-like events. The Doppler velocities combined with the SJ images suggest possible anti-parallel flows in finer loop strands. In the third group, interaction between two cool loop systems is observed. Evidence for magnetic reconnection between the two loop systems is reflected in the line profiles of explosive events, and a magnetic cancellation rate of $3\times10^{15}$ Mx/s observed in the corresponding area. The IRIS observations have thus opened a new window of opportunity for in-depth investigations of cool transition region loops. Further numerical experiments are crucial for understanding their physics and their role in the coronal heating processes.

The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1} M_\odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2 \times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.

We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.

Standing sausage modes are important in interpreting quasi-periodic pulsations in the lightcurves of solar flares. Their periods and damping times play an important role in seismologically diagnosing key parameters like the magnetic field strength in regions where flare energy is released. Usually such applications are based on theoretical results neglecting unresolved fine structures in magnetized loops. However, the existence of fine structuring is suggested on both theoretical and observational grounds. Adopting the framework of cold magnetohydrodynamics (MHD), we model coronal loops as magnetized cylinders with a transverse equilibrium density profile comprising a monolithic part and a modulation due to fine structuring in the form of concentric shells. The equation governing the transverse velocity perturbation is solved with an initial-value-problem approach, and the effects of fine structuring on the periods $P$ and damping times $\tau$ of global, leaky, standing sausage modes are examined. A parameter study shows that fine structuring, be it periodically or randomly distributed, brings changes of only a few percent to $P$ and $\tau$ when there are more than about ten shells. The monolithic part, its steepness in particular, plays a far more important role in determining $P$ and $\tau$. We conclude that when measured values of $P$ and $\tau$ of sausage modes are used for seismological purposes, it is justified to use theoretical results where the effects due to fine structuring are neglected.

We use properties of void populations identified in $N$-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \lcdm~cosmology. We analyze multiple simulation realizations, which were designed to mimic the expected number densities, volumes, and redshifts of the upcoming Euclid satellite and a lower-redshift ground-based counterpart survey, using the public {\tt VIDE} toolkit. We examine void abundances, ellipicities, radial density profiles, and radial velocity profiles at redshifts 1.0 and 0.43. We find that stronger f(R) coupling strengths eliminates small voids and produces voids up to $\sim 20\%$ larger in radius, leading to a significant tilt in the void number function. Additionally, under the influence of modified gravity, voids at all scales tend to be measurably emptier with correspondingly higher compensation walls. The velocity profiles reflect this, showing increased outflows inside voids and increased inflows outside voids. Using the void number function as an example, we forecast that future surveys can constrain the modified gravity coupling strength to $\sim 3 \times 10^{-5}$ using voids.

The critical densities and impact of forming \D resonances in neutron stars are investigated within an extended nonlinear relativistic mean-field (RMF) model. The critical densities for the formation of four different charge states of \D are found to depend differently on the separate kinetic and potential parts of nuclear symmetry energy, the first example of a microphysical property of neutron stars to do so. Moreover, they are sensitive to the in-medium Delta mass $m_{\Delta}$ and the completely unknown $\Delta$-$\rho$ coupling strength $g_{\rho\Delta}$. In the universal baryon-meson coupling scheme where the respective $\Delta$-meson and nucleon-meson coupling constants are assumed to be the same, the critical density for the first $\Delta^-(1232)$ to appear is found to be \rc=$(2.08\pm0.02)\rho_0$ using RMF model parameters consistent with current constraints on all seven macroscopic parameters usually used to characterize the equation of state (EoS) of isospin-asymmetric nuclear matter (ANM) at saturation density $\rho_0$. Moreover, the composition and the mass-radius relation of neutron stars are found to depend significantly on the values of the $g_{\rho\Delta}$ and $m_{\Delta}$.

The energy of a free gas of neutrons and protons is well known to be approximately isospin parabolic with a negligibly small quartic term of only $0.45$ MeV at the saturation density of nuclear matter $\rho_0=0.16/\rm{fm}^3$. Using an isospin-dependent single-nucleon momentum distribution including a high (low) momentum tail (depletion) with its shape parameters constrained by recent high-energy electron scattering and medium-energy nuclear photodisintegration experiments as well as the state-of-the-art calculations of the deuteron wave function and the equation of state of pure neutron matter near the unitary limit within several modern microscopic many-body theories, we show for the first time that the kinetic energy of interacting nucleons in neutron-rich nucleonic matter has a significant quartic term of $7.18\pm2.52\,\rm{MeV}$. Such a large quartic term has significant ramifications in determining the equation of state of neutron-rich nucleonic matter using both terrestrial and astrophysical observables.

We show that, in the presence of massive neutrinos, the Galileon gravity model provides a very good fit to the current CMB temperature, CMB lensing and BAO data. This model, which we dub ${\nu} \rm{Galileon}$, when assuming its stable attractor background solution, contains the same set of free parameters as $\Lambda\rm{CDM}$, although it leads to different expansion dynamics and nontrivial gravitational interactions. The data provide compelling evidence ($\gtrsim 6\sigma$) for nonzero neutrino masses, with $\Sigma m_\nu \gtrsim 0.4\ {\rm eV}$ at the $2\sigma$ level. Upcoming precision terrestrial measurements of the absolute neutrino mass scale therefore have the potential to test this model. We show that CMB lensing measurements at multipoles $l \lesssim 40$ will be able to discriminate between the ${\nu} \rm{Galileon}$ and $\Lambda\rm{CDM}$ models. Unlike $\Lambda\rm{CDM}$, the ${\nu} \rm{Galileon}$ model is consistent with local determinations of the Hubble parameter. The presence of massive neutrinos lowers the value of $\sigma_8$ substantially, despite of the enhanced gravitational strength on large scales. Unlike $\Lambda\rm{CDM}$, the ${\nu} \rm{Galileon}$ model predicts a negative ISW effect, which is difficult to reconcile with current observational limits.

In order to assess the ability of purely crust-driven glitch models to match the observed glitch activity in the Vela pulsar, we conduct a systematic analysis of the dependence of the fractional moment of inertia of the inner crustal neutrons on the stiffness of the nuclear symmetry energy at saturation density $L$. We take into account both crustal entrainment and the fact that only a fraction $Y_{\rm g}$ of the core neutrons may couple to the crust on the glitch-rise timescale. We use a set of consistently-generated crust and core compositions and equations-of-state which are fit to results of low-density pure neutron matter calculations. When entrainment is included at the level suggested by recent microscopic calculations and the core is fully coupled to the crust, the model is only able to account for the Vela glitch activity for a 1.4$M_{\odot}$ star if the equation of state is particularly stiff $L>100$ MeV. However, an uncertainty of about 10\% in the crust-core transition density and pressure allows for the Vela glitch activity to be marginally accounted for in the range $L\approx30-60$MeV consistent with a range of experimental results. Alternatively, only a small amount of core neutrons need be involved. If less than 50\% of the core neutrons are coupled to the crust during the glitch, we can also account for the Vela glitch activity using crustal neutrons alone for EOSs consistent with the inferred range of $L$. We also explore the possibility of Vela being a high-mass neutron star, and of crustal entrainment being reduced or enhanced relative to its currently predicted values.

A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy $L$ of such models, and constraints extracted on $L$ from confronting them with observations. We focus on six sets of observations and proposed explanations: (i) The cooling rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear pasta regions of the inner crust, (ii) the upper limit of the observed periods of young X-ray pulsars, confronting models of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer, (iii) glitches from the Vela pulsar, confronting the paradigm that they arise due to a sudden re-coupling of the crustal neutron superfluid to the crustal lattice after a period during which they were decoupled due to vortex pinning, (iv) The frequencies of quasi-periodic oscillations in the X-ray tail of light curves from giant flares from soft gamma-ray repeaters, confronting models of torsional crust oscillations, (v) the upper limit on the frequency to which millisecond pulsars can be spun-up due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold frequency determined in part by the viscous dissipation timescale at the crust-core boundary, and (vi) the observations of precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star just before merger.

We present a kinematic analysis of the globular cluster systems and diffuse stellar light of four intermediate luminosity (sub-$L^{\ast}$) early-type galaxies in the Virgo cluster based on Gemini/GMOS data. Our galaxy sample is fainter ($-23.8<M_K<-22.7$) than most previous studies, nearly doubling the number of galaxies in this magnitude range that now have GC kinematics. The data for the diffuse light extends to $4R_e$, and the data for the globular clusters reaches 8--$12R_e$. We find that the kinematics in these outer regions are all different despite the fact that these four galaxies have similar photometric properties, and are uniformly classified as "fast rotators" from their stellar kinematics within $1R_e$. The globular cluster systems exhibit a wide range of kinematic morphology. The rotation axis and amplitude can change between the inner and outer regions, including a case of counter-rotation. This difference shows the importance of wide-field kinematic studies, and shows that stellar and GC kinematics can change significantly as one moves beyond the inner regions of galaxies. Moreover, the kinematics of the globular cluster systems can differ from that of the stars, suggesting that the formation of the two populations are also distinct.

We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the clusters. We discuss a number of important issues that need to be taken into account, associated with the use of nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, the same as in $\Lambda$CDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the cluster radius. However, at distances $\sim \left[2-20\right] {\rm Mpc}/h$ from the cluster center, we find that the surrounding force profiles are enhanced by $\sim20-40\%$ in the Cubic Galileon model. This has an impact on dynamical mass estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can also be applied to the Cubic Galileon model.

We identify the coronal sources of the solar winds sampled by the ACE spacecraft during 1999-2008, and examine the in situ solar wind properties as a function of wind sources. The standard two-step mapping technique is adopted to establish the photospheric footpoints of the magnetic flux tubes along which the ACE winds flow. The footpoints are then placed in the context of EIT 284~\AA\ images and photospheric magnetograms, allowing us to categorize the sources into four groups: coronal holes (CHs), active regions (ARs), the quiet Sun (QS), and "Undefined". This practice also enables us to establish the response to solar activity of the fractions occupied by each kind of solar winds, and of their speeds and O$^{7+}$/O$^{6+}$ ratios measured in situ. We find that during the maximum phase, the majority of ACE winds originate from ARs. During the declining phase, CHs and ARs are equally important contributors to the ACE solar winds. The QS contribution increases with decreasing solar activity, and maximizes in the minimum phase when QS appear to be the primary supplier of the ACE winds. With decreasing activity, the winds from all sources tend to become cooler, as represented by the increasingly low O$^{7+}$/O$^{6+}$ ratios. On the other hand, during each activity phase, the AR winds tend to be the slowest and associated with the highest O$^{7+}$/O$^{6+}$ ratios, and the CH winds correspond to the other extreme, with the QS winds lying in between. Applying the same analysis method to the slow winds only, here defined as the winds with speeds lower than 500 km s$^{-1}$, we find basically the same overall behavior, as far as the contributions of individual groups of sources are concerned. This statistical study indicates that QS regions are an important source of the solar wind during the minimum phase.

A sample of 111 Fermi blazars each with a well-established radio core luminosity, broad-line luminosity, bolometric luminosity and black hole mass has been compiled from the literatures.We present a significant correlation between radio core and broad-line emission luminosities that supports a close link between accretion processes and relativistic jets. Analysis reveals a relationship of $\rm{LogL_{BLR}\sim(0.81\pm0.06)LogL_{R}^{C}}$ which is consistant with theoretical predicted coefficient and supports that blazar jets are powered by energy extraction from a rapidly spinning Kerr black hole through the magnetic field provided by the accretion disk. Through studying the correlation between the intrinsic bolometric luminosity and the black hole mass, we find a relationship of $\rm{{Log}\frac{L_{in}}{L_{\odot}}=(0.95\pm0.26){Log}\frac{M}{M_{\odot}}+(3.53\pm2.24)}$ which supports mass-luminosity relation for Fermi blazars derived in this work is a powerlaw relation similar to that for main-sequence stars. Finally, evolutionary sequence of blazars is discussed.

We compiled a sample of 24 flat-spectrum radio quasars (FSRQs), 21 BL Lacs, 13 FRI and 12 FRII radio galaxies to study the jet power and broad-line luminosity relation, and the jet power and black hole mass relation. Furthermore, we obtain the histograms of key parameters. Our main results are as the follows:(i) We found that the FSRQs are roughly separated from BL Lacs by the Ledlow-Owen FRI/FRII dividing line in the $\rm{\log{P_{jet}}-\log{M}}$ plane. This result supports the unification model that BL Lacs unify with FRI radio galaxies and FSRQs unify with FRII radio galaxies.(ii)We found that the coefficient of the best-fit linear regression equation of $\rm{\log{L_{BLR}}-\log{P_{jet}}}$ relation is very close to 1 for our sample. The correlation between broad line luminosity and jet power is significant which supports that jet power has a close link with accretion.

The neutron-proton effective mass splitting in neutron-rich nucleonic matter reflects the space-time nonlocality of the isovector nuclear interaction. It affects the neutron/proton ratio during the earlier evolution of the Universe, cooling of protoneutron stars, structure of rare isotopes and dynamics of heavy-ion collisions. While there is still no consensus on whether the neutron-proton effective mass splitting is negative, zero or positive and how it depends on the density as well as the isospin-asymmetry of the medium, significant progress has been made in recent yeas in addressing these issues. We first recall the connections among the neutron-proton effective mass splitting, the momentum dependence of the isovector potential and the density dependence of the symmetry energy. We then make a few observations about the progress in calculating the neutron-proton effective mass splitting using various nuclear many-body theories and its effects on the isospin-dependence of in-medium nucleon-nucleon cross sections. Perhaps, our most reliable knowledge so far about the neutron-proton effective mass splitting at saturation density of nuclear matter comes from optical model analyses of huge sets of nucleon-nucleus scattering data accumulated over the last five decades. The momentum dependence of the symmetry potential from these analyses provide a useful boundary condition at saturation density for calibrating nuclear many-body calculations. Several observables in heavy-ion collisions have been identified as sensitive probes of the neutron-proton effective mass splitting in dense neutron-rich matter based on transport model simulations. We review these observables and comment on the latest experimental findings.

We present a comprehensive derivation of linear perturbation equations for different matter species, including photons, baryons, cold dark matter, scalar fields, massless and massive neutrinos, in the presence of a generic conformal coupling. Starting from the Lagrangians, we show how the conformal transformation affects the dynamics. In particular, we discuss how to incorporate consistently the scalar coupling in the equations of the Boltzmann hierarchy for massive neutrinos and the subsequent fluid approximations. We use the recently proposed K-mouflage model as an example to demonstrate the numerical implementation of our linear perturbation equations. K-mouflage is a new mechanism to suppress the fifth force between matter particles induced by the scalar coupling, but in the linear regime the fifth force is unsuppressed and can change the clustering of different matter species in different ways. We show how the CMB, lensing potential and matter power spectra are affected by the fifth force, and find ranges of K-mouflage parameters whose effects could be seen observationally. We also find that the scalar coupling can have the nontrivial effect of shifting the amplitude of the power spectra of the lensing potential and density fluctuations in opposite directions, although both probe the overall clustering of matter. This paper can serve as a reference for those who work on generic coupled scalar field cosmology, or those who are interested in the cosmological behaviour of the K-mouflage model.

The nuclear symmetry energy (E_sym(\roh)) is a vital ingredient of our understanding of many processes, from heavy-ion collisions to neutron stars structure. While the total nuclear symmetry energy at nuclear saturation density (\rho_0) is relatively well determined, its value at supranuclear densities is not. The latter can be better constrained by separately examining its kinetic and potential terms and their density dependencies. The kinetic term of the symmetry energy, E_sym^kin(\rho_0), equals the difference in the per-nucleon kinetic energy between pure neutron matter (PNM) and symmetric nuclear matter (SNM), often calculated using a simple Fermi gas model. However, experiments show that tensor force induced short-range correlations (SRC) between proton-neutron pairs shift nucleons to high-momentum in SNM, but have almost no effect in PNM. We present an approximate analytical expression for E_sym^kin(\rho) of correlated nucleonic matter. In our model, E_sym^kin(\rho_0) = -10 MeV, which differs significantly from +12.5 MeV for the widely-used free Fermi gas model. This result is consistent with our analysis of recent data on the free proton-to-neutron ratios measured in intermediate energy nucleus-nucleus collisions as well as with microscopic many-body calculations, and previous phenomenological extractions. We then use our calculated E_sym^kin(\rho) in combination with the known total symmetry energy and its density dependence at saturation density to constrain the value and density dependence of the potential part and to extrapolate the total symmetry energy to supranuclear densities.

The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.

There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperons and new phases, the most uncertain part of the nucleonic EOS is currently the density dependence of nuclear symmetry energy. To provide information that may help break the EOS-gravity degeneracy, we investigate effects of symmetry energy within its uncertain range determined by terrestrial nuclear laboratory experiments on the gravitational binding energy and spacetime curvature of neutron stars within GR and the scalar-tensor (ST) theory of gravity. In particular, we focus on effects of the following parameters characterizing the EOS of neutron-rich nucleonic matter: (1) the incompressibility $K_0$ of symmetric nuclear matter, (2) the slope $L$ of symmetry energy at saturation density and (3) the high-density behavior of symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of symmetry energy leads to large changes in both the binding energy and curvature of neutron stars while effects of varying the more constrained $K_0$ are negligibly small. The difference in predictions using the GR and the ST theory appears only for massive neutron stars, and is significantly smaller than the differences resulting from variations in the symmetry energy. We conclude that within the ST theory of gravity, the EOS-gravity degeneracy has been broken by the recent relativistic pulsar measurements, and that measurements of neutron star properties sensitive to the compactness constrain mainly the density dependence of the symmetry energy.

It is known that two supermassive black holes (SMBHs) cannot merge in a spherical galaxy within a Hubble time; an emerging picture is that galaxy geometry, rotation, and large potential perturbations may usher the SMBH binary through the critical three-body scattering phase and ultimately drive the SMBH to coalesce. We explore the orbital content within an N-body model of a mildly- flattened, non-rotating, SMBH-embedded elliptical galaxy. When used as the foundation for a study on the SMBH binary coalescence, the black holes bypassed the binary stalling often seen within spherical galaxies and merged on Gyr timescales (Khan et al. 2013). Using both frequency-mapping and angular momentum criteria, we identify a wealth of resonant orbits in the axisymmetric model, including saucers, that are absent from an otherwise identical spherical system and that can potentially interact with the binary. We quantified the set of orbits that could be scattered by the SMBH binary, and found that the axisymmetric model contained nearly seven times the number of these potential loss cone orbits compared to our equivalent spherical model. In this flattened model, the mass of these orbits is roughly 3 times of that of the SMBH, which is consistent with what the SMBH binary needs to scatter to transition into the gravitational wave regime.

We model the linear and nonlinear growth of large scale structure in the Cubic Galileon gravity model, by running a suite of N-body cosmological simulations using the {\tt ECOSMOG} code. Our simulations include the Vainshtein screening effect, which reconciles the Cubic Galileon model with local tests of gravity. In the linear regime, the amplitude of the matter power spectrum increases by $\sim 25%$ with respect to the standard $\Lambda$CDM model today. The modified expansion rate accounts for $\sim 20%$ of this enhancement, while the fifth force is responsible for only $\sim 5%$. This is because the effective unscreened gravitational strength deviates from standard gravity only at late times, even though it can be twice as large today. In the nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{-1}$), the fifth force leads to only a modest increase ($\lesssim 8%$) in the clustering power on all scales due to the very efficient operation of the Vainshtein mechanism. Such a strong effect is typically not seen in other models with the same screening mechanism. The screening also results in the fifth force increasing the number density of halos by less than 10%, on all mass scales. Our results show that the screening does not ruin the validity of linear theory on large scales which anticipates very strong constraints from galaxy clustering data. We also show that, whilst the model gives an excellent match to CMB data on small angular scales ($l \gtrsim 50$), the predicted integrated Sachs-Wolf effect is in tension with Planck/WMAP results.

Numerical simulations in modified gravity have commonly been performed under the quasi-static approximation -- that is, by neglecting the effect of time derivatives in the equation of motion of the scalar field that governs the fifth force in a given modified gravity theory. To test the validity of this approximation, we analyse the case of $f(R)$ gravity beyond this quasi-static limit, by considering effects, if any, these terms have in the matter and velocity divergence cosmic fields. To this end, we use the adaptive mesh refinement code ECOSMOG to study three variants ($|f_{R}|= 10^{-4}[$F4$], 10^{-5}[$F5$]$ and $10^{-6}[$F6$]$) of the Hu-Sawicki $f(R)$ gravity model, each of which refers to a different magnitude for the scalar field that generates the fifth force. We find that for F4 and F5, which show stronger deviations from standard gravity, a low-resolution simulation is enough to conclude that time derivatives make a negligible contribution to the matter distribution. The F6 model shows a larger deviation from the quasi-static approximation, but one that diminishes when re-simulated at higher resolution. We therefore come to the conclusion that the quasi-static approximation is valid for the most practical applications in $f(R)$ cosmologies.

Magnetohydrodynamic waves interact with structured plasmas and reveal the internal magnetic and thermal structures therein, thereby having seismological applications in the solar atmosphere. We investigate the evolution of fast magnetoacoustic pulses in randomly structured plasmas, in the context of large-scale propagating waves in the solar atmosphere. We perform one dimensional numerical simulations of fast wave pulses propagating perpendicular to a constant magnetic field in a low-$\beta$ plasma with a random density profile across the field. Both linear and nonlinear regimes are considered. We study how the evolution of the pulse amplitude and width depends on their initial values and the parameters of the random structuring. A randomly structured plasma acts as a dispersive medium for a fast magnetoacoustic pulse, causing amplitude attenuation and broadening of the pulse width. After the passage of the main pulse, secondary propagating and standing fast waves appear in the plasma. Width evolution of both linear and nonlinear pulses can be well approximated by linear functions; however, narrow pulses may have zero or negative broadening. This arises because a narrow pulse is prone to splitting, while a broad pulse usually deviates less from their initial Gaussian shape and form ripple structures on top of the main pulse. A linear pulse decays at almost a constant rate, while a nonlinear pulse decays exponentially. A pulse interacts most efficiently with a random medium which has a correlation length of about half of its initial pulse width. The development of a detailed model of a fast MHD pulse propagating in highly structured medium substantiates the interpretation of EIT waves as fast magnetoacoustic waves. Evolution of a fast pulse provides us with a novel method to diagnose the sub-resolution filamentation of the solar atmosphere.

We analyse modelling techniques for the large-scale structure formed in scalar-tensor theories of constant Brans-Dicke parameter which match the concordance model background expansion history and produce a chameleon suppression of the gravitational modification in high-density regions. Thereby, we use a mass and environment dependent chameleon spherical collapse model, the Sheth-Tormen halo mass function and linear halo bias, the Navarro-Frenk-White halo density profile, and the halo model. Furthermore, using the spherical collapse model, we extrapolate a chameleon mass-concentration scaling relation from a LCDM prescription calibrated to N-body simulations. We also provide constraints on the model parameters to ensure viability on local scales. We test our description of the halo mass function and nonlinear matter power spectrum against the respective observables extracted from large-volume and high-resolution N-body simulations in the limiting case of f(R) gravity, corresponding to a vanishing Brans-Dicke parameter. We find good agreement between the two; the halo model provides a good qualitative description of the shape of the relative enhancement of the f(R) matter power spectrum with respect to LCDM caused by the extra attractive gravitational force but fails to recover the correct amplitude. Introducing an effective linear power spectrum in the computation of the two-halo term to account for an underestimation of the chameleon suppression at intermediate scales in our approach, we accurately reproduce the measurements from the N-body simulations.

We explore voids in dark matter and halo fields from simulations of $\Lambda$CDM and Hu-Sawicki $f(R)$ models. In $f(R)$ gravity, dark matter void abundances are greater than that of general relativity (GR). However, when using haloes to identify voids, the differences of void abundances become much smaller, but can still be told apart, in principle, at the 2, 6 and 14 $\sigma$ level for the $f(R)$ model parameter amplitudes of $|f_{R0}|=10^{-6}$, $10^{-5}$ and $10^{-4}$. In contrast, the abundance of large voids found using haloes in $f(R)$ gravity is lower than in GR. The more efficient halo formation in underdense regions makes $f(R)$ voids less empty of haloes. This counter intuitive result suggests that voids are not necessarily emptier in $f(R)$ if one looks at galaxies in voids. Indeed, the halo number density profiles of voids are not distinguishable from GR. However, the same $f(R)$ voids are more empty of dark matter. This can in principle be observed by weak gravitational lensing of voids, for which the combination of a spec-$z$ and a photo-$z$ survey over the same sky is necessary. For a volume of 1~(Gpc/$h$)$^3$, neglecting the lensing shape noise, $|f_{R0}|=10^{-5}$ and $10^{-4}$ may be distinguished from GR using the lensing tangential shear signal around voids by 4 and 8$\sigma$. The line-of-sight projection of large-scale structure is the main systematics that limits the significance of this signal, limiting the constraining power for $|f_{R0}|=10^{-6}$. The halo void abundance being smaller and the steepening of dark matter void profiles in $f(R)$ models are unique features that can be combined to break the degeneracy between $|f_{R0}|$ and $\sigma_8$. The outflow of mass from void centers and velocity dispersions are greater in $f(R)$. Model differences in velocity profiles imply potential powerful constraints of the model in phase space and in redshift space.

We investigate void properties in $f(R)$ models using N-body simulations, focusing on their differences from General Relativity (GR) and their detectability. In the Hu-Sawicki $f(R)$ modified gravity (MG) models, the halo number density profiles of voids are not distinguishable from GR. In contrast, the same $f(R)$ voids are more empty of dark matter, and their profiles are steeper. This can in principle be observed by weak gravitational lensing of voids, for which the combination of a spectroscopic redshift and a lensing photometric redshift survey over the same sky is required. Neglecting the lensing shape noise, the $f(R)$ model parameter amplitudes $|f_{R0}|=10^{-5}$ and $10^{-4}$ may be distinguished from GR using the lensing tangential shear signal around voids by 4 and 8$\sigma$ for a volume of 1~(Gpc/$h$)$^3$. The line-of-sight projection of large-scale structure is the main systematics that limits the significance of this signal for the near future wide angle and deep lensing surveys. For this reason, it is challenging to distinguish $|f_{R0}|=10^{-6}$ from GR. We expect that this can be overcome with larger volume. The halo void abundance being smaller and the steepening of dark matter void profiles in $f(R)$ models are unique features that can be combined to break the degeneracy between $|f_{R0}|$ and $\sigma_8$.

(Abridged) Despite the success of Cold Dark Matter (CDM) in explaining a wide range of observations, the microscopic nature of dark matter is still unknown. Attempts to detect WIMPs, the most commonly studied form, have not yet succeeded. Some predictions are even in apparent conflict with observations (e.g. cuspy-cored halos and `missing satellite problem'). Therefore, it is important to consider different forms of CDM. We focus on dark matter comprised of ultralight bosons that form a Bose-Einstein Condensate (BEC), described by a complex scalar field, and solve the Klein-Gordon and Einstein field equations to evolve the Friedmann-Robertson-Walker universe. We find, in addition to phases of radiation-, matter- and \Lambda-domination, an earlier phase of scalar-field-domination. Moreover, while WIMP CDM is non-relativistic at all times after it decouples, BEC scalar field dark matter (SFDM) is relativistic at early times, evolving from stiff to radiation-like, before it becomes non-relativistic (CDM-like) at late times. The timing of transitions between these phases yields fundamental constraints on SFDM model parameters, particle mass m and self-interaction coupling strength \lambda. We derive the range of particle parameters required to match observations of the evolving background universe, including the CMB and abundances of light elements produced by BBN, characterized by N_eff, the effective number of neutrino species, and the epoch of matter-radiation equality z_eq. This yields m >= 2.4*10^{-21} eV/c^2 and 9.5*10^{-19} eV^{-1}cm^3 <= \lambda/(mc^2)^2 <= 4*10^{-17} eV^{-1}cm^3. Our model accommodates current observations in which N_eff at BBN is higher than at z_eq, as probed by the CMB, otherwise unexplained by WIMP CDM. SFDM without self-interaction (`Fuzzy Dark Matter') cannot satisfy current BBN constraints within 68% confidence and is therefore disfavored.

Cosmic voids are becoming key players in testing the physics of our Universe. Here we concentrate on the abundances and the dynamics of voids as these are among the best candidates to provide information on cosmological parameters. Cai, Padilla \& Li (2014) use the abundance of voids to tell apart Hu \& Sawicki $f(R)$ models from General Relativity. An interesting result is that even though, as expected, voids in the dark matter field are emptier in $f(R)$ gravity due to the fifth force expelling away from the void centres, this result is reversed when haloes are used to find voids. The abundance of voids in this case becomes even lower in $f(R)$ compared to GR for large voids. Still, the differences are significant and this provides a way to tell apart these models. The velocity field differences between $f(R)$ and GR, on the other hand, are the same for halo voids and for dark matter voids. Paz et al. (2013), concentrate on the velocity profiles around voids. First they show the necessity of four parameters to describe the density profiles around voids given two distinct void populations, voids-in-voids and voids-in-clouds. This profile is used to predict peculiar velocities around voids, and the combination of the latter with void density profiles allows the construction of model void-galaxy cross-correlation functions with redshift space distortions. When these models are tuned to fit the measured correlation functions for voids and galaxies in the Sloan Digital Sky Survey, small voids are found to be of the void-in-cloud type, whereas larger ones are consistent with being void-in-void. This is a novel result that is obtained directly from redshift space data around voids. These profiles can be used to remove systematics on void-galaxy Alcock-Pacinsky tests coming from redshift-space distortions.

We reexamine the screening mechanism in $f(R)$ gravity using N-body simulations. By explicitly examining the relation between the extra scalar field $\delta f_R$ and the gravitational potential $\phi$ in the perturbed Universe, we find that the relation between these two fields plays an important role in understanding the screening mechanism. We show that the screening mechanism in $f(R)$ gravity depends mainly on the depth of the potential well, and find a useful condition for identifying unscreened halos in simulations. We also discuss the potential application of our results to real galaxy surveys.

The contribution of plumes to the solar wind has been subject to hot debate in the past decades. The EUV Imaging Spectrometer (EIS) on board Hinode provides a unique means to deduce outflow velocities at coronal heights via direct Doppler shift measurements of coronal emission lines. Such direct Doppler shift measurements were not possible with previous spectrometers. We measure the outflow velocity at coronal heights in several on-disk long-duration plumes, which are located in coronal holes and show significant blue shifts throughout the entire observational period. In one case, a plume is measured 4 hours apart. The deduced outflow velocities are consistent, suggesting that the flows are quasi-steady. Furthermore, we provide an outflow velocity profile along the plumes, finding that the velocity corrected for the line-of-sight effect can reach 10 km s$^{-1}$ at 1.02 $R_{\odot}$, 15 km s$^{-1}$ at 1.03 $R_{\odot}$, and 25 km s$^{-1}$ at 1.05 $R_{\odot}$. This clear signature of steady acceleration, combined with the fact that there is no significant blue shift at the base of plumes, provides an important constraint on plume models. At the height of 1.03 $R_{\odot}$, EIS also deduced a density of 1.3$\times10^{8}$ cm$^{-3}$, resulting in a proton flux of about 4.2$\times10^9$ cm$^{-2}$s$^{-1}$ scaled to 1AU, which is an order of magnitude higher than the proton input to a typical solar wind if a radial expansion is assumed. This suggests that, coronal hole plumes may be an important source of the solar wind.

We revisit the excursion set approach to calculate void abundances in chameleon-type modified gravity theories, which was previously studied by Clampitt, Cai and Li (2013). We focus on properly accounting for the void-in-cloud effect, i.e., the growth of those voids sitting in over-dense regions may be restricted by the evolution of their surroundings. This effect may change the distribution function of voids hence affect predictions on the differences between modified gravity and GR. We show that the thin-shell approximation usually used to calculate the fifth force is qualitatively good but quantitatively inaccurate. Therefore, it is necessary to numerically solve the fifth force in both over-dense and under-dense regions. We then generalise the Eulerian void assignment method of Paranjape, Lam and Sheth (2012) to our modified gravity model. We implement this method in our Monte Carlo simulations and compare its results with the original Lagrangian methods. We find that the abundances of small voids are significantly reduced in both modified gravity and GR due to the restriction of environments. However, the change in void abundances for the range of void radii of interest for both models is similar. Therefore, the difference between models remains similar to the results from the Lagrangian method, especially if correlated steps of the random walks are used. As Clampitt, Cai and Li (2013), we find that the void abundance is much more sensitive to modified gravity than halo abundances. Our method can then be a faster alternative to N-body simulations for studying the qualitative behaviour of a broad class of theories. We also discuss the limitations and other practical issues associated with its applications.

Models of modified gravity offer promising alternatives to the concordance $\Lambda$CDM cosmology to explain the late-time acceleration of the universe. A popular such model is $f(R)$ gravity, in which the Ricci scalar in the Einstein-Hilbert action is replaced by a general function of it. We study the $f(R)$ model of Hu & Sawicki (2007), which recovers standard General Relativity in high density regimes, while reproducing the desired late-time acceleration at cosmological scales. We run a suite of high resolution zoom simulations using the ECOSMOG code to examine the effect of $f(R)$ gravity on the properties of a halo that is analogous to the Virgo cluster. We show that the velocity dispersion profiles can potentially discriminate between $f(R)$ models and $\Lambda$CDM, and provide complementary analysis of lensing signal profiles to explore the possibility to further distinguish the different $f(R)$ models. Our results confirm the techniques explored by Cabre et al. (2012) to quantify the effect of environment in the behavior of $f(R)$ gravity, and we extend them to study halo satellites at various redshifts. We find that the modified gravity effects in our models are most observable at low redshifts, and that effects are generally stronger for satellites far from the center of the main halo. We show that the screening properties of halo satellites trace very well that of dark matter particles, which means that low-resolution simulations in which subhalos are not very well resolved can in principle be used to study satellite properties. We discuss observables, particularly for halo satellites, that can potentially be used to constrain the observational viability of $f(R)$ gravity.

We study the nonlinear growth of structure in nonlocal gravity models with the aid of N-body simulation and the spherical collapse and halo models. We focus on a model in which the inverse-squared of the d'Alembertian operator acts on the Ricci scalar in the action. For fixed cosmological parameters, this model differs from $\Lambda{\rm CDM}$ by having a lower late-time expansion rate and an enhanced and time-dependent gravitational strength ($\sim 6\%$ larger today). Compared to $\Lambda{\rm CDM}$ today, in the nonlocal model, massive haloes are slightly more abundant (by $\sim 10\%$ at $M \sim 10^{14} M_{\odot}/h$) and concentrated ($\approx 8\%$ enhancement over a range of mass scales), but their linear bias remains almost unchanged. We find that the Sheth-Tormen formalism describes the mass function and halo bias very well, with little need for recalibration of free parameters. The fitting of the halo concentrations is however essential to ensure the good performance of the halo model on small scales. For $k \gtrsim 1 h/{\rm Mpc}$, the amplitude of the nonlinear matter and velocity divergence power spectra exhibits a modest enhancement of $\sim 12\%$ to $15\%$, compared to $\Lambda{\rm CDM}$ today. This suggests that this model might only be distinguishable from $\Lambda{\rm CDM}$ by future observational missions. We point out that the absence of a screening mechanism may lead to tensions with Solar System tests due to local time variations of the gravitational strength, although this is subject to assumptions about the local time evolution of background averaged quantities.

We investigate the dependence of the Vainshtein screening mechanism on the cosmic web morphology of both dark matter particles and halos as determined by ORIGAMI. Unlike chameleon and symmetron screening, which come into effect in regions of high density, Vainshtein screening instead depends on the dimensionality of the system, and screened bodies can still feel external fields. ORIGAMI is well-suited to this problem because it defines morphologies according to the dimensionality of the collapsing structure and does not depend on a smoothing scale or density threshold parameter. We find that halo particles are screened while filament, wall, and void particles are unscreened, and this is independent of the particle density. However, after separating halos according to their large scale morphological environment, we find no difference in the screening properties of halos in filaments versus halos in clusters. We find that the fifth force enhancement of dark matter particles in halos is greatest well outside the virial radius. We confirm the theoretical expectation that even if the internal field is suppressed by the Vainshtein mechanism, the object still feels the fifth force generated by the external fields, by measuring peculiar velocities and velocity dispersions of halos. Finally, we investigate the morphology and gravity model dependence of halo spins, concentrations, and shapes.

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

Solar Optical Telescope onboard Hinode observed a sunspot (AR 11836) with two light bridges (LBs) on 31 Aug 2013. We analysed a 2-hour \ion{Ca}{2} H emission intensity data set and detected strong 5-min oscillation power on both LBs and in the inner penumbra. The time-distance plot reveals that 5-min oscillation phase does not vary significantly along the thin bridge, indicating that the oscillations are likely to originate from the underneath. The slit taken along the central axis of the wide light bridge exhibits a standing wave feature. However, at the centre of the wide bridge, the 5-min oscillation power is found to be stronger than at its sides. Moreover, the time-distance plot across the wide bridge exhibits a herringbone pattern that indicates a counter-stream of two running waves originated at the bridge sides. Thus, the 5-min oscillations on the wide bridge also resemble the properties of running penumbral waves. The 5-min oscillations are suppressed in the umbra, while the 3-min oscillations occupy all three cores of the sunspot's umbra, separated by the LBs. The 3-min oscillations were found to be in phase at both sides of the LBs. It may indicate that either LBs do not affect umbral oscillations, or umbral oscillations at different umbral cores share the same source. Also, it indicates that LBs are rather shallow objects situated in the upper part of the umbra. We found that umbral flashes follow the life cycles of umbral oscillations with much larger amplitudes. They cannot propagate across LBs. Umbral flashes dominate the 3-min oscillation power within each core, however, they do not disrupt the phase of umbral oscillation.

We use the latest CMB data from Planck, together with BAO measurements, to constrain the full parameter space of Galileon gravity. We constrain separately the three main branches of the theory known as the Cubic, Quartic and Quintic models, and find that all yield a very good fit to these data. Unlike in $\Lambda{\rm CDM}$, the Galileon model constraints are compatible with local determinations of the Hubble parameter and predict nonzero neutrino masses at over $5\sigma$ significance. We also identify that the low-$l$ part of the CMB lensing spectrum may be able to distinguish between $\Lambda{\rm CDM}$ and Galileon models. In the Cubic model, the lensing potential deepens at late times on sub-horizon scales, which is at odds with the current observational suggestion of a positive ISW effect. Compared to $\Lambda$CDM, the Quartic and Quintic models predict less ISW power in the low-$l$ region of the CMB temperature spectrum, and as such are slightly preferred by the Planck data. We illustrate that residual local modifications to gravity in the Quartic and Quintic models may render the Cubic model as the only branch of Galileon gravity that passes Solar System tests.

Magnetohydrodynamic waves are important for diagnosing the physical parameters of coronal plasmas. Field-aligned flows appear frequently in coronal loops.We examine the effects of transverse density and plasma flow structuring on standing sausage modes trapped in coronal loops, and examine their observational implications. We model coronal loops as straight cold cylinders with plasma flow embedded in a static corona. An eigen-value problem governing propagating sausage waves is formulated, its solutions used to construct standing modes. Two transverse profiles are distinguished, one being the generalized Epstein distribution (profile E) and the other (N) proposed recently in Nakariakov et al.(2012). A parameter study is performed on the dependence of the maximum period $P_\mathrm{max}$ and cutoff length-to-radius ratio $(L/a)_{\mathrm{cutoff}}$ in the trapped regime on the density parameters ($\rho_0/\rho_\infty$ and profile steepness $p$) and flow parameters (magnitude $U_0$ and profile steepness $u$). For either profile, introducing a flow reduces $P_\mathrm{max}$ relative to the static case. $P_\mathrm{max}$ depends sensitively on $p$ for profile N but is insensitive to $p$ for profile E. By far the most important effect a flow introduces is to reduce the capability for loops to trap standing sausage modes: $(L/a)_{\mathrm{cutoff}}$ may be substantially reduced in the case with flow relative to the static one. If the density distribution can be described by profile N, then measuring the sausage mode period can help deduce the density profile steepness. However, this practice is not feasible if profile E better describes the density distribution. Furthermore, even field-aligned flows with magnitudes substantially smaller than the ambient Alfv\'en speed can make coronal loops considerably less likely to support trapped standing sausage modes.

We present GPU implementations of two fast force calculation methods, based on series expansions of the Poisson equation. One is the Self-Consistent Field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other is the Multipole Expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field method and optimized for collisionless galactic dynamics, but while SCF is a "pure" expansion, MEX is an expansion in just the angular part; it is thus capable of capturing radial structure easily, where SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ~ 0.1 seconds (depending on expansion cutoff), making simulations with as many as $10^8$ particles fast on a comparatively small number of nodes.

The velocity field of dark matter and galaxies reflects the continued action of gravity throughout cosmic history. We show that the low-order moments of the pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws of gravity on cosmological scales. In particular, the projected line-of-sight galaxy pairwise velocity dispersion, $\sigma_{12}(r)$, is very sensitive to the presence of modified gravity. Using a set of high-resolution N-body simulations we compute the pairwise velocity distribution and its projected line-of-sight dispersion for a class of modified gravity theories: the chameleon \fR gravity and Galileon gravity (cubic and quartic). The velocities of dark matter halos with a wide range of masses would exhibit deviations from General Relativity at the $(5-10)\sigma$ level. We examine strategies for detecting these deviations in galaxy redshift and peculiar velocity surveys. If detected, this signature would be a "smoking gun" for modified gravity.

The heating of solar coronal loops is at the center of the problem of coronal heating. Given that the origin of the fast solar wind has been tracked down to atmospheric layers with transition region or even chromospheric temperatures, it is worthy attempting to address whether the mechanisms proposed to provide the basal heating of the solar wind apply to coronal loops as well. We extend the loop studies based on a classical parallel-cascade scenario originally proposed in the solar wind context by considering the effects of loop expansion, and perform a parametric study to directly contrast the computed loop densities and electron temperatures with those measured by TRACE and YOHKOH/SXT. This comparison yields that with the wave amplitudes observationally constrained by SUMER measurements, while the computed loops may account for a significant fraction of SXT loops, they seem too hot when compared with TRACE loops. Lowering the wave amplitudes does not solve this discrepancy, introducing magnetic twist will make the comparison even less desirable. We conclude that the nanoflare heating scenario better explains ultraviolet loops, while turbulence-based steady heating mechanisms may be at work in heating a fraction of soft X-ray loops.

According to the Hugenholtz-Van Hove theorem, the nuclear symmetry energy $S(\rho)$ and its slope $L(\rho)$ at arbitrary densities can be decomposed in terms of the density and momentum dependence of the single-nucleon potentials in isospin-asymmetric nuclear matter which are potentially accessible to experiment. We quantify the correlations between several well-known isovector observables and $L(\rho)$ to locate the density range in which each isovector observable is most sensitive to the density dependence of the $S(\rho)$. We then study the correlation coefficients between those isovector observables and all the components of the $L(\rho)$. The neutron skin thickness of $^{208}$Pb is found to be strongly correlated with the $L(\rho)$ at a subsaturation density of $\rho = 0.59 \rho_0$ through the density dependence of the first-order symmetry potential. Neutron star radii are found to be strongly correlated with the $L(\rho)$ over a wide range of supra-saturation densities mainly through both the density and momentum dependence of the first-order symmetry potential. Finally, we find that although the crust-core transition pressure has a complex correlation with the $L(\rho)$, it is strongly correlated with the momentum derivative of the first-order symmetry potential, and the density dependence of the second-order symmetry potential.

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

We study spherical collapse in the Quartic and Quintic Covariant Galileon gravity models within the framework of the excursion set formalism. We derive the nonlinear spherically symmetric equations in the quasi-static and weak-field limits, focusing on model parameters that fit current CMB, SNIa and BAO data. We demonstrate that the equations of the Quintic model do not admit physical solutions of the fifth force in high density regions, which prevents the study of structure formation in this model. For the Quartic model, we show that the effective gravitational strength deviates from unity at late times ($z \lesssim 1$), becoming larger if the density is low, but smaller if the density is high. This shows that the Vainshtein mechanism at high densities is not enough to screen all of the modifications of gravity. This makes halos that collapse at $z \lesssim 1$ feel an overall weaker gravity, which suppresses halo formation. However, the matter density in the Quartic model is higher than in standard $\Lambda$CDM, which boosts structure formation and dominates over the effect of the weaker gravity. In the Quartic model there is a significant overabundance of high-mass halos relative to $\Lambda$CDM. Dark matter halos are also less biased than in $\Lambda$CDM, with the difference increasing appreciably with halo mass. However, our results suggest that the bias may not be small enough to fully reconcile the predicted matter power spectrum with LRG clustering data.

A thorough understanding of many astrophysical phenomena associated with compact objects requires reliable knowledge about both the equation of state (EOS)of super-dense nuclear matter and the theory of strong-field gravity simultaneously because of the EOS-gravity degeneracy. Currently, variations of the neutron star (NS) mass-radius correlation from using alternative gravity theories are much larger than those from changing the NS matter EOS within known constraints. At least two independent observables are required to break the EOS-gravity degeneracy. Using model EOSs for hybrid stars and a Yukawa-type non-Newtonian gravity, we investigate both the mass-radius correlation and pulsating frequencies of NSs. While the maximum mass of NSs increases, the frequencies of the $f$, $p_1$, $p_2$, and $w_I$ pulsating modes are found to decrease with the increasing strength of the Yukawa-type non-Newtonian gravity, providing a useful reference for future determination simultaneously of both the strong-field gravity and the supranuclear EOS by combining data of x-ray and gravitational wave emissions of neutron stars.

We study the late-time Integrated Sachs-Wolfe (ISW) effect in $f(R)$ gravity using N-body simulations. In the $f(R)$ model under study, the linear growth rate is larger than that in general relativity (GR). This slows down the decay of the cosmic potential and induces a smaller ISW effect on large scales. Therefore, the $\dot\Phi$ (time derivative of the potential) power spectrum at $k<0.1h$/Mpc is suppressed relative to that in GR. In the non-linear regime, relatively rapid structure formation in $f(R)$ gravity boosts the non-linear ISW effect relative to GR, and the $\dot\Phi$ power spectrum at $k>0.1h$/Mpc is increased (100$\%$ greater on small scales at $z=0$). We explore the detectability of the ISW signal via stacking supercluster and supervoids. The differences in the corresponding ISW cold or hot spots are $\sim 20\%$ for structures of $\sim 100$Mpc/$h$. Such differences are greater for smaller structures, but the amplitude of the signal is lower. The high amplitude of ISW signal detected by Granett et al. can not explained in the $f(R)$ model. We find relatively big differences between $f(R)$ and GR in the transverse bulk motion of matter, and discuss its detectability via the relative frequency shifts of photons from multiple lensed images.

We study how the filamentary pattern of the cosmic web changes if the true gravity deviates from the general relativity (GR) on the large scale. The f(R) gravity whose strength is controlled to satisfy the current observational constraints on the cluster scale is adopted as our fiducial model and a large N-body simulation of high-resolution is utilized for this study. By applying the minimal spanning tree algorithm to the halo catalogs from the simulation at various epochs, we identify the main stems of the rich superclusters located in the most prominent filamentary section of the cosmic web and determine their spatial extents per member cluster as the degree of their straightness. It is found that the $f(R)$ gravity has an effect of significantly bending the superclusters and that the effect becomes stronger as the universe evolves. Even in the case where the deviation from GR is too small to be detectable by any other observables, the degree of the supercluster straightness exhibits conspicuous difference between the $f(R)$ and the GR models. Our results also imply that the supercluster straightness could be a useful discriminator of f(R) gravity from the coupled dark energy since it is shown to evolve differently between the two models. As a final conclusion, the degree of the straightness of the rich superclusters should provide a powerful cosmological test of large scale gravity.

In this presentation, we report for the first time the detailed temperature evolution process of the magnetic flux rope in a failed solar eruption. Occurred on January 05, 2013, the flux rope was impulsively accelerated to a speed of ~ 400 km/s in the first minute, then decelerated and came to a complete stop in two minutes. The failed eruption resulted in a large-size high-lying (~ 100 Mm above the surface) high-temperature "fire ball" sitting in the corona for more than two hours. The time evolution of the thermal structure of the flux rope was revealed through the differential emission measure analysis technique, which produced temperature maps using observations of the Atmospheric Imaging Assembly on board Solar Dynamic Observatory. The average temperature of the flux rope steadily increased from ~ 5 MK to ~ 10 MK during the first nine minutes of the evolution, which was much longer than the rise time (about three minutes) of the associated soft X-ray flare. We suggest that the flux rope be heated by the energy release of the continuing magnetic reconnection, different from the heating of the low-lying flare loops, which is mainly produced by the chromospheric plasma evaporation. The loop arcade overlying the flux rope was pushed up by ~ 10 Mm during the attempted eruption. The pattern of the velocity variation of the loop arcade strongly suggests that the failure of the eruption be caused by the strapping effect of the overlying loop arcade.

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

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

X-ray observations of the neutron star in the Cas A supernova remnant over the past decade suggest the star is undergoing a rapid drop in surface temperature of $\approx$ $2-5.5\%$. One explanation suggests the rapid cooling is triggered by the onset of neutron superfluidity in the core of the star, causing enhanced neutrino emission from neutron Cooper pair breaking and formation (PBF). Using consistent neutron star crust and core equations of state (EOSs) and compositions, we explore the sensitivity of this interpretation to the density dependence of the symmetry energy $L$ of the EOS used, and to the presence of enhanced neutrino cooling in the bubble phases of crustal "nuclear pasta". Modeling cooling over a conservative range of neutron star masses and envelope compositions, we find $L\lesssim70$ MeV, competitive with terrestrial experimental constraints and other astrophysical observations. For masses near the most likely mass of $M\gtrsim 1.65 M_{\odot}$, the constraint becomes more restrictive $35\lesssim L\lesssim 55$ MeV. The inclusion of the bubble cooling processes decreases the cooling rate of the star during the PBF phase, matching the observed rate only when $L\lesssim45$ MeV, taking all masses into consideration, corresponding to neutron star radii $\lesssim 11$km.

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion.

In this paper, we study the non-linear matter power spectrum in a specific family of $f(R)$ models that can reproduce the $\Lambda$CDM background expansion history, using high resolution $N$-body simulations based on the {\sc ecosmog} code. We measure the matter power spectrum in the range of $0.05h{\rm Mpc}^{-1}<k<10h{\rm Mpc}^{-1}$ from simulations for our $f(R)$ models and give theoretical explanations to their behaviour and evolution patterns. We also examine the chameleon mechanism for our models and find that it works throughout the cosmic history in dense regions, for our $f(R)$ models with $|f_{R0}|<10^{-4}$. On the other hand, for models with $|f_{R0}|>10^{-3}$, we find no chameleon screening in dense regions at late times ($z<3$), which means that those models could be ruled out due to the factor-of-$1/3$ enhancement to the strength of Newtonian gravity. We also give the best-fit parameters for a generalised PPF fitting formula which works well for the models studied here.

Current constraints on f(R) gravity from the large-scale structure are at the verge of penetrating into a region where the modified forces become nonlinearly suppressed. For a consistent treatment of observables at these scales, we study cluster quantities produced in chameleon and linearized Hu-Sawicki f(R) gravity dark matter N-body simulations. We find that the standard Navarro-Frenk-White halo density profile and the radial power law for the pseudo phase-space density provide equally good fits for f(R) clusters as they do in the Newtonian scenario. We give qualitative arguments for why this should be the case. For practical applications, we derive analytic relations, e.g., for the f(R) scalar field, the gravitational potential, and the velocity dispersion as seen within the virialized clusters. These functions are based on three degrees of freedom fitted to simulations, i.e., the characteristic density, scale, and velocity dispersion. We further analyze predictions for these fitting parameters from the gravitational collapse and the Jeans equation, which are found to agree well with the simulations. Our analytic results can be used to consistently constrain chameleon f(R) gravity with future observations on virialized cluster scales without the necessity of running a large number of simulations.

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

On cosmological scales, observations of the cluster abundance currently place the strongest constraints on f(R) gravity. These constraints lie in the large-field limit, where the modifications of general relativity can correctly be modeled by setting the Compton wavelength of the scalar field to its background value. These bounds are, however, at the verge of penetrating into a regime, where the modifications become nonlinearly suppressed due to the chameleon mechanism and cannot be described by this linearized approximation. For future constraints based on observations subjected to cluster abundance, it is therefore essential to consistently model the chameleon effect. We analyze descriptions of the halo mass function in chameleon f(R) gravity using a mass- and environment-dependent spherical collapse model in combination with excursion set theory and phenomenological fits to N-body simulations in the LCDM and f(R) gravity scenarios. Our halo mass functions consistently incorporate the chameleon suppression and cosmological parameter dependencies, improving upon previous formalisms and providing an important extension to N-body simulations for the application in consistent tests of gravity with observables sensitive to the abundance of clusters.

Infrared modifications of General Relativity (GR) can be revealed by comparing the mass of galaxy clusters estimated from weak lensing to that from infall kinematics. We measure the 2D galaxy velocity distribution in the cluster infall region by applying the galaxy infall kinematics (GIK) model developed by Zu and Weinberg (2013) to two suites of f(R) and Galileon modified gravity simulations. Despite having distinct screening mechanisms, namely, the Chameleon and the Vainshtein effects, the f(R) and Galileon clusters exhibit very similar deviations in their GIK profiles from GR, with ~ 100-200 k/s enhancement in the characteristic infall velocity at r=5 Mpc/h and 50-100 km/s broadening in the radial and tangential velocity dispersions across the entire infall region, for clusters with mass ~ 10^{14} Msol/h at z=0.25. These deviations are detectable via the GIK reconstruction of the redshift--space cluster-galaxy cross-correlation function, xi_cg^s(r_p,r_\pi), which shows ~ 1-2 Mpc/h increase in the characteristic line-of-sight distance r_\pi^c at r_p<6 Mpc/h from GR predictions. With overlapping deep imaging and large redshift surveys in the future, we expect that the GIK modelling of xi_cg^s, in combination with the stacked weak lensing measurements, will provide powerful diagnostics of modified gravity theories and the origin of cosmic acceleration.

We present the first constraints on the full parameter space of the Galileon modified gravity model, considering both the cosmological parameters and the coefficients which specify the additional terms in the Lagrangian due to the Galileon field, which we call the Galileon parameters. We use the latest cosmic microwave background measurements, along with distance measurements from supernovae and baryonic acoustic oscillations, performing a Monte Carlo Markov Chain exploration of the 9-dimensional parameter space. The integrated Sachs-Wolfe signal can be very different in Galileon models compared to standard gravity, making it essential to use the full CMB data rather than the CMB distance priors. We demonstrate that meaningful constraints are only possible in the Galileon parameter space after taking advantage of a scaling degeneracy. We find that the Galileon model can fit the WMAP 9-year results better than the standard \Lambda-Cold Dark Matter model, but gives a slightly worse fit overall once lower redshift distance measurements are included. The best-fitting cosmological parameters (e.g. matter density, scalar spectral index, fluctuation amplitude) can differ by more than 2\sigma\ in the Galileon model compared with \Lambda CDM. We highlight other potential constraints of the Galileon model using galaxy clustering and weak lensing measurements.

According to the Hugenholtz-Van Hove theorem, nuclear symmetry energy \esym and its slope \lr at an arbitrary density $\rho$ are determined by the nucleon isovector (symmetry) potential \usym and its momentum dependence $\frac{\partial U_{sym}}{\partial k}$. The latter determines uniquely the neutron-proton effective k-mass splitting $m^*_{n-p}(\rho,\delta)\equiv (m_{\rm n}^*-m_{\rm p}^*)/m$ in neutron-rich nucleonic matter of isospin asymmetry $\delta$. Using currently available constraints on the \es0 and \l0 at normal density $\rho_0$ of nuclear matter from 28 recent analyses of various terrestrial nuclear laboratory experiments and astrophysical observations, we try to infer the corresponding neutron-proton effective k-mass splitting $m^*_{n-p}(\rho_0,\delta)$. While the mean values of the $m^*_{n-p}(\rho_0,\delta)$ obtained from most of the studies are remarkably consistent with each other and scatter very closely around an empirical value of \emass$=0.27\cdot\delta$, it is currently not possible to scientifically state surely that the \emass is positive within the present knowledge of the uncertainties. Quantifying, better understanding and then further reducing the uncertainties using modern statistical and computational techniques in extracting the \es0 and \l0 from analyzing the experimental data are much needed.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

Fundamental standing modes and their overtones play an important role in coronal seismology. We examine how a significant field-aligned flow affects standing modes supported by coronal loops, modeled here as cold magnetic slabs. Of particular interest are the period ratios of the fundamental to its $(n-1)$-th overtone ($P_1/nP_n$) for both kink and sausage modes, and the threshold half-width-to-length ratio for sausage modes. For standing kink modes, the flow significantly reduces $P_1/nP_n$ in general, the effect being particularly strong for larger $n$ and when the density contrast $\rho_0/\rho_e$ between loops and their surroundings is weak. That said, even when $\rho_0/\rho_e$ approaches infinity, this effect is still substantial, reducing the minimal $P_1/nP_n$ by up to 13.7% (24.5%) for $n=2$ ($n=4$) relative to the static case, when the Alfv\'en Mach number $M_A$ reaches 0.8 where $M_A$ measures the loop flow speed in units of the internal Alfv\'en speed. For standing sausage modes, though not negligible, the flow effect in reducing $P_1/nP_n$ is not as strong. However, the threshold half-width-to-length ratio is considerably larger in the flowing case than its static counterpart. For $\rho_0/\rho_e$ in the range $[9, 1024]$ and $M_A$ in the range $[0, 0.5]$, an exhaustive parameter study yields that this threshold is well fitted by our Equation (23) which involves the two parameters in a simple way. This allows one to analytically constrain the combination $(\rho_0/\rho_e, M_A)$ for a loop with known width-to-length ratio when a standing sausage oscillation is identified therein. It also allows one to further examine the idea of partial sausage modes, and the flow is found to reduce significantly the spatial extent where partial modes are allowed.

Gravitational wave (GW) astronomy opens up an entirely new window on the Universe to probe the equations of state (EOS) of neutron-rich matter. With the advent of next generation GW detectors, measuring the gravitational radiation from coalescing binary neutron star systems, mountains on rotating neutron stars, and stellar oscillation modes may become possible in the near future. Using a set of model EOSs satisfying the latest constraints from terrestrial nuclear experiments, state of the art nuclear many-body calculations of the pure neutron matter EOS, and astrophysical observations consistently, we study various GW signatures of the high-density behavior of the nuclear symmetry energy, which is considered among the most uncertain properties of dense neutron-rich nucleonic matter. In particular, we find the tidal polarizability of neutron stars, potentially measurable in binary systems just prior to merger, is more sensitive to the high density component of the nuclear symmetry energy than the symmetry energy at nuclear saturation density. We also find that the upper limit on the GW strain amplitude from elliptically deformed stars is very sensitive to the density dependence of the symmetry energy. This suggests that future developments in modeling of the neutron star crust, and direct gravitational wave signals from accreting binaries will provide a wealth of information on the EOS of neutron-rich matter. We also review the sensitivity of the $r$-mode instability window to the density dependence of the symmetry energy. Whereas models with larger values of the density slope of the symmetry energy at saturation seem to be disfavored by the current observational data, within a simple $r$-mode model, we point out that a subsequent softer behavior of the symmetry energy at high densities (hinted at by recent observational interpretations) could rule them in.

There is growing interest in testing alternative gravity theories using the subtle gravitational redshifts in clusters of galaxies. However, current models all neglect a transverse Doppler redshift of similar magnitude, and some models are not self-consistent. An equilibrium model would fix the gravitational and transverse Doppler velocity shifts to be about 6sigma^2/c and 3sigma^2/2c in order to fit the observed velocity dispersion sigma self-consistently. This result comes from the Virial Theorem for a spherical isotropic cluster, and is insensitive to the theory of gravity. A gravitational redshift signal also does not directly distinguish between the Einsteinian and f(R) gravity theories, because each theory requires different dark halo mass function to keep the clusters in equilibrium. When this constraint is imposed, the gravitational redshift has no sensitivity to theory. Indeed our N-body simulations show that the halo mass function differs in f(R), and that the transverse Doppler effect is stronger than analytically predicted due to non-equilibrium.

We use a suite of high resolution state-of-the-art N-body Dark Matter simulations of chameleon $f(R)$ gravity to study the higher order volume averaged correlation functions $\bar{\xi_n}$ together with the hierarchical $n$-th order correlation amplitudes $S_n=\bar{\xi}_n/\bar{\xi}_2^{n-1}$ and density distribution functions (PDF). We show that under the non-linear modifications of gravity the hierarchical scaling of the reduced cumulants is preserved. This is however characterised by significant changes of both the $\bar{\xi_n}$ and $S_n$'s values and their scale dependence with respect to General Relativity gravity (GR). In addition, we measure a significant increase of the non linear $\sigma_8$ parameter reaching $14, 5$ and 0.5% in excess of the GR value for the three flavours of our $f(R)$ models. We further note that the values of the reduced cumulants up to order $n=9$ are significantly increased in $f(R)$ gravity for all our models at small scales $R\simlt 30\hmpc$. In contrast the values of the hierarchical amplitudes, $S_n$'s, are smaller in $f(R)$ indicating that the modified gravity density distribution functions are deviating from the GR case. Furthermore we find that the redshift evolution of relative deviations of the $f(R)$ hierarchical correlation amplitudes is fastest at high and moderate redshifts $1\leq z \leq4$. The growth of these deviations significantly slows down in the low redshift universe. We also compute the PDFs and show that for scales below $\sim 20\hmpc$ they are significantly shifted in $f(R)$ gravity towards the low densities. Finally we discuss the implications of our theoretical predictions for measurements of the hierarchical clustering in galaxy redshifts surveys, including the important problems of the galaxy biasing and redshifts space distortions.

We examine the electron-ion streaming instabilities driven by drift velocities of the order of the electron thermal velocity in a nonmagnetized plasma by one-dimensional electrostatic particle-in-cell code which adopts an ion-to-electron mass ratio of 1600. An initial state is set up where the ion bulk speed is zero while the electrons drift relative to ions, and where electrons are much hotter. We examine in detail four runs where drift velocity is systematically varied from lower than to larger than the electron thermal velocity. In all runs the Langmuir waves with Doppler-shifted frequencies dominate early on when streaming instabilities are too weak to discern. And then intense ion-acoustic waves or Buneman instabilities appear, which tend to be accompanied by localized electron and ion beams. Ion-acoustic modes and Buneman modes co-exist in the system when the initial drift velocity is just over the electron thermal speed. Beam modes are excited when the localized beams with large enough velocities appear. In the developed stage of instabilities, the direction in which density depressions propagate is always opposite to that of the localized ion beams. When the initial drift velocity is close to the electron thermal speed, categorizing the relevant instabilities is not easy, and one needs to examine in detail the wave dispersion diagrams at various stages of the evolution of the system.

There has been a lot of research interest in modified gravity theories which utilise the Vainshtein mechanism to recover standard general relativity in regions with high matter density, such as the Dvali-Gabadadze-Porrati and Galileon models. The strong nonlinearity in the field equations of these theories implies that accurate theoretical predictions could only be made using high-resolution cosmological simulations. Previously, such simulations were usually done on regular meshes, which limits both their performance and the accuracy. In this paper, we report the development of a new algorithm and code, based on ECOSMOG, that uses adaptive mesh refinements to improve the efficiency and precision in simulating the models with Vainshtein mechanism. We have made various code tests against the numerical reliability, and found consistency with previous simulations. We also studied the velocity field in the self-accelerating branch of the DGP model. The code, parallelised using MPI, is suitable for large cosmological simulations of Galileon-type modified gravity theories.

In this work we systematically study the linear and nonlinear structure formation in chameleon theories of modified gravity, using a generic parameterisation which describes a large class of models using only 4 parameters. For this we have modified the N-body simulation code ECOSMOG to perform a total of 65 simulations for different models and parameter values, including the default LCDM. These simulations enable us to explore a significant portion of the parameter space. We have studied the effects of modified gravity on the matter power spectrum and mass function, and found a rich and interesting phenomenology where the difference with the LCDM paradigm cannot be reproduced by a linear analysis even on scales as large as k~0.05h/Mpc, since the latter incorrectly assumes that the modification of gravity depends only on the background matter density. Our results show that the chameleon screening mechanism is significantly more efficient than other mechanisms such as the dilaton and symmetron, especially in high-density regions and at early times, and can serve as a guidance to determine the parts of the chameleon parameter space which are cosmologically interesting and thus merit further studies in the future.

The k-filtering technique and wave polarization analysis are applied to Cluster magnetic field data to study plasma turbulence at the scale of the ion gyroradius in the fast solar wind. Waves are found propagating in directions nearly perpendicular to the background magnetic field at such scales. The frequencies of these waves in the solar wind frame are much smaller than the proton gyro-frequency. After the wave vector ${\bf k}$ is determined at each spacecraft frequency $f_{sc}$, wave polarization property is analyzed in the plane perpendicular to ${\bf k}$. Magnetic fluctuations have $\delta B_\perp>\delta B_\parallel$ (here the $\parallel$ and $\perp$ refer to the background magnetic field ${\bf B}_0$). The wave magnetic field has right-handed polarization at propagation angles $\theta_{\bf kB}<90^\circ$ and $>90^\circ$. The magnetic field in the plane perpendicular to ${\bf B}_0$ however has no clear sense of a dominant polarization but local rotations. We discuss the merits and limitations of linear kinetic Aflv\'en waves (KAWs) and coherent Alfv\'en vortices in the interpretation of the data. We suggest that the fast solar wind turbulence may be populated with KAWs, small scale current sheets and Alfv\'en vortices at ion kinetic scales.

We investigate the behavior of the fifth force in voids in chameleon models using the spherical collapse method. Contrary to Newtonian gravity, we find the fifth force is repulsive in voids. The strength of the fifth force depends on the density inside and outside the void region as well as its radius. It can be many times larger than the Newtonian force and their ratio is in principle unbound. This is very different from the case in halos, where the fifth force is no more than 1/3 of gravity. The evolution of voids is governed by the Newtonian gravity, the effective dark energy force and the fifth force. While the first two forces are common in both LCDM and chameleon universes, the fifth force is unique to the latter. Driven by the outward-pointing fifth force, individual voids in chameleon models expand faster and grow larger than in a LCDM universe. The expansion velocity of the void shell can be 20% to 30% larger for voids of a few Mpc/h in radius, while their sizes can be larger by ~10%. These differences are smaller for larger voids of the same density. We compare void statistics using excursion set theory; for voids of the same size, their number density is larger in chameleon models. The fractional difference increases with void size. The chance of having voids of radius ~25 Mpc/h can be 2.5 times larger. This difference is about 10 times larger than that in the halo mass function. We find strong environmental dependence of void properties in chameleon models. The differences in size and expansion velocity with GR are both larger for small voids in high density regions. In general, the difference between chameleon models and LCDM in void properties (size, expansion velocity and distribution function) are larger than the corresponding quantities for halos. This suggests that voids might be better candidates than halos for testing gravity.

In the applications of solar magneto-seismology(SMS), employing the ratio of the period of the fundamental mode to twice the one of its first overtone, $P_1/2P_2$, plays an important role. We examine how field-aligned flows affect the dispersion properties, and hence the period ratios, of standing modes supported by magnetic slabs in the solar atmosphere. We numerically solve the dispersion relations and devise a graphic means to construct standing modes. For coronal slabs, we find that the flow effects are significant, for the fast kink and sausage modes alike. For the kink ones, they may reduce $P_1/2P_2$ by up to 23% compared with the static case, and the minimum allowed $P_1/2P_2$ can fall below the lower limit analytically derived for static slabs. For the sausage modes, while introducing the flow reduces $P_1/2P_2$ by typically $\lesssim 5$ % relative to the static case, it significantly increases the threshold aspect ratio only above which standing sausage modes can be supported, meaning that their detectability is restricted to even wider slabs. In the case of photospheric slabs, the flow effect is not as strong. However, standing modes are distinct from the coronal case in that standing kink modes show a $P_1/2P_2$ that deviates from unity even for a zero-width slab, while standing sausage modes no longer suffer from a threshold aspect ratio. We conclude that transverse structuring in plasma density and flow speed should be considered in seismological applications of multiple periodicities to solar atmospheric structures.

According to the Hugenholtz-Van Hove theorem, the symmetry energy and its density slope can be decomposed uniquely in terms of the single-nucleon potential in asymmetric nuclear matter which, at normal density, can be constrained by the nucleon optical model potential extracted from analyzing the nucleon-nucleus scattering data. To more accurately extract information on the symmetry energy and its density slope at normal density from neutron-nucleus scattering data, going beyond the well-known Lane potential, we include consistently the second order terms in isospin asymmetry in both the optical model potential and the symmetry energy decomposition. We find that the strength of the second-order symmetry potential $U_{sym,2}$ in asymmetric nuclear matter is significant compared to the first-order one U_{sym,1}$, especially at high nucleon momentum. While the $U_{sym,1}$ at normal density decreases with nucleon momentum, the $U_{sym,2}$ is found to have the opposite behavior. Moreover, we discuss effects of the $U_{sym,1}$ and $U_{sym,2}$ on determining the density dependence of the symmetry energy, and we find that the available neutron-nucleus scattering data favor a softer density dependence of the symmetry energy at normal density with a negligible contribution from the $U_{sym,2}$.

Effects of the non-Newtonian gravity on properties of finite nuclei are studied by consistently incorporating both the direct and exchange contribution of the Yukawa potential in the Hartree-Fock approach using a well-tested Skyrme force for the strong interaction. It is shown for the first time that the strength of the Yukawa term in the non-Newtonian gravity is limited to $\log(|\alpha|)<1.75/[\lambda(\rm fm)]^{0.54} + 33.6$ within the length scale of $\lambda=1-10$ fm in order for the calculated properties of finite nuclei not to be in conflict with accurate experimental data available.

We study the cosmology of Galileon modified gravity models in the linear perturbation regime. We derive the fully covariant and gauge invariant perturbed field equations using two different methods, which give consistent results, and solve them using a modified version of the {\tt CAMB} code. We find that, in addition to modifying the background expansion history and therefore shifting the positions of the acoustic peaks in the cosmic microwave background (CMB) power spectrum, the Galileon field can cluster strongly from early times, and causes the Weyl gravitational potential to grow, rather than decay, at late times. This leaves clear signatures in the low-$l$ CMB power spectrum through the modified integrated Sachs-Wolfe effect, strongly enhances the linear growth of matter density perturbations and makes distinctive predictions for other cosmological signals such as weak lensing and the power spectrum of density fluctuations. The quasi-static approximation is shown to work quite well from small to the near-horizon scales. We demonstrate that Galileon models display a rich phenomenology due to the large parameter space and the sensitive dependence of the model predictions on the Galileon parameters. Our results show that some Galileon models are already ruled out by present data and that future higher significance galaxy clustering, ISW and lensing measurements will place strong constraints on Galileon gravity.

We present a systematic survey the range of predictions of the neutron star inner crust composition, crust-core transition densities and pressures, and density range of the nuclear `pasta' phases at the bottom of the crust provided by the compressible liquid drop model in the light of current experimental and theoretical constraints on model parameters. Using a Skyrme-like model for nuclear matter, we construct baseline sequences of crust models by consistently varying the density dependence of the bulk symmetry energy at nuclear saturation density, $L$, under two conditions: (i) that the magnitude of the symmetry energy at saturation density $J$ is held constant, and (ii) $J$ correlates with $L$ under the constraint that the pure neutron matter (PNM) EoS satisfies the results of ab-initio calculations at low densities. Such baseline crust models facilitate consistent exploration of the $L$ dependence of crustal properties. The remaining surface energy and symmetric nuclear matter parameters are systematically varied around the baseline, and different functional forms of the PNM EoS at sub-saturation densities implemented, to estimate theoretical `error bars' for the baseline predictions. Inner crust composition and transition densities are shown to be most sensitive to the surface energy at very low proton fractions and to the behavior of the sub-saturation PNM EoS. Recent calculations of the energies of neutron drops suggest that the low-proton-fraction surface energy might be higher than predicted in Skyrme-like models, which our study suggests may result in a greatly reduced volume of pasta in the crust than conventionally predicted.

The modeling of many neutron star observables incorporates the microphysics of both the stellar crust and core, which is tied intimately to the properties of the nuclear matter equation of state (EoS). We explore the predictions of such models over the range of experimentally constrained nuclear matter parameters, focusing on the slope of the symmetry energy at nuclear saturation density $L$. We use a consistent model of the composition and EoS of neutron star crust and core matter to model the binding energy of pulsar B of the double pulsar system J0737-3039, the frequencies of torsional oscillations of the neutron star crust and the instability region for r-modes in the neutron star core damped by electron-electron viscosity at the crust-core interface. By confronting these models with observations, we illustrate the potential of astrophysical observables to offer constraints on poorly known nuclear matter parameters complementary to terrestrial experiments, and demonstrate that our models consistently predict $L<70$ MeV.

Significant progress has been made in recent years in constraining nuclear symmetry energy at and below the saturation density of nuclear matter using data from both terrestrial nuclear experiments and astrophysical observations. However, many interesting questions remain to be studied especially at supra-saturation densities. In this lecture note, after a brief summary of the currently available constraints on nuclear symmetry energy near the saturation density we first discuss the relationship between the symmetry energy and the isopin and momentum dependence of the single-nucleon potential in isospin-asymmetric nuclear medium. We then discuss several open issues regarding effects of the tensor force induced neutron-proton short-range correlation (SRC) on nuclear symmetry energy. Finally, as an example of the impacts of nuclear symmetry energy on properties of neutron stars and gravitational waves, we illustrate effects of the high-density symmetry energy on the tidal polarizability of neutron stars in coalescing binaries.

How the solar wind is accelerated to its supersonic speed is intimately related to how it is heated. Mechanisms based on ion-cyclotron resonance have been successful in explaining a large number of observations, those concerning the significant ion temperature anisotropy above coronal holes in particular. However, they suffer from the inconsistency with turbulence theory which says that the turbulent cascade in a low-beta medium like the solar corona should proceed in the perpendicular rather than the parallel direction, meaning that there is little energy in the ion gyro-frequency range for ions to absorb via ion-cyclotron resonance. Recently a mechanism based on the interaction between the solar wind particles and the anisotropic turbulence has been proposed, where the perpendicular proton energy addition is via the stochastic heating (Chandran et al. 2011). We extend this promising mechanism by properly accounting for the effect of proton temperature anisotropy on the propagation of Alfven waves, for the radiative losses of electron energy, and for the field line curvature that naturally accompanies solar winds in the corona. While this mechanism was shown in previous studies to apply to the polar fast solar wind, we demonstrate here for the first time that it applies also to the slow wind flowing along field lines bordering streamer helmets.

We investigate the effect of modified gravity on the specific angular momentum of galactic halos by analyzing the halo catalogs at z=0 from high-resolution N-body simulations for a f(R) gravity model that meets the solar-system constraint. It is shown that the galactic halos in the f(R) gravity model tend to acquire significantly higher specific angular momentum than those in the standard LCDM model. The largest difference in the specific angular momentum distribution between these two models occurs for the case of the isolated galactic halos with mass less than 10^{11}M_sun/h, which are likely least shielded by the chameleon screening mechanism. As the specific angular momentum of galactic halos is rather insensitive to the other cosmological parameters, it can in principle be an independent discriminator of modified gravity. We speculate a possibility of using the relative abundance of the low surface brightness galaxies (LSBGs) as a test of GR given that the formation of the LSBGs befalls in the fast spinning dark halos.

Using a set of model equations of state satisfying the latest constraints from both terrestrial nuclear experiments and astrophysical observations as well as state-of-the-art nuclear many-body calculations of the pure neutron matter equation of state, the tidal polarizability of canonical neutron stars in coalescing binaries is found to be a very sensitive probe of the high-density behavior of nuclear symmetry energy which is among the most uncertain properties of dense neutron-rich nucleonic matter. Moreover, it changes less than $\pm 10%$ by varying various properties of symmetric nuclear matter and symmetry energy around the saturation density within their respective ranges of remaining uncertainty.

We study the matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time

In this review, we will discuss the results of our recent work to study the general optimization of the pure isovector parameters of the popular relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) nuclear energy-density functionals (EDFs), using constraints on the pure neutron matter (PNM) equation of state (EoS) from recent {\sl ab initio} calculations. By using RMF and SHF parameterizations that give equivalent predictions for ground-state properties of doubly magic nuclei and properties of symmetric nuclear matter (SNM) and PNM, we found that such optimization leads to broadly consistent symmetry energy $J$ and its slope parameter $L$ at saturation density within a tight range of $\sigma(J) < 2$ MeV and $\sigma(L) < 6$ MeV. We demonstrate that a clear model dependence shows up (a) in the curvature parameter of the symmetry energy $K_{\rm sym}$, (b) the symmetry energy at supra-saturation densities, and (c) the radius of neutron stars.

We study the linear and nonlinear structure formation in the dilaton and symmetron models of modified gravity using a generic parameterisation which describes a large class of scenarios using only a few parameters, such as the coupling between the scalar field and the matter, and the range of the scalar force on very large scales. For this we have modified the N-body simulation code ECOSMOG, which is a variant of RAMSES working in modified gravity scenarios, to perform a set of 110 simulations for different models and parameter values, including the default LCDM. These simulations enable us to explore a large portion of the parameter space. We have studied the effects of modified gravity on the matter power spectrum and mass function, and found a rich and interesting phenomenology where the difference with the LCDM template cannot be reproduced by a linear analysis even on scales as large as k~0.05 h/Mpc. Our results show the full effect of screening on nonlinear structure formation and the associated deviation from LCDM. We also investigate how differences in the force mediated by the scalar field in modified gravity models lead to qualitatively different features for the nonlinear power spectrum and the halo mass function, and how varying the individual model parameters changes these observables. The differences are particularly large in the nonlinear power spectra whose shapes for f(R), dilaton and symmetron models vary greatly, and where the characteristic bump around 1 h/Mpc of f(R) models is preserved for symmetrons, whereas an increase on much smaller scales is particular to symmetrons. No bump is present for dilatons where a flattening of the power spectrum takes place on small scales. These deviations from LCDM and the differences between modified gravity models, such as dilatons and symmetrons, could be tested with future surveys.

We consider the effect of a scalar field degree of freedom on the dynamics of gravity from small to large scales. We show that the effects of modified gravity can be completely captured by the time variations of the scalar field mass and its coupling to matter. This leads to a parameterisation of modified gravity where local constraints are easy to analyse and large scale scale structure effects apparent.