BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kenneth Wood

Kenneth Wood appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-phastro-ph.GAastro-ph.SRastro-ph.COastro-ph.EPastro-ph.HEastro-ph.IMphysics.ao-phphysics.space-ph

Catalog footprint

What is connected

33works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

33 published item(s)

preprint2016arXiv

A Model for (Quasi-)Periodic Multi-wavelength Photometric Variability in Young Stellar Objects

We present radiation transfer models of rotating young stellar objects (YSOs) with hotspots in their atmospheres, inner disk warps and other 3-D effects in the nearby circumstellar environment. Our models are based on the geometry expected from the magneto-accretion theory, where material moving inward in the disk flows along magnetic field lines to the star and creates stellar hotspots upon impact. Due to rotation of the star and magnetosphere, the disk is variably illuminated. We compare our model light curves to data from the Spitzer YSOVAR project (Morales-Calderon et al. 2014, Cody et al. 2014) to determine if these processes can explain the variability observed at optical and mid-infrared wavelengths in young stars. We focus on those variables exhibiting "dipper" behavior that may be periodic, quasi-periodic, or aperiodic. We find that the stellar hotspot size and temperature affects the optical and near-infrared light curves, while the shape and vertical extent of the inner disk warp affects the mid-IR light curve variations. Clumpy disk distributions with non-uniform fractal density structure produce more stochastic light curves. We conclude that the magneto-accretion theory is consistent with certain aspects of the multi-wavelength photometric variability exhibited by low-mass YSOs. More detailed modeling of individual sources can be used to better determine the stellar hotspot and inner disk geometries of particular sources.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ionisation and discharge in cloud-forming atmospheres of brown dwarfs and extrasolar planets

Brown dwarfs and giant gas extrasolar planets have cold atmospheres with a rich chemical compositions from which mineral cloud particles form. Their properties, like particle sizes and material composition, vary with height, and the mineral cloud particles are charged due to triboelectric processes in such dynamic atmospheres. The dynamics of the atmospheric gas is driven by the irradiating host star and/or by the rotation of the objects that changes during its lifetime. Thermal gas ionisation in these ultra-cool but dense atmospheres allows electrostatic interactions and magnetic coupling of a substantial atmosphere volume. Combined with a strong magnetic field $\gg B_{\rm Earth}$, a chromosphere and aurorae might form as suggested by radio and X-ray observations of brown dwarfs. Non-equilibrium processes like cosmic ray ionisation and discharge processes in clouds will increase the local pool of free electrons in the gas. Cosmic rays and lighting discharges also alter the composition of the local atmospheric gas such that tracer molecules might be identified. Cosmic rays affect the atmosphere through air showers which was modelled with a 3D Monte Carlo radiative transfer code to be able to visualise their spacial extent. Given a certain degree of thermal ionisation of the atmospheric gas, we suggest that electron attachment to charge mineral cloud particles is too inefficient to cause an electrostatic disruption of the cloud particles. Cloud particles will therefore not be destroyed by Coulomb explosion for the local temperature in the collisional dominated brown dwarf and giant gas planet atmospheres. However, the cloud particles are destroyed electrostatically in regions with strong gas ionisation. The potential size of such cloud holes would, however, be too small and might occur too far inside the cloud to mimic the effect of, e.g., magnetic field induced star spots.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

CSI 2264: Simultaneous optical and infrared light curves of young disk-bearing stars in NGC 2264 with CoRoT and Spitzer-- evidence for multiple origins of variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$α$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Models of Diffuse Hα in the Interstellar Medium: The Relative Contributions from In Situ Ionisation and Dust Scattering

Using three dimensional Monte Carlo radiation transfer models of photoionisation and dust scattering, we explore different components of the widespread diffuse H$α$ emission observed in the interstellar medium of the Milky Way and other galaxies. We investigate the relative contributions of H$α$ from recombination emission in ionised gas and H$α$ that originates in H{\sc ii} regions near the Galactic midplane and scatters off high altitude dust in the diffuse interstellar medium. For the radiation transfer simulations we consider two geometries for the interstellar medium: a three dimensional fractal geometry that reproduces the average density structure inferred for hydrogen in the Milky Way, and a density structure from a magneto hydrodynamic simulation of a supernovae driven turbulent interstellar medium. Although some sight lines that are close to H{\sc ii} regions can be dominated by scattered light, overall we find that less than $\sim 20\%$ of the total H$α$ intensity in our simulations can be attributed to dust scattering. Our findings on the relative contribution of scattered H$α$ are consistent with previous observational and theoretical analyses. We also investigate the relative contributions of dust scattering and in situ ionisation of high density dust clouds in the diffuse gas. Dust scattering in these partially ionised clouds contribute $\sim 40\%$ to the total intensity of H$α$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Photoionisation and Heating of a Supernova Driven, Turbulent, Interstellar Medium

