Chiral selectivity vs. noise in spontaneous mirror symmetry breaking

Mirror symmetry breaking bifurcations, that occur in nonlinear chemical systems leading to final chiral states with very large enantiomeric excess, can be exploited as an eﬃcient chiral signal selector for even the smallest chiral polarizations. This eﬀect of the chiral polarization requires the system’s capacity for overcoming thermal noise, which is manifested as fluctuating reaction rate constants. Therefore, we investigate the chiral selectivity across a range of tiny parity-violating energy diﬀerences (PVED) in the presence of inevitable non-equilibrium temperature fluctuations. We use a stochastic diﬀerential equation simulation methodology (Ito process) that serves as a valuable tool in open systems for identifying the thresholds at which the chiral force induces chiral selectivity in the presence of non-equilibrium temperature fluctuations. This approach enables us to include and analyze chiral selectivity in the presence of other types of fluctuations, such as perturbations in the rate of fluid flow into and out of the reactor and in the clamped input concentrations. These concepts may be of practical interest ( i.e. , spontaneous deracemizations) but are also useful for a better understanding of the general principles governing the emergence of biological homochirality.


Introduction
Enantiomerism 1 refers to molecular geometrical structures that break parity invariance 2 The concept is related to chirality, which involves the presence of right-and left-handed objects.Chiral selectivity is the ability of a molecule or a system to recognize and differentiate between chiral forms or enantiomers.Biological homochirality, or one-handedness, refers to the uniform predominance of one and the same chiral sign of the building blocks of the functionalized biological polymers such as the L-amino acids in proteins and the D-sugars in DNA, supporting the basic functions of life such as catalysis, metabolism, self-reproduction and the ability for evolution. 3he homochirality of amino acids and sugar molecules in biosystems is a necessity for life, 4,5 and has its most striking experimental manifestation regarding the question of the origin of biological homochirality.
The homochiral property of life stands in striking contrast to the characteristics of artificial chemical processes that result in chiral species.In laboratory or industrial chemistry, when there is no external chiral polarization, the final output of a reaction is typically racemic, meaning it contains an equal proportion of both enantiomeric species.However, there are a few wellknown chemical experiments, such as the Soai reaction 6 and Viedma deracemization, 7 where the final chiral composition is often nearly homochiral rather than racemic.This phenomenon occurs in nonlinear chemical networks that exhibit enantioselective autocatalysis leading to a bifurcation when the racemic non-equilibrium stationary state (NESS) becomes unstable.In these cases, fluctuations drive the system stochastically towards one of the two possible enantiomorphic outputs.[10][11][12] However, due to fluctuations and in the absence of any chiral polarization, there is a stochastic distribution of chiral signs around the ideal racemic composition.It is important to emphasize that SMSB represents a genuine bifurcation scenario, which can be highly sensitive to even slight chiral biases. 13This sensitivity stands in contrast to the strong chiral polarizations that are required, for instance, in classical asymmetric synthesis. 149][20] Nevertheless, it is crucial to acknowledge that in real nonlinear chemical systems, any chiral polarization competes with the unavoidable thermal and environmental fluctuations.Therefore, it is likely that there exists a critical threshold of chiral polarization capable of surpassing the racemizing effects of these fluctuations.To estimate this threshold from first principles, we assume for an ideal solution with perfect mixing that, on the one hand, the changes in the reaction rate constants are due to the thermal fluctuations and, on the other hand, that there exists a tiny but permanent chiral polarization, namely the parity-violation energy difference (PVED) between the enantiomers. 21,22][25][26][27][28][29][30][31] Although PVED has not yet been detected in typical chemical experiments 29,30,[32][33][34] a logical connection exists between physical reality and biological chirality. 35,36However, the planned design and construction of ultrahigh resolution instruments for vibrational spectra, and the use of heavy atom chiral complexes, promises the possibility of successful PVED measurements being carried out in the near future.
In this paper, we follow a different approach, namely the use of SMSB [37][38][39][40] as a powerful amplifier of chiral polarizations.It is well-known that some enantioselective autocatalytic systems lead to outcomes yielding large deviations from the racemic composition.A stochastic distribution of final chiral signs between successive experiments is expected on the average, 41 although there was some initial doubt for both the Viedma deracemization 42 and the Soai reaction. 43Therefore, for any reasonable scenario of chemical evolution, and for a relatively large number of analogous systems, the ultimate net consequence should be the racemic composition.Nevertheless, since Kondepudi's pioneering work, 44 it has been assumed that weak neutral currents, extended over sufficiently long periods of time, during which there is a slow increase of the substrate concentration, could well lead to a deterministic and robust chiral bias in nonlinear chemical networks involving enantiomers, even in a randomly fluctuating environment. 45Motivated by this possibility, we developed a preliminary model 46 demonstrating that the very tiny 47 but permanent parity violating energy difference (PVED) 22 can be sufficient to deracemize a macroscopic and out of equilibrium chemical system, even in the presence of unavoidable and permanent thermal fluctuations.These calculations were based on nonlinear stochastic differential rate equations, where the reaction rate constants involving the enantiomers are sensitive to the thermal fluctuations, and as well to the PVED bias.It was shown that when the amplitude of the thermal fluctuations times the activation enthalpy is of the same order of magnitude as the PVED bias, a resilience of the parity-violation-induced chiral selectivity to nonequilibrium temperature fluctuations is observed.
Encouraged by this initial success, and being aware that there is always a risk that the influence of the fluctuations might not be taken into account sufficiently, we clarify here a number of important points not addressed or covered in our previous work.
The purpose of this paper is to focus on the details of the stochastic histories of the enantiomeric chemical compositions before and after the chiral symmetry breaking transition, showing the influence of the presence or absence of the PVED bias.Our findings demonstrate that in close proximity to the critical threshold, an increased number of simulation trials leads to a statistically more significant differentiation between the stochastic and deterministic outcomes.In the present study, we employ an improved stochastic simulation methodology that serves as a valuable tool for identifying the thresholds at which chiral forces, such as the PVED, induce chiral selectivity that remains resilient in the face of non-equilibrium temperature fluctuations in open systems.These findings shed light on the interplay between chiral forces and fluctuation effects, providing insights into the conditions under which chiral selectivity can be sustained in dynamic chemical systems.To accomplish this, we define and implement a different mathematical approach for the numerical simulation of stochastic differential rate equations which is based on Ito and Wiener processes, 48 and we include the complete untruncated thermal noise spectrum.This distinct approach will allow us to consider fluctuating boundary conditions, fluctuations in the system volume and also fluctuations in the input fluid flow rates to the reactor, and distinct input/output reactor configurations or reactor architectures.To this end, we employ a stochastic version of an open-flow, fully reversible, Frank model. 49We choose this dynamical model because it is now well known that its basic ingredients, namely enantioselective autocatalysis and mutual inhibition of opposite enantiomers, lie at the heart of the Soai reaction, the first experimentally proven absolute asymmetric synthesis. 6,50his paper is organized as follows.In Section 2 we review briefly the general expression for the amplitude of the rateconstant fluctuations, the chiral selectivity criterion and the open flow reaction scheme employed.We then present the set of coupled Ito stochastic differential equations.In the present approach, the entire range of temperature fluctuations are used in the numerical simulations, in contrast to the approximation used previously, and this is commented briefly in Section 3. The inclusion of chiral bias due to PVED in the vector Ito equations is presented in Section 4. The results of varying the relative magnitude of the temperature fluctuations with respect to the fixed PVED bias are presented in Section 5 for a range of bias values.Fluctuations in the boundary conditions associated with the system volume, the fluid flow rates, and the fixed input concentrations are treated in Section 6. Conclusions are given in Section 7.

