Low-lying electronic excitations and optical absorption spectra of the black dye sensitizer: a first-principles study

We report on ab-initio calculations of the electronic structure and optical absorption response of the black dye sensitizer in gas phase. We show that, despite the large size of this molecule, the second-order multiconfiguration quasi-degenerate perturbation theory (MC-QDPT) can be used to calculate vertical excitation energies, oscillator strengths and optical absorption spectra. The zeroth-order reference states entering perturbation calculations are complete active space (CAS) configuration interaction (CI) wave functions computed for 12 active electrons distributed in 12 active orbitals. We found that the CI approach is not enough for taking into account the strong dynamical correlation effects in this system. In fact, the excitation energies of the CAS-CI target states are strongly renormalized by the MC-QDPT calculations. In the calculated absorption spectra, the analysis of the perturbed wavefunctions revealed that the stronger absorption bands correspond to metal-to-ligand and ligand-to-ligand charge transfer processes. Comparison with independent time-dependent extension (TDDFT) calculations performed with different functionals shows that corrections to the long-range behavior of the functional is pivotal to achieve agreement with the MC-QDPT results.     

The dynamics of water evaporation from partially solvated cytochrome c in the gas phase

The study of evaporation of water from biological macromolecules is important for the understanding of electrospray mass spectrometry experiments. In electrospray ionization (ESI), electrically charged nanoscale droplets are formed from solutions of, for example, proteins. Then evaporation of the solvent leads to dry protein ions that can be analyzed in the mass spectrometer. In this work the dynamics of water evaporation from native cytochrome c covered by a monolayer of water is studied by molecular dynamics (MD) simulations at constant energy. A model of the initial conditions of the process is introduced. The temperature of the protein drops by about 100 K during the 400 picoseconds of the simulations. This sharp drop in temperature causes the water evaporation rate to decrease by about an order of magnitude, leaving the protein with 50% to 90% of the original water molecules, depending on the initial temperature of the simulation. The structural changes of the protein upon desolvation were considered through calculations of the radius of gyration and the root mean square (RMS) of the protein. A variation of 0.4 A in the radius of gyration, together with an RMS value of less than 3 A, indicates only minor changes in the overall shape of the protein structure. The water coordination number of the solvation shell is much smaller than that for bulk water. The mobility of water is high at the beginning of the simulations and drops as the simulation progresses and the temperature decreases. Incomplete desolvation of protein ions was also observed in recent experiments.     

Rotationally resolved electronic spectra of 9,10-dihydrophenanthrene. A "floppy" molecule in the gas phase

Rotationally resolved fluorescence excitation spectra of several bands in the S1<--S0 electronic spectrum of 9,10-dihydrophenanthrene (DHPH) have been observed and assigned. Each band was fit using rigid rotor Hamiltonians in both electronic states. Analyses of these data reveal that DHPH has a nonplanar configuration in its S0 state with a dihedral angle between the aromatic rings (phi) of approximately 21.5 degrees. The data also show that excitation of DHPH with UV light results in a more planar structure of the molecule in the electronically excited state, with phi approximately 8.5 degrees. Three prominent Franck-Condon progressions appear in the low resolution spectrum, all with fundamental frequencies lying below 300 cm(-1). Estimates of the potential energy surfaces along each of these coordinates have been obtained from analyses of the high resolution spectra. The remaining barrier to planarity in the S1 state is estimated to be approximately 2650 cm(-1) along the bridge deformation mode and is substantially reduced by excitation of the molecule along the (orthogonal) ring twisting coordinate.     

Low-lying electronic states of the Ti2 dimer: electronic absorption spectroscopy in rare gas matrices in concert with quantum chemical calculations

