Source author record

Piotr S. Żuchowski

Piotr S. Żuchowski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2020arXiv

Benchmarking the accuracy of seniority-zero wavefunction methods for non-covalent interactions

In this paper, we scrutinize the ability of seniority-zero wavefunction-based methods to model different types of non-covalent interactions, such as hydrogen bonds, dispersion, and mixed non-covalent interactions as well as prototypical model systems with various contributions of dynamic and static electron correlation effects. Specifically, we focus on the pair Coupled-Cluster Doubles (pCCD) ansatz combined with two different flavours of dynamic energy corrections, (i) based on a perturbation theory correction and (ii) on a linearized coupled-cluster ansatz on top of pCCD. We benchmark these approaches against the A24 data set [Řez{á}{č} and Hobza, J.~Chem.~Theory~Comput., 9, 2151-2155 (2013)] extrapolated to the basis set limit and some model non-covalent complexes that feature covalent bond breaking. By dissecting different types of interactions in the A24 data set within the Symmetry-Adapted Perturbation Theory (SAPT) framework, we demonstrate that pCCD can be classified as a dispersion-free method. Furthermore, we found that both flavours of post-pCCD approaches represent encouraging and computationally more efficient alternatives to standard electronic structure methods to model weakly-bound systems, resulting in small statistical errors. Finally, a linearized coupled cluster correction on top of pCCD proved to be most reliable for the majority of investigated systems, featuring smaller non-parallelity errors compared to perturbation-theory-based approaches.

preprint2019arXiv

Nonrigidity effects -- a missing puzzle piece in the description of low-energy anisotropic molecular collisions

Cold collisions serve as a very sensitive probe of the interaction potential. In the recent study of Klein et al. (Nature Phys. 13, 35-38 (2017)) the one-parameter scaling of the interaction potential was necessary to obtain agreement between theoretical and observed patterns of the orbiting resonances for excited metastable helium atoms colliding with hydrogen molecules. Here we show that the effect of nonrigidity of the H$_2$ molecule on the resonant structure, absent in the previous study, is critical to predict correct positions of the resonances in that case. We have complemented the theoretical description of the interaction potential and revised reaction rate coefficients by proper inclusion of the flexibility of the molecule. The calculated reaction rate coefficients are in remarkable agreement with the experimental data without empirical adjustment of the interaction potential. We have shown that even state-of-the-art calculations of the interaction energy cannot ensure agreement with the experiment if such an important physical effect as flexibility of the interacting molecule is neglected. Our findings about the significance of the nonrigidity effects can be especially crucial in cold chemistry, where the quantum nature of molecules is pronounced.

preprint2019arXiv

Product-state distribution after isotopic substitution in ultracold atom-molecule collisions

We show that products of the isotopic substitution reactions in experimentally accessible molecules such as NaK, RbCs, and SrF are cold according to their translational energy below hundreds of mK. For these chemical reactions, molecular products may occupy only the lowest rotational states. We also discuss the possibility of controlling the chemical reactions by the electric field in ultracold mixtures of molecules and atoms with low kinetic energy release, where one of the constituent atoms of colliding molecule is replaced by its isotope. This letter opens new avenues in investigating the branching ratios of chemical reactions in ultracold conditions.

preprint2019arXiv

Quantum chaos in Feshbach resonances of the ErYb system

We investigate ultracold magnetic-field-assisted collisions in the so far unexplored ErYb system. The nonsphericity of the Er atom leads to weakly anisotropic interactions that provide the mechanism for Feshbach resonances to emerge. The resonances are moderately sparsely distributed with a density of $0.1\,{\rm G}^{-1}-0.3\,{\rm G}^{-1}$ and exhibit chaotic statistics characterized by a Brody parameter $η\approx 0.5-0.7$. The chaotic behaviour of Feshbach resonances is accompanied by strong mixing of magnetic and rotational quantum numbers in near-threshold bound states. We predict the existence of broad resonances at fields $<300\,{\rm G}$ that may be useful for the precise control of scattering properties and magnetoassociation of ErYb molecules. The high number of bosonic Er-Yb isotopic combinations gives many opportunities for mass scaling of interactions. Uniquely, two isotopic combinations have nearly identical reduced masses (differing by less than $10^{-5}$ relative) that we expect to have strikingly similar Feshbach resonance spectra, which would make it possible to experimentally measure their sensitivity to hypothetical variations of proton-to-electron mass ratio.

preprint2015arXiv

Ab initio interaction potentials and scattering lengths for ultracold mixtures of metastable helium and alkali-metal atoms

