The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA.

Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $\sigma_{\chi A} \propto A^2 \mu_A^2 \sigma_{\chi N}\propto A^4 \sigma_{\chi N}$, where the $A^2$ comes from coherence and the $\mu_A^2\simeq A^2 m_N^2$ from kinematics for $m_\chi\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $\sigma_{\chi A} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $\sigma_{\chi N}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $\sigma_{\chi N}>10^{-32}\;\text{cm}^2$ with $m_\chi\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted.

The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to super-sample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0.