Circumstellar dust shells around longperiod variables
X. Dynamics of envelopes around standard luminous, Crich AGB stars
Technische Universität Berlin,
Zentrum für Astronomie und Astrophysik (ZAA), EW
81, Hardenbergstr. 36,
10623
Berlin,
Germany
email: dreyer@astro.physik.tuberlin.de
Received:
27
August
2010
Accepted:
15
October
2010
Context. Longperiod variables (LPVs) and Miras exhibit a pronounced variability in their luminosity with a more or less welldefined period, and they suffer large mass loss in the form of stellar winds. Owing to this extensive mass loss, they are surrounded by extended circumstellar dust shells (CDSs). The dynamics of these envelopes is the result of a complex interplay via an external excitation by the pulsating central star, dust formation, and radiative transfer.
Aims. Our study is aimed at an understanding of the dynamics of CDSs around carbonrich, standard luminous LPVs and Miras. These shells often show multiperiodicity with secondary periods as high as a few 10^{4} d superimposed on a main period that is in the range of approximately 10^{2}−10^{3} d. Such secondary periods may be caused at least in part by the presence of dust.
Methods. We consider an excitation of the CDSs either by a harmonic force, provided by the oscillation of the central star, or by a stochastic force with a continuous power spectrum. The resulting numerically computed dynamical behaviour of the shell is analysed with the help of Fourier analysis and stroboscopic maps.
Results. CDSs may be described as multioscillatory systems that are driven by the pulsating stars. A set of normal modes can be identified. The obtained periods of these modes are some 10^{3} d, which is a characteristic timescale for dust nucleation, growth, and elemental enrichment in the dust formation zone. Depending on the oscillation period and strength of the central star, the envelope reacts periodically, multi periodically, or irregularly.
Key words: stars: AGB and postAGB / circumstellar matter / stars: oscillations / chaos / hydrodynamics / methods: numerical
© ESO, 2010
1. Introduction
Longperiod variables (LPVs) and Miras are highly evolved stars located in the upper region of the asymptotic giant branch (AGB). They are characterised by low surface temperatures T < 3 × 10^{3} K and high luminosities of the order of 10^{4} L_{⊙}. They typically exhibit a pronounced variability of their intensity and suffer large mass loss in the form of a stellar wind with mass loss rates of up to 10^{4} M_{⊙} yr^{1}. Owing to their extensive mass loss, LPVs and Miras are surrounded by extended circumstellar envelopes that provide the necessary conditions for the formation of dust.
Dreyer et al. (2009, henceforth called Paper I), a previous paper of this series, studied the dynamics of circumstellar dust shells (CDSs) around high luminous LPVs and Miras. It turned out that for these kind of stars, the dynamics of the envelopes is clearly dominated by the interplay between dust formation and radiative transfer, i.e. the exterior κeffect (see Fleischer et al. 1995; Höfner et al. 1995). A regular oscillation of the shell can be solely explained by the formation of dust and the radiative pressure on it.
In this work we concentrate on the circumstellar envelopes of less luminous LPVs and Miras. For luminosities of less than 10^{4} L_{⊙}, additional input of energy and momentum is required to sustain the observed oscillations. This energy is provided by oscillations of the star’s deeper layers. Caused by the pulsational instability due to variations in the convective flux (cf. Fox & Wood 1982), hydrodynamic waves are generated in the dense stellar atmosphere and steepen into shock waves when travelling into the thinner circumstellar medium. The presence of such shock waves causes a levitation of the atmosphere producing favourable conditions for dust nucleation and growth. Once the first dust particles have formed, the momentum coupling between matter and radiation field is strongly enhanced, which again results in a massive stellar wind.
Although the shell’s oscillation has to be triggered by a stellar pulsation, the presence of dust influences the dynamics of the entire system. Even for a strictly harmonic excitation, the shell does not necessarily mimic the stellar pulsation, but instead exhibits its own normal modes (cf. Barthés & Tuchman 1994). Superimposed upon the excitation period that is in the range of approximately (10^{2}−10^{3}) d (Whitelock et al. 1991), the envelope typically shows oscillations with periods of up to a few 10^{4} d. This could explain, at least in part, the observed cycletocycle variation in the periods of some Miras (Eddington & Plakidis 1929). As pointed out by Percy & Colivas (1999), these fluctuations could come from multiperiodicity, with either a secondary period that is comparable to the primary one or a secondary period an order of magnitude longer than the primary one. In the latter case, the secondary oscillation may be explained by an oscillation of the circumstellar shell from the exterior κeffect.
In the context of dynamical systems, the CDS can be described as a multioscillatory system that is driven by the pulsating stellar atmosphere. The resulting dynamical behaviour depends on both the frequency and strength of the excitation, as well as, on the intrinsic dynamics of the shell. For a close examination of the possible normal modes of a CDS, which are controlled by the intrinsic timescales of the various physical and chemical processes, we consider response spectra of a circumstellar shell that is excited not only by a harmonic force but also by a stochastic force with a continuous power spectrum. The inherently nonlinear nature of the underlying equations introduces an element of randomness into the shell’s dynamics. This also becomes clear when the CDS is excited by a strictly harmonic force. For a sensitive range of frequencies and excitation strengths, the shell does not behave strictly periodically but shows signs of deterministic chaos.
This paper is organised as follows. Section 2 outlines our dynamical modelling method of a CDS, giving a special focus to the formulation of a purely harmonic and stochastic excitation of the CDS provided by a central star. Section 3 presents our direct methods as used to investigate the dynamical CDS behaviour. This section discusses the problems of creating and analysing discrete power spectra and deterministic maps. In Sect. 4, we apply the method to a reference model. The eigendynamics are first determined and then the complex interaction between the inner excitation by a pulsating star and the multioscillatory circumstellar shell is discussed, in this case depending on excitation frequency and strength. The paper finishes with some concluding remarks in Sect. 5.
2. Modelling of CDSs
The theoretical model of a CDS is conceived as the complete unique set of solutions of an appropriate system of fundamental dynamical equations, comprising the conservation laws, transport equations, and constitutive equations that describe both the local and global physical situation of the CDS and the various inherent processes, such as radiative transfer (Lucy 1971, 1976; Unno & Kondo 1976) or dust formation (Gail & Sedlmayr 1988; Gauger et al. 1990). These equations were listed and solved according to the procedure given in Fleischer et al. (1992). The CDSs are characterised by the fundamental parameters of their generating stars, which are mass M_{∗}, luminosity L_{∗}, effective temperature T_{∗}, and carbontooxygen abundance ratio C/O. As these parameters change only slightly with time because of stellar evolution, they are considered as given external parameters for modelling a circumstellar shell.
2.1. Monoperiodic excitation
The pulsation of an AGB star, more specifically the underlying inner κmechanism, is treated by the piston approximation adopted from the work of Wood (1979), Bowen (1988), and Bowen & Willson (1991). It is assumed that the radius of the inner boundary R_{in}(t) located a few scale heights below the stellar radius oscillates sinusoidally about an equilibrium position R_{0} with (1)This implies that, at the inner boundary the hydrodynamic velocity u_{in}, varies in a manner determined by (2)The velocity amplitude Δu, which simulates the strength of the interior pulsation of the star, as well as the pulsation period P, are input parameters.
2.2. Stochastic excitation
To examine the eigenfrequencies and the bandpass characteristics of circumstellar shells around standard and lowluminous AGB stars, we consider the excitation of the shell by external white noise as generated by, e.g., a convective atmosphere of the giants. In analogy to the piston approximation, the movement of the inner boundary R_{in}(t) is described by the following stochastic equations: (3)(4)The term Γ(t) describes the influence of the randomness or noise and is taken to be a Gaussian with a mean of zero: (5)Since we adopt white noise, Γ(t) is uncorrelated in time, which corresponds to the (auto)correlation function K(τ) being Dirac’s δfunction (6)where σ^{2} scales the intensity of the noise. According to the WienerKhinchin theorem, the spectral density H(f)^{2} of Γ(t) is given by the Fourier transform of its autocorrelation function, (7)and is therefore independent of the frequency f. The integral of the Gaussian white noise, i.e. the solution of Eq. (4), is known as the Wiener process. It is also Gaussian distributed, so the mean vanishes, and its variance increases linearly with time t. For a concrete numerical modelling of the CDS, only realisations of the stochastic process Eq. (3) where R_{in}(t) stays within predefined finite boundaries (to ensure that the inner boundary stays within the stellar radius R_{∗}) are considered. The strength of the noise σ is also an input parameter.
3. Diagnostic approach
The dynamics of (non)linear systems may be described by means of different characteristics. One of the classical ways to characterise the dynamics is to compute the power spectrum based on the Fourier transform. This provides the possibility of obtaining a deeper physical understanding of the real system dynamics by detecting various frequencies contained in a given signal. Another way is an examination by means of a deterministic approach, e.g., to create stroboscopic or specific Poincaré maps to find domains of stability in the phase space.
3.1. Fourier transform
Generally, the Fourier transform is simply a method of expressing a function in terms of the sum of its projections onto a set of basis functions. The standard basis functions used for Fourier transform are { sin(2πft), cos(2πft),f ∈ R } or, equivalently, { e^{2πift},f ∈ R } . It is the oscillatory frequency f that varies over the set of all real numbers to give an infinite collection of basis functions. We project a given time function h(t) onto our basis functions to get the Fourier amplitudes H(f) for each frequency f: (8)In general, H(f) will be complex. The norm of the amplitude, H(f) is called the Fourier spectrum, and the square H(f)^{2} is called the power spectrum of h(t).
When using a fast Fouriertransform (FFT)algorithm to compute the spectrum of a signal, we use a finite length of signal, consisting of N = 2^{γ},γ ∈ N discrete samples. However, the discretisation of a signal is not without problems. If the sampling rate f_{s} = 1/Δt is not high enough to sample the signal correctly, then a phenomenon called aliasing occurs; i.e., components of the signal at high frequencies are mistaken for components at lower frequencies.
3.2. Mapping
For dynamic systems having several degrees of freedom, it is not very practical to discuss the orbit in a multidimensional phase space. A more appropriate way is to study the intersections of orbit with a plane in phase space, which reduces the phase space dimension by one.
If the differential equations describing a system contain terms periodic in time (where P is the period), it is convenient to create stroboscopic maps showing the time trajectory for several time intervals Δt; i.e., the flow is sampled whenever t = nΔt for n ∈ N_{0}. This kind of map can be constructed for any temporal interval. Sampling by a fixed time interval Δt = P, i.e. by successive times when the orbit crosses the phase angle φ plane, i.e. t = k2π with k being an integer, results in a Poincaré map.
4. Application
The condition for a purely dustdriven wind by the exterior κmechanism is only properly fulfilled at the very end of the AGB evolution. If the stellar luminosity decreases, the dominance of the dust with respect to the CDS dynamics at first diminishes. Then, further reduction eventually leads to a vanishing outflow, if no additional energy input is provided. AGBs with low and standard luminosities generate neither selfinduced shocks and outflows, with selfmaintained oscillatory patterns caused solely by the dust formation (exterior κmechanism), nor a CDS. Generating a CDS requires input of momentum and energy via an additional mechanism in order to levitate material into the outer atmosphere, i.e. into the dust formation window. An example of such a mechanism can be the pulsation of the outer layers of a red giant. The pulsation provides the star with mechanical momentum for levitating the material in their atmosphere.
4.1. Eigendynamics of envelopes
To examine the eigendynamics of a CDS around a standard or low luminous AGB star (e.g. with the reference stellar parameters: L_{∗} = 9 × 10^{3} L_{⊙}, M_{∗} = 1 M_{⊙}, T_{∗} = 2,5 × 10^{3} K, and C/O = 1.75), the additional mechanism of pulsation is adopted to be irregular (stochastically determined), such as is generated by convection in a red giant atmosphere. The reason for this assumption is to avoid the dominance of the stellar pulsation on the CDS dynamics. The shell has been excited with different noise intensities σ = (0.1,0.5,1.0,2.0) km s^{1} and sequences from up to 24 models were calculated (cf. Table 1), each with the same set of fundamental stellar parameters, but for different realisations of the stochastic process as described in Sect. 2.2. The stochastic movement of the inner boundary R_{in}(t) and the corresponding velocity u_{in}(t) for an arbitrary stochastic realisation with intensity σ = 0.1 km s^{1} are displayed in Fig. 1.
Fig. 1 Stochastical movement of the inner boundary R_{in} about an equilibrium position R_{0} = 0.87 R_{∗} (upper panel) and associated velocity u_{in} (lower panel) for the reference model, stochastically excited with strength σ = 0.1 km s^{1}. 