We have obtained accurate ab initio quartet potentials for the diatomic metastable triplet helium + alkali-metal (Li, Na, K, Rb) systems, using all-electron restricted open-shell coupled cluster singles and doubles with noniterative triples corrections [CCSD(T)] calculations and accurate calculations of the long-range $C_6$ coefficients. These potentials provide accurate ab initio quartet scattering lengths, which for these many-electron systems is possible, because of the small reduced masses and shallow potentials that results in a small amount of bound states. Our results are relevant for ultracold metastable triplet helium + alkali-metal mixture experiments.

preprint2015arXiv

Magnetically tunable Feshbach resonances in Li+Er

We explore the magnetic Feshbach spectra of ultracold ground-state Li+Er systems. Our calculations predict many tunable resonances at fields below 1000 G that could be stably tuned in ultracold experiments. We show that Li+Er spectra are much less congested than those of systems involving heavier highly-magnetic atoms and exhibit non-chaotic properties. These features would facilitate identifying and addressing individual resonances. We derive a simple model for the mass-scaling shifting of low-field resonances that may simplify designing experiments with different Er bosonic isotopes. Our work establishes Li+Er as very promising systems for quantum simulation, precision measurements and the formation of polar paramagnetic molecules.

preprint2014arXiv

Ground and excited state properties of the polar and paramagnetic RbSr molecule: a comparative study

This paper deals with the electronic structure of RbSr, a molecule possessing both a permanent magnetic and electric dipole moment in its own frame allowing its manipulation with external fields. Two complementary ab-initio approaches are used for the ground and lowest excited states: first, an approach relying on optimized effective core potentials with core polarization potentials based on a full configuration interaction involving three valence electrons, and second, an approach using a small-size effective core potential with 19 correlated electrons in the framework of coupled-cluster theory. We have found excellent agreement between these two approaches for the ground state properties including the permanent dipole moment. We have focused on studies of excited states correlated to the two lowest asymptotes Rb(5p)+Sr(5s) and Rb(5s)+Sr(5s5p) relevant for ongoing experiments on quantum degenerate gases. We present also the Hund c) case potential curves obtained using atomic spin-orbit constants. These potential curves are an excellent starting point for experimental studies of molecular structure of RbSr using high-resolution spectroscopy.

preprint2013arXiv

First-principle interaction potentials for metastable He($^3$S) and Ne($^3$P) with closed-shell molecules. Application to Penning-Ionizing systems

We present new interaction potential curves, calculated from first-principle, for the He$(^3$S)$...$H$_2$ and He$(^3$S)$...$Ar systems, relevant in Penning's ionization experiments. Two different approaches were applied: supermolecular using coupled cluster theory (CC) and perturbational within symmetry-adapted perturbation theory (SAPT). Both methods gave consistent results and the potentials were used to determine the positions of shape resonances in low collision energy scattering regime. We found a good agreement with the most recent scattering experiment of Henson {\em et al.} [Science {\bf 338}, 234 (2012)]. In addition, we investigated two other dimers, composed of metastable Ne and ground state He and Ar atoms. For the Ne$(^3$P)$...$He system a good agreement between CC and SAPT approaches was obtained. The Ne$(^3$P)$...$Ar dimer was described only with SAPT, as CC gave divergent results. Ne$^*$ systems exhibit extremely small electronic orbital angular momentum anisotropy of the potentials. We attribute this effect to the screening of the open 2$p$ shell by the singly occupied 3$s$ shell.

preprint2013arXiv

Scattering lengths in isotopologues of the RbYb system

We model the binding energies of rovibrational levels of the RbYb molecule using experimental data from two-color photoassociation spectroscopy in mixtures of ultracold $^{87}$Rb with various Yb isotopes. The model uses a theoretical potential based on state-of-the-art \emph{ab initio} potentials, further improved by least-squares fitting to the experimental data. We have fixed the number of bound states supported by the potential curve, so that the model is mass scaled, that is, it accurately describes the bound state energies for all measured isotopic combinations. Such a model enables an accurate prediction of the s-wave scattering lengths of all isotopic combinations of the RbYb system. The reduced mass range is broad enough to cover the full scattering lengths range from $-\infty$ to $+\infty$. For example, the $^{87}$Rb$^{174}$Yb system is characterized by a large positive scattering length of $+880(120)$~a.u., while $^{87}$Rb$^{173}$Yb has $a=-626(88)$~a.u.. On the other hand $^{87}$Rb$^{170}$Yb has a very small scattering length of $-14.5(1.8)$~a.u. confirmed by the pair's extremely low thermalization rate. For isotopic combinations including $^{85}$Rb the variation of the interspecies scattering lengths is much smoother ranging from $+39.0(1.6)$ a.u. for $^{85}$Rb$^{176}$Yb to $+230(10)$ a.u. in the case of $^{85}$Rb$^{168}$Yb. Hyperfine corrections to these scattering lengths are also given. We further complement the fitted potential with interaction parameters calculated from alternative methods. The recommended value of the van der Waals coefficient is $C_6$=2837(13)~a.u. and is in agreement and more precise than the current state-of-the-art theoretical predictions (S. G. Porsev, M. S. Safronova, A. Derevianko, and C.W. Clark, arXiv:1307.2654 (2013)).

