Search

We demonstrate that contemporary computational resources allow for accurate theoretical studies of systems matching recent advances in experimental helicene chemistry. Concerned with first-principles calculations of carbohelicenes, our work surpasses CH[12] as the largest system investigated to date and unravels trends in the electronic structure of the low-lying states of the homologous series. Utilizing a highly efficient implementation of linear response algorithms, we present electronic circular dichroism (CD) spectra of carbohelicenes ranging from CH[5] to CH[30] at the level of Kohn–Sham density-functional theory. Our results for a systematic increase in system size show the emergence of new CD bands that subsequently rise to intensities dominating the spectrum. The spectral band progressions exhibit a periodicity directly linked to the number of overlapping layers of conjugation. While our findings rectify the current understanding of the electronic structure of carbohelicenes, they also serve as a general call for caution regarding the extrapolation of trends from small system ranges. [ABSTRACT FROM AUTHOR]

We report on accurate and efficient calculations of vibrationally resolved emission spectra for oligothiophenes from anharmonic vibrational configuration interaction wave-function calculations in reduced vibrational spaces. These reduced spaces are chosen based on the independent mode displaced harmonic oscillator model. Good agreement with experiment is obtained for all-trans oligothiophenes with two to five rings also when employing only a few active modes. Vibrational modes incorporating inter-ring carbon–carbon stretches and a ring breathing mode are found to be the main players in the vibrational progression for the emission from the first excited electronic state for all investigated oligothiophene derivatives. The presented framework is here illustrated for oligothiophenes, but we have made no underlying system-dependent assumptions and believe it to become a valuable tool for the rational design of fluorescence biomarkers. [ABSTRACT FROM AUTHOR]

DNA binding modes of the stereoisomeric rotamers of two dithenylethene derivatives (DTE1 and DTE2) representing candidate molecular photoswitches of great promise for photopharmacology and nanotechnology have been identified and characterized in terms of their binding energies and electronic circular dichroism (CD) responses. In the open form, two binding modes are identified namely minor-groove binding of the lowest-energy conformer with an anti-parallel arrangement of methyl groups and major-groove double-intercalation of the P-enantiomers of an intermediate-state rotamer. Only the latter binding mode is found to be enantiomerically selective and expected to have an overall negative linear dichroism (LD) as observed in the experiment for DTE1 (Angew. Chem., Int. Ed., 2013, 52, 4393). In the closed form, the most favorable binding mode is found to be minor groove binding. Also this binding mode is found to be enantiomerically selective and for DTE1, it is the M-enantiomer that binds the strongest, showing a positive theoretical signature CD band in the long wavelength region with origin in pyridinium ligands. The theoretical CD spectrum is found to be in good agreement with the experimental one, which provides an indirect evidence for a correct identification of the binding mode in the closed form. [ABSTRACT FROM AUTHOR]

Despite being very well established in the field of electro-optics, ferroelectric liquid crystals so far lacked interest from a ferroelectric device perspective due to a typically high operating temperature, a modest remnant polarization and/or poor polarization retention. Here, we experimentally demonstrate how simple structural modification of a prototypical ferroelectric liquid-crystal benzene-1,3,5-trisamide (BTA) – introduction of branched-tail substituents – results in materials with a wide operating temperature range and a data retention time of more than 10 years in thin-film solution-processed capacitor devices at room temperature. The observed differences between linear- and branched-tail compounds are analyzed using density functional theory (DFT) and molecular dynamics (MD) simulations. We conclude that morphological factors like improved packing quality and reduced disorder, rather than electrostatic interactions or intra/inter-columnar steric hindrance, underlay the superior properties of the branched-tailed BTAs. Synergistic effects upon blending of compounds with branched and linear side-chains can be used to further improve the materials’ characteristics. [ABSTRACT FROM AUTHOR]

