Limits and tradeoffs in the control of autocatalytic systems

Despite the complexity of biological networks, we find that certain common architectures govern network structures. These architectures impose fundamental constraints on system performance and create tradeoffs that the system must balance in the face of uncertainty in the environment. This means that while a system may be optimized for a specific function through evolution, the optimal achievable state must follow these constraints. One such constraining architecture is autocatalysis, as seen in many biological networks including glycolysis and ribosomal protein synthesis. Using a minimal model, we show that ATP autocatalysis in glycolysis imposes stability and performance constraints and that the experimentally well-studied glycolytic oscillations are in fact a consequence of a tradeoff between error minimization and stability. We also show that additional complexity in the network results in increased robustness. Ribosome synthesis is also autocatalytic where ribosomes must be used to make more ribosomal proteins. When ribosomes have higher protein content, the autocatalysis is increased. We show that this autocatalysis destabilizes the system, slows down response, and also constrains the system’s performance. On a larger scale, transcriptional regulation of whole organisms also follows architectural constraints and this can be seen in the differences between bacterial and yeast transcription networks. We show that the degree distributions of bacterial transcription network follow a power law distribution while the yeast network follows an exponential distribution. We then explored the evolutionary models that have previously been proposed and show that neither the preferential linking model nor the duplication-divergence model of network evolution generates the power-law, hierarchical structure found in bacteria. However, in real biological systems, the generation of new nodes occurs through both duplication and horizontal gene transfers, and we show that a biologically reasonable combination of the two mechanisms generates the desired network.

Experimental study on inertial effects in liquid-solid flows

This thesis presents experimental measurements of the rheological behavior of liquid-solid mixtures at moderate Reynolds (defined by the shear rate and particle diameter) and Stokes numbers, ranging from 3 ≤ Re ≤ 1.6 × 10³ and 0.4 ≤ St ≤ 195. The experiments use a specifically designed Couette cylindrical rheometer that allows for probing the transition from transporting a pure liquid to transporting a dense suspension of particles. Measurements of the shear stress are presented for a wide range of particle concentration (10 to 60% in volume) and for particle to fluid density ratio between 1 and 1.05. The effective relative viscosity exhibits a strong dependence on the solid fraction for all density ratios tested. For density ratio of 1 the effective viscosity increases with Stokes number (St) for volume fractions (φ) lower than 40% and becomes constant for higher φ. When the particles are denser than the liquid, the effective viscosity shows a stronger dependance on St. An analysis of the particle resuspension for the case with a density ratio of 1.05 is presented and used to predict the local volume fraction where the shear stress measurements take place. When the local volume fraction is considered, the effective viscosity for settling and no settling particles is consistent, indicating that the effective viscosity is independent of differences in density between the solid and liquid phase. Shear stress measurements of pure fluids (no particles) were performed using the same rheometer, and a deviation from laminar behavior is observed for gap Reynolds numbers above 4× 10³, indicating the presence of hydrodynamic instabilities associated with the rotation of the outer cylinder. The increase on the effective viscosity with Stokes numbers observed for mixtures with φ ≤ 30% appears to be affected by such hydrodynamic instabilities. The effective viscosity for the current experiments is considerably higher than the one reported in non-inertial suspensions.

PP backward elastic scattering between 0.70 and 2.37 GeV/c

̄pp backward elastic scattering has been measured for the cos θ_cm region between – 1.00 and – 0.88 and for the incident ̄p laboratory momentum region between 0.70 and 2.37 GeV/c. These measurements, done in intervals of approximately 0.1 GeV/c, have been performed at the Alternating Gradient Synchrotron at Brookhaven National Laboratory during the winter of 1968. The measured differential cross sections, binned in cos θ_cm intervals of 0.02, have statistical errors of about 10%. Backward dipping exists below 0.95 GeV/c and backward peaking above 0.95 GeV/c. The 180˚ differential cross section extrapolated from our data shows a sharp dip centered at 0.95 GeV/c and a broad hump centered near 1.4 GeV/c. Our data have been interpreted in terms of resonance effects and in terms of diffraction dominance effects.

Nuclear reactions leading to energy levels of Ar37 and K37