The Diffuse Ionised Gas (DIG) in galaxies traces photoionisation feedback from massive stars. Through three dimensional photoionisation simulations, we study the propagation of ionising photons, photoionisation heating and the resulting distribution of ionised and neutral gas within snapshots of magnetohydrodynamic simulations of a supernova driven turbulent interstellar medium. We also investigate the impact of non-photoionisation heating on observed optical emission line ratios. Inclusion of a heating term which scales less steeply with electron density than photoionisation is required to produce diagnostic emission line ratios similar to those observed with the Wisconsin Hα Mapper. Once such heating terms have been included, we are also able to produce temperatures similar to those inferred from observations of the DIG, with temperatures increasing to above 15000 K at heights |z| > 1 kpc. We find that ionising photons travel through low density regions close to the midplane of the simulations, while travelling through diffuse low density regions at large heights. The majority of photons travel small distances (< 100pc); however some travel kiloparsecs and ionise the DIG.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A systematic survey for eruptive young stellar objects using mid-infrared photometry

Accretion in young stellar objects (YSOs) is at least partially episodic, i.e. periods with high accretion rates ('bursts') are interspersed by quiescent phases. These bursts manifest themselves as eruptive variability. Here we present a systematic survey for eruptive YSOs aiming to constrain the frequency of accretion bursts. We compare mid-infrared photometry from Spitzer and WISE separated by ~5 yr for two samples of YSOs, in nearby star forming regions and in the Galactic plane, each comprising about 4000 young sources. All objects for which the brightness at 3.6 and 4.5 mu is increased by at least 1 mag between the two epochs may be eruptive variables and burst candidates. For these objects, we carry out follow-up observations in the near-infrared. We discover two new eruptive variables in the Galactic plane which could be FU Ori-type objects, with K-band amplitudes of more than 1.5 mag. One object known to undergo an accretion burst, V2492 Cyg, is recovered by our search as well. In addition, the young star ISO-Oph-50, previously suspected to be an eruptive object, is found to be better explained by a disk with varying circumstellar obscuration. In total, the number of burst events in a sample of 4000 YSOs is 1-4. Assuming that all YSOs undergo episodic accretion, this constraint can be used to show that phases of strong accretion (>10^-6 Msol/yr) occur in intervals of about 10^4 yr, most likely between 5000 and 50000 yr. This is consistent with the dynamical timescales for outflows, but not with the separations of emission knots in outflows, indicating that episodic accretion could either trigger or stop collimated large-scale outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three dimensional geometries and the analysis of H II regions

We compare emission line intensities from photoionization models of smooth and fractal shell geometries for low density H II regions, with particular focus on the low-ionization diagnostic diagram [N II]/H-alpha vs H-alpha. Building on previously published models and observations of Barnard's Loop, we show that the observed range of intensities and variations in the line intensity ratios may be reproduced with a three dimensional shell geometry. Our models adopt solar abundances throughout the model nebula, in contrast with previous one dimensional modeling which suggested the variations in line intensity ratios could only be reproduced if the heavy element abundances were increased by a factor of 1.4. For spatially resolved H II regions, the multiple sightlines that pierce and sample different ionization and temperature conditions within smooth and fractal shells produce a range of line intensities that are easily overlooked if only the total integrated intensities from the entire nebula model are computed. Our conclusion is that inference of H II region properties, such as elemental abundances, via photoionization models of one dimensional geometries must be treated with caution and further tested through three dimensional modeling.

0 citations0 reviewsNo code linkedOpen access
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
preprint2012arXiv

Probing the Structure and Kinematics of the transition Layer between the Magellanic Stream and the Halo in HI

