Molecular dynamics simulations applied to electric field induced second harmonic generation in dipolar chromophore solutions

Electric field induced second harmonic generation (EFISH) is an important experimental technique in extracting the first hyperpolarizability of an organic chromophore molecule. Such experiments are carried out in solutions with chromophore molecules dissolved in some common solvents. A known fact is that the first hyperpolarizabilities extracted from EFISH experiments are subject to the use of local field factors. In this work, we apply simulations to study the EFISH properties of chromophore solutions. By combining quantum chemistry calculations with the results derived from molecular dynamics simulations, we show how macroscopic EFISH properties can be modeled, using 4-(dimethylamino)-4'-nitroazobenzene dissolved in chloroform as a demonstration case. The focus of the study is on deriving accurate local field factors. We find that the local field approach applies very well to dipolar solutions, such as the one studied here, but that the local field factors derived are much smaller than the commonly used Onsager or Lorentz local field factors. Our study indicates that many of the reported first hyperpolarizabilities for dipolar molecules from EFISH experiments are most probably underestimated because the Onsager/Lorentz approach, commonly used in extracting the molecular first hyperpolarizability, neglects the effects of the shapes of dipolar chromophore molecules on the local field factors.     

Molecular theory on dielectric constant at interfaces: a molecular dynamics study of the water/vapor interface

Though the local dielectric constant at interfaces is an important phenomenological parameter in the analysis of surface spectroscopy, its microscopic definition has been uncertain. Here, we present a full molecular theory on the local field at interfaces with the help of molecular dynamics simulation, and thereby provide microscopic basis for the local dielectric constant so as to be consistent to the phenomenological three-layer model of interface systems. To demonstrate its performance, we applied the theory to the water/vapor interface, and obtained the local field properties near the interface where the simple dielectric model breaks down. Some computational issues pertinent to Ewald calculations of the dielectric properties are also discussed.     

First principles molecular dynamics simulation of liquid rubidium

We present a first principles molecular dynamics simulation of liquid rubidium. The atomic forces are obtained from a quantum mechanical calculation of its electronic structure within the local density approximation of the density functional formalism and using the pseudopotential plane-wave method. We compare our results with the structure and dynamics predicted by classical molecular dynamics simulations.     

Simulation of liquid imidazole using a high-rank quantum topological electrostatic potential

Rigid body molecular dynamics simulations were carried out on pure liquid imidazole at four different temperatures and at 1 atm. Imidazole, which is important both in life science and materials science, is one of the simplest molecules to possess both a lone pair and a π system. These two features are known to benefit from multipolar electrostatics. Here the electrostatic interaction is governed by atomic multipole moments obtained from topologically partitioned ab initio electron densities. The non-electrostatic terms are modeled with Lennard-Jones parameters adjusted to fit the experimental liquid density. All σ values are incrementally increased by one single scaling factor. We report on how the presence of multipolar electrostatics influences the local structure, dynamics and thermodynamics of the liquid compared to electrostatics by atomic point charges. The point charge force field exaggerates the number of π-stacked dimers in the liquid, and underestimates the number of hydrogen-bonded dimers. The effect of the temperature on the local structure of liquid imidazole was analysed using radial and spatial distribution functions.     

Solvent effects on electronic properties from Wannier functions in a dimethyl sulfoxide/water mixture

We present an efficient implementation for the calculation of maximally localized Wannier functions (MLWFs) during parallel Car-Parrinello molecular dynamics simulations. The implementation is based on a block Jacobi method. The calculation of MLWFs results in only a moderate (10%-20%) increase in computer time. Consequently it is possible to calculate MLWFs routinely during Car-Parrinello simulations. The Wannier functions are then applied to derive molecular dipole moments of dimethyl sulfoxide (DMSO) in gas phase and aqueous solution. We observe a large increase of the local dipole moment from 3.97 to 7.39 D. This large solvent effect is caused by strong hydrogen bonding at the DMSO oxygen atom and methyl groups. Decomposing the dipole moment into local contributions from the S-O bond and the methyl groups is used to understand the electrostatic response of DMSO in aqueous solution. A scheme is given to derive charges on individual atoms from the MLWFs using the D-RESP methodology. The charges also display large solvent effects and give insight into the transferability of recent force field models for DMSO.     

