Intermolecular vibrations and fast relaxations in supercooled ionic liquids

Short-time dynamics of ionic liquids has been investigated by low-frequency Raman spectroscopy (4 < omega < 100 cm(-1)) within the supercooled liquid range. Raman spectra are reported for ionic liquids with the same anion, bis(trifluoromethylsulfonyl)imide, and different cations: 1-butyl-3-methylimidazolium, 1-hexyl-3-methylimidazolium, 1-butyl-1-methylpiperidinium, trimethylbutylammonium, and tributylmethylammonium. It is shown that low-frequency Raman spectroscopy provides similar results as optical Kerr effect (OKE) spectroscopy, which has been used to study intermolecular vibrations in ionic liquids. The comparison of ionic liquids containing aromatic and non-aromatic cations identifies the characteristic feature in Raman spectra usually assigned to librational motion of the imidazolium ring. The strength of the fast relaxations (quasi-elastic scattering, QES) and the intermolecular vibrational contribution (boson peak) of ionic liquids with non-aromatic cations are significantly lower than imidazolium ionic liquids. A correlation length assigned to the boson peak vibrations was estimated from the frequency of the maximum of the boson peak and experimental data of sound velocity. The correlation length related to the boson peak (similar to 19 angstrom) does not change with the length of the alkyl chain in imidazolium cations, in contrast to the position of the first-sharp diffraction peak observed in neutron and X-ray scattering measurements of ionic liquids. The rate of change of the QES intensity in the supercooled liquid range is compared with data of excess entropy, free volume, and mean-squared displacement recently reported for ionic liquids. The temperature dependence of the QES intensity in ionic liquids illustrates relationships between short-time dynamics and long-time structural relaxation that have been proposed for glass-forming liquids. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3604533]

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

Simulations of charge transport in organic compounds

To aid the design of organic semiconductors, we study the charge transport properties of organic liquid crystals, i.e. hexabenzocoronene and carbazole macrocycle, and single crystals, i.e. rubrene, indolocarbazole and benzothiophene derivatives (BTBT, BBBT). The aim is to find structure-property relationships linking the chemical structure as well as the morphology with the bulk charge carrier mobility of the compounds. To this end, molecular dynamics (MD) simulations are performed yielding realistic equilibrated morphologies. Partial charges and molecular orbitals are calculated based on single molecules in vacuum using quantum chemical methods. The molecular orbitals are then mapped onto the molecular positions and orientations, which allows calculation of the transfer integrals between nearest neighbors using the molecular orbital overlap method. Thus we obtain realistic transfer integral distributions and their autocorrelations. In case of organic crystals the differences between two descriptions of charge transport, namely semi-classical dynamics (SCD) in the small polaron limit and kinetic Monte Carlo (KMC) based on Marcus rates, are studied. The liquid crystals are investigated solely in the hopping limit. To simulate the charge dynamics using KMC, the centers of mass of the molecules are mapped onto lattice sites and the transfer integrals are used to compute the hopping rates. In the small polaron limit, where the electronic wave function is spread over a limited number of neighboring molecules, the Schroedinger equation is solved numerically using a semi-classical approach. The results are compared for the different compounds and methods and, where available, with experimental data. The carbazole macrocycles form columnar structures arranged on a hexagonal lattice with side chains facing inwards, so columns can closely approach each other allowing inter-columnar and thus three-dimensional transport. When taking only intra-columnar transport into account, the mobility is orders of magnitude lower than in the three-dimensional case. BTBT is a promising material for solution-processed organic field-effect transistors. We are able to show that, on the time-scales of charge transport, static disorder due to slow side chain motions is the main factor determining the mobility. The resulting broad transfer integral distributions modify the connectivity of the system but sufficiently many fast percolation paths remain for the charges. Rubrene, indolocarbazole and BBBT are examples of crystals without significant static disorder. The high mobility of rubrene is explained by two main features: first, the shifted cofacial alignment of its molecules, and second, the high center of mass vibrational frequency. In comparsion to SCD, only KMC based on Marcus rates is capable of describing neighbors with low coupling and of taking static disorder into account three-dimensionally. Thus it is the method of choice for crystalline systems dominated by static disorder. However, it is inappropriate for the case of strong coupling and underestimates the mobility of well-ordered crystals. SCD, despite its one-dimensionality, is valuable for crystals with strong coupling and little disorder. It also allows correct treatment of dynamical effects, such as intermolecular vibrations of the molecules. Rate equations are incapable of this, because simulations are performed on static snapshots. We have thus shown strengths and weaknesses of two state of the art models used to study charge transport in organic compounds, partially developed a program to compute and visualize transfer integral distributions and other charge transport properties, and found structure-mobility relations for several promising organic semiconductors.

