Quantum Features of Anionic Species He*⁻ and He₂*⁻ in Small He(N) Clusters

We present variational calculations on systems containing a few boson helium atoms attached to electronically excited atomic and molecular helium anions He*⁻ and He₂*⁻ and characterize their structures and energetics. Previously reported high-level ab initio results [Huber, S. E.; Mauracher, A. Mol. Phys. 2014, 112, 794] to describe the interactions between excited (metastable) anions and a neutral He atom have been employed. For the case of the atomic species He*⁻, the corresponding interaction with He suggests large anharmonicity effects due to the presence of a deep well of ∼17,500 cm⁻¹ at short distances, together with a more external shallow secondary well of ∼4 cm⁻¹, both supporting bound levels. Moreover, when a sum of pairwise interactions is assumed to describe the full PES corresponding to the presence of several neutral He atoms, geometrical constraints already predict the complete solvation of the anionic impurity by six helium atoms, giving rise to a bipyramidal structure. In turn, for the anisotropic weak interaction He-He₂*⁻, where the anionic dimer is considered as a rigid rotor, the obtained structures show the tendency of the helium atoms to pack themselves together and largely far away from the dopant, thereby confirming the heliophobic character of He₂*⁻.

Reactive scattering calculations for (87)Rb+(87)RbHe→Rb2((3)Σ(u)(+),v)+He from ultralow to intermediate energies

We investigate atom-diatom reactive collisions, as a preliminary step,in order to assess the possibility of forming Rb(2) molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets [corrected]. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.

Computational investigations of the thermodynamic properties of size-selected water and Ar-water clusters: high-pressure transitions

Classical parallel-tempering Monte Carlo simulations in the isothermal-isobaric ensemble were carried out for the (H2O)20 and Ar(H2O)20 clusters, over a wide range of temperatures (30-1000 K) and pressures (3 kPa-10 GPa) in order to study their thermodynamic properties and structural changes. The TIP4P/ice water model is employed for the water-water interactions, while both semiempirical and ab initio-based potentials are used to model the interaction between the rare-gas atoms and the water molecules. Temperature-pressure phase diagrams for these cluster systems were constructed by employing a two-dimensional multiple-histogram method. Structural changes were detected by analyzing the heat capacity landscape and the Pearson correlation coefficient profile for the interaction energy and volume. Those at high pressure correspond to solid-to-solid transitions and are found to be related to clathrate-like cages around the Ar atom. It is also shown that the formation and thermodynamic stability of such structures are determined by the intermolecular interaction between the rare-gas atoms and the host water molecules.

Vibrational dynamics of the H5(+) and its isotopologues from multiconfiguration time-dependent Hartree calculations

Full-dimensional multiconfiguration time-dependent Hartree (MCTDH) computations are reported for the vibrational states of the H(5)(+) and its H(4)D(+), H(3)D(2)(+), H(2)D(3)(+), HD(4)(+), D(5)(+) isotopologues employing two recent analytical potential energy surfaces of Xie et al. [J. Chem. Phys. 122, 224307 (2005)] and Aguado et al. [J. Chem. Phys. 133, 024306 (2010)]. The potential energy operators are constructed using the n-mode representation adapted to a four-combined mode cluster expansion, including up to seven-dimensional grids, chosen adequately to take advantage in representing the MCTDH wavefunction. An error analysis is performed to quantify the convergence of the potential expansion to reproduce the reference surfaces at the energies of interest. An extensive analysis of the vibrational ground state properties of these isotopes and comparison with the reference diffusion Monte Carlo results by Acioli et al. [J. Chem. Phys. 128, 104318 (2008)] are presented. It is found that these systems are highly delocalized, interconverting between equivalent minima through rotation and internal proton transfer motions even at their vibrational ground state. Isotopic substitution affects the zero-point energy and structure, showing preference in the arrangements of the H and D within the mixed clusters, and the most stable conformers of each isotopomer are the ones with the H in the central position. Vibrational excited states are also computed and by comparing the energies and structures predicted from the two surfaces, the effect of the potential topology on them is discussed.

