Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “OPACITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Herbig stars' near-infrared excess: An origin in the protostellar disk's magnetically supported atmosphere

Young stars with masses 2-8 times solar, the Herbig Ae and Be stars, often show a near-infrared excess too large to explain with a hydrostatically supported circumstellar disk of gas and dust. At the same time, the accretion flow carrying the circumstellar gas to the star is thought to be driven by magnetorotational turbulence, which, according to numerical MHD modeling, yields an extended low-density atmosphere supported by the magnetic fields. We demonstrate that the base of the atmosphere can be optically thick to the starlight and that the parts lying near 1 AU are tall enough to double the fraction of the stellar luminosity reprocessed into the near-infrared. We generate synthetic spectral energy distributions (SEDs) using Monte Carlo radiative transfer calculations with opacities for submicron silicate and carbonaceous grains. The synthetic SEDs closely follow the median Herbig SED constructed recently by Mulders and Dominik and, in particular, match the large near-infrared flux, provided the grains have a mass fraction close to interstellar near the disk's inner rim.

79 ASTRONOMY AND ASTROPHYSICS↗

The emergent 1.1-1.7 μm spectrum of the exoplanet COROT-2B as measured using the Hubble space telescope

We have used Hubble/WFC3 and the G141 grism to measure the secondary eclipse of the transiting, very hot Jupiter CoRoT-2b in the 1.1-1.7 μm spectral region. We find an eclipse depth averaged over this band equal to 395{sub −45}{sup +69} parts per million, equivalent to a blackbody temperature of 1788 ± 18 K. We study and characterize several WFC3 instrumental effects, especially the 'hook' phenomenon described by Deming et al. We use data from several transiting exoplanet systems to find a quantitative relation between the amplitude of the hook and the exposure level of a given pixel. Although the uncertainties in this relation are too large to allow us to develop an empirical correction for our data, our study provides a useful guide for optimizing exposure levels in future WFC3 observations. We derive the planet's spectrum using a differential method. The planet-to-star contrast increases to longer wavelength within the WFC3 bandpass, but without water absorption or emission to a 3σ limit of 85 ppm. The slope of the WFC3 spectrum is significantly less than the slope of the best-fit blackbody. We compare all existing eclipse data for this planet to a blackbody spectrum, and to spectra from both solar abundance and carbon-rich (C/O = 1) models. A blackbody spectrum is an acceptable fit to the full data set. Extra continuous opacity due to clouds or haze, and flattened temperature profiles, are strong candidates to produce quasi-blackbody spectra, and to account for the amplitude of the optical eclipses. Our results show ambiguous evidence for a temperature inversion in this planet.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical evolution of viscous disks around be stars. II. Polarimetry

Be stars exhibit variability for a great number of observables. Putting the pieces of the disk dynamics together is not an easy task and requires arduous modeling before achieving a good fit to the observational data. In order to guide the modeling process and make it more efficient, it is very instructive to investigate reference dynamical cases. This paper focuses on continuum polarimetric quantities and is the second of a series that aims to demonstrate the capacity of deriving the dynamical history and fundamental parameters of a classical Be star through follow-up of various observables. After a detailed study of the different opacities at play in the formation of polarized spectra, we investigate predictions of polarimetric observables in the continuum for different dynamical scenarios. Our models are based on a coupling of a hydrodynamic viscous decretion simulations in a disk and a three-dimensional non-LTE radiative transfer code. Through introduction of the polarization color diagram (PCD), we show that certain combinations of polarimetric observables exhibit features that are characteristic of a mass-loss history. This diagram also enables estimates of fundamental parameters such as the inclination angle, disk density scale, and the α viscous diffusion parameter. We present the PCD as a powerful diagnosis tool to track the dynamical phases of a Be star, such as disk build-up, dissipation, periodic, and episodic outbursts. Finally, we confront our models with observations of four Be stars that exhibited long-term polarimetric activity.

79 ASTRONOMY AND ASTROPHYSICS↗

A monitoring campaign for Luhman 16AB. I. Detection of resolved near-infrared spectroscopic variability