Chiral selectivity criterion and a stochastic open-flow model
where DH ‡ is the activation enthalpy of the reaction.The (nonequilibrium) temperature fluctuations dT are distributed according to a Gaussian distribution with mean T and standard deviation s 51 where k B is the Boltzmann constant and C V the constant volume specific heat.This relation allows us to characterize the T-fluctuation term in eqn (1) in terms of the root-meansquare (rms) value h(dT) 2 i 1/2 as follows.We estimate the ratio in eqn (1) according to: ZðtÞ; so that eqn (1) can be expressed as where the amplitude x i of the fluctuation in the rate constant k i is given by and Z i (t) is a Gaussian white noise.As the transition enthalpy can only be positive DH ‡ i 4 0 we define a positive noise amplitude x i 4 0. We lose no generality in defining x in this way since the white noise factor it multiplies has both positive and negative fluctuations that increase or decrease the corresponding transition enthalpy factor.Both the amplitude x i and the noise Z i depend on the individual ith reaction; the forward and reverse steps are treated as independent reactions.Note moreover that x i in eqn (4) is a product of two independent contributions: the first one depends on the transition enthalpy divided by the mean temperature, whereas the second factor involves the relative rms temperature fluctuation.The former is a chemical property of each single reaction path, while the latter encodes the bulk thermodynamic properties of the system, via the specific heat C V at constant volume (a physical attribute).
Denote by DDG ‡ i 4 0 the minute energy difference between the two enantiomers due to their interaction with a chiral polarization, exerted for example by chiral media, or by a natural chiral physical force.Then, consider the case when g corresponds to the PVED bias in the ith rate constant.When this bias is greater than the rms temperature fluctuations times the change in enthalpy factor associated with that rate constant: then we would rightfully expect that the chiral bias g can overcome the thermal noise x i , and therefore select deterministically the final stable chiral outcome, provided the inequality g Z x i holds for all the reactions i in which the PVED bias intervenes.The relation eqn (5) implies a signal-to-noise (S/N) ratio g/x i Z 1, in which the PVED bias plays the role of the constant applied ''direct current'' or DC signal, and the temperature fluctuations comprise the ''alternating current'' or AC background noise.See Appendix A for derivation of particular cases of this chiral selectivity criterion for the kinetic model described below.We test this chiral selectivity criterion eqn (5) with the following model.It involves four species A, L, D, P and ten one-way reactions eqn ( 6)- (8), in addition to the five pseudoreactions or input/output flow terms eqn ( 9)-( 13), as defined by the following transformations: where k f = q/V, the volumetric flow rate q is in liters per second and V is in liters.Species A is an achiral precursor and P is an achiral product: the heterodimer LD = DL is the achiral (meso form) dimer of the two chiral enantiomers L and D.
The individual forward and reverse reaction rate constants in eqn ( 6) and ( 7) are constrained by the equality of the equilibrium constants for enantioselective autocatalysis and the direct production: that is, and this relates the differences in the free energies of activation of the forward and reverse reactions for autocatalysis and direct production: This relation follows from inserting the expression for the rate constants k into the above constraint, and using transition state theory (the Eyring-Polanyi equation for k).This is the most general statement that can be made regarding the activation energies.
The stochastic reaction rate equations that follow from the above transformations eqn ( 6)-( 13) and using eqn (3), are written as a set of coupled Ito stochastic differential equations, see ref. 48.The noise terms follow from substituting the expression for the fluctuating reaction rate constants k i ,

