A charge self-consistent LDA plus DMFT study of the spectral properties of hexagonal NiS

The electronic structure and spectral properties of hexagonal NiS have been studied in the high temperature paramagnetic phase and low temperature anti-ferromagnetic phase. The calculations have been performed using charge self-consistent density-functional theory in local density approximation combined with dynamical mean-field theory (LDA+DMFT). The photoemission spectra (PES) and optical properties have been computed and compared with the experimental data. Our results show that the dynamical correlation effects are important to understand the spectral and optical properties of NiS. These effects have been analyzed in detail by means of the computed real and imaginary part of the self-energy.

Phase diagram of the half-filled ionic Hubbard model

We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity using dynamical mean-field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered ionic potential Delta and the on-site Hubbard U. We find that for a finite Delta and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold U = U-AF via a first-order phase transition. For U smaller than U-AF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing U. We show that the results obtained within both methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature T-AF(U,Delta) or, equivalently, on decreasing U below U-AF(T,Delta)], within both methods, for weak to intermediate values of U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite temperature T, DMFT + CTQMC shows a second phase transition (not seen within DMFT + IPT) on increasing U beyond U-AF. At U-N > U-AF, when the Neel temperature T-N for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second-order transition. For U >> Delta, T-N similar to t(2)/U(1 - x(2)), where x = 2 Delta/U and thus T-N increases with increase in Delta/U. In the three-dimensional parameter space of (U/t, T/t, and Delta/t), as T increases, the surface of first-order transition at U-AF(T,Delta) and that of the second-order transition at U-N(T,Delta) approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.

Detection of a quantum particle on a lattice under repeated projective measurements

We consider a quantum particle, moving on a lattice with a tight-binding Hamiltonian, which is subjected to measurements to detect its arrival at a particular chosen set of sites. The projective measurements are made at regular time intervals tau, and we consider the evolution of the wave function until the time a detection occurs. We study the probabilities of its first detection at some time and, conversely, the probability of it not being detected (i.e., surviving) up to that time. We propose a general perturbative approach for understanding the dynamics which maps the evolution operator, which consists of unitary transformations followed by projections, to one described by a non-Hermitian Hamiltonian. For some examples of a particle moving on one-and two-dimensional lattices with one or more detection sites, we use this approach to find exact expressions for the survival probability and find excellent agreement with direct numerical results. A mean-field model with hopping between all pairs of sites and detection at one site is solved exactly. For the one-and two-dimensional systems, the survival probability is shown to have a power-law decay with time, where the power depends on the initial position of the particle. Finally, we show an interesting and nontrivial connection between the dynamics of the particle in our model and the evolution of a particle under a non-Hermitian Hamiltonian with a large absorbing potential at some sites.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Pattern-Formation During Self-organization in Three-Phase Eutectic Solidification

Eutectic growth is an interesting example for exploring the topic of pattern-formation in multi-phase systems, where the growth of the phases is coupled with the diffusive transport of one or more components in the melt. While in the case of binary alloys, the number of possibilities are limited (lamellae, rods, labyrinth etc.), their number rapidly increases with the number of components and phases. In this paper, we will investigate pattern formation during three-phase eutectic solidification using a state-of-the art phase-field method based on the grand-canonical density formulation. The major aim of the study is to highlight the role of two properties, which are the volume fraction of the solid phases and the solid-liquid interfacial energies, in the self-organization of the solid phases during directional growth. Thereafter, we will show representative phase-field simulations of a micro-structure in a real alloy (Ag-Al-Cu) using an asymmetric phase diagram as well as interfacial properties.

Evidence for 3D isotropic long range spin-spin interaction near the ferromagnetic transition in bulk and thin film SrRuO3

In the case of metallic ferromagnets there has always been a controversy, i.e. whether the magnetic interaction is itinerant or localized. For example SrRuO3 is known to be an itinerant ferromagnet where the spin-spin interaction is expected to be mean field in nature. However, it is reported to behave like Ising, Heisenberg or mean field by different groups. Despite several theoretical and experimental studies and the importance of strongly correlated systems, the experimental conclusion regarding the type of spin-spin interaction in SrRuO3 is lacking. To resolve this issue, we have investigated the critical behaviour in the vicinity of the paramagnetic-ferromagnetic phase transition using various techniques on polycrystalline as well as (001) oriented SrRuO3 films. Our analysis reveals that the application of a scaling law in the field-cooled magnetization data extracts the value of the critical exponent only when it is measured at H -> 0. To substantiate the actual nature without any ambiguity, the critical behavior is studied across the phase transition using the modified Arrott plot, Kouvel-Fisher plot and M-H isotherms. The critical analysis yields self-consistent beta, gamma and delta values and the spin interaction follows the long-range mean field model. Further the directional dependence of the critical exponent is studied in thin films and it reveals the isotropic nature. It is elucidated that the different experimental protocols followed by different groups are the reason for the ambiguity in determining the critical exponents in SrRuO3.

Evidence for 3D isotropic long range spin-spin interaction near the ferromagnetic transition in bulk and thin film SrRuO3