Open with DEXTER 
Fig. 2 Radial structure of stochastically excited CDSs with excitation strength σ = 0.1 km s^{1} (left) and σ = 1.0 km s^{1} (right). Upper panel: gas velocity u (solid red line), gas temperature T (dashed blue line). Second panel: mass density ρ (solid orange line), radiative acceleration on dust in units of the local gravitational deceleration α (dashed grey line). Third panel: stationary nucleation rate J_{∗} (solid light green line), number of dust grains per Hatom n_{d}/n_{H} (dashed dark green line). Lower panel: degree of condensation of carbon f_{cond} (solid black line), mean particle radius ⟨a_{m}⟩ (dashed cyan line). 

Open with DEXTER 
List of model calculations and resultant quantities for the stochastically excited reference model.
Figure 2 shows snapshots of the radial envelope structure from the stochastically excited model calculation at time t = 38 800 d for the intensity strengths σ = 0.1 km s^{1} (left) and σ = 1.0 km s^{1} (right) for various physical quantities describing a CDS. The higher input of mechanical energy in the model with σ = 1.0 km s^{1} can be seen immediately from the higher amplitude of the innermost shock compared to the model with σ = 0.1 km s^{1} (uppermost panels). Since the density scale height is the same for both models in the inner region, the higher energy input yields a smoother transition to the wind region (second panels), i.e. a more levitated atmosphere, hence more favourable conditions for dust formation. The outflow is initiated by the pulsation but is accelerated by the dust. The stellar radiation is blocked by the opacity of the instantly formed shell, thereby heating the material behind it, which is located at smaller radii. This backwarming effect prevents new dust from forming and demonstrates that formation and growth of grains also strongly influences the thermodynamical shell structure. As depicted in the lowermost left panel for σ = 0.1 km s^{1}, the radial structure is rather smooth with an almost constant degree of condensation f_{cond} ≈ 0.2; i.e., only 20% of the condensable carbon has actually condensed into solid particles.
This results in an optically thin, slowly dustdriven wind with a timeaveraged mass loss rate of ⟨Ṁ⟩^{1} ≈ 9.5 × 10^{8} M_{⊙} yr^{1} and terminal velocity of ⟨u_{∞}⟩ ≈ 3.6 km s^{1}. For σ = 1.0 km s^{1}, the wind becomes rather irregular, but shows the wellknown onionlike shell structure (lowermost right panel). The calculated mass loss rate is ⟨Ṁ⟩ = 4.5 × 10^{6} M_{⊙} yr^{1}, and the final outflow velocity shows a typical value of several tens of kilometres per second for an AGB star, namely ⟨u_{∞}⟩ ≈ 25.4 km s^{1}. For both excitation strengths the nucleation zone is localised around (2−3) R_{∗} (third panels).
These stochastically excited CDSs were subjected to a Fourier transform. Since the gas velocity u(r,t) of a CDS outflow is usually taken to be a tracer of (multi)periodicity, we focused on a detailed analysis of the velocity structure. Regularly sampled with an interval Δt over a time domain of t_{max} = 160 000 d, the timedependent velocity u(r,t) was transformed at each fixed radial position r for 1 R_{∗} ≤ r ≤ 5 R_{∗}, i.e. within the dustforming zone. The spatial resolution was chosen to be Δr = 0.25 R_{∗}. To identify whether spurious frequencies are generated either by interaction between the grid with numbers of gridpoints j_{max}, grid remapping rezoinc^{2} or by the choice of the sampling interval Δt, we studied several numerical parameter combinations. These studies suggest the following optimal numerical parameter set for our reference CDS: N = 512, j_{max} = 4098, rezoinc = 0.4 P, and Δt = 83 d, which yields a Nyquistfrequency of f_{c} = 6 × 10^{3} d^{1}. Averaging all stochastic realisations provides an overall mean power spectrum of the local expansion gas velocity u, which is presented in Fig. 3 for the same excitation strengths σ = (0.1,1.0) km s^{1} as discussed in Fig. 2. See also Table 1, which lists the model parameters of excitation strength σ, the number of ensemble members No., and the resulting quantities: ensembleaveraged CDSeigenmodes f_{κ}, timeaveraged mass loss rate ⟨Ṁ⟩, outflow velocity ⟨u_{∞}⟩, and dusttogas mass ratio ⟨ρ_{d}/ρ_{g}⟩ in the dust nucleation zone at r ≈ 2.5 R_{∗}.
Fig. 3 Normalised mean power spectra of the radial gas velocity u caused by stochastically excitation of intensity σ = 0.1 km s^{1} (left) and σ = 1.0 km s^{1} (right). 

