Naked Science Forum

Non Life Sciences => Chemistry => Topic started by: chiralSPO on 24/05/2022 18:32:34

Title: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 24/05/2022 18:32:34
This question is inspired by thinking about passive temperature control.
Please, don’t sidetrack this discussion with comments about alternative technologies, or putting on a sweater, or the pros/cons of a certain temperature, or anything about climate change. I want to keep the discussion focused on the chemistry/physics/engineering/math of optimizing a system as described in the next post (the remainder of this post is all background info pertaining to how I got to the question at hand—feel free to skip ahead)

A significant portion of our society’s energy usage goes towards heating and cooling our homes and places of business. There are many ways of passively controlling temperature inside building including insulation and reflective vs light-absorbing exteriors etc. But I want to focus this thread on materials/systems that have a mediating influence on temperature by virtue of their heat capacity (I believe “thermal mass” is the term used by architects).

Some substances have an extremely high heat capacity per unit volume. For example, liquid water; coming in at about 4.2 kJ•L–1 •K–1, with some minor variation across its liquid range. But this pales in comparison to the incredible latent heat of a phase change.

For example, t-butanol melts/freezes at 298 K (25 °C, or 77 °F), and has a latent heat of fusion of about 116 kJ•L–1 •K–1. Any room at equilibrium with a bottle of t-butanol would have a significant “thermal mass” crossing 298 K (25 °C, or 77 °F). A sufficiently large bottle of liquid t-butanol, with sufficiently fast heat exchange to the room would effectively prevent the room from dropping below 25 °C. Likewise, if the t-butanol were solid, the room could be held below this temperature. If you want a different temperature, you can pick a different substance, which will have a different melting point (and different latent heat of fusion)

The problem with this approach (ignoring issues specific to the substance of choice), is that it only works for a single temperature (very narrow temperature range). It probably wouldn’t be economically feasible to have such a large reservoir of this substance with heat exchangers etc. necessary to keep a room (or building) at a single temperature. And if the temperature at any point drifted far from the melting point, then the thermal mass of the substance would be very small compared to the latent heat of fusion.

One possible solution (no pun intended) would be to have a few different reservoirs containing different substances, with different critical temperatures. For example dmso melts at 19 °C (66 °F). So one could imagine a room that is kept between 19 °C and 25 °C with one reservoir of t-butanol and one of dmso. Or one could imagine using only one solvent, but with two (or more) reservoirs with different amounts of solute, resulting in slightly different melting points (by virtue of colligative properties).
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 24/05/2022 18:33:25
A more interesting (and potentially more useful) approach is to have an equilibrium mixture that has a broad working range, rather than a sharp phase transition.

Imagine a chemical reaction in equilibrium. For simplicity, we will just say there are two states, A and Z.

A 68468762664bf7f63435ea54ec87a726.gif Z

ΔGforward = ΔHforward – T• ΔSforward

For this reaction to work, entropy and enthalpy must be opposing each other ΔHforward and ΔSforward must both be positive (so as temperature increases, an endothermic reaction goes forward). Therefore, ΔHreverse and ΔSreverse must both be negative.

As the temperature changes, the equilibrium concentrations of the product/reactant states changes, absorbing heat as the temperature increases, and releasing heat as it decreases.

ΔG = –RTln([Z]/[A])
or
ΔH – T ΔS = –RTln([Z]/[A])

Some rearranging gives:
ΔH = T(ΔS – R ln([Z]/[A]))

In the real world, both ΔH and ΔS are also both functions of T, but let’s hold off on that for now, unless someone can provide a good reason why it must be considered. (I think it should be pretty subtle across the ranges of T we are trying to cover).

