Showing 1–46 of 46 results for author: Goldberg, M

The Prevalence of Resonance Among Young, Close-in Planets

Authors: Fei Dai, Max Goldberg, Konstantin Batygin, Jennifer van Saders, Eugene Chiang, Nick Choksi, Rixin Li, Erik A. Petigura, Gregory J. Gilbert, Sarah C. Millholland, Yuan-Zhe Dai, Luke Bouma, Lauren M. Weiss, Joshua N. Winn

Abstract: Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common… ▽ More Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common in younger systems. This prediction can now be tested, thanks to recent discoveries of young planets, particularly those in stellar clusters, by NASA's {\it TESS} mission. We divided the known planetary systems into three age groups: young ($<$100-Myr-old), adolescent (0.1-1-Gyr-old), and mature ($>1$-Gyr-old). The fraction of neighboring planet pairs having period ratios within a few percent of a first-order commensurability (e.g.~4:3, 3:2, or 2:1) is 70$\pm$15\% for young pairs, 24$\pm$8\% for adolescent pairs, and 15$\pm$2\% for mature pairs. The fraction of systems with at least one nearly commensurable pair (either first or second-order) is 86$\pm13$\% among young systems, 38$\pm12$\% for adolescent systems, and 23$\pm3$\% for mature systems. First-order commensurabilities prevail across all age groups, with an admixture of second-order commensurabilities. Commensurabilities are more common in systems with high planet multiplicity and low mutual inclinations. Observed period ratios often deviate from perfect commensurability by $\sim$1\% even among young planets, too large to be explained by resonant repulsion with equilibrium eccentricity tides. We also find that super-Earths in the radius gap ($1.5-1.9R_\oplus$) are less likely to be near-resonant (11.9$\pm2.0\%$) compared to Earth-sized planets ($R_p<1R_\oplus$; 25.3$\pm4.4\%$) or mini-Neptunes ($1.9R_\oplus \leq R_p<2.5R_\oplus$; 14.4$\pm1.8\%$). △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures, submitted to AAS

Breaking Giant Chains: Early-Stage Instabilities in Long-Period Giant Planet Systems

Authors: Vighnesh Nagpal, Max Goldberg, Konstantin Batygin

Abstract: Orbital evolution is a critical process that sculpts planetary systems, particularly during their early stages where planet-disk interactions are expected to lead to the formation of resonant chains. Despite the theoretically expected prominence of such configurations, they are scarcely observed among long-period giant exoplanets. This disparity suggests an evolutionary sequence wherein giant plan… ▽ More Orbital evolution is a critical process that sculpts planetary systems, particularly during their early stages where planet-disk interactions are expected to lead to the formation of resonant chains. Despite the theoretically expected prominence of such configurations, they are scarcely observed among long-period giant exoplanets. This disparity suggests an evolutionary sequence wherein giant planet systems originate in compact multiresonant configurations, but subsequently become unstable, eventually relaxing to wider orbits--a phenomenon mirrored in our own solar system's early history. In this work, we present a suite of N-body simulations that model the instability-driven evolution of giant planet systems, originating from resonant initial conditions, through phases of disk dispersal and beyond. By comparing the period ratio and normalized angular momentum distributions of our synthetic aggregate of systems with the observational census of long-period Jovian planets, we derive constraints on the expected rate of orbital migration, the efficiency of gas-driven eccentricity damping, and typical initial multiplicity. Our findings reveal a distinct inclination toward densely packed initial conditions, weak damping, and high giant planet multiplicities. Furthermore, our models indicate that resonant chain origins do not facilitate the formation of Hot Jupiters via the coplanar high-eccentricity pathway at rates high enough to explain their observed prevalence. △ Less

Submitted 3 August, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

Comments: ApJ published version. Changes to Hot Jupiter production calculation in Section 5; main conclusions remain unchanged

Journal ref: ApJ, 969, 2, 133 (2024)

Chaotic tides as a solution to the Hyperion problem

Authors: Max Goldberg, Konstantin Batygin

Abstract: The dynamics of the outer regular satellites of Saturn are driven primarily by the outward migration of Titan, but several independent constraints on Titan's migration are difficult to reconcile with the current resonant orbit of the small satellite Hyperion. We argue that Hyperion's rapid irregular tumbling greatly increases tidal dissipation with a steep dependence on orbital eccentricity. Reson… ▽ More The dynamics of the outer regular satellites of Saturn are driven primarily by the outward migration of Titan, but several independent constraints on Titan's migration are difficult to reconcile with the current resonant orbit of the small satellite Hyperion. We argue that Hyperion's rapid irregular tumbling greatly increases tidal dissipation with a steep dependence on orbital eccentricity. Resonant excitation from a migrating Titan is then balanced by damping in a feedback mechanism that maintains Hyperion's eccentricity without fine-tuning. The inferred tidal parameters of Hyperion are most consistent with rapid Titan migration enabled by a resonance lock with an internal mode of Saturn, but a scenario with only equilibrium dissipation in Saturn is also possible. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 11 pages, 4 figures, accepted to Icarus

A $5M_\text{Jup}$ Non-Transiting Coplanar Circumbinary Planet Around Kepler-1660AB

Authors: Max Goldberg, Daniel Fabrycky, David V. Martin, Simon Albrecht, Hans J. Deeg, Grzegorz Nowak

Abstract: Over a dozen transiting circumbinary planets have been discovered around eclipsing binaries. Transit detections are biased towards aligned planet and binary orbits, and indeed all of the known planets have mutual inclinations less than $4.5^{\circ}$. One path to discovering circumbinary planets with misaligned orbits is through eclipse timing variations (ETVs) of non-transiting planets. Borkovits… ▽ More Over a dozen transiting circumbinary planets have been discovered around eclipsing binaries. Transit detections are biased towards aligned planet and binary orbits, and indeed all of the known planets have mutual inclinations less than $4.5^{\circ}$. One path to discovering circumbinary planets with misaligned orbits is through eclipse timing variations (ETVs) of non-transiting planets. Borkovits et al. (2016) discovered ETVs on the 18.6 d binary Kepler-1660AB, indicative of a third body on a $\approx 236$ d period, with a misaligned orbit and a potentially planetary mass. Getley et al. (2017) agreed with the planetary hypothesis, arguing for a $7.7M_{\rm Jup}$ circumbinary planet on an orbit that is highly misaligned by $120^{\circ}$ with respect to the binary. In this paper, we obtain the first radial velocities of the binary. We combine these with an analysis of not only the ETVs but also the eclipse depth variations. We confirm the existence of a $239.5$ d circumbinary planet, but with a lower mass of $4.87M_{\rm Jup}$ and a coplanar orbit. The misaligned orbits proposed by previous authors are definitively ruled out by a lack of eclipse depth variations. Kepler-1660ABb is the first confirmed circumbinary planet found using ETVs around a main sequence binary. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: Resubmitted to MNRAS following positive referee report

Dynamics and Origins of the Near-Resonant Kepler Planets

Authors: Max Goldberg, Konstantin Batygin

Abstract: Short-period super-Earths and mini-Neptunes encircle more than $\sim50\%$ of Sun-like stars and are relatively amenable to direct observational characterization. Despite this, environments in which these planets accrete are difficult to probe directly. Nevertheless, pairs of planets that are close to orbital resonances provide a unique window into the inner regions of protoplanetary disks, as they… ▽ More Short-period super-Earths and mini-Neptunes encircle more than $\sim50\%$ of Sun-like stars and are relatively amenable to direct observational characterization. Despite this, environments in which these planets accrete are difficult to probe directly. Nevertheless, pairs of planets that are close to orbital resonances provide a unique window into the inner regions of protoplanetary disks, as they preserve the conditions of their formation, as well as the early evolution of their orbital architectures. In this work, we present a novel approach toward quantifying transit timing variations within multi-planetary systems and examine the near-resonant dynamics of over 100 planet pairs detected by \textit{Kepler}. Using an integrable model for first-order resonances, we find a clear transition from libration to circulation of the resonant angle at a period ratio of $\approx 0.6\%$ wide of exact resonance. The orbital properties of these systems indicate that they systematically lie far away from the resonant forced equilibrium. Cumulatively our modeling indicates that while orbital architectures shaped by strong disk damping or tidal dissipation are inconsistent with observations, a scenario where stochastic stirring by turbulent eddies augments the dissipative effects of protoplanetary disks reproduces several features of the data. △ Less

Submitted 3 May, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted to ApJ

Journal ref: The Astrophysical Journal, Vol. 948 (2023), 12

TOI-1136 is a Young, Coplanar, Aligned Planetary System in a Pristine Resonant Chain