Ferroelectrics find broad applications, e.g. in non-volatile memories, but the switching kinetics in real, disordered, materials is still incompletely understood. Here, we develop an electrostatic model to study ferroelectric switching using 3D Monte Carlo simulations. We apply this model to the prototypical small molecular ferroelectric trialkylbenzene-1,3,5-tricarboxamide (BTA) and find good agreement between the Monte Carlo simulations, experiments, and molecular dynamics studies. Since the model lacks any explicit steric effects, we conclude that these are of minor importance. While the material is shown to have a frustrated antiferroelectric ground state, it behaves as a normal ferroelectric under practical conditions due to the large energy barrier for switching that prevents the material from reaching its ground state after poling. We find that field-driven polarization reversal and spontaneous depolarization have orders of magnitude different switching kinetics. For the former, which determines the coercive field and is relevant for data writing, nucleation occurs at the electrodes, whereas for the latter, which governs data retention, nucleation occurs at disorder-induced defects. As a result, by reducing the disorder in the system, the polarization retention time can be increased dramatically while the coercive field remains unchanged. [ABSTRACT FROM AUTHOR]

A very stable binding site for the interaction between a pentameric oligothiophene and an amyloid-β(1–42) fibril has been identified by means of non-biased molecular dynamics simulations. In this site, the probe is locked in an all-trans conformation with a Coulombic binding energy of 1200 kJ mol−1 due to the interactions between the anionic carboxyl groups of the probe and the cationic ε-amino groups in the lysine side chain. Upon binding, the conformationally restricted probes show a pronounced increase in molecular planarity. This is in line with the observed changes in luminescence properties that serve as the foundation for their use as biomarkers. [ABSTRACT FROM AUTHOR]

QM/MM and QM/QM/MM protocols are applied to the ab initio study of the three linear birefringences Kerr, Cotton–Mouton, and Jones, as shown by acetonitrile in the gas and pure liquid phases. The relevant first-order properties as well as linear, quadratic, and cubic frequency-dependent response functions were computed using time-dependent Kohn–Sham density-functional theory with use of the standard CAM-B3LYP functional. In the liquid phase, a series of room temperature (293.15 K) molecular dynamics snapshots were selected, for which averaged values of the observables were obtained at an optical wavelength of 632.8 nm. The birefringences were computed for electric and magnetic induction fields corresponding to the laboratory setup previously employed by T. Roth and G. L. J. A. Rikken in Phys. Rev. Lett., 2000, 85, 4478. Under these conditions, acetonitrile is shown to exhibit a weak Jones response―in fact roughly 6.5 times smaller than the limit of detection of the apparatus employed in the measurements mentioned above. A comparison is made with the corresponding gas-phase results and an assessment is made of the index of measurability, estimating the degree of overlap of the three birefringences in actual measurements. For acetonitrile, it is shown that this index is a factor of 3.6 and 6.7 larger than that of methylcyclopentadienyl-Mn-tricarbonyl and cyclohexadienyl-Fe-tricarbonyl, respectively―two compounds reported in Phys. Rev. Lett., 2000, 85, 4478 to exhibit a strong Jones signal. [ABSTRACT FROM AUTHOR]

A computational protocol for magneto-chiral dichroism and magneto-chiral birefringence dispersion is presented within the framework of damped response theory, also known as complex polarization propagator theory, at the level of time-dependent Hartree–Fock and time-dependent density functional theory. Magneto-chiral dichroism and magneto-chiral birefringence spectra in the (resonant) frequency region below the first ionization threshold of R-methyloxirane and l-alanine are presented and compared with the corresponding results obtained for both the electronic circular dichroism and the magnetic circular dichroism. The additional information content yielded by the magneto-chiral phenomena, as well as their potential experimental detectability for the selected species, is discussed. [ABSTRACT FROM AUTHOR]