Quantum mechanical computations on very large molecular systems: The local self-consistent field method

Quantum chemical computations on a subset of a large molecule can be performed, at the neglect of diatomic differential overlap (NDDO) level, without further approximation provided that the atomic orbitals of the frontier atoms are replaced by parametrized orthogonal hybrid orbitals. The electrostatic interaction with the rest of the molecule, treated classically by the usual molecular mechanical approximations, is included into the self-consistent field (SCF) equations. The first and second derivatives of energy are obtained analytically, allowing the search for energy minima and transition states as well as the resolution of Newton equations in molecular dynamics simulations. The local self-consistent field (LSCF) method based on these approximations is tested by studying the intramolecular proton transfer in a Gly-Arg-Glu-Gly model tetrapeptide, which reveals an excellent agreement between a computation performed on the whole molecule and the results obtained by the present method, especially if the quantum subsystem includes the side chains and the peptidic unit in between. The merits of the LSCF method are examplified by a study of proton transfer in the Asp⁶⁹—Arg⁷¹ salt bridge in dihydrofolate reductase. Simulations of large systems, involving local changes of electronic structure, are therefore possible at a good degree of approximation by introducing a quantum chemical part in molecular dynamics studies. This methodology is expected to be very useful for reactivity studies in biomolecules or at the surface of covalent solids. © 1994 by John Wiley & Sons, Inc.     

Pressure calculation in hybrid particle-field simulations

In the framework of a recently developed scheme for a hybrid particle-field simulation techniques where self-consistent field (SCF) theory and particle models (molecular dynamics) are combined [J. Chem. Phys. 130, 214106 (2009)], we developed a general formulation for the calculation of instantaneous pressure and stress tensor. The expressions have been derived from statistical mechanical definition of the pressure starting from the expression for the free energy functional in the SCF theory. An implementation of the derived formulation suitable for hybrid particle-field molecular dynamics-self-consistent field simulations is described. A series of test simulations on model systems are reported comparing the calculated pressure with those obtained from standard molecular dynamics simulations based on pair potentials.     

Dielectric permittivity profiles of confined polar fluids

The dielectric response of a simple model of a polar fluid near neutral interfaces is examined by a combination of linear response theory and extensive molecular dynamics simulations. Fluctuation expressions for a local permittivity tensor epsilon(r) are derived for planar and spherical geometries, based on the assumption of a purely local relationship between polarization and electric field. While the longitudinal component of epsilon exhibits strong oscillations on the molecular scale near interfaces, the transverse component becomes ill defined and unphysical, indicating nonlocality in the dielectric response. Both components go over to the correct bulk permittivity beyond a few molecular diameters. Upon approaching interfaces from the bulk, the permittivity tends to increase, rather than decrease as commonly assumed, and this behavior is confirmed for a simple model of water near a hydrophobic surface. An unexpected finding of the present analysis is the formation of "electrostatic double layers" signaled by a dramatic overscreening of an externally applied field inside the polar fluid close to an interface. The local electric field is of opposite sign to the external field and of significantly larger amplitude within the first layer of polar molecules.     

Microstructure of neat alcohols: a molecular dynamics study

Neat methanol and tert-butanol are studied by molecular dynamics with the focus on the microstructure of these two alcohols. The site-site radial distribution functions, the corresponding structure factors, and an effective local one-body density function are shown to be the appropriate statistical quantities that point in a complementary manner towards the same microstructure for any given liquid. Methanol is found to be a weakly associated liquid forming various chainlike patterns (open and closed) while tert-butanol is almost entirely associated and forms micellelike primary pattern. The presence of stable local microheterogeneity within homogeneous disordered phase appears as a striking feature of these liquids. The absence of any such apparent clustering in water--a stronger hydrogen bonding liquid--through the same two statistical quantities is analyzed.     

