Probing Ion-Receptor Interactions in Halide Complexes of Octamethyl Calix[4]Pyrrole.

Ion receptors are molecular hosts that bind ionic guests, often with great selectivity. The interplay of solvation and ion binding in anion host-guest complexes in solution governs the binding efficiency and selectivity of such ion receptors. To gain molecular-level insight into the intrinsic binding properties of octamethyl calix[4]pyrrole (omC4P) host molecules with halide guest ions, we performed cryogenic ion vibrational spectroscopy (CIVS) of omC4P in complexes with fluoride, chloride, and bromide ions. We interpret the spectra using density functional theory, describing the infrared spectra of these complexes with both harmonic and anharmonic second-order vibrational perturbation theory (VPT2) calculations. The NH stretching modes of the pyrrole moieties serve as sensitive probes of the ion binding properties, as their frequencies encode the ion-receptor interactions. While scaled harmonic spectra reproduce the experimental NH stretching modes of the chloride and bromide complexes in broad strokes, the high proton affinity of fluoride introduces strong anharmonic effects. As a result, the spectrum of F-·omC4P is not even qualitatively captured by harmonic calculations, but it is recovered very well by VPT2 calculations. In addition, the VPT2 calculations recover the intricate coupling of the NH stretching modes with overtones and combination bands of CH stretching and NH bending modes and with low-frequency vibrations of the omC4P macrocycle, which are apparent for all of the halide ion complexes investigated here. A comparison of the CIVS spectra with infrared spectra of solutions of the same ion-receptor complexes in d3-acetonitrile and d6-acetone shows how ion solvation changes the ion-receptor interactions for the different halide ions.

Correlations between the Structures and Spectra of Protonated Water Clusters.

Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.

Infrared Spectrum of the Pyrene Anion in the CH Stretching Region.

In this work, we report the infrared spectrum of the pyrene anion, measured using messenger tagging with up to three Ar atoms. We assign the spectrum using density functional theory and vibrational perturbation theory. We discuss our results in the context of computed and experimental spectra from the literature as well as recent observations from astronomical sources, addressing the question of whether polycyclic aromatic hydrocarbon anions could contribute to the strong infrared emission bands at 3.29 µm from carbon-rich regions of space.

Infrared spectroscopy of the syn-methyl-substituted Criegee intermediate: A combined experimental and theoretical study.

An IR-vacuum ultraviolet (VUV) ion-dip spectroscopy method is utilized to examine the IR spectrum of acetaldehyde oxide (CH3CHOO) in the overtone CH stretch (2νCH) spectral region. IR activation creates a depletion of the ground state population that reduces the VUV photoionization signal on the parent mass channel. IR activation of the more stable and populated syn-CH3CHOO conformer results in rapid unimolecular decay to OH + vinoxy products and makes the most significant contribution to the observed spectrum. The resultant IR-VUV ion-dip spectrum of CH3CHOO is similar to that obtained previously for syn-CH3CHOO using IR action spectroscopy with UV laser-induced fluorescence detection of OH products. The prominent IR features at 5984 and 6081 cm-1 are also observed using UV + VUV photoionization of OH products. Complementary theoretical calculations utilizing a general implementation of second-order vibrational perturbation theory provide new insights on the vibrational transitions that give rise to the experimental spectrum in the overtone CH stretch region. The introduction of physically motivated small shifts of the harmonic frequencies yields remarkably improved agreement between experiment and theory in the overtone CH stretch region. The prominent features are assigned as highly mixed states with contributions from two quanta of CH stretch and/or a combination of CH stretch with an overtone in mode 4. The generality of this approach is demonstrated by applying it to three different levels of electronic structure theory/basis sets, all of which provide spectra that are virtually indistinguishable despite showing large deviations prior to introducing the shifts to the harmonic frequencies.

Exploring the Origins of the Intensity of the OH Stretch-HOH Bend Combination Band in Water.