Open with DEXTER 
To gain more insight, Fig. 4 illustrates the power spectra for some selected radial positions within the circumstellar envelope. For small excitation strength σ = 0.1 km s^{1}, the response spectra are dominated by three distinct peaks at f_{κ} = 0.06 × 10^{3} d^{1} ≈ (17 000 d)^{1}, at f_{κ} = 0.51 × 10^{3} d^{1} ≈ (1 960 d)^{1}, and at f_{κ} = 1.58 × 10^{3} d^{1} ≈ (630 d)^{1}. The amplitude of the third frequency reaches its maximum near the inner edge of the dust nucleation zone at r/R_{∗} = 1.75 and diminishes rapidly with increasing distance from the star. Beyond r/R_{∗} = 1.75, dust formation sets in and the first and second maxima at f_{κ} = 0.06 × 10^{3} d^{1} ≈ (17 000 d)^{1} and at f_{κ} = 0.51 × 10^{3} d^{1} ≈ (1 960 d)^{1} become visible. They eventually dominate the spectrum for r/R_{∗} > 1.75.
Fig. 4 Selected mean power spectrum profiles of the gas velocity u at different radii r for the noise intensities σ = 0.1 km s^{1} (upper panels) and σ = 1.0 km s^{1} (lower panels). 

Open with DEXTER 
For the higher excitation strength σ = 1.0 km s^{1}, the power spectrum looks much richer in detail than for σ = 0.1 km s^{1} with at least five distinct maxima. Even closer to the star where the medium is too hot to allow for effective dust formation, two strong peaks at f_{κ} = 1.13 × 10^{3} d^{1} ≈ (885 d)^{1} and f_{κ} = 1.50 × 10^{3} d^{1} ≈ (670 d)^{1} can be seen. Here, the dynamics of the gas are influenced by the backwarming of the dust and the dilution wave originating in the dust nucleation zone (cf. Paper I). With increasing distance from the star, additional maxima appear at f_{κ} = 0.06 × 10^{3} d^{1} ≈ (17 000 d)^{1}, f_{κ} = 2.13 × 10^{3} d^{1} ≈ (470 d)^{1}, and f_{κ} = 3.2 × 10^{3}d^{1} ≈ (310 d)^{1}.
Similar to the case of a selfexcited CDS (Paper I), which can be characterised by an eigenmode P_{κ}, these frequency maxima are characteristic of the dynamics of envelope around low and standard luminous LPVs and Miras. The result shows that a CDS does not simply mimic an external excitation but has its own eigenmodes, which are determined by the complex interplay between dust formation, growth, radiative transfer, and hydrodynamics. It has also become apparent that the response of the CDS strongly depends on the excitation strength. It is assumed that the maxima at f_{κ} = 1.58 × 10^{3} d^{1} ≈ (630 d)^{1} for σ = 0.1 km s^{1} and at f_{κ} = 1.50 × 10^{3}d^{1} ≈ (670 d)^{1} for σ = 1.0 km s^{1} belong to the same eigenmode, whose frequency is detuned towards lower values as the strength of the excitation increases.
4.2. Dynamics of periodically excited envelopes
This section examines the response of the CDS to an external excitation by a sinusoidally oscillating stellar atmosphere (cf. Sect. 2.1). Owing to the inherent nonlinearities of a dustforming system, the CDS does not generally oscillate with the applied frequency but usually exhibits a much more complex behaviour. Depending on the excitation frequency, it may exhibit an irregular or periodic behaviour. Since the dynamics in CDSs are strongly influenced by the presence of dust, we focus our investigation not only on the dustformation zone, but also on the physical quantities responsible for providing favourable conditions for dust formation, such as gas velocity u, density ρ, and temperature T.
4.2.1. Excitation frequency
Figure 5 shows the most prominent shell modes normalised by the eigenfrequency f_{CDS}/f_{κ} in the dust nucleation zone at r ≈ 2.5 R_{∗} as a function of the excitation period P in the case of our reference model. Once again, the perturbation force was taken as the sine of period P according to Eqs. (1) and (2), whereas the amplitude was held at Δu = 1.0 km s^{1}.
Fig. 5 Most dominant frequencies of the CDS f_{CDS} normalised by the eigenmode f_{κ} = 0.51 × 10^{3} d^{1} (triangles), periodically excited with various stellar pulsation periods P but fixed amplitude Δu = 1.0 km s^{1}, at the dust nucleation zone r ≈ 2.5 R_{∗}. The set of eigenfrequencies are indicated by circles. The eigenperiod is labelled with P_{κ}. The excitation frequencies lie on the solid line and their harmonics on the dashdotted lines. The grey line separates the irregular (left) from the pulsationdominated domain (right). 

