Treatment Of Anharmonicity Research Articles

Phonon collapse and anharmonic melting of the 3D charge-density wave in kagome metals

The charge-density wave (CDW) mechanism and resulting structure of the AV3Sb5 family of kagome metals has posed a puzzling challenge since their discovery four years ago. In fact, the lack of consensus on the origin and structure of the CDW hinders the understanding of the emerging phenomena. Here, by employing a non-perturbative treatment of anharmonicity from first-principles calculations, we reveal that the charge-density transition in CsV3Sb5 is driven by the large electron-phonon coupling of the material and that the melting of the CDW state is attributed to ionic entropy and lattice anharmonicity. The calculated transition temperature is in very good agreement with experiments, implying that soft mode physics are at the core of the charge-density wave transition. Contrary to the standard assumption associated with a pure kagome lattice, the CDW is essentially three-dimensional as it is triggered by an unstable phonon at the L point. The absence of involvement of phonons at the M point enables us to constrain the resulting symmetries to six possible space groups. The unusually large electron-phonon linewidth of the soft mode explains why inelastic scattering experiments did not observe any softened phonon. We foresee that large anharmonic effects are ubiquitous and could be fundamental to understand the observed phenomena also in other kagome families.

Twin peaks: Matrix isolation studies of H2S·amine complexes shedding light on fundamental S-H⋯N bonding.

Noncovalent bonding between atmospheric molecules is central to the formation of aerosol particles and cloud condensation nuclei and, consequently, radiative forcing. While our understanding of O-H⋯B interactions is well developed, S-H⋯B hydrogen bonding has received far less attention. Sulfur- and nitrogen-containing molecules, particularly amines, play a significant role in atmospheric chemistry, yet S-H⋯N interactions are not well understood at a fundamental level. To help characterize these systems, H2S and methyl-, ethyl-, n-propyl-, dimethyl-, and trimethylamine (MA, EA, n-PA, DMA, and TMA) have been investigated using matrix isolation Fourier transform infrared spectroscopy and high-level theoretical methods. Experiments showed that H2S forms hydrogen bonded complexes with each of the amines, with bond strengths following the trend MA ≈ EA ≈ n-PA < TMA ≤ DMA, in line with past experimental work on H2SO4·amine complexes. However, the calculated results indicated that the trend should be MA < DMA < TMA, in line with past theoretical work on H2SO4·amine complexes. Evidence of strong Fermi resonances indicated that anharmonicity may play a critical role in the stabilization of each complex. The theoretical results were able to replicate experiment only after binding energies were recalculated to include the anharmonic effects. In the case of H2SO4·amine complexes, our results suggest that the discrepancy between theory and experiment could be reconciled, given an appropriate treatment of anharmonicity.

The infrared absorption spectrum of phenylacetylene and its deuterated isotopologue in the mid- to far-IR.

Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, IR-UV ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectra of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm-1 range originate from CH(D) in-plane and out-of-plane wags and bends, as well as bending motions including the C≡C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm-1 for d1-phenylacetylene, a one-to-one match is seen when comparing computations and experiments with all features assigned to combination bands and overtones. The C≡C stretch observed in phenylacetylene is not observed in d1-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these results, it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general.

Thermoelectric properties of the bismuth oxychalcogenides Bi2SO2, Bi2SeO2 and Bi2TeO2

We present a detailed theoretical study of the thermoelectric properties of the bismuth oxychalcogenides Bi2ChO2 (Ch = S, Se, Te). The electrical transport is modelled using semi-classical Boltzmann transport theory with electronic structures from hybrid density-functional theory, including an approximate model for the electron lifetimes. The lattice thermal conductivity is calculated using first-principles phonon calculations with an explicit treatment of anharmonicity, yielding microscopic insight into how partial replacement of the chalcogen in the bismuth chalcogenides impacts the phonon transport. We find very good agreement between the predicted transport properties and a favourable cancellation of errors that allows for near-quantitative predictions of the thermoelectric figure of merit ZT. Our calculations suggest recent experiments on n-doped Bi2SeO2 have achieved close to the largest ZT possible in bulk materials, whereas the largest reported ZT for Bi2TeO2 could be improved sixfold by optimising the carrier concentration. We also predict that much larger ZT > 2.5, competitive with the benchmark thermoelectric SnSe, could be obtained for Bi2SO2 and Bi2SeO2 with heavy p-type doping. This study demonstrates the predictive power of this modelling approach for studying thermoelectrics and highlights several avenues for improving the performance of the Bi2ChO2.

Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules

Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.

