Showing 1–28 of 28 results for author: Wibking, B

On the limitations of H alpha luminosity as a star formation tracer in spatially resolved observations

Authors: Zipeng Hu, Benjamin D. Wibking, Mark R. Krumholz, Christoph Federrath

Abstract: This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scal… ▽ More This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scales from whole-galaxy to individual molecular clouds ($\lesssim 100$ pc). Our results reveal significant spatial mismatches between H$α$ and true YSO maps on sub-100 pc scales, primarily due to ionizing photon leakage, with a secondary contribution from young stars drifting away from their parent molecular clouds. On small scales these effect contribute significantly to the observed anti-correlation between gas and star formation, such that there is noticeably less anti-correlation if we replace an H$α$-based star formation map with a YSO-based one; this in turn implies that previous studies have underestimated the time it takes for young stars to disperse their parent molecular clouds. However, these effects are limited in dense regions with hydrogen columns $N_\mathrm{H} > 3 \times 10^{21}$ cm$^{-2}$, where the H$α$- and YSO-based SFR maps show better agreement. Based on this finding we propose a calibration model that can precisely measure the SFR of large molecular clouds (mean radius > 100 pc) with a combination of H$α$ and CO observations, which provides a foundation for future study of star formation processes in extragalactic molecular clouds. △ Less

Submitted 16 June, 2024; originally announced June 2024.

OGLE-2015-BLG-0845L: A low-mass M dwarf from the microlensing parallax and xallarap effects

Authors: Zhecheng Hu, Wei Zhu, Andrew Gould, Andrzej Udalski, Takahiro Sumi, Ping Chen, Sebastiano Calchi Novati, Jennifer C. Yee, Charles A. Beichman, Geoffery Bryden, Sean Carey, Michael Fausnaugh, B. Scott Gaudi, Calen B. Henderson, Yossi Shvartzvald, Benjamin Wibking, Przemek Mróz, Jan Skowron, Radoslaw Poleski, Michaeł K. Szymański, Igor Soszynśki, Paweł Pietrukowicz, Szymon Kozłowski, Krzysztof Ulaczyk, Krzysztof A. Rybicki , et al. (29 additional authors not shown)

Abstract: We present the analysis of the microlensing event OGLE-2015-BLG-0845, which was affected by both the microlensing parallax and xallarap effects. The former was detected via the simultaneous observations from the ground and Spitzer, and the latter was caused by the orbital motion of the source star in a relatively close binary. The combination of these two effects led to a direct mass measurement o… ▽ More We present the analysis of the microlensing event OGLE-2015-BLG-0845, which was affected by both the microlensing parallax and xallarap effects. The former was detected via the simultaneous observations from the ground and Spitzer, and the latter was caused by the orbital motion of the source star in a relatively close binary. The combination of these two effects led to a direct mass measurement of the lens object, revealing a low-mass ($0.14 \pm 0.05 M_{\odot}$) M-dwarf at the bulge distance ($7.6 \pm 1.0$ kpc). The source binary consists of a late F-type subgiant and a K-type dwarf of $\sim1.2 M_{\odot}$ and $\sim 0.9 M_{\odot}$, respectively, and the orbital period is $70 \pm 10$ days. OGLE-2015-BLG-0845 is the first single-lens event in which the lens mass is measured via the binarity of the source. Given the abundance of binary systems as potential microlensing sources, the xallarap effect may not be a rare phenomenon. Our work thus highlights the application of the xallarap effect in the mass determination of microlenses, and the same method can be used to identify isolated dark lenses. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 13 pages, 10 figures. Submitted to Monthly Notices of the Royal Astronomical Society

An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics

Authors: Chong-Chong He, Benjamin D. Wibking, Mark R. Krumholz

Abstract: Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable prop… ▽ More Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by orders of magnitude across the computational domain. Here we present a new implicit-explicit (IMEX) method for numerical RHD that has a number of desirable properties that have not previously been combined in a single method. Our scheme is based on moments and allows machine-precision conservation of energy and momentum, making it highly suitable for adaptive mesh refinement applications; it requires no more communication than hydrodynamics and includes no non-local iterative steps, making it highly suitable for massively parallel and GPU-based systems where communication is a bottleneck; and we show that it is asymptotically-accurate in the streaming, static diffusion, and dynamic diffusion limits, including in the so-called asymptotic diffusion regime where the computational grid does not resolve the photon mean free path. We implement our method in the GPU-accelerated RHD code QUOKKA and show that it passes a wide range of numerical tests. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 14 pages, 6 figures, submitted to MNRAS; Code available for download at https://github.com/quokka-astro/quokka

Quokka-based Understanding of Outflows (QED) I. Metal loading, phase structure, and convergence testing for Solar neighbourhood conditions

Authors: Aditi Vijayan, Mark R. Krumholz, Benjamin D. Wibking

Abstract: Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to res… ▽ More Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to resolve mass, momentum and energy exchanges between the different phases of the outflows. In this context, we present simulations of outflows, driven by SN feedback, conducted using \textsc{Quokka}, a new GPU-optimised AMR radiation-hydrodynamics code. This code allows us to reach combinations of resolution, simulation volume, and simulation duration larger than those that have previously been possible, and to resolve all gas phases from cold neutral medium, $T \sim 100$ K, to hot ionised gas, $T \gtrsim 10^7$ K. In this, a first of a series of papers exploring generation and evolution of multiphase outflows from a wide range of galactic environments and star formation rates, we quantify the extent of selective metal loading in Solar neighbourhood-like environments. We explain the selective metal loading we find as a result of the transport of metals within and between phases, a phenomenon we can study owing to the parsec-scale resolution that our simulations achieve. We also quantify the sensitivity of metal loading studies to numerical resolution, and present convergence criteria for future studies. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS. Full title="Quokka-based Understanding of Outflows Derived from Extensive, Repeated, Accurate, Thorough, Demanding, Expensive, Memory-consuming, Ongoing Numerical Simulations of Transport, Removal, Accretion, Nucleosynthesis, Deposition, and Uplifting of Metals (QUOD ERAT DEMONSTRANDUM, i.e, QED)"

