Critical clusters and efficient dynamics for frustrated spin models

A general method to find, in a systematic way, efficient Monte Carlo cluster dynamics among the avast class of dynamics introduced by Kandel et al. [Phys. Rev. Lett. 65, 941 (1990)] is proposed. The method is successfully applied to a class of frustrated two-dimensional Ising systems. In the case of the fully frustrated model, we also find the intriguing result that critical clusters consist of self-avoiding walk at the theta point.

Intramolecular coupling as a mechanism for a liquid-liquid phase transition

We study a model for water with a tunable intramolecular interaction Js, using mean-field theory and off-lattice Monte Carlo simulations. For all Js>~0, the model displays a temperature of maximum density. For a finite intramolecular interaction Js>0, our calculations support the presence of a liquid-liquid phase transition with a possible liquid-liquid critical point for water, likely preempted by inevitable freezing. For J=0, the liquid-liquid critical point disappears at T=0.

Precursor phenomena in frustrated systems

To understand the origin of the dynamical transition, between high-temperature exponential relaxation and low-temperature nonexponential relaxation, that occurs well above the static transition in glassy systems, a frustrated spin model, with and without disorder, is considered. The model has two phase transitions, the lower being a standard spin glass transition (in the presence of disorder) or fully frustrated Ising (in the absence of disorder), and the higher being a Potts transition. Monte Carlo results clarify that in the model with (or without) disorder the precursor phenomena are related to the Griffiths (or Potts) transition. The Griffiths transition is a vanishing transition which occurs above the Potts transition and is present only when disorder is present, while the Potts transition which signals the effect due to frustration is always present. These results suggest that precursor phenomena in frustrated systems are due either to disorder and/or to frustration, giving a consistent interpretation also for the limiting cases of Ising spin glass and of Ising fully frustrated model, where also the Potts transition is vanishing. This interpretation could play a relevant role in glassy systems beyond the spin systems case.

Essays on the treatment of cash flows under stochastic interest rates

Abstract This paper shows how to calculate recursively the moments of the accumulated and discounted value of cash flows when the instantaneous rates of return follow a conditional ARMA process with normally distributed innovations. We investigate various moment based approaches to approximate the distribution of the accumulated value of cash flows and we assess their performance through stochastic Monte-Carlo simulations. We discuss the potential use in insurance and especially in the context of Asset-Liability Management of pension funds.

Mean first-passage times for systems driven by gamma and McFadden dichotomous noise

We consider mean first-passage times (MFPTs) for systems driven by non-Markov gamma and McFadden dichotomous noises. A simplified derivation is given of the underlying integral equations and the theory for ordinary renewal processes is extended to modified and equilibrium renewal processes. The exact results are compared with the MFPT for Markov dichotomous noise and with the results of Monte Carlo simulations.

Damage spreading transition in glasses: A probe for the ruggedness of the configurational landscape

We consider damage spreading transitions in the framework of mode-coupling theory. This theory describes relaxation processes in glasses in the mean-field approximation which are known to be characterized by the presence of an exponentially large number of metastable states. For systems evolving under identical but arbitrarily correlated noises, we demonstrate that there exists a critical temperature T0 which separates two different dynamical regimes depending on whether damage spreads or not in the asymptotic long-time limit. This transition exists for generic noise correlations such that the zero damage solution is stable at high temperatures, being minimal for maximal noise correlations. Although this dynamical transition depends on the type of noise correlations, we show that the asymptotic damage has the good properties of a dynamical order parameter, such as (i) independence of the initial damage; (ii) independence of the class of initial condition; and (iii) stability of the transition in the presence of asymmetric interactions which violate detailed balance. For maximally correlated noises we suggest that damage spreading occurs due to the presence of a divergent number of saddle points (as well as metastable states) in the thermodynamic limit consequence of the ruggedness of the free-energy landscape which characterizes the glassy state. These results are then compared to extensive numerical simulations of a mean-field glass model (the Bernasconi model) with Monte Carlo heat-bath dynamics. The freedom of choosing arbitrary noise correlations for Langevin dynamics makes damage spreading an interesting tool to probe the ruggedness of the configurational landscape.

