Anti-symmetric Compton scattering in LiNiPO4: Towards a direct probe of the magneto-electric multipole moment

Background: Magnetoelectric multipoles, which break both space-inversion and time-reversal symmetries, play an important role in the magnetoelectric response of a material. Motivated by uncovering the underlying fundamental physics of the magnetoelectric multipoles and the possible technological applications of magnetoelectric materials, understanding as well as detecting such magnetoelectric multipoles has become an active area of research in condensed matter physics. Here we employ the well-established Compton scattering effect as a possible probe for the magnetoelectric toroidal moments in LiNiPO4. Methods: We employ combined theoretical and experimental techniques to compute as well as detect the antisymmetric Compton profile in LiNiPO4. For the theoretical investigation we use density functional theory to compute the anti-symmetric part of the Compton profile for the magnetic and structural ground state of LiNiPO4. For the experimental verification, we measure the Compton signals for a single magnetoelectric domain sample of LiNiPO4, and then again for the same sample with its magnetoelectric domain reversed. We then take the difference between these two measured signals to extract the antisymmetric Compton profile in LiNiPO4. Results: Our theoretical calculations indicate an antisymmetric Compton profile in the direction of the ty toroidal moment in momentum space, with the computed antisymmetric profile around four orders of magnitude smaller than the total profile. The difference signal that we measure is consistent with the computed profile, but of the same order of magnitude as the statistical errors and systematic uncertainties of the experiment. Conclusions: While the weak difference signal in the measurements prevents an unambiguous determination of the antisymmetric Compton profile in LiNiPO4, our results motivate further theoretical work to understand the factors that influence the size of the antisymmetric Compton profile, and to identify materials exhibiting larger effects.

A rush to explore protein–ligand electrostatic interaction energy with Charger

The mutual penetration of electron densities between two interacting molecules complicates the computation of an accurate electrostatic interaction energy based on a pseudo-atom representation of electron densities. The numerical exact potential and multipole moment (nEP/MM) method is time-consuming since it performs a 3D integration to obtain the electrostatic energy at short interaction distances. Nguyen et al. [(2018), Acta Cryst. A74, 524–536] recently reported a fully analytical computation of the electrostatic interaction energy (aEP/MM). This method performs much faster than nEP/MM (up to two orders of magnitude) and remains highly accurate. A new program library, Charger, contains an implementation of the aEP/MM method. Charger has been incorporated into the MoProViewer software. Benchmark tests on a series of small molecules containing only C, H, N and O atoms show the efficiency of Charger in terms of execution time and accuracy. Charger is also powerful in a study of electrostatic symbiosis between a protein and a ligand. It determines reliable protein–ligand interaction energies even when both contain S atoms. It easily estimates the individual contribution of every residue to the total protein–ligand electrostatic binding energy. Glutathione transferase (GST) in complex with a benzophenone ligand was studied due to the availability of both structural and thermodynamic data. The resulting analysis highlights not only the residues that stabilize the ligand but also those that hinder ligand binding from an electrostatic point of view. This offers new perspectives in the search for mutations to improve the interaction between the two partners. A proposed mutation would improve ligand binding to GST by removing an electrostatic obstacle, rather than by the traditional increase in the number of favourable contacts.

Spherical symmetry breaking in electric, magnetic and toroidal multipole moment radiations in spherical toroidal resonant cavities and optimum-efficiency antennas

This Letter reports the breaking of the spherical symmetry in the complete electromagnetic multipole expansion when its sources are distributed on spherical toroidal surfaces, identifying the specic geometrical and physical changes fromthe familiar case of sources on a spherical surface. In fact, for spherical toroids dened by concentric spherical rings and symmetric conical rings, the boundary conditions at the latter are not compatible in general with integer values for the orbital angular momentum label of the multipole moments: the polar angle eigenfunctions become Legendre functions of order λ and associativity m represented as innite series with a denite parity, and their complementary associated radial functions are spherical Bessel functions of the same order λ. Consequently, the corresponding multipole sources for the electric, magnetic and toroidal moments and their connections are identied within the Debye formalism, and theappropriate outgoing wave Green functions are constructed in the new basis of eigenfunctions of the Helmholtz equation. Our familiarity with the exact solutions, for the cases of the complete sphere and of cylindrical toroids, allow us to give a preliminary account of the electromagnetic elds for the spherical toroids via the integration of their sources and the Green function for resonant cavities and optimum effciency antennas.

