Next Article in Journal
The Investigation of Organic Binder Effect on Morphological Structure of Ceramic Membrane Support
Next Article in Special Issue
Mirror Symmetry Breaking in Liquids and Their Impact on the Development of Homochirality in Abiogenesis: Emerging Proto-RNA as Source of Biochirality?
Previous Article in Journal
Modified Evolutionary Algorithm and Chaotic Search for Bilevel Programming Problems
Previous Article in Special Issue
Chiral Interface of Amyloid Beta (Aβ): Relevance to Protein Aging, Aggregation and Neurodegeneration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spontaneous Chiral Symmetry Breaking and Entropy Production in a Closed System

1
Department of Chemistry, Wake Forest University, Winston-Salem, NC 27109, USA
2
BASF, 100 Park Avenue, Florham Park, NJ 07932, USA
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(5), 769; https://doi.org/10.3390/sym12050769
Submission received: 2 March 2020 / Revised: 13 April 2020 / Accepted: 14 April 2020 / Published: 6 May 2020
(This article belongs to the Special Issue Asymmetry in Biological Homochirality)

Abstract

:
In this short article, we present a study of theoretical model of a photochemically driven, closed chemical system in which spontaneous chiral symmetry breaking occurs. By making all the steps in the reaction elementary reaction steps, we obtained the rate of entropy production in the system and studied its behavior below and above the transition point. Our results show that the transition is similar to a second-order phase transition with rate of entropy production taking the place of entropy and the radiation intensity taking the place of the critical parameter: the steady-state entropy production, when plotted against the incident radiation intensity, has a change in its slope at the critical point. Above the critical intensity, the slope decreases, showing that asymmetric states have lower entropy than the symmetric state.

1. Introduction

Modern thermodynamics, formulated in the 20th century by Onsager [1], De Donder [2], Prigogine [3,4], and others, introduced a critical concept lacking in its classical formulation: rate of entropy change and its relationship to irreversible processes. Classical thermodynamics was concerned with functions of state, such as energy and entropy, and their change from one equilibrium state to another. Absent from this theory of states is consideration of the rates of processes. Changes in entropy for infinitely slow reversible processes are calculated using the relation dS = dQ/T, (in which dS is the change in entropy, T is the temperature in Kelvin, and dQ is the heat exchanged between a system and its exterior). However, for changes that take place in a finite time due to irreversible processes, the theory does not specify a way of calculating the entropy change; it is only stated that dS > dQ/T. Modern thermodynamics is a theory of processes in which thermodynamic forces and the flows they drive are identified and the rate of entropy production is expressed in terms of these thermodynamic forces and flows [5,6]. More specifically, the rate of entropy production per unit volume, σ, is expressed in terms of the forces and flows as
σ = d s d t = k F k J k
in which s is the entropy density, Fk are the thermodynamic forces and Jk are thermodynamic flows. A temperature gradient, for example, is the thermodynamic force, Fk, that drives thermodynamic flow, Jk, of heat current. The force that drives chemical reactions has been identified as affinity [2,5,6] and the corresponding flow is the rate of conversion form reactants to products. This flow is expressed as the time derivative dξ/dt (mol/s) of the extent of reaction ξ [5,6]. For an elementary chemical reaction step, the rate of entropy production can be written in terms of the forward reaction rate, Rf, and the reverse reaction rate, Rr [6]
σ = R ( R f R r ) ln ( R f R r )
in which R is the gas constant. The total entropy production for a sequence of reactions is the sum of entropy productions of each reaction (indexed by k) [6]
σ = R k ( R k f R k r ) ln ( R k f R k r )
We note that calculations of entropy production using the above formulas require that the reverse reaction rates, Rr, are non-zero. In considering kinetic equations of a chemical system, often the reverse rates have high reaction barriers and, correspondingly, very low reaction rates, and are assumed to be “zero” because they are negligible compared to the forward rates. Low reverse reaction rates keep the system from evolving to equilibrium state. In the model we will present below, the system is driven far from equilibrium by an inflow of radiation. In the absence of radiation, our system evolves to equilibrium—there are no very high reaction barriers to keep it from reaching its equilibrium state.
It is well known that nonequilibrium chemical systems can undergo spontaneous symmetry-breaking transitions to organized structures called dissipative structures [6,7]. The large class of dissipative structures includes spatial patterns, chemical clocks, and structures with chiral asymmetry. These structures arise when a nonequilibrium system becomes unstable and undergoes a transition to new organized state. The model we present below shows a chiral symmetry breaking transition to an asymmetric state as the intensity of radiation increases.
The study of spontaneous chiral symmetry breaking in chemical systems has a long history, starting in the 1950s with the first model by Frank [8]. In the 1970s, when the theory of dissipative structures was developed, instability, bifurcation, and spontaneous symmetry breaking in non-equilibrium systems became the foundation for the study of chiral asymmetry we see in nature [4,7,9,10,11]. A general theory of chiral symmetry breaking in chemical systems, independent of the details of the chemical kinetics, was formulated, and it was used to study the sensitivity of chiral symmetry breaking systems to small chiral influences [11,12,13]. Using this theory, it was possible to calculate the time scales needed for a chiral-symmetry-breaking chemical system to be influenced by the chiral asymmetry of electro-weak interaction in molecules [14], and it was found that this timescale is of the order to 104 years [15]. These results show that a lot of interesting and important general conclusions can be arrived at through using theoretical models [16,17]. Along these lines, we investigate the thermodynamic aspects of systems that spontaneously break chiral symmetry, using a model presented below.
Using nonequilibrium thermodynamics, we analyze the behavior of entropy production, σ, for a reaction scheme that consists of a photochemically driven closed system (that has no flow of matter). The incident radiation drives a generation–decomposition cycle of chiral molecules. At the critical intensity (above which the system breaks chiral symmetry), the slope of the rate of entropy production, σ, changes just as entropy does in a second-order phase transition. We have reported a similar result in a flow system with an inflow of reactants and outflow of products [18]. There is a basic difference, however: in our previous study, the slope of σ increased, but in our current study we find that the slope decreased, although in both cases the behavior of σ is similar to that in a second-order phase transition.

