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

Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems

The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.

A nuclear configuration interaction approach to study nuclear spin effects: an application to ortho- and para-3 He2 @C60

We introduce a non-orthogonal configuration interaction approach to investigate nuclear quantum effects on energies and densities of confined fermionic nuclei. The Hamiltonian employed draws parallels between confined systems and many-electron atoms, where effective non-Coulombic potentials represent the interactions of the trapped particles. One advantage of this method is its generality, as it offers the potential to study the nuclear quantum effects of various confined species affected by effective isotropic or anisotropic potentials. As a first application, we analyze the quantum states of two 3 He atoms encapsulated in C60 . At the Hartree-Fock level, we observe the breaking of spin and spatial symmetries. To ensure wavefunctions with the correct symmetries, we mix the broken-symmetry Hartree-Fock states within the non-orthogonal configuration interaction expansion. Our proposed approach predicts singly and triply degenerate ground states for the singlet (para-3 He2 @C60 ) and triplet (ortho-3 He2 @C60 ) nuclear spin configurations, respectively. The ortho-3 He2 @C60 ground state is 5.69 cm-1 higher in energy than the para-3 He2 @C60 ground state. The nuclear densities obtained for these states exhibit the icosahedral symmetry of the C60 embedding potential. Importantly, our calculated energies for the lowest 85 states are in close agreement with perturbation theory results based on a harmonic oscillator plus rigid rotor model of 3 He2 @C60 .

A nuclear spin and spatial symmetry-adapted full quantum method for light particles inside carbon nanotubes: clusters of 3He, 4He, and para-H2

A new nuclear spin and spatial symmetry-adapted full quantum method for light fermionic and bosonic particles under cylindrical carbon nanotube confinement.

Quantum effects on the stability of the He5 I2 van der Waals conformers

We present 15-dimensional quantum multiconfiguration time-dependent Hartree calculations of the vibrational levels of the He₅ I₂ van der Waals (vdW) complex employing an ab initio-based potential energy surface (PES). The energies and spatial features of such bound structures are analyzed, providing predictions on the structures and relative stabilities of its three lowest isomers. We found that the most stable isomer corresponds to all five He atoms encircling the I₂ molecule, indicating that in this case the anharmonic quantum effects do not stabilize the isomers involving a He atom in a linear configuration as reported previously for the smaller He_N I₂ systems. Such finding provides information on the overall structuring of the finite-size-solvent systems, highlighting the intriguing interplay between weak intermolecular interactions and quantum effects. © 2019 Wiley Periodicals, Inc.

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₂*⁻.

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.

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.

Ca impurity in small mixed (4)He-(3)He clusters

The structure of small mixed helium clusters doped with one calcium atom has been determined within the diffusion Monte Carlo framework. The results show that the calcium atom sits at the (4)He-(3)He interface. This is in agreement with previous studies, both experimental and theoretical, performed for large clusters. A comparison between the results obtained for the largest cluster we have considered for each isotope shows a clear tendency of the Ca atom to reside in a deep dimple at the surface of the cluster for (4)He clusters, and to become fully solvated for (3)He clusters. We have calculated the absorption spectrum of Ca around the 4s4p<--4s(2) transition and have found that it is blueshifted from that of the free-atom transition by an amount that depends on the size and composition of the cluster.

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.

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.

Solvent configuration-interaction calculations of intermolecular states in molecule-(atom)N clusters: application to Br2-4HeN

We describe variational calculations of J=0 intermolecular states in Br(2)-(4)He(N) clusters. The method employed is analogous to configuration-interaction calculations in electronic-structure work and relies on the ability to express the intermolecular Hamiltonian H(v) as a sum of one- and two-body terms. A basis set is built up from solutions to the Schrödinger equation in which only the one-body terms of H(v) are included. These configurations are products of N=1 eigenstates. The matrix of H(v) in a symmetry-adapted configuration basis is then computed, the two-body terms of H(v) serving to couple different configurations. This computation involves integrals of dimension five or less. Filter diagonalization is then used to obtain energies and eigenfunctions within a selected energy range. Results on clusters having N=2-5 are reported.

Isotopic replacement in ionic systems: The He2+4+He3→He3He+4+He4 reaction

Full quantum dynamics calculations have been carried out for the ionic reaction (4)He(2) (+)+(3)He and state-to-state reactive probabilities have been obtained using both time-dependent and time-independent approaches. An accurate ab initio potential-energy surface has been employed for the present quantum dynamics and the two sets of results are shown to be in agreement with each other. The results for zero total angular momentum suggest a marked presence of atom exchange (isotopic replacement) reaction with probabilities as high as 60%. The reaction probabilities are only weakly dependent on the initial vibrational state of the reactants, while they are slightly more sensitive to the degree of rotational excitation. A brief discussion of the results for selected higher total angular momentum values is also presented, while the l-shifting approximation [S. K. Gray et al., Phys. Chem. Chem. Phys. 1, 1141 (1999)] has been used to provide estimates of the total reaction rates for the title process. Such rates are found to be large enough to possibly become experimentally accessible.

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.