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

J. E. Bjorkman

J. E. Bjorkman 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.SRastro-ph.EPastro-ph.GA

Catalog footprint

What is connected

23works
4topics
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

23 published item(s)

preprint2015arXiv

Modeling the Complete Lightcurve of omega CMa

We have used the radiative transfer code HDUST to analyze and interpret the long-term photometric behavior of the Be star omega CMa, considering four complete cycles of disk formation and dissipation. This is the first time in which a full lightcurve of a Be star was investigated and modeled including both disk build-up and dissipation phases. Based on the quite good fit of the observed data we were able to derive the history of stellar mass decretion rates (including long- and short-term changes) during the disk formation and dissipation phases in all four cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multitechnique testing of the viscous decretion disk model I. The stable and tenuous disk of the late-type Be star $β$ CMi

The viscous decretion disk (VDD) model is able to explain most of the currently observable properties of the circumstellar disks of Be stars. However, more stringent tests, focusing on reproducing multitechnique observations of individual targets via physical modeling, are needed to study the predictions of the VDD model under specific circumstances. In the case of nearby, bright Be star $β$ CMi, these circumstances are a very stable low-density disk and a late-type (B8Ve) central star. The aim is to test the VDD model thoroughly, exploiting the full diagnostic potential of individual types of observations, in particular, to constrain the poorly known structure of the outer disk if possible, and to test truncation effects caused by a possible binary companion using radio observations. We use the Monte Carlo radiative transfer code HDUST to produce model observables, which we compare with a very large set of multitechnique and multiwavelength observations that include ultraviolet and optical spectra, photometry covering the interval between optical and radio wavelengths, optical polarimetry, and optical and near-IR (spectro)interferometry. Due to the absence of large scale variability, data from different epochs can be combined into a single dataset. A parametric VDD model with radial density exponent of $n$ = 3.5, which is the canonical value for isothermal flaring disks, is found to explain observables typically formed in the inner disk, while observables originating in the more extended parts favor a shallower, $n$ = 3.0, density falloff. Modeling of radio observations allowed for the first determination of the physical extent of a Be disk (35$^{+10}_{-5}$ stellar radii), which might be caused by a binary companion. Finally, polarization data allowed for an indirect measurement of the rotation rate of the star, which was found to be $W \gtrsim 0.98$, i.e., very close to critical.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-Dimensional Radiation Transfer in Young Stellar Objects

We have updated our publicly available dust radiative transfer code (HOCHUNK3D) to include new emission processes and various 3-D geometries appropriate for forming stars. The 3-D geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2-D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium. We include the option for simple power-law envelope geometry, which combined with fractal clumping, and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from PAHs and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: 1) Outflow cavities may be more clumpy than infalling envelopes. 2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband SED slope; and related to this, 3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as YSOs. 4) Our hydrostatic equilibrium models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamical Evolution of Viscous Disks around Be Stars. I: photometry

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard $α$-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk build-up as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a 3-D non-LTE radiative transfer code. We show that lightcurves for different wavelengths are specific of a mass loss history, inclination angle and $α$ viscosity parameter. The diagnostic potential of those lightcurves is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometry and velocity structure of HD 44179's bipolar jet

In this paper we analyse a set of 33 optical spectra, which were acquired with the ARCES echelle spectrograph (R = 38,000) on the 3.5-m telescope at the Apache Point Observatory. We examine the Hα profile in each of these observations in order to determine the geometry and velocity structure of the previously discovered bipolar jet, which originates from the secondary star of HD 44179 located at the centre of the Red Rectangle nebula. Using a 3D geometric model we are able to determine the orbital coverage during which the jet occults the primary star. During the occultation, part of the Hα line profile appears in absorption. The velocity structure of the jet was determined by modelling the absorption line profile using the Sobolev approximation for each orbital phase during which we have observations. The results indicate the presence of a wide angle jet, likely responsible for observed biconical structure of the outer nebula. Furthermore, we were able to determine a likely velocity structure and rule out several others. We find that the jet is comprised of low-density, high-velocity, central region and a higher-density, lower-velocity, conical shell.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing the properties of Be star discs with spectroastrometry and NLTE radiative transfer modelling: beta CMi

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R~80,000, sub milli-arcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with three-dimensional, NLTE radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws; Keperian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed SED, polarisation and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cyclic Variability of the Circumstellar Disc of the Be Star $ζ$ Tau. II. Testing the 2D Global Disc Oscillation Model