Simulations of lipid crystals: Characterization of potential energy functions and parameters for lecithin molecules

We evaluate an empirical potential energy function and associated parameters for classical molecular dynamics simulations of lecithins, a common class of lipid. The physical accuracy of the force field was tested through its application to molecular dynamics simulations of the known crystal structures of lipid molecules. Average atomic positions and molecular conformation are well maintained during the simulations despite considerable thermal motion. Calculated isotropic temperature factors correlate highly with those from experiment.     

Statistical mechanical theory for steady state systems. IV. Transition probability and simulation algorithm demonstrated for heat flow

Two microscopic transition theorems are given for the probability of nonequilibrium work performed on a subsystem of a thermal reservoir along the trajectory in phase space of the subsystem. The resultant transition probability is applied to the case of heat flow down an applied temperature gradient. A combined molecular dynamics and Monte Carlo algorithm is given for such a nonequilibrium steady state. Results obtained for the thermal conductivity are in good agreement with previous Green-Kubo and nonequilibrium molecular dynamics results.     

A comparative study of the properties of polar and nonpolar solvent/solute/polystyrene solutions in microwave fields via molecular dynamics

The influence of an applied microwave field on the dynamics of methylamine-dichloromethane (DCM) mixtures bound within atactic polystyrene (a-PS) over a range of polymer densities from 30 to 94 wt % polymer was examined using atomistic molecular dynamics simulations. This study is an extension of previous studies on methylamine transport in relatively polar polystyrene solutions of methanol and dimethylformamide [M. J. Purdue et al., J. Chem. Phys. 124, 204904 (2006)]. A direct comparison is made across the three types of polystyrene solutions. Consideration is given to both solvent and reagent transport within the polymer solutions under zero-field conditions and in an external electromagnetic field in the canonical ensemble (NVT) at 298.0 K. Various frequencies up to 10(4) GHz and a rms electric field intensity of 0.1 VA were applied. The simulation studies were validated by comparison of the simulated zero-field self-diffusion coefficients of DCM in a-PS with those obtained using pulsed-gradient spin-echo NMR spectrometry. Athermal effects of microwave fields on solute transport behavior within polymer solutions are discussed.     

Calculation of the microscopic and macroscopic linear and nonlinear optical properties of liquid acetonitrile. II. Local fields and linear and nonlinear susceptibilities in quadrupolar approximation

A discrete model based on the multipolar expansion including terms up to hexadecapoles was employed to describe the electrostatic interactions in liquid acetonitrile. Liquid structures obtained form molecular dynamics simulations with different classical, nonpolarizable potentials were used to analyze the electrostatic interactions. The computed average local field was employed for the determination of the environmental effects on the linear and nonlinear electrical molecular properties. Dipole-dipole interactions yield the dominant contribution to the local field, whereas higher multipolar contributions are small but not negligible. Using the effective in-phase properties, macroscopic linear and nonlinear susceptibilities of the liquid were computed. Depending on the partial charges describing the Coulomb interactions of the force field employed, either the linear properties (refractive index and dielectric constant) were reproduced in good agreement with experiment or the nonlinear properties [third-harmonic generation (THG) and electric field induced second-harmonic (EFISH) generation] and the bulk density but never both sets of properties together. It is concluded that the partial charges of the force fields investigated are not suitable for reliable dielectric properties. New methods are probably necessary for the determination of partial charges, which should take into account the collective and long-range nature of electrostatic interactions more precisely.     

Comparison of kinetic Monte Carlo and molecular dynamics simulations of diffusion in a model glass former