The S-1/S-2 exciton interaction in 2-pyridone center dot 6-methyl-2-pyridone: Davydov splitting, vibronic coupling, and vibronic quenching

The excitonic splitting between the S-1 and S-2 electronic states of the doubly hydrogen-bonded dimer 2-pyridone center dot 6-methyl-2-pyridone (2PY center dot 6M2PY) is studied in a supersonic jet, applying two-color resonant two-photon ionization (2C-R2PI), UV-UV depletion, and dispersed fluorescence spectroscopies. In contrast to the C-2h symmetric (2-pyridone) 2 homodimer, in which the S-1 <- S-0 transition is symmetry-forbidden but the S-2 <- S-0 transition is allowed, the symmetry-breaking by the additional methyl group in 2PY center dot 6M2PY leads to the appearance of both the S-1 and S-2 origins, which are separated by Delta(exp) = 154 cm(-1). When combined with the separation of the S-1 <- S-0 excitations of 6M2PY and 2PY, which is delta = 102 cm(-1), one obtains an S-1/S-2 exciton coupling matrix element of V-AB, el = 57 cm(-1) in a Frenkel-Davydov exciton model. The vibronic couplings in the S-1/S-2 <- S-0 spectrum of 2PY center dot 6M2PY are treated by the Fulton-Gouterman single-mode model. We consider independent couplings to the intramolecular 6a' vibration and to the intermolecular sigma' stretch, and obtain a semi-quantitative fit to the observed spectrum. The dimensionless excitonic couplings are C(6a') = 0.15 and C(sigma') = 0.05, which places this dimer in the weak-coupling limit. However, the S-1/S-2 state exciton splittings Delta(calc) calculated by the configuration interaction singles method (CIS), time-dependent Hartree-Fock (TD-HF), and approximate second-order coupled-cluster method (CC2) are between 1100 and 1450 cm(-1), or seven to nine times larger than observed. These huge errors result from the neglect of the coupling to the optically active intra-and intermolecular vibrations of the dimer, which lead to vibronic quenching of the purely electronic excitonic splitting. For 2PY center dot 6M2PY the electronic splitting is quenched by a factor of similar to 30 (i.e., the vibronic quenching factor is Gamma(exp) = 0.035), which brings the calculated splittings into close agreement with the experimentally observed value. The 2C-R2PI and fluorescence spectra of the tautomeric species 2-hydroxypyridine center dot 6-methyl-2-pyridone (2HP center dot 6M2PY) are also observed and assigned. (C) 2011 American Institute of Physics.

Intermolecular Clamping by Hydrogen Bonds: 2-Pyridone center dot NH3

A combined spectroscopic and ab initio theoretical study of the doubly hydrogen-bonded complex of 2-pyridone (2PY) with NH3 has been performed. The S-1 <- S-0 spectrum extends up to approximate to 1200 cm(-1) above the 0(0)(0) band, close to twice the range observed for 2PY. The S-1 state nonradiative decay for vibrations above approximate to 300 cm(-1) in the NH3 complex is dramatically slowed down relative to bare 2PY. Also, the Delta v=2,4,... overtone bands of the v(1)' and v(2)' out-of-plane vibrations that dominate the low-energy spectral region of 2PY are much weaker or missing for 2PY center dot NH3, which implies that the bridging (2PY)NH center dot center dot center dot NH3 and H2NH center dot center dot center dot O=C H-bonds clamp the 2PY at a planar geometry in the S-1 state. The mass-resolved UV vibronic spectra of jet-cooled 2PY center dot NH3 and its H/D mixed isotopomers are measured using two-color resonant two-photon ionization spectroscopy. The S-0 and S-1 equilibrium structures and normal-mode frequencies are calculated by density functional (B3LYP) and correlated ab initio methods (MP2 and approximate second-order coupled-cluster, CC2). The S-1 <- S-0 vibronic assignments are based on configuration interaction singles (CIS) and CC2 calculations. A doubly H-bonded bridged structure of C-S symmetry is predicted, in agreement with that of Held and Pratt [J. Am. Chem. Soc. 1993, 115, 9718]. While the B3LYP and MP2 calculated rotational constants are in very good agreement with experiment, the calculated H2NH center dot center dot center dot O=C H-bond distance is approximate to 0.7 angstrom shorter than that derived by Held and Pratt. On the other hand, this underlines their observation that ammonia can act as a strong H-bond donor when built into an H-bonded bridge. The CC2 calculations predict the H2NH center dot center dot center dot O distance to increase by 0.2 angstrom upon S-1 <- S-0 electronic excitation, while the (2PY)NH center dot center dot center dot NH3 H-bond remains nearly unchanged. Thus, the expansion of the doubly H-bonded bridge in the excited state is asymmetric and almost wholly due to the weakening of the interaction of ammonia with the keto acceptor group.