Aims. In this paper we model, in a self-consistent way, polarimetric, photometric, spectrophotometric and interferometric observations of the classical Be star $ζ$ Tauri. Our primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. We have carried out detailed three-dimensional, NLTE radiative transfer calculations using the radiative transfer code HDUST. For the input for the code we have used the most up-to-date research on Be stars to include a physically realistic description for the central star and the circumstellar disc. We adopt a rotationally deformed, gravity darkened central star, surrounded by a disc whose unperturbed state is given by a steady-state viscous decretion disc model. We further assume that disc is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disc model for $ζ$ Tauri and a rigorous solution of the radiative transfer, we have obtained a very good fit of the time-average properties of the disc. This provides strong theoretical evidence that the viscous decretion disc model is the mechanism responsible for disc formation. With the global oscillation model we have successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations of the H$α$ and Br$γ$ lines. This result convincingly demonstrates that the oscillation pattern in the disc is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Did Fomalhaut, HR 8799, and HL Tauri Form Planets via the Gravitational Instability? Placing Limits on the Required Disk Masses

Disk fragmentation resulting from the gravitational instability has been proposed as an efficient mechanism for forming giant planets. We use the planet Fomalhaut b, the triple-planetary system HR 8799, and the potential protoplanet associated with HL Tau to test the viability of this mechanism. We choose the above systems since they harbor planets with masses and orbital characteristics favored by the fragmentation mechanism. We do not claim that these planets must have formed as the result of fragmentation, rather the reverse: if planets can form from disk fragmentation, then these systems are consistent with what we should expect to see. We use the orbital characteristics of these recently discovered planets, along with a new technique to more accurately determine the disk cooling times, to place both lower and upper limits on the disk surface density--and thus mass--required to form these objects by disk fragmentation. Our cooling times are over an order of magnitude shorter than those of Rafikov (2005),which makes disk fragmentation more feasible for these objects. We find that the required mass interior to the planet's orbital radius is ~0.1 Msun for Fomalhaut b, the protoplanet orbiting HL Tau, and the outermost planet of HR 8799. The two inner planets of HR 8799 probably could not have formed in situ by disk fragmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Polarization and SEDs from Microlensing of Circumstellar Envelopes

Microlensing surveys have proven to be tremendously fruitful in providing valuable data products for many fields of astrophysics, from eclipse lightcurves for substellar candidates to limb darkening in stellar atmospheres. We report on a program of modeling observables from microlensing of circumstellar envelopes, particularly those of red giant stars that are the most likely to show finite source effects. We will summarize work for how polarization light curves can be used to infer envelope properties and will describe recent modeling of the time dependent spectral energy distributions (SEDs) for microlensing of dusty winds. One of the most exciting developments is the possibility of measuring variable polarization from microlensing in a suitable source using the RINGO polarimeter at La Palma. Also quite interesting is the possibility of probing a dusty wind using IRAC data for a suitable source in the event that Spitzer has a ``warm'' cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

CHARA Array K'-band Measurements of the Angular Dimensions of Be Star Disks

We present the first K'-band, long-baseline interferometric observations of the northern Be stars gamma Cas, phi Per, zeta Tau, and kappa Dra. The measurements were made with multiple telescope pairs of the CHARA Array interferometer, and in every case the observations indicate that the circumstellar disks of the targets are resolved. We fit the interferometric visibilities with predictions from a simple disk model that assumes an isothermal gas in Keplerian rotation. We derive fits of the four model parameters (disk base density, radial density exponent, disk normal inclination, and position angle) for each of the targets. The resulting densities are in broad agreement with prior studies of the IR excess flux and the resulting orientations generally agree with those from interferometric H-alpha and continuum polarimetric observations. We find that the angular size of the K' disk emission is smaller than that determined for the H-alpha emission, and we argue that the difference is the result of a larger H-alpha opacity and the relatively larger neutral hydrogen fraction with increasing disk radius. All the targets are known binaries with faint companions, and we find that companions appear to influence the interferometric visibilities in the cases of phi Per and kappa Dra. We also present contemporaneous observations of the H-alpha, H-gamma, and Br-gamma emission lines. Synthetic model profiles of these lines that are based on the same disk inclination and radial density exponent as derived from the CHARA Array observations match the observed emission line strength if the disk base density is reduced by approximately 1.7 dex.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Properties of the $δ$ Scorpii Circumstellar Disk from Continuum Modeling

