Molecular Dipole Polarizabilities Research Articles

Ambiguities in Decomposing Molecular Polarizability into Atomic Charge Flow and Induced Dipole Contributions.

The molecular dipole polarizability can be decomposed into components corresponding to the charge flow between atoms and changes in atomic dipole moments. Such decompositions are recognized to depend on how atoms are defined within a molecule, as, for example, by Hirshfeld, iterative Stockholder, or quantum topology partitioning of the electron density. For some of these, however, there are significant differences between the numerical results obtained by analytical response methods and finite field calculations. We show that this difference is due to analytical response methods accounting for (only) the change in electron density by a perturbation, while finite field methods may also include a component corresponding to a perturbation-dependent change in the definition of an atom within a molecule. For some atom-in-molecule definitions, such as the iterative Hirshfeld, iterative Stockholder, and quantum topology methods, the latter effect significantly increases the charge flow component. The decomposition of molecular polarizability into atomic charge flow and induced dipole components thus depends on whether the atom-in-molecule definition is taken to be perturbation-dependent.

Exploring nonlinear and linear optical properties of N, N’-bis(salicylidene)-o-phenylenediamine (salophen) through Z-scan technique and TD-DFT computational analysis

This investigation delves into the comprehensive analysis of the linear and nonlinear optical characteristics exhibited by N, N’-bis(salicylidene)-o-phenylenediamine (salophen) using a combination of Z-scan methodology and quantum chemistry calculations. The Z-scan technique facilitated meticulous computations of the third-order susceptibility, nonlinear absorption coefficient, and nonlinear refractive index of the specimen. Specifically, in the solvent DMSO, the assessed values for the nonlinear refractive index, nonlinear absorption coefficient, and third-order susceptibility were determined to be 0.035 × 10–10 cm2 W −1 , −0.024 × 10–5 cm W −1 , and 0.596 × 10–5 esu, respectively. Furthermore, quantum mechanical analyses were employed to meticulously calculate the molecular hyperpolarizabilities (β and γ), dipole moment (μ), and dipole polarizability (α) of salophen. This thorough exploration highlighted a notable congruence between the outcomes derived from experimental observations and those obtained through quantum mechanical simulations. The collective findings from both theoretical computations and experimental assessments distinctly showcase the robust nonlinear potential inherent in salophen. These insights suggest its promising suitability and potential as a viable candidate for applications in optical devices. The alignment between theoretical predictions and experimental results underlines the reliability and potential practicality of salophen in the realm of optical technology, emphasizing its significance as a potential material for advancing optical device functionalities.

Data-driven tailoring of molecular dipole polarizability and frontier orbital energies in chemical compound space.

Understanding correlations - or lack thereof - between molecular properties is crucial for enabling fast and accurate molecular design strategies. In this contribution, we explore the relation between two key quantities describing the electronic structure and chemical properties of molecular systems: the energy gap between the frontier orbitals and the dipole polarizability. Based on the recently introduced QM7-X dataset, augmented with accurate molecular polarizability calculations as well as analysis of functional group compositions, we show that polarizability and HOMO-LUMO gap are uncorrelated when considering sufficiently extended subsets of the chemical compound space. The relation between these two properties is further analyzed on specific examples of molecules with similar composition as well as homooligomers. Remarkably, the freedom brought by the lack of correlation between molecular polarizability and HOMO-LUMO gap enables the design of novel materials, as we demonstrate on the example of organic photodetector candidates.

Range-dependence of two-body intermolecular interactions and their energy components in molecular crystals.

Routinely assessing the stability of molecular crystals with high accuracy remains an open challenge in the computational sciences. The many-body expansion decomposes computation of the crystal lattice energy into an embarrassingly parallel collection of computations over molecular dimers, trimers, and so forth, making quantum chemistry techniques tractable for many crystals of small organic molecules. By examining the range-dependence of different types of energetic contributions to the crystal lattice energy, we can glean qualitative understanding of solid-state intermolecular interactions as well as practical, exploitable reductions in the number of computations required for accurate energies. Here, we assess the range-dependent character of two-body interactions of 24 small organic molecular crystals by using the physically interpretable components from symmetry-adapted perturbation theory (electrostatics, exchange-repulsion, induction/polarization, and London dispersion). We also examine correlations between the convergence rates of electrostatics and London dispersion terms with molecular dipole moments and polarizabilities, to provide guidance for estimating convergence rates in other molecular crystals.

