The Fujitsu Digital Annealer is designed to solve fully connected quadratic unconstrained binary optimization (QUBO) problems. It is implemented on application-specific CMOS hardware and currently ...solves problems of up to 1,024 variables. The Digital Annealer's algorithm is currently based on simulated annealing; however, it differs from it in its utilization of an efficient parallel-trial scheme and a dynamic escape mechanism. In addition, the Digital Annealer exploits the massive parallelization that custom application-specific CMOS hardware allows. We compare the performance of the Digital Annealer to simulated annealing and parallel tempering with isoenergetic cluster moves on two-dimensional and fully connected spin-glass problems with bimodal and Gaussian couplings. These represent the respective limits of sparse vs. dense problems, as well as high-degeneracy vs. low-degeneracy problems. Our results show that the Digital Annealer currently exhibits a time-to-solution speedup of roughly two orders of magnitude for fully connected spin-glass problems with bimodal or Gaussian couplings, over the single-core implementations of simulated annealing and parallel tempering Monte Carlo used in this study. The Digital Annealer does not appear to exhibit a speedup for sparse two-dimensional spin-glass problems, which we explain on theoretical grounds. We also benchmarked an early implementation of the Parallel Tempering Digital Annealer. Our results suggest an improved scaling over the other algorithms for fully connected problems of average difficulty with bimodal disorder. The next generation of the Digital Annealer is expected to be able to solve fully connected problems up to 8,192 variables in size. This would enable the study of fundamental physics problems and industrial applications that were previously inaccessible using standard computing hardware or special-purpose quantum annealing machines.
We present the largest submillimeter images that have been made of the extragalactic sky. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 660 deg2 with the PACS and ...SPIRE cameras in five photometric bands: 100, 160, 250, 350, and 500 m. In this paper we present the images from our two largest fields, which account for ∼75% of the survey. The first field is 180.1 deg2 in size, centered on the north Galactic pole (NGP), and the second is 317.6 deg2 in size, centered on the south Galactic pole. The NGP field serendipitously contains the Coma cluster. Over most (∼80%) of the images, the pixel noise, including both instrumental noise and confusion noise, is approximately 3.6, and 3.5 mJy pix−1 at 100 and 160 m, and 11.0, 11.1 and 12.3 mJy beam−1 at 250, 350 and 500 m, respectively, but reaches lower values in some parts of the images. If a matched filter is applied to optimize point-source detection, our total 1 map sensitivity is 5.7, 6.0, and 7.3 mJy at 250, 350, and 500 m, respectively. We describe the results of an investigation of the noise properties of the images. We make the most precise estimate of confusion in SPIRE maps to date, finding values of 3.12 0.07, 4.13 0.02, and 4.45 0.04 mJy beam−1 at 250, 350, and 500 m in our un-convolved maps. For PACS we find an estimate of the confusion noise in our fast-parallel observations of 4.23 and 4.62 mJy beam−1 at 100 and 160 m. Finally, we give recipes for using these images to carry out photometry, both for unresolved and extended sources.
Abstract
The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys ...reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a single Galaxy Sequence (GS) rather than a star-forming ‘main sequence’ and a separate region of ‘passive’ or ‘red-and-dead’ galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations – they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population.
Submillimeter galaxies (SMGs) at are luminous in the far-infrared, and have star formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true ...analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at -4. We present 70-160 m photometry and spectroscopy of the O iv26 m, Fe ii26 m, S iii33 m, Si ii34 m, O iii52 m, N iii57 m, and O i63 m fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected; we instead focus on stacking these spectra with observations of an additional 32 SMGs from the Herschel archive-representing a complete compilation of PACS spectroscopy of SMGs. We detect O i63 m, Si ii34 m, and N iii57 m at in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are , , and , respectively. Using the O iii52 m/N iii57 m emission line ratio, we show that SMGs have average gas-phase metallicities . By using PDR modeling and combining the new spectral measurements with integrated far-infrared fluxes and existing C ii158 m data, we show that SMGs have average gas densities, n, of and FUV field strengths, (in Habing units: ), consistent with both local ULIRGs and lower luminosity star-forming galaxies.
Star formation and accretion onto supermassive black holes in the nuclei of galaxies are the two most energetic processes in the universe, producing the bulk of the observed emission throughout its ...history. We simulated the luminosity functions of star-forming and active galaxies for spectral lines that are thought to be good spectroscopic tracers of either phenomenon, as a function of redshift. We focused on the infrared (IR) and submillimeter domains, where the effects of dust obscuration are minimal. Using three different and independent theoretical models for galaxy formation and evolution, constrained by multi-wavelength luminosity functions, we computed the number of star-forming and active galaxies per IR luminosity and redshift bin. We converted the continuum luminosity counts into spectral line counts using relationships that we calibrated on mid- and far-IR spectroscopic surveys of galaxies in the local universe. Our results demonstrate that future facilities optimized for survey-mode observations, i.e., the Space Infrared telescope for Cosmology and Astrophysics and the Cerro Chajnantor Atacama Telescope, will be able to observe thousands of z > 1 galaxies in key fine-structure lines, e.g., SiII, OI, OIII, CIII, in a half-square-degree survey, with 1 hr integration time per field of view. Fainter lines such as OIV, NeV, and H sub(2) (0-0)S1 will be observed in several tens of bright galaxies at 1 < z < 2, while diagnostic diagrams of active nucleus versus star formation activity will be feasible even for normal z ~ 1 galaxies. We discuss the new parameter space that these future telescopes will cover and that strongly motivates their construction.
Interstellar dust in galaxies can be traced either through its extinction effects on the star light or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of ...several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of 3. We investigate the dust energy balance for IC 4225 and NGC 5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC 4225 underestimates the observed fluxes by a factor of about 3. For NGC 5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.
The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty star-forming galaxies (DSFGs) with rising SPIRE flux densities (S500 > S350 > S250; the ...so-called "500 m risers") as an efficient way for identifying DSFGs at higher redshift (z > 4). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution Atacama Large Millimeter/submillimeter Array, Northern Extended Millimeter Array, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, ∼65% appear to be unlensed and ∼27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multiwavelength spectral energy distribution modeling to derive their physical properties and compare with the more numerous z ∼ 2 DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 × 1011 M , star formation rate (SFR) of 730 M yr−1, total dust luminosity of 9.0 × 1012 L , dust mass of 2.8 × 109 M , and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous, and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation.
Recently, there has been considerable interest in solving optimization problems by mapping these onto a binary representation, sparked mostly by the use of quantum annealing machines. Such binary ...representation is reminiscent of a discrete physical two-state system, such as the Ising model. As such, physics-inspired techniques—commonly used in fundamental physics studies—are ideally suited to solve optimization problems in a binary format. While binary representations can be often found for paradigmatic optimization problems, these typically result in k-local higher-order unconstrained binary optimization cost functions. In this work, we discuss the effects of locality reduction needed for the majority of the currently available quantum and quantum-inspired solvers that can only accommodate 2-local (quadratic) cost functions. General locality reduction approaches require the introduction of ancillary variables which cause an overhead over the native problem. Using a parallel tempering Monte Carlo solver on Microsoft Azure Quantum, as well as k-local binary problems with planted solutions, we show that post reduction to a corresponding 2-local representation the problems become considerably harder to solve. We further quantify the increase in computational hardness introduced by the reduction algorithm by measuring the variation of number of variables, statistics of the coefficient values, and the population annealing entropic family size. Our results demonstrate the importance of avoiding locality reduction when solving optimization problems.
We present results from a survey carried out by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) on a 9 deg2 field near the South Ecliptic Pole at 250, 350, and 500 Delta *mm. The ...median 1 Delta *s depths of the maps are 36.0, 26.4, and 18.4 mJy, respectively. We apply a statistical method to estimate submillimeter galaxy number counts and find that they are in agreement with other measurements made with the same instrument and with the more recent results from Herschel/SPIRE. Thanks to the large field observed, the new measurements give additional constraints on the bright end of the counts. We identify 132, 89, and 61 sources with S/N >=4 at 250, 350, 500 Delta *mm, respectively and provide a multi-wavelength combined catalog of 232 sources with a significance >=4 Delta *s in at least one BLAST band. The new BLAST maps and catalogs are available publicly at http://blastexperiment.info.