X-ray absorption spectra of carbon, silicon, germanium, and sulfur compounds have been investigated by means of damped four-component density functional response theory. It is demonstrated that a reliable description of relativistic effects is obtained at both K- and L-edges. Notably, an excellent agreement with experimental results is obtained for L2,3-spectra―with spin–orbit effects well accounted for―also in cases when the experimental intensity ratio deviates from the statistical one of 2 : 1. The theoretical results are consistent with calculations using standard response theory as well as recently reported real-time propagation methods in time-dependent density functional theory, and the virtues of different approaches are discussed. As compared to silane and silicon tetrachloride, an anomalous error in the absolute energy is reported for the L2,3-spectrum of silicon tetrafluoride, amounting to an additional spectral shift of ∼1 eV. This anomaly is also observed for other exchange–correlation functionals, but it is seen neither at other silicon edges nor at the carbon K-edge of fluorine derivatives of ethene. Considering the series of molecules SiH4−XFX with X = 1, 2, 3, 4, a gradual divergence from interpolated experimental ionization potentials is observed at the level of Kohn–Sham density functional theory (DFT), and to a smaller extent with the use of Hartree–Fock. This anomalous error is thus attributed partly to difficulties in correctly emulating the electronic structure effects imposed by the very electronegative fluorines, and partly due to inconsistencies in the spurious electron self-repulsion in DFT. Substitution with one, or possibly two, fluorine atoms is estimated to yield small enough errors to allow for reliable interpretations and predictions of L2,3-spectra of more complex and extended silicon-based systems. [ABSTRACT FROM AUTHOR]

A computational benchmark study on X-ray absorption spectra of water has been performed by means of transition-potential density functional theory (TP-DFT), damped time-dependent density functional theory (TDDFT), and damped coupled cluster (CC) linear response theory. For liquid water, using TDDFT with a tailored CAM-B3LYP functional and a polarizable embedding, we find that an embedding with over 2000 water molecules is required to fully converge spectral features for individual molecules, but a substantially smaller embedding can be used within averaging schemes. TP-DFT and TDDFT calculations on 100 MD structures demonstrate that TDDFT produces a spectrum with spectral features in good agreement with experiment, while it is more difficult to fully resolve the spectral features in the TP-DFT spectrum. Similar trends were also observed for calculations of bulk ice. In order to further establish the performance of these methods, small water clusters have been considered also at the CC2 and CCSD levels of theory. Issues regarding the basis set requirements for spectrum simulations of liquid water and the determination of gas-phase ionization potentials are also discussed. [ABSTRACT FROM AUTHOR]

Infrared (IR) absorption and vibrational Raman spectra of a family of branched oligothiophenes have been determined experimentally as well as theoretically. The molecular spectra have been compared to those of the linear analogues, with identification made of spectral features due to structural properties that are valued in organic solar cell applications. The theoretical spectra have been obtained through a newly developed method in which individual conformer spectra, calculated at the time-dependent DFT level in this work, are weighted by statistics extracted from classical molecular dynamics trajectories. The agreement with experiment for the resulting averaged spectra is at least as good as, and often better than, what is observed for Boltzmann-weighted spectra. As the weights are available before the costly step of spectrum calculation, the method has the additional advantage of enabling efficient approximations. For simulating the molecular dynamics of the studied α,β-linked thiophenes and 2-methylthiophenes, high quality parameters have been derived for the CHARMM force field. Furthermore, the temperature dependence of the IR and Raman spectra has been investigated, both experimentally and theoretically. [ABSTRACT FROM AUTHOR]

Based on resonant-convergent—and thereby complex—electric-dipole polarizabilities, formulae for differential cross sections in resonant Raman spectroscopy are presented. In absorptive regions of the spectrum, there are phase delays in the induced polarization that enter the expressions for the cross sections. The theory is exemplified by calculations employing the adiabatic time-dependent density functional theory approximation, with applications on the low-lying, dipole-allowed, electronic states of a sample of molecules, including hydrogen peroxide, pyridine, and trinitrotoluene. Results obtained with the Coulomb attenuated B3LYP exchange–correlation functional are found to be more accurate than those obtained with the conventional B3LYP functional—an observation that holds not only for trinitrotoluene with its transitions of charge-transfer character but also for the other cases. The qualitative features of the resonant Raman spectra for a given molecule vary strongly from one resonance wavelength to another, which is a fact that could further facilitate the use of this spectroscopy in applications of stand-off detection of gaseous samples in ultra-low concentrations. [ABSTRACT FROM AUTHOR]