Quantum-dynamics study of the H5+ cluster: full dimensional benchmark results on its vibrational states

A full-dimensional quantum dynamics study is carried out for the highly fluxional H(5)(+) cation on a recent reference potential energy surface by using the multi configuration time-dependent Hartree method. With five equivalent light atoms and shallow barriers between various low-lying stationary points on the surface, the spectroscopic characterization of H(5)(+) represents a huge challenge for accurate quantum dynamics simulations. The present calculation is the first such a study on this cation, which together with its isotope analogies are of primary importance in the interstellar chemistry. The vibrational ground state properties and several vibrationally excited states corresponding to low vibrational frequency motions, not yet directly observable by the experiment, are presented and analyzed.

Theoretical investigation of the He4Br2 conformers

Full dimensional quantum dynamics calculations of the three lowest isomers of the He(4)Br(2) van der Waals molecule in its ground electronic state are reported. The calculations are performed using the multiconfiguration time-dependent Hartree (MCTDH) method and a realistic potential form that includes the sum of three body ab initio coupled-cluster single double triple [CCSD(T)] He-Br(2) interactions plus the He-He and Br-Br interactions. This potential exhibits several multiple minima, with the three lowest ones lying very close in energy, just within 2 cm(-1). Such small differences are also found in the calculated binding energies of the three most stable conformers, indicating the floppiness of the system and, thus, the need of accurate potential forms and quantum full dynamics methods to treat this kind of complexes. The 12 dimensional results reported in this work present benchmark data and, thus, can serve to evaluate approximate methods aiming to describe higher order rare gas-dihalogen (N > 4) complexes. A comparison with previous studies using different potential forms and approaches to the energetics for the He(4)Br(2) cluster is also presented.

Spin-polarized Rb2 interacting with bosonic He atoms: potential energy surface and quantum structures of small clusters

A new full-dimension potential energy surface of the three-body He-Rb₂(³Σ(u)(+)) complex and a quantum study of small (⁴He)(N)-Rb₂(³Σ(u)(+)) clusters, 1 ≤ N ≤ 4, are presented. We have accurately fitted the ab initio points of the interaction to an analytical form and addressed the dopant's vibration, which is found to be negligible. A Variational approach and a Diffusion Monte Carlo technique have been applied to yield energy and geometric properties of the selected species. Our quantum structure calculations show a transition in the arrangements of the helium atoms from N = 2, where they tend to be separated across the diatomic bond, to N = 4, in which a closer packing of the rare gas particles is reached, guided by the dominance of the He-He potential over the weaker interaction of the latter adatoms with the doping dimer. The deepest well of the He-Rb₂ interaction is placed at the T-shape configuration, a feature which causes the dopant to be located as parallel to the helium "minidroplet". Our results are shown to agree with previous findings on this and on similar systems.

Weakly bound finite systems: (⁴He)N-Rb₂(³Σu), clustering structures from a quantum Monte Carlo approach

We report here ((4)He)(N)-Rb(2)((3)Σ(u)) complexes, 2 ≤ N ≤ 20, analysed through a quantum diffusion Monte Carlo stochastic approach. The calculations show that the spin stretched dimer molecule is bound outside the pure He sub-complex, due to the stronger He-He potential as compared with the He-Rb(2) interaction, while the rare gas atom moiety presents, in turn, a shell-like structure with ten He adatoms completing the first shell. Our results agree with previous findings on this and similarly weakly interacting systems.

Theoretical characterization of intermolecular vibrational states through the multi-configuration time dependent Hartree approach: the He2,3ICl clusters