Authors: Fei Dai, Kento Masuda, Corey Beard, Paul Robertson, Max Goldberg, Konstantin Batygin, Luke Bouma, Jack J. Lissauer, Emil Knudstrup, Simon Albrecht, Andrew W. Howard, Heather A. Knutson, Erik A. Petigura, Lauren M. Weiss, Howard Isaacson, Martti Holst Kristiansen, Hugh Osborn, Songhu Wang, Xian-Yu Wang, Aida Behmard, Michael Greklek-McKeon, Shreyas Vissapragada, Natalie M. Batalha, Casey L. Brinkman, Ashley Chontos , et al. (38 additional authors not shown)

Abstract: Convergent disk migration has long been suspected to be responsible for forming planetary systems with a chain of mean-motion resonances (MMR). Dynamical evolution over time could disrupt the delicate resonant configuration. We present TOI-1136, a 700-Myr-old G star hosting at least 6 transiting planets between $\sim$2 and 5 $R_\oplus$. The orbital period ratios deviate from exact commensurability… ▽ More Convergent disk migration has long been suspected to be responsible for forming planetary systems with a chain of mean-motion resonances (MMR). Dynamical evolution over time could disrupt the delicate resonant configuration. We present TOI-1136, a 700-Myr-old G star hosting at least 6 transiting planets between $\sim$2 and 5 $R_\oplus$. The orbital period ratios deviate from exact commensurability by only $10^{-4}$, smaller than the $\sim$\,$10^{-2}$ deviations seen in typical Kepler near-resonant systems. A transit-timing analysis measured the masses of the planets (3-8$M_\oplus$) and demonstrated that the planets in TOI-1136 are in true resonances with librating resonant angles. Based on a Rossiter-McLaughlin measurement of planet d, the star's rotation appears to be aligned with the planetary orbital planes. The well-aligned planetary system and the lack of detected binary companion together suggest that TOI-1136's resonant chain formed in an isolated, quiescent disk with no stellar fly-by, disk warp, or significant axial asymmetry. With period ratios near 3:2, 2:1, 3:2, 7:5, and 3:2, TOI-1136 is the first known resonant chain involving a second-order MMR (7:5) between two first-order MMR. The formation of the delicate 7:5 resonance places strong constraints on the system's migration history. Short-scale (starting from $\sim$0.1 AU) Type-I migration with an inner disk edge is most consistent with the formation of TOI-1136. A low disk surface density ($Σ_{\rm 1AU}\lesssim10^3$g~cm$^{-2}$; lower than the minimum-mass solar nebula) and the resultant slower migration rate likely facilitated the formation of the 7:5 second-order MMR. TOI-1136's deep resonance suggests that it has not undergone much resonant repulsion during its 700-Myr lifetime. One can rule out rapid tidal dissipation within a rocky planet b or obliquity tides within the largest planets d and f. △ Less

Submitted 14 November, 2022; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 48 pages, 23 figures, 8 tables. Accepted to AAS journals. Comments welcome!

A Criterion for the Stability of Planets in Chains of Resonances

Authors: Max Goldberg, Konstantin Batygin, Alessandro Morbidelli

Abstract: Uncovering the formation process that reproduces the distinct properties of compact super-Earth exoplanet systems is a major goal of planet formation theory. The most successful model argues that non-resonant systems begin as resonant chains of planets that later experience a dynamical instability. However, both the boundary of stability in resonant chains and the mechanism of the instability itse… ▽ More Uncovering the formation process that reproduces the distinct properties of compact super-Earth exoplanet systems is a major goal of planet formation theory. The most successful model argues that non-resonant systems begin as resonant chains of planets that later experience a dynamical instability. However, both the boundary of stability in resonant chains and the mechanism of the instability itself are poorly understood. Previous work postulated that a secondary resonance between the fastest libration frequency and a difference in synodic frequencies destabilizes the system. Here, we use that hypothesis to produce a simple and general criterion for resonant chain stability that depends only on planet orbital periods and masses. We show that the criterion accurately predicts the maximum mass of planets in synthetic resonant chains up to six planets. More complicated resonant chains produced in population synthesis simulations are found to be less stable than expected, although our criterion remains useful and superior to machine learning models. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 8 pages, 4 figures; accepted to Icarus

Cosmic flexion

Authors: Evan J. Arena, David M. Goldberg, David J. Bacon

Abstract: Cosmic flexion, like cosmic shear, is a correlation function whose signal originates from the large-scale structure of the universe. Building on the observational success of cosmic shear, along with the unprecedented quality of large-scale cosmological datasets, the time is ripe to explore the practical constraints from cosmic flexion. Unlike cosmic shear, which has a broad window function for pow… ▽ More Cosmic flexion, like cosmic shear, is a correlation function whose signal originates from the large-scale structure of the universe. Building on the observational success of cosmic shear, along with the unprecedented quality of large-scale cosmological datasets, the time is ripe to explore the practical constraints from cosmic flexion. Unlike cosmic shear, which has a broad window function for power, cosmic flexion is only measurable on small scales and therefore can uniquely place constraints on the small-scale matter power spectrum. Here, we present a full theoretical formalism for cosmic flexion, including both flexion-flexion and shear-flexion two-point correlations. We present forecasts for measuring cosmic flexion in the Dark Energy Survey (DES), a Stage III cosmological survey, and comment on the future prospects of measuring these cosmological flexion signals in the upcoming era of Stage IV experiments. △ Less

Submitted 29 June, 2022; v1 submitted 22 March, 2022; originally announced March 2022.

Comments: 21 pages, 9 figures; accepted for publication in Phys.Rev.D; revisions to Section IV. A and Appendix A, as well as other minor changes from submitted version

Report number: Phys. Rev. D 105, 123521 (2022)

Architectures of Compact Super-Earth Systems Shaped by Instabilities

Authors: Max Goldberg, Konstantin Batygin

Abstract: Compact non-resonant systems of sub-Jovian planets are the most common outcome of the planet formation process. Despite exhibiting broad overall diversity, these planets also display dramatic signatures of intra-system uniformity in their masses, radii, and orbital spacings. Although the details of their formation and early evolution are poorly known, sub-Jovian planets are expected to emerge from… ▽ More Compact non-resonant systems of sub-Jovian planets are the most common outcome of the planet formation process. Despite exhibiting broad overall diversity, these planets also display dramatic signatures of intra-system uniformity in their masses, radii, and orbital spacings. Although the details of their formation and early evolution are poorly known, sub-Jovian planets are expected to emerge from their natal nebulae as multi-resonant chains, owing to planet-disk interactions. Within the context of this scenario, the architectures of observed exoplanet systems can be broadly replicated if resonances are disrupted through post-nebular dynamical instabilities. Here, we generate an ad-hoc sample of resonant chains and use a suite of N-body simulations to show that instabilities can not only reproduce the observed period ratio distribution, but that the resulting collisions also modify the mass uniformity in a way that is consistent with the data. Furthermore, we demonstrate that primordial mass uniformity, motivated by the sample of resonant chains coupled with dynamical sculpting, naturally generates uniformity in orbital period spacing similar to what is observed. Finally, we find that almost all collisions lead to perfect mergers, but some form of post-instability damping is likely needed to fully account for the present-day dynamically cold architectures of sub-Jovian exoplanets. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 14 pages, 8 figures, accepted for publication in AJ

Bayes' Theorem, Inflation, and the Ekpyrotic Universe

Authors: Joseph Wraga, David M. Goldberg

Abstract: We calculate the Bayesian evidences for a class of Ekpyrotic universe models, and compare with a model of single field inflation with a Higgs-type potential. Combining parsimony and observational constraints, this gives us a systematic way to evaluate the degree to which Ekpyrotic models are constrained by CMB data from Planck. We integrate the equations of motion numerically to define a likelihoo… ▽ More We calculate the Bayesian evidences for a class of Ekpyrotic universe models, and compare with a model of single field inflation with a Higgs-type potential. Combining parsimony and observational constraints, this gives us a systematic way to evaluate the degree to which Ekpyrotic models are constrained by CMB data from Planck. We integrate the equations of motion numerically to define a likelihood using Planck 2018 data and sample this likelihood to obtain Bayesian evidences. Priors are justified and used to put Ekpyrotic models and inflation on equal footing. We find reasonable preference for one of the considered Ekpyrotic models over the others, but that even this one is disfavored compared with Higgs inflation. △ Less

Submitted 3 March, 2022; v1 submitted 21 January, 2022; originally announced January 2022.

Comments: 24 pages, 6 figures; Added references and corrected typos

A New Estimate of Galaxy Mass-to-Light Ratios from Flexion Lensing Statistics