Permutationally invariant polynomial representation of polarizability tensor surfaces for linear regression analysis.

A linearly parameterized functional form for a Cartesian representation of molecular dipole polarizability tensor surfaces (PTS) is described. The proposed expression for the PTS is a linearization of the recently reported power series ansatz of the original Applequist model, which by construction is non-linear in parameter space. This new approach possesses (i) a unique solution to the least-squares fitting problem; (ii) a low level of the computational complexity of the resulting linear regression procedure, comparable to those of the potential energy and dipole moment surfaces; and (iii) a competitive level of accuracy compared to the non-linear PTS model. Calculations of CH4 PTS, with polarizabilities fitted to 9000 training set points with the energies up to 14,000 cm-1 show an impressive level of accuracy of the linear PTS model obtained with ~1600 parameters: ~1% versus 0.3% RMSE for the non-linear vs. linear model on a test set of 1000 configurations.

Deep Neural Network Model to Predict the Electrostatic Parameters in the Polarizable Classical Drude Oscillator Force Field.

The Drude polarizable force field (FF) captures electronic polarization effects via auxiliary Drude particles that are attached to non-hydrogen atoms, distinguishing it from commonly used additive FFs that rely on fixed charges. The Drude FF currently includes parameters for biomolecules such as proteins, nucleic acids, lipids, and carbohydrates and small-molecule representative of those classes of molecules as well as a range of atomic ions. Extension of the Drude FF to novel small druglike molecules is challenging as it requires the assignment of partial charges, atomic polarizabilities, and Thole scaling factors. In the present article, deep neural network (DNN) models are trained on quantum mechanical (QM)-based partial charges and atomic polarizabilities along with Thole scale factors trained to target QM molecular dipole moments and polarizabilities. Training of the DNN model used a collection of 39 421 molecules with molecular weights up to 200 Da and containing H, C, N, O, P, S, F, Cl, Br, or I atoms. The DNN model utilizes bond connectivity, including 1,2, 1,3, 1,4, and 1,5 terms and distances of Drude FF atom types as the feature vector to build the model, allowing it to capture both local and nonlocal effects in the molecules. Novel methods have been developed to determine restrained electrostatic potential (RESP) charges on atoms and external points representing lone pairs and to determine Thole scale factors, which have no QM analogue. A penalty scheme is devised as a performance predictor of the trained model. Validation studies show that these DNN models can precisely predict molecular dipole and polarizabilities of Food and Drug Administration (FDA)-approved drugs compared to reference MP2 calculations. The availability of the DNN model allowing for the rapid estimation of the Drude electrostatic parameters will facilitate its applicability to a wider range of molecular species.

On the Cartesian Representation of the Molecular Polarizability Tensor Surface by Polynomial Fitting to Ab Initio Data.

We describe an approach to constructing an analytic Cartesian representation of the molecular dipole polarizability tensor surface in terms of polynomials in interatomic distances with a training set of ab initio data points obtained from a molecular dynamics (MD) simulation or by any other available means. The proposed formulation is based on a perturbation treatment of the unmodified point dipole polarizability model of Applequist [ J. Am. Chem. Soc. 1972, 94, 2952] and is shown here to be, by construction (i) free of short-range or other singularities or discontinuities, (ii) symmetric and translationally invariant, and (iii) nonreliant on a body-fixed coordinate system. Permutational invariance of like nuclei is demonstrated to be readily applicable, making this approach useful for highly fluxional and reactive systems. Derivation of the method is described in detail, adding brief didactic numerical examples of H2 and H2O and concluding with an MD simulation of the Raman spectrum of H5O2+ at 300 K with the polarizability tensor fitted to CCSD(T)/aug-cc-pVTZ data obtained using the HBB-4B potential [ J. Chem. Phys. 2005, 122, 044308].

Comparative analysis of nonrelativistic and relativistic calculations of electric dipole moments and polarizabilities of heteronuclear alkali-metal dimers

