Molecular simulation of static hyper-Rayleigh scattering: A calculation of the depolarization ratio and the local fields for liquid nitrobenzene

Modelling water diffusion in plasticizers: development and optimization of a force field for 2,4-dinitroethylbenzene and 2,4,6-trinitroethylbenzene

A classical all-atom force field has been developed for 2,4,6-trinitroethylbenzene and 2,4-dinitroethylbenzene and applied in molecular dynamics simulations of the two pure and two mixed plasticizer systems. Bonding parameters and partial charges were derived through electronic and geometry optimization of the single molecules. The other required parameters were derived from values already available in the literature for generic nitro aromatic compounds, which were adjusted to reproduce to a high level of accuracy the densities of 2,4-dinitroethylbenzene, 2,4,6-trinitroethylbenzene and the energetic plasticizers K10 and R8002. This force field has been applied to both K10 and R8002, which when used as plasticizers form an energetic binder with nitrocellulose. Nitrocellulose decomposes in storage, under varying conditions, but in particular where it may become increasingly dry. Following the derivation of the force field, we have therefore applied it to calculate water diffusion coefficients for each of the different materials at 298 K and 338 K, thereby providing a starting point for understanding water behaviour in a nitrocellulose binder.

A classical all-atom force field has been developed for the plasticizer molecules 2,4,6-trinitroethylbenzene and 2,4-dinitroethylbenzene which can be used to investigate properties and energetic output of nitrocellulose-based propellants and binders.

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.

The role of electronic properties to the mutagenic activity of 1,6- and 3,6-dinitrobenzo[a]pyrene isomers

Equilibrium geometries, infrared spectra, vertical first ionization potential (IP), electronic affinity (EA), dipole moment (mu) and electronic dipole polarizability (alpha) of 1,6- and 3,6-dinitrobenzo[a]pyrene isomers (1,6-DNBaP and 3,6-DNBaP) were evaluated by means of Density Functional Theory (DFT) and recent semiempirical PM6 method. Structural, energetic and vibrational properties of DNBaP isomers are substantially similar to each other. Calculated IP, EA and alpha values of these isomers are practically identical, while mu of 3,6-DNBaP (8.2 D at DFT level) is predicted to be ca. 4 times the value of 1,6-DNBaP isomer (1.9 D at DFT level), owing to favorable mutual orientation of the individual nitro group vectors. Higher direct-mutagenic activities of 3,6-DNBaP with respect to 1,6-DNBaP isomer by 1-2 orders of magnitude might be determined by its peculiar electronic charge distribution, which through stronger electrostatic and inductive interactions, can promote much more effectively binding to active-site of enzymes involved in mutagenic pathways. On the other hand, orientation of the nitro substituents relatively to the plane of the aromatic moiety, molecular sizes, as well as nitroreduction and oxidation reactions seem not to have a key role in the determination of the different mutagenic behaviour of these isomers.

Anomalous temperature behavior of nonlinear dielectric effect in supercooled nitrobenzene

Studies of the nonlinear dielectric effect (NDE), describing changes of the dielectric permittivity induced by a strong electric field, in normal and supercooled liquid nitrobenzene are presented. An unusual increase of the stationary NDE, portrayed by the critical-like relation approximately (T-T+)(-1), was obtained. Nitrobenzene samples solidified at TS approximately T+ +2 K , approximately 10 K below the reported melting temperature Tm approximately 278 K . The anomalous increase of the NDE coincided with a slow relaxation process, detected in time-resolved measurements. The results presented offer a reinterpretation of the classical description by Piekara and Piekara [A. Piekara and B. Piekara, C. R. Hebd. Seances Acad. Sci. 203, 852 (1936)] of a positive sign contribution to the NDE in liquids. It is suggested that the obtained anomaly may be associated with the appearance of local quasinematic structures. This is supported by a speculative link to a general model for liquid-liquid transitions [H. Tanaka, Phys. Rev. E 62, 6968 (2000)] and a phenomenological model originally developed for the self-focusing of laser beams [J. Hanus, Phys. Rev. 178, 420 (1969)]. The case of the isotropic-nematic transition in liquid crystalline materials is also recalled. The NDE results reported here are related neither to the glass transition phenomenon nor to the recently developed concept of a second liquid-liquid transition.

Refractive index of liquid water in different solvent models

We present a combined molecular dynamics/quantum chemical perturbation method for calculating the refractive index of liquid water at different temperatures. We compare results of this method with the refractive index obtained from other solvent models. The best agreement with the experimental refractive index of liquid water and its temperature dependence is obtained using correlated gas-phase polarizabilities in the classical Lorentz-Lorenz expression. Also, the iterative self-consistent reaction field approach in the semicontinuum implementation matches the experimental refractive index reasonably well.

Computer Simulation of the Linear and Nonlinear Optical Susceptibilities of p-Nitroaniline in Cyclohexane, 1,4-Dioxane, and Tetrahydrofuran in Quadrupolar Approximation. II. Local Field Effects and Optical Susceptibilitities