Benchmark, full-dimensional calculations on the ground and excited vibrational states for the tetra-, and penta-atomic weakly bound He(2,3)ICl complexes are reported. The representation of the potential energy surfaces includes three-body HeICl potentials parameterized to coupled-cluster singles, doubles, and perturbative triples ab initio data. These terms are important in accurately describing the interactions of such highly floppy systems. The corresponding 6D/9D computations are performed with the multi-configuration time dependent Hartree method, using natural potential fits, and a mode combination scheme to optimize the computational effort in the improved relaxation calculations. For these complexes several low-lying vibrational states are computed, and their binding energies and radial/angular probability density distributions are obtained. We found various isomers which are assigned to different structural models related with combinations of the triatomic isomers, like linear, T-shaped, and antilinear ones. Comparison of these results with recent experimental data is presented, and the quantitative deviations found with respect to the experiment are discussed.

Frequency domain description of Kohlrausch response through a pair of Havriliak-Negami-type functions: an analysis of functional proximity

An approximation to the Fourier transform (FT) of the Kohlrausch function (stretched exponential) with shape parameter 0 < β ≤ 1 is presented by using Havriliak-Negami-like functions. Mathematical expressions to fit their parameters α, γ, and τ, as functions of β (0 < β ≤ 1 and 1 < β < 2) are given, which allows a quick identification in the frequency domain of the corresponding shape factor β. Reconstruction via fast Fourier transform of frequency approximants to time domain are shown as good substitutes in short times though biased in long ones (increasing discrepancies as β → 1). The method is proposed as a template to commute time and frequency domains when analyzing complex data. Such a strategy facilitates intensive algorithmic search of parameters while adjusting the data of one or several Kohlrausch-Williams-Watts relaxations.

Full-dimensional multi configuration time dependent Hartree calculations of the ground and vibrationally excited states of He2,3Br2 clusters

Quantum dynamics calculations are reported for the tetra-, and penta-atomic van der Waals He(N)Br(2) complexes using the multiconfiguration time-dependent Hartree (MCTDH) method. The computations are carried out in satellite coordinates, and the kinetic energy operator in this set of coordinates is given. A scheme for the representation of the potential energy surface based on the sum of the three-body HeBr(2) interactions at CSSD(T) level plus the He-He interaction is employed. The potential surfaces show multiple close lying minima, and a quantum description of such highly floppy multiminima systems is presented. Benchmark, full-dimensional converged results on ground vibrational/zero-point energies are reported and compared with recent experimental data available for all these complexes, as well as with previous variational quantum calculations for the smaller HeBr(2) and He(2)Br(2) complexes on the same surface. Some low-lying vibrationally excited eigenstates are also computed by block improved relaxation calculations. The binding energies and the corresponding vibrationally averaged structures are determined for different conformers of these complexes. Their relative stability is discussed, and contributes to evaluate the importance of the multiple-minima topology of the underlying potential surface.

Quantum features of a barely bound molecular dopant: Cs2(3Σu) in bosonic helium droplets of variable size

We present in this work the study of small (4)He(N)-Cs(2)((3)Σ(u)) aggregates (2 ≤ N ≤ 30) through combined variational, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) calculations. The full surface is modeled as an addition of He-Cs(2) interactions and He-He potentials. Given the negligible strength and large range of the He-Cs(2) interaction as compared with the one for He-He, a propensity of the helium atoms to pack themselves together, leaving outside the molecular dopant is to be expected. DMC calculations determine the onset of helium gathering at N = 3. To analyze energetic and structural properties as a function of N, PIMC calculations with no bosonic exchange, i.e., Boltzmann statistics, at low temperatures are carried out. At T = 0.1 K, although acceptable one-particle He-Cs(2) distributions are obtained, two-particle He-He distributions are not well described, indicating that the proper symmetry should be taken into account. PIMC distributions at T = 1 K already compare well with DMC ones and show minor exchange effects, although binding energies are still far from having converged in terms of the number of quantum beads. As N increases, the He-He PIMC pair correlation function shows a clear tendency to coincide with the experimental boson-liquid helium one at that temperature. It supports the picture of a helium droplet which carries the molecular impurity on its surface, as found earlier for other triplet dimers.

Internal proton transfer and H2 rotations in the H5(+) cluster: a marked influence on its thermal equilibrium state

