-
Demystifying Spatial Confounding
Authors:
Emiko Dupont,
Isa Marques,
Thomas Kneib
Abstract:
Spatial confounding is a fundamental issue in regression models for spatially indexed data. It arises because spatial random effects, included to approximate unmeasured spatial variation, are typically not independent of the covariates in the model. This can lead to significant bias in covariate effect estimates. Despite extensive research, it is still a topic of much confusion with sometimes puzz…
▽ More
Spatial confounding is a fundamental issue in regression models for spatially indexed data. It arises because spatial random effects, included to approximate unmeasured spatial variation, are typically not independent of the covariates in the model. This can lead to significant bias in covariate effect estimates. Despite extensive research, it is still a topic of much confusion with sometimes puzzling and seemingly contradictory results. In this paper we develop a broad theoretical framework that brings mathematical clarity to the mechanisms of spatial confounding, providing explicit and interpretable analytical expressions for the resulting bias. From these, we see that it is a problem directly linked to spatial smoothing, and we can identify exactly how the features of the model and the data generation process affect the size and occurrence of bias. We also use our framework to understand and generalise some of the main results on spatial confounding in the past, including suggested methods for bias adjustment. Thus, our comprehensive and mathematically explicit approach clears up existing confusion and, indeed, demystifies the issue of spatial confounding.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
Bayesian spatial+: A joint model perspective
Authors:
Isa Marques,
Paul F. V. Wiemann
Abstract:
A common phenomenon in spatial regression models is spatial confounding. This phenomenon occurs when spatially indexed covariates modeling the mean of the response are correlated with a spatial effect included in the model. spatial+ Dupont et al. (2022) is a popular approach to reducing spatial confounding. spatial+ is a two-stage frequentist approach that explicitly models the spatial structure i…
▽ More
A common phenomenon in spatial regression models is spatial confounding. This phenomenon occurs when spatially indexed covariates modeling the mean of the response are correlated with a spatial effect included in the model. spatial+ Dupont et al. (2022) is a popular approach to reducing spatial confounding. spatial+ is a two-stage frequentist approach that explicitly models the spatial structure in the confounded covariate, removes it, and uses the corresponding residuals in the second stage. In a frequentist setting, there is no uncertainty propagation from the first stage estimation determining the residuals since only point estimates are used. Inference can also be cumbersome in a frequentist setting, and some of the gaps in the original approach can easily be remedied in a Bayesian framework. First, a Bayesian joint model can easily achieve uncertainty propagation from the first to the second stage of the model. In a Bayesian framework, we also have the tools to infer the model's parameters directly. Notably, another advantage of using a Bayesian framework we thoroughly explore is the ability to use prior information to impose restrictions on the spatial effects rather than applying them directly to their posterior. We build a joint prior for the smoothness of all spatial effects that simultaneously shrinks towards a high smoothness of the response and imposes that the spatial effect in the response is a smoother of the confounded covariates' spatial effect. This prevents the response from operating at a smaller scale than the covariate and can help to avoid situations where there is insufficient variation in the residuals resulting from the first stage model. We evaluate the performance of the Bayesian spatial+ via both simulated and real datasets.
△ Less
Submitted 11 September, 2023;
originally announced September 2023.
-
Roadmap for Optical Tweezers
Authors:
Giovanni Volpe,
Onofrio M. Maragò,
Halina Rubinzstein-Dunlop,
Giuseppe Pesce,
Alexander B. Stilgoe,
Giorgio Volpe,
Georgiy Tkachenko,
Viet Giang Truong,
Síle Nic Chormaic,
Fatemeh Kalantarifard,
Parviz Elahi,
Mikael Käll,
Agnese Callegari,
Manuel I. Marqués,
Antonio A. R. Neves,
Wendel L. Moreira,
Adriana Fontes,
Carlos L. Cesar,
Rosalba Saija,
Abir Saidi,
Paul Beck,
Jörg S. Eismann,
Peter Banzer,
Thales F. D. Fernandes,
Francesco Pedaci
, et al. (58 additional authors not shown)
Abstract:
Optical tweezers are tools made of light that enable contactless pushing, trapping, and manipulation of objects ranging from atoms to space light sails. Since the pioneering work by Arthur Ashkin in the 1970s, optical tweezers have evolved into sophisticated instruments and have been employed in a broad range of applications in life sciences, physics, and engineering. These include accurate force…
▽ More
Optical tweezers are tools made of light that enable contactless pushing, trapping, and manipulation of objects ranging from atoms to space light sails. Since the pioneering work by Arthur Ashkin in the 1970s, optical tweezers have evolved into sophisticated instruments and have been employed in a broad range of applications in life sciences, physics, and engineering. These include accurate force and torque measurement at the femtonewton level, microrheology of complex fluids, single micro- and nanoparticle spectroscopy, single-cell analysis, and statistical-physics experiments. This roadmap provides insights into current investigations involving optical forces and optical tweezers from their theoretical foundations to designs and setups. It also offers perspectives for applications to a wide range of research fields, from biophysics to space exploration.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
A multivariate Gaussian random field prior against spatial confounding
Authors:
Isa Marques,
Thomas Kneib,
Nadja Klein
Abstract:
Spatial models are used in a variety research areas, such as environmental sciences, epidemiology, or physics. A common phenomenon in many spatial regression models is spatial confounding. This phenomenon takes place when spatially indexed covariates modeling the mean of the response are correlated with the spatial random effect. As a result, estimates for regression coefficients of the covariates…
▽ More
Spatial models are used in a variety research areas, such as environmental sciences, epidemiology, or physics. A common phenomenon in many spatial regression models is spatial confounding. This phenomenon takes place when spatially indexed covariates modeling the mean of the response are correlated with the spatial random effect. As a result, estimates for regression coefficients of the covariates can be severely biased and interpretation of these is no longer valid. Recent literature has shown that typical solutions for reducing spatial confounding can lead to misleading and counterintuitive results. In this paper, we develop a computationally efficient spatial model in a Bayesian framework integrating novel prior structure that reduces spatial confounding. Starting from the univariate case, we extend our prior structure to case of multiple spatially confounded covariates. In a simulation study, we show that our novel model flexibly detects and reduces spatial confounding in spatial datasets, and it performs better than typically used methods such as restricted spatial regression. These results are promising for any applied researcher who wishes to interpret covariate effects in spatial regression models. As a real data illustration, we study the effect of elevation and temperature on the mean of daily precipitation in Germany.
△ Less
Submitted 7 June, 2021;
originally announced June 2021.
-
Multipole Engineering of Attractive-Repulsive and Bending Optical Forces
Authors:
Denis A. Kislov,
Egor A. Gurvitz,
Alexander A. Pavlov,
Dmitrii N. Redka,
Manuel I. Marques,
Pavel Ginzburg,
Alexander S. Shalin
Abstract:
Focused laser beams allow controlling mechanical motion of objects and can serve as a tool for assembling complex micro and nano structures in space. While in a vast majority of cases small particles experience attractive gradient forces and repulsive radiation pressure, introduction of additional degrees of freedom into optomechanical manipulation suggests approaching new capabilities. Here we an…
▽ More
Focused laser beams allow controlling mechanical motion of objects and can serve as a tool for assembling complex micro and nano structures in space. While in a vast majority of cases small particles experience attractive gradient forces and repulsive radiation pressure, introduction of additional degrees of freedom into optomechanical manipulation suggests approaching new capabilities. Here we analyze optical forces acting on a high refractive index silicon sphere in a focused Gaussian beam and reveal new regimes of particles anti-trapping. Multipolar analysis allows separating an optical force into interception and recoil components, which have a completely different physical nature resulting in different mechanical actions. In particular, interplaying interception radial forces and multipolar resonances within a particle can lead to either trapping or anti-trapping scenarios, depending of the overall system parameters. At the same time, the recoil force generates a significant azimuthal component along with an angular-dependent radial force. Those contribution enable enhancing either trapping or anti-trapping regimes and also introduce bending reactions. These effects are linked to the far-field multipole interference resulting and, specifically, to its asymmetric scattering diagrams. The later approach is extremely useful, as it allows assessing the nature of optomechanical motion by observing far-field patterns. Multipolar engineering of optical forces, being quite general approach, is not necessarily linked to simple spherical shapes and paves a way to new possibilities in microfluidic applications, including sorting and micro assembly of nontrivial volumetric geometries.
△ Less
Submitted 18 February, 2021;
originally announced February 2021.
-
Light induced "Mock Gravity" at the nanoscale
Authors:
Jorge Luis-Hita,
Manuel I. Marques,
Rafael Delgado-Buscalioni,
Nuno de Sousa,
Luis S. Froufe-Perez,
Frank Scheffold,
Juan Jose Saenz
Abstract:
The origin of long-range attractive interactions has fascinated scientist along centuries. The remarkable Fatio-LeSage's corpuscular theory, introduced as early as in 1690 and generalized to electromagnetic waves by Lorentz, proposed that, due to their mutual shadowing, two absorbing particles in an isotropic radiation field experience an attractive force which follows a gravity-like inverse squar…
▽ More
The origin of long-range attractive interactions has fascinated scientist along centuries. The remarkable Fatio-LeSage's corpuscular theory, introduced as early as in 1690 and generalized to electromagnetic waves by Lorentz, proposed that, due to their mutual shadowing, two absorbing particles in an isotropic radiation field experience an attractive force which follows a gravity-like inverse square distance law. Similar "Mock Gravity" interactions were later introduced by Spitzer and Gamow in the context of Galaxy formation but their actual relevance in Cosmology has never been unambiguously established. Here we predict the existence of Mock-Gravity, inverse square distance, attractive forces between two identical molecules or nanoparticles in a quasi monochromatic isotropic random light field, whenever the light frequency is tuned to an absorption line such that the real part of the particle's electric polarizability is zero, i.e. at the so-called Froehlich resonance. These interactions are scale independent, holding for both near and far-field separation distances.
△ Less
Submitted 15 February, 2018;
originally announced February 2018.
-
Strong collective currents of gold nanoparticles in an optical vortex lattice
Authors:
R. Delgado-Buscalioni,
M. Meléndez,
J. Luis-Hita,
M. I. Marqués,
J. J. Sáenz
Abstract:
Self-organized collective behaviour of active units is inspiring new designs of artificial swarms of micron-sized objects. However, active control at the nanoscale remains elusive. We have accurately solved the collective optofluidic dynamics of gold plasmonic 50 nm-radius nanoparticles moving in aqueous solution under a nonconservative optical vortex lattice. Above a critical field intensity and…
▽ More
Self-organized collective behaviour of active units is inspiring new designs of artificial swarms of micron-sized objects. However, active control at the nanoscale remains elusive. We have accurately solved the collective optofluidic dynamics of gold plasmonic 50 nm-radius nanoparticles moving in aqueous solution under a nonconservative optical vortex lattice. Above a critical field intensity and concentration, the interplay between optical forces, thermal fluctuations and hydrodynamic pairing leads to a spontaneous transition towards synchronised motion exhibiting a rich assortment of collective dynamics. The critical phenomenon creates strong unidirectional currents of nanoparticles, reaching speeds of centimeters per second along the diagonals of the lattice.
△ Less
Submitted 29 September, 2017; v1 submitted 13 September, 2017;
originally announced September 2017.
-
Effect of long range spatial correlations on the lifetime statistics of an emitter in a two-dimensional disordered lattice
Authors:
N. de Sousa,
J. J. Sáenz,
A. García-Martín,
L. S. Froufe-Pérez,
M. I. Marqués
Abstract:
The effect of spatial correlations on the Purcell effect in a bidimensional dispersion of resonant nanoparticles is analyzed. We perform extensive calculations on the fluorescence decay rate of a point emitter embedded in a system of nanoparticles statistically distributed according to a simple 2D lattice-gas model near the critical point. For short-range correlations (high temperature thermalizat…
▽ More
The effect of spatial correlations on the Purcell effect in a bidimensional dispersion of resonant nanoparticles is analyzed. We perform extensive calculations on the fluorescence decay rate of a point emitter embedded in a system of nanoparticles statistically distributed according to a simple 2D lattice-gas model near the critical point. For short-range correlations (high temperature thermalization) the Purcell factors present a long-tailed statistic which evolves towards a bimodal distribution when approaching the critical point where the spatial correlation length diverges. Our results suggest long-range correlations as a possible origin of the large fluctuations of experimental decay rates in disordered metal films.
△ Less
Submitted 10 June, 2014; v1 submitted 16 January, 2014;
originally announced January 2014.
-
A model for vortex formation in magnetic nanodots
Authors:
S. A. Leonel,
I. A. Marques,
P. Z. Coura,
B. V. Costa
Abstract:
We use Monte Carlo simulation to study the vortex nucleation on magnetic nanodots at low temperature. In our simulations, we have considered a simple microscopic two-dimensional anisotropic Heisenberg model with term to describe the anisotropy due to the presence of the nanodot edge. We have considered the thickness of the edge, which was not considered in previous works, introducing a term that c…
▽ More
We use Monte Carlo simulation to study the vortex nucleation on magnetic nanodots at low temperature. In our simulations, we have considered a simple microscopic two-dimensional anisotropic Heisenberg model with term to describe the anisotropy due to the presence of the nanodot edge. We have considered the thickness of the edge, which was not considered in previous works, introducing a term that controls the energy associated to the edge. Our results clearly show that the thickness of the edge has a considerable influence in the vortex nucleation on magnetic nanodots. We have obtained the hysteresis curve for several values of the surface anisotropy and skin depth parameter ($ξ$). The results are in excellent agreement with experimental data.
△ Less
Submitted 21 October, 2010;
originally announced October 2010.
-
A high pressure calorimetric experiment to validate the liquid-liquid critical point hypothesis in water
Authors:
Manuel I. Marques
Abstract:
An experimental proposal to test the existence of a liquid-liquid critical point in water, based on high pressure calorimetric measurements, is presented on this paper. Considering the existence of an intramolecular correlation in the water molecule we show how the response of the specific heat at high pressure is different depending on the existence, or not, of the second critical point. If the…
▽ More
An experimental proposal to test the existence of a liquid-liquid critical point in water, based on high pressure calorimetric measurements, is presented on this paper. Considering the existence of an intramolecular correlation in the water molecule we show how the response of the specific heat at high pressure is different depending on the existence, or not, of the second critical point. If the liquid-liquid critical point hypothesis is true there must be a maximum in the specific heat at some temperature $T>T_{H}$ for any pressure $P>P_{c}$ (being $T_{H}$ the homogeneous nucleation temperature and $P_{c}$ the pressure of the second critical point). This maximum does not appear for the singularity free scenario.
△ Less
Submitted 23 January, 2007;
originally announced January 2007.
-
Relationship between Fragility, Diffusive Directions and Energy Barriers in a Supercooled Liquid
Authors:
Manuel I. Marques,
H. Eugene Stanley
Abstract:
An analysis of diffusion in a supercooled liquid based solely in the density of diffusive directions and the value of energy barriers shows how the potential energy landscape (PEL) approach is capable of explaining the $α$ and $β$ relaxations and the fragility of a glassy system. We find that the $β$ relaxation is directly related to the search for diffusive directions. Our analysis shows how in…
▽ More
An analysis of diffusion in a supercooled liquid based solely in the density of diffusive directions and the value of energy barriers shows how the potential energy landscape (PEL) approach is capable of explaining the $α$ and $β$ relaxations and the fragility of a glassy system. We find that the $β$ relaxation is directly related to the search for diffusive directions. Our analysis shows how in strong liquids diffusion is mainly energy activated, and how in fragile liquids the diffusion is governed by the density of diffusive directions. We describe the fragile-to-strong crossover as a change in the topography of the PEL sampled by the system at a certain crossover temperature $T_\times$.
△ Less
Submitted 28 March, 2006;
originally announced March 2006.
-
Behavior of the local mode potential in BaTiO3 studied by Effective Hamiltonian numerical simulations
Authors:
Manuel I. Marques,
Julio A. Gonzalo
Abstract:
The behavior of the local mode potential energy in the unit cell of BaTiO3 is studied by effective Hamiltonian numerical simulations. The study is focused on the evolution of the potential minimum felt by the Ti ion. At the paraelectric phase, the single instantaneous minimum is located at one of eight off-centered regions at the body diagonals of the unit cell. The study shows graphically how t…
▽ More
The behavior of the local mode potential energy in the unit cell of BaTiO3 is studied by effective Hamiltonian numerical simulations. The study is focused on the evolution of the potential minimum felt by the Ti ion. At the paraelectric phase, the single instantaneous minimum is located at one of eight off-centered regions at the body diagonals of the unit cell. The study shows graphically how the soft mode behavior of the perovskite comes from the average response of the Ti ion to the potential changes in the unit cell, showing that the single minimum of the time-averaged potential is centered at the paraelectric phase and off-centered at the ferroelectric phase.
△ Less
Submitted 16 March, 2006;
originally announced March 2006.
-
Poling effect on distribution of quenched random fields in a uniaxial relaxor ferroelectric
Authors:
Manuel I. Marques,
Carmen Arago
Abstract:
The frequency dependence of the dielectric permitivity's maximum has been studied for poled and unpoled doped relaxor strontium barium niobate $Sr_{0.61}Ba_{0.39}Nb_{2}O_{6}:Cr^{3+}$ (SBN-61:Cr). In both cases the maximum found is broad and the frequency dispersion is strong. The present view of random fields compensation in the unpoled sample is not suitable for explaining this experimental res…
▽ More
The frequency dependence of the dielectric permitivity's maximum has been studied for poled and unpoled doped relaxor strontium barium niobate $Sr_{0.61}Ba_{0.39}Nb_{2}O_{6}:Cr^{3+}$ (SBN-61:Cr). In both cases the maximum found is broad and the frequency dispersion is strong. The present view of random fields compensation in the unpoled sample is not suitable for explaining this experimental result. We propose a new mechanism where the dispersion of quenched random electric fields, affecting the nanodomains, is minimized after poling. We test our proposal by numerical simulations on a random field Ising model. Results obtained are in agreement with the polarization's measurements presented by Granzow et al. [Phys. Rev. Lett {\bf 92}, 065701 (2004)].
△ Less
Submitted 14 March, 2006;
originally announced March 2006.
-
The quantum paraelectric behavior of SrTiO_{3} revisited: relevance of the structural phase transition temperature
Authors:
Manuel I. Marques,
Carmen Arago,
Julio A. Gonzalo
Abstract:
It has been known for a long time that the low temperature behavior shown by the dielectric constant of quantum paraelectric $SrTiO_{3}$ can not be fitted properly by Barrett's formula using a single zero point energy or saturation temperature ($T_{1}$). As it was originally shown [K. A. Müller and H. Burkard, Phys. Rev. B {\bf 19}, 3593 (1979)] a crossover between two different saturation tempe…
▽ More
It has been known for a long time that the low temperature behavior shown by the dielectric constant of quantum paraelectric $SrTiO_{3}$ can not be fitted properly by Barrett's formula using a single zero point energy or saturation temperature ($T_{1}$). As it was originally shown [K. A. Müller and H. Burkard, Phys. Rev. B {\bf 19}, 3593 (1979)] a crossover between two different saturation temperatures ($T_{1l}$=77.8K and $T_{1h}$=80K) at $T\sim10K$ is needed to explain the low and high temperature behavior of the dielectric constant. However, the physical reason for the crossover between these two particular values of the saturation temperature at $T\sim10K$ is unknown. In this work we show that the crossover between these two values of the saturation temperature at $T\sim10K$ can be taken as a direct consequence of (i) the quantum distribution of frequencies $g(Ω)\proptoΩ^{2}$ associated with the complete set of low-lying modes and (ii) the existence of a definite maximum phonon frequency given by the structural transition critical temperature $T_{tr}$.
△ Less
Submitted 13 March, 2006;
originally announced March 2006.
-
Test of cold denaturation mechanism for proteins as a function of water's structure
Authors:
Manuel I. Marques
Abstract:
In a recent paper [PRL 91, 138103 (2003)] a new mechanism to explain the cold denaturation of proteins, based on the loss of local low-density water structure, has been proposed. In the present paper this mechanism is tested by means of full atom numerical simulations. In good agreement with this proposal, cold denaturation resulting in the unfolded state was found at the High Density Liquid (HD…
▽ More
In a recent paper [PRL 91, 138103 (2003)] a new mechanism to explain the cold denaturation of proteins, based on the loss of local low-density water structure, has been proposed. In the present paper this mechanism is tested by means of full atom numerical simulations. In good agreement with this proposal, cold denaturation resulting in the unfolded state was found at the High Density Liquid (HDL) state of water, at which the amount of open tetragonal hydrogen bonds decreases at cooling.
△ Less
Submitted 9 March, 2006;
originally announced March 2006.
-
Extinction of BKT transition by nonmagnetic disorder in planar-symmetry spin models
Authors:
G. M. Wysin,
A. R. Pereira,
I. A. Marques,
S. A. Leonel,
P. Z. Coura
Abstract:
The Berezinskii-Kosterlitz-Thouless (BKT) transition in two-dimensional planar rotator and XY models on a square lattice, diluted by randomly placed vacancies, is studied here using hybrid Monte Carlo simulations that combine single spin flip, cluster and over-relaxation techniques. The transition temperature $T_c$ is determined as a function of vacancy density $ρ_{vac}$ by calculations of the h…
▽ More
The Berezinskii-Kosterlitz-Thouless (BKT) transition in two-dimensional planar rotator and XY models on a square lattice, diluted by randomly placed vacancies, is studied here using hybrid Monte Carlo simulations that combine single spin flip, cluster and over-relaxation techniques. The transition temperature $T_c$ is determined as a function of vacancy density $ρ_{vac}$ by calculations of the helicity modulus and the by finite-size scaling of the in-plane magnetic susceptibility. The results for $T_c$ are consistent with those from the much less precise fourth-order cumulant of Binder. $T_c$ is found to decrease monotonically with increasing $ρ_{vac}$, and falls to zero close to the square lattice percolation limit, $ρ_{vac}\approx 0.41$ . The result is physically reasonable: the long-range orientational order of the low-temperature phase cannot be maintained in the absence of sufficient spin interactions across the lattice.
△ Less
Submitted 6 April, 2005;
originally announced April 2005.
-
A possible mechanism for cold denaturation of proteins at high pressure
Authors:
Manuel I. Marques,
Jose M. Borreguero,
H. Eugene Stanley,
Nikolay V. Dokholyan
Abstract:
We study cold denaturation of proteins at high pressures. Using multicanonical Monte Carlo simulations of a model protein in a water bath, we investigate the effect of water density fluctuations on protein stability. We find that above the pressure where water freezes to the dense ice phase ($\approx2$ kbar), the mechanism for cold denaturation with decreasing temperature is the loss of local lo…
▽ More
We study cold denaturation of proteins at high pressures. Using multicanonical Monte Carlo simulations of a model protein in a water bath, we investigate the effect of water density fluctuations on protein stability. We find that above the pressure where water freezes to the dense ice phase ($\approx2$ kbar), the mechanism for cold denaturation with decreasing temperature is the loss of local low-density water structure. We find our results in agreement with data of bovine pancreatic ribonuclease A.
△ Less
Submitted 11 June, 2003; v1 submitted 30 December, 2002;
originally announced December 2002.
-
Accurate Monte Carlo critical exponents for Ising lattices
Authors:
Jorge Garcia,
Julio A. Gonzalo,
Manuel I. Marques
Abstract:
A careful Monte Carlo investigation of the phase transition very close to the critical point (T -> Tc, H -> 0) in relatively large d = 3, s = 1/2 Ising lattices did produce critical exponents beta = 0.3126(4) =~ 5/16, delta^{-1} = 0.1997(4) =~ 1/5 and gamma_{3D} = 1.253(4) =~ 5/4. Our results indicate that, within experimental error, they are given by simple fractions corresponding to the linear…
▽ More
A careful Monte Carlo investigation of the phase transition very close to the critical point (T -> Tc, H -> 0) in relatively large d = 3, s = 1/2 Ising lattices did produce critical exponents beta = 0.3126(4) =~ 5/16, delta^{-1} = 0.1997(4) =~ 1/5 and gamma_{3D} = 1.253(4) =~ 5/4. Our results indicate that, within experimental error, they are given by simple fractions corresponding to the linear interpolations between the respective two-dimensional (Onsager) and four-dimensional (mean field) critical exponents. An analysis of our inverse susceptibility data chi^{-1}(T) vs. /T - Tc/ shows that these data lead to a value of gamma_{3D} compatible with gamma' = gamma and Tc = 4.51152(12), while gamma values obtained recently by high and low temperature series expansions and renormalization group methods are not.
△ Less
Submitted 14 November, 2002;
originally announced November 2002.
-
Intra-molecular coupling as a mechanism for a liquid-liquid phase transition
Authors:
Giancarlo Franzese,
Manuel I. Marques,
H. Eugene Stanley
Abstract:
We study a model for water with a tunable intra-molecular interaction $J_σ$, using mean field theory and off-lattice Monte Carlo simulations. For all $J_σ\geq 0$, the model displays a temperature of maximum density.For a finite intra-molecular interaction $J_σ> 0$,our calculations support the presence of a liquid-liquid phase transition with a possible liquid-liquid critical point for water, lik…
▽ More
We study a model for water with a tunable intra-molecular interaction $J_σ$, using mean field theory and off-lattice Monte Carlo simulations. For all $J_σ\geq 0$, the model displays a temperature of maximum density.For a finite intra-molecular interaction $J_σ> 0$,our calculations support the presence of a liquid-liquid phase transition with a possible liquid-liquid critical point for water, likely pre-empted by inevitable freezing. For J=0 the liquid-liquid critical point disappears at T=0.
△ Less
Submitted 29 October, 2002; v1 submitted 18 December, 2001;
originally announced December 2001.
-
Dynamic Scaling in Diluted Systems Phase Transitions: Deactivation trough Thermal Dilution
Authors:
M. I. Marques,
J. A. Gonzalo
Abstract:
Activated scaling is confirmed to hold in transverse field induced phase transitions of randomly diluted Ising systems. Quantum Monte Carlo calculations have been made not just at the percolation threshold but well bellow and above it including the Griffiths-McCoy phase. A novel deactivation phenomena in the Griffiths-McCoy phase is observed using a thermal (in contrast to random) dilution of th…
▽ More
Activated scaling is confirmed to hold in transverse field induced phase transitions of randomly diluted Ising systems. Quantum Monte Carlo calculations have been made not just at the percolation threshold but well bellow and above it including the Griffiths-McCoy phase. A novel deactivation phenomena in the Griffiths-McCoy phase is observed using a thermal (in contrast to random) dilution of the system.
△ Less
Submitted 1 September, 2000;
originally announced September 2000.
-
Evolution of the universality class in slightly diluted (1>p>0.8) Ising systems
Authors:
M. I. Marques,
J. A. Gonzalo
Abstract:
The crossover of a pure (undiluted) Ising system (spin per site probability p=1) to a diluted Ising system (spin per site probability p<0.8) is studied by means of Monte Carlo calculations with p ranging between 1 and 0.8 at intervals of 0.025. The evolution of the self averaging is analyzed by direct determination of the normalized square widths for magnetization and susceptibility as a functio…
▽ More
The crossover of a pure (undiluted) Ising system (spin per site probability p=1) to a diluted Ising system (spin per site probability p<0.8) is studied by means of Monte Carlo calculations with p ranging between 1 and 0.8 at intervals of 0.025. The evolution of the self averaging is analyzed by direct determination of the normalized square widths for magnetization and susceptibility as a function of p. We find a monotonous and smooth evolution from the pure to the randomly diluted universality class. The p-dependent transition is found to be independent of the size (L). This property is very convenient for extrapolation towards the randomly diluted universality class avoiding complications resulting from finite size effects.
△ Less
Submitted 14 March, 2000;
originally announced March 2000.
-
Revisiting "swings" in the crossover features of Ising thin films near Tc(D)
Authors:
M. I. Marques,
J. A. Gonzalo,
J. Romero,
L. F. Fonseca
Abstract:
"Swing" effects at the onset of crossover towards two dimensional behavior in thin Ising films are investigated close to Tc(D) by means of Monte Carlo calculations. We find that the effect is extremely large for the specific heat effective critical exponent, in comparison with the "swing" already noted by Capehart and Fisher for the susceptibility. These effects change considerably the system's…
▽ More
"Swing" effects at the onset of crossover towards two dimensional behavior in thin Ising films are investigated close to Tc(D) by means of Monte Carlo calculations. We find that the effect is extremely large for the specific heat effective critical exponent, in comparison with the "swing" already noted by Capehart and Fisher for the susceptibility. These effects change considerably the system's evolution with thickness (D) from two-dimensional to three-dimensional behavior, forcing the effective exponents to pass near characteristic Tri Critical Point (TCP) values.
△ Less
Submitted 17 February, 2000;
originally announced February 2000.
-
Universality Class of Thermally Diluted Ising Systems at Criticality
Authors:
M. I. Marques,
J. A. Gonzalo,
J. Iniguez
Abstract:
The universality class of thermally diluted Ising systems, in which the realization of the disposition of magnetic atoms and vacancies is taken from the local distribution of spins in the pure original Ising model at criticality, is investigated by finite size scaling techniques using the Monte Carlo method. We find that the critical temperature, the critical exponents and therefore the universa…
▽ More
The universality class of thermally diluted Ising systems, in which the realization of the disposition of magnetic atoms and vacancies is taken from the local distribution of spins in the pure original Ising model at criticality, is investigated by finite size scaling techniques using the Monte Carlo method. We find that the critical temperature, the critical exponents and therefore the universality class of these thermally diluted Ising systems depart markedly from the ones of short range correlated disordered systems. Our results agree fairly well with theoretical predictions previously made by Weinrib and Halperin for systems with long range correlated disorder.
△ Less
Submitted 26 January, 2000;
originally announced January 2000.
-
Possible evidence for a change in cosmic equation of state at decoupling
Authors:
J. A. Gonzalo,
M. I. Marques
Abstract:
Indirect but significative evidence (from fusion research data, and from the known 4He abundance) relevant to primordial nucleosynthesis is examined. It is shown that this evidence supports a change in cosmic equation of state at decoupling, compatible with Einstein's cosmological equations, from ROm T^4 aprox. cte (ROm=ROr) prior to decoupling (t<taf), to ROmT^3 aprox. cte after (t>taf), being…
▽ More
Indirect but significative evidence (from fusion research data, and from the known 4He abundance) relevant to primordial nucleosynthesis is examined. It is shown that this evidence supports a change in cosmic equation of state at decoupling, compatible with Einstein's cosmological equations, from ROm T^4 aprox. cte (ROm=ROr) prior to decoupling (t<taf), to ROmT^3 aprox. cte after (t>taf), being taf the atom formation time. Implications at the baryon threshold time (temperature) and beyond are discussed.
△ Less
Submitted 19 January, 2000; v1 submitted 13 May, 1999;
originally announced May 1999.
-
Self-averaging of random and thermally disordered diluted Ising systems
Authors:
M. I. Marques,
J. A. Gonzalo
Abstract:
Self-averaging of singular thermodynamic quantities at criticality for randomly and thermally diluted three dimensional Ising systems has been studied by the Monte Carlo approach. Substantially improved self-averaging is obtained for critically clustered (critically thermally diluted) vacancy distributions in comparison with the observed self-averaging for purely random diluted distributions. Cr…
▽ More
Self-averaging of singular thermodynamic quantities at criticality for randomly and thermally diluted three dimensional Ising systems has been studied by the Monte Carlo approach. Substantially improved self-averaging is obtained for critically clustered (critically thermally diluted) vacancy distributions in comparison with the observed self-averaging for purely random diluted distributions. Critically thermal dilution, leading to maximum relative self-averaging, corresponds to the case when the characteristic vacancy ordering temperature is made equal to the magnetic critical temperature for the pure 3D Ising systems. For the case of a high ordering temperature, the self-averaging obtained is comparable to that in a randomly diluted system.
△ Less
Submitted 22 April, 1999;
originally announced April 1999.