Resampling images in Fourier domain

When simulating sky images, one often takes a galaxy image $F(x)$ defined by a set of pixelized samples and an interpolation kernel, and then wants to produce a new sampled image representing this galaxy as it would appear with a different point-spread function, a rotation, shearing, or magnification, and/or a different pixel scale. These operations are sometimes only possible, or most efficiently executed, as resamplings of the Fourier transform $\tilde F(u)$ of the image onto a $u$-space grid that differs from the one produced by a discrete Fourier transform (DFT) of the samples. In some applications it is essential that the resampled image be accurate to better than 1 part in $10^3$, so in this paper we first use standard Fourier techniques to show that Fourier-domain interpolation with a wrapped sinc function yields the exact value of $\tilde F(u)$ in terms of the input samples and kernel. This operation scales with image dimension as $N^4$ and can be prohibitively slow, so we next investigate the errors accrued from approximating the sinc function with a compact kernel. We show that these approximations produce a multiplicative error plus a pair of ghost images (in each dimension) in the simulated image. Standard Lanczos or cubic interpolators, when applied in Fourier domain, produce unacceptable artifacts. We find that errors $<1$ part in $10^3$ can be obtained by (1) 4-fold zero-padding of the original image before executing the $x\rightarrow u$ DFT, followed by (2) resampling to the desired $u$ grid using a 6-point, piecewise-quintic interpolant that we design expressly to minimize the ghosts, then (3) executing the DFT back to $x$ domain.Comment: Typographical and one algebraic correction, to appear in PASP March 201

Optimizing weak lensing mass estimates for cluster profile uncertainty

Weak lensing measurements of cluster masses are necessary for calibrating mass-observable relations (MORs) to investigate the growth of structure and the properties of dark energy. However, the measured cluster shear signal varies at fixed mass M_200m due to inherent ellipticity of background galaxies, intervening structures along the line of sight, and variations in the cluster structure due to scatter in concentrations, asphericity and substructure. We use N-body simulated halos to derive and evaluate a weak lensing circular aperture mass measurement M_ap that minimizes the mass estimate variance <(M_ap - M_200m)^2> in the presence of all these forms of variability. Depending on halo mass and observational conditions, the resulting mass estimator improves on M_ap filters optimized for circular NFW-profile clusters in the presence of uncorrelated large scale structure (LSS) about as much as the latter improve on an estimator that only minimizes the influence of shape noise. Optimizing for uncorrelated LSS while ignoring the variation of internal cluster structure puts too much weight on the profile near the cores of halos, and under some circumstances can even be worse than not accounting for LSS at all. We briefly discuss the impact of variability in cluster structure and correlated structures on the design and performance of weak lensing surveys intended to calibrate cluster MORs.Comment: 11 pages, 5 figures; accepted by MNRA

Characterization and correction of charge-induced pixel shifts in DECam

Interaction of charges in CCDs with the already accumulated charge distribution causes both a flux dependence of the point-spread function (an increase of observed size with flux, also known as the brighter/fatter effect) and pixel-to-pixel correlations of the Poissonian noise in flat fields. We describe these effects in the Dark Energy Camera (DECam) with charge dependent shifts of effective pixel borders, i.e. the Antilogus et al. (2014) model, which we fit to measurements of flat-field Poissonian noise correlations. The latter fall off approximately as a power-law r^-2.5 with pixel separation r, are isotropic except for an asymmetry in the direct neighbors along rows and columns, are stable in time, and are weakly dependent on wavelength. They show variations from chip to chip at the 20% level that correlate with the silicon resistivity. The charge shifts predicted by the model cause biased shape measurements, primarily due to their effect on bright stars, at levels exceeding weak lensing science requirements. We measure the flux dependence of star images and show that the effect can be mitigated by applying the reverse charge shifts at the pixel level during image processing. Differences in stellar size, however, remain significant due to residuals at larger distance from the centroid.Comment: typo and formatting fixes, matches version published in JINS

Phase clustering in complex networks of delay-coupled oscillatos

We study the clusterization of phase oscillators coupled with delay in complex networks. For the case of difussive oscillators, we formulate the equations relating the topology of the network and the phases and frequencies of the oscillators (functional response). We solve them exactly in directed networks for the case of perfect synchronization. We also compare the reliability of the solution of the linear system for non-linear couplings. Taking advantage of the form of the solution, we propose a frequency adaptation rule to achieve perfect synchronization. We also propose a mean-field theory for uncorrelated random networks that proves to be pretty accurate to predict phase synchronization in real topologies, as for example the C.elegans or the Autonomous Systems connectivity.Peer reviewe

Synchronous Rotation in the (136199) Eris–Dysnomia System

We combine photometry of Eris from a 6 month campaign on the Palomar 60 inch telescope in 2015, a 1 month Hubble Space Telescope WFC3 campaign in 2018, and Dark Energy Survey data spanning 2013–2018 to determine a light curve of definitive period 15.771 ± 0.008 days (1σ formal uncertainties), with nearly sinusoidal shape and peak-to-peak flux variation of 3%. This is consistent at part-per-thousand precision with the P = 15.785 90 ± 0.00005 day sidereal period of Dysnomia's orbit around Eris, strengthening the recent detection of synchronous rotation of Eris by Szakáts et al. with independent data. Photometry from Gaia are consistent with the same light curve. We detect a slope of 0.05 ± 0.01 mag per degree of Eris's brightness with respect to illumination phase averaged across g, V, and r bands, intermediate between Pluto's and Charon's values. Variations of 0.3 mag are detected in Dysnomia's brightness, plausibly consistent with a double-peaked light curve at the synchronous period. The synchronous rotation of Eris is consistent with simple tidal models initiated with a giant-impact origin of the binary, but is difficult to reconcile with gravitational capture of Dysnomia by Eris. The high albedo contrast between Eris and Dysnomia remains unexplained in the giant-impact scenario

