Showing 1–50 of 112 results for author: McMillan, S

Massive Star Cluster Formation II. Runaway Stars as Fossils of Sub-Cluster Mergers

Authors: Brooke Polak, Mordecai-Mark Mac Low, Ralf S. Klessen, Simon Portegies Zwart, Eric P. Andersson, Sabrina M. Appel, Claude Cournoyer-Cloutier, Simon C. O. Glover, Stephen L. W. McMillan

Abstract: Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway star… ▽ More Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the sub-cluster ejection scenario (SCES), where an infalling sub-cluster of stars is ejected out of the cluster by a slingshot interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of a young massive cluster forming from a $10^6\rm~M_\odot$ gas cloud using the Torch framework. This star cluster forms hierarchically through a sequence of sub-cluster mergers, determined by the initial turbulent, spherical initial conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming sub-clusters fall into the central potential, where they are ejected on a slingshot trajectory, forming groups of runaway stars that are distributed highly anisotropically. Runaways formed by the same SCES share similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling sub-clusters were tightened through ejection. This allows for another velocity kick via subsequent BSS, which is a promising avenue for producing hypervelocity stars unbound to the Galaxy. The SCES occurs when sub-cluster formation is resolved. We expect non-spherical initial gas distributions to increase runaway star numbers. The observation of groups of runaway stars formed via SCES thus reveals the assembly history of their natal clusters. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 11 pages, 10 figures. Submitted to A&A. Comments welcome!

Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars

Authors: Brooke Polak, Mordecai-Mark Mac Low, Ralf S. Klessen, Jia Wei Teh, Claude Cournoyer-Cloutier, Eric P. Andersson, Sabrina M. Appel, Aaron Tran, Sean C. Lewis, Maite J. C. Wilhelm, Simon Portegies Zwart, Simon C. O. Glover, Long Wang, Stephen L. W. McMillan

Abstract: The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg… ▽ More The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Submitted to A&A

Early Evolution and 3D Structure of Embedded Star Clusters

Authors: Claude Cournoyer-Cloutier, Alison Sills, William E. Harris, Sabrina M. Appel, Sean C. Lewis, Brooke Polak, Aaron Tran, Martijn J. C. Wilhelm, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Simon Portegies Zwart

Abstract: We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at… ▽ More We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at 0.0683 pc resolution for $\sim$ 2 Myr after the onset of star formation, with virial parameters $α$ = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial $α$ having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over $\sim$ 0.01 Myr timescales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses $\lesssim$ 1000 M$_{\odot}$ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 15 pages, 10 figures, to be published in MNRAS

Early-Forming Massive Stars Suppress Star Formation and Hierarchical Cluster Assembly

Authors: Sean C. Lewis, Stephen L. W. McMillan, Mordecai-Mark Mac Low, Claude Cournoyer-Cloutier, Brooke Polak, Martijn J. C. Wilhelm, Aaron Tran, Alison Sills, Simon Portegies Zwart, Ralf S. Klessen, Joshua E. Wall

Abstract: Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation… ▽ More Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suite \texttt{Torch}, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer, $N$-body dynamics, and stellar feedback to model four initially identical $10^4$ M$_\odot$ giant molecular clouds with a Gaussian density profile peaking at $521.5 \mbox{ cm}^{-3}$. Using the \texttt{Torch} software suite through the \texttt{AMUSE} framework we modify three of the models to ensure that the first star that forms is very massive (50, 70, 100 M$_\odot$). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star forming region. The star formation rate is suppressed, reducing the total mass of stars formed. Our fiducial control model without an early massive star has a larger star formation rate and total efficiency by up to a factor of three and a higher average star formation efficiency per free-fall time by up to a factor of seven. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster. △ Less

Submitted 28 February, 2023; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 17 pages, 7 figures. Published in ApJ

Journal ref: ApJ 944 211 (2023)

Compositional Convection in the Deep Interior of Uranus

Authors: Dustin J. Hill, Krista M. Soderlund, Stephen L. W. McMillan

Abstract: Uranus and Neptune share properties that are distinct from the other giant planets in the solar system, but they are also distinct from one another, particularly in their relative internal heat flux. Not only does Neptune emit about ten times the amount of heat that emitted by Uranus, the relative amount of emitted heat to the energy they absorb from the sun also differs greatly, being comparable… ▽ More Uranus and Neptune share properties that are distinct from the other giant planets in the solar system, but they are also distinct from one another, particularly in their relative internal heat flux. Not only does Neptune emit about ten times the amount of heat that emitted by Uranus, the relative amount of emitted heat to the energy they absorb from the sun also differs greatly, being comparable at Uranus and the largest of all giant planets at Neptune. As a result, it is questionable whether thermal convection occurs within the interior of Uranus. However, the presence of an intrinsic magnetic field implies that interior fluid motions must exist. Here, we consider compositional convection driven by the release of hydrogen associated with the formation of large organic networks or diamond precipitation in the deep interior. We test this hypotheses using a set of numerical rotating convection models where the convective driving is varied between thermal and compositional sources and is sufficiently vigorous to not be strongly constrained by rotation. In most cases, we find ice-giant-like zonal flows develop, with three bands characterized by a retrograde equatorial jet and prograde jets at higher latitudes. Large-scale circulation cells also develop and lead to heat and mass fluxes that tend to exhibit local maxima along the equatorial plane. This similarity between convective flows driven by thermal and compositional buoyancy therefore predict Uranus and Neptune to have similar interior dynamics despite Uranus' minimal internal heat flow and may thus explain why both ice giants have comparable magnetic fields. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 30 Pages, 15 figures

Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity

Authors: Simon F. Portegies Zwart, Tjarda C. N. Boekholt, Emiel Por, Adrian S. Hamers, Steve L. W. McMillan

Abstract: In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, bu… ▽ More In self-gravitating $N$-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with $N$ up to $N \sim 30$, but that for larger systems, the growth rate has a weaker scaling with $N$. For concentrated systems, however, the rate of exponential growth continues to scale with $N$. In relativistic self-gravitating systems, the rate of growth is almost independent of $N$. This effect, however, is only noticeable when the system's mean velocity approaches the speed of light to within three orders of magnitude. The chaotic behavior of systems with more than a dozen bodies for the usually adopted approximation of only solving the pairwise interactions in the Einstein-Infeld-Hoffmann equation of motion is qualitatively different than when the interaction terms (or cross terms) are taken into account. This result provides a strong motivation for follow-up studies on the microscopic effect of general relativity on orbital chaos, and on the influence of higher-order cross-terms in the Taylor-series expansion of the Einstein-Infeld-Hoffmann equations of motion. △ Less

Submitted 17 November, 2021; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: A&A in press

Journal ref: A&A 659, A86 (2022)

Trigonometric parallaxes of two T dwarfs with Keck and ShaneAO astrometry

Authors: Jayke Nguyen, S. Mark Ammons, Kaitlin Dennison, E. Victor Garcia, Jessica R. Lu, Stephen McMillan, Maissa Salama

Abstract: We present trigonometric parallax and proper motion measurements for two T-type brown dwarfs. We derive our measurements from infrared laser guide star adaptive optics observations spanning five years from the ShaneAO/SHARCS and NIRC2/medium-cam instruments on the Shane and Keck telescopes, respectively. To improve our astrometric precision, we measure and apply a distortion correction to our fiel… ▽ More We present trigonometric parallax and proper motion measurements for two T-type brown dwarfs. We derive our measurements from infrared laser guide star adaptive optics observations spanning five years from the ShaneAO/SHARCS and NIRC2/medium-cam instruments on the Shane and Keck telescopes, respectively. To improve our astrometric precision, we measure and apply a distortion correction to our fields for both instruments. We also transform the Keck and ShaneAO astrometric reference frames onto the ICRS using five-parameter parallax and proper motion solutions for background reference stars from Gaia DR2. Fitting for parallax and proper motion, we measure parallaxes of $73.5\pm9.2$ mas and $70.1\pm6.7$ mas for WISEJ19010703+47181688 (WISE1901) and WISEJ21543294+59421370 (WISE2154), respectively. We utilize Monte Carlo methods to estimate the error in our sparse field methods, taking into account overfitting and differential atmospheric refraction. Comparing to previous measurements in the literature, all of our parallax and proper motion values fall within $2σ$ of the published measurements, and 4 of 6 measurements are within $1σ$. These data are among the first parallax measurements of these T dwarfs and serve as precise measurements for calibrating stellar formation models. These two objects are the first results of an ongoing survey of T dwarfs with Keck/NIRC2 and the Shane Adaptive Optics system at Lick Observatory. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 18 Pages, 19 Figures. Accepted by PASP

