Effective pair potentials and structure of water

Classical Monte Carlo calculations have been performed in order to investigate the ability of the TIP4P, SPC, and SPCE water models to reproduce the structural features of liquid water. The simulations were carried out in the NPT ensemble at 4 thermodynamic conditions. The results are compared with recent neutron diffraction data. Essentially, the three models capture equally well the thermodynamic and structural features of water. Although they were parametrized to reproduce the water properties at ambient conditions, the agreement with the experimental pair correlation functions was even better at supercritical conditions. This is because the effective pair potentials have some difficulty to reproduce cooperative interactions, like hydrogen bonds. These interactions are less effective at supercritical conditions, where the liquid behaves roughly like a gas.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Finite-temperature correlations and density profiles of an inhomogeneous interacting one-dimensional Bose gas

We calculate the density profiles and density correlation functions of the one-dimensional Bose gas in a harmonic trap, using the exact finite-temperature solutions for the uniform case, and applying a local density approximation. The results are valid for a trapping potential that is slowly varying relative to a correlation length. They allow a direct experimental test of the transition from the weak-coupling Gross-Pitaevskii regime to the strong-coupling, fermionic Tonks-Girardeau regime. We also calculate the average two-particle correlation which characterizes the bulk properties of the sample, and find that it can be well approximated by the value of the local pair correlation in the trap center.

Domain growth in binary mixtures at low temperatures

We have studied domain growth during spinodal decomposition at low temperatures. We have performed a numerical integration of the deterministic time-dependent Ginzburg-Landau equation with a variable, concentration-dependent diffusion coefficient. The form of the pair-correlation function and the structure function are independent of temperature but the dynamics is slower at low temperature. A crossover between interfacial diffusion and bulk diffusion mechanisms is observed in the behavior of the characteristic domain size. This effect is explained theoretically in terms of an equation of motion for the interface.

Ordered structures in rotating ultracold Bose gases

Two-dimentional systems of trapped samples of few cold bosonic atoms submitted to strong rotation around the perpendicular axis may be realized in optical lattices and microtraps. We investigate theoretically the evolution of ground state structures of such systems as the rotational frequency Omega increases. Various kinds of ordered structures are observed. In some cases, hidden interference patterns exhibit themselves only in the pair correlation function; in some other cases explicit broken-symmetry structures appear that modulate the density. For N < 10 atoms, the standard scenario, valid for large sytems is absent, and is only gradually recovered as N increases. On the one hand, the Laughlin state in the strong rotational regime contains ordered structures much more similar to a Wigner molecule than to a fermionic quantum liquid. On the other hand, in the weak rotational regime, the possibility to obtain equilibrium states, whose density reveals an array of vortices, is restricted to the vicinity of some critical values of the rotational frequency Omega.

Statistical dynamics of dislocation systems: The influence of dislocation-dislocation correlations.

During plastic deformation of crystalline materials, the collective dynamics of interacting dislocations gives rise to various patterning phenomena. A crucial and still open question is whether the long range dislocation-dislocation interactions which do not have an intrinsic range can lead to spatial patterns which may exhibit well-defined characteristic scales. It is demonstrated for a general model of two-dimensional dislocation systems that spontaneously emerging dislocation pair correlations introduce a length scale which is proportional to the mean dislocation spacing. General properties of the pair correlation functions are derived, and explicit calculations are performed for a simple special case, viz pair correlations in single-glide dislocation dynamics. It is shown that in this case the dislocation system exhibits a patterning instability leading to the formation of walls normal to the glide plane. The results are discussed in terms of their general implications for dislocation patterning.

Long-range-interactions induced ordered structures in deposition processes

We present a new model of sequential adsorption in which the adsorbing particles experience dipolar interactions. We show that in the presence of these long-range interactions, highly ordered structures in the adsorbed layer may be induced at low temperatures. The new phenomenology is manifest through significant variations of the pair correlation function and the jamming limit, with respect to the case of noninteracting particles. Our study could be relevant in understanding the adsorption of magnetic colloidal particles in the presence of a magnetic field.

Influence of hydrodynamic interactions on the adsoption process of large particles

We have studied the adsorption process of non-Brownian particles on a line. Our work differs from previously proposed models in that we have incorporated hydrodynamic interactions between the incoming particles and the preadsorbed particles as well as the surface. We then numerically analyze the effect of these interactions on quantities related to the adsorption process. Comparing our model to the ballistic deposition model (BM) shows a significant discrepancy in the pair correlation function. These results can explain some differences between recent experiments and BM predictions. Finally, the limitations of the applicability of BM are addressed.

Quantum Hall Phases of two-component bosons

The recent production of synthetic magnetic fields acting on electroneutral particles, such as atoms or photons, has boosted interest in the quantum Hall physics of bosons. Adding pseudospin 1/2 to the bosons greatly enriches the scenario, as it allows them to form an interacting integer quantum Hall (IQH) phase with no fermionic counterpart. Here we show that, for a small two-component Bose gas on a disk, the complete strongly correlated regime, extending from the integer phase at filling factor ν = 2 to the Halperin phase at filling factor ν = 2 / 3, is well described by composite fermionization of the bosons. Moreover we study the edge excitations of the IQH state, which, in agreement with expectations from topological field theory, are found to consist of forward-moving charge excitations and backward-moving spin excitations. Finally, we demonstrate how pair-correlation functions allow one to experimentally distinguish the IQH state from competing states, such as non-Abelian spin singlet (NASS) states.

