The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

This is the second paper in a series where we propose a method of indirectly measuring large scale structure using information from small scale perturbations. The idea is to build a quadratic estimator from small scale modes that provides a map of structure on large scales. We demonstrated in the first paper that the quadratic estimator works well on a dark-matter-only N-body simulation at a snapshot of $z = 0$. Here we generalize the theory to the case of a light-cone with halo bias and redshift space distortions taken into consideration. We successfully apply the generalized version of the quadratic estimator to a light cone halo catalog of an N-body simulation of size $\sim5.6\,(h^{-1}\,\rm Gpc)^3$. The most distant point in the light cone is at a redshift of $1.4$, which indicates that we might be able to apply our method to next generation galaxy surveys.

While magnetic fields are important in contemporary star formation, their role in primordial star formation is unknown. Magnetic fields of order 10^-16 G are produced by the Biermann battery due to the curved shocks and turbulence associated with the infall of gas into the dark matter minihalos that are the sites of formation of the first stars. These fields are rapidly amplified by a small-scale dynamo until they saturate at or near equipartition with the turbulence in the central region of the gas. Analytic results are given for the outcome of the dynamo, including the effect of compression in the collapsing gas. The mass-to-flux ratio in this gas is 2-3 times the critical value, comparable to that in contemporary star formation. Predictions of the outcomes of simulations using smooth particle hydrodynamics (SPH) and grid-based adaptive mesh refinement (AMR) are given. Because the numerical viscosity and resistivity for the standard resolution of 64 cells per Jeans length are several orders of magnitude greater than the physical values, dynamically significant magnetic fields affect a much smaller fraction of the mass in simulations than in reality. An appendix gives an analytic treatment of free-fall collapse, including that in a constant density background. Another appendix presents a new method of estimating the numerical viscosity; results are given for both SPH and grid-based codes.

We consider the high spin expansion for the null geodesics in the Kerr spacetime. We expand the null geodesic equation successively to higher orders in deviation from extremity. Via the method of matched asymptotic expansion, the radial integrals are obtained analytically. It turns out that the analytic expressions are very sensitive to the value of the shifted Carter constant $q$. We show that for a large $q$, the analytic expressions can be used to study observational electromagnetic signatures for astrophysical black holes like M87*. However, for a small $q$, the high spin expansion method can only be applied to (near-) extreme black holes.

We explore the prospects of detecting of Galactic double white dwarf (DWD) binaries with the space-based gravitational wave (GW) observatory TianQin. In this work we analyse both a sample of currently known and a realistic synthetic population of DWDs to assess the number of guaranteed detections and the full capacity of the mission. We find that TianQin can detect 12 out of $\sim100$ known DWDs; GW signals of these binaries can be modelled in detail ahead of the mission launch, and therefore they can be used as verification sources. Besides we estimate that TianQin has potential to detect as many as $10^4$ DWDs in the Milky Way. TianQin is expected to measure their orbital periods and amplitudes with accuracy of $\sim10^{-7}$ and $\sim0.2$ respectively, and to localize on the sky a large fraction (39%) of the detected population to better than 1 deg$^2$. We conclude that TianQin has the potential to significantly advance our knowledge on Galactic DWDs by increasing the sample up to 2 orders of magnitude, and will allow their multi-messenger studies in combination with electromagnetic telescopes. We also test the possibilities of different configuration of TianQin: 1) the same mission with a different orientation, 2) two perpendicular constellations combined into a network, and 3) the combination of the network with the ESA-lead Laser Interferometer Space Antenna. We find that the network of detectors boosts the accuracy on the measurement of source parameters by $1-2$ orders of magnitude, with the improvement on sky localization being the most significant.

