Hirshfeld-based intrinsic polarizability density representations as a tool to analyze molecular polarizability

Origin-Independent Dynamic Polarizability Density from Coupled Cluster Response Theory

The calculation of the origin-independent density of the dynamic electric dipole polarizability, previously presented for uncorrelated and density functional theory (DFT)-based methods, has been developed and implemented at the coupled cluster singles and doubles (CCSD) level of theory. A pointwise analysis of polarizability densities calculated for a number of molecules at Hartree-Fock (HF) and CCSD clearly shows that the electron correlation effect is much larger than one would argue considering the integrated dipole electric polarizability alone. Large error compensations occur during the integration process, which hide fairly large deviations mainly located in the internuclear regions. The same is observed when calculated CCSD and B3LYP polarizability densities are compared, with the remarkable feature that positive/negative deviations between CCSD and HF reverse sign, becoming negative/positive when comparing CCSD to B3LYP.

Electronic current densities and origin-independent property densities induced by optical fields

The interaction of a molecule with optical fields is customarily interpreted by means of induced time-dependent electric polarizabilities, magnetizabilities and mixed electric-magnetic polarizabilities. In general, these properties can be rationalized by integrals of density functions formulated in terms of induced charge and current densities. In this perspective, we focus on what has been done so far at the theoretical level, and on what can be expected to be unveiled from the topological study of suitable density functions, endowed with the fundamental requirement of origin invariance. Densities characterized by such a property can be integrated all over the configuration space to obtain electric dipole polarizability and optical rotatory power. Corresponding maps visualize domains mainly involved in the molecular response. The diagonal components of origin-independent density tensor functions that, on integration, yield corresponding electric dipole polarizability tensor of benzene, naphthalene, phenanthrene and ovalene, have been computed, confirming the ubiquitous presence of counter-polarization regions in the proximity of the atomic nuclei. They are associated with toroidal electron currents, induced by time derivative of the electric field of impinging radiation. Electron (de)localization in these systems is readily observed and estimated. The optical rotation density of the carbonyl chromophore is studied in detail. Its essential feature is the separation in quadrants of alternating sign of density about the CO bond. The presence of an extrachromophoric perturbation determines asymmetry in the extension of the quadrant distribution, thus causing optical rotation.

Benchmark of a functional-group database for distributed polarizability and dipole moment in biomolecules

The extraction of functional-group properties in condensed phases is very useful for predicting material behaviors, including those of biomaterials. For this reason, computational approaches based on partitioning schemes have been developed aiming at rapidly and accurately estimating properties from chemically meaningful building blocks. A comprehensive database of group polarizabilities and dipole moments is useful not only to predict the optical properties of biomacromolecules but also to improve molecular force fields focused on simulating biochemical processes. In this work we benchmark a database of distributed polarizabilities and dipole moments for functional groups extracted from a series of polypeptides. This allows reconstruction of a variety of relevant chemical environments. The accuracy of our database was tested to predict the electro-optical properties of larger peptides and also simpler amino acids for which density functional theory calculations at the M06-HF/aug-cc-pVDZ level of theory was chosen as the reference. This approach is reasonably accurate for the diagonal components of the polarizability tensor, with errors not larger than 15-20%. The anisotropy of the polarizability is predicted with smaller efficacy though. Solvent effects were included explicitly by surrounding the database entries by a box of water molecules whose distribution was optimized using the CHARMM force field.

Origin-Independent Densities of Static and Dynamic Molecular Polarizabilities

The notion of the electric dipole polarizability density function of atoms and molecules has been considered. The current density induced by the time derivative of the electric field of monochromatic light allows for a new definition of the electric dipole polarizability density, which is translationally invariant. This translational invariance provides the physical meaning that was lacking in previous defined polarizability densities. The new polarizability density has been implemented at the TD-DFT level of theory. The origin independence has been proven in silico to hold regardless of the basis set size. Two emblematic molecules, i.e., CO and N₂, which in many respects display similar electric response, have been studied in detail. The substantial differences, which have been highlighted in the topology of the parallel and perpendicular polarizability density tensor components of CO and N₂, are grossly hidden by compensation, when integration is carried out to get the molecular properties.

Accurate Atom-Dipole Interaction Model for Prediction of Electro-optical Properties: From van der Waals Aggregates to Covalently Bonded Clusters

This work aims at the accurate estimation of the electro-optical properties of atoms and functional groups in organic crystals. A better understanding of the nature of building blocks and the way they interact with each other enables more efficient prediction of self-assembly, and thus physical properties in condensed matter. We propose a modified version of an atom-dipole interaction model that is based on atomic dipole moments calculated from the quantum theory of atoms in molecules. The method is very reliable for the prediction of various optical and electric properties in diverse chemical environments, ranging from hydrocarbon molecules bonded by dispersive interactions to polar rings connected by hydrogen bonds, or even polymeric structures whose monomers are covalently linked. Distributed polarizabilities and electrostatic potentials are compared to those obtained using a complete quantum-mechanical approach on finite-size aggregates. Our electrostatic approximation recovers isotropic polarizabilities with an accuracy of ca. 5 au and electrostatic potentials of ca. 0.05 au, even in the worst-case scenario in which polarization and charge-transfer effects are large. Functional groups are highly exportable, estimating the properties of small peptides and polyaromatics with a maximum deviation as low as ca. 15%.

