Showing 1–48 of 48 results for author: Bassett, R

Uncertainty Quantification Using Ensemble Learning and Monte Carlo Sampling for Performance Prediction and Monitoring in Cell Culture Processes

Authors: Thanh Tung Khuat, Robert Bassett, Ellen Otte, Bogdan Gabrys

Abstract: Biopharmaceutical products, particularly monoclonal antibodies (mAbs), have gained prominence in the pharmaceutical market due to their high specificity and efficacy. As these products are projected to constitute a substantial portion of global pharmaceutical sales, the application of machine learning models in mAb development and manufacturing is gaining momentum. This paper addresses the critica… ▽ More Biopharmaceutical products, particularly monoclonal antibodies (mAbs), have gained prominence in the pharmaceutical market due to their high specificity and efficacy. As these products are projected to constitute a substantial portion of global pharmaceutical sales, the application of machine learning models in mAb development and manufacturing is gaining momentum. This paper addresses the critical need for uncertainty quantification in machine learning predictions, particularly in scenarios with limited training data. Leveraging ensemble learning and Monte Carlo simulations, our proposed method generates additional input samples to enhance the robustness of the model in small training datasets. We evaluate the efficacy of our approach through two case studies: predicting antibody concentrations in advance and real-time monitoring of glucose concentrations during bioreactor runs using Raman spectra data. Our findings demonstrate the effectiveness of the proposed method in estimating the uncertainty levels associated with process performance predictions and facilitating real-time decision-making in biopharmaceutical manufacturing. This contribution not only introduces a novel approach for uncertainty quantification but also provides insights into overcoming challenges posed by small training datasets in bioprocess development. The evaluation demonstrates the effectiveness of our method in addressing key challenges related to uncertainty estimation within upstream cell cultivation, illustrating its potential impact on enhancing process control and product quality in the dynamic field of biopharmaceuticals. △ Less

Submitted 3 September, 2024; originally announced September 2024.

ACM Class: A.1; J.3; I.2.0; I.2.6; I.2.m; I.5.4

Optimal Design of Resolvent Splitting Algorithms

Authors: Robert L Bassett, Peter Barkley

Abstract: In this paper, we introduce a novel semidefinite programming framework for designing custom frugal resolvent splitting algorithms which find a zero in the sum of n monotone operators. This framework features a number of design choices which facilitate creating resolvent splitting algorithms with specific communication structure. We illustrate these design choices using a variety of constraint sets… ▽ More In this paper, we introduce a novel semidefinite programming framework for designing custom frugal resolvent splitting algorithms which find a zero in the sum of n monotone operators. This framework features a number of design choices which facilitate creating resolvent splitting algorithms with specific communication structure. We illustrate these design choices using a variety of constraint sets and objective functions, as well as the use of a mixed-integer SDP to minimize time per iteration or required number of communications between the resolvents which define the splitting. Using the Performance Estimation Problem (PEP) framework, we provide parameter selections, such as step size, which for high dimensional problems provide optimal contraction factors in the algorithms. Among the algorithm design choices we introduce, we provide a characterization of algorithm designs which provide minimal convergence time given structural properties of the monotone operators, resolvent computation times, and communication latencies. △ Less

Submitted 3 September, 2024; v1 submitted 22 July, 2024; originally announced July 2024.

MSC Class: 49M37; 65K05

Galaxy 3D Shape Recovery using Mixture Density Network

Authors: Suk Yee Yong, K. E. Harborne, Caroline Foster, Robert Bassett, Gregory B. Poole, Mitchell Cavanagh

Abstract: Since the turn of the century, astronomers have been exploiting the rich information afforded by combining stellar kinematic maps and imaging in an attempt to recover the intrinsic, three-dimensional (3D) shape of a galaxy. A common intrinsic shape recovery method relies on an expected monotonic relationship between the intrinsic misalignment of the kinematic and morphological axes and the triaxia… ▽ More Since the turn of the century, astronomers have been exploiting the rich information afforded by combining stellar kinematic maps and imaging in an attempt to recover the intrinsic, three-dimensional (3D) shape of a galaxy. A common intrinsic shape recovery method relies on an expected monotonic relationship between the intrinsic misalignment of the kinematic and morphological axes and the triaxiality parameter. Recent studies have, however, cast doubt about underlying assumptions relating shape and intrinsic kinematic misalignment. In this work, we aim to recover the 3D shape of individual galaxies using their projected stellar kinematic and flux distributions using a supervised machine learning approach with mixture density network (MDN). Using a mock dataset of the EAGLE hydrodynamical cosmological simulation, we train the MDN model for a carefully selected set of common kinematic and photometric parameters. Compared to previous methods, we demonstrate potential improvements achieved with the MDN model to retrieve the 3D galaxy shape along with the uncertainties, especially for prolate and triaxial systems. We make specific recommendations for recovering galaxy intrinsic shapes relevant for current and future integral field spectroscopic galaxy surveys. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: Accepted for publication in PASA. 18 pages, 12 figures, 2 tables

Journal ref: Publ. Astron. Soc. Aust. 41 (2024) e033

Adversarial Perturbations of Physical Signals

Authors: Robert L. Bassett, Austin Van Dellen, Anthony P. Austin

Abstract: We investigate the vulnerability of computer-vision-based signal classifiers to adversarial perturbations of their inputs, where the signals and perturbations are subject to physical constraints. We consider a scenario in which a source and interferer emit signals that propagate as waves to a detector, which attempts to classify the source by analyzing the spectrogram of the signal it receives usi… ▽ More We investigate the vulnerability of computer-vision-based signal classifiers to adversarial perturbations of their inputs, where the signals and perturbations are subject to physical constraints. We consider a scenario in which a source and interferer emit signals that propagate as waves to a detector, which attempts to classify the source by analyzing the spectrogram of the signal it receives using a pre-trained neural network. By solving PDE-constrained optimization problems, we construct interfering signals that cause the detector to misclassify the source even though the perturbations to the spectrogram of the received signal are nearly imperceptible. Though such problems can have millions of decision variables, we introduce methods to solve them efficiently. Our experiments demonstrate that one can compute effective and physically realizable adversarial perturbations for a variety of machine learning models under various physical conditions. △ Less

Submitted 26 February, 2024; originally announced February 2024.

MSC Class: 65K05; 90C30; 65M60

Applications of Machine Learning in Biopharmaceutical Process Development and Manufacturing: Current Trends, Challenges, and Opportunities

Authors: Thanh Tung Khuat, Robert Bassett, Ellen Otte, Alistair Grevis-James, Bogdan Gabrys

Abstract: While machine learning (ML) has made significant contributions to the biopharmaceutical field, its applications are still in the early stages in terms of providing direct support for quality-by-design based development and manufacturing of biopharmaceuticals, hindering the enormous potential for bioprocesses automation from their development to manufacturing. However, the adoption of ML-based mode… ▽ More While machine learning (ML) has made significant contributions to the biopharmaceutical field, its applications are still in the early stages in terms of providing direct support for quality-by-design based development and manufacturing of biopharmaceuticals, hindering the enormous potential for bioprocesses automation from their development to manufacturing. However, the adoption of ML-based models instead of conventional multivariate data analysis methods is significantly increasing due to the accumulation of large-scale production data. This trend is primarily driven by the real-time monitoring of process variables and quality attributes of biopharmaceutical products through the implementation of advanced process analytical technologies. Given the complexity and multidimensionality of a bioproduct design, bioprocess development, and product manufacturing data, ML-based approaches are increasingly being employed to achieve accurate, flexible, and high-performing predictive models to address the problems of analytics, monitoring, and control within the biopharma field. This paper aims to provide a comprehensive review of the current applications of ML solutions in a bioproduct design, monitoring, control, and optimisation of upstream, downstream, and product formulation processes. Finally, this paper thoroughly discusses the main challenges related to the bioprocesses themselves, process data, and the use of machine learning models in biopharmaceutical process development and manufacturing. Moreover, it offers further insights into the adoption of innovative machine learning methods and novel trends in the development of new digital biopharma solutions. △ Less

Submitted 15 October, 2023; originally announced October 2023.

Comments: 155 pages

ACM Class: A.1; J.3; I.2.0; I.2.6; I.2.m; I.5.4

Signal Subspace Methods Which Are Robust to Impulsive Noise

Authors: Robert L. Bassett, Micah Y. Oh

Abstract: We consider the problem of estimating a signal subspace in the presence of interference that contaminates some proportion of the received observations. Our emphasis is on detecting the contaminated observations so that the signal subspace can be estimated with the contaminated observations discarded. To this end, we employ a signal model which explicitly includes an interference term that is disti… ▽ More We consider the problem of estimating a signal subspace in the presence of interference that contaminates some proportion of the received observations. Our emphasis is on detecting the contaminated observations so that the signal subspace can be estimated with the contaminated observations discarded. To this end, we employ a signal model which explicitly includes an interference term that is distinct from environmental noise. To detect when the interference term is nonzero, we estimate the interference term using an optimization problem with a sparsity-inducing group SLOPE penalty which accounts for simultaneous sparsity across all channels of the multichannel signal. We propose an iterative algorithm which efficiently computes the observations estimated to contain interference. Theoretical support for the accuracy of our interference estimator is provided by bounding its false discovery rate, the expected proportion of uncontaminated observations among those estimated to be contaminated. Finally, we demonstrate the empirical performance of our contributions in a number of simulated experiments. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Exponential Family Trend Filtering on Lattices