Vortex state and effect of anisotropy in sub-100 nm magnetic nanodots

Magnetic properties of Fe nanodots are simulated using a scaling technique and Monte Carlo method, in good agreement with experimental results. For the 20-nm-thick dots with diameters larger than 60¿nm, the magnetization reversal via vortex state is observed. The role of magnetic interaction between dots in arrays in the reversal process is studied as a function of nanometric center-to-center distance. When this distance is more than twice the dot diameter, the interaction can be neglected and the magnetic properties of the entire array are determined by the magnetic configuration of the individual dots. The effect of crystalline anisotropy on the vortex state is investigated. For arrays of noninteracting dots, the anisotropy strongly affects the vortex nucleation field and coercivity, and only slightly affects the vortex annihilation field

Effective squark/chargino/neutralino couplings: MadGraph implementation

We have included the effective description of squark interactions with charginos/neutralinos in the MadGraph MSSM model. This effective description includes the effective Yukawa couplings, and another logarithmic term which encodes the supersymmetry-breaking. We have performed an extensive test of our implementation analyzing the results of the partial decay widths of squarks into charginos and neutralinos obtained by using FeynArts/FormCalc programs and the new model file in MadGraph. We present results for the cross-section of top-squark production decaying into charginos and neutralinos.

DnaSP, DNA polymorphism analyses by the coalescent and other methods

DnaSP is a software package for the analysis of DNA polymorphism data. Present version introduces several new modules and features which, among other options allow: (1) handling big data sets (~5 Mb per sequence); (2) conducting a large number of coalescent-based tests by Monte Carlo computer simulations; (3) extensive analyses of the genetic differentiation and gene flow among populations; (4) analysing the evolutionary pattern of preferred and unpreferred codons; (5) generating graphical outputs for an easy visualization of results. Availability: The software package, including complete documentation and examples, is freely available to academic users from: http://www.ub.es/dnasp

Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered.

Linear least squares estimation of the first order moving average parameter

We propose an iterative procedure to minimize the sum of squares function which avoids the nonlinear nature of estimating the first order moving average parameter and provides a closed form of the estimator. The asymptotic properties of the method are discussed and the consistency of the linear least squares estimator is proved for the invertible case. We perform various Monte Carlo experiments in order to compare the sample properties of the linear least squares estimator with its nonlinear counterpart for the conditional and unconditional cases. Some examples are also discussed

Morphological changes in convection-diffusion-limited deposition

The effect of hydrodynamic flow upon diffusion-limited deposition on a line is investigated using a Monte Carlo model. The growth process is governed by the convection and diffusion field. The convective diffusion field is simulated by the biased-random walker resulting from a superimposed drift that represents the convective flow. The development of distinct morphologies is found with varying direction and strength of drift. By introducing a horizontal drift parallel to the deposition plate, the diffusion-limited deposit changes into a single needle inclined to the plate. The width of the needle decreases with increasing strength of drift. The angle between the needle and the plate is about 45° at high flow rate. In the presence of an inclined drift to the plate, the convection-diffusion-limited deposit leads to the formation of a characteristic columnar morphology. In the limiting case where the convection dominates, the deposition process is equivalent to ballistic deposition onto an inclined surface.

Diffusion of passive scalars under stochastic convection

The diffusion of passive scalars convected by turbulent flows is addressed here. A practical procedure to obtain stochastic velocity fields with well¿defined energy spectrum functions is also presented. Analytical results are derived, based on the use of stochastic differential equations, where the basic hypothesis involved refers to a rapidly decaying turbulence. These predictions are favorable compared with direct computer simulations of stochastic differential equations containing multiplicative space¿time correlated noise.