A clear transition of galaxy quenching is identified in the multi-parameter space of stellar mass ($M_*$), bulge to total mass ratio ($B/T_{\rm m}$), halo mass ($M_{\rm h}$) and halo-centric distance ($r/r_{180}$). For given halo mass, the characteristic stellar mass ($M_{*, \rm ch}$) for the transition is about one-fifth of that of the corresponding central galaxy, and almost independent of $B/T_{\rm m}$. Once $B/T_{\rm m}$ is fixed, the quenched fraction of galaxies with $M_*<M_{*, \rm ch}$ increases with $M_{\rm h}$, but decreases with $M_*$ in the inner part of halos ($r/r_{180}<0.5$). In the outer part ($r/r_{180}>0.5$), the trend with $M_{\rm h}$ remains but the correlation with $M_*$ is absent or becomes positive. For galaxies above $M_{\rm *, ch}$ and with $B/T_{\rm m}$ fixed, the quenched fraction increases with $M_{\rm *}$, but depends only weakly on $M_{\rm h}$ in both the inner and outer regions. Moreover, at fixed $B/T_{\rm m}$ and $M_*$, the quenched fraction increases with the decreases of $r/r_{180}$ for galaxies with $M_*<M_{*, \rm ch}$, while remaining constant with $r/r_{180}$ for galaxies with $M_*>M_{*, \rm ch}$. Our finding provides a physically-motivated way to classify galaxies in halos into two classes based on their quenching properties: an `upper class' with $M_*>M_{\rm *,ch}$ and a `lower class' with $M_*<M_{\rm *,ch}$. Environmental quenching is important only for `lower class' galaxies, while internal quenching plays the dominating role for the `upper class'.

Since the first detection in 2015, ten black hole binary merger events have been formally reported. These gravitational wave signals have significantly enhanced our understanding of the black hole astrophysics. In general, the properties of a black hole are comprehensively described by its mass, spin, and charge. The third parameter is often ignored because of the very low value expected in the realistic astrophysical environment. In this work, we constrain the amount of charge in a way which is equivalent to the parametrized post-Einsteinian framework by treating the charge as a small perturbation with Bayesian method, and we find that the current limit on the charge is similar to the one obtained by the Fisher information matrix method in previous works. Then we develop a zeroth order post-Newtonian waveform for charged binary black hole inspirals, and apply this charged waveform to the binary black hole merger events observed by LIGO-Virgo in their first two runs. A Bayesian model selection is performed among the post-Newtonian waveform, charged post-Newtonian waveform, and full inspiral-merger-ringdown waveform. Remarkably, we show that charged post-Newtonian waveform matches the LIGO-Virgo data better compared with non-charged post-Newtonian waveform, mainly due to the influence in the late inspiral period, and the full inspiral-merger-ringdown outrun both by far. The inclusion of charge in the waveform might shift the estimated parameters significantly. Finally, we obtain certain constraints on the binary black hole charges and discuss the drawbacks. Our work calls for the development of a full inspiral-merger-ringdown waveform for charged binary black hole merger, so that more meaningful constraints could be obtained.

Large scale density modes are difficult to measure because they are sensitive to systematic observational errors in galaxy surveys, but we can study them indirectly by observing their impact on small scale perturbations. Cosmological perturbation theory predicts that second-order density inhomogeneities are a convolution of a short- and a long-wavelength mode. This arises physically because small scale structures grow at different rates depending on the large scale environment in which they reside. This induces an off-diagonal term in the two-point statistics in Fourier space that we use as the basis for a quadratic estimator for the large scale field. We demonstrate that this quadratic estimator works well on an N-body simulation of size (2.5 h^{-1} Gpc)^3. In particular, the quadratic estimator successfully reconstructs the long-wavelength modes using only small-scale information. This opens up novel opportunities to study structure on the largest observable scales.

We present rotation curve fits to 175 late-type galaxies from the Spitzer Photometry & Accurate Rotation Curves (SPARC) database using seven dark matter (DM) halo profiles: pseudo-isothermal (pISO), Burkert, Navarro-Frenk-White (NFW), Einasto, Di Cintio (2014, DC14), coreNFW, and a new semi-empirical profile named Lucky13. We marginalize over stellar mass-to-light ratio, galaxy distance, disk inclination, halo concentration and halo mass (and an additional shape parameter for Einasto) using a Markov Chain Monte Carlo method. We find that cored halo models such as the DC14 and Burkert profiles generally provide better fits to rotation curves than the cuspy NFW profile. The stellar mass-halo mass relation from abundance matching is recovered by all halo profiles once imposed as a Bayesian prior, whereas the halo mass-concentration relation is not reproduced in detail by any halo model. We provide an extensive set of figures as well as best-fit parameters in machine-readable tables to facilitate model comparison and the exploration of DM halo properties.

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $\mu m$ continuum emissions. The angles between outflow axes and corresponding filaments ($\gamma_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large scale (>~$10^3$ AU) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of $\mu_\phi$=2, 8, and 32. The initial separations of protostars in multiples depends on the global magnetic field strength, with stronger magnetic fields (e.g., $\mu_\phi$=2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4-0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin mis-alignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to <~ $10^2$ AU separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.

Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.

We present an empirical method to measure the halo mass function (HMF) of galaxies. We determine the relation between the \hi\ line-width from single-dish observations and the dark matter halo mass ($M_{200}$) inferred from rotation curve fits in the SPARC database, then we apply this relation to galaxies from the \hi\ Parkes All Sky Survey (HIPASS) to derive the HMF. This empirical HMF is well fit by a Schecther function, and matches that expected in $\Lambda$CDM over the range $10^{10.5} < M_{200} < 10^{12}\;\mathrm{M}_{\odot}$. More massive halos must be poor in neutral gas to maintain consistency with the power law predicted by $\Lambda$CDM. We detect no discrepancy at low masses. The lowest halo mass probed by HIPASS, however, is just greater than the mass scale where the Local Group missing satellite problem sets in. The integrated mass density associated with the dark matter halos of \hi-detected galaxies sums to $\Omega_{\rm m,gal} \approx 0.03$ over the probed mass range.

Galaxies are observed to obey a strict set of dynamical scaling relations. We review these relations for rotationally supported disk galaxies spanning many decades in mass, surface brightness, and gas content. The behavior of these widely varied systems can be summarized with a handful of empirical laws connected by a common acceleration scale.

We present the B-fields mapped in IRDC G34.43+0.24 using 850\,$\mu$m polarized dust emission observed with the POL-2 instrument at JCMT. We examine the magnetic field geometries and strengths in the northern, central, and southern regions of the filament. The overall field geometry is ordered and aligned closely perpendicular to the filament's main axis, particularly in regions containing the central clumps MM1 and MM2, whereas MM3 in the north has field orientations aligned with its major axis. The overall field orientations are uniform at large (POL-2 at 14$\arcsec$ and SHARP at 10$\arcsec$) to small scales (TADPOL at 2.5$\arcsec$ and SMA at 1.5$\arcsec$) in the MM1 and MM2 regions. SHARP/CSO observations in MM3 at 350\,$\mu$m from Tang et al. show a similar trend as seen in our POL-2 observations. TADPOL observations demonstrate a well-defined field geometry in MM1/MM2 consistent with MHD simulations of accreting filaments. We obtained a plane-of-sky magnetic field strength of 470$\pm$190\,$\mu$G, 100$\pm$40\,$\mu$G, and 60$\pm$34\,$\mu$G in the central, northern and southern regions of G34, respectively, using the updated Davis-Chandrasekhar-Fermi relation. The estimated value of field strength, combined with column density and velocity dispersion values available in the literature, suggests G34 to be marginally critical with criticality parameter $\rm \lambda$ values 0.8$\pm$0.4, 1.1$\pm$0.8, and 0.9$\pm$0.5 in the central, northern, and southern regions, respectively. The turbulent motions in G34 are sub-Alfv\'{e}nic with Alfv\'{e}nic Mach numbers of 0.34$\pm$0.13, 0.53$\pm$0.30, and 0.49$\pm$0.26 in the three regions. The observed aligned B-fields in G34.43+0.24 are consistent with theoretical models suggesting that B-fields play an important role in guiding the contraction of the cloud driven by gravity.

The surface of last scattering of the photons in the cosmic microwave background is not a spherical shell. Apart from its finite width, each photon experiences a different gravitational potential along its journey to us, leading to different travel times in different directions. Since all photons were released at the same cosmic time, the photons with longer travel times started farther away from us than those with shorter times. Thus, the surface of last scattering is corrugated, a deformed spherical shell. We present an estimator quadratic in the temperature and polarization fields that could provide a map of the time delays as a function of position on the sky. The signal to noise of this map could exceed unity for the dipole, thereby providing a rare insight into the universe on the largest observable scales.

Context. Radio occultation is a technique used to study planetary atmospheres by means of the refraction and absorption of a spacecraft carrier signal through the atmosphere of the celestial body of interest, as detected from a ground station on Earth. This technique is usually employed by the deep space tracking and communication facilities (e.g., NASA's Deep Space Network (DSN), ESA's Estrack). Aims. We want to characterize the capabilities of the Planetary Radio Interferometry and Doppler Experiment (PRIDE) technique for radio occultation experiments, using radio telescopes equipped with Very Long Baseline Interferometry (VLBI) instrumentation. Methods. We conducted a test with ESA's Venus Express (VEX), to evaluate the performance of the PRIDE technique for this particular application. We explain in detail the data processing pipeline of radio occultation experiments with PRIDE, based on the collection of so-called open-loop Doppler data with VLBI stations, and perform an error propagation analysis of the technique. Results. With the VEX test case and the corresponding error analysis, we have demonstrated that the PRIDE setup and processing pipeline is suited for radio occultation experiments of planetary bodies. The noise budget of the open-loop Doppler data collected with PRIDE indicated that the uncertainties in the derived density and temperature profiles remain within the range of uncertainties reported in previous Venus' studies. Open-loop Doppler data can probe deeper layers of thick atmospheres, such as that of Venus, when compared to closed-loop Doppler data. Furthermore, PRIDE through the VLBI networks around the world, provides a wide coverage and range of large antenna dishes, that can be used for this type of experiments.

We perform ideal MHD high resolution AMR simulations with driven turbulence and self-gravity and find that long filamentary molecular clouds are formed at the converging locations of large-scale turbulence flows and the filaments are bounded by gravity. The magnetic field helps shape and reinforce the long filamentary structures. The main filamentary cloud has a length of ~4.4 pc. Instead of a monolithic cylindrical structure, the main cloud is shown to be a collection of fiber/web-like sub-structures similar to filamentary clouds such as L1495. Unless the line-of-sight is close to the mean field direction, the large-scale magnetic field and striations in the simulation are found roughly perpendicular to the long axis of the main cloud, similar to 1495. This provides strong support for a large-scale moderately strong magnetic field surrounding L1495. We find that the projection effect from observations can lead to incorrect interpretations of the true three-dimensional physical shape, size, and velocity structure of the clouds. Helical magnetic field structures found around filamentary clouds that are interpreted from Zeeman observations can be explained by a simple bending of the magnetic field that pierces through the cloud. We demonstrate that two dark clouds form a T-shape configuration which are strikingly similar to the Infrared dark cloud SDC13 leading to the interpretation that SDC13 results from a collision of two long filamentary clouds. We show that a moderately strong magnetic field (M_A ~ 1) is crucial for maintaining a long and slender filamentary cloud for a long period of time ~0.5 million years.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $\mu$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex.

We study the scaling relations between dark matter (DM) haloes and galaxy discs using 175 galaxies from the SPARC database. We explore two cosmologically motivated DM halo profiles: the Einasto profile from DM-only simulations and the DC14 profile from hydrodynamic simulations. We fit the observed rotation curves using a Markov Chain Monte Carlo method and break the disc-halo degeneracy using near-infrared photometry and $\Lambda$CDM-motivated priors. We find that the characteristic volume density $\rho_{\rm s}$ of DM haloes is nearly constant over $\sim$5 decades in galaxy luminosity. The scale radius $r_s$ and the characteristic surface density $\rho_s\cdot r_s$, instead, correlate with galaxy luminosity. These scaling relations provide an empirical benchmark to cosmological simulations of galaxy formation.

We constrain the mass--richness scaling relation of redMaPPer galaxy clusters identified in the Dark Energy Survey Year 1 data using weak gravitational lensing. We split clusters into $4\times3$ bins of richness $\lambda$ and redshift $z$ for $\lambda\geq20$ and $0.2 \leq z \leq 0.65$ and measure the mean masses of these bins using their stacked weak lensing signal. By modeling the scaling relation as $\langle M_{\rm 200m}|\lambda,z\rangle = M_0 (\lambda/40)^F ((1+z)/1.35)^G$, we constrain the normalization of the scaling relation at the 5.0 per cent level as $M_0 = [3.081 \pm 0.075 ({\rm stat}) \pm 0.133 ({\rm sys})] \cdot 10^{14}\ {\rm M}_\odot$ at $\lambda=40$ and $z=0.35$. The richness scaling index is constrained to be $F=1.356 \pm 0.051\ ({\rm stat})\pm 0.008\ ({\rm sys})$ and the redshift scaling index $G=-0.30\pm 0.30\ ({\rm stat})\pm 0.06\ ({\rm sys})$. These are the tightest measurements of the normalization and richness scaling index made to date. We use a semi-analytic covariance matrix to characterize the statistical errors in the recovered weak lensing profiles. Our analysis accounts for the following sources of systematic error: shear and photometric redshift errors, cluster miscentering, cluster member dilution of the source sample, systematic uncertainties in the modeling of the halo--mass correlation function, halo triaxiality, and projection effects. We discuss prospects for reducing this systematic error budget, which dominates the uncertainty on $M_0$. Our result is in excellent agreement with, but has significantly smaller uncertainties than, previous measurements in the literature, and augurs well for the power of the DES cluster survey as a tool for precision cosmology and upcoming galaxy surveys such as LSST, Euclid and WFIRST.

There is undetermined potential function $V(\phi)$ in the action of mimetic gravity which should be resolved through physical means. In general relativity(GR), the static spherically symmetric(SSS) solution to the Einstein equation is a benchmark and its deformation also plays a crucial role in mimetic gravity. The equation of motion is provided with high nonlinearity, but we can reduce primal nonlinearity to a frequent Riccati form in the SSS case of mimetic gravity. In other words, we obtain an expression of solution to the functional differential equation of motion with any potential function. Remarkably, we proved rigorously that there is a zero point of first order for the metric function $\beta(r)$ if another metric function $\alpha(r)$ possesses a pole of first order within mimetic gravity. The zero point theorem may be regarded as the generalization of Birkhoff theorem $\alpha\beta=1$ in GR. As a corollary, we show that there is a modified black hole solution for any given $V(\phi)$, which can pass the test of solar system. As another corollary, the zero point theorem provides a dynamical mechanism for the maximum size of galaxies. Especially, there are two analytic solutions which provide good fits to the rotation curves of galaxies without the demand for particle dark matter.

Galaxies follow a tight radial acceleration relation (RAR): the acceleration observed at every radius correlates with that expected from the distribution of baryons. We use the Markov Chain Monte Carlo method to fit the mean RAR to 175 individual galaxies in the SPARC database, marginalizing over stellar mass-to-light ratio ($\Upsilon_{\star}$), galaxy distance, and disk inclination. Acceptable fits with astrophysically reasonable parameters are found for the vast majority of galaxies. The residuals around these fits have an rms scatter of only 0.057 dex ($\sim$13$\%$). This is in agreement with the predictions of modified Newtonian dynamics (MOND). We further consider a generalized version of the RAR that, unlike MOND, permits galaxy-to-galaxy variation in the critical acceleration scale. The fits are not improved with this additional freedom: there is no credible indication of variation in the critical acceleration scale. The data are consistent with the action of a single effective force law. The apparent universality of the acceleration scale and the small residual scatter are key to understanding galaxies.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 \textit{Planck} Galactic Cold Clumps (PGCCs) in the second quadrant of the Milky Way. Using the $^{13}$CO and C$^{18}$O $J = 1 - 0$ lines, and 850\,$\mu$m continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the $^{13}$CO line and 850\,$\mu$m continuum maps, respectively. We make use of the Bayesian Distance Calculator and derived the distances of all 64 PGCCs. We found that in general, the mass-size plane follows a relation of $m\sim r^{1.67}$. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass-size relation suggests that the PGCCs in our sample have a low core formation efficiency ($\sim$3.0\%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850\,$\mu$m cores are more turbulent, more optically thick, and denser than the $^{13}$CO clumps for star formation candidates, suggesting that the 850\,$\mu$m cores are likely more appropriate future star-formation candidates than the $^{13}$CO clumps.

Stellar feedback from high-mass stars (e.g., H{\sc ii} regions) can strongly influence the surrounding interstellar medium and regulate star formation. Our new ALMA observations reveal sequential high-mass star formation taking place within one sub-virial filamentary clump (the G9.62 clump) in the G9.62+0.19 complex. The 12 dense cores (MM 1-12) detected by ALMA are at very different evolutionary stages, from starless core phase to UC H{\sc ii} region phase. Three dense cores (MM6, MM7/G, MM8/F) are associated with outflows. The mass-velocity diagrams of outflows associated with MM7/G and MM8/F can be well fitted with broken power laws. The mass-velocity diagram of SiO outflow associated with MM8/F breaks much earlier than other outflow tracers (e.g., CO, SO, CS, HCN), suggesting that SiO traces newly shocked gas, while the other molecular lines (e.g., CO, SO, CS, HCN) mainly trace the ambient gas continuously entrained by outflow jets. Five cores (MM1, MM3, MM5, MM9, MM10) are massive starless core candidates whose masses are estimated to be larger than 25 M$_{\sun}$, assuming a dust temperature of $\leq$ 20 K. The shocks from the expanding H{\sc ii} regions ("B" \& "C") to the west may have great impact on the G9.62 clump through compressing it into a filament and inducing core collapse successively, leading to sequential star formation. Our findings suggest that stellar feedback from H{\sc ii} regions may enhance the star formation efficiency and suppress the low-mass star formation in adjacent pre-existing massive clumps.

Star formation in a filamentary infrared dark cloud (IRDC) is simulated over a dynamic range of 4.2 pc to 28 au for a period of $3.5\times 10^5$ yr, including magnetic fields and both radiative and outflow feedback from the protostars. At the end of the simulation, the star formation efficiency is 4.3 per cent and the star formation rate per free fall time is $\epsilon_{\rm ff}\simeq 0.04$, within the range of observed values (Krumholz et al. 2012a). The total stellar mass increases as $\sim\,t^2$, whereas the number of protostars increases as $\sim\,t^{1.5}$. We find that the density profile around most of the simulated protostars is $\sim\,\rho\propto r^{-1.5}$, as predicted by Murray & Chang (2015). At the end of the simulation, the protostellar mass function approaches the Chabrier (2005) stellar initial mass function. We infer that the time to form a star of median mass $0.2\,M_\odot$ is about $1.4\times 10^5$~yr from the median mass accretion rate. We find good agreement among the protostellar luminosities observed in the large sample of Dunham et al. (2013), our simulation, and a theoretical estimate, and conclude that the classical protostellar luminosity problem Kenyon et al. (1990) is resolved. The multiplicity of the stellar systems in the simulation agrees to within a factor 2 of observations of Class I young stellar objects; most of the simulated multiple systems are unbound. Bipolar protostellar outflows are launched using a sub-grid model, and extend up to 1 pc from their host star. The mass-velocity relation of the simulated outflows is consistent with both observation and theory.

In this paper, we perform a cosmological model-independent test of the cosmic distance-duality relation (CDDR) in terms of the ratio of angular diameter distance (ADD) $D=D_{\rm A}^{\rm sl}/D_{\rm A}^{\,\rm s}$ from strong gravitational lensing (SGL) and the ratio of luminosity distance (LD) $D^\ast=D_{\rm L}^{\,\rm l}/D_{\rm L}^{\,\rm s}$ obtained from the joint of type Ia supernovae (SNIa) Union2.1 compilation and the latest Gamma-Ray Bursts (GRBs) data, where the superscripts s and l correspond to the redshifts $z_{\,\rm s}$ and $z_{\,\rm l}$ at the source and lens from SGL samples. The purpose of combining GRB data with SNIa compilation is to test CDDR in a wider redshift range. The LD associated with the redshits of the observed ADD, is obtained through two cosmological model-independent methods, namely, method A: binning the SNIa+GRBs data, and method B: reconstructing the function of DL by combining the Crossing Statistic with the smoothing method. We find that CDDR is compatible with the observations at $1\sigma$ confidence level for the power law model which is assumed to describe the mass distribution of lensing systems with method B in a wider redshift range.

We present a detailed study of the Vaidya solution and its generalization in de Rham-Gabadadze-Tolley (dRGT) theory. Since the diffeomorphism invariance can be restored with the St\"{u}ckelberg fields $\phi^a$ introduced, there is a new invariant $I^{ab}=g^{\mu \nu}\partial_\mu \phi^a\partial_\nu \phi^b$ in the massive gravity, which adds to the ones usually encountered in general relativity. There is no conventional Vaidya solution if we choose unitary gauge. In this paper, we obtain three types of self-consistent ansatz with some nonunitary gauge, and find accordingly the Vaidya, generalized Vaidya and furry Vaidya solution. As by-products, we obtain a series of furry black hole. The Vaidya solution and its generalization in dRGT massive gravity describe the black holes with a variable horizon.

We present a detailed study of the static spherically symmetric solutions in de Rham-Gabadadze-Tolley (dRGT) theory. Since the diffeomorphism invariance can be restored by introducing the St\"{u}ckelberg fields $\phi^a$, there is new invariant $I^{ab}=g^{\mu\nu}\partial_{\mu}\phi^a\partial_\nu\phi^b$ in the massive gravity, which adds to the ones usually encountered in general relativity (GR). In the unitary gauge $\phi^a=x^\mu\delta_\mu^a$, any inverse metric $g^{\mu\nu}$ that has divergence including the coordinate singularity in GR would exhibit a singularity in the invariant $I^{ab}$. Therefore, there is no conventional Schwarzschild metric if we choose unitary gauge. In this paper, we obtain a self-consistent static spherically symmetric ansatz in the nonunitary gauge. Under this ansatz, we find that there are seven solutions including the Schwarzschild solution, Reissner-Nordstr\"{o}m solution and five other solutions. These solutions may possess an event horizon depending upon the physical parameters (Schwarzschild radius $r_s$, scalar charge $S$ and/or electric charge $Q$). If these solutions possess an event horizon, we show that the singularity of $I^{ab}$ is absent at the horizon. Therefore, these solutions may become candidates for black holes in dRGT.

The intermittent dissipation of interstellar turbulence is an important energy source in the diffuse ISM. Though on average smaller than the heating rates due to cosmic rays and the photoelectric effect on dust grains, the turbulent cascade can channel large amounts of energy into a relatively small fraction of the gas that consequently undergoes significant heating and chemical enrichment. In particular, this mechanism has been proposed as a solution to the long-standing problem of the high abundance of CH+ along diffuse molecular sight lines, which steady-state, low temperature models under-produce by over an order of magnitude. While much work has been done on the structure and chemistry of these small-scale dissipation zones, comparatively little attention has been paid to relating these zones to the properties of the large-scale turbulence. In this paper, we attempt to bridge this gap by estimating the temperature and CH+ column density along diffuse molecular sight-lines by post-processing 3-dimensional MHD turbulence simulations. Assuming reasonable values for the cloud density (30 / cm^3), size (20 pc), and velocity dispersion (2.3 km / s), we find that our computed abundances compare well with CH+ column density observations, as well as with observations of emission lines from rotationally excited H2 molecules.

The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density$^{0.65}$ can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand $M_\odot$ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $M_{A,0}= 1$) and one with a weak initial field ($M_{A,0}=10$). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium is $<B_{tot}> \sim 42 n_{H,4}^{0.65} \mu$G. Furthermore, the median value of the ratio of the line-of-sight field to density$^{0.65}$ in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong and weak field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70% of the line-of-sight value. This is larger than the 50% usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (~ factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al. (2009), which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas.

We report the first results of DarkSide-50, a direct search for dark matter operating in the underground Laboratori Nazionali del Gran Sasso (LNGS) and searching for the rare nuclear recoils possibly induced by weakly interacting massive particles (WIMPs). The dark matter detector is a Liquid Argon Time Projection Chamber with a (46.4+-0.7) kg active mass, operated inside a 30 t organic liquid scintillator neutron veto, which is in turn installed at the center of a 1 kt water Cherenkov veto for the residual flux of cosmic rays. We report here the null results of a dark matter search for a (1422+-67) kg d exposure with an atmospheric argon fill. This is the most sensitive dark matter search performed with an argon target, corresponding to a 90% CL upper limit on the WIMP-nucleon spin-independent cross section of 6.1x10^-44 cm^2 for a WIMP mass of 100 GeV/c^2.

It is generally inferred from astronomical measurements that Dark Matter (DM) comprises approximately 27\% of the energy-density of the universe. If DM is a subatomic particle, a possible candidate is a Weakly Interacting Massive Particle (WIMP), and the DarkSide-50 (DS) experiment is a direct search for evidence of WIMP-nuclear collisions. DS is located underground at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, and consists of three active, embedded components; an outer water veto (CTF), a liquid scintillator veto (LSV), and a liquid argon (LAr) time projection chamber (TPC). This paper describes the data acquisition and electronic systems of the DS detectors, designed to detect the residual ionization from such collisions.

As part of the DarkSide program of direct dark matter searches using liquid argon TPCs, a prototype detector with an active volume containing 10 kg of liquid argon, DarkSide-10, was built and operated underground in the Gran Sasso National Laboratory in Italy. A critically important parameter for such devices is the scintillation light yield, as photon statistics limits the rejection of electron-recoil backgrounds by pulse shape discrimination. We have measured the light yield of DarkSide-10 using the readily-identifiable full-absorption peaks from gamma ray sources combined with single-photoelectron calibrations using low-occupancy laser pulses. For gamma lines of energies in the range 122-1275 keV, we get consistent light yields averaging 8.887+-0.003(stat)+-0.444(sys) p.e./keVee. With additional purification, the light yield measured at 511 keV increased to 9.142+-0.006(stat) p.e./keVee.

The temperature of the gas in molecular clouds is a key determinant of the characteristic mass of star formation. Ambipolar diffusion (AD) is considered one of the most important heating mechanisms in weakly ionized molecular clouds. In this work, we study the AD heating rate using 2-fluid turbulence simulations and compare it with the overall heating rate due to turbulent dissipation. We find that for observed molecular clouds, which typically have Alfven Mach numbers of ~1 (Crutcher 1999) and AD Reynolds numbers of ~20 (McKee et al. 2010), about 70% of the total turbulent dissipation is in the form of AD heating. AD has an important effect on the length scale where energy is dissipated: when AD heating is strong, most of the energy in the cascade is removed by ion-neutral drift, with a comparatively small amount of energy making it down to small scales. We derive a relation for the AD heating rate that describes the results of our simulations to within a factor of two. Turbulent dissipation, including AD heating, is generally less important that cosmic-ray heating in molecular clouds, although there is substantial scatter in both.

In this paper, we propose a new class of parametrization of the equation of state of dark energy. In contrast with the famous CPL parametrization, these new parametrization of the equation of state does not divergent during the evolution of the Universe even in the future. Also, we perform a observational constraint on two simplest dark energy models belonging to this new class of parametrization, by using the Markov Chain Monte Carlo (MCMC) method and the combined latest observational data from the type Ia supernova compilations including Union2(557), cosmic microwave background, and baryon acoustic oscillation.

Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing positivity of the density and pressure. We have implemented a robust combination of reconstruction schemes, Riemann solvers, limiters, and Constrained Transport EMF averaging schemes that can meet this challenge, and using this strategy, we have developed a new Adaptive Mesh Refinement (AMR) MHD module of the ORION2 code. We investigate the effects of AMR on several statistical properties of a turbulent ideal MHD system with a thermal Mach number of 10 and a plasma $\beta_0$ of 0.1 as initial conditions; our code is shown to be stable for simulations with higher Mach numbers ($M_rms = 17.3$) and smaller plasma beta ($\beta_0 = 0.0067$) as well. Our results show that the quality of the turbulence simulation is generally related to the volume-averaged refinement. Our AMR simulations show that the turbulent dissipation coefficient for supersonic MHD turbulence is about 0.5, in agreement with unigrid simulations.

Ambipolar diffusion (AD) is believed to be a crucial process for redistributing magnetic flux in the dense molecular gas that occurs in regions of star formation. We carry out numerical simulations of this process in regions of low ionization using the heavy ion approximation. The simulations are for regions of strong field (plasma \beta=0.1) and mildly supersonic turbulence (M=3, corresponding to an Alfven mach number of 0.67). The velocity power spectrum of the neutral gas changes from an Iroshnikov-Kraichnan spectrum in the case of ideal MHD to a Burgers spectrum in the case of a shock-dominated hydrodynamic system. The magnetic power spectrum shows a similar behavior. We use a 1D radiative transfer code to post-process our simulation results; the simulated emission from the CS J=2-1 and H13CO+ J=1-0 lines shows that the effects of AD are observable in principle. Linewidths of ions are observed to be less than those of neutrals, and we confirm previous suggestions that this is due to AD. We show that AD is unlikely to affect the Chandrasekhar-Fermi method for inferring field strengths unless the AD is stronger than generally observed. Finally, we present the first fully 3D study of the enhancement of AD by turbulence, finding that AD is accelerated by factor 2-4.5 for non self-gravitating systems with the level of turbulence we consider.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

Ambipolar diffusion is important in redistributing magnetic flux and in damping Alfven waves in molecular clouds. The importance of ambipolar diffusion on a length scale $\ell$ is governed by the ambipolar diffusion Reynolds number, $\rad=\ell/\lad$, where $\lad$ is the characteristic length scale for ambipolar diffusion. The logarithmic mean of the AD Reynolds number in a sample of 15 molecular clumps with measured magnetic fields (Crutcher 1999) is 17, comparable to the theoretically expected value. We identify several regimes of ambipolar diffusion in a turbulent medium, depending on the ratio of the flow time to collision times between ions and neutrals; the clumps observed by Crutcher (1999) are all in the standard regime of ambipolar diffusion, in which the neutrals and ions are coupled over a flow time. We have carried out two-fluid simulations of ambipolar diffusion in isothermal, turbulent boxes for a range of values of $\rad$. The mean Mach numbers were fixed at $\calm=3$ and $\ma=0.67$; self-gravity was not included. We study the properties of overdensities--i.e., clumps--in the simulation and show that the slope of the higher-mass portion of the clump mass spectrum increases as $\rad$ decreases, which is qualitatively consistent with Padoan et al. (2007)'s finding that the mass spectrum in hydrodynamic turbulence is significantly steeper than in ideal MHD turbulence. For a value of $\rad$ similar to the observed value, we find a slope that is consistent with that of the high-mass end of the Initial Mass Function for stars. However, the value we find for the spectral index in our ideal MHD simulation differs from theirs, presumably because our simulations have different initial conditions. This suggests that the mass spectrum of the clumps in the Padoan et al. (2007) turbulent fragmentation model for the IMF depends on the environment, which would conflict with evidence ...

Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse through the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256^3 and 512^3 simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li, McKee, & Klein (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low order statistics of supersonic turbulence, neglecting gravity (e.g. gravitational fragmentation, accretion and merging). In this work the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~1000^3 computational zones) and three different codes (Enzo, Zeus and the Stagger Code). The model predicts a power law distribution for large masses, related to the turbulence energy power spectrum slope, and the shock jump conditions. This power law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 solar masses, as long as the column density of such clouds is not much in excess of 10^22 cm^-2.

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

The currently most viable methods to estimate magnetic field strengths in molecular cloud cores are Zeeman measurements and the Chandrasekhar-Fermi (CF) method. The CF-method estimates magnetic field strengths from polarimetry and relies on equipartition between the turbulent kinetic and turbulent magnetic energy in the observed region. Thus, its application to objects not dominated by hydromagnetic turbulence is questionable. We calibrate the CF-method and Zeeman measurements against numerical models of self-gravitating molecular cloud cores. We find -- in agreement with previous results -- that the CF-method is on the average accurate up to a factor of 2, while a single estimate can be off by a factor of 7. The CF-method is surprisingly robust even when applied to regions not dominated by hydromagnetic turbulence and in the face of small-number statistics. Zeeman measurements can systematically overestimate the field strength in high-density regions, if the field is predominantly oriented along the line of sight, while they tend to underestimate the field in other cases.

We use ZEUS-MP to perform high resolution, three-dimensional, super-Alfvenic turbulent simulations in order to investigate the role of magnetic fields in self-gravitating core formation within turbulent molecular clouds. Statistical properties of our super-Alfvenic model without gravity agree with previous similar studies. Including self-gravity, our models give the following results. They are consistent with the turbulent fragmentation prediction of the core mass distribution of Padoan & Nordlund. They also confirm that local gravitational collapse is not prevented by magnetohydrodynamic waves driven by turbulent flows, even when the turbulent Jeans mass exceeds the mass in the simulation volume. Comparison of results between 256^3 and 512^3 zone simulations reveals convergence in the collapse rate. Analysis of self-gravitating cores formed in the simulation shows that: (1) All cores formed are magnetically supercritical by at least an order of magnitude. (2) A power law relation between central magnetic field strength and density B_c propto rho_c^{1/2} is observed despite the cores being strongly supercritical. (3) Specific angular momentum j propto R^{3/2} for cores with radius R. (4) Most cores are prolate and triaxial in shape, in agreement with the results of Gammie et al.

After a couple of years of quiescence, the soft gamma repeater SGR 1900+14 suddenly reactivated on 18 April 2001, with the emission of a very intense, long and modulated flare, only second in intensity and duration to the 27 August 1998 giant flare. BeppoSAX caught the large flare with its Gamma Ray Burst Monitor and with one of the Wide Field Cameras. The Wide Field Cameras also detected X-ray bursting activity shortly before the giant flare. A target of opportunity observation was started only 8 hours after the large flare with the Narrow Field Instruments, composed of two 60-ks long pointings. These two observations show an X-ray afterglow of the persistent SGR 1900+14 source, decaying with time according to a power law of index -0.6.

The dynamical state of star-forming molecular clouds cannot be understood without determining the structure and strength of their magnetic fields. Measurements of polarized far-infrared radiation from thermally aligned dust grains are used to map the orientation of the field and estimate its strength, but the accuracy of the results has remained in doubt. In order to assess the reliability of this method, we apply it to simulated far-infrared polarization maps derived from three-dimensional simulations of supersonic magnetohydrodynamical turbulence, and compare the estimated values to the known magnetic field strengths in the simulations. We investigate the effects of limited telescope resolution and self-gravity on the structure of the maps. Limited observational resolution affects the field structure such that small scale variations can be completely suppressed, thus giving the impression of a very homogeneous field. The Chandrasekhar-Fermi method of estimating the mean magnetic field in a turbulent medium is tested, and we suggest an extension to measure the rms field. Both methods yield results within a factor of 2 for field strengths typical of molecular clouds, with the modified version returning more reliable estimates for slightly weaker fields. However, neither method alone works well for very weak fields, missing them by a factor of up to 150. Taking the geometric mean of both methods estimates even the weakest fields accurately within a factor of 2.5. Limited telescope resolution leads to a systematic overestimation of the field strengths for all methods. We discuss the effects responsible for this overestimation and show how to extract information on the underlying (turbulent) power spectrum.

We present a 2 - 10 keV ASCA observation of the field around the soft gamma repeater SGR1627-41. A quiescent X-ray source was detected in this observation whose position was consistent both with that of a recently discovered BeppoSAX X-ray source and with the Interplanetary Network localization for this SGR. In 2 - 10 keV X-rays, the spectrum of the X-ray source may be fit equally well by a power law, blackbody, or bremsstrahlung function, with unabsorbed flux ~5 x 10^-12 erg cm^-2 s^-1. We do not confirm a continuation of a fading trend in the flux, and we find no evidence for periodicity, both noted in the earlier BeppoSAXobservations.

We have used the ROSAT High Resolution Imager to search for quiescent X-ray counterparts to four gamma-ray bursts which were localized to small (< 10 sq. arcmin.) error boxes with the Interplanetary Network. The observations took place years after the bursts, and the effective exposure times for each target varied from ~16 - 23 ks. We have not found any X-ray sources inside any of the error boxes. The 0.1 - 2.4 keV 3 sigma flux upper limits range from around 5 x 10^-14 erg cm^-2 s^-1 to 6 x 10^-13 erg cm^-2 s^-1 depending on the burst and the assumed shape of the quiescent spectrum. We consider four types of X-ray emitting galaxies (normal, AGN, faint, and star-forming) and use the flux upper limits to constrain their redshifts. We then use the GRB fluences to constrain the total energies of the bursts.

We present a 2 - 10 keV ASCA observation of the field around the soft gamma repeater SGR1900+14. One quiescent X-ray source was detected in this observation, and it was in the SGR error box. In 2 - 10 keV X-rays, its spectrum may be fit by a power law with index -2.2, and its unabsorbed flux is 9.6 x 10^-12 erg cm^-2 s^-1. We also find a clear 5.16 s period. The properties of the three well-studied soft gamma repeaters are remarkably similar to one another, and provide evidence that all of them are associated with young, strongly magnetized neutron stars in supernova remnants.