A non-local problem for the Fokker-Planck equation related to the Becker-D\"{o}ring model

This paper concerns a Fokker-Planck equation on the positive real line modeling nucleation and growth of clusters. The main feature of the equation is the dependence of the driving vector field and boundary condition on a non-local order parameter related to the excess mass of the system. The first main result concerns the well-posedness and regularity of the Cauchy problem. The well-posedness is based on a fixed point argument, and the regularity on Schauder estimates. The first a priori estimates yield H\"older regularity of the non-local order parameter, which is improved by an iteration argument. The asymptotic behavior of solutions depends on some order parameter $\rho$ depending on the initial data. The system shows different behavior depending on a value $\rho_s>0$, determined from the potentials and diffusion coefficient. For $\rho \leq \rho_s$, there exists an equilibrium solution $c^{\text{eq}}_{(\rho)}$. If $\rho\le\rho_s$ the solution converges strongly to $c^{\text{eq}}_{(\rho)}$, while if $\rho > \rho_s$ the solution converges weakly to $c^{\text{eq}}_{(\rho_s)}$. The excess $\rho - \rho_s$ gets lost due to the formation of larger and larger clusters. In this regard, the model behaves similarly to the classical Becker-D\"oring equation. The system possesses a free energy, strictly decreasing along the evolution, which establishes the long time behavior. In the subcritical case $\rho<\rho_s$ the entropy method, based on suitable weighted logarithmic Sobolev inequalities and interpolation estimates, is used to obtain explicit convergence rates to the equilibrium solution. The close connection of the presented model and the Becker-D\"oring model is outlined by a family of discrete Fokker-Planck type equations interpolating between both of them. This family of models possesses a gradient flow structure, emphasizing their commonality.Comment: Minor revised version accepted for publication in Discrete & Continuous Dynamical Systems -

A benign juvenile environment reduces the strength of antagonistic pleiotropy and genetic variation in the rate of senescence

The environment can play an important role in the evolution of senescence because the optimal allocation between somatic maintenance and reproduction depends on external factors influencing life expectancy. The aims of this study were to experimentally test whether environmental conditions during early life can shape senescence schedules, and if so, to examine whether variation among individuals or genotypes with respect to the degree of ageing differs across environments. We tested life-history plasticity and quantified genetic effects on the pattern of senescence across different environments within a reaction norm framework by using an experiment on the three-spined stickleback (Gasterosteus aculeatus, Linnaeus) in which F1 families originating from a wild annual population experienced different temperature regimes. Male sticklebacks that had experienced a more benign environment earlier in life subsequently reduced their investment in carotenoid-based sexual signals early in the breeding season, and consequently senesced at a slower rate later in the season, compared to those that had developed under harsher conditions. This plasticity of ageing was genetically determined. Both antagonistic pleiotropy and genetic variation in the rate of senescence were evident only in the individuals raised in the harsher environment. The experimental demonstration of genotype-by-environment interactions influencing the rate of reproductive senescence provides interesting insights into the role of the environment in the evolution of life histories. The results suggest that benign conditions weaken the scope for senescence to evolve and that the dependence on the environment may maintain genetic variation under selection

Formation of Kuiper Belt Binaries

The discovery that a substantial fraction of Kuiper Belt objects (KBOs) exists in binaries with wide separations and roughly equal masses, has motivated a variety of new theories explaining their formation. Goldreich et al. (2002) proposed two formation scenarios: In the first, a transient binary is formed, which becomes bound with the aid of dynamical friction from the sea of small bodies (L^2s mechanism); in the second, a binary is formed by three body gravitational deflection (L^3 mechanism). Here, we accurately calculate the L^2s and L^3 formation rates for sub-Hill velocities. While the L^2s formation rate is close to previous order of magnitude estimates, the L^3 formation rate is about a factor of 4 smaller. For sub-Hill KBO velocities (v << v_H) the ratio of the L^3 to the L^2s formation rate is 0.05 (v/v_H) independent of the small bodies' velocity dispersion, their surface density or their mutual collisions. For Super-Hill velocities (v >> v_H) the L^3 mechanism dominates over the L^2s mechanism. Binary formation via the L^3 mechanism competes with binary destruction by passing bodies. Given sufficient time, a statistical equilibrium abundance of binaries forms. We show that the frequency of long-lived transient binaries drops exponentially with the system's lifetime and that such transient binaries are not important for binary formation via the L^3 mechanism, contrary to Lee et al. (2007). For the L^2s mechanism we find that the typical time, transient binaries must last, to form Kuiper Belt binaries (KBBs) for a given strength of dynamical friction, D, increases only logarithmically with D. Longevity of transient binaries only becomes important for very weak dynamical friction (i.e. D \lesssim 0.002) and is most likely not crucial for KBB formation.Comment: 20 pages, 3 figures, Accepted for publication in ApJ, correction of minor typo

Constraints on the Growth and Spin of the Supermassive Black Hole in M32 From High Cadence Visible Light Observations

