Quantum rigid dipole in a permanent electric field. II. Model of librational motion in liquid water and ice Ih: Preliminary results

A theoretical one-body model of librational motion in liquid water and ice Ih is proposed within the quantum mechanical treatment of a rigid and pointlike dipole in a permanent electric field and the concept of the internal electric field justified for V-structure of liquid water.     

Exchange-hole dipole moment and the dispersion interaction: high-order dispersion coefficients

In recent publications [A. D. Becke and E. R. Johnson, J. Chem. Phys. 122, 154104 (2005); E. R. Johnson and A. D. Becke 123, 024101 (2005)] we have demonstrated that the position-dependent dipole moment of the exchange hole can be used to generate dispersion interactions between closed-shell systems. Remarkably accurate C6 coefficients and intermolecular potential-energy surfaces can be obtained from Hartree-Fock occupied orbitals and polarizability data alone. In the present work, our model is extended to predict C8 and C10 coefficients as well. These higher-order coefficients are obtained as easily as C6 and with comparable accuracy.     

The polarization and dipole moment of lattice water in hydrogen bonded crystals

Dipole moments and charge distributions are calculated for the water molecules in four hydrogen bonded crystals, α and β oxalic dihydrates and the hexagonal and cubic ices.     

Study of the molecular configuration and the dipole moment in fluorinated liquid crystals

We have investigated the relationship between the molecular configuration and dipole moment of some fluorinated liquid crystals (LCs). The geometries of the molecules were preliminarily optimized at empirical AM1 and then were further optimized at B3LYP/6‐31G(d) level. The dipole moment has been calculated. It is strongly influenced by the position and number of fluorine substituents in the benzene ring of the molecule. The polarizability, mean polarizabilities, and anisotropic polarizability of the phenylbicyclohexane (PBC) fluorine substituents are also given and discussed. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Elastic properties of bent-core nematic liquid crystals: the role of the bend angle

ABSTRACT

The nematic phase formed from bent-core liquid crystals has been the focus of intensive research for more than a decade. With the potential of biaxial nematic phase and other interesting features, such as high Kerr constants, large flexoelectric coefficients and anomalous elastic constants, these bent-core materials have been in the limelight of research. This paper presents a mini-review of the interesting elastic behaviour reported in various bent-core compounds. The review further focusses on two different types of bent-core compounds as exemplars: one derived from an oxadiazole and other a thiadiazole, highlighting the importance of bend angle in bent-core compounds. The universality of the unique behaviour of bend elastic constants via molecular field theory and atomistic modelling has also been summarised.     

The permanent electric dipole moment and hyperfine interaction in ruthenium monoflouride (RuF)

Ruthenium monofluoride, RuF, has been detected using low-resolution laser-induced fluorescence (LIF) in the visible and near infrared spectral regions. A visible band, designated as [18.2]5.5-X 4Phi(9/2), has been recorded field-free and in the presence of a static electric field using high-resolution LIF spectroscopy. The r0 internuclear distances for the [18.2]5.5 and X 4Phi(9/2) states were determined to be 1.911 and 1.916 A, respectively. The vibrational interval DeltaG(1/2) of 534(15) cm-1 for the X 4Phi(9/2) state was determined from the analysis of the dispersed LIF. The Stark shifts of the visible band were analyzed to produce permanent electric dipole moments of 1.97(8) and 5.34(7) D for the [18.2]5.5 and X 4Phi(9/2), states, respectively. The fluorine magnetic hyperfine structure associated with spectral features was analyzed. The hyperfine structure and dipole moments are interpreted using a molecular-orbital correlation model and compared with FeF and other ruthenium-containing molecules.     

Quantum effects in liquid water and ice: model dependence