This is the second part of a study to elucidate the local field effects on the nonlinear optical properties of p-nitroaniline (pNA) in three solvents of different multipolar character, that is, cyclohexane (CH), 1,4-dioxane (DI), and tetrahydrofuran (THF), employing a discrete description of the solutions. By the use of liquid structure information from molecular dynamics simulations and molecular properties computed by high-level ab initio methods, the local field and local field gradients on p-nitroaniline and the solvent molecules are computed in quadrupolar approximation. To validate the simulations and the induction model, static and dynamic (non)linear properties of the pure solvents are also computed. With the exception of the static dielectric constant of pure THF, a good agreement between computed and experimental refractive indices, dielectric constants, and third harmonic generation signals is obtained for the solvents. For the solutions, it is found that multipole moments up to two orders higher than quadrupole have a negligible influence on the local fields on pNA, if a simple distribution model is employed for the electric properties of pNA. Quadrupole effects are found to be nonnegligible in all three solvents but are especially pronounced in the 1,4-dioxane solvent, in which the local fields are similar to those in THF, although the dielectric constant of DI is 2.2 and that of the simulated THF is 5.4. The electric-field-induced second harmonic generation (EFISH) signal and the hyper-Rayleigh scattering signal of pNA in the solutions computed with the local field are in good to fair agreement with available experimental results. This confirms the effect of the "dioxane anomaly" also on nonlinear optical properties. Predictions based on an ellipsoidal Onsager model as applied by experimentalists are in very good agreement with the discrete model predictions. This is in contrast to a recent discrete reaction field calculation of pNA in 1,4-dioxane, which found that the predicted first hyperpolarizability of pNA deviated strongly from the predictions obtained using Onsager-Lorentz local field factors.

Molecular Dynamics Study of 2-Nitrophenyl Octyl Ether and Nitrobenzene

The pure organic liquids nitrobenzene (NB) and 2-nitrophenyl octyl ether (NPOE) have been studied by means of molecular dynamics simulations. Both solvents are extremely important in various interfacial processes, mainly connected with ion transfer taking place across the interface with water. Thermodynamic (mass density, enthalpy of vaporization, isothermal compressibility, dipole moment) and dynamic (viscosities and self-diffusion coefficients) properties of both liquids have been calculated and are in very good agreement with the experimental data. In the case of NB, several potentials have been tested and the obtained results compared and discussed. In most cases, the OPLS all-atom potential gives results that are in better agreement with available experimental values. Atomic radial distribution functions, dihedral and angle distributions, as well as dipole-orientation correlation functions are used to probe the structure and interactions of the bulk molecules of both organic solvents. These were seen to be very similar in terms of structure and thermodynamics, but quite distinct in terms of dynamic behavior, with NPOE showing a much slower dynamic response than NB. A simulation study of the simple Cl- and K+ ions dissolved in both solvents has been also undertaken, revealing details about the diffusion and solvation mechanisms of these ions. It was found that in both liquids the positive potassium ion is solvated by the negative end of the molecular dipole, whereas the negative chloride ion is solvated by the positive end of the dipole.

Ferroelectric domains in nitrobenzene-nitromethane solutions measured by hyper-Rayleigh scattering

Hyper-Rayleigh scattering (HRS) spectra were measured for liquid solutions of C6H5NO2 and CH3NO2 at T=300 K. The depolarized HRS spectra at small frequency shift are dominated by two components due to reorientation of the nitrobenzene molecules. One is a Lorentzian with spectral width nu1=0.16-0.45 cm(-1) and corresponding orientation relaxation time tau=33-12 ps. The second component is a narrow spike with spectral width <2 MHz and corresponding relaxation time tau>80 ns, attributed to HRS from slowly relaxing ferroelectric domains. The dipole order parameter g0=0.053+/-0.005, saturation parameter p=0.9+/-0.1, and volume V=20+/-6 nm3 for these domains in nitromethane were determined from measurements of the nitrobenzene-concentration dependence of the intensity ratio for these two spectral components. Orientation of the 230 nitromethane molecules within each domain is inhomogenous but highly ordered.

Optimization of Nonlinear Optical Properties by Substituent Position, Geometry and Symmetry of the Molecule: An ab Initio Study

Static polarizability and first- and second-order hyperpolarizability tensors are computed at the correlated level for a series of para-nitroaniline derivatives. The importance of including electron correlation effects in the determination of equilibrium structure and the molecular properties is investigated. A qualitative description of the substitution effects, planarity, and symmetric effect of the molecule on the molecular susceptibility is discussed.

Study of intermolecular interactions in liquid nitrobenzene by depolarized hyper-Rayleigh scattering

We employed depolarized hyper-Rayleigh scattering (HRS) to investigate the intermolecular interactions in liquid nitrobenzene. By comparing the depolarization ratios of the second-harmonic scattered light from neat nitrobenzene and mixtures of nitrobenzene and methanol of varying mixing ratios, we demonstrated the existence of a coherent component of HRS in liquid nitrobenzene. The coherent component was found to essentially disappear at a sufficiently high dilution of the nitrobenzene liquid. We also observed that both localized orientational correlation and delocalized libron excitation contribute to coherent HRS in liquid nitrobenzene. The delocalized contribution to coherent HRS was found to diminish much more readily with the introduction of interstitial foreign molecules than the localized contribution.