Theoretical Kinetic Study on Hydrogen Abstraction Reactions from n-Pentane by NO2.

The interaction of fuel with NOx chemistry is important for the construction of the reaction mechanism and engine application. In this work, the reaction pathways of nC5H12 + NO2 were studied by high-level electronic structure calculations (DLPNO-CCSD(T)-F12/cc-pVTZ-F12//B2PLYPD3/cc-pVTZ). The rate constants were calculated by using the multistructural canonical transition-state theory with the Eckart tunneling method (TST/MS-T/ET). The studied condition is in a wide temperature range of 298-2400 K. The influence of MS-T anharmonicity and tunneling effect will be clarified for these site-specific H-abstraction pathways. The result reflects the large deviation introduced by the treatment of MS-T anharmonicity, especially at a high temperature. For the same type of reactions, the rate constants of H-abstraction both occurring at the secondary carbon are not almost identical. The branching ratios show that abstraction from the secondary site forming cis-HONO (R2c) contributes 36-78% to nC5H12 consumption in the temperature range of 298-2400 K. The current results show that the multistructural torsional anharmonicity has a crucial influence on the accurate estimation of branching ratios.

Lithium partition and isotopic fractionation for cation exchange in clay: Insights from molecular simulations

Chemical weathering is one of the most important processes on the Earth’s surface, and lithium isotopes have been used as a promising tracer of silicate weathering, as the secondary clay minerals precipitated during silicate weathering are thought to preferentially incorporate the lighter 6Li isotope. However, to utilize lithium isotopes as a geochemical tracer, it is essential to quantitatively understand the partitioning of Li in clay minerals, including structurally bound Li and exchangeable interlayer ions, as well as the associated isotope fractionation. With classical molecular dynamic simulations and thermodynamic calculations, the cation exchange between Na-/K-montmorillonite and solvated Li in environmental solution are disclosed under different aqueous activity conditions. The Li+ ion, whose hydration ability is stronger than Na+ and K+ ions, is easier to swell the montmorillonite and shows less affinity for the clay. And the calculated cationic selectivity coefficients between Na+/K+ and Li+ ions are Kc (Li → Na) = 2.0, Kc (Li → K) = 10.0 for Arizona-type montmorillonite and Kc (Li → Na) = 1.0, Kc (Li → K) = 4.0 for Wyoming-type montmorillonite, which are comparable to reported experiments. The cationic partition of Li in the clay interlayer decreases when the hydration states of montmorillonite transform from triple-layer to bilayer and then to monolayer. Based on the density functional theory, the Li isotope fractionations for Li in clay interlayer sites and structural octahedral sites are further derived by thermodynamic calculations with path integral molecular dynamic simulation, which deals with the quantum anharmonicity. The behaviors of Li isotope in the bilayer hydration state closely resembles that in the aqueous solution, and only a slight Li isotope fractionation is found between the monolayer hydrated octahedral-substituted clay and solution phase (Δ7Limonolayer-aq = −2.0‰). While Li isotope fractionation between aqueous solution and mica structural sites is significant (Δ7Listructural-aq = −10.4‰) and shows an obvious quantum anharmonicity, which indicates the importance of the anharmonic treatment for Li isotope fractionation calculations. The newly calculated Li equilibrium isotope fractionation between clay and aqueous solution deviates from the experimentally determined fractionation factor, which can be attributed to the kinetic effects. Moreover, on the basis of the surface kinetic reaction model, the Li isotope fractionation becomes more pronounced with higher clay growth rate. By integrating the cationic selectivity coefficient and Li isotope fractionation, it is possible to constrain the fluid composition during Li-clay deposit and diagenesis, and to evaluate the intensity of chemical weathering processes as well as the related carbon cycle.

First-Principles Thermodynamics of CsSnI3.

CsSnI3 is a promising ecofriendly solution for energy harvesting technologies. It exists at room temperature in either a black perovskite polymorph or a yellow 1D double-chain, which irreversibly deteriorates in the air. In this work, we unveil the relative thermodynamic stability between the two structures with a first-principles sampling of the CsSnI3 finite-temperature phase diagram, discovering how it is driven by anomalously large quantum and anharmonic ionic fluctuations. Thanks to a comprehensive treatment of anharmonicity, the simulations deliver a remarkable agreement with known experimental data for the transition temperatures of the orthorhombic, rhombohedral, and cubic perovskite structures and the thermal expansion coefficient. We disclose how the perovskite polymorphs are the ground state above 270 K and discover an abnormal decrease in heat capacity upon heating in the cubic black perovskite. Our results also significantly downplay the Cs+ rattling modes' contribution to mechanical instability. The remarkable agreement with experiments validates our methodology, which can be systematically applied to all metal halides.

Anharmonic Coupling of Stretching Vibrations in Ice: A Periodic VSCF and VCI Description.

The anharmonicity of O-H stretching vibrations of water ice is characterized by use of a periodic implementation of the vibrational self-consistent field (VSCF) and vibrational configuration interaction (VCI) methods, which take phonon-phonon couplings explicitly into account through numerical evaluation of high-order terms of the nuclear potential. The low-temperature, proton-ordered phase of water ice (namely, ice XI) is investigated. The net effect of a coupled anharmonic treatment of stretching modes is not just a rigid blue-shift of the respective harmonic spectral frequencies but rather a complex change of their relative spectral positions, which cannot be captured by simple scaling strategies based on harmonic calculations. The adopted techniques allow for a hierarchical treatment of anharmonic terms of the nuclear potential, which is key to an effective identification of leading factors. We show that the anharmonic independent-mode approximation─only describing the "intrinsic anharmonicity" of the O-H stretches─is unable to capture the correct physics, and that couplings among O-H stretches must be described. Inspection of harmonic normal coordinates allows identification of specific features of the O-H stretching motions which most likely enable strong mode-mode couplings. Finally, by coupling O-H stretches to all other possible modes of ice XI (THz collective vibrations, molecular librations, bendings), we identify specific types of motion which significantly affect O-H stretching states: in particular, molecular librations are found to affect the stretching states more than molecular bendings.

Electronic and mechanical anharmonicities in the vibrational spectra of the H-bonded, cryogenically cooled X- · HOCl (X=Cl, Br, I) complexes: Characterization of the strong anionic H-bond to an acidic OH group.

We report vibrational spectra of the H2-tagged, cryogenically cooled X- · HOCl (X = Cl, Br, and I) ion-molecule complexes and analyze the resulting band patterns with electronic structure calculations and an anharmonic theoretical treatment of nuclear motions on extended potential energy surfaces. The complexes are formed by "ligand exchange" reactions of X- · (H2O)n clusters with HOCl molecules at low pressure (∼10-2mbar) in a radio frequency ion guide. The spectra generally feature many bands in addition to the fundamentals expected at the double harmonic level. These "extra bands" appear in patterns that are similar to those displayed by the X- · HOD analogs, where they are assigned to excitations of nominally IR forbidden overtones and combination bands. The interactions driving these features include mechanical and electronic anharmonicities. Particularly intense bands are observed for the v = 0 → 2 transitions of the out-of-plane bending soft modes of the HOCl molecule relative to the ions. These involve displacements that act to break the strong H-bond to the ion, which give rise to large quadratic dependences of the electric dipoles (electronic anharmonicities) that drive the transition moments for the overtone bands. On the other hand, overtone bands arising from the intramolecular OH bending modes of HOCl are traced to mechanical anharmonic coupling with the v = 1 level of the OH stretch (Fermi resonances). These interactions are similar in strength to those reported earlier for the X- · HOD complexes.

Astrochemical significance and spectroscopy of tetratomic [H, P, S, O

Context. Phosphorus is integral to life on Earth, and its role in the chemistry of the interstellar medium is highly debated and unknown. Only a handful of phosphorus-bearing species have been detected thus far, with the most recent confirmed detection taking place in 2014. The simultaneous detection of molecules such as PO, SH, and OH indicate the possibility of reactive intermediate species existing in the interstellar medium and circumstellar envelopes of evolved stars. To explore this possibility, the [H, P, S, O] tetratomic isomer family was characterized using high level ab initio methods. Aims. The aim of this study is to provide rotational, vibrational, and electronic spectroscopic data to drive experimental and observational detection of new phosphorus and sulfur-bearing molecules. Additionally, chemical pathways are explored to explain possible reservoirs for the as of yet undetected PH and PS diatomic molecules. Methods. Coupled cluster quantum chemistry methods were used to calculate the equilibrium electronic structure followed by the anharmonic treatment of the cubic and quartic force fields to obtain accurate rotational and vibrational data. Møller–Plesset perturbation theory in conjunction with coupled cluster methods were used to explore bimolecular reaction pathways. Multi-reference methods were then used to characterize the photochemical pathways of the excited electronic states and simulate the electronic absorption spectrum. Results. The reaction between detected molecules SH and PO is highly exothermic and forms the HSPO isomer. Deeply submerged transition state barriers allow for facile isomerization to other isomers, especially HOPS. The dominant photochemical process predicted for HOPS is dissociation to form OH + PS, while that of HSPO is a combination of photodissociation to form H + SPO and SH + PO, depending on the wavelength of light absorbed. If PH and PS are formed in the early outflows from evolved stars, bimolecular reactions may act as a reservoir and partially account for their lack of detection to date. The electronic absorption spectrum is predicted to be congested in the 175–200 nm region for both HOPS and HSPO. Differentiating peaks exist &gt;400 nm, which can be used for spectral assignment. Vibrationally corrected rotational constants and anharmonic vibrational frequencies were calculated to assist in the laboratory and observation identification of the most stable molecules. The PO stretch is predicted to be the most intense vibrational mode in both HOPS isomers, and a frequency difference of 20 cm−1 may prove to help differentiate the conformers in an experimental spectrum.

Tailored anharmonic-harmonic vibrational profiles for fluorescent biomarkers.

We propose a hybrid anharmonic-harmonic scheme for vibrational broadenings, which embeds a reduced-space vibrational configuration interaction (VCI) anharmonic wave function treatment in the independent-mode displaced harmonic oscillator (IMDHO) model. The resulting systematically-improvable VCI-in-IMDHO model allows including the vibronic effects of all vibrational degrees of freedom, while focusing the effort on the important degrees of freedom with minimal extra computational effort compared to a reduced-space VCI treatment. We show for oligothiophene examples that the VCI-in-IMDHO approach can yield accurate vibrational profiles employing smaller vibrational spaces in the VCI part than the reduced-space VCI approach. By this, the VCI-in-IMDHO model enables accurate calculation of vibrational profiles of common fluorescent dyes with more than 100 vibrational degrees of freedom. We illustrate this for three examples of fluorescent biomarkers of current interest. These are the oligothiophene-based fluorescent dye called HS84, 1,4-diphenylbutadiene, and an anthracene diimide. For all examples, we assess the impact of the anharmonic treatment on the vibrational broadening, which we find to be more pronounced for the intensities than for the peak positions.

Infrared Spectroscopy of Jet-cooled “GrandPAHs” in the 3–100 μm Region

Abstract Although large polycyclic aromatic hydrocarbons (PAHs) are likely to be responsible for IR emission of gaseous and dusty regions, their neutral experimental high-resolution gas-phase IR spectra—needed to construct accurate astronomical models—have so far remained out of reach because of their nonvolatility. Applying laser desorption to overcome this problem, we report here the first IR spectra of the jet-cooled large PAHs coronene (C24H12), peropyrene (C26H14), ovalene (C32H14), and hexa(peri)benzocoronene (C42H18) in the 3–100 μm region. Apart from providing experimental spectra that can be compared directly to astronomical data, such IR spectra are crucial for assessing the accuracy of theoretically predicted spectra used to interpret interstellar IR emission. Here we use the experimental spectra to evaluate the performance of conventional calculations using the harmonic approximation, as well as calculations with an anharmonic (GVPT2) treatment. The harmonic prediction agrees well with the experiment between 100 and 1000 cm−1 (100 and 10 μm) but shows significant shortcomings in the combination band (1600–2000 cm−1, 6.25–5 μm) and CH-stretch (2950–3150 cm−1, 3.4–3.17 μm) regions. Especially the CH-stretch region is known to be dominated by the effects of anharmonicity, and we find that large PAHs are no exception. However, for the CH out-of-plane region (667–1000 cm−1, 15–10 μm) the anharmonic treatment that significantly improves the predicted spectra for small PAHs leads to large and unrealistic frequency shifts, and intensity changes for large PAHs, thereby rendering the default results unreliable. A detailed analysis of the results of the anharmonic treatment suggests a possible route for improvement, although the underlying cause for the large deviations remains a challenge for theory.

Meta‐analysis of uniform scaling factors for harmonic frequency calculations

AbstractVibrational frequency calculations performed under the harmonic approximation are widespread across chemistry. However, it is well‐known that the calculated harmonic frequencies tend to systematically overestimate experimental fundamental frequencies; a limitation commonly overcome with multiplicative scaling factors. In practice, multiplicative scaling factors are derived for each individual model chemistry choice (i.e., a level of theory and basis set pair), where performance is judged by, for example, the root‐mean square error (RMSE) between the predicted scaled and experimental frequencies. However, despite the overwhelming number of scaling factors reported in the literature and model chemistry approximations available, there is little guidance for users on appropriate model chemistry choices for harmonic frequency calculations. Here, we compile and analyze the data for 1495 scaling factors calculated using 141 levels of theory and 109 basis sets. Our meta‐analysis of this data shows that scaling factors and RMSE approach convergence with only hybrid functionals and double‐zeta basis sets, with anharmonicity error already dominating model chemistry errors. Noting inconsistent data and the lack of independent testing, we can nevertheless conclude that a minimum error of 25 cm−1—arising from insufficiently accurate treatment of anharmonicity—is persistent regardless of the model chemistry choice. Based on the data we compiled and cautioning the need for a future systematic benchmarking study, we recommend ωB97X‐D/def2‐TZVP for most applications and B2PLYP/def2‐TZVPD for superior intensity predictions. With a smaller benchmark set, direct comparison prefers ωB97X‐D/6‐31G* to B3LYP/6‐31G*.This article is categorized under: Electronic Structure Theory &gt; Density Functional Theory Theoretical and Physical Chemistry &gt; Spectroscopy

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. The thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 $\phantom{\rule{0.16em}{0ex}}{\mathrm{Wm}}^{\ensuremath{-}1}\phantom{\rule{0.16em}{0ex}}{\mathrm{K}}^{\ensuremath{-}1}$ are estimated at ambient and extreme conditions ($T=4000\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ and $P=135\phantom{\rule{0.16em}{0ex}}\mathrm{GPa}$), respectively.

Anharmonicity measure for materials

Theoretical frameworks used to qualitatively and quantitatively describe nuclear dynamics in solids are often based on the harmonic approximation. However, this approximation is known to become inaccurate or to break down completely in many modern functional materials. Interestingly, there is no reliable measure to quantify anharmonicity so far. Thus, a systematic classification of materials in terms of anharmonicity and a benchmark of methodologies that may be appropriate for different strengths of anharmonicity is currently impossible. In this work, we derive and discuss a statistical measure that reliably classifies compounds across temperature regimes and material classes by their "degree of anharmonicity". This enables us to distinguish "harmonic" materials, for which anharmonic effects constitute a small perturbation on top of the harmonic approximation, from strongly "anharmonic" materials, for which anharmonic effects become significant or even dominant and the treatment of anharmonicity in terms of perturbation theory is more than questionable. We show that the analysis of this measure in real and reciprocal space is able to shed light on the underlying microscopic mechanisms, even at conditions close to, e.g., phase transitions or defect formation. Eventually, we demonstrate that the developed approach is computationally efficient and enables rapid high-throughput searches by scanning over a set of several hundred binary solids. The results show that strong anharmonic effects beyond the perturbative limit are not only active in complex materials or close to phase transitions, but already at moderate temperatures in simple binary compounds.

The first microsolvation step for furans: New experiments and benchmarking strategies

The site-specific first microsolvation step of furan and some of its derivatives with methanol is explored to benchmark the ability of quantum-chemical methods to describe the structure, energetics, and vibrational spectrum at low temperature. Infrared and microwave spectra in supersonic jet expansions are used to quantify the docking preference and some relevant quantum states of the model complexes. Microwave spectroscopy strictly rules out in-plane docking of methanol as opposed to the top coordination of the aromatic ring. Contrasting comparison strategies, which emphasize either the experimental or the theoretical input, are explored. Within the harmonic approximation, only a few composite computational approaches are able to achieve a satisfactory performance. Deuteration experiments suggest that the harmonic treatment itself is largely justified for the zero-point energy, likely and by design due to the systematic cancellation of important anharmonic contributions between the docking variants. Therefore, discrepancies between experiment and theory for the isomer abundance are tentatively assigned to electronic structure deficiencies, but uncertainties remain on the nuclear dynamics side. Attempts to include anharmonic contributions indicate that for systems of this size, a uniform treatment of anharmonicity with systematically improved performance is not yet in sight.

Computational analysis of the proton-bound acetonitrile dimer, (ACN)2 H.

In atmospheric pressure ionization mass spectrometry the theoretical thermodynamic treatment of proton-bound cluster stabilities helps us to understand the prevailing chemical processes. However, such calculations are rather challenging because low-barrier internal rotations and strong anharmonicity of the hydrogen bonds cause the breakdown of the usually applied harmonic approximation. Even the implemented anharmonic treatment in standard ab initio software failed in the case of (ACN)2 H+ . For a case study of the proton-bound acetonitrile dimer, (ACN)2 H+ , we scan the potential energy surface (PES) for the internal rotation and the proton movement in all three spatial directions. We correct the partition functions by treating the internal rotation as a free rotor and by solving the nuclear Schrödinger equation explicitly for the proton movement. An additional PES scan for the dissociation surface further improves the understanding of the cluster behavior. The internal rotation is essentially barrier free (V0 = 2.6 × 10-6 eV) and the proton's movement between the two nitrogen atoms follows a quartic rather than quadratic potential. As a figure of merit we calculate the free dissociation enthalpy of the dimer. Our description significantly improves the standard results from about 118.3 kJ/mol to 99.6 kJ/mol, compared with the experimentally determined value of 92.2 kJ/mol. The dissociation surface reveals strong crosstalk between modes and is essentially responsible for the observed errors. The presented corrections to the partition functions significantly improve their accuracy and are rather easy to implement. In addition, this work stresses the importance of alternative theoretical methods for proton-bound cluster systems besides the standard harmonic approximations.

Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene.

Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S1 ← S0 and S2 ← S0 absorption and S2 → S0 emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.

Gas-Phase Vibrational Spectroscopy of the Hydrocarbon Cations l-C3H+, HC3H+, and c-C3H2+: Structures, Isomers, and the Influence of Ne-Tagging.

We report the first gas-phase vibrational spectra of the hydrocarbon ions C3H+ and C3H2+. The ions were produced by electron impact ionization of allene. Vibrational spectra of the mass-selected ions tagged with Ne were recorded using infrared predissociation spectroscopy in a cryogenic ion trap instrument using the intense and widely tunable radiation of a free electron laser. Comparison of high-level quantum chemical calculations and resonant depletion measurements revealed that the C3H+ ion is exclusively formed in its most stable linear isomeric form, whereas two isomers were observed for C3H2+. Bands of the energetically favored cyclic c-C3H2+ are in excellent agreement with calculated anharmonic frequencies, whereas for the linear open-shell HCCCH+ (2Πg) a detailed theoretical description of the spectrum remains challenging because of Renner–Teller and spin–orbit interactions. Good agreement between theory and experiment, however, is observed for the frequencies of the stretching modes for which an anharmonic treatment was possible. In the case of linear l-C3H+, small but non-negligible effects of the attached Ne on the ion fundamental band positions and the overall spectrum were found.