Authors: Joseph Michael Fabritius II, David M Goldberg

Abstract: We perform a flexion based weak gravitational analysis of the first two Hubble Frontier Field clusters: Abell 2744 and MACS 0416. A parametric method for using projected flexion signals as a probe of cluster member mass is described in detail. The normalization and slope of a $L-θ_E$ (as a proxy for $L-σ$) scaling relation in each cluster is determined using measured flexion signals. A parallel fi… ▽ More We perform a flexion based weak gravitational analysis of the first two Hubble Frontier Field clusters: Abell 2744 and MACS 0416. A parametric method for using projected flexion signals as a probe of cluster member mass is described in detail. The normalization and slope of a $L-θ_E$ (as a proxy for $L-σ$) scaling relation in each cluster is determined using measured flexion signals. A parallel field analysis is undertaken concurrently to provide a baseline measure of method effectiveness. We find an agreement in the Faber-Jackson slope $\ell$ associated with galaxy age and morphology for both clusters, as well as theoretical distinction in the cluster normalization mass. △ Less

Submitted 27 July, 2022; v1 submitted 11 August, 2021; originally announced August 2021.

A Tidal Origin for a 3-body Resonance in Kepler-221

Authors: Max Goldberg, Konstantin Batygin

Abstract: Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here we show that the four-planet system… ▽ More Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here we show that the four-planet system Kepler-221, not previously recognized to have active orbital resonances, has a three-body commensurability relation unique within the Kepler sample. Using a suite of numerical experiments as well as a perturbative analysis, we demonstrate that this system likely began as a resonant chain and proceeded to undergo large-scale divergence away from resonance, under the action of tidal dissipation. Our results further indicate that obliquity tides, driven by a secular spin-orbit resonance and mutual inclination, are an excellent candidate for driving this orbital divergence, and that the high tidal luminosity may also explain the anomalous size of planet b, which lies within the Fulton radius gap. △ Less

Submitted 16 May, 2021; originally announced May 2021.

Comments: 15 pages, 9 figures, accepted for publication in AJ

Flat Maps that improve on the Winkel Tripel

Authors: J. Richard Gott III, David M. Goldberg, Robert J. Vanderbei

Abstract: Goldberg & Gott (2008) developed six error measures to rate flat map projections on their verisimilitude to the sphere: Isotropy, Area, Flexion, Skewness, Distances, and Boundary Cuts. The first two depend on the metric of the projection, the next two on its first derivatives. By these criteria, the Winkel Tripel (used by National Geographic for world maps) was the best scoring of all the known pr… ▽ More Goldberg & Gott (2008) developed six error measures to rate flat map projections on their verisimilitude to the sphere: Isotropy, Area, Flexion, Skewness, Distances, and Boundary Cuts. The first two depend on the metric of the projection, the next two on its first derivatives. By these criteria, the Winkel Tripel (used by National Geographic for world maps) was the best scoring of all the known projections with a sum of squares of the six errors of 4.563, normalized relative to the Equirectangular in each error term. We present here a useful Gott-Wagner variant with a slightly better error score of only 4.497. We also present a radically new class of flat double-sided maps (like phonograph records) which have correct topology and vastly improved error scores: 0.881 for the azimuthal equidistant version. We believe it is the most accurate flat map of Earth yet. We also show maps of other solar system objects and sky maps. △ Less

Submitted 15 February, 2021; originally announced February 2021.

Shape, Color and Distance in Weak Gravitational Flexion

Authors: Joseph M. Fabritius II, Evan J. Arena, David M. Goldberg

Abstract: Canonically, elliptical galaxies might be expected to have a perfect rotational symmetry, making them ideal targets for flexion studies - however, this assumption hasn't been tested. We have undertaken an analysis of low and high redshift galaxy catalogs of known morphological type with a new gravitational lensing code, Lenser. Using color measurements in the u-r band and fit Sersic index values,… ▽ More Canonically, elliptical galaxies might be expected to have a perfect rotational symmetry, making them ideal targets for flexion studies - however, this assumption hasn't been tested. We have undertaken an analysis of low and high redshift galaxy catalogs of known morphological type with a new gravitational lensing code, Lenser. Using color measurements in the u-r band and fit Sersic index values, objects with characteristics consistent with early-type galaxies are found to have a lower intrinsic scatter in flexion signal than late-type galaxies. We find this measured flexion noise can be reduced by more than a factor of two at both low and high redshift. △ Less

Submitted 10 January, 2021; v1 submitted 5 June, 2020; originally announced June 2020.

Comments: 7 pages, 13 figures, submitted to MNRAS

Prospects for Refining Kepler TTV Masses using TESS Observations

Authors: Max Goldberg, Sam Hadden, Matthew J. Payne, Matthew J. Holman

Abstract: In this paper we investigate systems previously identified to exhibit transit timing variations (TTVs) in Kepler data, with the goal of predicting the expected improvements to the mass and eccentricity constraints that will arise from combining Kepler data with future data from the TESS mission. We advocate for the use of the Kullback-Leibler (KL) divergence as a means to quantify improvements in… ▽ More In this paper we investigate systems previously identified to exhibit transit timing variations (TTVs) in Kepler data, with the goal of predicting the expected improvements to the mass and eccentricity constraints that will arise from combining Kepler data with future data from the TESS mission. We advocate for the use of the Kullback-Leibler (KL) divergence as a means to quantify improvements in the measured constraints. Compared to the original Kepler data, the TESS data will have a lower signal-to-noise ratio, rendering some of the planetary transits undetectable, and lowering the accuracy with which the transit mid-time can be estimated. Despite these difficulties, out of the 55 systems (containing 143 planets) investigated, we predict that the collection of short-cadence data by TESS will be of significant value (i.e. it will improve the mass uncertainty such that the KL divergence is > 0.1) for approximately 6 - 14 planets during the nominal mission, with the range primarily driven by the uncertain precision with which transit mid-times will be recovered from TESS data. In an extended mission this would increase to a total of approximately 12 - 25 planets. △ Less

Submitted 5 October, 2018; originally announced October 2018.

Comments: 11 pages, 8 figures, submitted to AAS journals

Flexion in Abell 2744

Authors: Justin P. Bird, David M. Goldberg

Abstract: We present the first flexion-focused gravitational lensing analysis of the first of the strong-lensing "cosmic telescope" galaxy clusters, observed as part of the Hubble Frontier Fields initiative. Using HST observations of Abell 2744 (z = 0.308), we apply a modified Analytic Image Model (AIM) technique to measure source galaxy flexion and shear values at a final number density of 82 arcmin… ▽ More We present the first flexion-focused gravitational lensing analysis of the first of the strong-lensing "cosmic telescope" galaxy clusters, observed as part of the Hubble Frontier Fields initiative. Using HST observations of Abell 2744 (z = 0.308), we apply a modified Analytic Image Model (AIM) technique to measure source galaxy flexion and shear values at a final number density of 82 arcmin$^{-2}$. By using flexion data alone we are able to identify the primary mass structure aligned along the heart of the cluster in addition to a major substructure peak offset 1.43' from the cluster core. We generate two types of nonparametric reconstructions: a flexion aperture mass map, which identifies the central potential and substructure peak with mass signal-to-noise of 3.5$σ$ and 2.3$σ$ respectively; and a convergence map derived directly from the smoothed flexion field. For the primary peak we find a mass of $1.93\times10^{14}\,h^{-1}\,M_{\odot}$ within a 45" (145h$^{-1}$ kpc) aperture, and for the western substructure we find a mass of $7.12\times10^{13}\,h^{-1}\,M_{\odot}$ within a 25" (80h$^{-1}$ kpc) aperture. The associated peak velocity dispersions were determined to be $σ_v$ = 1630 km/s and $σ_v$ = 766 km/s, respectively, by fitting nonsingular isothermal sphere profiles to the flexion data. Additionally, we use simultaneous shear measurements to independently reconstruct the broader cluster mass structure, and find that it is unable to reproduce the small-scale structure associated with the flexion reconstructions. Finally, we perform the same analysis on the Abell 2744 parallel sky field, and find no strong phantom signals in the noise reconstructions. △ Less

Submitted 23 November, 2016; originally announced November 2016.

Comments: 12 pages, 5 figures. Submitted to MNRAS

The Galaxy Cluster Concentration-Mass Scaling Relation

Authors: Austen M. Groener, David M. Goldberg, Mauro Sereno

