Draft version September 1, 2018 Preprint typeset using LAT X style AASTeX6 v. 1.0 E 1 THE EFFECT OF ATMOSPHERIC COOLING ON THE VERTICAL VELOCITY DISPERSION AND DENSITY DISTRIBUTION OF BROWN DWARFS Russell E. Ryan Jr.1, Paul A. Thorman2, Sarah J. Schmidt3, Seth H. Cohen4, Nimish P. Hathi1,5, Benne W. Holwerda6,7, Jonathan I. Lunine8, Nor Pirzkal1, Rogier A. Windhorst4, Erick Young9 1Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218 2Haverford College, Departments of Physics and Astronomy, 370 Lancaster Ave., Haverford, PA 19041 3Leibniz-Institut f¨urAstrophysik | Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany 4Arizona State University, School of Earth and Space Exploration, P.O. Box 871404, Tempe, AZ 85281-1404 5 Aix Marseille Universit´, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 13388, Marseille, France 6Leiden Observatory, Leiden University, Niels Bohrweg 2, NL-2333 CA Leiden, the Netherlands 7University of Louisville, Department of Physics and Astronomy, Louisville, KY 40292 8Department of Astronomy and Carl Sagan Institute, Space Sciences Building Cornell University, Ithaca, NY 14853 9NASA Ames SOFIA Science Center, N211, Mountain View, CA 94043 ABSTRACT We present a Monte Carlo simulation designed to predict the vertical velocity dispersion of brown dwarfs in the Milky Way. We show that since these stars are constantly cooling, the velocity dispersion has a noticeable trend with spectral type. With realistic assumptions for the initial-mass function, star- formation history, and the cooling models, we show that the velocity dispersion is roughly consistent with what is observed for M dwarfs, decreases to cooler spectral types, and increases again for the coolest types in our study (∼T9). We predict a minimum in the velocity dispersions for L/T transition objects, however the detailed properties of the minimum predominately depend on the star-formation history. Since this trend is due to brown dwarf cooling, we expect the velocity dispersion as a function of spectral type should deviate from constancy around the hydrogen-burning limit. We convert from velocity dispersion to vertical scale height using standard disk models, and present similar trends in disk thickness as a function of spectral type. We suggest that future, wide-field photometric and/or spectroscopic missions may collect sizable samples of distant (∼1 kpc) of dwarfs that span the hydrogen-burning limit. As such, we speculate that such observations may provide a unique way of constraining the average spectral type of hydrogen-burning. Keywords: Galaxy: disk | Galaxy: structure | brown dwarfs | stars: low mass 1. INTRODUCTION Disk kinematics are powerful probes of the structure Brown dwarfs are substellar objects that lack the nec- and evolution of the Milky Way, but also provide con- essary mass to sustain hydrogen fusion in their cores straints on the process of star formation. It is com- (e.g. Kumar 1962; Hayashi & Nakano 1963). Hence, monly accepted that stars are formed in giant molec- ular clouds that occupy a thin Galactic disk, then dif- arXiv:1708.02591v1 [astro-ph.SR] 8 Aug 2017 these objects are constantly cooling throughout their lifetimes (e.g. Burrows, Hubbard, & Lunine 1989), which fuse over time by successive interactions with the disk leads to the known mass-age degeneracy: a brown dwarf material (e.g. Spitzer & Schwarzschild 1951). This so- detected with a specific temperature could have a range called disk heating can be used as a Galactic chronome- of ages, depending on its mass (e.g. Baraffe & Chabrier ter, which makes stellar kinematics a weak tracer of stel- 1996; Burrows et al. 1997). lar population age (e.g. Wielen 1977). Nearby popula- tions of early-M dwarfs are observed to have kinematic ages of ∼ 3 Gyr (e.g. Reid et al. 2002; Reiners & Basri 1 Support for program #13266 was provided by NASA through 2009), but extending this chronology to the substellar a grant from the Space Telescope Science Institute, which is op- population is not trivial, particularly if their formation erated by the Associations of Universities for Research in Astron- omy, Incorporated, under NASA contract NAS 5-26555. process differs from that of main-sequence stars. Luh- man(2012) review the proposed formation scenarios for [email protected] brown dwarfs and the various effects on their initial kine- 2 Ryan et al. matics and present-day spatial distributions. Early for- To characterize the effects of brown dwarf cooling on mation models speculated that brown dwarfs are pref- their Galactic properties, we develop a simple Monte erentially ejected from their birthsites (e.g. Reipurth & Carlo simulation that is similar to several other studies Clarke 2001), which should enhance their initial veloc- (e.g. Reid et al. 1999; Burgasser 2004; Day-Jones et al. ities and imply kinematic ages older than the main- 2013; Burgasser et al. 2015). We start with the usual sequence population. Most kinematic and spatial distri- assumption that the creation function is separable in butions of cluster brown dwarfs seem to rule out the ejec- mass and time: tion scenario, and favor models that have initial kine- C(m; t) dm dt = φ(m) (t) dm dt; (1) matics comparable to main-sequence stars (e.g. White & Basri 2003; Joergens 2006; Luhman et al. 2006; Parker where φ(m) is the IMF and (t) is the star-formation et al. 2011). Moreover disk brown dwarfs are reported history (SFH). Although these two distributions are of- with a range of kinematic ages, in same cases younger ten the focus of many studies, we only wish to adopt than M-dwarfs (e.g. Zapatero Osorio et al. 2007; Faherty plausible functional forms and reasonable range of pa- et al. 2009; Schmidt et al. 2010) or older (e.g. Seifahrt et rameters to set bounds on the effects of atmospheric al. 2009; Blake, Charbonneau, & White 2010; Burgasser cooling. Therefore we use a powerlaw IMF et al. 2015). Whether the lack of empirical consensus φ(m; α) dm / m−α dm; (2) stems from some observational effect (such as small sam- ple size) or points to new astrophysics (such as forma- for 0:0005 ≤ m ≤ 0:1 M , and we consider models with tion mechanisms) is not readily clear. However brown α 2 [0; 0:5; 1:0; 1:5], but discuss the role of these mass dwarf kinematics provide important insights into their limits in AppendixA. For these mass limits, our simu- formation, cooling, and disk heating. lation is valid for spectral types ∼M8{T9. The effect of atmospheric cooling on populations of We adopt a gamma-distribution form for the SFH: brown dwarfs has been discussed in great detail through- t β t dt out the literature, but there are three works that are (t; β; τ) dt / exp − (3) particularly relevant to our study. (1) Burgasser(2004) τ τ τ present a series of Monte Carlo simulations to predict for 0 ≤ t ≤ t0, where t0 is the age of the Milky Way, the present-day mass function (PDMF) of brown dwarfs. which we take to be 12 Gyr. To avoid confusion, we Using various assumptions for the initial mass function use t to refer to the time since the formation of the (IMF), birth rate, and cooling models (Burrows et al. Milky Way disk and a as the age of a brown dwarf. 1997; Baraffe et al. 2003), they confirm the relatively These two quantities are related as: a = t0 − t. We low number density of mid-L dwarfs, even for shallow consider models with τ 2 [0:3; 0:5; 0:7; 0:9; 1:1] Gyr and IMFs. (2) Deacon & Hambly(2006) develop a similar β 2 [2:5; 3:5; 4:5; 5:5; 6:5] (an exponential model being calculation to model the present-day luminosity func- the special case of β =0). The parameters β and τ rep- tion (PDLF), which have been used to rule out ex- resent the power-law slope at early times (t τ) and treme (halo-type) birth rates (Day-Jones et al. 2013) the exponential decay at late times (t∼τ), respectively. and confirm the low space density of L/T transition This parameterization of the SFH has a maximum star- dwarfs (Marocco et al. 2015). The increasing number formation rate at tmax =τβ, and it is qualitatively sim- of late-T field dwarfs has been seen in several studies ilar to that proposed by Snaith et al.(2015). Also, this (e.g. Metchev et al. 2008; Kirkpatrick et al. 2012; Burn- SFH is consistent with that of the typical Milky Way- ingham et al. 2013), which is seemingly a pile up of cool like galaxy (Papovich et al. 2015), which we estimate dwarfs. (3) Burgasser et al.(2015) propose a two-phase from their work will have (β; τ)≈(4:5; 0:7 Gyr). star-formation history to explain the old kinematic ages With these two distribution functions, the creation of (6:5 ± 0:4 Gyr) of L dwarfs in their sample of 85 M8{L6 simulated brown dwarfs is parameterized by three tun- dwarfs. Therefore the aim of our work is to build on able values (α; β; τ), which we refer to as the galaxy such kinematic models and predict the vertical number model. For each galaxy model, we draw 109 random density of brown dwarfs in the Galactic disk. masses and formation times (or ages) from Equations2 We organize this paper as follows: in Section2 we and3, respectively. describe our Monte Carlo simulation, in Section3 we We adopt the cooling models of Burrows et al.(1997) present our estimates for the disk thickness, in Section4 for solar metallicity to describe the effective tempera- we discuss our results in the context of what can be ture as a function of stellar mass and age.