This paper explores the influence of choice of potential model on the quantum effects observed in liquid water and ice. This study utilizes standard rigid models and a more formal context for the rigid-body centroid molecular dynamics methodology used to perform the quantum simulations is provided. Quantum and classical molecular dynamics simulations are carried out for liquid water and ice Ih at 298 and 220 K, respectively, with the simple point charge/extended and TIP4P-Ew water models. The results obtained for equilibrium and dynamical properties are compared with those recently reported on TIP4P [L. Hernandez de la Pena and P. G. Kusalik, J. Chem. Phys. 121, 5992 (2004); L. Hernandez de la Pena et al., J. Chem. Phys 123, 144506 (2005)]. For the liquid, an energy shift of about 8% and an average molecular uncertainty of about 11 degrees were found independently of the water model. The self-diffusion coefficient consistently increases by more than 50% when going from the classical to the quantum system and quantum dynamics are found to reproduce the experimental isotopic shifts with the models examined. The ice results compare remarkably well with those previously reported for the TIP4P water model; they confirm that quantum effects are considerable and that the quantum mechanical uncertainty and the energy shifts due to quantization are smaller in ice than in liquid water. The relevance of these findings in the context of the construction of water models is briefly discussed.     

A global, high accuracy ab initio dipole moment surface for the electronic ground state of the water molecule

A highly accurate, global dipole moment surface (DMS) is calculated for the water molecule using ab initio quantum chemistry methods. The new surface is named LTP2011 and is based on all-electron, internally contracted multireference configuration interaction, including size-extensivity corrections in the aug-cc-pCV6Z basis set. Dipoles are computed as energy derivatives and small corrections due to relativistic effects included. The LTP2011 DMS uses an appropriate functional form that guarantees qualitatively correct behaviour even for most high energies configuration (up to about 60,000 cm(-1)), including, in particular, along the water dissociation channel. Comparisons with high precision experimental data show agreement within 1% for medium-strength lines. The new DMS and all the ab initio data are made available in the supplementary material.     

The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface

In this work we present an implementation for the calculation of the melting point of ice I(h) from direct coexistence of the solid-liquid interface. We use molecular dynamics simulations of boxes containing liquid water and ice in contact. The implementation is based on the analysis of the evolution of the total energy along NpT simulations at different temperatures. We report the calculation of the melting point of ice I(h) at 1 bar for seven water models: SPC/E, TIP4P, TIP4P-Ew, TIP4P/ice, TIP4P/2005, TIP5P, and TIP5P-E. The results for the melting temperature from the direct coexistence simulations of this work are in agreement (within the statistical uncertainty) with those obtained previously by us from free energy calculations. By taking into account the results of this work and those of our free energy calculations, recommended values of the melting point of ice I(h) at 1 bar for the above mentioned water models are provided.     

Fifth-order two-dimensional Raman spectroscopy of liquid water, crystalline ice Ih and amorphous ices: sensitivity to anharmonic dynamics and local hydrogen bond network structure

The theoretical study of off-resonant fifth-order two-dimensional (2D)-Raman spectroscopy is made to analyze the intermolecular dynamics of liquid and solid water. The 2D-Raman spectroscopy is susceptible to the nonlinear anharmonic dynamics and local hydrogen bond structure in water. It is found that the distinct 2D-Raman response appears as the negative signal near the t(2) axis. The origin of this negative signal for t(2)<15 fs is from the nonlinear polarizability in the librational motions, whereas that for 30 fs相似文献   

The melting point of hexagonal ice (Ih) is strongly dependent on the quadrupole of the water models

It is shown that quadrupolar interactions play a determinant role in the melting temperatures of common water models and that there is a simple relationship between the strength of the quadrupolar forces and the position of the negative charge; our conclusion is that acceptable predictions for the melting temperature can only be obtained when the negatively charged site is shifted from the oxygen atom towards the hydrogens.     

On the origin dependence of the angle made by the electric and magnetic vibrational transition dipole moment vectors

The concept of robustness of rotational strengths of vibrational modes in a VCD spectrum has been introduced as an aid in assignment of the absolute configuration with the help of the VCD spectrum. The criteria for robustness have been based on the distribution around 90° of the angles ξ(i) between electric and magnetic transition dipoles of all the modes i of a molecule. The angles ξ(i) (not, of course, the rotational strengths) are, however, dependent on the choice of origin. The derived criteria are for the center of mass chosen as the origin of the coordinate system. We stress in this note that application of the derived criteria assumes that excessive translation of the coordinate origin is not applied. Although the ξ(i) angles are not very sensitive to the position of the origin, very small displacements (a few ?) are not a problem, excessive translation of the origin does have considerable effect on the ξ(i) angles. In this note we quantify this effect and demonstrate how the distribution of ξ(i) angles is affected. Although it is possible to recalibrate the robustness criteria for the angles for a specific (large) displacement, we recommend that such displacement simply be avoided. It is to be noted that some modeling software does yield output with excessively displaced coordinate origin; this should be checked and corrected.     