We present 1-second cadence observations of M32 (NGC221) with the CHIMERA instrument at the Hale 200-inch telescope of the Palomar Observatory. Using field stars as a baseline for relative photometry, we are able to construct a light curve of the nucleus in the g-prime and r-prime band with 1sigma=36 milli-mag photometric stability. We derive a temporal power spectrum for the nucleus and find no evidence for a time-variable signal above the noise as would be expected if the nuclear black hole were accreting gas. Thus, we are unable to constrain the spin of the black hole although future work will use this powerful instrument to target more actively accreting black holes. Given the black hole mass of (2.5+/-0.5)*10^6 Msun inferred from stellar kinematics, the absence of a contribution from a nuclear time-variable signal places an upper limit on the accretion rate which is 4.6*10^{-8} of the Eddington rate, a factor of two more stringent than past upper limits from HST. The low mass of the black hole despite the high stellar density suggests that the gas liberated by stellar interactions was primarily at early cosmic times when the low-mass black hole had a small Eddington luminosity. This is at least partly driven by a top-heavy stellar initial mass function at early cosmic times which is an efficient producer of stellar mass black holes. The implication is that supermassive black holes likely arise from seeds formed through the coalescence of 3-100 Msun mass black holes that then accrete gas produced through stellar interaction processes.Comment: 8 pages, 3 figures, submitted to the Astrophysical Journal, comments welcom

Analytic study of the urn model for separation of sand

We present an analytic study of the urn model for separation of sand recently introduced by Lipowski and Droz (Phys. Rev. E 65, 031307 (2002)). We solve analytically the master equation and the first-passage problem. The analytic results confirm the numerical results obtained by Lipowski and Droz. We find that the stationary probability distribution and the shortest one among the characteristic times are governed by the same free energy. We also analytically derive the form of the critical probability distribution on the critical line, which supports their results obtained by numerically calculating Binder cumulants (cond-mat/0201472).Comment: 6 pages including 3 figures, RevTe

Out of equilibrium dynamics of coherent non-abelian gauge fields

We study out-of-equilibrium dynamics of intense non-abelian gauge fields. Generalizing the well-known Nielsen-Olesen instabilities for constant initial color-magnetic fields, we investigate the impact of temporal modulations and fluctuations in the initial conditions. This leads to a remarkable coexistence of the original Nielsen-Olesen instability and the subdominant phenomenon of parametric resonance. Taking into account that the fields may be correlated only over a limited transverse size, we model characteristic aspects of the dynamics of color flux tubes relevant in the context of heavy-ion collisions.Comment: 12 pages, 10 figures; PRD version, minor change

Analytic study of the three-urn model for separation of sand

We present an analytic study of the three-urn model for separation of sand. We solve analytically the master equation and the first-passage problem. We find that the stationary probability distribution obeys the detailed balance and is governed by the {\it free energy}. We find that the characteristic lifetime of a cluster diverges algebraically with exponent 1/3 at the limit of stability.Comment: 5pages, 4 figures include

Jet pumps for thermoacoustic applications: design guidelines based on a numerical parameter study

The oscillatory flow through tapered cylindrical tube sections (jet pumps) is characterized by a numerical parameter study. The shape of a jet pump results in asymmetric hydrodynamic end effects which cause a time-averaged pressure drop to occur under oscillatory flow conditions. Hence, jet pumps are used as streaming suppressors in closed-loop thermoacoustic devices. A two-dimensional axisymmetric computational fluid dynamics model is used to calculate the performance of a large number of conical jet pump geometries in terms of time-averaged pressure drop and acoustic power dissipation. The investigated geometrical parameters include the jet pump length, taper angle, waist diameter and waist curvature. In correspondence with previous work, four flow regimes are observed which characterize the jet pump performance and dimensionless parameters are introduced to scale the performance of the various jet pump geometries. The simulation results are compared to an existing quasi-steady theory and it is shown that this theory is only applicable in a small operation region. Based on the scaling parameters, an optimum operation region is defined and design guidelines are proposed which can be directly used for future jet pump design.Comment: The following article has been accepted by the Journal of the Acoustical Society of America. After it is published, it will be found at http://scitation.aip.org/JAS

Breakdown of Burton-Prim-Slichter approach and lateral solute segregation in radially converging flows

A theoretical study is presented of the effect of a radially converging melt flow, which is directed away from the solidification front, on the radial solute segregation in simple solidification models. We show that the classical Burton-Prim-Slichter (BPS) solution describing the effect of a diverging flow on the solute incorporation into the solidifying material breaks down for the flows converging along the solidification front. The breakdown is caused by a divergence of the integral defining the effective boundary layer thickness which is the basic concept of the BPS theory. Although such a divergence can formally be avoided by restricting the axial extension of the melt to a layer of finite height, radially uniform solute distributions are possible only for weak melt flows with an axial velocity away from the solidification front comparable to the growth rate. There is a critical melt velocity for each growth rate at which the solution passes through a singularity and becomes physically inconsistent for stronger melt flows. To resolve these inconsistencies we consider a solidification front presented by a disk of finite radius $R_0$ subject to a strong converging melt flow and obtain an analytic solution showing that the radial solute concentration depends on the radius $r$ as $\sim\ln^{1/3}(R_0/r)$ and $\sim\ln(R_0/r)$ close to the rim and at large distances from it. The logarithmic increase of concentration is limited in the vicinity of the symmetry axis by the diffusion becoming effective at a distance comparable to the characteristic thickness of the solute boundary layer. The converging flow causes a solute pile-up forming a logarithmic concentration peak at the symmetry axis which might be an undesirable feature for crystal growth processes.Comment: 15 pages, 5 figure