We report resolved near-infrared spectroscopic monitoring of the nearby L dwarf/T dwarf binary WISE J104915.57–531906.1AB (Luhman 16AB), as part of a broader campaign to characterize the spectral energy distribution and temporal variability of this system. A continuous 45 minute sequence of low-resolution IRTF/SpeX data spanning 0.8-2.4 μm were obtained, concurrent with combined-light optical photometry with ESO/TRAPPIST. Our spectral observations confirm the flux reversal of this binary, and we detect a wavelength-dependent decline in the relative spectral fluxes of the two components coincident with a decline in the combined-light optical brightness of the system over the course of the observation. These data are successfully modeled as a combination of achromatic (brightness) and chromatic (color) variability in the T0.5 Luhman 16B, consistent with variations in overall cloud opacity; and no significant variability was found in L7.5 Luhman 16A, consistent with recent resolved photometric monitoring. We estimate a peak-to-peak amplitude of 13.5% at 1.25 μm over the full light curve. Using a simple two-spot brightness temperature model for Luhman 16B, we infer an average cold covering fraction of ≈30%-55%, varying by 15%-30% over a rotation period assuming a ≈200-400 K difference between hot and cold regions. We interpret these variations as changes in the covering fraction of a high cloud deck and corresponding 'holes' which expose deeper, hotter cloud layers, although other physical interpretations are possible. A Rhines scale interpretation for the size of the variable features explains an apparent correlation between period and amplitude for Luhman 16B and the variable T dwarfs SIMP 0136+0933 and 2MASS J2139+0220, and predicts relatively fast winds (1-3 km s{sup –1}) for Luhman 16B consistent with light curve evolution on an advective time scale (1-3 rotation periods). The strong variability observed in this flux reversal brown dwarf pair supports the model of a patchy disruption of the mineral cloud layer as a universal feature of the L dwarf/T dwarf transition.

79 ASTRONOMY AND ASTROPHYSICS↗

Temperature fluctuations as a source of brown dwarf variability

A number of brown dwarfs are now known to be variable with observed amplitudes as large as 10%-30% at some wavelengths. While spatial inhomogeneities in cloud coverage and thickness are likely responsible for much of the observed variability, it is possible that some of the variations arise from atmospheric temperature fluctuations instead of, or in addition to, clouds. To better understand the role that thermal variability might play we present a case study of brown dwarf variability using a newly developed one-dimensional, time-stepping model of atmospheric thermal structure. We focus on the effects of thermal perturbations, intentionally simplifying the problem through omission of clouds and atmospheric circulation. Model results demonstrate that thermal perturbations occurring deep in the atmosphere (at pressures greater than 10 bar) of a model T-dwarf can be communicated to the upper atmosphere through radiative heating via the windows in near-infrared water opacity. The response time depends on where in the atmosphere a thermal perturbation is introduced. We show that, for certain periodic perturbations, the emission spectrum can have complex time- and wavelength-dependent behaviors, including phase shifts in times of maximum flux observed at different wavelengths. Since different wavelengths probe different levels in the atmosphere, these variations track a wavelength-dependent set of radiative exchanges happening between different atmospheric levels as a perturbation evolves in time. We conclude that thermal—as well as cloud—fluctuations must be considered as possible contributors to the observed brown dwarf variability.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiation magnetohydrodynamic simulations of protostellar collapse: Low-metallicity environments

Among many physical processes involved in star formation, radiation transfer is one of the key processes because it dominantly controls the thermodynamics. Because metallicities control opacities, they are one of the important environmental parameters that affect star formation processes. In this work, I investigate protostellar collapse in solar-metallicity and low-metallicity (Z = 0.1 Z {sub ☉}) environments using three-dimensional radiation hydrodynamic and magnetohydrodynamic simulations. Because radiation cooling in high-density gas is more effective in low-metallicity environments, first cores are colder and have lower entropies. As a result, first cores are smaller, less massive, and have shorter lifetimes in low-metallicity clouds. Therefore, first cores would be less likely to be found in low-metallicity star forming clouds. This also implies that first cores tend to be more gravitationally unstable and susceptible to fragmentation. The evolution and structure of protostellar cores formed after the second collapse weakly depend on metallicities in the spherical and magnetized models, despite the large difference in the metallicities. Because this is due to the change of the heat capacity by dissociation and ionization of hydrogen, it is a general consequence of the second collapse as long as the effects of radiation cooling are not very large during the second collapse. On the other hand, the effects of different metallicities are more significant in the rotating models without magnetic fields, because they evolve slower than other models and therefore are more affected by radiation cooling.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultraviolet observations of Super-Chandrasekhar mass type Ia supernova candidates with swift UVOT