2. Materials and Methods

All reactions in our scheme are reversible so that the system can reach a state of chemical equilibrium. At equilibrium, the rate of entropy production is zero. The reaction scheme of this photochemically driven closed system is shown in Table 1. It is assumed to take place in a homogeneous aqueous phase in which radiation is incident. For a photochemical reaction, the intensity of a narrow band of wavelength is the relevant intensity, therefore it will have a low value compared to a typical intensity of blackbody radiation that includes all wavelengths—for example, radiation from the sun. An achiral molecule, T, is photochemically activated to an excited species or a more reactive isomer, TE, as shown in reaction (R1); II is the intensity of radiation of the exciting wavelength. TE can radiate its excitation energy and return to the unexcited state, T (R1a), through various processes. During the transition of TE to T, the emitted radiation undergoes scattering and thermalizes to the temperature of the system. (R1a) represents all the processes that keep TE and T at equilibrium (including thermal radiation [6]) in the absence of external radiation. In writing the reaction rates for T and TE, we only need to include additional term in the forward reaction that includes II and the all reverse reactions, TE to T, combined into one reaction rate. The excited species TE can also react with an achiral molecule S to form a chiral species, XL or XD (R2 and R3). The set of reactions R4a and R4b are elementary steps of an autocatalytic reaction for XL with the intermediate species SL; similar reaction steps result in the autocatalytic production of XD, as shown in R5a and R5b. Reaction steps R6, R7 and R8 show reactions through which the species decompose into starting material S and T. The scheme explicitly has all the steps needed for the system to reach chemical equilibrium in the absence of radiation.
The above model is a variation of the models in our earlier studies [11,12,15] which are modifications of the original model of Frank [8]. The modifications allow us to analyze non-equilibrium symmetry breaking and rate of entropy production. Models such as this are used to extract general properties that are not model dependent. Examples of such properties are the qualitative behavior of steady-state rate of entropy production as a function of a parameter that drives the system away from equilibrium (such as the incident radiation intensity II in the above model). The difference in concentration between enantiomers of a chiral species as a function of a parameter, such as the intensity II, is generally parabolic, as predicted by bifurcation theory based on the symmetry group (mirror symmetry in this case) of the system, regardless of the details of the chemical reactions that break chiral symmetry.
Though there is currently no known reaction that has all the properties in the above model, the reaction has no steps that are implausible. For example, reactions (R1), (R1a), (R2) and (R3) comprise a photoaddition reaction that produces a chiral compound. An example is the following reaction series [19,20]:
(Ph)2 C = CH(R1) + hv → [(Ph)2 C = CH(R1)]* (R9)
[(Ph)2 C = CH(R1)]* + (R2)OH → [(Ph)2 HC − CH(R1)(OR2)]L(R10)
[(Ph)2 C = CH(R1)]* + (R2)OH → [(Ph)2 HC − CH(R1)(OR2)]D (R11)
in which Ph is the phenyl group and R1 = CH3 or C2H5 and R2 = CH3, C2H5 or C3H7. In the reaction (R9), a photon is absorbed by the electrons in the C=C double bond and the molecule transitions to a reactive excited state [(Ph)2 C = CH(R1)]*. In the addition reaction shown in (R10) and (R11), the excited molecule reacts with an alcohol, (R2) OH, and produces a chiral compound (Ph)2 HC − CH(R1)(OR2) in the L and D enantiomeric forms. In this compound, the carbon shown in boldface is a chiral carbon (its tetrahedral bonds to four different groups makes it so). Other examples of photoaddition reactions producing chiral products from achiral reactants can be found in [19]. We note that TE need not be an excited state; it could be a different, more reactive isomer of the T [19].
The reactions (R4a)–(R5b) are steps leading to chiral autocatalysis. This involves the formation of a chiral complex of S and X in their enantiomeric forms. Examples of chiral complexes resulting in reactions with a high degree of chiral selectivity have been known for a long time [21,22]. In an article published in 1984, we noted some mechanisms that are based on chiral ligands in rhodium phosphine catalysts which could lead to chiral autocatalysis [12,22]. To date, there are several chirally autocatalytic reactions that have been experimentally studied. Chiral symmetry breaking was noticed and systematically studied first in NaClO3 in 1990 [23], and in 1995 chiral autocatalysis and amplification of small initial enantiomeric excess was reported in inorganic reactions involving cobalt complexes [24] and in organic reactions involving alkylation of aldehydes [25]. Since then, these and closely related systems have been extensively experimentally studied and the mechanisms of chiral autocatalysis have been investigated [26,27,28,29,30,31]. A variant of chiral symmetry breaking in stirred crystallization was reported in 2005, and it too has been studied extensively [32,33]. The mechanisms of chiral autocatalysis vary in these systems: in crystallization, it is secondary nucleation, in the organic and inorganic reactions, cluster/complex formation seems to be involved. Reactions (R4a) and (R5a) may be thought of as a simple form of chiral complex formation. As was noted in a review [27], the exact details of chiral catalysis are not of significance for the general properties symmetry-breaking bifurcation and thermodynamic properties of such systems. In particular, properties such as phase-transitions-like behavior we present here are quite independent of the details of chemical kinetics. Examples of reaction (R6), the dimer formation of enantiomers, are also known; in fact, such dimerization of certain chiral catalysts leads to asymmetric amplification [34].
For the above theoretical model (R1)–(R8), the corresponding forward and reverse rate for each reaction are written as follows, in which concentrations are shown explicitly as functions of time:
R1f = (k1f + k1 II)T[t], R1r = k1r TE[t]
R2f = k2f S[t] TE[t], R2r = k2r XL[t]
R3f = k3f S[t] TE[t], R3r = k3r XD[t]
R4af = k4af S[t] XL[t], R4ar = k4ar SL[t]
R4bf = k4bf SL[t] TE[t], R4br = k4br (XL[t])2
R5af = k5af S[t] XD[t], R5ar = k5ar SD[t]
R5bf = k5bf SD[t] TE[t], R5br = k5br (XL[t])2
R6f = k6f XL[t]XD[t], R6r = k6r P[t]W[t]
R7f = k7f P[t], R7r = k7r (S[t])2
R8f = k8f W[t], R8r = k8r (T[t])2
In these equations, the rate constants are written as k1f, k1r etc., and the concentration are written as T[t], S[t], etc. In terms of these forward and reverse rates, the rate equations for the concentrations can be written as:
d T[t]/dt = −R1f + R1r + 2R8f − 2R8r
d TE[t]/dt = −R1f + R1r − R2f + R2r − R3f + R3r − R4bf + R4br − R5bf + R5br
d S[t]/dt = −R2f + R2r − R3f + R3r –R4af + R4ar − R5af + R5ar + 2R7f − 2R7r
d SL[t]/dt = R4af − R4ar − R4bf + R4br
d SD[t]/dt = R5af − R5ar − R5bf + R5br
d XL[t]/dt = R2f − R2r − R4af + R4ar + 2R4bf − 2R4br − R6f +
d XD[t]/dt = R3f − R3r − R5af + R5ar + 2R5bf − 2R5br − R6f + R6r
d P[t]/dt = R6f − R6r − R7f + R7r
d W[t]/dt = R6f − R6r − R8f + R8r
This set of coupled non-linear equations were solved numerically using Mathematica NDSolve. NDSolve is a Mathematica command that has the following structure: NDSolve[{Equations}, {yi},{t, tmin, tmax}], in which “Equations” are the differentials equations for the set of functions {yi} with t as the independent variable; numerical solutions are obtained in the range tmin, to tmax. More details can be found in the online documentation that comes with Mathematica. The rate constants used for the numerical solutions are summarized in Figure 1. In assigning values to rate constants, there are consistency conditions that must be met. For example, since reactions (R4a) and (R4b) together are equivalent to (R2), the products of the equilibrium constants of R4a and R4b must equal the equilibrium constant of R2. This gives us the following condition for the rate constants:
(k4af/k4ar)(k4bf/k4br) = k2f/k2r
Numerical values were assigned to rate constants so as to fulfill these requirements. The units were chosen so that all concentrations are in mol/L. Assigned numerical values are such that the concentrations of the reactants have realistic values. Symmetric (racemic) and asymmetric states are parametrized by α = (XL − XD), in the symmetric state α = 0 and in the asymmetric state α ≠ 0.

