Any symmetry reduces a second-order differential equation to a first integral: variational symmetries of the action (exemplified by central field dynamics) lead to conservation laws, but symmetries of only the equations of motion (exemplified by scale-invariant hydrostatics) yield first-order {\em non-conservation laws} between invariants. We obtain these non-conservation laws by extending Noether's Theorem to non-variational symmetries and present an innovative variational formulation of spherical adiabatic hydrostatics. For the scale-invariant case, this novel synthesis of group theory, hydrostatics, and astrophysics allows us to recover all the known properties of polytropes and define a {\em core radius}, inside which polytropes of index $n$ share a common core mass density structure, and outside of which their envelopes differ. The Emden solutions (regular solutions of the Lane-Emden equation) are obtained, along with useful approximations. An appendix discusses the $n=3$ polytrope in order to emphasize how the same mechanical structure allows different thermal structures in relativistic degenerate white dwarfs and zero age main sequence stars.

Any symmetry reduces a second-order differential equation to a first-order equation: variational symmetries of the action (exemplified by central field dynamics) lead to conservation laws, but symmetries of only the equations of motion (exemplified by scale-invariant hydrostatics), yield first-order {\em non-conservation laws} between invariants. We obtain these conservation laws by extending Noether's Theorem to non-variational symmetries, and present a variational formulation of spherical adiabatic hydrostatics. For scale-invariant hydrostatics, we directly recover all the published properties of polytropes and define a {\em core radius}, a new measure of mass concentration in polytropes of index n. The Emden solutions (regular solutions of the Lane-Emden equation) are finally obtained, along with useful approximations. An appendix discusses the special n=3 polytrope, emphasizing how the same mechanical structure allows different {\em thermostatic} structures in relativistic degenerate white dwarfs and and zero age main sequence stars.

Because scaling symmetries of the Euler-Lagrange equations are generally not variational symmetries of the action, they do not lead to conservation laws. Instead, an extension of Noether's theorem reduces the equations of motion to evolutionary laws that prove useful, even if the transformations are not generalized symmetries of the equations of motion. In the case of scaling, symmetry leads to a scaling evolutionary law, a first-order equation in terms of scale invariants, linearly relating kinematic and dynamic degrees of freedom. This scaling evolutionary law appears in dynamical and in static systems. Applied to dynamical central-force systems, the scaling evolutionary equation leads to generalized virial laws, which linearly connect the kinetic and potential energies. Applied to barotropic hydrostatic spheres, the scaling evolutionary equation linearly connects the gravitational and internal energy densities. This implies well-known properties of polytropes, describing degenerate stars and chemically homogeneous non-degenerate stellar cores.

In scale invariant hydrostatic barotropes, the radial evolutionary equation linearly relates the local gravitational and internal energies. From this first-order equation, directly follow all the properties of polytropes and the important mass-radius relation. Quadrature then leads to the regular Lane-Emden functions and their Picard and Pade approximations, which are useful wherever stars are approximately or exactly polytropic. We illustrate this particularly for the n=3 regular polytrope and obtain analytic approximations to the solution of the Lane-Emden equation, valid over the bulk of relativistic degenerate stars (massive white dwarfs) and chemically homogeneous stars in radiative equilibrium (ZAMS stars).

Scaling symmetries of the Euler-Lagrange equations are generally not variational symmetries of the action and do not lead to conservation laws. Nevertheless, by an extension of Noether's theorem, scaling symmetries lead to useful {\em nonconservation} laws, which still reduce the Euler-Lagrange equations to first order in terms of scale invariants. We illustrate scaling symmetry dynamically and statically. Applied dynamically to systems of bodies interacting via central forces, the nonconservation law is Lagrange's identity, leading to generalized virial laws. Applied to self-gravitating spheres in hydrostatic equilibrium, the nonconservation law leads to well-known properties of polytropes describing degenerate stars and chemically homogeneous nondegenerate stellar cores.

We review current observations of the homogeneous cosmological expansion which, because they measure only kinematic variables, cannot determine the dynamics driving the recent accelerated expansion. The minimal fit to the data, the flat $\Lambda CDM$ model, consisting of cold dark matter and a cosmological constant, interprets $4\Lambda$ geometrically as a classical spacetime curvature constant of nature, avoiding any reference to quantum vacuum energy. (The observed Uehling and Casimir effects measure forces due to QED vacuum polarization, but not any quantum material vacuum energies.) An Extended Anthropic Principle, that Dark Energy and Dark Gravity be indistinguishable, selects out flat $\Lambda CDM$. Prospective cosmic shear and galaxy clustering observations of the growth of fluctuations are intended to test whether the 'dark energy' driving the recent cosmological acceleration is static or moderately dynamic. Even if dynamic, observational differences between an additional negative-pressure material component within general relativity (Dark Energy) and low-curvature modifications of general relativity (Dark Gravity) will be extremely small.

