% aa.dem % AA vers. 4.01, LaTeX class for Astronomy & Astrophysics % demonstration file % (c) Springer-Verlag HD %----------------------------------------------------------------------- % %\documentclass[referee]{aa} % for a referee version % \documentclass{aa} % \begin{document} \thesaurus{06 % A&A Section 6: Form. struct. and evolut. of stars (03.11.1; % Cosmogony, 16.06.1; % Planets and satellites: general, 19.06.1; % Solar system: general, 19.37.1; % Stars: formation of, 19.53.1; % Stars: oscillations of, 19.63.1)} % Stars: structure of. % \title{Hydrodynamics of giant planet formation} \subtitle{I. Overviewing the $\kappa$-mechanism} \author{G. Wuchterl % \inst{1} % \and % C. Ptolemy\inst{2}\fnmsep\thanks{Just to show the usage % of the elements in the author field} } \offprints{G. Wuchterl} \institute{Institute for Astronomy (IfA), University of Vienna, T\"urkenschanzstrasse 17, A-1180 Vienna\\ email: wuchterl@amok.ast.univie.ac.at % \and % University of Alexandria, Department of Geography\\ % email: c.ptolemy@hipparch.uheaven.space % \thanks{The university of heaven temporarily does not % accept e-mails} } \date{Received September 15, 1996; accepted March 16, 1997} \maketitle \begin{abstract} To investigate the physical nature of the `nuc\-leated instability' of proto giant planets (Mizuno \cite{mizuno}), the stability of layers in static, radiative gas spheres is analysed on the basis of Baker's \cite{baker} standard one-zone model. It is shown that stability depends only upon the equations of state, the opacities and the local thermodynamic state in the layer. Stability and instability can therefore be expressed in the form of stability equations of state which are universal for a given composition. The stability equations of state are calculated for solar composition and are displayed in the domain $-14 \leq \lg \rho / \mathrm{[g\, cm^{-3}]} \leq 0 $, $ 8.8 \leq \lg e / \mathrm{[erg\, g^{-1}]} \leq 17.7$. These displays may be used to determine the one-zone stability of layers in stellar or planetary structure models by directly reading off the value of the stability equations for the thermodynamic state of these layers, specified by state quantities as density $\rho$, temperature $T$ or specific internal energy $e$. Regions of instability in the $(\rho,e)$-plane are described and related to the underlying microphysical processes. Vibrational instability is found to be a common phenomenon at temperatures lower than the second He ionisation zone. The $\kappa$-mechanism is widespread under `cool' conditions. \keywords{giant planet formation -- $\kappa$-mechanism -- stability of gas spheres } \end{abstract} % %________________________________________________________________ \section{Introduction} In the \emph{nucleated instability\/} (also called core instability) hypothesis of giant planet formation, a critical mass for static core envelope protoplanets has been found. Mizuno (\cite{mizuno}) determined the critical mass of the core to be about $12 \,M_\oplus$ ($M_\oplus=5.975 \, 10^{27}\,\mathrm{g}$ is the Earth mass), which is independent of the outer boundary conditions and therefore independent of the location in the solar nebula. This critical value for the core mass corresponds closely to the cores of today's giant planets. Although no hydrodynamical study has been available many workers conjectured that a collapse or rapid contraction will ensue after accumulating the critical mass. The main motivation for this article is to investigate the stability of the static envelope at the critical mass. With this aim the local, linear stability of static radiative gas spheres is investigated on the basis of Baker's (\cite{baker}) standard one-zone model. The nonlinear, hydrodynamic evolution of the protogiant planet beyond the critical mass, as calculated by Wuchterl (\cite{wuchterl}), will be described in a forthcoming article. The fact that Wuchterl (\cite{wuchterl}) found the excitation of hydrodynamical waves in his models raises considerable interest on the transition from static to dynamic evolutionary phases of the protogiant planet at the critical mass. The waves play a crucial role in the development of the so-called nucleated instability in the nucleated instability hypothesis. They lead to the formation of shock waves and massive outflow phenomena. The protoplanet evolves into a new quasi-equilibrium structure with a \emph{pulsating\/} envelope, after the mass loss phase has declined. Phenomena similar to the ones described above for giant planet formation have been found in hydrodynamical models concerning star formation where protostellar cores explode (Tscharnuter \cite{tscharnuter}, Balluch \cite{balluch}), whereas earlier studies found quasi-steady collapse flows. The similarities in the (micro)physics, i.e., constitutive relations of protostellar cores and protogiant planets serve as a further motivation for this study. %__________________________________________________________________ \section{Baker's standard one-zone model} % Two column figure (place early!) %______________________________________________ Gamma_1 (lg rho, lg e) \begin{figure*} % \vspace{4cm} \rule{0.4pt}{4cm}% line thickness, height of picture % \caption{Adiabatic exponent $\Gamma_1$} \hfill \parbox[b]{5cm}{\caption[]{Adiabatic exponent $\Gamma_1$. $\Gamma_1$ is plotted as a function of $\lg$ internal energy $\mathrm{[erg\,g^{-1}]}$ and $\lg$ density $\mathrm{[g\,cm^{-3}]}$ }}% \label{FigGam}% \end{figure*} % In this section the one-zone model of Baker (\cite{baker}), originally used to study the Cephe{\"{\i}}d pulsation mechanism, will be briefly reviewed. The resulting stability criteria will be rewritten in terms of local state variables, local timescales and constitutive relations. Baker (\cite{baker}) investigates the stability of thin layers in self-gravitating, spherical gas clouds with the following properties: \begin{itemize} \item hydrostatic equilibrium, \item thermal equilibrium, \item energy transport by grey radiation diffusion. \end{itemize} For the one-zone-model Baker obtains necessary conditions for dynamical, secular and vibrational (or pulsational) stability [Eqs.\ (34a,\,b,\,c) in Baker \cite{baker}]. Using Baker's notation: \[ \begin{array}{lp{0.8\linewidth}} M_{\mathrm{r}} & mass internal to the radius $r$ \\ m & mass of the zone \\ r_0 & unperturbed zone radius \\ \rho_0 & unperturbed density in the zone \\ T_0 & unperturbed temperature in the zone \\ L_{r0} & unperturbed luminosity \\ E_{\mathrm{th}} & thermal energy of the zone \end{array} \] \noindent and with the definitions of the \emph{local cooling time\/} (see Fig.~\ref{FigGam}) \begin{equation} \tau_{\mathrm{co}} = \frac{E_{\mathrm{th}}}{L_{r0}} \,, \end{equation} and the \emph{local free-fall time} \begin{equation} \tau_{\mathrm{ff}} = \sqrt{ \frac{3 \pi}{32 G} \frac{4\pi r_0^3}{3 M_{\mathrm{r}}} }\,, \end{equation} Baker's $K$ and $\sigma_0$ have the following form: \begin{eqnarray} \sigma_0 & = & \frac{\pi}{\sqrt{8}} \frac{1}{ \tau_{\mathrm{ff}}} \\ K & = & \frac{\sqrt{32}}{\pi} \frac{1}{\delta} \frac{ \tau_{\mathrm{ff}} } { \tau_{\mathrm{co}} }\,; \end{eqnarray} where $ E_{\mathrm{th}} \approx m (P_0/{\rho_0})$ has been used and \begin{equation} \begin{array}{l} \delta = - \left( \frac{ \partial \ln \rho }{ \partial \ln T } \right)_P \\[1cm] e=mc^2 \end{array} \end{equation} is a thermodynamical quantity which is of order $1$ and equal to $1$ for nonreacting mixtures of classical perfect gases. The physical meaning of $ \sigma_0 $ and $K$ is clearly visible in the equations above. $\sigma_0$ represents a frequency of the order one per free-fall time. $K$ is proportional to the ratio of the free-fall time and the cooling time. Substituting into Baker's criteria, using thermodynamic identities and definitions of thermodynamic quantities, \begin{displaymath} \Gamma_1 = \left( \frac{ \partial \ln P}{ \partial\ln \rho} \right)_{S} \, , \; \chi^{}_\rho = \left( \frac{ \partial \ln P}{ \partial\ln \rho} \right)_{T} \, , \; \kappa^{}_{P} = \left( \frac{ \partial \ln \kappa}{ \partial\ln P} \right)_{T} \end{displaymath} \begin{displaymath} \nabla_{\mathrm{ad}} = \left( \frac{ \partial \ln T} { \partial\ln P} \right)_{S} \, , \; \chi^{}_T = \left( \frac{ \partial \ln P} { \partial\ln T} \right)_{\rho} \, , \; \kappa^{}_{T} = \left( \frac{ \partial \ln \kappa} { \partial\ln T} \right)_{T} \end{displaymath} one obtains, after some pages of algebra, the conditions for \emph{stability\/} given below: \begin{eqnarray} \frac{\pi^2}{8} \frac{1}{\tau_{\mathrm{ff}}^2} ( 3 \Gamma_1 - 4 ) & > & 0 \label{ZSDynSta} \\ \frac{\pi^2}{\tau_{\mathrm{co}} \tau_{\mathrm{ff}}^2} \Gamma_1 \nabla_{\mathrm{ad}} \left[ \frac{ 1- 3/4 \chi^{}_\rho }{ \chi^{}_T } ( \kappa^{}_T - 4 ) + \kappa^{}_P + 1 \right] & > & 0 \label{ZSSecSta} \\ \frac{\pi^2}{4} \frac{3}{\tau_{ \mathrm{co} } \tau_{ \mathrm{ff} }^2 } \Gamma_1^2 \, \nabla_{\mathrm{ad}} \left[ 4 \nabla_{\mathrm{ad}} - ( \nabla_{\mathrm{ad}} \kappa^{}_T + \kappa^{}_P ) - \frac{4}{3 \Gamma_1} \right] & > & 0 \label{ZSVibSta} \end{eqnarray} % For a physical discussion of the stability criteria see Baker (\cite{baker}) or Cox (\cite{cox}). We observe that these criteria for dynamical, secular and vibrational stability, respectively, can be factorized into \begin{enumerate} \item a factor containing local timescales only, \item a factor containing only constitutive relations and their derivatives. \end{enumerate} The first factors, depending on only timescales, are positive by definition. The signs of the left hand sides of the inequalities~(\ref{ZSDynSta}), (\ref{ZSSecSta}) and (\ref{ZSVibSta}) therefore depend exclusively on the second factors containing the constitutive relations. Since they depend only on state variables, the stability criteria themselves are \emph{ functions of the thermodynamic state in the local zone}. The one-zone stability can therefore be determined from a simple equation of state, given for example, as a function of density and temperature. Once the microphysics, i.e.\ the thermodynamics and opacities (see Table~\ref{KapSou}), are specified (in practice % by specifying a chemical composition) the one-zone stability can be inferred if the thermodynamic state is specified. The zone -- or in other words the layer -- will be stable or unstable in whatever object it is imbedded as long as it satisfies the one-zone-model assumptions. Only the specific growth rates (depending upon the time scales) will be different for layers in different objects. %__________________________________________________ One column table \begin{table} \caption[]{Opacity sources} \label{KapSou} \[ \begin{array}{p{0.5\linewidth}l} \hline \noalign{\smallskip} Source & T / {[\mathrm{K}]} \\ \noalign{\smallskip} \hline \noalign{\smallskip} Yorke 1979, Yorke 1980a & \leq 1700^{\mathrm{a}} \\ % Yorke 1979, Yorke 1980a & \leq 1700 \\ Kr\"ugel 1971 & 1700 \leq T \leq 5000 \\ Cox \& Stewart 1969 & 5000 \leq \\ \noalign{\smallskip} \hline \end{array} \] \begin{list}{}{} \item[$^{\mathrm{a}}$] This is footnote a \end{list} \end{table} % % %___________________________________ Two column table (place early!) \begin{table*} \caption[]{Regions of secular instability} \label{TabSecInst} \vspace{4cm} \end{table*} We will now write down the sign (and therefore stability) determining parts of the left-hand sides of the inequalities (\ref{ZSDynSta}), (\ref{ZSSecSta}) and (\ref{ZSVibSta}) and thereby obtain \emph{stability equations of state}. The sign determining part of inequality~(\ref{ZSDynSta}) is $3\Gamma_1 - 4$ and it reduces to the criterion for dynamical stability \begin{equation} \Gamma_1 > \frac{4}{3}\,. \end{equation} Stability of the thermodynamical equilibrium demands \begin{equation} \chi^{}_\rho > 0, \;\; c_v > 0\, , \end{equation} and \begin{equation} \chi^{}_T > 0 \end{equation} holds for a wide range of physical situations. With \begin{eqnarray} \Gamma_3 - 1 = \frac{P}{\rho T} \frac{\chi^{}_T}{c_v}&>&0\\ \Gamma_1 = \chi_\rho^{} + \chi_T^{} (\Gamma_3 -1)&>&0\\ \nabla_{\mathrm{ad}} = \frac{\Gamma_3 - 1}{\Gamma_1} &>&0 \end{eqnarray} we find the sign determining terms in inequalities~(\ref{ZSSecSta}) and (\ref{ZSVibSta}) respectively and obtain the following form of the criteria for dynamical, secular and vibrational \emph{stability}, respectively: \begin{eqnarray} 3 \Gamma_1 - 4 =: S_{\mathrm{dyn}} > & 0 & \label{DynSta} \\ % \frac{ 1- 3/4 \chi^{}_\rho }{ \chi^{}_T } ( \kappa^{}_T - 4 ) + \kappa^{}_P + 1 =: S_{\mathrm{sec}} > & 0 & \label{SecSta} \\ % 4 \nabla_{\mathrm{ad}} - (\nabla_{\mathrm{ad}} \kappa^{}_T + \kappa^{}_P) - \frac{4}{3 \Gamma_1} =: S_{\mathrm{vib}} > & 0\,.& \label{VibSta} \end{eqnarray} The constitutive relations are to be evaluated for the unperturbed thermodynamic state (say $(\rho_0, T_0)$) of the zone. We see that the one-zone stability of the layer depends only on the constitutive relations $\Gamma_1$, $\nabla_{\mathrm{ad}}$, $\chi_T^{},\,\chi_\rho^{}$, $\kappa_P^{},\,\kappa_T^{}$. These depend only on the unperturbed thermodynamical state of the layer. Therefore the above relations define the one-zone-stability equations of state $S_{\mathrm{dyn}},\,S_{\mathrm{sec}}$ and $S_{\mathrm{vib}}$. See Fig.~\ref{FigVibStab} for a picture of $S_{\mathrm{vib}}$. Regions of secular instability are listed in Table~\ref{TabSecInst}. % % One column figure %----------------------------------------------------------- S_vib \begin{figure} \vspace{5cm} \caption[]{Vibrational stability equation of state $S_{\mathrm{vib}}(\lg e, \lg \rho)$. $>0$ means vibrational stability } \label{FigVibStab} \end{figure} % %______________________________________________________________ \section{Conclusions} \begin{enumerate} \item The conditions for the stability of static, radiative layers in gas spheres, as described by Baker's (\cite{baker}) standard one-zone model, can be expressed as stability equations of state. These stability equations of state depend only on the local thermodynamic state of the layer. \item If the constitutive relations -- equations of state and Rosseland mean opacities -- are specified, the stability equations of state can be evaluated without specifying properties of the layer. \item For solar composition gas the $\kappa$-mechanism is working in the regions of the ice and dust features in the opacities, the $\mathrm{H}_2$ dissociation and the combined H, first He ionization zone, as indicated by vibrational instability. These regions of instability are much larger in extent and degree of instability than the second He ionization zone that drives the Cephe{\"\i}d pulsations. \end{enumerate} \begin{acknowledgements} Part of this work was supported by the German \emph{Deut\-sche For\-schungs\-ge\-mein\-schaft, DFG\/} project number Ts~17/2--1. \end{acknowledgements} \begin{thebibliography}{} \bibitem[1966]{baker} Baker N., 1966, in: Stellar Evolution, eds.\ R. F. Stein, A. G. W. Cameron, Plenum, New York, p.\ 333 \bibitem[1988]{balluch} Balluch M., 1988, A\&A 200, 58 \bibitem[1980]{cox} Cox J. P., 1980, Theory of Stellar Pulsation, Princeton University Press, Princeton, p.\ 165 \bibitem[1969]{cox69} Cox A. N., Stewart J. N., 1969, Academia Nauk, Scientific Information 15, 1 \bibitem[1971]{kruegel} Kr\"ugel E., 1971, Der Rosselandsche Mittelwert bei tiefen Temperaturen, Diplom--Thesis, Univ.\ G\"ottingen \bibitem[1980]{mizuno} Mizuno H., 1980, Prog. Theor. Phys. 64, 544 \newpage % if processed with CM fonts, balances the references columns \bibitem[1987]{tscharnuter} Tscharnuter W. M., 1987, A\&A 188, 55 \bibitem[1989]{wuchterl} Wuchterl G., 1989, Zur Entstehung der Gasplaneten.\ Ku\-gel\-sym\-me\-tri\-sche Gas\-str\"o\-mun\-gen auf Pro\-to\-pla\-ne\-ten, Dissertation, Univ.\ Wien \bibitem[1979]{yorke79} Yorke H. W., 1979, A\&A 80, 215 \bibitem[1980a]{yorke80a} Yorke H. W., 1980a, A\&A 86, 286 \end{thebibliography} \end{document} Examples for figures using psfig and epsf respectively %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Figure up to 8.8cm (use the vertical space and the horizontal space to % place the figure in the right place, top and left of page) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \vspace{0cm} \hspace{0cm}\psfig{figure=994f9.ps,width=8.8cm} \vspace{0cm} \vspace{0cm} \hspace{0cm}\epsfxsize=8.8cm \epsfbox{file.ps} \vspace{0cm} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 2 figures side by side %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \vspace{0cm} \hbox{\hspace{0cm}\psfig{figure=994f9.ps,width=14.8cm}\hspace{0cm} \psfig{figure=994f9.ps,width=14.8cm}} \vspace{0cm} \vspace{0cm} \hbox{\hspace{0cm}\epsfxsize=7.5cm \epsfbox{file.ps} \epsfxsize=7.5cm \epsfbox{file.ps}} \vspace{0cm} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Figure with caption on the right side (0.5cm place between the figure % and the legend) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \vbox{\psfig{figure=file.ps,width=12.0cm}\vspace{-3cm}} \hfill\parbox[b]{5.5cm}{\caption[]{}} \vbox{\epsfxsize=12cm \epsfbox{file.ps}\vspace{-3cm}} \hfill\parbox[b]{5.5cm}{\caption[]{}} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % working with epsf.sty and rotate.sty to rotate figures (the figure % will be rotated out of the paper thats why you have to use big % hspaces and vspaces to put the figure in the right position) % but first try it with the psfig-macro (angle-command) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \psfig{figure=file.ps,width=8.8cm,angle=-90} \vbox{\vspace{-5.2cm}\hbox{\hspace{8.5cm}\epsfxsize=5.9cm \rotate[l]{\epsfbox{file.ps}}}} \vspace{5.4cm} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % to cut somthing from the top, bottom, right and left of the figure % change the Bounding Box of the figure with the folowing parameters % or use the command clip= (in epsf: \epsfclipon) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \psfig{figure=file.ps,width=8.8cm,bbllx=20pt,bblly=20pt,% bburx=365pt,bbury=567pt} \psfig{figure=file.ps,width=8.8cm,clip=} \epsfxsize=8.8cm \epsfbox[20 20 300 300]{aa2283.f1}