3. Results

The rate equations were first solved for an equilibrium state where the incident radiation intensity, II, was set to 0. The initial concentrations of species S and T were set to 0.01 M and the initial concentrations of all other species were set to 0.0 M. Under these conditions, the system evolves to its racemic equilibrium state, in which α = 0.
The simulation code was run for sufficient time (about 104 s) to ensure the concentrations of all species have reached a steady state, which is the equilibrium state. At t = 104 s, the concentrations at equilibrium were: S = T = 8.478 × 10−3 M, TE = 8.477 × 10−6 M, SL = SD = 3.047 × 10−6 M, XL = XD = 3.593 × 10−5 M, and P = W = 7.188 × 10−4 M. The conversion of the initial species S and T compared to other species was rather small for the numerical values of the rate constants shown in Figure 1. By choosing a different set of rate constants, the conversion could be increased. The numerical values confirm that complete symmetry of the system was maintained when no incident radiation is present. Figure 2 shows the time evolution of the chiral species SL, SD, XL, XD from t = 0 s, to t = 1000 s. As the system evolved to its equilibrium state, the entropy production σ was monitored; initially, it took a nonzero value but, as expected, its value decreased to zero at the equilibrium state.
We then used these values for a symmetric equilibrium state as initial values for a system subject to a radiation input, i.e., II > 0. This radiation input serves as a means to push the system away from thermodynamic equilibrium. A very small excess, about 0.1% (3 × 10−8 M) of XL was introduced into the system as a “random fluctuation”. If the system has the mechanism to break chiral symmetry, it will have a critical value IIC. At values of II < IIC, this excess 0.1% of XL will decrease and the system will again evolve into a steady state where XL = XD; at values of II > IIC, the excess will increase and lead to a steady state in which XL > XD. In a real system, this slight perturbation may be due to a random fluctuation such as a local excess of one enantiomer that may then be amplified, resulting in a state of broken symmetry. The overall behavior of the system is summarized in Figure 3.
It was found that this system indeed has the mechanism needed for breaking chiral symmetry. At values of II < 0.004, the system evolved to a symmetric state corresponding to XL = XD and α = 0. Figure 4 shows the time evolution of the chiral species when II = 0.0030. A steady state is reached in approximately 600 s.
For values of II > 0.004, the small excess of 0.1% of XL in the initial concentration increased, thus driving all chiral species to an asymmetric state. A time evolution of the chiral species in an asymmetric state when II = 0.008 is shown in Figure 5.
As shown, steady state is reached about 7000 s. The time taken to reach steady state, the relaxation time, depends on the value of II and on the amount of initial excess (0.1% of XL). As is well known in the study of stability of steady states [6,7,9], the initial exponential growth of the small enantiomeric excess depends on the eigenvalue of the unstable mode of the linearized equations derived from the set (14)−(22) around the initial state. These linearized equations are obtained by assuming a small perturbation of the concentrations, δCk, from the initial steady state. This leads to a set of linear equations of the type dδCk /dt = Σl LklδCl [6,7,9]. The eigenvalues of Lkl with positive real parts are the exponents that determine the initial rate of growth of the enantiomeric excess. However, the later growth and leveling off at the steady state depend on the kinetics and rate constants. In general, near the critical point IIC, the relaxation time is long, the so-called “critical slowing”, but as the value of II increases, the growth rate becomes faster and the relaxation time decreases.
To determine the relation between the steady-state value of α on II, the reaction was run for various values of II, from 0.0035 to 0.0045, and the corresponding steady-state values of α were obtained. As noted above, the time it takes for the system to reach an asymmetric steady state depends on the value of II; close to the critical value IIC (0.004 in this case), the relaxation to asymmetric steady state is slower and it becomes faster as the value of II increases. The exact quantitative relationship between relaxation time and II depends on the kinetics and rate constants and not of significance to the current study. In these runs, to obtain both positive and negative branches of α, the initial condition with a 0.1% excess of XD was also included. The dependence of α on II is shown in Figure 6, demonstrating the typical bifurcation of asymmetric states above the critical value IIC = 0.004. As is expected, in a chiral symmetry breaking transition, the values of α above the critical point are parabolic.
With these results, we now turn to the rate of entropy production σ. As stated above, initially the system is in the state of equilibrium, with the intensity of radiation II = 0 and σ = 0. Then the intensity of the radiation II is increased to a non-zero value and the rate of entropy production σ is monitored. Initially, σ sharply increases and, as the system reaches its steady-state, corresponding to the set value of II, the entropy production also reaches its steady state value. Figure 7 shows the evolution of σ from its value when t = 0, to its final steady-state value when II = 0.008. The final steady state value of σ is small compared to its initial value, but it is nonzero.
We would like to note that the results shown above do not depend crucially on the particular values assigned to rate constants. Whether symmetry breaking occurs or does not depends mostly on the mechanism of reactions in the model, not on a narrow range of values of rate constants. In general, qualitative properties change drastically due to small changes in rate constants only in singular cases. The presented model is not singular. In our study, we have tried a range of rate constants and observed symmetry breaking. A typical value is presented in this article.
Our objective is to study the behavior of σ as the intensity of radiation II moves form a value below to a value above the critical value IIC. In a previous study [18] of an open system, σ behaved as entropy does in a second order phase transition: its slope is discontinuous at the transition point. In addition, we noted that its slope increased above the critical point. This behavior was consistent with the Maximum Entropy Production (MEP) hypothesis [35,36,37,38,39,40] that states that non-equilibrium steady states maximize the rate of entropy production. In other words, the breaking of symmetry was a reflection of the general tendency of non-equilibrium systems and, from this point of view, it was only to be expected. This would imply that biochemical asymmetry is a consequence of MEP. To date, there is no general proof for MEP; indeed, in our own investigation we found that MEP was valid for some systems but not all. MEP, in general, is a controversial hypothesis and its applications to various complex systems have been questioned [41]. Hence, we investigated MEP in the context of chiral symmetry breaking.
Figure 8 shows the behavior of entropy production in the closed system we present here. It shows a change in the slope at the transition point, IIC, as in our previous study. From the point of view of symmetry breaking transition, we see that entropy production behaves in a way that is similar to that of entropy in a second-order phase transition. However, unlike the previous result, the slope decreases above the critical point. As indicated in Figure 8, when II > IIC, if initially we set the system in a non-equilibrium symmetric state, which is an extrapolation of the symmetric state below IIC, the system evolves to an asymmetric state in which σ is lower. The extrapolation of the symmetric state beyond the critical value IIC is possible on a computer because, without a small perturbation in XL or XD, or other chiral species, the system stays in an unstable symmetric steady state. With a small perturbation, it evolves to the stable asymmetric state. The time evolution of the system from this initial state to its final asymmetric steady state results in the decrease in σ, thus indicating that the stable asymmetric state is associated with a lower value of σ compared with that of a symmetric state. Thus, we see that the entropy production in a closed system is not consistent with MEP, because MEP would predict higher values for σ in the asymmetric state.