The sub-critical illusion: synthetic Zeeman effect observations from galactic zoom-in simulations

Authors: Zipeng Hu, Benjamin D. Wibking, Mark R. Krumholz

Abstract: Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this pap… ▽ More Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this paper, we carry out simulations of the formation of magnetised molecular clouds, zooming in from an entire galaxy to sub-pc scales, which we post-process to produce synthetic HI and OH Zeeman measurements. The mass-to-flux ratios we recover from the simulated observations show a transition in magnetic criticality that closely matches observations, but we find that the gravitational-magnetic energy ratios on corresponding scales are mostly super-critical, even in the HI regime. We conclude that HI clouds in the process of assembling to form molecular clouds are already super-critical even before H_2 forms, and that the apparent transition from sub- to super-criticality between HI and H_2 is primarily an illusion created by chemical and excitation biases affecting the Zeeman measurements. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 12 pages, 10 figures, submitted to MNRAS

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov… ▽ More We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used -- but less accurate -- closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 7 pages, 2 figures, to appear in the proceedings of IAUS362: The predictive power of computational astrophysics as discovery tool (D. Bisikalo, T. Hanawa, C. Boily, J. Stone, eds.). arXiv admin note: substantial text overlap with arXiv:2202.08778

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing.… ▽ More We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS

Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs

Authors: Benjamin D. Wibking, Mark R. Krumholz

Abstract: We present Quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics code optimised for graphics processing units (GPUs). Quokka solves the equations of hydrodynamics with the piecewise parabolic method (PPM) in a method-of-lines formulation, and handles radiative transfer via the variable Eddington tensor (VET) radiation moment equations with a loca… ▽ More We present Quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics code optimised for graphics processing units (GPUs). Quokka solves the equations of hydrodynamics with the piecewise parabolic method (PPM) in a method-of-lines formulation, and handles radiative transfer via the variable Eddington tensor (VET) radiation moment equations with a local closure. We use the AMReX library to handle the adaptive mesh management. In order to maximise GPU performance, we combine explicit-in-time evolution of the radiation moment equations with the reduced speed-of-light approximation. We show results for a wide range of test problems for hydrodynamics, radiation, and coupled radiation hydrodynamics. On uniform grids in 3D on a single GPU, our code achieves > 250 million hydrodynamic updates per second and almost 40 million radiation hydrodynamic updates per second. For radiation hydrodynamics problems on uniform grids in 3D, our code scales from 4 GPUs to 256 GPUs with an efficiency of 76 per cent. The code is publicly released under an open-source license on GitHub. △ Less

Submitted 30 January, 2022; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: Revised version re-submitted to MNRAS, added new performance measurements with A100 GPUs

Exploiting Non-linear Scales in Galaxy-Galaxy Lensing and Galaxy Clustering: A Forecast for the Dark Energy Survey

Authors: Andrés N. Salcedo, David H. Weinberg, Hao-Yi Wu, Benjamin D. Wibking

Abstract: The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a powerful probe of low redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we extend the N-body and halo occupation distribution (HOD) emulator method of arxiv:1907.06293 to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding paramet… ▽ More The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a powerful probe of low redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we extend the N-body and halo occupation distribution (HOD) emulator method of arxiv:1907.06293 to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias $A_{lens}$. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile $ΔΣ(r_p)$ and the projected galaxy correlation function $w_{p,gg}(r_p)$ in the final (Year 6) DES data set over scales $r_p=0.3-30h^{-1}$ Mpc. For a $3\%$ prior on $A_{lens}$ we forecast precisions of $1.9\%$, $2.0\%$, and $1.9\%$ on $Ω_m$, $σ_8$, and $S_8 \equiv σ_8Ω_m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as $A_{lens}$ and a point-mass contribution to $ΔΣ$. Adding scales $r_p=0.3-3h^{-1}$ Mpc improves the $S_8$ precision by a factor of $\sim1.6$ relative to a large scale ($3.0-30.0h^{-1}$ Mpc) analysis, equivalent to increasing the survey area by a factor of ${\sim}2.6$. Sharpening the $A_{lens}$ prior to $1\%$ further improves the $S_8$ precision by a factor of $1.7$ (to $1.1\%$), and it amplifies the gain from including non-linear scales. Our emulator achieves percent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: 17 pages, 7 figures, to be submitted to MNRAS

The global structure of magnetic fields and gas in simulated Milky Way-analogue galaxies

Authors: Benjamin D. Wibking, Mark R. Krumholz

Abstract: We simulate an isolated, magnetised Milky Way-like disc galaxy using a self-consistent model of unresolved star formation and feedback, evolving the system until it reaches statistical steady state. We show that the quasi-steady-state structure is distinctly layered in galactocentric height $z$, with an innermost region having comparable gas and magnetic pressures (plasma beta $β\sim 1$), an outer… ▽ More We simulate an isolated, magnetised Milky Way-like disc galaxy using a self-consistent model of unresolved star formation and feedback, evolving the system until it reaches statistical steady state. We show that the quasi-steady-state structure is distinctly layered in galactocentric height $z$, with an innermost region having comparable gas and magnetic pressures (plasma beta $β\sim 1$), an outermost region having dominant gas pressures ($β\gg 1$), and an intermediate region between $300$ pc $\lesssim |z| \lesssim 3$ kpc that is dynamically dominated by magnetic fields ($β\ll 1$). We find field strengths, gas surface densities, and star formation rates that agree well with those observed both in the Galactic centre and in the Solar neighbourhood. The most significant dynamical effect of magnetic fields on the global properties of the disc is a reduction of the star formation rate by a factor of 1.5-2 with respect to an unmagnetised control simulation. At fixed star formation rate, there is no significant difference in the mass outflow rates or profiles between the magnetised and non-magnetised simulations. Our results for the global structure of the magnetic field have significant implications for models of cosmic ray-driven winds and cosmic-ray propagation in the Galaxy, and can be tested against observations with the forthcoming Square Kilometre Array and other facilities. Finally, we report the discovery of a physical error in the implementation of neutral gas heating and cooling in the popular GIZMO code, which may lead to qualitatively incorrect phase structures if not corrected. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 12 pages + appendix, 7 figures. Simulated with 100 percent renewable electricity

