I've been working on trying to calculate atmospheres for a given equilibrium temperature $T_{eq}$ planetary mass $M_{planet}$ and radius $r_{planet}$, how would you calculate the maximum possible atmospheric pressure for that planet? I've found a link relating atmospheric pressure and planetary mass together, but I'm having trouble bridging the connection.
For ideal gases, $v_{esc} =6v_{therm} $, and $v_{therm}^2 = \frac{3RT}{m_{mol}}$, but I'm not understanding how to tie it all together.
I'm assuming spherical cows, and that forces besides simple atmospheric losses are not at play. I know that atmospheric losses can be a function of time, and am assuming a geologic timescale >3bn years.
Another page I found that I could not quite make sense of: How tall do atmosphere retaining walls on rotating space habitats need to be?
Edit:
I've been researching quite a bit since I've asked this question, and I've had some forward progress. I had looked at the barometric formula, but there was two different problems I couldn't resolve. The first is that that integral resolves to the same thing when $z=0$, which is $P_0 = \frac{G_\oplus M_{atm}}{A_\oplus}$. With $G_\oplus$ as the specific gravity for that body ($G_\oplus=9.81$ here on Earth). The other one is that the behavior is divergent, with $lim_{P_0\rightarrow \infty}=\infty$ which isn't matching my qualitative understanding of the problem.
So, I took my understanding of the problem, and tried to approach the problem under these assumptions:
- The maximum stable surface pressure would be where the outermost molecules will need 0 additional energy to escape orbit.
- For a gas volume at the surface, that volume would need to expend all its kinetic energy in order to reach escape velocity.
- A mass of air's total kinetic energy is a function of its temperature.
- I'm working with an ideal gas, and its behavior will be consistent.
My rephrased question became: "At that surface temperature, what pressure would have the same energy on the surface as an air mass at escape velocity?" And I started with 3 different formulas as givens: $v_a^2 = \frac{3RT}{M}\quad PV=nRT \quad v_{esc}^2=\frac{2G\oplus M_\oplus}{R_\oplus}$, with $R$ as the Boltzmann constant, $T$ as temperature, $n$ as number of atmospheric mols, $P$ as pressure and $V$ as volume.
I started by applying Boyle's law: $P_1 V_1=nRT_1 \quad P_2 V_2 = nRT_2$. Now, since we are exhausting all of the energy in an atmospheric volume, we can algebraically say $(P_1-P_2)V=nR(T_1-T_2)$. The volume stays constant, and the number of mols stays constant, which makes the only free variables pressure and temperature. My algebra (which doesn't really matter, since it led to an incorrect answer) led me to $\frac{2G_\oplus M_\oplus M}{3R_\oplus}=P_0$, where M was the molar mass. It felt really good, but when I started to test it with known planets, I was off by several orders of magnitude. My suspicion was that some nuance of fluids is at play, which is out of my field.
Back to the drawing board. After some searching, I came across a masters paper that was asking basically the same question. This is ultimately what I'm after, and I completely understand that the atmospheric mass flux, $\Phi$ is a function of time, and that $\frac{d\Phi}{dt}\neq 0$, but I can't identify the function of time. He mentions $F_{stel}(t)$ in relation to (12)
, but since that function does not have a time variable, I'm at a loss.
In case that paper gets deleted, its title is "Atmospheric Evolution on Low-gravity Waterworlds" by Constantin Arnscheidt. The function of import (I think) is: $$ F_{stel} = \frac{\big(\frac{r_{LW}}{r_s}\big)^2 F_{out}+(g_sr_s+L)\Phi}{\big(\frac{r_{SW}}{r_s}\big)^2 \frac{1}{4}(1-A)} $$