Authors: Veeranjaneyulu Sadhanala, Robert Bassett, James Sharpnack, Daniel J. McDonald

Abstract: Trend filtering is a modern approach to nonparametric regression that is more adaptive to local smoothness than splines or similar basis procedures. Existing analyses of trend filtering focus on estimating a function corrupted by homoskedastic Gaussian noise, but our work extends this technique to general exponential family distributions. This extension is motivated by the need to study massive, g… ▽ More Trend filtering is a modern approach to nonparametric regression that is more adaptive to local smoothness than splines or similar basis procedures. Existing analyses of trend filtering focus on estimating a function corrupted by homoskedastic Gaussian noise, but our work extends this technique to general exponential family distributions. This extension is motivated by the need to study massive, gridded climate data derived from polar-orbiting satellites. We present algorithms tailored to large problems, theoretical results for general exponential family likelihoods, and principled methods for tuning parameter selection without excess computation. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 53 pages; 6 figures; 3 tables

A Cautionary Tale of LyC Escape Fraction Estimates from High Redshift Galaxies

Authors: R. Bassett, E. V. Ryan-Weber, J. Cooke, U. Mestric, L. J. Prichard, M. Rafelski, I. Iwata, M. Sawicki, S. Gwyn, S. Arnouts

Abstract: Measuring the escape fraction, $f_{\rm esc}$, of ionizing, Lyman Continuum (LyC) radiation is key to our understanding of the process of cosmic reionization. In this paper we provide a methodology for recovering the posterior probability distribution of the LyC escape fraction, $f_{\rm esc}^{\rm PDF}$, considering both the observational uncertainties and ensembles of simulated transmission functio… ▽ More Measuring the escape fraction, $f_{\rm esc}$, of ionizing, Lyman Continuum (LyC) radiation is key to our understanding of the process of cosmic reionization. In this paper we provide a methodology for recovering the posterior probability distribution of the LyC escape fraction, $f_{\rm esc}^{\rm PDF}$, considering both the observational uncertainties and ensembles of simulated transmission functions through the intergalactic medium (IGM). We present an example of this method applied to a VUDS galaxy at $z=3.64$ and find $f_{\rm esc}^{\rm PDF}$ = 0.51$^{+0.33}_{-0.34}$ and compare this to the values computed assuming averaged IGM transmission with and without consideration of detection bias along average sightlines yielding $f_{\rm esc}^{\langle T \rangle}$ = 1.40$^{+0.80}_{-0.42}$, and $f_{\rm esc}^{\rm bias}$ = 0.82$^{+0.33}_{-0.16}$. Our results highlight the limitations of methods assuming average, smooth transmission functions. We also present MOSFIRE data for a sample of seven LyC candidates selected based on photometric redshifts at $z > 3.4$, but find that all seven have overestimated photometric redshifts by $Δz \sim 0.2$ making them unsuitable for LyC measurements. This results likely due to a bias induced by our selection criteria. △ Less

Submitted 6 February, 2022; originally announced February 2022.

Comments: 12 pages, 9 figues, 2 tables, accepted for publication in MNRAS

Lyman Continuum Galaxy Candidates in COSMOS

Authors: Laura J. Prichard, Marc Rafelski, Jeff Cooke, Uros Mestric, Robert Bassett, Emma V. Ryan-Weber, Ben Sunnquist, Anahita Alavi, Nimish Hathi, Xin Wang, Mitchell Revalski, Varun Bajaj, John M. O'Meara, Lee Spitler

Abstract: Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$,… ▽ More Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$, $\sim29\text{--}30$~mag depth at 5$σ$), we attempt to identify a sample of Lyman continuum galaxies (LCGs). These are individual sources that emit ionizing flux below the Lyman break ($<912~\textÅ$). This population would allow us to constrain cosmic reionization parameters such as the number density and escape fraction ($f_{\rm esc}$) of ionizing sources. We compile a comprehensive parent sample that does not rely on the Lyman-break technique for redshifts. We present three new spectroscopic candidates at $z\sim3.7\text{--}4.4$, and 32 new photometric candidates. The high-resolution multi-band HST imaging and new Keck/Low Resolution Imaging Spectrometer (LRIS) redshifts make these promising spectroscopic LCG candidates. Using both a traditional and probabilistic approach, we find the most likely $f_{\rm esc}$ values for the three spectroscopic LCG candidates are $>100\%$, and therefore not physical. We are unable to confirm the true nature of these sources with the best available imaging and direct blue Keck/LRIS spectroscopy. More spectra, especially from the new class of 30 m telescopes, will be required to build a statistical sample of LCGs to place firm observational constraints on cosmic reionization. △ Less

Submitted 14 February, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ. 38 pages (28 in main body), 14 figures

Upper limits on the escape fraction of ionizing radiation from galaxies at $2\lesssim z < 6$

Authors: U. Meštrić, E. V. Ryan-Weber, J. Cooke, R. Bassett, L. J. Prichard, M. Rafelski

Abstract: In this work, we investigate upper limits on the global escape fraction of ionizing photons ($f_{\rm esc/global}^{\rm abs}$) from a sample of galaxies probed for Lyman-continuum (LyC) emission characterized as non-LyC and LyC leakers. We present a sample of 9 clean non-contaminated (by low redshift interlopers, CCD problems and internal reflections of the instrument) galaxies which do not show sig… ▽ More In this work, we investigate upper limits on the global escape fraction of ionizing photons ($f_{\rm esc/global}^{\rm abs}$) from a sample of galaxies probed for Lyman-continuum (LyC) emission characterized as non-LyC and LyC leakers. We present a sample of 9 clean non-contaminated (by low redshift interlopers, CCD problems and internal reflections of the instrument) galaxies which do not show significant ($>$ $3σ$) LyC flux between 880Å $<λ_{rest}<$ 910Å. The 9 galaxy stacked spectrum reveals no significant LyC flux with an upper limit of $f_{\rm esc}^{\rm abs} \leq 0.06$. In the next step of our analysis, we join all estimates of $f_{\rm esc}^{\rm abs}$ upper limits derived from different samples of $2\lesssim z < 6$ galaxies from the literature reported in last $\sim$20 years and include the sample presented in this work. We find the $f_{\rm esc}^{\rm abs}$ upper limit $\leq$ 0.084 for the galaxies recognized as non-LyC leakers. After including all known detections from literature $f_{\rm esc/global}^{\rm abs}$ upper limit $\leq$ 0.088 for all galaxies examined for LyC flux. Furthermore, $f_{\rm esc}^{\rm abs}$ upper limits for different groups of galaxies indicate that the strongest LyC emitters could be galaxies classified as Lyman alpha emitters. We also discuss the possible existence of a correlation among the observed flux density ratio $(F_ν^{LyC}/F_ν^{UV})_{\rm obs}$ and Lyman alpha equivalent width EW(Ly$α)$, where we confirm the existence of moderately significant correlation among galaxies classified as non-LyC leakers. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 17 pages, 5 Figures

DeepFake Detection with Inconsistent Head Poses: Reproducibility and Analysis

Authors: Kevin Lutz, Robert Bassett

Abstract: Applications of deep learning to synthetic media generation allow the creation of convincing forgeries, called DeepFakes, with limited technical expertise. DeepFake detection is an increasingly active research area. In this paper, we analyze an existing DeepFake detection technique based on head pose estimation, which can be applied when fake images are generated with an autoencoder-based face swa… ▽ More Applications of deep learning to synthetic media generation allow the creation of convincing forgeries, called DeepFakes, with limited technical expertise. DeepFake detection is an increasingly active research area. In this paper, we analyze an existing DeepFake detection technique based on head pose estimation, which can be applied when fake images are generated with an autoencoder-based face swap. Existing literature suggests that this method is an effective DeepFake detector, and its motivating principles are attractively simple. With an eye towards using these principles to develop new DeepFake detectors, we conduct a reproducibility study of the existing method. We conclude that its merits are dramatically overstated, despite its celebrated status. By investigating this discrepancy we uncover a number of important and generalizable insights related to facial landmark detection, identity-agnostic head pose estimation, and algorithmic bias in DeepFake detectors. Our results correct the current literature's perception of state of the art performance for DeepFake detection. △ Less

Submitted 28 August, 2021; originally announced August 2021.

The Maximal Eigengap Estimator for Acoustic Vector-Sensor Processing

Authors: Robert Bassett, Jacob Foster, Kay L. Gemba, Paul Leary, Kevin B. Smith

Abstract: This paper introduces the maximal eigengap estimator for finding the direction of arrival of a wideband acoustic signal using a single vector-sensor. We show that in this setting narrowband cross-spectral density matrices can be combined in an optimal weighting that approximately maximizes signal-to-noise ratio across a wide frequency band. The signal subspace resulting from this optimal combinati… ▽ More This paper introduces the maximal eigengap estimator for finding the direction of arrival of a wideband acoustic signal using a single vector-sensor. We show that in this setting narrowband cross-spectral density matrices can be combined in an optimal weighting that approximately maximizes signal-to-noise ratio across a wide frequency band. The signal subspace resulting from this optimal combination of narrowband power matrices defines the maximal eigengap estimator. We discuss the advantages of the maximal eigengap estimator over competing methods, and demonstrate its utility in a real-data application using signals collected in 2019 from an acoustic vector-sensor deployed in the Monterey Bay. △ Less