Abstract: Scaling relations of clusters have made them particularly important cosmological probes of structure formation. In this work, we present a comprehensive study of the relation between two profile observables, concentration ($\mathrm{c_{vir}}$) and mass ($\mathrm{M_{vir}}$). We have collected the largest known sample of measurements from the literature which make use of one or more of the following… ▽ More Scaling relations of clusters have made them particularly important cosmological probes of structure formation. In this work, we present a comprehensive study of the relation between two profile observables, concentration ($\mathrm{c_{vir}}$) and mass ($\mathrm{M_{vir}}$). We have collected the largest known sample of measurements from the literature which make use of one or more of the following reconstruction techniques: Weak gravitational lensing (WL), strong gravitational lensing (SL), Weak+Strong Lensing (WL+SL), the Caustic Method (CM), Line-of-sight Velocity Dispersion (LOSVD), and X-ray. We find that the concentration-mass (c-M) relation is highly variable depending upon the reconstruction technique used. We also find concentrations derived from dark matter only simulations (at approximately $\mathrm{M_{vir} \sim 10^{14} M_{\odot}}$) to be inconsistent with the WL and WL+SL relations at the $\mathrm{1σ}$ level, even after the projection of triaxial halos is taken into account. However, to fully determine consistency between simulations and observations, a volume-limited sample of clusters is required, as selection effects become increasingly more important in answering this. Interestingly, we also find evidence for a steeper WL+SL relation as compared to WL alone, a result which could perhaps be caused by the varying shape of cluster isodensities, though most likely reflects differences in selection effects caused by these two techniques. Lastly, we compare concentration and mass measurements of individual clusters made using more than one technique, highlighting the magnitude of the potential bias which could exist in such observational samples. △ Less

Submitted 20 November, 2015; v1 submitted 7 October, 2015; originally announced October 2015.

Comments: 34 pages, 7 figures, 4 tables

Journal ref: MNRAS (January 01, 2016) Volume 455, Issue 1: 892-919

Shape Profiles and Orientation Bias for Weak and Strong Lensing Cluster Halos

Authors: A. M. Groener, D. M. Goldberg

Abstract: We study the intrinsic shape and alignment of isodensities of galaxy cluster halos extracted from the MultiDark MDR1 cosmological simulation. We find that the simulated halos, are extremely prolate on small scales, and increasingly spherical on larger ones. Due to this trend, analytical projection along the line of sight produces an overestimate of the concentration index as a decreasing function… ▽ More We study the intrinsic shape and alignment of isodensities of galaxy cluster halos extracted from the MultiDark MDR1 cosmological simulation. We find that the simulated halos, are extremely prolate on small scales, and increasingly spherical on larger ones. Due to this trend, analytical projection along the line of sight produces an overestimate of the concentration index as a decreasing function of radius, which we quantify by using both the intrinsic distribution of 3D concentrations ($c_{200}$) and isodensity shape on weak and strong lensing scales. We find this difference to be $\sim 18\%$ ($\sim 9\%$) for low (medium) mass cluster halos with intrinsically low concentrations ($c_{200}=1-3$), while we find virtually no difference for halos with intrinsically high concentrations. Isodensities are found to be fairly well-aligned throughout the entirety of the radial scale of each halo population. However, major axes of individual halos have been found to deviate by as much as $\sim 30^{\circ}$. We also present a value-added catalog of our analysis results, which we have made publicly available to download. △ Less

Submitted 22 July, 2014; v1 submitted 8 May, 2014; originally announced May 2014.

Comments: 9 pages, 6 figures, 3 tables

Journal ref: The Astrophysical Journal, 2014, 795, 153

Mass Reconstruction using Particle Based Lensing II: Quantifying substructure with Strong+Weak lensing and X-rays

Authors: Sanghamitra Deb, Andrea Morandi, Kristian Pedersen, Signe Riemer-Sorensen, David M. Goldberg, Hakon Dahle

Abstract: We report a mass reconstruction of A1689 using Particle Based Lensing (PBL), a new technique for Strong+Weak lensing that allows a variable resolution depending on the data density and the signal-to-noise. Using PBL we also calculate the covariance matrix for the resulting mass map. The reconstruction of A1689 shows a secondary mass peak in the north-east direction confirming previous optical obse… ▽ More We report a mass reconstruction of A1689 using Particle Based Lensing (PBL), a new technique for Strong+Weak lensing that allows a variable resolution depending on the data density and the signal-to-noise. Using PBL we also calculate the covariance matrix for the resulting mass map. The reconstruction of A1689 shows a secondary mass peak in the north-east direction confirming previous optical observations. This indicates that the central region of the cluster is still undergoing a weak merger. We have used this mass map to measure power ratios of the dark matter distribution and compared it to the X-ray distribution. We find that the power in the X-ray distribution is lower suggesting a smoother and rounder gas distribution compared to the dark matter distribution. We fitted an NFW profile to the profile derived from the mass map and we find that the lensing mass within 1 Mpc is $1.5\pm0.33\times 10^{15}M_\odot$. This is higher than the X-ray mass. △ Less

Submitted 17 January, 2012; originally announced January 2012.

Comments: submitted to APJ. 21 pages, 6 figures

Report number: ANL-HEP-PR-12-3

New Constraints on the Complex Mass Substructure in Abell 1689 from Gravitational Flexion

Authors: Adrienne Leonard, Lindsay J. King, David M. Goldberg

Abstract: In a recent publication, the flexion aperture mass statistic was found to provide a robust and effective method by which substructure in galaxy clusters might be mapped. Moreover, we suggested that the masses and mass profile of structures might be constrained using this method. In this paper, we apply the flexion aperture mass technique to HST ACS images of Abell 1689. We demonstrate that the fle… ▽ More In a recent publication, the flexion aperture mass statistic was found to provide a robust and effective method by which substructure in galaxy clusters might be mapped. Moreover, we suggested that the masses and mass profile of structures might be constrained using this method. In this paper, we apply the flexion aperture mass technique to HST ACS images of Abell 1689. We demonstrate that the flexion aperture mass statistic is sensitive to small-scale structures in the central region of the cluster. While the central potential is not constrained by our method, due largely to missing data in the central 0.5$^\prime$ of the cluster, we are able to place constraints on the masses and mass profiles of prominent substructures. We identify 4 separate mass peaks, and use the peak aperture mass signal and zero signal radius in each case to constrain the masses and mass profiles of these substructures. The three most massive peaks exhibit complex small-scale structure, and the masses indicated by the flexion aperture mass statistic suggest that these three peaks represent the dominant substructure component of the cluster ($\sim 7\times 10^{14}h^{-1}M_\odot$). Their complex structure indicates that the cluster -- far from being relaxed -- may have recently undergone a merger. The smaller, subsidiary peak is located coincident with a group of galaxies within the cluster, with mass $\sim 1\times10^{14}h^{-1}M_\odot$. These results are in excellent agreement with previous substructure studies of this cluster. △ Less

Submitted 7 December, 2010; v1 submitted 6 September, 2010; originally announced September 2010.

Comments: 18 pages, 10 figures, MNRAS accepted (7 Dec 2010)

Analyzing the Flux Anomalies of the Large-Separation Lensed Quasar SDSS J1029+2623

Authors: Rachael M. Kratzer, Gordon T. Richards, David M. Goldberg, Masamune Oguri, Christopher S. Kochanek, Jacqueline A. Hodge, Robert H. Becker, Naohisa Inada

Abstract: Using a high resolution radio image, we successfully resolve the two fold image components B and C of the quasar lens system SDSS J1029+2623. The flux anomalies associated with these two components in the optical regime persist, albeit less strongly, in our radio observations, suggesting that the cluster must be modeled by something more than a single central potential. We argue that placing subst… ▽ More Using a high resolution radio image, we successfully resolve the two fold image components B and C of the quasar lens system SDSS J1029+2623. The flux anomalies associated with these two components in the optical regime persist, albeit less strongly, in our radio observations, suggesting that the cluster must be modeled by something more than a single central potential. We argue that placing substructure close to one of the components can account for a flux anomaly with negligible changes in the component positions. Our best fit model has a substructure mass of ~10^8 solar masses up to the mass-sheet degeneracy, located roughly 0.1 arcsecs West and 0.1 arcsecs North of component B. We demonstrate that a positional offset between the centers of the source components can explain the differences between the optical and radio flux ratios. △ Less

Submitted 6 January, 2011; v1 submitted 13 August, 2010; originally announced August 2010.

Comments: 5 pages, 3 figures, 2 tables, accepted to ApJ Letters

Measuring dark matter ellipticity of Abell 901/902 using Particle Based Lensing

Authors: Sanghamitra Deb, David M. Goldberg, Catherine Heymans, Andrea Morandi