N-H center dot center dot center dot pi hydrogen-bonding and large-amplitude tipping vibrations in jet-cooled pyrrole-benzene

The N-H center dot center dot center dot pi hydrogen bond is an important intermolecular interaction in many biological systems. We have investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet cooled complex of pyrrole with benzene and benzene-d(6) (Pyr center dot Bz, Pyr center dot Bz-d(6)). DFT-D density functional, SCS-MP2 and SCS-CC2 calculations predict a T-shaped and (almost) C(s) symmetric structure with an N-H center dot center dot center dot pi hydrogen bond to the benzene ring. The pyrrole is tipped by omega(S(0)) = +/- 13 degrees relative to the surface normal of Bz. The N center dot center dot center dot ring distance is 3.13 angstrom. In the S(1) excited state, SCS-CC2 calculations predict an increased tipping angle omega(S(1)) = +/- 21 degrees. The IR depletion spectra support the T-shaped geometry: The NH stretch is redshifted by -59 cm(-1), relative to the "free" NH stretch of pyrrole at 3531 cm(-1), indicating a moderately strong N-H center dot center dot center dot pi interaction. The interaction is weaker than in the (Pyr)(2) dimer, where the NH donor shift is -87 cm(-1) [Dauster et al., Phys. Chem. Chem. Phys., 2008, 10, 2827]. The IR C-H stretch frequencies and intensities of the Bz subunit are very similar to those of the acceptor in the (Bz)(2) dimer, confirming that Bz acts as the acceptor. While the S(1) <- S(0) electronic origin of Bz is forbidden and is not observable in the gas-phase, the UV spectrum of Pyr center dot Bz in the same region exhibits a weak 0(0)(0) band that is red-shifted by 58 cm(-1) relative to that of Bz (38 086 cm(-1)). The origin appears due to symmetry-breaking of the p-electron system of Bz by the asymmetric pyrrole NH center dot center dot center dot pi hydrogen bond. This contrasts with (Bz)(2), which does not exhibit a 0(0)(0) band. The Bz moiety in Pyr center dot Bz exhibits a 6a(0)(1) band at 0(0)(0) + 518 cm(-1) that is about 20x more intense than the origin band. The symmetry breaking by the NH center dot center dot center dot pi hydrogen bond splits the degeneracy of the v(6)(e(2g)) vibration, giving rise to 6a' and 6b' sub-bands that are spaced by similar to 6 cm(-1). Both the 0(0)(0) and 6(0)(1) bands of Pyr center dot Bz carry a progression in the low-frequency (10 cm(-1)) excited-state tipping vibration omega', in agreement with the change of the omega tipping angle predicted by SCS-MP2 and SCS-CC2 calculations.

Intermolecular and Intramolecular Charge Transfer in Functional Molecular Materials

