BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Kenneth Wood

Kenneth Wood contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-phastro-ph.GAastro-ph.COastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2012arXiv

A self-consistent model of Galactic stellar and dust infrared emission and the abundance of polycyclic aromatic hydrocarbons

We present a self-consistent three-dimensional Monte-Carlo radiative transfer model of the stellar and dust emission in the Milky-Way, and have computed synthetic observations of the 3.6 to 100 microns emission in the Galactic mid-plane. In order to compare the model to observations, we use the GLIMPSE, MIPSGAL, and IRAS surveys to construct total emission spectra, as well as longitude and latitude profiles for the emission. The distribution of stars and dust is taken from the SKY model, and the dust emissivities includes an approximation of the emission from polycyclic aromatic hydrocarbons in addition to thermal emission. The model emission is in broad agreement with the observations, but a few modifications are needed to obtain a good fit. Firstly, by adjusting the model to include two major and two minor spiral arms rather than four equal spiral arms, the fit to the longitude profiles for |l|>30 degrees can be improved. Secondly, introducing a deficit in the dust distribution in the inner Galaxy results in a better fit to the shape of the IRAS longitude profiles at 60 and 100 microns. With these modifications, the model fits the observed profiles well, although it systematically under-estimates the 5.8 and 8.0 microns fluxes. One way to resolve this discrepancy is to increase the abundance of PAH molecules by 50% compared to the original model, although we note that changes to the dust distribution or radiation field may provide alternative solutions. Finally, we use the model to quantify which stellar populations contribute the most to the heating of different dust types, and which stellar populations and dust types contribute the most to the emission at different wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

MHD simulations of a supernova-driven ISM and the warm ionized medium using a positivity preserving ideal MHD scheme

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The 3-Dimensional Distribution of Dust in NGC 891

We produce three-dimensional Monte-Carlo radiative transfer models of the edge-on spiral galaxy NGC 891, a fast-rotating galaxy thought to be an analogue to the Milky Way. The models contain realistic spiral arms and a fractal distribution of clumpy dust. We fit our models to Hubble Space Telescope images corresponding to the B and I bands, using shapelet analysis and a genetic algorithm to generate 30 statistically best-fitting models. These models have a strong preference for spirality and clumpiness, with average face-on attenuation decreasing from 0.24(0.16) to 0.03(0.03) mag in the B(I) band between 0.5 and 2 radial scale-lengths. Most of the attenuation comes from small high-density clumps with low (<10%) filling factors. The fraction of dust in clumps is broadly consistent with results from fitting NGC 891&#39;s spectral energy distribution. Because of scattering effects and the intermixed nature of the dust and starlight, attenuation is smaller and less wavelength-dependent than the integrated dust column-density. Our clumpy models typically have higher attenuation at low inclinations than previous radiative transfer models using smooth distributions of stars and dust, but similar attenuation at inclinations above 70 degrees. At all inclinations most clumpy models have less attenuation than expected from previous estimates based on minimizing scatter in the Tully-Fisher relation. Mass-to-light ratios are higher and the intrinsic scatter in the Tully-Fisher relation is larger than previously expected for galaxies similar to NGC 891. The attenuation curve changes as a function of inclination, with R_(B,B-I)=A_(B)/E(B-I) increasing by ~0.75 from face-on to near-edge-on orientations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Evolution of brown dwarf disks: A Spitzer survey in Upper Scorpius

