Because ,
Using and the given photon density in the deuterium equilibrium abundance formula gives
Although the exponential favors deuterium once , the very small factor expresses the enormous excess of photons. The high-energy tail therefore photodissociates newly formed deuterium until the temperature is well below ; this is the deuterium bottleneck.
A larger increases the deuterium fraction at fixed temperature, so the bottleneck ends at a higher temperature and earlier time. With weak decoupling unchanged, fewer free neutrons then beta-decay before nuclear reactions bind them. Since primordial helium production is neutron-limited, the baryon-density effect on primordial helium makes larger than the standard value.
The number of momentum states in a shell is
With this gives
At zero temperature all states up to are occupied, hence
The energy is . Since a nonrelativistic ideal gas in dimensions obeys ,
Since ,
Using and gives
This has the Sackur-Tetrode form with the effective volume , but its additive constant is rather than . The extra unit arises because the atmosphere's vertical extent grows with temperature and its mean gravitational energy per atom is .
Increasing makes smaller at every fixed temperature, so freeze-out occurs earlier, at a higher . The equilibrium value is then closer to one. Although some neutrons subsequently decay, the final neutron fraction is larger, and
increases, approaching unity if the surviving neutron and proton abundances become comparable.
Such a universe begins with less hydrogen and more helium. It has less fuel for long-lived hydrogen-burning stars; stellar evolution would tend to proceed through hotter, faster helium-burning stages, reducing typical stellar lifetimes. The shorter stable-energy window, altered heavy-element production, and lower primordial hydrogen abundance would make familiar water-rich, slowly evolving habitats less likely, although a detailed conclusion would require a stellar-evolution model.
At equal pressure and temperature, ideal gases have equal molecular number density. The mean air molecular mass is , while helium has mass , so
For , remains stationary. At large , evaluating the free energy near a positive makes the negative cubic energy dominate the bounded entropy cost, so some has lower free energy than zero. At high temperature entropy uniquely favors zero. Continuity in temperature therefore forces a transition at an intermediate .
Conservation of thermal energy and Fourier's law give the heat equations
The density times the specific heat capacity, , is the volumetric heat capacity, and is the thermal conductivity.
The condition holds the left exterior face at a fixed zero reference temperature; does the same for the right exterior face. The condition expresses perfect thermal contact, with no temperature jump at the interface. Finally,
expresses continuity of the signed heat flux in the positive direction: the interface neither stores nor produces heat.
For , substitution of the separated solutions into the heat equations gives
The zero exterior boundary values are incorporated by
To obtain a eigenvalue equation for heat conduction in two joined slabs, write and . Continuity of temperature and heat flux at zero gives the homogeneous linear system
A nonzero amplitude pair exists precisely when its determinant vanishes, equivalently
If , division by gives the printed transcendental equation
The printed necessity claim needs the qualification that both tangents are defined. If just one of vanishes, . If both cosines vanish, however, and legitimate normal modes remain. Their temperature at the thermal interface need not be zero; the temperature condition fixes the amplitude ratio, while each heat flux is zero there.
For a concrete counterexample, take identical slabs: , , , and . Then
with each formula restricted to its own slab, satisfies both heat equations and all four boundary conditions. But , so neither tangent is defined. This directly disproves the unqualified printed assertion. The determinant condition retains these missing normal modes.
Finally, must be checked before dividing by it. Then , , and interface continuity gives , while heat flux continuity gives . Positivity of the thermal conductivities forces . Thus In the determinant formula, its apparent zero at is spurious because the chosen sine basis degenerates there.
Choose the positive weights
These are the volumetric heat capacities. The heat equations and integration by parts give
The exterior boundary terms vanish because the temperatures there are zero. At the thermal interface, the remaining term is
by continuity of temperature and heat flux. Hence the weighted energy estimate is
Both weights are positive, so . Zero initial temperatures give , forcing for all . Each temperature therefore vanishes throughout its slab. Applying this argument to a difference also proves uniqueness for arbitrary common initial data in this thermal interface problem.
Assume the rod has constant cross-sectional area , no internal source of heat, no bulk motion, and no lateral heat flux. Assume also that each section has a uniform temperature, so that depends only on , and that Fourier's law applies with the stated constant thermal conductivity . The longitudinal heat flux density is then .
By conservation of energy, the rate of increase of internal energy in a slice equals incoming minus outgoing heat:
Divide by and take the limit . For a sufficiently differentiable temperature field this gives , and hence
Here is mass density, is specific heat capacity, and is thermal diffusivity.
A steady state has , so . The fixed temperature gives , and thermal insulation gives . Thus
Subtracting this steady state produces the homogeneous mixed boundary conditions
The partial differential equation and boundary conditions apply for . The initial condition is attained as . Its slope at is , so it does not satisfy the insulated-end condition at the instant . The resulting heat equation solution is smooth for and approaches the initial temperature continuously; its spatial derivative need not be continuous at .
For separation of variables, put and set . The boundary conditions are , . Multiplying by and using integration by parts gives
A nonzero solution therefore has : the zero case forces constant, and then forces . Writing , the eigenfunctions and their temporal factors are
These sine functions are orthogonal on , with squared norm . Expanding the initial condition and multiplying by gives, using the permitted interchange of the sum and integral,
Here and . Thus the Fourier series solution is
The exponential decay permits differentiation of the series for , and each term satisfies both boundary conditions. At it is the mixed-end Fourier expansion of .
The gained heat is the increase in internal energy. Since the paper specifies no cross-sectional area, first give the answer per unit cross-sectional area; multiply by for a rod of cross-sectional area . The initial temperature has integral , whereas . At this gives
Numerically . This is positive and less than the eventual gain , as expected from relaxation towards the uniform steady state.
The specific heat capacity is the energy needed per unit mass per unit increase in temperature. With constant, a small element of volume has mass and thermal energy , measured relative to the reference at which . Summing over gives
Choosing a different reference temperature only adds a time-independent constant to this energy, so it does not affect the heat balance with volumetric heat generation.
Suppose and are two sufficiently regular steady state solutions with the same prescribed temperature on the boundary. Their difference satisfies
The identical PDE source terms cancel. Apply the product rule and the divergence theorem to ; this is the weighted form of Green's first identity:
Since , the nonnegative integrand must vanish, so . Hence is constant on each connected component of , and its zero Dirichlet boundary condition makes each constant zero. Therefore : the solution, if it exists, is unique. This proves Dirichlet uniqueness for variable-conductivity heat conduction; it does not require a constant thermal conductivity.
Let be the constant cross-sectional area, and suppose that the temperature is uniform on each cross-section. Consider a slice . Its stored thermal energy is . Since the lateral surface is insulated and there is no internal heat generation, conservation of energy gives
where is the heat flux in the positive direction. By Fourier's law, , so
The slice is arbitrary; therefore the local heat equation is
The constant is the thermal diffusivity. The cancellation of uses a rod of uniform cross-section; the one-dimensional model assumes no transverse temperature variation.
Subtract the constant steady state profile: set . Then
The last equation follows from Fourier's law and zero heat flux at the insulated end. For a nonzero separated solution , separation of variables gives
To see that for a nonzero eigenfunction, multiply the spatial equation by and use integration by parts:
Both boundary terms vanish. If the right side were zero, would be constant and would make it zero. Thus write with . Solving gives , and requires . The mixed Dirichlet-Neumann modes on an interval therefore have
The eigenvalues of are ; are their positive square roots. By superposition, with coefficients for which the resulting series and its required derivatives converge at positive times,
The eigenfunctions are orthogonal on and have squared norm . Expanding the initial excess temperature gives
Thus the requested temperature is
At the insulated endpoint cancels the alternating coefficient sign, leaving
Every summand is positive. The endpoint temperature never drops below ; it decreases towards as and starts at .
For completeness, the Sturm-Liouville eigenfunction expansion gives the initial profile in the sense. Here the coefficients are absolutely summable, so the Weierstrass M-test also gives uniform convergence on , and the initial values are recovered everywhere. For positive times the exponential factors justify all required termwise derivatives, giving a classical solution. There is, however, an initial incompatibility with a Neumann boundary condition at : the initial slope is , while the insulated Neumann boundary condition requires zero slope for . The temperature is continuous to the initial edge, but cannot extend continuously to that corner. One must not differentiate the initial series term by term at .
If for a twice continuously differentiable potential , then and . Equality of mixed partial derivatives gives the necessary condition for an exact differential:
This condition alone does not ensure global exactness on a domain with holes.
For the Helmholtz free energy , use the stated identity for internal energy and the product rule:
Thus and . Equality of mixed partial derivatives yields the Maxwell relation
The temperature and volume are the independent variables here, with composition and amount of gas fixed.
Next express and substitute its total differential into :
Matching the coefficient of gives
This is the thermodynamic internal pressure identity: an equation of state determines how internal energy changes under an isothermal change of volume.
For the Dieterici equation, set while differentiating at fixed . Then
Substitute into the thermodynamic internal pressure identity. The term proportional to cancels , leaving
The physical gas domain has and .
For the ideal gas instead, gives . Therefore is independent of volume at fixed amount of gas. Integrating the remaining condition , with constant heat capacity at constant volume, gives
where is independent of both and and fixes the reference for internal energy.
Insert the linear initial profile into the transformed integral from (b):
The constant part integrates to . The other part is an absolutely integrable odd function, so its integral vanishes. Hence
As a direct check, and . Thus a linear temperature profile is stationary on the whole line. More generally, every affine function is preserved by this heat-kernel solution.
Along a change in volume at fixed temperature, the first law of thermodynamics gives
Use the Helmholtz Maxwell relation derived in (a) to replace the entropy derivative. This proves
For the stated Van der Waals equation, temperature-independent give . Hence . Integrating with respect to yields
The dilute-limit condition holds at every fixed and forces , so
The result is the internal energy of a van der Waals gas in the fixed-amount convention of the question. Its volume derivative is positive even though the attractive energy correction itself is negative.