Hydrogen bond and dipole moment in sub- and supercritical water close to the saturation curve

The Car-Parrinello molecular dynamics method was used to determine the distributions of water dipole moments under normal conditions, at the critical point, and in six thermodynamic states in sub- and supercritical phase diagram regions. Dipole moment changes along the saturation curve, the 650 K isotherm, and the 30 MPa isobar were analyzed.     

Predicting the infrared transition intensities in the Ar-HF complex: the key role of the dipole moment surface accuracy

The method proposed earlier for the generation of the full-dimensional energy surface for van der Waals complexes [P. Jankowski, J. Chem. Phys. 121, 1655 (2004)] is used to obtain a fulldimensional dipole moment surface for the atom-diatom complex in calculations based on the coupled-cluster with single, double, and noniterative triple excitation approach and the aug-cc-pVQZ basis sets. This surface has been employed to calculate transition intensities of the infrared spectra of Ar-HF. Special attention has been paid to study the problem of relative intensities of the different bands which have not been properly predicted within the long-range models of the dipole moment [A. E. Thornley and J. M. Hutson, J. Chem. Phys. 101, 5578 (1994)]. The intensities calculated with the present dipole moment surface agree very well with the experimental data, which indicate that the short-range interactions significantly affect the dipole moment surface and the calculated intensities. To investigate the role of the accuracy of the dipole moment surface on infrared transition intensities in atom-diatom complexes, four models of increasing complexity are studied. Their performance is shown to strongly depend on the region of the interaction energy surface probed by the initial and final states of the individual transitions.     

On the trapping of Bjerrum defects in ice Ih: the case of the molecular vacancy

We present a density-function theory (DFT) study of Bjerrum-defect trapping centers involving the molecular vacancy in ice Ih. As a first step, we compute the intrinsic migration barrier to D-defect motion using the nudged elastic band (NEB) method and find them to be of the same order of magnitude as the energy barriers involving intrinsic L-defect motion. This finding suggests that intrinsic mobility factors cannot explain the experimentally observed inactivity of D defects, supporting the idea that D defects are trapped at other lattice-defect sites. Next we study the defect complexes formed by the combination of isolated D and L defects with a molecular vacancy. The corresponding geometries show that the formation of these aggregates significantly reduces elastic distortions that are present in isolated Bjerrum defects. An analysis of the energetics involved in the formation of both defect complexes reveals a significant binding energy, indicating that the molecular vacancy represents a strong trapping center for Bjerrum defects. On the other hand, the fact that there is no difference between the absolute values of the binding energies for both D and L defects suggests that the vacancy affects both species of Bjerrum defects in a similar fashion, possibly ruling out the vacancy trapping centers as an explanation for the experimentally observed inactivity of D defects.     

Pretransitional behavior in the isotropic phase of a nematic liquid crystal with the transverse permanent dipole moment

The results presented give the evidence for the quasicritical, pretransitional behavior of dielectric properties in the isotropic phase of a rodlike nematic liquid crystal with the transverse permanent dipole moment. Studies were conducted in 2-cyano-4-pentylbiphenyl 4-(trans-4-pentylcyclohexyl) benzoate, focusing on the static-and ionic-dominated low-frequency (LF) regions. For the static dielectric permittivity [epsilon(')(T)] the application of the derivative analysis revealed the pretransitional anomaly associated with the specific heat exponent alpha approximately 0.5. For the LF domain the contribution to epsilon(')(T) from residual ionic impurities follows a linear temperature dependence on approaching the isotropic-nematic (I-N) transition. This dependence and pretransitional anomalies of electric conductivity and dielectric modulus can be associated with the influence of prenematic fluctuations. "Linear" dielectric studies were supported by the static nonlinear dielectric effect measurements, which delivered reliable estimations of the temperature of the hypothetical continuous phase transition T(*) and the discontinuity of the I-N transition DeltaT approximately 1.7 K.