Identification of adsorbed molecules via STM tip manipulation: CO, H₂O, and O₂ on TiO₂ anatase (101)

While Scanning Tunneling Microscopy (STM) has evolved as an ideal tool to study surface chemistry at the atomic scale, the identification of adsorbed species is often not straightforward. This paper describes a way to reliably identify H2O, CO and O2 on the TiO2 anatase (101) surface with STM. These molecules are of a key importance in the surface chemistry of this and many other (photo-) catalytic materials. They exhibit a wide variety of contrasts in STM images, depending on the tip condition. With clean, metallic tips the molecules appear very similar, i.e., as bright, dimer-like features located in the proximity of surface Ti5c atoms. However, each species exhibits a specific response to the electric field applied by the STM tip. It is shown that this tip–adsorbate interaction can be used to reliably ascertain the identity of such species. The tip–adsorbate interactions, together with comparison of experimental and calculated STM images, are used to analyse and revisit the assignments of molecular adsorbed species reported in recent studies.

Distributed Atomic Polarizabilities of Amino Acids and their Hydrogen-Bonded Aggregates

With the purpose of rational design of optical materials, distributed atomic polarizabilities of amino acid molecules and their hydrogen-bonded aggregates are calculated in order to identify the most efficient functional groups, able to buildup larger electric susceptibilities in crystals. Moreover, we carefully analyze how the atomic polarizabilities depend on the one-electron basis set or the many-electron Hamiltonian, including both wave function and density functional theory methods. This is useful for selecting the level of theory that best combines high accuracy and low computational costs, very important in particular when using the cluster method to estimate susceptibilities of molecular-based materials.

Interactions of Polyvinylpyrrolidone with Chlorin e6-Based Photosensitizers Studied by NMR and Electronic Absorption Spectroscopy

Polyvinylpyrrolidone (PVP) can act as potential drug delivery vehicle for porphyrin-based photosensitizers in photodynamic therapy (PDT) to enhance their stability and prevent porphyrin self-association. In the present study the interactions of PVP (MW 10 kDa) were probed with five different derivatives of chlorin e6 (CE6) bearing either one of the amino acids serine, lysine, tyrosine or arginine, or monoamino-hexanoic acid as substituent. All derivatives of CE6 (xCE) formed aggregates of a similar structure in aqueous buffer in the millimolar range. In the presence of PVP monomerization of all xCE aggregates could be proved by 1H NMR spectroscopy. xCE-PVP complex formation was confirmed by 1H NMR T2 relaxation and diffusion ordered spectroscopy (DOSY). 1H1H-NOESY data suggested that the xCE uptake into the PVP polymer matrix is governed by hydrophobic interactions. UV–vis absorption and fluorescence emission bands of xCE in the micromolar range revealed characteristic PVP-induced bathochromic shifts. The presented data point out the potential of PVP as carrier system for amphiphilic derivatives of chlorin e6. The capacity of PVP to monomerize xCE aggregates may enhance their efficiency as possible photosensitizers in PDT.

Simulating Energy Transfer and Upconversion in β-NaYF 4 : Yb 3+ , Tm 3+

The design of upconversion phosphors with higher quantum yield requires a deeper understanding of the detailed energy transfer and upconversion processes between active ions inside the material. Rate equations can model those processes by describing the populations of the energy levels of the ions as a function of time. However, this model presents some drawbacks: energy migration is assumed to be infinitely fast, it does not determine the detailed interaction mechanism (multipolar or exchange), and it only provides the macroscopic averaged parameters of interaction. Hence, a rate equation model with the same parameters cannot correctly predict the time evolution of upconverted emission and power dependence under a wide range of concentrations of active ions. We present a model that combines information about the host material lattice, the concentration of active ions, and a microscopic rate equation system. The extent of energy migration is correctly taken into account because the energy transfer processes are described on the level of the individual ions. This model predicts the decay curves, concentration, and excitation power dependences of the emission. This detailed information can be used to predict the optimal concentration that results in the maximum upconverted emission.

Probing the Structure, Pseudorotation, and Radial Vibrations of Cyclopentane by Femtosecond Rotational Raman Coherence Spectroscopy

