Molecular dynamics simulation of the dielectric constant of water: the effect of bond flexibility

The role of bond flexibility on the dielectric constant of water is investigated via molecular dynamics simulations using a flexible intermolecular potential SPC/Fw [Y. Wu, H. L. Tepper, and G. A. Voth, J. Chem. Phys. 128, 024503 (2006)]. Dielectric constants and densities are reported for the liquid phase at temperatures of 298.15 K and 473.15 K and the supercritical phase at 673.15 K for pressures between 0.1 MPa and 200 MPa. Comparison with both experimental data and other rigid bond intermolecular potentials indicates that introducing bond flexibility significantly improves the prediction of both dielectric constants and pressure-temperature-density behavior. In some cases, the predicted densities and dielectric constants almost exactly coincide with experimental data. The results are analyzed in terms of dipole moments, quadrupole moments, and equilibrium bond angles and lengths. It appears that bond flexibility allows the molecular dipole and quadrupole moment to change with the thermodynamic state point, and thereby mimic the change of the intermolecular interactions in response to the local environment.     

Simulations of high-dielectric Stockmayer fluids in hyperspherical geometry

The static dielectric properties of Stockmayer fluids are investigated in the hyperspherical geometry, S(3). Different methods of obtaining the static dielectric constant ε(r) are compared. Tested methods include the evaluation of the Kirkwood factor, fluctuations of the total dipole moment, and a two-center potential correlation formula to obtain the dielectric constant through effective interactions. With no coupling to the "surrounding," the different methods give consistent estimates of the dielectric constant. Adding a coupling to the surrounding gives large size dependencies and the two-center potential correlation formula breaks down. For low dipole moments, there is a good agreement in the dielectric constant with previous studies.     

The role of intermolecular hydrogen bond on dielectric properties in hydrogen-bonded material 5-bromo-9-hydroxyphenalenone: theoretical investigation

Dielectric properties of the hydrogen-bonded material, 5-bromo-9-hydroxyphenalenone (C(13)H(7)O(2)Br; BrHPLN), are investigated theoretically by means of electronic structure calculations and Monte Carlo simulations. The density functional calculations of BrHPLN crystals have revealed that the polarization per one molecule can be about 1.7 times larger than that of the isolated monomer. It is also found that there exists significant electron density (0.01 e bohr(-3)) in an intermolecular C-H···O region, which, together with the interatomic distances of 2.39 ? for H···O and 3.34 ? for C···O, suggests the existence of intermolecular weak hydrogen bonding that may enhance the molecular polarization. The induced polarization effects in various intermolecular configurations are evaluated with the Fragment Molecular Orbital method. In addition to the π-π stacking interactions, two types of "in plane" intermolecular weak hydrogen-bonding configurations are found to affect the molecular dipole moment most significantly. These effects are efficiently included in a Monte Carlo simulation method in terms of "dipole corrections" as functions of both the intermolecular arrangements and the intramolecular proton configurations. The application to the dielectric phase transition of a BrHPLN crystal shows that the dipole corrections almost double the transition temperature, toward better agreement with experiments, and qualitatively affect the temperature dependence of the dielectric constant. Discussions are given to support that the results will remain adequate and consistent even after explicit inclusion of the quantum tunneling effects.     

Dipole moment and dielectric relaxation of poly(dimethylsilazane) oligomers

The polarity of polysilazane chains is studied using as models hexamethyldisilazane (HMDSZ) and an oligomer of amine-terminated poly(1,1-dimethylsilazane) (PDMSZ). Solvent effects and the contribution of the atomic polarization to the dipole moment of both PDMSZ and HMDSZ are discussed. A qualitative interpretation of the dipole moment of the oligomer suggests that gauche states about the skeletal bonds may be disfavored with respect to the alternative trans states. The dielectric spectrum of PDMSZ presents a glass-rubber relaxation followed by two weak secondary absorptions. The experimental results show that the relaxation strength Δε of PDMSZ is lower than that corresponding to poly(1,1-dimethylsiloxane) (PDMS). The analysis of the glass-rubber relaxation by using the coupling scheme suggests that strong intermolecular interactions occur in the chains that presumably arise from hydrogen bonding interactions between amine groups. © 1993 John Wiley & Sons, Inc.     

The flexible, polarizable, thole-type interaction potential for water (TTM2-F) revisited