Classical and path integral Monte Carlo (CMC, PIMC) "on the fly" calculations are carried out to investigate anharmonic quantum effects on the thermal equilibrium structure of the H5(+) cluster. The idea to follow in our computations is based on using a combination of the above-mentioned nuclear classical and quantum statistical methods, and first-principles density functional (DFT) electronic structure calculations. The interaction energies are computed within the DFT framework using the B3(H) hybrid functional, specially designed for hydrogen-only systems. The global minimum of the potential is predicted to be a nonplanar configuration of C(2v) symmetry, while the next three low-lying stationary points on the surface correspond to extremely low-energy barriers for the internal proton transfer and to the rotation of the H2 molecules, around the C2 axis of H5(+), connecting the symmetric C(2v) minima in the planar and nonplanar orientations. On the basis of full-dimensional converged PIMC calculations, results on the quantum vibrational zero-point energy (ZPE) and state of H5(+) are reported at a low temperature of 10 K, and the influence of the above-mentioned topological features of the surface on its probability distributions is clearly demonstrated.

Toward a realistic density functional theory potential energy surface for the H5+ cluster

The potential energy surface of H(5)(+) is characterized using density functional theory. The hypersurface is evaluated at selected configurations employing different functionals, and compared with results obtained from ab initio CCSD(T) calculations. The lowest ten stationary points (minima and saddle-points) on the surface are located, and the features of the short-, intermediate-, and long-range intermolecular interactions are also investigated. A detailed analysis of the surface's topology, and comparisons with extensive CCSD(T) results, as well as a recent ab initio analytical surface, shows that density functional theory calculations using the B3(H) functional represent very well all aspects studied on the H(5)(+) potential. These include the tiny energy difference between the minimum at 1-C(2v) configuration and the 2-D(2d) one corresponding to the transition state for the proton transfer between the two equivalent C(2v) minima, and also the correct asymptotic behavior of the long-range interactions. The calculated binding energy and dissociation enthalpies compare very well with previous benchmark coupled-cluster ab initio data, and with experimental data available. Based on these results the use of such approach to perform first-principles molecular dynamics simulations could provide reliable information regarding the dynamics of protonated hydrogen clusters.

Exact, Born-Oppenheimer, and quantum-chemistry-like calculations in helium clusters doped with light molecules: The He2N2(X) system

Helium clusters doped with diatomic molecules, He(N)-BC, have been recently studied by means of a quantum-chemistry-like approach. The model treats He atoms as "electrons" and dopants as "nuclei" in standard electronic structure calculations. Due to the large mass difference between He atoms and electrons, and to the replacement of Coulomb interactions by intermolecular potentials, it is worth assessing up to what extent are the approximations involved in this model, i.e., decoupling of the BC rotation from the He-atom orbital angular momenta and Born-Oppenheimer separation of the BC stretch versus the He motions, accurate enough. These issues have been previously tackled elsewhere for the (4)He(2)-Br(2)(X) system, which contains a heavy dopant [Roncero et al., Int. J. Quantum Chem. 107, 2756 (2007)]. Here, we consider a similar cluster but with a much lighter dopant such as N(2)(X). Although the model does not provide the correct energy levels for the cluster, positions and intensities of the main detectable lines of the vibrotational Raman spectrum at low temperature are accurately reproduced.

Ab initiovibrational predissociation dynamics of He–I2(B) complex

