21,206 research outputs found

    Geometric approach to Fletcher's ideal penalty function

    Get PDF
    Original article can be found at: www.springerlink.com Copyright Springer. [Originally produced as UH Technical Report 280, 1993]In this note, we derive a geometric formulation of an ideal penalty function for equality constrained problems. This differentiable penalty function requires no parameter estimation or adjustment, has numerical conditioning similar to that of the target function from which it is constructed, and also has the desirable property that the strict second-order constrained minima of the target function are precisely those strict second-order unconstrained minima of the penalty function which satisfy the constraints. Such a penalty function can be used to establish termination properties for algorithms which avoid ill-conditioned steps. Numerical values for the penalty function and its derivatives can be calculated efficiently using automatic differentiation techniques.Peer reviewe

    Magnetic fields and spiral arms in the galaxy M51

    Full text link
    (Abridged) We use new multi-wavelength radio observations, made with the VLA and Effelsberg telescopes, to study the magnetic field of the nearby galaxy M51 on scales from 200\pc to several \kpc. Interferometric and single dish data are combined to obtain new maps at \wwav{3}{6} in total and polarized emission, and earlier \wav{20} data are re-reduced. We compare the spatial distribution of the radio emission with observations of the neutral gas, derive radio spectral index and Faraday depolarization maps, and model the large-scale variation in Faraday rotation in order to deduce the structure of the regular magnetic field. We find that the \wav{20} emission from the disc is severely depolarized and that a dominating fraction of the observed polarized emission at \wav{6} must be due to anisotropic small-scale magnetic fields. Taking this into account, we derive two components for the regular magnetic field in this galaxy: the disc is dominated by a combination of azimuthal modes, m=0+2m=0+2, but in the halo only an m=1m=1 mode is required to fit the observations. We disuss how the observed arm-interarm contrast in radio intensities can be reconciled with evidence for strong gas compression in the spiral shocks. The average arm--interam contrast, representative of the radii r>2\kpc where the spiral arms are broader, is not compatible with straightforward compression: lower arm--interarm contrasts than expected may be due to resolution effects and \emph{decompression} of the magnetic field as it leaves the arms. We suggest a simple method to estimate the turbulent scale in the magneto-ionic medium from the dependence of the standard deviation of the observed Faraday rotation measure on resolution. We thus obtain an estimate of 50\pc for the size of the turbulent eddies.Comment: 21 pages, 18 figures (some at lower resolution than submitted version), accepted for publication in MNRA

    Spin Exchange Rates in Electron-Hydrogen Collisions

    Get PDF
    The spin temperature of neutral hydrogen, which determines the 21 cm optical depth and brightness temperature, is set by the competition between radiative and collisional processes. In the high-redshift intergalactic medium, the dominant collisions are typically those between hydrogen atoms. However, collisions with electrons couple much more efficiently to the spin state of hydrogen than do collisions with other hydrogen atoms and thus become important once the ionized fraction exceeds ~1%. Here we compute the rate at which electron-hydrogen collisions change the hydrogen spin. Previous calculations included only S-wave scattering and ignored resonances near the n=2 threshold. We provide accurate results, including all partial wave terms through the F-wave, for the de-excitation rate at temperatures T_K < 15,000 K; beyond that point, excitation to n>=2 hydrogen levels becomes significant. Accurate electron-hydrogen collision rates at higher temperatures are not necessary, because collisional excitation in this regime inevitably produces Lyman-alpha photons, which in turn dominate spin exchange when T_K > 6200 K even in the absence of radiative sources. Our rates differ from previous calculations by several percent over the temperature range of interest. We also consider some simple astrophysical examples where our spin de-excitation rates are useful.Comment: submitted to MNRAS, 9 pages, 5 figure

    A theoretical investigation of the effect of proliferation & adhesion on monoclonal conversion in the colonic crypt

    Get PDF
    The surface epithelium lining the intestinal tract renews itself rapidly by a coordinated programme of cell proliferation, migration and differentiation events that is initiated in the crypts of Lieberkühn. It is generally believed that colorectal cancer arises due to mutations that disrupt the normal cellular dynamics of the crypts. Using a spatially structured cell-based model of a colonic crypt, we investigate the likelihood that the progeny of a mutated cell will dominate, or be sloughed out of, a crypt. Our approach is to perform multiple simulations, varying the spatial location of the initial mutation, and the proliferative and adhesive properties of the mutant cells, to obtain statistical distributions for the probability of their domination. Our simulations lead us to make a number of predictions. The process of monoclonal conversion always occurs, and does not require that the cell which initially gave rise to the population remains in the crypt. Mutations occurring more than one to two cells from the base of the crypt are unlikely to become the dominant clone. The probability of a mutant clone persisting in the crypt is sensitive to dysregulation of adhesion. By comparing simulation results with those from a simple one-dimensional stochastic model of population dynamics at the base of the crypt, we infer that this sensitivity is due to direct competition between wild-type and mutant cells at the base of the crypt. We also predict that increases in the extent of the spatial domain in which the mutant cells proliferate can give rise to counter-intuitive, non-linear changes to the probability of their fixation, due to effects that cannot be captured in simpler models

    A theoretical investigation of the effect of proliferation and\ud adhesion on monoclonal conversion in the colonic crypt

    Get PDF
    Colorectal cancers are initiated by the accumulation of mutations in the colonic epithelium. Using a spatially structured cell-based model of a colonic crypt, we investigate the likelihood that the progeny of a mutated cell will dominate, or be sloughed out of, a crypt. Our approach is to perform multiple simulations, varying the spatial location of the initial mutation, and its proliferative and adhesive properties, to obtain statistical distributions for the probability of domination. Our simulations lead us to make a number of predictions. The process of monoclonal conversion always occurs, and does not require that the cell which initially gave rise to the population remains in the crypt. Mutations occurring more than one to two cells from the base of the crypt are unlikely to become the dominant clone. The probability of a mutant clone persisting in the crypt is sensitive to dysregulation of adhesion, and comparison with a one-dimensional model suggests that this is caused by competition directly at the base of the crypt.\ud We also predict that increases in the extent of the spatial domain in which the mutant cells proliferate cause counter-intuitive non-linear changes to the probability of its fixation, due to effects that cannot be captured in simpler models

    Predicting Proton-Air Cross Sections at sqrt s ~30 TeV, using Accelerator and Cosmic Ray Data

    Full text link
    We use the high energy predictions of a QCD-inspired parameterization of all accelerator data on forward proton-proton and antiproton-proton scattering amplitudes, along with Glauber theory, to predict proton-air cross sections at energies near \sqrt s \approx 30 TeV. The parameterization of the proton-proton cross section incorporates analyticity and unitarity, and demands that the asymptotic proton is a black disk of soft partons. By comparing with the p-air cosmic ray measurements, our analysis results in a constraint on the inclusive particle production cross section.Comment: 9 pages, Revtex, uses epsfig.sty, 5 postscript figures. Minor text revisions. Systematic errors in k included, procedure for extracting k clarified. Previously undefined symbols now define

    The supernova-regulated ISM. I. The multi-phase structure

    Get PDF
    We simulate the multi-phase interstellar medium randomly heated and stirred by supernovae, with gravity, differential rotation and other parameters of the solar neighbourhood. Here we describe in detail both numerical and physical aspects of the model, including injection of thermal and kinetic energy by SN explosions, radiative cooling, photoelectric heating and various transport processes. With 3D domain extending 1 kpc^2 horizontally and 2 kpc vertically, the model routinely spans gas number densities 10^-5 - 10^2 cm^-3, temperatures 10-10^8 K, local velocities up to 10^3 km s^-1 (with Mach number up to 25). The thermal structure of the modelled ISM is classified by inspection of the joint probability density of the gas number density and temperature. We confirm that most of the complexity can be captured in terms of just three phases, separated by temperature borderlines at about 10^3 K and 5x10^5 K. The probability distribution of gas density within each phase is approximately lognormal. We clarify the connection between the fractional volume of a phase and its various proxies, and derive an exact relation between the fractional volume and the filling factors defined in terms of the volume and probabilistic averages. These results are discussed in both observational and computational contexts. The correlation scale of the random flows is calculated from the velocity autocorrelation function; it is of order 100 pc and tends to grow with distance from the mid-plane. We use two distinct parameterizations of radiative cooling to show that the multi-phase structure of the gas is robust, as it does not depend significantly on this choice.Comment: 28 pages, 22 figures and 8 table

    Transverse Magnetoresistance of GaAs/AlGaAs Heterojunctions in the Presence of Parallel Magnetic Fields

    Full text link
    We have calculated the resistivity of a GaAs\slash AlGaAs heterojunction in the presence of both an in--plane magnetic field and a weak perpendicular component using a semiclassical Boltzmann transport theory. These calculations take into account fully the distortion of the Fermi contour which is induced by the parallel magnetic field. The scattering of electrons is assumed to be due to remote ionized impurities. A positive magnetoresistance is found as a function of the perpendicular component, in good qualitative agreement with experimental observations. The main source of this effect is the strong variation of the electronic scattering rate around the Fermi contour which is associated with the variation in the mean distance of the electronic states from the remote impurities. The magnitude of the positive magnetoresistance is strongly correlated with the residual acceptor impurity density in the GaAs layer. The carrier lifetime anisotropy also leads to an observable anisotropy in the resistivity with respect to the angle between the current and the direction of the in--plane magnetic field.Comment: uuencoded file containing a 26 page RevTex file and 14 postscript figures. Submitted to Phys. Rev.

    The Elusive p-air Cross Section

    Full text link
    For the \pbar p and pppp systems, we have used all of the extensive data of the Particle Data Group[K. Hagiwara {\em et al.} (Particle Data Group), Phys. Rev. D 66, 010001 (2002).]. We then subject these data to a screening process, the ``Sieve'' algorithm[M. M. Block, physics/0506010.], in order to eliminate ``outliers'' that can skew a χ2\chi^2 fit. With the ``Sieve'' algorithm, a robust fit using a Lorentzian distribution is first made to all of the data to sieve out abnormally high \delchi, the individual ith^{\rm th} point's contribution to the total χ2\chi^2. The χ2\chi^2 fits are then made to the sieved data. We demonstrate that we cleanly discriminate between asymptotic lns\ln s and ln2s\ln^2s behavior of total hadronic cross sections when we require that these amplitudes {\em also} describe, on average, low energy data dominated by resonances. We simultaneously fit real analytic amplitudes to the ``sieved'' high energy measurements of pˉp\bar p p and pppp total cross sections and ρ\rho-values for s6\sqrt s\ge 6 GeV, while requiring that their asymptotic fits smoothly join the the σpˉp\sigma_{\bar p p} and σpp\sigma_{pp} total cross sections at s=\sqrt s=4.0 GeV--again {\em both} in magnitude and slope. Our results strongly favor a high energy ln2s\ln^2s fit, basically excluding a lns\ln s fit. Finally, we make a screened Glauber fit for the p-air cross section, using as input our precisely-determined pppp cross sections at cosmic ray energies.Comment: 15 pages, 6 figures, 2 table,Paper delivered at c2cr2005 Conference, Prague, September 7-13, 2005. Fig. 2 was missing from V1. V3 fixes all figure
    corecore