Open with DEXTER 
The excitation frequency f_{p} = 1/P itself and some of its harmonics are recognisable for all periods P, whereas the determined eigenmodes f_{κ} are only visible for small excitation periods P ≲ 1 250 d.

Irregular domain. For small and intermediate excitation periods P ≲ 1 250 d, the power spectra show a rather extended continuum with a number of distinct maxima. Some of these could be identified as eigenmodes of the CDS or as their harmonics. Also, the subharmonics f_{p}/2 and f_{p}/3 of the excitation frequency f_{p} = 1/P show up. As the pulsation period P increases, the excitation period and its harmonics become clearly recognisable. In contrast, the amplitude of the envelope’s eigenmode mostly decreases as P increases. Within this domain a kind of resonance phenomenon between excitation and envelope period occurs. The envelope reacts multiperiodically for a sensitive range of excitation period P and strength Δu.

Pulsationdominated domain. For longer periods P ≳ 1 250 d, the CDS becomes finally enslaved by the excitation force. The most prominent feature of the spectrum is the pulsation mode supplemented by a set of its harmonics.
List of resultant quantities from a CDS periodically disturbed.
Fig. 6 Power spectra of the local expansion velocity u (left panels) and corresponding maps of the (u,ρ,T) phase space (green) and projection onto (u,ρ)plane (yellow) (right panels) in the dust nucleation zone for different monoperiodic excitation periods P = 716 d (upper), 450 d (middle), and 376 d (lower), respectively. The stroboscopic maps were obtained by sampling [u(n0.02P), ρ(n0.02P), T(n0.02P)] (green) and the Poincaré maps by [u(nP), ρ(nP), T(nP)] (black) for 1 ≤ n ≤ ⌊ t_{max}/P ⌋ ,n ∈ N. The excitation frequencies are highlighted in the spectra. 