Report number: LLNL-JRNL-821003

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

Authors: Claude Cournoyer-Cloutier, Aaron Tran, Sean Lewis, Joshua E. Wall, William E. Harris, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Simon Portegies Zwart, Alison Sills

Abstract: The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled M… ▽ More The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 15 pages, 14 figures, submitted to MNRAS and edited to address positive referee's report

Modelling of the Effects of Stellar Feedback during Star Cluster Formation Using a Hybrid Gas and N-Body Method

Authors: Joshua E. Wall, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Ralf S. Klessen, Simon Portegies Zwart, Andrew Pellegrino

Abstract: Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of… ▽ More Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds and supernovae in FLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ~1000 stars formed from clouds with a mass range of 10^3-10^5 M_sun. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H II region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations. △ Less

Submitted 10 October, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted to ApJ. 29 pages, 18 figures. Source code at https://bitbucket.org/torch-sf/torch/ and documentation at https://torch-sf.bitbucket.io/

TYCHO -- Realistically Simulating Exoplanets within Stellar Clusters I: Improving the Monte Carlo Approach

Authors: Joseph P. Glaser, Stephen L. W. McMillan, Aaron M. Geller, Jonathan D. Thornton, Mark R. Giovinazzi

Abstract: To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses $>0.5M_\odot$ (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additi… ▽ More To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses $>0.5M_\odot$ (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale $N$-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, TYCHO, which employs a variety of different computational techniques, including multiple $N$-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1$\%$ of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which perturb long-period planets lead to 38.3$\%$ of solar-system-like planetary systems becoming long-term unstable. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: 16 Pages, 9 Figures, 4 Tables; Submitted to ApJ

Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics

Authors: Simon Portegies Zwart, Inti Pelupessy, Carmen Martinez-Barbosa, Arjen van Elteren, Steve McMillan

Abstract: Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method… ▽ More Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations of Bridge include higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processes. Such augmentation allows for including additional processes in a classic Newtonian N-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset. The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism. ... abridged ... △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: Accepted for publication in Communications in Nonlinear Science and Numerical Simulation (CNSNS) The associated software is part of the AMUSE framework and can be downloaded from http:www.amusecode.org

A satellite orbit drift in binary near-Earth asteroids (66391) 1999 KW4 and (88710) 2001 SL9 -- Indication of the BYORP effect

Authors: P. Scheirich, P. Pravec, P. Kušnirák, K. Hornoch, J. McMahon, D. J. Scheeres, D. Čapek, D. P. Pray, H. Kučáková, A. Galád, J. Vraštil, Yu. N. Krugly, N. Moskovitz, L. D. Avner, B. Skiff, R. S. McMillan, J. A. Larsen, M. J. Brucker, A. F. Tubbiolo, W. R. Cooney, J. Gross, D. Terrell, O. Burkhonov, K. E. Ergashev, Sh. A. Ehgamberdiev , et al. (12 additional authors not shown)

Abstract: We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 taken from 2000 to 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of… ▽ More We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 taken from 2000 to 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 +/- 0.16 deg/yr2 (all quoted uncertainties are 3sigma). This means that the semimajor axis of the mutual orbit of the components of this binary system increases in time with a mean rate of 1.2 +/- 0.3 cm/yr. The data for 2001 SL9 are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 +/- 0.2 and 5.2 +/- 0.2 deg/yr2, implying that the semimajor axis of the mutual orbit of the components decreases in time with a mean rate of -2.8 +/- 0.2 or -5.1 +/- 0.2 cm/yr for the two solutions, respectively. The expanding orbit of 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides. In the case of 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit. △ Less

Submitted 4 January, 2021; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 44 pages, 17 figures. Accepted for publication in Icarus journal

Survivability of planetary systems in young and dense star clusters

Authors: Arjen van Elteren, Simon Portegies Zwart, Inti Pelupessy, Maxwell Cai, Steve McMillan

Abstract: We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced… ▽ More We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced. Four, five, or size planets per star were introduced in orbit around the 500 solar-like stars with a maximum orbital separation of 400au. Our study focuses on the production of free-floating planets. A total of 357 become unbound from a total of 2522 planets in the initial conditions of the simulation. Of these, 281 leave the cluster within the crossing timescale of the star cluster; the others remain bound to the cluster as free-floating intra-cluster planets. Five of these free-floating intra-cluster planets are captured at a later time by another star. The two main mechanisms by which planets are lost from their host star, ejection upon a strong encounter with another star or internal planetary scattering, drive the evaporation independent of planet mass of orbital separation at birth. The effect of small perturbations due to slow changes in the cluster potential are important for the evolution of planetary systems. In addition, the probability of a star to lose a planet is independent of the planet mass and independent of its initial orbital separation. As a consequence, the mass distribution of free-floating planets indistinguishable from the mass distribution of planets bound to their host star. △ Less

Submitted 11 March, 2019; v1 submitted 12 February, 2019; originally announced February 2019.

Comments: A&A in press

Journal ref: A&A 624, A120 (2019)

Collisional N-Body Dynamics Coupled to Self-Gravitating Magnetohydrodynamics Reveals Dynamical Binary Formation

Authors: Joshua E. Wall, Stephen L. W. McMillan, Mordecai-Mark Mac Low, Ralf S. Klessen, Simon Portegies Zwart

Abstract: We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for… ▽ More We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for resolving binary evolution (multiples). This combination yields unique star formation simulations that allow us to study binaries formed dynamically from interactions with both other stars and dense, magnetized gas subject to stellar feedback during the birth and early evolution of stellar clusters. We find that for massive stars, our simulations are consistent with the observed dynamical binary fractions and mass ratios. However, our binary fraction drops well below observed values for lower mass stars, presumably due to unincluded binary formation during initial star formation. Further, we observe a build up of binaries near the hard-soft boundary that may be an important mechanism driving early cluster contraction. △ Less

Submitted 4 January, 2019; originally announced January 2019.

Comments: 14 pages, 12 figures. Submitted to ApJ

Mathematical Foundations of the GraphBLAS

Authors: Jeremy Kepner, Peter Aaltonen, David Bader, Aydın Buluc, Franz Franchetti, John Gilbert, Dylan Hutchison, Manoj Kumar, Andrew Lumsdaine, Henning Meyerhenke, Scott McMillan, Jose Moreira, John D. Owens, Carl Yang, Marcin Zalewski, Timothy Mattson

Abstract: The GraphBLAS standard (GraphBlas.org) is being developed to bring the potential of matrix based graph algorithms to the broadest possible audience. Mathematically the Graph- BLAS defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the mathematics of th… ▽ More The GraphBLAS standard (GraphBlas.org) is being developed to bring the potential of matrix based graph algorithms to the broadest possible audience. Mathematically the Graph- BLAS defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the mathematics of the GraphBLAS. Graphs represent connections between vertices with edges. Matrices can represent a wide range of graphs using adjacency matrices or incidence matrices. Adjacency matrices are often easier to analyze while incidence matrices are often better for representing data. Fortunately, the two are easily connected by matrix mul- tiplication. A key feature of matrix mathematics is that a very small number of matrix operations can be used to manipulate a very wide range of graphs. This composability of small number of operations is the foundation of the GraphBLAS. A standard such as the GraphBLAS can only be effective if it has low performance overhead. Performance measurements of prototype GraphBLAS implementations indicate that the overhead is low. △ Less

Submitted 13 July, 2016; v1 submitted 18 June, 2016; originally announced June 2016.

Comments: 9 pages; 11 figures; accepted to IEEE High Performance Extreme Computing (HPEC) conference 2016. arXiv admin note: text overlap with arXiv:1504.01039

Evolution of Binary Stars in Multiple-Population Globular Clusters - II. Compact Binaries

Authors: Jongsuk Hong, Enrico Vesperini, Antonio Sollima, Stephen L. W. McMillan, Franca D'Antona, Annibale D'Ercole

Abstract: We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases mo… ▽ More We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases more rapidly than that of the FG population. The difference between the FG and SG binary fraction is qualitatively similar to but quantitatively smaller than that found for wider binaries in our previous investigations.The evolution of the radial variation of the binary fraction is driven by the interplay between binary segregation, ionization and ejection. Ionization and ejection counteract in part the effects of mass segregation but for compact binaries the effects of segregation dominate and the inner binary fraction increases during the cluster evolution. We explore the variation of the difference between the FG and the SG binary fraction with the distance from the cluster centre and its dependence on the binary binding energy and cluster structural parameters. The difference between the binary fraction in the FG and the SG populations found in our simulations is consistent with the results of observational studies finding a smaller binary fraction in the SG population. △ Less

Submitted 4 April, 2016; originally announced April 2016.

Comments: 9 pages, 12 figures. Accepted for publication in MNRAS

The NEOWISE-Discovered Comet Population and the CO+CO2 production rates

Authors: James M. Bauer, Rachel Stevenson, Emily Kramer, A. K. Mainzer, Tommy Grav, Joseph R. Masiero, Yan R. Fernández, Roc M. Cutri, John W. Dailey, Frank J. Masci, Karen J. Meech, Russel Walker, C. M. Lisse, Paul R. Weissman, Carrie R. Nugent, Sarah Sonnett, Nathan Blair, Andrew Lucas, Robert S. McMillan, Edward L. Wright, the WISE, NEOWISE Teams

Abstract: The 163 comets observed during the WISE/NEOWISE prime mission represent the largest infrared survey to date of comets, providing constraints on dust, nucleus sizes, and CO+CO2 production. We present detailed analyses of the WISE/NEOWISE comet discoveries, and discuss observations of the active comets showing 4.6 $μ$m band excess. We find a possible relation between dust and CO+CO2 production, as w… ▽ More The 163 comets observed during the WISE/NEOWISE prime mission represent the largest infrared survey to date of comets, providing constraints on dust, nucleus sizes, and CO+CO2 production. We present detailed analyses of the WISE/NEOWISE comet discoveries, and discuss observations of the active comets showing 4.6 $μ$m band excess. We find a possible relation between dust and CO+CO2 production, as well as possible differences in the sizes of long and short period comet nuclei. △ Less

Submitted 7 October, 2015; v1 submitted 28 September, 2015; originally announced September 2015.

Comments: 67 pages, 10 figures, 6 Tables; Accepted by The Astrophysical Journal

Evolution of Binary Stars in Multiple-Population Globular Clusters

Authors: Jongsuk Hong, Enrico Vesperini, Antonio Sollima, Stephen L. W. McMillan, Franca D'Antona, Annibale D'Ercole

Abstract: The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and… ▽ More The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions. △ Less

Submitted 6 March, 2015; originally announced March 2015.

Comments: 10 pages, 12 figures. Accepted for publication in MNRAS

Kinematical fingerprints of star cluster early dynamical evolution

Authors: Enrico Vesperini, Anna Lisa Varri, Stephen L. W. McMillan, Stephen E. Zepf

Abstract: We study the effects of the external tidal field on the violent relaxation phase of star clusters dynamical evolution, with particular attention to the kinematical properties of the equilibrium configurations emerging at the end of this phase.We show that star clusters undergoing the process of violent relaxation in the tidal field of their host galaxy can acquire significant internal differential… ▽ More We study the effects of the external tidal field on the violent relaxation phase of star clusters dynamical evolution, with particular attention to the kinematical properties of the equilibrium configurations emerging at the end of this phase.We show that star clusters undergoing the process of violent relaxation in the tidal field of their host galaxy can acquire significant internal differential rotation and are characterized by a distinctive radial variation of the velocity anisotropy. These kinematical properties are the result of the symmetry breaking introduced by the external tidal field in the collapse phase and of the action of the Coriolis force on the orbit of the stars. The resulting equilibrium configurations are characterized by differential rotation, with a peak located between one and two half-mass radii. As for the anisotropy, similar to clusters evolving in isolation, the systems explored in this Letter are characterized by an inner isotropic core, followed by a region of increasing radial anisotropy. However for systems evolving in an external tidal field the degree of radial anisotropy reaches a maximum in the cluster intermediate regions and then progressively decreases, with the cluster outermost regions being characterized by isotropy or a mild tangential anisotropy. Young or old but less-relaxed dynamically young star clusters may keep memory of these kinematical fingerprints of their early dynamical evolution. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: 5 pages 4 figures. Accepted for publication in MNRAS Letters

Dynamical Processes in Globular Clusters

Authors: Stephen L. W. McMillan

Abstract: Globular clusters are among the most congested stellar systems in the Universe. Internal dynamical evolution drives them toward states of high central density, while simultaneously concentrating the most massive stars and binary systems in their cores. As a result, these clusters are expected to be sites of frequent close encounters and physical collisions between stars and binaries, making them e… ▽ More Globular clusters are among the most congested stellar systems in the Universe. Internal dynamical evolution drives them toward states of high central density, while simultaneously concentrating the most massive stars and binary systems in their cores. As a result, these clusters are expected to be sites of frequent close encounters and physical collisions between stars and binaries, making them efficient factories for the production of interesting and observable astrophysical exotica. I describe some elements of the competition among stellar dynamics, stellar evolution, and other processes that control globular cluster dynamics, with particular emphasis on pathways that may lead to the formation of blue stragglers. △ Less

Submitted 13 June, 2014; originally announced June 2014.

Comments: Chapter 10, in Ecology of Blue Straggler Stars, H.M.J. Boffin, G. Carraro & G. Beccari (Eds), Astrophysics and Space Science Library, Springer

Evolution of star clusters in a cosmological tidal field

Authors: Steven Rieder, Tomoaki Ishiyama, Paul Langelaan, Junichiro Makino, Stephen L. W. McMillan, Simon Portegies Zwart

Abstract: We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializ… ▽ More We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo. △ Less

Submitted 27 September, 2013; originally announced September 2013.

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

The Astrophysical Multipurpose Software Environment

Authors: F. I. Pelupessy, A. van Elteren, N. de Vries, S. L. W. McMillan, N. Drost, S. F. Portegies Zwart

Abstract: We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guaran… ▽ More We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guarantees that the whole framework is well-suited for distributed computation. It includes facilities for unit handling and data storage. Currently it includes codes for gravitational dynamics, stellar evolution, hydrodynamics and radiative transfer. Within each domain the interfaces to the codes are as similar as possible. We describe the design and implementation of AMUSE, as well as the main components and community codes currently supported and we discuss the code interactions facilitated by the framework. Additionally, we demonstrate how AMUSE can be used to resolve complex astrophysical problems by presenting example applications. △ Less

Submitted 11 July, 2013; originally announced July 2013.

Comments: 23 pages, 25 figures, accepted for A&A

Simulating Star Clusters with the AMUSE Software Framework: I. Dependence of Cluster Lifetimes on Model Assumptions and Cluster Dissolution Modes

Authors: Alfred J. Whitehead, Stephen L. W. McMillan, Enrico Vesperini, Simon Portegies Zwart

Abstract: We perform a series of simulations of evolving star clusters using AMUSE (the Astrophysical Multipurpose Software Environment), a new community-based multi-physics simulation package, and compare our results to existing work. These simulations model a star cluster beginning with a King model distribution and a selection of power-law initial mass functions, and contain a tidal cut-off. They are evo… ▽ More We perform a series of simulations of evolving star clusters using AMUSE (the Astrophysical Multipurpose Software Environment), a new community-based multi-physics simulation package, and compare our results to existing work. These simulations model a star cluster beginning with a King model distribution and a selection of power-law initial mass functions, and contain a tidal cut-off. They are evolved using collisional stellar dynamics and include mass loss due to stellar evolution. After determining that the differences between AMUSE results and prior publications are understood, we explored the variation in cluster lifetimes due to the random realization noise introduced by transforming a King model to specific initial conditions. This random realization noise can affect the lifetime of a simulated star cluster by up to 30%. Two modes of star cluster dissolution were identified: a mass evolution curve that contains a run-away cluster dissolution with a sudden loss of mass, and a dissolution mode that does not contain this feature. We refer to these dissolution modes as "dynamical" and "relaxation" dominated respectively. For Salpeter-like initial mass functions, we determined the boundary between these two modes in terms of the dynamical and relaxation time scales. △ Less

Submitted 16 October, 2013; v1 submitted 11 March, 2013; originally announced March 2013.

Comments: 12 pages, 10 figures, accepted by ApJ 2013-10-07

The Pan-STARRS Moving Object Processing System

Authors: Larry Denneau, Robert Jedicke, Tommy Grav, Mikael Granvik, Jeremy Kubica, Andrea Milani, Peter Veres, Richard Wainscoat, Daniel Chang, Francesco Pierfederici, N. Kaiser, K. C. Chambers, J. N. Heasley, Eugene. A. Magnier, P. A. Price, Jonathan Myers, Jan Kleyna, Henry Hsieh, Davide Farnocchia, Chris Waters, W. H. Sweeney, Denver Green, Bryce Bolin, W. S. Burgett, J. S. Morgan , et al. (19 additional authors not shown)

Abstract: We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan… ▽ More We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development. △ Less

Submitted 28 February, 2013; originally announced February 2013.

Comments: 57 Pages, 26 Figures, 13 Tables

The MICA Experiment: Astrophysics in Virtual Worlds

Authors: S. G. Djorgovski, Piet Hut, Rob Knop, Giuseppe Longo, Steve McMillan, Enrico Vesperini, Ciro Donalek, Matthew Graham, Asish Mahabal, Franz Sauer, Charles White, Crista Lopes

Abstract: We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. W… ▽ More We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted. △ Less

Submitted 28 January, 2013; originally announced January 2013.

Comments: 10 pages, 8 figures; invited paper for the refereed proc. of the SLACTIONS 2012 conference

Journal ref: Ref. proc. of the SLACTIONS 2012 International Research Conference on Virtual Worlds,ceds. L. Morgado, Y. Sivan, A.M. Maia, et al., pp. 49-58 (2012). Vila Real, PT: UTAD, ISBN 978-989-704-102-0, avail. at http://hdl.handle.net/10348/2305

Dynamical Evolution and Spatial Mixing of Multiple Population Globular Clusters

Authors: Enrico Vesperini, Stephen L. W. McMillan, Francesca D'Antona, Annibale D'Ercole

Abstract: In this paper we study the long-term dynamical evolution of multiple-population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars.In previous studies we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we… ▽ More In this paper we study the long-term dynamical evolution of multiple-population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars.In previous studies we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we explore the time scales and the dynamics of the spatial mixing of the FG and the SG populations and their dependence on the SG initial concentration.Our simulations show that, as the evolution proceeds, the radial profile of the SG/FG number ratio, NSG/NFG, is characterized by three regions: 1) a flat inner part; 2) a declining part in which FG stars are increasingly dominant; and 3) an outer region where the NSG/NFG profile flattens again (the NSG/NFG profile may rise slightly again in the outermost cluster regions). The radial variation of NSG/NFG implies that the fraction of SG stars determined by observations covering a limited range of radial distances is not, in general, equal to the SG global fraction, (NSG/NFG)glob. The distance at which NSG/NFG equals (NSG/NFG)glob is approximately between 1 and 2 cluster half-mass radii. The results of our simulations suggest that in many Galactic globular clusters the SG should still be more spatially concentrated than the FG.[abridged] △ Less

