Numerical evaluation of orientation averages and its application to molecular physics

Abstract

In molecular physics, it is often necessary to average over the orientation of molecules when calculating observables, in particular when modelling experiments in the liquid or gas phase. Evaluated in terms of Euler angles, this is closely related to integration over two- or three-dimensional unit spheres, a common problem discussed in numerical analysis. The computational cost of the integration depends significantly on the quadrature method, making the selection of an appropriate method crucial for the feasibility of simulations. After reviewing several classes of spherical quadrature methods in terms of their efficiency and error distribution, we derive guidelines for choosing the best quadrature method for orientation averages and illustrate these with three examples from chiral molecule physics. While Gauss quadratures allow for achieving numerically exact integration for a wide range of applications, other methods offer advantages in specific circumstances. Our guidelines can also by applied to higher-dimensional spherical domains and other geometries. We also present a Python package providing a flexible interface to a variety of quadrature methods.

References in corpus (22)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Measurement of areas on a sphere using Fibonacci and latitude-longitude lattices
• Locally and globally chiral fields for ultimate control of chiral light matter interaction
• Quaternions in molecular modeling
• MOLSCAT: a program for non-reactive quantum scattering calculations on atomic and molecular collisions
• Selective Orientation of Chiral Molecules by Laser Fields with Twisted Polarization
• Observation of enhanced chiral asymmetries in the inner-shell photoionization of uniaxially oriented methyloxirane enantiomers
• Interatomic potentials, electric properties, and spectroscopy of the ground and excited states of the Rb_2 molecule: Ab initio calculations and effect of a non-resonant field
• Theoretical description of circular dichroism in photoelectron angular distributions of randomly oriented chiral molecules after multi-photon photoionization
• Quantum control of photoelectron circular dichroism
• Controlling the sense of molecular rotation: classical vs quantum analysis
• Orienting Asymmetric Molecules by Laser Fields with Skewed Polarization
• Semiclassical catastrophes and accumulative angular squeezing of a kicked quantum rotor
• Photoelectron circular dichroism in the multiphoton ionization by short laser pulses: III. Photoionization of fenchone in different regimes
• Photoelectron circular dichroism of O 1$s$-photoelectrons of uniaxially oriented trifluoromethyloxirane: Energy dependence and sensitivity to molecular configuration
• Photoelectron circular dichroism in the multiphoton ionization by short laser pulses: II. Three- and four-photon ionization of fenchone and camphor
• Detecting handedness of spatially oriented molecules by Coulomb explosion imaging
• Photoelectron circular dichroism with Lissajous-type bichromatic fields: One- vs two-photon ionization of chiral molecules
• Near optimal angular quadratures for polarised radiative transfer
• Perfect control of photoelectron anisotropy for randomly oriented ensembles of molecules by XUV REMPI and polarization shaping
• $SO(3)$ quadratures in angular-momentum projection