Among Type Ia supernovae (SNe Ia), a class of overluminous objects exist whose ejecta mass is inferred to be larger than the canonical Chandrasekhar mass. We present and discuss the UV/optical photometric light curves, colors, absolute magnitudes, and spectra of three candidate Super-Chandrasekhar mass SNe—2009dc, 2011aa, and 2012dn—observed with the Swift Ultraviolet/Optical Telescope. The light curves are at the broad end for SNe Ia, with the light curves of SN 2011aa being among the broadest ever observed. We find all three to have very blue colors which may provide a means of excluding these overluminous SNe from cosmological analysis, though there is some overlap with the bluest of 'normal' SNe Ia. All three are overluminous in their UV absolute magnitudes compared to normal and broad SNe Ia, but SNe 2011aa and 2012dn are not optically overluminous compared to normal SNe Ia. The integrated luminosity curves of SNe 2011aa and 2012dn in the UVOT range (1600-6000 Å) are only half as bright as SN 2009dc, implying a smaller {sup 56}Ni yield. While it is not enough to strongly affect the bolometric flux, the early time mid-UV flux makes a significant contribution at early times. The strong spectral features in the mid-UV spectra of SNe 2009dc and 2012dn suggest a higher temperature and lower opacity to be the cause of the UV excess rather than a hot, smooth blackbody from shock interaction. Further work is needed to determine the ejecta and {sup 56}Ni masses of SNe 2011aa and 2012dn and to fully explain their high UV luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗

Searching for circumplanetary disks around LkCa 15

We present Karl G. Jansky Very Large Array (VLA) observations of the 7 mm continuum emission from the disk surrounding the young star LkCa 15. The observations achieve an angular resolution of 70 mas and spatially resolve the circumstellar emission on a spatial scale of 9 AU. The continuum emission traces a dusty annulus of 45 AU in radius that is consistent with the dust morphology observed at shorter wavelengths. The VLA observations also reveal a compact source at the center of the disk, possibly due to thermal emission from hot dust or ionized gas located within a few AU from the central star. No emission is observed between the star and the dusty ring and, in particular, at the position of the candidate protoplanet LkCa 15 b. By comparing the observations with theoretical models for circumplanetary disk emission, we find that if LkCa 15 b is a massive planet (>5 M{sub J} ) accreting at a rate greater than 10{sup –6} M{sub J} yr{sup –1}, then its circumplanetary disk is less massive than 0.1 M{sub J} , or smaller than 0.4 Hill radii. Similar constraints are derived for any possible circumplanetary disk orbiting within 45 AU from the central star. The mass estimates are uncertain by at least one order of magnitude due to the uncertainties on the mass opacity. Future ALMA observations of this system might be able to detect circumplanetary disks down to a mass of 5 × 10{sup –4} M{sub J} and as small as 0.2 AU, providing crucial constraints on the presence of giant planets in the act of forming around this young star.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconstructing the stellar mass distributions of galaxies using S{sup 4}G IRAC 3.6 and 4.5 μm images. II. The conversion from light to mass

We present a new approach for estimating the 3.6 μm stellar mass-to-light (M/L) ratio Y{sub 3.6} in terms of the [3.6]-[4.5] colors of old stellar populations. Our approach avoids several of the largest sources of uncertainty in existing techniques using population synthesis models. By focusing on mid-IR wavelengths, we gain a virtually dust extinction-free tracer of the old stars, avoiding the need to adopt a dust model to correctly interpret optical or optical/near-IR colors normally leveraged to assign the mass-to-light ratio Y. By calibrating a new relation between near-IR and mid-IR colors of giant stars observed in GLIMPSE we also avoid the discrepancies in model predictions for the [3.6]-[4.5] colors of old stellar populations due to uncertainties in the molecular line opacities assumed in template spectra. We find that the [3.6]-[4.5] color, which is driven primarily by metallicity, provides a tight constraint on Y{sub 3.6}, which varies intrinsically less than at optical wavelengths. The uncertainty on Y{sub 3.6} of ∼0.07 dex due to unconstrained age variations marks a significant improvement on existing techniques for estimating the stellar M/L with shorter wavelength data. A single Y{sub 3.6} = 0.6 (assuming a Chabrier initial mass function (IMF)), independent of [3.6]-[4.5] color, is also feasible because it can be applied simultaneously to old, metal-rich and young, metal-poor populations, and still with comparable (or better) accuracy (∼0.1 dex) than alternatives. We expect our Y{sub 3.6} to be optimal for mapping the stellar mass distributions in S{sup 4}G galaxies, for which we have developed an independent component analysis technique to first isolate the old stellar light at 3.6 μm from nonstellar emission (e.g., hot dust and the 3.3 polycyclic aromatic hydrocarbon feature). Our estimate can also be used to determine the fractional contribution of nonstellar emission to global (rest-frame) 3.6 μm fluxes, e.g., in WISE imaging, and establishes a reliable basis for exploring variations in the stellar IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