The homogeneous expansion history H(z) of our universe measures only kinematic variables, but cannot fix the underlying dynamics driving the recent acceleration: cosmographic measurements of the homogeneous universe, are consistent with either a static finely-tuned cosmological constant or a dynamic `dark energy' mechanism, which may be material Dark Energy or modified gravity (Dark Gravity). Resolving the composition and foreground noise to reduce their large systematic errors.dynamics of either kind of `dark energy', will require complementing the homogeneous expansion observations with observations of the growth of cosmological fluctuations. Because the 'dark energy' evolution is at least quasi-static, any dynamical effects on the fluctuation growth function g(z) will be minimal. They will be best studied in the weak lensing convergence of light from galaxies at 0<z<5, from neutral hydrogen at 6<z<20, and ultimately from the CMB last scattering surface at z=1089. Galaxy clustering also measures g(z), but requires large corrections for baryonic composition and foreground noise. Projected observations potentially distinguish static from dynamic `dark energy', but distinguishing dynamic Dark Energy from Dark Gravity will require a weak lensing shear survey more ambitious than any now projected. Low-curvature modifications of Einstein gravity are also, in principle, observable in the solar system or in isolated galaxy clusters. The cosmological constant can only be fine-tuned, at present. The Cosmological Coincidence Problem, that we live when the ordinary matter density approximates the 'gravitational vacuum energy', on the other hand, is a material problem, calling for an understanding of the observers' role in cosmology.

We review the evidence for recently accelerating cosmological expansion or "dark energy", either a negative pressure constituent in General Relativity (Dark Energy) or modified gravity (Dark Gravity), without any Dark Energy constituent. If constituent Dark Energy does not exist, so that our universe is now dominated by pressure-free matter, Einstein gravity must be modified at low curvature. The vacuum symmetry of any Robertson-Walker universe then characterizes Dark Gravity as low- or high-curvature modifications of Einstein gravity. The dynamics of either kind of "dark energy" cannot be derived from the homogeneous expansion alone, but requires also observing the growth of inhomogeneities. Present and projected observations are all consistent with a small fine tuned cosmological constant, with nearly static Dark Energy, or with gravity modified at cosmological scales. The growth of cosmological fluctuations will potentially distinguish static "dark energy" from dynamic "dark energy" with equation of state $w(z)$ either changing rapidly or tracking the background matter. But to cosmologically distinguish $\Lambda$CDM from modified gravity will require a weak lensing shear survey more ambitious than any now projected. Dvali-Gabadadze-Porrati modifications of Einstein gravity may also be detected in refined bservations in the solar system or at the intermediate Vainstein scale. Dark Energy's epicyclic character, failure to explain the original Cosmic Coincidence ("Why now?") without fine tuning, inaccessibility to laboratory or solar system tests, along with braneworld theories, now motivate future precision solar system, Vainstein-scale and cosmological-scale studies of Dark Gravity.

Good tracking requires that the quintessence energy fraction slowly increase while the roll $\lambda\equiv -d\ln V/\varkappa d\phi$ slowly decreases, but is not yet truly slow-rolling. The supernova bound on the present quintessence equation of state requires either (1) a cosmological constant or other fine-tuned "crawling quintessence" or (2) "roll-over quintessence" that tracked until recently, but now became slow rolling, because of a sharp increase in potential curvature. Thus, fine-tuning is required by constant equation of state and inverse power potentials, but can be avoided by the SUGRA and Skordis-Albrecht potentials and other good trackers, provided quintessence energy domination and slow roll {\em both} began only recently. This makes the time in which we live special in {\em two} respects.

Good tracking requires significant quintessence energy fraction, even in the past, but a potential that is not yet truly slow-rolling. The supernova bound on cosmic acceleration excludes constant equation of state and inverse power potentials, but allows the SUGRA potential and other good trackers, in which quintessence energy domination and kinetic energy suppression {\em both} began only recently. This makes the time in which we live special in {\em two} respects.