We present a revision of the flexible, polarizable, Thole-type interaction potential for water [J. Chem. Phys.2002, 116, 5115], which allows for condensed-phase simulations. The revised version (TTM2.1-F) of the potential correctly describes the individual water molecular dipole moment and alleviates problems arising at short intermolecular separations that can be sampled in the course of molecular dynamics and Monte Carlo simulations of condensed environments. Furthermore, its parallel implementation under periodic boundary conditions enables the efficient calculation of the macroscopic structural and thermodynamic properties of liquid water, as its performance scales superlinearly with up to a number of 64 processors for a simulation box of 512 molecules. We report the radial distribution functions, average energy, internal geometry, and dipole moment in the liquid as well as the density, dielectric constant, and self-diffusion coefficient at T = 300 K from (NVT) and (NPT) classical molecular dynamics simulations by using the revised version of the potential.     

Polarizabilities and dipole moments of benzaldehyde, benzoic acid and oxalic acid in polar and nonpolar solvents

The purpose of this report is to calculate the orientation polarizability of benzaldehyde, benzoic acid and oxalic acid in polar and nonpolar solvents. The calculations are based on the knowledge of permanent dipole moment of the solutions. Other important physical quantities such as refractive index, density, specific volume, dielectric constant, molar polarization and molar refractivity are also calculated. Dipole moments of the solutions are calculated by using measured dielectric constants of the solutions. The dielectric constant measurements were made at 100 kHz. Relationships between the polarizability and concentration, specific volume, dielectric constant and dipole moment of the solutions are suggested.     

Ordering of a monotropic nematic liquid crystal in dielectric medium at phase transition temperature: A computational analysis

Ordering in a monotropic nematic liquid crystal, 4’-n-butyl-4-cyanobiphenyl, was examined by the methods of quantum mechanics and intermolecular forces. The atomic charges and dipole moment at each atomic center were evaluated by the CNDO/2 method. The configurational energy was computed using the Rayleigh-Schrödinger perturbation method. The total interaction energies obtained by these computations were used to calculate the probability of each configuration in a dielectric medium (noninteracting and nonmessogenic solvent, benzene) at the phase transition temperature using the Maxwell-Boltzmann formula. In the dielectric medium, the energies/probabilities are redistributed, and there is a considerable rise in the probability of interactions, although the order of preference remains the same. An attempt was made to develop a new and interesting model of nematogen in a dielectric medium. A theoretical support is offered to the experimental observations.     

Intermolecular vibrational relaxation in liquids

The influence of intermolecular vibrational relaxation on dipole moment correlation functions, as obtained from IR band shapes, is discussed. It is explicitly shown that vibrational relaxation due to intermolecular interactions depends on the reorientational behaviour of the molecules in the liquid.Therefore, an a priori separation of the dipole moment correlation function into independent reorientational and vibrational factors is not generally possible. The implications for various procedures used to “correct” Raman and IR band shapes for vibrational relaxation are discussed.The expression derived for the intermolecular vibrational relaxation is used to calculate theoretically the effect of transition dipole-transition dipole coupling on dipole moment correlation functions.Experimental data obtained from isotopic dilution measurements support the interpretation of the isotopic dilution effect in terms of the transition dipole-transition dipole coupling.     

Space partitioning and the effects of molecular proximity on electrostatic moments of the crystalline formamide molecule

Application of two different spatial partitioning techniques to determine the dipole moment of the formamide molecule in the crystal indicates approximate cancellation of the effects of molecular truncation introduced by space partitioning and intermolecular interactions, leading to a solid state dipole moment which is in reasonable agreement with the gas phase value.     

Competition between dipolar and steric interactions in swallow-tailed compounds

A new class of swallow-tailed mesogen with the longitudinal dipole moment in the direction of the non-branched 'head' is presented. Mixtures of these compounds with 'conventional' swallow-tailed substances were investigated by polarization microscopy and dielectric measurements. In the mixtures, dipolar and steric interactions favour opposite arrangements of the molecules in the short range order.     

Polarizability and Kerr constant of proteins by boundary element methods

A precise implementation of the boundary element method has been applied to the computation of the polarizability and the Kerr constant of eight soluble proteins. The method is demonstrated to be accurate and precise by comparison with analytical values for spheroids. Two different integral equations have been solved: (1) an exact equation with explicit dielectric constant dependence, and (2) an exact equation for a metallic body. The dielectric dependence for the metallic body case is built in with a generalization of the ellipsoid formula. Both methods agree quantitatively with each other for low relative dielectric constants. A full tensor expression for the Kerr constant yields perfect agreement with experiment for some proteins and badly under reports for the rest. The protein structure is obtained from a crystallographic database and is assumed rigid throughout the TEB measurement. Electrolyte effects are neglected. The Kerr constant is dominated by the protein dipole moment and is quite sensitive to the orientation of the dipole moment relative to the principal axes of the polarizability tensor. Several possible reasons for the large discrepancy between some computed and measured values are discussed.     