A three-dimensional view of the thermal structure in a super-penumbral canopy

We investigate the three-dimensional (3D) thermal topology in a super-penumbral canopy of an active region (AR). We derive temperature stratifications in the AR by an inversion of the Ca II IR line at 854.2 nm, assuming local thermal equilibrium. We find that about half of the radially oriented fibrils in the super-penumbral canopy form short, low-lying (h < 1 Mm) loops in the 3D temperature cube. These closed loops connect from bright grains in or close to the penumbra to the photosphere a few Mms away from the sunspot. The other half of the fibrils monotonically rise with distance from the sunspot. Many of the fibrils show a central dark core and two lateral brightenings in line-core intensity images. The corresponding velocity image shows fibrils that are as wide as the fibrils seen in intensity without a lateral substructure. Additionally, we study a feature from a different class of structures without prominent mass flows. Its 3D topology is formed by two parallel, closed loops that connect patches of opposite polarity. We present evidence that the inverse Evershed flow into the sunspot in the lower chromosphere is the consequence of siphon flows along short loops that connect photospheric foot points. The dark-cored structure of the chromospheric fibrils cannot have a convective origin because of their location above regular granulation. The dark core most likely results from an opacity difference between the central axis and the lateral edges caused by the significant flow speed along the fibrils.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure and dynamics of the accretion process and wind in TW Hya

Time-domain spectroscopy of the classical accreting T Tauri star, TW Hya, covering a decade and spanning the far UV to the near-infrared spectral regions can identify the radiation sources, the atmospheric structure produced by accretion, and properties of the stellar wind. On timescales from days to years, substantial changes occur in emission line profiles and line strengths. Our extensive time-domain spectroscopy suggests that the broad near-IR, optical, and far-uv emission lines, centered on the star, originate in a turbulent post-shock region and can undergo scattering by the overlying stellar wind as well as some absorption from infalling material. Stable absorption features appear in Hα, apparently caused by an accreting column silhouetted in the stellar wind. Inflow of material onto the star is revealed by the near-IR He I 10830 Å line, and its free-fall velocity correlates inversely with the strength of the post-shock emission, consistent with a dipole accretion model. However, the predictions of hydrogen line profiles based on accretion stream models are not well-matched by these observations. Evidence of an accelerating warm to hot stellar wind is shown by the near-IR He I line, and emission profiles of C II, C III, C IV, N V, and O VI. The outflow of material changes substantially in both speed and opacity in the yearly sampling of the near-IR He I line over a decade. Terminal outflow velocities that range from 200 km s{sup –1} to almost 400 km s{sup –1} in He I appear to be directly related to the amount of post-shock emission, giving evidence for an accretion-driven stellar wind. Calculations of the emission from realistic post-shock regions are needed.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating three-dimensional nonthermal high-energy photon emission in colliding-wind binaries

Massive stars in binary systems have long been regarded as potential sources of high-energy γ rays. The emission is principally thought to arise in the region where the stellar winds collide and accelerate relativistic particles which subsequently emit γ rays. On the basis of a three-dimensional distribution function of high-energy particles in the wind collision region—as obtained by a numerical hydrodynamics and particle transport model—we present the computation of the three-dimensional nonthermal photon emission for a given line of sight. Anisotropic inverse Compton emission is modeled using the target radiation field of both stars. Photons from relativistic bremsstrahlung and neutral pion decay are computed on the basis of local wind plasma densities. We also consider photon-photon opacity effects due to the dense radiation fields of the stars. Results are shown for different stellar separations of a given binary system comprising of a B star and a Wolf-Rayet star. The influence of orbital orientation with respect to the line of sight is also studied by using different orbital viewing angles. For the chosen electron-proton injection ratio of 10{sup –2}, we present the ensuing photon emission in terms of two-dimensional projections maps, spectral energy distributions, and integrated photon flux values in various energy bands. Here, we find a transition from hadron-dominated to lepton-dominated high-energy emission with increasing stellar separations. In addition, we confirm findings from previous analytic modeling that the spectral energy distribution varies significantly with orbital orientation.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining the location of gamma-ray flares in luminous blazars