Abstract: We present a non-parametric measure of the ellipticity and the alignment of the dark matter halos in Abell 901/902 supercluster. This super-cluster is a system of four separate peaks in a $0.5^{\circ}\times0.5^{\circ}$ field of view. We map the mass distribution of each individual peak using an improved version of Particle Based Lensing (PBL) and measure the ellipticity of the dark matter halos… ▽ More We present a non-parametric measure of the ellipticity and the alignment of the dark matter halos in Abell 901/902 supercluster. This super-cluster is a system of four separate peaks in a $0.5^{\circ}\times0.5^{\circ}$ field of view. We map the mass distribution of each individual peak using an improved version of Particle Based Lensing (PBL) and measure the ellipticity of the dark matter halos associated with two of the peaks directly from the mass map and by fitting them to a singular isothermal ellipse. The parametric and non-parametric measurements are consistent for A901b while the position angle for the Southwest Group is different for the two techniques. We account for this discrepancy to substructure present in the Southwest Peak. We estimate an axis ratio of $0.37\pm 0.1$ for A901b and $0.54^{+0.08}_{-0.09}$ for the Southwest Group. △ Less

Submitted 21 December, 2009; originally announced December 2009.

Comments: submitted to APJ, 25 pages, 13 figures

Journal ref: Astrophys.J.721:124-136,2010

Fold Lens Flux Anomalies: A Geometric Approach

Authors: David M. Goldberg, Mary K. Chessey, Wendy B. Harris, Gordon T. Richards

Abstract: We develop a new approach for studying flux anomalies in quadruply-imaged fold lens systems. We show that in the absence of substructure, microlensing, or differential absorption, the expected flux ratios of a fold pair can be tightly constrained using only geometric arguments. We apply this technique to 11 known quadruple lens systems in the radio and infrared, and compare our estimates to the… ▽ More We develop a new approach for studying flux anomalies in quadruply-imaged fold lens systems. We show that in the absence of substructure, microlensing, or differential absorption, the expected flux ratios of a fold pair can be tightly constrained using only geometric arguments. We apply this technique to 11 known quadruple lens systems in the radio and infrared, and compare our estimates to the Monte Carlo based results of Keeton, Gaudi, and Petters. We show that a robust estimate for a flux ratio from a smoothly varying potential can be found, and at long wavelengths those lenses deviating from from this ratio almost certainly contain significant substructure. △ Less

Submitted 5 May, 2010; v1 submitted 4 December, 2009; originally announced December 2009.

Comments: 16 pages, including 8 figures

Journal ref: Astrophys.J.715:793,2010; Astrophys.J.715:793-802,2010

Beowulf Analysis Symbolic INterface BASIN: Interactive Parallel Data Analysis for Everyone

Authors: Enrico Vesperini, David M. Goldberg, Stephen L. W. McMillan, James Dura, Douglas Jones

Abstract: The advent of affordable parallel computers such as Beowulf PC clusters and, more recently, of multi-core PCs has been highly beneficial for a large number of scientists and smaller institutions that might not otherwise have access to substantial computing facilities. However, there has not been an analogous progress in the development and dissemination of parallel software: scientists need the… ▽ More The advent of affordable parallel computers such as Beowulf PC clusters and, more recently, of multi-core PCs has been highly beneficial for a large number of scientists and smaller institutions that might not otherwise have access to substantial computing facilities. However, there has not been an analogous progress in the development and dissemination of parallel software: scientists need the expertise to develop parallel codes and have to invest a significant amount of time in the development of tools even for the most common data analysis tasks. We describe the Beowulf Analysis Symbolic INterface (BASIN) a multi-user parallel data analysis and visualization framework. BASIN is aimed at providing scientists with a suite of parallel libraries for astrophysical data analysis along with general tools for data distribution and parallel operations on distributed data to allow them to easily develop new parallel libraries for their specific tasks. △ Less

Submitted 29 April, 2008; originally announced April 2008.

Comments: 12 pages, 3 figures, Submitted to Computing in Science and Engineering-special issue on Computational Astrophysics

Reconstruction of Cluster Masses using Particle Based Lensing I: Application to Weak Lensing

Authors: Sanghamitra Deb, David M. Goldberg, Vede J. Ramdass

Abstract: We present Particle-Based Lensing (PBL), a new technique for gravitational lensing mass reconstructions of galaxy clusters. Traditionally, most methods have employed either a finite inversion or gridding to turn observational lensed galaxy ellipticities into an estimate of the surface mass density of a galaxy cluster. We approach the problem from a different perspective, motivated by the success… ▽ More We present Particle-Based Lensing (PBL), a new technique for gravitational lensing mass reconstructions of galaxy clusters. Traditionally, most methods have employed either a finite inversion or gridding to turn observational lensed galaxy ellipticities into an estimate of the surface mass density of a galaxy cluster. We approach the problem from a different perspective, motivated by the success of multi-scale analysis in smoothed particle hydrodynamics. In PBL, we treat each of the lensed galaxies as a particle and then reconstruct the potential by smoothing over a local kernel with variable smoothing scale. In this way, we can tune a reconstruction to produce constant signal-noise throughout, and maximally exploit regions of high information density. PBL is designed to include all lensing observables, including multiple image positions and fluxes from strong lensing, as well as weak lensing signals including shear and flexion. In this paper, however, we describe a shear-only reconstruction, and apply the method to several test cases, including simulated lensing clusters, as well as the well-studied ``Bullet Cluster'' (1E0657-56). In the former cases, we show that PBL is better able to identify cusps and substructures than are grid-based reconstructions, and in the latter case, we show that PBL is able to identify substructure in the Bullet Cluster without even exploiting strong lensing measurements. We also make our codes publicly available. △ Less

Submitted 13 June, 2008; v1 submitted 31 January, 2008; originally announced February 2008.

Comments: Accepted for publication in ApJ; Codes available at http://www.physics.drexel.edu/~deb/PBL.htm ; 12 pages,9 figures, section 3 shortened

Weak lensing ellipticities in a strong lensing regime

Authors: Richard Massey, David M. Goldberg

Abstract: It is now routine to measure the weak gravitational lensing shear signal from the mean ellipticity of distant galaxies. However, conversion between ellipticity and shear assumes local linearity of the lensing potential (ie that the spatial derivatives of the shear are small), and this condition is not satisfied in some of the most interesting regions of the sky. We extend a derivation of lensing… ▽ More It is now routine to measure the weak gravitational lensing shear signal from the mean ellipticity of distant galaxies. However, conversion between ellipticity and shear assumes local linearity of the lensing potential (ie that the spatial derivatives of the shear are small), and this condition is not satisfied in some of the most interesting regions of the sky. We extend a derivation of lensing equations to include higher order terms, and assess the level of biases introduced by assuming that first-order weak lensing theory holds in a relatively strong shear regime. We find that, even in a worst-case scenario, a fully linear analysis is accurate to within 1% outside ~1.07 times the Einstein radius of a lens, by deriving an analytic function that can be used to estimate the applicability of any first-order analysis. The effect is too small to explain the discrepancy between weak- and strong-lensing estimates of the mass of the bullet cluster, and should not impact cluster surveys for the forseeable future. In fact, it means that arclets can be used to measure shears closer to a cluster core than has been generally appreciated. However, at the level of accuracy demanded by future lensing surveys, this bias ought to be considered in measurements of the inner slope of cluster mass distributions and the small-scale end of the mass power spectrum. Both of these are central in determining the relationship between baryonic and dark matter. △ Less

Submitted 5 January, 2009; v1 submitted 10 September, 2007; originally announced September 2007.

Comments: 5 pages, 1 figure. Matches published version

Journal ref: Astrophys.J.673:L111-L114,2008

Gravitational Shear, Flexion and Strong Lensing in Abell 1689

Authors: Adrienne Leonard, David M. Goldberg, Jason L. Haaga, Richard Massey

Abstract: We present a gravitational lensing analysis of the galaxy cluster Abell 1689, incorporating measurements of the weak shear, flexion, and strong lensing induced in background galaxies. This is the first time that a shapelet technique has been used to reconstruct the distribution of mass in this cluster, and the first time that a flexion signal has been measured using cluster members as lenses. Fr… ▽ More We present a gravitational lensing analysis of the galaxy cluster Abell 1689, incorporating measurements of the weak shear, flexion, and strong lensing induced in background galaxies. This is the first time that a shapelet technique has been used to reconstruct the distribution of mass in this cluster, and the first time that a flexion signal has been measured using cluster members as lenses. From weak shear measurements alone, we generate a non-parametric mass reconstruction, which shows significant substructure corresponding to groups of galaxies within the cluster. Additionally, our galaxy-galaxy flexion signal demonstrates that the cluster galaxies can be well-fit by a singular isothermal sphere model with a characteristic velocity dispersion of $σ_v = 295\pm 40 km/s $. We identify a major, distinct dark matter clump, offset by 40$h^{-1}$kpc from the central cluster members, which was not apparent from shear measurements alone. This secondary clump is present in a parametric mass reconstruction using flexion data alone, and its existence is suggested in a non-parametric reconstruction of the cluster using a combination of strong and weak lensing. As found in previous studies, the mass profile obtained by combining weak and strong lensing data shows a much steeper profile than that obtained from only weak lensing data. △ Less