It is demonstrated that time-dependent density functional theory (DFT) calculations can accurately predict changes in near-UV electronic circular dichroism (ECD) spectra of DNA as the structure is altered from the linear (free) B-DNA form to the supercoiled N-DNA form found in nucleosome core particles. At the DFT/B3LYP level of theory, the ECD signal response is reduced by a factor of 6.7 in going from the B-DNA to the N-DNA form, and it is illustrated how more than 90% of the individual base-pair dimers contribute to this strong hypochromic effect. Of the several inter-base pair parameters, an increase in twist angles is identified as to strongly contribute to a reduced ellipticity. The present work provides first evidence that first-principles calculations can elucidate changes in DNA dichroism due to the supramolecular organization of the nucleoprotein particle and associates these changes with the local structural features of nucleosomal DNA.

The Rayleigh light scattering properties of (H2SO4)a(NH3)b and (H2SO4)a((CH3)2NH)b atmospheric molecular clusters have been investigated using a response theory approach. Using density functional theory the molecular structures and stepwise formation free energies of clusters with a and b up to 4 have been re-investigated. The Rayleigh scattering intensities are calculated from the dipole polarizability tensor α using the CAM-B3LYP functional by applying linear response methods. The intrinsic scattering properties of (H2SO4)a(NH3)b and (H2SO4)a((CH3)2NH)b indicate that amine containing clusters scatter light significantly more efficiently then their ammonia containing counterparts. Using the Atmospheric Cluster Dynamics Code (ACDC) the steady state cluster concentrations are estimated and the effective scattering is calculated. The effective scattering is shown to be highly dependent on the estimated concentrations and indicates that there exist competitive pathways, such as nucleation and coagulation, which influence the cluster distributions. The frequency dependence of the scattering is found to depend on the cluster composition and show increased responses when clusters contain more bases than acid molecules. Based on structures obtained using semi-empirical molecular dynamics simulations the Rayleigh scattering properties of clusters with up to 20 acid-base pairs are evaluated. This study represents the first step towards gaining a fundamental understanding of the scattering properties of small atmospheric clusters in the ambient atmosphere.

The molecular recognition of human telomeric G-quadruplexes by a novel cationic π-extended Ni(II)-porphyrin (Ni(II)-TImidP4) is studied in aqueous solutions via (chir)optical spectroscopy, Fluorescence Resonance Energy Transfer (FRET) melting assay, and computational molecular modeling. The results are systematically compared with the recognition by a conventional meso-substituted Ni(II)-porphyrin (Ni(II)-TMPyP4), which allows us to pinpoint the differences in binding modes depending on the G-quadruplex topology. Importantly, FRET melting assays show the higher selectivity of Ni(II)-TImidP4 towards human telomeric G4 than that of Ni(II)-TMPyP4.

The Rayleigh and hyper Rayleigh scattering properties of the binary (H2SO4)(H2O)n and ternary (H2SO4)(NH3)(H2O)n clusters are investigated using a quantum mechanical response theory approach. The molecular Rayleigh scattering intensities are expressed using the dipole polarizability α and hyperpolarizability β tensors. Using density functional theory, we elucidate the effect of cluster morphology on the scattering properties using a combinatorial sampling approach. We find that the Rayleigh scattering intensity depends quadratically on the number of water molecules in the cluster and that a single ammonia molecule is able to induce a high anisotropy, which further increases the scattering intensity. The hyper Rayleigh scattering activities are found to be extremely low. This study presents the first attempt to map the scattering of atmospheric molecular clusters using a bottom-up approach.