Structure, thermodynamics, and liquid-vapor equilibrium of ethanol from molecular-dynamics simulations using nonadditive interactions

We present a molecular-dynamics simulation study of the bulk and liquid-vapor interfacial properties of ethanol using a polarizable force field based on the fluctuating charge (FQ) formalism, as well as the nonpolarizable CHARMM22 force field. Both models are competitive with respect to the prediction of ambient liquid properties such as liquid density, enthalpy of vaporization, dielectric constant, and self-diffusion constants. The polarizable model predicts an average condensed-phase dipole moment of 2.2 D associated with an induced liquid-phase dipole moment of 0.6 D; though qualitatively in agreement with earlier nonadditive models as well as recent Car-Parinello calculations, the current FQ model underestimates the condensed-phase dipole moment. In terms of liquid structure, both models are in agreement with recent neutron-diffraction results of liquid ethanol structure, although the polarizable model predicts the hydroxyl-hydrogen-hydroxyl-hydrogen structure factor in closer agreement with the experimental data. In terms of interfacial properties, both models predict ambient surface tension to within 4% of the experimental value of 22.8 dyncm, while overestimating the surface excess entropy by almost a factor of 2. Both models display the characteristic preferential orientation of interfacial molecules. The polarizable model allows for a monotonic variation of the average molecular dipole moment from the bulk value to that of the vapor phase. Consequently, there is a dramatic difference in the surface potential predicted by the polarizable and nonpolarizable models. The polarizable model estimates a surface potential of -209+/-3 mV, while the nonpolarizable model yields a value of -944+/-10 mV. Finally, based on the vapor-liquid equilibrium simulation data from several temperatures, we estimate the critical properties of both models. As observed with other FQ models for associating fluids (such as water and methanol), and counter to what one would anticipate by modeling more physically the electrostatic response to local environment, the current FQ model underestimates the critical temperature and overestimates the critical density of ethanol; moreover, the FQ model is, in this respect, equivalent to the underlying fixed-charge model. These results further suggest the need to revisit polarizable models in terms of quantitative vapor-liquid equilibrium prediction.     

The inclusion of electron correlation in intermolecular potentials: applications to the formamide dimer and liquid formamide

A test of the quality of the electrostatic properties and polarizabilities used in the nonempirical molecular orbital (NEMO) potential is carried out for formamide by calculating the molecular dipole moment and polarizability at the second-order M?ller–Plesset (MP2) level of theory. The molecular dipole moment is 11% lower at the MP2 level than at the Hartree–Fock (HF) level, whereas the isotropic part of the polarizability is increased by 36% by adding electron correlation and using a considerably larger basis set. The atomic charges, dipole moments and polarizabilities obtained at the HF level are rescaled to get the correct molecular properties at the MP2 level. The potential minimum for the cyclic dimer of formamide is −17.50 kcal/mol with the MP2-scaled properties and is significantly lower than other potentials give. Two intermolecular potentials are constructed and used in subsequent molecular dynamics simulations: one with the regular NEMO potential and the other with the rescaled MP2 properties. A damping of the electrostatic field at short intermolecular distances is included in the present NEMO model. The average energies for liquid formamide are lower for the MP2-scaled model and are in good agreement with experimental results. The lowering of the simulation energy for the MP2-scaled potential indicates the strong dispersive interactions in liquid formamide. Received: 20 March 2000 / Accepted: 18 April 2000 / Published online: 18 August 2000     

Dielectric response of one-dimensional polar chains

We propose a theory for the dielectric constant of materials made of parallel infinite one-dimensional chains of dipoles. Each dipole is allowed to rotate in three dimensions. Monte Carlo simulations show that the Kirkwood factor of the chain grows with increasing dipole moment much faster than in the case of three-dimensional polar fluids. With increasing dipole moment or cooling the one-dimensional chain undergoes a continuous order-disorder transition to the ferroelectric phase, in which the dielectric constant is limited by the size of ferroelectric domains along the chain.     

The Influence of Large Pendent Groups on Chain Anisotropy and Electrical Energy Loss of Polyimides at High Frequency through All-Atomic Molecular Simulation