PCCP Paper
This journal is © the Owner Societies 2023 eqn (3), into the differential reaction rate equations that follow from the above transformations.The result is then expressed as Ito stochastic differential equations, using the definition of Wiener processes to write the noise terms as indicated below.
We can separate each equation into its deterministic (pdt) and random (pdW j (t)) parts in order to cast the equations as a vector Ito process, 48 using its attendant standard notation, as follows: There are ten independent one-way transformations in eqn ( 6)-( 8), hence these stochastic equations will depend on ten independent Wiener processes {W i } i=1

10
, where dW i (t) = Z i (t)dt; 48 one such Wiener process (and Gaussian white noise Z i ) for each one of the ten independent one-way reactions.In the above, we have taken all x i = x to be equal, a reasonable first approximation; see eqn (4).We comment on the scaling.The reaction rates (the coefficients of the dt) scale as the inverse volume BV À1 .The noise terms (see coefficients of the dW) are given by the product of x times the reaction rates, and the noise amplitude x scales as the inverse square root of the heat capacity at constant volume: x B (1/C V ) 1/2 .The overall product, which gives the noise term thus scales as BV À1 (1/C V ) 1/2 .The constant volume heat capacity scales with the volume C V B V. Thus the overall noise terms scale as BV À3/2 .Thus relative to the deterministic contributions (the coefficients of dt) to the rate equation which scale as BV À1 , the noise terms (coefficients of the dW) therefore scale as BV À1/2 .
We remark that an attempt to estimate the individual noise amplitudes x i would make sense if one were to consider using the rather more involved and phenomenologically realistic kinetic models of say, the Soai reaction, 52 and provided of course that one could reasonably estimate the values of the transition enthalpies DH ‡ i of the individual reactions belonging to those kinetic and thermodynamic models.In that case, one could use the van't Hoff equation 53 for this purpose, and would therefore need to know the equilibrium constant of a given reaction for at least two different temperatures, to be able to calculate the transition enthalpy of that reaction.