We analyze the molecular electric dipole moments (PDMs) and static electric dipole polarizabilities of heteronuclear alkali dimers in their ground states by employing coupled-cluster theory, both in the non-relativistic and four-component relativistic frameworks. The roles of electron correlations as well as relativistic effects are demonstrated by studying them at different levels of theory, followed by a comprehensive treatment of error estimates. We compare our obtained values with the previous non-relativistic calculations, some of which include lower-order relativistic corrections, as well as with the experimental values, wherever available. We find that the PDMs are very sensitive to relativistic effects, as compared to polarizabilities; this aspect can explain the long-standing question on the difference between experimental values and theoretical results for LiNa. We show that consideration of relativistic values of PDMs improves significantly the isotropic Van der Waals $C_6$ coefficients of the investigated alkali dimers over the previously reported non-relativistic calculations. The dependence of dipole polarizabilities on molecular volume is also illustrated.

Ab initio calculations of permanent dipole moments and dipole polarizabilities of alkaline-earth monofluorides

The ground-state permanent dipole moments (PDMs) and molecular dipole polarizabilities (DPs) of open-shell alkaline-earth monofluorides, and atomic DPs of alkaline-earth- and fluorine atoms are reported at the Kramers-restricted configuration interaction level of theory limited to single and double excitations (KRCISD), using the finite-field approach. Sufficiently large basis sets such as quadrupole-zeta (QZ) and augmented-QZ basis sets together with the generalized active space technique is employed to carry out the field-dependent energy calculations at the KRCISD level. The PDMs and the components of DPs are extracted from the linear- and quadratic fit of energies against perturbative electric field, respectively. Accuracy of the present calculations for the electronic properties is examined by comparison with the measurements and calculations where ever available.

Accurate molecular polarizabilities with coupled cluster theory and machine learning

The molecular dipole polarizability describes the tendency of a molecule to change its dipole moment in response to an applied electric field. This quantity governs key intra- and intermolecular interactions, such as induction and dispersion; plays a vital role in determining the spectroscopic signatures of molecules; and is an essential ingredient in polarizable force fields. Compared with other ground-state properties, an accurate prediction of the molecular polarizability is considerably more difficult, as this response quantity is quite sensitive to the underlying electronic structure description. In this work, we present highly accurate quantum mechanical calculations of the static dipole polarizability tensors of 7,211 small organic molecules computed using linear response coupled cluster singles and doubles theory (LR-CCSD). Using a symmetry-adapted machine-learning approach, we demonstrate that it is possible to predict the LR-CCSD molecular polarizabilities of these small molecules with an error that is an order of magnitude smaller than that of hybrid density functional theory (DFT) at a negligible computational cost. The resultant model is robust and transferable, yielding molecular polarizabilities for a diverse set of 52 larger molecules (including challenging conjugated systems, carbohydrates, small drugs, amino acids, nucleobases, and hydrocarbon isomers) at an accuracy that exceeds that of hybrid DFT. The atom-centered decomposition implicit in our machine-learning approach offers some insight into the shortcomings of DFT in the prediction of this fundamental quantity of interest.

On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions.

In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C 6 dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H2O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α iso) and nearly 50% to the anisotropic (α aniso) polarizability values. In doing so, we find that α iso = 9.8307 bohr3, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.

Effect of the Molecular Polarizability of SAMs on the Work Function Modification of Gold: Closed‐ versus Open‐Shell Donor–Acceptor SAMs

AbstractCharge injection barriers at metal/organic interfaces can be tuned by modifying the work function of metallic electrodes using self‐assembled monolayers (SAMs) of polar molecules. An interesting example of polar molecules is offered by donor–acceptor (D–A) dyads based on ferrocene (Fc) as electron‐donor unit and either a polychlorotriphenylmethyl radical or a polychlorotriphenylmethane as electron‐acceptor units, connected by a π‐conjugated vinylene bridge. The D–A radical exhibits high chemical and thermal stability and presents different electronic, optical, and magnetic properties with respect to the closed‐shell form. The magnitude of the shift in the charge injection barriers for these two D–A systems is estimated by means of surface potential measurements performed by Kelvin probe force microscopy. The experimental data are compared with density functional theory calculations, which evidence the importance of the molecular dipole moments and polarizabilities to understand the experimental values. In order to achieve high work function shifts of metals upon SAM formation, the molecules forming the SAM have to exhibit both a high permanent dipole moment and a low polarizability along the direction normal to the substrate. In presence of polarizable molecules, the work function shifts can be enhanced by reducing the intermolecular interactions; by using mixed SAMs with active molecules embedded into a passive matrix.