Submitted 7 July, 2021; v1 submitted 21 April, 2021; originally announced April 2021.

CRKSPH-compatible discretization of the SUPG and SAAF transport equations

Authors: Brody R. Bassett, J. Michael Owen

Abstract: The self-adjoint angular flux and streamline-upwind Petrov-Galerkin transport equations are discretized using reproducing kernels with the collocation method to produce a discretization that is compatible with conservative reproducing kernel smoothed particle hydrodynamics. A novel second derivative is derived for the diffusion-like term in the self-adjoint angular flux equation. The resulting equ… ▽ More The self-adjoint angular flux and streamline-upwind Petrov-Galerkin transport equations are discretized using reproducing kernels with the collocation method to produce a discretization that is compatible with conservative reproducing kernel smoothed particle hydrodynamics. A novel second derivative is derived for the diffusion-like term in the self-adjoint angular flux equation. The resulting equations involve only evaluations of kernels and physical data at the nodal centers. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: Submitted to the International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C) 2021

Report number: LLNL-CONF-819823

Systematic difference between ionized and molecular gas velocity dispersion in $z\sim1-2$ disks and local analogues

Authors: M. Girard, D. B. Fisher, A. D. Bolatto, R. Abraham, R. Bassett, K. Glazebrook, R. Herrera-Camus, E. Jiménez, L. Lenkić, D. Obreschkow

Abstract: We compare the molecular and ionized gas velocity dispersion of 9 nearby turbulent disks, analogues to high-redshift galaxies, from the DYNAMO sample using new ALMA and GMOS/Gemini observations. We combine our sample with 12 galaxies at $z\sim $0.5-2.5 from the literature. We find that the resolved velocity dispersion is systematically lower by a factor $2.45\pm0.38$ for the molecular gas compared… ▽ More We compare the molecular and ionized gas velocity dispersion of 9 nearby turbulent disks, analogues to high-redshift galaxies, from the DYNAMO sample using new ALMA and GMOS/Gemini observations. We combine our sample with 12 galaxies at $z\sim $0.5-2.5 from the literature. We find that the resolved velocity dispersion is systematically lower by a factor $2.45\pm0.38$ for the molecular gas compared to the ionized gas, after correcting for thermal broadening. This offset is constant within the galaxy disks and indicates the co-existence of a thin molecular and thick ionized gas disks. This result has a direct impact on the Toomre $Q$ and pressure derived in galaxies. We obtain pressures $\sim0.22$ dex lower on average when using the molecular gas velocity dispersion, $σ_{0,mol}$. We find that $σ_{0,mol}$ increases with gas fraction and star formation rate. We also obtain an increase with redshift and show that the EAGLE and FIRE simulations overall overestimate $σ_{0,mol}$ at high redshift. Our results suggest that efforts to compare the kinematics of gas using ionized gas as a proxy for the total gas may overestimate the velocity dispersion by a significant amount in galaxies at the peak of cosmic star formation. When using the molecular gas as a tracer, our sample is not consistent with predictions from constant efficiency star formation models, even when including transport as a source of turbulence. Feedback models with variable star formation efficiency, $ε_{ff}$, and/or feedback efficiency, $p_*/m_*$, better predict our observations. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in ApJ. 19 pages, 6 figures

IGM Transmission Bias for $z$ $\geq$ 2.9 Lyman Continuum Detected Galaxies

Authors: R. Bassett, E. V. Ryan-Weber, J. Cooke, U. Meštrić, K. Kakiichi, L. Prichard, M. Rafelski

Abstract: Understanding the relationship between the underlying escape fraction of Lyman continuum (LyC) photons ($f_{\rm esc}$) emitted by galaxies and measuring the distribution of observed $f_{\rm esc}$ values at high redshift is fundamental to the interpretation of the reionization process. In this paper we perform a statistical exploration of the attenuation of LyC photons by neutral hydrogen in the in… ▽ More Understanding the relationship between the underlying escape fraction of Lyman continuum (LyC) photons ($f_{\rm esc}$) emitted by galaxies and measuring the distribution of observed $f_{\rm esc}$ values at high redshift is fundamental to the interpretation of the reionization process. In this paper we perform a statistical exploration of the attenuation of LyC photons by neutral hydrogen in the intergalactic medium using ensembles of simulated transmission functions. We show that LyC detected galaxies are more likely to be found in sightlines with higher-than-average transmission of LyC photons. This means that adopting a mean transmission at a given redshift leads to an overestimate of the true $f_{\rm esc}$ for LyC detected galaxies. We note, however, that mean values are appropriate for $f_{\rm esc}$ estimates of larger parent samples that include LyC non-detected galaxies. We quantify this IGM transmission bias for LyC detections in photometric and spectroscopic surveys in the literature and show that the bias is stronger for both shallower observations and for fainter parent samples (i.e. Lyman $α$ emitters versus Lyman break galaxies). We also explore the effects of varying the underlying probability distribution function (PDF) of $f_{\rm esc}$ on recovered values, showing that the underlying $f_{\rm esc}$ PDF may depend on sample selection by comparing with observational surveys. This work represents a first step in improved interpretation of LyC detections in the context of understanding $f_{\rm esc}$ from high redshift galaxies. △ Less

Submitted 3 January, 2021; originally announced January 2021.

Comments: 20 pages, 17 figures, resubmitted to MNRAS after minor revision

The MAGPI Survey -- science goals, design, observing strategy, early results and theoretical framework

Authors: C. Foster, J. T. Mendel, C. D. P. Lagos, E. Wisnioski, T. Yuan, F. D'Eugenio, T. M. Barone, K. E. Harborne, S. P. Vaughan, F. Schulze, R. -S. Remus, A. Gupta, F. Collacchioni, D. J. Khim, P. Taylor, R. Bassett, S. M. Croom, R. M. McDermid, A. Poci, A. J. Battisti, J. Bland-Hawthorn, S. Bellstedt, M. Colless, L. J. M. Davies, C. Derkenne , et al. (18 additional authors not shown)

Abstract: We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive… ▽ More We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive optics aided MUSE observations of fields selected from the GAMA survey, providing a wealth of publicly available ancillary multi-wavelength data. With these data, MAGPI will map the kinematic and chemical properties of stars and ionised gas for a sample of 60 massive (> 7 x 10^10 M_Sun) central galaxies at 0.25 < z <0.35 in a representative range of environments (isolated, groups and clusters). The spatial resolution delivered by MUSE with Ground Layer Adaptive Optics (GLAO, 0.6-0.8 arcsec FWHM) will facilitate a direct comparison with Integral Field Spectroscopy surveys of the nearby Universe, such as SAMI and MaNGA, and at higher redshifts using adaptive optics, e.g. SINS. In addition to the primary (central) galaxy sample, MAGPI will deliver resolved and unresolved spectra for as many as 150 satellite galaxies at 0.25 < z <0.35, as well as hundreds of emission-line sources at z < 6. This paper outlines the science goals, survey design, and observing strategy of MAGPI. We also present a first look at the MAGPI data, and the theoretical framework to which MAGPI data will be compared using the current generation of cosmological hydrodynamical simulations including EAGLE, Magneticum, HORIZON-AGN, and Illustris-TNG. Our results show that cosmological hydrodynamical simulations make discrepant predictions in the spatially resolved properties of galaxies at z ~ 0.3. MAGPI observations will place new constraints and allow for tangible improvements in galaxy formation theory. △ Less

Submitted 14 June, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 31 pages, 13 figures, PASA accepted

Color and Edge-Aware Adversarial Image Perturbations

Authors: Robert Bassett, Mitchell Graves, Patrick Reilly

Abstract: Adversarial perturbation of images, in which a source image is deliberately modified with the intent of causing a classifier to misclassify the image, provides important insight into the robustness of image classifiers. In this work we develop two new methods for constructing adversarial perturbations, both of which are motivated by minimizing human ability to detect changes between the perturbed… ▽ More Adversarial perturbation of images, in which a source image is deliberately modified with the intent of causing a classifier to misclassify the image, provides important insight into the robustness of image classifiers. In this work we develop two new methods for constructing adversarial perturbations, both of which are motivated by minimizing human ability to detect changes between the perturbed and source image. The first of these, the Edge-Aware method, reduces the magnitude of perturbations permitted in smooth regions of an image where changes are more easily detected. Our second method, the Color-Aware method, performs the perturbation in a color space which accurately captures human ability to distinguish differences in colors, thus reducing the perceived change. The Color-Aware and Edge-Aware methods can also be implemented simultaneously, resulting in image perturbations which account for both human color perception and sensitivity to changes in homogeneous regions. Because Edge-Aware and Color-Aware modifications exist for many image perturbations techniques, we also focus on computation to demonstrate their potential for use within more complex perturbation schemes. We empirically demonstrate that the Color-Aware and Edge-Aware perturbations we consider effectively cause misclassification, are less distinguishable to human perception, and are as easy to compute as the most efficient image perturbation techniques. Code and demo available at https://github.com/rbassett3/Color-and-Edge-Aware-Perturbations △ Less

Submitted 22 March, 2021; v1 submitted 27 August, 2020; originally announced August 2020.

MSC Class: 68T45; 62F35

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

Authors: Brody R. Bassett, J. Michael Owen

Abstract: The time-dependent radiation transport equation is discretized using the meshless-local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The… ▽ More The time-dependent radiation transport equation is discretized using the meshless-local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. Results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels. △ Less