Cosmology with galaxy cluster weak lensing: statistical limits and experimental design

Authors: Hao-Yi Wu, David H. Weinberg, Andrés N. Salcedo, Benjamin D. Wibking

Abstract: We forecast constraints on the amplitude of matter clustering sigma_8(z) achievable with the combination of cluster weak lensing and number counts, in current and next-generation weak lensing surveys. We advocate an approach, analogous to galaxy-galaxy lensing, in which the observables in each redshift bin are the mean number counts and the mean weak lensing profile of clusters above a mass proxy… ▽ More We forecast constraints on the amplitude of matter clustering sigma_8(z) achievable with the combination of cluster weak lensing and number counts, in current and next-generation weak lensing surveys. We advocate an approach, analogous to galaxy-galaxy lensing, in which the observables in each redshift bin are the mean number counts and the mean weak lensing profile of clusters above a mass proxy threshold. The primary astrophysical nuisance parameter is the logarithmic scatter between the mass proxy and true mass near the threshold. For surveys similar to the Dark Energy Survey (DES), the Roman Space Telescope High Latitude Survey (HLS), and the Rubin Observatory Legacy Survey of Space and Time (LSST), we forecast aggregate precision on sigma_8 of 0.26%, 0.24%, and 0.10%, respectively, if the mass-observable scatter has an external prior better than 0.01. These constraints would be degraded by about 20% for a 0.05 prior on scatter in the case of DES or HLS and for a 0.016 prior for LSST. A one-month observing program with Roman Space Telescope targeting approximately 2500 massive clusters could achieve a 0.5% constraint on sigma_8(z=0.7) on its own, or a ~0.33% constraint in combination with the HLS. Realizing the constraining power of clusters requires accurate knowledge of the mass-observable relation and stringent control of systematics. We provide analytic approximations to our numerical results that allow easy scaling to other survey assumptions or other methods of cluster mass estimation. △ Less

Submitted 26 March, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: 23 pages, 16 figures; replaced to match published version

Covariance matrices for galaxy cluster weak lensing: from virial regime to uncorrelated large-scale structure

Authors: Hao-Yi Wu, David H. Weinberg, Andrés N. Salcedo, Benjamin D. Wibking, Ying Zu

Abstract: Next-generation optical imaging surveys will revolutionise the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological paramet… ▽ More Next-generation optical imaging surveys will revolutionise the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological parameters. We use a combination of analytical calculations and high-resolution N-body simulations to derive accurate covariance matrices that span from the virial regime to linear scales of the cluster-matter cross-correlation. We validate this calculation using a public ray-tracing lensing simulation and provide a software package for calculating covariance matrices for a wide range of cluster and source sample choices. We discuss the relative importance of shape noise and density fluctuations, the impact of radial bin size, and the impact of off-diagonal elements. For a weak lensing source density 10 per square arcmin, shape noise typically dominates the variance on comoving scales less than 5 Mpc/h. However, for 60 per square arcmin, potentially achievable with future weak lensing experiments, density fluctuations typically dominate the variance at scales greater than 1 Mpc/h and remain comparable to shape noise on smaller scales. △ Less

Submitted 15 October, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 23 pages, 14 figures; replaced to match published version; software at https://github.com/hywu/cluster-lensing-cov

Cosmology with galaxy-galaxy lensing on non-perturbative scales: Emulation method and application to BOSS LOWZ

Authors: Benjamin D. Wibking, David H. Weinberg, Andrés N. Salcedo, Hao-Yi Wu, Sukhdeep Singh, Sergio Rodríguez-Torres, Lehman H. Garrison, Daniel J. Eisenstein

Abstract: We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \… ▽ More We describe our nonlinear emulation (i.e., interpolation) framework that combines the halo occupation distribution (HOD) galaxy bias model with $N$-body simulations of nonlinear structure formation, designed to accurately predict the projected clustering and galaxy-galaxy lensing signals from luminous red galaxies (LRGs) in the redshift range $0.16 < z < 0.36$ on comoving scales $0.6 < r_p < 30$ \hMpc. The interpolation accuracy is $\lesssim 1-2$ per cent across the entire physically plausible range of parameters for all scales considered. We correctly recover the true value of the cosmological parameter $S_8 = ({σ_8}/{0.8228}) ({Ω_{\text{m}}}/{0.3107})^{0.6}$ from mock measurements produced via subhalo abundance matching (SHAM)-based lightcones designed to approximately match the properties of the SDSS LOWZ galaxy sample. Applying our model to Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 14 (DR14) LOWZ galaxy clustering and galaxy-shear cross-correlation measurements made with Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) imaging, we perform a prototype cosmological analysis marginalizing over $w$CDM cosmological parameters and galaxy HOD parameters. We obtain a 4.4 per cent measurement of $S_8 = 0.847 \pm 0.037$, in $3.5σ$ tension with the Planck cosmological results of $1.00 \pm 0.02$. We discuss the possibility of underestimated systematic uncertainties or astrophysical effects that could explain this discrepancy. △ Less

Submitted 3 December, 2019; v1 submitted 14 July, 2019; originally announced July 2019.

Comments: Updated to match accepted version

Cosmology with Stacked Cluster Weak Lensing and Cluster-Galaxy Cross-Correlations

Authors: Andrés N. Salcedo, Benjamin D. Wibking, David H. Weinberg, Hao-Yi Wu, Douglas Ferrer, Daniel Eisenstein, Philip Pinto