Submitted 11 December, 2012; originally announced December 2012.

Comments: 10 pages 9 figures. Accepted for publication in MNRAS

Rapid Mass Segregation in Massive Star Clusters

Authors: Stephen McMillan, Enrico Vesperini, Nicholas Kruczek

Abstract: Several dynamical scenarios have been proposed that can lead to prompt mass segregation on the crossing time scale of a young cluster. They generally rely on cool and/or clumpy initial conditions, and are most relevant to small systems. As a counterpoint, we present a novel dynamical mechanism that can operate in relatively large, homogeneous, cool or cold systems. This mechanism may be important… ▽ More Several dynamical scenarios have been proposed that can lead to prompt mass segregation on the crossing time scale of a young cluster. They generally rely on cool and/or clumpy initial conditions, and are most relevant to small systems. As a counterpoint, we present a novel dynamical mechanism that can operate in relatively large, homogeneous, cool or cold systems. This mechanism may be important in understanding the assembly of large mass-segregated clusters from smaller clumps. △ Less

Submitted 30 October, 2012; originally announced October 2012.

Comments: 2 pages, 2 figures; Invited talk presented as Special Session 1 "Origin and Complexity of Massive Star Clusters" at IAU General Assembly XXVIII, Beijing 2012. To appear in Highlights of Astronomy