While the intensity of the OH stretching fundamental transition is strongly correlated to hydrogen-bond strength, the intensity of the corresponding transition to the state with one quantum of excitation in both the OH stretching and HOH bending vibrations in the same water molecule shows a much weaker sensitivity to the hydrogen-bonding environment. The origins of this difference are explored through analyses of the contributions of terms in the expansion of the dipole moment to the calculated intensity. It is found that the leading contribution to the stretch-bend intensity involves the second derivative of the dipole moment with respect to the OH bond length and HOH angle. While this is not surprising, the insensitivity of this derivative to the hydrogen-bonding environment is unexpected. Possible contributions of mode mixing are also explored. While mode mixing leads to splittings of the energies of nearly degenerate excited states, it does not result in significant changes in the sum of the intensities of these transitions. Analysis of changes in the partial charges on the hydrogen atoms upon displacement of the HOH angles shows that these charges generally increase with increasing HOH angle. This effect is partially canceled by a decrease in the charge of the hydrogen atom when a hydrogen bond is broken. The extent of this cancellation increases with the hydrogen bond strength, which is reflected in the observed insensitivity of the intensity of the stretch-bend transition to hydrogen-bond strength.

Isotope Effects on Ground and Excited States of Ethyl Cation, H+(C2H4).

The structure and spectra of ethyl cation, H+(C2H4), and its deuterated analogues are investigated using diffusion Monte Carlo (DMC). These calculations all show that the ground state wave function for H+(C2H4) is localized near the minimum energy configuration in which the excess proton is in a bridging configuration, although the amplitude of the vibrational motions of the bridging proton is large. Deuteration of the bridging proton reduces the amplitude of this motion, while deuteration of only the ethylenic hydrogen atoms in H+(C2D4) has little effect on the amplitude of the motion of the bridging proton. Excited states that are accessed by spectroscopically observed transitions in H+(C2H4) are calculated using fixed-node DMC. The calculated and measured frequencies for the states with one quantum of excitation in the ethylenic CH stretching vibrations show good agreement. We also explore the excited state with one quantum of excitation in the proton transfer vibration of the bridging proton and obtain a frequency of 616 cm-1 for H+(C2H4). This frequency increases to 629 cm-1 in H+(C2D4). Deuteration decreases this frequency to 491 and 495 cm-1 in D+(C2H4) and D+(C2D4), respectively. The effects of partial deuteration on the frequencies of the CH stretching vibrations, and the corresponding probability amplitudes are also explored. Finally, we report the vibrationally averaged rotational constants for the four isotopologues of ethyl cation considered in this study.

High-Resolution Photoelectron Spectroscopy of Vibrationally Excited Vinoxide Anions.

High-resolution photoelectron spectra of vibrationally pre-excited vinoxide anions (CH2CHO-) are reported using the recently developed IR-cryo-SEVI technique. This method is combined with a newly developed implementation of vibrational perturbation theory that can readily identify relevant anharmonic couplings among nearly degenerate vibrational states. IR-cryo-SEVI spectra are obtained by resonant infrared excitation of vinoxide anions via the fundamental C-O (ν4, 1566 cm-1) or isolated C-H (ν3, 2540 cm-1) stretching vibrations prior to photodetachment. Excitation of the ν4 mode leads to a well-resolved photoelectron spectrum that is in excellent agreement with a harmonic Franck-Condon simulation. Excitation of the higher-energy ν3 mode results in a more complicated spectrum that requires consideration of the calculated anharmonic resonances in both the anion and the neutral. From this analysis, information about the zeroth-order states that contribute to the nominal ν3 wave function in the anion is obtained. In the neutral, we observe anharmonic splitting of the ν3 fundamental into a polyad feature with peaks at 2737(22), 2 835(18), and 2910(12) cm-1, for which only the center frequency has been previously reported. Overall, 9 of the 12 fundamental frequencies of the vinoxy radical are extracted from the IR-cryo-SEVI and ground-state cryo-SEVI spectra, most of which are consistent with previous measurements. However, we provide a new estimate of the ν5 (CH2 scissoring) fundamental frequency at 1395(11) cm-1 and attribute the discrepancy with previously reported values to a Fermi resonance with the 2ν11 overtone (CH2 wagging).

Exploring the Origins of Spectral Signatures of Strong Hydrogen Bonding in Protonated Water Clusters.