We can imagine the enthalpy of the reaction (ΔH) is what we are able to select. The more positive it is, the more energy is absorbed per molecule that is transformed from A to Z (and more released when going from Z to A). But if the equilibrium lies too far to one side, we will be limited in how much energy can be handled in one direction. (ie if 99.999% is in state A at a given temperature, then only a tiny amount of energy can possibly be released when the temperature falls). If we want the system to be able to both absorb and release significant amount of heat at a comfortable temperature, the larger the ΔH term is, the larger the ΔS term must be as well (so that T• ΔS is close to ΔH when [A] is close to [Z].

So the question becomes. For a given high and low temperatures (Thigh and Tlow), what is the optimal choice of ΔH and ΔS, such that the system is most stable between those to temperatures?[/i]

My algebra and calculus is a bit rusty. So, while my intuition indicates that there will be a very easy way to express these trade-offs mathematically, and find the optimum values, I have not been able to convince myself that the answers I am finding are correct.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Halc on 24/05/2022 19:19:40
But this pales in comparison to the incredible latent heat of a phase change.
I don't get it. Suppose I have a material that melts/freezes at room temperature. This only works once, and then it's done. Say I want to heat my building in the winter.  I have liquid 'stuff' that freezes as the room temp drops just below where I want it, so it keeps the room warm until it's entirely frozen. Now what? How am I going to get it into liquid state again? I have to turn the heater on and it has all the much more work to do since it has to melt all this nice stuff on top of actually heating the place. It seems I've saved no energy at all, so I'm not sure what you're getting at.
Heating/cooling is all about insulation, not thermal capacity. The more thermal energy that passes from the hot side to the cold side, the more energy it takes to put it back.

Industry, the primary consumer of resources, seems not to care. In the middle of winter I watched the power consumed by the air conditioners in the computer lab. All it needed was a fresh air fan on the roof since it was well below freezing outside, and there they are pumping heat out of the lab to the radiator on the roof, and not even into the heating system keeping the offices warm.
Another building (built for IBM) had the heater break down on an August day.  We had the doors/windows open and still had to wear winter coats because there was no heat to mix with the cold system. Temp was set just like water in houses: by mixing just the right amount of hot and cold, and not just turning off the whole system when it was cool enough. Apparently the utility bill was of no concern.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 24/05/2022 19:47:57
I don't get it. Suppose I have a material that melts/freezes at room temperature. This only works once, and then it's done. Say I want to heat my building in the winter.  I have liquid 'stuff' that freezes as the room temp drops just below where I want it, so it keeps the room warm until it's entirely frozen. Now what? How am I going to get it into liquid state again? I have to turn the heater on and it has all the much more work to do since it has to melt all this nice stuff on top of actually heating the place. It seems I've saved no energy at all, so I'm not sure what you're getting at.
Heating/cooling is all about insulation, not thermal capacity. The more thermal energy that passes from the hot side to the cold side, the more energy it takes to put it back.
Sure, the use of a phase-change material would not solve this particular problem. Imagine a climate (or time of year), when the temperature changes throughout the day, oscillating between a low of 15 °C and a high of 30 °C (≈59 °F to ≈86 °F), and you would rather keep the house at 22 °C (≈72.5 °F) with minimal energy usage.
Industry, the primary consumer of resources, seems not to care. In the middle of winter I watched the power consumed by the air conditioners in the computer lab. All it needed was a fresh air fan on the roof since it was well below freezing outside, and there they are pumping heat out of the lab to the radiator on the roof, and not even into the heating system keeping the offices warm.
Another building (built for IBM) had the heater break down on an August day.  We had the doors/windows open and still had to wear winter coats because there was no heat to mix with the cold system. Temp was set just like water in houses: by mixing just the right amount of hot and cold, and not just turning off the whole system when it was cool enough. Apparently the utility bill was of no concern.
They should care. But, I would like to remind folks that the point of this thread is not to answer problems of economics, politics, or even practicality. I want to know how to express the relationships between ΔH and ΔS and working temperature range of such an equilibrium. Simple thermodynamics/algebra/calculus.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: evan_au on 24/05/2022 23:11:59
Quote from: Halc
I don't get it.
Another on the sideline of practicality: It could make sense in a desert, where it gets well above 30C almost every day (including winter), and well below 15C almost every night (even in summer).

You would need two tubs of your phase change material which can be selectively exposed to outside air or inside air.
- During the day, you expose one tub to the outside air to heat it, while the other (cold) tub is exposed to inside air to keep the house cool.
- Reverse it overnight: Expose the hot tub to the inside to heat the house, while the other tub is exposed to the outside air to cool it, ready for the next day.

Effectively, you want a system that will produce a 12-hour phase change* in the outside temperature cycle.
     * This is a temporal phase change, achieved by a physical phase change

For it to work well, your house would need to be well-insulated, and you would need fans to circulate air past the tubs of hot/cold materials.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 25/05/2022 03:19:48
Effectively, you want a system that will produce a 12-hour phase change* in the outside temperature cycle.
     * This is a temporal phase change, achieved by a physical phase change
I see what you did there :-)
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: paul cotter on 25/05/2022 12:40:33
I understand the original idea as using an arbitrarily large quantity of phase change material to store energy and release it only at a transition temperature. I see no problem in principle other than the vast quantity of material needed. If going for a chemical change process things become more difficult as a lot of processes that are thermodynamically favourable do not proceed for kinetic reasons. With exothermic processes there is a risk of positive feedback leading to a runaway.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 25/05/2022 13:17:42
You would need two tubs of your phase change material which can be selectively exposed to outside air or inside air.
Just make the walls out of tanks of the material.

The idea works, but the amount of stuff you need is huge unless you have really good insulation.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 25/05/2022 23:03:16
I understand the original idea as using an arbitrarily large quantity of phase change material to store energy and release it only at a transition temperature. I see no problem in principle other than the vast quantity of material needed. If going for a chemical change process things become more difficult as a lot of processes that are thermodynamically favourable do not proceed for kinetic reasons. With exothermic processes there is a risk of positive feedback leading to a runaway.

Thankfully, any chemical equilibrium that satisfies the reversibility requirements stipulated are not going to be at risk of thermal runaway.

Also, while it is true that favorable thermodynamics does not require practical kinetics, for now, let's leave that wrinkle out. I am envisioning very simple, fast, and reversible reactions such as proton transfer reactions—the equilibrium of protonation states of weak acids and bases, and their conjugate bases and acids.

A-H + B 68468762664bf7f63435ea54ec87a726.gif A + B-H+

One can think of controlling the ΔH by selecting A-H and B with appropriate relative proton affinities, and controlling the ΔS by selecting the ratio of A-H and B. In the extreme, we can think about one of them essentially being the solvent, for example a weak acid in water, (A-H +H2O 68468762664bf7f63435ea54ec87a726.gif A + H3O+) where the hotter it is, the more the acid dissociates, pushing the reaction to the right).

But what I am most interested in at the moment, is just the math involved in figuring out how to maximize the ability of the solution to absorb/release energy over a desired temperature range.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 25/05/2022 23:55:43
Hi.

This question is inspired
   Well, it is quite a good idea.

ΔG = –RTln([Z]/[A])
    Could you clarify this please?   I'm not sure what your ΔG is,  is it actually ΔG° ?   Are  [Z] and [A] concentrations at equillibirum only?   i.d.k.

This is the conventional equation:
    ΔG =   RT   Ln (Q/K) 
