Electric potential invariants and ions-in-molecules effective potentials for molecular Rydberg states.

The dependence of multipole moments and polarizabilities on external fields appears in many applications including biomolecular molecular mechanics, optical non-linearity, nanomaterial calculations, and the perturbation of spectroscopic signatures in atomic clocks. Over a wide range of distances, distributed multipole and polarizability potentials can be applied to obtain the variation of atom-centered atoms-in-molecules electric properties like bonding-quenched polarizability. For cylindrically symmetric charge distributions, we examine single-center and atom-centered effective polarization potentials in a non-relativistic approximation for Rydberg states. For ions, the multipole expansion is strongly origin-dependent, but we note that origin-independent invariants can be defined. The several families of invariants correspond to optimized representations differing by origin and number of terms. Among them, a representation at the center of dipole polarizability optimizes the accuracy of the potential with terms through 1/r4. We formulate the single-center expansion in terms of polarization-modified effective multipole moments, defining a form related to the source-multipole expansion of Brink and Satchler. Atom-centered potentials are an origin independent alternative but are limited both by the properties allowed at each center and by the neglected effects like bond polarizability and charge flow. To enable comparisons between single-center effective potentials in Cartesian or spherical form and two-center effective potentials with differing levels of mutual induction between atomic centers, we give analytical expressions for the bond-length and origin-dependence of multipole and polarizability terms projected in the multipole and polarizability expansion of Buckingham. The atom-centered potentials can then be used with experimental data and ab initio calculations to estimate atoms-in-molecules properties. Some results are given for BaF+ and HF showing the utility and limitations of the approach. More detailed results on X 1Σ+ CaF+ are published separately.