Nuclear response for the Skyrme effective interaction with zero-range tensor terms. II. Sum rules and instabilities

The formalism of linear response theory for Skyrme forces including tensor terms presented in article [1] is generalized for the case of a Skyrme energy density functional in infinite matter. We also present analytical results for the odd-power sum rules, with particular attention to the inverse energy weighted sum rule, $M_{-1}$, as a tool to detect instabilities in Skyrme functionals.Comment: Submitted to Phys. Rev.

Understanding fragility in supercooled Lennard-Jones mixtures. II. Potential energy surface

We numerically investigated the connection between isobaric fragility and the properties of high-order stationary points of the potential energy surface in different supercooled Lennard-Jones mixtures. The increase of effective activation energies upon supercooling appears to be driven by the increase of average potential energy barriers measured by the energy dependence of the fraction of unstable modes. Such an increase is sharper, the more fragile is the mixture. Correlations between fragility and other properties of high-order stationary points, including the vibrational density of states and the localization features of unstable modes, are also discussed.Comment: 13 pages, 13 figures, minor revisions, one figure adde

Pairing correlations of cold fermionic gases at overflow from a narrow to a wide harmonic trap

Within the context of Hartree-Fock-Bogoliubov theory, we study the behavior of superfluid Fermi systems when they pass from a small to a large container. Such systems can be now realized thanks to recent progress in experimental techniques. It will allow to better understand pairing properties at overflow and in general in rapidly varying external potentials

Nuclear response for the Skyrme effective interaction with zero-range tensor terms. III. Neutron matter and neutrino propagation

The formalism of the linear response for the Skyrme energy density functional including tensor terms derived in articles [1,2] for nuclear matter is applied here to the case of pure neutron matter. As in article [2] we present analytical results for the response function in all channels, the Landau parameters and the odd-power sum rules. Special emphasis is given to the inverse energy weighted sum rule because it can be used to detect non physical instabilities. Typical examples are discussed and numerical results shown. Moreover, as a direct application, neutrino propagation in neutron matter is investigated through its neutrino mean free path at zero temperature. This quantity turns out to be very sensitive to the tensor terms of the Skyrme energy density functional

Dynamics of conduction blocks in a model of paced cardiac tissue

We study numerically the dynamics of conduction blocks using a detailed electrophysiological model. We find that this dynamics depends critically on the size of the paced region. Small pacing regions lead to stationary conduction blocks while larger pacing regions can lead to conduction blocks that travel periodically towards the pacing region. We show that this size-dependence dynamics can lead to a novel arrhythmogenic mechanism. Furthermore, we show that the essential phenomena can be captured in a much simpler coupled-map model.Comment: 8 pages 6 figure

Nuclear energy density optimization: Shell structure

Nuclear density functional theory is the only microscopical theory that can be applied throughout the entire nuclear landscape. Its key ingredient is the energy density functional. In this work, we propose a new parameterization UNEDF2 of the Skyrme energy density functional. The functional optimization is carried out using the POUNDerS optimization algorithm within the framework of the Skyrme Hartree-Fock-Bogoliubov theory. Compared to the previous parameterization UNEDF1, restrictions on the tensor term of the energy density have been lifted, yielding a very general form of the energy density functional up to second order in derivatives of the one-body density matrix. In order to impose constraints on all the parameters of the functional, selected data on single-particle splittings in spherical doubly-magic nuclei have been included into the experimental dataset. The agreement with both bulk and spectroscopic nuclear properties achieved by the resulting UNEDF2 parameterization is comparable with UNEDF1. While there is a small improvement on single-particle spectra and binding energies of closed shell nuclei, the reproduction of fission barriers and fission isomer excitation energies has degraded. As compared to previous UNEDF parameterizations, the parameter confidence interval for UNEDF2 is narrower. In particular, our results overlap well with those obtained in previous systematic studies of the spin-orbit and tensor terms. UNEDF2 can be viewed as an all-around Skyrme EDF that performs reasonably well for both global nuclear properties and shell structure. However, after adding new data aiming to better constrain the nuclear functional, its quality has improved only marginally. These results suggest that the standard Skyrme energy density has reached its limits and significant changes to the form of the functional are needed.Comment: 18 pages, 13 figures, 12 tables; resubmitted for publication to Phys. Rev. C after second review by refere