The effects of anharmonicity on the spectral features of strong ionic hydrogen bonds are explored through reduced dimensional studies of the couplings between the hydrogen bonding OH and the donor-acceptor OO stretching vibrations in protonated water clusters with 2-4 water molecules. Specifically, this study focuses on how the anharmonicities and couplings in these ions are reflected in the vibrational spectra by exploring the intensities of the transitions to states with excitation in both the OH and the OO stretching vibrations and changes in the frequency of the OO stretching vibration when the OH stretching vibration is excited. These questions are addressed through the application of several approximate treatments that are based on an adiabatic separation of the high-frequency OH and low-frequency OO stretching vibrations as well as low-order expansions of the potential and dipole surfaces. While an adiabatic approximation captures most of the trends found in the spectra and from an analysis of the two-dimensional model, a vibrational Franck-Condon approach fails to capture the intensities of these transitions. Of the terms in the expansion of the dipole moment function, those that are proportional to ΔrOH and ΔrOH2 are found to provide the largest contributions to the calculated intensities of the transitions involving excitation of both the OH and the OO stretches. This leads to the conclusion that the intensities of these transitions encode information about the frequency and anharmonicity of the OH stretching vibration and how they are affected by changes in the OO distance. The anharmonicity of the potential also leads to changes in the OO stretching frequency with excitation of the OH stretching vibration. The direction of this change in frequency encodes additional information about the strength of the ionic hydrogen bond.

Exploring Expansions of the Potential and Dipole Surfaces Used for Vibrational Perturbation Theory.

A scheme for evaluating expansions of the potential and dipole moment surfaces for vibrational perturbation theory is described. The approach is based on numerical differentiation of the Hessian in the coordinates of interest. It is shown that performing these calculations in internal coordinates generates expansions that are transferable among isotopologues of the molecule of interest. Additionally, re-expressing the expansion of the potential in terms of functions of the internal coordinates, for example, cosines of angles or exponential functions of the bond length displacements, provides expansions that can be used for higher-order perturbation theory calculations. The approach is explored and the results are discussed for water, HOD, ammonia, isomers of HNO3, and halogenated methane.

Vibrational Signatures of HNO3 Acidity When Complexed with Microhydrated Alkali Metal Ions, M+·(HNO3)(H2O)n=5 (M = Li, K, Na, Rb, Cs), at 20 K.

The speciation of strong acids like HNO3 under conditions of restricted hydration is an important factor in the rates of chemical reactions at the air-water interface. Here, we explore the trade-offs at play when HNO3 is attached to alkali ions (Li+-Cs+) with four water molecules in their primary hydration shells. This is achieved by analyzing the vibrational spectra of the M+·(HNO3)(H2O)5 clusters cooled to about 20 K in a cryogenic photofragmentation mass spectrometer. The local acidity of the acidic OH group is estimated by the extent of the red shift in its stretching frequency when attached to a single water molecule. The persistence of this local structural motif (HNO3-H2O) in all of these alkali metal clusters enables us to determine the competition between the effect of the direct complexation of the acid with the cation, which acts to enhance acidity, and the role of the water network in the first hydration shell around the ions, which acts to counter (screen) the intrinsic effect of the ion. Analysis of the vibrational features associated with the acid molecule, as well as those of the water network, reveals how cooperative interactions in the microhydration regime conspire to effectively offset the intrinsic enhancement of HNO3 acidity afforded by attachment to the smaller cations.

Fast Near Ab Initio Potential Energy Surfaces Using Machine Learning.

A machine-learning based approach for evaluating potential energies for quantum mechanical studies of properties of the ground and excited vibrational states of small molecules is developed. This approach uses the molecular-orbital-based machine learning (MOB-ML) method to generate electronic energies with the accuracy of CCSD(T) calculations at the same cost as a Hartree-Fock calculation. To further reduce the computational cost of the potential energy evaluations without sacrificing the CCSD(T) level accuracy, GPU-accelerated Neural Network Potential Energy Surfaces (NN-PES) are trained to geometries and energies that are collected from small-scale Diffusion Monte Carlo (DMC) simulations, which are run using energies evaluated using the MOB-ML model. The combined NN+(MOB-ML) approach is used in variational calculations of the ground and low-lying vibrational excited states of water and in DMC calculations of the ground states of water, CH5+, and its deuterated analogues. For both of these molecules, comparisons are made to the results obtained using potentials that were fit to much larger sets of electronic energies than were required to train the MOB-ML models. The NN+(MOB-ML) approach is also used to obtain a potential surface for C2H5+, which is a carbocation with a nonclassical equilibrium structure for which there is currently no available potential surface. This potential is used to explore the CH stretching vibrations, focusing on those of the bridging hydrogen atom. For both CH5+ and C2H5+ the MOB-ML model is trained using geometries that were sampled from an AIMD trajectory, which was run at 350 K. By comparison, the structures sampled in the ground state calculations can have energies that are as much as ten times larger than those used to train the MOB-ML model. For water a higher temperature AIMD trajectory is needed to obtain accurate results due to the smaller thermal energy. A second MOB-ML model for C2H5+ was developed with additional higher energy structures in the training set. The two models are found to provide nearly identical descriptions of the ground state of C2H5+.