Submitted 5 August, 2020; originally announced August 2020.

Report number: LLNL-JRNL-813299

One-Step Estimation With Scaled Proximal Methods

Authors: Robert Bassett, Julio Deride

Abstract: We study statistical estimators computed using iterative optimization methods that are not run until completion. Classical results on maximum likelihood estimators (MLEs) assert that a one-step estimator (OSE), in which a single Newton-Raphson iteration is performed from a starting point with certain properties, is asymptotically equivalent to the MLE. We further develop these early-stopping resul… ▽ More We study statistical estimators computed using iterative optimization methods that are not run until completion. Classical results on maximum likelihood estimators (MLEs) assert that a one-step estimator (OSE), in which a single Newton-Raphson iteration is performed from a starting point with certain properties, is asymptotically equivalent to the MLE. We further develop these early-stopping results by deriving properties of one-step estimators defined by a single iteration of scaled proximal methods. Our main results show the asymptotic equivalence of the likelihood-based estimator and various one-step estimators defined by scaled proximal methods. By interpreting OSEs as the last of a sequence of iterates, our results provide insight on scaling numerical tolerance with sample size. Our setting contains scaled proximal gradient descent applied to certain composite models as a special case, making our results applicable to many problems of practical interest. Additionally, our results provide support for the utility of the scaled Moreau envelope as a statistical smoother by interpreting scaled proximal descent as a quasi-Newton method applied to the scaled Moreau envelope. △ Less

Submitted 24 June, 2021; v1 submitted 29 April, 2020; originally announced April 2020.

MSC Class: 62F12; 65K10; 90C53

Outside the Lyman-break box: detecting Lyman continuum emitters at $3.5<z<5.1$ with CLAUDS

Authors: Uroš Meštrić, E. V. Ryan-Weber, J. Cooke, R. Bassett, M. Sawicki, A. L. Faisst, K. Kakiichi, A. K. Inoue, M. Rafelski, L. J. Prichard, S. Arnouts, T. Moutard, J. Coupon, A. Golob, S. Gwyn

Abstract: Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Comp… ▽ More Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$α)\leq50$Å and we find no correlation or anti-correlation between EW (Ly$α$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 22 pages, 12 figures, accepted for publication 30/3/2020

Efficient smoothed particle radiation hydrodynamics II: Radiation hydrodynamics

Authors: Brody R. Bassett, J. Michael Owen, Thomas A. Brunner

Abstract: The radiation hydrodynamics equations for smoothed particle hydrodynamics are derived by operator splitting the radiation and hydrodynamics terms, including necessary terms for material motion, and discretizing each of the sets of equations separately in time and space. The implicit radiative transfer discussed in the first paper of this series is coupled to explicit smoothed particle hydrodynamic… ▽ More The radiation hydrodynamics equations for smoothed particle hydrodynamics are derived by operator splitting the radiation and hydrodynamics terms, including necessary terms for material motion, and discretizing each of the sets of equations separately in time and space. The implicit radiative transfer discussed in the first paper of this series is coupled to explicit smoothed particle hydrodynamics. The result is a multi-material meshless radiation hydrodynamics code with arbitrary opacities and equations of state that performs well for problems with significant material motion. The code converges with second-order accuracy in space and first-order accuracy in time to the semianalytic solution for the Lowrie radiative shock problem and has competitive performance compared to a mesh-based radiation hydrodynamics code for a multi-material problem in two dimensions and an ablation problem inspired by inertial confinement fusion in two and three dimensions. △ Less

Submitted 16 November, 2020; v1 submitted 30 January, 2020; originally announced January 2020.

Report number: LLNL-JRNL-799714

Efficient smoothed particle radiation hydrodynamics I: Thermal radiative transfer

Authors: Brody R. Bassett, J. Michael Owen, Thomas A. Brunner

Abstract: This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled… ▽ More This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled material and radiation energy equations are solved using an inexact Newton iteration scheme based on nonlinear elimination, which reduces the number of Newton iterations needed to converge within each time step. During each Newton iteration, the radiation diffusion equation is solved using Krylov iterative methods with a multigrid preconditioner, which abstracts and optimizes much of the communication when running in parallel. The code is verified for an infinite medium problem, a one-dimensional Marshak wave, and a two and three-dimensional manufactured problem, and exhibits first-order convergence in time and second-order convergence in space. For these problems, the number of iterations needed to converge the inexact Newton scheme and the diffusion equation are independent of the number of spatial points and the number of processors. △ Less

Submitted 16 November, 2020; v1 submitted 30 January, 2020; originally announced January 2020.

Report number: LLNL-JRNL-799713

The intricate link between galaxy dynamics and intrinsic shape (or why so-called prolate rotation is a misnomer)

Authors: Caroline Foster, Robert Bassett

Abstract: Many recent integral integral field spectroscopy (IFS) survey teams have used stellar kinematic maps combined with imaging to statistically infer the underlying distributions of galaxy intrinsic shapes. With now several IFS samples at our disposal, the method, which was originally proposed by M. Franx and collaborators in 1991, is gaining in popularity, having been so far applied to ATLAS3D, SAMI,… ▽ More Many recent integral integral field spectroscopy (IFS) survey teams have used stellar kinematic maps combined with imaging to statistically infer the underlying distributions of galaxy intrinsic shapes. With now several IFS samples at our disposal, the method, which was originally proposed by M. Franx and collaborators in 1991, is gaining in popularity, having been so far applied to ATLAS3D, SAMI, MANGA and MASSIVE. We present results showing that a commonly assumed relationship between dynamical and intrinsic shape alignment does not hold in Illustris, affecting our ability to recover accurate intrinsic shape distributions. A further implication is that so-called "prolate rotation", where the bulk of stars in prolate galaxies are thought to rotate around the projected major axis, is a misnomer. △ Less

Submitted 22 August, 2019; originally announced August 2019.

Comments: 4 pages, 1 figure, IAU353 Proceedings

How many ways to color the map of America?

Authors: Rebekah Bassett, Jennifer Canizales, Jasbir S. Chahal, Thomas Fackrell, Vanessa Rico

Abstract: Although the Four Color Conjecture originated in cartography, surprisingly, there is nothing in the literature on the number of ways to color an actual geographic map with four or fewer colors. In this paper, we compute these numbers, with exponentially increasing order of difficulty, for the maps of Canada, France, and the USA. Our attempts to compute the latter two lead to some new results on th… ▽ More Although the Four Color Conjecture originated in cartography, surprisingly, there is nothing in the literature on the number of ways to color an actual geographic map with four or fewer colors. In this paper, we compute these numbers, with exponentially increasing order of difficulty, for the maps of Canada, France, and the USA. Our attempts to compute the latter two lead to some new results on the chromatic polynomial of graphs. △ Less

Submitted 15 August, 2019; originally announced August 2019.

Prospects for Recovering Galaxy Intrinsic Shapes from Projected Quantities

Authors: Robert Bassett, Caroline Foster

Abstract: The distribution of three dimensional intrinsic galaxy shapes has been a longstanding open question. The difficulty stems from projection effects meaning one must rely on statistical methods applied to galaxy samples to infer intrinsic shape distributions. Theoretical work using analytical galaxy potentials suggests a relationship between galaxy intrinsic shape (as defined by its "triaxiality", in… ▽ More The distribution of three dimensional intrinsic galaxy shapes has been a longstanding open question. The difficulty stems from projection effects meaning one must rely on statistical methods applied to galaxy samples to infer intrinsic shape distributions. Theoretical work using analytical galaxy potentials suggests a relationship between galaxy intrinsic shape (as defined by its "triaxiality", in practice a proxy for how prolate a galaxy is) and the intrinsic misalignment angle between kinematic and morphological axes ($Ψ_{\rm int}$). This relationship reduces the number of unknowns, providing more reliable inferred intrinsic shape distributions than methods using photometry alone. Here we explore the connection between intrinsic shape and stellar kinematics using cosmological hydrodynamical simulations from the Illustris project. The strongest relationship we find is that galaxy intrinsic flattening is correlated with specific angular momentum (j) with high j galaxies being flatter than galaxies with low specific angular momentum. Our analysis shows that, although the majority of kinematically misaligned galaxies exhibit prolate shapes, examples of kinematically aligned prolate galaxies are also present. Clearly a direct correspondence between prolate shape and minor-axis rotation (often referred to as "prolate rotation") is not present in Illustris. Thus, we demonstrate that the assumption of a simple relationship between $Ψ_{\rm int}$ and intrinsic shape commonly employed in shape recovery studies is not valid for Illustris galaxies. We suggest improvements on the method as well as some alternative methods for future work in this area. △ Less

Submitted 8 May, 2019; originally announced May 2019.

Comments: 18 pages, 9 figures, accepted for publication in MNRAS

On the Lack of Correlation Between [OIII]/[OII] and Lyman-Continuum Escape Fraction

Authors: R. Bassett, E. V. Ryan-Weber, J. Cooke, C. G. Diaz, T. Nanayakkara, T. -T. Yuan, L. R. Spitler, U. Meštrić, T. Garel, M. Sawicki, S. Gwyn, A. Golob

Abstract: We present the first results of our pilot study of 8 photometrically selected Lyman continuum (LyC) emitting galaxy candidates from the COSMOS field and focus on their optical emission line ratios. Observations were performed in the H and K bands using the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) instrument at the Keck Observatory, targeting the [OII], H$β$, and [OIII] emissio… ▽ More We present the first results of our pilot study of 8 photometrically selected Lyman continuum (LyC) emitting galaxy candidates from the COSMOS field and focus on their optical emission line ratios. Observations were performed in the H and K bands using the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) instrument at the Keck Observatory, targeting the [OII], H$β$, and [OIII] emission lines. We find that photometrically selected LyC emitting galaxy candidates have high ionization parameters, based on their high [OIII]/[OII] ratios (O32), with an average ratio for our sample of 2.5$\pm$0.2. Preliminary results of our companion Low Resolution Imaging Spectrometer (LRIS) observations, targeting LyC and Ly$α$, show that those galaxies with the largest O32 are typically found to also be Ly$α$ emitters. High O32 galaxies are also found to have tentative non-zero LyC escape fractions ($f_{esc}(LyC)$) based on $u$ band photometric detections. These results are consistent with samples of highly ionized galaxies, including confirmed LyC emitting galaxies from the literature. We also perform a detailed comparison between the observed emission line ratios and simulated line ratios from density bounded H$_{\textrm{II}}$ regions modeled using the photoionization code MAPPINGS V. Estimates of $f_{esc}(LyC)$ for our sample fall in the range from 0.0-0.23 and suggest possible tension with published correlations between O32 and $f_{esc}(LyC)$, adding weight to dichotomy of arguments in the literature. We highlight the possible effects of clumpy geometry and mergers that may account for such tension. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 21 pages, 11 figures, 3 tables, accepted for publication in MNRAS