This thesis is devoted to the investigation of inter and intramolecular charge transfer (CT) in molecular functional materials and specifically organic dyes and CT crystals. An integrated approach encompassing quantum-chemical calculations, semiempirical tools, theoretical models and spectroscopic measurements is applied to understand structure-property relationships governing the low-energy physics of these materials. Four main topics were addressed: 1) Spectral properties of organic dyes. Charge-transfer dyes are constituted by electron donor (D) and electron acceptor (A) units linked through bridge(s) to form molecules with different symmetry and dimensionality. Their low-energy physics is governed by the charge resonance between D and A groups and is effectively described by a family of parametric Hamiltonians known as essential-state models. These models account for few electronic states, corresponding to the main resonance structures of the relevant dye, leading to a simple picture that is completed introducing the coupling of the electronic system to molecular vibrations, treated in a non-adiabatic way, and an effective classical coordinate, describing polar solvation. In this work a specific essential-state model was proposed and parametrized for the dye Brilliant Green. The central issue in this work has been the definition of the diabatic states, a not trivial task for a multi-branched chromophore. In a second effort, we have used essential-state models for the description of the early-stage dynamics of excited states after ultrafast excitation. Crucial to this work is the fully non-adiabatic treatment of the coupled electronic and vibrational motion, allowing for a reliable description of the dynamics of systems showing a multistable, broken-symmetry excited state. 2) Mixed-stack CT salts. Mixed-stack (MS) CT crystals are an interesting class of multifunctional molecular materials, where D and A molecules arrange themselves to form stacks, leading to delocalized electrons in one dimension. The interplay between the intermolecular CT, electrostatic interactions, lattice phonons and molecular vibrations leads to intriguing physical properties that include (photoinduced) phase transitions, multistability, antiferromagnetism, ferroelectricity and potential multiferroicity. The standard microscopic model to describe this family of materials is the Modified Hubbard model accounting for electron-phonon coupling (Peierls coupling), electron-molecular vibrations coupling (Holstein coupling) and electrostatic interactions. We adopt and validate a method, based on DFT calculations on dimeric DA structures, to extract relevant model parameters. The approach offers a powerful tool to shed light on the complex physics of MS-CT salts. 3) Charge transfer in organic radical dipolar dyes. In collaboration with the group of Prof. Jaume Veciana (ICMAB- Barcellona), we have studied spectral properties of a special class of CT dyes with D-bridge-A structure where the acceptor group is a stable radical (of the perchlorotriphenylmethyl, PTM, family), leading to an open-shell CT dyes. These materials are of interest since they associate the electronic and optical properties of CT dyes with magnetic properties from the unpaired electron. The first effort was devoted to the parametrization of the relevant essential-state model. Two strategies were adopted, one based on the calculation of the low-energy spectral properties, the other based on the variation of ground state properties with an applied electric field. 4) The spectral properties of organic nanoparticles based on radical species are investigated in collaboration with Dr. I. Ratera (ICMAB- Barcellona). Intriguing spectroscopic behavior was observed pointing to the presence of excimer states. In an attempt to rationalize these findings, extensive calculations (TD-DFT and ZINDO) were performed. The results for the isolated dyes are validated against experimental spectra in solution. To address intermolecular interactions we studied dimeric structures in the gas phase, but the preliminary results obtained do not support excimer formation.

EVALUATING VIBRATIONS ON A SMALL-SCALE MODEL OF A TIMBER FOOTBRIDGE

In Brazil, the study of pedestrian-induced vibration on footbridges has been undertaken since the early 1990s, for concrete and steel footbridges. However, there are no recorded studies of this kind for timber footbridges. Brazilian code ABNT NBR 7190 (1997) gives design requirements only for static loads in the case of timber footbridges, without considering the serviceability limit state from pedestrian-induced vibrations. The aim of this work is to perform a theoretical dynamic, numerical and experimental analysis on simply-supported timber footbridges, by using a small-scale model developed from a 24 m span and 2 m width timber footbridge, with two main timber beams. Span and width were scaled down (1:4) to 6 m e 0.5 in, respectively. Among the conclusions reached herein, it is emphasized that the Euler-Bernoulli beam theory is suitable for calculating the vertical and lateral first natural frequencies in simply-supported timber footbridges; however, special attention should be given to the evaluation of lateral bending stiffness, as it leads to conservative values.

Intermolecular energy transfer and photostability of luminescence-tuneable multicolour PMMA films doped with lanthanide-beta-diketonate complexes