Absorption spectra were measured for Ti2 in Ne and Ar matrices. The spectra give evidence for several electronic transitions in the region between 4000 and 10 000 cm(-1) and provide important information about some excited electronic states of Ti2 in proximity to the ground state. The vibrational fine structure measured for these transitions allowed to calculate the force constants and the anharmonicity of the potential energy curves of the excited states, and to estimate changes in the internuclear Ti-Ti distances relative to the electronic ground state. The quantum chemical studies confirm the previously suggested (3)Delta(g) state as the ground state of Ti2. The equilibrium bond distance is calculated to be 195.4 pm. The calculated harmonic frequency of 432 cm(-1) is in good agreement with the experimental value of 407.0 cm(-1). With the aid of the calculations it was possible to assign the experimentally observed transitions in the region between 4000 and 10 000 cm(-1) to the 1 (3)Pi(u)<--(3)Delta(g), 1 (3)Phi(u)<--(3)Delta(g), 2 (3)Pi(u)<--(3)Delta(g), 2 (3)Phi(u)<--(3)Delta(g), and (3)Delta(u)<--(3)Delta(g) excitations (in the order of increasing energy). The calculated relative energies and harmonic frequencies are in pleasing agreement with the experimentally obtained values, with deviations of less than 5% and 2%, respectively. The bond distances estimated on the basis of the experimental spectra tally satisfactorily with the predictions of our calculations.     

Comparative characteristics of a tryptophan molecule in the gas phase and water

The L and D isomers of the tryptophan (Trp) molecule and the (Trp)⁺ cation in the gas phase and water are calculated at the DFT level to reveal the effect of water considered in the dielectric continuum approximation on the electronic characteristics of the molecule. The distribution of effective atomic charges and bond lengths enables the prediction of the most probable parts of the chemical bond cleavage during the fragmentation of the molecule under the ionizing particle flux. These data are supplemented with a calculation of fragmentation energies. Zwitterionic structures characterized by the appearance of considerable dipole moments and a change in their orientation with respect to the ground state are distinguished among the possible isomeric forms in water.     

Resonance energy transfer from a fluorescent dye molecule to plasmon and electron-hole excitations of a metal nanoparticle

The authors study the distance dependence of the rate of electronic excitation energy transfer from a dye molecule to a metal nanoparticle. Using the spherical jellium model, they evaluate the rates corresponding to the excitation of l=1, 2, and 3 modes of the nanoparticle. The calculation takes into account both the electron-hole pair and the plasmon excitations of the nanoparticle. The rate follows conventional R(-6) dependence at large distances while small deviations from this behavior are observed at shorter distances. Within the framework of the jellium model, it is not possible to attribute the experimentally observed d(-4) dependence of the rate to energy transfer to plasmons or electron-hole pair excitations.     

tert-Butylphosphonic acid: from the bulk to the gas phase

The structure of tert-butylphosphonic acid in the solid, in solution, and in the gas phase was studied by single-crystal X-ray diffraction, (1)H and (31)P NMR spectroscopic studies in solution, solid-state (31)P NMR spectroscopy, and electrospray ionization mass spectrometry. In addition, density functional theory (DFT) calculations at the B3LYP/6-31G*, B3LYP/6-31+G*, and B3LYP/6-311+G* level of theory for a large number of H-bonded aggregates of the type (tBuPO(3)H(2))(n) (C(n), P(n); n=1-7) support the experimental work. Crystallization of tBuPO(3)H(2) from polar solvents such as CH(3)CN or THF gives the H-bonded one-dimensional polymer 2, whereas crystallization from the less polar solvent CDCl(3) favors the formation of the H-bonded cluster (tBuPO(3)H(2))(6).CDCl(3) (1). In CDCl(3) the hexamer (tBuPO(3)H(2))(6) (C(6)) is replaced by smaller aggregates down to the monomer with decreasing concentration. DFT calculations and natural bond orbital (NBO) analyses for the clusters C(1)-C(7) and the linear arrays P(1)-P(7) reveal the hexamer C(6) to be the energetically favored structure resulting from cooperative strengthening of the hydrogen bonds in the H-bonded framework. However, the average hydrogen bond strengths calculated for C(6) and P(2) do not differ significantly (42-43 kJ mol(-1)). The average distances r(O.O), r(Obond;H), r(Pdbond;O), and r(Pbond;OH) in C(1)-C(7) and P(1)-P(7) are closely related to the hydrogen bond strength. Electrospray ionization mass spectrometry shows the presence of different anionic species of the type [(tBuPO(3)H(2))(n)-H](-) (A(1)-A(7), n=1-7) depending on the instrumental conditions. DFT calculations at the B3LYP/6-31G* level of theory were carried out for A(1)-A(6). We suggest the dimer [(tBuPO(3)H(2))(2)-H](-) (A(2)) and the trimer [(tBuPO(3)H(2))(3)-H](-) (A(3)) are the energetically favored anionic structures. A hydrogen bond energy of approximately 83 kJ mol(-1) was calculated for A(2). Electrospray ionization mass spectrometry is not suitable to study the assembling process of neutral H-bonded tert-butylphosphonic acid since the removal of a proton from the neutral aggregates has a large influence on the hydrogen bond strength and the cluster structure.     