Fused Density Estimation: Theory and Methods

Authors: Robert Bassett, James Sharpnack

Abstract: In this paper we introduce a method for nonparametric density estimation on geometric networks. We define fused density estimators as solutions to a total variation regularized maximum-likelihood density estimation problem. We provide theoretical support for fused density estimation by proving that the squared Hellinger rate of convergence for the estimator achieves the minimax bound over univaria… ▽ More In this paper we introduce a method for nonparametric density estimation on geometric networks. We define fused density estimators as solutions to a total variation regularized maximum-likelihood density estimation problem. We provide theoretical support for fused density estimation by proving that the squared Hellinger rate of convergence for the estimator achieves the minimax bound over univariate densities of log-bounded variation. We reduce the original variational formulation in order to transform it into a tractable, finite-dimensional quadratic program. Because random variables on geometric networks are simple generalizations of the univariate case, this method also provides a useful tool for univariate density estimation. Lastly, we apply this method and assess its performance on examples in the univariate and geometric network setting. We compare the performance of different optimization techniques to solve the problem, and use these results to inform recommendations for the computation of fused density estimators. △ Less

Submitted 4 December, 2018; v1 submitted 8 May, 2018; originally announced May 2018.

MSC Class: 62G07; 62G20; 90C35; 49M37

Robust Cross-correlation-based Measurement of Clump Sizes in Galaxies

Authors: Kamran Ali, Danail Obreschkow, David B. Fisher, Karl Glazebrook, Ivana Damjanov, Roberto G. Abraham, Robert Bassett

Abstract: Stars form in molecular complexes that are visible as giant clouds ($\sim 10^{5-6} \mathrm{M}_\odot$) in nearby galaxies and as giant clumps ($\sim 10^{8-9}\mathrm{M}_\odot$) in galaxies at redshifts $z\approx1$$-$$3$. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution… ▽ More Stars form in molecular complexes that are visible as giant clouds ($\sim 10^{5-6} \mathrm{M}_\odot$) in nearby galaxies and as giant clumps ($\sim 10^{8-9}\mathrm{M}_\odot$) in galaxies at redshifts $z\approx1$$-$$3$. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g.\ H$α$) can be recovered statistically using the two-point correlation function (2PCF), if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at $\sim20\%$ accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical H$α$-clump size measurements agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps. △ Less

Submitted 8 January, 2018; originally announced January 2018.

Comments: Published in ApJ on 2017 August 9. 13 pages, 8 figures

Journal ref: The Astrophysical Journal, Volume 845, Issue 1, article id. 37, 12 pp. (2017)

The connection between the peaks in velocity dispersion and star-forming clumps of turbulent galaxies

Authors: P. Oliva-Altamirano, D. Fisher, K. Glazebrook, E. Wisnioski, G. Bekiaris, R. Bassett, D. Obreschkow, R. Abraham

Abstract: We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable… ▽ More We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable PSF and spatial sampling, which we account for in our analysis. We present gas velocity dispersion maps corrected for these effects, and confirm that DYNAMO galaxies do have high gas velocity dispersion ($σ=40-80$\kms), even at high spatial sampling. We find statistically significant structure in 6 out of 7 galaxies. The most common distance between the peaks in velocity dispersion and emission line peaks is $\sim0.5$~kpc, we note this is very similar to the average size of a clump measured with HST H$α$ maps. This could suggest that the peaks in velocity dispersion in clumpy galaxies likely arise due to some interaction between the clump and the surrounding ISM of the galaxy, though our observations cannot distinguish between outflows, inflows or velocity shear. Observations covering a wider area of the galaxies will be needed to confirm this result. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in MNRAS

The Formation of S0 Galaxies with Counter-Rotating Neutral and Molecular Hydrogen

Authors: Robert Bassett, Kenji Bekki, Luca Cortese, Warrick J. Couch

Abstract: The observation of counter rotation in galaxies (i.e. gas that rotates in the opposite direction to the stellar component or two co-spatial stellar populations with opposite rotation) is becoming more commonplace with modern integral field spectroscopic surveys. In this paper we explore the emergence of counter-rotation (both stellar and gaseous) in S0 galaxies from smoothed-particle hydrodynamics… ▽ More The observation of counter rotation in galaxies (i.e. gas that rotates in the opposite direction to the stellar component or two co-spatial stellar populations with opposite rotation) is becoming more commonplace with modern integral field spectroscopic surveys. In this paper we explore the emergence of counter-rotation (both stellar and gaseous) in S0 galaxies from smoothed-particle hydrodynamics simulations of 1/10 mass ratio minor mergers between a 10^10.8 solar mass disk galaxy with a bulge-to-total ration of 0.17 and a gas rich companion (gas-to-stellar mass fraction of 5.0). These simulations include a self-consistent treatment of gas dynamics, star formation, the production/destruction of H2 and dust, and the time evolution of the interstellar radiation field. We explore the effect of retrograde versus prograde obits, gas and bulge mass fractions of the primary galaxy, and orbital parameters of the companion. The key requirement for producing counter rotation in stars or gas in a merger remnant is a retrograde primary, while the relative spin of the companion affects only the radial extent of the accreted gas. We also find that including a significant amount of gas in the primary can prevent the emergence of counter-rotating gas, although accreted stars retain counter-rotation. Bulge mass and orbit have a secondary effect, generally influencing the final distribution of accreted stars and gas within the framework outlined above. In addition to our primary focus of counter-rotating components in galaxies, we also make some predictions regarding the SFRs, H2 distributions, and dust in minor merger remnants. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 16 pages, 11 figues, 3 tables, accepted for publication in MNRAS

The SAMI Galaxy Survey: Kinematics of Dusty Early-Type Galaxies

Authors: Robert Bassett, K. Bekki, L. Cortese, W. J. Couch, A. E. Sansom, J. van de Sande, J. J. Bryant, C. Foster, S. M. Croom, S. Brough, S. M. Sweet, A. M. Medling, M. S. Owers, S. P. Driver, L. J. M. Davies, O. I. Wong, B. A. Groves, J. Bland-Hawthorn, S. N. Richards, M. Goodwin, I. S. Konstantopoulos, J. S. Lawrence

Abstract: Recently, large samples of visually classified early-type galaxies (ETGs) containing dust have been identified using space-based infrared observations with the Herschel Space Telescope. The presence of large quantities of dust in massive ETGs is peculiar as X-ray halos of these galaxies are expected to destroy dust in 10 Myr (or less). This has sparked a debate regarding the origin of the dust: is… ▽ More Recently, large samples of visually classified early-type galaxies (ETGs) containing dust have been identified using space-based infrared observations with the Herschel Space Telescope. The presence of large quantities of dust in massive ETGs is peculiar as X-ray halos of these galaxies are expected to destroy dust in 10 Myr (or less). This has sparked a debate regarding the origin of the dust: is it internally produced by asymptotic giant branch (AGB) stars, or is it accreted externally through mergers? We examine the 2D stellar and ionised gas kinematics of dusty ETGs using IFS observations from the SAMI galaxy survey, and integrated star-formation rates, stellar masses, and dust masses from the GAMA survey. Only 8% (4/49) of visually-classified ETGs are kinematically consistent with being dispersion-supported systems. These "dispersion-dominated galaxies" exhibit discrepancies between stellar and ionised gas kinematics, either offsets in the kinematic position angle or large differences in the rotational velocity, and are outliers in star-formation rate at a fixed dust mass compared to normal star-forming galaxies. These properties are suggestive of recent merger activity. The remaining 90% of dusty ETGs have low velocity dispersions and/or large circular velocities, typical of "rotation-dominated galaxies". These results, along with the general evidence of published works on X-ray emission in ETGs, suggest that they are unlikely to host hot, X-ray gas consistent with their low stellar mass when compared to dispersion-dominated galaxies. This means dust will be long lived and thus these galaxies do not require external scenarios for the origin of their dust content. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