Submitted 7 June, 2007; v1 submitted 9 February, 2007; originally announced February 2007.

Comments: 17 pages, 12 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.666:51-63,2007

Flexion and Skewness in Map Projections of the Earth

Authors: David M. Goldberg, J. Richard Gott III

Abstract: Tissot indicatrices have provided visual measures of local area and isotropy distortions. Here we show how large scale distortions of flexion (bending) and skewness (lopsidedness) can be measured. Area and isotropy distortions depend on the map projection metric, flexion and skewness, which manifest themselves on continental scales, depend on the first derivatives of the metric. We introduce new… ▽ More Tissot indicatrices have provided visual measures of local area and isotropy distortions. Here we show how large scale distortions of flexion (bending) and skewness (lopsidedness) can be measured. Area and isotropy distortions depend on the map projection metric, flexion and skewness, which manifest themselves on continental scales, depend on the first derivatives of the metric. We introduce new indicatrices that show not only area and isotropy distortions but flexion and skewness as well. We present a table showing error measures for area, isotropy, flexion, skewness, distances, and boundary cuts allowing us to compare different world map projections. We find that the Winkel-Tripel projection (already adopted for world maps by the National Geographic), has low distortion on most measures and excellent quality overall. △ Less

Submitted 12 April, 2007; v1 submitted 23 August, 2006; originally announced August 2006.

Comments: 31 pages, including 27 postscript figures. Accepted to Cartographica. Detailed discussion and code at http://www.physics.drexel.edu/~goldberg/projections

Measuring Flexion

Authors: David M. Goldberg, Adrienne Leonard

Abstract: We describe practical approaches to measuring flexion in observed galaxies. In particular, we look at the issues involved in using the Shapelets and HOLICs techniques as means of extracting 2nd order lensing information. We also develop an extension of HOLICs to estimate flexion in the presence of noise, and with a nearly isotropic PSF. We test both approaches in simple simulated lenses as well… ▽ More We describe practical approaches to measuring flexion in observed galaxies. In particular, we look at the issues involved in using the Shapelets and HOLICs techniques as means of extracting 2nd order lensing information. We also develop an extension of HOLICs to estimate flexion in the presence of noise, and with a nearly isotropic PSF. We test both approaches in simple simulated lenses as well as a sample of possible background sources from ACS observations of A1689. We find that because noise is weighted differently in shapelets and HOLICs approaches, that the correlation between measurements of the same object is somewhat diminished, but produce similar scatter due to measurement noise. △ Less

Submitted 23 April, 2007; v1 submitted 26 July, 2006; originally announced July 2006.

Comments: 10 pages, including 13 postscript figures

Journal ref: Astrophys.J.660:1003-1015,2007

Weak Gravitational Flexion

Authors: D. J. Bacon, D. M. Goldberg, B. T. P. Rowe, A. N. Taylor

Abstract: Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy-galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe t… ▽ More Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy-galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy-galaxy lensing, examining isothermal sphere and Navarro, Frenk & White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in 2-D projection on the sky, and in 3-D in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence-flexion cross-correlations, and find that the signal is an effective probe of structure on non-linear scales. △ Less

Submitted 21 April, 2005; originally announced April 2005.

Comments: 17 pages, including 12 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 365 (2006) 414-428

Mapping the Cosmic Web with the Sloan Digital Sky Survey

Authors: Michael S. Vogeley, Fiona Hoyle, Randall R. Rojas, David M. Goldberg

Abstract: Wide-angle, moderately deep redshift surveys such as that conducted as part of the Sloan Digital Sky Survey (SDSS) allow study of the relationship between the structural elements of the large-scale distribution of galaxies -- including groups, cluster, superclusters, and voids -- and the dependence of galaxy formation and evolution on these enviroments. We present a progress report on mapping ef… ▽ More Wide-angle, moderately deep redshift surveys such as that conducted as part of the Sloan Digital Sky Survey (SDSS) allow study of the relationship between the structural elements of the large-scale distribution of galaxies -- including groups, cluster, superclusters, and voids -- and the dependence of galaxy formation and evolution on these enviroments. We present a progress report on mapping efforts with the SDSS and discuss recently constructed catalogs of clusters, voids, and void galaxies, and evidence for a 420Mpc/h supercluster or ``Great Wall.'' Analysis of multi-band photometry and moderate-resolution spectroscopy from the SDSS reveals environmental dependence of the star formation history of galaxies that extends over more than a factor of 100 in density, from clusters all the way to the deep interiors of voids. On average, galaxies in the rarified environments of voids exhibit bluer colors, higher specific star formation rates, lower dust content, and more disk-like morphology than objects in denser regions. This trend persists in comparisons of samples in low vs. high-density regions with similar luminosity and morphology, thus this dependence is not simply an extension of the morphology-density relation. Large-scale modulation of the halo mass function and the temperature of the intergalactic medium might explain this dependence of galaxy evolution on the large-scale environment. △ Less

Submitted 31 August, 2004; originally announced August 2004.

Comments: 7 pages, invited review to appear in Outskirts of Galaxy Clusters: intense life in the suburbs, Proceedings IAU COlloquium No. 195, 2004, A. Diaferio, ed

The Mass Function of Void Galaxies in the SDSS Data Release 2

Authors: David M. Goldberg, Timothy D. Jones, Fiona Hoyle, Randall R. Rojas, Michael S. Vogeley, Michael R. Blanton

Abstract: We estimate the Mass Function of void galaxies in the second public data release of the Sloan Digital Sky Survey from a sample of 1000 galaxies with local density contrasts of delta_v < -0.6. The galaxy sample is split into ellipticals and spirals using a color-Sersic index criteria. We estimate the virial masses of ellipticals using the measured spectral line-widths along with the observed size… ▽ More We estimate the Mass Function of void galaxies in the second public data release of the Sloan Digital Sky Survey from a sample of 1000 galaxies with local density contrasts of delta_v < -0.6. The galaxy sample is split into ellipticals and spirals using a color-Sersic index criteria. We estimate the virial masses of ellipticals using the measured spectral line-widths along with the observed size. Projection effects and uncertainties in halo properties make mass estimates of spirals more difficult. We use an inversion of the Tully-Fisher relation to estimate the isothermal rotational velocity, and introduce a scaling factor to estimate the halo extent. We then fit the measured mass function against a theoretical Press-Schechter model, and find that the distribution of galaxies in voids appears to be nearly unbiased compared to the mass. △ Less

Submitted 11 November, 2004; v1 submitted 23 June, 2004; originally announced June 2004.

Comments: 9 pages, 5 postscript figures. Accepted to Astrophysical Journal; Revisions in response to referee comments

Journal ref: Astrophys.J. 621 (2005) 643-650

Galaxy-Galaxy Flexion: Weak Lensing to Second Order

Authors: David M. Goldberg, David J. Bacon

Abstract: In this paper, we develop a new gravitational lensing inversion technique. While traditional approaches assume that the lensing field varies little across a galaxy image, we note that this variation in the field can give rise to a ``Flexion'' or bending of a galaxy image, which may then be used to detect a lensing signal with increased signal to noise. Since the significance of the Flexion signa… ▽ More In this paper, we develop a new gravitational lensing inversion technique. While traditional approaches assume that the lensing field varies little across a galaxy image, we note that this variation in the field can give rise to a ``Flexion'' or bending of a galaxy image, which may then be used to detect a lensing signal with increased signal to noise. Since the significance of the Flexion signal increases on small scales, this is ideally suited to galaxy-galaxy lensing. We develop an inversion technique based on the ``Shapelets'' formalism of Refregier (2003). We then demonstrate the proof of this concept by measuring a Flexion signal in the Deep Lens Survey. Assuming an intrinsically isothermal distribution, we find from the Flexion signal alone a velocity width of v_c=221\pm 12 km/s for lens galaxies of r < 21.5, subject to uncertainties in the intrinsic Flexion distribution. △ Less

Submitted 15 October, 2004; v1 submitted 16 June, 2004; originally announced June 2004.

Comments: 11 pages, Latex, 4 figures. Accepted by ApJ, changes include revision of errors from previous draft

