Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone.

We investigate the effects of stellar evolution and dust on measurements of stellar velocity dispersion in mergers of disk galaxies. $N$-body simulations and radiative transfer analysis software are used to obtain mass-weighted and flux-weighted measurements of stellar velocity dispersion. We find that the distribution of dust with respect to the distribution of young stars in such systems is more important than the total degree of attenuation. The presence of dust typically causes flux-weighted measurements of stellar velocity dispersion to be elevated with respect to mass-weighted measurements because dust preferentially obscures young stars, which tend to be dynamically cooler than older stellar populations in such systems. In exceptional situations, in which young stars are not preferentially obscured by dust, flux-weighted velocity dispersion measurements tend to be negatively offset with respect to mass-weighted measurements because the dynamically cool young stellar populations are more luminous, per unit mass, than older stellar populations. Our findings provide a context for comparing observationally-obtained measurements of velocity dispersion with measurements of velocity dispersion obtained from galaxy merger simulations.

We introduce NebulOS, a Big Data platform that allows a cluster of Linux machines to be treated as a single computer. With NebulOS, the process of writing a massively parallel program for a datacenter is no more complicated than writing a Python script for a desktop computer. The platform enables most pre-existing data analysis software to be used, as scale, in a datacenter without modification. The shallow learning curve and compatibility with existing software greatly reduces the time required to develop distributed data analysis pipelines. The platform is built upon industry-standard, open-source Big Data technologies, from which it inherits several fault tolerance features. NebulOS enhances these technologies by adding an intuitive user interface, automated task monitoring, and other usability features. We present a summary of the architecture, provide usage examples, and discuss the system's performance scaling.

Highly accurate models of the galaxy population over cosmological volumes are necessary in order to predict the performance of upcoming cosmological missions. We present a data-driven model of the galaxy population constrained by deep 0.1-8 $\rm \mu m$ imaging and spectroscopic data in the COSMOS survey, with the immediate goal of simulating the spectroscopic redshift performance of the proposed SPHEREx mission. SPHEREx will obtain over the full-sky $R\sim41$ spectrophotometry at moderate spatial resolution ($\sim6"$) over the wavelength range 0.75-4.18 $\rm \mu m$ and $R\sim135$ over the wavelength range 4.18-5 $\rm \mu m$. We show that our simulation accurately reproduces a range of known galaxy properties, encapsulating the full complexity of the galaxy population and enables realistic, full end-to-end simulations to predict mission performance. Finally, we discuss potential applications of the simulation framework to future cosmology missions and give a description of released data products.

In order to better understand stellar dynamics in merging systems, such as NGC 6240, we examine the evolution of central stellar velocity dispersion (\sig) in dissipative galaxy mergers using a suite of binary disk merger simulations that include feedback from stellar formation and active galactic nuclei (AGNs). We find that $\sigma_*$ undergoes the same general stages of evolution that were observed in our previous dissipationless simulations: coherent oscillation, then phase mixing, followed by dynamical equilibrium. We also find that measurements of $\sigma_*$ that are based only upon the youngest stars in simulations consistently yield lower values than measurements based upon the total stellar population. This finding appears to be consistent with the so-called "$\sigma_*$ discrepancy," observed in real galaxies. We note that quasar-level AGN activity is much more likely to occur when $\sigma_*$ is near its equilibrium value rather than during periods of extreme $\sigma_*$. Finally, we provide estimates of the scatter inherent in measuring $\sigma_*$ in ongoing mergers.

Using N-body simulations, we studied the detailed evolution of central stellar velocity dispersion, {\sigma}, during dissipationless binary mergers of galaxies. Stellar velocity dispersion was measured using the common mass-weighting method as well as a flux-weighting method designed to simulate the technique used by observers. A toy model for dust attenuation was introduced in order to study the effect of dust attenuation on measurements of {\sigma}. We found that there are three principal stages in the evolution of {\sigma} in such mergers: oscillation, phase mixing, and dynamical equilibrium. During the oscillation stage, {\sigma} undergoes damped oscillations of increasing frequency. The oscillation stage is followed by a phase mixing stage during which the amplitude of the variations in {\sigma} is smaller and more chaotic than in the oscillation stage. Upon reaching dynamical equilibrium, {\sigma} assumes a stable value. We used our data regarding the evolution of {\sigma} during mergers to characterize the scatter inherent in making measurements of {\sigma} in non-quiescent systems. In particular, we found that {\sigma} does not fall below 70% nor exceed 200% of its final, quiescent value during a merger and that a random measurement of {\sigma} in such a system is much more likely to fall near the equilibrium value than near an extremum. Our toy model of dust attenuation suggested that dust can systematically reduce observational measurements of {\sigma} and increase the scatter in {\sigma} measurements.