Self-Similar Polytropic Champagne Flows in H II Regions

We explore large-scale hydrodynamics of H II regions for various self-similar shock flows of a polytropic gas cloud under self-gravity and with quasi-spherical symmetry. We formulate cloud dynamics by invoking specific entropy conservation along streamlines and obtain global self-similar "champagne flows" for a conventional polytropic gas with shocks as a subclass. Molecular cloud cores are ionized and heated to high temperatures after the onset of nuclear burning of a central protostar. We model subsequent evolutionary processes in several ways and construct possible self-similar shock flow solutions. We may neglect the mass and gravity of the central protostar. The ionization and heating of the surrounding medium drive outflows in the inner cloud core and a shock travels outwards, leading to the so-called "champagne phase" with an expanding outer cloud envelope. Complementarily, we also consider the expansion of a central cavity around the centre. As the inner cloud expands plausibly due to powerful stellar winds, a cavity (i.e., `void' or `bubble') can be created around the centre, and when the cavity becomes sufficiently large, one may neglect the gravity of the central protostar. We thus present self-similar shock solutions for "champagne flows" with an expanding central void. We compare our solutions with isothermal solutions and find that the generalization to the polytropic regime brings about significant differences of the gas dynamics, especially for cases of n<1, where n is a key scaling index in the self-similar transformation. We also compare our global polytropic self-similar solutions with numerical simulations on the expansion of H II regions.Comment: 17 pages, 10 figures, accepted for publication by MNRA

General Polytropic Magnetofluid under Self-Gravity: Voids and Shocks

We study the self-similar magnetohydrodynamics (MHD) of a quasi-spherical expanding void (viz. cavity or bubble) in the centre of a self-gravitating gas sphere with a general polytropic equation of state. We show various analytic asymptotic solutions near the void boundary in different parameter regimes and obtain the corresponding void solutions by extensive numerical explorations. We find novel void solutions of zero density on the void boundary. These new void solutions exist only in a general polytropic gas and feature shell-type density profiles. These void solutions, if not encountering the magnetosonic critical curve (MCC), generally approach the asymptotic expansion solution far from the central void with a velocity proportional to radial distance. We identify and examine free-expansion solutions, Einstein-de Sitter expansion solutions, and thermal-expansion solutions in three different parameter regimes. Under certain conditions, void solutions may cross the MCC either smoothly or by MHD shocks, and then merge into asymptotic solutions with finite velocity and density far from the centre. Our general polytropic MHD void solutions provide physical insight for void evolution, and may have astrophysical applications such as massive star collapses and explosions, shell-type supernova remnants and hot bubbles in the interstellar and intergalactic media, and planetary nebulae.Comment: 21 pages, 15 figures, accepted for publication on New Astronom

Robust control strategies for multi-inventory systems with average flow constraints.

In this paper, we consider multi-inventory systems in the presence of uncertain demand. We assume that (i) demand is unknown but bounded in an assigned compact set and (ii) the control inputs (controlled flows) are subject to assigned constraints. Given a long-term average demand, we select a nominal flow that feeds such a demand. In this context, we are interested in a control strategy that meets at each time all possible current demands and achieves the nominal flow in the average. We provide necessary and sufficient conditions for such a strategy to exist and we characterize the set of achievable flows. Such conditions are based on linear programming and thus they are constructive. In the special case of a static flow (i.e. a system with 0-capacity buffers) we show that the strategy must be affine. The dynamic problem can be solved by a linear-saturated control strategy (inspired by the previous one). We provide numerical analysis and illustrative examples

Measurement of CP Asymmetries and Branching Fractions in Charmless Two-Body B-Meson Decays to Pions and Kaons

We present improved measurements of CP-violation parameters in the decays $B^0 \to \pi^+ \pi^-$, $B^0 \to K^+ \pi^-$, and $B^0 \to \pi^0 \pi^0$, and of the branching fractions for $B^0 \to \pi^0 \pi^0$ and $B^0 \to K^0 \pi^0$. The results are obtained with the full data set collected at the $\Upsilon(4S)$ resonance by the BABAR experiment at the PEP-II asymmetric-energy $B$ factory at the SLAC National Accelerator Laboratory, corresponding to $467 \pm 5$ million $B\bar B$ pairs. We find the CP-violation parameter values and branching fractions $S_{\pi^+\pi^-} = -0.68 \pm 0.10 \pm 0.03, C_{\pi^+\pi^-} = -0.25 \pm 0.08 \pm 0.02, A_{K^-\pi^+} = -0.107 \pm 0.016 ^{+0.006}_{-0.004}, C_{\pi^0\pi^0} = -0.43 \pm 0.26 \pm 0.05, Br(B^0 \to \pi^0 \pi^0) = (1.83 \pm 0.21 \pm 0.13) \times 10^{-6}, Br(B^0 \to K^0 \pi^0) = (10.1 \pm 0.6 \pm 0.4) \times 10^{-6},$ where in each case, the first uncertainties are statistical and the second are systematic. We observe CP violation with a significance of 6.7 standard deviations for $B^0 \to\pi^+\pi^-$ and 6.1 standard deviations for $B^0 \to K^+ \pi^-$, including systematic uncertainties. Constraints on the Unitarity Triangle angle $\alpha$ are determined from the isospin relations among the $B \to \pi\pi$ rates and asymmetries. Considering only the solution preferred by the Standard Model, we find $\alpha$ to be in the range $[71^\circ,109^\circ]$ at the 68% confidence level.Comment: 18 pages, 11 postscript figures, submitted to Phys. Rev.

Measurement of the Time-Dependent CP Asymmetry of Partially Reconstructed B0->D*+D*- Decays

We present a new measurement of the time-dependent CP asymmetry of B0->D*+D*- decays using (471+-5) million BBbar pairs collected with the BaBar detector at the PEP-II B Factory at the SLAC National Accelerator Laboratory. Using the technique of partial reconstruction, we measure the time-dependent CP asymmetry parameters S=-0.34+-0.12+-0.05$ and C=+0.15+-0.09+-0.04. Using the value for the CP-odd fraction R_perp=0.158+-0.028+-0.006, previously measured by BaBar with fully reconstructed B0->D*+D*- events, we extract the CP-even components S+=-0.49+-0.18+-0.07+-0.04 and C+=+0.15+-0.09+-0.04. In each case, the first uncertainty is statistical and the second is systematic; the third uncertainty on S+ is the contribution from the uncertainty on R_perp. The measured value of the CP-even component S+ is consistent with the value of sin(2Beta) measured in b->(ccbar)s transitions, and with the Standard Model expectation of small penguin contributions.Comment: 17 pages, 7 figures, submitted to Physical Review

Search for lepton-number violating processes in B+ -> h- l+ l+ decays

We have searched for the lepton-number violating processes B+ -> h- l+ l+ with h- = K-/pi- and l+ = e+/mu+, using a sample of 471+/-3 million BBbar events collected with the BaBar detector at the PEP-II e+e- collider at the SLAC National Accelerator Laboratory. We find no evidence for these decays and place 90% confidence level upper limits on their branching fractions Br(B+ -> pi- e+ e+) K- e+ e+) pi- mu+ mu+) K- mu+ mu+) < 6.7 x 10^{-8}.Comment: 8 pages, 4 postscript figures, submitted to Phys. Rev. D. R