One- and Two-Component Organogels Containing Cyanostilbene without any Auxiliary Substituents

Pyridyl acrylonitrile without traditional auxiliary groups form stable organogels in ethanol. The addition of a second non-gelating cyanostilbene component results in a more stable two-component gel. Single crystal X-ray data reveal the influence of C-H⋅ ⋅ ⋅N, C-H⋅ ⋅ ⋅π, and π-π interactions in the formation of organogels. The morphology of the xerogels was studied by using SEM, which showed the self-assembly of molecules to fibers and sheet-like structures, and phase differences upon the gel formation and the structural phase characterization was measured using powder XRD. Exposure of the organogels to acidic (TFA) vapors results in distinct color changes and loss of gelation properties, thus highlighting the potential of these gels in sensing. The results represent a rare example of two-component organogels using two different cyanostilbene units and show that functional two-component organogels can be formed by utilizing the synergistic effects of the individual components.

Potential Application of h-BNC Structures in SERS and SEHRS Spectroscopies: A Theoretical Perspective

In this work, the electronic and optical properties of hybrid boron-nitrogen-carbon structures (h-BNCs) with embedded graphene nanodisks are investigated. Their molecular affinity is explored using pyridine as model system and comparing the results with the corresponding isolated graphene nanodisks. Time-dependent density functional theory (TDDFT) analysis of the electronic excited states was performed in the complexes in order to characterize possible surface and charge transfer resonances in the UV region. Static and dynamic (hyper)polarizabilities were calculated with coupled-perturbed Kohn-Sham theory (CPKS) and the linear and nonlinear optical responses of the complexes were analyzed in detail using laser excitation wavelengths available for (Hyper)Raman experiments and near-to-resonance excitation wavelengths. Enhancement factors around 10³ and 10⁸ were found for the polarizability and first order hyperpolarizability, respectively. The quantum chemical simulations performed in this work point out that nanographenes embedded within hybrid h-BNC structures may serve as good platforms for enhancing the (Hyper)Raman activity of organic molecules immobilized on their surfaces and for being employed as substrates in surface enhanced (Hyper)Raman scattering (SERS and SEHRS). Besides the better selectivity and improved signal-to-noise ratio of pristine graphene with respect to metallic surfaces, the confinement of the optical response in these hybrid h-BNC systems leads to strong localized surface resonances in the UV region. Matching these resonances with laser excitation wavelengths would solve the problem of the small enhancement factors reported in Raman experiments using pristine graphene. This may be achieved by tuning the size/shape of the embedded nanographene structure.

Partition of optical properties into orbital contributions

A new tool to analyze the response property through the partition of nonlinear optical properties in terms of orbital contributions (PNOC), valuable in the assessment of the electronic structure methods in the NLOPs computations, is presented.

A Computational Study of the Interaction and Polarization Effects of Complexes Involving Molecular Graphene and C60 or a Nucleobases

A systematic analysis of the molecular structure, energetics, electronic (hyper)polarizabilities and their interaction-induced counterparts of C60 with a series of molecular graphene (MG) models, CmHn, where m = 24, 84, 114, 222, 366, 546 and n = 12, 24, 30, 42, 54, 66, was performed. All the reported data were computed by employing density functional theory and a series of basis sets. The main goal of the study is to investigate how alteration of the size of the MG model affects the strength of the interaction, charge rearrangement, and polarization and interaction-induced polarization of the complex, C60-MG. A Hirshfeld-based scheme has been employed in order to provide information on the intrinsic polarizability density representations of the reported complexes. It was found that the interaction energy increases approaching a limit of -26.98 kcal/mol for m = 366 and 546; the polarizability and second hyperpolarizability increase with increasing the size of MG. An opposite trend was observed for the dipole moment. Interestingly, the variation of the first hyperpolarizability is relatively small with m. Since polarizability is a key factor for the stability of molecular graphene with nucleobases (NB), a study of the magnitude of the interaction-induced polarizability of C84H24-NB complexes is also reported, aiming to reveal changes of its magnitude with the type of NB. The binding strength of C84H24-NB complexes is also computed and found to be in agreement with available theoretical and experimental data. The interaction involved in C60 B12N12H24-NB complexes has also been considered, featuring the effect of contamination on the binding strength between MG and NBs.

Establishing the pivotal role of local aromaticity in the electronic properties of boron-nitride graphene lateral hybrids

Within an attempt to unravel the conundrum of irregular bandgap variations in hybrids of white-graphene (hBN) and graphene (G) observed in both experiment and theory, strong proofs about the decisive role of aromaticity in their electronic properties are brought to light.