Where ΔG = Gibbs free energy change for the system, (in the forward direction and at the specified concentrations).
Q = quotient of concentrations of products / reactants =  [Z] / [A]
K = chemical equilibrium constant =  Quotient as above but AT EQUILIBRIUM.
   
    Just to clarify this,  this ΔG is a function of 3 variables:   The temperature, T,  and the concentrations [Z] and [A].

[reference:  https://chem.libretexts.org/Courses/Grand_Rapids_Community_College/CHM_120_-_Survey_of_General_Chemistry/7%3A_Equilibrium_and_Thermodynamics/7.11%3A_Gibbs_Free_Energy_and_Equilibrium ]

   There seems to be a K missing in your expression,    much as if you were assuming K always = 1.   
This could be enough to stop your idea working completely.    If  K = 1 always, then the net reaction never shifts forward or backward -  the equillibrium point remains with equal concentrations of products and reactants [Z] = [A]  regardless of what happens.   In this way it won't respond to changes in temperature at all.   
    To re-phrase this  K ≠1.  It is essential that K = K(T) = some function of temperature.

     Using conventional theory,   it seems that we can approximate  K(T) = equillibrium constant at temperature T  as
K(T)     ≈   e -(ΔG°/ RT )

    This quantity, ΔG°  is not a function of the concentrations of the products and reactants.   At most it is a function of the temperature, T, but more usually the temperature and pressure are also assumed to be standard temp. and pressure.   Since you're interested in changes occurring around room temp. and pressure, it shouldn't be a problem to assume  ΔG°  is just a constant  which you can find in a book for the reaction A → Z.

   Anyway, re-arranging that equation we obtain:   ΔG°  = -RT Ln (K)  =   -RT Ln ([Z]/[A])     where  [Z] and [A] are now only to be taken as the concentrations at equillibrium.     That might have been the equation you were suggesting in your original post.  It matters a lot because, if that was what you were doing,  then when you re-arranged it to find ΔH  I don't think it was the ΔH that you were actually hoping or thinking you'd find.

   Summary:  Sorry that was confusing.  I'm confused and just trying to match up your notation with that used in some other texts on the subject.  I need you to check or explain what it was you were hoping to suggest with your formula  ΔG = -RT Ln ([Z] / [A])  . 

Best Wishes.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: alancalverd on 26/05/2022 06:56:43
The idea is not entirely impractical.

As BC says, it does depend on having really good insulation but if you know your maximum insolation, minimum outdoor temperature, conductivity of your exterior cladding, and internal heat load (say 100W per person plus whatever other stuff you have in the building) you can calculate the mass of phase-change material required to keep the interior wall close to the melting point provided the night temperature is below your target for long enough.

And a practical thought: whilst heat storage systems have been proposed and implemented around melting, it's worth noting that the latent heat of vaporisation of water is about 7 times that of fusion, and the transition temperature can be relatively easily controlled by changing the pressure in the containment vessel. Worth looking at some more volatile materials?   
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 26/05/2022 08:54:32
Trouble with vapourisations is the very low density of vapours- you need big tanks.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: alancalverd on 26/05/2022 10:54:42
But you can start with a tank of liquid and a fairly small head space. LV decreases with increasing pressure but may remain practicable up to 1000 bar or so depending on the critical point of the liquid.

Or consider a sublimating solid?
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 26/05/2022 16:54:23
  Could you clarify this please?   I'm not sure what your ΔG is,  is it actually ΔG° ?   Are  [Z] and [A] concentrations at equillibirum only?   i.d.k.

Thank you for engaging my intended line of inquiry and trying to clarify points of confusion :-)

I am definitely not talking about ΔG° because I am trying to compare different temperatures and concentrations.

I think that [A] and [Z] are referring only to equilibrium concentrations (at each of the temperatures under consideration), so I am thinking about K, rather than Q. In the simplest case, A 68468762664bf7f63435ea54ec87a726.gif Z,  Keq = [Z]/[A], though for an acid/base case (which I mentioned in another post), B + A-H+ 68468762664bf7f63435ea54ec87a726.gif B-H + A, it could be Keq = ([B][A-H+])/([B-H][A]). Another option would be a dissociation, A-B 68468762664bf7f63435ea54ec87a726.gif A+ + B, where Keq = [A+][B]/[A-B]
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 26/05/2022 19:43:30
But you can start with a tank of liquid and a fairly small head space. LV decreases with increasing pressure but may remain practicable up to 1000 bar or so depending on the critical point of the liquid.

And you can watch the heat of vaporisation crash practically to zero (to exactly zero at  T Crit.)

https://en.wikipedia.org/wiki/Enthalpy_of_vaporization
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: alancalverd on 26/05/2022 22:33:12
All control systems have limits, which is why we need to use numbers to design them.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 26/05/2022 23:44:36
All control systems have limits, which is why we need to use numbers to design them.
Sometimes the number you need is "about a hundred or a thousand" which is the right ballpark for how much denser a solid is, compared to the vapour near atmospheric pressure.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 27/05/2022 04:40:39
Hi again.

   I've had another look at the mathematics and tried a few approaches.  I've got something that looks reasonable but it's taking too long to write it down and mark up the equations with LaTex.    I'll try and get it done but I'm off to the big city tomorrow and I know I won't be writing much here for at least a day.

    The final results I can write down.   They are dissapointing and probably only tell you what you should have already been able to guess.   On the other hand... they make sense and that suggests I haven't gone too far wrong.   Glass half-full or half-empty.

1.   You want to have as large an amount of Substance A (or equivalently of substance Z) as possible.    i.e.  Whatever you do, do lots of it -  make the biggest system you can.

2.    You want a reaction with as large a ΔH  as you can get.   Infinity would be great but otherwise just as large as you can get.

3.    What about ΔS?    If you want maximum response from the system at room temperature   (say about 300 Kelvin).   Then you will need ΔS  ≈  ΔH / 300.       More generally the system will act most effectively as a buffer for the room temperature when  ΔS  =   (whatever your value of ΔH is)   divided by (whatever value of T you're trying to hold the room at).   Changing the value of ΔS from that value has a rapid and detrimental effect:   The system becomes almost non-changing (and hence unable to buffer the temperature) at your desired room temperature.   Instead it works well as a buffer at a different temperature   (which will be precisely T = ΔH / ΔS ).

    You could have guessed these couldn't you?   The mathematics will show it but it's going to be a day or more.

Best Wishes.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 27/05/2022 09:00:21
While it's true that ΔS  ≈  ΔH / 300 is pretty much the criterion for optimal regulation at 300K, it is just restating the requirement that the melting point is near 300.

And tables of melting points are easier to find than tables of ΔS

It's also worth thinking about the "units".
As Alan didn't notice,  ΔH is typically given per mole or per unit mass.
On a practical basis you might want to look at the value per cubic metre.
And then you need to look at   ΔH per £.

 

Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: alancalverd on 27/05/2022 10:39:35
Sometimes the number you need is "about a hundred or a thousand" which is the right ballpark for how much denser a solid is, compared to the vapour near atmospheric pressure.
Which is why I suggested using a liquid. Or a sublimating solid.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 27/05/2022 13:36:43
Sometimes the number you need is "about a hundred or a thousand" which is the right ballpark for how much denser a solid is, compared to the vapour near atmospheric pressure.
Which is why I suggested using a liquid. Or a sublimating solid.
Ok, because the vapour takes up an impractical amount of space you suggested turning a solid or liquid into the (impractical)  vapour.
Were you aiming to fail?
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 30/05/2022 15:25:22
Hi again.

I'll try and get it done but I'm off to the big city tomorrow and I know I won't be writing much here for at least a day.
    I'm really sorry but I'm not going to get this done any time soon.

   I'm already on pages of work.   Therefore no-one will want to read it, not even @chiralSPO .
I do not have the time to complete.  I'm way behind on several other important tasks.  The whole thing is now begining to cause me some stress and that's not how a forum should be. 
   I'm very sorry @chiralSPO  but you should not wait for a detailed response from me, it would safer to assume it won't happen.

Apologies,  Eternal Student.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 31/05/2022 20:49:04
Thanks to Eternal Student and Bored chemist for trying to answer the questions I posed.

Hi again.

I'll try and get it done but I'm off to the big city tomorrow and I know I won't be writing much here for at least a day.
    I'm really sorry but I'm not going to get this done any time soon.

   I'm already on pages of work.   Therefore no-one will want to read it, not even @chiralSPO .
I do not have the time to complete.  I'm way behind on several other important tasks.  The whole thing is now begining to cause me some stress and that's not how a forum should be. 
   I'm very sorry @chiralSPO  but you should not wait for a detailed response from me, it would safer to assume it won't happen.

Apologies,  Eternal Student.

Thanks for trying. I do not mean to give anybody homework to stress over.

I too have gotten bogged down in the math, which is why I reached out.

I have gotten somewhere:

Assuming a two state model (A 68468762664bf7f63435ea54ec87a726.gif Z), the equilibrium ratio of concentrations is:

[Z]T/[A]T = e–RT/(ΔΗ–TΔS)

Because we want ΔH–TΔS to be 0 at comfortable temperatures, that restricts the range of ΔH as a function of ΔS (or vice versa). The bigger ΔH is, the more thermal energy can be absorbed per mole of transformed molecule as T increases (and vice versa). But due to the exponential nature of the function providing [Z]/[A], there will be only a very narrow range in which there is an appreciable amount of both A and Z.

Because one is turning into the other, at a given temperature, T1, [Z]T1 = 1–[A]T1
KT1= (1–[A]T1)/[A]T1
And then at a different temperature, T2, KT2= (x+1–[A]T1)/([A]T1–x), where x is how many moles is getting converted.

amount of energy absorbed when going from T1 to T2 = ΔH•x

And then I get lost in the algebraic rearrangements when trying to get ΔH•x as a function of T1 and T2. (which then means I cannot use calc to optimize it).


I may keep banging my head against this wall for a while (hints or answers would still be appreciated). But I think the long and short of it is: there isn't really much advantage to trying for a gradual transition--it doesn't change the overall amount of energy that can be absorbed (ΔH). So, as Bored chemist says, maximizing ΔH/£ (or ΔH/$) for something with a comfortable equilibrium temperature is probably the way to go, and then just use as much as needed to keep the temp where you want it.

 :)
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 31/05/2022 22:02:50
It's worth remembering that, for decades, we used steam radiators The surface temperature of them was near 100 C.
But we used some sort of thermostat to maintain the temperature near 20C .
You can get away with the equilibrium temperature being a little "off" the temperature you want to maintain.

Having said that, using something that melts near 21C is a good starting point..
And if you have that, you automatically have the right value of delta S for the material's delta H.

So, as you say  what you want is something cheap  that "melts" near 21C and has a large delta H per $
I think soft paraffin wax is likely to be a good candidate.
The other thing you need to consider is thermal conductivity. Wax isn't so good at that.
 
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 01/06/2022 02:37:00
Hi again.

   @chiralSPO  - this is what I had started writing but stopped.   It's rough, not spell-checked and just not finished.   However, if it helps you,  then you can see it.     You'll realise why I stopped.... it's too long and I started by explaining things too simply, as if other users might also be reading it.

    Actually it just is too long....  the forum won't allow me to post anything with too many characters.....   I'll split it into two chunks.   I'll also put it under a "spoiler" so it won't clog up all the space for people who are reading this thread.
- - -  - - - - - - - - -

Chunk 1
Spoiler: show

But what I am most interested in at the moment, is just the math involved in figuring out how to maximize the ability of the solution to absorb/release energy over a desired temperature range.
    Let's minimise peoples disinterest and just do something for the simple situation.
   
You have the simple chemical reaction forming an equilibrium   A 68468762664bf7f63435ea54ec87a726.gif Z

   Since the reaction is 1 molecule  to/from  1 other molecule,  we have

NA   + NZ  =  N   = a constant at all times 
[Eqn 1]
   with  NA = number of A molecules;   NZ = no. of Z molecules;   N = total number of molecules.

We're going to tackle this situation with the equilibrium constant.  I've considered alternative strategies but this one is reasonably easy, makes good progress and produces answers that seem to agree with what you might have expected.

K  =   [Z] / [A]   =  ratio of concentrations at equilibrium.   
[Eqn 2]

      The way you've described it we can assume the reactants and products are in the same container,  possibly in some solvent....  Whatever details apply - they slosh around in or contribute to some fluid of volume V.   The volume might change slightly as products and reactants come and go, or temperature causes expansion and contraction, that won't matter.  All that matters is that at any instant, the products and reactants are in the same volume.   So we get some cancellation and can replace concentrations with numbers of molecules in the equilibrium constant.

     Hence,    K  =  NZ  / NA.     
[Eqn 3]

Combine [Eqn 1]  and [Eqn 3]  to obtain,

    NA  =  a867215751a34630a3d87f730593a0da.gif   
[Eqn 4]
That's a decent start,  we know the number of A molecules that should be present at any time.   Note that the equilibrium constant K  is a not a constant, it changes depending on various parameters.

   The balance of the reaction equilibrium will be entirely determined by K.   We take this expression for K
    b4e3f1d47db6a5dfe91396a94dadda96.gif