The SAMI Galaxy Survey: mass as the driver of the kinematic morphology - density relation in clusters

Authors: Sarah Brough, Jesse van de Sande, Matt S. Owers, Francesco d'Eugenio, Rob Sharp, Luca Cortese, Nicholas Scott, Scott M. Croom, Rob Bassett, Kenji Bekki, Joss Bland-Hawthorn, Julia J. Bryant, Roger Davies, Michael J. Drinkwater, Simon P. Driver, Caroline Foster, Gregory Goldstein, A. R. Lopez-Sanchez, Anne M. Medling, Sarah M. Sweet, Dan S. Taranu, Chiara Tonini, Sukyoung K. Yi, Michael Goodwin, J. S. Lawrence , et al. (1 additional authors not shown)

Abstract: We examine the kinematic morphology of early-type galaxies (ETGs) in eight galaxy clusters in the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The clusters cover a mass range of 14.2<log(M_200/M_odot)<15.2 and we measure spatially-resolved stellar kinematics for 315 member galaxies with stellar masses 10.0<log(M_*/M_odot)<11.7 within 1R_200 of the cluster centers. We c… ▽ More We examine the kinematic morphology of early-type galaxies (ETGs) in eight galaxy clusters in the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The clusters cover a mass range of 14.2<log(M_200/M_odot)<15.2 and we measure spatially-resolved stellar kinematics for 315 member galaxies with stellar masses 10.0<log(M_*/M_odot)<11.7 within 1R_200 of the cluster centers. We calculate the spin parameter, lambda_R and use that to classify the kinematic morphology of the galaxies as fast or slow rotators. The total fraction of slow rotators in the early-type galaxy population, F_SR=0.14+/-0.02 and does not depend on host cluster mass. Across the eight clusters, the fraction of slow rotators increases with increasing local overdensity. We also find that the slow-rotator fraction increases at small clustercentric radii (R_cl<0.3R_200), and note that there is also an increase in slow-rotator fraction at R_cl~0.6R_200. The slow rotators at these larger radii reside in cluster substructure. We find the strongest increase in slow-rotator fraction occurs with increasing stellar mass. After accounting for the strong correlation with stellar mass, we find no significant relationship between spin parameter and local overdensity in the cluster environment. We conclude that the primary driver for the kinematic morphology--density relationship in galaxy clusters is the changing distribution of galaxy stellar mass with local environment. The presence of slow rotators in substructure suggests that the cluster kinematic morphology--density relationship is a result of mass segregation of slow-rotating galaxies forming in groups that later merge with clusters and sink to the cluster center via dynamical friction. △ Less

Submitted 15 June, 2017; v1 submitted 4 April, 2017; originally announced April 2017.

Comments: Accepted by ApJ 15 June 2017

Maximum a Posteriori Estimators as a Limit of Bayes Estimators

Authors: Robert Bassett, Julio Deride

Abstract: Maximum a posteriori and Bayes estimators are two common methods of point estimation in Bayesian Statistics. It is commonly accepted that maximum a posteriori estimators are a limiting case of Bayes estimators with 0-1 loss. In this paper, we provide a counterexample which shows that in general this claim is false. We then correct the claim that by providing a level-set condition for posterior den… ▽ More Maximum a posteriori and Bayes estimators are two common methods of point estimation in Bayesian Statistics. It is commonly accepted that maximum a posteriori estimators are a limiting case of Bayes estimators with 0-1 loss. In this paper, we provide a counterexample which shows that in general this claim is false. We then correct the claim that by providing a level-set condition for posterior densities such that the result holds. Since both estimators are defined in terms of optimization problems, the tools of variational analysis find a natural application to Bayesian point estimation. △ Less

Submitted 21 January, 2018; v1 submitted 17 November, 2016; originally announced November 2016.

MSC Class: 62C10; 62F10; 62F15; 65K10

Dust Attenuation in Clumpy, Star-Forming Galaxies at 0.07 < z < 0.14

Authors: Robert Bassett, Karl Glazebrook, David B. Fisher, Emily Wisnioski, Ivana Damjanov, Roberto Abraham, Danail Obreschkow, Andrew W. Green, Elisabete da Cunha, Peter J. McGregor

Abstract: Dust attenuation in galaxies has been extensively studied nearby, however, there are still many unknowns regarding attenuation in distant galaxies. We contribute to this effort using observations of star-forming galaxies in the redshift range z = 0.05-0.15 from the DYNAMO survey. Highly star-forming DYNAMO galaxies share many similar attributes to clumpy, star-forming galaxies at high redshift. Co… ▽ More Dust attenuation in galaxies has been extensively studied nearby, however, there are still many unknowns regarding attenuation in distant galaxies. We contribute to this effort using observations of star-forming galaxies in the redshift range z = 0.05-0.15 from the DYNAMO survey. Highly star-forming DYNAMO galaxies share many similar attributes to clumpy, star-forming galaxies at high redshift. Considering integrated Sloan Digital Sky Survey observations, trends between attenuation and other galaxy properties for DYNAMO galaxies are well matched to star-forming galaxies at high redshift. Integrated gas attenuations of DYNAMO galaxies are 0.2-2.0 mags in the V-band, and the ratio of stellar E(B-V) and gas E(B-V) is 0.78-0.08 (compared to 0.44 at low redshift). Four highly star-forming DYNAMO galaxies were observed at H-alpha using the Hubble Space Telescope and at Pa-alpha using integral field spectroscopy at Keck. The latter achieve similar resolution (~0.8-1 kpc) to our HST imaging using adaptive optics, providing resolved observations of gas attenuations of these galaxies on sub-kpc scales. We find < 1.0 mag of variation in attenuation (at H-alpha) from clump to clump, with no evidence of highly attenuated star formation. Attenuations are in the range 0.3-2.2 mags in the V band, consistent with attenuations of low redshift star-forming galaxies. The small spatial variation on attenuation suggests that a majority of the star-formation activity in these four galaxies occurs in relatively unobscured regions and, thus, star-formation is well characterised by our H-alpha observations. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: 16 pages, 16 figures, accepted for publication in MNRAS

DYNAMO-HST Survey: Clumps in Nearby Massive Turbulent Disks and the Effects of Clump Clustering on Kiloparsec Scale Measurements of Clumps

Authors: David B. Fisher, Karl Glazebrook, Ivana Damjanov, Roberto G. Abraham, Danail Obreschkow, Emily Wisnioski, Robert Bassett, Andy Green, Peter McGregor

Abstract: We present $\sim$100 pc resolution Hubble Space Telescope H$α$ images of 10 galaxies from the DYnamics of Newly-Assembled Massive Objects (DYNAMO) survey of low-$z$ turbulent disk galaxies, and use these to undertake the first detailed systematic study of the effects of resolution and clump clustering on observations of clumps in turbulent disks. In the DYNAMO-{\em HST} sample we measure clump dia… ▽ More We present $\sim$100 pc resolution Hubble Space Telescope H$α$ images of 10 galaxies from the DYnamics of Newly-Assembled Massive Objects (DYNAMO) survey of low-$z$ turbulent disk galaxies, and use these to undertake the first detailed systematic study of the effects of resolution and clump clustering on observations of clumps in turbulent disks. In the DYNAMO-{\em HST} sample we measure clump diameters spanning the range $d_{clump} \sim 100-800$~pc, and individual clump star formation rates as high as $\sim5$~M$_{\odot}$~yr$^{-1}$. DYNAMO clumps have very high SFR surface densities, $Σ_{SFR}\sim 15$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$, $\sim100\times$ higher than in H{\sc ii} regions of nearby spirals. Indeed, SFR surface density provides a simple dividing line between massive star forming clumps and local star forming regions, where massive star forming clumps have $Σ_{SFR}> 0.5$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$. When degraded to match the observations of galaxies in $z\sim 1-3$ surveys, DYNAMO galaxies are similar in morphology and measured clump properties to clumpy galaxies observed in the high-$z$ Universe. Emission peaks in the simulated high-redshift maps typically correspond to multiple clumps in full resolution images. This clustering of clumps systematically increases the apparent size and SFR of clumps in 1~kpc resolution maps, and decreases the measured SFR surface density of clumps by as much as a factor of 20$\times$. From these results we can infer that clump clustering is likely to strongly effect the measured properties of clumps in high-$z$ galaxies, which commonly have kiloparsec scale resolution. △ Less

Submitted 29 August, 2016; originally announced August 2016.

Comments: Accepted to MNRAS

A Human Computer Interaction Solution for Radiology Reporting: Evaluation of the Factors of Variation

Authors: Tejaswini Ganapathi, David Vining, Roland Bassett, Naveen Garg, Mia Markey

Abstract: The purpose of this research is to evaluate the human and technical factors required to create a human-computer interface (HCI) for a structured reporting solution based on eye-gaze and speech signals. Gaze and speech signals from radiologists acquired during simulated image interpretation and dictation sessions were analyzed to determine a) variation of temporal relationship between eye gaze and… ▽ More The purpose of this research is to evaluate the human and technical factors required to create a human-computer interface (HCI) for a structured reporting solution based on eye-gaze and speech signals. Gaze and speech signals from radiologists acquired during simulated image interpretation and dictation sessions were analyzed to determine a) variation of temporal relationship between eye gaze and speech in a dictation environment, and b) variation in eye movements for a particular image interpretation task among radiologists. Knowledge of these factors provides information regarding the complexity of the image interpretation or dictation task, and provides information that can be used to design a HCI for use in diagnostic radiology. Our ultimate goal is to use these data to create an HCI to automate the generation of a particular type of structured radiology report. Our data indicate that the a) temporal relationships between eye gaze and speech and b) scan paths substantially vary among radiologists, thus implying that an HCI system based on eye gaze and speech for automating the capture of data for structured reporting processes should be customized for each user. The image resolution and layout, image content, and order of targets during an image interpretation session are not relevant factors to consider when designing an HCI. Our findings can be applied to the design of other HCI solutions for radiological applications that involve visual inspection and verbal descriptions of image findings. △ Less