A flexible approach to vibrational perturbation theory using sparse matrix methods.

A sparse linear algebra based implementation of Rayleigh-Schrödinger vibrational perturbation theory is presented. This implementation allows for flexibility in the coordinates used to expand the vibrational Hamiltonian as well as the order to which the perturbation theory is performed. It also provides a powerful tool for investigating the origin of spectral intensity and transition frequencies. Specifically, this flexibility allows for the analysis of which terms in the expansions of the Hamiltonian and dipole surface lead to the largest corrections to the energies and transition intensities, and how these conclusions depend on the coordinates used for these expansions. Comparisons of corrections to transition frequencies are reported for the Morse oscillator when the potential is expanded in Δr and Morse coordinates as well as for water, water dimer, and peroxynitrous acid when the molecular Hamiltonians and dipole surfaces are expanded in Cartesian displacement coordinates and in the displacements of the bond-angle-dihedral internal coordinates. Further comparisons of the corrections to the transitions moments are made for H2O and (H2O)2. It is found that while the transition frequencies and intensities are independent of coordinate choice, a good choice of coordinates leads to a cleaner interpretation of the origins of the anharmonicities in these systems.

A wave function correction-based approach to the identification of resonances for vibrational perturbation theory.

An approach for identifying resonances in vibrational perturbation theory calculations is introduced. This approach makes use of the corrections to the wave functions that are obtained from non-degenerate perturbation theory calculations to identify spaces of states that must be treated with degenerate perturbation theory. Pairs of states are considered to be in resonance if the magnitude of expansion coefficients in the corrections to the wave functions in the non-degenerate perturbation theory calculation is greater than a specified threshold, χmax. This approach is applied to calculations of the vibrational spectra of CH4, H2CO, HNO3, and cc-HOONO. The question of how the identified resonances depend on the value of χmax and how the choice of the resonance spaces affects the calculated vibrational spectrum is further explored for H2CO. The approach is also compared to the Martin test [J. M. L. Martin et al., J. Chem. Phys. 103, 2589-2602 (1995)] for calculations of the vibrational spectra of H2CO and cc-HOONO.

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-2 mbar) 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.

Infrared spectroscopic signature of a hydroperoxyalkyl radical (•QOOH).

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH2(CH3)2COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm-1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. The distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

GPU-Accelerated Neural Network Potential Energy Surfaces for Diffusion Monte Carlo.

Diffusion Monte Carlo (DMC) provides a powerful method for understanding the vibrational landscape of molecules that are not well-described by conventional methods. The most computationally demanding step of these calculations is the evaluation of the potential energy. In this work, a general approach is developed in which a neural network potential energy surface is trained by using data generated from a small-scale DMC calculation. Once trained, the neural network can be evaluated by using highly parallelizable calls to a graphics processing unit (GPU). The power of this approach is demonstrated for DMC simulations on H2O, CH5+, and (H2O)2. The need to include permutation symmetry in the neural network potentials is explored and incorporated into the molecular descriptors of CH5+ and (H2O)2. It is shown that the zero-point energies and wave functions obtained by using the neural network potentials are nearly identical to the results obtained when using the potential energy surfaces that were used to train the neural networks at a substantial savings in the computational requirements of the simulations.

Using Diffusion Monte Carlo Wave Functions to Analyze the Vibrational Spectra of H7O3+ and H9O4.

An approach for evaluating spectra from ground state probability amplitudes (GSPA) obtained from diffusion Monte Carlo (DMC) simulations is extended to improve the description of excited state energies and allow for coupling among vibrational excited states. This approach is applied to studies of the protonated water trimer and tetramer, and their deuterated analogs. These ions provide models for solvated hydronium, and analysis of these spectra provides insights into spectral signatures of proton transfer in aqueous environments. In this approach, we obtain a separable set of internal coordinates from the DMC ground state probability amplitude. A basis is then developed from products of the DMC ground state wave function and low-order polynomials in these internal coordinates. This approach provides a compact basis in which the Hamiltonian and dipole moment matrix are evaluated and used to obtain the spectrum. The resulting spectra are in good agreement with experiment and in many cases provide comparable agreement to the results obtained using much larger basis sets. In addition, the compact basis allows for interpretation of the spectral features and how they evolve with cluster size and deuteration.