Locating the gamma-ray emission sites in blazar jets is a long standing and highly controversial issue. We jointly investigate several constraints on the distance scale r and Lorentz factor Γ of the gamma-ray emitting regions in luminous blazars (primarily flat spectrum radio quasars). Working in the framework of one-zone external radiation Comptonization models, we perform a parameter space study for several representative cases of actual gamma-ray flares in their multiwavelength context. We find a particularly useful combination of three constraints: from an upper limit on the collimation parameter Γθ ≲ 1, from an upper limit on the synchrotron self-Compton (SSC) luminosity L{sub SSC} ≲ L{sub X}, and from an upper limit on the efficient cooling photon energy E{sub cool,obs} ≲ 100 MeV. These three constraints are particularly strong for sources with low accretion disk luminosity L{sub d}. The commonly used intrinsic pair-production opacity constraint on Γ is usually much weaker than the SSC constraint. The SSC and cooling constraints provide a robust lower limit on the collimation parameter Γθ ≳ 0.1-0.7. Typical values of r corresponding to moderate values of Γ ∼ 20 are in the range 0.1-1 pc, and are determined primarily by the observed variability timescale t{sub var,obs}. Alternative scenarios motivated by the observed gamma-ray/millimeter connection, in which gamma-ray flares of t{sub var,obs} ∼ a few days are located at r ∼ 10 pc, are in conflict with both the SSC and cooling constraints. Moreover, we use a simple light travel time argument to point out that the gamma-ray/millimeter connection does not provide a significant constraint on the location of gamma-ray flares. We argue that spine-sheath models of the jet structure do not offer a plausible alternative to external radiation fields at large distances; however, an extended broad-line region is an idea worth exploring. We propose that the most definite additional constraint could be provided by determination of the synchrotron self-absorption frequency for correlated synchrotron and gamma-ray flares.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic pair creation transparency in gamma-ray pulsars

Magnetic pair creation, γ → e {sup +} e {sup –}, has been at the core of radio pulsar paradigms and central to polar cap models of gamma-ray pulsars for over three decades. The Fermi gamma-ray pulsar population now exceeds 140 sources and has defined an important part of Fermi's science legacy, providing rich information for the interpretation of young energetic pulsars and old millisecond pulsars. Among the population characteristics well established is the common occurrence of exponential turnovers in their spectra in the 1-10 GeV range. These turnovers are too gradual to arise from magnetic pair creation in the strong magnetic fields of pulsar inner magnetospheres. By demanding insignificant photon attenuation precipitated by such single-photon pair creation, the energies of these turnovers for Fermi pulsars can be used to compute lower bounds for the typical altitude of GeV band emission. This paper explores such pair transparency constraints below the turnover energy and updates earlier altitude bound determinations that have been deployed in various Fermi pulsar papers. For low altitude emission locales, general relativistic influences are found to be important, increasing cumulative opacity, shortening the photon attenuation lengths, and also reducing the maximum energy that permits escape of photons from a neutron star magnetosphere. Rotational aberration influences are also explored, and are found to be small at low altitudes, except near the magnetic pole. The analysis presented in this paper clearly demonstrates that including near-threshold physics in the pair creation rate is essential to deriving accurate attenuation lengths and escape energies. The altitude bounds are typically in the range of 2-7 stellar radii for the young Fermi pulsar population, and provide key information on the emission altitude in radio quiet pulsars that do not possess double-peaked pulse profiles. The bound for the Crab pulsar is at a much higher altitude, with the putative detection by MAGIC out to 350-400 GeV implying a lower bound of 310 km to the emission region, i.e., approximately 20% of the light cylinder radius. These results are also extended to the super-critical field domain, where it is found that emission in neurton stars originating below around 10 stellar radii will not appear in the Fermi-LAT band.

79 ASTRONOMY AND ASTROPHYSICS↗

Optical observations of the type Ic supernova 2007gr in NGC 1058

