Temperature and anharmonic effects on the infrared absorption spectrum from a quantum statistical approach: application to naphthalene

Aromatic and Acetylenic C-H or C-D Stretching Bands Anharmonicity Detection of Phenylacetylene by UV Laser-Induced Vibrational Emission

The anharmonic infrared (IR) emission spectra of phenylacetylene C₆H₅CCH and an isotopologue C₆H₅CCD induced by 193 nm UV excitation have been investigated in the gas phase. The study has been operated with a homemade IR spectrometer enabling to record time- and wavelength-resolved spectra between 2.5 and 4.5 μm, emitted all along the collisional cooling. The analysis is supported by a kinetic Monte Carlo simulation in the vibrational harmonic approximation. For both species, the anharmonic shifts of the acetylenic C-H or C-D stretching modes and the aromatic C-H stretching modes are studied for band positions and bandwidths in terms of the internal energy. For C₆H₅CCD, the internal energy dependence of the emission intensity band ratio is investigated and rationalized. This work demonstrates the potential of time-resolved IR emission spectroscopy to explore anharmonicity of astrophysically relevant molecules.

Anharmonicity and the IR Emission Spectrum of Neutral Interstellar PAH Molecules

The characteristics of the CH stretching and out-of-plane bending modes in polycyclic aromatic hydrocarbon molecules are investigated using anharmonic density functional theory (DFT) coupled to a vibrational second-order perturbation treatment taking resonance effects into account. The results are used to calculate the infrared emission spectrum of vibrationally excited species in the collision-less environment of interstellar space. This model follows the energy cascade as the molecules relax after the absorption of a UV photon in order to calculate the detailed profiles of the infrared bands. The results are validated against elegant laboratory spectra of polycyclic aromatic hydrocarbon absorption and emission spectra obtained in molecular beams. The factors which influence the peak position, spectral detail, and relative strength of the CH stretching and out-of-plane bending modes are investigated, and detailed profiles for these modes are derived. These are compared to observations of astronomical objects in space, and the implications for our understanding of the characteristics of the molecular inventory of space are assessed.

Anharmonic Infrared Spectra of Thermally Excited Pyrene (C16H10): A Combined View of DFT-Based GVPT2 with AnharmonicCaOs, and Approximate DFT Molecular dynamics with DemonNano

The study of the Aromatic Infrared Bands (AIBs) in astronomical environments has opened interesting spectroscopic questions on the effect of anharmonicity on the infrared (IR) spectrum of hot polycyclic aromatic hydrocarbons (PAHs) and related species in isolated conditions. The forthcoming James Webb Space Telescope will unveil unprecedented spatial and spectral details in the AIB spectrum; significant advancement is thus necessary now to model the infrared emission of PAHs, their presumed carriers, with enough detail to exploit the information content of the AIBs. This requires including anharmonicity in such models, and to do so systematically for all species included, requiring a difficult compromise between accuracy and efficiency. We performed a benchmark study to compare the performances of two methods in calculating anharmonic spectra, comparing them to available experimental data. One is a full quantum method, AnharmoniCaOs, relying on an ab initio potential, and the other relies on Molecular Dynamics simulations using a Density Functional based Tight Binding potential. The first one is found to be very accurate and detailed, but it becomes computationally very expensive for increasing temperature; the second is faster and can be used for arbitrarily high temperatures, but is less accurate. Still, its results can be used to model the evolution with temperature of isolated bands. We propose a new recipe to model anharmonic AIB emission using minimal assumptions on the general behaviour of band positions and widths with temperature, which can be defined by a small number of empirical parameters. Modelling accuracy will depend critically on these empirical parameters, allowing for an incremental improvement in model results, as better estimates become gradually available.

Experimental Approach to the Study of Anharmonicity in the Infrared Spectrum of Pyrene from 14 to 723 K

Quantifying the effect of anharmonicity on the infrared spectrum of large molecules such as polycyclic aromatic hydrocarbons (PAHs) at high temperatures is the focus of a number of theoretical and experimental studies, many of them motivated by astrophysical applications. We recorded the IR spectrum of pyrene C16H10 microcrystals embedded in KBr pellets over a wide range of temperatures (14-723 K) and studied the evolution of band positions, widths, and integrated intensities with temperature. We identified jumps for some of the spectral characteristics of some bands in the 423-473 K range. These were attributed to a change of phase from crystal to molten in condensed pyrene, which appears to affect more strongly bands involving large CH motions. Empirical anharmonicity factors that quantify the linear evolution of band positions and widths with temperature for values larger than ∼150-250 K, depending on the band, were retrieved from both phases and averaged to provide recommended values for these anharmonicity factors. The derived values were found to be consistent with available gas phase data. We conclude about the relevance of the methodology to produce data that can be compared with calculated anharmonic IR spectra and provide input for models that simulate the IR emission of astro-PAHs.