Physical Parameters of Asteroids Estimated from the WISE 3 Band Data and NEOWISE Post-Cryogenic Survey

Authors: A. Mainzer, T. Grav, J. Masiero, J. Bauer, R. M. Cutri, R. S. McMillan, C. R. Nugent, D. Tholen, R. Walker, E. L. Wright

Abstract: Enhancements to the science data processing pipeline of NASA's Wide-field Infrared Explorer (WISE) mission, collectively known as NEOWISE, resulted in the detection of $>$158,000 minor planets in four infrared wavelengths during the fully cryogenic portion of the mission. Following the depletion of its cryogen, NASA's Planetary Science Directorate funded a four month extension to complete the surv… ▽ More Enhancements to the science data processing pipeline of NASA's Wide-field Infrared Explorer (WISE) mission, collectively known as NEOWISE, resulted in the detection of $>$158,000 minor planets in four infrared wavelengths during the fully cryogenic portion of the mission. Following the depletion of its cryogen, NASA's Planetary Science Directorate funded a four month extension to complete the survey of the inner edge of the Main Asteroid Belt and to detect and discover near-Earth objects (NEOs). This extended survey phase, known as the NEOWISE Post-Cryogenic Survey, resulted in the detection of $\sim$6500 large Main Belt asteroids and 88 NEOs in its 3.4 and 4.6 $μ$m channels. During the Post-Cryogenic Survey, NEOWISE discovered and detected a number of asteroids co-orbital with the Earth and Mars, including the first known Earth Trojan. We present preliminary thermal fits for these and other NEOs detected during the 3-Band Cryogenic and Post-Cryogenic Surveys. △ Less

Submitted 1 October, 2012; originally announced October 2012.

Comments: ApJL accepted

Characterizing Subpopulations within the Near Earth Objects with NEOWISE: Preliminary Results

Authors: A. Mainzer, T. Grav, J. Masiero, J. Bauer, R. S. McMillan, J. Giorgini, T. Spahr, R. M. Cutri, D. J. Tholen, R. Jedicke, R. Walker, E. Wright, C. R. Nugent

Abstract: We present the preliminary results of an analysis of the sub-populations within the near-Earth asteroids, including the Atens, Apollos, Amors, and those that are considered potentially hazardous using data from the Wide-field Infrared Survey Explorer (WISE). In order to extrapolate the sample of objects detected by WISE to the greater population, we determined the survey biases for asteroids detec… ▽ More We present the preliminary results of an analysis of the sub-populations within the near-Earth asteroids, including the Atens, Apollos, Amors, and those that are considered potentially hazardous using data from the Wide-field Infrared Survey Explorer (WISE). In order to extrapolate the sample of objects detected by WISE to the greater population, we determined the survey biases for asteroids detected by the project's automated moving object processing system (known as NEOWISE) as a function of diameter, visible albedo, and orbital elements. Using this technique, we are able to place constraints on the number of potentially hazardous asteroids (PHAs) larger than 100 m and find that there are $\sim4700\pm1450$ such objects. As expected, the Atens, Apollos, and Amors are revealed by WISE to have somewhat different albedo distributions, with the Atens being brighter than the Amors. The cumulative size distributions of the various near-Earth object (NEO) subgroups vary slightly between 100 m and 1 km. A comparison of the observed orbital elements of the various sub-populations of the NEOs with the current best model is shown. △ Less

Submitted 16 May, 2012; originally announced May 2012.

Comments: ApJ accepted

Multi-physics simulations using a hierarchical interchangeable software interface

Authors: Simon Portegies Zwart, Steve McMillan, Arjen van Elteren, Inti Pelupessy, Nathan de Vries

Abstract: We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery f… ▽ More We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery from errors. The framework can be efficiently implemented and incurs an acceptable overhead. In practice, we measure the time spent in the framework to be less than 1% of the wall-clock time. Due to the unified structure of the interface, incorporating multiple modules addressing the same physics in different ways is relatively straightforward. Different modules may be advanced serially or in parallel. Despite initial concerns, we have encountered relatively few problems with this strict separation between modules, and the results of our simulations are consistent with earlier results using more traditional monolithic approaches. This framework provides a platform to combine existing simulation codes or develop new physical solver codes within a rich "ecosystem" of interchangeable modules. △ Less

Submitted 16 September, 2012; v1 submitted 24 April, 2012; originally announced April 2012.

Comments: 36 pages, 12 figures and source code (accepted for publication in Computer Physics Communications)

A revised asteroid polarization-albedo relationship using WISE/NEOWISE data

Authors: Joseph R. Masiero, A. K. Mainzer, T. Grav, J. M. Bauer, E. L. Wright, R. S. McMillan, D. J. Tholen, A. W. Blain

Abstract: We present a reanalysis of the relationship between asteroid albedo and polarization properties using the albedos derived from the Wide-field Infrared Survey Explorer. We find that the function that best describes this relation is a three-dimensional linear fit in the space of log(albedo)-log(polarization slope)-log(minimum polarization). When projected to two dimensions the parameters of the fit… ▽ More We present a reanalysis of the relationship between asteroid albedo and polarization properties using the albedos derived from the Wide-field Infrared Survey Explorer. We find that the function that best describes this relation is a three-dimensional linear fit in the space of log(albedo)-log(polarization slope)-log(minimum polarization). When projected to two dimensions the parameters of the fit are consistent with those found in previous work. We also define p* as the quantity of maximal polarization variation when compared with albedo and present the best fitting albedo-p* relation. Some asteroid taxonomic types stand out in this three-dimensional space, notably the E, B, and M Tholen types, while others cluster in clumps coincident with the S- and C-complex bodies. We note that both low albedo and small (D<30 km) asteroids are under-represented in the polarimetric sample, and we encourage future polarimetric surveys to focus on these bodies. △ Less