Geometric structure of molecule of N-nitropyrrolidine in the gas phase

M. V. Lomonosov Moscow State University. Translated from Zhurnal Strukturnoi Khimii, Vol. 33, No. 1, pp. 46–50, January–February, 1992.     

Cluster expansion of electronic excitations: application to fcc Ni-Al alloys

The cluster expansion method is applied to electronic excitations and a set of effective cluster densities of states (ECDOS) is defined, analogous to effective cluster interactions (ECIs). The ECDOSs are used to generate alloy thermodynamic properties as well as the equation of state (EOS) of electronic excitations for the fcc Ni-Al systems. When parent clusters have a small size, the convergence of the expansion is not so good but the electronic density of state (DOS) is well reproduced. However, the integrals of the DOS such as the cluster expanded free energy, entropy, and internal energy associated with electronic excitations are well described at the level of the tetrahedron-octahedron cluster approximation, indicating that the ECDOS is applicable to produce electronic ECIs for cluster variation method (CVM) or Monte Carlo calculations. On the other hand, the Gruneisen parameter, calculated with first-principles methods, is no longer a constant and implies that the whole DOS profile should be considered for EOS of electronic excitations, where ECDOS adapts very well for disordered alloys and solid solutions.     

Generation and orientation of organoxenon molecule H-Xe-CCH in the gas phase

We report on the first observation of the organoxenon HXeCCH molecule in the gas phase. This molecule has been prepared in a molecular beam experiment by 193 nm photolysis of an acetylene molecule on Xe(n) clusters (n approximately 390). Subsequently the molecule has been oriented via the pseudo-first-order Stark effect in a strong electric field of the polarized laser light combined with the weak electrostatic field in the extraction region of a time-of-flight spectrometer. The experimental evidence for the oriented molecule has been provided by measurements of its photodissociation. For comparison, photolysis of C(2)H(2) on Ar(n) clusters (n approximately 280) has been measured. Here the analogous rare gas molecule HArCCH could not be generated. The interpretation of our experimental findings has been supported by ab initio calculations. In addition, the experiment together with the calculations reveals information on the photochemistry of the HXeCCH molecule. The 193 nm radiation excites the molecule predominantly into the 2 (1)Sigma(+) state, which cannot dissociate the Xe-H bond directly, but the system evolves along the Xe-C coordinate to a conical intersection of a slightly nonlinear configuration with the dissociative 1 (1)Pi state, which then dissociates the Xe-H bond.     

Phase control of the competition between electronic transitions in a solvated laser dye

Excited state population can be manipulated by resonant chirped laser pulses through pump–dump processes. We investigate these processes in the laser dye LD690 as a function of wavelength by monitoring the saturated absorption of chirped ultrafast pulses. The resulting nonlinear absorption spectrum becomes increasingly complex as the pulse is tuned to shorter wavelengths. However, fluorescence measurements indicate that the excited state population depends weakly on chirp when the pump wavelength is far from the lowest order electronic transition. Using a learning algorithm and closed-loop control, we find nonlinear chirp parameters that optimize features in the transmission spectrum. The results are discussed in terms of competition between excited state absorption and stimulated resonant Raman scattering.     

