As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.

Although photometric redshifts (photo-z's) are crucial ingredients for current and upcoming large-scale surveys, the high-quality spectroscopic redshifts currently available to train, validate, and test them are substantially non-representative in both magnitude and color. We investigate the nature and structure of this bias by tracking how objects from a heterogeneous training sample contribute to photo-z predictions as a function of magnitude and color, and illustrate that the underlying redshift distribution at fixed color can evolve strongly as a function of magnitude. We then test the robustness of the galaxy-galaxy lensing signal in 120 deg$^2$ of HSC-SSP DR1 data to spectroscopic completeness and photo-z biases, and find that their impacts are sub-dominant to current statistical uncertainties. Our methodology provides a framework to investigate how spectroscopic incompleteness can impact photo-z-based weak lensing predictions in future surveys such as LSST and WFIRST.

We use the UniverseMachine to analyze the source of scatter between the central galaxy mass, the total stellar mass in the halo, and the dark matter halo mass. We also propose a new halo mass estimator, the cen+N mass: the sum of the stellar mass of the central and the N most massive satellites. We show that, when real space positions are perfectly known, the cen+N mass has scatter competitive with that of richness-based estimators. However, in redshift space, the cen+N mass suffers less from projection effects in the UniverseMachine model. The cen+N mass is therefore a viable low scatter halo mass estimator, and should be considered an important tool to constrain cosmology with upcoming spectroscopic data from DESI. We analyze the scatter in stellar mass at fixed halo mass and show that the total stellar mass in a halo is uncorrelated with secondary halo properties, but that the central stellar mass is a function of both halo mass and halo age. This is because central galaxies in older halos have had more time to grow via accretion. If the UniverseMachine model is correct, accurate galaxy-halo modeling of mass selected samples therefore needs to consider halo age in addition to mass.

Using deep images from the Hyper Suprime-Cam (HSC) survey and taking advantage of its unprecedented weak lensing capabilities, we reveal a remarkably tight connection between the stellar mass distribution of massive central galaxies and their host dark matter halo mass. Massive galaxies with more extended stellar mass distributions tend to live in more massive dark matter haloes. We explain this connection with a phenomenological model that assumes, (1) a tight relation between the halo mass and the total stellar content in the halo, (2) that the fraction of in-situ and ex-situ mass at $r<10$ kpc depends on halo mass. This model provides an excellent description of the stellar mass functions (SMF) of total stellar mass ($M_{\star}^{\rm Max}$) and stellar mass within inner 10 kpc ($M_{\star}^{10}$) and also reproduces the HSC weak lensing signals of massive galaxies with different stellar mass distributions. The best-fit model shows that halo mass varies significantly at fixed total stellar mass (as much as 0.4 dex) with a clear dependence on $M_{\star}^{10}$. Our two-parameter $M_{\star}^{\rm Max}$-$M_{\star}^{10}$ description provides a more accurate picture of the galaxy-halo connection at the high-mass end than the simple stellar-halo mass relation (SHMR) and opens a new window to connect the assembly history of halos with those of central galaxies. The model also predicts that the ex-situ component dominates the mass profiles of galaxies at $r< 10$ kpc for $\log M_{\star} \ge 11.7$). The code used for this paper is available online: https://github.com/dr-guangtou/asap

Protostellar outflows provide a means to probe the accretion process of forming stars and their ability to inject energy into their surroundings. However, conclusions based on outflow observations depend upon the degree of accuracy with which their properties can be estimated. We examine the quality of Atacama Large Millimeter/submillimeter Array (ALMA) observations of protostellar outflows by producing synthetic $^{12}$CO(1-0) and $^{13}$CO(1-0) observations of numerical simulations. We use various ALMA configurations, observational parameters, and outflow inclinations to assess how accurately different assumptions and setups can recover underlying properties. We find that more compact arrays and longer observing times can improve the mass and momentum recovery by a factor of two. During the first $\sim$0.3 Myr of evolution, $^{12}$CO(1-0) is optically thick, even for velocities $|v|\ge 1$ km s$^{-1}$, and outflow mass is severely underestimated without an optical depth correction. Likewise, $^{13}$CO(1-0) is optically thick during the first $\simeq 0.1$ Myr. However, underestimation due to shorter observing time, missing flux, and optical depth are partially offset by the assumption of local thermodynamic equilibrium and higher excitation temperatures. Overall, we expect that full ALMA $^{13}$CO(1-0) observations of protostellar sources within 500 pc with observing times $\gtrsim 1$ hrs and assumed excitation temperatures of $T<20$K will reliably measure mass and line-of-sight momentum to within 20\%.

