Dynamics towards the Feigenbaum attractor

We expose at a previously unknown level of detail the features of the dynamics of trajectories that either evolve towards the Feigenbaum attractor or are captured by its matching repellor. Amongst these features are the following: i) The set of preimages of the attractor and of the repellor are embedded (dense) into each other. ii) The preimage layout is obtained as the limiting form of the rank structure of the fractal boundaries between attractor and repellor positions for the family of supercycle attractors. iii) The joint set of preimages for each case form an infinite number of families of well-defined phase-space gaps in the attractor or in the repellor. iv) The gaps in each of these families can be ordered with decreasing width in accord to power laws and are seen to appear sequentially in the dynamics generated by uniform distributions of initial conditions. v) The power law with log-periodic modulation associated to the rate of approach of trajectories towards the attractor (and to the repellor) is explained in terms of the progression of gap formation. vi) The relationship between the law of rate of convergence to the attractor and the inexhaustible hierarchy feature of the preimage structure is elucidated.Comment: 8 pages, 12 figure

Rheology of a sonofluidized granular packing

We report experimental measurements on the rheology of a dry granular material under a weak level of vibration generated by sound injection. First, we measure the drag force exerted on a wire moving in the bulk. We show that when the driving vibration energy is increased, the effective rheology changes drastically: going from a non-linear dynamical friction behavior - weakly increasing with the velocity- up to a linear force-velocity regime. We present a simple heuristic model to account for the vanishing of the stress dynamical threshold at a finite vibration intensity and the onset of a linear force-velocity behavior. Second, we measure the drag force on spherical intruders when the dragging velocity, the vibration energy, and the diameters are varied. We evidence a so-called ''geometrical hardening'' effect for smaller size intruders and a logarithmic hardening effect for the velocity dependence. We show that this last effect is only weakly dependent on the vibration intensity.Comment: Accepted to be published in EPJE. v3: Includes changes suggested by referee

Labyrinthine pathways towards supercycle attractors in unimodal maps

We uncover previously unknown properties of the family of periodic superstable cycles in unimodal maps characterized each by a Lyapunov exponent that diverges to minus infinity. Amongst the main novel properties are the following: i) The basins of attraction for the phases of the cycles develop fractal boundaries of increasing complexity as the period-doubling structure advances towards the transition to chaos. ii) The fractal boundaries, formed by the preimages of the repellor, display hierarchical structures organized according to exponential clusterings that manifest in the dynamics as sensitivity to the final state and transient chaos. iii) There is a functional composition renormalization group (RG) fixed-point map associated to the family of supercycles. iv) This map is given in closed form by the same kind of $q$-exponential function found for both the pitchfork and tangent bifurcation attractors. v) There is a final stage ultra-fast dynamics towards the attractor with a sensitivity to initial conditions that decreases as an exponential of an exponential of time.Comment: 8 pages, 13 figure

Entropies for severely contracted configuration space

We demonstrate that dual entropy expressions of the Tsallis type apply naturally to statistical-mechanical systems that experience an exceptional contraction of their configuration space. The entropic index $\alpha>1$ describes the contraction process, while the dual index $\alpha ^{\prime }=2-\alpha<1$ defines the contraction dimension at which extensivity is restored. We study this circumstance along the three routes to chaos in low-dimensional nonlinear maps where the attractors at the transitions, between regular and chaotic behavior, drive phase-space contraction for ensembles of trajectories. We illustrate this circumstance for properties of systems that find descriptions in terms of nonlinear maps. These are size-rank functions, urbanization and similar processes, and settings where frequency locking takes place

Application of the gradient method to Hartree-Fock-Bogoliubov theory

A computer code is presented for solving the equations of Hartree-Fock-Bogoliubov (HFB) theory by the gradient method, motivated by the need for efficient and robust codes to calculate the configurations required by extensions of HFB such as the generator coordinate method. The code is organized with a separation between the parts that are specific to the details of the Hamiltonian and the parts that are generic to the gradient method. This permits total flexibility in choosing the symmetries to be imposed on the HFB solutions. The code solves for both even and odd particle number ground states, the choice determined by the input data stream. Application is made to the nuclei in the $sd$-shell using the USDB shell-model Hamiltonian.Comment: 20 pages, 5 figures, 3 table

Pairing gaps in Hartree-Fock Bogoliubov theory with the Gogny D1S interaction

As part of a program to study odd-A nuclei in the Hartree-Fock-Bogoliubov (HFB) theory, we have developed a new calculational tool to find the HFB minima of odd-A nuclei based on the gradient method and using interactions of Gogny's form. The HFB minimization includes both time-even and time-odd fields in the energy functional, avoiding the commonly used "filling approximation". Here we apply the method to calculate neutron pairing gaps in some representative isotope chains of spherical and deformed nuclei, namely the Z=8,50 and 82 spherical chains and the Z=62 and 92 deformed chains. We find that the gradient method is quite robust, permitting us to carry out systematic surveys involving many nuclei. We find that the time-odd field does not have large effect on the pairing gaps calculated with the Gogny D1S interaction. Typically, adding the T-odd field as a perturbation increases the pairing gap by ~100 keV, but the re-minimization brings the gap back down. This outcome is very similar to results reported for the Skyrme family of nuclear energy density functionals. Comparing the calculated gaps with the experimental ones, we find that the theoretical errors have both signs implying that the D1S interaction has a reasonable overall strength. However, we find some systematic deficiencies comparing spherical and deformed chains and comparing the lighter chains with the heavier ones. The gaps for heavy spherical nuclei are too high, while those for deformed nuclei tend to be too low. The calculated gaps of spherical nuclei show hardly any A-dependence, contrary to the data. Inclusion of the T-odd component of the interaction does not change these qualitative findings

Incidence of qq-statistics in rank distributions

We show that size-rank distributions with power-law decay (often only over a limited extent) observed in a vast number of instances in a widespread family of systems obey Tsallis statistics. The theoretical framework for these distributions is analogous to that of a nonlinear iterated map near a tangent bifurcation for which the Lyapunov exponent is negligible or vanishes. The relevant statistical-mechanical expressions associated with these distributions are derived from a maximum entropy principle with the use of two different constraints, and the resulting duality of entropy indexes is seen to portray physically relevant information. While the value of the index $\alpha$ fixes the distribution's power-law exponent, that for the dual index $2-\alpha$ ensures the extensivity of the deformed entropy.Comment: Santa Fe Institute working paper: http://www.santafe.edu/media/workingpapers/14-07-024.pdf. see: http://www.pnas.org/content/early/2014/09/03/1412093111.full.pdf+htm