preprint2011arXiv

Interaction between LiH molecule and Li atom from state-of-the-art electronic structure calculations

State-of-the-art ab initio techniques have been applied to compute the potential energy surface for the lithium atom interacting with the lithium hydride molecule in the Born-Oppenheimer approximation. The interaction potential was obtained using a combination of the explicitly correlated unrestricted coupled-cluster method with single, double, and noniterative triple excitations [UCCSD(T)-F12] for the core-core and core-valence correlation and full configuration interaction for the valence-valence correlation. The potential energy surface has a global minimum 8743 cm^{-1} deep if the Li-H bond length is held fixed at the monomer equilibrium distance or 8825 cm^{-1} deep if it is allowed to vary. In order to evaluate the performance of the conventional CCSD(T) approach, calculations were carried out using correlation-consistent polarized valence X-tuple-zeta bases, with X ranging from 2 to 5. The contribution beyond the CCSD(T)-F12 model, obtained from full configuration interaction (FCI) calculations for the valence-valence correlation, was shown to be very small. At linear LiH-Li geometries the ground-state potential shows an avoided crossing with an ion-pair potential.Using both adiabatic and diabatic pictures we analyse the interaction between the two potential energy surfaces and its possible impact on the collisional dynamics. When the LiH bond is allowed to vary, a seam of conical intersections appears at C_{2v} geometries. At the linear LiH-Li geometry, the conical intersection is at a Li-H distance which is only slightly larger than the monomer equilibrium distance, but for nonlinear geometries it quickly shifts to Li-H distances that are well outside the classical turning points of the ground-state potential of LiH. Finally, the reaction channels for the exchange and insertion reactions are also analyzed.

preprint2011arXiv

Prospects for sympathetic cooling of molecules in electrostatic, ac and microwave traps

We consider how trapped molecules can be sympathetically cooled by ultracold atoms. As a prototypical system, we study LiH molecules co-trapped with ultracold Li atoms. We calculate the elastic and inelastic collision cross sections of LiH + Li with the molecules initially in the ground state and in the first rotationally excited state. We then use these cross sections to simulate sympathetic cooling in a static electric trap, an ac electric trap, and a microwave trap. In the static trap we find that inelastic losses are too great for cooling to be feasible for this system. The ac and microwave traps confine ground-state molecules, and so inelastic losses are suppressed. However, collisions in the ac trap can take molecules from stable trajectories to unstable ones and so sympathetic cooling is accompanied by trap loss. In the microwave trap there are no such losses and sympathetic cooling should be possible.

preprint2011arXiv

The prospects of sympathetic cooling of NH molecules with Li atoms

We calculate the quartet potential energy surface for Li+NH and use it to calculate elastic and spin-relaxation cross sections for collisions in magnetically trappable spin-stretched states. The potential is strongly anisotropic but spin-relaxation collisions are still suppressed by centrifugal barriers when both species are in spin-stretched states. In the ultracold regime, both the elastic and inelastic cross sections fluctuate dramatically as the potential is varied because of Feshbach resonances. The potential-dependence is considerably reduced at higher energies. The major effect of using an unconverged basis set in the scattering calculations is to shift the resonances without changing their general behaviour. We have calculated the ratio of elastic and spin-relaxation cross sections, as a function of collision energy and magnetic field, for a variety of potential energy surfaces. Most of the surfaces produce ratios that are favorable for sympathetic cooling, at temperatures below about 20 mK.

preprint2010arXiv

A DFT Approach to Non-Covalent Interactions via Monomer Polarization and Pauli Blockade

We propose a "DFT+dispersion" treatment which avoids double counting of dispersion terms by deriving the dispersion-free density functional theory (DFT) interaction energy and combining it with DFT-based dispersion. The formalism involves self-consistent polarization of DFT monomers restrained by the exclusion principle via the Pauli blockade technique. Any exchange-correlation potential can be used within monomers, but only the exchange operates between them. The applications to rare-gas dimers, ion-rare gas interactions and hydrogen bonds demonstrate that the interaction energies agree with benchmark values.

preprint2010arXiv

Cold collisions of N atoms and NH molecules in magnetic fields

We calculate the interaction potential between N atoms and NH molecules and use it to investigate cold and ultracold collisions important for sympathetic cooling. The ratio of elastic to inelastic cross sections is large over a wide range of collision energy and magnetic field for most isotopic combinations, so that sympathetic cooling of NH molecules by N atoms is a good prospect. However, there are important effects due to a p-wave shape resonance that may inhibit cooling in some cases. We show that scaling the reduced mass used in the collision is approximately equivalent to scaling the interaction potential. We then explore the dependence of the scattering properties on the reduced mass and explain the resonant effects observed using angular-momentum-insensitive quantum defect theory