Femtosecond time-resolved Raman rotational coherence spectroscopy (RCS) is employed to determine accurate rotational, vibration–rotation coupling constants, and centrifugal distortion constants of cyclopentane (C⁵H¹⁰). Its lowest-frequency vibration is a pseudorotating ring deformation that interconverts 10 permutationally distinct but energetically degenerate “twist” minima interspersed by 10 “bent” conformers. While the individual twist and bent structures are polar asymmetric tops, the pseudorotation is fast on the time scale of external rotation, rendering cyclopentane a fluxionally nonpolar symmetric top molecule. The pseudorotational level pattern corresponds to a one-dimensional internal rotor with a pseudorotation constant Bps ≈ 2.8 cm⁻¹. The pseudorotational levels are significantly populated up to l = ± 13 at 298 K; <10% of the molecules are in the l = 0 level. The next-higher vibration is the “radial” ν²³ ring deformation mode at 273 cm⁻¹, which is far above the pseudorotational fundamental. Femtosecond Raman RCS measurements were performed in a gas cell at T = 293 K and in a pulsed supersonic jet at T ≈ 90 K. The jet cooling reduces the pseudorotational distribution to l < ±8 and eliminates the population of ν²³, allowing one to determine the rotational constant as A0 = B0 = 6484.930(11) MHz. This value is ∼300 times more precise than the previous value. The fit of the RCS transients reveals that the rotation–pseudorotation coupling constant αe,psB = −0.00070(1) MHz is diminutive, implying that excitation of the pseudorotation has virtually no effect on the B0 rotational constant of cyclopentane. The smallness of αe,psB can be realized when comparing to the vibration–rotation coupling constant of the ν²³ vibration, αe,23B = −9.547(1) MHz, which is about 10⁴ times larger.

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

Modeling the Histidine–Phenylalanine Interaction: The NH···π Hydrogen Bond of Imidazole·Benzene

NH···π hydrogen bonds occur frequently between the amino acid side groups in proteins and peptides. Data-mining studies of protein crystals find that ~80% of the T-shaped histidine···aromatic contacts are CH···π, and only ~20% are NH···π interactions. We investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet-cooled imidazole·benzene (Im·Bz) complex as a model for the NH···π interaction between histidine and phenylalanine. Ground- and excited-state dispersion-corrected density functional calculations and correlated methods (SCS-MP2 and SCS-CC2) predict that Im·Bz has a Cs-symmetric T-shaped minimum-energy structure with an NH···π hydrogen bond to the Bz ring; the NH bond is tilted 12° away from the Bz C₆ axis. IR depletion spectra support the T-shaped geometry: The NH stretch vibrational fundamental is red shifted by −73 cm⁻¹ relative to that of bare imidazole at 3518 cm⁻¹, indicating a moderately strong NH···π interaction. While the Sₒ(A1g) → S₁(B₂u) origin of benzene at 38 086 cm⁻¹ is forbidden in the gas phase, Im·Bz exhibits a moderately intense Sₒ → S₁ origin, which appears via the D₆h → Cs symmetry lowering of Bz by its interaction with imidazole. The NH···π ground-state hydrogen bond is strong, De=22.7 kJ/mol (1899 cm⁻¹). The combination of gas-phase UV and IR spectra confirms the theoretical predictions that the optimum Im·Bz geometry is T shaped and NH···π hydrogen bonded. We find no experimental evidence for a CH···π hydrogen-bonded ground-state isomer of Im·Bz. The optimum NH···π geometry of the Im·Bz complex is very different from the majority of the histidine·aromatic contact geometries found in protein database analyses, implying that the CH···π contacts observed in these searches do not arise from favorable binding interactions but merely from protein side-chain folding and crystal-packing constraints. The UV and IR spectra of the imidazole·(benzene)₂ cluster are observed via fragmentation into the Im·Bz+ mass channel. The spectra of Im·Bz and Im·Bz₂ are cleanly separable by IR hole burning. The UV spectrum of Im·Bz₂ exhibits two 000 bands corresponding to the Sₒ → S₁ excitations of the two inequivalent benzenes, which are symmetrically shifted by −86/+88 cm⁻¹ relative to the 000 band of benzene.

Carnot Cycles: Traditional and Stefan Boltzmann, and the First Planck Argument

The Stefan Boltzmann equation is obtained using a non-traditional Carnot Engine. In addition, the original Planck argument for radiation density is given.

MO visualization of Pi Orbitals

In the spirit of trying to convert people to understanding atomic orbitals centered elsewhere than the origin, we continue the discussion of visualizing molecular orbitals, so called LCAO-MO, using various plotting tricks in Maple.

Lee O Case's "Elements of the Phase Rule", Chapter 2

The second chapter of the Case text on the phase rule, this time dealing with two component systems, added to by the author, is presented here.

Guessing Harmonic Oscillator Wavefunctions Using Maple

The guessing of eigenfunctions is not trivial at higher quantum numbers, no matter what the system being considered. Instead of guessing, one can employ a symbolic calculus progam (Maple in this case) to aid in the reasoning process.

Parameterizing a Simple Potential Energy Function Appropriate for Vibrational Frequencies Using Molecular Dynamics

A Maple scheme for quickly parameterizing vibrational potential energy functions is presented. As an example, the potential energy function's parameters for the vibrational motions in H_2O_2 are obtained assuming the simplest potential energy function. This paper was originally written as a research paper, but rejected by the referees. It is therefore being edited into an ``educational'' paper for student usage.