Density Matrix Functional Calculations for Matter in Strong Magnetic Fields: I. Atomic Properties

We report on a numerical study of the density matrix functional introduced by Lieb, Solovej and Yngvason for the investigation of heavy atoms in high magnetic fields. This functional describes {\em exactly} the quantum mechanical ground state of atoms and ions in the limit when the nuclear charge $Z$ and the electron number $N$ tend to infinity with $N/Z$ fixed, and the magnetic field $B$ tends to infinity in such a way that $B/Z^{4/3}\to\infty$. We have calculated electronic density profiles and ground state energies for values of the parameters that prevail on neutron star surfaces and compared them with results obtained by other methods. For iron at $B=10^{12}$ G the ground state energy differs by less than 2 \% from the Hartree-Fock value. We have also studied the maximal negative ionization of heavy atoms in this model at various field strengths. In contrast to Thomas-Fermi type theories atoms can bind excess negative charge in the density matrix model. For iron at $B=10^{12}$ G the maximal excess charge in this model corresponds to about one electron.Comment: Revtex, 13 pages with 6 eps figures include

Structural Invariance of Sunspot Umbrae Over the Solar Cycle: 1993-2004

Measurements of maximum magnetic flux, minimum intensity, and size are presented for 12 967 sunspot umbrae detected on the NASA/NSO spectromagnetograms between 1993 and 2004 to study umbral structure and strength during the solar cycle. The umbrae are selected using an automated thresholding technique. Measured umbral intensities are first corrected for a confirming observation of umbral limb-darkening. Log-normal fits to the observed size distribution confirm that the size spectrum shape does not vary with time. The intensity-magnetic flux relationship is found to be steady over the solar cycle. The dependence of umbral size on the magnetic flux and minimum intensity are also independent of cycle phase and give linear and quadratic relations, respectively. While the large sample size does show a low amplitude oscillation in the mean minimum intensity and maximum magnetic flux correlated with the solar cycle, this can be explained in terms of variations in the mean umbral size. These size variations, however, are small and do not substantiate a meaningful change in the size spectrum of the umbrae generated by the Sun. Thus, in contrast to previous reports, the observations suggest the equilibrium structure, as testified by the invariant size-magnetic field relationship, as well as the mean size (i.e. strength) of sunspot umbrae do not significantly depend on solar cycle phase.Comment: 17 pages, 6 figures. Published in Solar Physic

Baxterization, dynamical systems, and the symmetries of integrability

