N-vector spin models on the sc and the bcc lattices: a study of the critical behavior of the susceptibility and of the correlation length by high temperature series extended to order beta^{21}

High temperature expansions for the free energy, the susceptibility and the second correlation moment of the classical N-vector model [also known as the O(N) symmetric classical spin Heisenberg model or as the lattice O(N) nonlinear sigma model] on the sc and the bcc lattices are extended to order beta^{21} for arbitrary N. The series for the second field derivative of the susceptibility is extended to order beta^{17}. An analysis of the newly computed series for the susceptibility and the (second moment) correlation length yields updated estimates of the critical parameters for various values of the spin dimensionality N, including N=0 [the self-avoiding walk model], N=1 [the Ising spin 1/2 model], N=2 [the XY model], N=3 [the Heisenberg model]. For all values of N, we confirm a good agreement with the present renormalization group estimates. A study of the series for the other observables will appear in a forthcoming paper.Comment: Revised version to appear in Phys. Rev. B Sept. 1997. Revisions include an improved series analysis biased with perturbative values of the scaling correction exponents computed by A. I. Sokolov. Added a reference to estimates of exponents for the Ising Model. Abridged text of 19 pages, latex, no figures, no tables of series coefficient

Post-Newtonian SPH calculations of binary neutron star coalescence. I. Method and first results

We present the first results from our Post-Newtonian (PN) Smoothed Particle Hydrodynamics (SPH) code, which has been used to study the coalescence of binary neutron star (NS) systems. The Lagrangian particle-based code incorporates consistently all lowest-order (1PN) relativistic effects, as well as gravitational radiation reaction, the lowest-order dissipative term in general relativity. We test our code on sequences of single NS models of varying compactness, and we discuss ways to make PN simulations more relevant to realistic NS models. We also present a PN SPH relaxation procedure for constructing equilibrium models of synchronized binaries, and we use these equilibrium models as initial conditions for our dynamical calculations of binary coalescence. Though unphysical, since tidal synchronization is not expected in NS binaries, these initial conditions allow us to compare our PN work with previous Newtonian results. We compare calculations with and without 1PN effects, for NS with stiff equations of state, modeled as polytropes with $\Gamma=3$. We find that 1PN effects can play a major role in the coalescence, accelerating the final inspiral and causing a significant misalignment in the binary just prior to final merging. In addition, the character of the gravitational wave signal is altered dramatically, showing strong modulation of the exponentially decaying waveform near the end of the merger. We also discuss briefly the implications of our results for models of gamma-ray bursts at cosmological distances.Comment: RevTeX, 37 pages, 17 figures, to appear in Phys. Rev. D, minor corrections onl

Standard random walks and trapping on the Koch network with scale-free behavior and small-world effect

A vast variety of real-life networks display the ubiquitous presence of scale-free phenomenon and small-world effect, both of which play a significant role in the dynamical processes running on networks. Although various dynamical processes have been investigated in scale-free small-world networks, analytical research about random walks on such networks is much less. In this paper, we will study analytically the scaling of the mean first-passage time (MFPT) for random walks on scale-free small-world networks. To this end, we first map the classical Koch fractal to a network, called Koch network. According to this proposed mapping, we present an iterative algorithm for generating the Koch network, based on which we derive closed-form expressions for the relevant topological features, such as degree distribution, clustering coefficient, average path length, and degree correlations. The obtained solutions show that the Koch network exhibits scale-free behavior and small-world effect. Then, we investigate the standard random walks and trapping issue on the Koch network. Through the recurrence relations derived from the structure of the Koch network, we obtain the exact scaling for the MFPT. We show that in the infinite network order limit, the MFPT grows linearly with the number of all nodes in the network. The obtained analytical results are corroborated by direct extensive numerical calculations. In addition, we also determine the scaling efficiency exponents characterizing random walks on the Koch network.Comment: 12 pages, 8 figures. Definitive version published in Physical Review

Renormalization group approach to an Abelian sandpile model on planar lattices

