Random path sampling applied to vector models of magnetism

The recently reported Monte Carlo Random Path Sampling method (RPS) is here improved and its application is expanded to the study of the 2D and 3D Ising and discrete Heisenberg models. The methodology was implemented to allow use in both CPU-based high-performance computing infrastructures (C/MPI) and GPU-based (CUDA) parallel computation, with significant computational performance gains. Convergence is discussed, both in terms of free energy and magnetization dependence on field/temperature. From the calculated magnetization-energy joint density of states, fast calculations of field and temperature dependent thermodynamic properties are performed, including the effects of anisotropy on coercivity, and the magnetocaloric effect. The emergence of first-order magneto-volume transitions in the compressible Ising model is interpreted using the Landau theory of phase transitions. Using metallic Gadolinium as a real-world example, the possibility of using RPS as a tool for computational magnetic materials design is discussed. Experimental magnetic and structural properties of a Gadolinium single crystal are compared to RPS-based calculations using microscopic parameters obtained from Density Functional Theory.

Observations of turbulence within the surf and swash zone of a field-scale sandy laboratory beach

Current coastal-evolution models generally lack the ability to accurately predict bed level change in shallow (<~2 m) water, which is, at least partly, due to the preclusion of the effect of surface-induced turbulence on sand suspension and transport. As a first step to remedy this situation, we investigated the vertical structure of turbulence in the surf and swash zone using measurements collected under random shoaling and plunging waves on a steep (initially 1:15) field-scale sandy laboratory beach. Seaward of the swash zone, turbulence was measured with a vertical array of three Acoustic Doppler Velocimeters (ADVs), while in the swash zone two vertically spaced acoustic doppler velocimeter profilers (Vectrino profilers) were applied. The vertical turbulence structure evolves from bottom-dominated to approximately vertically uniform with an increase in the fraction of breaking waves to ~ 50%. In the swash zone, the turbulence is predominantly bottom-induced during the backwash and shows a homogeneous turbulence profile during uprush. We further find that the instantaneous turbulence kinetic energy is phase-coupled with the short-wave orbital motion under the plunging breakers, with higher levels shortly after the reversal from offshore to onshore motion (i.e. wavefront).

Undersampled Critical Branching Processes on Small-World and Random Networks Fail to Reproduce the Statistics of Spike Avalanches

The power-law size distributions obtained experimentally for neuronal avalanches are an important evidence of criticality in the brain. This evidence is supported by the fact that a critical branching process exhibits the same exponent t~3=2. Models at criticality have been employed to mimic avalanche propagation and explain the statistics observed experimentally. However, a crucial aspect of neuronal recordings has been almost completely neglected in the models: undersampling. While in a typical multielectrode array hundreds of neurons are recorded, in the same area of neuronal tissue tens of thousands of neurons can be found. Here we investigate the consequences of undersampling in models with three different topologies (two-dimensional, small-world and random network) and three different dynamical regimes (subcritical, critical and supercritical). We found that undersampling modifies avalanche size distributions, extinguishing the power laws observed in critical systems. Distributions from subcritical systems are also modified, but the shape of the undersampled distributions is more similar to that of a fully sampled system. Undersampled supercritical systems can recover the general characteristics of the fully sampled version, provided that enough neurons are measured. Undersampling in two-dimensional and small-world networks leads to similar effects, while the random network is insensitive to sampling density due to the lack of a well-defined neighborhood. We conjecture that neuronal avalanches recorded from local field potentials avoid undersampling effects due to the nature of this signal, but the same does not hold for spike avalanches. We conclude that undersampled branching-process-like models in these topologies fail to reproduce the statistics of spike avalanches.

Random matrix theory eta aplikazioa kaos kuantikoan

[eus] Gradu amaierako lan honetan ausazko matrizeen teoriari, RMT-ri, buruzko sarrera orokor bat egiten da ondoren aplikazio fisiko bat emateko. Teoriaren aplikazioa egiteko Kaos kuantikoa deritzon fisikaren arloa erabiliko da. Lehenik eta behin, RMT-ren kontzeptu batzuk azalduko dira helburutzat lehen auzokideen distantziaren distribuzioaren espresio lortzea izanik. Izan ere, distribuzio honek erakutsiko baititu Kaosak kuantikoki uzten dituen aztarnak. Bigarren kapituluan, aplikazio fisikoa azalduko da. Lehenengo Kaosean RMT nola aplikatzen den ikusiko da, ondoren adibide batzuen bidez argituz, eremu magnetiko batean dagoen hidrogeno atomoa eta billar kuantikoak izenarekin ezagutzen diren sistemak, batik bat.

Conditional mutual inclusive information enables accurate quantification of associations in gene regulatory networks

International audience

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.