The optical properties, including UV-vis spectra and resonance Raman profiles, of pristine and defected single-walled carbon nanotubes (SWCNTs) are computed using state-of-the-art time-dependent density functional theory (TDDFT) as implemented using the Liouville-Lanczos approach to linear-response TDDFT. The CNT defects were of the form of Stone-Wales and diatom-vacancies. Our results are in very good agreement with experimental results where defects were introduced into a part of defect-free CNTs. In particular, we show that the first and second π-π* excitation energies are barely shifted due to the defects and associated with a relatively small reduction in the maxima of the absorption bands. In contrast, the resonance Raman spectra show close to an order of magnitude reduction in intensities, offering a means to distinguish between pristine and defected SWCNTs even at low defect concentrations.

Theoretical analysis using quadratic response theory within the time-dependent density functional theory (TDDFT) formalism shows that the dermacozines, a group of phenazine-based compounds isolated from cultures of Dermacoccus abyssi found in the Mariana Trench, possess large first hyperpolarisability (β) values at common incident laser wavelengths that are highly sensitive to the degree and type of substitution of the core structure. The phenazine moiety is a versatile and tunable chromophore for non-linear optics and this work serves to highlight the potential that (marine) natural products, even those found in the darkest places on the planet, may have for aiding developments in optical materials design.

The development of electronic response theory in quantum chemistry has been reviewed, starting from the early 1970's and reaching the current state-of-the-art. The general theory has been applied to the calculation of a large number of spectroscopic parameters over the years, and it has been implemented for the majority of standard electronic structure methods. Two formulations of response theory, the Ehrenfest expectation value and the quasi-energy derivative formulation, have turned into leading alternatives for the derivation of computationally tractable expressions of response functions, and they are here reviewed with an attempt to, as far as possible, leave out technical details. A set of four steps are identified as common in derivations of response functions, and the two formulations are compared along this series of steps. Particular emphasis is given to the situation when the oscillation of the weak external electromagnetic field is in resonance with a transition frequency of the system. The formation of physically sound response functions in resonance regions of the spectrum is discussed in light of the causality condition and the Kramers-Kronig relations, and it is achieved in wave function theory by means of the introduction of relaxation parameters in a manner that mimics what one sees in density matrix theory. As a working example, equations are illustrated by their application to a two-state model for para-nitroaniline including the ground and the lowest charge-transfer state in the electric dipole approximation.

CHARMM force field parameter values for a class of oligothiophene derivatives have been derived with reference to density functional theory/B3LYP potential energy surfaces. The force field parametrization of these luminescent conjugated polyelectrolytes includes the electronic ground state as well as the strongly light absorbing first excited state. In conjunction with quantum chemical response theory calculations of transition state properties, a molecular dynamical model of the Stokes shift is obtained. The theoretical model is benchmarked against experimental data recorded at room temperature which refer to sodium salts of p-HTAA and p-FTAA with distilled water as a solvent. For p-HTAA the theoretically predicted Stokes shift of 112 nm is in good agreement with the experimental result of 124 nm, given the approximations about exciton localization that were introduced to obtain a force field for the excited state., (This journal is © the Owner Societies 2011)

The relativistic four-component static exchange approach for calculation of near-edge X-ray absorption spectra has been reviewed. Application of the method is made to the Au(111) interface and the adsorption of methanethiol by a study of the near sulfur L-edge spectrum. The binding energies of the sulfur 2p(1/2) and 2p(3/2) sublevels in methanethiol are determined to be split by 1.2 eV due to spin-orbit coupling, and the binding energy of the 2p(3/2) shell is lowered from 169.2 eV for the isolated system to 167.4 and 166.7-166.8 eV for methanethiol in mono- and di-coordinated adsorption sites, respectively (with reference to vacuum). In the near L-edge X-ray absorption fine structure spectrum only the sigma*(S-C) peak at 166 eV remains intact by surface adsorption, whereas transitions of predominantly Rydberg character are largely quenched in the surface spectra. The sigma*(S-H) peak of methanethiol is replaced by low-lying, isolated, sigma*(S-Au) peak(s), where the number of peaks in the latter category and their splittings are characteristic of the local bonding situation of the sulfur.