Quantitative estimates of relationships between geomagnetic activity and equatorial spread-F as determined by TID occurrence levels

Using a world-wide set of stations for 15 years, quantitative estimates of changes to equatorial spread-F (ESF) occurrence rates obtained from ionogram scalings, have been determined for a range of geomagnetic activity (GA) levels, as well as for four different levels of solar activity. Average occurrence rates were used as a reference. The percentage changes vary significantly depending on these subdivisions. For example for very high GA the inverse association is recorded by a change of -33% for R-z greater than or equal to 150, and -10% for R-z < 50. Using data for 9 years for the equatorial station, Huancayo, these measurements of ESF which indicate the presence of TIDs, have also been investigated by somewhat similar analyses. Additional parameters were used which involved the local times of GA, with the ESF being examined separately for occurrence pre-midnight (PM) and after-midnight (AM). Again the negative changes were most pronounced for high GA in R-z-max years (-21%). This result is for PM ESF for GA at a local time of 1700. There were increased ESF levels (+31%) for AM ESF in R-z-min years for high GA around 2300 LT. This additional knowledge of the influence of GA on ESF occurrence involving not only percentage changes, but these values for a range of parameter levels, may be useful if ever short-term forecasts are needed. There is some discussion on comparisons which can be made between ESF results obtained by coherent scatter from incoherent-scatter equipment and those obtained by ionosondes

Scaling Behavior of the Landau Gauge Overlap Quark Propagator