The Phoenix stream : a cold stream in the southern hemisphere

We report the discovery of a stellar stream in the Dark Energy Survey Year 1 (Y1A1) data. The discovery was made through simple color–magnitude filters and visual inspection of the Y1A1 data. We refer to this new object as the Phoenix stream, after its resident constellation. After subtraction of the background stellar population we detect a clear signal of a simple stellar population. By fitting the ridge line of the stream in color–magnitude space, we find that a stellar population with age τ=11.5±0.5 Gyr and [Fe/H]<−1.6, located 17.5±0.9 kpc from the Sun, gives an adequate description of the stream stellar population. The stream is detected over an extension of 8°.1 (2.5 kpc) and has a width of ∼54 pc assuming a Gaussian profile, indicating that a globular cluster (GC) is a probable progenitor. There is no known GC within 5 kpc that is compatible with being the progenitor of the stream, assuming that the stream traces its orbit. We examined overdensities (ODs) along the stream, however, no obvious counterpart-bound stellar system is visible in the coadded images. We also find ODs along the stream that appear to be symmetrically distributed—consistent with the epicyclic OD scenario for the formation of cold streams—as well as a misalignment between the northern and southern part of stream. Despite the close proximity we find no evidence that this stream and the halo cluster NGC 1261 have a common accretion origin linked to the recently found EriPhe OD

Trans-Neptunian objects found in the first four years of the Dark Energy Survey

We present a catalog of 316 trans-Neptunian bodies (TNOs) detected from the first four seasons ("Y4" data) of the Dark Energy Survey (DES). The survey covers a contiguous 5000 deg(2) of the southern sky in the grizY optical/NIR filter set, with a typical TNO in this part of the sky being targeted by 25-30 Y4 exposures. This paper focuses on the methods used to detect these objects from the 60,000 Y4 exposures, a process made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. Newly developed techniques include: transient/moving object detection by comparison of single-epoch catalogs to catalogs of "stacked" images; quantified astrometric error from atmospheric turbulence; new software for detecting TNO linkages in a temporally sparse transient catalog, and for estimating the rate of spurious linkages; use of faint stars to determine the detection efficiency versus magnitude in all exposures. Final validation of the reality of linked orbits uses a new "sub-threshold confirmation" test, wherein we demand the object be detectable in a stack of the exposures in which the orbit indicates an object should be present, but was not individually detected. This catalog contains all validated TNOs which were detected on >= 6 unique nights in the Y4 data, and is complete to r less than or similar to 23.3 mag with virtually no dependence on orbital properties for bound TNOs at distance 30 au d 0.3 mag more depth, and arcs of >4 yr for nearly all detections.This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Searching for dark matter annihilation in recently discovered Milky Way satellites with Fermi-LAT

We search for excess γ-ray emission coincident with the positions of confirmed and candidate Milky Way satellite galaxies using six years of data from the Fermi Large Area Telescope (LAT). Our sample of 45 stellar systems includes 28 kinematically confirmed dark-matter-dominated dwarf spheroidal galaxies (dSphs) and 17 recently discovered systems that have photometric characteristics consistent with the population of known dSphs. For each of these targets, the relative predicted γ-ray flux due to dark matter annihilation is taken from kinematic analysis if available, and estimated from a distance-based scaling relation otherwise, assuming that the stellar systems are DM-dominated dSphs. LAT data coincident with four of the newly discovered targets show a slight preference (each ~2σ local) for γ-ray emission in excess of the background. However, the ensemble of derived γ-ray flux upper limits for individual targets is consistent with the expectation from analyzing random blank-sky regions, and a combined analysis of the population of stellar systems yields no globally significant excess (global significance 1 TeV and mDM,t+t-> 70 GeV) and weakening by a factor of ~1.5 at lower masses relative to previously observed limits

Redshift distributions of galaxies in the Dark Energy Survey Science Verification shear catalogue and implications for weak lensing

We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based photometric redshift methods—ANNZ2, BPZ calibrated against BCC-Ufig simulations, SKYNET, and TPZ—are analyzed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-z’s. From the galaxies in the DES SV shear catalogue, which have mean redshift 0.72 0.01 over the range 0.3 < z < 1.3, we construct three tomographic bins with means of z ¼ f0.45; 0.67; 1.00g. These bins each have systematic uncertainties δz ≲ 0.05 in the mean of the fiducial SKYNET photo-z nðzÞ. We propagate the errors in the redshift distributions through to their impact on cosmological parameters estimated with cosmic shear, and find that they cause shifts in the value of σ8 of approximately 3%. This shift is within the one sigma statistical errors on σ8 for the DES SV shear catalogue. We further study the potential impact of systematic differences on the critical surface density, Σcrit, finding levels of bias safely less than the statistical power of DES SV data. We recommend a final Gaussian prior for the photo-z bias in the mean of nðzÞ of width 0.05 for each of the three tomographic bins, and show that this is a sufficient bias model for the corresponding cosmology analysis

The Dark Energy Survey : more than dark energy – an overview

This overview paper describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4 m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion, the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterize dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large-scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper, we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from ‘Science Verification’, and from the first, second and third seasons of observations), what DES can tell us about the Solar system, the Milky Way, galaxy evolution, quasars and other topics. In addition, we show that if the cosmological model is assumed to be +cold dark matter, then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 trans-Neptunian objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed)