Can quantum mechanics predict the likelihood of chemical reactions in mixtures?

[Big-Daddy]

How do I write the non-approximated Schrodinger equation Hamiltonian for a mixture containing 25% by partial pressure of H₂ gas and 75% by partial pressure of He gas, at 100 KPa pressure and 298 K?

 

[Jano L.]

If you want exact Schroedinger equation, you have to take all charged point-like particles - nuclei and electrons - as basic constituents. The number of molecules is irrelevant.


In order to write a meaningful Schroedinger equation, you have to have the particles enclosed in some finite volume. You can model this situation by a potential well described by some function ##V(x,y,z)##. You can model the molecules as a system of charged particles with total potential energy ##U## given by the Coulomb electrostatic energy. The Hamiltonian is

$$
\hat H = \sum_a \frac{p_a^2}{2m_a} + U + \sum_a q_a V(\mathbf r_a).
$$

The index ##a## runs over all nuclei and electrons. The masses and charges are arbitrary and you can choose them such that they correspond to your problem.

Then the Schroedinger equation for the function ##\psi(\mathbf r_1, \mathbf r_2, ...)## is

$$
\partial_t \psi = \frac{1}{i\hbar} \hat{H} \psi.
$$

 

[Big-Daddy]

What does [itex]p_a[/itex] represent, and don't we have to specify the function [itex]V(\mathbf r_a)[/itex]?

Shouldn't there be temperature and pressure dependence? Shouldn't there be dependence on how much of each gas is present in the mixture?

 

[Jorriss]

What does [itex]p_a[/itex] represent, and don't we have to specify the function [itex]V(\mathbf r_a)[/itex]?

[itex]p_a[/itex] is the momentum of particle a.

Shouldn't there be temperature and pressure dependence? Shouldn't there be dependence on how much of each gas is present in the mixture?

Temperature and pressure are statistical mechanical quantities. I suspect, though one may correct me on this, that pressure can still be reasonably interpreted as the negative derivative of the energy with respect to volume?

"How much" of each gas is present will correspond to different solutions. He is more stable than hydrogen and will correspond to lower energy solutions. I don't really know if that is the right way to think about it though - saying the system is 3/4ths He and 1/4th H2 feels like a classical way of thinking about it.

 

[Jano L.]

Shouldn't there be temperature and pressure dependence? Shouldn't there be dependence on how much of each gas is present in the mixture?

These are the quantities that restrict the set of applicable wave functions for your situation. But these restrictions can be applied only after the possible wave functions, solving the Schroedinger equation, are found. The equation itself remains unchanged.

 

[Big-Daddy]

I see. So once all of the solutions are found for the equation (each solution corresponding to one wave-function), how do I then apply a temperature/pressure condition or a mixture composition condition (in this case, 25% H2 and 75% He) to find which wave-function will be representative of the sample under those conditions?

I'm guessing all wave-functions which remain when the set of all solutions is restricted firstly by the temperature/pressure condition and secondly by the mixture composition condition will be the possible isomers under those conditions. And of these, whichever isomers correspond to a lower energy wave-functions will be more prevalent than others under those conditions.

 

[Big-Daddy]

And can one go directly from a molecular wave-function to the molecular orbital set for that molecule?

 

[Jorriss]

And can one go directly from a molecular wave-function to the molecular orbital set for that molecule?

There's no reason to get molecular orbitals if you have the actual wavefunction. MO's are an approximation on the order of hartree fock theory.

 

[Big-Daddy]

Ok then. And let's say I want to know if the H₂ in this mixture (i.e. mixed in with some He) will react with O₂ under certain statistical mechanical conditions (i.e. temperature and pressure). How in principle do we write the equations to figure this out? Solving them as you've said is an entirely different matter.

Can we specify a product and get the rate of the reaction we've written? Or, can we get the expected and likeliest product out of the two reactants, if they will indeed react under the given conditions?

 

[Jorriss]

So you're moving out from quantum mechanics then. If you want to know IF they'll react, regardless of the rate, you'd want to compare the free energies. The free energy difference can give the relative proportion of reactants and products.

Rates are a lot more complicated. If you're comfortable starting to discuss approximations I can tell you a few to look into.

 

[Big-Daddy]

So you're moving out from quantum mechanics then. If you want to know IF they'll react, regardless of the rate, you'd want to compare the free energies. The free energy difference can give the relative proportion of reactants and products.

Ah no. I'm well aware of thermodynamic/kinetic methods. I was just wondering if quantum mechanics will tell us the answers. And using free energies answers only the question: will these reactants react to produce these products? It does not answer the more general question: will these reactants react at all, and if so, what are the likeliest products for their reaction.