4. Concluding Remarks

Our results show several aspects of chiral symmetry breaking. First, they show that Frank’s original model can be modified with several additional steps—all of which are possible elementary chemical reaction steps—to demonstrate that spontaneous chiral symmetry breaking can occur in a photochemically driven closed system. Our model is motivated by the fact that the earth is a essentially a closed system (except for a very small influx in interstellar matter such as meteorites) and the evolution of life was driven by incident solar radiation. The incident radiation drives a cycle of generation and decomposition of chiral compounds that, at a sufficiently high intensity of radiation, makes a transition to a state of broken chiral symmetry. This example demonstrates that life on earth could have evolved under such conditions of prebiotic molecular chiral asymmetry.
From a thermodynamic viewpoint, this study also confirms that chiral symmetry breaking transitions are similar to second-order phase transitions, with entropy production taking the place of entropy. We note that the general qualitative features of bifurcation of asymmetric states and change in the behavior entropy production at the critical point, shown in Figure 6 and Figure 8, are a consequence of the two-fold symmetry of the system, not the particularities of the chemical reactions in the model. In fact, using the system’s symmetry group, it is possible to derive the following generic equation for the evolution of α near the critical point: dα/dt = −Aα3 + Bα + C, in which A, B and C are functions of the kinetic constants of the chiral symmetry breaking chemical reactions [11,12,13,15]. Finally, we show that the behavior of entropy production in this chiral-symmetry-breaking system is not consistent with the MEP hypothesis, because MEP implies that the state of broken symmetry will have a higher rate of entropy production compared to a symmetric state, but we find the opposite to be the case in this model. In a system with an inflow of reactants and outflow of products, however, the entropy production was higher in the asymmetric state. This indicates that MEP is valid for a certain class of systems, but what this class is has not yet been clearly identified.