Within present constraints on the observed smooth energy and its equation of state parameter, it is important to find out whether the smooth energy is static (cosmological constant) or dynamic (quintessence). The most dynamical quintessence fields observationally allowed are now still fast-rolling and no longer satisfy the tracker approximation if the equation of state parameter varies moderately with cosmic scale. We are optimistic about distinguishing between a cosmological constant and appreciably dynamic quintessence, by measuring average values for the effective equation of state parameter. However, reconstructing the quintessence potential from observations of any scale dependence appears problematic in the near future. For our flat universe, at present dominated by smooth energy in the form of either a cosmological constant (LCDM) or quintessence (QCDM), we calculate the asymptotic collapsed mass fraction to be maximal at the observed smooth energy/matter ratio. Identifying this collapsed fraction as a conditional probability for habitable galaxies, we infer that the prior distribution is flat. Interpreting this prior as a distribution over theories, rather than as a distribution over unobservable subuniverses, leads us to heuristic predictions about the class of future quantum cosmology theories and the static or quasi-static nature of the smooth energy.

This paper proposes two phenomenological quintessence potentials with parameters fitted to the presently observed ratio $\rat$ of smooth energy to clustered mass and limits on the equation of state parameter $w_Q(a)=P/\rhoq$, for which the quintessence potential is now slow- or fast-rolling. These two potentials are intended to illustrate rather different quintessence phenomenology since tracking and into the near future. Neither can be a fundamental potential derivable from present-day string theory.

The four equations of stellar structure are reformulated as two alternate pairs of variational principles. Different thermodynamic representations lead to the same hydromechanical equations, but the thermal equations require, not the entropy, but the temperature as the thermal field variable. Our treatment emphasizes the hydrostatic energy and the entropy production rate of luminosity produced and transported. The conceptual and calculational advantages of integral over differential formulations of stellar structure are discussed along with the difficulties in describing stellar chemical evolution by variational principles.

All stars exhibit universal central behavior in terms of new homology variables (u,w). In terms of these variables, we obtain simple analytic fits to numerical standard solar models for the core and radiative zones of the ZAMS and present Suns, with a few global parameters. With these analytic fits, different theoretical models of the solar core, neutrino fluxes, and helioseismic observations can be parametrized and compared.

Precise numerical standard solar models (SSMs) now agree with one another and with helioseismological observations in the convective and outer radiative zones. Nevertheless these models obscure how luminosity, neutrino production and g-mode core helioseismology depend on such inputs as opacity and nuclear cross sections. Although the Sun is not homologous, its inner core by itself is chemically evolved and almost homologous, because of its compactness, radiative energy transport, and ppI-dominated luminosity production. We apply luminosity-fixed homology transformations to the core to estimate theoretical uncertainties in the SSM and to obtain a broad class of non-SSMs, parametrized by central temperature and density and purely radiative energy transport in the core.

For a flat universe presently dominated by smooth energy, either cosmological constant (LCDM) or quintessence (QCDM), we calculate the asymptotic collapsed mass fraction as function of the present ratio of smooth energy to matter energy $\mathcal R_0$. Identifying the normalized collapsed fraction as a conditional probability for habitable galaxies, we observe that the observed present ratio $\mathcal R_0 \sim 2$ is likely in LCDM, but more likely in QCDM. Inverse application of Bayes' Theorem makes the Anthropic Principle a predictive scientific principle: the data implies that the prior probability for $\mathcal R_0$ must be essentially flat over the anthropically allowed range. Interpreting this prior as a distribution over {\em theories} lets us predict that any future theory of initial conditions must be indifferent to $\mathcal R_0$. This application of the Anthropic Principle does not demand the existence of other universes.

For a universe presently dominated by static or dynamic vacuum energy, cosmological constant (LCDM) or quintessence (QCDM), we calculate the asymptotic collapsed mass fraction as function of the present ratio of vacuum energy to clustered mass, \Omega_{Q0}/\Omega_{m0}. Identifying these collapsed fractions as anthropic probabilities, we find the value \Omega_{Q0}/\Omega_{m0} \sim 2 to be reasonably likely in LCDM, and very likely in QCDM.

We review the cosmological evidence for a low matter density universe and a cosmological constant or dynamical vacuum energy and address the cosmolog$ coincidence problem: why is the matter density about one-half the vacuum energy {\em now}. This is reasonble, following the anthropic argument of Efstathiou and of Martel, Schapiro & Weinberg.

Direct observations of the supposedly universal primordial deuterium abundance imply a relatively large baryon density $\Omega_B= (0.019-0.030)h^{-2}$ (95% C.L.). On the other hand, concordance between the previously accepted $^4 He$ and $^7 Li$ abundances and standard Big Bang Nucleosynthesis requires the thrice smaller value $\Omega_B=(0.005-0.010)h^{-2}$ (95% C.L.). For each $\Omega_B$, we use X-ray and Sunyaev-Zeldovich observations of the baryon fraction $f_B$ in rich clusters of galaxies, in order to obtain limits on the total mass density $\Omega_{cl}$ in clusters of galaxies. The higher-$\Omega_B$ values are consistent with clusters being a fair sample of the Universe, and then imply $\Omega_{m} =(0.3-0.9)$, a medium or critical density Universe. Said otherwise, the observed limits $f_B > 0.1, ~ \Omega_m > 0.3$ imply $\Omega_B > 0.03$. If the newer$^4 He$ abundance observations are accepted, this is consistent with standard BBN.

