Computational challenges in atomic, molecular and optical physics

Six challenges are discussed. These are the laser-driven helium atom; the laser-driven hydrogen molecule and hydrogen molecular ion: electron scattering (with ionization) from one-electron atoms; the vibrational and rotational structure of molecules such as H-3(+) and water at their dissociation limits; laser- heated clusters; and quantum degeneracy and Bose-Einstein condensation. The first four concern fundamental few-body systems where use of high-performance computing (HPC) is currently making possible accurate modelling from first principles. This leads to reliable predictions and support for laboratory experiment as well as true understanding of the dynamics. Important aspects of these challenges addressable only via a terascale facility are set out. Such a facility makes the last two challenges in the above list meaningfully accessible for the first time, and the scientific interest together with the prospective role for HPC in these is emphasized.

Relation between effective atomic charge and chemical shifts in X-ray absorption spectra of transition-metal compounds

Chemical shifts of K absorption discontinuities, Delta E, of several manganese, iron and cobalt oxides with the metal in the formal oxidation states between +2 and +4, have been measured. These data, together with data in the literature on other compounds of these metals, can be fitted into the expression Delta E=aq+bq2, where q is the effective atomic charge on the metal. Theoretical considerations also support this functional relationship between Delta E and q.

PAMOP: Petascale atomic, molecular and optical collision calculations

Petaflop architectures are currently being utilized efficiently to perform large scale computations in Atomic, Molecular and Optical Collisions. We solve the Schr\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach. We briefly outline the parallel methodology used and implemented for the current suite of Breit-Pauli and DARC codes. In this report, various examples are shown from our theoretical results compared with experimental results obtained from Synchrotron Radiation facilities where the Cray architecture at HLRS is playing an integral part in our computational projects.

Intra- vs intermolecular photoinduced electron transfer reactions of a macrocyclic donor-acceptor dyad

The synthesis, structural characterization, and photophysical behavior of a 14-membered tetraazamacrocycle with pendant 4-dimethylaminobenzyl (DMAB) and 9-anthracenylmethyl groups is reported (L-3, 6-((9-anthracenylmethyl)amino)-trans-6,13-dimethyl-13-((4-dimethylaminobenzyl)amino)-1,4,8,11-tetraaza-cyclotetradecane). In its free base form, this compound displays rapid intramolecular photoinduced electron transfer (PET) quenching of the anthracene emission, with both the secondary amines and the DMAB group capable of acting as electron donors. When complexed with Zn(II), the characteristic fluorescence of the anthracene chromophore is restored as the former of these pathways is deactivated by coordination. Importantly, it is shown that the DMAB group, which remains uncoordinated and PET active, acts only very weakly to quench emission, by comparison to the behavior of a model Zn complex lacking the pendant DMAB group, [ZnL2](2+) (Chart 1). By contrast, Stern-Volmer analysis of intermolecular quenching of [ZnL2](2+) by N,N-dimethylaniline (DMA) has shown that this reaction is diffusion limited. Hence, the pivotal role of the bridge in influencing intramolecular PET is highlighted.

Classical and quantum dynamics of a model for atomic-molecular Bose-Einstein condensates

We study a model for a two-mode atomic-molecular Bose-Einstein condensate. Starting with a classical analysis we determine the phase space fixed points of the system. It is found that bifurcations of the fixed points naturally separate the coupling parameter space into four regions. The different regions give rise to qualitatively different dynamics. We then show that this classification holds true for the quantum dynamics.

Electronic structure and chemical bonding in W(2) molecule

The electronic structure of the lowest-lying electronic states of W(2) were investigated at the CASPT2 level. The ground state is a X(1)Sigma(+)(g) state, followed by the a(3)Delta(u), b(3)Sigma(+)(u) and A(1)Delta(u) electronic states. Seven low-lying Omega-states were computed: (1)0(g)(+), (2)3(u), (3)2(u), (4)1(u), (5)0(u)(-), (6)1(u), and (7)2(u), with the ground state corresponding to the (1)0(g)(+)(X(1)Sigma(+)(g)) state. Comparison with the other VIB transition metal group dimers indicates a common pattern of electronic structure and spectroscopic properties. (C) 2010 Elsevier B.V. All rights reserved.

Electronic Structure and Chemical Bonding in the Lowest Electronic States of TcN

Multiconfiguration second-order perturbation theory, with the inclusion of relativistic effects and spin-orbit Coupling, was employed to investigate the nature of the ground and low-lying Lambda-S and Omega states of the TcN molecule. Spectroscopic constants, effective bond order, and potential energy curves for 13 low-lying Lambda-S states and 5 Omega states are given, The computed ground state of TcN is of Omega = 3 symmetry (R(e) = 1.605 angstrom and omega(e) = 1085 cm(-1)), originating mainly from the (3)Delta Lambda-S ground state. This result is contrasted with the nature of the ground state for other VIIB transtion-metal mononitrides, including X(3)Sigma(-) symmetry for MnN and Omega = 0(+) symmetry for ReN, derived also from a X(3)Sigma(-) state.

Electronic structure and chemical bonding in the ground states of Tc(2) and Re(2)

