Physical conditions in PDRs revealed by IGRINS ro-vibrational H2 observations

We compare H2 high ro-vibrational observations of 5 PDRs (S140, IC63, Horsehead Nebula, NGC 2023 and Orion Bar) obtained with the IGRINS spectrograph to PDR models produced with the Meudon PDR code, in order to estimate the physical conditions and to constrain the physical and chemical processes that govern PDRs. We use newly flux-calibrated and extinction-corrected IGRINS H2 observations covering the 1.45-2.45 um range at a resolution of 45000, and adjust Meudon PDR models using both a simple chi2 minimization approach and a robust method using Bayesian inversion followed by posterior exploration via an advanced MCMC technique. The PDR models use specific incident FUV spectra for each PDRs based on theoretical stellar spectra from the Pollux database. The Meudon PDR code is able to reproduce more than 85% of the observed H2 ro-vibrational line intensities within a factor of two. We find that (1) a realistic modeling of the incident FUV field (both in term of geometry and of spectral shape) and (2) the inclusion of recent data on collisional de-excitation rates for high vibrational levels of H2 are key to this success. H2 ro-vibrational emission lines are found to provide good constraints on both the thermal pressure and the incident FUV field strength, G0, although with a non-negligible remaining degeneracy in most PDRs. Additional constraints, such as the spatial scales of the PDR derived from ALMA or JWST observations, are found to be able to lift this degeneracy. For low excitation PDRs, the observations provide evidence that nascent H2 molecules formed on grains have relatively low rotational energy and high vibrational energy, as predicted by theoretical and experimental studies of the surface formation of H2.