Constraints on big bang nucleosynthesis (BBN) and on cosmological parameters from conflicting deuterium observations in different high red-shift QSO systems are discussed. The high deuterium observations by Carswell et al., Songaila et al., and Rugers & Hogan is consistent with 4He and 7Li observations and Standard BBN (N_nu =3) and allows N_nu \leq 3.6 at 95% C.L., but is inconsistent with local observations of D and 3He in the context of conventional theories of stellar and Galactic evolution. In contrast, the low deuterium observations by Tytler, Fan & Burles and Burles & Tytler are consistent with the constraints from local Galactic observations, but require N_nu = 1.9 +/- 0.3 at 68% C.L., excluding Standard BBN at 99.9% C.L., unless the systematic uncertainties in the 4He observations have been underestimated by a large amount. The high and low primordial deuterium abundances imply, respectively, Omega_B h^2 = 0.005 - 0.01 and Omega_B h^2 = 0.02 - 0.03 at 95% C.L. When combined with the high baryon fraction inferred from x-ray observations of rich clusters, the corresponding total mass densities (for 50 \le H_0 \le 90) are Omega_M = 0.05 - 0.20 and Omega_M = 0.2 - 0.7, respectively (95\% C.L.) The range of Omega_M corresponding to high D is in conflict with dynamical constraints (Omega_M > 0.2 - 0.3) and with the shape parameter constraint (Gamma = Omega_M h = 0.25 +/- 0.05) from large scale structure formation in CDM and Lambda-CDM models.

A new evaluation of the constraint on the number of light neutrino species (N_nu) from big bang nucleosynthesis suggests a discrepancy between the predicted light element abundances and those inferred from observations, unless the inferred primordial 4He abundance has been underestimated by 0.014 +/- 0.004 (1 sigma) or less than 10% (95%C.L.) of 3He survives stellar processing. With the quoted systematic errors in the observed abundances and a conservative chemical evolution parameterization, the best fit to the combined data is N_nu = 2.1 +/- 0.3 (1 sigma) and the upper limit is N_nu < 2.6 (95% C.L.). The data are inconsistent with the Standard Model (N_nu = 3) at the 98.6% C.L.

We consider the most general solar model, using the neutrino fluxes as free parameters constrained only by the solar luminosity, and show that the combined solar neutrino data exclude any astrophysical solution at 98\% C.L.\ Our best fit to the $^7$Be and $^8$B fluxes is respectively $<$7\% and 37$\pm$4\% of the standard solar model prediction, but only with a large $\chi^2$ (5.6 for 1 d.f.). This best fit to the fluxes contradicts explicit nonstandard solar models, which generally reduce the $^8$B flux more than the $^7$Be. Those models are well parameterized by a single parameter, the central temperature.

Constraints on the core temperature (T_c) of the Sun and on neutrino- oscillation parameters are obtained from the existing solar neutrino data, including the recent GALLEX and Kamiokande III results. (1) A purely astrophysical solution to the solar neutrino problem is strongly disfavored by the data: the best fit in a cooler Sun model requires an 8% reduction in T_c, but the chi-sqaured test rejects this hypothesis at 99.99% C.L., suggesting new neutrino physics. (2) Assuming the Standard Solar Model (SSM) and MSW oscillations, the MSW parameters are constrained to two small regions: one in the non-adiabatic region and the other in the large-mixing region. The non-adiabatic solution gives a considerably better fit. For nu_e oscillations into sterile neutrinos, the allowed region (90%) is constrained to non- adiabatic oscillations. As long as the SSM is assumed, the neutrino mixing angles are at least four times larger, or considerably smaller, than the corresponding quark mixing angles. (3) Allowing both MSW oscillations and a non-standard core temperature, a) the experiments determine the core temperature at the 5% level, yielding a value consistent with the SSM prediction. b) When T_c is used as a free parameter, the allowed MSW region is broadened: a cooler Sun (T_c=0.95) allows mass and mixing implied by the SUSY SO(10) GUT, while a warmer Sun (T_c=1.05) allows parameter space suggested by intermediate-scale SO(10) GUTs. Superstring-inspired models are consistent with all solutions. (4) From the narrowed parameter space, we predict the neutrino spectral shape which should be observed in SNO. Throughout the calculation we use the latest Bahcall-Pinsonneault SSM, and include nuclear and astrophysical uncertainties in a simplified, but physically transparent way.