Submitted 3 August, 2016; v1 submitted 22 July, 2016; originally announced July 2016.

Comments: 17 pages, 7 figures, 1 table, In preparation for submission

Log-Concave Duality in Estimation and Control

Authors: Robert Bassett, Michael Casey, Roger J-B Wets

Abstract: In this paper we generalize the estimation-control duality that exists in the linear-quadratic-Gaussian setting. We extend this duality to maximum a posteriori estimation of the system's state, where the measurement and dynamical system noise are independent log-concave random variables. More generally, we show that a problem which induces a convex penalty on noise terms will have a dual control p… ▽ More In this paper we generalize the estimation-control duality that exists in the linear-quadratic-Gaussian setting. We extend this duality to maximum a posteriori estimation of the system's state, where the measurement and dynamical system noise are independent log-concave random variables. More generally, we show that a problem which induces a convex penalty on noise terms will have a dual control problem. We provide conditions for strong duality to hold, and then prove relaxed conditions for the piecewise linear-quadratic case. The results have applications in estimation problems with nonsmooth densities, such as log-concave maximum likelihood densities. We conclude with an example reconstructing optimal estimates from solutions to the dual control problem, which has implications for sharing solution methods between the two types of problems. △ Less

Submitted 8 July, 2016; originally announced July 2016.

The Spitzer-HETDEX Exploratory Large-Area Survey

Authors: Casey Papovich, H. V. Shipley, N. Mehrtens, C. Lanham, M. Lacy, R. Ciardullo, S. L. Finkelstein, R. Bassett, P Behroozi, G. A. Blanc, R. S. de Jong, D. L. DePoy, N. Drory, E. Gawiser, K. Gebhardt, C. Gronwall, G. J. Hill, U. Hopp, S. Jogee, L. Kawinwanichakij, J. L. Marshall, E. McLinden, E. Mentuch Cooper, R. S. Somerville, M. Steinmetz , et al. (5 additional authors not shown)