Open with DEXTER 
Figure 6 demonstrates the basic behaviour of the envelope for some selected monoperiodic excitation periods P = 376 d, 450 d, and 716 d chosen from the irregular shell dynamics regime. The left panels display the power spectra of the local expansion velocity u by scaling to the overall strongest mode. Finally, the right panels present the stroboscopic maps in the (u, ρ, T)phase space, supplemented by a projection onto the (u,ρ)plane. Both were sampled at an interval Δt = 0.02 P. They also show the corresponding state of the system at a constant phase angle φ, i.e. Δt = P (Poincaré map) in the dust nucleation zone at r ≈ 2.5 R_{∗}.
As can be seen in the uppermost and lowermost spectra, the CDS reacts with the same dominant mode f = 1.34 × 10^{3} d^{1} = (716 d)^{1}, even though it is excited with different periods P = 716 d and P = 376 d. This is similar to the resonance phenomenon (iP = jP_{CDS}) studied in Paper I for highluminous Miras and LPVs. However, for the case of standard luminous LPVs and Miras studied in this paper we always found j = 1, which means that the period of the CDS P_{CDS} appears as an integer multiple of the excitation period P (iP = P_{CDS}). In the context of circumstellar envelopes, this phenomenon is often referred to as multiperiodicity (cf. Winters et al. 1994; Fleischer et al. 1995; Höfner et al. 1996). As indicated by the number of the black clusters in the stroboscopic maps, the system reacts monoperiodically (single black cluster, upper right panel, i.e. i = 1) for an excitation with P = 716 d, whereas the system’s response is doubleperiodic (two black features in the lowermost plot on right, i.e. i = 2) for an excitation with P = 376 d. In other words, the system returns to the same state after a time interval Δt = iP. The corresponding curves in the (u,ρ,T)hyper plane are restricted and closed. Nevertheless, some parts of the trajectory are passed through extremely rapidly, so that these parts are only depicted by a few points on the corresponding stroboscopic map.
Fig. 7 Power spectrum and stroboscopic maps for a CDS excited with a period P = 376 d and strength Δu = 9 km s^{1}. The plot parameters for the maps are the same as in Fig. 6. 