The Magellanic Stream (MS) is a nearby laboratory for studying the fate of cool gas streams injected into a gaseous galactic halo. We investigate properties of the boundary layer between the cool MS gas and the hot Milky Way halo with 21 cm HI observations of a relatively isolated cloud having circular projection in the northern MS. Through averaging and modeling techniques, our observations obtained with the Robert C. Byrd Green Bank Telescope (GBT), reach unprecedented 3σ sensitivity of ~10^17/cm^2, while retaining the telescope's 9.1' resolution in the essential radial dimension. We find an envelope of diffuse neutral gas with FWHM of 60 km/s, associated in velocity with the cloud core having FWHM of 20 km/s, extending to 3.5 times the core radius with a neutral mass seven times that of the core. We show that the envelope is too extended to represent a conduction-dominated layer between the core and the halo. Its observed properties are better explained by a turbulent mixing layer driven by hydrodynamic instabilities. The fortuitous alignment of the NGC 7469 background source near the cloud center allows us to combine UV absorption and HI emission data to determine a core temperature of 8350 +/- 50 K. We show that the HI column density and size of the core can be reproduced when a slightly larger cloud is exposed to Galactic and extragalactic background ionizing radiation. Cooling in the large diffuse turbulent mixing layer envelope extends the cloud lifetime by at least a factor of two relative to a simple hydrodynamic ablation case, suggesting that the cloud is likely to reach the Milky Way disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The 3-Dimensional Structure of NGC 891 and M51

We investigate the three-dimensional structure of the nearby edge-on spiral galaxy NGC 891 using 3D Monte Carlo radiative transfer models, with realistic spiral structure and fractally clumped dust. Using the spiral and clumpiness parameters found from recently completed scattered light models we produce lower resolution SED models which reproduce the global UV-to-FIR SED of NGC 891. Our models contain a color gradient across the major axis of the galaxy - similar to what is seen in images of the NGC 891. With minor adjustment our SED models are able to match the majority of M51's SED, a similar galaxy at a near face-on different inclination.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Infrared Spectroscopy of Halos of Edge-on Galaxies

We present a study of ionized gas, PAHs, and molecular hydrogen emission in the halos of three edge-on galaxies, NGC 891, NGC 5775 and NGC 3044, based on 10-20 micron Spitzer Space Telescope spectra. The [Ne III]/[Ne II] ratio, an excellent measure of radiation hardness, rises with z in the halo of NGC 891. It is also higher in the halo of NGC 5775 than in the disk. NGC 3044 presents a more confusing situation. To explain the [Ne III]/[Ne II] as well as optical line ratio behavior in NGC 891, we carry out a simple exploration of parameter space with CLOUDY, which indicates a large increase in radiation temperature with height. Illustrative examples of physical models using a Monte Carlo radiative transfer code show that the rising neon ratio may be explained by adding a vertically extended, hot stellar source to a thin disk of massive stars. However, several other sources of hard spectra may be relevant. PAH features have scale heights of 430--530 pc in NGC 891 and 720--1080 pc in NGC 5775, suggesting they can be transported by disk-halo flows. Within NGC 891 and NGC 5775, scale heights are similar for all PAHs. For NGC 891, the scale heights exceed that of 8 micron emission, indicating a transition from more ionized to more neutral PAHs with height. Most PAH equivalent widths are higher in the halos. Molecular hydrogen 17.03 micron emission with scale heights of 550-580 pc in NGC 891 and 850 pc in NGC 5775 suggests a molecular component in a surprisingly thick layer.

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'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
preprint2011arXiv

The Infrared Properties of Super Star Clusters: Predictions from Three-Dimensional Radiative Transfer Models

With high-resolution infrared data becoming available that can probe the formation of high-mass stellar clusters for the first time, models that make testable predictions of these objects are necessary. We utilize a three-dimensional radiative transfer code, including a hierarchically clumped medium, to study the earliest stages of super star cluster evolution. We explore a range of parameter space in geometric sequences that mimic the evolution of an embedded super star cluster. The inclusion of a hierarchically clumped medium can make the envelope porous, in accordance with previous models and supporting observational evidence. The infrared luminosity inferred from observations can differ by a factor of two from the true value in the clumpiest envelopes depending on the viewing angle. The infrared spectral energy distribution also varies with viewing angle for clumpy envelopes, creating a range in possible observable infrared colors and magnitudes, silicate feature depths and dust continua. General observable features of cluster evolution differ between envelopes that are relatively opaque or transparent to mid-infrared photons. The [70]-[160] color can be used to determine star formation efficiency; the Spitzer IRAC/MIPS [8.0]-[24] color is able to constrain Rin and Rout values; and the IRAC [3.6]-[5.8] color is sensitive to the fraction of the dust distributed in clumps. Finally, in a comparison of these models to data of ultracompact HII regions, we find good agreement, suggesting that these models are physically relevant, and will provide useful diagnostic ability for datasets of resolved, embedded SSCs with the advent of high-resolution infrared telescopes like JWST.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Stability of Low Surface Brightness Disks Based on Multi-Wavelength Modeling