Optimized geometry, spectroscopic characterization and nonlinear optical properties of carbazole picrate: a density functional theory study

The molecular modeling of carbazole picrate (CP) was carried out by using B3LYP and HSEH1PBE levels of density functional theory and 6-311++G(d,p) basis set by means of Gaussian 09 revision D.01 program. These methods have been used to determine the optimized molecular geometries, vibrational frequencies, electronic transitions and bonding features of the title compound. The computed small energy gap between HOMO and LUMO energies shows that the charge transfer occurs within the investigated compound. Additionally, the intensive interactions characterized by high stabilization energies were the powerful indicators of intra- and intermolecular charge transfer interactions. The obtained molecular dipole moment (μ), polarizability ( $$\left\langle \alpha \right\rangle$$ ) and hyperpolarizability ( $$\left\langle \beta \right\rangle$$ ) indicates that CP exhibits considerable nonlinear optical characteristic. The theoretical structural parameters such as bond lengths and bond angles are in a good agreement with the experimental values of the title compound. Additionally, the hydrogen bonding interactions were visualized via molecular electrostatic potential surface.

Analytical Insight into the Effect of Electric Field on Molecular Properties of Homonuclear Diatomic Molecules

Background: The diatomic molecules represent one of the most active classes of materials, and they have been widely used as active materials for important applications. Objective: In this study, the effect of external electric field on molecular properties of four homonuclear diatomic molecules H2, C2, Cl2, and Br2 has been analyzed. Methods: All the calculations are based on the density functional theory (DFT) at the B3LYP/6- 311++G(3df,3pd) level through the Gaussian 09W program package. Results: The studied molecular properties include total energy, density of state (DOS) analysis, frontier molecular orbital energies HOMO and LUMO, energy gap, absorption threshold, internuclear separation, harmonic vibrational frequency, bond force constant, molecular dipole moment, polarizability, chemical potential, chemical hardness and electrophilicity. Conclusion: It has been observed that the interaction of molecules with external electric fields can be used to control and improve the basic properties of homonuclear diatomic molecules for their exploitation in some significant applications in spintronics. Keywords: B3LYP, DFT, diatomic molecules, energy gap, external field.

Modeling molecular boiling points using computed interaction energies.

The noncovalent van der Waals interactions between molecules in liquids are typically described in textbooks as occurring between the total molecular dipoles (permanent, induced, or transient) of the molecules. This notion was tested by examining the boiling points of 67 halogenated hydrocarbon liquids using quantum chemically calculated molecular dipole moments, ionization potentials, and polarizabilities obtained from semi-empirical (AM1 and PM3) and ab initio Hartree-Fock [HF 6-31G(d), HF 6-311G(d,p)], and density functional theory [B3LYP/6-311G(d,p)] methods. The calculated interaction energies and an empirical measure of hydrogen bonding were employed to model the boiling points of the halocarbons. It was found that only terms related to London dispersion energies and hydrogen bonding proved significant in the regression analyses, and the performances of the models generally improved at higher levels of quantum chemical computation. An empirical estimate for the molecular polarizabilities was also tested, and the best models for the boiling points were obtained using either this empirical polarizability itself or the polarizabilities calculated at the B3LYP/6-311G(d,p) level, along with the hydrogen-bonding parameter. The results suggest that the cohesive forces are more appropriately described as resulting from highly localized interactions rather than interactions between the global molecular dipoles.

Effect of an external electric field on the C[sbnd]N cleavage reactions in nitromethane and triaminotrinitrobenzene

Density functional theory calculations were performed to provide an explicit physical understanding on how an external electric field affects molecular reactivity by taking energetic molecules nitromethane and triaminotrinitrobenzene as examples. The activation energies of CN cleavage, which is fixed to parallel or antiparallel to the applied field, in these molecules decrease when the bonds are either shortened or lengthened by the applied field. Our calculations reveal that the coupling of molecular dipole moment and/or polarizability with the applied field has a stabilization effect on the activated bond. Such effect increases during the CN bond cleavage, resulting in a lower activation energy than that in absence of the applied field. Moreover, the reduction in bond activation energy is favored in conjugated structures.

