Response functions of cold neutron matter: density fluctuations

We compute the finite temperature density response function of nonrelativistic cold fermions with an isotropic condensate. The pair-breaking contribution to the response function is evaluated in the limit of small three-momentum transfers q within an effective theory which exploits series expansion in powers of small q/p_F, where p_F is the Fermi momentum. The leading order O(q^2) contribution is universal and depends only on two fundamental scales, the Fermi energy and the pairing gap. The particle-hole Landau Fermi-liquid interaction contributes first at the next-to-leading-order O(q^4). The scattering contribution to the polarization tensor is nonperturbative (in the above sense) and is evaluated numerically. The spectral functions of density fluctuations are constructed and the relevance of the q^2 scaling for the pair-breaking neutrino emission from neutron stars is discussed.Comment: v2: 11 pages, 4 figures, matches published version

Cheap Newton steps for optimal control problems: automatic differentiation and Pantoja's algorithm

Original article can be found at: http://www.informaworld.com/smpp/title~content=t713645924~db=all Copyright Taylor and Francis / Informa.In this paper we discuss Pantoja's construction of the Newton direction for discrete time optimal control problems. We show that automatic differentiation (AD) techniques can be used to calculate the Newton direction accurately, without requiring extensive re-writing of user code, and at a surprisingly low computational cost: for an N-step problem with p control variables and q state variables at each step, the worst case cost is 6(p + q + 1) times the computational cost of a single target function evaluation, independent of N, together with at most p3/3 + p2(q + 1) + 2p(q + 1)2 + (q + l)3, i.e. less than (p + q + l)3, floating point multiply-and-add operations per time step. These costs may be considerably reduced if there is significant structural sparsity in the problem dynamics. The systematic use of checkpointing roughly doubles the operation counts, but reduces the total space cost to the order of 4pN floating point stores. A naive approach to finding the Newton step would require the solution of an Np Np system of equations together with a number of function evaluations proportional to Np, so this approach to Pantoja's construction is extremely attractive, especially if q is very small relative to N. Straightforward modifications of the AD algorithms proposed here can be used to implement other discrete time optimal control solution techniques, such as differential dynamic programming (DDP), which use state-control feedback. The same techniques also can be used to determine with certainty, at the cost of a single Newton direction calculation, whether or not the Hessian of the target function is sufficiently positive definite at a point of interest. This allows computationally cheap post-hoc verification that a second-order minimum has been reached to a given accuracy, regardless of what method has been used to obtain it.Peer reviewe

A Dynamically Diluted Alignment Model Reveals the Impact of Cell Turnover on the Plasticity of Tissue Polarity Patterns

The polarisation of cells and tissues is fundamental for tissue morphogenesis during biological development and regeneration. A deeper understanding of biological polarity pattern formation can be gained from the consideration of pattern reorganisation in response to an opposing instructive cue, which we here consider by example of experimentally inducible body axis inversions in planarian flatworms. Our dynamically diluted alignment model represents three processes: entrainment of cell polarity by a global signal, local cell-cell coupling aligning polarity among neighbours and cell turnover inserting initially unpolarised cells. We show that a persistent global orienting signal determines the final mean polarity orientation in this stochastic model. Combining numerical and analytical approaches, we find that neighbour coupling retards polarity pattern reorganisation, whereas cell turnover accelerates it. We derive a formula for an effective neighbour coupling strength integrating both effects and find that the time of polarity reorganisation depends linearly on this effective parameter and no abrupt transitions are observed. This allows to determine neighbour coupling strengths from experimental observations. Our model is related to a dynamic $8$-Potts model with annealed site-dilution and makes testable predictions regarding the polarisation of dynamic systems, such as the planarian epithelium.Comment: Preprint as prior to first submission to Journal of the Royal Society Interface. 25 pages, 6 figures, plus supplement (18 pages, contains 1 table and 7 figures). A supplementary movie is available from https://dx.doi.org/10.6084/m9.figshare.c388781

Marx's reproduction schemes and the Keynesian multiplier: a reply to Sardoni

In a recent contribution to this journal, C. Sardoni takes issue with the identification by Trigg, in a 2006 publication, of a role for the Keynesian investment multiplier in Marx's schemes of reproduction. Indirectly, Sardoni also expresses his disagreement with Hartwig (by attributing one of his statements to Trigg). We appreciate the opportunity to defend our view against Sardoni's critique. This reply shows that a bridging point between Marx and Keynes can be established without recourse to microfoundations. As suggested by both Trigg, in 2006, and Hartwig, in 2004, the well known Harrod-Domar model of economic growth provides an interpretation of Marx's reproduction schemes that has the Keynesian multiplier as a constituent element. This note will further explore the assumptions underlying the interface between Marx and Keynes, in response to the challenging questions raised in Sardoni's contributio