Theoretical and experimental results for the surface core-level binding energy, BE, shifts, SCLS, for MgO(100) are presented and the anomalous O(1s) SCLS is interpreted in terms of the surface electronic structure. While the Mg(2p) surface BE shifts to a higher value than bulk by ≈1 eV as expected from the different surface and bulk Madelung potentials, the O(1s) SCLS is almost 0 rather than ≈-1 eV, expected from the Madelung potentials. The distortion of the surface atoms from the spherical symmetry of the bulk Mg and O atoms is examined by a novel theoretical procedure. The anomalous O SCLS is shown to arise from the increase of the effective size of surface O anions.

The design of accurate helium-solute interaction potentials for the simulation of chemically complex molecules solvated in superfluid helium has long been a cumbersome task due to the rather weak but strongly anisotropic nature of the interactions. We show that this challenge can be met by using a combination of an effective pair potential for the He-He interactions and a flexible high-dimensional neural network potential (NNP) for describing the complex interaction between helium and the solute in a pairwise additive manner. This approach yields an excellent agreement with a mean absolute deviation as small as 0.04 kJ mol for the interaction energy between helium and both hydronium and Zundel cations compared with coupled cluster reference calculations with an energetically converged basis set. The construction and improvement of the potential can be performed in a highly automated way, which opens the door for applications to a variety of reactive molecules to study the effect of solvation on the solute as well as the solute-induced structuring of the solvent. Furthermore, we show that this NNP approach yields very convincing agreement with the coupled cluster reference for properties like many-body spatial and radial distribution functions. This holds for the microsolvation of the protonated water monomer and dimer by a few helium atoms up to their solvation in bulk helium as obtained from path integral simulations at about 1 K.

Superfluid helium has not only fascinated scientists for centuries but is also the ideal matrix for the investigation of chemical systems under ultra-cold conditions in helium nanodroplet isolation experiments. Together with related experimental techniques such as helium tagging photodissociation spectroscopy, these methods have provided unique insights into many interesting systems. Complemented by theoretical work, they were additionally able to greatly expand our general understanding of manifestations of superfluid behavior in finite sized clusters and their response to molecular impurities. However, most theoretical studies up to now have not included the reactivity and flexibility of molecular systems embedded in helium. In this perspective, the theoretical foundation of simulating fluxional molecules and reactive complexes in superfluid helium is presented in detail. Special emphasis is put on recent developments for the converged description of both the molecular interactions and the quantum nature of the nuclei at ultra-low temperatures. As a first step, our hybrid path integral molecular dynamics/bosonic path integral Monte Carlo method is reviewed. Subsequently, methods for efficient path integral sampling tailored for this hybrid coupling scheme are discussed while also introducing new developments to enhance the accurate incorporation of the solute⋯solvent coupling. Finally, highly accurate descriptions of the interactions in solute⋯helium systems using machine learning techniques are addressed. Our current automated and adaptive fitting procedures to parameterize high-dimensional neural network potentials for both the full-dimensional potential energy surface of solutes and the solute⋯solvent interaction potentials are concisely presented. They are demonstrated to faithfully represent many-body potential functions able to describe chemically complex and reactive solutes in helium environments seamlessly from one He atom up to bulk helium at the accuracy level of coupled cluster electronic structure calculations. Together, these advances allow for converged quantum simulations of fluxional and reactive solutes in superfluid helium under cryogenic conditions.

Path integral methods provide a rigorous and systematically convergent framework to include the quantum mechanical nature of atomic nuclei in the evaluation of the equilibrium properties of molecules, liquids or solids at finite temperature. Such nuclear quantum effects often are already significant for light nuclei at room temperature, but become crucial at cryogenic temperatures such as those provided by superfluid helium as a solvent. Unfortunately, the cost of converged path integral simulations increases significantly upon lowering the temperature so that the computational burden of simulating matter at the typical superfluid helium temperatures becomes prohibitive. Here we investigate how accelerated path integral techniques based on colored noise generalized Langevin equations, in particular the so called PIGLET variant, perform in this extreme quantum regime using as an example the quasi rigid methane molecule and its highly fluxional protonated cousin, CH + 5 . We show that the PIGLET technique gives a speedup of two orders of magnitude in the evaluation of structural observables and quantum kinetic energy at ultralow temperatures. Moreover, we computed the spatial spread of the quantum nuclei in CH 4 to illustrate the limits of using such colored noise thermostats close to the many body quantum ground state.

Messenger infrared predissociation spectroscopy is an powerful technique to investigate elusive molecules. Despite accumulating evidence that He-tagging has a minor effect on structure due to its weak interactions, its impact on intramolecular large-amplitude motion is much less clear. Using one of the most floppy molecules known, protonated methane, we investigate to what extent helium tags perturb its salient fluxionality, underlying full hydrogen scrambling. Bosonic path integral quantum simulations show that the impact of attaching up to four He messengers to CH is negligible compared to the non-tagged bare species. This non-invasive nature is unrelated to the bosonic character of He. Yet, hydrogen scrambling has a pronounced effect on the microsolvation shell pattern of CH ⋅He that emerges as a result of tagging. These findings should provide impetus to explore the power of He-tagging, particularly in the realm of weakly bound intermolecular complexes.