We present extensive optical observations of the normal Type Ic supernova (SN) 2007gr, spanning from about one week before maximum light to more than one year thereafter. The optical light and color curves of SN 2007gr are very similar to those of the broad-lined Type Ic SN 2002ap, but the spectra show remarkable differences. The optical spectra of SN 2007gr are characterized by unusually narrow lines, prominent carbon lines, and slow evolution of the line velocity after maximum light. The earliest spectrum (taken at t = –8 days) shows a possible signature of helium (He I λ5876 at a velocity of ∼19,000 km s{sup –1}). Moreover, the larger intensity ratio of the [O I] λ6300 and λ6364 lines inferred from the early nebular spectra implies a lower opacity of the ejecta shortly after the explosion. These results indicate that SN 2007gr perhaps underwent a less energetic explosion of a smaller-mass Wolf-Rayet star (∼8-9 M{sub ☉}) in a binary system, as favored by an analysis of the progenitor environment through pre-explosion and post-explosion Hubble Space Telescope images. In the nebular spectra, asymmetric double-peaked profiles can be seen in the [O I] λ6300 and Mg I] λ4571 lines. We suggest that the two peaks are contributed by the blueshifted and rest-frame components. The similarity in velocity structure and the different evolution of the strength of the two components favor an aspherical explosion with the ejecta distributed in a torus or disk-like geometry, but inside the ejecta the O and Mg have different distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the Alfvénic nature of the interstellar medium: Velocity anisotropy revisited

The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel and Lazarian method to estimate the Alfvén Mach number using the structure function anisotropy in velocity centroid data from Position-Position-Velocity maps. We utilize three-dimensional magnetohydrodynamic simulations of fully developed turbulence, with a large range of sonic and Alfvénic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfvén Mach number dependency found in Esquivel and Lazarian might change when taking the second moment of the Position-Position-Velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e., Alfvén Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line of sight is parallel to up to ≈45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvénic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques.

79 ASTRONOMY AND ASTROPHYSICS↗

New evidence for a substellar luminosity problem: Dynamical mass for the brown dwarf binary Gl 417BC

We present new evidence for a problem with cooling rates predicted by substellar evolutionary models that implies that model-derived masses in the literature for brown dwarfs and directly imaged planets may be too high. Based on our dynamical mass for Gl 417BC (L4.5+L6) and a gyrochronology system age from its young, solar-type host star, commonly used models predict luminosities 0.2-0.4 dex lower than we observe. This corroborates a similar luminosity-age discrepancy identified in our previous work on the L4+L4 binary HD 130948BC, which coincidentally has nearly identical component masses (≈50-55 M{sub Jup}) and age (≈800 Myr) as Gl 417BC. Such a luminosity offset would cause systematic errors of 15%-25% in model-derived masses at this age. After comparing different models, including cloudless models that should not be appropriate for mid-L dwarfs like Gl 417BC and HD 130948BC but actually match their luminosities better, we speculate the observed overluminosity could be caused by opacity holes (i.e., patchy clouds) in these objects. Moreover, from hybrid substellar evolutionary models that account for cloud disappearance, we infer the corresponding phase of overluminosity may extend from a few hundred million years up to a few gigayears and cause masses to be overestimated by up to 25%, even well after clouds disappear from view entirely. Thus, the range of ages and spectral types affected by this potential systematic shift in luminosity evolution would encompass most known directly imaged gas-giants and field brown dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

Brown dwarf disks with ALMA

We present Atacama Large Millimeter/submillimeter Array continuum and spectral line data at 0.89 mm and 3.2 mm for three disks surrounding young brown dwarfs and very low mass stars in the Taurus star forming region. Dust thermal emission is detected and spatially resolved for all the three disks, while CO(J = 3-2) emission is seen in two disks. We analyze the continuum visibilities and constrain the disks' physical structure in dust. The results of our analysis show that the disks are relatively large; the smallest one has an outer radius of about 70 AU. The inferred disk radii, radial profiles of the dust surface density, and disk to central object mass ratios lie within the ranges found for disks around more massive young stars. We derive from our observations the wavelength dependence of the millimeter dust opacity. In all the three disks, data are consistent with the presence of grains with at least millimeter sizes, as also found for disks around young stars, and confirm that the early stages of the solid growth toward planetesimals occur also around very low-mass objects. We discuss the implications of our findings on models of solids evolution in protoplanetary disks, the main mechanisms proposed for the formation of brown dwarfs and very low-mass stars, as well as the potential of finding rocky and giant planets around very low-mass objects.

79 ASTRONOMY AND ASTROPHYSICS↗