Femtosecond electronic relaxation of excited metalloporphyrins in the gas phase

A systematic study of the ultrafast decay of metalloporphyrins containing various transition metals with partially filled 3d shells and zinc (3d filled) is reported here after excitation in the second excited state of the system (Soret band). Both time-of-flight mass spectrometry and velocity map imaging have been used for detection. A general biexponential decay with a short time constant tau1 approximately 100 fs is observed for the transition metal porphyrins, followed by a tau2 approximately 1 ps time decay. This evolution is interpreted as a porphyrin-to-metal charge transfer, tau1, followed by a back transfer, tau2, which leads to an excited state (d,d*) localized on the metal. These conclusions stem from the different behaviors of zinc and the transition metal porphyrins. A porphyrin-to-metal charge transfer model is chosen to describe the relaxation mechanism, based upon the fact that transition metalloporphyrins can accept electrons on the metal site, in contrast to zinc porphyrins.     

Rotationally resolved electronic spectroscopy of tryptophol in the gas phase

High resolution S1-S0 fluorescence excitation spectra of tryptophol have been observed in the collision-free environment of a supersonic beam. Each origin band has been assigned to a unique conformer based on its observed rotational constants. Unlike its close relative tryptamine, which exhibits seven distinguishable conformers under similar conditions, tryptophol exhibits only four (GPy-in, GPh-in, and two anti structures). Possible reasons for this difference in behavior are discussed.     

Low-lying electronic states of HNCS and its ions: a CASSCF/CASPT2 study

Complete active space self-consistent-field (CASSCF) and multiconfigurational second-order perturbation theory (CASPT2) calculations in conjunction with the ANO-L basis set were performed to investigate systematically the low-lying electronic states of HNCS and its ions in C _s symmetry. Our highly accurate calculation indicated that theoretically determined geometric parameters and harmonic vibrational frequencies for the ground-state X ¹A′ are in good agreement with observed experimental data. The geometry of triplet HNCS is clearly favored C ₁ symmetry, and the relative energy is predicted to be 3.000 eV (69.2 kcal/mol). The vertical transition energies for the selected excited states of HNCS were calculated at CASSCF/CASPT2/ANO-L level of theory based on CASSCF optimized geometry. Except for a few linear states of X ²Π (1²A′, 1²A″), 1⁴Σ⁻ (1⁴A″), and 1²Σ⁺ (3²A′) states of HNCS⁺, our results confirmed that the majority of excited states are twisted trans-bend structures. The existence of bound excited anion states has been found for the first time in HNCS⁻. A more elaborate examination of ionization potential of HNCS (AIP, VIP) than previous reports has been presented.     

Phosphorous disulfide and its ions in the electronic ground state: a gas phase theoretical study

Ab initio molecular orbital theory and density functional theory calculations were performed on the electronic ground states of the open-shell PS₂ molecule and its singly charged ions. A comparison of the optimized molecular structures indicates as the stepwise one-electron reduction of the PS₂⁺ ion, to yield PS₂ and PS₂⁻, provokes a symmetric elongation of both PS bonds along with a bending of its linear equilibrium geometry. The ionization potential (IP), adiabatic electron affinity (EA_ad), and atomization energy (AE) of the open-shell PS₂ molecule were calculated at different levels of theory. The following values were obtained at the more realistic UMP4SDTQ/6-311+G(3df)//UHF/6-311+G(3df) level of theory: IP=8.32 eV, EA_ad=3.03 eV and AE=12.40 eV. At the same level of theory, the calculated vertical detachment energy (VDE) of the PS₂⁻ anion is 3.22 eV. The donor–acceptor complexes formed in the gas-phase upon interaction of either one or two ammonia molecules with PS₂⁺ were also investigated. The calculated gas-phase binding energies indicate that the formation of the bis-adduct is favored over that of the mono-adduct by a binding energy gain of about 20 kcal/mol.