We have carried out a Spitzer survey for brown dwarf (BD) disks in the ~5 Myr old Upper Scorpius (UpSco) star forming region, using IRS spectroscopy from 8 to 12μm and MIPS photometry at 24μm. Our sample consists of 35 confirmed very low mass members of UpSco. Thirteen objects in this sample show clear excess flux at 24μm, explained by dust emission from a circum-sub-stellar disk. Objects without excess emission either have no disks at all or disks with inner opacity holes of at least ~5 AU radii. Our disk frequency of 37\pm 9% is higher than what has been derived previously for K0-M5 stars in the same region (on a 1.8 sigma confidence level), suggesting a mass-dependent disk lifetime in UpSco. The clear distinction between objects with and without disks as well as the lack of transition objects shows that disk dissipation inside 5 AU occurs rapidly, probably on timescales of <~10^5 years. For the objects with disks, most SEDs are uniformly flat with flux levels of a few mJy, well modeled as emission from dusty disks affected by dust settling to the midplane, which also provides indirect evidence for grain growth. The silicate feature around 10μm is either absent or weak in our SEDs, arguing for a lack of hot, small dust grains. Compared with younger objects in Taurus, BD disks in UpSco show less flaring. Taken together, these results clearly demonstrate that we see disks in an advanced evolutionary state: Dust settling and grain growth are ubiquituous in circum-sub-stellar disks at ages of 5 Myr, arguing for planet forming processes in BD disks. For almost all our targets, results from high-resolution spectroscopy and high-spatial resolution imaging have been published before, thus providing a large sample of BDs for which information about disks, accretion, and binarity is available. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Interpreting Spectral Energy Distributions from Young Stellar Objects. I. A grid of 200,000 YSO model SEDs

We present a grid of radiation transfer models of axisymmetric young stellar objects (YSOs), covering a wide range of stellar masses (from 0.1Msun to 50Msun) and evolutionary stages (from the early envelope infall stage to the late disk-only stage). The grid consists of 20,000 YSO models, with spectral energy distributions (SEDs) and polarization spectra computed at ten viewing angles for each model, resulting in a total of 200,000 SEDs. [...]. These models are publicly available on a dedicated WWW server: http://www.astro.wisc.edu/protostars/ . In this paper we summarize the main features of our models, as well as the range of parameters explored. [...]. We examine the dependence of the spectral indices of the model SEDs on envelope accretion rate and disk mass. In addition, we show variations of spectral indices with stellar temperature, disk inner radius, and disk flaring power for a subset of disk-only models. We also examine how changing the wavelength range of data used to calculate spectral indices affects their values. We show sample color-color plots of the entire grid as well as simulated clusters at various distances with typical {\it Spitzer Space Telescope} sensitivities. We find that young embedded sources generally occupy a large region of color-color space due to inclination and stellar temperature effects. Disk sources occupy a smaller region of color-color space, but overlap substantially with the region occupied by embedded sources, especially in the near- and mid-IR. We identify regions in color-color space where our models indicate that only sources at a given evolutionary stage should lie. [...].

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Interpreting Spectral Energy Distributions from Young Stellar Objects. II. Fitting observed SEDs using a large grid of pre-computed models

We present a method to analyze the spectral energy distributions (SEDs) of young stellar objects (YSOs). Our approach is to fit data with pre-computed 2-D radiation transfer models spanning a large region of parameter space. This allows us to determine not only a single set of physical parameter values but the entire range of values consistent with the multi-wavelength observations of a given source. In this way we hope to avoid any over-interpretation when modeling a set of data. We have constructed spectral energy distributions from optical to sub-mm wavelengths, including new Spitzer IRAC and MIPS photometry, for 30 young and spatially resolved sources in the Taurus-Auriga star-forming region. We demonstrate fitting model SEDs to these sources, and find that we correctly identify the evolutionary stage and physical parameters found from previous independent studies, such as disk mass, disk accretion rate, and stellar temperature. We also explore how fluxes at various wavelengths help to constrain physical parameters, and show examples of degeneracies that can occur when fitting SEDs. A web-based version of this tool is available to the community at http://www.astro.wisc.edu/protostars .

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Estimating the Porosity of the Interstellar Medium from Three-Dimensional Photoionization Modeling of H II Regions