Structural and spectroscopic investigation on antioxidant dipeptide, l-Methionyl-l-Serine: A combined experimental and DFT study

The focus of this study is to determine the conformational, structural and vibrational properties of Methionyl-Serine dipeptide (l-Methionyl-l-Serine, Met-Ser), a biological active molecule. To investigate their energetically preferred conformations, molecular mechanics methods were utilized to determine the optimal conformations of the 3402 different dihedral angle values of the backbone and side chains. It was found that the extended (e) backbone shape in the LB conformational range was the most stable l-Methionyl-l-Serine dipeptide conformation, with 3.12 kcal/mol of energy. Density Functional Theory (DFT) was used to determine the optimized geometry, the vibrational wavenumbers and modes of the title dipeptide values, with 6-31G (d,p) and 6-311++G (d,p) basis sets. The potential energy distribution data was used to carry out the assignment of the bands. In addition, the vibrational spectra of the most stable conformer and its dimer form were determined and the obtained results were compared with the experimental IR and Raman spectra in the solid phase. To determine the presence of intramolecular charge transfer, molecular dipole moment, polarizability and hyperpolarizability, the Natural Bond Orbital (NBO), HOMO–LUMO calculations, the linear polarizability (α) and the first order hyperpolarizability (β0) value analyses of the investigated molecule were carried out using the DFT with the B3LYP/6-31++G(d,p) basis set. This study aims to determine a relatively stable conformation of antioxidant dipeptide and to investigate the molecular geometry, molecular vibrations and hydrogen bonding interactions between monomeric and dimeric forms of Methiony-Serine dipeptide.

Direct Reduction of Alkyl Monohalides at Silver in Dimethylformamide: Effects of Position and Identity of the Halogen

AbstractElectrochemical reduction of primary (1‐halodecane), secondary (2‐halohexane and cyclohexyl halide), and tertiary halides (tert‐butyl halide) at a silver cathode in dimethylformamide (DMF) containing 0.050 M tetramethylammonium perchlorate (TMAP) was investigated with the aid of cyclic voltammetry and controlled‐potential (bulk) electrolysis for iodides, bromides, and chlorides. Selected reductions were probed in dried DMF–TMAP and DMF containing tetra‐n‐hexylammonium perchlorate (THAP). Cyclic voltammograms reveal that the number of cathodic peaks and their potentials are highly dependent on the identity and position of the halogen as well as the composition of the solvent–electrolyte. Intermediates arising from bulk electrolysis of these halides at silver undergo both radical and carbanion reactions. Electrolysis products depend on the identity and position of the halogen, amount of residual water, and size of the electrolyte cation in the solvent–electrolyte. Theoretical calculations of molecular dipole moments and polarizabilities were performed and compared to experimental observations.

Hyperpolarizability calculation and kinetic effect of impurities on LVP

In order to investigate microscopic second order nonlinear optical properties of L-valinium picrate (LVP) crystals, the molecular dipole moment, polarizability, and first hyperpolarizability are estimated by density functional theory (DFT). Apart from that, a quantitative study using Atomic force microscopy (AFM) combined with crystal growth experiments on the effects of impurities (magnesium and copper) on LVP is presented. Two-dimensional (2D) nucleation growth mechanism and typical surface morphologies are described and discussed. In the presence of impurities, the decreasing of the size of 2D islands is probably incurred by the reducing of the step kinetic coefficient. From the supersaturation dependence of the normal growth rate, the effective step edge free energies are evaluated as γ0=2.90kT, γMg=2.76kT, γCu=2.95kT.

Polarizability of neutral copper clusters.

Neutral copper clusters are characterized through their molecular structures, binding energy and electric dipole polarizability. It is shown that the mean and anisotropy of polarizability tensor are useful properties in the characterization and rationalization of reactivity and growth patterns in copper clusters. We also found a relationship between softness per atom and cubic root polarizability per atom which can be useful to get global softness in copper clusters.