It is reported in this work the preparation, characterisation and photoluminescence study of poly(methylmethacrylate) (PMMA) thin films co-doped with [Eu(tta)(3)(H(2)O)(2)] and [Tb(acac)(3)(H(2)O)(3)] complexes. Both the composition and excitation wavelength may be tailored to fine-tune the emission properties of these Ln(3+)-beta-diketonate doped polymer films, exhibiting green and red primary colours, as well as intermediate colours. In addition to the ligand-Ln(3+) intramolecular energy transfer, it is observed an unprecedented intermolecular energy transfer process from the (5)D(4) emitting level of the Tb(3+) ion to the excited triplet state T(1) of the tta ligand coordinated to the Eu(3+) ion. The PMMA polymer matrix acts as a co-sensitizer and enhances the overall luminescence intensity of the polymer films. Furthermore, it provides considerable UV protection for the luminescent species and improves the photostability of the doped system.

Modeling and Analysis of Piezoelectric Energy Harvesting From Aeroelastic Vibrations Using the Doublet-Lattice Method

Multifunctional structures are pointed out as an important technology for the design of aircraft with volume, mass, and energy source limitations such as unmanned air vehicles (UAVs) and micro air vehicles (MAVs). In addition to its primary function of bearing aerodynamic loads, the wing/spar structure of an UAV or a MAV with embedded piezoceramics can provide an extra electrical energy source based on the concept of vibration energy harvesting to power small and wireless electronic components. Aeroelastic vibrations of a lifting surface can be converted into electricity using piezoelectric transduction. In this paper, frequency-domain piezoaeroelastic modeling and analysis of a canti-levered platelike wing with embedded piezoceramics is presented for energy harvesting. The electromechanical finite-element plate model is based on the thin-plate (Kirchhoff) assumptions while the unsteady aerodynamic model uses the doublet-lattice method. The electromechanical and aerodynamic models are combined to obtain the piezoaeroelastic equations, which are solved using a p-k scheme that accounts for the electromechanical coupling. The evolution of the aerodynamic damping and the frequency of each mode are obtained with changing airflow speed for a given electrical circuit. Expressions for piezoaeroelastically coupled frequency response functions (voltage, current, and electrical power as well the vibratory motion) are also defined by combining flow excitation with harmonic base excitation. Hence, piezoaeroelastic evolution can be investigated in frequency domain for different airflow speeds and electrical boundary conditions. [DOI:10.1115/1.4002785]

A Short Note on Transverse Vibrations of a Shaft Carrying Two (or One) Disk Excited by a Nonideal Motor

The paper presents a number of numerical simulations of the transverse vibrations of two (or one) imbalanced rotors forced by an electric motor with limited power supply, during the passage through of the two resonance zones (increasing and decreasing input voltages). The predominant presence of the Sommerfeld effect. when the rotational velocity of the motor is captured, in the second resonance frequency is demonstrated. We have shown that the hysteretic jump phenomenon exists in a rotor system with two (or one) disks, and with this, we have shown that a torque is influenced by the dynamical behavior of die rotor [DOI: 10.1115/1.3007979]

Thoracic vibrations in stingless bees (Melipona seminigra): resonances of the thorax influence vibrations associated with flight but not those associated with sound production

Bees generate thoracic vibrations with their indirect flight muscles in various behavioural contexts. The main frequency component of non-flight vibrations, during which the wings are usually folded over the abdomen, is higher than that of thoracic vibrations that drive the wing movements for flight. So far, this has been concluded from an increase in natural frequency of the oscillating system in association with the wing adduction. In the present study, we measured the thoracic oscillations in stingless bees during stationary flight and during two types of non-flight behaviour, annoyance buzzing and forager communication, using laser vibrometry. As expected, the flight vibrations met all tested assumptions for resonant oscillations: slow build-up and decay of amplitude; increased frequency following reduction of the inertial load; and decreased frequency following an increase of the mass of the oscillating system. Resonances, however, do not play a significant role in the generation of non-flight vibrations. The strong decrease in main frequency at the end of the pulses indicates that these were driven at a frequency higher than the natural frequency of the system. Despite significant differences regarding the main frequency components and their oscillation amplitudes, the mechanism of generation is apparently similar in annoyance buzzing and forager vibrations. Both types of non-flight vibration induced oscillations of the wings and the legs in a similar way. Since these body parts transform thoracic oscillations into airborne sounds and substrate vibrations, annoyance buzzing can also be used to study mechanisms of signal generation and transmission potentially relevant in forager communication under controlled conditions.