Theoretical study of the finite temperature spectroscopy in van der Waals clusters. III Solvated Chromophore as an effective diatomics

The absorption spectroscopy of calcium-doped argon clusters is described in terms of an effective diatomics molecule Ca-(Ar_n), in the framework of semiclassical vertical transitions. We show how, upon choosing a suitable reaction coordinate, the effective finite-temperature equilibrium properties can be obtained for the ground- and excited-surfaces from the potential of mean force (PMF). An extension of the recent multiple range random-walk method is used to calculate the PMF over continuous intervals of distances. The absorption spectra calculated using this single-coordinate description are found to be in good agreement with the spectra obtained from high-statistics Monte Carlo data, in various situations. For CaAr$_{13}$, we compare the performances of two different choices of the reaction coordinate. For CaAr_37, the method is seen to be accurate enough to distinguish between different low-energy structures. Finally, the idea of casting the initial many-body problem into a single degree of freedom problem is tested on the spectroscopy of calcium in bulk solid argon.Comment: 8 pages, 9 figure

Theoretical study of finite temperature spectroscopy in van der Waals clusters. II Time-dependent absorption spectra

Using approximate partition functions and a master equation approach, we investigate the statistical relaxation toward equilibrium in selected CaAr$_n$ clusters. The Gaussian theory of absorption (previous article) is employed to calculate the average photoabsorption intensity associated with the 4s^2-> 4s^14p^1 transition of calcium as a function of time during relaxation. In CaAr_6 and CaAr_10 simple relaxation is observed with a single time scale. CaAr_13 exhibits much slower dynamics and the relaxation occurs over two distinct time scales. CaAr_37 shows much slower relaxation with multiple transients, reminiscent of glassy behavior due to competition between different low-energy structures. We interpret these results in terms of the underlying potential energy surfaces for these clusters.Comment: 10 pages, 9 figure

Mechanical, Electrical, and Magnetic Properties of Ni Nanocontacts

The dynamic deformation upon stretching of Ni nanowires as those formed with mechanically controllable break junctions or with a scanning tunneling microscope is studied both experimentally and theoretically. Molecular dynamics simulations of the breaking process are performed. In addition, and in order to compare with experiments, we also compute the transport properties in the last stages before failure using the first-principles implementation of Landauer's formalism included in our transport package ALACANT.Comment: 5 pages, 6 figure

Formation of a Metallic Contact: Jump to Contact Revisited

The transition from tunneling to metallic contact between two surfaces does not always involve a jump, but can be smooth. We have observed that the configuration and material composition of the electrodes before contact largely determines the presence or absence of a jump. Moreover, when jumps are found preferential values of conductance have been identified. Through combination of experiments, molecular dynamics, and first-principles transport calculations these conductance values are identified with atomic contacts of either monomers, dimers or double-bond contacts.Comment: 4 pages, 5 figure

Stacked clusters of polycyclic aromatic hydrocarbon molecules

Clusters of polycyclic aromatic hydrocarbon (PAH) molecules are modelled using explicit all-atom potentials using a rigid body approximation. The PAH's considered range from pyrene (C10H8) to circumcoronene (C54H18), and clusters containing between 2 and 32 molecules are investigated. In addition to the usual repulsion-dispersion interactions, electrostatic point-charge interactions are incorporated, as obtained from density functional theory calculations. The general electrostatic distribution in neutral or singly charged PAH's is reproduced well using a fluctuating charges analysis, which provides an adequate description of the multipolar distribution. Global optimization is performed using a variety of methods, including basin-hopping and parallel tempering Monte Carlo. We find evidence that stacking the PAH molecules generally yields the most stable motif. A structural transition between one-dimensional stacks and three-dimensional shapes built from mutiple stacks is observed at larger sizes, and the threshold for this transition increases with the size of the monomer. Larger aggregates seem to evolve toward the packing observed for benzene in bulk.Difficulties met in optimizing these clusters are analysed in terms of the strong anisotropy of the molecules. We also discuss segregation in heterogeneous clusters and vibrational properties in the context of astrophysical observations.Comment: 12 pages, 7 figure

Electrostatic potential variations on stellarator magnetic surfaces in low collisionality regimes

The component of the neoclassical electrostatic potential that is non-constant on the magnetic surface, that we denote by $\tilde\varphi$, can affect radial transport of highly charged impurities, and this has motivated its inclusion in some modern neoclassical codes. The number of neoclassical simulations in which $\tilde\varphi$ is calculated is still scarce, partly because they are usually demanding in terms of computational resources, especially at low collisionality. In this paper the size, the scaling with collisionality and with aspect ratio, and the structure of $\tilde\varphi$ on the magnetic surface are analytically derived in the $1/\nu$, $\sqrt{\nu}$ and superbanana-plateau regimes of stellarators close to omnigeneity; i. e. stellarators that have been optimized for neoclassical transport. It is found that the largest $\tilde\varphi$ that the neoclassical equations admit scales linearly with the inverse aspect ratio and with the size of the deviation from omnigeneity. Using a model for a perturbed omnigeneous configuration, the analytical results are verified and illustrated with calculations by the code KNOSOS. The techniques, results and numerical tools employed in this paper can be applied to neoclassical transport problems in tokamaks with broken axisymmetry.Comment: 30 pages, 12 figures, 1 table. Published versio

Numerical Implementation of Gradient Algorithms

A numerical method for computational implementation of gradient dynamical systems is presented. The method is based upon the development of geometric integration numerical methods, which aim at preserving the dynamical properties of the original ordinary differential equation under discretization. In particular, the proposed method belongs to the class of discrete gradients methods, which substitute the gradient of the continuous equation with a discrete gradient, leading to a map that possesses the same Lyapunov function of the dynamical system, thus preserving the qualitative properties regardless of the step size. In this work, we apply a discrete gradient method to the implementation of Hopfield neural networks. Contrary to most geometric integration methods, the proposed algorithm can be rewritten in explicit form, which considerably improves its performance and stability. Simulation results show that the preservation of the Lyapunov function leads to an improved performance, compared to the conventional discretization.Spanish Government project no. TIN2010-16556 Junta de Andalucía project no. P08-TIC-04026 Agencia Española de Cooperación Internacional para el Desarrollo project no. A2/038418/1

Theoretical study of finite temperature spectroscopy in van der Waals clusters. I. Probing phase changes in CaAr_n

The photoabsorption spectra of calcium-doped argon clusters CaAr_n are investigated at thermal equilibrium using a variety of theoretical and numerical tools. The influence of temperature on the absorption spectra is estimated using the quantum superposition method for a variety of cluster sizes in the range 6<=n<=146. At the harmonic level of approximation, the absorption intensity is calculated through an extension of the Gaussian theory by Wadi and Pollak [J. Chem. Phys. vol 110, 11890 (1999)]. This theory is tested on simple, few-atom systems in both the classical and quantum regimes for which highly accurate Monte Carlo data can be obtained. By incorporating quantum anharmonic corrections to the partition functions and respective weights of the isomers, we show that the superposition method can correctly describe the finite-temperature spectroscopic properties of CaAr_n systems. The use of the absorption spectrum as a possible probe of isomerization or phase changes in the argon cluster is discussed at the light of finite-size effects.Comment: 17 pages, 9 figure