A numerical study of Penrose-like inequalities in a family of axially symmetric initial data

Our current picture of black hole gravitational collapse relies on two assumptions: i) the resulting singularity is hidden behind an event horizon -- weak cosmic censorship conjecture -- and ii) spacetime eventually settles down to a stationarity state. In this setting, it follows that the minimal area containing an apparent horizon is bound by the square of the total ADM mass (Penrose inequality conjecture). Following Dain et al. 2002, we construct numerically a family of axisymmetric initial data with one or several marginally trapped surfaces. Penrose and related geometric inequalities are discused for these data. As a by-product, it is shown how Penrose inequality can be used as a diagnosis for an apparent horizon finder numerical routine.Comment: Contribution to the "Encuentros Relativistas Espanoles - Spanish Relativity Meeting ERE07" Proceedings, Tenerife, Spain (September 2007

Survival probability of a diffusing particle in the presence of Poisson-distributed mobile traps

The problem of a diffusing particle moving among diffusing traps is analyzed in general space dimension d. We consider the case where the traps are initially randomly distributed in space, with uniform density rho, and derive upper and lower bounds for the probability Q(t) (averaged over all particle and trap trajectories) that the particle survives up to time t. We show that, for 1<=d<2, the bounds converge asymptotically to give $Q(t) \sim exp(-\lambda_d t^{d/2})$ where $\lambda_d = (2/\pi d) sin(\pi d/2) (4\pi D)^{d/2} \rho$ and D is the diffusion constant of the traps, and that $Q(t) \sim exp(- 4\pi\rho D t/ln t)$ for d=2. For d>2 bounds can still be derived, but they no longer converge for large t. For 1<=d<=2, these asymptotic form are independent of the diffusion constant of the particle. The results are compared with simulation results obtained using a new algorithm [V. Mehra and P. Grassberger, Phys. Rev. E v65 050101 (2002)] which is described in detail. Deviations from the predicted asymptotic forms are found to be large even for very small values of Q(t), indicating slowly decaying corrections whose form is consistent with the bounds. We also present results in d=1 for the case where the trap densities on either side of the particle are different. For this case we can still obtain exact bounds but they no longer converge.Comment: 13 pages, RevTeX4, 6 figures. Figures and references updated; equations corrected; discussion clarifie

Ground state non-universality in the random field Ising model

Two attractive and often used ideas, namely universality and the concept of a zero temperature fixed point, are violated in the infinite-range random-field Ising model. In the ground state we show that the exponents can depend continuously on the disorder and so are non-universal. However, we also show that at finite temperature the thermal order parameter exponent one half is restored so that temperature is a relevant variable. The broader implications of these results are discussed.Comment: 4 pages 2 figures, corrected prefactors caused by a missing factor of two in Eq. 2., added a paragraph in conclusions for clarit

Calculation of ground states of four-dimensional +or- J Ising spin glasses

Ground states of four-dimensional (d=4) EA Ising spin glasses are calculated for sizes up to 7x7x7x7 using a combination of a genetic algorithm and cluster-exact approximation. The ground-state energy of the infinite system is extrapolated as e_0=-2.095(1). The ground-state stiffness (or domain wall) energy D is calculated. A D~L^{\Theta} behavior with \Theta=0.65(4) is found which confirms that the d=4 model has an equilibrium spin-glass-paramagnet transition for non-zero T_c.Comment: 5 pages, 3 figures, 31 references, revtex; update of reference

Cluster persistence in one-dimensional diffusion--limited cluster--cluster aggregation

The persistence probability, $P_C(t)$, of a cluster to remain unaggregated is studied in cluster-cluster aggregation, when the diffusion coefficient of a cluster depends on its size $s$ as $D(s) \sim s^\gamma$. In the mean-field the problem maps to the survival of three annihilating random walkers with time-dependent noise correlations. For $\gamma \ge 0$ the motion of persistent clusters becomes asymptotically irrelevant and the mean-field theory provides a correct description. For $\gamma < 0$ the spatial fluctuations remain relevant and the persistence probability is overestimated by the random walk theory. The decay of persistence determines the small size tail of the cluster size distribution. For $0 < \gamma < 2$ the distribution is flat and, surprisingly, independent of $\gamma$.Comment: 11 pages, 6 figures, RevTeX4, submitted to Phys. Rev.

Hydrodynamics of domain growth in nematic liquid crystals

We study the growth of aligned domains in nematic liquid crystals. Results are obtained solving the Beris-Edwards equations of motion using the lattice Boltzmann approach. Spatial anisotropy in the domain growth is shown to be a consequence of the flow induced by the changing order parameter field (backflow). The generalization of the results to the growth of a cylindrical domain, which involves the dynamics of a defect ring, is discussed.Comment: 12 revtex-style pages, including 12 figures; small changes before publicatio

The three-dimensional random field Ising magnet: interfaces, scaling, and the nature of states

The nature of the zero temperature ordering transition in the 3D Gaussian random field Ising magnet is studied numerically, aided by scaling analyses. In the ferromagnetic phase the scaling of the roughness of the domain walls, $w\sim L^\zeta$, is consistent with the theoretical prediction $\zeta = 2/3$. As the randomness is increased through the transition, the probability distribution of the interfacial tension of domain walls scales as for a single second order transition. At the critical point, the fractal dimensions of domain walls and the fractal dimension of the outer surface of spin clusters are investigated: there are at least two distinct physically important fractal dimensions. These dimensions are argued to be related to combinations of the energy scaling exponent, $\theta$, which determines the violation of hyperscaling, the correlation length exponent $\nu$, and the magnetization exponent $\beta$. The value $\beta = 0.017\pm 0.005$ is derived from the magnetization: this estimate is supported by the study of the spin cluster size distribution at criticality. The variation of configurations in the interior of a sample with boundary conditions is consistent with the hypothesis that there is a single transition separating the disordered phase with one ground state from the ordered phase with two ground states. The array of results are shown to be consistent with a scaling picture and a geometric description of the influence of boundary conditions on the spins. The details of the algorithm used and its implementation are also described.Comment: 32 pp., 2 columns, 32 figure

Fluctuation-dissipation relations in the non-equilibrium critical dynamics of Ising models

We investigate the relation between two-time, multi-spin, correlation and response functions in the non-equilibrium critical dynamics of Ising models in d=1 and d=2 spatial dimensions. In these non-equilibrium situations, the fluctuation-dissipation theorem (FDT) is not satisfied. We find FDT `violations' qualitatively similar to those reported in various glassy materials, but quantitatively dependent on the chosen observable, in contrast to the results obtained in infinite-range glass models. Nevertheless, all FDT violations can be understood by considering separately the contributions from large wavevectors, which are at quasi-equilibrium and obey FDT, and from small wavevectors where a generalized FDT holds with a non-trivial limit fluctuation-dissipation ratio X. In d=1, we get X = 1/2 for spin observables, which measure the orientation of domains, while X = 0 for observables that are sensitive to the domain-wall motion. Numerical simulations in d=2 reveal a unique X = 0.34 for all observables. Measurement protocols for X are discussed in detail. Our results suggest that the definition of an effective temperature Teff = T / X for large length scales is generically possible in non-equilibrium critical dynamics.Comment: 26 pages, 10 figure

Scaling properties in off equilibrium dynamical processes

In the present paper, we analyze the consequences of scaling hypotheses on dynamic functions, as two times correlations $C(t,t')$. We show, under general conditions, that $C(t,t')$ must obey the following scaling behavior $C(t,t') = \phi_1(t)^{f(\beta)}{\cal{S}}(\beta)$, where the scaling variable is $\beta=\beta(\phi_1(t')/\phi_1(t))$ and $\phi_1(t')$, $\phi_1(t)$ two undetermined functions. The presence of a non constant exponent $f(\beta)$ signals the appearance of multiscaling properties in the dynamics.Comment: 6 pages, no figure

Domain growth and aging scaling in coarsening disordered systems

Using extensive Monte Carlo simulations we study aging properties of two disordered systems quenched below their critical point, namely the two-dimensional random-bond Ising model and the three-dimensional Edwards-Anderson Ising spin glass with a bimodal distribution of the coupling constants. We study the two-times autocorrelation and space-time correlation functions and show that in both systems a simple aging scenario prevails in terms of the scaling variable $L(t)/L(s)$, where $L$ is the time-dependent correlation length, whereas $s$ is the waiting time and $t$ is the observation time. The investigation of the space-time correlation function for the random-bond Ising model allows us to address some issues related to superuniversality.Comment: 8 pages, 9 figures, to appear in European Physical Journal