Bridge-Mediated RET between Two Chiral Molecules

Molecular quantum electrodynamics (QED) theory is employed to calculate the rate of resonance energy transfer (RET) between a donor, D, described by an electric dipole and quadrupole, and magnetic dipole coupling, and an identical acceptor molecule, A, that is mediated by a third body, T, which is otherwise inert. A single virtual photon propagates between D and T, and between T and A. Time-dependent perturbation theory is used to compute the matrix element, from which the transfer rate is evaluated using the Fermi golden rule. This extends previous studies that were limited to the electric dipole approximation only and admits the possibility of the exchange of excitation between a chiral emitter and absorber. Rate terms are computed for specific pure and mixed multipole-dependent contributions of D and A for both an oriented arrangement of the three particles and for the freely tumbling situation. Mixed multipole moment contributions, such as those involving electric–magnetic dipole or electric dipole–quadrupole coupling at one center, do not survive random orientational averaging. Interestingly, the mixed electric–magnetic dipole D and A rate term is non-vanishing and discriminatory, exhibiting a dependence on the chirality of the emitter and absorber, and is entirely retarded. It vanishes, however, if D and A are oriented perpendicularly to one another. Near- and far-zone asymptotes of isotropic contributions to the rate are also evaluated, demonstrating radiationless short-range transfer and inverse-square radiative exchange at very large separations.

Multipole-moment effects in ion–molecule reactions at low temperatures: part I – ion-dipole enhancement of the rate coefficients of the He+ + NH3 and He+ + ND3 reactions at collisional energies Ecoll/kB near 0 K

The energy dependence of the rates of the reactions between He+ and ammonia (NY3, Y = {H,D}), forming NY2+, Y and He as well as NY+, Y2 and He has been measured at low collision energies near 0 K.

Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. III. Application to crystal structures via the Ewald and direct summation methods

The previously reported exact potential and multipole moment (EP/MM) method for fast and accurate evaluation of the intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density [Volkov, Koritsanszky & Coppens (2004). Chem. Phys. Lett. 391, 170–175; Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524–536; Nguyen & Volkov (2019). Acta Cryst. A75, 448–464] is extended to the calculation of electrostatic interaction energies in molecular crystals using two newly developed implementations: (i) the Ewald summation (ES), which includes interactions up to the hexadecapolar level and the EP correction to account for short-range electron-density penetration effects, and (ii) the enhanced EP/MM-based direct summation (DS), which at sufficiently large intermolecular separations replaces the atomic multipole moment approximation to the electrostatic energy with that based on the molecular multipole moments. As in the previous study [Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524–536], the EP electron repulsion integral is evaluated analytically using the Löwdin α-function approach. The resulting techniques, incorporated in the XDPROP module of the software package XD2016, have been tested on several small-molecule crystal systems (benzene, L-dopa, paracetamol, amino acids etc.) and the crystal structure of a 181-atom decapeptide molecule (Z = 4) using electron densities constructed via the University at Buffalo Aspherical Pseudoatom Databank [Volkov, Li, Koritsanszky & Coppens (2004). J. Phys. Chem. A, 108, 4283–4300]. Using a 2015 2.8 GHz Intel Xeon E3-1505M v5 computer processor, a 64-bit implementation of the Löwdin α-function and one of the higher optimization levels in the GNU Fortran compiler, the ES method evaluates the electrostatic interaction energy with a numerical precision of at least 10−5 kJ mol−1 in under 6 s for any of the tested small-molecule crystal structures, and in 48.5 s for the decapeptide structure. The DS approach is competitive in terms of precision and speed with the ES technique only for crystal structures of small molecules that do not carry a large molecular dipole moment. The electron-density penetration effects, correctly accounted for by the two described methods, contribute 28–64% to the total electrostatic interaction energy in the examined systems, and thus cannot be neglected.