We resolve the `baxterization' problem with the help of the automorphism group of the Yang-Baxter (resp. star-triangle, tetrahedron, \dots) equations. This infinite group of symmetries is realized as a non-linear (birational) Coxeter group acting on matrices, and exists as such, {\em beyond the narrow context of strict integrability}. It yields among other things an unexpected elliptic parametrization of the non-integrable sixteen-vertex model. It provides us with a class of discrete dynamical systems, and we address some related problems, such as characterizing the complexity of iterations.Comment: 25 pages, Latex file (epsf style). WARNING: Postscript figures are BIG (600kB compressed, 4.3MB uncompressed). If necessary request hardcopy to [email protected] and give your postal mail addres

In Situ Measurements of Interstellar Dust

We present the mass distribution of interstellar grains measured in situ by the Galileo and Ulysses spaceprobes as cumulative flux. The derived in situ mass distribution per logarithmic size interval is compared to the distribution determined by fitting extinction measurements. Large grains measured in situ contribute significantly to the overall mass of dust in the local interstellar cloud. The problem of a dust-to-gas mass ratio that contradicts cosmic abundances is discussed.Comment: 4 pages and two figure

On the Potts model partition function in an external field

We study the partition function of Potts model in an external (magnetic) field, and its connections with the zero-field Potts model partition function. Using a deletion-contraction formulation for the partition function Z for this model, we show that it can be expanded in terms of the zero-field partition function. We also show that Z can be written as a sum over the spanning trees, and the spanning forests, of a graph G. Our results extend to Z the well-known spanning tree expansion for the zero-field partition function that arises though its connections with the Tutte polynomial

Thermal Density Functional Theory in Context

This chapter introduces thermal density functional theory, starting from the ground-state theory and assuming a background in quantum mechanics and statistical mechanics. We review the foundations of density functional theory (DFT) by illustrating some of its key reformulations. The basics of DFT for thermal ensembles are explained in this context, as are tools useful for analysis and development of approximations. We close by discussing some key ideas relating thermal DFT and the ground state. This review emphasizes thermal DFT's strengths as a consistent and general framework.Comment: Submitted to Spring Verlag as chapter in "Computational Challenges in Warm Dense Matter", F. Graziani et al. ed

A mathematical framework for critical transitions: normal forms, variance and applications

Critical transitions occur in a wide variety of applications including mathematical biology, climate change, human physiology and economics. Therefore it is highly desirable to find early-warning signs. We show that it is possible to classify critical transitions by using bifurcation theory and normal forms in the singular limit. Based on this elementary classification, we analyze stochastic fluctuations and calculate scaling laws of the variance of stochastic sample paths near critical transitions for fast subsystem bifurcations up to codimension two. The theory is applied to several models: the Stommel-Cessi box model for the thermohaline circulation from geoscience, an epidemic-spreading model on an adaptive network, an activator-inhibitor switch from systems biology, a predator-prey system from ecology and to the Euler buckling problem from classical mechanics. For the Stommel-Cessi model we compare different detrending techniques to calculate early-warning signs. In the epidemics model we show that link densities could be better variables for prediction than population densities. The activator-inhibitor switch demonstrates effects in three time-scale systems and points out that excitable cells and molecular units have information for subthreshold prediction. In the predator-prey model explosive population growth near a codimension two bifurcation is investigated and we show that early-warnings from normal forms can be misleading in this context. In the biomechanical model we demonstrate that early-warning signs for buckling depend crucially on the control strategy near the instability which illustrates the effect of multiplicative noise.Comment: minor corrections to previous versio

Study of the q^2-Dependence of B --> pi ell nu and B --> rho(omega)ell nu Decay and Extraction of |V_ub|

We report on determinations of |Vub| resulting from studies of the branching fraction and q^2 distributions in exclusive semileptonic B decays that proceed via the b->u transition. Our data set consists of the 9.7x10^6 BBbar meson pairs collected at the Y(4S) resonance with the CLEO II detector. We measure B(B0 -> pi- l+ nu) = (1.33 +- 0.18 +- 0.11 +- 0.01 +- 0.07)x10^{-4} and B(B0 -> rho- l+ nu) = (2.17 +- 0.34 +0.47/-0.54 +- 0.41 +- 0.01)x10^{-4}, where the errors are statistical, experimental systematic, systematic due to residual form-factor uncertainties in the signal, and systematic due to residual form-factor uncertainties in the cross-feed modes, respectively. We also find B(B+ -> eta l+ nu) = (0.84 +- 0.31 +- 0.16 +- 0.09)x10^{-4}, consistent with what is expected from the B -> pi l nu mode and quark model symmetries. We extract |Vub| using Light-Cone Sum Rules (LCSR) for 0<= q^2<16 GeV^2 and Lattice QCD (LQCD) for 16 GeV^2 <= q^2 < q^2_max. Combining both intervals yields |Vub| = (3.24 +- 0.22 +- 0.13 +0.55/-0.39 +- 0.09)x10^{-3}$ for pi l nu, and |Vub| = (3.00 +- 0.21 +0.29/-0.35 +0.49/-0.38 +-0.28)x10^{-3} for rho l nu, where the errors are statistical, experimental systematic, theoretical, and signal form-factor shape, respectively. Our combined value from both decay modes is |Vub| = (3.17 +- 0.17 +0.16/-0.17 +0.53/-0.39 +-0.03)x10^{-3}.Comment: 45 pages postscript, also available through http://w4.lns.cornell.edu/public/CLNS, submitted to PR

Measurement of the Charge Asymmetry in B→K∗(892)±π∓B\to K^* (892)^{\pm}\pi^{\mp}

We report on a search for a CP-violating asymmetry in the charmless hadronic decay B -> K*(892)+- pi-+, using 9.12 fb^-1 of integrated luminosity produced at \sqrt{s}=10.58 GeV and collected with the CLEO detector. We find A_{CP}(B -> K*(892)+- pi-+) = 0.26+0.33-0.34(stat.)+0.10-0.08(syst.), giving an allowed interval of [-0.31,0.78] at the 90% confidence level.Comment: 7 pages postscript, also available through http://w4.lns.cornell.edu/public/CLNS, submitted to PR

Observation of a Narrow Resonance of Mass 2.46 GeV/c^2 Decaying to D_s^*+ pi^0 and Confirmation of the D_sJ^* (2317) State

Using 13.5 inverse fb of e+e- annihilation data collected with the CLEO II detector we have observed a narrow resonance in the Ds*+pi0 final state, with a mass near 2.46 GeV. The search for such a state was motivated by the recent discovery by the BaBar Collaboration of a narrow state at 2.32 GeV, the DsJ*(2317)+ that decays to Ds+pi0. Reconstructing the Ds+pi0 and Ds*+pi0 final states in CLEO data, we observe peaks in both of the corresponding reconstructed mass difference distributions, dM(Dspi0)=M(Dspi0)-M(Ds) and dM(Ds*pi0)=M(Ds*pi0)-M(Ds*), both of them at values near 350 MeV. We interpret these peaks as signatures of two distinct states, the DsJ*(2317)+ plus a new state, designated as the DsJ(2463)+. Because of the similar dM values, each of these states represents a source of background for the other if photons are lost, ignored or added. A quantitative accounting of these reflections confirms that both states exist. We have measured the mean mass differences = 350.0 +/- 1.2 [stat] +/- 1.0 [syst] MeV for the DsJ*(2317) state, and = 351.2 +/- 1.7 [stat] +/- 1.0 [syst] MeV for the new DsJ(2463)+ state. We have also searched, but find no evidence, for decays of the two states via the channels Ds*+gamma, Ds+gamma, and Ds+pi+pi-. The observations of the two states at 2.32 and 2.46 GeV, in the Ds+pi0 and Ds*+pi0 decay channels respectively, are consistent with their interpretations as (c anti-strange) mesons with orbital angular momentum L=1, and spin-parities of 0+ and 1+.Comment: 16 pages postscript, also available through http://w4.lns.cornell.edu/public/CLNS, version to be published in Physical Review D; minor modifications and fixes to typographical errors, plus an added section on production properties. The main results are unchanged; they supersede those reported in hep-ex/030501