Thermal noise statistics
A major important difference between the stochastic differential equations used in ref. 46, and the vector Ito process used here, has to do with the treatment of the noise.For purposes of the numerical simulations, the noise was generated as continuous random functions, the latter obtained by interpolating sequences of random numbers generated from a normal distribution.These continuous random functions were then included in the differential equations.We also truncated the spectrum of the relative temperature fluctuations dTðtÞ T so as to lie within one standard deviation of the mean temperature T.
In the present paper, we employ vector Ito process to model the stochastic equations and use independent Wiener processes for the noise terms, so no such restriction is imposed on the noise.Whereas the former truncation accounted for approximately 67% of all the fluctuations, here all the fluctuations are included, see Fig. 1.The ''missin'' fluctuations, those beyond one standard deviation, make up for the remaining 33% of the noise spectrum.Recall that the probability distribution for nonequilibrium temperature fluctuations dT is given by ref. 46 where j is the normal frequency function 54 and the standard deviation s = dT rms = h(dT) Using the full range of fluctuations will have an impact on the ability of the PVED bias g to overcome the thermal noise and so select a final chiral sign (see Section 5), but is not as drastic as one might suppose, since the very largest fluctuations, which are those lying beyond three standard deviations, 3s, from the mean, have exponentially vanishing probabilities, and so are extremely rare events, see Fig. 1.Moreover, as white Gaussian noise includes all frequencies, and with equal weighting (its frequency spectrum is flat), there always exists arbitrarily large frequencies o such that oDt 4 1 is satisfied, where Dt is the smallest time step used in the numerical solution of the stochastic differential equations.Compliance with this inequality implies that statistical randomness between successive simulation time steps is assured.

