-
Discontinuous Galerkin finite element method for the continuum radiative transfer problem inside axis-symmetric circumstellar envelopes
Authors:
J. Perdigon,
M. Faurobert,
G. Niccolini
Abstract:
The study of the continuum radiative transfer problem inside circumstellar envelopes is both a theoretical and numerical challenge, especially in the frequency-dependent and multi-dimensional case. While approximate methods are easier to handle numerically, they often fail to accurately describe the radiation field inside complex geometries. For these cases, it is necessary to directly solve numer…
▽ More
The study of the continuum radiative transfer problem inside circumstellar envelopes is both a theoretical and numerical challenge, especially in the frequency-dependent and multi-dimensional case. While approximate methods are easier to handle numerically, they often fail to accurately describe the radiation field inside complex geometries. For these cases, it is necessary to directly solve numerically the radiative transfer equation. We investigate the accuracy of a discontinuous Galerkin finite element method (DGFEM hereafter) applied to the frequency-dependent two dimensional radiative transfer equation, coupled with the radiative equilibrium equation, inside axis-symmetric circumstellar envelopes. The DGFEM is a variant of finite element methods. It employs discontinuous elements and flux integrals along their boundaries, ensuring local conservation. However, as opposed to the classical finite-element methods, the solution is discontinuous across element edges. We implemented the method in a code and tested its accuracy by comparing our results with the benchmarks from the literature. For all the tested cases, the temperature profile agrees within one percent. Additionally, the emerging spectral energy distributions (SEDs) and images, obtained subsequently by ray-tracing techniques from the DGFEM solution, agree on average within $5~\mathrm{\%}$ and $10~\mathrm{\%}$, respectively. We show that the DGFEM can accurately describe the temperature profile inside axis-symmetric circumstellar envelopes. Consecutively the emerging SEDs and images are also well reproduced. The discontinuous Galerkin finite element method provides an alternative method (other than Monte Carlo methods for instance) for solving the radiative transfer equation, and could be used in cases that are more difficult to handle with the other methods.
△ Less
Submitted 11 October, 2024; v1 submitted 27 August, 2024;
originally announced August 2024.
-
New boundary conditions for the approximate flux-limited diffusion radiative transfer in circumstellar environments: test case study for spherically symmetric envelopes
Authors:
Jérémy Perdigon,
Gilles Niccolini,
Marianne Faurobert
Abstract:
In order to constrain the models describing circumstellar environments, it is necessary to solve the radiative transfer equation in the presence of absorption and scattering, coupled with the equation for radiative equilibrium. However, solving this problem requires much CPU time, which makes the use of automatic minimisation procedures to characterise these environments challenging. The use of ap…
▽ More
In order to constrain the models describing circumstellar environments, it is necessary to solve the radiative transfer equation in the presence of absorption and scattering, coupled with the equation for radiative equilibrium. However, solving this problem requires much CPU time, which makes the use of automatic minimisation procedures to characterise these environments challenging. The use of approximate methods is then of primary interest. One promising candidate method is the flux-limited diffusion (FLD), which recasts the radiative transfer problem into a non-linear diffusion equation. One important aspect for the accuracy of the method lies in the implementation of appropriate boundary conditions (BCs). We present new BCs for the FLD approximation in circumstellar environments that we apply here to spherically symmetric envelopes. At the inner boundary, the entering flux (coming from the star and from the envelope itself) may be written in the FLD formalism and provides us with an adequate BC. At the free outer boundary, we used the FLD formalism to constrain the ratio of the mean radiation intensity over the emerging flux. In both cases we derived non-linear mixed BCs relating the surface values of the mean specific intensity and its gradient. We implemented these conditions and compared the results with previous benchmarks and the results of a Monte Carlo radiative transfer code. A comparison with results derived from BCs that were previously proposed in other contexts is presented as well. For all the tested cases, the average relative difference with the benchmark results is below 2\% for the temperature profile and below 6\% for the corresponding spectral energy distribution or the emerging flux. We point out that the FLD method together with the new outer BC also allows us to derive an approximation for the emerging flux.
△ Less
Submitted 29 July, 2021;
originally announced July 2021.
-
Mid-infrared circumstellar emission of the long-period Cepheid l Carinae resolved with VLTI/MATISSE
Authors:
V. Hocdé,
N. Nardetto,
A. Matter,
E. Lagadec,
A. Mérand,
P. Cruzalèbes,
A. Meilland,
F. Millour,
B. Lopez,
P. Berio,
G. Weigelt,
R. Petrov,
J. W. Isbell,
W. Jaffe,
P. Kervella,
A. Glindemann,
M. Schöller,
F. Allouche,
A. Gallenne,
A. Domiciano de Souza,
G. Niccolini,
E. Kokoulina,
J. Varga,
S. Lagarde,
J. -C. Augereau
, et al. (129 additional authors not shown)
Abstract:
The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. W…
▽ More
The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars.
△ Less
Submitted 31 March, 2021;
originally announced March 2021.
-
A thin shell of ionized gas explaining the IR excess of classical Cepheids
Authors:
V. Hocdé,
N. Nardetto,
E. Lagadec,
G. Niccolini,
A. Domiciano de Souza,
A. Mérand,
P. Kervella,
A. Gallenne,
M. Marengo,
B. Trahin,
W. Gieren,
G. Pietrzynski,
S. Borgniet,
L. Breuval,
B. Javanmardi
Abstract:
Despite observational evidences, InfraRed (IR) excess of classical Cepheids are seldom studied and poorly understood, but probably induces systematics on the Period-Luminosity (PL) relation used in the calibration of the extragalactic distance scale. This study aims to understand the physical origin of the IR excess found in the spectral energy distribution (SED) of 5 Cepheids : RS Pup (P=41.46d),…
▽ More
Despite observational evidences, InfraRed (IR) excess of classical Cepheids are seldom studied and poorly understood, but probably induces systematics on the Period-Luminosity (PL) relation used in the calibration of the extragalactic distance scale. This study aims to understand the physical origin of the IR excess found in the spectral energy distribution (SED) of 5 Cepheids : RS Pup (P=41.46d), zeta Gem (P=10.15d), eta Aql (P=7.18d), V Cen (P=5.49d) and SU Cyg (P=3.85d). A time series of atmospheric models along the pulsation cycle are fitted to a compilation of data, including optical and near-IR photometry, Spitzer spectra (secured at a specific phase), interferometric angular diameters, effective temperature and radial velocity measurements. Herschel images in two bands are also analyzed qualitatively. In this fitting process, based on the SPIPS algorithm, a residual is found in the SED, whatever the pulsation phase, and for wavelengths larger than about $1.2μ$m, which corresponds to the so-determined infrared excess of Cepheids. This IR excess is then corrected from interstellar medium absorption in order to infer or not the presence of dust shells, and is finally used in order to fit a model of a shell of ionized gas. For all Cepheids, we find a continuum IR excess increasing up to about -0.1 magnitudes at 30$μ$m, which cannot be explained by a hot or cold dust model of CircumStellar Environment (CSE). We show, for the first time, that the IR excess of Cepheids can be explained by free-free emission from a thin shell of ionized gas, with a thickness of about 15% of the star radius, a mass of $10^{-9}-10^{-7}$ Msol and a temperature ranging from 3500 to 4500K. This result has to be tested with interferometers operating in visible, in the mid-IR or in the radio domain. The impact of such CSEs of ionized gas on the PL relation of Cepheids needs also more investigations.
△ Less
Submitted 26 September, 2019;
originally announced September 2019.
-
ALMA imaging of the nascent planetary nebula IRAS 15103-5754
Authors:
Jose F. Gomez,
Gilles Niccolini,
Olga Suarez,
Luis F. Miranda,
J. Ricardo Rizzo,
Lucero Uscanga,
James A. Green,
Itziar de Gregorio-Monsalvo
Abstract:
We present continuum and molecular line (CO, C$^{18}$O, HCO$^+$) observations carried out with the Atacama Large Millimeter/submillimeter Array toward the "water fountain" star IRAS 15103-5754, an object that could be the youngest PN known. We detect two continuum sources, separated by $0.39\pm 0.03$ arcsec. The emission from the brighter source seems to arise mainly from ionized gas, thus confirm…
▽ More
We present continuum and molecular line (CO, C$^{18}$O, HCO$^+$) observations carried out with the Atacama Large Millimeter/submillimeter Array toward the "water fountain" star IRAS 15103-5754, an object that could be the youngest PN known. We detect two continuum sources, separated by $0.39\pm 0.03$ arcsec. The emission from the brighter source seems to arise mainly from ionized gas, thus confirming the PN nature of the object. The molecular line emission is dominated by a circumstellar torus with a diameter of $\simeq 0.6$ arcsec (2000 au) and expanding at $\simeq 23$ km s$^{-1}$. We see at least two gas outflows. The highest-velocity outflow (deprojected velocities up to 250 km s$^{-1}$), traced by the CO lines, shows a biconical morphology, whose axis is misaligned $\simeq 14^\circ$ with respect to the symmetry axis of the torus, and with a different central velocity (by $\simeq 8$ km s$^{-1}$). An additional high-density outflow (traced by HCO$^+$) is oriented nearly perpendicular to the torus. We speculate that IRAS 15103-5754 was a triple stellar system that went through a common envelope phase, and one of the components was ejected in this process. A subsequent low-collimation wind from the remaining binary stripped out gas from the torus, creating the conical outflow. The high velocity of the outflow suggests that the momentum transfer from the wind was extremely efficient, or that we are witnessing a very energetic mass-loss event.
△ Less
Submitted 22 August, 2018; v1 submitted 10 August, 2018;
originally announced August 2018.
-
Grown-up stars physics with MATISSE
Authors:
Florentin Millour,
J Hron,
A Chiavassa,
G Weigelt,
A Soulain,
Z Khorrami,
A Meilland,
N Nardetto,
C Paladini,
A Domiciano de Souza,
G Niccolini,
K. -H Hofmann,
D Schertl,
P Stee,
P Bendjoya,
F Thévenin,
F Vakili,
P Berio,
T Lanz,
A Matter,
P Cruzalèbes,
R Petrov,
B Lopez
Abstract:
MATISSE represents a great opportunity to image the environment around massive and evolved stars. This will allow one to put constraints on the circumstellar structure, on the mass ejection of dust and its reorganization , and on the dust-nature and formation processes. MATISSE measurements will often be pivotal for the understanding of large multiwavelength datasets on the same targets collected…
▽ More
MATISSE represents a great opportunity to image the environment around massive and evolved stars. This will allow one to put constraints on the circumstellar structure, on the mass ejection of dust and its reorganization , and on the dust-nature and formation processes. MATISSE measurements will often be pivotal for the understanding of large multiwavelength datasets on the same targets collected through many high-angular resolution facilities at ESO like sub-millimeter interferometry (ALMA), near-infrared adaptive optics (NACO, SPHERE), interferometry (PIONIER, GRAVITY), spectroscopy (CRIRES), and mid-infrared imaging (VISIR). Among main sequence and evolved stars, several cases of interest have been identified that we describe in this paper.
△ Less
Submitted 8 July, 2016;
originally announced July 2016.
-
Stars with the B[e] phenomenon seen by long baseline interferometry
Authors:
Marcelo Borges Fernandes,
Olivier Chesneau,
Michaela Kraus,
Lydia Cidale,
Anthony Meilland,
Philippe Bendjoya,
Armando Domiciano de Souza,
Gilles Niccolini,
Ileana Andruchow,
Samer Kanaan,
Philippe Stee,
Florentin Millour,
Alain Spang,
Michel Cure
Abstract:
Thanks to the high spatial resolution provided by long baseline interferometry, it is possible to understand the complex circumstellar geometry around stars with the B[e] phenomenon. These stars are composed by objects in different evolutionary stages, like high- and low-mass evolved stars, intermediate-mass pre-main sequence stars and symbiotic objects. However, up to now more than 50% of the con…
▽ More
Thanks to the high spatial resolution provided by long baseline interferometry, it is possible to understand the complex circumstellar geometry around stars with the B[e] phenomenon. These stars are composed by objects in different evolutionary stages, like high- and low-mass evolved stars, intermediate-mass pre-main sequence stars and symbiotic objects. However, up to now more than 50% of the confirmed B[e] stars are not well classified, being called unclassified B[e] stars. From instruments like VLTI/AMBER and VLTI/MIDI, we have identified the presence of gaseous and dusty circumstellar disks, which have provided us with some hints related to the nature of these objects. Here, we show our results for two galactic stars with the B[e] phenomenon, HD 50138 and CPD-529243, based on interferometric measurements.
△ Less
Submitted 15 October, 2012; v1 submitted 9 October, 2012;
originally announced October 2012.
-
The galactic unclassified B[e] star HD 50138. II. Interferometric constraints on the close circumstellar environment
Authors:
Marcelo Borges Fernandes,
Anthony Meilland,
Philippe Bendjoya,
Armando Domiciano de Souza,
Gilles Niccolini,
Olivier Chesneau,
Florentin Millour,
Alain Spang,
Philippe Stee,
Michaela Kraus
Abstract:
HD 50138 is a southern star that presents the B[e] phenomenon, but its evolutionary stage is still not well known. This object presents spectral variability, which can be explained by outbursts or shell phases and spectropolarimetric observations have shown the presence of a non-spherically symmetric circumstellar environment that is responsible for the B[e] phenomenon. Based on recent optical lon…
▽ More
HD 50138 is a southern star that presents the B[e] phenomenon, but its evolutionary stage is still not well known. This object presents spectral variability, which can be explained by outbursts or shell phases and spectropolarimetric observations have shown the presence of a non-spherically symmetric circumstellar environment that is responsible for the B[e] phenomenon. Based on recent optical long baseline interferometric observations from the VLTI/MIDI and VLTI/AMBER, and also from the Keck segment-tilting experiment, we study the structure of the circumstellar environment of HD 50138, through a geometrical analytical modeling, also using the recent LITpro software and considering a large space of parameters. We resolve and describe its circumstellar geometry for the first time in detail. The presence of a dusty circumstellar disk with an orientation onto the sky-plane of 71+-7 degrees, which is perpendicular to the polarimetric measurements from the literature, was derived. We also derived that HD 50138 is seen under an intermediate angle related to the line of sight, 56+-4 degrees. In addition, the structure of the disk and the flux contributions of the gas and dust components is discussed.
△ Less
Submitted 25 January, 2011;
originally announced January 2011.
-
Fast ray-tracing algorithm for circumstellar structures (FRACS) I. Algorithm description and parameter-space study for mid-IR interferometry of B[e] stars
Authors:
Gilles Niccolini,
Philippe Bendjoya,
Armando Domiciano De Souza
Abstract:
The physical interpretation of spectro-interferometric data is strongly model-dependent. On one hand, models involving elaborate radiative transfer solvers are too time consuming in general to perform an automatic fitting procedure and derive astrophysical quantities and their related errors. On the other hand, using simple geometrical models does not give sufficient insights into the physics of t…
▽ More
The physical interpretation of spectro-interferometric data is strongly model-dependent. On one hand, models involving elaborate radiative transfer solvers are too time consuming in general to perform an automatic fitting procedure and derive astrophysical quantities and their related errors. On the other hand, using simple geometrical models does not give sufficient insights into the physics of the object. We propose to stand in between these two extreme approaches by using a physical but still simple parameterised model for the object under consideration. Based on this philosophy, we developed a numerical tool optimised for mid-infrared (mid-IR) interferometry, the fast ray-tracing algorithm for circumstellar structures (FRACS) which can be used as a stand-alone model, or as an aid for a more advanced physical description or even for elaborating observation strategies. FRACS is based on the ray-tracing technique without scattering, but supplemented with the use of quadtree meshes and the full symmetries of the axisymmetrical problem to significantly decrease the necessary computing time to obtain e.g. monochromatic images and visibilities. We applied FRACS in a theoretical study of the dusty circumstellar environments (CSEs) of B[e] supergiants (sgB[e]) in order to determine which information (physical parameters) can be retrieved from present mid-IR interferometry (flux and visibility). From a set of selected dusty CSE models typical of sgB[e] stars we show that together with the geometrical parameters (position angle, inclination, inner radius), the temperature structure (inner dust temperature and gradient) can be well constrained by the mid-IR data alone. Our results also indicate that the determination of the parameters characterising the CSE density structure is more challenging but, in some cases, upper limits as well as correlations on the parameters characterising the mass loss can be obtained. Good constraints for the sgB[e] central continuum emission (central star and inner gas emissions) can be obtained whenever its contribution to the total mid-IR flux is only as high as a few percents. Ray-tracing parameterised models such as FRACS are thus well adapted to prepare and/or interpret long wavelengths (from mid-IR to radio) observations at present (e.g. VLTI/MIDI) and near-future (e.g. VLTI/MATISSE, ALMA) interferometers.
△ Less
Submitted 28 September, 2010;
originally announced September 2010.
-
Fast ray-tracing algorithm for circumstellar structures (FRACS). II. Disc parameters of the B[e] supergiant CPD-57° 2874 from VLTI/MIDI data
Authors:
Armando Domiciano De Souza,
Philippe Bendjoya,
Gilles Niccolini,
Olivier Chesneau,
Marcelo Borges Fernandes,
A. C. Carciofi,
A. Spang,
Philippe Stee,
Thomas Driebe
Abstract:
B[e] supergiants are luminous, massive post-main sequence stars exhibiting non-spherical winds, forbidden lines, and hot dust in a disc-like structure. The physical properties of their rich and complex circumstellar environment (CSE) are not well understood, partly because these CSE cannot be easily resolved at the large distances found for B[e] supergiants (typically $\ga 1$~kpc). From mid-IR spe…
▽ More
B[e] supergiants are luminous, massive post-main sequence stars exhibiting non-spherical winds, forbidden lines, and hot dust in a disc-like structure. The physical properties of their rich and complex circumstellar environment (CSE) are not well understood, partly because these CSE cannot be easily resolved at the large distances found for B[e] supergiants (typically $\ga 1$~kpc). From mid-IR spectro-interferometric observations obtained with VLTI/MIDI we seek to resolve and study the CSE of the Galactic B[e] supergiant CPD-57\degr\,2874. For a physical interpretation of the observables (visibilities and spectrum) we use our ray-tracing radiative transfer code (FRACS), which is optimised for thermal spectro-interferometric observations. Thanks to the short computing time required by FRACS ($<10$~s per monochromatic model), best-fit parameters and uncertainties for several physical quantities of CPD-57\degr\,2874 were obtained, such as inner dust radius, relative flux contribution of the central source and of the dusty CSE, dust temperature profile, and disc inclination. The analysis of VLTI/MIDI data with FRACS allowed one of the first direct determinations of physical parameters of the dusty CSE of a B[e] supergiant based on interferometric data and using a full model-fitting approach. In a larger context, the study of B[e] supergiants is important for a deeper understanding of the complex structure and evolution of hot, massive stars.
△ Less
Submitted 23 September, 2010; v1 submitted 17 September, 2010;
originally announced September 2010.
-
Searching for sub-stellar companion into the LkCa15 proto-planetary disk
Authors:
M. Bonavita,
G. Chauvin,
A. Boccaletti,
V. Pietu,
P. Baudoz,
J. L. Beuzit,
A. Dutrey,
S. Guilloteau,
A. M. Lagrange,
D. Mouillet,
G. Niccolini
Abstract:
Recent sub-millimetric observations at the Plateau de Bure interferometer evidenced a cavity at ~ 46 AU in radius into the proto-planetary disk around the T Tauri star LkCa15 (V1079 Tau), located in the Taurus molecular cloud. Additional Spitzer observations have corroborated this result possibly explained by the presence of a massive (>= 5 MJup) planetary mass, a brown dwarf or a low mass star co…
▽ More
Recent sub-millimetric observations at the Plateau de Bure interferometer evidenced a cavity at ~ 46 AU in radius into the proto-planetary disk around the T Tauri star LkCa15 (V1079 Tau), located in the Taurus molecular cloud. Additional Spitzer observations have corroborated this result possibly explained by the presence of a massive (>= 5 MJup) planetary mass, a brown dwarf or a low mass star companion at about 30 AU from the star. We used the most recent developments of high angular resolution and high contrast imaging to search directly for the existence of this putative companion, and to bring new constraints on its physical and orbital properties. The NACO adaptive optics instrument at VLT was used to observe LkCa15 using a four quadrant phase mask coronagraph to access small angular separations at relatively high contrast. A reference star at the same parallactic angle was carefully observed to optimize the quasi-static speckles subtraction (limiting our sensitivity at less than 1.0). Although we do not report any positive detection of a faint companion that would be responsible for the observed gap in LkCa15's disk (25-30 AU), our detection limits start constraining its probable mass, semi-major axis and eccentricity. Using evolutionary model predictions, Monte Carlo simulations exclude the presence of low eccentric companions with masses M >= 6 M Jup and orbiting at a >= 100 AU with significant level of confidence. For closer orbits, brown dwarf companions can be rejected with a detection probability of 90% down to 80 AU (at 80% down to 60 AU). Our detection limits do not access the star environment close enough to fully exclude the presence of a brown dwarf or a massive planet within the disk inner activity (i.e at less than 30 AU). Only, further and higher contrast observations should unveil the existence of this putative companion inside the LkCa15 disk.
△ Less
Submitted 20 April, 2010;
originally announced April 2010.
-
Dust cloud formation in stellar environments. II. Two-dimensional models for structure formation around AGB stars
Authors:
P. Woitke,
G. Niccolini
Abstract:
This paper reports on computational evidence for the formation of cloud-like dust structures around C-rich AGB stars. This spatio-temporal structure formation process is caused by a radiative/thermal instability of dust forming gases as identified by Woitke et al.(2000). Our 2D (axisymmetric) models combine a time-dependent description of the dust formation process according to Gail & Sedlmayr (…
▽ More
This paper reports on computational evidence for the formation of cloud-like dust structures around C-rich AGB stars. This spatio-temporal structure formation process is caused by a radiative/thermal instability of dust forming gases as identified by Woitke et al.(2000). Our 2D (axisymmetric) models combine a time-dependent description of the dust formation process according to Gail & Sedlmayr (1988) with detailed, frequency-dependent continuum radiative transfer by means of a Monte Carlo method (Niccolini et al.2003) in an otherwise static medium (v=0). These models show that the formation of dust behind already condensed regions, which shield the stellar radiation field, is strongly favoured. In the shadow of these clouds, the temperature decreases by several hundred Kelvin which triggers the subsequent formation of dust and ensures its thermal stability. Considering an initially dust-free gas with small density inhomogeneities, we find that finger-like dust structures develop which are cooler than the surroundings and point towards the centre of the radiant emission, similar to the cometary knots observed in planetary nebulae and star formation regions. Compared to a spherical symmetric reference model, the clumpy dust distribution has little effect on the spectral energy distribution, but dominates the optical appearance in near IR monochromatic images.
△ Less
Submitted 28 April, 2004;
originally announced April 2004.
-
The 2D Continuum Radiative Transfer Problem: Benchmark Results for Disk Configurations
Authors:
I. Pascucci,
S. Wolf,
J. Steinacker,
C. P. Dullemond,
Th. Henning,
G. Niccolini,
P. Woitke,
B. Lopez
Abstract:
We present benchmark problems and solutions for the continuum radiative transfer (RT) in a 2D disk configuration. The reliability of three Monte-Carlo and two grid-based codes is tested by comparing their results for a set of well-defined cases which differ for optical depth and viewing angle. For all the configurations, the overall shape of the resulting temperature and spectral energy distribu…
▽ More
We present benchmark problems and solutions for the continuum radiative transfer (RT) in a 2D disk configuration. The reliability of three Monte-Carlo and two grid-based codes is tested by comparing their results for a set of well-defined cases which differ for optical depth and viewing angle. For all the configurations, the overall shape of the resulting temperature and spectral energy distribution is well reproduced. The solutions we provide can be used for the verification of other RT codes.We also point out the advantages and disadvantages of the various numerical techniques applied to solve the RT problem.
△ Less
Submitted 15 February, 2004;
originally announced February 2004.
-
Are dust shell models well-suited to explain interferometric data of late-type stars in the near-infrared?
Authors:
P. Schuller,
P. Salomé,
G. Perrin,
B. Mennesson,
G. Niccolini,
P. de Laverny,
S. Ridgway,
V. Coudé du Foresto,
W. A. Traub,
.
Abstract:
Recently available near-infrared interferometric data on late-type stars show a strong increase of diameter for asymptotic giant branch (AGB) stars between the K (2.0 - 2.4 μm) and L (3.4 - 4.1 μm) bands. Aiming at an explanation of these findings, we chose the objects αOrionis (Betelgeuse), SW Virginis, and R Leonis, which are of different spectral types and stages of evolution, and which are s…
▽ More
Recently available near-infrared interferometric data on late-type stars show a strong increase of diameter for asymptotic giant branch (AGB) stars between the K (2.0 - 2.4 μm) and L (3.4 - 4.1 μm) bands. Aiming at an explanation of these findings, we chose the objects αOrionis (Betelgeuse), SW Virginis, and R Leonis, which are of different spectral types and stages of evolution, and which are surrounded by circumstellar envelopes with different optical thicknesses. For these stars, we compared observations with spherically symmetric dust shell models. Photometric and 11 μm interferometric data were also taken into account to further constrain the models. -- [...] -- We conclude that AGB models comprising a photosphere and dust shell, although consistent with SED data and also interferometric data in K and at 11 μm, cannot explain the visibility data in L; an additional source of model opacity, possibly related to a gas component, is needed in L to be consistent with the visibility data.
△ Less
Submitted 22 January, 2004;
originally announced January 2004.