Submitted 10 February, 2012; originally announced February 2012.

Comments: 16 pages, Accepted to ApJ

PSDF: Particle Stream Data Format for N-Body Simulations

Authors: Will M. Farr, Jeff Ames, Piet Hut, Junichiro Makino, Steve McMillan, Takayuki Muranushi, Koichi Nakamura, Keigo Nitadori, Simon Portegies Zwart

Abstract: We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same… ▽ More We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same approach as XML but with a simpler syntax. Together with a specification of PSDF, we provide background and motivation, as well as specific examples in a variety of computer languages. We also offer a web site from which these examples can be retrieved, in order to make it easy to augment existing codes in order to give them the option to produce PSDF output. △ Less

Submitted 9 January, 2012; originally announced January 2012.

Comments: 5 pages; submitted to New Astronomy

Journal ref: New Astronomy 17 (2012), pp. 520--523

WISE/NEOWISE observations of Active Bodies in the Main Belt

Authors: James M. Bauer, A. K. Mainzer, Tommy Grav, Russell G. Walker, Joseph R. Masiero, Erin K. Blauvelt, Robert S. McMillan, Yan R. Fernández, Karen J. Meech, Carey M. Lisse, Roc M. Cutri, John W. Dailey, David J. Tholen, Timm Riesen, Laurie Urban, Alain Khayat, George Pearman, James V. Scotti, Emily Kramer, De'Andre Cherry, Thomas Gautier, Stephanie Gomillion, Jessica Watkins, Edward L. Wright, the WISE Team

Abstract: We report results based on mid-infrared photometry of 5 active main belt objects (AMBOs) detected by the Wide-field Infrared Survey Explorer (WISE) spacecraft. Four of these bodies, P/2010 R2 (La Sagra), 133P/Elst-Pizarro, (596) Scheila, and 176P/LINEAR, showed no signs of activity at the time of the observations, allowing the WISE detections to place firm constraints on their diameters and albedo… ▽ More We report results based on mid-infrared photometry of 5 active main belt objects (AMBOs) detected by the Wide-field Infrared Survey Explorer (WISE) spacecraft. Four of these bodies, P/2010 R2 (La Sagra), 133P/Elst-Pizarro, (596) Scheila, and 176P/LINEAR, showed no signs of activity at the time of the observations, allowing the WISE detections to place firm constraints on their diameters and albedos. Geometric albedos were in the range of a few percent, and on the order of other measured comet nuclei. P/2010 A2 was observed on April 2-3, 2010, three months after its peak activity. Photometry of the coma at 12 and 22 μm combined with ground-based visible-wavelength measurements provides constraints on the dust particle mass distribution (PMD), dlogn/dlogm, yielding power-law slope values of α = -0.5 +/- 0.1. This PMD is considerably more shallow than that found for other comets, in particular inbound particle fluence during the Stardust encounter of comet 81P/Wild 2. It is similar to the PMD seen for 9P/Tempel 1 in the immediate aftermath of the Deep Impact experiment. Upper limits for CO2 & CO production are also provided for each AMBO and compared with revised production numbers for WISE observations of 103P/Hartley 2. △ Less

Submitted 13 December, 2011; originally announced December 2011.

Comments: 32 Pages, including 5 Figures

Simulations of Dense Stellar Systems with the AMUSE Software Toolkit

Authors: Stephen McMillan, Simon Portegies Zwart, Arjen van Elteren, Alfred Whitehead

Abstract: We describe AMUSE, the Astrophysical Multipurpose Software Environment, a programming framework designed to manage multi-scale, multi-physics simulations in a hierarchical, extensible, and internally consistent way. Constructed as a collection of individual modules, AMUSE allows computational tools for different physical domains to be easily combined into a single task. It facilitates the coupling… ▽ More We describe AMUSE, the Astrophysical Multipurpose Software Environment, a programming framework designed to manage multi-scale, multi-physics simulations in a hierarchical, extensible, and internally consistent way. Constructed as a collection of individual modules, AMUSE allows computational tools for different physical domains to be easily combined into a single task. It facilitates the coupling of modules written in different languages by providing inter-language tools and a standard programming interface that represents a balance between generality and computational efficiency. The framework currently incorporates the domains of stellar dynamics, stellar evolution, gas dynamics, and radiative transfer. We present some applications of the framework and outline plans for future development of the package. △ Less

Submitted 16 November, 2011; originally announced November 2011.

Comments: 8 pages, 4 figures; Invited talk; to appear in the Proceedings of "Advances in Computational Astrophysics: methods, tools and outcomes," Cefalù (Sicily, Italy), June 13-17, 2011

NEOWISE Studies of Asteroids with Sloan Photometry: Preliminary Results

Authors: A. Mainzer, J. Masiero, T. Grav, J. Bauer, D. J. Tholen, R. S. McMillan, E. Wright, T. Spahr, R. M. Cutri, R. Walker, W. Mo, J. Watkins, E. Hand, C. Maleszewski

Abstract: We have combined the NEOWISE and Sloan Digital Sky Survey data to study the albedos of 24,353 asteroids with candidate taxonomic classifications derived using Sloan photometry. We find a wide range of moderate to high albedos for candidate S-type asteroids that are analogous to the S-complex defined by previous spectrophotometrically-based taxonomic systems. The candidate C-type asteroids, while g… ▽ More We have combined the NEOWISE and Sloan Digital Sky Survey data to study the albedos of 24,353 asteroids with candidate taxonomic classifications derived using Sloan photometry. We find a wide range of moderate to high albedos for candidate S-type asteroids that are analogous to the S-complex defined by previous spectrophotometrically-based taxonomic systems. The candidate C-type asteroids, while generally very dark, have a tail of higher albedos that overlaps the S types. The albedo distribution for asteroids with a photometrically derived Q classification is extremely similar to those of the S types. Asteroids with similar colors to (4) Vesta have higher albedos than the S types, and most have orbital elements similar to known Vesta family members. Finally, we show that the relative reflectance at 3.4 and 4.6 $μ$m is higher for D-type asteroids and suggest that their red visible and near-infrared spectral slope extends out to these wavelengths. Understanding the relationship between size, albedo, and taxonomic classification is complicated by the fact that the objects with classifications were selected from the visible/near-infrared Sloan Moving Object Catalog, which is biased against fainter asteroids, including those with lower albedos. △ Less

Submitted 22 October, 2011; originally announced October 2011.

Comments: ApJ accepted

Multi-physics simulations using a hierarchical interchangeable software interface

Authors: Simon Portegies Zwart, Stephen McMillan, Inti Pelupessy, Arjen van Elteren

Abstract: We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified software interface. Our framework is intrinsically parallel, and conveniently separates all components in memory. It performs unit conversion between different modules automatically and defines common data structures to communicate across different codes. We use the framework to… ▽ More We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified software interface. Our framework is intrinsically parallel, and conveniently separates all components in memory. It performs unit conversion between different modules automatically and defines common data structures to communicate across different codes. We use the framework to simulate embedded star clusters. For this purpose we couple solvers for gravitational dynamics, stellar evolution and hydrodynamics to self consistently resolve the dynamical evolution simultaneousy with the internal nuclear evolution of the stars and the hydrodynamic response of the gas. We find, in contrast to earlier studies, that the survival of a young star cluster depends only weakly on the efficiency of star formation. The main reason for this weak dependency is the asymmetric expulsion of the embedding gas from the cluster. △ Less

Submitted 12 October, 2011; originally announced October 2011.

Comments: To appear in "Advances in Computational Astrophysics: methods, tools and outcomes", ASP Conference Series, R. Capuzzo-Dolcetta, M. Limongi and A. Tornambè eds

WISE/NEOWISE Observations of the Hilda Population: Preliminary Results

Authors: T. Grav, A. K. Mainzer, J. Bauer, J. Masiero, T. Spahr, R. S. McMillan, R. Walker, R. Cutri, E. Wright, P. R. Eisenhardt, E. Blauvelt, E. DeBaun, D. Elsbury, T. Gautier IV, S. Gomillion, E. Hand, A. Wilkins

