Energy levels, radiative rates, and electron impact excitation rates for transitions in Li-like ions with 21≤Z≤28

We report calculations of energy levels, radiative rates, and electron impact excitation rates for transitions in Li-like ions with 21≤Z≤28. The General-Purpose Relativistic Atomic Structure Package is adopted for calculating energy levels and radiative rates, while for determining the collision strengths and subsequently the excitation rates, the Dirac Atomic R-matrix Code is used. Oscillator strengths, radiative rates, and line strengths are listed for all E1, E2, M1, and M2 transitions among the lowest 24 levels of the Li-like ions considered. Collision strengths have been averaged over a Maxwellian velocity distribution, and the effective collision strengths obtained are given over a wide temperature range up to 10^7.8 K. Additionally, lifetimes are listed for all calculated levels of these ions. Finally, extensive comparisons are made with results available in the literature, as well as with our analogous calculations for all parameters with the Flexible Atomic Code, in order to assess the accuracy of the results.

Total cross sections for transfer ionization in fast ion-helium collisions

The effects of electron correlation and second-order terms on theoretical total cross sections of transfer ionization in collisions of the helium atom with fast H+, He2+ and Li3+ ions are studied and reported. The total cross sections are calculated using highly correlated wavefunctions with expansion of the transition amplitude in the Born series through the second order. The results of these calculations are in sensible agreement with experimental data.

Entanglement between two superconducting qubits via interaction with nonclassical radiation

We propose a scheme to physically interface superconducting nanocircuits and quantum optics. We address the transfer of quantum information between systems having different physical natures and defined in Hilbert spaces of different dimensions. In particular, we investigate the transfer of the entanglement initially in a nonclassical state of an infinite dimensional system to a pair of superconducting charge qubits. This setup is able to drive an initially separable state of the qubits into an almost pure, highly entangled state suitable for quantum information processing.

A discrete time-dependent method for metastable atoms and molecules in intense fields

The full-dimensional time-dependent Schrodinger equation for the electronic dynamics of single-electron systems in intense external fields is solved directly using a discrete method. Our approach combines the finite-difference and Lagrange mesh methods. The method is applied to calculate the quasienergies and ionization probabilities of atomic and molecular systems in intense static and dynamic electric fields. The gauge invariance and accuracy of the method is established. Applications to multiphoton ionization of positronium, the hydrogen atom and the hydrogen molecular ion are presented. At very high laser intensity, above the saturation threshold, we extend the method using a scaling technique to estimate the quasienergies of metastable states of the hydrogen molecular ion. The results are in good agreement with recent experiments. (C) 2004 American Institute of Physics.

Energy dependence for pair production with electron capture in relativistic heavy-ion collisions

In 'Charge transfer from the negative-energy continuum: alternative mechanism for pair production in relativistic atomic collisions', Eichler (1995 Phys. Rev. Lett. 75 3653) proposes an alternative mechanism for capture by pair production, and from it derives an analytic expression for the total cross section with a surprisingly strong energy dependence. We show that, in fact, there is no alternative mechanism; rather the above mechanism may be more transparently viewed as an ionization-like transition in one centre with inclusion of continuum distortion by the second centre. We further show that to Centre the initial and final states on the target and projectile leads to confusion in the momentum transfer vectors, and hence, respectively that the alleged high-energy behaviour is erroneous.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Concertedness and solvent effects in multiple proton transfer reactions: The formic acid dimer in solution

The issue of multiple proton transfer (PT) reactions in solution is addressed by performing molecular dynamics simulations for a formic acid dimer embedded in a water cluster. The reactant species is treated quantum mechanically, within a density functional approach, while the solvent is represented by a classical model. By constraining different distances within the dimer we analyze the PT process in a variety of situations representative of more complex environments. Free energy profiles are presented, and analyzed in terms of typical solvated configurations extracted from the simulations. A decrease in the PT barrier height upon solvation is rationalized in terms of a transition state which is more polarized than the stable states. The dynamics of the double PT process is studied in a low-barrier case and correlated with solvent polarization fluctuations. Cooperative effects in the motion of the two protons are observed in two different situations: when the solvent polarization does not favor the transfer of one of the two protons and when the motion of the two protons is not synchronized. This body of observations is correlated with local structural and dynamical properties of the solvent in the vicinity of the reactant. (C) 2000 American Institute of Physics. [S0021-9606(00)51121-0].

Molecular conduction: Do time-dependent simulations tell you more than the Landauer approach?

A dynamical method for simulating steady-state conduction in atomic and molecular wires is presented which is both computationally and conceptually simple. The method is tested by calculating the current-voltage spectrum of a simple diatomic molecular junction, for which the static Landauer approach produces multiple steady-state solutions. The dynamical method quantitatively reproduces the static results and provides information on the stability of the different solutions. (c) 2006 American Institute of Physics.

Glory and thresholds effects in H+D-2 reactive angular scattering

We analyse H + D-2 reactive angular scattering using the S- matrix elements obtained by Aoiz et al. and Althorpe et al. Enhancement of small angle scattering in the v' = 3

DFT studies of long-chain FAMEs: theoretical justification for determining chain length and unsaturation from experimental Raman spectra

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of a series of long-chain, saturated fatty acid methyl esters (FAMEs) with the formula CH2nO2 (n = 5-21) and two series of unsaturated FAMEs. The calculations showed that the lowest energy conformer within the saturated FAMEs is the simple (all-trans) structure and, in general, it was possible to reproduce experimental data using calculations on only the all-trans conformer. The only exception was C6H12O2, where a second low-lying conformer had to be included in order to correctly simulate the experimental Raman spectrum. The objective of the work was to provide theoretical justification for the methods that are commonly used to determine the properties of the fats and oils, such as chain length and degree of unsaturation, from experimental Raman data. Here it is shown that the calculations reproduce the trends and calibration curves that are found experimentally and also allow the reasons for the failure of what would appear to be rational measurements to be understood. This work shows that although the assumption that each FAME can simply be treated as a collection of functional groups can be justified in some cases, many of the vibrational modes are complex motions of large sections of the molecules and thus would not be expected to show simple linear trends with changes in structure, such as increasing chain length and/or unsaturation. Simple linear trends obtained from experimental data may thus arise from cancellation of opposing effects, rather than reflecting an underlying simplicity.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

Time delay between singly and doubly ionizing wavepackets in laser-driven helium

We present calculations of the time delay between single and double ionization of helium, obtained from full-dimensionality numerical integrations of the helium-laser Schroedinger equation. The notion of a quantum mechanical time delay is defined in terms of the interval between correlated bursts of single and double ionization. Calculations are performed at 390 and 780 nm in laser intensities that range from 2 X 10^14 to 14 X 10^14 W /cm^2. We find results consistent with the rescattering model of double ionization but supporting its classical interpretation only at 780 nm.