Abstract: Cluster weak lensing is a sensitive probe of cosmology, particularly the amplitude of matter clustering $σ_8$ and matter density parameter $Ω_m$. The main nuisance parameter in a cluster weak lensing cosmological analysis is the scatter between the true halo mass and the relevant cluster observable, denoted $σ_{\ln Mc}$. We show that combining the cluster weak lensing observable $ΔΣ$ with the proj… ▽ More Cluster weak lensing is a sensitive probe of cosmology, particularly the amplitude of matter clustering $σ_8$ and matter density parameter $Ω_m$. The main nuisance parameter in a cluster weak lensing cosmological analysis is the scatter between the true halo mass and the relevant cluster observable, denoted $σ_{\ln Mc}$. We show that combining the cluster weak lensing observable $ΔΣ$ with the projected cluster-galaxy cross-correlation function $w_{p,cg}$ and galaxy auto-correlation function $w_{p,gg}$ can break the degeneracy between $σ_8$ and $σ_{\ln Mc}$ to achieve tight, percent-level constraints on $σ_8$. Using a grid of cosmological N-body simulations, we compute derivatives of $ΔΣ$, $w_{p,cg}$, and $w_{p,gg}$ with respect to $σ_8$, $Ω_m$, $σ_{\ln Mc}$ and halo occupation distribution (HOD) parameters describing the galaxy population. We also compute covariance matrices motivated by the properties of the Dark Energy Suvery (DES) cluster and weak lensing survey and the BOSS CMASS galaxy redshift survey. For our fiducial scenario combining $ΔΣ$, $w_{p,cg}$, and $w_{p,gg}$ measured over $0.3-30.0 \; h^{-1} \; \mathrm{Mpc}$, for clusters at $z=0.35-0.55$ above a mass threshold $M_c\approx 2\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$, we forecast a $1.4\%$ constraint on $σ_8$ while marginalizing over $σ_{\ln Mc}$ and all HOD parameters. Reducing the mass threshold to $1\times 10^{14} \; h^{-1} \; \mathrm{M_{\odot}}$ and adding a $z=0.15-0.35$ redshift bin sharpens this constraint to $0.8\%$. The small scale $(r_p < 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``mass function'' and large scale $(r_p > 3.0 \; h^{-1} \; \mathrm{Mpc})$ ``halo-mass cross-correlation'' regimes of $ΔΣ$ have comparable constraining power, allowing internal consistency tests from such an analysis. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 23 pages, 11 figures, 1 appendix (with extra figure)

Radiation pressure in galactic disks: stability, turbulence, and winds in the single-scattering limit

Authors: Benjamin D. Wibking, Todd A. Thompson, Mark R. Krumholz

Abstract: The radiation force on dust grains may be dynamically important in driving turbulence and outflows in rapidly star-forming galaxies. Recent studies focus on the highly optically-thick limit relevant to the densest ultra-luminous galaxies and super star clusters, where reprocessed infrared photons provide the dominant source of electromagnetic momentum. However, even among starburst galaxies, the g… ▽ More The radiation force on dust grains may be dynamically important in driving turbulence and outflows in rapidly star-forming galaxies. Recent studies focus on the highly optically-thick limit relevant to the densest ultra-luminous galaxies and super star clusters, where reprocessed infrared photons provide the dominant source of electromagnetic momentum. However, even among starburst galaxies, the great majority instead lie in the so-called "single-scattering" limit, where the system is optically-thick to the incident starlight, but optically-thin to the re-radiated infrared. In this paper we present a stability analysis and multidimensional radiation-hydrodynamic simulations exploring the stability and dynamics of isothermal dusty gas columns in this regime. We describe our algorithm for full angle-dependent radiation transport based on the discontinuous Galerkin finite element method. For a range of near-Eddington fluxes, we show that the medium is unstable, producing convective-like motions in a turbulent atmosphere with a scale height significantly inflated compared to the gas pressure scale height and mass-weighted turbulent energy densities of $\sim 0.01-0.1$ of the midplane radiation energy density, corresponding to mass-weighted velocity dispersions of Mach number $\sim 0.5-2$. Extrapolation of our results to optical depths of $10^3$ implies maximum turbulent Mach numbers of $\sim20$. Comparing our results to galaxy-averaged observations, and subject to the approximations of our calculations, we find that radiation pressure does not contribute significantly to the effective supersonic pressure support in star-forming disks, which in general are substantially sub-Eddington. We further examine the time-averaged vertical density profiles in dynamical equilibrium and comment on implications for radiation-pressure-driven galactic winds. △ Less

Submitted 4 January, 2018; originally announced January 2018.

Comments: 19 pages, 15 figures, submitted to MNRAS. Comments welcome

Emulating galaxy clustering and galaxy-galaxy lensing into the deeply nonlinear regime: methodology, information, and forecasts

Authors: Benjamin D. Wibking, Andrés N. Salcedo, David H. Weinberg, Lehman H. Garrison, Douglas Ferrer, Jeremy Tinker, Daniel Eisenstein, Marc Metchnik, Philip Pinto

Abstract: The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy… ▽ More The combination of galaxy-galaxy lensing (GGL) with galaxy clustering is one of the most promising routes to determining the amplitude of matter clustering at low redshifts. We show that extending clustering+GGL analyses from the linear regime down to $\sim 0.5 \, h^{-1}$ Mpc scales increases their constraining power considerably, even after marginalizing over a flexible model of non-linear galaxy bias. Using a grid of cosmological N-body simulations, we construct a Taylor-expansion emulator that predicts the galaxy autocorrelation $ξ_{\text{gg}}(r)$ and galaxy-matter cross-correlation $ξ_{\text{gm}}(r)$ as a function of $σ_8$, $Ω_m$, and halo occupation distribution (HOD) parameters, which are allowed to vary with large scale environment to represent possible effects of galaxy assembly bias. We present forecasts for a fiducial case that corresponds to BOSS LOWZ galaxy clustering and SDSS-depth weak lensing (effective source density $\sim 0.3$ arcmin$^{-2}$). Using tangential shear and projected correlation function measurements over $0.5 \leq r_p \leq 30 \, h^{-1}$ Mpc yields a 1.8% constraint on the parameter combination $σ_8Ω_m^{0.58}$, a factor of two better than a constraint that excludes non-linear scales ($r_p > 2 \, h^{-1}$ Mpc, $4 \, h^{-1}$ Mpc for $γ_t,w_p$). Much of this improvement comes from the non-linear clustering information, which breaks degeneracies among HOD parameters that would otherwise degrade the inference of matter clustering from GGL. Increasing the effective source density to $3$ arcmin$^{-2}$ sharpens the constraint on $σ_8Ω_m^{0.58}$ by a further factor of two. With robust modeling into the non-linear regime, low-redshift measurements of matter clustering at the 1-percent level with clustering+GGL alone are well within reach of current data sets such as those provided by the Dark Energy Survey. △ Less

Submitted 20 September, 2017; originally announced September 2017.

Comments: 23 pages, 7 figures

Towards Accurate Modelling of Galaxy Clustering on Small Scales: Testing the Standard $Λ\mathrm{CDM}$ + Halo Model

Authors: Manodeep Sinha, Andreas A. Berlind, Cameron K. McBride, Roman Scoccimarro, Jennifer A. Piscionere, Benjamin D. Wibking

Abstract: Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In th… ▽ More Interpreting the small-scale clustering of galaxies with halo models can elucidate the connection between galaxies and dark matter halos. Unfortunately, the modelling is typically not sufficiently accurate for ruling out models statistically. It is thus difficult to use the information encoded in small scales to test cosmological models or probe subtle features of the galaxy-halo connection. In this paper, we attempt to push halo modelling into the "accurate" regime with a fully numerical mock-based methodology and careful treatment of statistical and systematic errors. With our forward-modelling approach, we can incorporate clustering statistics beyond the traditional two-point statistics. We use this modelling methodology to test the standard $Λ\mathrm{CDM}$ + halo model against the clustering of SDSS DR7 galaxies. Specifically, we use the projected correlation function, group multiplicity function and galaxy number density as constraints. We find that while the model fits each statistic separately, it struggles to fit them simultaneously. Adding group statistics leads to a more stringent test of the model and significantly tighter constraints on model parameters. We explore the impact of varying the adopted halo definition and cosmological model and find that changing the cosmology makes a significant difference. The most successful model we tried (Planck cosmology with Mvir halos) matches the clustering of low luminosity galaxies, but exhibits a 2.3$σ$ tension with the clustering of luminous galaxies, thus providing evidence that the "standard" halo model needs to be extended. This work opens the door to adding interesting freedom to the halo model and including additional clustering statistics as constraints. △ Less

Submitted 6 August, 2018; v1 submitted 16 August, 2017; originally announced August 2017.

Comments: Replaced to match the published version

Journal ref: http://adsabs.harvard.edu/abs/2018MNRAS.478.1042S

Ground-based parallax confirmed by Spitzer: binary microlensing event MOA-2015-BLG-020

Authors: Tianshu Wang, Wei Zhu, Shude Mao, I. A. Bond, A. Gould, A. Udalski, T. Sumi, V. Bozza, C. Ranc, A. Cassan, J. C. Yee, C. Han, F. Abe, Y. Asakura, R. Barry, D. P. Bennett, A. Bhattacharya, M. Donachie, P. Evans, A. Fukui, Y. Hirao, Y. Itow, K. Kawasaki, N. Koshimoto, M. C. A. Li , et al. (61 additional authors not shown)

Abstract: We present the analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameter… ▽ More We present the analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameters constrained by the ground-based data are confirmed by the Spitzer observations through the satellite parallax method. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined. It is found that the binary lens is composed of two dwarf stars with masses $M_1 = 0.606 \pm 0.028M_\odot$ and $M_2 = 0.125 \pm 0.006M_\odot$ in the Galactic disk. Assuming the source star is at the same distance as the bulge red clump stars, we find the lens is at a distance $D_L = 2.44 \pm 0.10 kpc$. In the end, we provide a summary and short discussion of all published microlensing events in which the annual parallax effect is confirmed by other independent observations. △ Less

Submitted 26 July, 2017; v1 submitted 24 July, 2017; originally announced July 2017.

Comments: 16 pages, 5 figures

OGLE-2015-BLG-1482L: the first isolated low-mass microlens in the Galactic bulge

Authors: S. -J. Chung, W. Zhu, A. Udalski, C. -U. Lee, Y. -H. Ryu, Y. K. Jung, I. -G. Shin, J. C. Yee, K. -H. Hwang, A. Gould, M. Albrow, S. -M. Cha, C. Han, D. -J. Kim, H. -W. Kim, S. -L. Kim, Y. -H. Kim, Y. Lee, B. -G. Park, R. W. Pogge, R. Poleski, P. Mróz, P. Pietrukowicz, J. Skowron, M. K. Szymański , et al. (13 additional authors not shown)

Abstract: We analyze the single microlensing event OGLE-2015-BLG-1482 simultaneously observed from two ground-based surveys and from \textit{Spitzer}. The \textit{Spitzer} data exhibit finite-source effects due to the passage of the lens close to or directly over the surface of the source star as seen from \textit{Spitzer}. Such finite-source effects generally yield measurements of the angular Einstein radi… ▽ More We analyze the single microlensing event OGLE-2015-BLG-1482 simultaneously observed from two ground-based surveys and from \textit{Spitzer}. The \textit{Spitzer} data exhibit finite-source effects due to the passage of the lens close to or directly over the surface of the source star as seen from \textit{Spitzer}. Such finite-source effects generally yield measurements of the angular Einstein radius, which when combined with the microlens parallax derived from a comparison between the ground-based and the \textit{Spitzer} light curves, yields the lens mass and lens-source relative parallax. From this analysis, we find that the lens of OGLE-2015-BLG-1482 is a very low-mass star with the mass $0.10 \pm 0.02 \ M_\odot$ or a brown dwarf with the mass $55\pm 9 \ M_{J}$, which are respectively located at $D_{\rm LS} = 0.80 \pm 0.19\ \textrm{kpc}$ and $ D_{\rm LS} = 0.54 \pm 0.08\ \textrm{kpc}$, and thus it is the first isolated low-mass microlens that has been decisively located in the Galactic bulge. The fundamental reason for the degeneracy is that the finite-source effect is seen only in a single data point from \textit{Spitzer} and this single data point gives rise to two solutions for $ρ$. Because the $ρ$ degeneracy can be resolved only by relatively high cadence observations around the peak, while the \textit{Spitzer} cadence is typically $\sim 1\,{\rm day}^{-1}$, we expect that events for which the finite-source effect is seen only in the \textit{Spitzer} data may frequently exhibit this $ρ$ degeneracy. For OGLE-2015-BLG-1482, the relative proper motion of the lens and source for the low-mass star is $μ_{\rm rel} = 9.0 \pm 1.9\ \textrm{mas yr$^{-1}$}$, while for the brown dwarf it is $5.5 \pm 0.5\ \textrm{mas yr$^{-1}$}$. Hence, the degeneracy can be resolved within $\sim 10\ \rm yrs$ from direct lens imaging by using next-generation instruments with high spatial resolution. △ Less

Submitted 21 March, 2017; v1 submitted 17 March, 2017; originally announced March 2017.

Comments: 18 pages, 4 figures, accepted for publication in ApJ

Toward a Galactic Distribution of Planets. I. Methodology & Planet Sensitivities of the 2015 High-Cadence Spitzer Microlens Sample

Authors: Wei Zhu, A. Udalski, S. Calchi Novati, S. -J. Chung, Y. K. Jung, Y. -H. Ryu, I. -G. Shin, A. Gould, C. -U. Lee, M. D. Albrow, J. C. Yee, C. Han, K. -H. Hwang, S. -M. Cha, D. -J. Kim, H. -W. Kim, S. -L. Kim, Y. -H. Kim, Y. Lee, B. -G. Park, R. Pogge, R. Poleski, J. Skowron, P. Mroz, M. K. Szymanski , et al. (13 additional authors not shown)

Abstract: We analyze an ensemble of microlensing events from the 2015 Spitzer microlensing campaign, all of which were densely monitored by ground-based high-cadence survey teams. The simultaneous observations from Spitzer and the ground yield measurements of the microlensing parallax vector $π_{\rm E}$, from which compact constraints on the microlens properties are derived, including $\lesssim$25\% uncerta… ▽ More We analyze an ensemble of microlensing events from the 2015 Spitzer microlensing campaign, all of which were densely monitored by ground-based high-cadence survey teams. The simultaneous observations from Spitzer and the ground yield measurements of the microlensing parallax vector $π_{\rm E}$, from which compact constraints on the microlens properties are derived, including $\lesssim$25\% uncertainties on the lens mass and distance. With the current sample, we demonstrate that the majority of microlenses are indeed in the mass range of M dwarfs. The planet sensitivities of all 41 events in the sample are calculated, from which we provide constraints on the planet distribution function. In particular, assuming a planet distribution function that is uniform in $\log{q}$, where $q$ is the planet-to-star mass ratio, we find a $95\%$ upper limit on the fraction of stars that host typical microlensing planets of 49\%, which is consistent with previous studies. Based on this planet-free sample, we develop the methodology to statistically study the Galactic distribution of planets using microlensing parallax measurements. Under the assumption that the planet distributions are the same in the bulge as in the disk, we predict that $\sim$1/3 of all planet detections from the microlensing campaigns with Spitzer should be in the bulge. This prediction will be tested with a much larger sample, and deviations from it can be used to constrain the abundance of planets in the bulge relative to the disk. △ Less

Submitted 3 November, 2017; v1 submitted 18 January, 2017; originally announced January 2017.

Comments: published on Astronomical Journal

Journal ref: 2017, AJ, 154, 210

OGLE-2015-BLG-0196: Ground-based Gravitational Microlens Parallax Confirmed By Space-Based Observation

Authors: C. Han, A. Udalski, A. Gould, Wei Zhu, M. K. Szymański, I. Soszyński, J. Skowron, P. Mróz, R. Poleski, P. Pietrukowicz, S. Kozłowski, K. Ulaczyk, M. Pawlak, J. C. Yee, C. Beichman, S. Calchi Novati, S. Carey, C. Bryden, M. Fausnaugh, B. S. Gaudi, Calen B. Henderson, Y. Shvartzvald, B. Wibking

Abstract: In this paper, we present the analysis of the binary gravitational microlensing event OGLE-2015-BLG-0196. The event lasted for almost a year and the light curve exhibited significant deviations from the lensing model based on the rectilinear lens-source relative motion, enabling us to measure the microlens parallax. The ground-based microlens parallax is confirmed by the data obtained from space-b… ▽ More In this paper, we present the analysis of the binary gravitational microlensing event OGLE-2015-BLG-0196. The event lasted for almost a year and the light curve exhibited significant deviations from the lensing model based on the rectilinear lens-source relative motion, enabling us to measure the microlens parallax. The ground-based microlens parallax is confirmed by the data obtained from space-based microlens observations using the {\it Spitzer} telescope. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined up to the two-fold degeneracy: $u_0<0$ and $u_0>0$ solutions caused by the well-known "ecliptic" degeneracy. It is found that the binary lens is composed of two M dwarf stars with similar masses $M_1=0.38\pm 0.04\ M_\odot$ ($0.50\pm 0.05\ M_\odot)$ and $M_2=0.38\pm 0.04\ M_\odot$ ($0.55\pm 0.06\ M_\odot$) and the distance to the lens is $D_{\rm L}=2.77\pm 0.23$ kpc ($3.30\pm 0.29$ kpc). Here the physical parameters out and in the parenthesis are for the $u_0<0$ and $u_0>0$ solutions, respectively. △ Less

Submitted 8 November, 2016; originally announced November 2016.

Comments: 7 pages, 3 tables, 6 figures

OGLE-2015-BLG-0479LA,B: Binary Gravitational Microlens Characterized by Simultaneous Ground-based and Space-based Observation

Authors: C. Han, A. Udalski, A. Gould, Wei Zhu, R. A. Street, J. C. Yee, C. Beichman, C. Bryden, S. Calchi Novati, S. Carey, M. Fausnaugh, B. S. Gaudi, Calen B. Henderson, Y. Shvartzvald, B. Wibking, M. K. Szymański, I. Soszyński, J. Skowron, P. Mróz, R. Poleski, P. Pietrukowicz, S. Kozłowski, K. Ulaczyk, Ł. Wyrzykowski, M. Pawlak , et al. (38 additional authors not shown)

Abstract: We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by paralla… ▽ More We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by parallax effects. From modeling the light curves, we measure the space-based microlens parallax. Combined with the angular Einstein radius measured by analyzing the caustic crossings, we determine the mass and distance of the lens. We find that the lens is a binary composed of two G-type stars with masses $\sim 1.0\ M_\odot$ and $\sim 0.9\ M_\odot$ located at a distance $\sim 3$ kpc. In addition, we are able to constrain the complete orbital parameters of the lens thanks to the precise measurement of the microlens parallax derived from the joint analysis. In contrast to the binary event OGLE-2014-BLG-1050, which was also observed by {\it Spitzer}, we find that the interpretation of OGLE-2015-BLG-0479 does not suffer from the degeneracy between $(\pm,\pm)$ and $(\pm,\mp)$ solutions, confirming that the four-fold parallax degeneracy in single-lens events collapses into the two-fold degeneracy for the general case of binary-lens events. The location of the blend in the color-magnitude diagram is consistent with the lens properties, suggesting that the blend is the lens itself. The blend is bright enough for spectroscopy and thus this possibility can be checked from future follow-up observations. △ Less

Submitted 30 June, 2016; originally announced June 2016.

Comments: 9 pages, 4 figures

First simultaneous microlensing observations by two space telescopes: $Spitzer$ & $Swift$ reveal a brown dwarf in event OGLE-2015-BLG-1319

Authors: Y. Shvartzvald, Z. Li, A. Udalski, A. Gould, T. Sumi, R. A. Street, S. Calchi Novati, M. Hundertmark, V. Bozza, C. Beichman, G. Bryden, S. Carey, J. Drummond, M. Fausnaugh, B. S. Gaudi, C. B. Henderson, T. G. Tan, B. Wibking, R. W. Pogge, J. C. Yee, W. Zhu, Y. Tsapras, E. Bachelet, M. Dominik, D. M. Bramich , et al. (68 additional authors not shown)

Abstract: Simultaneous observations of microlensing events from multiple locations allow for the breaking of degeneracies between the physical properties of the lensing system, specifically by exploring different regions of the lens plane and by directly measuring the "microlens parallax". We report the discovery of a 30-55$M_J$ brown dwarf orbiting a K dwarf in microlensing event OGLE-2015-BLG-1319. The sy… ▽ More Simultaneous observations of microlensing events from multiple locations allow for the breaking of degeneracies between the physical properties of the lensing system, specifically by exploring different regions of the lens plane and by directly measuring the "microlens parallax". We report the discovery of a 30-55$M_J$ brown dwarf orbiting a K dwarf in microlensing event OGLE-2015-BLG-1319. The system is located at a distance of $\sim$5 kpc toward the Galactic bulge. The event was observed by several ground-based groups as well as by $Spitzer$ and $Swift$, allowing the measurement of the physical properties. However, the event is still subject to an 8-fold degeneracy, in particular the well-known close-wide degeneracy, and thus the projected separation between the two lens components is either $\sim$0.25 AU or $\sim$45 AU. This is the first microlensing event observed by $Swift$, with the UVOT camera. We study the region of microlensing parameter space to which $Swift$ is sensitive, finding that while for this event $Swift$ could not measure the microlens parallax with respect to ground-based observations, it can be important for other events. Specifically, for detecting nearby brown dwarfs and free-floating planets in high magnification events. △ Less

Submitted 7 June, 2016; originally announced June 2016.

Comments: 28 pages, 6 figures, 2 tables. Submitted to ApJ

Spitzer Observations of OGLE-2015-BLG-1212 Reveal a New Path to Breaking Strong Microlens Degeneracies

Authors: V. Bozza, Y. Shvartzvald, A. Udalski, S. Calchi Novati, I. A. Bond, C. Han, M. Hundertmark, R. Poleski, M. Pawlak, M. K. Szymański, J. Skowron, P. Mróz, S. Kozłowski, Ł. Wyrzykowski, P. Pietrukowicz, I. Soszyński, K. Ulaczyk, C. Beichman, G. Bryden, S. Carey, M. Fausnaugh, B. S. Gaudi, A. Gould, C. B. Henderson, R. W. Pogge , et al. (76 additional authors not shown)

Abstract: Spitzer microlensing parallax observations of OGLE-2015-BLG-1212 decisively breaks a degeneracy between planetary and binary solutions that is somewhat ambiguous when only ground-based data are considered. Only eight viable models survive out of an initial set of 32 local minima in the parameter space. These models clearly indicate that the lens is a stellar binary system possibly located within t… ▽ More Spitzer microlensing parallax observations of OGLE-2015-BLG-1212 decisively breaks a degeneracy between planetary and binary solutions that is somewhat ambiguous when only ground-based data are considered. Only eight viable models survive out of an initial set of 32 local minima in the parameter space. These models clearly indicate that the lens is a stellar binary system possibly located within the bulge of our Galaxy, ruling out the planetary alternative. We argue that several types of discrete degeneracies can be broken via such space-based parallax observations. △ Less

Submitted 10 February, 2016; v1 submitted 7 January, 2016; originally announced January 2016.

Comments: 27 pages, 6 figures, accepted by ApJ

Spitzer Microlensing Program as a Probe for Globular Cluster Planets. Analysis of OGLE-2015-BLG-0448

Authors: Radosław Poleski, Wei Zhu, Grant W. Christie, Andrzej Udalski, Andrew Gould, Etienne Bachelet, Jesper Skottfelt, Sebastiano Calchi Novati, M. K. Szymański, I. Soszyński, G. Pietrzyński, Ł. Wyrzykowski, K. Ulaczyk, P. Pietrukowicz, Szymon Kozłowski, J. Skowron, P. Mróz, M. Pawlak, C. Beichman, G. Bryden, S. Carey, M. Fausnaugh, B. S. Gaudi, C. B. Henderson, R. W. Pogge , et al. (61 additional authors not shown)

Abstract: The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($μ_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the sou… ▽ More The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($μ_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using OGLE photometry is consistent with the value found based on the light curve displacement between Earth and Spitzer. △ Less

Submitted 29 October, 2019; v1 submitted 28 December, 2015; originally announced December 2015.

Comments: ApJ submitted, 12 pages, 8 figures, 2 tabels

Mass Measurements of Isolated Objects from Space-based Microlensing

Authors: Wei Zhu, S. Calchi Novati, A. Gould, A. Udalski, C. Han, Y. Shvartzvald, C. Ranc, U. G. Jorgensen, R. Poleski, V. Bozza, C. Beichman, G. Bryden, S. Carey, B. S. Gaudi, C. B. Henderson, R. W. Pogge, I. Porritt, B. Wibking, J. C. Yee, M. Pawlak, M. K. Szymanski, J. Skowron, P. Mroz, S. Kozlowski, L. Wyrzykowski , et al. (56 additional authors not shown)

Abstract: We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$… ▽ More We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$ $M_{\rm J}$ brown dwarf at $5.9\pm1.0$ kpc. The lens of of the second event, OGLE-2015-BLG-0763, is a $0.50\pm0.04$ $M_\odot$ star at $6.9\pm1.0$ kpc. We show that the probability to definitively measure the mass of isolated microlenses is dramatically increased once simultaneous ground- and space-based observations are conducted. △ Less

Submitted 21 April, 2016; v1 submitted 7 October, 2015; originally announced October 2015.

Comments: 10 papers, 4 figures, 2 tables; ApJ in press

Spitzer IRAC Photometry for Time Series in Crowded Fields

Authors: S. Calchi Novati, A. Gould, J. C. Yee, C. Beichman, G. Bryden, S. Carey, M. Fausnaugh, B. S. Gaudi, C. B. Henderson, R. W. Pogge, Y. Shvartzvald, B. Wibking, W. Zhu, A. Udalski, R. Poleski, M. Pawlak, M. K. Szymański, J. Skowron, P. Mróz, S. Kozłowski, Ł. Wyrzykowski, P. Pietrukowicz, G. Pietrzyński, I. Soszyński, K. Ulaczyk

Abstract: We develop a new photometry algorithm that is optimized for $Spitzer$ time series in crowded fields and that is particularly adapted to faint and/or heavily blended targets. We apply this to the 170 targets from the 2015 $Spitzer$ microlensing campaign and present the results of three variants of this algorithm in an online catalog. We present detailed accounts of the application of this algorithm… ▽ More We develop a new photometry algorithm that is optimized for $Spitzer$ time series in crowded fields and that is particularly adapted to faint and/or heavily blended targets. We apply this to the 170 targets from the 2015 $Spitzer$ microlensing campaign and present the results of three variants of this algorithm in an online catalog. We present detailed accounts of the application of this algorithm to two difficult cases, one very faint and the other very crowded. Several of $Spitzer$'s instrumental characteristics that drive the specific features of this algorithm are shared by $Kepler$ and $WFIRST$, implying that these features may prove to be a useful starting point for algorithms designed for microlensing campaigns by these other missions. △ Less

Submitted 13 October, 2015; v1 submitted 31 August, 2015; originally announced September 2015.

Comments: accepted for publication in The Astrophysical Journal Supplement, online catalog available at http://www.astronomy.ohio-state.edu/Spitzer2015/

Journal ref: ApJ, 814, 92 (2015)

Spitzer Microlens Measurement of a Massive Remnant in a Well-Separated Binary

Authors: Y. Shvartzvald, A. Udalski, A. Gould, C. Han, V. Bozza, M. Friedmann, M. Hundertmark, C. Beichman, G. Bryden, S. Calchi Novati, S. Carey, M. Fausnaugh, B. S. Gaudi, C. B. Henderson, T. Kerr, R. W. Pogge, W. Varricatt, B. Wibking, J. C. Yee, W. Zhu, R. Poleski, M. Pawlak, M. K. Szymański, J. Skowron, P. Mróz , et al. (41 additional authors not shown)

Abstract: We report the detection and mass measurement of a binary lens OGLE-2015-BLG-1285La,b, with the more massive component having $M_1>1.35\,M_\odot$ (80% probability). A main-sequence star in this mass range is ruled out by limits on blue light, meaning that a primary in this mass range must be a neutron star or black hole. The system has a projected separation $r_\perp= 6.1\pm 0.4\,{\rm AU}$ and lies… ▽ More We report the detection and mass measurement of a binary lens OGLE-2015-BLG-1285La,b, with the more massive component having $M_1>1.35\,M_\odot$ (80% probability). A main-sequence star in this mass range is ruled out by limits on blue light, meaning that a primary in this mass range must be a neutron star or black hole. The system has a projected separation $r_\perp= 6.1\pm 0.4\,{\rm AU}$ and lies in the Galactic bulge. These measurements are based on the "microlens parallax" effect, i.e., comparing the microlensing light curve as seen from $Spitzer$, which lay at $1.25\,{\rm AU}$ projected from Earth, to the light curves from four ground-based surveys, three in the optical and one in the near infrared. Future adaptive optics imaging of the companion by 30m class telescopes will yield a much more accurate measurement of the primary mass. This discovery both opens the path and defines the challenges to detecting and characterizing black holes and neutron stars in wide binaries, with either dark or luminous companions. In particular, we discuss lessons that can be applied to future $Spitzer$ and $Kepler$ K2 microlensing parallax observations. △ Less

Submitted 18 December, 2015; v1 submitted 26 August, 2015; originally announced August 2015.

Comments: ApJ. accepted, 34 pages, 7 figures, 2 tables