Open with DEXTER 
Parameter study for various excitation strength Δu with fixed period P = 376 d.
For P = 450 d, as can be seen from the power spectra in the middleleft panel, the excitation period does not dominate the envelope. The amplitude of the piston frequency f_{p} = 2.22 × 10^{3} d^{1} is weak, and the frequency spectrum shows up with distinct maxima at f = 0.44 × 10^{3} d^{1} ≈ (2300 d)^{1}, which is a characteristical timescale of the CDS and f = 1.35 × 10^{3} d^{1} ≈ (740 d)^{1}, which seems to correspond to the somewhatshifted feature P = 716 d providing the resonance case for this model. This is also consistent with the stroboscopic map given in the middleright panel. In contrast to P = 716 d and 376 d, the system is not limited to a closed trajectory in the (u, ρ, T)space but rather completely fills out a continuous part of the phase space. However, the dynamics of the envelope once more seem to synchronise with the excitation period to some extent. For a constant phase angle φ of the stellar pulsation, i.e. Δt = P, the system stays inside welldefined stripes that cut through the entire phasespace area filled out by the system.
4.2.2. Excitation strength
This section studies the abovementioned reference model with regard to the response of the CDS for a sinusoidally stellar pulsation with period P = 376 d, i.e. a typical oscillation period of LPVs and Miras, over a sequence of varying strengths, i.e. amplitudes Δu. The amplitude starts with a value of Δu = 1.0 km s^{1}, so the disturbance is subsonic, and is then increased up to Δu = 9.0 km s^{1} i.e. into the supersonice regime. Table 3 lists the results of these calculations.
The excitation frequency f_{p} = 2.66 × 10^{3} d^{1} = (376 d)^{1} itself is only noticeable for a low excitation strength of Δu = 1.0 km s^{1}. Because weak in power, it vanishes with increasing excitation strengths. The inertial envelope is not able to follow the short excitation period; in fact, this mode at no time does it dominate the envelope dynamics. Actually, the first subharmonic of the excitation frequency f_{p}/2 = 1.33 × 10^{3} d^{1} ≈ (752 d)^{1} dominates the dynamics. It can be regarded as a detuned envelopeeigenmode f_{κ} = 1.50 × 10^{3} d ≈ (670 d)^{1}. Subsequently, the eigenmode at f_{κ} = 0.50 × 10^{3} d^{1} ≈ (2000 d)^{1} takes over the dominance of the shell dynamics. In general, the CDSeigenmodes are detuned towards lower frequencies as already shown in the case of a CDS around a highluminous LPV (cf. Paper I). The stellar pulsation leads to a levitation of the atmosphere; i.e., it lifts the material out of the gravitational field into the dustforming zone. With increasing strength of the pulsation amplitude Δu, the additional input of energy and momentum leads to an enhanced density in the dustforming region and supports further dust formation rather than accelerating the wind. Consequently, the stellar wind becomes far more massive, and the mass loss rate increases from ⟨Ṁ⟩ = 1.0 × 10^{5} M_{⊙} yr^{1} for Δu = 1 km s^{1} to ⟨Ṁ⟩ = 7.4 × 10^{5} M_{⊙} yr^{1} for Δu = 9 km s^{1} at nearly constant velocity. As the wind becomes more intense, the periodic depletion in dustforming material also increases, and it takes more time to reenrich the nucleation zone. The result is a longer shell period as seen in Table 3.
Figure 7 depicts the spectrum and maps of the reference CDS excited with period P = 376 d and amplitude Δu = 9 km s^{1}. With increasing amplitude the CDS behaviour changes from multiperiodic to irregular, as can be seen in the stroboscopic maps (cf. Fig. 6, lowermost panels, right and Fig. 7, right). The double periodicity breaks up. As can be seen in the power spectrum, the timescale of the stellar pulsation (excitation modes) imprints their dynamics onto the envelope close to the star. With onset of the dust formation at r ≈ 2.5 R_{∗} a continuous power spectrum appears rather than distinct modes in terms of the radial excitation frequency.
5. Summary and outlook
In this article we have investigated the dynamics of envelopes around standard and lowluminous, carbonrich AGB stars. We analysed the numerical solutions of the coupled equation system of hydrodynamics, thermodynamics, dust nucleation and growth, and radiative transfer by means of power spectra and stroboscopic maps of quantities, such as radial outflow velocity, temperature, and density.
As shown for the case of CDSs around highly luminous AGBs, the presence of dust influences the dynamical and energetic behaviour of the entire system. Firstly, it amplifies the momentum coupling between matter and radiation field, and secondly, it enforces an additional dynamical behaviour (exterior κmechanism) besides the pulsation to the dust envelope. The consequence of the coupling is the development of an envelope dynamics, which is in fact excited by the momentum and energy input by the stellar pulsation, but is operating on the timescales of the physical processes that determining the envelope. These dynamics can be characterised by a set of autonomous frequencies that are not identical to the frequency of the pulsation.
For CDSs around low and standard luminous AGBs, such distinct frequencies are not obvious. To study the shell eigendynamics, the system has to be supplied with energy and momentum. This we have done by introducing an additional force, regarded stochastic, e.g. provided by convection in the giant’s atmosphere. The goal is to supply enough energy to introduce a CDS generating mass loss, but also to minimise the influence on the CDS’s eigendynamics.
When applying this method with different stochastic noise intensities to a reference CDS, the power spectra of the radial velocity show up as a set of eigenmodes of about 10^{3} d, which represents the dynamical timescale of dust formation and element enhancement in the dust formation zone. To study the interaction between the envelopes’s eigendynamics and stellar pulsation, a parameter study of various stellar excitation periods and strengths was carried out. Depending on the stimulation, the system reacts periodically, multiperiodically, or irregularly. For short excitation periods, the considered power spectra suggest that the CDS models tend to be semiregular or even chaotic. The timescale of the excitation is much smaller than the timescale
of the various physical processes (in particular the timescale of dust formation) in the envelope, so that the shell is not able to follow the excitation. For long periods the system is dominated by the external excitation. In contrast to envelopes around high luminous AGBs, no eigenmodedominated domain could be found. In a closer examination of excitation periods from the irregular domain, we found the same resonance phenomenon (iP = jP_{CDS}) as in the case of envelopes around highluminous LPVs and Miras. However, it turns out that the integer j always equal 1 for standard luminous AGBs.
Combined with the analysis method of shells around highly luminous AGBs already presented in a previous article, we are now able to start a systematic stellar parameter study. In order to compare the results with observations, this planned study will also analyse synthetic lightcurves and spectra.
For details of grid redistribution see Fleischer (1994).
References
 Barthés, D., & Tuchman, Y. 1994, A&A, 289, 429 [NASA ADS] (In the text)
 Bowen, G. H. 1988, ApJ, 329, 299 [NASA ADS] [CrossRef] (In the text)
 Bowen, G. H., & Willson, L. A. 1991, ApJ, 375, L53 [NASA ADS] [CrossRef] (In the text)
 Eddington, A. S., & Plakidis, L. 1929, MNRAS, 90, 65 [NASA ADS] [CrossRef] (In the text)
 Dreyer, C., Hegmann, M., & Sedlmayr, E. 2009, A&A, 499, 765 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Fleischer, A. J. 1994, Ph.D. Thesis, Technische Universität Berlin, FRG (In the text)
 Fleischer, A. J., Gauger, A., & Sedlmayr, E. 1992, A&A, 266, 321 [NASA ADS] (In the text)
 Fleischer, A. J., Gauger, A., & Sedlmayr, E. 1995, A&A, 297, 543 [NASA ADS] (In the text)
 Fox, M. W., & Wood, P. R. 1982, ApJ, 259, 198 [NASA ADS] [CrossRef] (In the text)
 Gail, H. P., & Sedlmayr, E. 1988, A&A, 206, 153 [NASA ADS] (In the text)
 Gauger, A., Gail, H. P., & Sedlmayr, E. 1990, A&A, 235, 345 [NASA ADS] (In the text)
 Höfner, S., Feuchtinger, M. U., & Dorfi, E. A. 1995, A&A, 297, 815 [NASA ADS] (In the text)
 Höfner, S., Fleischer, A. J., Gauger, A., et al. 1996, A&A, 314, 204 [NASA ADS] (In the text)
 Whitelock, P. A., Feast, M. W., & Catchpole, R. M. 1991, MNRAS, 248, 276 [NASA ADS] [CrossRef] (In the text)
 Lucy, L. B. 1971, ApJ, 163, 95 [NASA ADS] [CrossRef] (In the text)
 Lucy, L. B. 1976, ApJ, 205, 482 [NASA ADS] [CrossRef] (In the text)
 Percy, J. R., & Colivas, T. 1999, PASP, 111, 94 [NASA ADS] [CrossRef] (In the text)
 Unno, W., & Kondo, M. 1976, PASJ, 28, 347 [NASA ADS] (In the text)
 Winters, J. M., Fleischer, A. J., Gauger, A., & Sedlmayr, E. 1994, A&A, 290, 623 [NASA ADS] (In the text)
 Wood, P. R. 1979, ApJ, 227, 220 [NASA ADS] [CrossRef] (In the text)
