So, habitable zone calculations are a pain. For a start, most habitable zone calculations make some key assumptions:
- The orbiting planet(s) is/are Earth-like, or at least similar to a more welcoming Venus or Mars.
- The orbits remain fully inside the habitable zone.
- The planets don't have any freaky axial tilts.
- The luminosity of the star remains constant (many calculations do look at how the habitable zone changes over a stars life, but variable stars could have much shorter oscillations in luminosity). Essentially, we want to orbit well-behaved stars.
- We want liquid water on the surface.
These assumptions do not cover all of the possible scenarios in which life could arise. For instance, they ignore the possibility of life on moons orbiting gas giants, where tidal forces could provide heat (hello, Europa and Enceladus!). They also imply that life must be carbon-based, using water as a solvent. Essentially, the term "circumstellar habitable zone" should really be "circumstellar this-seems-about-right-for-Earth-and-humans-to-live-don't-you-think zone".
The boundaries are also highly dependent on climate models, as we saw earlier - the section on Wikipedia detailing various predicted Solar System habitable zones should convince you of this. coming soonChoices of the four constants for $S_{eff}$, for instance, have dramatic effects on a planet, changing it from a Venusian hell to a cold Martian twin. ButModels from first-principles need to take into account, for instance, the greenhouse effect (radiative forcing, anyone?).
So, here's how to determine the habitable zone, in a nutshell:
- Choose the properties of your star at a given time - essentially, luminosity.
- Choose the physical properties of the sort of planet you want, early in its life. These include atmospheric composition, mass and radius (maybe), albedo, etc.
- Create models of the evolution of the planet depending on the incident stellar flux.
- Determine the range of fluxes in which such worlds can be habitable.
- Calculate the radii at which the stellar flux will take these values.
All of this, for the best models, is extraordinarily complicated. I thinkdon't know how to do most of it. However, we can look at one type of case which is really simple: the above answered your questionidealized greenhouse model. A simple derivation can be found here.
Let $T_s$ be the surface temperature and $T_a$ be the atmospheric temperature (assuming that both are roughly uniform across the planet). The planet itself has albedo $A$, and its atmosphere has an absorption constant $f$, which depends on its composition. The incident flux is $F_s$. The energy balance equation for the planet is
$$\frac{F_s(1-A)}{4}=(1-f)\sigma T_s^4+f\sigma T_a^4\tag{Planet}$$
The equation for the atmosphere is
$$f\sigma T_s^4=2f\sigma T_a^4\tag{Atmosphere}$$
Putting these two together yields
$$\frac{F_s(1-A)}{4}=(1-f)\sigma T_s^4+\frac{f}{2}\sigma T_s^4=\left(1-\frac{f}{2}\right)\sigma T_s^4$$
Rearranging, we get
$$T_s=\left[\frac{F_s(1-A)}{4\sigma\left(1-\frac{f}{2}\right)}\right]^{\frac{1}{4}}\quad\text{or}\quad F_s=\frac{4\sigma T_s^4\left(1-\frac{f}{2}\right)}{1-A}$$
The latter is probably more helpful to us if we want to find the boundaries of the habitable zone, although because it is so idealized, it still doesn't take into account more complicated effects like radiative forcing. I should also note that setting $f=0$ makes $T_s$ simple the planet's effective temperature. However, it is almost never true that $f=0$; on Earth, $f\simeq0.77$.
I suspect I may have gone into a little more detail than you needed. As I said at the beginning, you really only need those two equations to figure out the rough boundaries of the habitable zone. Still, I hope that the rest of this answer was a little - dare I make this pun - illuminating.