Journal ref: Astrophys.J. 619 (2005) 741-748

Simulating Voids

Authors: David M. Goldberg, Michael S. Vogeley

Abstract: We present a novel method for simulation of the interior of large cosmic voids, suitable for study of the formation and evolution of objects lying within such regions. Following Birkhoff's theorem, void regions dynamically evolve as universes with cosmological parameters that depend on the underdensity of the void. We derive the values of $Ω_M$, $Ω_Λ$, and $H_0$ that describe this evolution. We… ▽ More We present a novel method for simulation of the interior of large cosmic voids, suitable for study of the formation and evolution of objects lying within such regions. Following Birkhoff's theorem, void regions dynamically evolve as universes with cosmological parameters that depend on the underdensity of the void. We derive the values of $Ω_M$, $Ω_Λ$, and $H_0$ that describe this evolution. We examine how the growth rate of structure and scale factor in a void differ from the background universe. Together with a prescription for the power spectrum of fluctuations, these equations provide the initial conditions for running specialized void simulations. The increased efficiency of such simulations, in comparison with general-purpose simulations, allows an improvement of upwards of twenty in the mass resolution. As a sanity check, we run a moderate resolution simulation ($N=128^3$ particles) and confirm that the resulting mass function of void halos is consistent with other theoretical and numerical models. △ Less

Submitted 12 December, 2003; v1 submitted 9 July, 2003; originally announced July 2003.

Comments: Accepted to ApJ., Dec. 11, 2003. 9 pages, 3 postscript figures. Changes made to respond to referee's comments

Journal ref: Astrophys.J.605:1-6,2004

The Galaxy Octopole Moment as a Probe of Weak Lensing Shear Fields

Authors: D. M. Goldberg, P. Natarajan

Abstract: In this paper, we introduce the octopole moment of the light distribution in galaxies as a probe of the weak lensing shear field. While traditional ellipticity estimates of the local shear derived from the quadrupole moment are limited by the width of the intrinsic ellipticity distribution of background galaxies, the dispersion in the intrinsic octopole distribution is expected to be much smalle… ▽ More In this paper, we introduce the octopole moment of the light distribution in galaxies as a probe of the weak lensing shear field. While traditional ellipticity estimates of the local shear derived from the quadrupole moment are limited by the width of the intrinsic ellipticity distribution of background galaxies, the dispersion in the intrinsic octopole distribution is expected to be much smaller, implying that the signal from this higher order moment is ultimately limited by measurement noise, and not by intrinsic scatter. We present the computation of the octopole moment and show that current observations are at the regime where the octopole estimates will soon be able to contribute to the overall accuracy of the estimates of local shear fields. Therefore, the prospects for this estimator from future datasets like the Advanced Camera for Survey and the Next Generation Space Telescope are very promising. △ Less

Submitted 10 July, 2001; originally announced July 2001.

Comments: 9 pages, 2 PostScript figures; Submitted to Astrophysical Journal

Journal ref: Astrophys.J.564:65,2002

Using Perturbative Least Action to Reconstruct the Local Group

Authors: David M. Goldberg

Abstract: In this paper, we apply the Perturbative Least Action Method to model the Local Group of galaxies for various Cosmological models. We show that though the galaxy masses are theoretically good discriminators of $Ω_M$ given some observed MW and M31 separation and radial velocity, current estimates of the masses are insufficient to make any cosmological claims. We then discuss additional complicati… ▽ More In this paper, we apply the Perturbative Least Action Method to model the Local Group of galaxies for various Cosmological models. We show that though the galaxy masses are theoretically good discriminators of $Ω_M$ given some observed MW and M31 separation and radial velocity, current estimates of the masses are insufficient to make any cosmological claims. We then discuss additional complications to expand this analysis. △ Less

Submitted 6 September, 2000; v1 submitted 4 September, 2000; originally announced September 2000.

Comments: 18 Pages LaTeX, including 5 postscript figures. Submitted to Astrophysical Journal

Journal ref: Astrophys. J.550:87,2001

Using Perturbative Least Action to Reconstruct Redshift Space Distortions

Authors: David M. Goldberg

Abstract: In this paper, we present a redshift space reconstruction scheme which is analogous to and extends the Perturbative Least Action (PLA) method described by Goldberg & Spergel (2000). We first show that this scheme is effective in reconstructing even nonlinear observations. We then suggest that by varying the cosmology to minimize the quadrupole moment of a reconstructed density field, it may be p… ▽ More In this paper, we present a redshift space reconstruction scheme which is analogous to and extends the Perturbative Least Action (PLA) method described by Goldberg & Spergel (2000). We first show that this scheme is effective in reconstructing even nonlinear observations. We then suggest that by varying the cosmology to minimize the quadrupole moment of a reconstructed density field, it may be possible to lower the errorbars on the redshift distortion parameter, $β$ as well as to break the degeneracy between the linear bias parameter, $b$, and $Ω_M$. Finally, we discuss how PLA might be applied to realistic redshift surveys. △ Less

Submitted 17 August, 2000; originally announced August 2000.

Comments: 34 Pages LaTeX, including 10 postscript figures. Submitted to Astrophysical Journal

Journal ref: Astrophys.J.552:413,2001

Using Perturbative Least Action to Recover Cosmological Initial Conditions

Authors: David M. Goldberg, David N. Spergel

Abstract: We introduce a new method for generating initial conditions consistent with highly nonlinear observations of density and velocity fields. Using a variant of the Least Action method, called Perturbative Least Action (PLA), we show that it is possible to generate several different sets of initial conditions, each of which will satisfy a set of highly nonlinear observational constraints at the pres… ▽ More We introduce a new method for generating initial conditions consistent with highly nonlinear observations of density and velocity fields. Using a variant of the Least Action method, called Perturbative Least Action (PLA), we show that it is possible to generate several different sets of initial conditions, each of which will satisfy a set of highly nonlinear observational constraints at the present day. We then discuss a code written to test and apply this method and present the results of several simulations. △ Less

Submitted 17 April, 2000; v1 submitted 19 December, 1999; originally announced December 1999.

Comments: 24 pages, 6 postscript figures. Accepted for publication in Astrophysical Journal

Journal ref: Astrophys.J.544:21,2000

Using Perturbative Least Action to Run N Body Simulations Back in Time

Authors: D. M. Goldberg, D. N. Spergel

Abstract: In this report, we present a new method for reconstructing N body initial conditions from a proscribed final density field. This method, Perturbative Least Action (PLA) is similar to traditional least action approaches, except that orbits of particles are found as expansions around previously determined and physically motivated orbits. In this report, we present a new method for reconstructing N body initial conditions from a proscribed final density field. This method, Perturbative Least Action (PLA) is similar to traditional least action approaches, except that orbits of particles are found as expansions around previously determined and physically motivated orbits. △ Less

Submitted 2 September, 1999; originally announced September 1999.

Comments: 4 pages, 1 figure, Cosmic Flows Workshop, Victoria BC, July 1999, S. Courteau, M. Strauss, and J. Willick eds

A Comparison of Simple Mass Estimators for Galaxy Clusters

Authors: Tereasa G. Brainerd, Candace Oaxaca Wright, David M. Goldberg, Jens Verner Villumsen

Abstract: High-resolution N-body simulations are used to investigate systematic trends in the mass profiles and total masses of clusters as derived from 3 simple estimators: (1) the weak gravitational lensing shear field under the assumption of an isothermal cluster potential, (2) the dynamical mass obtained from the measured velocity dispersion under the assumption of an isothermal cluster potential, and… ▽ More High-resolution N-body simulations are used to investigate systematic trends in the mass profiles and total masses of clusters as derived from 3 simple estimators: (1) the weak gravitational lensing shear field under the assumption of an isothermal cluster potential, (2) the dynamical mass obtained from the measured velocity dispersion under the assumption of an isothermal cluster potential, and (3) the classical virial estimator. The clusters consist of order 2.5e+05 particles of mass $m_p \simeq 10^{10} \Msun$, have triaxial mass distributions, and significant substructure exists within their virial radii. Not surprisingly, the level of agreement between the mass profiles obtained from the various estimators and the actual mass profiles is found to be scale-dependent. The virial estimator yields a good measurement of the total cluster mass, though it is systematically underestimated by of order 10%. This result suggests that, at least in the limit of ideal data, the virial estimator is quite robust to deviations from pure spherical symmetry and the presence of substructure. The dynamical mass estimate based upon a measurement of the cluster velocity dispersion and an assumption of an isothermal potential yields a poor measurement of the total mass. The weak lensing estimate yields a very good measurement of the total mass, provided the mean shear used to determine the equivalent cluster velocity dispersion is computed from an average of the lensing signal over the entire cluster (i.e. the mean shear is computed interior to the virial radius). [abridged] △ Less