We present optical $WBVR$ and infrared $JHKL$ photometric observations of the Be binary system $δ$ Sco, obtained in 2000--2005, mid-infrared (10 and $18 μ$m) photometry and optical ($λλ$ 3200--10500 Å) spectropolarimetry obtained in 2001. Our optical photometry confirms the results of much more frequent visual monitoring of $δ$ Sco. In 2005, we detected a significant decrease in the object's brightness, both in optical and near-infrared brightness, which is associated with a continuous rise in the hydrogen line strenghts. We discuss possible causes for this phenomenon, which is difficult to explain in view of current models of Be star disks. The 2001 spectral energy distribution and polarization are succesfully modeled with a three-dimensional non-LTE Monte Carlo code which produces a self-consistent determination of the hydrogen level populations, electron temperature, and gas density for hot star disks. Our disk model is hydrostatically supported in the vertical direction and radially controlled by viscosity. Such a disk model has, essentially, only two free parameters, viz., the equatorial mass loss rate and the disk outer radius. We find that the primary companion is surrounded by a small (7 $R_\star$), geometrically-thin disk, which is highly non-isothermal and fully ionized. Our model requires an average equatorial mass loss rate of $1.5\times 10^{-9} M_{\sun}$ yr$^{-1}$.

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

NLTE Monte Carlo Radiative Transfer: I. The Thermal Properties of Keplerian Disks around Classical Be Stars

We present a 3-D NLTE Monte Carlo radiative transfer code that we use to study the temperature and ionization structure of Keplerian disks around Classical Be stars. The method we employ is largely similiar to the Monte Carlo transition probability method developed by Lucy. Here we present a simplification of his method that avoids the use of the macro atom concept. Our investigations of the temperature structure of Be star disks show that the disk temperature behavior is a hybrid between the behavior of Young Stellar Object (YSO) disks and Hot Star winds. The optically thick inner parts of Be star disks have temperatures that are similar to YSO disks, while the optically thin outer parts are like stellar winds. Thus, the temperature at the disk midplane initially drops, reaching a minimum at 3--5 stellar radii, after which it rises back to the optically thin radiative equilibrium temperature at large distances. On the other hand, the optically thin upper layers of the disk are approximately isothermal -- a behavior that is analogous to the hot upper layers of YSO disks. We also find that the disks are fully ionized, as expected, but there is an ionization minimum in the vicinity of the temperature minimum. Finally, we find that, despite the complex temperature structure, the infrared excess is well-approximated by an equivalent isothermal disk model whose temperature is about 60% of the stellar temperature. This is largely because, at long wavelengths, the effective photosphere of the disk is located in its isothermal regions.

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
preprint2005arXiv

The polarisation signature from microlensing of circumstellar envelopes in caustic crossing events

The high gradient of magnification across the source during both small impact parameter events and caustic crossings offers a unique opportunity for determining the surface brightness profile of the source. Furthermore, models indicate that these events can also provide an appreciable polarisation signal -- arising from differential magnification across the otherwise symmetric source. In a previous paper, polarimetric variations from point lensing of a circumstellar envelope were considered, as would be suitable for an extended envelope around a red giant. In this work we examine the polarisation in the context of caustic crossing events, the scenario which represents the most easily accessible situation for actually observing a polarisation signal in Galactic microlensing.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Observations and modeling of the inner disk region of T Tauri stars