One important step in the renormalization group (RG) approach to a lattice sandpile model is the exact enumeration of all possible toppling processes of sandpile dynamics inside a cell for RG transformations. Here we propose a computer algorithm to carry out such exact enumeration for cells of planar lattices in RG approach to Bak-Tang-Wiesenfeld sandpile model [Phys. Rev. Lett. {\bf 59}, 381 (1987)] and consider both the reduced-high RG equations proposed by Pietronero, Vespignani, and Zapperi (PVZ) [Phys. Rev. Lett. {\bf 72}, 1690 (1994)] and the real-height RG equations proposed by Ivashkevich [Phys. Rev. Lett. {\bf 76}, 3368 (1996)]. Using this algorithm we are able to carry out RG transformations more quickly with large cell size, e.g. $3 \times 3$ cell for the square (sq) lattice in PVZ RG equations, which is the largest cell size at the present, and find some mistakes in a previous paper [Phys. Rev. E {\bf 51}, 1711 (1995)]. For sq and plane triangular (pt) lattices, we obtain the only attractive fixed point for each lattice and calculate the avalanche exponent $\tau$ and the dynamical exponent $z$. Our results suggest that the increase of the cell size in the PVZ RG transformation does not lead to more accurate results. The implication of such result is discussed.Comment: 29 pages, 6 figure

Enhanced heat capacity and a new temperature instability in superfluid He-4 in the presence of a constant heat flux near T-lambda

We present the first experimental evidence that the heat capacity of superfluid 4He, at temperatures very close to the lambda point Tλ, is enhanced by a constant heat flux Q. The heat capacity at constant Q, CQ, is predicted to diverge at a temperature Tc(Q)<Tλ at which superflow becomes unstable. In agreement with previous measurements, we find that dissipation enters our cell at a temperature, TDAS(Q), below the theoretical value, Tc(Q). We argue that TDAS(Q) can be accounted for by a temperature instability at the cell wall, and is therefore distinct from Tc(Q). The excess heat capacity we measure has the predicted scaling behavior as a function of T and Q, but it is much larger than predicted by current theory

Monte Carlo calculations of energy depositions and radiation transport. Volume 1 - Validation of COHORT codes

Monte Carlo codes for IBM 7090 digital computer to calculate radiation heating in propellant tanks, and radiation environment about nuclear rocket stag

Cyclic Distributed Space–Time Codes for Wireless Relay Networks With No Channel Information

In this paper, we present a coding strategy for half duplex wireless relay networks, where we assume no channel knowledge at any of the transmitter, receiver, or relays. The coding scheme uses distributed space–time coding, that is, the relay nodes cooperate to encode the transmitted signal so that the receiver senses a space–time codeword. It is inspired by noncoherent differential techniques. The proposed strategy is available for any number of relays nodes. It is analyzed, and shown to yield a diversity linear in the number of relays. We also study the resistance of the scheme to relay node failures, and show that a network with R relay nodes and d of them down behaves, as far as diversity is concerned, as a network with R-d nodes. Finally, our construction can be easily generalized to the case where the transmitter and receiver nodes have several antennas

Growing dust grains in protoplanetary discs - I. Radial drift with toy growth models

In a series of papers, we present a comprehensive analytic study of the global motion of growing dust grains in protoplanetary discs, addressing both the radial drift and the vertical settling of the particles. Here we study how the radial drift of dust particles is affected by grain growth. In a first step, toy models in which grain growth can either be constant, accelerate or decelerate are introduced. The equations of motion are analytically integrable and therefore the grains dynamics is easy to understand. The radial motion of growing grains is governed by the relative efficiency of the growth and migration processes which is expressed by the dimensionless parameter Lambda, as well as the exponents for the gas surface density and temperature profiles, denoted p and q respectively. When Lambda is of order unity, growth and migration are strongly coupled, providing the most efficient radial drift. For the toy models considered, grains pile up when -p+q+1/2<0. Importantly, we show the existence of a second process which can help discs to retain their solid materials. For accelerating growth, grains end up their migration at a finite radius, thus avoiding being accreted onto the central star.Comment: 12 pages, 9 figures. Accepted for publication in MNRAS. v2: typos correcte

Spacetime as a topological insulator: Mechanism for the origin of the fermion generations

We suggest a mechanism whereby the three generations of quarks and leptons correspond to surface modes in a five-dimensional theory. These modes arise from a nonlinear fermion dispersion relation in the extra dimension, much in the same manner as fermion surface modes in a topological insulator or lattice implementation of domain wall fermions. We also show that the topological properties can persist in a deconstructed version of the model in four dimensions.Comment: Substantially revised version, to appear in Phys. Rev. Let