Including PVED bias
To include a chiral bias g due to a small energy difference between enantiomers, in the above stochastic differential rate equations, we make the following substitutions in the first two equations eqn ( 14) and ( 15): and in only those terms that contain [L], but in no term containing [D]: This chiral bias favors the D-enantiomer over the L (changing the sign from Àg to +g in the above favors instead the L-enantiomer over the D).This leads to the modified stochastic equations for [A] and [L]: The remaining two stochastic equations eqn ( 16) and ( 17) are unchanged.We can then carry out simulations of the four coupled stochastic differential equations eqn ( 16), ( 17), ( 22) and ( 23) to assess the ability of a range of PVED biases g to select deterministically final stable scalemic states in the presence of random temperature fluctuations x, according to eqn (5).The simulations of the set of Ito stochastic differential equations have been performed using the Euler-Maruyama method, as defined and implemented in Mathematica v13.
5 Fluctuating reaction rate constants: stochastic vs. deterministic regimes The model parameters are taken as follows: and the flow rate k f = 10 À3 .From eqn ( 6) and ( 7) these imply an equilibrium constant K eq = 10 5 .We initiate all the simulations on the idealized racemic composition: [L] 0 = [D] 0 = 10 À3 , with [P] 0 = 10 À6 and [A] 0 = [A] in .These values locate the system on the unstable racemic branch in the absence of chiral bias.We first generate multiple runs of the system subject to the temperature fluctuations alone (no PVED bias: g = 0), and for strict racemic initial conditions.This yields a set of statistically independent random histories, a typical example of which is shown in Fig. 2.After an induction period the system bifurcates, and the times at which the concentrations of the two enantiomers separate from each other is spread out in time, this dispersion is due to the fluctuations in the rate constants.We can map these bifurcations to histograms of the final enantiomeric excesses ee = ([L] À [D])/([L] + [D]) by taking time-slices through the bifurcation curves (at the final times at t = 5000 s) in these sets of histories.Note, for each individual run, there is a majority (on the upper branch) and a minority (on the lower branch) enantiomer, which can be either the L or the D enantiomer, respectively, and we calculate the corresponding ee for this run.As each outcome is essentially homochiral (already for t = 5000 s, there are seven orders of magnitude difference in the majority vs. minority enantiomer concentrations) this single run either yields ee D À1 (when the D is the majority enantiomer) or else ee D +1, (when the L is the majority enantiomer).This procedure is then repeated for the entire set of all n runs, and so we calculate P n j ee j ðþÞ and P n j ee j ðÀÞ, where ee j (AE) = AE1, where +1 holds for the majority D-enantiomer, and À1 holds for the majority L-enantiomer, respectively.We plot the results in histograms to show the distribution of the set of all outcomes.The data slice is taken at t = 5000 s.Next, we include the PVED bias g 4 0 in the stochastic equations and make multiple runs, again with racemic initial conditions.For each value of the PVED bias chosen, we test the chiral selectivity against various noise amplitudes x, for noise greater than, equal to, and less than the bias: x 4 g, x = g, x o g.We map the outcomes to histograms of final enantiomeric excesses ee as outlined above.
Here in contrast, after the induction period the times at which the concentrations of two enantiomers separate from each other is sharply focused in time, see Fig. 3.The PVED bias has overcome the thermal fluctuations in the rate constants.See Fig. 4 for characteristic examples sampled over a wide range of g and also Table 1 for a guide to the characterization of the distinct outcomes.Very recent calculations of PVED for cinnabar crystals, 55 lead to values as large as g = 2.2 Â 10 À12 , which we include in our analysis below.Some comments regarding these results, with respect to the those obtained in ref. 46, are in order.There is now a qualitative and quantitative pattern of statistical uniformity in the outcomes of the numerical results with respect to the range of g bias values tested in that: (i) all the g's tested here select the final state deterministically when the noise level is an order of magnitude smaller than the bias: x o g (see Panel (a) in Fig. 4), (ii) for noise of the same order of magnitude as the bias x = g, the selection is still biased strongly, but is slightly skewed (see Panel (b) in Fig. 4), and for noise an order of magnitude larger than the bias x 4 g, the outcomes are bimodal, no statistically significant selection of a final scalemic state is possible, see Panel (c) in Fig. 4.Here we solve numerically for the four coupled Ito stochastic differential equations eqn ( 16), ( 17), ( 22) and (23).The main mathematical difference between this approach and the previous one is now the random T-fluctuations are taken from a complete normal distribution, and are not restricted to lie within (AEs) one standard deviation of the mean (see discussion in Section 3 on this point).That is, all the temperature fluctuations beyond one standard deviation from the mean are now included.This translates into a diminished chiral selectivity when the fluctuations are greater than the bias x 4 g, which is natural and is to be expected.The largest temperature fluctuations (meaning those that lie beyond three standard deviations from the mean) diminish the chiral selectivity of the PVED bias, but they are also less likely than the smaller fluctuations that lie within one standard deviation of the mean temperature, see Fig. 1.Further details of the individual numerical simulations are compiled in Table 2.For a given chiral bias g 4 0, noise amplitude x 4 0 and number of statistically independent runs n, we obtain the total numbers of final homochiral L and D outcomes, where L + D = n.
From this we calculate the net enantiomeric excess heei n = (L À D)/(L + D) averaged over the n runs, and also the fraction of homochiral D-outcomes p D = D/(L + D).We compare the latter to the 99% confidence interval (C.I.) for expecting a racemic outcome, i.e., when p D = 1/2.This C.I. eqn (24) depends on p, the confidence level and also on n.Note, in the so-called ''severe'' region of parameter space, for which the strength of the chiral bias is less in magnitude than the noise amplitude the runs, qualitatively similar to that in Fig. 4(c).See text for how these outcomes are mapped to histograms and see Table 1 for their interpretation.16), ( 17), ( 22) and ( 23) with bias g = 10 À18 with x = 10 À19 .The bias g reduces the spread in the bifurcation times.Compare to Fig. 2. Computing the enantiomeric excess ee at time slice t = 5000 s leads to a unimodal histogram for all the runs, qualitatively similar to that in Fig. 4(a).See text for how these outcomes are mapped to histograms and see Table 1 for their interpretation.g o x, we have also further increased the number of runs n from 100 to 500 and then to 1000 in order to ascertain whether the final stationary outcomes can be regarded as either racemic or else biased.We use the 99% confidence interval (C.I.) for expecting a racemic outcome as criterion: where p is the proportion of interest (e.g., for racemic proportion p = 1/2), n is the sample size, a the desired confidence, and z 1À a 2 the Z-value for the desired confidence, z 1À a 2 ¼ 2:57 for 99% confidence. 56These confidence interval bounds are computed and listed in the right-hand column of Table 2.The table is based on time-slices through the bifurcation curves at the final times at t = 5000 s.When the results in Fig. 4 are plotted against the logarithms of the bias g and of the noise x, a simple pattern emerges as shown in Fig. 5. Indeed, there is a line, or separatrix, that divides the region of deterministic chiral outcomes from the region of stochastic ones, for the range of bias and noises considered.This is a linear relation between log(g) and log(x).We also see the effect of using the full noise noise spectrum here in contrast to the simulations of ref. 46.Also, for noise levels x = 10 g an order of magnitude greater that the bias, increasing the number of runs n narrows the confidence interval CI eqn (24), so what was a racemic outcome for low n can become a chiral one for larger n, since now the fraction p D lies outside the 99% CI, see Fig. 5.
Table 1 Characterization of the various outcomes: random, biased or deterministic, in terms of relationships between the sums of the enantiomeric excesses of the individual homochiral runs ee j (AE) = AE1.A statistically more precise discrimination between the random or biased outcomes is provided by estimating the confidence intervals (C.I.) for having a racemic outcome, as discussed in the text; see also Tables 2 and 3 16), ( 17), ( 22) and ( 23).consequence of lowering the equilibrium constant by a factor of 10 2 , by choosing the rate constants such that k d = 10 À6 and k a = 1 keeping all other rate constants and parameters fixed, then K eq = 10 3 .Then, mirror symmetry is still broken by temperature fluctuations, but now the reactions are less exergonic, or more reversible, than before.The results are tabulated in Table 3.For the same values of chiral bias g, noise amplitude x and number of runs n, the chiral selectivity of the bias is enhanced, for a lower value of the equilibrium constant, with respect to the results in Table 2.This indicates that when the reactions are more reversible (i.e., are less exergonic) the greater is the chiral selectivity of the bias g with respect to a given noise level x (see Table 3).We again take time-slices through the bifurcation curves at the final times at t = 5000 s.The relationship between the number of experiments needed to yield an adequate CI in order to be able to distinguish between stochastic or deterministic outputs is a point that cannot be detected using common chemical kinetic methods.Despite the fact that the ratio g/x = 1 defines the transition from stochastic to deterministic regimes (see Tables 2 and 3 and Fig. 4), the CI (24) indicates the efficiency of the system towards chiral selectivity when the ratio between chiral induction and noise amplitude changes are 41; a large number of experiments, or in a origin of life scenario a large set of separated similar systems, is necessary to detect the chiral selection effect in SMSB.However, when g/x c 1 only a small number of experiments, or a few separated systems, will already lead to the same final chiral sign selected by g.