The properties of the momentum space quark propagator in Landau gauge are examined for the overlap quark action in quenched lattice QCD. Numerical calculations are done on three lattices with different lattice spacings and similar physical volumes to explore the approach of the quark propagator towards the continuum limit. We have calculated the nonperturbative momentum-dependent wavefunction renormalization function $Z(p^2)$ and the nonperturbative mass function $M(p^2)$ for a variety of bare quark masses and extrapolate to the chiral limit. We find the behavior of $Z(p^2)$ and $M(p^2)$ are in good agreement for the two finer lattices in the chiral limit. The quark condensate is also calculated.Comment: 3 pages, Lattice2003(Chiral fermions

Scaling behavior of quark propagator in full QCD

We study the scaling behavior of the quark propagator on two lattices with similar physical volume in Landau gauge with 2+1 flavors of dynamical quarks in order to test whether we are close to the continuum limit for these lattices. We use configurations generated with an improved staggered (``Asqtad'') action by the MILC collaboration. The calculations are performed on $28^3\times 96$ lattices with lattice spacing $a = 0.09$ fm and on $20^3\times 64$ lattices with lattice spacing $a = 0.12$ fm. We calculate the quark mass function, $M(q^2)$, and the wave-function renormalization function, $Z(q^2)$, for a variety of bare quark masses. Comparing the behavior of these functions on the two sets of lattices we find that both $Z(q^2)$ and $M(q^2)$ show little sensitivity to the ultraviolet cutoff.Comment: 6 pages, 5 figure

Computational science and re-discovery: open-source implementations of ellipsoidal harmonics for problems in potential theory

We present two open-source (BSD) implementations of ellipsoidal harmonic expansions for solving problems of potential theory using separation of variables. Ellipsoidal harmonics are used surprisingly infrequently, considering their substantial value for problems ranging in scale from molecules to the entire solar system. In this article, we suggest two possible reasons for the paucity relative to spherical harmonics. The first is essentially historical---ellipsoidal harmonics developed during the late 19th century and early 20th, when it was found that only the lowest-order harmonics are expressible in closed form. Each higher-order term requires the solution of an eigenvalue problem, and tedious manual computation seems to have discouraged applications and theoretical studies. The second explanation is practical: even with modern computers and accurate eigenvalue algorithms, expansions in ellipsoidal harmonics are significantly more challenging to compute than those in Cartesian or spherical coordinates. The present implementations reduce the "barrier to entry" by providing an easy and free way for the community to begin using ellipsoidal harmonics in actual research. We demonstrate our implementation using the specific and physiologically crucial problem of how charged proteins interact with their environment, and ask: what other analytical tools await re-discovery in an era of inexpensive computation?Comment: 25 pages, 3 figure

An Exactly Conservative Integrator for the n-Body Problem

The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the computed value of the total energy and angular momentum. Even symplectic integration schemes exactly conserve only an approximate Hamiltonian. We present an algorithm that conserves the true Hamiltonian and the total angular momentum to machine precision. It is derived by applying conventional discretizations in a new space obtained by transformation of the dependent variables. We develop the method first for the restricted circular three-body problem, then for the general two-dimensional three-body problem, and finally for the planar n-body problem. Jacobi coordinates are used to reduce the two-dimensional n-body problem to an (n-1)-body problem that incorporates the constant linear momentum and center of mass constraints. For a four-body choreography, we find that a larger time step can be used with our conservative algorithm than with symplectic and conventional integrators.Comment: 17 pages, 3 figures; to appear in J. Phys. A.: Math. Ge

Casimir-Polder force density between an atom and a conducting wall

In this paper we calculate the Casimir-Polder force density (force per unit area acting on the elements of the surface) on a metallic plate placed in front of a neutral atom. To obtain the force density we use the quantum operator associated to the electromagnetic stress tensor. We explicitly show that the integral of this force density over the plate reproduces the total force acting on the plate. This result shows that, although the force is obtained as a sum of surface element-atom contributions, the stress-tensor method includes also nonadditive components of Casimir-Polder forces in the evaluation of the force acting on a macroscopic object.Comment: 5 page

Dynamic stereo microscopy for studying particle sedimentation

We demonstrate a new method for measuring the sedimentation of a single colloidal bead by using a combination of optical tweezers and a stereo microscope based on a spatial light modulator. We use optical tweezers to raise a micron-sized silica bead to a fixed height and then release it to observe its 3D motion while it sediments under gravity. This experimental procedure provides two independent measurements of bead diameter and a measure of Faxén’s correction, where the motion changes due to presence of the boundary

Space and Time pattern of mid-velocity IMF emission in peripheral heavy-ion collisions at Fermi energies

The emission pattern in the V_perp - V_par plane of Intermediate Mass Fragments with Z=3-7 (IMF) has been studied in the collision 116Sn + 93Nb at 29.5 AMeV as a function of the Total Kinetic Energy Loss of the reaction. This pattern shows that for peripheral reactions most of IMF's are emitted at mid-velocity. Coulomb trajectory calculations demonstrate that these IMF's are produced in the early stages of the reaction and shed light on geometrical details of these emissions, suggesting that the IMF's originate both from the neck and the surface of the interacting nuclei.Comment: 4 pages, 3 figures, RevTex 3.1, submitted to Phys. Rev. Letter

Z-dependent Barriers in Multifragmentation from Poissonian Reducibility and Thermal Scaling

We explore the natural limit of binomial reducibility in nuclear multifragmentation by constructing excitation functions for intermediate mass fragments (IMF) of a given element Z. The resulting multiplicity distributions for each window of transverse energy are Poissonian. Thermal scaling is observed in the linear Arrhenius plots made from the average multiplicity of each element. ``Emission barriers'' are extracted from the slopes of the Arrhenius plots and their possible origin is discussed.Comment: 15 pages including 4 .ps figures. Submitted to Phys. Rev. Letters. Also available at http://csa5.lbl.gov/moretto

A non-perturbative renormalization group study of the stochastic Navier--Stokes equation

We study the renormalization group flow of the average action of the stochastic Navier--Stokes equation with power-law forcing. Using Galilean invariance we introduce a non-perturbative approximation adapted to the zero frequency sector of the theory in the parametric range of the H\"older exponent $4-2\,\varepsilon$ of the forcing where real-space local interactions are relevant. In any spatial dimension $d$, we observe the convergence of the resulting renormalization group flow to a unique fixed point which yields a kinetic energy spectrum scaling in agreement with canonical dimension analysis. Kolmogorov's -5/3 law is, thus, recovered for $\varepsilon=2$ as also predicted by perturbative renormalization. At variance with the perturbative prediction, the -5/3 law emerges in the presence of a \emph{saturation} in the $\varepsilon$-dependence of the scaling dimension of the eddy diffusivity at $\varepsilon=3/2$ when, according to perturbative renormalization, the velocity field becomes infra-red relevant.Comment: RevTeX, 18 pages, 5 figures. Minor changes and new discussion