Author Contributions

Conceptualization, D.K.; methodology, D.K. and Z.M.; software and Mathematica Code, D.K. & Z.M. All authors have read and agreed to the published version of the manuscript.

Funding

Wake Forest University, Winston Salem, NC 27109, USA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Onsager, L. Reciprocal Relations in Irreversible Processes. Phys. Rev. 1931, 37, 405–426. [Google Scholar] [CrossRef]
  2. De Donder, T.; Van Rysselberghe, P. Affinity; Stanford University Press: Menlo Park, CA, USA, 1936. [Google Scholar]
  3. Prigogine, I. Etude Thermodynamique des Phénomènes Irréversibles; Editions Desoer: Liège, Belgium, 1947. [Google Scholar]
  4. Prigogine, I. Introduction to Thermodynamics of Irreversible Processes; John Wiley: London, UK, 1967. [Google Scholar]
  5. De Groot, S.R.; Mazur, P. Non-Equilibrium Thermodynamics; North Holland Press: Amsterdam, The Netherlands, 1969. [Google Scholar]
  6. Kondepudi, D.K.; Prigogine, I. Modern Thermodynamics: From Heat Engines to Dissipative Structures, 2nd ed.; John Wiley: New York, NY, USA, 2015. [Google Scholar]
  7. Nicolis, G.; Prigogine, I. Self-Organization in Non-equilibrium Systems; Wiley: New York, NY, USA, 1977. [Google Scholar]
  8. Frank, F.C. Spontaneous asymmetric synthesis. Biochem. Biophys. Acta 1953, 11, 459–463. [Google Scholar] [CrossRef]
  9. Haken, H. Synergetics: An Introduction; Springer: Berlin, Germany, 1977. [Google Scholar]
  10. Kondepudi, D.K.; Hegstrom, R.A. The Handedness of the Universe. Sci. Am. 1990, 262, 108–115. [Google Scholar]
  11. Kondepudi, D.K.; Nelson, G.W. Chiral Symmetry Breaking in Nonequilibrium Systems. Phys. Rev. Lett. 1983, 50, 1023–1026. [Google Scholar] [CrossRef]
  12. Kondepudi, D.K.; Nelson, G. Chiral Symmetry Breaking states and Their Sensitivity in Nonequilibrium Systems. Physica A 1984, 125, 465–496. [Google Scholar] [CrossRef]
  13. Kondepudi, D.K. State Selection Dynamics in Symmetry Breaking Transitions. In Noise in Nonlinear Dynamical Systems; Moss, F., McClintock, P.V., Eds.; Cambridge University Press: Cambridge, UK, 1989; Volume 2, pp. 251–270. [Google Scholar]
  14. Hegstrom, R.A.; Rein, D.W.; Sanders, P.G.H. Calculation of the parity nonconserving energy difference between mirror-image molecules. J. Chem. Phys. 1980, 73, 2329–2341. [Google Scholar] [CrossRef]
  15. Kondepudi, D.K.; Nelson, G.W. Weak Neutral Currents and the Origin of Biomolecular Chirality. Nature 1985, 314, 438–441. [Google Scholar] [CrossRef]
  16. Decker, P. The origin of molecular asymmetry through the amplification of “stochastic information” (noise) in bioids, open systems which can exist in several steady states. J. Mol. Evol. 1974, 4, 49–65. [Google Scholar] [CrossRef]
  17. Plasson, R.; Kondepudi, D.K.; Bersini, H.; Commeyras, A.; Asakura, K. Emergence of homochirality in far-for-equilibrium systems: Mechanisms and role in prebiotic chemistry. Chirality 2007, 19, 589–600. [Google Scholar] [CrossRef]
  18. Kondepudi, D.K.; Kapcha, L. Entropy production in chiral symmetry breaking transitions. Chirality 2008, 20, 524–528. [Google Scholar] [CrossRef]
  19. Inoue, Y. Enantiodifferentiating Photosensitized Reactions. In Chiral Photochemistry; CRC Press: Boca Raton, FL, USA, 2004; Chapter 4. [Google Scholar]
  20. Inoue, Y.; Ramamurthy, V. (Eds.) Chiral Photochemistry; Marcel Decker: New York, NY, USA, 2004. [Google Scholar]
  21. Kagan, H.B. Chiral Ligands for Asymmetric Catalysis. In Asymmetric Synthesis: Chiral Catalysis; Morrison, J.D., Ed.; Academic Press: New York, NY, USA, 1985; Chapter 1; Volume 5, pp. 1–39. [Google Scholar]
  22. Fryzuk, M.D.; Bosnich, B. Asymmetric synthesis. An asymmetric homogeneous hydrogenation catalyst which breeds its own chirality. J. Am. Chem. Soc. 1978, 100, 5491–5494. [Google Scholar] [CrossRef]
  23. Kondepudi, D.K.; Kaufmann, R.; Singh, N. Chiral Symmetry Breaking in Sodium Chlorate Crystallization. Science 1990, 250, 975–976. [Google Scholar] [CrossRef] [PubMed]
  24. Asakura, K.; Kobayashi, K.; Mizusawa, Y.; Ozawa, T.; Osanai, S.; Yoshikawa, S. Generation of an optically active octahedral cobalt complex by a chiral autocatalysis. Physica D 1995, 84, 72–78. [Google Scholar] [CrossRef]
  25. Soai, K.; Shibata, T.; Morioka, H.; Choji, K. Asymmetric autocatalysis and amplification of enantiomeric excess of a chiral molecule. Nature 1995, 378, 767–768. [Google Scholar] [CrossRef]
  26. Kondepudi, D.K.; Asakura, K.; Martin, R. Mechanism of Chiral Asymmetry Generation by Chiral Autocatalysis in the Preparation of Chiral Octahedral Cobalt Complex. Chirality 1998, 10, 343–348. [Google Scholar]
  27. Kondepudi, D.K.; Asakura, K. Chiral Autocatalysis, Spontaneous Symmetry Breaking, and Stochastic Behavior. Acc. Chem. Res. 2001, 34, 946–954. [Google Scholar] [CrossRef]
  28. Soai, K.; Shibata, T.; Sato, I. Enantioselective Automultiplication of Chiral Molecules by Asymmetric Autocatalysis. Acc. Chem. Res. 2000, 33, 382–390. [Google Scholar] [CrossRef]
  29. Gehring, T.; Busch, M.; Schlageter, M.; Weingand, D. A concise summary of experimental facts about the Soai reaction. Chirality 2010, 22 (Suppl. 1), E173–E182. [Google Scholar] [CrossRef]
  30. Soai, K.; Kawasaki, T.; Matsumoto, A. Asymmetric Autocatalysis of Pyrimidyl Alkanol and Its Application to the Study on the Origin of Homochirality. Acc. Chem. Res. 2014, 47, 3643–3654. [Google Scholar] [CrossRef]
  31. Blackmond, D.G. Autocatalytic Models for the Origin of Biological Homochirality. Chem. Rev. 2019. [Google Scholar] [CrossRef]
  32. Viedma, C. Chiral Symmetry Breaking During Crystallization: Complete Chiral Purity Induced by Nonlinear Autocatalysis and Recycling. Phys. Rev. Lett. 2005, 94, 065504–065506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Sögütoglu, L.C.; Steendam, R.R.; Meekes, H.; Vlieg, E.; Rutjes, F.P. Viedma ripening: A reliable crystallisation method to reach single chirality. Chem. Soc. Rev. 2015, 44, 6723–6732. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Girard, C.; Kagan, H. Nonlinear Effects in Asymmetric Synthesis and Stereoselective Reactions: Ten Years of Investigation. Angew. Chem. Int. Ed. 1998, 37, 2922–2959. [Google Scholar] [CrossRef]
  35. Swenson, R. Emergent attractors and the law of maximum entropy production: Foundations to a theory of general evolution. Syst. Res. 1989, 6, 187–197. [Google Scholar] [CrossRef]
  36. Swenson, R.; Turvey, M.T. Thermodynamic Reasons for Perception—Action Cycles. Ecol. Psychol. 1991, 3, 317–348. [Google Scholar] [CrossRef]
  37. Martyushev, L.M.; Seleznev, V.D. Maximum entropy production principle in physics, chemistry and biology. Phys. Rep. 2006, 426, 1–45. [Google Scholar] [CrossRef]
  38. Maysman, F.J.R.; Bruers, S. A thermodynamic perspective on food webs: Quantifying entropy production within detrital-based ecosystems. J. Theor. Biol. 2007, 249, 124–139. [Google Scholar] [CrossRef]
  39. Kleidon, A. Nonequilibrium thermodynamics and maximum entropy production in the Earth system. Naturwissenschaften 2009, 96, 653–677. [Google Scholar] [CrossRef] [Green Version]
  40. Niven, R.K. Steady state of a dissipative flow-controlled system and the maximuentropy production principle. Phys. Rev. E 2009, 80, 021113. [Google Scholar] [CrossRef] [Green Version]
  41. Ross, J.; Corlan, A.D.; Muller, S.C. Proposed Principle of Maximum Local Entropy Production. J. Phys. Chem. B 2012, 116, 7858–7865. [Google Scholar] [CrossRef]