[Eqn 5]
A reference for this result was given in a previous post.   Here's a slightly different one:  https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Supplemental_Modules_(Physical_and_Theoretical_Chemistry)/Thermodynamics/Fundamentals_of_Thermodynamics/Free_Energy_and_Equilibrium
     This is where we will make contact with your (@chiralSPO) ideas of varying the enthalpy and entropy.  We haven't forgotten about it.  All the information and contact we shall need will come through the equilibrium constant.  That constant tells us how the equilibrium shifts. 
     We're going to need a moment to explain this formula [Eqn 5].  As already mentioned in an earlier post, some texts would have used ΔG° in place of ΔG in that formula and everyone has a slightly different view on what ΔG° is supposed to be.   In the worst case, some texts insist that ΔG° has been determined at a standard reference temperature, in that situation ΔG° is nothing more than a fixed value you can find in a book, a constant regardless of whatever else is happening.
    Now we MUST HAVE a quantity,  ΔG,   which depends on temperature.   Consider for a moment that it doesn't but instead we are using a quantity which I'll call ΔG" and which is just a constant regardless of whatever else is happening.  An equilibrium with K described by [Eqn 3] achieves a perfect balance with an equal amount of product and reactant   if and only if     K = 1.     Next use [Eqn 5] with ΔG = ΔG" = a fixed constant:
    K =  e (-ΔG"/RT).   
[Eqn 6]
So   K = 1   if and only if    one of the following holds:   (1)  Temperature T →∞  ,  or else  (2) ΔG" = 0.   We're not interested in infinite temperatures, so we would want ΔG" = 0.   However ΔG" was assumed to be a constant, so using  [Eqn 6] again,  we have   K = e0 =  1   regardless of temperature.    The equilibrium never changes, it doesn't respond to temperature at all.   We've discussed this before - it's no use.
   The approximation where a fixed standard value ΔG" was used in [Eqn 5] is acceptable for many situations especially where you are working around standard temperatures AND  ΔG" is not approximately 0.   However, that is not our situation.  We are very likely to want an equilibrium operating around the perfect balance point with K ≈ 1 and hence ΔG" ≈ 0.    We must have a better approximation for the equilibrium constant.
 
   Here's [Eqn 5] again:
b4e3f1d47db6a5dfe91396a94dadda96.gif
[Eqn 5]
That's a perfectly good expression where ΔG  is allowed to be determined by temperature.   Many texts will still call the value that will be inserted into that formula ΔG° but it is fully assumed to be dependant on temperature.  It is just the free energy change that occurs for the forward reaction  A → Z  at a given temperature.   (Further discussion of this removed, this post is already too long).
Anyway, this ΔG is calculated as follows:
     ΔG =  ΔH - TΔS
[Eqn 7]

   As @chiralSPO  mentioned,   ΔH   and  ΔS  do have some temperature dependence but over the range of temperatures we are considering it's not important.   ΔH ≈ ΔH° = constant   and  ΔS ≈ ΔS° = constant.     We obtain,
    ΔG°  =   ΔG°(T) =     ΔH° - TΔS°
[Eqn 8]
    Note that this ΔG°  does retain some temperature dependence and it is sufficient for our purposes,   the quantity T is a variable and appears as the multiple of ΔS°.   (It's only the more extreme definitions of ΔG° where the temperature would have been locked at one value and become what I previously described as ΔG" ).

   I appreciate we might lose some people here.  It would help to get a diagram to show what we're doing:

Sketch to show  ΔG°   and  K  as functions of T  given by Eqn 8 and Eqn 5

* graph2.png (19.6 kB . 1495x935 - viewed 2326 times)

Features of the sketch:    K = 0 at T=0 kelvin.      K  → finite value = eΔS°/R  as T →∞.   It has a point of inflection around  T=ΔH° / 2R.     
   ΔG°  is just a straight line,  with  ΔG° = 0  at  temp. T0 = ΔH°/ΔS°   (easily obtained from Eqn 8 ).
   The point of inflection for K(T) is not too important but is easily obtained from the expression   K = e-ΔG°/RT  =  (eΔS°/R) .  e-ΔH°/RT    =    (constant) . e-(ΔH°/RT)   and differentiating twice with respect to  T.   
     The most important characteristic is that K passes through 1  at T = T0  this is where the equilibrium would have an equal amount of product and reactant.


---End of Chunk 1 ---
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 01/06/2022 02:39:40
Chunk 2
Spoiler: show

So the question becomes. For a given high and low temperatures (Thigh and Tlow), what is the optimal choice of ΔH and ΔS, such that the system is most stable between those to temperatures?
    Stability is not something you've decleared very well.   I'm going to interpret the notion of "stability" as implying that the reaction would absorb (or release) as much heat as possible to counter a temperature change in the room.
   Let's just spend a moment to see that this is a sensible approach:

Diagram 2:   How the room loses / gains thermal energy

   Thermal Energy change
   in the Room,  cΔT               ----> heat lost to the outside world, ΔE   (due to... poor insulation etc.)
       /\
        |
        | 
  Heat released by chemical equillibrium system, ΔU


    Under the situation shown in the above diagram we have     
Change in thermal energy of the room  =   c.ΔT  =    ΔU - ΔE                 
[eqn 10]
  with c = heat capacity of the room;   ΔE = energy transferred out of the room;   ΔU = heat released by the chemical system into the room.
   Note the sign conventions we have established:    Room cooling: ΔT<0   ;  ΔU >0 the chemical system will release heat (thermal energy) into the room  ;  ΔE > 0  heat is lost to the outside world.    Let's not spend too long discussing this, it should be apparent that ΔE can be considered as the sum of natural heat exchange to the world outside the room + energy exchanged due to heating or cooling systems like radiators or air-conditioning units that a person has inside the room.  A negative value indicates a net heat input to the room (e.g. from the house central heating system).
    Note also that the chemical system is considered as part of the contents of the room.   In the most degenerate situation the room is only the chemical system but usually the room is the chemical system + the rest of the room (like walls and furniture) such that the chemical system is just a minor contribution to the mass and heat capacity for the room.   (LATE EDITING:  Stuff removed - the above is just conservation of energy ideas and doesn't need further comment).

   The following results follow from [Eqn 10]
(i) ΔU ≠ ΔE  iff the room changes temperature.
(ii)  .. result removed in late editing... serious error here, sorry.   The following result still holds but for different reasons...
   Therefore,  |ΔU| ≤ |ΔE| always and equality holds if and only if ΔT =0.   
[Result 11]


   We can now see that no chemcal equillibrium system A↔Z   exists which would perfectly maintain a fixed temperature in the room.    Assume some energy ΔE≠0 is removed from (or added to) the room.  For ΔT = 0, we must have ΔU = ΔE  but  the chemical equillibrium would not change at all if the temperature didn't change, so we would have ΔU = 0.    Therefore, when some energy enters or leaves the room we are always in the situation with |ΔU| < |ΔE| and the best we can hope for is that the chemical equillibrium opposes  ΔE as best it can (so that the temperature change will be less severe than if the chemical equillibrium system wasn't there).

    I would suggest that the best interpretation of "maximum stability" for our temperature buffering system will be trying to keep |ΔT| as low as possible for a given amount of ΔE (energy transferred to/from something outside the room).  For our purposes, the quantity ΔE is beyond our control or influence, that requires someone to insulate the room and turning on/off the central heating system of the house etc.

    From [Eqn 10] ; [Results 11] and recognising that under our sign conventions ΔU opposes the sign of ΔT we have,
   |ΔT| =  5c8a579dc12db701a9a839324727b38e.gif
[Equation 12]
So, that minimising the size of the temperature change for the room due to a given Energy transfer into / out of the room, requires the size of |ΔU| to be MAXIMISED.

     CLAIM  that ΔU = heat released from the chemical system = function of the current Temperature T, the parameters ΔH° and ΔS° of the reaction A→B   and, of course, of the temperature change ΔT that occurrs.  PROOF: We have the following,
   ΔU =  ΔNA . ΔH° / (Avogadro)   
[Equation 13]
 
   Where, ΔNA = the change in NA (change in the number of molecules of A);     ΔH° is the specific enthalpy of the reaction change (heat change for 1 mol of A being converted to Z);     (Avogadro) = Avogadros number.   The signs work (I hope),  recall ΔNA is negative, NA decreases, when molecules change from A to Z   and  ΔH° = negative when  A→Z is exothermic,  making  ΔU = heat released to the room = positive.
    This is taking too long.... let's just speed this up.... we can obtain an expression for ΔNA   from  [Equation 4] in terms of the change in the equillibrium constant K.   We have a thermodynamic expression for K from [equation 5].   That establishes the "CLAIM" that ΔU is a function of the things mentioned above.  The dependance on ΔT (the change in temperature) came from ΔNA but that will be clearer in a moment...

    During the actual operation of the chemical equillibrium device,  you cannot change the parameters ΔH° or ΔS°  of the reaction or add extra molecules of A,  or do very much to change the way the device can abosrb or release heat.   The ony changes that occurr are temperature changes in the room.  So in calculating the quantity ΔNA during the actual operation of the device,   we can treat all other parameters as fixed constants.   This saves a mess writing derivatives as partial derivatives,  for example  K = K(T) = a function only of Temperature during actual operation,   even though  K has dependance on ΔH and ΔS  looking at [eqn 5].   This should also help to explain why finding partial derivatives w.r.t.  variables like  ΔH°   is not something we will do.   Just to be clear on this:  It might seem that since we are trying to maximise ΔU  (amount of heat released by the chemical system) with a carefull selection of parameters like ΔH° for the reaction, then we would want to find derivatives w.r.t. ΔH° and treat the problem as a typical maximisation problem.   However, that is not what we will do.   Instead, we treat the reaction parameters as a fixed constant (don't specify the actual value - just recognise that it is a constant) and see how the system will naturally behave in operation.   After that, we are free to see how those results change or what they will mean when the reaction parameters like ΔH° and ΔS° are given a particular numerical value.   This method seems to work better than just directly trying to find maxima from derivatives w.r.t.  reaction parameters.   (For example, there doesn't seem to be any local maximum.   This might be discussed later -  the larger ΔH° can be the better, the progression is endless and no local maxima exists  BUT also  you just don't get much useful information about the system anyway).

   Let's see [Eqn 13] again:
δU =  δNA . ΔH° / (Avogadro)   
[Equation 13' ]
     We'll drop the Avogadros number (as if ΔH° was stated per molecule not per mole), it only scales the value of ΔU,  the maximum won't be affected without it anyway.   Note that we are now considering only small changes, as marked with δ symbols -   the quantity  ΔH°  remains with a big delta because it is just a reaction parameter, a constant and not something that changes during the operation of the device.  We'll also just write ΔH instead of ΔH°  because it's shorter and should be clear enough that this is a fixed constant, a reaction parameter not something that changes dynamically during operation.
   From [Eqn 4] we can differentiate, noting that temperature is the only variable that changes during operation of the device:
   δNA   =    468a0c3e5b20bec18b49e354e8569341.gif
K and the derivative dK/dT can be found from [Eqn 5] and [Eqn 7]
  b4e3f1d47db6a5dfe91396a94dadda96.gif  =  c1ad89a28eacefddbf42653d9a4e2e73.gif

   b32640fae9cdd0cc5c8e0c33070b3eef.gif

[Check the signs work here:  If ΔH >0,  then A→Z is endothermic.  K = [Z]/[A]  and we expect K to increase when T increases.  dK/dT >0  is what we have.  Next up δNA has a negative sign out at the front, so it's negative and yes we have A→Z so the number of molecules of A would decrease.   This all looks OK to me. ]

  Substitute this into [Eqn 13' ] and we have:

δU =   9b01b4c318b24e8cd8e8b4287531115e.gif   
[Eqn 14]
   where  ζ =  eΔS/R  =  a convenient constant to hold all the entropy information.   Note that a lot of the quantities appearing in {Eqn 14] are constants,    T is the only variable.

Recall that we wish to maximise δU  given a fixed transfer of energy δE from the room.  Straight away from [Eqn 14]  we can see that  if N (the total number of molecules of A and Z) is increased then δU is increased.   There is also a term on the numerator  ~  (ΔH)2   and although ΔH does appear elsewhere its effect is most noticeable here. *LATE EDITING:  Reasonably true for small to medium ΔH.   However, when ΔH becomes very large, then the exponential function e-ΔH will dominate and bring the numerator towards 0.   This effect is more important than I first thought because ΔH usually is large - hundreds of thousands in standard units of Joules/mol. * 
   We can tidy up our expression a little.   We are trying to maximise δU  for a fixed value δE of energy transferred out of (or into) the room.   However, we can imagine the maximisation requirement differently.   Assume that a small change in temperature (δT) of the room has happened.   We would like to have this cause a release (or absorption) of energy by the chemical system  (δU)  that is as large as possible because that would have forced δE to be as large as possible.   [Eqn 12]  can be used to formalise this condition and establish that it is equivalent.  The post is too long, I'm skipping that.  We will just re-state the maximisation problem as required:   To buffer the temperature of the room as much as possible,  we wish to maximise δU  for a given fixed change in temperature δT.    Anyway, with that modification we can divide by δT  and obtain the following result (where δT→0 to form the derivative on the R.H.S.)

3c6c95c1fbe188a08312bad221eb9acc.gif   
[Eqn 15]
     Our objective is to maximise this quantity.

.....  I stopped writing here.... there are some derivatives and some graphical plots to show what is happening that i would have inlcuded.... but by now   no-one is reading this   and the whole set of expressions have become sufficiently messy,   no-one wants to see them anymore.   They belong on a computer system like SymbolLab or Mathematica... or whatever your preferred mathematical manipulation software might be.   They do not belong on a text document for people to read or differentiate by hand.

    Minor Notes:    We were really maximising  |δU|  ,  the size of a change.  As it happens  δU  will oppose  δT in sign,   which means that we are often actually trying to minimise  δU  because  its negative and a lower negative has a greater size.   I started modifying everything to |δU|  in the early part of the post... but time limits cut in... I never carried it through to the end of the post.
    Here are some plots of what we get from [Eqn 15].

Parameters   N=1  (arbitrary scalar anyway).         ΔH = 10 000   J/mol       ΔS = 33.33  J/mol/K    (deliberately chosen so that  ΔH/ΔS ≈ 300 kelvin ≈ approx. room temperature  and ΔG would become 0 at that temp).

* Plot-300.JPG (66.42 kB . 1363x418 - viewed 2329 times)

   Note that although I would have expected the minima   (maxima of |dU/dT| ) to be  around  T = 300 k,  it was actually a bit below that.


Paremeters as before,   ΔH raised to 20 000 J/mol.    ΔS = 33.33.     Expected  T = 600 kelvin but the maxima is still slightly off.

* plot-600.JPG (82.3 kB . 1401x535 - viewed 2340 times)


Anyway, the general  behaviour of the chemical system seems to be reasonable.....   it always has maximum temperature buffering ability at a particular value of temperature.   Changing the values of ΔH  and  ΔS   move that  maxima   up or down.     You can probably also see why you can't get clear results by taking something like [Eqn 15] and directly differentiating w.r.t.  reaction parameters like ΔH.....      You can have a system that is actually tuned to buffer around a temperature of  100 kelvin *   instead of  300 kelvin   BUT   if  the value of N (the number of molecules in the system) is large enough and also  ΔH   is large eneough then it would still be a numerically better buffer at room temperature than a smaller and more sensibly ΔH valued system.    * LATE EDITING in the previous sentence:  It's more noticeable where the system is tuned to buffer below room temperature not above.  See the actual plots with given parameters.  It is, in my opinion, only by analysing the situation with ΔH and ΔS as fixed parameters and recognising that T is the only thing that can chnage in operation of the device where the natural behaviour of the device is observed.   For a realistic situation, where you can't have an infintely sized system or infinitely high values of ΔH,  you would want to work with the natural behaviour of the system and set it up with ΔH and ΔS  parameters that put it's maximum buffering ability right at the room temperature you wish to try and maintain.


Best Wishes.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: chiralSPO on 01/06/2022 12:30:01
thank you! it will take some time to read, but it looks like you made more headway than I did, so this should be very informative!
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 01/06/2022 12:40:47
You may find this useful.
https://en.wikipedia.org/wiki/Trouton%27s_rule
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 04/06/2022 12:09:23
Hi.

   Come on then @chiralSPO  let's have an update on how you're doing.

[Z]T/[A]T = e–RT/(ΔΗ–TΔS)
    Presumably you noticed that you had the exponent upside down.     e-(ΔH-TΔS)/RT    should have been there.

   However, overall the approach using  the equilibrium constant, K, is what we both went for.

    I think I kept mentioning that ΔH should be as large as possible  BUT   it's actually a bit indirect....  it's the size of ΔS that makes the more obvious difference  however you need  ΔH  to be increasing with it just to keep the maximum response of the chemical equilibrium in the right sort of position (at the right temperature).

    Anyway... by now I was thinking you might have looked at the last formula,  I think it was labelled [Eqn 15] and decided if I made some errors along the way to getting there.   If it's alright, you can continue with it, that's what I would have done.   That's a thing you can examine directly and optimise it by varying the ΔH and ΔS parameters of the reaction.

   Anyway, I don't think you need your two temperatures  T1   and T2    just a single figure  T  where you hope to keep the temperature of the room.

   You can develop the idea and see what happens with more complicated reactions involving multiple particle species and having an equilibrium constant involving products of concentrations raised to some powers.  I'm not saying I want to do it.... the algebra is only going to get more messy... just that it could be done.

Best Wishes.
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Bored chemist on 04/06/2022 16:04:02
  I think I kept mentioning that ΔH should be as large as possible  BUT   it's actually a bit indirect....  it's the size of ΔS that makes the more obvious difference 
You are talking as if they are independent properties, but once you decide on the temperature you want, there's not much choice.
At equilibrium delta G is zero,
so T Delta S = delta H

What you want to do is store heat and that's delta H so you want that to be big (for a given mass, volume, price or whatever).
Title: Re: How can I find the optimum ΔH and ΔS for passive T control?
Post by: Eternal Student on 04/06/2022 18:02:01
Hi.

You are talking as if they are independent properties, but once you decide on the temperature you want, there's not much choice.
    I completely agree.   That was the essence of it. 

    There is one, if quite silly, choice you can make that does allow you to have ΔS and ΔH values chosen independently.    You can develop the chemical equilibrium system to have ΔG of the reaction = 0 at some other temperature (not the desired room temperature) and just accept the consequences.   You'd have the equilibrium constant k >> 1   or  else  k <<1   but never the less,  it would be some equilibrium with  0 < k < ∞.   By Le Chatelier's principle alone if you want to keep it simple, the system will oppose a change in temperature.  It can't perform as well as when the system had k≈1 but it will function to some extent.   So you can compensate by just making the system massive,  i.e. have many molecules of the stuff in the reaction.   
    To phrase this another way, a massive but badly set-up equilibrium system can still oppose temperature changes and help to hold a room at a given temperature.  If it is massive enough then it will be better than a smaller system set up to have ΔG ≈ 0 right at approx. room temperature.
    Other than cost of production, the massive but badly set up system may have other undesirable properties:   It may not help the person in the room to establish the preferred room temperature as much as the finely tuned and well set up system would.   It just maintains it once you have got the room to the preferred temperature initially.   Specifically, the system becomes less responsive on one side of the preferred temperature than on the other side of that temperature.  Let's say it was set up with ΔH and ΔS values to be an optimal temperature buffer at 0 deg. C then it's effectiveness falls off as you move further away from 0 deg C.   So that the person can push the room temperature above the preferred temperature by accident just by leaving their heating system on for an extra half-hour.  Meanwhile, if the equilibrium system is set up to have maximum temperature bufferring ability precisely at 22 deg C then it's ability to absorb / release additional heat remains more consistent on both sides of the preferred temperature.  Therefore it's more forgiving of a careless person who left the heating on for an extra half hour by mistake - they won't raise the room temperature that much higher than the preferred value.

Best Wishes.