Abstract: We present the preliminary analysis of 1023 known asteroids in the Hilda region of the Solar System observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). The sizes of the Hildas observed range from $\sim 3 - 200$km. We find no size - albedo dependency as reported by other projects. The albedos of our sample are low, with a weighted mean value $p_V = 0.055\pm0.018$, f… ▽ More We present the preliminary analysis of 1023 known asteroids in the Hilda region of the Solar System observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). The sizes of the Hildas observed range from $\sim 3 - 200$km. We find no size - albedo dependency as reported by other projects. The albedos of our sample are low, with a weighted mean value $p_V = 0.055\pm0.018$, for all sizes sampled by the NEOWISE survey. We observed a significant fraction of the objects in the two known collisional families in the Hilda population. It is found that the Hilda collisional family is brighter, with weighted mean albedo of $p_V = 0.061\pm0.011$, than the general population and dominated by D-type asteroids, while the Schubart collisional family is darker, with weighted mean albedo of ($p_V = 0.039\pm0.013$). Using the reflected sunlight in the two shortest WISE bandpasses we are able to derive a method for taxonomic classification of $\sim 10%$ of the Hildas detected in the NEOWISE survey. For the Hildas with diameter larger than 30km there are $67^{+7}_{-15}%$ D-type asteroids and $26^{+17}_{-5}%$ C-/P-type asteroids (with the majority of these being P-types). △ Less

Submitted 3 October, 2011; originally announced October 2011.

Comments: Accepted for publication in Astrophysical Journal. Electronic table to be published on the publishers website

Journal ref: ApJ (2012), v744, p197

WISE/NEOWISE Observations of the Jovian Trojans: Preliminary Results

Authors: T. Grav, A. K. Mainzer, J. Bauer, J. Masiero, T. Spahr, R. S. McMillan, R. Walker, R. Cutri, E. Wright, E. Blauvelt, E. DeBaun, D. Elsbury, T. Gautier IV, S. Gomillion, E. Hand, A. Wilkins

Abstract: We present the preliminary analysis of over 1739 known and 349 candidate Jovian Trojans observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). With this survey the available diameters, albedos and beaming parameters for the Jovian Trojans have been increased by more than an order of magnitude compared to previous surveys. We find that the Jovian Trojan population is v… ▽ More We present the preliminary analysis of over 1739 known and 349 candidate Jovian Trojans observed by the NEOWISE component of the Wide-field Infrared Survey Explorer (WISE). With this survey the available diameters, albedos and beaming parameters for the Jovian Trojans have been increased by more than an order of magnitude compared to previous surveys. We find that the Jovian Trojan population is very homogenous for sizes larger than $\sim10$km (close to the detection limit of WISE for these objects). The observed sample consists almost exclusively of low albedo objects, having a mean albedo value of $0.07\pm0.03$. The beaming parameter was also derived for a large fraction of the observed sample, and it is also very homogenous with an observed mean value of $0.88\pm0.13$. Preliminary debiasing of the survey shows our observed sample is consistent with the leading cloud containing more objects than the trailing cloud. We estimate the fraction to be N(leading)/N(trailing) $\sim 1.4 \pm 0.2$, lower than the $1.6 \pm 0.1$ value derived by others. △ Less

Submitted 3 October, 2011; originally announced October 2011.

Comments: Accepted for publication in Astrophysical Journal. Electronic table will be available at the publishers website

Journal ref: ApJ (2011), v742, p40

NEOWISE Studies of Spectrophotometrically Classified Asteroids: Preliminary Results

Authors: A. Mainzer, T. Grav, J. Masiero, J. Bauer, E. Hand, D. Tholen, R. S. McMillan, T. Spahr, R. M. Cutri, E. Wright, J. Watkins, W. Mo, C. Maleszewski

Abstract: The NEOWISE dataset offers the opportunity to study the variations in albedo for asteroid classification schemes based on visible and near-infrared observations for a large sample of minor planets. We have determined the albedos for nearly 1900 asteroids classified by the Tholen, Bus and Bus-DeMeo taxonomic classification schemes. We find that the S-complex spans a broad range of bright albedos, p… ▽ More The NEOWISE dataset offers the opportunity to study the variations in albedo for asteroid classification schemes based on visible and near-infrared observations for a large sample of minor planets. We have determined the albedos for nearly 1900 asteroids classified by the Tholen, Bus and Bus-DeMeo taxonomic classification schemes. We find that the S-complex spans a broad range of bright albedos, partially overlapping the low albedo C-complex at small sizes. As expected, the X-complex covers a wide range of albedos. The multi-wavelength infrared coverage provided by NEOWISE allows determination of the reflectivity at 3.4 and 4.6 $μ$m relative to the visible albedo. The direct computation of the reflectivity at 3.4 and 4.6 $μ$m enables a new means of comparing the various taxonomic classes. Although C, B, D and T asteroids all have similarly low visible albedos, the D and T types can be distinguished from the C and B types by examining their relative reflectance at 3.4 and 4.6 $μ$m. All of the albedo distributions are strongly affected by selection biases against small, low albedo objects, as all objects selected for taxonomic classification were chosen according to their visible light brightness. Due to these strong selection biases, we are unable to determine whether or not there are correlations between size, albedo and space weathering. We argue that the current set of classified asteroids makes any such correlations difficult to verify. A sample of taxonomically classified asteroids drawn without significant albedo bias is needed in order to perform such an analysis. △ Less

Submitted 29 September, 2011; originally announced September 2011.

Comments: Accepted to ApJ

NEOWISE Observations of Near-Earth Objects: Preliminary Results

Authors: A. Mainzer, T. Grav, J. Bauer, J. Masiero, R. S. McMillan, R. M. Cutri, R. Walker, E. Wright, P. Eisenhardt, D. J. Tholen, T. Spahr, R. Jedicke, L. Denneau, E. DeBaun, D. Elsbury, T. Gautier, S. Gomillion, E. Hand, W. Mo, J. Watkins, A. Wilkins, G. L. Bryngelson, A. Del Pino Molina, S. Desai, M. Go'mez Camus , et al. (12 additional authors not shown)

Abstract: With the NEOWISE portion of the \emph{Wide-field Infrared Survey Explorer} (WISE) project, we have carried out a highly uniform survey of the near-Earth object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 $μ$m, allowing us to refine estimates of their numbers, sizes, and albedos. The NEOWISE survey detected NEOs the same way whether they were previously known or not, subje… ▽ More With the NEOWISE portion of the \emph{Wide-field Infrared Survey Explorer} (WISE) project, we have carried out a highly uniform survey of the near-Earth object (NEO) population at thermal infrared wavelengths ranging from 3 to 22 $μ$m, allowing us to refine estimates of their numbers, sizes, and albedos. The NEOWISE survey detected NEOs the same way whether they were previously known or not, subject to the availability of ground-based follow-up observations, resulting in the discovery of more than 130 new NEOs. The survey's uniformity in sensitivity, observing cadence, and image quality have permitted extrapolation of the 428 near-Earth asteroids (NEAs) detected by NEOWISE during the fully cryogenic portion of the WISE mission to the larger population. We find that there are 981$\pm$19 NEAs larger than 1 km and 20,500$\pm$3000 NEAs larger than 100 m. We show that the Spaceguard goal of detecting 90% of all 1 km NEAs has been met, and that the cumulative size distribution is best represented by a broken power law with a slope of 1.32$\pm$0.14 below 1.5 km. This power law slope produces $\sim13,200\pm$1,900 NEAs with $D>$140 m. Although previous studies predict another break in the cumulative size distribution below $D\sim$50-100 m, resulting in an increase in the number of NEOs in this size range and smaller, we did not detect enough objects to comment on this increase. The overall number for the NEA population between 100-1000 m are lower than previous estimates. The numbers of near-Earth comets will be the subject of future work. △ Less

Submitted 29 September, 2011; originally announced September 2011.

Comments: Accepted to ApJ

Main Belt Asteroids with WISE/NEOWISE I: Preliminary Albedos and Diameters

Authors: Joseph R. Masiero, A. K. Mainzer, T. Grav, J. M. Bauer, R. M. Cutri, J. Dailey, P. R. M. Eisenhardt, R. S. McMillan, T. B. Spahr, M. F. Skrutskie, D. Tholen, R. G. Walker, E. L. Wright, E. DeBaun, D. Elsbury, T. Gautier IV, S. Gomillion, A. Wilkins