Polyimide is a potential material for high-performance printed circuit boards because of its chemical stability and excellent thermal and mechanical properties. Flexible printed circuit boards must have a low static dielectric constant and dielectric loss to reduce signal loss in high-speed communication devices. Engineering the molecular structure of polyimides with large pendant groups is a strategy to reduce their dielectric constant. However, there is no systematic study on how the large pendant groups influence electrical energy loss. We integrated all-atomic molecular dynamics and semi-empirical quantum mechanical calculations to examine the influence of pendant groups on polymer chain anisotropy and electrical energy loss at high frequencies. We analyzed the radius of gyration, relative shape anisotropy, dipole moment, and degree of polarization of the selected polyimides (TPAHF, TmBPHF, TpBPHF, MPDA, TriPMPDA, m-PDA, and m-TFPDA). The simulation results show that anisotropy perpendicular to chain direction and local chain rigidity correlate to electrical energy loss rather than dipole moment magnitudes. Polyimides with anisotropic pendant groups and significant local chain rigidity reduce electrical energy loss. The degree of polarization correlated well with the dielectric loss with a moderate computational cost, and difficulties in directly calculating the dielectric loss were circumvented.     

Ordering in higher homologous series of <Emphasis Type="Italic">p-n</Emphasis>-alkylbenzoic acids having eight alkyl chain carbon atoms: a computational analysis

A computational analysis of ordering in the higher homologous series of p-n-alkylbenzoic acids (nBAC) having 8(8BAC) carbon atoms in the alkyl chain was performed based on quantum mechanics and intermolecular forces. The atomic charge and dipole moment at each atomic centre were evaluated using the all-valence electron CNDO/2 method. The modified Rayleigh-Schrödinger perturbation theory and the multi-center-multipole expansion method were employed to evaluate long-range intermolecular interactions, while a 6-exp potential function was assumed for short-range interactions. The total interaction energy values obtained in these computations were used as input for calculating the probability of each configuration in a noninteracting and nonmesogenic solvent (i.e., benzene) at room temperature (300 K) using the Maxwell-Boltzmann formula. The molecular parameters of 8BAC, including the total energy, binding energy, and total dipole moment, were compared with those of 7BAC and 9BAC. The present article offers a theoretical support for the experimental findings.     

Study of structural and dynamic properties of liquid phenyltrimethoxysilane

Herein, we present a combined experimental and computational study of liquid phenyltrimethoxysilane. A femtosecond time-resolved optical Kerr effect experiment has been performed to study the rotational diffusion of the molecule. A new all-atoms molecular model of the compound, based on the OPLS force field, has been developed to reproduce the rotational diffusion time constant and other physical and dynamic properties available in the literature. The density obtained from the simulations is 1074 ± 4 kg m(-3), which is within 1% of the experimental value of 1062 kg m(-3). The viscosity from the simulations is 1.6 ± 0.1 mPa s while the experimental value is 2.1 mPa s. The average bulk dipole moment of 1.8 ± 0.5 Debye obtained from the simulation matches the experimental value of 1.77 Debye. The average relative dielectric constant from the simulations is 3.86 ± 0.04, which is within 13% of the experimental value (4.4). The rotational diffusion time of the dipole moment obtained from the simulations is 20.39 ± 0.06 ps, which is in excellent agreement with the experimental value of 20 ± 1 ps obtained from our measurements. The new model has also been used to calculate structural and dynamic properties of the molecule not yet determined experimentally.     

Thermodynamical study of the thermoelectric effect for magnesium silicide

The thermoelectric effect of magnesium silicide is studied by using a thermodynamical method in the presence of an electric field. The thermoelectric potential is evaluated from the partial derivative of free energy with respect to charge in which the free energy is calculated at the B3LYP/6-31G(d,p) level of density functional theory. This free energy is also utilized to determine the average dipole moment from which the polarizability, alpha; molar polarization, Psi; and dielectric constant can be computed. The present calculation for the dielectric constant (approximately 24-20) is in very good agreement with the experimental value (20). This accurate dielectric constant can be used to derive the relation of the thermoelectric potential with respect to temperature, from which the thermoelectric power or the Seebeck coefficients are calculated. The present result shows good agreement with experiment measurement for the Seebeck coefficients. In comparison, that calculation from the energy band structure theory is far off from the experimental values.     

Quantum-Mechanical and Solvatochromic Characterization of Quercetin

ABSTRACT

Quercetin was subjected to a quantum-mechanical analysis and energetic and electro-optic parameters were evaluated. A solvatochromic study was used to estimate the dipole moment of quercetin in the excited state and to determine the angle between the dipole moments of this molecule in the electronic states responsible for the visible absorption band. The contribution of each type of intermolecular interactions to the total spectral shift experimentally recorded in quercetin solution was also established in this study.