Quantum Hall Phases of two-component bosons

The recent production of synthetic magnetic fields acting on electroneutral particles, such as atoms or photons, has boosted interest in the quantum Hall physics of bosons. Adding pseudospin 1/2 to the bosons greatly enriches the scenario, as it allows them to form an interacting integer quantum Hall (IQH) phase with no fermionic counterpart. Here we show that, for a small two-component Bose gas on a disk, the complete strongly correlated regime, extending from the integer phase at filling factor ν = 2 to the Halperin phase at filling factor ν = 2 / 3, is well described by composite fermionization of the bosons. Moreover we study the edge excitations of the IQH state, which, in agreement with expectations from topological field theory, are found to consist of forward-moving charge excitations and backward-moving spin excitations. Finally, we demonstrate how pair-correlation functions allow one to experimentally distinguish the IQH state from competing states, such as non-Abelian spin singlet (NASS) states.

PREDICTING DIFFUSIVITIES IN DENSE FLUIDS

In this work the Weeks-Chandler-Andersen (WCA) perturbation theory coupled with the Enskog’s solution of the Boltzmann equation for dense hard-sphere fluids is employed for estimating diffusion coefficients in compressed pure liquids and fluids and dense fluid mixtures. The effect of density correction on the estimation of diffusivities is analyzed using the Carnahan-Starling pair correlation function and the correlation of Speedy and Harris which have been proposed as models of self-diffusion coefficient of hard-sphere fluids. The approach presented here is based on the smooth hard-sphere theory without any binary adjustable parameters and can be readily used for estimating diffusivities in multicomponent fluid mixtures. It is shown that the correlated and the predicted diffusivities are in good agreement with the experimental data and much better than estimates of Wilke-Chang equation.

Statistical dynamics of dislocation systems: The influence of dislocation-dislocation correlations.

During plastic deformation of crystalline materials, the collective dynamics of interacting dislocations gives rise to various patterning phenomena. A crucial and still open question is whether the long range dislocation-dislocation interactions which do not have an intrinsic range can lead to spatial patterns which may exhibit well-defined characteristic scales. It is demonstrated for a general model of two-dimensional dislocation systems that spontaneously emerging dislocation pair correlations introduce a length scale which is proportional to the mean dislocation spacing. General properties of the pair correlation functions are derived, and explicit calculations are performed for a simple special case, viz pair correlations in single-glide dislocation dynamics. It is shown that in this case the dislocation system exhibits a patterning instability leading to the formation of walls normal to the glide plane. The results are discussed in terms of their general implications for dislocation patterning.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Processo de difusão com agregação e reorganização espontânea em uma rede 2D

Difusive processes are extremely common in Nature. Many complex systems, such as microbial colonies, colloidal aggregates, difusion of fluids, and migration of populations, involve a large number of similar units that form fractal structures. A new model of difusive agregation was proposed recently by Filoche and Sapoval [68]. Based on their work, we develop a model called Difusion with Aggregation and Spontaneous Reorganization . This model consists of a set of particles with excluded volume interactions, which perform random walks on a square lattice. Initially, the lattice is occupied with a density p = N/L2 of particles occupying distinct, randomly chosen positions. One of the particles is selected at random as the active particle. This particle executes a random walk until it visits a site occupied by another particle, j. When this happens, the active particle is rejected back to its previous position (neighboring particle j), and a new active particle is selected at random from the set of N particles. Following an initial transient, the system attains a stationary regime. In this work we study the stationary regime, focusing on scaling properties of the particle distribution, as characterized by the pair correlation function ø(r). The latter is calculated by averaging over a long sequence of configurations generated in the stationary regime, using systems of size 50, 75, 100, 150, . . . , 700. The pair correlation function exhibits distinct behaviors in three diferent density ranges, which we term subcritical, critical, and supercritical. We show that in the subcritical regime, the particle distribution is characterized by a fractal dimension. We also analyze the decay of temporal correlations

Sol phase and sol-gel transition in SnO2 colloidal suspensions

The effect of concentration on the structure of SnO2 colloids in aqueous suspension, on their spatial correlation and on the gelation process was studied by small angle x-ray scattering (SAXS). The shape of the experimental SAXS curves varies with suspension concentration. For diluted suspensions ([SnO2] less than or equal to 0.13 mol L-1), SAXS results indicate the presence of colloidal fractal aggregates with an internal correlation length xi congruent to 20 Angstrom, without any noticeable spatial correlation between them. This suggests that the aggregates are spatially arranged without any significant interaction like in ideal gas structures. For higher concentrations ([SnO2] = 0.16, 0.32, and 0.64 mol L-1), the colloidal aggregates are larger (xi = 24 Angstrom) and exhibit a certain degree of spatial correlation between them. The pair correlation function corresponding to the sol with the highest concentration (0.92 mol L-1) reveals a rather strong short range order between aggregates, characteristic of a fluid-like structure, with an average nearest-neighbor distance between aggregates d(1) = 125 Angstrom and an average second-neighbor distance d(2) = 283 Angstrom. The pair distribution function remains essentially invariant during the sol-gel transition, suggesting that gelation involves the formation of a few points of connection between the aggregates resulting in a gel network constituted by essentially linear chains of clusters..