We apply our three dimensional photoionization code to model Wisconsin H-alpha Mapper observations of the H II region surrounding the O9.5V star Zeta Oph. Our models investigate the porosity of the interstellar medium around zeta Oph and the effects of 3D densities on the H-alpha surface brightness and variation in the [N II]/H-alpha line ratio. The Zeta Oph H II region has a well characterized ionizing source, so it is an excellent starting point for 3D models of diffuse ionized gas. We investigate various hierarchically clumped density structures, varying the overall smoothness within the clumping algorithm. By simulating the observations, we can estimate the porosity of the medium in the vicinity of Zeta Oph and find that within the context of our hierarchically clumped models, around 50% to 80% of the volume is occupied by clumps surrounded by a low density smooth medium. We also conclude that in order for O stars to ionize the diffuse Warm Ionized Medium, the O star environment must be more porous than that surrounding Zeta Oph, with clumps occupying less than one half of the interstellar volume. Our clumpy models have irregular boundaries, similar to observed H II regions. However, in observed H II regions it is difficult to identify the precise location of the boundary because of the foreground and/or background emission from the widespread Warm Ionized Medium. This complicates the interpretation of the predicted rapid rise of some emission line ratios near the edge of uniform density H II regions and combined with the three dimensional clumpy nature of the interstellar medium may explain the apparent lack of distinctive emission line ratios near H I -- H II interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Modeling the Photopolarimetric Variability of AA Tau

We present Monte Carlo scattered light models of a warped disk that reproduce the observed photopolarimetric variability of the classical T Tauri star AA Tau. For a system inclination of 75 deg and using an analytic description for a warped inner disk, we find that the shape and amplitude of the photopolarimetric variability are reproduced with a warp which occults the star, located at 0.07 au, amplitude 0.016 au, extending over radial and azimuthal ranges 0.0084 au and 145 deg. We also show a time sequence of high spatial resolution scattered light images, showing a dark shadow cast by the warp sweeping round the disk. Using a modified smooth particle hydrodynamics code, we find that a stellar dipole magnetic field of strength 5.2 kG, inclined at 30 deg to the stellar rotation axis can reproduce the required disk warping to explain AA Tau&#39;s photopolarimetric variability.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Observational Properties of Proto-planetary Disk Gaps

We study the effects of an annular gap induced by an embedded proto-planet on disk scattered light images and the infrared spectral energy distribution. We find that the outer edge of a gap is brighter in the scattered light images than a similar location in a gap-free disk. The stellar radiation that would have been scattered by material within in the gap is instead scattered by the disk wall at the outer edge of the gap, producing a bright ring surrounding the dark gap in the images. Given sufficient resolution, such gaps can be detected by the presence of this bright ring in scattered light images. A gap in a disk also changes the shape of the SED. Radiation that would have been absorbed by material in the gap is instead reprocessed by the outer gap wall. This leads to a decrease in the SED at wavelengths corresponding to the temperature at the radius of the missing gap material and a corresponding flux increase at longer wavelengths corresponding to the temperature of the outer wall. We note that, unlike an inner hole in the disk, the presence of an annular gap does not change the bolometric IR flux; it simply redistributes the radiation, previously produced by material within the gap, to longer wavelengths. This implies that the changes in the SED generally will be smaller for gaps than holes. Although it will be difficult on the basis of the SED alone to distinguish between the presence of a gap and other physical effects, the level of changes can be sufficiently large to be measurable with current instruments (e.g., Spitzer).

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Structure of Brown Dwarf Circumstellar Disks

We present synthetic spectra for circumstellar disks that are heated by radiation from a central brown dwarf. Under the assumption of vertical hydrostatic equilibrium, our models yield scaleheights for brown dwarf disks in excess of three times those derived for classical T Tauri (CTTS) disks. If the near-IR excess emission observed from brown dwarfs is indeed due to circumstellar disks, then the large scaleheights we find could have a significant impact on the optical and near-IR detectability of such systems. Our radiation transfer calculations show that such highly flared disks around brown dwarfs will result in a large fraction of obscured sources due to extinction of direct starlight by the disk over a wide range of sightlines. The obscured fraction for a &#39;typical&#39; CTTS is less than 20%. We show that the obscured fraction for brown dwarfs may be double that for CTTS, but this depends on stellar and disk mass. We also comment on possible confusion in identifying brown dwarfs via color-magnitude diagrams: edge-on CTTS display similar colors and magnitudes as a face-on brown dwarf plus disk systems.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

2-D Radiative Transfer in Protostellar Envelopes: I. Effects of Geometry on Class I Sources