Measurement of B(B-->X_s {\gamma}), the B-->X_s {\gamma} photon energy spectrum, and the direct CP asymmetry in B-->X_{s+d} {\gamma} decays

The photon spectrum in B --> X_s {\gamma} decay, where X_s is any strange hadronic state, is studied using a data sample of (382.8\pm 4.2) \times 10^6 e^+ e^- --> \Upsilon(4S) --> BBbar events collected by the BABAR experiment at the PEP-II collider. The spectrum is used to measure the branching fraction B(B --> X_s \gamma) = (3.21 \pm 0.15 \pm 0.29 \pm 0.08)\times 10^{-4} and the first, second, and third moments = 2.267 \pm 0.019 \pm 0.032 \pm 0.003 GeV,, )^2> = 0.0484 \pm 0.0053 \pm 0.0077 \pm 0.0005 GeV^2, and )^3> = -0.0048 \pm 0.0011 \pm 0.0011 \pm 0.0004 GeV^3, for the range E_\gamma > 1.8 GeV, where E_{\gamma} is the photon energy in the B-meson rest frame. Results are also presented for narrower E_{\gamma} ranges. In addition, the direct CP asymmetry A_{CP}(B --> X_{s+d} \gamma) is measured to be 0.057 \pm 0.063. The spectrum itself is also unfolded to the B-meson rest frame; that is the frame in which theoretical predictions for its shape are made.Comment: 37 pages, 19 postscript figures, submitted to Phys. Rev. D. No analysis or results have changed from previous version. Some changes to improve clarity based on interactions with Phys. Rev. D referees, including one new Figure (Fig. 13), and some minor wording/punctuation/spelling mistakes fixe