Multiconfiguration second-order perturbation theory, including relativistic effects and spin-orbit coupling, has been employed to investigate the nature of the chemical bonding in the ground state of Tc(2) and Re(2). The Tc(2) ground state is found to be a 0(g)(+) state, with an effective bond order (EBO) of 4.4, and a dissociation energy of 3.25 eV. The Re(2) ground state is a 1(g) state, with EBO = 4.3. Almost degenerate to it, is a 0(g)(+) state (T(e) = 77 cm(-1)), with EBO = 4.1. Experimental evidence also indicates that the ground state is of 1(g) nature. The dissociation energy is computed to be 5.0 eV in agreement with an experimental estimate of 4 +/- 1 eV.

Electronic Structure and Chemical Bonding in the Ground and Low-Lying Electronic States of Ta(2)

The electronic structure and chemical bonding of the ground and low-lying Lambda - S and Omega states of Ta(2) were investigated at the multiconfiguration second-order perturbation theory (CASSCF//CASPT2) level. The ground state of Ta(2) is computed to be a X(3)Sigma(-)(g) state (R(e) = 2.120 angstrom, omega(e) = 323 cm(-1), and D(e) = 4.65 eV), with two low-lying singlet states close to it (a(1) Sigma(+)(g) : T(e) = 409 cm(-1), R(e) = 2.131 angstrom, and omega(e) = 313 cm(-1); b(1) Gamma(g): T(e) = 1, 038 cm(-1), R(e) = 2.127 angstrom, and omega(e) = 316 cm(-1)). These electronic states are derived from the same electronic configuration: vertical bar 13 sigma(2)(g)14 sigma(2)(g)7 delta(2)(g)13 pi(4)(u)>. The effective bond order of the X(3) Sigma(-)(g) state is 4.52, which indicates that the Ta atoms are bound by a quintuple chemical bond. The a(1) Sigma(+)(g) state interacts strongly with the X(3)Sigma(-)(g) g ground state by a second-order spin-orbit interaction, giving rise to the (1)0(g)(+) (ground state) (dominated by the X(3)Sigma(-)(g) Lambda - S ground state) and (9)0(g)(+) (dominated by the a(1) Sigma(+)(g) Lambda - S state) Omega states. These results are in line with those reported for the group 5B homonuclear transition metal diatomics. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1306-1315, 2011

Molecular and Electronic Structure of the Peptide Subunit of Geobacter sulfurreducens Conductive Phi from First Principles

The respiration of metal oxides by the bacterium Geobacter sulfurreducens requires the assembly of a small peptide (the GS pilin) into conductive filaments termed pili. We gained insights into the contribution of the GS pilin to the pilus conductivity by developing a homology model and performing molecular dynamics simulations of the pilin peptide in vacuo and in solution. The results were consistent with a predominantly helical peptide containing the conserved a-helix region required for pilin assembly but carrying a short carboxy-terminal random-coiled segment rather than the large globular head of other bacterial pilins. The electronic structure of the pain was also explored from first principles and revealed a biphasic charge distribution along the pilin and a low electronic HOMO-LUMO gap, even in a wet environment. The low electronic band gap was the result of strong electrostatic fields generated by the alignment of the peptide bond dipoles in the pilin's alpha-helix and by charges from ions in solution and amino acids in the protein. The electronic structure also revealed some level of orbital delocalization in regions of the pilin containing aromatic amino acids and in spatial regions of high resonance where the HOMO and LUMO states are, which could provide an optimal environment for the hopping of electrons under thermal fluctuations. Hence, the structural and electronic features of the pilin revealed in these studies support the notion of a pilin peptide environment optimized for electron conduction.

Three-dimensional solitons in coupled atomic-molecular Bose-Einstein condensates

We present a theoretical analysis of three-dimensional (3D) matter-wave solitons and their stability properties in coupled atomic and molecular Bose-Einstein condensates (BECs). The soliton solutions to the mean-field equations are obtained in an approximate analytical form by means of a variational approach. We investigate soliton stability within the parameter space described by the atom-molecule conversion coupling, the atom-atom s-wave scattering, and the bare formation energy of the molecular species. In terms of ordinary optics, this is analogous to the process of sub- or second-harmonic generation in a quadratic nonlinear medium modified by a cubic nonlinearity, together with a phase mismatch term between the fields. While the possibility of formation of multidimensional spatiotemporal solitons in pure quadratic media has been theoretically demonstrated previously, here we extend this prediction to matter-wave interactions in BEC systems where higher-order nonlinear processes due to interparticle collisions are unavoidable and may not be neglected. The stability of the solitons predicted for repulsive atom-atom interactions is investigated by direct numerical simulations of the equations of motion in a full 3D lattice. Our analysis also leads to a possible technique for demonstrating the ground state of the Schrodinger-Newton and related equations that describe Bose-Einstein condensates with nonlocal interparticle forces.

Ionization of atomic hydrogen and He+ by slow antiprotons

We study the ionization of H(1s), He+(1s) and He+(2s) by antiprotons in the energy range from 0.1 to 500 keV. We adopt a semiclassical single centre close-coupling approach in which the wavefunction for the electron is expanded in a B-spline basis centred on the nucleus of the atom/ion. Comparison is made with existing theoretical calculations and available experimental data. The results are encouraging.