Are dipolar liquids ferroelectric? Simulation studies

Dielectric continuum model examination of real-space electrostatic treatments

Electrostatic interaction is long ranged; thus, the accurate calculation is not an easy task in molecular dynamics or Monte Carlo simulations. Though the rigorous Ewald method based on the reciprocal space has been established, real-space treatments have recently become an attractive alternative because of the efficient calculation. However, the construction is not yet completed and is now a challenging subject. In an earlier theoretical study, Neumann and Steinhauser employed the Onsager dielectric continuum model to explain how simple real-space cutoff produces artificial dipolar orientation. In the present study, we employ this continuum model to explore the fundamental properties of the recently developed real-space treatments of three shifting schemes. The result of the distance-dependent Kirkwood function G_K(R) showed that the simple bare cutoff produces a well-known hole-shaped artifact, whereas the shift treatments do not. Two-dimensional mapping of electric field well explained how these shift treatments remove the hole-shaped artifact. Still, the shift treatments are not sufficient because they do not produce a flat G_K(R) profile unlike ideal no-cutoff treatment. To test the continuum model results, we also performed Monte Carlo simulations of dipolar particles. The results found that the continuum model could predict the qualitative tendency as to whether each electrostatic treatment produces the hole-shaped artifact of G_K(R) or not. We expect that the present study using the continuum model offers a stringent criterion to judge whether the primitive electrostatic behavior is correctly described or not, which will be useful for future construction of electrostatic treatments.

Polarization of acetonitrile under thermal fields via non-equilibrium molecular dynamics simulations

We show that thermal gradients polarize liquid and supercritical acetonitrile. The polarization results in a stationary electrostatic potential that builds up between hot and cold regions. The strength of the field increases with the static dielectric constant or with decreasing temperature. At near standard conditions, the thermal polarization coefficient is ∼-0.6 mV/K, making it possible to induce significant electrostatic fields, ∼10³ V/m, with thermal gradients ∼1 K/μm. At supercritical conditions, ∼600 K and 0.249 g/cm³ (the critical isochore), the electrostatic field is of the same order, despite the low dielectric constant of the fluid. In this case, the electrostatic field is determined by the enhanced rotational diffusion of the molecules and stronger cross-coupling between heat and polarization fluxes. We show that the coupling between the heat and polarization fluxes influences the thermal conductivity of acetonitrile, which becomes a worse heat conductor. For the thermodynamic states investigated in this work, the thermal polarization effect leads to a ∼2%-5% reduction in thermal conductivity.

What Does Second-Harmonic Scattering Measure in Diluted Electrolytes?

We derive a theoretical expression of the second harmonic scattering signal in diluted electrolytes compared with bulk water. We show that the enhancement of the signal with respect to pure water observed recently for electrolytes at very low dilution in the micromolar range is a mere manifestation of the Debye screening that makes the infinite-range dipole-dipole solvent correlations in 1/ r³ disappear as soon as the ionic concentration becomes finite. In q space, this translates into a correlation function having a well known singular behavior around q = 0, which drives the observed ionic effects. We find that the signal is independent of the ion-induced long-range behavior of the function ⟨cos ϕ( r)⟩ that has been recently discussed. We find also that the enhancement depends on the experimental geometry and occurs only for in-plane polarization detection, as observed experimentally. On the contrary, the measured isotope effect between light and heavy water cannot be fully explained.

Second-Harmonic Scattering as a Probe of Structural Correlations in Liquids

Second-harmonic scattering experiments of water and other bulk molecular liquids have long been assumed to be insensitive to interactions between the molecules. The measured intensity is generally thought to arise from incoherent scattering due to individual molecules. We introduce a method to compute the second-harmonic scattering pattern of molecular liquids directly from atomistic computer simulations, which takes into account the coherent terms. We apply this approach to large-scale molecular dynamics simulations of liquid water, where we show that nanosecond second-harmonic scattering experiments contain a coherent contribution arising from radial and angular correlations on a length scale of ≲1 nm, much shorter than had been recently hypothesized ( Shelton , D. P. J. Chem. Phys. 2014 , 141 ). By combining structural correlations from simulations with experimental data ( Shelton , D. P. J. Chem. Phys. 2014 , 141 ), we can also extract an effective molecular hyperpolarizability in the liquid phase. This work demonstrates that second-harmonic scattering experiments and atomistic simulations can be used in synergy to investigate the structure of complex liquids, solutions, and biomembranes, including the intrinsic intermolecular correlations.

Size effect and stability of polarized fluid phases

The existence of a ferroelectric fluid phase for systems of 1000-2000 dipolar hard or soft spheres is well established by numerical simulations. Theoretical approaches proposed to determine the stability of such a phase are either in qualitative agreement with the simulation results or disagree with them. Experimental results for systems of molecules or particles with large electric or magnetic dipole moments are also inconclusive. As a contribution to the question of existence and stability of a fluid ferroelectric phase this simulation work considers system sizes of the order of 10 000 particles, thus an order of magnitude larger than those used in previous studies. It shows that although ferroelectricity is not affected by an increase of system size, different spatial arrangements of the dipolar hard spheres in such a phase are possible whose free energies seem to differ only marginally.

Thermodynamics of ferrofluids in applied magnetic fields