In the case of metallic ferromagnets there has always been a controversy, i.e. whether the magnetic interaction is itinerant or localized. For example SrRuO3 is known to be an itinerant ferromagnet where the spin-spin interaction is expected to be mean field in nature. However, it is reported to behave like Ising, Heisenberg or mean field by different groups. Despite several theoretical and experimental studies and the importance of strongly correlated systems, the experimental conclusion regarding the type of spin-spin interaction in SrRuO3 is lacking. To resolve this issue, we have investigated the critical behaviour in the vicinity of the paramagnetic-ferromagnetic phase transition using various techniques on polycrystalline as well as (001) oriented SrRuO3 films. Our analysis reveals that the application of a scaling law in the field-cooled magnetization data extracts the value of the critical exponent only when it is measured at H -> 0. To substantiate the actual nature without any ambiguity, the critical behavior is studied across the phase transition using the modified Arrott plot, Kouvel-Fisher plot and M-H isotherms. The critical analysis yields self-consistent beta, gamma and delta values and the spin interaction follows the long-range mean field model. Further the directional dependence of the critical exponent is studied in thin films and it reveals the isotropic nature. It is elucidated that the different experimental protocols followed by different groups are the reason for the ambiguity in determining the critical exponents in SrRuO3.

Complex patterns arise through spontaneous symmetry breaking in dense homogeneous networks of neural oscillators

There has been much interest in understanding collective dynamics in networks of brain regions due to their role in behavior and cognitive function. Here we show that a simple, homogeneous system of densely connected oscillators, representing the aggregate activity of local brain regions, can exhibit a rich variety of dynamical patterns emerging via spontaneous breaking of permutation or translational symmetries. Upon removing just a few connections, we observe a striking departure from the mean-field limit in terms of the collective dynamics, which implies that the sparsity of these networks may have very important consequences. Our results suggest that the origins of some of the complicated activity patterns seen in the brain may be understood even with simple connection topologies.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).

Weibull modulus for diverse strength due to sample-specificity

By sample specificity it is meant that specimens with the same nominal material parameters and tested under the same environmental conditions may exhibit different behavior with diversified strength. Such an effect has been widely observed in the testing of material failure and is usually attributed to the heterogeneity of material at the mesoscopic level. The degree with which mesoscopic heterogeneity affects macroscopic failure is still not clear. Recently, the problem has been examined by making use of statistical ensemble evolution of dynamical system and the mesoscopic stress re-distribution model (SRD). Sample specificity was observed for non-global mean stress field models, such as the duster mean field model, stress concentration at tip of microdamage, etc. Certain heterogeneity of microdamage could be sensitive to particular SRD leading to domino type of coalescence. Such an effect could start from the microdamage heterogeneity and then be magnified to other scale levels. This trans-scale sensitivity is the origin of sample specificity. The sample specificity leads to a failure probability Phi (N) with a transitional region 0 < (N) < 1, so that globally stable and catastrophic modes could co-exist. It is found that the scatter in strength can fit the Weibull distribution very well. Hence, the Weibull modulus is indicative of sample specificity. Numerical results obtained from the SRD for different non-global mean stress fields show that Weibull modulus increases with increasing sample span and influence region of microdamage.

A backward particle interpretation of Feynman-Kac formulae

We design a particle interpretation of Feynman-Kac measures on path spaces based on a backward Markovian representation combined with a traditional mean field particle interpretation of the flow of their final time marginals. In contrast to traditional genealogical tree based models, these new particle algorithms can be used to compute normalized additive functionals "on-the-fly" as well as their limiting occupation measures with a given precision degree that does not depend on the final time horizon. We provide uniform convergence results with respect to the time horizon parameter as well as functional central limit theorems and exponential concentration estimates. Our results have important consequences for online parameter estimation for non-linear non-Gaussian state-space models. We show how the forward filtering backward smoothing estimates of additive functionals can be computed using a forward only recursion.

Relationship between strain localization and catastrophic rupture

In order to explore a prior warning to catastrophic rupture of heterogeneous media, like rocks, the present study investigates the relationship between surface strain localization and catastrophic rupture. Instrumented observations on the evolution of surface strain field and the catastrophic rupture of a rock under uniaxial compression were carried out. It is found that the evolution of surface strain field displays two phases: at the early stage, the strain field keeps nearly uniform with weak fluctuations increasing slowly; but at the stage prior to catastrophic rupture, a certain accelerating localization develops and a localized zone emerges. Based on the measurements, an analysis was performed with local mean-field approximation. More importantly, it is found that the scale of localized zone is closely related to the catastrophic rupture strain and the rupture strain can be calculated in accord with the local-mean-field model satisfactorily. This provides a possible clue to the forecast of catastrophic rupture. (c) 2007 Elsevier Ltd. All rights reserved.

Elastic and statistically brittle (ESB) model, damage localization and catastrophic rupture

In this paper, an elastic and statistically brittle (ESB) model is applied to the process of damage evolution induced catastrophic rupture and the influence of localization and softening on catastrophic rupture is discussed. According to the analysis, the uncertainty of catastrophic rupture should be attributed to the unknown scale of localized zone. Based on the elastic and statistically brittle model but local mean field approximation, the relation between the scale of localized zone and catastrophic rupture is obtained and then justified with experiments. These results can not only give a deeper understanding of the mechanism governing catastrophic rupture, but also provide a possible tool to foresee the occurrence of catastrophic rupture.