Improved Limits on B0B^{0} decays to invisible (+γ)(+\gamma) final states

We establish improved upper limits on branching fractions for B0 decays to final States 10 where the decay products are purely invisible (i.e., no observable final state particles) and for final states where the only visible product is a photon. Within the Standard Model, these decays have branching fractions that are below the current experimental sensitivity, but various models of physics beyond the Standard Model predict significant contributions for these channels. Using 471 million BB pairs collected at the Y(4S) resonance by the BABAR experiment at the PEP-II e+e- storage ring at the SLAC National Accelerator Laboratory, we establish upper limits at the 90% confidence level of 2.4x10^-5 for the branching fraction of B0-->Invisible and 1.7x10^-5 for the branching fraction of B0-->Invisible+gammaComment: 8 pages, 3 postscript figures, submitted to Phys. Rev. D (Rapid Communications

Precise Measurement of the e+ e- --> pi+ pi- (gamma) Cross Section with the Initial-State Radiation Method at BABAR

A precise measurement of the cross section of the process $e^+e^-\to\pi^+\pi^-(\gamma)$ from threshold to an energy of 3GeV is obtained with the initial-state radiation (ISR) method using 232fb$^{-1}$ of data collected with the BaBar detector at $e^+e^-$ center-of-mass energies near 10.6GeV. The ISR luminosity is determined from a study of the leptonic process $e^+e^-\to\mu^+\mu^-(\gamma)\gamma_{\rm ISR}$, which is found to agree with the next-to-leading-order QED prediction to within 1.1%. The cross section for the process $e^+e^-\to\pi^+\pi^-(\gamma)$ is obtained with a systematic uncertainty of 0.5% in the dominant $\rho$ resonance region. The leading-order hadronic contribution to the muon magnetic anomaly calculated using the measured $\pi\pi$ cross section from threshold to 1.8GeV is $(514.1 \pm 2.2({\rm stat}) \pm 3.1({\rm syst}))\times 10^{-10}$.Comment: 58 pages, 56 figures, to be submitted to Phys. Rev.

Measurement of Branching Fractions and Rate Asymmetries in the Rare Decays B -> K(*) l+ l-

In a sample of 471 million BB events collected with the BABAR detector at the PEP-II e+e- collider we study the rare decays B -> K(*) l+ l-, where l+ l- is either e+e- or mu+mu-. We report results on partial branching fractions and isospin asymmetries in seven bins of di-lepton mass-squared. We further present CP and lepton-flavor asymmetries for di-lepton masses below and above the J/psi resonance. We find no evidence for CP or lepton-flavor violation. The partial branching fractions and isospin asymmetries are consistent with the Standard Model predictions and with results from other experiments.Comment: 16 pages, 14 figures, accepted by Phys. Rev.