We present 2-D radiation transfer models of Class I Protostars and show the effect of including more realistic geometries on the resulting spectral energy distributions and images. We begin with a rotationally flattened infalling envelope as our comparison model, and add a flared disk and bipolar cavity. The disk affects the spectral energy distribution most strongly at edge-on inclinations, causing a broad dip at about 10 um (independent of the silicate feature) due to high extinction and low scattering albedo in this wavelength region. The bipolar cavities allow more direct stellar+disk radiation to emerge into polar directions, and more scattering radiation to emerge into all directions. The wavelength-integrated flux, often interpreted as luminosity, varies with viewing angle, with pole-on viewing angles seeing 2-4 times as much flux as edge-on, depending on geometry. Thus, observational estimates of luminosity should take into account the inclination of a source. The envelopes with cavities are significantly bluer in near-IR and mid-IR color-color plots than those without cavities. Using 1-D models to interpret Class I sources with bipolar cavities would lead to an underestimate of envelope mass and an overestimate of the implied evolutionary state. We compute images at near-, mid-, and far-IR wavelengths. We find that the mid-IR colors and images are sensitive to scattering albedo, and that the flared disk shadows the midplane on large size scales at all wavelengths plotted. Finally, our models produce polarization spectra which can be used to diagnose dust properties, such as albedo variations due to grain growth. Our results of polarization across the 3.1 um ice feature agree well with observations for ice mantles covering 5% of the radius of the grains.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

2-D Radiative Transfer in Protostellar Envelopes: II. An Evolutionary Sequence

We present model spectral energy distributions, colors, polarization, and images for an evolutionary sequence of a low-mass protostar from the early collapse stage (Class 0) to the remnant disk stage (Class III). We find a substantial overlap in colors and SEDs between protostars embedded in envelopes (Class 0-I) and T Tauri disks (Class II), especially at mid-IR wavelengths. Edge-on Class I-II sources show double-peaked spectral energy distributions, with a short-wavelength hump due to scattered light and the long-wavelength hump due to thermal emission. These are the bluest sources in mid-IR color-color diagrams. Since Class 0 and I sources are diffuse, the size of the aperture over which fluxes are integrated has a substantial effect on the computed colors, with larger aperture results showing significantly bluer colors. This causes overlap in color-color diagrams between all evolutionary states, especially in the mid-IR. However the near-IR polarization of the Class 0 sources is much higher than the Class I-II sources, providing a means to separate these evolutionary states. We varied the grain properties in the circumstellar envelope, allowing for larger grains in the disk midplane and smaller in the envelope. We find that grain growth in disks of Class I sources can be detected at wavelengths greater than 100 $μ$m. Our image calculations predict that the diffuse emission from edge-on Class I and II sources should be detectable in the mid-IR with the Space Infrared Telescope Facility (SIRTF) in nearby star forming regions (out to several hundred parsecs).

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Constraints on a planetary origin for the gap in GM Aurigae&#39;s protoplanetary disc

The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer

We describe a general radiative equilibrium and temperature correction procedure for use in Monte Carlo radiation transfer codes with sources of temperature-independent opacity, such as astrophysical dust. The technique utilizes the fact that Monte Carlo simulations track individual photon packets, so we may easily determine where their energy is absorbed. When a packet is absorbed, it heats a particular cell within the envelope, raising its temperature. To enforce radiative equilibrium, the absorbed packet is immediately re-emitted. To correct the cell temperature, the frequency of the re-emitted packet is chosen so that it corrects the temperature of the spectrum previously emitted by the cell. The re-emitted packet then continues being scattered, absorbed, and re-emitted until it finally escapes from the envelope. As the simulation runs, the envelope heats up, and the emergent spectral energy distribution (SED) relaxes to its equilibrium value, without iteration. This implies that the equilibrium temperature calculation requires no more computation time than the SED calculation of an equivalent pure scattering model with fixed temperature. In addition to avoiding iteration, our method conserves energy exactly, because all injected photon packets eventually escape. Furthermore, individual packets transport energy across the entire system because they are never destroyed. This long-range communication, coupled with the lack of iteration, implies that our method does not suffer the convergence problems commonly associated with lambda-iteration. To verify our temperature correction procedure, we compare our results to standard benchmark tests, and finally we present the results of simulations for two-dimensional axisymmetric density structures.

0 citations0 reviewsNo code linkedOpen access