Figure 1. The figure shows numerical values assigned to rate constants to obtain numerical solution for the rate equations of the system. Units of volume are assumed to be liters (L) and concentrations mol/L (M). Units of variable parameter II may be thought of as W/m2. If II is thought of as a radiation from the sun, its blackbody temperature is very high compared to the temperature of the system. All rate constants are assumed to have the appropriate units, though not written explicitly.
Figure 1. The figure shows numerical values assigned to rate constants to obtain numerical solution for the rate equations of the system. Units of volume are assumed to be liters (L) and concentrations mol/L (M). Units of variable parameter II may be thought of as W/m2. If II is thought of as a radiation from the sun, its blackbody temperature is very high compared to the temperature of the system. All rate constants are assumed to have the appropriate units, though not written explicitly.
Symmetry 12 00769 g001
Figure 2. Time evolution of chiral species to their equilibrium state in the model reaction scheme when II = 0 and S = T = 0.01 M at t = 0. Concentrations of all other species was set to zero at t = 0. Due to the symmetry of the system, at equilibrium, SL = SD (dashed line) and XL = XD (solid line), so the curves for the two enantiomers overlap.
Figure 2. Time evolution of chiral species to their equilibrium state in the model reaction scheme when II = 0 and S = T = 0.01 M at t = 0. Concentrations of all other species was set to zero at t = 0. Due to the symmetry of the system, at equilibrium, SL = SD (dashed line) and XL = XD (solid line), so the curves for the two enantiomers overlap.
Symmetry 12 00769 g002
Figure 3. Schematic of reaction system. Radiation (shown as hν and denoted as II in the reaction scheme)) is incident on a closed chemical system. The radiation drives a generation–decomposition cycle of enantiomeric species X and other compounds, as shown in A. When II > IIC, the system evolves to one of two asymmetric states, B or C. In state B, the amount of XL > XD, and in state C, XD > XL. The asymmetry is parametrized by α = (XL − XD).
Figure 3. Schematic of reaction system. Radiation (shown as hν and denoted as II in the reaction scheme)) is incident on a closed chemical system. The radiation drives a generation–decomposition cycle of enantiomeric species X and other compounds, as shown in A. When II > IIC, the system evolves to one of two asymmetric states, B or C. In state B, the amount of XL > XD, and in state C, XD > XL. The asymmetry is parametrized by α = (XL − XD).
Symmetry 12 00769 g003
Figure 4. Time evolution of chiral species using equilibrium values for initial conditions when II = 0.003. The solid line represents XL and XD and the dashed SL and SD. Since symmetry is not broken, the amounts of each chiral species are exactly equal, hence the overlapping curves (α = 0).
Figure 4. Time evolution of chiral species using equilibrium values for initial conditions when II = 0.003. The solid line represents XL and XD and the dashed SL and SD. Since symmetry is not broken, the amounts of each chiral species are exactly equal, hence the overlapping curves (α = 0).
Symmetry 12 00769 g004
Figure 5. Time evolution of chiral species in an asymmetric state. Here, II = 0.008 and an asymmetric steady state is reached in about 7000 s. The blue solid and dashed lines represent XD and SD, respectively, and the red solid and dashed lines represent XL and SL, respectively. The black line represents α, which takes on a nonzero value once symmetry is broken.
Figure 5. Time evolution of chiral species in an asymmetric state. Here, II = 0.008 and an asymmetric steady state is reached in about 7000 s. The blue solid and dashed lines represent XD and SD, respectively, and the red solid and dashed lines represent XL and SL, respectively. The black line represents α, which takes on a nonzero value once symmetry is broken.
Symmetry 12 00769 g005
Figure 6. Dependence of α on II. Units of α are M and II are Wm−2. Steady-state values of α are plotted as a function of II. When II > IIC, α takes a positive (XL> XD, shown in green X) or negative (XL< XD, shown in red X) value, depending on the random perturbation that drives the system away from the unstable racemic state α = 0. The blue Xs show the symmetric branch which is unstable above the critical point. In the region II > IIC, α increases in a characteristically parabolic manner.
Figure 6. Dependence of α on II. Units of α are M and II are Wm−2. Steady-state values of α are plotted as a function of II. When II > IIC, α takes a positive (XL> XD, shown in green X) or negative (XL< XD, shown in red X) value, depending on the random perturbation that drives the system away from the unstable racemic state α = 0. The blue Xs show the symmetric branch which is unstable above the critical point. In the region II > IIC, α increases in a characteristically parabolic manner.
Symmetry 12 00769 g006
Figure 7. Rate of entropy production in an asymmetric state when II = 0.008. Units of σ are JK−1L−1s−1. As in Figure 3, steady state is reached at about 7000 s. Although it may appear that σ is approaching 0, it is not so; σ maintains a nonzero value at steady state after symmetry is broken.
Figure 7. Rate of entropy production in an asymmetric state when II = 0.008. Units of σ are JK−1L−1s−1. As in Figure 3, steady state is reached at about 7000 s. Although it may appear that σ is approaching 0, it is not so; σ maintains a nonzero value at steady state after symmetry is broken.
Symmetry 12 00769 g007
Figure 8. Rate of entropy production σ and α as a function of II. The critical value IIC = 0.0040. Beyond this point, α takes on nonzero values and there is a decrease in the slope of σ. Units of II are Wm−2. The dashed arrow indicates a transition from an unstable state, where α = 0, to a stable asymmetric state where α is nonzero. In the region II > IIC rate of entropy production of asymmetric state is lower than that of the symmetric state. The solid arrow shows that transition from an unstable symmetric state to a stable asymmetric state results in the lowering of σ, the rate of entropy production.
Figure 8. Rate of entropy production σ and α as a function of II. The critical value IIC = 0.0040. Beyond this point, α takes on nonzero values and there is a decrease in the slope of σ. Units of II are Wm−2. The dashed arrow indicates a transition from an unstable state, where α = 0, to a stable asymmetric state where α is nonzero. In the region II > IIC rate of entropy production of asymmetric state is lower than that of the symmetric state. The solid arrow shows that transition from an unstable symmetric state to a stable asymmetric state results in the lowering of σ, the rate of entropy production.
Symmetry 12 00769 g008
Table 1. A model photochemically driven reaction scheme in a closed system. The table lists all the elementary steps in the model. The excitation of T to the state TE drives the reaction that generate chiral species XL and XD. Reactions R4a, R4b, R5a, R5b are the autocatalytic steps for XL and XD. Autocatalysis and reaction R6 result in spontaneous breaking of chiral symmetry when the intensity of radiation, II, is above a critical intensity, IIC.
Table 1. A model photochemically driven reaction scheme in a closed system. The table lists all the elementary steps in the model. The excitation of T to the state TE drives the reaction that generate chiral species XL and XD. Reactions R4a, R4b, R5a, R5b are the autocatalytic steps for XL and XD. Autocatalysis and reaction R6 result in spontaneous breaking of chiral symmetry when the intensity of radiation, II, is above a critical intensity, IIC.
Chemical ReactionsNumber
S + II ⇄ TE(R1)
T ⇄ TE(R1a)
S + TE ⇄ XL(R2)
S + TE ⇄ XD(R3)
S + XL ⇄ SL (R4a)
SL + TE ⇄ 2XL(R4b)
S + XD ⇄ SD(R5a)
SD + TE ⇄ 2XD(R5b)
XL + XD ⇄ P + W(R6)
P ⇄ 2S (R7)
W ⇄ 2T (R8)

Share and Cite

MDPI and ACS Style

Kondepudi, D.; Mundy, Z. Spontaneous Chiral Symmetry Breaking and Entropy Production in a Closed System. Symmetry 2020, 12, 769. https://doi.org/10.3390/sym12050769

AMA Style

Kondepudi D, Mundy Z. Spontaneous Chiral Symmetry Breaking and Entropy Production in a Closed System. Symmetry. 2020; 12(5):769. https://doi.org/10.3390/sym12050769

Chicago/Turabian Style

Kondepudi, Dilip, and Zachary Mundy. 2020. "Spontaneous Chiral Symmetry Breaking and Entropy Production in a Closed System" Symmetry 12, no. 5: 769. https://doi.org/10.3390/sym12050769

APA Style

Kondepudi, D., & Mundy, Z. (2020). Spontaneous Chiral Symmetry Breaking and Entropy Production in a Closed System. Symmetry, 12(5), 769. https://doi.org/10.3390/sym12050769

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop