Three approaches for determining the thermodynamic stability of irreversible processes are described in generalized formulations. The simplest is the Gibbs-Duhem theory, specialized to irreversible trajectories, which uses the concept of virtual displacement in the reverse direction. Its only drawback is that even a trajectory leading to an explosion is identified as a thermodynamically stable motion. In the second approach, we use a thermodynamic Lyapunov function and its time rate from the Lyapunov thermodynamic stability theory (LTS, previously known as CTTSIP). In doing so, we demonstrate that the second differential of entropy, a frequently used Lyapunov function, is useful only for investigating the stability of equilibrium states. Nonequilibrium steady states do not qualify. Without using explicit perturbation coordinates, we further identify asymptotic thermodynamic stability and thermodynamic stability under constantly acting disturbances of unperturbed trajectories as well as of nonequilibrium steady states. The third approach is also based on the Lyapunov function from LTS, but here we additionally use the rates of perturbation coordinates, based on the Gibbs relations and without using their explicit expressions, to identify not only asymptotic thermodynamic stability but also thermodynamic stability under constantly acting disturbances. Only those trajectories leading to an infinite rate of entropy production (unstable states) are excluded from this conclusion. Finally, we use these findings to formulate the Fourth Law of thermodynamics based on the thermodynamic stability. It is a comprehensive statement covering all nonequilibrium trajectories, close to as well as far from equilibrium. Unlike previous suggested \"fourth laws\", this one meets the same level of generality that is associated with the original zeroth to third laws. The above is illustrated using the Schlögl reaction with its multiple steady states in certain regions of operation.

To explore abiotic theories related to the origin of biomolecular homochirality, we analyze two entirely reversible kinetic models composed of an enantioselective autocatalysis with limited stereoselectivity that is coupled to an enantiomeric mutual inhibition (Frank-like models). The two models differ in their autocatalytic steps in respect to the formation of monomer species in one model and of dimer species in the other. While fully reversible and running in a closed system, spontaneous mirror symmetry breaking (SMSB) gives rise to transient chiral excursions, even when starting from a strictly achiral situation. Before the SMSB, the two models differ in the main dissipative processes. At the SMSB, the entropy production rate reaches its maximum in both models. Here it is the enantioselective autocatalysis with retention of the winner enantiomer that dominates. During the terminal phase, the enantioselective autocatalysis with inversion prevails, while the entropy production rate vanishes, thus fulfilling the conditions of microscopic reversibility. SMSB does not occur if the autocatalytic rate constant is too strong or too weak. However, when the autocatalysis is relatively weak, the temporary chiral excursions last for long periods of time and could be the starting point of a cascade of asymmetric reactions. The realism of such Frank-like models is discussed from the viewpoint of their relevance to prebiotic chemistry.