To investigate the structure and composition of the dusty interstellar medium (ISM) of low surface brightness (LSB) disk galaxies, we have used multiwavelength photometry to construct spectral energy distributions for three low-mass, edge-on LSB galaxies. We use Monte Carlo radiation transfer codes that include the effects of transiently heated small grains and polycyclic aromatic hydrocarbon molecules to model and interpret the data. We find that unlike the high surface brightness galaxies previously modeled, the dust disks appear to have scale heights equal to or exceeding their stellar scale heights. This result supports the findings of previous studies that low mass disk galaxies have dust scale heights comparable to their stellar scale heights and suggests that the cold ISM of low mass, LSB disk galaxies may be stable against fragmentation and gravitational collapse. This may help to explain the lack of observed dust lanes in edge-on LSB galaxies and their low current star formation rates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The standard model of low-mass star formation applied to massive stars: multi-wavelength modelling of IRAS 20126+4104

In order to investigate whether massive stars form similarly to their low-mass counterparts, we have used the standard envelope plus disc geometry successfully applied to low-mass protostars to model the near-IR to sub-millimetre SED and several mid-IR images of the embedded massive star IRAS20126+4104. We have used a Monte Carlo radiative transfer dust code to model the continuum absorption, emission and scattering through two azimuthally symmetric dust geometries, the first consisting of a rotationally flattened envelope with outflow cavities, and the second which also includes a flared accretion disc. Our results show that the envelope plus disc model reproduces the observed SED and images more accurately than the model without a disc, although the latter model more closely reproduces the morphology of the mid-IR emission within a radius of 1.1" or ~1800au. We have put forward several possible causes of this discontinuity, including inner truncation of the disc due to stellar irradiation, or precession of the outflow cavity. Our best fitting envelope plus disc model has a disc radius of 9200 au. We find that it is unlikely that the outer regions of such a disc would be in hydrostatic or centrifugal equilibrium, however we calculate that the temperatures within the disc would keep it stable to fragmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Photoionization of High Altitude Gas in a Supernova-Driven Turbulent Interstellar Medium

We investigate models for the photoionization of the widespread diffuse ionized gas in galaxies. In particular we address the long standing question of the penetration of Lyman continuum photons from sources close to the galactic midplane to large heights in the galactic halo. We find that recent hydrodynamical simulations of a supernova-driven interstellar medium have low density paths and voids that allow for ionizing photons from midplane OB stars to reach and ionize gas many kiloparsecs above the midplane. We find ionizing fluxes throughout our simulation grids are larger than predicted by one dimensional slab models, thus allowing for photoionization by O stars of low altitude neutral clouds in the Galaxy that are also detected in Halpha. In previous studies of such clouds the photoionization scenario had been rejected and the Halpha had been attributed to enhanced cosmic ray ionization or scattered light from midplane H II regions. We do find that the emission measure distributions in our simulations are wider than those derived from Halpha observations in the Milky Way. In addition, the horizontally averaged height dependence of the gas density in the hydrodynamical models is lower than inferred in the Galaxy. These discrepancies are likely due to the absence of magnetic fields in the hydrodynamic simulations and we discuss how magnetohydrodynamic effects may reconcile models and observations. Nevertheless, we anticipate that the inclusion of magnetic fields in the dynamical simulations will not alter our primary finding that midplane OB stars are capable of producing high altitude diffuse ionized gas in a realistic three-dimensional interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hot spots and a clumpy disk: Variability of brown dwarfs and stars in the young Sigma Ori cluster