All Tables
List of model calculations and resultant quantities for the stochastically excited reference model.
All Figures
Fig. 1 Stochastical movement of the inner boundary R_{in} about an equilibrium position R_{0} = 0.87 R_{∗} (upper panel) and associated velocity u_{in} (lower panel) for the reference model, stochastically excited with strength σ = 0.1 km s^{1}. 

Open with DEXTER  
In the text 
Fig. 2 Radial structure of stochastically excited CDSs with excitation strength σ = 0.1 km s^{1} (left) and σ = 1.0 km s^{1} (right). Upper panel: gas velocity u (solid red line), gas temperature T (dashed blue line). Second panel: mass density ρ (solid orange line), radiative acceleration on dust in units of the local gravitational deceleration α (dashed grey line). Third panel: stationary nucleation rate J_{∗} (solid light green line), number of dust grains per Hatom n_{d}/n_{H} (dashed dark green line). Lower panel: degree of condensation of carbon f_{cond} (solid black line), mean particle radius ⟨a_{m}⟩ (dashed cyan line). 

Open with DEXTER  
In the text 
Fig. 3 Normalised mean power spectra of the radial gas velocity u caused by stochastically excitation of intensity σ = 0.1 km s^{1} (left) and σ = 1.0 km s^{1} (right). 