Anharmonic vibrational spectroscopy of polycyclic aromatic hydrocarbons (PAHs)

While powerful techniques exist to accurately account for anharmonicity in vibrational molecular spectroscopy, they are computationally very expensive and cannot be routinely employed for large species and/or at non-zero vibrational temperatures. Motivated by the study of Polycyclic Aromatic Hydrocarbon (PAH) emission in space, we developed a new code, which takes into account all modes and can describe all infrared transitions including bands becoming active due to resonances as well as overtone, combination, and difference bands. In this article, we describe the methodology that was implemented and discuss how the main difficulties were overcome, so as to keep the problem tractable. Benchmarking with high-level calculations was performed on a small molecule. We carried out specific convergence tests on two prototypical PAHs, pyrene (C₁₆H₁₀) and coronene (C₂₄H₁₂), aiming at optimising tunable parameters to achieve both acceptable accuracy and computational costs for this class of molecules. We then report the results obtained at 0 K for pyrene and coronene, comparing the calculated spectra with available experimental data. The theoretical band positions were found to be significantly improved compared to harmonic density functional theory calculations. The band intensities are in reasonable agreement with experiments, the main limitation being the accuracy of the underlying calculations of the quartic force field. This is a first step toward calculating moderately high-temperature spectra of PAHs and other similarly rigid molecules using Monte Carlo sampling.

Hydrogenation of Furfural with Nickel Nanoparticles Stabilized on Nitrogen-Rich Carbon Core-Shell and Its Transformations for the Synthesis of γ-Valerolactone in Aqueous Conditions

In this article, we report the synthesis of nitrogen-rich carbon layer-encapsulated Ni(0) nanoparticles as a core-shell structure (Ni@N/C-g-800) for the catalytic hydrogenation of furfural to furfuryl alcohol. The nickel nanoparticles were stabilized by the nitrogen-rich graphitic framework, which formed during the agitation of nickel acetate-impregnated cucurbit[6]uril surface in a reducing atmosphere. Furthermore, the catalyst was characterized using various physicochemical methods such as powder X-ray diffraction, Raman, field emission-scanning electron microscopy, high-resolution transmission electron microscopy, X-ray photoelectron spectroscopy, Brunauer-Emmett-Teller surface area, CO₂-temperature-programmed desorption, inductive coupled plasma, and CHN analyses. The nitrogen-rich environment of the solid support with metallic Ni nanoparticles was found to be active and selective for the catalytic hydrogenation of furfural with molecular H₂ in an aqueous medium at 100 °C. To understand the reaction mechanism, the diffuse reflectance infrared Fourier transform study was performed, which revealed that the C═O bond is activated in the presence of a catalyst. In addition, we have extended our methodology toward the synthesis of "levulinic acid" and "γ-valerolactone", by successive hydrolysis and hydrogenation of furfuryl alcohol and levulinic acid, respectively, in an aqueous medium. Moreover, the heterogeneous catalysts used in all of the three consecutive steps help in recovery and recycling of the catalyst and easy separation of products.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Nonlinear effects in infrared action spectroscopy of silicon and vanadium oxide clusters: experiment and kinetic modeling

Nonlinear effects in infrared action spectroscopy are experimentally quantified and successfully modeled for different inorganic clusters.

Accurate modeling of infrared multiple photon dissociation spectra: the dynamical role of anharmonicities

The dynamical response of a molecular system to a macropulse typically produced by a free-electron laser is theoretically modeled over experimentally long times, within a realistic kinetic Monte Carlo framework that incorporates absorption, stimulated emission, spontaneous emission, and dissociation events. The simulation relies on an anharmonic potential energy surface obtained from quantum chemistry calculations. Application to cationic naphthalene yields a better agreement with measurements than the anharmonic linear absorption spectrum, thus emphasizing the importance of specific dynamical effects on the spectral properties.

Probing the spin multiplicity of gas-phase polycyclic aromatic hydrocarbons through their infrared emission spectrum: a theoretical study