preprint2010arXiv

Derivation of the Supermolecular Interaction Energy from the Monomer Densities in the Density Functional Theory

The density functional theory (DFT) interaction energy of a dimer is rigorously derived from the monomer densities. To this end, the supermolecular energy bifunctional is formulated in terms of mutually orthogonal sets of orbitals of the constituent monomers. The orthogonality condition is preserved in the solution of the Kohn-Sham equations through the Pauli blockade method. Numerical implementation of the method provides interaction energies which agree with those obtained from standard supermolecular calculations within less than 0.1% error for three example functionals: Slater-Dirac, PBE0 and B3LYP, and for two model van der Waals dimers: Ne2 and (C2H4)2, and two model H-bond complexes: (HF)2 and (NH3)2.

preprint2009arXiv

Ab initio potential energy surfaces for NH-NH with analytical long range

We present four-dimensional ab initio potential energy surfaces for the three spin states of the NH-NH complex. The potentials are partially based on the work of Dhont et al. [J. Chem. Phys. 123, 184302 (2005)]. The surface for the quintet state is obtained at the RCCSD(T)/aug-cc-pVTZ level of theory and the energy diferences with the singlet and triplet states are calculated at the CASPTn/aug-cc-pVTZ (n = 2; 3) level of theory. The ab initio potentials are fitted to coupled spherical harmonics in the angular coordinates, and the long range is further expanded as a power series in 1/R. The RCCSD(T) potential is corrected for a size-consistency error prior to fitting. The long-range coeficients obtained from the fit are found to be in good agreement with perturbation theory calculations.

preprint2009arXiv

Low-energy collisions of NH3 and ND3 with ultracold Rb atoms

We carry out quantum inelastic scattering calculations of collisions of Rb atoms with inverting NH3 and ND3 molecules in the energy range between 0 and 100 cm-1, which are important for experiments using velocity-controlled molecular beams to probe scattering resonances. We focus on molecules initially in the upper level of the ammonia inversion doublet, which are low-field-seeking and can be controlled in a Stark decelerator. We calculate the integral elastic and state-to-state inelastic cross sections in the coupled states approximation. We demonstrate the presence of both shape and Feshbach resonances in the elastic and inelastic cross sections at low collision energies and discuss their origin in terms of the bound states of Rb--ND3 complex. We also consider elastic and inelastic cross sections in the ultracold regime, using close-coupling calculations, in order to assess the viability of sympathetic cooling of ND3 by Rb. The inelastic cross section for relaxation to the lower level of the inversion doublet is smaller than expected for such a strongly coupled system, but is still likely to be too large to allow sympathetic cooling for ND3 in low-field-seeking states. However, there is a good prospect that sympathetic cooling will be possible for molecules in high-field-seeking states, even when the collision partner is a magnetically trapped atom in a low-field-seeking state.

preprint2009arXiv

Prospects for sympathetic cooling of polar molecules: NH with alkali-metal and alkaline-earth atoms - a new hope

We explore the potential energy surfaces for NH molecules interacting with alkali-metal and alkaline-earth atoms using highly correlated ab-initio electronic structure calculations. The surfaces for interaction with alkali-metal atoms have deep wells dominated by covalent forces. The resulting strong anisotropies will produce strongly inelastic collisions. The surfaces for interaction with alkaline-earth atoms have shallower wells that are dominated by induction and dispersion forces. For Be and Mg the anisotropy is small compared to the rotational constant of NH, so that collisions will be relatively weakly inelastic. Be and Mg are thus promising coolants for sympathetic cooling of NH to the ultracold regime.

preprint2008arXiv

The prospects for producing ultracold NH$_3$ molecules by sympathetic cooling: a survey of interaction potentials

We investigate the possibility of producing ultracold NH3 molecules by sympathetic cooling in a bath of ultracold atoms. We consider the interactions of NH3 with alkali-metal and alkaline-earth atoms, and with Xe, using ab initio coupled-cluster calculations. For Rb-NH3 and Xe-NH3 we develop full potential energy surfaces, while for the other systems we characterize the stationary points (global and local minima and saddle points). We also calculate isotropic and anisotropic Van der Waals C6 coefficients for all the systems. The potential energy surfaces for interaction of NH3 with alkali-metal and alkaline-earth atoms all show deep potential wells and strong anisotropies. The well depths vary from 887 1/cm for Mg-NH3 to 5104 1/cm for Li-NH3. This suggests that all these systems will exhibit strong inelasticity whenever inelastic collisions are energetically allowed and that sympathetic cooling will work only when both the atoms and the molecules are already in their lowest internal states. Xe-NH3 is more weakly bound and less anisotropic.