In view of recent interest in the Cl³⁷ (ʋ solar’^e-)Ar³⁷ reaction cross section, information on some aspects of mass 37 nuclei has been obtained using the K³⁹ (d, ∝) Ar³⁷ and Cl³⁵ (He³, p) Ar³⁷ reactions. Ar³⁷ levels have been found at 0, 1.41, 1.62, 2.22, 2.50, 2.80, 3.17, 3.27, 3.53, 3.61, 3.71, (3.75), (3.90), 3.94, 4.02, (4.21), 4.28, 4.32, 4.40, 4.45, 4.58, 4.63, 4.74, 4.89, 4.98, 5.05, 5.10, 5.13, 5.21, 5.35, 5.41, 5.44, 5.54, 5.58, 5.67, 5.77, and 5.85 MeV (the underlined values correspond to previously tabulated levels). The nuclear temperature calculated from the Ar³⁷ level density is 1.4 MeV. Angular distributions of the lowest six levels with the K³⁹ (d, ∝) Ar³⁷ reaction at E_d = 10 MeV indicate a dominant direct interaction mechanism and the inapplicability of the 2I + 1 rule of the statistical model. Comparison of the spectra obtained with the K³⁹ (d, ∝) Ar³⁷ and Cl³⁵ (He³, p) Ar³⁷ reactions leads to the suggestion that the 5.13-MeV level is the T = 3/2 Cl³⁷ ground state analog. The ground state Q-value of the Ca⁴⁰ (p, ∝) K³⁷ reaction has been measured: -5179 ± 9 keV. This value implies a K³⁷ mass excess of -24804 ± 10 keV. Description of a NMR magnetometer and a sixteen-detector array used in conjunction with a 61-cm double-focusing magnetic spectrometer are included in appendices.

Inelastic K+p reactions at incident momenta from 1.37 GeV/c to 2.17 GeV/c

A number of recent experiments have suggested the possibility of a highly inelastic resonance in K⁺p scattering. To study the inelastic K⁺p reactions, a 400 K exposure has been taken at the L.R.L. 25 inch bubble chamber. The data are spread over seven K⁺ momenta between 1.37 and 2.17 GeV/c.

Cross-sections have been measured for the reaction K⁺p → pK°π+ which is dominated by the quasi-two body channels K∆ and K*N. Both these channels are strongly peripheral, as at other momenta. The decay of the ∆ is in good agreement with the predictions of the rho-photon analogy of Stodolsky and Sakurai. The data on the K*p channel show evidence of both pseudo scalar and vector exchange.

Cross-sections for the final state pK⁺π+π- shows a strong contribution from the quasi-two body channel K*∆. This reaction is also very peripheral even at threshold. The decay angular distributions indicate the reaction is dominated as at higher momenta by a pion exchange mechanism. The data are also in good agreement with the quark model predictions of Bialas and Zalewski for the K* and ∆ decay.

CP and the three pion decay of the K°

The time distribution of the decays of an initially pure K° beam into π+π-π° has been analyzed to determine the complex parameter W (also known as Ƞ^+-° and (x + iy)). The K° beam was produced in a brass target by the interactions of a 2.85 GeV/c π- beam which was generated on an internal target in the Lawrence Radiation Laboratory (LRL) Bevatron. The counters and hodoscopes in the apparatus selected for events with a neutral (K°) produced in the brass target, two charged secondaries passing through a magnet spectrometer and a ɣ-ray shower in a shower hodoscope.

From the 275K apparatus triggers, 148 K → π+π-π° events were isolated. The presence of a ɣ-ray shower in the optical shower chambers and a two-prong vee in the optical spark chambers were devices used to isolate the events. The backgrounds were further reduced by reconstructing the momenta of the two charged secondaries and applying kinematic constraints.

The best fit to the final sample of 148 events distributed between .3 and 7.0 K_S lifetimes gives:

ReW = -.05 ±.17

ImW = +.39 +.35/-.37

This result is consistent with both CPT invariance (ReW = 0) and CP invariance (W = 0). Backgrounds are estimated to be less than 10% and systematic effects have also been estimated to be negligible.

An analysis of the present data on CP violation in this decay mode and other K° decay modes has estimated the phase of ɛ to be 45.3 ± 2.3 degrees. This result is consistent with the super weak theories of CP violation which predicts the phase of ɛ to be 43°. This estimate is in turn used to predict the phase of Ƞ°° to be 48.0 ± 7.9 degrees. This is a substantial improvement on presently available measurements. The largest error in this analysis comes from the present limits on W from the world average of recent experiments. The K → πuʋ mode produces the next largest error. Therefore further experimentation in these modes would be useful.