Coupling of torsion and OH-stretching in tert-butyl hydroperoxide. II. The OH-stretching fundamental and overtone spectra.

The vibrational spectra of gas phase tert-butyl hydroperoxide have been recorded in the OH-stretching fundamental and overtone regions (ΔvOH = 1-5) at room temperature using conventional Fourier transform infrared (ΔvOH = 1-3) and cavity ring-down (ΔvOH = 4-5) spectroscopy. In hydroperoxides, the OH-stretching and COOH torsion vibrations are strongly coupled. The double-well nature of the COOH torsion potential leads to tunneling splitting of the energy levels and, combined with the low frequency of the torsional vibration, results in spectra in the OH-stretching regions with multiple vibrational transitions. In each of the OH-stretching regions, both an OH-stretching and a stretch-torsion combination feature are observed, and we show direct evidence for the tunneling splitting in the OH-stretching fundamental region. We have developed two complementary vibrational models to describe the spectra of the OH-stretching regions, a reaction path model and a reduced dimensional local mode model, both of which describe the features of the vibrational spectra well. We also explore the torsional dependence of the OH-stretching transition dipole moment and show that a Franck-Condon treatment fails to capture the intensity in the region of the stretch-torsion combination features. The accuracy of the Franck-Condon treatment of these features improves with increasing ΔvOH.

Coupling of torsion and OH-stretching in tert-butyl hydroperoxide. I. The cold and warm first OH-stretching overtone spectrum.

The infrared (IR) spectrum of tert-butyl hydroperoxide (TBHP) in the region of the first OH-stretching overtone has been observed under jet-cooled and thermal (300 K, 3 Torr) conditions at ∼7017 cm-1. The jet-cooled spectrum is recorded by IR multiphoton excitation with UV laser-induced fluorescence detection of OH radical products, while direct IR absorption is utilized under thermal conditions. Prior spectroscopic studies of TBHP and other hydroperoxides have shown that the OH-stretch and XOOH (X = H or C) torsion vibrations are strongly coupled, resulting in a double well potential associated with the torsional motion about the OO bond that is different for each of the OH-stretching vibrational states. A low barrier between the wells on the torsional potential results in tunneling split energy levels, which leads to four distinct transitions associated with excitation of the coupled OH-stretch-torsion states. In order to interpret the experimental results, two theoretical models are used that include the OH-stretch-torsion coupling in TBHP. Both methods are utilized to compute the vibrational transitions associated with the coupled OH-stretch-torsion states of TBHP, revealing the underlying transitions that compose the experimentally observed features. A comparison between theory and experiment illustrates the necessity for treatments that include OH-stretch and COOH torsion in order to unravel the spectral features observed in the first OH-stretching overtone region of TBHP.

Guided Diffusion Monte Carlo: A Method for Studying Molecules and Ions That Display Large Amplitude Vibrational Motions.

Diffusion Monte Carlo provides an effective and efficient approach for calculating ground state properties of molecular systems based on potential energy surfaces. The approach has been shown to require increasingly large ensembles when intra- and intermolecular vibrations are weakly coupled. We recently proposed a guided variant of diffusion Monte Carlo to address these challenges for water clusters [Lee, V. G. M.; McCoy, A. B. J. Phys. Chem. A 2019, 123, 8063-8070]. In the present study, we extend this approach and apply it to more strongly bound molecular ions, specifically CH5+ and H+(H2O)n=1-4. For the protonated water systems, we show that the guided DMC approach that was developed for studies of (H2O)n can be used to describe the OH stretches and HOH bends in the solvating water molecules, as well as the free OH stretches in the hydronium core. For the hydrogen bonded OH stretches in the H3O+ core of H+(H2O)n and the CH stretches in CH5+, we develop adaptive guiding functions based on the instantaneous structure of the ion of interest. Using these guiding functions, we demonstrate that we are able to obtain converged zero-point energies and ground state wave functions using ensemble sizes that are as small as 10% the size that is needed to obtain similar accuracy from unguided calculations.