Correlations between 1-10 Hz quasi-periodic oscillations (QPOs) and spectral power law index have been reported for black hole (BH) candidate sources and one neutron star source, 4U 1728-34. An examination of QPO frequency and index relationships in Sco X-1 is reported herein. We discovered that Sco X-1, representing Z-source groups, can be adequately modeled by a simple two-component model of Compton up-scattering with a soft photon electron temperature of about 0.4 keV, plus an Iron K-line. The results show a strong correlation between spectral power law index and kHz QPOs. Because Sco X-1 radiates near the Eddington limit, one can infer that the geometrical configuration of the Compton cloud (CC) is quasi-spherical because of high radiation pressure in the CC. Thus, we conclude that the high Thomson optical depth of the Compton cloud, in the range of ~5-6 from the best-fit model parameters, is consistent with the neutron star's surface being obscured by material. Moreover, a spin frequency of Sco X-1 is likely suppressed due to photon scattering off CC electrons. Additionally, we demonstrate how the power spectrum evolves when Sco X-1 transitions from the horizontal branch to the normal branch.

We present a method for determining the B-field around neutron stars based on observed kHz and viscous QPO frequencies used in combination with the best-fit optical depth and temperature of a Comptonization model. In the framework of the transition layer QPO model, we analyze magnetoacoustic wave formation in the layer between a neutron star surface and the inner edge of a Keplerian disk. We derive formulas for the magnetoacoustic wave frequencies for different regimes of radial transition layer oscillations. We demonstrate that our model can use the QPO as a new kind of probe to determine the magnetic field strengths for 4U 1728-42, GX 340+0, and Sco X-1 in the zone where the QPOs occur. Observations indicate that the dependence of the viscous frequency on the Keplerian frequency is closely related to the inferred dependence of the magnetoacoustic wave frequency on the Keplerian frequency for a dipole magnetic field. The magnetoacoustic wave dependence is based on a single parameter, the magnetic moment of the star as estimated from the field strength in the transition layer. The best-fit magnetic moment parameter is about (0.5-1)x 10^{25} G cm^3 for all studied sources. From observational data, the magnetic fields within distances less 20 km from neutron star for all three sources are strongly constrained to be dipole fields with the strengths 10^{7-8} G on the neutron star surface.

We present a comprehensive classification of all observed QPOs within the framework of the transition layer model using a large set of Rossi X-ray Timing Explorer (RXTE) data for Sco X-1. The model assumes an optically thin material along the observer's line of sight in the horizontal branch and an increasingly optically thick material while in the other two branches that is consistent with X-ray and radio observations and the disk transition layer model of QPOs. We identify the ~ 6 Hz frequencies in the normal branch as acoustic oscillations of a spherical shell around the neutron star (NS) that is formed after radiation pressure near the Eddington accretion rate destroys the disk. The size of the shell is on the order of one NS radii from the NS. We also estimate the upper limit of Sco X-1's magnetic field to be 0.7 x 10^6 G at about one NS radii above the NS surface while in the horizontal X-ray branch.

The radio components associated with the LMXB Sco X-1 have been monitored with extensive VLBI imaging at 1.7 and 5.0 GHz over four years, including a 56-hour continuous VLBI observation in 1999 June. We often detected one strong and one weak compact radio component, moving in opposite directions from the radio core. We suggest that the moving components are radio lobes generated by the disruption of energy flow in a twin-beam from the binary system. The average lifetime of a lobe-pair, the space motion of the lobes and the measured energy flow in the beam are discussed in arXiv:astro-ph/0104325. The lobe has a flux density that is variable over a time-scale of one hour, a measured minimum size of 1 mas (2.8 au), and is extended perpendicular to its motion. This short electron radiative lifetime may be caused by synchrotron losses if the lobe magnetic field is 300 G, or by adiabatic expansion of the electrons as soon as they are produced at the working surface. The lobes also show periods of slow expansion and a steepening radio spectrum, perhaps related to the characteristics of the beam energy flow. The radio morphology for Sco X-1 is more simple than for most other Galactic jet sources. The lobes of Sco X-1 are similar to hot-spots found in many extragalactic double sources. Scaling the phenomena observed in Sco X-1 to extragalactic sources implies hot-spot variability time-scales of 10^4 yr and hot-spot lifetimes of 10^5 yr. The recurrent formation of lobes in Sco X-1 probably does not occur for extragalactic radio sources.

The evolution of the radio emission from Sco X-1 is determined from a 56-hour continuous VLBI observation and from shorter observations over a four-year period. The radio source consists of a variable core near the binary, and two variable compact radio lobes which form near the core, move diametrically outward, then fade away. Subsequently, a new lobe-pair form near the core and the behavior repeats. The differences in the radio properties of the two lobes are consistent with the delay and Doppler-boosting associated with an average space velocity of 0.45c at 44 deg to the line of sight. Four lobe speeds, between 0.32c and 0.57c, were measured for several lobe-pairs on different days. The speed during each epoch remained constant over many hours. The direction of motion of the lobes over all epochs remained constant to a few degrees. Two core flares are contemporaneous with two lobe flares after removal of the delay associated with an energy burst moving with speed >0.95c in a twin-beam from the core to each lobe. This is the first direct measurement of the speed of energy flow within an astrophysical jet. The similarity of the core and lobe flares suggests that the twin-beam flow is symmetric and that the core is located near the base of the beam.