*
Cococubed.com

Abundance variable
second derivatives

Home

Astronomy research
Software instruments
Presentations
Illustrations
cococubed YouTube
Bicycle adventures
Public Outreach
Education materials
     Solar Systems Astronomy
     Energy in Everyday Life
     Geometry of Art and Nature
     Calculus
2023 ASU Solar Systems Astronomy
2023 ASU Energy in Everyday Life

AAS Journals
2023 AAS YouTube
2023 AAS Peer Review Workshops
2023 MESA VI
2023 MESA Marketplace
2023 MESA Classroom
2023 Neutrino Emission from Stars
2023 White Dwarfs & 12C(α,γ)16O
2022 Earendel, A Highly Magnified Star
2022 Black Hole Mass Spectrum
2022 MESA in Don't Look Up
2021 Bill Paxton, Tinsley Prize


Contact: F.X.Timmes
my one page vitae,
full vitae,
research statement, and
teaching statement.

$ \def\drvop#1{{\frac{d}{d{#1}}}} \def\drvf#1#2{{\frac{d{#1}}{d{#2}}}} \def\ddrvf#1#2{{\frac{d^2{#1}}{d{#2}^2}}} \def\partop#1{{\frac{\partial}{\partial {#1}}}} \def\ppartop#1{{\frac{\partial^2}{\partial {#1}^2}}} \def\partf#1#2{{\frac{\partial{#1}}{\partial{#2}}}} \def\ppartf#1#2{{\frac{\partial^2{#1}}{\partial{#2}^2}}} \def\mpartf#1#2#3{{\frac{\partial^2{#1}}{\partial{#2} \ {\partial{#3}}}}} $ A pdf of this note is avaliable.

One may ask why second derivatives are needed. If the equations being evolved contains derivative quantities, for example the $\partial e / \partial Y_i$ ``chemical potential'' term from the first law of thermodynamics, and if an implicit time integration is desirable, for example the system is stiff, then the Jacobian matrix will contain terms such as $\partial^2 e / \partial Y_i^2$.



Its been previously shown that the average of any quantity $\overline{\beta}$ by the number density $n_i$ weighted average \begin{equation} \overline{\beta} = \frac{\sum \beta_i Y_i}{\sum Y_i} \ , \label{eq:betabar} \end{equation} whose first partial derivative with respect to abundance $Y_i$ is \begin{equation} \frac{ \partial \overline{\beta}}{\partial Y_i} = \frac{\beta_i}{\sum Y_i} - \frac{\sum \beta_i Y_i}{\left ( \sum Y_i \right )^2} = \frac{\beta_i}{\sum Y_i} - \frac{\overline{\beta}}{\sum Y_i} = \frac{\beta_i - \overline{\beta}}{\sum Y_i} = \overline{\rm{A}} \ ( \beta_i - \overline{\beta} ) \ . \end{equation} The second partial derivative with respect to abundance $Y_i$ is then \begin{align} \ppartf{\overline{\beta}}{Y_i} & = \partop{Y_i} \left [ \frac{\beta_i}{\sum Y_i} - \frac{\sum \beta_i Y_i}{\left ( \sum Y_i \right )^2} \right ] \notag \\[8pt] & = -\frac{\beta_i}{(\sum Y_i)^2} - \frac{\beta_i}{(\sum Y_i)^2} + 2 \frac{\sum \beta_i Y_i}{\left ( \sum Y_i \right )^3} \notag \\[8pt] & = 2 \left ( \frac{\overline{\beta}}{(\sum Y_i)^2} - \frac{\beta_i}{(\sum Y_i)^2} \right ) \notag \\[8pt] & = 2 \overline{\rm{A}}^2 \ ( \overline{\beta} - \beta_i ) \notag \\[8pt] & = 2 \overline{\rm{A}} \ \frac{ \partial \overline{\beta}}{\partial Y_i} \ , \label{eq:azbar2nd} \end{align} which is a handy expression. Explicitly, \begin{align} \ppartf{\overline{{\rm A}}}{Y_i} & = 2 \overline{\rm{A}} \ \frac{ \partial \overline{{\rm A}}}{\partial Y_i} = - 2 \overline{\rm{A}}^3 \notag \\[8pt] \ppartf{\overline{{\rm Z}}}{Y_i} & = 2 \overline{\rm{A}} \ \frac{ \partial \overline{{\rm Z}}}{\partial Y_i} \end{align}


It's worth considering the general case for second full derivative as its not common. The differential operator \begin{equation} d = dx \partop{x} + dy \partop{y} \end{equation} applied to $f$ gives \begin{equation} df = dx \partf{f}{x} + dy \partf{f}{y} \end{equation} The second differential operator \begin{equation} d^2 = \left ( dx \partop{x} + dy \partop{y}\right ) \left ( dx \partop{x} + dy \partop{y}\right ) \end{equation} applied to $f$ gives \begin{align} d^2f & = \left ( dx \partop{x} + dy \partop{y}\right ) \left ( dx \partop{x} + dy \partop{y}\right ) f \notag \\[8pt] & = \left ( d^2x \ppartop{x} + d^2y \ppartop{y} + dx \ dy \partop{x} \ \partop{y} + dy \ dx \partop{y} \ \partop{x} \right ) f \end{align} If the partial derivatives commute such that \begin{equation} \mpartf{f}{x}{y} = \mpartf{f}{y}{x} \ , \end{equation} then \begin{equation} d^2f = d^2x \ \ppartf{f}{x} + d^2y \ \ppartf{f}{y} + 2 \ dx \ dy \mpartf{f}{x}{y} \ , \end{equation} and for an arbitrary quantity $z$ \begin{equation} \ddrvf{f}{z} = \ddrvf{x}{z} \ \ppartf{f}{x} + \ddrvf{y}{z} \ \ppartf{f}{y} + 2 \ \drvf{x}{z} \ \drvf{y}{z} \mpartf{f}{x}{y} \ . \label{eq:2ndfull} \end{equation}


For the case of composition variables, for an arbitray quantity $\alpha$, applying equation ($\ref{eq:2ndfull}$) yields \begin{equation} \ddrvf{\alpha}{Y_i} = \ddrvf{\overline{\rm{Z}}}{Y_i} \ \ppartf{\alpha}{\overline{\rm{Z}}} + \ddrvf{\overline{\rm{A}}}{Y_i} \ \ppartf{\alpha}{\overline{\rm{A}}} + 2 \ \drvf{\overline{\rm{Z}}}{Y_i} \ \drvf{\overline{\rm{A}}}{Y_i} \ \mpartf{\alpha}{\overline{\rm{Z}}}{\overline{\rm{A}}} \ . \end{equation} One assumes all partials of $\alpha$ with respect to $\overline{\rm{A}}$ and $\overline{\rm{Z}}$ are available from the physics is at hand (e.g., from an eos). The second partials of $\overline{\rm{A}}$ and $\overline{\rm{Z}}$ are given by equation ($\ref{eq:azbar2nd}$), and the first partials have been given previously.