We present results from kinetic Monte Carlo (KMC) simulations of diffusion in a model glass former. We find that the diffusion constants obtained from KMC simulations have Arrhenius temperature dependence, while the correct behavior, obtained from molecular dynamics simulations, can be super-Arrhenius. We conclude that the discrepancy is due to undersampling of higher-lying local minima in the KMC runs. We suggest that the relevant connectivity of minima on the potential energy surface is proportional to the energy density of the local minima, which determines the "inherent structure entropy." The changing connectivity with potential energy may produce a correlation between dynamics and thermodynamics.     

Local density augmentation and dynamic properties of hydrogen-and non-hydrogen-bonded supercritical fluids: a molecular dynamics study

The local density inhomogeneities in neat supercritical fluids were investigated via canonical molecular dynamics simulations. The selected systems under investigation were the polar and hydrogen-bonded fluid methanol as well as the quadrupolar non-hydrogen-bonded carbon dioxide one. Effective local densities, local density augmentation, and enhancement factors were calculated at state points along an isotherm close to the critical temperature of each system (T(r)=1.03). The results obtained reveal strong influence of the polarity and hydrogen bonding upon the intensity of the local density augmentation. It is found that this effect is sufficiently larger in the case of the polar and associated methanol in comparison to those predicted for carbon dioxide. For both fluids the local density augmentation values are maximized in the bulk density region near 0.7rho(c), a result that is in agreement with experiment. In addition, the local density dynamics of each fluid were investigated in terms of the appropriate time correlation functions. The behavior of these functions reveals that the bulk density dependence of the local density reorganization times is very sensitive to the specific intermolecular interactions and to the size of the local region. Also, the estimated local density reorganization time as a function of bulk density of each fluid was further analyzed and successfully related to two different time-scale relaxation mechanisms. Finally, the results obtained indicate a possible relationship between the single-molecule reorientational dynamics and the local density reorganization ones.     

Spontaneous Polarization of Thin Water Films in Nonwettable Capillaries by the Data of Numerical Experiment

Abstract–Water films in the flat capillaries with ideal walls were simulated by molecular dynamics method. The profiles of local density and distribution functions of the dipole moment orientations, the directions of valence and H-H bonds were obtained. Profiles of local electric field and electric potential were calculated. Potential jumps at the water–vapor interface were estimated.     

Electromagnetic model and calculations of the surface-enhanced Raman-shifted emission from Langmuir-Blodgett films on metal nanostructures

We present a theoretical study of the electromagnetic contribution to surface-enhanced Raman scattering (SERS) from a Langmuir-Blodgett film close to a metal surface. This macroscopic dipolar model fully accounts for the Raman-shifted emission so that meaningful SERS (electromagnetic) enhancement factors that do not depend only on the local electromagnetic field enhancement at the pump frequency are defined. For a plane metal surface, analytical SERS enhancement factors that are consistent for all pump beam polarization and molecular orientation are obtained. In order to investigate SERS on complex nanostructured metal surfaces, we introduce this model into the formally exact, Green's theorem surface integral equation formulation of the scattered electromagnetic field. This formulation is thus employed to calculate numerically the near-field and far-field emissions at the Raman-shifted frequency for very rough, random nanostructured surfaces, with emphasis on the impact of collective processes for varying pump frequency and Raman shift. Our results reveal that the widely used |E|4 approximation tends to overestimate average SERS enhancement factors.     

Bidirectional molecular dynamics: Interpretation in terms of a modern formulation of classical mechanics

The validity and applicability of bidirectional molecular dynamics is shown in terms of modern classical mechanics. A simple interpretation of bidirectional molecular dynamics is given. This interpretation justifies an easy approach to start from an equilibrium configuration, to perform simulations in parallel in “forward” and “backward” time direction, and to combine the two obtained trajectories into a single one for proper evaluation of statistical quantities. Practical results, obtained for liquid ammonia, are presented. © 1996 by John Wiley & Sons, Inc.