Abstract: We present initial results from the Wide-field Infrared Survey Explorer (WISE), a four-band all-sky thermal infrared survey that produces data well suited to measuring the physical properties of asteroids, and the NEOWISE enhancement to the WISE mission allowing for detailed study of Solar system objects. Using a NEATM thermal model fitting routine we compute diameters for over 100,000 Main Belt a… ▽ More We present initial results from the Wide-field Infrared Survey Explorer (WISE), a four-band all-sky thermal infrared survey that produces data well suited to measuring the physical properties of asteroids, and the NEOWISE enhancement to the WISE mission allowing for detailed study of Solar system objects. Using a NEATM thermal model fitting routine we compute diameters for over 100,000 Main Belt asteroids from their IR thermal flux, with errors better than 10%. We then incorporate literature values of visible measurements (in the form of the H absolute magnitude) to determine albedos. Using these data we investigate the albedo and diameter distributions of the Main Belt. As observed previously, we find a change in the average albedo when comparing the inner, middle, and outer portions of the Main Belt. We also confirm that the albedo distribution of each region is strongly bimodal. We observe groupings of objects with similar albedos in regions of the Main Belt associated with dynamical breakup families. Asteroid families typically show a characteristic albedo for all members, but there are notable exceptions to this. This paper is the first look at the Main Belt asteroids in the WISE data, and only represents the preliminary, observed raw size and albedo distributions for the populations considered. These distributions are subject to survey biases inherent to the NEOWISE dataset and cannot yet be interpreted as describing the true populations; the debiased size and albedo distributions will be the subject of the next paper in this series. △ Less

Submitted 19 September, 2011; originally announced September 2011.

Comments: Accepted to ApJ. Online table to also appear on the publisher's website

WISE/NEOWISE observations of comet 103P/Hartley 2

Authors: James M. Bauer, Russell G. Walker, A. K. Mainzer, Joseph R. Masiero, Tommy Grav, John W. Dailey, Robert S. McMillan, Carey M. Lisse, Yan R. Fernández, Karen J. Meech, Jana Pittichova, Erin K. Blauvelt, Frank J. Masci, Michael F. A'Hearn, Roc M. Cutri, James V. Scotti, David J. Tholen, Emily DeBaun, Ashlee Wilkins, Emma Hand, Edward L. Wright, the WISE Collaboration

Abstract: We report results based on mid-infrared photometry of comet 103P/Hartley 2 taken during May 4-13, 2010 (when the comet was at a heliocentric distance of 2.3 AU, and an observer distance of 2.0 AU) by the Wide-field Infrared Survey Explorer (Wright et al. 2010). Photometry of the coma at 22 microns and data from the University of Hawaii 2.2-m telescope obtained on May 22, 2010 provide constraints o… ▽ More We report results based on mid-infrared photometry of comet 103P/Hartley 2 taken during May 4-13, 2010 (when the comet was at a heliocentric distance of 2.3 AU, and an observer distance of 2.0 AU) by the Wide-field Infrared Survey Explorer (Wright et al. 2010). Photometry of the coma at 22 microns and data from the University of Hawaii 2.2-m telescope obtained on May 22, 2010 provide constraints on the dust particle size distribution, dlogn/dlogm, yielding power-law slope values of alpha = -0.97 +/- 0.10, steeper than that found for the inbound particle fluence during the Stardust encounter of comet 81P/Wild 2 (Green et al. 2004). The extracted nucleus signal at 12 microns is consistent with a body of average spherical radius of 0.6 +/- 0.2 km (one standard deviation), assuming a beaming parameter of 1.2. The 4.6 micron-band signal in excess of dust and nucleus reflected and thermal contributions may be attributed to carbon monoxide or carbon dioxide emission lines and provides limits and estimates of species production. Derived carbon dioxide coma production rates are 3.5(+/- 0.9) \times 10^24 molecules per second. Analyses of the trail signal present in the stacked image with an effective exposure time of 158.4 seconds yields optical-depth values near 9 x 10^-10 at a delta mean anomaly of 0.2 deg trailing the comet nucleus, in both 12 and 22 μm bands. A minimum chi-squared analysis of the dust trail position yields a beta-parameter value of 1.0 x10^-4, consistent with a derived mean trail-grain diameter of 1.1/ρ cm for grains of ρ g/cm^3 density. This leads to a total detected trail mass of at least 4 x 10^10 ρ kg. △ Less

Submitted 13 July, 2011; originally announced July 2011.

Comments: Accepted to ApJ

Thermal Model Calibration for Minor Planets Observed with WISE/NEOWISE: Comparison with IRAS

Authors: A. Mainzer, T. Grav, J. Masiero, J. Bauer, E. Wright, R. M. Cutri, R. Walker, R. S. McMillan

Abstract: With thermal infrared observations detected by the NEOWISE project, we have measured diameters for 1742 minor planets that were also observed by the Infrared Astronomical Satellite (IRAS). We have compared the diameters and albedo derived by applying a spherical thermal model to the objects detected by NEOWISE and find that they are in generally good agreement with the IRAS values. We have shown t… ▽ More With thermal infrared observations detected by the NEOWISE project, we have measured diameters for 1742 minor planets that were also observed by the Infrared Astronomical Satellite (IRAS). We have compared the diameters and albedo derived by applying a spherical thermal model to the objects detected by NEOWISE and find that they are in generally good agreement with the IRAS values. We have shown that diameters computed from NEOWISE data are often less systematically biased than those found with IRAS. This demonstrates that the NEOWISE dataset can provide accurate physical parameters for the >157,000 minor planets that were detected by NEOWISE. △ Less

Submitted 28 June, 2011; originally announced June 2011.

Comments: ApJ Letters accepted

Binary Star Disruption in Globular Clusters with Multiple Stellar Populations

Authors: Enrico Vesperini, Stephen L. W. McMillan, Francesca D'Antona, Annibale D'Ercole

Abstract: The discovery of multiple stellar populations in globular clusters raises fundamental questions concerning the formation and dynamical history of these systems. In a previous study aimed at exploring the formation of second-generation (SG) stars from the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution of the cluster, we showed that SG stars are expected to form in… ▽ More The discovery of multiple stellar populations in globular clusters raises fundamental questions concerning the formation and dynamical history of these systems. In a previous study aimed at exploring the formation of second-generation (SG) stars from the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution of the cluster, we showed that SG stars are expected to form in a dense subsystem concentrated in the inner regions of the FG cluster. In this paper we explore the implications of the structural properties of multiple-population clusters, and in particular the presence of the inner SG subsystem, for the disruption of binary stars. We quantify the enhancement of the binary disruption rate due to the presence of the central SG subsystem for a number of different initial conditions. Our calculations show that SG binaries, which are assumed to be more concentrated in the cluster inner regions, are disrupted at a substantially larger rate than FG binaries. Assuming a similar initial fraction of FG and SG binaries, our dynamical study indicates that the SG population is now expected to contain a significantly smaller binary fraction than the FG population. △ Less

Submitted 3 June, 2011; originally announced June 2011.

Comments: 6 pages, 6 figures. Accepted for publication in MNRAS

Thermal Model Calibration for Minor Planets Observed with WISE/NEOWISE

Authors: A. Mainzer, T. Grav, J. Masiero, J. Bauer, E. Wright, R. M. Cutri, R. S. McMillan, M. Cohen, M. Ressler, P. Eisenhardt

Abstract: With the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), we have observed over 157,000 minor planets (Mainzer et al. 2011). Included in these are a number of near-Earth objects, Main Belt Asteroids, and irregular satellites which have well-measured physical properties via radar, occultation and in situ imaging. We have used these objects to validate models of thermal models using t… ▽ More With the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), we have observed over 157,000 minor planets (Mainzer et al. 2011). Included in these are a number of near-Earth objects, Main Belt Asteroids, and irregular satellites which have well-measured physical properties via radar, occultation and in situ imaging. We have used these objects to validate models of thermal models using the WISE measurements, as well as the color corrections derived in Wright et al. (2010) for the four WISE bandpasses as a function of effective temperature. We have used 50 objects with diameters measured by radar, occultation or in situ imaging to characterize the systematic errors implicit in using the WISE data with a faceted spherical NEATM model to compute diameters and albedos. By using the previously measured diameters and H magnitudes with a spherical NEATM model, we compute the predicted fluxes after applying the color corrections given in Wright et al. (2010) in the WISE bands and compare them to the measured magnitudes. We find minimum systematic flux errors of 5-10%, yielding minimum relative diameter and albedo errors of ~10% and ~20%, respectively. Visible albedos for the objects are computed and compared to the albedos at 3.4 and 4.6 microns, which contain a mix of reflected sunlight and thermal emission for most asteroids. We derive a linear relationship between subsolar temperature and effective temperature, which allows the color corrections given in Wright et al. (2010) to be used for asteroids by computing only subsolar temperature instead of a faceted thermal model. The thermal models derived in this paper are not intended to supplant previous measurements made using radar or spacecraft imaging; rather, we have used them to characterize the errors that should be expected when computing diameters and albedos of WISE asteroids using a spherical NEATM model. △ Less