We present observations of four T Tauri stars using long baseline infrared interferometry from the Palomar Testbed Interferometer. The target sources, T Tau N, SU Aur, RY Tau and DR Tau, are all known to be surrounded by dusty circumstellar disks. The observations directly trace the inner regions (< 1 AU) of the disk and can be used to constrain the physical properties of this material. For three of the sources observed, the infrared emission is clearly resolved. We first use geometric models to characterize the emission region size, which ranges from 0.04 to 0.3 AU in radius. We then use Monte Carlo radiation transfer models of accretion disks to jointly model the spectral energy distribution and the interferometric observations with disk models including accretion and scattering. With these models, we are able to reproduce the data set with extended emission arising from structures larger than 10 milliarcseconds contributing less than 6% of the K band emission, consistent with there being little or no envelope remaining for these Class II sources (d log(lambda*F_lambda)/d log(lambda) ~ -2 to 0 in the infrared). The radiation transfer models have inner radii for the dust similar to the geometric models; however, for RY Tau emission from gas within the inner dust radius contributes significantly to the model flux and visibility at infrared wavelengths. The main conclusion of our modeling is that emission from inner gas disks (between the magnetic truncation radius and the dust destruction radius) can be a significant component in the inner disk flux for sources with large inner dust radii.

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
preprint2003arXiv

Effects of Grain Size on the Spectral Energy Distribution of Dusty Circumstellar Envelopes

We study the effects of dust grain size on the spectral energy distribution (SED) of spherical circumstellar envelopes. Based on the self-similarity relations of dusty SEDs derived by Ivezic & Elitzur (1997), we expect an approximate invariance of the IR SED for models with different grain sizes. Approximate invariance follows from the fact that differently sized grains have similar optical properties at long wavelengths where the dust reprocesses the starlight. In this paper, we discuss what are the physical requirements on the model parameters to maintain the approximate invariance of the IR SED. Single grain size models are studied for a wide range of grain sizes in three optical depth regimes. In this study, we find limits for the cases where the IR SED is and is not capable of conveying information about grain sizes, and to what extent it does so. We find that approximate invariance occurs for a much larger range of grain sizes than previously believed, and, when approximate invariance holds, the SED is controlled mainly by one parameter, the reprocessing optical depth, a quantity that measures the fraction of starlight that is absorbed by the dust grains. Models with a grain size distribution are studied as well. For these models, we find that, in many instances, the concept of approximate invariance may be extended from the IR SED to all wavelengths. This means that, for a wide range of optical depths, models with different grain size distributions will produce very similar SEDs and, hence, the reprocessing optical depth is the only quantity that can be unambiguously obtained from the SED. The observational consequences of this result are discussed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Spatial study with the VLT of a new resolved edge-on circumstellar dust disk discovered at the periphery of the rho Ophiuchi dark cloud

We report the discovery in NIR with SofI at the NTT of a resolved circumstellar dust disk around a 2MASS source at the periphery of the rho Ophiuchi dark cloud. We present follow-up observations in J, H, and Ks-band obtained with ISAAC at the VLT, under 0.4"-seeing conditions, which unveil a dark dust lane oriented East-West between two characteristic northern and southern reflection nebulae. This new circumstellar dust disk has a radius of 2.15" (300 AU at 140 pc), and a width of 1.2" (170 AU at 140 pc). Thanks to its location at the periphery of the dense cores, it suffers small foreground visual extinction (Av=2.1\pm2.6 mag). Although this disk is seen close to edge-on, the two reflection nebulae display very different colors. We introduce a new NIR data visualization called ``Pixel NIR Color Mapping'', which allows to visualize directly the NIR colors of the nebula pixels. Thanks to this method we identify a ridge, 0.3" (40 AU at 140 pc) to the north of the dark lane and parallel to it, which displays a NIR color excess. This ridge corresponds to an unusual increase of brightness from J to Ks, which is also visible in the NTT observation obtained 130 days before the VLT one. We also find that the northern nebula shows ~3 mag more extinction than the southern nebula. We compute axisymmetric disk models to reproduce the VLT scattered light images and the spectral energy distribution from optical to NIR. Our best model, with a disk inclination i=86\pm1 deg, correctly reproduces the extension of the southern reflection nebula, but it is not able to reproduce either the observed NIR color excess in the northern nebula or the extinction difference between the two reflection nebulae. We discuss the possible origin of the peculiar asymmetrical NIR color properties of this object.

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