Abstract: We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectr… ▽ More We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$α$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $μ$Jy at both 3.6 and 4.5 micron (1$σ$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication. △ Less

Submitted 17 March, 2016; originally announced March 2016.

Comments: Accepted for publication in ApJS. 30 pages, emulateapj format. Many figures. All data products, including images and catalogs available http://irsa.ipac.caltech.edu/data/SPITZER/SHELA

Multistage Portfolio Optimization: A Duality Result in Conic Market Models

Authors: Robert Bassett, Khoa Le

Abstract: We prove a general duality result for multi-stage portfolio optimization problems in markets with proportional transaction costs. The financial market is described by Kabanov's model of foreign exchange markets over a finite probability space and finite-horizon discrete time steps. This framework allows us to compare vector-valued portfolios under a partial ordering, so that our model does not req… ▽ More We prove a general duality result for multi-stage portfolio optimization problems in markets with proportional transaction costs. The financial market is described by Kabanov's model of foreign exchange markets over a finite probability space and finite-horizon discrete time steps. This framework allows us to compare vector-valued portfolios under a partial ordering, so that our model does not require liquidation into some numeraire at terminal time. We embed the vector-valued portfolio problem into the set-optimization framework, and generate a problem dual to portfolio optimization. Using recent results in the development of set optimization, we then show that a strong duality relationship holds between the problems. △ Less

Submitted 22 January, 2016; v1 submitted 4 January, 2016; originally announced January 2016.

Low Angular Momentum in Clumpy, Turbulent Disk Galaxies

Authors: Danail Obreschkow, Karl Glazebrook, Robert Bassett, David B. Fisher, Roberto G. Abraham, Emily Wisnioski, Andrew W. Green, Peter J. McGregor, Ivana Damjanov, Attila Popping, Inger Jorgensen

Abstract: We measure the stellar specific angular momentum jstar=Jstar/Mstar in four nearby (z~0.1) disk galaxies that have stellar masses Mstar near the break M* of the galaxy mass function, but look like typical star-forming disks at z~2 in terms of their low stability (Q~1), clumpiness, high ionized gas dispersion (40-50 km/s), high molecular gas fraction (20-30%) and rapid star formation (~20 Msun/yr).… ▽ More We measure the stellar specific angular momentum jstar=Jstar/Mstar in four nearby (z~0.1) disk galaxies that have stellar masses Mstar near the break M* of the galaxy mass function, but look like typical star-forming disks at z~2 in terms of their low stability (Q~1), clumpiness, high ionized gas dispersion (40-50 km/s), high molecular gas fraction (20-30%) and rapid star formation (~20 Msun/yr). Combining high-resolution (Keck-OSIRIS) and large-radius (Gemini-GMOS) spectroscopic maps, only available at low z, we discover that these targets have about three times less stellar angular momentum than typical local spiral galaxies of equal stellar mass and bulge fraction. Theoretical considerations show that this deficiency in angular momentum is the main cause of their low stability, while the high gas fraction plays a complementary role. Interestingly, the low jstar values of our targets are similar to those expected in the M*-population at higher z from the approximate theoretical scaling jstar~(1+z)^(-1/2) at fixed Mstar. This suggests that a change in angular momentum, driven by cosmic expansion, is the main cause for the remarkable difference between clumpy M*-disks at high z (which likely evolve into early-type galaxies) and mass-matched local spirals. △ Less

Submitted 16 November, 2015; v1 submitted 19 August, 2015; originally announced August 2015.

Comments: 4 Figures (including one interactive 3D figure), 1 Table

The host galaxies of X-ray selected Active Galactic Nuclei to z=2.5: Structure, star-formation and their relationships from CANDELS and Herschel/PACS

Authors: D. J. Rosario, D. H. McIntosh, A. van der Wel, J. Kartaltepe, P. Lang, P. Santini, S. Wuyts, D. Lutz, M. Rafelski, C. Villforth, D. M. Alexander, F. E. Bauer, E. F. Bell, S. Berta, W. N. Brandt, C. J. Conselice, A. Dekel, S. M. Faber, H. C. Ferguson, R. Genzel, N. A. Grogin, D. D. Kocevski, A. M. Koekemoer, D. C. Koo, J. M. Lotz , et al. (38 additional authors not shown)

Abstract: We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of s… ▽ More We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5. △ Less

Submitted 17 September, 2014; originally announced September 2014.

Comments: 25 pages (double-column), 17 Figures, under review at A&A

Journal ref: A&A 573, A85 (2015)

Extreme gas fractions in clumpy, turbulent disk galaxies at z~0.1

Authors: David B. Fisher, Karl Glazebrook, Alberto Bolatto, Danail Obreschkow, Erin Mentuch-Cooper, Emily Wisnioski, Robert BAssett, Roberto G. Abraham, Ivana Damjanov, Andy Green, Peter McGregor

Abstract: In this letter we report the discovery of CO fluxes, suggesting very high gas fractions in three disk galaxies seen in the nearby Universe (z ~ 0.1). These galaxies were investigated as part of the DYnamics of Newly Assembled Massive Objects (DYNAMO) survey. High-resolution Hubble Space Telescope imaging of these objects reveals the presence of large star forming clumps in the bodies of the galaxi… ▽ More In this letter we report the discovery of CO fluxes, suggesting very high gas fractions in three disk galaxies seen in the nearby Universe (z ~ 0.1). These galaxies were investigated as part of the DYnamics of Newly Assembled Massive Objects (DYNAMO) survey. High-resolution Hubble Space Telescope imaging of these objects reveals the presence of large star forming clumps in the bodies of the galaxies, while spatially resolved spectroscopy of redshifted Halpha reveals the presence of high dispersion rotating disks. The internal dynamical state of these galaxies resembles that of disk systems seen at much higher redshifts (1 < z < 3). Using CO(1-0) observations made with the Plateau de Bure Interferometer, we find gas fractions of 20-30% and depletion times of tdep ~ 0.5 Gyr (assuming a Milky Way-like CO conversion factor). These properties are unlike those expected for low- redshift galaxies of comparable specific star formation rate, but they are normal for their high-z counterparts. DYNAMO galaxies break the degeneracy between gas fraction and redshift, and we show that the depletion time per specific star formation rate for galaxies is closely tied to gas fraction, independent of redshift. We also show that the gas dynamics of two of our local targets corresponds to those expected from unstable disks, again resembling the dynamics of high-z disks. These results provide evidence that DYNAMO galaxies are local analogues to the clumpy, turbulent disks, which are often found at high redshift. △ Less

Submitted 28 May, 2014; originally announced May 2014.

Comments: Accepted to ApJ Letters

DYNAMO II: Coupled Stellar and Ionized Gas Kinematics in Two Low Redshift Clumpy Disks

Authors: Robert Bassett, Karl Glazebrook, David B. Fisher, Andrew W. Green, Emily Wisnioski, Danail Obreschkow, Erin Mentuch Cooper, Roberto G. Abraham, Ivana Damjanov, Peter J. McGregor

Abstract: We study the spatially resolved stellar kinematics of two star-forming galaxies at z = 0.1 from the larger DYnamics of Newly Assembled Massive Objects (DYNAMO) sample. These galaxies, which have been characterized by high levels of star formation and large ionized gas velocity dispersions, are considered possible analogs to high-redshift clumpy disks. They were observed using the GMOS instrument i… ▽ More We study the spatially resolved stellar kinematics of two star-forming galaxies at z = 0.1 from the larger DYnamics of Newly Assembled Massive Objects (DYNAMO) sample. These galaxies, which have been characterized by high levels of star formation and large ionized gas velocity dispersions, are considered possible analogs to high-redshift clumpy disks. They were observed using the GMOS instrument in integral field spectroscopy (IFS) mode at the Gemini Observatory with high spectral resolution (R=5400, equivalent to 24 km/s at the observed wavelengths) and 6 hour exposure times in order to measure the resolved stellar kinematics via absorption lines. We also obtain higher-quality emission line kinematics than previous observations. The spatial resolution (1.2 kpc) is sufficient to show that the ionized gas in these galaxies (as traced by H-beta emission) is morphologically irregular, forming multiple giant clumps while stellar continuum light is smooth and well described by an exponential profile. Clumpy gas morphologies observed in IFS data are confirmed by complementary narrow band H-alpha imaging from the Hubble Space Telescope. Morphological differences between the stars and ionized gas are not reflected dynamically as stellar kinematics are found the be closely coupled to the kinematics of the ionized gas: both components are smoothly rotating with large velocity dispersions (~40 km/s) suggesting that the high gas dispersions are not primarily driven by star-formation feedback. In addition, the stellar population ages of these galaxies are estimated to be quite young (60-500 Myr). The large velocity dispersions measured for these young stars suggest that we are seeing the formation of thick disks and/or stellar bulges in support of recent models which produce these from clumpy galaxies at high redshift. △ Less

Submitted 26 May, 2014; originally announced May 2014.

Comments: 18 pages, 11 figures, 4 tables. Accepted for publication, MNRAS

A Quantitative Doignon-Bell-Scarf Theorem

Authors: Iskander Aliev, Robert Bassett, Jesus A. De Loera, Quentin Louveaux

Abstract: The famous Doignon-Bell-Scarf Theorem is a Helly-type result about the existence of integer solutions on systems of linear inequalities. The purpose of this paper is to present the following quantitative generalization: Given an integer $k$, we prove that there exists a constant $c(n,k)$, depending only on the dimension $n$ and $k$, such that if a polyhedron ${x: Ax \leq b}$ contains exactly k int… ▽ More The famous Doignon-Bell-Scarf Theorem is a Helly-type result about the existence of integer solutions on systems of linear inequalities. The purpose of this paper is to present the following quantitative generalization: Given an integer $k$, we prove that there exists a constant $c(n,k)$, depending only on the dimension $n$ and $k$, such that if a polyhedron ${x: Ax \leq b}$ contains exactly k integer solutions, then there exists a subset of the rows, of cardinality no more than $c(n,k)$, defining a polyhedron that contains exactly the same $k$ integer points. In this case $c(n,0) = 2^n$ is the original case of Doignon-Bell-Scarf for infeasible systems of inequalities. We work on both upper and lower bounds for the constant $c(n,k)$ and discuss some consequences, including a Clarkson-style algorithm to find the $l$-th best solution of an integer program with respect to the ordering induced by the objective function. △ Less

Submitted 7 September, 2015; v1 submitted 10 May, 2014; originally announced May 2014.

CANDELS Visual Classifications: Scheme, Data Release, and First Results

Authors: Jeyhan S. Kartaltepe, Mark Mozena, Dale Kocevski, Daniel H. McIntosh, Jennifer Lotz, Eric F. Bell, Sandra Faber, Henry Ferguson, David Koo, Robert Bassett, Maksym Bernyk, Kirsten Blancato, Frederic Bournaud, Paolo Cassata, Marco Castellano, Edmond Cheung, Christopher J. Conselice, Darren Croton, Tomas Dahlen, Duilia F. de Mello, Laura DeGroot, Jennifer Donley, Javiera Guedes, Norman Grogin, Nimish Hathi , et al. (27 additional authors not shown)

Abstract: We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wid… ▽ More We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band. △ Less

Submitted 9 October, 2015; v1 submitted 10 January, 2014; originally announced January 2014.

Comments: 15 pages, 14 figures, Accepted for publication in The Astrophysical Journal Supplement Series For access to data tables prior to publication contact jeyhan@astro.rit.edu

CANDELS Observations of the Environmental Dependence of the Color-Mass-Morphology Relation at z = 1.6

Authors: Robert Bassett, Casey Papovich, Jennifer M. Lotz, Eric F. Bell, Steven L. Finkelstein, Jeffrey A. Newman, Kim-Vy Tran, Omar Almaini, Caterina Lani, Michael Cooper, Darren Croton, Avishai Dekel, Henry C. Ferguson, Dale D. Kocevski, Anton M. Koekemoer, David C. Koo, Elizabeth J. McGrath, Daniel H. McIntosh, Risa H. Wechsler

Abstract: We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and Sérsic index, n. In both the cluster and field, approximately half of the bulge-dominated… ▽ More We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and Sérsic index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (n > 2) reside on the red sequence of the color-magnitude diagram, and most disk-dominated galaxies (n < 2) have colors expected for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U - B colors and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster have larger median effective radii and smaller Sérsic indices compared to the field with a significance of 2?. These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < Rproj < 1.5 Mpc, which have low Sérsic indices and possibly larger effective radii, more consistent with star-forming galaxies at this epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are consistent with models where gas accretion onto these galaxies is suppressed from processes associated with the cluster environment. △ Less

Submitted 2 May, 2013; originally announced May 2013.

Comments: ApJ accepted, 19 pages, 10 figures

CANDELS Observations of the Structural Properties and Evolution of Galaxies in a Cluster at z=1.62

Authors: Casey Papovich, R. Bassett, J. M. Lotz, A. van der Wel, K. -V. Tran, S. L. Finkelstein, E. F. Bell, C. J. Conselice, A. Dekel, J. S. Dunlop, Yicheng Guo, S. M. Faber, D. Farrah, H. C. Ferguson, K. D. Finkelstein, B. Häußler, D. D. Kocevski, A. Koekemoer, D. C. Koo, R. J. Mclure, E. J. McGrath, D. H. McIntosh, I. Momcheva, J. A. Newman, G. Rudnick , et al. (3 additional authors not shown)

Abstract: We discuss the structural and morphological properties of galaxies in a z=1.62 proto-cluster using near-IR imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: galaxies with colors of quiescent stellar populations generally have morphologies consis… ▽ More We discuss the structural and morphological properties of galaxies in a z=1.62 proto-cluster using near-IR imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: galaxies with colors of quiescent stellar populations generally have morphologies consistent with spheroids, and galaxies with colors consistent with ongoing star formation have disk-like and irregular morphologies. The size distribution of the quiescent cluster galaxies shows a deficit of compact (< 1kpc), massive galaxies compared to CANDELS field galaxies at z=1.6. As a result the cluster quiescent galaxies have larger average effective sizes compared to field galaxies at fixed mass at greater than 90% significance. Combined with data from the literature, the size evolution of quiescent cluster galaxies is relatively slow from z~1.6 to the present, growing as (1+z)^(-0.6+/-0.1). If this result is generalizable, then it implies that physical processes associated with the denser cluster region seems to have caused accelerated size growth in quiescent galaxies prior to z=1.6 and slower subsequent growth at z<1.6 compared to galaxies in the lower density field. The quiescent cluster galaxies at z=1.6 have higher ellipticities compared to lower redshift samples at fixed mass, and their surface-brightness profiles suggest that they contain extended stellar disks. We argue the cluster galaxies require dissipationless (i.e., gas-poor or "dry") mergers to reorganize the disk material and to match the relations for ellipticity, stellar mass, size, and color of early-type galaxies in z<1 clusters. △ Less

Submitted 29 February, 2012; v1 submitted 17 October, 2011; originally announced October 2011.

Comments: Accepted for publication in ApJ. 14 pages in emulateapj format. Replacement includes improvements from referee report, and updates and additions to references

CANDELS: The Evolution of Galaxy Rest-Frame Ultraviolet Colors from z = 8 to 4

Authors: Steven L. Finkelstein, Casey Papovich, Brett Salmon, Kristian Finlator, Mark Dickinson, Henry C. Ferguson, Mauro Giavalisco, Anton M. Koekemoer, Naveen A. Reddy, Robert Bassett, Christopher J. Conselice, James S. Dunlop, S. M. Faber, Norman A. Grogin, Nimish P. Hathi, Dale D. Kocevski, Kamson Lai, Kyoung-Soo Lee, Ross J. McLure, Bahram Mobasher, Jeffrey A. Newman

Abstract: We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synth… ▽ More We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8. △ Less

Submitted 8 August, 2012; v1 submitted 17 October, 2011; originally announced October 2011.

Comments: Accepted to the Astrophysical Journal. 19 pages, 8 figures and 5 tables. The complete data for Table 3 will be available in the online version of the Astrophysical Journal, and a version is included in the source files for this submission