Submitted 5 May, 2011; originally announced May 2011.

Comments: ApJ accepted

Stellar Astrophysics with a Dispersed Fourier Transform Spectrograph. II. Orbits of Double-lined Spectroscopic Binaries

Authors: Bradford B. Behr, Andrew T. Cenko, Arsen R. Hajian, Robert S. McMillan, Marc Murison, Jeff Meade, Robert Hindsley

Abstract: We present orbital parameters for six double-lined spectroscopic binaries (iota Pegasi, omega Draconis, 12 Bootis, V1143 Cygni, beta Aurigae, and Mizar A) and two double-lined triple star systems (kappa Pegasi and eta Virginis). The orbital fits are based upon high-precision radial velocity observations made with a dispersed Fourier Transform Spectrograph, or dFTS, a new instrument which combines… ▽ More We present orbital parameters for six double-lined spectroscopic binaries (iota Pegasi, omega Draconis, 12 Bootis, V1143 Cygni, beta Aurigae, and Mizar A) and two double-lined triple star systems (kappa Pegasi and eta Virginis). The orbital fits are based upon high-precision radial velocity observations made with a dispersed Fourier Transform Spectrograph, or dFTS, a new instrument which combines interferometric and dispersive elements. For some of the double-lined binaries with known inclination angles, the quality of our RV data permits us to determine the masses M_1 and M_2 of the stellar components with relative errors as small as 0.2%. △ Less

Submitted 7 April, 2011; originally announced April 2011.

Comments: 41 pages, 8 figures, accepted by AJ

Preliminary Results from NEOWISE: An Enhancement to the Wide-field Infrared Survey Explorer for Solar System Science

Authors: A. Mainzer, J. Bauer, T. Grav, J. Masiero, R. M. Cutri, J. Dailey, P. Eisenhardt, R. S. McMillan, E. Wright, R. Walker, R. Jedicke, T. Spahr, D. Tholen, R. Alles, R. Beck, H. Brandenburg, T. Conrow, T. Evans, J. Fowler, T. Jarrett, K. Marsh, F. Masci, H. McCallon, S. Wheelock, M. Wittman , et al. (10 additional authors not shown)

Abstract: The \emph{Wide-field Infrared Survey Explorer} has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the \emph{Infrared Astronomical Satellite} and the \emph{Cosmic Background Explorer}. NASA's Planetary Science Division has funded an enhancement to the \WISE\ data processing system called "NEOWISE" that allo… ▽ More The \emph{Wide-field Infrared Survey Explorer} has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the \emph{Infrared Astronomical Satellite} and the \emph{Cosmic Background Explorer}. NASA's Planetary Science Division has funded an enhancement to the \WISE\ data processing system called "NEOWISE" that allows detection and archiving of moving objects found in the \WISE\ data. NEOWISE has mined the \WISE\ images for a wide array of small bodies in our Solar System, including Near-Earth Objects (NEOs), Main Belt asteroids, comets, Trojans, and Centaurs. By the end of survey operations in February 2011, NEOWISE identified over 157,000 asteroids, including more than 500 NEOs and $\sim$120 comets. The NEOWISE dataset will enable a panoply of new scientific investigations. △ Less

Submitted 9 February, 2011; originally announced February 2011.

Comments: ApJ accepted

The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos

Authors: Tomoaki Ishiyama, Steven Rieder, Junichiro Makino, Simon Portegies Zwart, Derek Groen, Keigo Nitadori, Cees de Laat, Stephen McMillan, Kei Hiraki, Stefan Harfst

Abstract: We present the results of the "Cosmogrid" cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10… ▽ More We present the results of the "Cosmogrid" cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10^7 Msun. We have analyzed the spherically averaged density profiles of the three most massive halos which are of galaxy group size and contain at least 170 million particles. The slopes of these density profiles become shallower than -1 at the inner most radius. We also find a clear correlation of halo concentration with mass. The mass dependence of the concentration parameter cannot be expressed by a single power law, however a simple model based on the Press-Schechter theory proposed by Navarro et al. gives reasonable agreement with this dependence. The spin parameter does not show a correlation with the halo mass. The probability distribution functions for both concentration and spin are well fitted by the log-normal distribution for halos with the masses larger than ~10^8 Msun. The subhalo abundance depends on the halo mass. Galaxy-sized halos have 50% more subhalos than ~10^{11} Msun halos have. △ Less

Submitted 8 April, 2013; v1 submitted 10 January, 2011; originally announced January 2011.

Comments: 15 pages, 18 figures, accepted by ApJ

Journal ref: 2013, ApJ, 767, 146

The Wide-field Infrared Survey Explorer (WISE): Mission Description and Initial On-orbit Performance

Authors: Edward L. Wright, Peter R. M. Eisenhardt, Amy Mainzer, Michael E. Ressler, Roc M. Cutri, Thomas Jarrett, J. Davy Kirkpatrick, Deborah Padgett, Robert S. McMillan, Michael Skrutskie, S. A. Stanford, Martin Cohen, Russell G. Walker, John C. Mather, David Leisawitz, Thomas N. Gautier III, Ian McLean, Dominic Benford, Carol J. Lonsdale, Andrew Blain, Bryan Mendez, William R. Irace, Valerie Duval, Fengchuan Liu, Don Royer , et al. (13 additional authors not shown)

Abstract: The all sky surveys done by the Palomar Observatory Schmidt, the European Southern Observatory Schmidt, and the United Kingdom Schmidt, the InfraRed Astronomical Satellite and the 2 Micron All Sky Survey have proven to be extremely useful tools for astronomy with value that lasts for decades. The Wide-field Infrared Survey Explorer is mapping the whole sky following its launch on 14 December 2009.… ▽ More The all sky surveys done by the Palomar Observatory Schmidt, the European Southern Observatory Schmidt, and the United Kingdom Schmidt, the InfraRed Astronomical Satellite and the 2 Micron All Sky Survey have proven to be extremely useful tools for astronomy with value that lasts for decades. The Wide-field Infrared Survey Explorer is mapping the whole sky following its launch on 14 December 2009. WISE began surveying the sky on 14 Jan 2010 and completed its first full coverage of the sky on July 17. The survey will continue to cover the sky a second time until the cryogen is exhausted (anticipated in November 2010). WISE is achieving 5 sigma point source sensitivities better than 0.08, 0.11, 1 and 6 mJy in unconfused regions on the ecliptic in bands centered at wavelengths of 3.4, 4.6, 12 and 22 microns. Sensitivity improves toward the ecliptic poles due to denser coverage and lower zodiacal background. The angular resolution is 6.1, 6.4, 6.5 and 12.0 arc-seconds at 3.4, 4.6, 12 and 22 microns, and the astrometric precision for high SNR sources is better than 0.15 arc-seconds. △ Less

Submitted 9 November, 2010; v1 submitted 30 July, 2010; originally announced August 2010.

Comments: 22 pages with 19 included figures. Updated to better match the accepted version in the AJ

Journal ref: Edward L. Wright et al 2010 The Astronomical Journal 140 1868

The fraction of globular cluster second-generation stars in the Galactic halo

Authors: Enrico Vesperini, Stephen L. W. McMillan, Francesca D'Antona, Annibale D'Ercole

Abstract: Many observational studies have revealed the presence of multiple stellar generations in Galactic globular clusters. These studies suggest that second-generation stars make up a significant fraction of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 to 80 per cent in individual clusters. In this Letter we carry out hydrodynamical simulations to expl… ▽ More Many observational studies have revealed the presence of multiple stellar generations in Galactic globular clusters. These studies suggest that second-generation stars make up a significant fraction of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 to 80 per cent in individual clusters. In this Letter we carry out hydrodynamical simulations to explore the dependence of the mass of second-generation stars on the initial mass and structural parameters and stellar initial mass function (IMF) of the parent cluster. We then use the results of these simulations to estimate the fraction, f_{SG,H}, of the mass of the Galactic stellar halo composed of second-generation stars that originated in globular clusters. We study the dependence of f_{SG,H} on the parameters of the initial mass function of the Galactic globular cluster system. For a broad range of initial conditions, we find that the fraction of mass of the Galactic stellar halo in second-generation stars is always small, f_{SG,H}<4-6 per cent for a Kroupa-1993 IMF and f_{SG,H}<7-9 per cent for a Kroupa-2001 IMF. △ Less

Submitted 9 July, 2010; originally announced July 2010.

Comments: 5 pages, 3 figures, Accepted for publication in The Astrophysical Journal, Letters