Indeterminacy of Spatiotemporal Cardiac Alternans

Cardiac alternans, a beat-to-beat alternation in action potential duration (at the cellular level) or in ECG morphology (at the whole heart level), is a marker of ventricular fibrillation, a fatal heart rhythm that kills hundreds of thousands of people in the US each year. Investigating cardiac alternans may lead to a better understanding of the mechanisms of cardiac arrhythmias and eventually better algorithms for the prediction and prevention of such dreadful diseases. In paced cardiac tissue, alternans develops under increasingly shorter pacing period. Existing experimental and theoretical studies adopt the assumption that alternans in homogeneous cardiac tissue is exclusively determined by the pacing period. In contrast, we find that, when calcium-driven alternans develops in cardiac fibers, it may take different spatiotemporal patterns depending on the pacing history. Because there coexist multiple alternans solutions for a given pacing period, the alternans pattern on a fiber becomes unpredictable. Using numerical simulation and theoretical analysis, we show that the coexistence of multiple alternans patterns is induced by the interaction between electrotonic coupling and an instability in calcium cycling.Comment: 20 pages, 10 figures, to be published in Phys. Rev.

Effects of patch size and number within a simple model of patchy colloids

We report on a computer simulation and integral equation study of a simple model of patchy spheres, each of whose surfaces is decorated with two opposite attractive caps, as a function of the fraction $\chi$ of covered attractive surface. The simple model explored --- the two-patch Kern-Frenkel model --- interpolates between a square-well and a hard-sphere potential on changing the coverage $\chi$. We show that integral equation theory provides quantitative predictions in the entire explored region of temperatures and densities from the square-well limit $\chi = 1.0$ down to $\chi \approx 0.6$. For smaller $\chi$, good numerical convergence of the equations is achieved only at temperatures larger than the gas-liquid critical point, where however integral equation theory provides a complete description of the angular dependence. These results are contrasted with those for the one-patch case. We investigate the remaining region of coverage via numerical simulation and show how the gas-liquid critical point moves to smaller densities and temperatures on decreasing $\chi$. Below $\chi \approx 0.3$, crystallization prevents the possibility of observing the evolution of the line of critical points, providing the angular analog of the disappearance of the liquid as an equilibrium phase on decreasing the range for spherical potentials. Finally, we show that the stable ordered phase evolves on decreasing $\chi$ from a three-dimensional crystal of interconnected planes to a two-dimensional independent-planes structure to a one-dimensional fluid of chains when the one-bond-per-patch limit is eventually reached.Comment: 26 pages, 11 figures, J. Chem. Phys. in pres

Partial wave decomposition of the N3LO equation of state

By means of a partial wave decomposition, we separate their contributions to the equation of state (EoS) of symmetric nuclear matter for the N3LO pseudo-potential. In particular, we show that although both the tensor and the spin-orbit terms do not contribute to the EoS, they give a non-vanishing contribution to the separate (JLS) channels

Partial-wave decomposition of the finite-range effective tensor interaction

We perform a detailed analysis of the properties of the finite-range tensor term associated with the Gogny and M3Y effective interactions. In particular, by using a partial wave decomposition of the equation of state of symmetric nuclear matter, we show how we can extract their tensor parameters directly from microscopic results based on bare nucleon-nucleon interactions. Furthermore, we show that the zero-range limit of both finite-range interactions has the form of the N3LO Skyrme pseudo-potential, which thus constitutes a reliable approximation in the density range relevant for finite nuclei. Finally, we use Brueckner-Hartree-Fock results to fix the tensor parameters for the three effective interactions