Role of fluctuations in the boundary conditions: volume, flow rates and input concentrations
Fluctuations in the input and output rate dk f (see the pseudoreactions in eqn ( 9)-( 13)) can arise due to either fluctuations in Fig. 5 Left: Influence of the calculation technique on the system output.Open circles or squares: truncated fluctuations, filled circles and squares: unrestricted fluctuations.Red round symbols: racemic output, green square symbols: chiral output.For an equivalent number of trials the stochastic Ito process is less selective than the calculations carried out in ref. 46.This can be seen at log(noise) = À11 and log(bias) = À12.For n = 100, Ito gives a racemic output (red filled dot) while truncated calculation gives a chiral output (open green square, n = 120).The same effect is also seen at log(noise) = À16 and log(bias) = À17 where truncated calculation gives a chiral output (open green square) while a red filled dot (racemic output) is expected from the Ito process.The dashed line (log(noise) = log(bias) + 0.5) is parallel to the first diagonal.It appears as a sensitive zone (separatrix between racemic and chiral domains from the Ito process) where the noise is about three times larger than the bias.Right: Statistical analysis of a series of parallel runs using a 99% confidence interval (CI) for expecting a racemic output in the presence of a noise 10 times greater than the bias.Results show that a high run number helps to unveil the tiny (but permanent) effect of the chiral bias which, otherwise would remain undetected (100 vs. 1000 runs).It is also shown that a higher bias is needed to overcome the randomizing effect of a greater noise level (bias 410 À14 at 500 runs).
Table 3 Results of simulations for an equilibrium constant K eq = 10 3 , two orders of magnitude smaller than that corresponding to Table 2. g: the chiral bias in the rate constants; x: the amplitude of the temperature fluctuations, n: the number of statistically independent runs; L and D: the total number of homochiral outcomes for the L-and D-enantiomers, respectively; heei n = (L À D)/(L + D): the net enantiomeric excess averaged over the n runs; p D = D/(L + D) is the fraction of homochiral D-outcomes for n.Final column gives the 99% confidence interval (C.I.) for expecting a racemic outcome that is, for p D = 1/2.Note there is always a chiral output.Compare to Table 2  the fluid flow rate itself dq and or to fluctuations in the system volume dV: ) These are independent from the internal temperature fluctuations.We may include them in the above Ito formulation by making the following substitution in all the terms involving the flow rate: where Z f is a Gaussian white noise, with an amplitude x f and introducing the single Wiener process dW f [t] = Z f (t)dt for the fluctuating term.Then using eqn (27) in eqn ( 16), ( 17), ( 22) and ( 23) leads to the following replacements, where X = A in , A, L, D, P: Whereas the rate constant fluctuations dk i (due to the temperature fluctuations) for the ten individual one-way reactions eqn ( 6)-( 8) are statistically independent, the open-flow terms eqn ( 9)-( 13) will fluctuate coherently and in unison.This is because the input and output channels depend on one and the same flow rate constant k f .In particular, this means the fluctuations in the enantiomeric concentrations [L] and [D] due to the random flow rate dq and or volume fluctuations dV are racemic, and also of first order.So, no mirror symmetry breaking is expected due to these types of racemic, or ''diagonal'', fluctuations.We have confirmed this with simulations.
The outcome is similar to the Panel (a) in Fig. 4: a deterministic result.By way of example, even for very large amplitude fluctuations in the flow rate k f of the order of x f = 0.16 (eqn ( 27)) and for a minuscule chiral bias g = 10 À18 , and in the presence of T-fluctuations, the chiral bias still selects deterministically for g 4 x, while there is a statistically significant bias when x = 10 À18 and then no significant bias for x = 10 À17 , which is an order magnitude greater that the fixed chiral bias.Indeed, the qualitative pattern is that of the Panels (a) to (c) in Fig. 4. So, it is statistically indistinguishable from the case where the flow rate does not fluctuate.This is a consequence of the input and output flow scheme or flow architecture, for which all species concentrations are removed at the same rate in one unique output channel and where only the achiral substrate A is input to the reactor, see Fig. 6.
A distinct and realistic open-flow set-up from the chemical point of view, is to arrange for both enantiomers to flow into the reaction tank in a common shared channel and with an average racemic composition [L] in = [D] in .We can arrange for the achiral species A flow to in on a separate channel, see Fig. 6.There are thus two input channels, and a single common output channel carrying all four species A, L, D, P. In this case, the pseudo-reactions corresponding to the input flows are as follows: where [X] in denotes the fixed on average, or clamped, input concentration, and we define the individual input flow rates as follows: where the fluxes f X,in satisfy the constant volume constraint f out ¼ P j¼A;L;D f j;in . 57We consider fluctuations on the inflows of the enantiomers, see Fig. 6.From eqn (30), ( 32) and ( 33 as follows using substitutions analogous to eqn (28): is sufficiently large.The constant volume constraint is This can be satisfied by taking the flux fractions to be f for the two input channels, this implies k fA = k fL = k fD = k f /3.Examples of the competition between chiral bias and noise on the enantiomer input channel are shown below in the Fig. 7-9.
For the range of PVED biases g considered above, and for fixed thermal noise levels x 4 0, we can estimate the maximum amplitudes of the enantiomeric input line noises z 4 0 that can be tolerated before washing out the PVED selectivity.Deterministic outcomes can tolerate an external noise level z up to approximately 10 3 greater than the chiral bias g, and in the presence of temperature fluctuations, whereas the chiral bias gets obliterated for external noise levels on the order of approximately 10 5 greater than the chiral bias, see Fig. 7-9.

Concluding remarks
In our open flow stochastic Frank model, subject to a chiral bias g and temperature fluctuations x, the ratio g/x = 1 (see Section 4), defines the transition from stochastic to deterministic regimes.This condition represents a critical threshold for chiral sign selection in SMSB.When the ratio of chiral bias to temperature fluctuations is greater than unity g/x 4 1,  a calculated confidence interval is used to indicate the efficiency of the system towards chiral selectivity.A larger number n of experiments is necessary to detect statistically the chiral selection effect in SMSB.This relationship, depending on the system parameters, between the number of experiments needed to distinguish between stochastic or deterministic outputs, is a point that cannot be detected using common chemical kinetic methods.Moreover, since each simulation describes the behavior of a very small volume element, the confidence interval gives a description of what occurs in a macroscopic volume resulting from addition of large numbers of small volume elements.Therefore, larger sets of small volume reactors lead to a greater confidence for the distinction between stochastic or deterministic outcomes.This result is interesting because the chemical systems description of abiotic evolution are those based on the cooperative interaction between many micro systems (e.g.ref. 58).As expected, when g/x c 1 already a small number of experiments will lead to the g-selected final chiral sign.A large non-equilibrium system containing many interacting chiral species may behave as a collection of many microsystems.In this regard, a recent work shows a phase transition towards a homochiral state is likely to occur as the number of chiral species becomes large, see ref. 59.Since it is expected that larger amplitude chiral fluctuations in real systems (imperfect mixing, limited diffusion rates. ..), occur distinctly from those originating from the temperature fluctuations in the values of the rate constants, the experimental detection of a possible PVED effect calls for a difficult and challenging experimental procedure able to eliminate or suppress the sources of the chiral fluctuations that are not due the temperature fluctuations in thermally isolated systems.
Our results show that for the range of PVED biases g considered above, and for fixed thermal noise levels x 4 0, there exists a maximum tolerable amplitude of the enantiomeric input line noises x.In the presence of temperature fluctuations, deterministic outcomes can tolerate an external noise level x of up to approximately 10 greater that the chiral bias g.However, the PVED selectivity is washed out for external noise levels x about 10 2 times greater than the chiral bias.These refer to the thermal noise levels that can be tolerated by the chiral bias.By contrast, estimates concerning the chiral noise levels (chiral fluctuations on the ideal racemic composition that feeds into the reactor, see right hand side of Fig. 6), that can be tolerated by the chiral bias, and in the presence of fixed level of thermal noise, are notably more generous.
Results on permanent physical chiral forces overcoming chiral fluctuations are widely reported, and even in the case where the chiral force is at the order of cryptochirality. 41ryptochirality refers to minute perturbations of the chiral symmetry which remain undetected (hence, hidden chirality) in classical physical measurements or in enantioselective reactions, but which are amplified in the bifurcation in these polarized (chiral biased) spontaneous mirror symmetry breakings.Note this definition encompasses not only PVED, but also other experimentally undetectable sources of chiral polarizations, be they of a structural, an environmental or a statistical nature.This means that the presence of a chiral polarization below the experimental limits for detecting chirality by physical methods can determine the output chiral sign.6][17] Note, however that contrary to such cryptochirality, whose sign can be random, the effect of PVED must always be present, and precisely because it adds to, or subtracts from, biases exerted by other chiral forces, its detection could be indirectly made by means of an accurate statistical analysis of multiple outputs of cryptochiral enantiomorphic polarizations.A detailed discussion of cryptochirality can be found in ref. 41.
Finally, it is worthwhile pointing out that this fact is contrary to the older opinions and models which supposed that biological homochirality could arise from stochasticity.Here we demonstrate that low exergonic enantioselective autocatalysis can lead to SMSB in a deterministic way.

