Matter in the Universe is distributed in a cosmic web, with a filament of matter typically connecting neighboring galaxy pairs, separated by tens of millions of light-years. Filaments are known to influence the spins of galaxies around them, but it remains unknown whether filaments themselves spin. Here, we measure the angular momentum around intergalactic filaments in a dark-matter cosmological simulation, finding that matter generally rotates substantially around them. Our results suggest that spin is a fundamental property of filaments, part of a full picture of structure formation. Filaments are likely to be the longest coherently rotating objects in the Universe. This finding helps to explain how galaxies spin in the cosmic web, and could provide a mechanism to generate intergalactic magnetic fields.

The first charge-coupled device photometric light curves in B, V, R${_c}$, and I${_c}$ bands of the short-period contact binary star OQ Cam are presented. Through analyzing the light curves with the Wilson-Devinney code, it is found that OQ Cam is a W-type shallow contact binary star with a contact degree of f =20.2%. Based on the O-C analyzing, the orbit period displays a long term increasing with a rate of $dP/dt=4.40\times10{^{-7}} day \cdot year{^{-1}}$. The increasing in orbit period can be interpreted by mass transfer from the less massive component to the more massive one. As the orbit period increasing, OQ Cam would evolve from the present shallow contact configuration to a none contact stage. So it may be a potential candidate to confirm the thermal relaxation oscillation theory.

To effectively exploit large-scale structure surveys, we depend on accurate and reliable predictions of non-linear cosmological structure formation. Tools for efficient and comprehensive computational modelling are therefore essential to perform cosmological parameter inference analyses. We present the public software package ReACT, demonstrating its capability for the fast and accurate calculation of non-linear power spectra from non-standard physics. We showcase ReACT through a series of forecast constraints on the DGP and $f(R)$ gravity models, adopting LSST-like cosmic shear power spectra. Accurate non-linear modelling with ReACT more than doubles LSST's constraining power on the $f(R)$ parameter, in contrast to an analysis that is limited to the quasi-linear regime. We find that ReACT is sufficiently robust for the inference of consistent constraints on theories beyond $\Lambda$CDM for current and ongoing surveys. With further improvement, particularly in terms of the accuracy of the non-linear $\Lambda$CDM power spectrum, ReACT can, in principle, meet the accuracy requirements for future surveys such as Euclid and LSST.

In this paper, we present the first light curve synthesis and orbital period change analysis of nine contact binaries around the short period limit. It is found that all these systems are W-subtype contact binaries. One of them is a medium contact system while the others are shallow contact ones. Four of them manifest obvious O'Connell effect explained by a dark spot or hot spot on one of the component stars. Third light was detected in three systems. By investigating orbital period variations, we found that four of the targets display a secular period decrease while the others exhibit a long-term period increase. The secular period decrease is more likely caused by angular momentum loss while the long-term period increase is due to mass transfer from the less massive component to the more massive one. Based on the statistic of 19 ultrashort period contact binaries with known orbital period changes, we found that seven of them display long-term decrease (three of them also exhibit cyclic variations), ten of them manifest long-term increase while two of them only show cyclic variation and that most of them are shallow contact binaries supporting the long timescale angular momentum loss theory suggested by Stepien. For the three deep contact systems, we found that they are probably triple systems. The tertiary companion plays an essential role during their formation and evolution.

We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the Kilo-Degree Survey, the Red Cluster Sequence Lensing Survey and the Canada-France-Hawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of $3 - 5\, h^{-1}\text{Mpc}$ is chosen as a proxy for the location of filaments. We measure the average weak lensing shear around $\sim$11,000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we report a $3.4\,\sigma$ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from $N$-body simulations, the average density at the centre of these filamentary structures is found to be $15 \pm 4$ times the critical density.

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian streaming model to calculate the quasi non-linear redshift space matter power spectrum multipoles from a non-linear prescription for the matter power spectrum, which remains agnostic about gravity or dark energy. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model. When comparing the redshift space multipoles for halos, we find that a linear bias implementation is insufficient to be competitive with the TNS-based model used in the BOSS clustering analysis, which however has additional degrees of freedom. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

By analyzing two sets of complete $BVR_cI_c$ light curves of V342 UMa and three sets of complete $BVR_cI_c$ light curves of V509 Cam, we determined that the two systems are both W-subtype contact binaries and that V342 UMa shows a shallow contact configuration, while V509 Cam exhibits a medium contact configuration. Since both of them are totally eclipsing binaries, the physical parameters derived only by the photometric light curves are reliable.Meanwhile, the period changes of the two targets were analyzed based on all available eclipsing times. We discovered that V342 UMa shows long-term period decrease with a rate of $-1.02(\pm0.54)\times10^{-7}$ days/year and that V509 Cam exhibits long-term period increase with a rate of $3.96(\pm0.90)\times10^{-8}$ days/year. Both the conservative mass transfer and AML via magnetic stellar winds can interpret the long-term period decrease of V342 UMa. The long-term period increase of V509 Cam can be explained by mass transfer from the less massive star to the more massive one. The absolute parameters of the two binaries were estimated according to the Gaia distances and our derived photometric solution results. This method can be extended to other contact binaries without radial velocities but with reliable photometric solutions. The evolutionary states of them were discussed, we found that they reveal identical properties of other W-subtype contact systems.

The period distribution of contact binaries exhibits a very sharp short period cut-off at 0.22 days. In order to provide valuable information on this short period limit, we observed ten totally eclipsing contact binaries with orbital periods near this cut-off. By detailed analysis using W-D code, we determined that two of these systems are A-subtype contact binaries while the others are W-subtype contact binaries and all the targets show shallow contact configurations. Half of the targets exhibit stellar spot activity and four of them have third light. A statistical work on well studied USPCBs was carried out, and physical parameters of fifty-five USPCBs were obtained. Some common properties for these systems were derived. Due to the study of the period-color diagram of USPCBs, we found that the period-color relation of USPCBs is different from other W UMa type contact binaries and USPCBs are metal poor old stars. In addition, the evolutionary states of these systems were discussed by constructing the color-density diagram. We derived that the evolutionary states of the two components of USPCBs show identical characteristics of other contact binaries despite slower evolutionary status caused by smaller mass. We suggested that both the fully convective limit claimed by Rucinski (1992) and the AML theory recommended by Stepien (2006) can produce the short period cut-off and a tertiary companion is very important during the formation of the short period contact binaries by removing angular momentum from the host eclipsing pairs under certain circumstances.

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

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a subhalo abundance matching technique coupled to merger trees extracted from a N-body simulation. Our approach consist of forcing the model to match the observed SMF at redshift $z \sim 2.3$, and let it evolve down to $z \sim 0.3$ according to a $\tau$ model, an exponentially declining functional form which describes the star formation rate decay of both satellite and central galaxies. In this study, we use three different sets of SMFs: ZFOURGE data from Tomczak et al.; UltraVISTA data from Ilbert et al. and COSMOS data from Davidzon et al. We also build a mock survey combining UltraVISTA with ZFOURGE. Our modelling of quenching timescales is consistent with the evolution of the SMF down to $z \sim 0.3$, with different accuracy depending on the particular survey used for calibration. We tested our model against the observed SMFs at low redshift and it predicts residuals (observation versus model) within $1\sigma$ observed scatter along most of the stellar mass range investigated, and with mean residuals below 0.1 dex in the range $\sim [10^{8.7}-10^{11.7}] M_{\odot}$. We then compare the SFR-$M_*$ relation predicted by the model with the observed one at different redshifts. The predicted SFR-$M_*$ relation underpredicts the median SFR at fixed stellar mass relative to observations at all redshifts. Nevertheless, the shapes are consistent with the observed relations up to intermediate-mass galaxies, followed by a rapid decline for massive galaxies.

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

The radial number density and flattening of the Milky Way's stellar halo is measured with $\mathrm{5351}$ metal-poor ([Fe/H]$<-1$) K giants from LAMOST DR3, using a nonparametric method which is model independent and largely avoids the influence of halo substucture. The number density profile is well described by a single power law with index $5.03^{+0.64}_{-0.64}$, and flattening that varies with radius. The stellar halo traced by LAMOST K giants is more flattened at smaller radii, and becomes nearly spherical at larger radii. The flattening, $q$, is about 0.64, 0.8, 0.96 at $r=15$, 20 and 30 kpc (where $r=\sqrt{R^2+\left[Z/q\left(r\right)\right]^2}$), respectively. Moreover, the leading arm of the Sagittarius dwarf galaxy tidal stream in the north, and the trailing arm in the south, are significant in the residual map of density distribution. In addition, an unknown overdensity is identified in the residual map at (R,Z)=(30,15) kpc.

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a Subhalo Abundance Matching technique coupled to merger trees extracted from a N-body simulation. Our approach, which considers both galaxy mergers and stellar stripping, is to force the model to match the observed SMF at redshift $z>2$, and let it evolve down to the present time according to the observed (SFR-$M_*$) relation. In this study, we use two different sets of SMFs and two SFR-$M_*$ relations: a simple power law and a relation with a mass-dependent slope. Our analysis shows that the evolution of the SMF is more consistent with a SFR-$M_*$ relation with a mass-dependent slope, in agreement with predictions from other models of galaxy evolution and recent observations. In order to fully and realistically describe the evolution of the SMF, both mergers and stellar stripping must be considered, and we find that both have almost equal effects on the evolution of SMF at the massive end. Taking into account the systematic uncertainties in the observed data, the high-mass end of the SMF obtained by considering stellar stripping results in good agreement with recent observational data from the Sloan Digital Sky Survey (SDSS). At $\log M_* < 11.2$, our prediction at z=0.1 is close to \citet{li-white09} data, but the high-mass end ($\log M_* > 11.2$) is in better agreement with \citet{dsouza15} data which account for more massive galaxies.

Measurement of the local dark matter density plays an important role in both Galactic dynamics and dark matter direct detection experiments. However, the estimated values from previous works are far from agreeing with each other. In this work, we provide a well-defined observed sample with 1427 G \& K type main-sequence stars from the LAMOST spectroscopic survey, taking into account selection effects, volume completeness, and the stellar populations. We apply a vertical Jeans equation method containing a single exponential stellar disk, a razor thin gas disk, and a constant dark matter density distribution to the sample, and obtain a total surface mass density of $\rm {78.7 ^{+3.9}_{-4.7}\ M_{\odot}\ pc^{-2}}$ up to 1 kpc and a local dark matter density of $0.0159^{+0.0047}_{-0.0057}\,\rm M_{\odot}\,\rm pc^{-3}$. We find that the sampling density (i.e. number of stars per unit volume) of the spectroscopic data contributes to about two-thirds of the uncertainty in the estimated values. We discuss the effect of the tilt term in the Jeans equation and find it has little impact on our measurement. Other issues, such as a non-equilibrium component due to perturbations and contamination by the thick disk population, are also discussed.

Recent observations reveal that magnetic turbulence in the nearly colisionless solar wind plasma extends to scales smaller than the plasma microscales, such as ion gyroradius and ion inertial length. Measured breaks in the spectra of magnetic and density fluctuations at high frequencies are thought to be related to the transition from large-scale hydromagnetic to small-scale kinetic turbulence. The scales of such transitions and the responsible physical mechanisms are not well understood however. In the present work we emphasize the crucial role of the plasma parameters in the transition to kinetic turbulence, such as the ion and electron plasma beta, the electron to ion temperature ratio, the degree of obliquity of turbulent fluctuations. We then propose an explanation for the spectral breaks reported in recent observations.

We use about 15,000 F/G nearby dwarf stars selected from the LAMOST pilot survey to map the U-V velocity distribution in the solar neighbourhood. An extreme deconvolution algorithm is applied to reconstruct an empirical multi-Gaussian model. In addition to the well known substructures, e.g., Sirius, Coma Berenices, Hyades-Pleiades over-densities, several new substructures are unveiled. A ripple-like structure from (U, V) = (-120, -5) to (103, -32)km/s is clearly seen in the U-V distribution. This structure seems associated with resonance induced by the Galactic bar, since it is extended in U while having a small dispersion in V at the same time. A ridge structure between (U, V) = (-60, 40) and (-15, 15) km/s is also found. Although similar substructures have been seen in the Hipparcos data, their origin is still unclear. Another compact over-density is seen at (U, V) = (-102, -24). With this large data sample, we find that the substructure located at V~70 km/s and the Arcturus group are essentially parallel in V, which may indicate that they originate from an unrelaxed disk component perturbed by the rotating bar.

Recent theoretical studies argue that the rate of stochastic ion heating in low-frequency Alfv\'en-wave turbulence is given by $Q_\perp = c_1 [(\delta u)^3 /\rho] \exp(-c_2/\epsilon)$, where $\delta u$ is the rms turbulent velocity at the scale of the ion gyroradius $\rho$, $\epsilon = \delta u/v_{\perp \rm i}$, $v_{\perp \rm i}$ is the perpendicular ion thermal speed, and $c_1$ and $c_2$ are dimensionless constants. We test this theoretical result by numerically simulating test particles interacting with strong reduced magnetohydrodynamic (RMHD) turbulence. The heating rates in our simulations are well fit by this formula. The best-fit values of $c_1$ are $\sim 1$. The best-fit values of $c_2$ decrease (i.e., stochastic heating becomes more effective) as the grid size and Reynolds number of the RMHD simulations increase. As an example, in a $1024^2 \times 256$ RMHD simulation with a dissipation wavenumber of order the inverse ion gyroradius, we find $c_2 = 0.21$. We show that stochastic heating is significantly stronger in strong RMHD turbulence than in a field of randomly phased Alfv\'en waves with the same power spectrum, because coherent structures in strong RMHD turbulence increase orbit stochasticity in the regions where ions are heated most strongly. We find that $c_1$ increases by a factor of $\sim 3$ while $c_2$ changes very little as the ion thermal speed increases from values $\ll v_{\rm A}$ to values $\sim v_{\rm A}$, where $v_{\rm A}$ is the Alfv\'en speed. We discuss the importance of these results for perpendicular ion heating in the solar wind.

The nature of subproton scale fluctuations in the solar wind is an open question, partly because two similar types of electromagnetic turbulence can occur: kinetic Alfven turbulence and whistler turbulence. These two possibilities, however, have one key qualitative difference: whistler turbulence, unlike kinetic Alfven turbulence, has negligible power in density fluctuations. In this Letter, we present new observational data, as well as analytical and numerical results, to investigate this difference. The results show, for the first time, that the fluctuations well below the proton scale are predominantly kinetic Alfven turbulence, and, if present at all, the whistler fluctuations make up only a small fraction of the total energy.

One proposed mechanism for heating the solar wind, from close to the sun to beyond 10 AU, invokes low-frequency, oblique, Alfven-wave turbulence. Because small-scale oblique Alfven waves (kinetic Alfven waves) are compressive, the measured density fluctuations in the solar wind place an upper limit on the amplitude of kinetic Alfven waves and hence an upper limit on the rate at which the solar wind can be heated by low-frequency, Alfvenic turbulence. We evaluate this upper limit for both coronal holes at 5 solar radii and in the near-Earth solar wind. At both radii, the upper limit we find is consistent with models in which the solar wind is heated by low-frequency Alfvenic turbulence. At 1 AU, the upper limit on the turbulent heating rate derived from the measured density fluctuations is within a factor of 2 of the measured solar wind heating rate. Thus if low-frequency Alfvenic turbulence contributes to heating the near-Earth solar wind, kinetic Alfven waves must be one of the dominant sources of solar wind density fluctuations at frequencies of order 1 Hz. We also present a simple argument for why density fluctuation measurements do appear to rule out models in which the solar wind is heated by non-turbulent high-frequency waves ``sweeping'' through the ion-cyclotron resonance, but are compatible with heating by low-frequency Alfvenic turbulence.