Open with DEXTER  
In the text 
Fig. 4 Selected mean power spectrum profiles of the gas velocity u at different radii r for the noise intensities σ = 0.1 km s^{1} (upper panels) and σ = 1.0 km s^{1} (lower panels). 

Open with DEXTER  
In the text 
Fig. 5 Most dominant frequencies of the CDS f_{CDS} normalised by the eigenmode f_{κ} = 0.51 × 10^{3} d^{1} (triangles), periodically excited with various stellar pulsation periods P but fixed amplitude Δu = 1.0 km s^{1}, at the dust nucleation zone r ≈ 2.5 R_{∗}. The set of eigenfrequencies are indicated by circles. The eigenperiod is labelled with P_{κ}. The excitation frequencies lie on the solid line and their harmonics on the dashdotted lines. The grey line separates the irregular (left) from the pulsationdominated domain (right). 

Open with DEXTER  
In the text 
Fig. 6 Power spectra of the local expansion velocity u (left panels) and corresponding maps of the (u,ρ,T) phase space (green) and projection onto (u,ρ)plane (yellow) (right panels) in the dust nucleation zone for different monoperiodic excitation periods P = 716 d (upper), 450 d (middle), and 376 d (lower), respectively. The stroboscopic maps were obtained by sampling [u(n0.02P), ρ(n0.02P), T(n0.02P)] (green) and the Poincaré maps by [u(nP), ρ(nP), T(nP)] (black) for 1 ≤ n ≤ ⌊ t_{max}/P ⌋ ,n ∈ N. The excitation frequencies are highlighted in the spectra. 

Open with DEXTER  
In the text 
Fig. 7 Power spectrum and stroboscopic maps for a CDS excited with a period P = 376 d and strength Δu = 9 km s^{1}. The plot parameters for the maps are the same as in Fig. 6. 

Open with DEXTER  
In the text 