• Methyl cyanide is an important trace molecule in space, especially in star-forming regions where it is one of the more common molecules used to derive kinetic temperatures. We want to obtain accurate spectroscopic parameters of minor isotopologs of methyl cyanide in their lowest excited $v_8 = 1$ vibrational states to support astronomical observations, in particular, with interferometers such as ALMA. The laboratory rotational spectrum of methyl cyanide in natural isotopic composition has been recorded from the millimeter to the terahertz regions. Transitions with good signal-to-noise ratios could be identified for the three isotopic species CH$_3^{13}$CN, $^{13}$CH$_3$CN, and CH$_3$C(15)N up to about 1.2 THz ($J'' \le 66$). Accurate spectroscopic parameters were obtained for all three species. The present data were already instrumental in identifying $v_8 = 1$ lines of methyl cyanide with one $^{13}$C in IRAM 30 m and ALMA data toward Sagittarius B2(N).
  • Rotational transitions of $iso$-propyl cyanide, (CH$_3$)$_2$CHCN, also known as $iso$-butyronitrile, were recorded using long-path absorption spectroscopy in selected regions between 37 and 600 GHz. Further measurements were carried out between 6 and 20 GHz employing Fourier transform microwave (FTMW) spectroscopy on a pulsed molecular supersonic jet. The observed transitions reach $J$ and $K_a$ quantum numbers of 103 and 59, respectively, and yield accurate rotational constants as well as distortion parameters up to eighth order. The $^{14}$N nuclear hyperfine splitting was resolved in particular by FTMW spectroscopy yielding spin-rotation parameters as well as very accurate quadrupole coupling terms. In addition, Stark effect measurements were carried out in the microwave region to obtain a largely revised $c$-dipole moment component and to improve the $a$-component. The hyperfine coupling and dipole moment values are compared with values for related molecules both from experiment and from quantum chemical calculations.
  • Methyl cyanide is an important trace molecule in star-forming regions. It is one of the more common molecules used to derive kinetic temperatures in such sources. As preparatory work for Herschel, SOFIA, and in particular ALMA we want to improve the rest frequencies of the main as well as minor isotopologs of methyl cyanide. The laboratory rotational spectrum of methyl cyanide in natural isotopic composition has been recorded up to 1.63 THz. Transitions with good signal-to-noise ratio could be identified for CH$_3$CN, $^{13}$CH$_3$CN, CH$_3^{13}$CN, CH$_3$C$^{15}$N, CH$_2$DCN, and $^{13}$CH$_3^{13}$CN in their ground vibrational states up to about 1.2 THz. The main isotopic species could be identified even in the highest frequency spectral recordings around 1.6 THz. The highest $J'$ quantum numbers included in the fit are 64 for $^{13}$CH$_3^{13}$CN and 89 for the main isotopic species. Greatly improved spectroscopic parameters have been obtained by fitting the present data together with previously reported transition frequencies. The present data will be helpful to identify isotopologs of methyl cyanide in the higher frequency bands of instruments such as the recently launched Herschel satellite, the upcoming airplane mission SOFIA or the radio telescope array ALMA.
  • We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the HLLD approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibited a convergent picture of the early nonlinear growth of the MRI in a global model, and indicated that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.
  • The rotational spectrum of the formaldehyde isotopologue H$_2$C$^{17}$O was investigated between 0.56 and 1.50 THz using a sample of natural isotopic composition. In addition, transition frequencies were determined for H$_2$C$^{18}$O and H$_2$C$^{16}$O between 1.37 and 1.50 THz. The data were combined with critically evaluated literature data to derive improved sets of spectroscopic parameters which include $^{17}$O or H nuclear hyperfine structure parameters.
  • We develop a time-ordered data simulator and map-maker for the proposed PIXIE Fourier transform spectrometer and use them to investigate the impact of polarization leakage, imperfect collimation, elliptical beams, sub-pixel effects, correlated noise and spectrometer mirror jitter on the PIXIE data analysis. We find that PIXIE is robust to all of these effects, with the exception of mirror jitter which could become the dominant source of noise in the experiment if the jitter is not kept significantly below $0.1\mu m\sqrt{s}$. Source code is available at https://github.com/amaurea/pixie.
  • K. C. Chambers, E. A. Magnier, N. Metcalfe, H. A. Flewelling, M. E. Huber, C. Z. Waters, L. Denneau, P. W. Draper, D. Farrow, D. P. Finkbeiner, C. Holmberg, J. Koppenhoefer, P. A. Price, A. Rest, R. P. Saglia, E. F. Schlafly, S. J. Smartt, W. Sweeney, R. J. Wainscoat, W. S. Burgett, S. Chastel, T. Grav, J. N. Heasley, K. W. Hodapp, R. Jedicke, N. Kaiser, R.-P. Kudritzki, G. A. Luppino, R. H. Lupton, D. G. Monet, J. S. Morgan, P. M. Onaka, B. Shiao, C. W. Stubbs, J. L. Tonry, R. White, E. Bañados, E. F. Bell, R. Bender, E. J. Bernard, M. Boegner, F. Boffi, M. T. Botticella, A. Calamida, S. Casertano, W.-P. Chen, X. Chen, S. Cole, N. Deacon, C. Frenk, A. Fitzsimmons, S. Gezari, V. Gibbs, C. Goessl, T. Goggia, R. Gourgue, B. Goldman, P. Grant, E. K. Grebel, N.C. Hambly, G. Hasinger, A. F. Heavens, T. M. Heckman, R. Henderson, T. Henning, M. Holman, U. Hopp, W.-H. Ip, S. Isani, M. Jackson, C. D. Keyes, A. M. Koekemoer, R. Kotak, D. Le, D. Liska, K. S. Long, J. R. Lucey, M. Liu, N. F. Martin, G. Masci, B. McLean, E. Mindel, P. Misra, E. Morganson, D. N. A. Murphy, A. Obaika, G. Narayan, M. A. Nieto-Santisteban, P. Norberg, J. A. Peacock, E. A. Pier, M. Postman, N. Primak, C. Rae, A. Rai, A. Riess, A. Riffeser, H. W. Rix, S. Röser, R. Russel, L. Rutz, E. Schilbach, A. S. B. Schultz, D. Scolnic, L. Strolger, A. Szalay, S. Seitz, E. Small, K. W. Smith, D. R. Soderblom, P. Taylor, R. Thomson, A. N. Taylor, A. R. Thakar, J. Thiel, D. Thilker, D. Unger, Y. Urata, J. Valenti, J. Wagner, T. Walder, F. Walter, S. P. Watters, S. Werner, W. M. Wood-Vasey, R. Wyse
    Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.
  • This paper describes the organization of the database and the catalog data products from the Pan-STARRS1 $3\pi$ Steradian Survey. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI. The database is described in detail, including the construction of the database, the provenance of the data, the schema, and how the database tables are related. Examples of queries for a range of science goals are included. The catalog data products are available in the form of an SQL-based relational database from MAST, the Mikulski Archive for Space Telescopes at STScI.
  • The Pan-STARRS1 Science Consortium have carried out a set of imaging surveys using the 1.4 giga-pixel GPC1 camera on the PS1 telescope. As this camera is composed of many individual electronic readouts, and covers a very large field of view, great care was taken to ensure that the many instrumental effects were corrected to produce the most uniform detector response possible. We present the image detrending steps used as part of the processing of the data contained within the public release of the Pan-STARRS1 Data Release 1 (DR1). In addition to the single image processing, the methods used to transform the 375,573 individual exposures into a common sky-oriented grid are discussed, as well as those used to produce both the image stack and difference combination products.
  • We present the details of the photometric and astrometric calibration of the Pan-STARRS1 $3\pi$ Survey. The photometric goals were to reduce the systematic effects introduced by the camera and detectors, and to place all of the observations onto a photometric system with consistent zero points over the entire area surveyed, the ~30,000 square degrees north of $\delta$ = -30 degrees. The astrometric calibration compensates for similar systematic effects so that positions, proper motions, and parallaxes are reliable as well. The Pan-STARRS Data Release 2 (DR2) astrometry is tied to the Gaia DR1 release.
  • The Pan-STARRS Data Processing System is responsible for the steps needed to downloaded, archive, and process all images obtained by the Pan-STARRS telescopes, including real-time detection of transient sources such as supernovae and moving objects including potentially hazardous asteroids. With a nightly data volume of up to 4 terabytes and an archive of over 4 petabytes of raw imagery, Pan-STARRS is solidly in the realm of Big Data astronomy. The full data processing system consists of several subsystems covering the wide range of necessary capabilities. This article describes the Image Processing Pipeline and its connections to both the summit data systems and the outward-facing systems downstream. The latter include the Moving Object Processing System (MOPS) & the public database: the Published Science Products Subsystem (PSPS).
  • Over 3 billion astronomical objects have been detected in the more than 22 million orthogonal transfer CCD images obtained as part of the Pan-STARRS1 $3\pi$ survey. Over 85 billion instances of those objects have been automatically detected and characterized by the Pan-STARRS Image Processing Pipeline photometry software, psphot. This fast, automatic, and reliable software was developed for the Pan-STARRS project, but is easily adaptable to images from other telescopes. We describe the analysis of the astronomical objects by psphot in general as well as for the specific case of the 3rd processing version used for the first two public releases of the Pan-STARRS $3\pi$ survey data, DR1 & DR2.
  • (Abr.) Laser guide stars employed at astronomical observatories provide artificial wavefront reference sources to help correct (in part) the impact of atmospheric turbulence on astrophysical observations. Following the recent commissioning of the 4 Laser Guide Star Facility (4LGSF) on UT4 at the VLT, we characterize the spectral signature of the uplink beams from the 22W lasers to assess the impact of laser scattering from the 4LGSF on science observations. We use the MUSE optical integral field spectrograph to acquire spectra at a resolution of R~3000 of the uplink laser beams over the wavelength range of 4750\AA\ to 9350\AA. We report the first detection of laser-induced Raman scattering by N2, O2, CO2, H2O and (tentatively) CH4 molecules in the atmosphere above the astronomical observatory of Cerro Paranal. In particular, our observations reveal the characteristic spectral signature of laser photons -- but 480\AA\ to 2210\AA\ redder than the original laser wavelength of 5889.959\AA\ -- landing on the 8.2m primary mirror of UT4 after being Raman-scattered on their way up to the sodium layer. Laser-induced Raman scattering is not unique to the observatory of Cerro Paranal, but common to any astronomical telescope employing a laser-guide-star (LGS) system. It is thus essential for any optical spectrograph coupled to a LGS system to handle thoroughly the possibility of a Raman spectral contamination via a proper baffling of the instrument and suitable calibrations procedures. These considerations are particularly applicable for the HARMONI optical spectrograph on the upcoming Extremely Large Telescope. At sites hosting multiple telescopes, laser collision prediction tools also ought to account for the presence of Raman emission from the uplink laser beam(s) to avoid the unintentional contamination of observations acquired with telescopes in the vicinity of a LGS system.
  • Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.
  • Meaningful laws of nature must be independent of the units employed to measure the variables. The principle of similitude (Rayleigh 1915) or dimensional homogeneity, states that only commensurable quantities (ones having the same dimension) may be compared, therefore, meaningful laws of nature must be homogeneous equations in their various units of measurement, a result which was formalized in the $\rm \Pi$ theorem (Vaschy 1892; Buckingham 1914). However, most relations in allometry do not satisfy this basic requirement, including the `3/4 Law' (Kleiber 1932) that relates the basal metabolic rate and body mass, which it is sometimes claimed to be the most fundamental biological rate (Brown et al. 2004) and the closest to a law in life sciences (West \& Brown 2004). Using the $\rm \Pi$ theorem, here we show that it is possible to construct a unique homogeneous equation for the metabolic rates, in agreement with data in the literature. We find that the variations in the dependence of the metabolic rates on body mass are secondary, coming from variations in the allometric dependence of the heart frequencies. This includes not only different classes of animals (mammals, birds, invertebrates) but also different exercise conditions (basal and maximal). Our results demonstrate that most of the differences found in the allometric exponents (White et al. 2007) are due to compare incommensurable quantities and that our dimensionally homogenous formula, unify these differences into a single formulation. We discuss the ecological implications of this new formulation in the context of the Malthusian's, Fenchel's and the total energy consumed in a lifespan relations.
  • Machine learning techniques offer a precious tool box for use within astronomy to solve problems involving so-called big data. They provide a means to make accurate predictions about a particular system without prior knowledge of the underlying physical processes of the data. In this article, and the companion papers of this series, we present the set of Generalized Linear Models (GLMs) as a fast alternative method for tackling general astronomical problems, including the ones related to the machine learning paradigm. To demonstrate the applicability of GLMs to inherently positive and continuous physical observables, we explore their use in estimating the photometric redshifts of galaxies from their multi-wavelength photometry. Using the gamma family with a log link function we predict redshifts from the PHoto-z Accuracy Testing simulated catalogue and a subset of the Sloan Digital Sky Survey from Data Release 10. We obtain fits that result in catastrophic outlier rates as low as ~1% for simulated and ~2% for real data. Moreover, we can easily obtain such levels of precision within a matter of seconds on a normal desktop computer and with training sets that contain merely thousands of galaxies. Our software is made publicly available as an user-friendly package developed in Python, R and via an interactive web application (https://cosmostatisticsinitiative.shinyapps.io/CosmoPhotoz). This software allows users to apply a set of GLMs to their own photometric catalogues and generates publication quality plots with minimum effort from the user. By facilitating their ease of use to the astronomical community, this paper series aims to make GLMs widely known and to encourage their implementation in future large-scale projects, such as the Large Synoptic Survey Telescope.
  • Inferring model parameters from experimental data is a grand challenge in many sciences, including cosmology. This often relies critically on high fidelity numerical simulations, which are prohibitively computationally expensive. The application of deep learning techniques to generative modeling is renewing interest in using high dimensional density estimators as computationally inexpensive emulators of fully-fledged simulations. These generative models have the potential to make a dramatic shift in the field of scientific simulations, but for that shift to happen we need to study the performance of such generators in the precision regime needed for science applications. To this end, in this work we apply Generative Adversarial Networks to the problem of generating weak lensing convergence maps. We show that our generator network produces maps that are described by, with high statistical confidence, the same summary statistics as the fully simulated maps.
  • During its first observing run, in late 2015, the advanced LIGO facility announced 3 gravitational wave (GW) triggers to electromagnetic follow-up partners. Two of these have since been confirmed as being of astrophysical origin: both are binary black hole mergers at ~500 Mpc; the other trigger was later found not to be astrophysical. In this paper we report on the Swift follow up observations of the second and third triggers, including details of 21 X-ray sources detected; none of which can be associated with the GW event. We also consider the challenges that the next GW observing run will bring as the sensitivity and hence typical distance of GW events will increase. We discuss how to effectively use galaxy catalogues to prioritise areas for follow up, especially in the presence of distance estimates from the GW data. We also consider two galaxy catalogues and suggest that the high completeness at larger distances of the 2MASS Photometric Redshift Catalogue (2MPZ) makes it very well suited to optimise Swift follow-up observations.
  • This is a follow-up sensitivity study on r-mode gravitational wave signals from newborn neutron stars illustrating the applicability of machine learning algorithms for the detection of long-lived gravitational-wave transients. In this sensitivity study we examine three machine learning algorithms (MLAs): artificial neural networks (ANNs), support vector machines (SVMs) and constrained subspace classifiers (CSCs). The objective of this study is to compare the detection efficiency that MLAs can achieve with the efficiency of conventional detection algorithms discussed in an earlier paper. Comparisons are made using 2 distinct r-mode waveforms. For the training of the MLAs we assumed that some information about the distance to the source is given so that the training was performed over distance ranges not wider than half an order of magnitude. The results of this study suggest that machine learning algorithms are suitable for the detection of long-lived gravitational-wave transients and that when assuming knowledge of the distance to the source, MLAs are at least as efficient as conventional methods.
  • ANTARES Collaboration: A. Albert, M. André, M. Anghinolfi, G. Anton, M. Ardid, J.-J. Aubert, T. Avgitas, B. Baret, J. Barrios-Martí, S. Basa, B. Belhorma, V. Bertin, S. Biagi, R. Bormuth, S. Bourret, M.C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, A. Capone, L. Caramete, J. Carr, S. Celli, R. Cherkaoui El Moursli, T. Chiarusi, M. Circella, J.A.B. Coelho, A. Coleiro, R. Coniglione, H. Costantini, P. Coyle, A. Creusot, A. F. Díaz, A. Deschamps, G. De Bonis, C. Distefano, I. Di Palma, A. Domi, C. Donzaud, D. Dornic, D. Drouhin, T. Eberl, I. El Bojaddaini, N. El Khayati, D. Elsässer, A. Enzenhöfer, A. Ettahiri, F. Fassi, I. Felis, L.A. Fusco, S. Galatà, P. Gay, V. Giordano, H. Glotin, T. Grégoire, R. Gracia Ruiz, K. Graf, S. Hallmann, H. van Haren, A.J. Heijboer, Y. Hello, J.J. Hernández-Rey, J. Hößl, J. Hofestädt, C. Hugon, G. Illuminati, C.W. James, M. de Jong, M. Jongen, M. Kadler, O. Kalekin, U. Katz, D. Kießling, A. Kouchner, M. Kreter, I. Kreykenbohm, V. Kulikovskiy, C. Lachaud, R. Lahmann, D. Lefèvre, E. Leonora, M. Lotze, S. Loucatos, M. Marcelin, A. Margiotta, A. Marinelli, J.A. Martínez-Mora, R. Mele, K. Melis, T. Michael, P. Migliozzi, A. Moussa, S. Navas, E. Nezri, M. Organokov, G.E. Păvălaş, C. Pellegrino, C. Perrina, P. Piattelli, V. Popa, T. Pradier, L. Quinn, C. Racca, G. Riccobene, A. Sánchez-Losa, M. Saldaña, I. Salvadori, D. F. E. Samtleben, M. Sanguineti, P. Sapienza, F. Schüssler, C. Sieger, M. Spurio, Th. Stolarczyk, M. Taiuti, Y. Tayalati, A. Trovato, D. Turpin, C. Tönnis, B. Vallage, V. Van Elewyck, F. Versari, D. Vivolo, A. Vizzoca, J. Wilms, J.D. Zornoza, J. Zúñiga
    Dec. 5, 2018 hep-ex, astro-ph.IM, astro-ph.HE
    A search for cosmic neutrino sources using the data collected with the ANTARES neutrino telescope between early 2007 and the end of 2015 is performed. For the first time, all neutrino interactions --charged and neutral current interactions of all flavours-- are considered in a search for point-like sources with the ANTARES detector. In previous analyses, only muon neutrino charged current interactions were used. This is achieved by using a novel reconstruction algorithm for shower-like events in addition to the standard muon track reconstruction. The shower channel contributes about 23\% of all signal events for an $E^{-2}$ energy spectrum. No significant excess over background is found. The most signal-like cluster of events is located at $(\alpha,\delta) = (343.8^\circ, 23.5^\circ)$ with a significance of $1.9\sigma$. The neutrino flux sensitivity of the search is about $E^2 d\varPhi/dE = 6\cdot10^{-9} GeV cm^{-2} s^{-1}$ for declinations from $-90^\circ$ up to $-42^\circ$, and below $10^{-8} GeV cm^{-2} s^{-1}$ for declinations up to $5^{\circ}$. The directions of 106 source candidates and of 13 muon track events from the IceCube HESE sample are investigated for a possible neutrino signal and upper limits on the signal flux are determined.
  • We have carried a detailed analysis on the impact of cosmological redshift in the non-parametric approach to automated galaxy morphology classification. We artificially redshifted each galaxy from the EFIGI 4458 sample (re-centred at $z \sim 0$) simulating SDSS, DES, LSST and HST instruments set-ups over the range $0 < z < 1.5$. We then traced how the morphometry is degraded in each $z$ using MORFOMETRYKA. In the process, we re-sampled all catalogues to several resolutions and to a diverse signal-to-noise range, allowing us to understand the impact of image sampling and noise on our measurements separately. We summarize by exploring the impact of these effects on our capacity to perform automated galaxy supervised morphological classification by investigating the degradation of our classifier's metrics as a function of redshift for each instrument. The overall conclusion is that we can make reliable classification with MORFOMETRYKA for $z < 0.2$ with SDSS, for $z < 0.5$ with DES, for $z < 0.8$ with LSST and for at least $z < 1.5$ with HST.
  • We present in this paper the general formalism and data processing steps used in the MATISSE data reduction software, as it has been developed by the MATISSE consortium. The MATISSE instrument is the mid-infrared new generation interferometric instrument of the Very Large Telescope Interferometer (VLTI). It is a 2-in-1 instrument with 2 cryostats and 2 detectors: one 2k x 2k Rockwell Hawaii 2RG detector for L\&M-bands, and one 1k x 1k Raytheon Aquarius detector for N-band, both read at high framerates, up to 30 frames per second. MATISSE is undergoing its first tests in laboratory today.
  • The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or even years to complete, while an orbit-sampling Monte Carlo approach cannot adequately model the dynamics in a dense cluster core, particularly in the presence of many black holes. We have developed a new technique combining the precision of a direct N-body integration with the speed of a Monte Carlo approach. Our Rapid And Precisely Integrated Dynamics code, the RAPID code, statistically models interactions between neighboring stars and stellar binaries while integrating directly the orbits of stars or black holes in the cluster core. This allows us to accurately simulate the dynamics of the black holes in a realistic globular cluster environment without the burdensome N^2 scaling of a full N-body integration. We compare RAPID models of idealized globular clusters to identical models from the direct N-body and Monte Carlo methods. Our tests show that RAPID can reproduce the half-mass radii, core radii, black hole ejection rates, and binary properties of the direct N-body models far more accurately than a standard Monte Carlo integration while remaining significantly faster than a full N-body integration. With this technique, it will be possible to create more realistic models of Milky Way globular clusters with sufficient rapidity to explore the full parameter space of dense stellar clusters.
  • We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting 'preferred directions' on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common 'fully-thermal' (energy-dump) or 'fully-kinetic' (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution >100 solar masses, they diverge by orders-of-magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution. However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.
  • The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.