The thermodynamic properties of ferrofluids in applied magnetic fields are examined using theory and computer simulation. The dipolar hard sphere model is used. The second and third virial coefficients (B(2) and B(3)) are evaluated as functions of the dipolar coupling constant λ, and the Langevin parameter α. The formula for B(3) for a system in an applied field is different from that in the zero-field case, and a derivation is presented. The formulas are compared to results from Mayer-sampling calculations, and the trends with increasing λ and α are examined. Very good agreement between theory and computation is demonstrated for the realistic values λ≤2. The analytical formulas for the virial coefficients are incorporated in to various forms of virial expansion, designed to minimize the effects of truncation. The theoretical results for the equation of state are compared against results from Monte Carlo simulations. In all cases, the so-called logarithmic free energy theory is seen to be superior. In this theory, the virial expansion of the Helmholtz free energy is re-summed in to a logarithmic function. Its success is due to the approximate representation of high-order terms in the virial expansion, while retaining the exact low-concentration behavior. The theory also yields the magnetization, and a comparison with simulation results and a competing modified mean-field theory shows excellent agreement. Finally, the putative field-dependent critical parameters for the condensation transition are obtained and compared against existing simulation results for the Stockmayer fluid. Dipolar hard spheres do not undergo the transition, but the presence of isotropic attractions, as in the Stockmayer fluid, gives rise to condensation even in zero field. A comparison of the relative changes in critical parameters with increasing field strength shows excellent agreement between theory and simulation, showing that the theoretical treatment of the dipolar interactions is robust.

Chain formation and aging process in biocompatible polydisperse ferrofluids: experimental investigation and Monte Carlo simulations

We review the use of Monte Carlo simulations in the description of magnetic nanoparticles dispersed in a liquid carrier. Our main focus is the use of theory and simulation as tools for the description of the properties of ferrofluids. In particular, we report on the influence of polydispersity and short-range interaction on the self-organization of nanoparticles. Such contributions are shown to be extremely important for systems characterized by particles with diameters smaller than 10nm. A new 3D polydisperse Monte Carlo implementation for biocompatible magnetic colloids is proposed. As an example, theoretical and simulation results for an ionic-surfacted ferrofluid dispersed in a NaCl solution are directly compared to experimental data (transmission electron microscopy - TEM, magneto-transmissivity, and electron magnetic resonance - EMR). Our combined theoretical and experimental results suggest that during the aging process two possible mechanisms are likely to be observed: the nanoparticle's grafting decreases due to aggregate formation and the Hamaker constant increases due to oxidation. In addition, we also briefly discuss theoretical agglomerate formation models and compare them to experimental data.

Hyper-Rayleigh scattering from correlated molecules

The polarization dependence of hyper-Rayleigh scattering has been calculated for spherical domains of orientation correlated molecules. Distributions with radial or azimuthal mean polar orientation of the molecules are found that give results consistent with experimental observations, and expressions for the polarization ratios in terms of the product of correlation strength and correlated domain size are derived for these distributions. Assuming a plausible correlation strength, it is estimated that the correlated domain size in typical polar liquids is of order 100 molecular diameters.

Simulation study of the polarizable Stockmayer fluid in an external field

Gas-liquid phase coexistence curves of the polarizable Stockmayer fluid in external electric fields are computed using molecular dynamics computer simulation. We study in particular the critical-point shift dependence on polarizability and external electric field distinguishing the cases of fixed charge density and fixed potential. The results are compared to a previously developed mean-field theory for the polarizable Stockmayer fluid in an external field. We also investigate the behavior of the isochoric heat capacity near gas-liquid criticality via finite-size scaling depending on polarizability and external field.

Dielectric constants of simple liquids: stockmayer and ellipsoidal fluids

Coarse-grained models of molecular interactions are of interest because they convey the essence of molecular interactions in simple and easy to understand form. However, coarse-grained models fail to adequately predict some material properties, such as the failure of the Stockmayer model to reproduce the dielectric behavior of highly polar liquids. We examine the behavior of the Stockmayer fluid over a range of dipole densities that covers known organic solvents, as well as that of an ellipsoidal Stockmayer-like fluid, using NVT rigid-body Monte Carlo simulations. Both fluids are examined under different electrostatic boundary conditions and ensemble sizes. While the Stockmayer model predicts that liquids of similar dipole density to acetonitrile would be ferroelectric and have a dielectric constant far higher than shown by experiment, the ellipsoidal model provides a better accounting of dielectric behavior. This result bodes well for the use of coarse-grained solvent models for large-scale simulations.

Dipolar particles in an external field: Molecular dynamics simulation and mean field theory

Using molecular dynamics computer simulation we compute gas-liquid phase coexistence curves for the Stockmayer fluid in an external electric field. We observe a field-induced shift of the critical temperature DeltaTc. The sign of DeltaTc depends on whether the potential or the surface charge density is held constant assuming that the dielectric material fills the space between capacitor plates. Our own as well as previous literature data for DeltaTc are compared to and interpreted in terms of a simple mean field theory. Despite considerable errors in the simulation results, we find consistency between the simulation results obtained by different groups including our own and the mean field description. The latter ties the sign of DeltaTc to the outside constraints via the electric field dependence of the orientation part of the mean field free energy.

Phase behavior of the Stockmayer fluid via molecular dynamics simulation

The gas-isotropic liquid-nematic liquid phase behavior of the Stockmayer fluid is studied using molecular dynamics simulation together with a mean field lattice model. We obtain coexistence curves of the Stockmayer fluid over a wide range of dipole strengths, temperatures, and densities, including the transition from the isotropic liquid to the ferroelectric liquid. In our simulations we do not observe the disappearance of the isotropic gas-isotropic liquid coexistence at high dipole strength contrary to earlier findings based on Monte Carlo techniques. Even though the formation of reversible dipole chains strongly affects the location of the critical point, it does not lead to its disappearance. These results are supported by a mean field lattice model which yields good qualitative, and in parts quantitative, agreement with our simulations. In addition, we also investigate the gas-isotropic liquid phase behavior for different polarizabilities.