Simulation of the Zero Temperature Behavior of a 3-Dimensional Elastic Medium

We have performed numerical simulation of a 3-dimensional elastic medium, with scalar displacements, subject to quenched disorder. We applied an efficient combinatorial optimization algorithm to generate exact ground states for an interface representation. Our results indicate that this Bragg glass is characterized by power law divergences in the structure factor $S(k)\sim A k^{-3}$. We have found numerically consistent values of the coefficient $A$ for two lattice discretizations of the medium, supporting universality for $A$ in the isotropic systems considered here. We also examine the response of the ground state to the change in boundary conditions that corresponds to introducing a single dislocation loop encircling the system. Our results indicate that the domain walls formed by this change are highly convoluted, with a fractal dimension $d_f=2.60(5)$. We also discuss the implications of the domain wall energetics for the stability of the Bragg glass phase. As in other disordered systems, perturbations of relative strength $\delta$ introduce a new length scale $L^* \sim \delta^{-1/\zeta}$ beyond which the perturbed ground state becomes uncorrelated with the reference (unperturbed) ground state. We have performed scaling analysis of the response of the ground state to the perturbations and obtain $\zeta = 0.385(40)$. This value is consistent with the scaling relation $\zeta=d_f/2- \theta$, where $\theta$ characterizes the scaling of the energy fluctuations of low energy excitations.Comment: 20 pages, 13 figure

Critical slowing down in polynomial time algorithms

Combinatorial optimization algorithms which compute exact ground state configurations in disordered magnets are seen to exhibit critical slowing down at zero temperature phase transitions. Using arguments based on the physical picture of the model, including vanishing stiffness on scales beyond the correlation length and the ground state degeneracy, the number of operations carried out by one such algorithm, the push-relabel algorithm for the random field Ising model, can be estimated. Some scaling can also be predicted for the 2D spin glass.Comment: 4 pp., 3 fig

Computational Complexity of Determining the Barriers to Interface Motion in Random Systems

The low-temperature driven or thermally activated motion of several condensed matter systems is often modeled by the dynamics of interfaces (co-dimension-1 elastic manifolds) subject to a random potential. Two characteristic quantitative features of the energy landscape of such a many-degree-of-freedom system are the ground-state energy and the magnitude of the energy barriers between given configurations. While the numerical determination of the former can be accomplished in time polynomial in the system size, it is shown here that the problem of determining the latter quantity is NP-complete. Exact computation of barriers is therefore (almost certainly) much more difficult than determining the exact ground states of interfaces.Comment: 8 pages, figures included, to appear in Phys. Rev.

On The Existence Of Anisotropic Cosmological Models In Higher-Order Theories Of Gravity

We investigate the behaviour on approach to the initial singularity in higher-order extensions of general relativity by finding exact cosmological solutions for a wide class of models in which the Lagrangian is allowed to depend nonlinearly upon the three possible linear and quadratic scalars built from the Riemann tensor; $R$, $R_{ab}R^{ab}$ and $R_{abcd}R^{abcd}$. We present new anisotropic vacuum solutions analagous to the Kasner solutions of general relativity and extend previous results to a much wider range of fourth order theories of gravity. We discuss the implications of these results for the behaviour of the more general anisotropic Bianchi type VIII and IX cosmologies as the initial singularity is approached. Furthermore, we also consider the existence conditions for some other simple anisotropic Bianchi I vacuum solutions in which the expansion in each direction is of exponential, rather than power-law behaviour and their relevance for cosmic ``no-hair'' theorems.Comment: 24 pages, submitted to CQ

Hysteretic dynamics of domain walls at finite temperatures

Theory of domain wall motion in a random medium is extended to the case when the driving field is below the zero-temperature depinning threshold and the creep of the domain wall is induced by thermal fluctuations. Subject to an ac drive, the domain wall starts to move when the driving force exceeds an effective threshold which is temperature and frequency-dependent. Similarly to the case of zero-temperature, the hysteresis loop displays three dynamical phase transitions at increasing ac field amplitude $h_0$. The phase diagram in the 3-d space of temperature, driving force amplitude and frequency is investigated.Comment: 4 pages, 2 figure