The anharmonic infrared emission spectrum following an optical excitation has been calculated for a variety of polycyclic aromatic hydrocarbon molecules in their ground singlet electronic state or in their triplet state. The computational protocol relies on second-order perturbation theory and involves a quartic vibrational Hamiltonian, the vibrational quantum numbers being sampled according to a Monte Carlo procedure. In the case of neutral naphthalene, the IR spectrum obtained in the (ground) singlet state differs significantly from the spectrum in the triplet state, especially for out-of-plane CH bending modes. Although not as prominent, spectral differences in larger molecules are still observable.

A new Monte Carlo method for getting the density of states of atomic cluster systems

A novel Monte Carlo flat histogram algorithm is proposed to get the classical density of states in terms of the potential energy, g(E(p)), for systems with continuous variables such as atomic clusters. It aims at avoiding the long iterative process of the Wang-Landau method and controlling carefully the convergence, but keeping the ability to overcome energy barriers. Our algorithm is based on a preliminary mapping in a series of points (called a σ-mapping), obtained by a two-parameter local probing of g(E(p)), and it converges in only two subsequent reweighting iterations on large intervals. The method is illustrated on the model system of a 432 atom cluster bound by a Rydberg type potential. Convergence properties are first examined in detail, particularly in the phase transition zone. We get g(E(p)) varying by a factor 10(3700) over the energy range [0.01 < E(p) < 6000 eV], covered by only eight overlapping intervals. Canonical quantities are derived, such as the internal energy U(T) and the heat capacity C(V)(T). This reveals the solid to liquid phase transition, lying in our conditions at the triple point. This phase transition is further studied by computing a Lindemann-Berry index, the atomic cluster density n(r), and the pressure, demonstrating the progressive surface melting at this triple point. Some limited results are also given for 1224 and 4044 atom clusters.

Temperature effects on the rovibrational spectra of pyrene-based PAHs

Absorption infrared spectra have been computed for a variety of polycyclic aromatic hydrocarbon molecules of the pyrene family, taking into account anharmonicity and temperature effects, rovibrational quantization, and couplings. The energy levels are described by a second-order perturbative expansion of the rovibrational Hamiltonian in the vibrational and rotational quantum numbers, as relevant for a symmetric-top molecule, with ingredients obtained from quantum chemistry calculations. Multicanonical Monte Carlo simulations are carried out to compute bidimensional IR intensity histograms as a function of total energy and vibrational frequency, which then provide the absorption spectrum at arbitrary temperatures via a Laplace transformation. The main spectral features analyzed for neutral, anionic, and cationic pyrene indicate a strong dependence on temperature, in agreement with existing laboratory experiments, and a significant contribution of rotational degrees of freedom to the overall broadenings. The spectral shifts and broadenings reveal some sensitivity of anharmonicities to the charge and protonation states and, in the case of protonated pyrene and pyrenyl cation, on possible isomers and between aromatic and aliphatic C-H bands. Implications of the present work to the general issue of interstellar emission features are discussed.

Molecular dynamics simulations on [FePAH]+ π-complexes of astrophysical interest: anharmonic infrared spectroscopy

In this article, classical Born-Oppenheimer molecular dynamics (MD) simulations in the microcanonical ensemble are performed on neutral and cationic polycyclic aromatic hydrocarbon (PAH) species, focusing on [FePAH](+)π-complexes. Their anharmonic mid-infrared (mid-IR) spectra in the classical approximation are derived. This approach allows us to describe the influence of the energy of a system on its IR spectrum in terms of band-shifts and broadenings. The MD simulations are performed on a potential energy surface (PES) described at the self-consistent-charge density functional tight-binding level of theory. The PES is benchmarked on DFT calculations, showing the validity of the approach for complexes of Fe(+) with PAHs larger than coronene (C(24)H(12)) that are of astrophysical interest. MD simulations at high temperature show the occurrence of the diffusion of the Fe cation on the surface of the PAH. It proceeds through the edge of the carbon skeleton which is the lowest energy pathway presenting barriers smaller than 1 eV. Although only qualitative information on the band broadenings can be obtained, we show that the dependence of the computed positions of the main bands of [C(24)H(12)](0/+)and [FeC(24)H(12)](+)π-complexes on temperature can be fit by linear laws. The spectral trends determined for [FeC(24)H(12)](+) are compared to those of N-substituted [C(23)NH(12)](+)and [SiC(24)H(12)](+)π-complexes of astrophysical interest.