The properties of accretion disks around stars and brown dwarfs in the SOri cluster (age 3 Myr) are studied based on NIR time series photometry supported by MIR spectral energy distributions. We monitor ~30 young low-mass sources over 8 nights in the J- and K-band using the duPont telescope at Las Campanas. We find three objects showing variability with J-band amplitudes >0.5 mag; five additional objects exhibit low-level variations. All three highly variable sources have been previously identified as highly variable; thus we establish the long-term nature of their flux changes. The lightcurves contain periodic components with timescales of ~0.5-8 days, but have additional irregular variations superimposed -- the characteristic behaviour for classical T Tauri stars. Based on the colour variability, we conclude that hot spots are the dominant cause of the variations in two objects, including one likely brown dwarf, with spot temperatures in the range of 6000-7000 K. For the third one (#2), a brown dwarf or very low mass star, inhomogenities at the inner edge of the disk are the likely origin of the variability. Based on mid-infrared data from Spitzer, we confirm that the three highly variable sources are surrounded by circum-(sub)-stellar disks. They show typical SEDs for T Tauri-like objects. Using SED models we infer an enhanced scaleheight in the disk for the object #2, which favours the detection of disk inhomogenities in lightcurves and is thus consistent with the information from variability. In the SOri cluster, about every fifth accreting low-mass object shows persistent high-level photometric variability. We demonstrate that estimates for fundamental parameters in such objects can be significantly improved by determining the extent and origin of the variations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Ionized Gas Towards Molecular Clumps: Physical Properties of Massive Star Forming Regions

We have conducted a search for ionized gas at 3.6 cm, using the Very Large Array, towards 31 Galactic intermediate- and high-mass clumps detected in previous millimeter continuum observations. In the 10 observed fields, 35 HII regions are identified, of which 20 are newly discovered. Many of the HII regions are multiply peaked indicating the presence of a cluster of massive stars. We find that the ionized gas tends to be associated towards the millimeter clumps; of the 31 millimeter clumps observed, 9 of these appear to be physically related to ionized gas, and a further 6 have ionized gas emission within 1'. For clumps with associated ionized gas, the combined mass of the ionizing massive stars is compared to the clump masses to provide an estimate of the instantaneous star formation efficiency. These values range from a few percent to 25%, and have an average of 7 +/- 8%. We also find a correlation between the clump mass and the mass of the ionizing massive stars within it, which is consistent with a power law. This result is comparable to the prediction of star formation by competitive accretion that a power law relationship exists between the mass of the most massive star in a cluster and the total mass of the remaining stars.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

IRS Characterization of a Debris Disk around an M-type star in NGC2547

We present 5 to 15 micron Spitzer Infrared Spectrograph (IRS) low resolution spectral data of a candidate debris disk around an M4.5 star identified as a likely member of the ~40 Myr old cluster NGC2547. The IRS spectrum shows a silicate emission feature, indicating the presence of warm, small, (sub)micron-sized dust grains in the disk. Of the fifteen previously known candidate debris disks around M-type stars, the one we discuss in this paper is the first to have an observed mid-infrared spectrum and is also the first to have measured silicate emission. We combined the IRS data with ancillary data (optical, JHKs, and Spitzer InfraRed Array Camera and 24 micron data) to build the spectral energy distribution (SED) of the source. Monte Carlo radiation transfer modeling of the SED characterized the dust disk as being very flat (h100=2AU) and extending inward within at least 0.13AU of the central star. Our analysis shows that the disk is collisionally dominated and is likely a debris disk.

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'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 'typical' 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'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
preprint1999arXiv

Ultraviolet Imaging Polarimetry of the Large Magellanic Cloud. II. Models

Motivated by new sounding-rocket wide-field polarimetric images of the Large Magellanic Cloud, we have used a three-dimensional Monte Carlo radiation transfer code to investigate the escape of near-ultraviolet photons from young stellar associations embedded within a disk of dusty material (i.e. a galaxy). As photons propagate through the disk, they may be scattered or absorbed by dust. Scattered photons are polarized and tracked until they escape to be observed; absorbed photons heat the dust, which radiates isotropically in the far-infrared, where the galaxy is optically thin. The code produces four output images: near- UV and far-IR flux, and near-UV images in the linear Stokes parameters Q and U. From these images we construct simulated UV polarization maps of the LMC. We use these maps to place constraints on the star + dust geometry of the LMC and the optical properties of its dust grains. By tuning the model input parameters to produce maps that match the observed polarization maps, we derive information about the inclination of the LMC disk to the plane of the sky, and about the scattering phase function g. We compute a grid of models with i = 28 deg., 36 deg., and 45 deg., and g = 0.64, 0.70, 0.77, 0.83, and 0.90. The model which best reproduces the observed polarization maps has i = 36 +2/-5 degrees and g ~0.7. Because of the low signal-to-noise in the data, we cannot place firm constraints on the value of g. The highly inclined models do not match the observed centro-symmetric polarization patterns around bright OB associations, or the distribution of polarization values. Our models approximately reproduce the observed ultraviolet photopolarimetry of the western side of the LMC; however, the output images depend on many input parameters and are nonunique.

0 citations0 reviewsNo code linkedOpen access