Quasi-static cracks and minimal energy surfaces

We compare the roughness of minimal energy(ME) surfaces and scalar ``quasi-static'' fracture surfaces(SQF). Two dimensional ME and SQF surfaces have the same roughness scaling, w sim L^zeta (L is system size) with zeta = 2/3. The 3-d ME and SQF results at strong disorder are consistent with the random-bond Ising exponent zeta (d >= 3) approx 0.21(5-d) (d is bulk dimension). However 3-d SQF surfaces are rougher than ME ones due to a larger prefactor. ME surfaces undergo a ``weakly rough'' to ``algebraically rough'' transition in 3-d, suggesting a similar behavior in fracture.Comment: 7 pages, aps.sty-latex, 7 figure

Scaling of interfaces in brittle fracture and perfect plasticity

The roughness properties of two-dimensional fracture surfaces as created by the slow failure of random fuse networks are considered and compared to yield surfaces of perfect plasticity with similar disorder. By studying systems up to a linear size L=350 it is found that in the cases studied the fracture surfaces exhibit self-affine scaling with a roughness exponent close to 2/3, which is asymptotically exactly true for plasticity though finite-size effects are evident for both. The overlap of yield or minimum energy and fracture surfaces with exactly the same disorder configuration is shown to be a decreasing function of the system size and to be of a rather large magnitude for all cases studied. The typical ``overlap cluster'' length between pairs of such interfaces converges to a constant with $L$ increasing.Comment: Accepted for publication in Phys. Rev.

Breakdown of Simple Scaling in Abelian Sandpile Models in One Dimension

We study the abelian sandpile model on decorated one dimensional chains. We determine the structure and the asymptotic form of distribution of avalanche-sizes in these models, and show that these differ qualitatively from the behavior on a simple linear chain. We find that the probability distribution of the total number of topplings $s$ on a finite system of size $L$ is not described by a simple finite size scaling form, but by a linear combination of two simple scaling forms $Prob_L(s) = 1/L f_1(s/L) + 1/L^2 f_2(s/L^2)$, for large $L$, where $f_1$ and $f_2$ are some scaling functions of one argument.Comment: 10 pages, revtex, figures include

Interface Motion in Random Media at Finite Temperature

We have studied numerically the dynamics of a driven elastic interface in a random medium, focusing on the thermal rounding of the depinning transition and on the behavior in the $T=0$ pinned phase. Thermal effects are quantitatively more important than expected from simple dimensional estimates. For sufficient low temperature the creep velocity at a driving force equal to the $T=0$ depinning force exhibits a power-law dependence on $T$, in agreement with earlier theoretical and numerical predictions for CDW's. We have also examined the dynamics in the $T=0$ pinned phase resulting from slowly increasing the driving force towards threshold. The distribution of avalanche sizes $S_\|$ decays as $S_\|^{-1-\kappa}$, with $\kappa = 0.05\pm 0.05$, in agreement with recent theoretical predictions.Comment: harvmac.tex, 30 pages, including 9 figures, available upon request. SU-rm-94073

Ab-initio self-consistent Gorkov-Green's function calculations of semi-magic nuclei - I. Formalism at second order with a two-nucleon interaction

An ab-initio calculation scheme for finite nuclei based on self-consistent Green's functions in the Gorkov formalism is developed. It aims at describing properties of doubly-magic and semi-magic nuclei employing state-of-the-art microscopic nuclear interactions and explicitly treating pairing correlations through the breaking of U(1) symmetry associated with particle number conservation. The present paper introduces the formalism, necessary to undertake applications at (self-consistent) second-order using two-nucleon interactions, in a detailed and self-contained fashion. First applications of such a scheme will be reported soon in a forthcoming publication. Future works will extend the present scheme to include three-nucleon interactions and implement more advanced truncation schemes.Comment: 38 page