Fig. 1
Fig. 1 The standard normal frequency function j, eqn (18).The approximation used in ref. 46 corresponds to relative fluctuations dT dT rms

Fig. 2
Fig. 2 Random outcomes.Example of n = 20 independent runs of the stochastic equations eqn (14)-(17) without PVED bias (g = 0); with x = 10 À19 .The symmetry breaking bifurcation is spread out in time by the noise, but all individual histories converge to the same asymptotic state.With no bias, both the upper and lower sets of branches are equally populated by L and D outcomes on average.Computing the enantiomeric excess ee for each individual run at time slice t = 5000 s leads to a bimodal histogram for all the runs, qualitatively similar to that in Fig.4(c).See text for how these outcomes are mapped to histograms and see Table1for their interpretation.

Fig. 3
Fig. 3 Deterministic outcomes.Example of 20 independent runs of the stochastic equations eqn (16), (17), (22) and (23) with bias g = 10 À18 with x = 10 À19 .The bias g reduces the spread in the bifurcation times.Compare to Fig.2.Computing the enantiomeric excess ee at time slice t = 5000 s leads to a unimodal histogram for all the runs, qualitatively similar to that in Fig.4(a).See text for how these outcomes are mapped to histograms and see Table1for their interpretation.
for two independent Wiener processes W Lin , W Din .z 4 0 represents the amplitude of noise in the enantiomer concentrations[L] in and[D] in .Of course, we choose racemic input conditions: k fL = k fD ; equal input flux fractions for both enantiomers, see eqn (32) and(33), and [L] in = [D] in , so that any deviations from mirror symmetry coming from these terms eqn (34) and (35) are thus due strictly to the independent noises on the common L,D-input channel, as quantified via z 4 0. Since W Lin , W Din are independent Wiener processes, the fluctuations in the monomial terms k fL[L] in and k fD[D] in are therefore uncorrelated.So, they can give rise to net chiral fluctuations which can degrade the chiral selectivity of the PVED bias g provided the input noise amplitude z 4 0