Three-dimensional quantum mechanical calculations on the vibrational predissociation dynamics of HeI2 B state complex are performed using a potential energy surface accurately fitted to unrestricted open-shell coupled cluster ab initio data, further enabling extrapolation for large I2 bond lengths. A Lanczos iterative method with an optimized complex absorbing potential is used to determine energies and lifetimes of the vibrationally predissociating He,I2(B,v') complex for v'<or=26 of I2 vibrational excitations. The calculated predissociating state energies agree with recent experimental results within 0.5 cm(-1). This excellent agreement is remarkable since no adjustment was made with respect to the experiments. The present ab initio approach, however, shows its limitations in the fact that the computed lifetimes that are highly sensitive to subtle details of the potential energy surfaces such as anisotropy, are a factor of 1.5 larger than the available experimental data.

Ab initio potential energy surface and spectrum of the B(Π3) state of the HeI2 complex

The three-dimensional interaction potential for I2(B 3Pi0u+)+He is computed using accurate ab initio methods and a large basis set. Scalar relativistic effects are accounted for by large-core relativistic pseudopotentials for the iodine atoms. Using multireference configuration interaction calculations with subsequent treatment of spin-orbit coupling, it is shown for linear and perpendicular structures of the complex that the interaction potential for I2(B 3Pi0u+)+He is very well approximated by the average of the 3A' and 3A" interaction potentials obtained without spin-orbit coupling. The three-dimensional 3A' and 3A" interaction potentials are computed at the unrestricted open-shell coupled-cluster level of theory using large basis sets. Bound state calculations based on the averaged surface are carried out and binding energies, vibrationally averaged structures, and frequencies are determined. These results are found to be in excellent accord with recent experimental measurements from laser-induced fluorescence and action spectra of HeI2. Furthermore, in combination with a recent X-state potential, the spectral blueshift is obtained and compared with available experimental values.

A full-configuration interaction “nuclear orbital” method to study doped HeN3 clusters (N⩽4)

An efficient full configuration interaction (FCI) treatment, based on the Jacobi-Davidson algorithm, is developed in order to study small doped (3)He(N) clusters. The state of each He atom in a given cluster is described by a set of wave-functions which by extention of the quantum-chemistry notation are caller here "nuclear orbitals". The FCI treatment is applied to the calculation of binding energies and helium natural orbitals of (3)He(N)...Br(2)(X) complexes. In agreement with our previous calculations using a Hartree-Fock approach [Phys. Rev. Lett. 93, 053401 (2004)], in which the He-He interaction is modified at small distances to account for short-range correlation effects, the lowest-energy states of each multiplet are found to be very close in energy. The natural orbital analysis, in turn, indicates the adequacy of the "nuclear orbital" approach in these systems.

A theoretical study of He2ICl van der Waals cluster

The structure, energetics, and dynamics of He2ICl complex in its ground state are studied by means of ab initio electronic structure and quantum-mechanical calculations. Interaction energies for selected He2ICl configurations are calculated at the coupled-cluster [CCSD(T)] level of theory using a large-core pseudopotential for the I atom and the aug-cc-pVTZ and aug-cc-pV5Z basis sets for the Cl and He atoms, respectively. The surface is characterized around its lower five minima and the minimum energy pathways through them. The global minimum of the potential corresponds to a "police-nightstick (1)" configuration, the second one to a linear, the next one to tetrahedral configuration, and the following two to "bifork" and "police-nightstick (2)" structures, with well depths of -99.12, -97.42, -88.32, -85.84, and -78.54 cm(-1), respectively. An analytical form based on the sum of the three-body parametrized HeICl interactions plus the He-He interaction is found to represent very well the tetra-atomic CSSD(T) results. The present potential expression is employed to perform variational five-dimensional quantum-mechanical calculations to study the vibrational bound states of the van der Waals He2ICl complex. Results for total angular momentum J = 0 provide the binding energy D0 and the corresponding vibrationally averaged structure for different isomers of the cluster. Comparison of these results with recent experimental observations further justifies the potential used in this work.

Vibrational quenching of CO2(010) by collisions with O(P3) at thermal energies: A quantum-mechanical study

The CO(2)(010)-O((3)P) vibrational energy transfer (VET) efficiency is a key input to aeronomical models of the energy budget of the upper atmospheres of Earth, Venus, and Mars. This work addresses the physical mechanisms responsible for the high efficiency of the VET process at the thermal energies existing in the terrestrial upper atmosphere (150 K</=T</=550 K). We present a quantum-mechanical study of the process within a reduced-dimensionality approach. In this model, all the particles remain along a plane and the O((3)P) atom collides along the C(2v) symmetry axis of CO(2), which can present bending oscillations around the linear arrangement, while the stretching C-O coordinates are kept fixed at their equilibrium values. Two kinds of scattering calculations are performed on high-quality ab initio potential energy surfaces (PESs). In the first approach, the calculations are carried out separately for each one of the three PESs correlating to O((3)P). In the second approach, nonadiabatic effects induced by spin-orbit couplings (SOC) are also accounted for. The results presented here provide an explanation to some of the questions raised by the experiments and aeronomical observations. At thermal energies, nonadiabatic transitions induced by SOC play a key role in causing large VET efficiencies, the process being highly sensitive to the initial fine-structure level of oxygen. At higher energies, the two above-mentioned approaches tend to coincide towards an impulsive Landau-Teller mechanism of the vibrational to translational (V-T) energy transfer.

Rovibrational Structures in Floppy Triatomics: Distributed Gaussian Functions Treatment for the Ne2H- System

The full sequence of the bound states for a very floppy triatomic complex, Ne2H- in its ground electronic state, are initially computed for the rotationless situation and employing a variational approach that expands the total nuclear wave function over a large set of symmetry-adapted, distributed Gaussian functions and employs accurate atom-atom potential energy data. The results are tested for numerical convergence, compared with the behavior of both its diatomic fragments, Ne2 and NeH-, and further compared with the results for the Ne3 case. The computational analysis is extended to the production of the rotational constants for the very nonclassical ground state vibrational configuration by making use of the previous findings. The method is shown to provide us with several illuminating details on the nanoscopic internal dynamics of this very weakly bound quantum aggregate.

Br2(X) Microsolvation in Helium Clusters: Effect of the Interaction on the Quantum Solvent Density Distribution

The Born-Oppenheimer potential energy surface for the Br2(X) molecule interacting with a varying number of 4He bosons is constructed following two different schemes which employ either a full ab initio evaluation of the Br2-He interaction forces or an estimate of the latter through an empirical model. Both descriptions are employed by carrying out diffusion Monte Carlo (DMC) calculations of the ground-state energies and quantum wavefunctions for Br2-(He)n clusters with n up to 24. The results clearly indicate, for both interactions, the occurrence of the full solvation of the molecular dopant within the quantum bosonic "solvent" but also show differences between the two models in terms of the expected density distributions of the surrounding particles within the shorter-range region that makes up the clusters with smaller n values. Our calculations also show that such differences become insignificant for the larger 4He clusters surrounding the Br2 molecule, where density profiles and bulk behaviour are chiefly driven by the solvent structure, once n values reach the region of 15-20 adatoms.

Bound-state energies in argon trimers via a variational expansion: The effects from many-body corrections

In this paper we study the bound-state energies and geometries of Ar(3) for J=0, using the distributed Gaussian functions method that provides a configurational description of the different structures contributing to these states. Atom-atom potentials are employed and three-body long-range effects are also included in the computational treatment by adding to the sum of potentials the Axilrod-Teller triple-dipole correction for the whole rotationless energy spectrum. An estimate of the total number of bound states for the Ar trimer is given. With respect to previous calculations, limited to the lower-lying states, our results show slightly larger nonadditive effects and are further able to predict the full range of the bound spectrum. Changes on the geometries of a large part of the vibrationally excited states of Ar(3) when the Axilrod-Teller term is included in the molecular potential are found by the present study.

The binding of He4 and He3 to a hydrogen molecule: A computational study for pH2 and oH2

One bound state for the very weakly interacting complex of H(2) with one He atom is computed using a recently obtained, very accurate potential energy surface [A. I. Boothroyd and P. G. Martin, J. Chem. Phys. 119, 3187 (2003)] and employing a discrete variable representation expansion. The spatial features of this very diffused, "halo" state are analyzed and discussed for various isotopic variants of the title system.

A complete configurational study for the bound states of Ne trimers

The structural properties and the energetics of the ground and the excited bound states of Ne(3) for zero total angular momentum are examined using different modelings for the two-body interactions. We employ a method consisting of a variational approach with a distributed Gaussian functions (DGF) basis set expansion. We discuss at length the advantages and possible limitations of such an approach, comparing it to other methods which have been applied in the literature to the same system. The DGF method turns out to be very accurate in giving us the bound states energetics and also provides in a natural way a convincing pictorial description of all the states, including those with dominant linear configurations. Additional bound states are found for the Ne(3) system with respect to those indicated in previous works and we suggest a "stabilization" procedure that can be used to assess the truly bound nature of a state. Some considerations on the relative reliability of the examined two-body interactions are also reported.

Ab initio calculations, potential representation and vibrational dynamics of He2Br2 van der Waals complex

An intermolecular potential energy surface for He(2)Br(2) complex in the ground state is calculated at the levels of fourth-order (MP4) Moller-Plesset and coupled-cluster [CCSD(T)] approximations, using large-core pseudopotential for Br atoms and the aug-cc-pV5Z basis set for He. The surface is characterized by three minima and the minimum energy pathways through them. The global minimum corresponds to a linear He-Br(2)-He configuration, while the two other ones to "police-nightstick" and tetrahedral structures. The corresponding well depths are -90.39/-89.18, -81.23/-80.78 and -74.40/-74.02 cm(-1), respectively, at MP4/CCSD(T) levels of theory. It is found that results obtained by summing three-body parametrized HeBr(2) interactions and the He-He interaction are in very good accord with the corresponding MP4/CSSD(T) configuration energies of the He(2)Br(2). Variational calculations using a sum of three-body interactions are presented to study the bound states of the vdW He(2)Br(2) complex. The binding energy D(0) and the corresponding vibrationally averaged structure are determined for different isomers of the cluster and their comparison with the available experimental data is discussed.

Raman spectra of (He)N-Br2(X) clusters: The role of boson/fermion statistics in a quantum solvent

The aim of this paper is to elucidate the role played by the bosonic/fermionic character of N He atoms solvating a Br2(X) molecule. To this end, an adiabatic model in the molecular stretching coordinate is assumed and the ground energy levels of the complexes are searched by means of Hartree (or Hartree-Fock) Quantum Chemistry calculations for 4He (or 3He) solvent atoms. Simulations of vib-rotational Raman spectra point at the spin multiplicity as the main feature responsible for the drastic difference in the rotational structures of molecules embedded in boson or fermion helium drops as already observed by the experiments of Grebenev et al. [S. Grebenev, J. P. Toennies, and A. F. Vilesov, Science 279 (1998) 2083].

Role of boson-fermion statistics on the Raman spectra of Br2(X) in helium clusters

The role played by the bosonic or fermionic character of He atoms surrounding a Br2(X) molecule is analyzed through vibrotational Raman spectra simulations. Quantum chemistry-type calculations reveal the spin multiplicity to be chiefly responsible for the drastic difference observed by Grebenev et al. [Science 279, 2083 (1998)]] in the rotational structure of molecules embedded in helium droplets.

The open-shell interaction of He with the B 3Πu(0+) state of Br2: Anab initiostudy and its comparison with a diatomics-in-molecule perturbation model

The interaction of He with Br2 in electronically excited B 3Piu state is investigated using spin-unrestricted single and double coupled-cluster approach with noniterative perturbative treatment of triple excitations. Internal electrons of the Br atom are described by effective core pseudopotentials. The validity of this approach is analyzed by comparing the lowest 2Sigma+ and 2Pi electronic states of the HeBr molecule with those obtained in all electron calculations [J. Chem. Phys. 115, 10438 (2001)]. In this context, we examine the performance of different basis sets and saturation with bond functions. The comparison of theoretical blue-shifts with the experiment provides confidence about the present ab initio calculations. In addition, He-Br results of ab initio calculations at the same level are used to obtain approximate He-Br2 (3Piu) interactions in the framework of the diatomics-in-molecule first order perturbation theory (IDIM-PT1) [J. Chem. Phys. 104, 9913 (1996)]. Overall, the IDIM-PT1 model results show a good agreement with the ab initio ones, being the main difference the sensitivity to the elongation of the Br-Br bond.