Submitted 3 March, 1999; originally announced March 1999.

Comments: Accepted for publication in The Astrophysical Journal. Complete paper, including 3 large colour figures can also be obtained from http://bu-ast.bu.edu/~brainerd/preprints/

The Microwave Background Bispectrum, Paper I: Basic Formalism

Authors: David N. Spergel, David M. Goldberg

Abstract: In this paper, we discuss the potential importance of measuring the CMB anisotropy bispectrum. We develop a formalism for computing the bispectrum and for measuring it from microwave background maps. As an example, we compute the bispectrum resulting from the 2nd order Rees-Sciama effect, and find that is undetectable with current and upcoming missions. In this paper, we discuss the potential importance of measuring the CMB anisotropy bispectrum. We develop a formalism for computing the bispectrum and for measuring it from microwave background maps. As an example, we compute the bispectrum resulting from the 2nd order Rees-Sciama effect, and find that is undetectable with current and upcoming missions. △ Less

Submitted 4 January, 1999; v1 submitted 16 November, 1998; originally announced November 1998.

Comments: 18 Pages, 3 Postscript Figures; Minor changes in response to referee's report

Journal ref: Phys. Rev. D 59, 103001 (1999)

The Microwave Background Bispectrum, Paper II: A Probe of the Low Redshift Universe

Authors: David M. Goldberg, David N. Spergel

Abstract: Gravitational fluctuations along the line-of-sight from the surface of last scatter to the observer distort the microwave background in several related ways: The fluctuations deflect the photon path (gravitational lensing), the decay of the gravitational potential produces additional fluctuations (ISW effect) and scattering off of hot gas in clusters produce additional fluctuations (Sunyaev-Zel'… ▽ More Gravitational fluctuations along the line-of-sight from the surface of last scatter to the observer distort the microwave background in several related ways: The fluctuations deflect the photon path (gravitational lensing), the decay of the gravitational potential produces additional fluctuations (ISW effect) and scattering off of hot gas in clusters produce additional fluctuations (Sunyaev-Zel'dovich effect). Even if the initial fluctuations generated at the surface of last scatter were Gaussian, the combination of these effects produce non-Gaussian features in the microwave sky. We discuss the microwave bispectrum as a tool for measuring a studying this signal. For MAP, we estimate that these measurements will enable us to determine the fraction of ionized gas and to probe the time evolution of the gravitational potential. △ Less

Submitted 4 January, 1999; v1 submitted 16 November, 1998; originally announced November 1998.

Comments: 15 Pages, 5 Postscript Figures; Minor Changes in Response to Referee's report

Journal ref: Phys.Rev.D59:103002,1999

Astrometric Shifts in the OGLE-1 Microlensing Events

Authors: David M. Goldberg, Przemyslaw Wozniak

Abstract: We measure the astrometric shifts of the light centroids for the microlensing events in the OGLE-1 database. Of the 15 consistently detected events, 7 had significant shifts which we were able to measure with a fair degree of confidence. Those events with large shifts are also expected to be unresolved with a background ``blend'', and thus, we suggest that we have identified events in the OGLE-1… ▽ More We measure the astrometric shifts of the light centroids for the microlensing events in the OGLE-1 database. Of the 15 consistently detected events, 7 had significant shifts which we were able to measure with a fair degree of confidence. Those events with large shifts are also expected to be unresolved with a background ``blend'', and thus, we suggest that we have identified events in the OGLE-1 catalog which are strongly blended. Though we concentrate on the OGLE-1 database, and use the DoPHOT package in order to perform our analysis, we suggest that this shift is a generic effect, and should be observable in any crowded field. △ Less

Submitted 18 December, 1997; originally announced December 1997.

Comments: 18 pages, including 6 postscript figures

Journal ref: Acta Astron. 48 (1998) 19

Using Astrometry to Deblend Microlensing Events

Authors: David M. Goldberg

Abstract: We discuss the prospect of deblending microlensing events by observing astrometric shifts of the lensed stars. Since microlensing searches are generally performed in very crowded fields, it is expected that stars will be confusion limited rather than limited by photon statistics. By performing simulations of events in crowded fields, we find that if we assume a dark lens and that the lensed star… ▽ More We discuss the prospect of deblending microlensing events by observing astrometric shifts of the lensed stars. Since microlensing searches are generally performed in very crowded fields, it is expected that stars will be confusion limited rather than limited by photon statistics. By performing simulations of events in crowded fields, we find that if we assume a dark lens and that the lensed star obeys a power law luminosity function, $n(L)\propto L^{-β}$, over half the simulated events show a measurable astrometric shift. Our simulations included 20000 stars in a $256\times 256$ Nyquist sampled CCD frame. For $β=2$, we found that 58% of the events were significantly blended $(F_{\ast}/F_{tot}\leq 0.9)$, and of those, 73% had a large astrometric shift $(\geq 0.5 pixels)$. Likewise, for $β=3$, we found that 85% of the events were significantly blended, and that 85% of those had large shifts. Moreover, the shift is weakly correlated to the degree of blending, suggesting that it may be possible not only to detect the existence of a blend, but also to deblend events statistically using shift information. △ Less

Submitted 19 August, 1997; originally announced August 1997.

Comments: 24 pages, 7 postscript Figures

Journal ref: Astrophys.J.498:156,1998

Determination of the Baryon Density from Large Scale Galaxy Redshift Surveys

Authors: David M. Goldberg, Michael A. Strauss

Abstract: We estimate the degree to which the baryon density, $Ω_{b}$, can be determined from the galaxy power spectrum measured from large scale galaxy redshift surveys, and in particular, the Sloan Digital Sky Survey. A high baryon density will cause wiggles to appear in the power spectrum, which should be observable at the current epoch. We assume linear theory on scales $\geq 20h^{-1}Mpc$ and do not i… ▽ More We estimate the degree to which the baryon density, $Ω_{b}$, can be determined from the galaxy power spectrum measured from large scale galaxy redshift surveys, and in particular, the Sloan Digital Sky Survey. A high baryon density will cause wiggles to appear in the power spectrum, which should be observable at the current epoch. We assume linear theory on scales $\geq 20h^{-1}Mpc$ and do not include the effects of redshift distortions, evolution, or biasing. With an optimum estimate of $P(k)$ to $k\sim 2π/(20 h^{-1} Mpc)$, the $1 σ$ uncertainties in $Ω_{b}$ are roughly 0.07 and 0.016 in flat and open ($Ω_{0}=0.3$) cosmological models, respectively. This result suggests that it should be possible to test for consistency with big bang nucleosynthesis estimates of $Ω_{b}$ if we live in an open universe. △ Less

Submitted 17 July, 1997; originally announced July 1997.

Comments: 23 Pages, 10 Postscript figures

Journal ref: Astrophys.J.495:29,1998

High-Resolution Simulations of Cluster Formation

Authors: Tereasa G. Brainerd, David M. Goldberg, Jens Verner Villumsen

Abstract: The formation history of rich clusters is investigated using a hybrid N-body simulation in which high spatial and mass resolution can be achieved self-consistently within a small region of a very large volume. The evolution of three massive clusters is studied via mass accretion, spherically-averaged density profiles, three-dimensional and projected shapes, and degree of substructure. Each clust… ▽ More The formation history of rich clusters is investigated using a hybrid N-body simulation in which high spatial and mass resolution can be achieved self-consistently within a small region of a very large volume. The evolution of three massive clusters is studied via mass accretion, spherically-averaged density profiles, three-dimensional and projected shapes, and degree of substructure. Each cluster consists of order $4\times 10^5$ particles at the present epoch and in the case that rich cluster evolution is well-described by a 1-parameter family, the simulations have sufficient resolution to demonstrate this. At $z=0$ the clusters have similar masses, $M(r \le 1.5h^{-1} Mpc) \sim 2\times 10^{15} h^{-1} M_\odot$, and similar spherically-averaged density profiles, however markedly different formation histories are observed. No single, dominant pattern is apparent in the time variation of the mass accretion rate, the cluster shape, or the degree of substructure. Although not a statistically large sample, these objects suggest that the detailed formation history of rich clusters cannot be characterized by a simple 1-parameter family. These results suggest that the use of observations of rich clusters over a wide range of redshifts to constrain cosmological parameters may not be entirely straightforward. △ Less

Submitted 16 June, 1997; originally announced June 1997.

Comments: 28 pages, 8 figures, plain TeX

Report number: BU97-TGB-01